*HEADING
 NAFEMS TEST R0031/1, laminated strip under threepoint bending [SC6R]
 Modelled by stacking composite continuum shell elements
*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=SC6R,ELSET=ELALL
     1,   1,   2, 101, 12001, 12002, 12101
     2, 101,   2, 102, 12101, 12002, 12102
 12001,12001,12002,12101,13001,13002,13101
 12002,12101,12002,12102,13101,13002,13102
 13001,13001,13002,13101,16001,16002,16101
 13002,13101,13002,13102,16101,16002,16102
*ELGEN,ELSET=ELFACE_BOT 
 1, 10,1,2,2,100,100
 2, 10,1,2,2,100,100
*ELGEN,ELSET=ELCORE
 12001, 10,1,2,2,100,100
 12002, 10,1,2,2,100,100
*ELGEN,ELSET=ELFACE_TOP
 13001, 10,1,2,2,100,100
 13002, 10,1,2,2,100,100
*ELSET, ELSET=ELALL
 ELFACE_BOT,
 ELCORE,
 ELFACE_TOP,
*ELSET, ELSET=MID_FACE
 120, 12120, 13120, 20, 12020, 13020
*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