*HEADING
ONE DIMENSIONAL COUPLED DEMAND WETTABILITY, WITHOUT GEL, QUADRATIC ELEMENTS
*** UNITS: M, SEC, NEWTON
*NODE,NSET=ALLN
*NGEN,NSET=BOT
*NGEN,NSET=TOP
*NFILL,NSET=ALLN
BOT,TOP,20,5
*NSET,NSET=LHS,GENERATE
*NSET,NSET=RHS,GENERATE
*NSET,NSET=PORN,GENERATE
*ELEMENT,TYPE=CPE8RP,ELSET=BLOCK
*ELGEN,ELSET=BLOCK
*ELSET,ELSET=OUTE
***********************************************
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
.02,
*MATERIAL,NAME=CORE
*ELASTIC
10000.,0.
*POROUS BULK MODULI
*PERMEABILITY,SPECIFIC=10000.
3.7E-4,
*SORPTION
*SORPTION,TYPE=EXSORPTION
***********************************************
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,.05
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,-10000.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,5.
*BOUNDARY
BOT,2
ALLN,1
*RESTART,WRITE,FREQUENCY=10
***********************************************
*Elset,Elset=ODB_ESET
BLOCK,
*Nset, Nset=ODB_NSET
ALLN, BOT, TOP, ALLN
LHS, RHS, PORN,
***QA sets
*nset, nset=qa_test_lhs
 lhs,
*elset, elset=qa_test_oute
 oute,
***********************************************
*STEP
*SOILS,CONSOLIDATION
*DLOAD
10,P3,-500.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=5,NSET=LHS
U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,ELSET=OUTE
S,E
SAT,POR,VOIDR
*OUTPUT,FIELD, Frequency=100
*NODE OUTPUT,NSET=LHS
 U,POR
*ELEMENT OUTPUT,ELSET=OUTE
 S,SAT,VOIDR,GELVR
*OUTPUT,HISTORY,Frequency=5
*NODE OUTPUT,NSET=LHS
 U,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=OUTE
 S22,POR
**
*OUTPUT,FIELD, time interval=250
*NODE OUTPUT,NSET=qa_test_LHS
 U,POR
*ELEMENT OUTPUT,ELSET=qa_test_OUTE
 S,SAT,VOIDR,GELVR
*OUTPUT,HISTORY,time interval=250
*NODE OUTPUT,NSET=qa_test_LHS
 U2,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=qa_test_OUTE
 S22,POR
*END STEP
***********************************************
*STEP,INC=100,NLGEOM
*SOILS,CONSOLIDATION,UTOL=10000.
*BOUNDARY
1,8,,-100.
3,8,,-100.
*NODE PRINT,FREQUENCY=5,NSET=LHS
U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,ELSET=OUTE
S,E
SAT,POR,VOIDR
*OUTPUT,FIELD, Frequency=100
*NODE OUTPUT,NSET=LHS
 U,POR
*ELEMENT OUTPUT,ELSET=OUTE
 S,SAT,VOIDR,GELVR
*OUTPUT,HISTORY,Frequency=5
*NODE OUTPUT,NSET=LHS
 U,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=OUTE
 S22,POR
**
*OUTPUT,FIELD, time interval=250
*NODE OUTPUT,NSET=qa_test_LHS
 U,POR
*ELEMENT OUTPUT,ELSET=qa_test_OUTE
 S,SAT,VOIDR,GELVR
*OUTPUT,HISTORY,time interval=250
*NODE OUTPUT,NSET=qa_test_LHS
 U2,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=qa_test_OUTE
 S22,POR
*END STEP