*HEADING
 MLTF ELBOW: IN-PLANE BENDING. SYM.   
 HALF S8R5 MODEL,LAMDA=.167
*NSET,NSET=N100
100, 
*NODE
 100,0.,24.,-24.
 101,0.,31.808,-24.
 113,0.,16.192,-24.
 900,0.,24.0,0.
 901,0.,31.808,0.
 913,0.,16.192,0.
2100,0.,0.,24.
2101,0.,0.,31.808
2113,0.,0.,16.192
4500,0.,-72.,24.
4501,0.,-72.,31.808
4513,0.,-72.,16.192
*NGEN,NSET=PYZ
100,900,100
101,901,100
113,913,100
2100,4500,100
2101,4501,100
2113,4513,100
*NGEN,LINE=C,NSET=PYZ
 900,2100,100,0,0.,0.,0.,1.
 901,2101,100,0,0.,0.,0.,1.
 913,2113,100,0,0.,0.,0.,1.
*NGEN,LINE=C,NSET=LOADEND
101,113,1,0,0.,24.,-24.,0.,0.,-1.
*NGEN,LINE=C ,NSET=END2
901,913,1,0,0.,24.,  0.0,0.,0.,-1.
*NGEN,LINE=C,NSET=END3
2101,2113,1,0,0.,0.,24., 0.,1.,0.
*NGEN,LINE=C ,NSET=FIXEDEND
4501,4513,1,0,0.,-72.,24., 0.,1.,0.
*NFILL
LOADEND,END2,8,100
END3,FIXEDEND,24,100
*NGEN,LINE=C
1001,1013,1,1000,,,,0.,.130526192,-.9914448614
1101,1113,1,1100,,,,0.,.258819045,-.965925826
1201,1213,1,1200,,,,0.,.382683432,-.923879532
1301,1313,1,1300,,,,0.,.5,-.866025404
1401,1413,1,1400,,,,0.,.60876143,-.79335334
1501,1513,1,1500,,,,0.,.707106781,-.707106781
1601,1613,1,1600,,,,0.,.79335334,-.60876143
1701,1713,1,1700,,,,0.,.866025404,-.5
1801,1813,1,1800,,,,0.,.923879532,-.382683432
1901,1913,1,1900,,,,0.,.965925826,-.258819045
2001,2013,1,2000,,,,0.,.9914448614,-.130526192
*ELEMENT,TYPE=S8R5
101,101,301,303,103,201,302,203,102
1101,2101,2301,2303,2103,2201,2302,2203,2102
*ELGEN,ELSET=LEG
101,6,2,1,4,200,100
1101,6,2,1,12,200,100
*ELEMENT,TYPE=S8R5
501,901,1101,1103,903,1001,1102,1003,902
*ELGEN,ELSET=LBOW
501,6,2,1,6,200,100
*ELSET,ELSET=ALL
 LEG,LBOW
*MATERIAL ,NAME=PIPE
*ELASTIC
28.1E6,.2642
*PLASTIC
39440.,0.
50170.,.00473
54950.,.01264
58540.,.02836
61520.,.0491
76520.,.105
*SHELL SECTION,ELSET=LEG,MATERIAL=PIPE
.37,7
*SHELL SECTION,ELSET=LBOW,MATERIAL=PIPE
.41,7
*BOUNDARY
PYZ,1
PYZ,5,6
FIXEDEND,1,6
*NSET,NSET=COUPLED_END,GENERATE
101,113,1
*KINEMATIC COUPLING, REF NODE = 100
COUPLED_END,
*NSET,NSET=NOUT
N100,LOADEND,END2,END3,FIXEDEND
*RESTART,WRITE,FREQUENCY=5
** *STEP,INC=1,NLGEOM
**   APPLY PRESSURE
** *STATIC
** *DLOAD,OP=NEW
**  ALL,P,500.
** *CLOAD,OP=NEW,FOLLOWER
**  100,3,-4.7882E4
** *EL PRINT,FREQUENCY=0
**  PE
** *NODE PRINT,NSET=N100
**  U
**  RF
** *NODE FILE,NSET=NOUT
**  U,RF
** *END STEP
*STEP,NLGEOM
HOLD PRESSURE CONSTANT, APPLY BENDING ACTION
*STATIC,RIKS
 .05, 1.,.01,.5,2.0,100,4,-.2
*BOUNDARY
100,4,4,-.2
*EL PRINT,FREQUENCY=0
PE, 
*NODE PRINT,NSET=N100
U, 
RF, 
*NODE FILE,NSET=NOUT
U,RF
*END STEP