*HEADING
 RIKS ANALYSIS OF CYLINDER BUCKLING
 (1% thickness,*IMPERFECTION)
*NODE
*NGEN,LINE=C,NSET=LOADB
*NGEN,LINE=C,NSET=BND3
*IMPERFECTION, FILE=bucklecylshell_s4r5_n1, STEP=1
*NGEN,NSET=BND1
*NGEN,NSET=REST
*NGEN,NSET=BND2
*NSET,NSET=ALL
 BND1,BND2,BND3,REST
*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,ELSET=EALL,MATERIAL=MAT
 .25,
*MATERIAL,NAME=MAT
*ELASTIC
 30.E6,.3
*RESTART,WRITE
**
*STEP,NLGEOM,INC=60
*STATIC,RIKS
 .1,1.,,,,14,3,.6
*CLOAD
 LDB1,3, 49969.31158
 LDB2,3, 99938.62316
*EL PRINT, FREQUENCY=99
 S,E,
*EL FILE,FREQUENCY=99
 S,E,
*NODE PRINT,NSET=LDBFIL,FREQUENCY=2,GLOBAL=NO
 U,
 CF,
 RF,
*NODE FILE,FREQUENCY=1,NSET=LDBFIL,GLOBAL=NO
 U,CF,RF
*MONITOR,NODE=14,DOF=3
*END STEP