*HEADING
ONE DIMENSIONAL DESATURATION PROBLEM, COUPLED
*** UNITS: M, TON, SEC, KN
*NODE,NSET=ALLN
*NGEN,NSET=BOT
*NGEN,NSET=TOP
*NFILL,NSET=ALLN
BOT,TOP,20,5
*NSET,NSET=LHS,GENERATE
*NSET,NSET=RHS,GENERATE
*ELEMENT,TYPE=CPE8RP,ELSET=BLOCK
*ELGEN,ELSET=BLOCK
*ELSET,ELSET=OUTE
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
*MATERIAL,NAME=CORE
*ELASTIC
1.3E3,0.
*DENSITY
1.5,
*POROUS BULK MODULI
*PERMEABILITY,SPECIFIC=10.
4.5E-6,
*PERMEABILITY,TYPE=SATURATION
.666666,.85
1.,1.
*SORPTION
*SORPTION,TYPE=EXSORPTION
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,1.
*NSET,NSET=PORN,GENERATE
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,0.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,.4235
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
BLOCK,0.,1.,-17.9750615,0.,0.,0.
*EQUATION
*BOUNDARY
ALLN,1
BOT,2
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=1
*GEOSTATIC
1.E-6,1.E-6
*DLOAD
BLOCK,GRAV,10.,0.,-1.,0.
*BOUNDARY
PORN,8,,0.
*NODE PRINT,FREQUENCY=5,NSET=LHS
U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,ELSET=OUTE
S,E
SAT,POR,VOIDR
*NODE FILE,FREQUENCY=10,NSET=LHS
U,RF,POR,RVT
*EL FILE,FREQUENCY=10,ELSET=OUTE
S,E
SAT,POR,VOIDR
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,Frequency=1
*NODE OUTPUT,NSET=LHS
U,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=OUTE
S,E
SAT,POR,VOIDR
*END STEP
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=10.
*BOUNDARY,OP=NEW
1,8,,0.
ALLN,1,,0.
BOT,2,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*END STEP