*HEADING
 GEOTECHNICAL ANALYSIS WITH CPE4P,CPE4RP,CPE4PH and CPE4RPH ELEMENTS 
 MATERIALS : CLAY PLASTICITY
 COMPUTATION OF INITIAL STATE WITH ZERO (SMALL) DISPLACEMENTS
**
**   .................... PART DEFINITION
**
*PART,NAME=PART4
  *NODE,NSET=SOIL
   1, 0.0, 0.0
   2, 1.0, 0.0
   3, 1.0, 1.0
   4, 0.0, 1.0
   **
   5, 2.0, 0.0
   6, 3.0, 0.0
   7, 3.0, 1.0
   8, 2.0, 1.0
   **
   9, 4.0, 0.0
  10, 5.0, 0.0
  11, 5.0, 1.0
  12, 4.0, 1.0
   **
  13, 6.0, 0.0
  14, 7.0, 0.0
  15, 7.0, 1.0
  16, 6.0, 1.0
  *NSET,NSET=TOP
   3,4,7,8,11,12,15,16
  *NSET,NSET=BOTTOM
   1,2,5,6, 9,10,13,14
  *NSET,NSET=LEFT
   1,4,5,8, 9,12,13,16
  *ELEMENT,TYPE=CPE4P,ELSET=SOIL
   1, 1,2,3,4
  *ELEMENT,TYPE=CPE4RP,ELSET=SOIL
   2, 5,6,7,8
**  *ELEMENT,TYPE=CPE4P,ELSET=SOIL
  *ELEMENT,TYPE=CPE4P,ELSET=SOIL
   3, 9,10,11,12
**  *ELEMENT,TYPE=CPE4RPH,ELSET=SOIL
  *ELEMENT,TYPE=CPE4RP,ELSET=SOIL
   4, 13,14,15,16
  *ELSET,ELSET=EHG
   2,4
*END PART
**
**   .................... ASSEMBLY
**
*ASSEMBLY,NAME=GEOASSEMBLY
  **
  **   ... INSTANCES
  **
  *INSTANCE,NAME=ICLAYPLA,PART=PART4
    *SOLID SECTION,MATERIAL=MATCLAYPLA, ELSET=SOIL, CONTROLS=SC
  *END INSTANCE
  *INSTANCE,NAME=IMC,PART=PART4
    *SOLID SECTION,MATERIAL=MATMC, ELSET=SOIL, CONTROLS=SC
  *END INSTANCE
*END ASSEMBLY
**
**   .................... SETS
**
*NSET,NSET=TOP
 ICLAYPLA.TOP,IMC.TOP
*NSET,NSET=BOTTOM
 ICLAYPLA.BOTTOM,IMC.BOTTOM
*NSET,NSET=LEFT
 ICLAYPLA.LEFT,IMC.LEFT
*NSET,NSET=SOIL
 ICLAYPLA.SOIL,IMC.SOIL
*ELSET,ELSET=SOIL
 ICLAYPLA.SOIL,IMC.SOIL
*ELSET,ELSET=EHG
 ICLAYPLA.EHG,IMC.EHG
**
**   .................... MATERIALS
**
**   ... Clay Plasticity
**
*MATERIAL,NAME=MATCLAYPLA
*elastic
 1000.
*clay plasticity,hardening=tab
,1.0,0.5,1.0
*clay hardening
 200.0 ,  0.000
1000.0,   2.0
*PERMEABILITY,SPECIFIC=10.0
 1.0,0.0
**
**   ... Mohr  Coulomb
**
*MATERIAL,NAME=MATMC
*ELASTIC
 1000.0,.3
*MOHR COULOMB
 40.,35.0
*MOHR COULOMB HARDENING
  20.0, 0.
  30.0, 0.02
  40.0, 0.063333
  50.0, 0.11
  50.0, 1.
*PERMEABILITY,SPECIFIC=10.0
 1.0,0.0
**
*SECTION CONTROLS, NAME=SC, HOURGLASS=ENHANCED
**
**   .................... INITIAL CONDITIONS
**
*INITIAL CONDITIONS,TYPE=RATIO
 SOIL,1.08,0.,1.08,1.
*INITIAL CONDITIONS,TYPE=STRESS
 SOIL,-10.0,-10.0,-10.0
** IMC.SOIL,-30.0,-30.0,-30.0
**
**   .................... ZERO BOUNDARY CONDITIONS
**
*BOUNDARY
 LEFT,   1,1
 BOTTOM, 2,2
**
**   .................... STEP
**
*STEP,UNSYM=YES,NLGEOM=YES,INC=500
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC,UTOL
 0.01,1.0
*DLOAD
 SOIL,P2,150.
*BOUNDARY
 BOTTOM, 8,8
 TOP, 8,8
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT
 S,E,PE,PEEQ
*NODE OUTPUT
 U,POR
*END STEP