*Heading
** Job name: Job-1 Model name: desaturation_cax6mp_deform
**
** PARTS
**
*Part, name=Part-1
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=PART-1-1, part=Part-1
*Node
      1,           0.,           0.
      2,          0.1,           0.
      3,          0.1,           1.
      4,           0.,           1.
      5,          0.1,          0.1
      6,          0.1,          0.2
      7,          0.1,          0.3
      8,          0.1,          0.4
      9,          0.1,          0.5
     10,          0.1,          0.6
     11,          0.1,          0.7
     12,          0.1,          0.8
     13,          0.1,          0.9
     14,           0.,          0.9
     15,           0.,          0.8
     16,           0.,          0.7
     17,           0.,          0.6
     18,           0.,          0.5
     19,           0.,          0.4
     20,           0.,          0.3
     21,           0.,          0.2
     22,           0.,          0.1
     23,         0.05,         0.05
     24,         0.05,           0.
     25,          0.1,         0.05
     26,         0.05,         0.15
     27,         0.05,          0.2
     28,           0.,         0.15
     29,         0.05,          0.3
     30,         0.05,         0.35
     31,           0.,         0.35
     32,         0.05,         0.45
     33,          0.1,         0.45
     34,         0.05,          0.5
     35,         0.05,          0.6
     36,         0.05,         0.65
     37,           0.,         0.65
     38,         0.05,          0.7
     39,         0.05,         0.75
     40,           0.,         0.75
     41,          0.1,         0.75
     42,         0.05,          0.8
     43,          0.1,         0.85
     44,         0.05,         0.85
     45,         0.05,          0.9
     46,           0.,         0.85
     47,         0.05,         0.95
     48,           0.,         0.95
     49,          0.1,         0.95
     50,         0.05,           1.
     51,          0.1,         0.65
     52,         0.05,         0.55
     53,          0.1,         0.55
     54,           0.,         0.55
     55,         0.05,          0.4
     56,           0.,         0.45
     57,          0.1,         0.35
     58,         0.05,         0.25
     59,          0.1,         0.25
     60,           0.,         0.25
     61,         0.05,          0.1
     62,          0.1,         0.15
     63,           0.,         0.05
*Element, type=CAX6MP
 1,  5,  1,  2, 23, 24, 25
 2, 22,  6, 21, 26, 27, 28
 3, 20,  7, 19, 29, 30, 31
 4, 18,  8,  9, 32, 33, 34
 5, 17, 10, 16, 35, 36, 37
 6, 16, 11, 15, 38, 39, 40
 7, 15, 11, 12, 39, 41, 42
 8, 15, 12, 13, 42, 43, 44
 9, 15, 13, 14, 44, 45, 46
10, 14, 13,  4, 45, 47, 48
11, 13,  3,  4, 49, 50, 47
12, 16, 10, 11, 36, 51, 38
13, 17,  9, 10, 52, 53, 35
14, 18,  9, 17, 34, 52, 54
15, 19,  8, 18, 55, 32, 56
16,  8, 19,  7, 55, 30, 57
17,  7, 20,  6, 29, 58, 59
18, 20, 21,  6, 60, 27, 58
19,  6, 22,  5, 26, 61, 62
20, 22,  1,  5, 63, 23, 61
** Region: (Section-1-BLOCK:Picked)
*Elset, elset=_I1, internal, generate
  1,  20,   1
** Section: Section-1-BLOCK
*Solid Section, elset=_I1, material=CORE
1.,
*End Instance
*Nset, nset=ALLN, instance=PART-1-1, generate
  1,  63,   1
*Elset, elset=ALLN, instance=PART-1-1, generate
  1,  20,   1
*Nset, nset=BOT, instance=PART-1-1
  1,  2, 24
*Elset, elset=BOT, instance=PART-1-1
 1,
*Nset, nset=LHS, instance=PART-1-1
  1,  4, 14, 15, 16, 17, 18, 19, 20, 21, 22, 28, 31, 37, 40, 46
 48, 54, 56, 60, 63
*Elset, elset=LHS, instance=PART-1-1
  2,  3,  5,  6,  9, 10, 14, 15, 18, 20
*Nset, nset=OUTE, instance=PART-1-1, generate
  1,  63,   1
*Elset, elset=OUTE, instance=PART-1-1, generate
  1,  20,   1
*Nset, nset=POREP, instance=PART-1-1, generate
  1,  63,   1
*Elset, elset=POREP, instance=PART-1-1, generate
  1,  20,   1
*Nset, nset=RHS, instance=PART-1-1
  2,  3,  5,  6,  7,  8,  9, 10, 11, 12, 13, 25, 33, 41, 43, 49
 51, 53, 57, 59, 62
*Elset, elset=RHS, instance=PART-1-1
  1,  4,  7,  8, 11, 12, 13, 16, 17, 19
*Nset, nset=TOP, instance=PART-1-1
  3,  4, 50
*Elset, elset=TOP, instance=PART-1-1
 11,
*Nset, nset=_Set22, internal, instance=PART-1-1
  1,  2, 24
*Elset, elset=_Set22, internal, instance=PART-1-1
 1,
*Nset, nset=_Set23, internal, instance=PART-1-1, generate
  1,  63,   1
*Elset, elset=_Set23, internal, instance=PART-1-1, generate
  1,  20,   1
*End Assembly
** 
** MATERIALS
** 
*Material, name=CORE
*ELASTIC
1.3E3,0.
*DENSITY
1.5, 
*POROUS BULK MODULI
 ,2.E6
*PERMEABILITY,SPECIFIC=10.
4.5E-6, 
*PERMEABILITY,TYPE=SATURATION
.666666,.85
1.,1.
*SORPTION
-100.,.8
-10.,.85
-9.5,.9
-9.,.922
-8.,.947
-7.,.961
-6.,.973
-4.,.988
-2.,.999
0.,1.
*SORPTION,TYPE=EXSORPTION
-100.,.8
-10.,.85
-9.5,.9
-9.,.922
-8.,.947
-7.,.961
-6.,.973
-4.,.988
-2.,.999
0.,1.
** 
** BOUNDARY CONDITIONS
** 
** Name: Disp-BC-1 Type: Displacement/Rotation
*Boundary
ALLN, 1, 1
** Name: Disp-BC-2 Type: Displacement/Rotation
*Boundary
BOT, 2, 2
*initial conditions, type=saturation
ALLN,1.
*initial conditions,type=pore pressure
porep,0.
*initial conditions,type=stress,geostatic
ALLN,0.,1.,-17.9750615,0.,0.,0.
*initial conditions,type=ratio
ALLN,0.4235
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1
*Geostatic
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-3 Type: Pore pressure
*Boundary
_Set22, 8, 8, 0.
** 
** LOADS
** 
** Name: Load-1   Type: Gravity
*Dload
_Set23, GRAV, 10., 0., -1.
** 
** OUTPUT REQUESTS
** 
*Restart, write, frequency=10
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field,frequency=10
*Node Output
POR,U
*Element Output
S,SAT,VOIDR 
** 
** HISTORY OUTPUT: H-Output-1
** 
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=999999
*End Step
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=20.
1.,1000000.
*BOUNDARY
_set22,8,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT 
,1.,
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
,1.,
*END STEP