*HEADING
 TRANSIENT INTERNAL PRESSURE ON VISCOEL. CYLINDER -- CAX8R
 VISCOELASTIC MATERIAL IN ELASTIC CASE
 INTERNAL PRESSURE THAT IS FUNCTION OF TIME
 UNITS: SI
*RESTART,WRITE
*NODE,NSET=INNER
1,10.,0.
2,10.,1.
3,10.,2.
*NODE,NSET=OUTER
401,50.,0.
402,50.,1.
403,50.,2.
*NODE,NSET=CASE
421,50.5,0.
422,50.5,1.
423,50.5,2.
*NFILL,NSET=CYLINDER,BIAS=0.9,TWO STEP
INNER,OUTER,40,10
*NFILL,NSET=CASE
OUTER,CASE,2,10
*NSET,NSET=ALL
CYLINDER,CASE
*ELEMENT,TYPE=CAX8R
1,1,21,23,3,11,22,13,2
*ELGEN,ELSET=CYLINDER
1,20,20,1
*ELEMENT,TYPE=CAX8R,ELSET=CASE
101,401,421,423,403,411,422,413,402
*ELSET,ELSET=ELPRINT
1,21
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=VISCO
*SOLID SECTION,ELSET=CASE,MATERIAL=ELASTIC
*MATERIAL,NAME=ELASTIC
*ELASTIC,TYPE=ISOTROPIC
 200000, 3.00E-01
*MATERIAL,NAME=VISCO
*ELASTIC,MODULI=INSTANTANEOUS
 1651.59, 0.4999
*VISCOELASTIC,TIME=PRONY
 0.1986, ,0.281E-7
 0.1828, ,0.281E-5
 0.1388, ,0.281E-3
 0.2499, ,0.281E-1
 0.1703, ,0.281E1
 0.0593, ,0.281E3
*BOUNDARY
ALL,2
*STEP,INC=50
*VISCO,CETOL=7.E-3
0.01,0.5
*DLOAD
1,P4NU
*NODE FILE,NSET=INNER,FREQUENCY=1
U,RF
*EL FILE,ELSET=ELPRINT,FREQUENCY=1
S,E,CE
LOADS,
*NODE PRINT, NSET=INNER
U,RF
*EL PRINT,ELSET=ELPRINT
S, 
E, 
EE, 
CE, 
LOADS,
*END STEP