*HEADING
 GEOTECHNICAL ANALYSIS WITH CPE6MP ELEMENTS AND DIFFERENT MATERIALS: 
 COMPUTATION OF INITIAL STATE WITH ZERO (SMALL) DISPLACEMENTS
**
**   .................... PART DEFINITION
**
*PART,NAME=PART2
  *NODE,NSET=SOIL
   1, 
   2, 1.0
   3, 1.0,1.0
   4, 0.0,1.0
   5, 0.5,0.0
   6, 1.0,0.5
   7, 0.5,1.0
   8, 0.0,0.5
   9, 0.5,0.5
  *NSET,NSET=TOP
   3,4,7
  *NSET,NSET=BOTTOM
   1,2,5
  *NSET,NSET=LEFT
   1,4,8
  *ELEMENT,TYPE=CAX6MP,ELSET=SOIL
   1, 1,2,4,5,9,8
   2, 2,3,4,6,7,9
  *ELSET,ELSET=EDLOAD
   2,
*END PART
*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, 0.5,0.0
   6, 1.0,0.5
   7, 0.5,1.0
   8, 0.0,0.5
   9, 0.5,0.5
   **
   10, 2.0,0.0
   20, 3.0,0.0
   30, 3.0,1.0
   40, 2.0,1.0
   50, 2.5,0.0
   60, 3.0,0.5
   70, 2.5,1.0
   80, 2.0,0.5
   90, 2.5,0.5
  *NSET,NSET=TOP
   3,4,7, 30,40,70
  *NSET,NSET=BOTTOM
   1,2,5, 10,20,50
  *NSET,NSET=LEFT
   1,4,8, 10,40,80
  *ELEMENT,TYPE=CAX6MP,ELSET=SOIL
   1, 1,2,4,5,9,8
   2, 2,3,4,6,7,9
  *ELEMENT,TYPE=CAX6MPH,ELSET=SOIL
   3, 10,20,40,50,90,80
   4, 20,30,40,60,70,90
  *ELSET,ELSET=EDLOAD
   2,4
*END PART
**
**   .................... ASSEMBLY
**
*ASSEMBLY,NAME=GEOASSEMBLY
**
**   ... INSTANCES
**
*INSTANCE,NAME=IELA,PART=PART4
 0,0,0
 *SOLID SECTION,MATERIAL=MATELA, ELSET=SOIL
*END INSTANCE
*INSTANCE,NAME=IPOR,PART=PART2
 0,3,0
 *SOLID SECTION,MATERIAL=MATPOR, ELSET=SOIL
*END INSTANCE
*END ASSEMBLY
**
**   .................... SETS
**
*NSET,NSET=TOP
 IELA.TOP,IPOR.TOP
*NSET,NSET=BOTTOM
 IELA.BOTTOM, IPOR.BOTTOM
*NSET,NSET=LEFT
 IELA.LEFT,IPOR.LEFT
*NSET,NSET=SOIL
 IELA.SOIL,IPOR.SOIL
*ELSET,ELSET=SOIL
 IELA.SOIL,IPOR.SOIL
*ELSET,ELSET=EDLOAD
 IELA.EDLOAD,IPOR.EDLOAD
**
**   .................... MATERIALS
**
**   ... Elastic 
**
*MATERIAL,NAME=MATELA
*ELASTIC
 1000.0,0.3
*PERMEABILITY,SPECIFIC=10.0
 1.0,0.0
**
**   ... Porous Elastic
**
*MATERIAL,NAME=MATPOR
*POROUS ELASTIC
 .026,.3,100.0
*PERMEABILITY,SPECIFIC=10.0
 1.0,0.0
**
**   .................... INITIAL CONDITIONS
**
*INITIAL CONDITIONS,TYPE=RATIO
 SOIL,1.08,0.,1.08,1.
**
**   .................... ZERO BOUNDARY CONDITIONS
**
*BOUNDARY
 LEFT,   1,1
 BOTTOM, 2,2
**
**   .................... STEP-1
**
*STEP,UNSYM=YES,INC=1000
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC,UTOL
 1.0,1.0
*DLOAD
 EDLOAD,P1,100.
*BOUNDARY
 BOTTOM, 8,8
 TOP,    8,8, 100.0
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
 S,E,PE
*NODE OUTPUT
 U,POR
*END STEP
**
**   .................... STEP-2
**
*STEP,UNSYM=YES,INC=1000
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC
 1.0,1.0
*DLOAD
 EDLOAD,P1,100.
*BOUNDARY
 BOTTOM, 8,8
 TOP,    8,8, 100.0
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
 S,E,PE
*NODE OUTPUT
 U,POR
*END STEP