*HEADING
: two shell layers, element[S4R], multiple *SUBMODEL cards
*NODE
1,    0.
4,   10.
5,    0.,3.
8,   10.,3.
13,    6.666667
14,   10.
17,    6.666667,3.
18,   10.,3.
*NGEN,NSET=NTOP
1,4,1
5,8,1
*NSET,NSET=NBOT
13,14,17,18
*NSET,NSET=BC
1,5
*NSET,NSET=NALL
NTOP,NBOT
*NSET,NSET=LOADUP
 4,8
*NSET,NSET=LOADDOWN
 14,18
*ELEMENT,TYPE=S4R,ELSET=EALL
1,1,2,6,5
4,2,13,17,6
5,13,14,18,17
*ELGEN,ELSET=EALL
1,3,1,1
*ELSET,ELSET=ECOMMON
1,
*ELSET,ELSET=ETOP
2,3
*ELSET,ELSET=EBOT
4,5
*ELSET,ELSET=EALL
ETOP,EBOT
*SHELL SECTION,ELSET=ECOMMON,MATERIAL=MAT1
0.001, 
*SHELL SECTION,ELSET=ETOP,MATERIAL=MAT1
0.001, 
*SHELL SECTION,ELSET=EBOT,MATERIAL=MAT1
0.001, 
*MATERIAL,NAME=MAT1
*ELASTIC
1.E6,0.3
*BOUNDARY
BC,ENCASTRE
*STEP
*STATIC
0.5,1.
*BOUNDARY
 LOADUP,3,3,1.0
 LOADDOWN,3,3,-1.0
*EL FILE,FREQ=999
S,E,
*NODE FILE,FREQ=1
U,
*END STEP