*HEADING
 SHELL TO SOLID SUBMODEL                       
*PREPRINT,ECHO=NO
*RESTART,WRITE
*SYSTEM
0.0,0.0,0.0, 1.0, 0.0, 0.0
0.0, 1.0, 0.466308 
*NODE
10000, 0.0, 24.0, 0.0
  1, 0.0, 15.9882, 0.0
 21, 0.0, 32.0118, 0.0
2647, 0.0, 16.3976, 0.0
2667, 0.0, 31.6024, 0.0
10001,
*NGEN,NSET=A1,LINE=C
1,21,1,10000,,,,0.0,0.0,1.0
*NGEN,NSET=A2,LINE=C
2647,2667,1,10000,,,,0.0,0.0,1.0
*SYSTEM
*NFILL,NSET=SNODE_1
A1,A2,6,441
*NCOPY,CHANGE NUMBER=21,OLD SET=SNODE_1,NEW SET=ELBOW,SHIFT,MULTIPLE=20
0.0,0.0,0.0
0.0,0.0,0.0, 1.0,0.0,0.0, 2.0
*NSET,NSET=YZ_PLANE,GEN
22,400,21
463,841,21
904,1282,21
1345,1723,21
1786,2164,21
2227,2605,21
2668,3046,21
42,420,21
483,861,21
924,1302,21
1365,1743,21
1806,2184,21
2247,2625,21
2688,3066,21
3129,3507,21
*NSET,NSET=SNODE_2,GEN
421,441
862,882
1303,1323
1744,1764
2185,2205
2626,2646
3067,3087
*NSET,NSET=FILOUT,GEN
1324,1344
1744,1764
**
** C3D20 and C3D8 elements can also be used with this element definition
**
*ELEMENT, TYPE=C3D27R
1, 883, 885, 927, 925, 1, 3, 45, 43, 884, 906, 926, 904, 2, 24, 44, 
  22, 442, 444, 486, 484, 464, 905, 23, 443, 465, 485, 463
*ELGEN,ELSET=EALL
1,10,42,10,10,2,1,3,882,100
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT1
*ELSET,ELSET=SURF1,GEN
241,250
*ELSET,ELSET=SURF2,GEN
141,150
*ELSET,ELSET=SURF3,GEN
41,50
*surface,NAME=SIDE
SURF1,S5
SURF2,S5
SURF3,S5
*Elset, ELSET=ELOAD, GEN
201, 300
*Surface, Name=ELOAD
ELOAD, S1
*MATERIAL,NAME=MAT1
*ELASTIC,TYPE=ISO
28.1E6,0.2642
*PLASTIC
39440.0,0.0
50170.0,0.00473
54950.0,0.01264
58540.0,0.02836
61520.0,0.0491
76520.0,0.105
*SUBMODEL,EXT=0.05,ABS=0.06,SHELL TO SOLID,SHELL THICKNESS=0.41
SNODE_1,SNODE_2
***STEP,NLGEOM
** APPLY PRESSURE
***STATIC
**0.05,1.0
***BOUNDARY,SUBMODEL,STEP=1
**SNODE_1,0.005
**SNODE_2,0.005
***BOUNDARY
**YZ_PLANE,1,1
***DLOAD
**ELOAD,P1,500.0
***EL PRINT,F=0
***NODE PRINT,F=0
***END STEP
**
**
*INTEGRATED OUTPUT SECTION,NAME=IOS-1,SURFACE=SIDE,
REF NODE=10001,POSITION=CENTER,REFNODEMOTION=AVERAGETRANS,PROJECT=YES
**
*STEP,NLGEOM
 HOLD PRESSURE CONSTANT, APPLY BENDING ACTION
*STATIC
0.1, 0.963 
*BOUNDARY,SUBMODEL,STEP=1
SNODE_1,0.005
SNODE_2,0.005
*BOUNDARY
YZ_PLANE,1,1
*EL PRINT,F=4
S, 
SINV, 
E, 
*OUTPUT,HISTORY,OP=ADD,FREQ=4
*INTEGRATED OUTPUT,SECTION=IOS-1,VAR=PRESELECT
*SECTION PRINT,NAME=SID,SURFACE=SIDE,AXES=LOCAL,F=4
*SECTION FILE,NAME=SID,SURFACE=SIDE,AXES=LOCAL,F=999
*EL PRINT,F=999
*NODE PRINT,NSET=SNODE_1,F=4
U, 
*NODE PRINT,NSET=SNODE_2,F=4
U, 
*NODE PRINT,NSET=FILOUT
U, 
*NODE PRINT,F=999
*NODE FILE,NSET=FILOUT,F=20
U,RF
*END STEP