*HEADING
ONE DIMENSIONAL DESATURATION PROBLEM, UNCOUPLED
*** UNITS: M, TON, SEC, KN
*** A DUMMY ORIENTATION IS INTRODUCED AS A FURTHER TEST
**DEFORMABLE TO RIGID ELEMENT TESTING (CPE8RP)
*NODE,NSET=ALLN
1,0.,0.
3,.1,0.
101,0.,1.
103,.1,1.
*NGEN,NSET=BOT
1,3,1
*NGEN,NSET=TOP
101,103,1
*NFILL,NSET=ALLN
BOT,TOP,20,5
*NSET,NSET=LHS,GENERATE
1,101,5
*NSET,NSET=RHS,GENERATE
3,103,5
*NODE, NSET=NREF
99999, 0.0, 0.0
*ELEMENT,TYPE=CPE8RP,ELSET=BLOCK
1,1,3,13,11,2,8,12,6
*ELGEN,ELSET=BLOCK
1,10,10,1
*ELSET,ELSET=OUTE
1,3,5,7,9
*ORIENTATION,NAME=ROT
0.,1.,0.,-1.,0.,0.
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE,ORIENTATION=ROT
*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.
*RIGID BODY, ELSET=BLOCK, REF NODE=99999
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,1.
*NSET,NSET=PORN,GENERATE
1,101,10
3,103,10
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,0.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,.4235
*EQUATION
2, 
3,8,1.,1,8,-1.
*BOUNDARY
99999, 1, 2
99999, 6
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=10.
20.,50000.
*BOUNDARY
1,8,,0.
*DLOAD
BLOCK,GRAV,10.,0.,-1.,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=5,NSET=LHS
U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,ELSET=OUTE
S,E
SAT,POR,VOIDR
*NODE FILE,FREQUENCY=10,NSET=LHS
U,RF,POR,RVT
*EL FILE,FREQUENCY=10,ELSET=OUTE
S,E
SAT,POR,VOIDR
*OUTPUT,FIELD,VARIABLE=PRESELECT,FREQUENCY=10
*OUTPUT,HISTORY,Frequency=1
*NODE OUTPUT,NSET=LHS
U,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=OUTE
S,E
SAT,POR,VOIDR
*END STEP