*HEADING
 EIGENVALUE BUCKLING OF CYLINDER UNDER AXIAL LOAD (N=4,ALPHA=PI/8) 
*NODE
 10,  100.,0.,0.
 18,  92.39,0.,38.27
 410, 100.,400.,0.
 418, 92.39,400.,38.27
*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,10
*NGEN,NSET=REST
 11,411,10
 12,412,10
 13,413,10
 14,414,10
 15,415,10
 16,416,10
 17,417,10
*NGEN,NSET=BND2
 18,418,10
*NSET,NSET=ALL
 BND1,BND2,BND3,REST
*NSET,NSET=LDB1
 10,18
*NSET,NSET=LDB4
 11,13,15,17
*NSET,NSET=LDB2
 12,14,16
*****
*ELEMENT,TYPE=STRI65
    1,    10,    30,    32,    20,    31,    21
50001,    10,    32,    12,    21,    22,    11
*ELGEN,ELSET=EALL
    1,     4,     2,     2,    20,    20,     8
50001,     4,     2,     2,    20,    20,     8
*ELSET,ELSET=EPRINT
 1,2,3,4
*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
*****
**********
*STEP
*BUCKLE
1, , 10,  40
*CLOAD
 LDB1,3,500.0
 LDB2,3,1000.0
 LDB4,3,2000.0
*BOUNDARY
 BND1,YSYMM
 BND2,YSYMM
 BND3,ZSYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
*BOUNDARY,LOAD CASE=2, OP=NEW
 BND1,YSYMM
 BND2,YASYMM
 BND3,ZSYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
*MODAL FILE
*END STEP