To be easier for hand modifications, I used very course mesh (two quad4 elements per width).
Code: Select all
*Node, NSET=Nall
1, 6000, 0, 2500
2, 6000, 0, 0
3, 6000, 2000, 2500
4, 6000, 2000, 0
5, 0, 0, 2500
6, 0, 0, 0
7, 0, 2000, 2500
8, 0, 2000, 0
9, 6000, 0, 833.3333333333
10, 6000, 0, 1666.666666667
11, 6000, 1000, 2500
12, 6000, 2000, 833.3333333333
13, 6000, 2000, 1666.666666667
14, 6000, 1000, 0
15, 0, 0, 833.3333333333
16, 0, 0, 1666.666666667
17, 0, 1000, 2500
18, 0, 2000, 833.3333333333
19, 0, 2000, 1666.666666667
20, 0, 1000, 0
21, 1000, 2000, 2500
22, 2000, 2000, 2500
23, 3000, 2000, 2500
24, 4000, 2000, 2500
25, 5000, 2000, 2500
26, 1000, 0, 2500
27, 2000, 0, 2500
28, 3000, 0, 2500
29, 4000, 0, 2500
30, 5000, 0, 2500
31, 6000, 994.5004138628, 831.4934419314
32, 6000, 996.8858840134, 1665.089340799
33, 0, 994.5004138628, 831.4934419314
34, 0, 996.8858840134, 1665.089340799
35, 4999.983071052, 999.9830710519, 2500
36, 996.0531366691, 996.0531366691, 2500
37, 1998.547335024, 998.547335024, 2500
38, 2999.640259853, 999.6402598526, 2500
39, 3999.92550759, 999.9255075899, 2500
** DUPLICATED NODES
** 1 -> 40, 3 -> 41, 5 -> 42, 7 -> 43, 11 -> 44, 17 -> 45
40, 6000, 0, 2500
41, 6000, 2000, 2500
42, 0, 0, 2500
43, 0, 2000, 2500
44, 6000, 1000, 2500
45, 0, 1000, 2500
** Face elements
*Element, TYPE=S4, ELSET=Efaces
33, 31, 32, 10, 9
34, 32, 11, 1, 10
35, 32, 31, 12, 13
36, 11, 32, 13, 3
37, 2, 14, 31, 9
38, 12, 31, 14, 4
39, 33, 15, 16, 34
40, 34, 16, 5, 17
41, 34, 19, 18, 33
42, 17, 7, 19, 34
43, 6, 15, 33, 20
44, 18, 8, 20, 33
**45, 11, 3, 25, 35 REWRITTEN ELEMENTS
45, 44, 41, 25, 35
**46, 30, 1, 11, 35
46, 30, 40, 44, 35
47, 23, 22, 37, 38
48, 35, 25, 24, 39
49, 24, 23, 38, 39
50, 37, 27, 28, 38
51, 29, 30, 35, 39
52, 38, 28, 29, 39
**53, 17, 36, 21, 7
53, 45, 36, 21, 43
54, 37, 36, 26, 27
55, 21, 36, 37, 22
**56, 36, 17, 5, 26
56, 36, 45, 42, 26
** BINDING DOF OF DUPLICATED NODES
** NODES 1 AND 40
*EQUATION
2
1, 1, 1, 40, 1, -1
2
1, 2, 1, 40, 2, -1
2
1, 3, 1, 40, 3, -1
** NODES 3 AND 41
2
3, 1, 1, 41, 1, -1
2
3, 2, 1, 41, 2, -1
2
3, 3, 1, 41, 3, -1
** NODES 7 AND 43
2
7, 1, 1, 42, 1, -1
2
7, 2, 1, 42, 2, -1
2
7, 3, 1, 42, 3, -1
** NODES 11 AND 44
2
11, 1, 1, 44, 1, -1
2
11, 2, 1, 44, 2, -1
2
11, 3, 1, 44, 3, -1
** NODES 17 AND 45
2
17, 1, 1, 45, 1, -1
2
17, 2, 1, 45, 2, -1
2
17, 3, 1, 45, 3, -1
I tried also variant with *MPC cards. It triggered nonlinear calculation and gave me no results, spooles.out file saying: "matrix found to be singular". Just instead of *EQUATIONS, there were: