*HEADING
 RIKS ANALYSIS OF CYLINDER BUCKLING        
*NODE,SYSTEM=C,NSET=ALL,INPUT=CYLIN.NOD
*NSET,NSET=BND1,GENERATE
 10,410,10
*NSET,NSET=BND2,GENERATE
 18,418,10
*NSET,NSET=BND3,GENERATE
 410,418
*NSET,NSET=LOADB,GENERATE
 10,18
*NSET,NSET=LDB1
 10,18
*NSET,NSET=LDB2,GENERATE
 11,17
*NSET,NSET=LDBFIL
 LDB1,LDB2
*TRANSFORM,TYPE=C,NSET=ALL
 0.,0.,0., 0.,0.,1.
**********
*BOUNDARY
 BND1,YSYMM
 BND2,YSYMM
 BND3,ZSYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
**********
*ELEMENT,TYPE=S4R5
 1, 10,50,51,11
*ELGEN,ELSET=EALL
 1, 8,1,1, 10,40,10
*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