*HEADING
 GEOTECHNICAL ANALYSIS WITH C3D20P,C3D20RP,C3D20PH and C3D20RPH 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, 0.5,0.0
  10, 1.0,0.5	
  11, 0.5,1.0
  12, 0.0,0.5
   **
  13, 0.5,0.0,1.0
  14, 1.0,0.5,1.0
  15, 0.5,1.0,1.0
  16, 0.0,0.5,1.0
   **
   17, 0.0,0.0,0.5
   18, 1.0,0.0,0.5
   19, 1.0,1.0,0.5
   20, 0.0,1.0,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
   **
   109, 2.5,0.0
   110, 3.0,0.5	
   111, 2.5,1.0
   112, 2.0,0.5
   **
   113, 2.5,0.0,1.0
   114, 3.0,0.5,1.0
   115, 2.5,1.0,1.0
   116, 2.0,0.5,1.0
   **
   117, 2.0,0.0,0.5
   118, 3.0,0.0,0.5
   119, 3.0,1.0,0.5
   120, 2.0,1.0,0.5
  *NSET,NSET=NXZ
   1,2,5,6,9,13,17,18, 101,102,105,106,109,113,117,118
  *NSET,NSET=NXY
   1,2,3,4,9,10,11,12, 101,102,103,104,109,110,111,112
  *NSET,NSET=NXYTOP
   5,6,7,8,13,14,15,16, 105,106,107,108,113,114,115,116
  *NSET,NSET=NYZ
   1,4,5,8, 12,16,17,20, 101,104,105,108, 112,116,117,120
  *ELEMENT,TYPE=C3D20P,ELSET=SOIL
   1, 1,2,3,4, 5,6,7,8, 9,10,11,12, 13,14,15,
      16,17,18,19,20
  *ELEMENT,TYPE=C3D20RP,ELSET=SOIL
   2, 101,102,103,104, 105,106,107,108, 109,110,111,112, 113,114,115,
      116,117,118,119,120
*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, 0.5,0.0
  10, 1.0,0.5	
  11, 0.5,1.0
  12, 0.0,0.5
   **
  13, 0.5,0.0,1.0
  14, 1.0,0.5,1.0
  15, 0.5,1.0,1.0
  16, 0.0,0.5,1.0
   **
   17, 0.0,0.0,0.5
   18, 1.0,0.0,0.5
   19, 1.0,1.0,0.5
   20, 0.0,1.0,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
   **
   109, 2.5,0.0
   110, 3.0,0.5	
   111, 2.5,1.0
   112, 2.0,0.5
   **
   113, 2.5,0.0,1.0
   114, 3.0,0.5,1.0
   115, 2.5,1.0,1.0
   116, 2.0,0.5,1.0
   **
   117, 2.0,0.0,0.5
   118, 3.0,0.0,0.5
   119, 3.0,1.0,0.5
   120, 2.0,1.0,0.5
   ****************
   201, 4.0,0.0
   202, 5.0,0.0
   203, 5.0,1.0
   204, 4.0,1.0
   **
   205, 4.0,0.0,1.0
   206, 5.0,0.0,1.0
   207, 5.0,1.0,1.0
   208, 4.0,1.0,1.0
   **
   209, 4.5,0.0
   210, 5.0,0.5	
   211, 4.5,1.0
   212, 4.0,0.5
   **
   213, 4.5,0.0,1.0
   214, 5.0,0.5,1.0
   215, 4.5,1.0,1.0
   216, 4.0,0.5,1.0
   **
   217, 4.0,0.0,0.5
   218, 5.0,0.0,0.5
   219, 5.0,1.0,0.5
   220, 4.0,1.0,0.5
   ********
   ****************
   301, 6.0,0.0
   302, 7.0,0.0
   303, 7.0,1.0
   304, 6.0,1.0
   **
   305, 6.0,0.0,1.0
   306, 7.0,0.0,1.0
   307, 7.0,1.0,1.0
   308, 6.0,1.0,1.0
   **
   309, 6.5,0.0
   310, 7.0,0.5	
   311, 6.5,1.0
   312, 6.0,0.5
   **
   313, 6.5,0.0,1.0
   314, 7.0,0.5,1.0
   315, 6.5,1.0,1.0
   316, 6.0,0.5,1.0
   **
   317, 6.0,0.0,0.5
   318, 7.0,0.0,0.5
   319, 7.0,1.0,0.5
   320, 6.0,1.0,0.5
  *NSET,NSET=NXZ
   1,2,5,6,9,13,17,18, 101,102,105,106,109,113,117,118,
   201,202,205,206,209,213,217,218, 301,302,305,306,309,313,317,318
  *NSET,NSET=NXY
   1,2,3,4,9,10,11,12, 101,102,103,104,109,110,111,112,
   201,202,203,204,209,210,211,212, 301,302,303,304,309,310,311,312
  *NSET,NSET=NXYTOP
   5,6,7,8,13,14,15,16, 105,106,107,108,113,114,115,116,
   205,206,207,208,213,214,215,216, 305,306,307,308,313,314,315,316
  *NSET,NSET=NYZ
   1,4,5,8, 12,16,17,20, 101,104,105,108, 112,116,117,120,
   201,204,205,208, 212,216,217,220, 301,304,305,308, 312,316,317,320
  *ELEMENT,TYPE=C3D20P,ELSET=SOIL
   1, 1,2,3,4, 5,6,7,8, 9,10,11,12, 13,14,15,
      16,17,18,19,20
  *ELEMENT,TYPE=C3D20RP,ELSET=SOIL
   2, 101,102,103,104, 105,106,107,108, 109,110,111,112, 113,114,115,
      116,117,118,119,120
  *ELEMENT,TYPE=C3D20PH,ELSET=SOIL
   3, 201,202,203,204, 205,206,207,208, 209,210,211,212, 213,214,215,
      216,217,218,219,220
  *ELEMENT,TYPE=C3D20RPH,ELSET=SOIL
   4, 301,302,303,304, 305,306,307,308, 309,310,311,312, 313,314,315,
      316,317,318,319,320
*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=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.
**
**   .................... ZERO BOUNDARY CONDITIONS
**
*BOUNDARY
 NXY, 3,3
 NXZ, 2,2
 NYZ, 1,1
**
**   .................... STEP
**
*STEP,UNSYM=YES
  GEOSTATIC INITIAL STRESS STATE
*GEOSTATIC,UTOL
*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