*HEADING
 GEOTECHNICAL ANALYSIS WITH COH3D6P ELEMENTS AND DIFFERENT MATERIALS: 
 COMPUTATION OF INITIAL STATE WITH ZERO (SMALL) DISPLACEMENTS
**
**   .................... PART DEFINITION
**
*PART,NAME=PART1
  *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, 0.0,0.0,0.5
  10, 1.0,0.0,0.5
  11, 1.0,1.0,0.5
  12, 0.0,1.0,0.5
  *NSET,NSET=NXZ
   1,2,5,6,9,10
  *NSET,NSET=NXY
   1,2,3,4,
  *NSET,NSET=NXYTOP
   5,6,7,8
  *NSET,NSET=NYZ
   1,4,5,8,9,12
  *ELSET,ELSET=EDLOAD
   1,2,
  *ELEMENT,TYPE=COH3D6P,ELSET=SOIL
   1, 1,2,4, 5,6,8, 9,10,12
   2, 2,3,4, 6,7,8, 10,11,12
*END PART
**
**   .................... ASSEMBLY
**
*ASSEMBLY,NAME=GEOASSEMBLY
**
**   ... INSTANCES
**
*INSTANCE,NAME=IELA,PART=PART1
 0,0,0
 *COHESIVE SECTION, THICKNESS=SPECIFIED, MATERIAL=MATELA, ELSET=SOIL, RESPONSE=CONTINUU
*END INSTANCE
*INSTANCE,NAME=IPOR,PART=PART1
 3,0,0
** *COHESIVE SECTION, THICKNESS=SPECIFIED, MATERIAL=MATPOR, ELSET=SOIL, RESPONSE=CONTINUU
 *COHESIVE SECTION, THICKNESS=SPECIFIED, MATERIAL=MATELA, ELSET=SOIL, RESPONSE=CONTINUU
*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
*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.
*INITIAL CONDITIONS,TYPE=STRESS
 SOIL,-10.,-10.,-10.
**
**   .................... ZERO BOUNDARY CONDITIONS
**
*BOUNDARY
 NXY, 3,3
 NXZ, 2,2
 NYZ, 1,1
**
**   .................... STEP
**
*STEP,UNSYM=YES
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC,UTOL
*DLOAD
 EDLOAD,P2,100.
*BOUNDARY
 NXY, 8,8
 NXYTOP, 8,8, 100.0
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
 S,E,PE,P
*NODE OUTPUT
 U,POR
*END STEP