*HEADING
 GEOTECHNICAL ANALYSIS WITH C3D8P,C3D8RP,C3D8PH and C3D8RPH ELEMENTS 
 AND DIFFERENT MATERIALS: COMPUTATION OF INITIAL STATE WITH 
 ZERO (SMALL) DISPLACEMENTS
**
**   .................... PART DEFINITION
**
*PART,NAME=PART2
  *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.0, 0.0, 1.0
   6, 1.0, 0.0, 1.0
   7, 1.0, 1.0, 1.0
   8, 0.0, 1.0, 1.0
   **
   9, 2.0, 0.0
  10, 3.0, 0.0
  11, 3.0, 1.0
  12, 2.0, 1.0
  13, 2.0, 0.0, 1.0
  14, 3.0, 0.0, 1.0
  15, 3.0, 1.0, 1.0
  16, 2.0, 1.0, 1.0
  *NSET,NSET=NXZ
   1,2,5,6, 9,10,13,14
  *NSET,NSET=NXY
   1,2,3,4, 9,10,11,12
  *NSET,NSET=NXYTOP
   5,6,7,8, 13,14,15,16
  *NSET,NSET=NYZ
   1,4,5,8, 9,12,13,16
  *ELEMENT,TYPE=C3D8P,ELSET=SOIL
   1, 1,2,3,4, 5,6,7,8
  *ELEMENT,TYPE=C3D8RP,ELSET=SOIL
   2, 9,10,11,12, 13,14,15,16
  *ELSET,ELSET=EHG
   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.0, 0.0, 1.0
   6, 1.0, 0.0, 1.0
   7, 1.0, 1.0, 1.0
   8, 0.0, 1.0, 1.0
   **
   9, 2.0, 0.0
  10, 3.0, 0.0
  11, 3.0, 1.0
  12, 2.0, 1.0
  13, 2.0, 0.0, 1.0
  14, 3.0, 0.0, 1.0
  15, 3.0, 1.0, 1.0
  16, 2.0, 1.0, 1.0
   **
  17, 4.0, 0.0
  18, 5.0, 0.0
  19, 5.0, 1.0
  20, 4.0, 1.0
  21, 4.0, 0.0, 1.0
  22, 5.0, 0.0, 1.0
  23, 5.0, 1.0, 1.0
  24, 4.0, 1.0, 1.0
   **
  25, 6.0, 0.0
  26, 7.0, 0.0
  27, 7.0, 1.0
  28, 6.0, 1.0
  29, 6.0, 0.0, 1.0
  30, 7.0, 0.0, 1.0
  31, 7.0, 1.0, 1.0
  32, 6.0, 1.0, 1.0
  *NSET,NSET=NXZ
   1,2,5,6, 9,10,13,14, 17,18,21,22, 25,26,29,30 
  *NSET,NSET=NXY
   1,2,3,4, 9,10,11,12, 17,18,19,20, 25,26,27,28
  *NSET,NSET=NXYTOP
   5,6,7,8, 13,14,15,16, 21,22,23,24, 29,30,31,32
  *NSET,NSET=NYZ
   1,4,5,8, 9,12,13,16, 17,20,21,24, 25,28,29,32
  *ELEMENT,TYPE=C3D8P,ELSET=SOIL
   1, 1,2,3,4, 5,6,7,8
  *ELEMENT,TYPE=C3D8RP,ELSET=SOIL
   2, 9,10,11,12, 13,14,15,16
  *ELEMENT,TYPE=C3D8PH,ELSET=SOIL
   3, 17,18,19,20, 21,22,23,24
  *ELEMENT,TYPE=C3D8RPH,ELSET=SOIL
   4, 25,26,27,28, 29,30,31,32
  *ELSET,ELSET=EHG
   2,4
*END PART
**
**   .................... ASSEMBLY
**
*ASSEMBLY,NAME=GEOASSEMBLY
**
**   ... INSTANCES
**
*INSTANCE,NAME=IELA,PART=PART4
 0,0,0
 *SOLID SECTION,MATERIAL=MATELA, ELSET=SOIL, CONTROLS=SC
*END INSTANCE
*INSTANCE,NAME=IPOR,PART=PART2
 0,3,0
 *SOLID SECTION,MATERIAL=MATPOR, ELSET=SOIL, CONTROLS=SC
*END INSTANCE
*END ASSEMBLY
**
**   .................... SETS
**
*NSET,NSET=NXY
 IELA.NXY,IPOR.NXY
*NSET,NSET=NXYTOP
 IELA.NXYTOP,IPOR.NXYTOP
*NSET,NSET=NXZ
 IELA.NXZ,IPOR.NXZ
*NSET,NSET=NYZ
 IELA.NYZ,IPOR.NYZ
*NSET,NSET=SOIL
 IELA.SOIL,IPOR.SOIL
*ELSET,ELSET=SOIL
 IELA.SOIL,IPOR.SOIL
**
**   .................... 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.
***INITIAL CONDITIONS,TYPE=STRESS
** SOIL,-10.,-10.,-10.
**
*SECTION CONTROLS, NAME=SC, HOURGLASS=ENHANCED
**
**   .................... ZERO BOUNDARY CONDITIONS
**
*BOUNDARY
 NXY, 3,3
 NXZ, 2,2
 NYZ, 1,1
**
**   .................... STEP
**
*STEP,UNSYM=YES
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC,UTOL
 0.5,1.0
*DLOAD
 SOIL,P4,100.
*BOUNDARY
 NXY, 8,8
 NXYTOP, 8,8, 100.0
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
 S,E,PE
*NODE OUTPUT
 U,POR
*END STEP