*HEADING
test for viscoelasticity with anisotropic elasticity
The elasticity properties are isotropic but defined using 
TYPE=ENGINEERING CONSTANTS or TYPE=LAMINA
*NODE, NSET=C3D8R
1, 0., 0.,
2,10., 0.,
3,10.,10.,
4, 0.,10.,
5, 0., 0., 10.
6,10., 0., 10.
7,10.,10., 10.
8, 0.,10., 10.
101, 0.,-20.,
102,10.,-20.,
103,10.,-10.,
104, 0.,-10.,
105, 0.,-20., 10.
106,10.,-20., 10.
107,10.,-10., 10.
108, 0.,-10., 10.
*ELEMENT, TYPE=C3D8R, ELSET=C3D8R
1, 1,2,3,4,5,6,7,8
101, 101,102,103,104,105,106,107,108
*SOLID SECTION, ELSET=C3D8R, MATERIAL=POLY, CONTROLS=HOURGLASS, ORIENTATION=ORI
*NODE, NSET=CPE4R
11,20., 0.,
12,30., 0.,
13,30.,10.,
14,20.,10.,
111,20.,-20.,
112,30.,-20.,
113,30.,-10.,
114,20.,-10.,
*ELEMENT, TYPE=CPE4R, ELSET=CPE4R
11,11,12,13,14
111,111,112,113,114
*SOLID SECTION, ELSET=CPE4R, MATERIAL=POLY, CONTROLS=HOURGLASS, ORIENTATION=ORI
*NODE, NSET=CPS4R
21,40., 0.,
22,50., 0.,
23,50.,10.,
24,40.,10.,
121,40.,-20.,
122,50.,-20.,
123,50.,-10.,
124,40.,-10.,
*ELEMENT, TYPE=CPS4R, ELSET=CPS4R
21,21,22,23,24
121,121,122,123,124
*SOLID SECTION, ELSET=CPS4R, MATERIAL=POLY, CONTROLS=HOURGLASS, ORIENTATION=ORI
*NODE, NSET=CPS4R-1
221,100., 0.,
222,110., 0.,
223,110.,10.,
224,100.,10.,
321,100.,-20.,
322,110.,-20.,
323,110.,-10.,
324,100.,-10.,
*ELEMENT, TYPE=CPS4R, ELSET=CPS4R-1
221,221,222,223,224
321,321,322,323,324
*SOLID SECTION, ELSET=CPS4R-1, MATERIAL=POLY-LAMINA, CONTROLS=HOURGLASS, ORIENTATION=ORI
*NODE, NSET=S4R
31,60., 0.,
32,70., 0.,
33,70.,10.,
34,60.,10.,
131,60.,-20.,
132,70.,-20.,
133,70.,-10.,
134,60.,-10.,
*ELEMENT, TYPE=S4R, ELSET=S4R
31,31,32,33,34
131,131,132,133,134
*SHELL SECTION, ELSET=S4R, MATERIAL=POLY, CONTROLS=HOURGLASS, ORIENTATION=ORI
1.0, 3
*NODE, NSET=S4R-1
231,120., 0.,
232,130., 0.,
233,130.,10.,
234,120.,10.,
331,120.,-20.,
332,130.,-20.,
333,130.,-10.,
334,120.,-10.,
*ELEMENT, TYPE=S4R, ELSET=S4R-1
231,231,232,233,234
331,331,332,333,334
*SHELL SECTION, ELSET=S4R-1, MATERIAL=POLY-LAMINA, CONTROLS=HOURGLASS, ORIENTATION=ORI
1.0, 3
*NODE, NSET=M3D4R
41,80., 0.,
42,90., 0.,
43,90.,10.,
44,80.,10.,
141,80.,-20.,
142,90.,-20.,
143,90.,-10.,
144,80.,-10.,
*ELEMENT, TYPE=M3D4R, ELSET=M3D4R
41,41,42,43,44
141,141,142,143,144
*MEMBRANE SECTION, ELSET=M3D4R, MATERIAL=POLY, CONTROLS=HOURGLASS
1.0,
*NODE, NSET=M3D4R-1
241,140., 0.,
242,150., 0.,
243,150.,10.,
244,140.,10.,
341,140.,-20.,
342,150.,-20.,
343,150.,-10.,
344,140.,-10.,
*ELEMENT, TYPE=M3D4R, ELSET=M3D4R-1
241,241,242,243,244
341,341,342,343,344
*MEMBRANE SECTION, ELSET=M3D4R-1, MATERIAL=POLY-LAMINA, CONTROLS=HOURGLASS
1.0,
*SECTION CONTROLS, NAME=HOURGLASS, HOURGLASS=STIFFNESS
*ORIENTATION, NAME=ORI
1.,0.,0.,0.,1.,0.
3,0.
*MATERIAL,NAME=POLY
*DENSITY
1.0e-5,
*ELASTIC,MODULI=LONG TERM, TYPE=ENGINEERING CONSTANTS
30,30,30,0.4,0.4,0.4,10.71428571,10.71428571
10.71428571
*VISCOELASTIC,TIME=PRONY
0.25,0.05,5.
0.25,0.15,10.
*TRS, DEFINITION=WLF
 5.25, 0.163277, 269.23
