*HEADING
ONE DIMENSIONAL WICKING PROBLEM, COUPLED      
*** UNITS: M, TON, SEC, KN
*NODE,NSET=ALLN
1,0.,0.,0.
3,.1,0.,0.
101,0.,0.,1.
103,.1,0.,1.
1001,0.,-1.,0.
1003,.1,-1.,0.
1101,0.,-1.,1.
1103,.1,-1.,1.
*NGEN,NSET=BOTB
1,3,2
*NGEN,NSET=TOPB
101,103,2
*NFILL,NSET=ALLN
BOTB,TOPB,10,10
*NGEN,NSET=BOTF
1001,1003,2
*NGEN,NSET=TOPF
1101,1103,2
*NFILL,NSET=ALLN
BOTF,TOPF,10,10
*NSET,NSET=LHS,GENERATE
1,101,10
1001,1101,10
*NSET,NSET=POREP,GENERATE
1, 101, 10
3, 103, 10
1001, 1101, 10
1003, 1103, 10
*ELEMENT,TYPE=C3D8PH,ELSET=BLOCK
1,1,3,13,11,1001,1003,1013,1011
*ELGEN,ELSET=BLOCK
1,10,10,1
*ELSET,ELSET=OUTE
1,3,5,7,9
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
*MATERIAL,NAME=CORE
*ELASTIC
50.,0.
*DENSITY
.1, 
*POROUS BULK MODULI
 ,2.E6
*PERMEABILITY,SPECIFIC=10.
3.7E-4, 
*SORPTION
-100.,.04
-10.,.05
-4.5,.1
-3.5,.18
-2.,.45
-1.,.91
0.,1.
*SORPTION,TYPE=EXSORPTION
-100.,.09
-10.,.1
-8.,.11
-6.,.18
-4.5,.33
-3.,.79
-2.,.91
0.,1.
*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
2, 
3,8,1.,1,8,-1.
2, 
1003,8,1.,1,8,-1.
2, 
1001,8,1.,1,8,-1.
*BOUNDARY
ALLN,1
ALLN,2
BOTB,3
BOTF,3
*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
*NODE FILE,FREQUENCY=10,NSET=LHS
U,RF,POR,RVT
*EL FILE,FREQUENCY=10,ELSET=OUTE
S,E
SAT,POR,VOIDR
*END STEP
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=20.
1.,1000000.
*BOUNDARY
1,8,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT 
,1.,
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
,1.,
*END STEP