*HEADING
  COMPOSITE CYLINDRICAL PANEL WITH CIRCULAR HOLE--S4
**
**  FROM STUDY IN THESIS OF G.M. STANLEY
**
*NODE
*NGEN,NSET=OUTSIDE
**
** DEFINE CIRCLE
**
*SYSTEM
*NODE
*NGEN,LINE=C,NSET=HOLE
*SYSTEM
*NFILL,NSET=ALL
HOLE,OUTSIDE,20,100
*NSET,NSET=XPARB,GENERATE
*NSET,NSET=XPART,GENERATE
*NSET,NSET=YPAR,GENERATE
*NSET,NSET=PROFILE,GENERATE
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
*ELEMENT,TYPE=S4
*ELGEN,ELSET=ALL
*ELEMENT,TYPE=S4
*ELGEN,ELSET=ALL
*ELSET,ELSET=PROFILE,GENERATE
*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
*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
*END STEP