*HEADING
 EIGENVALUE BUCKLING OF CYLINDER UNDER AXIAL LOAD (N=7,ALPHA=PI/14)
*NODE
*NGEN,LINE=C,NSET=LOADB
*NGEN,LINE=C,NSET=BND3
*NGEN,NSET=BND1
*NGEN,NSET=REST
*NGEN,NSET=BND2
*NSET,NSET=ALL
 BND1,BND2,BND3,REST
*NSET,NSET=LDB1
*NSET,NSET=LDB2
*****
*ELEMENT,TYPE=S4R5
*ELGEN,ELSET=EALL
*ELSET,ELSET=EPRINT
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT
 .25,
*TRANSFORM,TYPE=C,NSET=ALL
*MATERIAL,NAME=MAT
*ELASTIC
 30.E6,.3
**********
*STEP
*BUCKLE
1, , 10, 100
*CLOAD
 LDB1,3,1500.
 LDB2,3,3000.
*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,ZASYMM
 LOADB,1,2
 LOADB,4
 LOADB,6
*MODAL FILE
*END STEP