*HEADING
*PARAMETER
E0 = 30000.0
nu = 0.3
e0 = 0.666
Kf = 2.0E6
Ks = 1.0E15
rhof = 1000.0
rhos = 2650.0
n0 = e0/(1.0+e0)
rho = (1.0-n0)*rhos + n0 * rhof
ko=5e-14
**
sigy0 = 1000
sigy1 = 1200
*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
*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
*POROUS BULK MODULI
 <Ks>,<Kf>
*PERMEABILITY, SPECIFIC=1.0e-6
<ko>
**----------------------------
*PORE FLUID PRESSURE, FIELD=2
**----------------------------
*INITIAL CONDITIONS,TYPE=RATIO
NALL,<e0>
*BOUNDARY
BOT,YSYMM
AXIS,XSYMM
SIDE,1,1
*AMPLITUDE, NAME=AMP-1
0.0, 0.0, 500.0, 1.0, 1000.0, 1.0
*INCLUDE, INPUT=amp_por_crushablefoam.inp
*STEP,inc=500,NLGEOM,AMPLITUDE=RAMP
*STATIC
10,1000.,,10.
*CONTROLS,PARAMETERS=FIELD
 , , , , ,1.e-2 
*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
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
PEEQ,S,POR
*OUTPUT,HISTORY,VAR=PRESELECT
*NODE OUTPUT,NSET=BOT
RF
*END STEP