*HEADING
  COMPOSITE CYLINDRICAL PANEL WITH CIRCULAR HOLE-
  S9R5   
**  FROM STUDY IN THESIS OF G.M. STANLEY
*RESTART,WRITE,FREQUENCY=1
*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.
*ELEMENT,TYPE=S9R5
1,1,201,203,3,101,202,103,2,102
*ELGEN,ELSET=ALL
1,31,2,1,10,200,100
*ELEMENT,TYPE=S9R5
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,COMPOSITE,
ORIENTATION=SECORI
.0056,,LAMINA, 45.
.0056,,LAMINA,-45.
.0056,,LAMINA, 90.
.0056,,LAMINA,  0.
.0056,,LAMINA,  0.
.0056,,LAMINA, 90.
.0056,,LAMINA,-45.
.0056,,LAMINA, 45.
*** CENTER LINE
.0056,,LAMINA, 45.
.0056,,LAMINA,-45.
.0056,,LAMINA, 90.
.0056,,LAMINA,  0.
.0056,,LAMINA,  0.
.0056,,LAMINA, 90.
.0056,,LAMINA,-45.
.0056,,LAMINA, 45.
*MATERIAL,NAME=LAMINA
*ELASTIC,TYPE=LAMINA
19.6E6,  1.89E6, .38, .93E6, .93E6, .63E6
*ORIENTATION,SYSTEM=CYLINDRICAL,NAME=SECORI
0.,0.,0., 0.,0., 1.
1, 0.
*STEP
*STATIC
*BOUNDARY
XPARB,1,6
XPART,1,2
XPART,4,6
YPAR,1,2
XPART,3,,-.0316
*EL PRINT,ELSET=PROFILE,POSITION=AVERAGED AT NODES
SF, 
*EL FILE,ELSET=PROFILE,POSITION=AVERAGED AT NODES
SF, 
*NODE PRINT,NSET=XPART,TOTALS=YES
U, 
RF, 
*NODE PRINT,NSET=XPARB,TOTALS=YES
U, 
RF, 
*NODE FILE,NSET=XPART
U,RF
*OUTPUT,FIELD,VAR=PRESELECT
*ELEMENT OUTPUT,ELSET=PROFILE,POSITION=AVERAGED AT NODES
SF, 
*OUTPUT,HIST
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP