*HEADING
  this input file is modified from mgrpro3vlp.inp
  Multiple Load case test for hyperelasticity, element type - continuum-C3D8
: RATE POROUS MISES PLASTICITY, UNIAXIAL TENSION TEST, C3D8
WITH LINEAR PERTURBATIONS
*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.
*ELEMENT,TYPE=C3D8,ELSET=ALLE
1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
200.E3,.3
*PLASTIC
200.,0.
220.,.0009
220.,.0029
*PLASTIC,RATE=0.010
238.07, 0.0000
261.88, 0.0009
261.88, 0.0029
*PLASTIC,RATE=0.100
260.34, 0.0000
286.38, 0.0009
286.38, 0.0029
*PLASTIC,RATE=1.000
295.64, 0.0000
325.20, 0.0009
325.20, 0.0029
*PLASTIC,RATE=10.000
351.57, 0.0000
386.73, 0.0009
386.73, 0.0029
*PLASTIC,RATE=100.000
440.22, 0.0000
484.25, 0.0009
484.25, 0.0029
*PLASTIC,RATE=1000.000
580.73, 0.0000
638.80, 0.0009
638.80, 0.0029
*POROUS METAL PLASTICITY,RELATIVE DENSITY=0.95
**q1, q2, q3
1., 1., 1.
*STEP,INC=10
STEP 1 - NONLINEAR
*STATIC,DIRECT
1.E-3,10.E-3
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,.002
5,3,,.002
6,3,,.002
8,3,,.002
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
VVF, 
*NODE PRINT
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,PE,VVF
*END STEP
*STEP,PERTURBATION,INC=1
STEP 2 - LINEAR PERTURBATION
*STATIC,DIRECT
.1E-3,.1E-3
*LOAD CASE,name=PULL
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,.0002
5,3,,.0002
6,3,,.0002
8,3,,.0002
*END LOAD CASE
*LOAD CASE,NAME=PUSH
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,-.0002
5,3,,-.0002
6,3,,-.0002
8,3,,-.0002
*END LOAD CASE
*END STEP
*STEP,INC=10
STEP 4 - REST OF NONLINEAR
*STATIC,DIRECT
1.E-3,10.E-3
*BOUNDARY
7,3,,.004
5,3,,.004
6,3,,.004
8,3,,.004
*END STEP
*STEP,PERTURBATION,INC=1
STEP 5 - LINEAR PERTURBATION
*STATIC,DIRECT
.1E-3,.1E-3
*LOAD CASE,NAME=PULL-1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,.0004
5,3,,.0004
6,3,,.0004
8,3,,.0004
*END LOAD CASE
*LOAD CASE,NAME=PUSH-2
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,-.0004
5,3,,-.0004
6,3,,-.0004
8,3,,-.0004
*END LOAD CASE
*EL PRINT
S, 
SINV, 
E, 
VVF, 
*NODE PRINT
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,VVF
*END STEP