*MATERIAL,NAME=POLY-LAMINA
*DENSITY
1.0e-5,
*ELASTIC,MODULI=LONG TERM, TYPE=LAMINA
30,30,0.499,10.00667111,10.00667111,10.00667111
*VISCOELASTIC,TIME=PRONY
0.25,0.05,5.
0.25,0.15,10.
*TRS, DEFINITION=WLF
 5.25, 0.163277, 269.23
*NSET, NSET=BOTTOM
101,102,105,106,111,112,121,122,131,132,141,142,
321,322,331,332,341,342
*NSET, NSET=TOP
103,104,107,108,113,114,123,124,133,134,143,144,
323,324,333,334,343,344,
*BOUNDARY
BOTTOM,1,6
TOP,2,6
*BOUNDARY
1,1,3
2,2,3
3,3,3
4,1,1
4,3,3
5,1,2
6,2,2
8,1,1
11,1,2
12,2,2
14,1,1
21,1,2
22,2,2
24,1,1
31,1,2
32,2,2
34,1,1
41,1,2
42,2,2
44,1,1
221,1,2
222,2,2
224,1,1
231,1,2
232,2,2
234,1,1
241,1,2
242,2,2
244,1,1
S4R,4,6
S4R-1,4,6
*NSET, NSET=NALL
C3D8R,CPE4R,CPS4R,S4R,M3D4R,CPS4R-1,S4R-1,M3D4R-1
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL, 615.0
*AMPLITUDE, NAME=RELAXATION, TIME=TOTAL TIME
0., 0., 2.0, 0.2, 20.0, 0.2
*STEP
*DYNAMIC, EXPLICIT, IMPROVED DT METHOD=NO
 ,5.
*BOUNDARY, AMPLITUDE=RELAXATION
3,2,2,1.0
4,2,2,1.0
7,2,2,1.0
8,2,2,1.0
13,2,2,1.0
14,2,2,1.0
23,2,2,1.0
24,2,2,1.0
33,2,2,1.0
34,2,2,1.0
43,2,2,1.0
44,2,2,1.0
223,2,2,1.0
224,2,2,1.0
233,2,2,1.0
234,2,2,1.0
243,2,2,1.0
244,2,2,1.0
TOP,1,1,1.0
*TEMPERATURE
NALL, 615.0
*RESTART, WRITE, NUMBER=10
*NSET, NSET=QA_TEST_NSET
C3D8R,CPE4R,CPS4R,S4R,M3D4R,CPS4R-1,S4R-1,M3D4R-1
*ELSET, ELSET=QA_TEST_ELSET
C3D8R,CPE4R,CPS4R,S4R,M3D4R,CPS4R-1,S4R-1,M3D4R-1
*OUTPUT, FIELD
*NODE OUTPUT
U,V,RF
*ELEMENT OUTPUT
S,LE
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=QA_TEST_ELSET
S,LE
*ENERGY OUTPUT
ALLAE,ALLWK,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD,ETOTAL
*ENERGY OUTPUT, ELSET=C3D8R
ALLAE,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD
*ENERGY OUTPUT, ELSET=CPE4R
ALLAE,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD
*ENERGY OUTPUT, ELSET=CPS4R
ALLAE,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD
*ENERGY OUTPUT, ELSET=S4R
ALLAE,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD
*ENERGY OUTPUT, ELSET=M3D4R
ALLAE,ALLKE,ALLIE,ALLSE,ALLCD,ALLVD
*END STEP