*HEADING
*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>
*EXPANSION, USER, TYPE=ISO, PROPERTIES=3
 <alpha>, <Ks>, <cte>
**-----------------------------------------
***PORE FLUID PRESSURE, USER, PROPERTIES=0
*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
*AMPLITUDE, NAME=AMP-1
0.0, 0.0, 500.0, 1.0, 1000.0, 1.0
*INCLUDE, INPUT=amp_por_elas.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