*HEADING
 TRANSIENT THERMAL LOADING OF VISCOELASTIC SLAB -- CPE8R 
*NODE
1,0.,0.
3,0.2,0.
7,0.5,0.
13,0.8,0.
19,0.95,0.
21,1.,0.
201,0.,1.
203,0.2,1.
207,0.5,1.
213,0.8,1.
219,0.95,1.
221,1.,1.
*NGEN,NSET=BOT
1,3
3,7
7,13
13,21
*NGEN,NSET=TOP
201,203
203,207
207,213
213,221
*NFILL
BOT,TOP,2,100
*NSET,NSET=RIGHT
3,103,203
*NSET,NSET=LEFT
1,101,201
*ELEMENT,TYPE=CPE8R,ELSET=SOLID
1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=SOLID
1,10,2,1
*ELSET,ELSET=LEFT
1,
*ELSET,ELSET=RIGHT
10,
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT1
*MATERIAL,NAME=MAT1
*ELASTIC,MODULI=INSTANTANEOUS
10000.,0.4833
*VISCOELASTIC,TIME=PRONY
0.901001,0.0,0.99
*TRS,DEFINITION=WLF
70.,4.92,215.
*EXPANSION
1.E-5, 
*BOUNDARY
BOT,2
TOP,2
LEFT,1
*INITIAL CONDITIONS,TYPE=TEMPERATURE,FILE=viscoslabthermload_heat,STEP=1,INC=0
*RESTART,WRITE,FREQUENCY=1
*STEP,INC=50
*VISCO,CETOL=2.E-3
5.E-4,6.
*TEMPERATURE,FILE=viscoslabthermload_heat,ESTEP=2
*OUTPUT,FIELD,VARIABLE=PRESELECT,FREQUENCY=1
*OUTPUT,FIELD,OP=ADD,FREQUENCY=1
*ELEMENT OUTPUT
E,
*OUTPUT,HISTORY,FREQUENCY=1
*ELEMENT OUTPUT,ELSET=LEFT
S,E
*ELEMENT OUTPUT,ELSET=RIGHT
S,E
*EL PRINT,FREQUENCY=10
COORD, EE11, E11, CE11, S22, TEMP
*NODE PRINT,FREQUENCY=0
*FILE FORMAT,ZERO INCREMENT
*END STEP