*HEADING
 GEOTECHNICAL ANALYSIS WITH C3D10MP and C3D10MPH  ELEMENTS AND DIFFERENT 
 MATERIALS:  COMPUTATION OF INITIAL STATE WITH ZERO (SMALL) DISPLACEMENTS
**
**   .................... PART DEFINITION
**
*PART,NAME=PART2
  *NODE,NSET=SOIL
   1, 2.0, 0.0
   2, 3.0, 0.0
   3, 3.0, 1.0
   4, 2.0, 1.0
   **
   5, 2.0, 0.0, 1.0
   6, 3.0, 0.0, 1.0
   7, 3.0, 1.0, 1.0
   8, 2.0, 1.0, 1.0
   **
   24, 2.5, 0.5
   26, 3.0, 0.0, 0.5
   48, 2.0, 1.0, 0.5
   68, 2.5, 0.5, 1.0
   **
   27, 3.0, 0.5, 0.5
   47, 2.5, 1.0, 0.5
   67, 3.0, 0.5, 1.0
   78, 2.5, 1.0, 1.0
   **
   23, 3.0, 0.5
   34, 2.5, 1.0
   37, 3.0, 1.0, 0.5
   **
   25, 2.5, 0.0, 0.5
   45, 2.0, 0.5, 0.5
   56, 2.5, 0.0, 1.0
   58, 2.0, 0.5, 1.0
   ***
   12, 2.5, 0.0
   14, 2.0, 0.5
   15, 2.0, 0.0, 0.5
   **
   28, 2.5, 0.5, 0.5
   *******************************************
   101, 2.0, 0.0
   102, 3.0, 0.0
   103, 3.0, 1.0
   104, 2.0, 1.0
   **
   105, 2.0, 0.0, 1.0
   106, 3.0, 0.0, 1.0
   107, 3.0, 1.0, 1.0
   108, 2.0, 1.0, 1.0
   **
   124, 2.5, 0.5
   126, 3.0, 0.0, 0.5
   148, 2.0, 1.0, 0.5
   168, 2.5, 0.5, 1.0
   **
   127, 3.0, 0.5, 0.5
   147, 2.5, 1.0, 0.5
   167, 3.0, 0.5, 1.0
   178, 2.5, 1.0, 1.0
   **
   123, 3.0, 0.5
   134, 2.5, 1.0
   137, 3.0, 1.0, 0.5
   **
   125, 2.5, 0.0, 0.5
   145, 2.0, 0.5, 0.5
   156, 2.5, 0.0, 1.0
   158, 2.0, 0.5, 1.0
   ***
   112, 2.5, 0.0
   114, 2.0, 0.5
   115, 2.0, 0.0, 0.5
   **
   128, 2.5, 0.5, 0.5
  *NSET,NSET=NXZ
   1,2,5,6, 12,15,25,26,56, 101,102,105,106, 112,115,125,126,156
  *NSET,NSET=NXY
   101,102,103,104, 112,114,123,124,134, 101,102,103,104, 112,114,123,124,134
  *NSET,NSET=NXYTOP
   105,106,107,108, 156,158,167,168,178, 105,106,107,108, 156,158,167,168,178
  *NSET,NSET=NYZ
   1,4,5,8, 14,15,45,48,58, 101,104,105,108, 114,115,145,148,158
  *ELSET,ELSET=EDLOADP2
   1,101
  *ELSET,ELSET=EDLOADP4
   3,103
  *ELEMENT,TYPE=C3D10MP,ELSET=SOIL
   1, 2,3,4,7, 23,34,24, 27,37,47
   2, 2,4,8,7, 24,48,28, 27,47,78
   3, 2,8,6,7, 28,68,26, 27,78,67
   **
   4, 1,2,4,5, 12,24,14, 15,25,45
   5, 8,4,2,5, 48,24,28, 58,45,25
   6, 8,2,6,5, 28,26,68, 58,25,56
  *ELEMENT,TYPE=C3D10MPH,ELSET=SOIL
   101, 102,103,104,107, 123,134,124, 127,137,147
   102, 102,104,108,107, 124,148,128, 127,147,178
   103, 102,108,106,107, 128,168,126, 127,178,167
   **
   104, 101,102,104,105, 112,124,114, 115,125,145
   105, 108,104,102,105, 148,124,128, 158,145,125
   106, 108,102,106,105, 128,126,168, 158,125,156
*END PART
**
**   .................... ASSEMBLY
**
*ASSEMBLY,NAME=GEOASSEMBLY
**
**   ... INSTANCES
**
*INSTANCE,NAME=IELA,PART=PART2
 0,0,0
 *SOLID SECTION,MATERIAL=MATELA, ELSET=SOIL
*END INSTANCE
*INSTANCE,NAME=IPOR,PART=PART2
 3,0,0
** *SOLID SECTION,MATERIAL=MATPOR, ELSET=SOIL
 *SOLID SECTION,MATERIAL=MATELA, ELSET=SOIL
*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=EDLOADP2
 IELA.EDLOADP2,IPOR.EDLOADP2
*ELSET,ELSET=EDLOADP4
 IELA.EDLOADP4,IPOR.EDLOADP4
**
**   .................... 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
 EDLOADP2,P2,100.
 EDLOADP4,P4,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