*HEADING
  NONLINEAR POSTBUCKLING ANALYSIS--10%  IMPERFECTION  S4  
  WITH CONTRIBUTIONS FROM FIRST FOUR MODES
  COMPOSITE CYLINDRICAL PANEL WITH CIRCULAR HOLE
**
**  FROM STUDY IN THESIS OF G.M. STANLEY
**
*NODE,INPUT=laminpanel_s4r5_node.inp
*NSET,NSET=OUTSIDE,GENERATE
2001,2017
2017,2033
2033,2049
2049,2065
*NSET,NSET=HOLE,GENERATE
1,17,1
17,33,1
33,49,1
49,65,1
*NSET,NSET=XPARB,GENERATE
2001,2017
*NSET,NSET=XPART,GENERATE
2033,2049
*NSET,NSET=YPAR,GENERATE
2049,2065
2017,2033
*NSET,NSET=PROFILE,GENERATE
57,2057,100
25,2025,100
*EQUATION
2, 
2033,3,1.,2041,3,-1.
2, 
2035,3,1.,2041,3,-1.
2, 
2037,3,1.,2041,3,-1.
2, 
2039,3,1.,2041,3,-1.
2, 
2043,3,1.,2041,3,-1.
2, 
2045,3,1.,2041,3,-1.
2, 
2047,3,1.,2041,3,-1.
2, 
2049,3,1.,2041,3,-1.
*NSET,NSET=MAIN
2041, 
*ELEMENT,TYPE=S4
1,1,201,203,3
*ELGEN,ELSET=ALL
1,31,2,1,10,200,100
*ELEMENT,TYPE=S4
32,63,263,201,1
*ELGEN,ELSET=ALL
32,10,200,100
*ELSET,ELSET=PROFILE,GENERATE
12, 912,100
13, 913,100
28, 928,100
29, 929,100
*SHELL GENERAL SECTION,ELSET=ALL
7.90239D5,2.51367D5,7.90239D5,-3.08578D-6,-7.94285D-5,2.69436D5,0.,0.
0., 4.92719D+02,0.,0.,0.,1.91513D+02,5.17951D+02,0.
0.,0.,  1.89245D+01 , 1.89245D+01 ,2.03602D+02
*TRANSVERSE SHEAR STIFFNESS
49573. ,    52967.    , 2.3382
*BOUNDARY
XPARB,1,6
XPART,1,2
XPART,4,6
YPAR,1,2
*RESTART,WRITE,FREQUENCY=5
*STEP,NLGEOM,INC=15
*STATIC,RIKS
1.,1.,,,,2041,3,-0.08
*CLOAD
2041,3,-10000.
*MONITOR,NODE=2041,DOF=3
*EL PRINT,ELSET=PROFILE,POSITION=AVERAGED AT NODES,FREQUENCY=5
SF, 
*NODE PRINT,NSET=XPART,TOTALS=YES,FREQUENCY=5
U, 
RF, 
*NODE PRINT,NSET=XPARB,TOTALS=YES,FREQUENCY=5
U, 
RF, 
*NODE FILE,FREQUENCY=1,NSET=MAIN
U,CF
*END STEP