*HEADING
 predefined pore fluid pressure (via field variable)
*PREPRINT, MODEL=YES
*PARAMETER
E0 = 30000.0
nu = 0.3
G0 = E0/(2.0*(1.0+nu))
K0 = E0/(3.0*(1.0-2.0*nu))
e0 = 0.666
Kf = 2.0E6
Ks = 1.0E15
**Ks = 1.0E9
rhof = 1000.0
rhos = 2650.0
ko = 5e-14
specWeight = 1.0e-6
conduct=ko/specWeight
n0 = e0/(1.0+e0)
MassScaling = 1.e+3
rho = ((1.0-n0)*rhos + n0 * rhof)*MassScaling
alpha = 1.0 -K0/Ks
Qinv = (alpha-1.0)/Ks + (1.0-n0)/Ks + n0/Kf
Qinv0 = alpha/Ks 
Qinv1 = (alpha-1.0)/Ks + 1.0/Kf
C = Qinv/rho 
C0 = Qinv0/rho 
C1 = Qinv1/rho 
expan0 = - 1.0/(3.0*Ks)
cte = 1E-5
**
Keff = K0 + alpha*alpha/Qinv
tscale = sqrt(K0/Keff)
*NODE
1,0.0,2.5
11,0.0,2.0
31,0.0,0.0
201,3.175,2.5
211,3.175,2.0
231,3.175,0.0
*NGEN
1,11,1
11,31,1
201,211,1
211,231,1
*NSET,GENERATE,NSET=AXIS
1,31,1
*NSET,GENERATE,NSET=SIDE
201,231,1
*NFILL,NSET=NALL
AXIS,SIDE,2,100
*NSET,GENERATE,NSET=BOT
31,231,100
*NSET,GENERATE,NSET=TOP
1,201,100
*ELEMENT,TYPE=CAX4R
1,3,203,201,1
*ELGEN,ELSET=EALL
1,15,2,1
*ELSET,ELSET=OUTPUT
1,2
*SOLID SECTION, ELSET=EALL, MATERIAL=TISSUE
*MATERIAL, NAME=TISSUE
*DENSITY
 <rho>
*ELASTIC
 <E0>,<nu>
*CRUSHABLE FOAM, HARDENING=VOLUMETRIC
 1.1, 0.1
*CRUSHABLE FOAM HARDENING
 2200.000,  0.0
 2465.121,  0.1
 2729.380,  0.2
 2990.192,  0.3
 3245.491,  0.4
 3493.453,  0.5
 3732.598,  0.6
 3961.709,  0.7
 4180.111,  0.8
 4387.210,  0.9
 4582.682,  1.0
 4938.430,  1.2
 5248.363,  1.4
 5515.284,  1.6
 5743.073,  1.8
 5935.948,  2.0
 6293.596,  2.5
 6519.947,  3.0
 6833.387,  5.0
 6883.322, 10.0
**-------------------------------------------
**---*EXPANSION, USER, TYPE=ISO, PROPERTIES=3
**--- <alpha>, <Ks>, <cte>
**----------------------------------------
*PORE FLUID PRESSURE, FIELD=2
**----------------------------------------
**-----*USER DEFINED FIELD, PROPERTIES=2
**----- <alpha>, <Ks>
*DEPVAR
4,
1,Porosity, Porosity
2,VOIDR, Void Ratio
3,POR, Pore pressure
4,EVOLR, Volumetric Strain Rate
*INITIAL CONDITIONS,TYPE=SOLUTION
EALL,<n0>,<e0>
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=1
NALL,<n0>
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=2
NALL,0.
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,0.
*BOUNDARY
BOT,YSYMM
AXIS,XSYMM
SIDE,1,1
TOP,11,11
**BOT,8,8
**SIDE,8,8
*AMPLITUDE, NAME=AMP-1
0.0, 0.0, 500.0, 1.0, 1000.0, 1.0
*INCLUDE, INPUT=amp_por_crushablefoam.inp
*STEP
*DYNAMIC, EXPLICIT, scale factor=<tscale>
 ,1000.
** ,50.
*TEMPERATURE
 NALL, 0.
*FIELD, VARIABLE=2, AMPLITUDE=N1_POR
 1, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N3_POR
 3, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N5_POR
 5, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N7_POR
 7, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N9_POR
 9, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N11_POR
 11, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N13_POR
 13, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N15_POR
 15, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N17_POR
 17, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N19_POR
 19, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N21_POR
 21, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N23_POR
 23, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N25_POR
 25, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N27_POR
 27, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N29_POR
 29, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N31_POR
 31, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N201_POR
 201, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N203_POR
 203, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N205_POR
 205, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N207_POR
 207, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N209_POR
 209, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N211_POR
 211, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N213_POR
 213, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N215_POR
 215, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N217_POR
 217, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N219_POR
 219, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N221_POR
 221, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N223_POR
 223, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N225_POR
 225, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N227_POR
 227, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N229_POR
 229, 1.
*FIELD, VARIABLE=2, AMPLITUDE=N231_POR
 231, 1.
*BOUNDARY, AMPLITUDE=AMP-1
TOP,2,2,-0.125
*OUTPUT,FIELD,VAR=PRESELECT
*ELEMENT OUTPUT
SDV, FV, THE, TEMP
*OUTPUT,HISTORY,VAR=PRESELECT
*NODE OUTPUT,NSET=BOT
RF
*ELEMENT OUTPUT,ELSET=OUTPUT
TEMP, SDV1, FV2
*END STEP