*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
Keff = E0/3.0/(1.0-2.0*nu)
alpha = 1.0-Keff/Ks
AlphaExpan = alpha/3.0/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
*SOLID SECTION, ELSET=EALL, MATERIAL=TISSUE
*MATERIAL, NAME=TISSUE
***DENSITY
** <rho>
*ELASTIC
<E0>,<nu>
*POROUS BULK MODULI
 <Ks>,<Kf>
*PERMEABILITY, SPECIFIC=1.0e-6
<ko>
*********************************
** predefined pore fluid pressure
*********************************
*PORE FLUID PRESSURE, FIELD=1
*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_elas.inp
*STEP,inc=500,NLGEOM,AMPLITUDE=RAMP
*STATIC
10, 1000., , 10.
*FIELD, VARIABLE=1, AMPLITUDE=N1_POR
 1, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N3_POR
 3, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N5_POR
 5, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N7_POR
 7, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N9_POR
 9, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N11_POR
 11, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N13_POR
 13, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N15_POR
 15, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N17_POR
 17, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N19_POR
 19, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N21_POR
 21, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N23_POR
 23, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N25_POR
 25, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N27_POR
 27, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N29_POR
 29, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N31_POR
 31, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N201_POR
 201, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N203_POR
 203, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N205_POR
 205, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N207_POR
 207, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N209_POR
 209, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N211_POR
 211, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N213_POR
 213, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N215_POR
 215, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N217_POR
 217, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N219_POR
 219, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N221_POR
 221, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N223_POR
 223, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N225_POR
 225, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N227_POR
 227, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N229_POR
 229, 1.
*FIELD, VARIABLE=1, AMPLITUDE=N231_POR
 231, 1.
*BOUNDARY, AMPLITUDE=AMP-1
TOP,2,2,-0.125
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S, POR
*NODE OUTPUT
 RF
*END STEP