*HEADING
 NAFEMS TEST R0031/1, laminated strip under threepoint bending [SC6R]
 Modelled with 1 composite element through the thickness
*PREPRINT, MODEL=YES
*NODE
     1,  0., 0., -0.5
    11, 25., 0., -0.5
   201,  0., 5., -0.5
   211, 25., 5., -0.5
 10001,  0., 0., +0.5
 10011, 25., 0., +0.5
 10201,  0., 5., +0.5
 10211, 25., 5., +0.5
*NGEN,NSET=NFRONT
 1,11
 10001,10011
*NGEN,NSET=NBACK
 201,211
 10201, 10211
*NFILL,NSET=NALL
 NFRONT,NBACK,2,100
*NSET,NSET=NRT
    11,  111,  211
 10011,10111,10211
*NSET,NSET=SUP
     5,  105,   205
 10005,10105, 10205
*ELEMENT,TYPE=SC6R,ELSET=ELALL
 1,   1,   2, 101, 10001, 10002, 10101
 2, 101,   2, 102, 10101, 10002, 10102
*ELGEN,ELSET=ELALL
 1, 10,1,2,2,100,100
 2, 10,1,2,2,100,100
*SHELL SECTION,ELSET=ELALL,COMPOSITE
 .1,3,MAT,0.
 .1,3,MAT,90.
 .1,3,MAT,0.
 .4,3,MAT,90.
 .1,3,MAT,0.
 .1,3,MAT,90.
 .1,3,MAT,0.
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 1.E5,5.E3,5.E3,.4,.3,.3,3.E3,2.E3,
 2.E3, 
*BOUNDARY
 NBACK,2
 NRT,1
 NRT,2
** PIN SUPPORTS AT A AND B 
*EQUATION
 2, 
   5, 3, 1.0, 10005, 3, 1.0
 2, 
 105, 3, 1.0, 10105, 3, 1.0
 2, 
 205, 3, 1.0, 10205, 3, 1.0
*STEP
*STATIC
 1.,1.
*CLOAD
 10011,3,-6.25
 10111,3,-12.5
 10211,3,-6.25
*NODE PRINT
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 1, 2, 3, 4, 5, 6, 19, 20, 21
 S,
 TSHR,
*NODE FILE 
 U,
*EL FILE
 1, 2, 3, 4, 5, 6, 19, 20, 21
 S,
 TSHR,
*END STEP