*HEADING
: NAFEMS TEST 1, laminated strip under threepoint bending [C3D20]
*NODE
1,
11,25.
201,0.,5.
211,25.,5.
*NGEN,NSET=NFRONT
1,11
*NGEN,NSET=NBACK
201,211
*NFILL,NSET=NALL0
NFRONT,NBACK,2,100
*NCOPY,CHANGE NUMBER=1000,OLD SET=NALL0,SHIFT,NEW SET=NALL1
0.,0.,.05
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=2000,OLD SET=NALL0,SHIFT,NEW SET=NALL2
0.,0.,.1
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=3000,OLD SET=NALL0,SHIFT,NEW SET=NALL3
0.,0.,.15
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=4000,OLD SET=NALL0,SHIFT,NEW SET=NALL4
0.,0.,.2
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=5000,OLD SET=NALL0,SHIFT,NEW SET=NALL5
0.,0.,.25
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=6000,OLD SET=NALL0,SHIFT,NEW SET=NALL6
0.,0.,.3
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=7000,OLD SET=NALL0,SHIFT,NEW SET=NALL7
0.,0.,.5
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=8000,OLD SET=NALL0,SHIFT,NEW SET=NALL8
0.,0.,.7
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=9000,OLD SET=NALL0,SHIFT,NEW SET=NALL9
0.,0.,.75
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=10000,OLD SET=NALL0,SHIFT,NEW SET=NALL10
0.,0.,.8
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=11000,OLD SET=NALL0,SHIFT,NEW SET=NALL11
0.,0.,.85
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=12000,OLD SET=NALL0,SHIFT,NEW SET=NALL12
0.,0.,.9
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=13000,OLD SET=NALL0,SHIFT,NEW SET=NALL13
0.,0.,.95
0.,0.,0.,0.,0.,0.,0.
*NCOPY,CHANGE NUMBER=14000,OLD SET=NALL0,SHIFT,NEW SET=NALL14
0.,0.,1.
0.,0.,0.,0.,0.,0.,0.
*NSET,NSET=NYZ,GENERATE
11,14011,1000
111,14111,1000
211,14211,1000
*NSET,NSET=SUP
5,105,205
*NSET,NSET=NXZ,GENERATE
201,14201,1000
202,14202,1000
203,14203,1000
204,14204,1000
205,14205,1000
206,14206,1000
207,14207,1000
208,14208,1000
209,14209,1000
210,14210,1000
211,14211,1000
*ELEMENT,TYPE=C3D20,ELSET=ELALL0
1,1,3,203,201,7001,7003,7203,7201,2,103,202,101,
7002,7103,7202,7101,5001,5003,5203,5201
*ELEMENT,TYPE=C3D20R,ELSET=ELALL90
7001,7001,7003,7203,7201,14001,14003,14203,14201,7002,7103,7202,7101,
14002,14103,14202,14101,10001,10003,10203,10201
*ELGEN,ELSET=ELALLBOT
1,5,2,2
*ELGEN,ELSET=ELALLTOP
7001,5,2,2
*ELSET,ELSET=OUTPUT
9,7009
*SOLID SECTION,ELSET=ELALLBOT,COMPOSITE
.1,3,MAT,ORIENT0
.1,3,MAT,ORIENT90
.1,3,MAT,ORIENT0
.2,3,MAT,ORIENT90
*SOLID SECTION,ELSET=ELALLTOP,COMPOSITE
.2,3,MAT,ORIENT90
.1,3,MAT,ORIENT0
.1,3,MAT,ORIENT90
.1,3,MAT,ORIENT0
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ENGINEERING CONSTANTS
1.E5,5.E3,5.E3,.4,.3,.3,3.E3,2.E3,
2.E3, 
*ORIENTATION,NAME=ORIENT0
1.,0.,0.,0.,1.,0.
*ORIENTATION,NAME=ORIENT90
0.,1.,0.,-1.,0.,0.
*BOUNDARY
NXZ,2
NXZ,4
NXZ,6
NYZ,1
NYZ,5,6
SUP,3
**
*STEP
*STATIC
1.,1.
*CLOAD
14011,3,-4.6665
14111,3,-16.6665
14211,3,-4.6665
*EL FILE,ELSET=OUTPUT,POSITION=NODES
1,2,3,4,5,6,7,8,9,10,11,12
S, 
*END STEP