*HEADING
 TRANSIENT INTERNAL PRESSURE ON VISCOEL. CYLINDER -- CPE4I
 VISCOEL. MATER. IN ELASTIC CASE UNDER INTERNAL PRESSURE (FUNC. OF TIME)
 UNITS: SI
*RESTART,WRITE
*NODE,NSET=INNERX
*NODE,NSET=OUTERX
*NODE,NSET=CASEX
*NFILL,NSET=CENTER,BIAS=0.9,TWO STEP
INNERX,OUTERX,20,20
*NFILL,NSET=CCASE
OUTERX,CASEX,1,20
*NCOPY, CHANGE NUMBER=1, OLD SET=CENTER, SHIFT, NEW SET=SIDES
*NCOPY, CHANGE NUMBER=-1, OLD SET=CENTER, SHIFT, NEW SET=SIDES
*NCOPY, CHANGE NUMBER=1, OLD SET=CCASE, SHIFT, NEW SET=SIDCAS
*NCOPY, CHANGE NUMBER=-1, OLD SET=CCASE, SHIFT, NEW SET=SIDCAS
*NSET,NSET=INNER
*NSET,NSET=OUTER
*NSET,NSET=CASE
*NSET,NSET=ALL
CENTER,CCASE,SIDES,SIDCAS
*TRANSFORM,NSET=ALL,TYPE=C
*ELEMENT,TYPE=CPE4I
*ELGEN,ELSET=CYLINDER
*ELEMENT,TYPE=CPE4I,ELSET=CASE
*ELSET,ELSET=ELPRINT
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=VISCO,ORIENTATION=CYL
*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
*ORIENTATION,SYSTEM=CYLINDRICAL,NAME=CYL
*BOUNDARY
SIDES,YSYMM
SIDCAS,YSYMM
*DIAGNOSTICS,NONHYBRID=WARNING
*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
*NODE PRINT, NSET=ALL
U,RF
*EL PRINT,ELSET=ELPRINT
S,
E,
EE,
CE,
*END STEP