*HEADING
:  LARGE STRAIN VISCO - VOLUMETRIC RELAXATION CPS4
*RESTART,WRITE
*NODE,nset=nall
1,
2,10.
3,10.,5.,
4,0.,5.,
*ELEMENT,TYPE=CPS4,ELSET=ONE
1,1,2,3,4
*SOLID SECTION,ELSET=ONE,MATERIAL=SIMPLE
1., 
*surface,NAME=SIDE
1,S3
*MATERIAL,NAME=SIMPLE
** *HYPERELASTIC,N=1,MODULI=INSTANTANEOUS
*HYPERELASTIC,N=1,MODULI=long term
 8.,2.,0.2
** 4.,1.,0.2  --> shouldn't change G if g's=0
**8.,2.,0.1
** *VISCOELASTIC,TIME=PRONY
** *VISCOELASTIC,FREQUENCY=PRONY,datat2f=yes
** 0.,0.5,3.
*VISCOELASTIC,FREQUENCY=tabular
*include,input=mvhcdo2zzz.inp
*density
1.0,
*BOUNDARY
1,1,2
2,2,2
4,1,1
*STEP,INC=500
*STEADY STATE DYNAMIC,DIRECT
.1,100.,5
*CLOAD
3,1,-28.
2,1,-28.
4,2,-14.
3,2,-14.
*NODE PRINT,NSET=nall
U,
RF,
PU,
PRF,
*EL PRINT
S,
E,
*ENERGY PRINT
*NODE FILE,NSET=nall
U,
RF,
PU,
PRF,
*EL FILE
S,
E,
*ENERGY FILE
*OUTPUT,FIELD
*NODE OUTPUT,NSET=nall
U,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=nall
U,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*OUTPUT,HISTORY
*NODE OUTPUT, Nset=nall
RF,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*OUTPUT,HISTORY
*ELEMENT OUTPUT, Elset=ONE
S,
E,
*OUTPUT,HISTORY
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP