*HEADING
o DSA ANALYSIS FOR SHELLS
o TOTAL DSA WITH TWO STEPS
*PARAMETER
PHI=45.0
T=0.5
E=2.e05
X=0.0
*DESIGN PARAMETER
PHI,T,E,X
*NODE
** S4R
 1,  0.0, 190.0, 0.0
 2,  0.0, 200.0, 0.0
 3, 10.0, 190.0, 0.0
 4, 10.0, 200.0, 0.0
 5, 20.0, 190.0, 0.0
 6, 20.0, 200.0, 0.0
** S4
13,  0.0, 150.0, 0.0
14,  0.0, 160.0, 0.0
15, 10.0, 150.0, 0.0
16, 10.0, 160.0, 0.0
17, 20.0, 150.0, 0.0
18, 20.0, 160.0, 0.0
** S3R
19,  0.0, 130.0, 0.0
20,  0.0, 140.0, 0.0
21, 10.0, 130.0, 0.0
22, 10.0, 140.0, 0.0
23, 20.0, 130.0, 0.0
24, 20.0, 140.0, 0.0
** STRI3
25,  0.0, 110.0, 0.0
26,  0.0, 120.0, 0.0
27, 10.0, 110.0, 0.0
28, 10.0, 120.0, 0.0
29, 20.0, 110.0, 0.0
30, 20.0, 120.0, 0.0
** S8R
31,  0.0,  90.0, 0.0
32,  0.0,  95.0, 0.0
33,  0.0, 100.0, 0.0
34,  5.0,  90.0, 0.0
35,  5.0, 100.0, 0.0
36, 10.0,  90.0, 0.0
37, 10.0,  95.0, 0.0
38, 10.0, 100.0, 0.0
39, 15.0,  90.0, 0.0
40, 15.0, 100.0, 0.0
41, 20.0,  90.0, 0.0
42, 20.0,  95.0, 0.0
43, 20.0, 100.0, 0.0
*PARAMETER SHAPE VARIATION, PARAMETER=X
3, 1.0, 0.0, 0.0
4, 1.0, 0.0, 0.0
5, 2.0, 0.0, 0.0
6, 2.0, 0.0, 0.0
15,1.0, 0.0, 0.0
16,1.0, 0.0, 0.0
17,2.0, 0.0, 0.0
18,2.0, 0.0, 0.0
21,1.0, 0.0, 0.0
22,1.0, 0.0, 0.0
23,2.0, 0.0, 0.0
24,2.0, 0.0, 0.0
27,1.0, 0.0, 0.0
28,1.0, 0.0, 0.0
29,2.0, 0.0, 0.0
30,2.0, 0.0, 0.0
36,1.0, 0.0, 0.0
37,1.0, 0.0, 0.0
38,1.0, 0.0, 0.0
39,1.5, 0.0, 0.0
40,1.5, 0.0, 0.0
41,2.0, 0.0, 0.0
42,2.0, 0.0, 0.0
43,2.0, 0.0, 0.0
*ELEMENT, TYPE=S4R, ELSET=S4R
1, 1, 3, 4, 2
2, 3, 5, 6, 4
*ELEMENT, TYPE=S4, ELSET=S4
5, 13, 15, 16, 14
6, 15, 17, 18, 16
*ELEMENT, TYPE=S3R, ELSET=S3R
7, 19, 22, 20
8, 19, 21, 22
9, 21, 24, 22
10,21, 23, 24
*ELEMENT, TYPE=STRI3, ELSET=STRI3
11,25, 28, 26
12,25, 27, 28
13,27, 30, 28
14,27, 29, 30
*ELEMENT, TYPE=S8R, ELSET=S8R
15, 31, 36, 38, 33, 34, 37, 35, 32
16, 36, 41, 43, 38, 39, 42, 40, 37
*ELSET, ELSET=EALL
S4R,S4,S3R,STRI3,S8R
*SHELL SECTION, COMPOSITE, ELSET=EALL
0.5, 3, MAT1, <PHI>
<T>, 3, MAT2, 90.0
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ENGINEERING CONSTANTS
<E>, 1000., 1000., 0.0, 0.0, 0.0, 1e05, 500.
500.,
*DENSITY
0.001,
*MATERIAL, NAME=MAT2
*ELASTIC
2E05, 0.3
*DENSITY
0.01,
*NSET, NSET=FIXED
1, 2,   13, 14, 19, 20, 25, 26,
31, 32, 33
*NSET, NSET=LOAD
5,6,17,18,23,24,29,30,
41,43
*BOUNDARY
FIXED, ENCASTRE
*DSA CONTROLS, FORMULATION=TOTAL
*STEP,NLGEOM,DSA
*STATIC
0.3, 1.0
*CLOAD
LOAD, 3, 20.0
*DESIGN RESPONSE
*NODE RESPONSE
U,
RF,
*ELEMENT RESPONSE
S,
LE,
*END STEP
*STEP,DSA
*FREQUENCY,EIGENSOLVER=LANCZOS,SIM=NO
10,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*DESIGN RESPONSE
*END STEP