*HEADING
  1-D Consolidation under surcharge
  Modified Cam-Clay Model
  Units :  Length - m  Force - kN  Time - s
*ELEMENT, TYPE=CPE8P, ELSET=ONE
 1,  1,  3,  23,  21,  2,  13,  22,  11
*ELGEN, ELSET=SOIL
 1, 20, 20
*ELSET, ELSET=TOP
  20
*ELEMENT, TYPE=CPE8, ELSET=MAST
 21, 401,  403,  423, 421, 402, 413, 422, 411
*ELGEN, ELSET=SURCH
21, 10, 20
*********************
*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.
601,  0.0,  15.
602,  0.5,  15.
603,  1.0,  15.
**********************
*NGEN, NSET=NALL
  1,  601, 10
  2,  602, 20
  3,  603, 10
*NSET, NSET=NALLC, GENERATE
  1,  401, 10
  3,  403, 10
**---------------------
*NSET, NSET=BASE
 1, 2, 3
**---------------------
*NSET, NSET=TOP
 401, 402, 403
*NSET, NSET=TOPC
 401, 403
**----------------------------------------
*SOLID SECTION, MATERIAL=A1, ELSET=SOIL
*SOLID SECTION, MATERIAL=EMBANK, ELSET=SURCH
*MATERIAL, NAME=A1
*POROUS ELASTIC
  0.1,  0.3
***CLAY PLASTICITY, DEPENDENCIES=1
**  0.5,  1.,  29.,  1.,  1., ,  0.
**  0.5,  1.,   5.,  1.,  1., , 10.
*CLAY PLASTICITY, INTERCEPT=3.27726
  0.5,  1.,  ,  1.,  1.,,
*PERMEABILITY, SPECIFIC=10
 0.11574E-7
**=========================================
**       density  0.939  at y =  0. (base)
**                0.829  ay y = 10. (top)
**=========================================
*DENSITY, DEPENDENCIES=1
  0.939, 0.,  0.
  0.829, 0., 10.
** 
*USER OUTPUT VARIABLES
  10
**
*MATERIAL, NAME=EMBANK
*ELASTIC
  1.E5,  0.25
*DENSITY
  2.0
*INITIAL CONDITION, TYPE=STRESS, GEOSTATIC, USER
*INITIAL CONDITION, TYPE=RATIO, USER
**
**INITIAL CONDITIONS, TYPE=FIELD, USER
**
**INITIAL CONDITION, TYPE=STRESS, GEOSTATIC
**SOIL, -60., 0., -10., 10., 0.6, 0.6
*INITIAL CONDITION, TYPE=PORE PRESSURE
NALLC, 100., 0.,  0., 10.
**INITIAL CONDITION, TYPE=RATIO
**NALLC, 1.277, 0.,  2.157, 10.
**
*BOUNDARY
BASE, 1,2
NALL, 1
 401, 8
 403, 8
**--------------------------------
*RESTART, WRITE, FREQUENCY=20
*STEP,NLGEOM
 ....   GEOSTATIC STEP
*GEOSTATIC
*FIELD,VARIABLE=1, USER
NALL
*DLOAD
 20, P3, 10.
*DLOAD
SOIL, GRAV,  10., 0., -1., 0.
*MODEL CHANGE, REMOVE
SURCH
**-------------------------------
*EL PRINT, 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, INC=100
**-------------------------------
 SET UP INITIAL PORE PRESSURES 
**-------------------------------
*SOILS,CONSOLIDATION, UTOL=100
5.0, 5000.0, 5.E-5, 
*DLOAD
 30, P3, 10.
*MODEL CHANGE, ADD
SURCH
*DLOAD
SURCH, GRAV, 10,  0, -1, 0
**
*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
**-------------------------------
** step 3 - CONSOLIDATION
**-------------------------------
**
*STEP,NLGEOM, AMPLITUDE=RAMP, INC=200
    consolidation
**
*SOILS,CONSOLIDATION, UTOL=10
 1.E-8, 5.0E10,
**
*CONTROLS, PARAMETERS=FIELD, FIELD=DISPLACEMENT
  ,0.7,
*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