*HEADING
 NAFEMS TEST R0031/1, laminated strip under threepoint bending [SC6R]
 Modelled by stacking homogeneous elements 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.4
 10011, 25., 0., -0.4
 10201,  0., 5., -0.4
 10211, 25., 5., -0.4
 11001,  0., 0., -0.3
 11011, 25., 0., -0.3
 11201,  0., 5., -0.3
 11211, 25., 5., -0.3
 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
 14001,  0., 0., +0.3
 14011, 25., 0., +0.3
 14201,  0., 5., +0.3
 14211, 25., 5., +0.3
 15001,  0., 0., +0.4
 15011, 25., 0., +0.4
 15201,  0., 5., +0.4
 15211, 25., 5., +0.4
 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
 10001,10011
 11001,11011
 12001,12011
 13001,13011
 14001,14011
 15001,15011
 16001,16011
*NGEN,NSET=NBACK
 201,211
 10201, 10211
 11201, 11211
 12201, 12211
 13201, 13211
 14201, 14211
 15201, 15211
 16201, 16211
*NFILL,NSET=NALL
 NFRONT,NBACK,2,100
*NSET,NSET=NRT
    11,  111,  211
 10011,10111,10211 
 11011,11111,11211
 12011,12111,12211
 13011,13111,13211
 14011,14111,14211
 15011,15111,15211
 16011,16111,16211
*NSET,NSET=SUP
     5,  105,   205
 10005,10105, 10205
 11005,11105, 11205
 12005,12105, 12205
 13005,13105, 13205
 14005,14105, 14205
 15005,15105, 15205
 16005,16105, 16205
*ELEMENT,TYPE=SC6R,ELSET=ELALL
     1,   1,   2, 101, 10001, 10002, 10101
     2, 101,   2, 102, 10101, 10002, 10102
 10001,10001,10002,10101,11001,11002,11101
 10002,10101,10002,10102,11101,11002,11102
 11001,11001,11002,11101,12001,12002,12101
 11002,11101,11002,11102,12101,12002,12102
 12001,12001,12002,12101,13001,13002,13101
 12002,12101,12002,12102,13101,13002,13102
 13001,13001,13002,13101,14001,14002,14101
 13002,13101,13002,13102,14101,14002,14102
 14001,14001,14002,14101,15001,15002,15101
 14002,14101,14002,14102,15101,15002,15102
 15001,15001,15002,15101,16001,16002,16101
 15002,15101,15002,15102,16101,16002,16102
*ELGEN,ELSET=ELFACE_0 
     1, 10,1,2,2,100,100
     2, 10,1,2,2,100,100
 11001, 10,1,2,2,100,100
 11002, 10,1,2,2,100,100
 13001, 10,1,2,2,100,100
 13002, 10,1,2,2,100,100
 15001, 10,1,2,2,100,100
 15002, 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_90
 10001, 10,1,2,2,100,100
 10002, 10,1,2,2,100,100
 14001, 10,1,2,2,100,100
 14002, 10,1,2,2,100,100
*ELSET, ELSET=ELALL
ELFACE_0,
ELCORE,
ELFACE_90,
*ELSET, ELSET=MID_FACE
 120, 10120, 11120, 12120, 13120, 14120, 15120, 
  20, 10020, 11020, 12020, 13020, 14020, 15020, 
*SURFACE, NAME=SURF
 MID_FACE, S4
*COUPLING, CONSTRAINT NAME=DCOUP, REF=999999999, SURFACE=SURF
*DISTRIBUTING  
*SHELL SECTION,ELSET=ELFACE_0,COMPOSITE
 .1,3,MAT,0.
*TRANSVERSE SHEAR STIFFNESS
 1.25E3,1.25E3
*SHELL SECTION,ELSET=ELCORE,COMPOSITE
 .4,3,MAT,90.
*TRANSVERSE SHEAR STIFFNESS
 1.25E3,1.25E3
*SHELL SECTION,ELSET=ELFACE_90,COMPOSITE 
 .1,3,MAT,90.
*TRANSVERSE SHEAR STIFFNESS
 1.25E3,1.25E3
*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
** 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
 S,
 TSHR,
*END STEP