*HEADING
ONE DIMENSIONAL WICKING PROBLEM, COUPLED
*** UNITS: M, TON, SEC, KN
*NODE,NSET=ALLN
*NGEN,NSET=BOTB
*NGEN,NSET=TOPB
*NFILL,NSET=ALLN
 BOTB,TOPB,10,10
*NGEN,NSET=BOTF
*NGEN,NSET=TOPF
*NFILL,NSET=ALLN
 BOTF,TOPF,10,10
*NSET,NSET=LHS,GENERATE
*NSET,NSET=POREP,GENERATE
*ELEMENT,TYPE=C3D6Ph,ELSET=BLOCK
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
*MATERIAL,NAME=CORE
*ELASTIC
 50.,0.
*DENSITY
 .1,
*POROUS BULK MODULI
*PERMEABILITY,SPECIFIC=10.
 3.7E-4,
*SORPTION
*SORPTION,TYPE=EXSORPTION
*INITIAL CONDITIONS,TYPE=SATURATION
 ALLN,.05
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
 POREP, -22.0, 1.0, -12.0, 0.0
*INITIAL CONDITIONS,TYPE=RATIO
 ALLN,5.
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
 BLOCK,-1.1,1.,-2.016666667,0.,0.,0.
*EQUATION
*BOUNDARY
 ALLN,1
 ALLN,2
 BOTB,3
 BOTF,3
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=1
*GEOSTATIC
 1.E-6,1.E-6
*BOUNDARY
 1,8,,-12.
*DLOAD
 BLOCK,GRAV,10.,0.,0.,-1.
*NODE PRINT,FREQUENCY=5,NSET=LHS
 U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,POSITION=AVERAGED AT NODES
 S,
 E,
 SAT,POR,VOIDR
*OUTPUT, FIELD, frequency=10
*ELEMENT OUTPUT
 S,POR,SAT
*NODE OUTPUT
 U,POR
*END STEP
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=20.
*BOUNDARY
 1,8,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT
,1.,
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
,1.,
*END STEP