*HEADING
 SHELL TO SOLID SUBMODEL
*PREPRINT,ECHO=NO
*RESTART,WRITE
*SYSTEM
*NODE
*NGEN,NSET=A1,LINE=C
*NGEN,NSET=A2,LINE=C
*SYSTEM
*NFILL,NSET=SNODE_1
A1,A2,6,441
*NCOPY,CHANGE NUMBER=21,OLD SET=SNODE_1,NEW SET=ELBOW,SHIFT,MULTIPLE=20
*NSET,NSET=YZ_PLANE,GEN
*NSET,NSET=SNODE_2,GEN
*NSET,NSET=FILOUT,GEN
**
** C3D20 and C3D8 elements can also be used with this element definition
**
*ELEMENT, TYPE=C3D27R
*ELGEN,ELSET=EALL
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT1
*ELSET,ELSET=SURF1,GEN
*ELSET,ELSET=SURF2,GEN
*ELSET,ELSET=SURF3,GEN
*surface,NAME=SIDE
SURF1,S5
SURF2,S5
SURF3,S5
*Elset, ELSET=ELOAD, GEN
*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