*HEADING
TWO DIMENSIONAL CYCLIC DEMAND WETTABILITY, WITH GEL
*** UNITS: M, SEC, NEWTON
*PREPRINT,MODEL=YES
*NODE,NSET=ALLN
*NSET,NSET=BOT,GEN
*NSET,NSET=TOP,GEN
*NSET,NSET=OUTN
*NSET,NSET=QA_TEST_NODE
*ELEMENT,TYPE=CPE4RP,ELSET=BLOCK
*ELGEN,ELSET=BLOCK
*ELSET,ELSET=OUTE
*ELSET,ELSET=QA_TEST_ELEM
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
.02,
*MATERIAL,NAME=CORE
*PERMEABILITY,SPECIFIC=10000.
3.7E-4,
*SORPTION
*SORPTION,TYPE=EXSORPTION
*elastic
100000.0, 0.0
*POROUS BULK MODULI
*GEL
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,.05
*NSET,NSET=PORN,GENERATE
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,-10000.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,5.
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=100,AMPLITUDE=STEP
*SOILS,CONSOLIDATION,UTOL=10000.
*BOUNDARY
ALLN,PINNED
1,8,,-100.
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
 .001,,,1.E-6
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=10,NSET=PORN
POR,RVF,RVT
*EL PRINT,FREQUENCY=10
SAT,POR,VOIDR,GELVR
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,FIELD
*NODE OUTPUT
POR,U
*ELEMENT OUTPUT
POR,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=OUTN
 POR,RVF,RVT
*ELEMENT OUTPUT,ELSET=OUTE
 SAT,POR,VOIDR,GELVR
*OUTPUT,HISTORY,FREQ=5
*NODE OUTPUT,NSET=QA_TEST_NODE
 POR,RVF,RVT
*ELEMENT OUTPUT,ELSET=QA_TEST_ELEM
 SAT,POR,VOIDR,GELVR
*END STEP
*STEP,INC=100,AMPLITUDE=RAMP
*SOILS,CONSOLIDATION,UTOL=10000.
*BOUNDARY
ALLN,PINNED
1,8,,-10000.
*END STEP
*STEP,INC=100,AMPLITUDE=RAMP
*SOILS,CONSOLIDATION,UTOL=10000.
*BOUNDARY
ALLN,PINNED
1,8,,-100.
*END STEP