*HEADING
   THREE LAYER COMPOSITE PLATE:  PAGANO (1969) S4 ELEMENTS
*RESTART, WRITE, FREQ=999
*NODE
1, 0.,0.,0
11,10.,0.,0.
101,0.,1.0,0.
111,10.,1.0,0.
*NGEN,NSET=EDGE
1,11
101,111
*NSET,NSET=LEFT
1,101
*NSET,NSET=RIGHT
11,111
*ELEMENT,TYPE=S4
1,1,2,102,101
*ELGEN,ELSET=PLATE
1,10,1,1
*SHELL SECTION,ELSET=PLATE,COMPOSITE
0.83333,5,LAMINA,PL0
0.83333,5,LAMINA,PL90
0.83333,5,LAMINA,PL0
*MATERIAL,NAME=LAMINA
*ELASTIC,TYPE=LAMINA
25.E6,1.E6,.25,0.5E6,0.5E6,0.2E6
*FAIL STRESS
30.,-12.,.5,-1.5,1.0,0.
*FAIL STRAIN
17.E-2,-7.E-2,5.E-2,-1.3E-2,11.E-2
*ELSET,ELSET=FMEAS
5,
*ORIENTATION,NAME=PL0
1,0,0,0,1,0
3,0
*ORIENTATION,NAME=PL90
1,0,0,0,1,0
3,90
*BOUNDARY
LEFT,1,3
RIGHT,2,3
EDGE,2
EDGE,4
EDGE,6
**
*STEP,PERTURBATION
*STATIC
 1.0, 
*DLOAD
PLATE,PNU
*EL PRINT,POSITION=AVERAGED AT NODES
1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
S, 
TSHR, 
*EL PRINT
 SF, 
 SE, 
*EL FILE
1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
S, 
TSHR, 
CFAILURE,
*EL FILE
SF,SE
*NODE FILE
U, 
*NODE PRINT
U, 
*EL PRINT,ELSET=FMEAS,POSITION=AVERAGED AT NODES
1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
CFAILURE,
*OUTPUT,HISTORY,VAR=PRESELECT
*OUTPUT,FIELD,VAR=PRESELECT
*ELEMENT OUTPUT,POSITION=AVERAGED AT NODES
1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
S, 
TSHR, 
*END STEP