*HEADING
 RIKS ANALYSIS OF CYLINDER BUCKLING
*NODE,SYSTEM=C,NSET=ALL,INPUT=CYLIN.NOD
*NSET,NSET=BND1,GENERATE
*NSET,NSET=BND2,GENERATE
*NSET,NSET=BND3,GENERATE
*NSET,NSET=LOADB,GENERATE
*NSET,NSET=LDB1
*NSET,NSET=LDB2,GENERATE
*NSET,NSET=LDBFIL
 LDB1,LDB2
*TRANSFORM,TYPE=C,NSET=ALL
**********
*BOUNDARY
 BND1,YSYMM
 BND2,YSYMM
 BND3,ZSYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
**********
*ELEMENT,TYPE=S4R5
*ELGEN,ELSET=EALL
*SHELL SECTION,MATERIAL=MAT,ELSET=EALL
 .25,3
*MATERIAL,NAME=MAT
*ELASTIC
 30.E6,.3
*****
**********
*RESTART,WRITE,FREQUENCY=12
*STEP,NLGEOM,INC=60
*STATIC,RIKS
 .1,1.,,,,14,3,.6
*CLOAD
 LDB1,3, 49969.31158
 LDB2,3, 99938.62316
*NODE PRINT,NSET=BND3,FREQUENCY=1
U,
RF,
*NODE PRINT,NSET=LOADB,FREQUENCY=1
U,
RF,
*NODE PRINT,FREQUENCY=10
U,
RF,
*EL PRINT,FREQUENCY=99
 S,E,
*EL FILE,FREQUENCY=99
 S,E,
*NODE FILE,FREQUENCY=1,NSET=LDBFIL
U,CF,RF
*MONITOR,NODE=14,DOF=3
*END STEP