*HEADING
  This input file is modified from mgrooo2hut.inp
  Multiple Load case test for Porous metal plasticity, element type - continuum-CAX4
: UNIAXIAL TENSION, DISPL CONTROL, 2D, CAX4
*NODE
 1, 0. , 0.
 2, 1. , 0.
 3, 1. , 1.
 4, 0. , 1.
*ELEMENT,TYPE=CAX4,ELSET=ONE
    1,    1, 2, 3, 4
*BOUNDARY
 1 , 2 , , 0.
 2 , 2 , , 0.
*NSET,NSET=TOP
 3 , 4
*SOLID SECTION,MATERIAL=BLOCK,ELSET=ONE
*MATERIAL,NAME=BLOCK
*ELASTIC
300., .3
*PLASTIC
 1.0000000, 0.000
 1.7411011, 0.050
 2.7276924, 0.500
 2.9950454, 0.800
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.95
**q1, q2 , q3 [above line D_0 = 1 - f_0 = 0.95]
 1. , 1. , 1.
*STEP,INC=400,NLGEOM
*STATIC
 .01 , 1. , , .01
*BOUNDARY
TOP , 2 , 2 , 0.5
*EL PRINT, FREQUENCY=10
S, 
E, 
PE, 
RD,VVF,VVFG
*NODE PRINT,FREQUENCY=10
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,
PE,
E,
RD,VVF,VVFG
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
U,
RF,
*END STEP


*STEP,PERTURBATION
*STATIC
 1,1
*LOAD CASE,name=PULL
*BOUNDARY
TOP , 2 , 2 , 0.5
*END LOAD CASE
*LOAD CASE,name=CLOAD
*CLOAD
TOP,2,20000
*END LOAD CASE
*EL PRINT, FREQUENCY=10
S, 
E, 
PE, 
RD,VVF,VVFG
*NODE PRINT,FREQUENCY=10
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,
PE,
E,
RD,VVF,VVFG
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
U,
RF,
*END STEP