*HEADING
  1-d consolidation of a column of soil
*preprint,echo=yes,model=yes,history=yes
**
*NODE
  1,  0.0,   0.
  2,  0.5,   0.
  3,  1.0,   0.
401,  0.0,  10.
402,  0.5,  10.
403,  1.0,  10.
*********************
*NGEN, NSET=NALL
  1,  401, 10
  2,  402, 20
  3,  403, 10
*NSET, NSET=NALLC, GENERATE
  1,  401, 20
  3,  403, 20
**---------------------
*NSET, NSET=BASE
 1, 2, 3
*NSET, NSET=BASEC
 1, 3
**---------------------
*NSET, NSET=TOP
 401, 402, 403
*NSET, NSET=TOPC
 401, 403
**---------------------
*NSET, NSET=VERBND, GENERATE
   1, 401, 10
   3, 403, 10
**----------------------------------------
*ELEMENT, TYPE=CPE8P, ELSET=ONE
 1,  1,  3,  23,  21,  2,  13,  22,  11
*ELGEN, ELSET=SOIL
 1, 20, 20
*ELSET, ELSET=TOP
  20
**----------------------------------------
*SOLID SECTION, MATERIAL=A1, ELSET=SOIL
*MATERIAL, NAME=A1
*POROUS ELASTIC
** kappa v
  0.1,  0.3
**CLAY PLASTICITY, DEPENDENCIES=1
** Lam  M    ao   beta  K  T  f1
**  0.5,  1.,  29.,  1.,  1., ,  0.
**  0.5,  1.,   5.,  1.,  1., , 10.
*******************************************
*CLAY PLASTICITY, INTERCEPT=3.2773
**  Lam  M   ao   beta  K  T 
   0.5,  1.,   ,  1.,   1., ,
*******************************************
*PERMEABILITY, SPECIFIC=10
 0.11574E-7,
**
*DENSITY
  0.88
*INITIAL CONDITION, TYPE=RATIO
NALL, 1.277, 0., 2.157, 10.
*INITIAL CONDITION, TYPE=STRESS, GEOSTATIC
SOIL, -60., 0., -10., 10., 0.6, 0.6
*INITIAL CONDITION, TYPE=SOLUTION, USER
*INITIAL CONDITION, TYPE=FIELD
NALL
*BOUNDARY
BASE, 1,2
VERBND, 1
 401, 8
 403, 8
**--------------------------------
*RESTART, WRITE, FREQUENCY=20
*STEP,NLGEOM
 ....   GEOSTATIC STEP
*GEOSTATIC
*FIELD,VARIABLE=1, USER
**NALL
**
** surcharge load
**
*DLOAD
 20, P3, 10.
**
** Gravity load
**--------------------------------------
*DLOAD
SOIL, GRAV,  10., 0., -1., 0.
**--------------------------------------
*EL PRINT, POSITION=CENTROID, ELSET=SOIL
S, 
*EL FILE, ELSET=SOIL
S, 
*NODE PRINT,NSET=NALL
U, POR
*NODE FILE,NSET=NALL
U, POR
*MONITOR, NODE=401, DOF=2
*END STEP
**--------------------------------------
**
**  step 2 : load application + 100 KPA (Total=110)
**
*STEP,NLGEOM, AMPLITUDE=RAMP
**-------------------------------
** SET UP INITIAL PORE PRESSURES 
**-------------------------------
*SOILS,CONSOLIDATION
50.0, 500.0
*DLOAD
 20, P3, 110.
*NODE PRINT, FREQUENCY=5,NSET=NALL
U,RF,POR,RVT
*EL PRINT, FREQUENCY=25, POSITION=CENTROID
 S, E
*EL FILE, FREQUENCY=25, POSITION=CENTROID
 S, E
*NODE PRINT, NSET=TOPC
 U, POR
**
*END STEP
**-------------------------------
** CONSOLIDATION
**-------------------------------
**
*STEP,NLGEOM,INC=200
    consolidation
**
*SOILS,CONSOLIDATION, UTOL=10
 1, 5.0E8,  1.E-5, 
**
*BOUNDARY
TOPC, 8, 8, 0.
**
*NODE FILE, NSET=NALLC, FREQUENCY=20
  U, POR
*NODE PRINT, NSET=NALLC, FREQUENCY=20
  U, POR
*EL PRINT, FREQUENCY=25, POSITION=CENTROID, ELSET=SOIL
 S, E, VOIDR
*EL FILE, FREQUENCY=25, POSITION=CENTROID, ELSET=SOIL
 S, E, VOIDR
*MONITOR, NODE=401, DOF=2
*END STEP