*HEADING
  COMPOSITE CYLINDRICAL PANEL WITH CIRCULAR HOLE--S4 
**
**  FROM STUDY IN THESIS OF G.M. STANLEY
**
*RESTART,WRITE,FREQUENCY=10
*NODE
2001,15.,     0.,  0.
2017,15.,  14.556, 0.
2033,15.,  14.556, 14.
2049,15.,     0. , 14.
2065,15.,     0.,  0.
*NGEN,NSET=OUTSIDE
2001,2017
2017,2033
2033,2049
2049,2065
**
** DEFINE CIRCLE
**
*SYSTEM
15.,7.27802298,7.0
*NODE
9999, 0.,0.
1,0., -.70710678, -.70710678
17,0., .70710678,-.70710678
33,0., .70710678, .70710678
49,0., -.70710678, .70710678
65,0., -.70710678, -.70710678
*NGEN,LINE=C,NSET=HOLE
1,17,1,9999
17,33,1,9999
33,49,1,9999
49,65,1,9999
*SYSTEM
0.,0.,0.
*NFILL,NSET=ALL
HOLE,OUTSIDE,20,100
*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
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
0.,0.,0.,   0.,0.,1.
15.,0.,0.
1.,3.8197307  ,1.
*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.
*ELEMENT,TYPE=S4
1,1,201,203,3,101,202,103,2,102
*ELGEN,ELSET=ALL
1,31,2,1,10,200,100
*ELEMENT,TYPE=S4
32,63,263,201,1,163,264,101,64,164
*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
*STEP
*BUCKLE
5,,50,20
*CLOAD
2041,3,-1000.
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*MODAL FILE
*OUTPUT,HISTORY,FREQUENCY=1
*MODAL OUTPUT
*END STEP