*HEADING
 NAFEMS TEST R0031/1, laminated strip under threepoint bending [SC8R]
 Modelled by stacking composite continuum shell
*PREPRINT, MODEL=YES
*NODE
     1,  0., 0., -0.5
    11, 25., 0., -0.5
   201,  0., 5., -0.5
   211, 25., 5., -0.5
 12001,  0., 0., -0.2
 12011, 25., 0., -0.2
 12201,  0., 5., -0.2
 12211, 25., 5., -0.2
 13001,  0., 0., +0.2
 13011, 25., 0., +0.2
 13201,  0., 5., +0.2
 13211, 25., 5., +0.2
 16001,  0., 0., +0.5
 16011, 25., 0., +0.5
 16201,  0., 5., +0.5
 16211, 25., 5., +0.5
*NODE
 999999999, 25.0, 2.5, 0.0
*NGEN,NSET=NFRONT
 1,11
 12001,12011
 13001,13011
 16001,16011
*NGEN,NSET=NBACK
 201,211
 12201, 12211
 13201, 13211
 16201, 16211
*NFILL,NSET=NALL
 NFRONT,NBACK,2,100
*NSET,NSET=NRT
    11,  111,  211
 12011,12111,12211
 13011,13111,13211 
 16011,16111,16211
*NSET,NSET=SUP
     5,  105,   205
 12005,12105, 12205
 13005,13105, 13205
 16005,16105, 16205
*ELEMENT,TYPE=SC8R,ELSET=ELALL
     1,    1,    2,  102,  101,12001,12002,12102,12101
 12001,12001,12002,12102,12101,13001,13002,13102,13101
 13001,13001,13002,13102,13101,16001,16002,16102,16101
*ELGEN,ELSET=ELFACE_BOT 
    1, 10,1,1,2,100,100
*ELGEN,ELSET=ELCORE
12001, 10,1,1,2,100,100
*ELGEN,ELSET=ELFACE_TOP
13001, 10,1,1,2,100,100
*ELSET, ELSET=ELALL
ELFACE_BOT,
ELCORE,
ELFACE_TOP,
*ELSET, ELSET=MID_FACE
 110, 12110, 13110, 10, 12010, 13010
*SURFACE, NAME=SURF
 MID_FACE, S4
*COUPLING, CONSTRAINT NAME=DCOUP, REF=999999999, SURFACE=SURF
*DISTRIBUTING  
*SHELL SECTION,ELSET=ELFACE_BOT,COMPOSITE
 .1,3,MAT,0.
 .1,3,MAT,90.
 .1,3,MAT,0.
*SHELL SECTION,ELSET=ELCORE,COMPOSITE
 .4,3,MAT,90.
*SHELL SECTION,ELSET=ELFACE_TOP,COMPOSITE 
 .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, 16005, 3, 1.0
 2, 
 105, 3, 1.0, 16105, 3, 1.0
 2, 
 205, 3, 1.0, 16205, 3, 1.0
*STEP
*STATIC
 1.,1.
*CLOAD
 999999999, 3, -25.0
*NODE PRINT
 U,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 1, 2 ,3, 4, 5, 6
 S,
 TSHR,
*END STEP