*HEADING
: Two-layer viscoplasticity, Time Viscous Law, Uniaxial Tension
: and relaxation, C3D8R
: CETOL is simply used as a mechanism to limit total number of
: increments (actual value makes very small difference in results)
*NODE,NSET=ALLN
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*NSET,NSET=PULL
5,6,7,8
*ELEMENT,TYPE=C3D8R,ELSET=ALLE
1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 10.e3, 0.25
 alle, 20.0E6, 0.3
*MATERIAL,NAME=ALLE
*ELASTIC
 distElastic,
*PLASTIC
200.
*VISCOUS,LAW=TIME
1.0e-6,1.,0.,0.2,0.
1.0e-6,1.,0.,0.25,20.
1.0e-6,1.,0.,0.3,40.
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALLN,20.
*STEP, INC=300, AMP=RAMP
*VISCO, CETOL=1.0e-4
0.01,2.
*BOUNDARY
PULL,3,3,2.e-4
*EL FILE,FREQ=5
S,
PS,
VS,
PE,
E,
VE,
EE,
*ENERGY FILE,FREQ=5
*NODE FILE, FREQ=5
U,RF
*EL PRINT
S,
E,
PE,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=ALLE
S
E
*OUTPUT,FIELD,FREQUENCY=0
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY
PULL,3,3,1.
*EL FILE
S,
PS,
VS,
PE,
E,
VE,
EE,
*ENERGY FILE
*END STEP
*STEP, INC=300, AMP=RAMP
*VISCO, CETOL=1.0e-4
0.01,2.
*EL FILE,FREQ=5
S,
PS,
VS,
PE,
E,
VE,
EE,
*ENERGY FILE,FREQ=5
*NODE FILE, FREQ=5
U,RF
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=ALLE
S
E
*END STEP