*HEADING
 EIGENVALUE BUCKLING OF CYLINDER UNDER AXIAL LOAD  
 (N=1,ALPHA=PI/4)  
*NODE
 10, 100.,0.,0.
 18, 70.71,0.,70.71
 410,100.,400.,0.
 418,70.71,400.,70.71
*NGEN,LINE=C,NSET=LOADB
  10,18,1,,0.,0.,0.
*NGEN,LINE=C,NSET=BND3
 410,418,1, ,0.,400.,0.
*NGEN,NSET=BND1
 10,410,40
*NGEN,NSET=REST
 11,411,40
 12,412,40
 13,413,40
 14,414,40
 15,415,40
 16,416,40
 17,417,40
*NGEN,NSET=BND2
 18,418,40
*NSET,NSET=ALL
 BND1,BND2,BND3,REST
*NSET,NSET=LDB1
 10,18
*NSET,NSET=LDB2,GENERATE
 11,17
*ELEMENT,TYPE=S4R5
 1, 10,50,51,11
*ELGEN,ELSET=EALL
 1,8,1,1,10,40,10
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT
 .25, 
*TRANSFORM,TYPE=C,NSET=ALL
 0.,0.,0., 0.,1.,0.
*MATERIAL,NAME=MAT
*ELASTIC
 30.E6,.3
*RESTART,WRITE
**
*STEP
*BUCKLE
1,,20, 80
*BOUNDARY
 BND1,YSYMM
 BND2,YSYMM
 BND3,ZSYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
*CLOAD
 LDB1,3,1500.
 LDB2,3,3000.
*NODE PRINT, GLOBAL=YES
U, 
*NODE FILE, LAST MODE=1, GLOBAL=YES
U, 
*END STEP