*HEADING
: RATE DEPENDENT POROUS 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
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 10.e3, 0.25
 alle, 200.0e3, 0.3
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
distElastic,
*PLASTIC
200.,0.
220.,.0009
220.,.0029
*RATE DEPENDENT,TYPE=POWER LAW
40.,5.
*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
*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
*EL PRINT
S, 
SINV, 
E, 
VVF, 
*NODE PRINT
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,VVF
*END STEP
*STEP,INC=10
STEP 3 - 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 4 - LINEAR PERTURBATION
*STATIC,DIRECT
.1E-3,.1E-3
*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
*EL PRINT
S, 
SINV, 
E, 
VVF, 
*NODE PRINT
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,VVF
*END STEP