*HEADING
 TRANSIENT INTERNAL PRESSURE ON VISCOEL. CYLINDER -- CAX4I
 VISCOEL. MATER. IN ELASTIC CASE UNDER INTERNAL PRESSURE (FUNC. OF TIME)
 UNITS: SI
 DYNAMIC ANALYSIS
*NODE,NSET=INNER
*NODE,NSET=OUTER
*NODE,NSET=CASE
*NFILL,NSET=CYLINDER,BIAS=0.9,TWO STEP
INNER,OUTER,20,20
*NFILL,NSET=CASE
OUTER,CASE,1,20
*NSET,NSET=ALL
CYLINDER,CASE
*ELEMENT,TYPE=CAX4I
*ELGEN,ELSET=CYLINDER
*ELEMENT,TYPE=CAX4I,ELSET=CASE
*ELSET,ELSET=ELPRINT
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=VISCO
*SOLID SECTION,ELSET=CASE,MATERIAL=ELASTIC
*MATERIAL,NAME=ELASTIC
*DENSITY
7.8E-6,
*ELASTIC,TYPE=ISOTROPIC
 200000, 3.00E-01
*MATERIAL,NAME=VISCO
*DENSITY
1.755E-6,
*HYPERELASTIC,MODULI=INSTANTANEOUS
 275.247,0.,7.E-7
*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
*DIAGNOSTICS,NONHYBRID=WARNING
*STEP,NLGEOM,INC=500
*DYNAMIC,HAFTOL=1000.
0.0001,0.05,,0.002
*DLOAD
1,P4NU
*NODE FILE,NSET=INNER,FREQUENCY=3
U,RF
*EL FILE,ELSET=ELPRINT,FREQUENCY=3
S,E,CE
LOADS,
*NODE PRINT, NSET=INNER
U,RF
*EL PRINT,ELSET=ELPRINT
S,E
CE,
LOADS,
*END STEP