*HEADING TWO DIMENSIONAL CYCLIC DEMAND WETTABILITY, WITH GEL *** UNITS: M, SEC, NEWTON *PREPRINT,MODEL=YES *NODE,NSET=ALLN 1 , 0.0000 , 0.0000 3 , 5.08000E-03, 0.0000 5 , 1.01600E-02, 0.0000 7 , 1.52400E-02, 0.0000 9 , 2.03200E-02, 0.0000 11 , 2.54000E-02, 0.0000 13 , 3.04800E-02, 0.0000 15 , 3.55600E-02, 0.0000 17 , 4.06400E-02, 0.0000 19 , 4.57200E-02, 0.0000 21 , 5.08000E-02, 0.0000 81 , 0.0000 , 5.08000E-03 83 , 5.08000E-03, 5.08000E-03 85 , 1.01600E-02, 5.08000E-03 87 , 1.52400E-02, 5.08000E-03 89 , 2.03200E-02, 5.08000E-03 91 , 2.54000E-02, 5.08000E-03 93 , 3.04800E-02, 5.08000E-03 95 , 3.55600E-02, 5.08000E-03 97 , 4.06400E-02, 5.08000E-03 99 , 4.57200E-02, 5.08000E-03 101 , 5.08000E-02, 5.08000E-03 161 , 0.0000 , 1.01600E-02 163 , 5.08000E-03, 1.01600E-02 165 , 1.01600E-02, 1.01600E-02 167 , 1.52400E-02, 1.01600E-02 169 , 2.03200E-02, 1.01600E-02 171 , 2.54000E-02, 1.01600E-02 173 , 3.04800E-02, 1.01600E-02 175 , 3.55600E-02, 1.01600E-02 177 , 4.06400E-02, 1.01600E-02 179 , 4.57200E-02, 1.01600E-02 181 , 5.08000E-02, 1.01600E-02 241 , 0.0000 , 1.52400E-02 243 , 5.08000E-03, 1.52400E-02 245 , 1.01600E-02, 1.52400E-02 247 , 1.52400E-02, 1.52400E-02 249 , 2.03200E-02, 1.52400E-02 251 , 2.54000E-02, 1.52400E-02 253 , 3.04800E-02, 1.52400E-02 255 , 3.55600E-02, 1.52400E-02 257 , 4.06400E-02, 1.52400E-02 259 , 4.57200E-02, 1.52400E-02 261 , 5.08000E-02, 1.52400E-02 321 , 0.0000 , 2.03200E-02 323 , 5.08000E-03, 2.03200E-02 325 , 1.01600E-02, 2.03200E-02 327 , 1.52400E-02, 2.03200E-02 329 , 2.03200E-02, 2.03200E-02 331 , 2.54000E-02, 2.03200E-02 333 , 3.04800E-02, 2.03200E-02 335 , 3.55600E-02, 2.03200E-02 337 , 4.06400E-02, 2.03200E-02 339 , 4.57200E-02, 2.03200E-02 341 , 5.08000E-02, 2.03200E-02 401 , 0.0000 , 2.54000E-02 403 , 5.08000E-03, 2.54000E-02 405 , 1.01600E-02, 2.54000E-02 407 , 1.52400E-02, 2.54000E-02 409 , 2.03200E-02, 2.54000E-02 411 , 2.54000E-02, 2.54000E-02 413 , 3.04800E-02, 2.54000E-02 415 , 3.55600E-02, 2.54000E-02 417 , 4.06400E-02, 2.54000E-02 419 , 4.57200E-02, 2.54000E-02 421 , 5.08000E-02, 2.54000E-02 481 , 0.0000 , 3.04800E-02 483 , 5.08000E-03, 3.04800E-02 485 , 1.01600E-02, 3.04800E-02 487 , 1.52400E-02, 3.04800E-02 489 , 2.03200E-02, 3.04800E-02 491 , 2.54000E-02, 3.04800E-02 493 , 3.04800E-02, 3.04800E-02 495 , 3.55600E-02, 3.04800E-02 497 , 4.06400E-02, 3.04800E-02 499 , 4.57200E-02, 3.04800E-02 501 , 5.08000E-02, 3.04800E-02 561 , 0.0000 , 3.55600E-02 563 , 5.08000E-03, 3.55600E-02 565 , 1.01600E-02, 3.55600E-02 567 , 1.52400E-02, 3.55600E-02 569 , 2.03200E-02, 3.55600E-02 571 , 2.54000E-02, 3.55600E-02 573 , 3.04800E-02, 3.55600E-02 575 , 3.55600E-02, 3.55600E-02 577 , 4.06400E-02, 3.55600E-02 579 , 4.57200E-02, 3.55600E-02 581 , 5.08000E-02, 3.55600E-02 641 , 0.0000 , 4.06400E-02 643 , 5.08000E-03, 4.06400E-02 645 , 1.01600E-02, 4.06400E-02 647 , 1.52400E-02, 4.06400E-02 649 , 2.03200E-02, 4.06400E-02 651 , 2.54000E-02, 4.06400E-02 653 , 3.04800E-02, 4.06400E-02 655 , 3.55600E-02, 4.06400E-02 657 , 4.06400E-02, 4.06400E-02 659 , 4.57200E-02, 4.06400E-02 661 , 5.08000E-02, 4.06400E-02 721 , 0.0000 , 4.57200E-02 723 , 5.08000E-03, 4.57200E-02 725 , 1.01600E-02, 4.57200E-02 727 , 1.52400E-02, 4.57200E-02 729 , 2.03200E-02, 4.57200E-02 731 , 2.54000E-02, 4.57200E-02 733 , 3.04800E-02, 4.57200E-02 735 , 3.55600E-02, 4.57200E-02 737 , 4.06400E-02, 4.57200E-02 739 , 4.57200E-02, 4.57200E-02 741 , 5.08000E-02, 4.57200E-02 801 , 0.0000 , 5.08000E-02 803 , 5.08000E-03, 5.08000E-02 805 , 1.01600E-02, 5.08000E-02 807 , 1.52400E-02, 5.08000E-02 809 , 2.03200E-02, 5.08000E-02 811 , 2.54000E-02, 5.08000E-02 813 , 3.04800E-02, 5.08000E-02 815 , 3.55600E-02, 5.08000E-02 817 , 4.06400E-02, 5.08000E-02 819 , 4.57200E-02, 5.08000E-02 821 , 5.08000E-02, 5.08000E-02 *NSET,NSET=BOT,GEN 1,21,1 *NSET,NSET=TOP,GEN 801,821,1 *NSET,NSET=OUTN 1,165,329,493,657,821 *NSET,NSET=QA_TEST_NODE 1, 493, 821 *ELEMENT,TYPE=CPE4RP,ELSET=BLOCK 1,1,3,83,81 *ELGEN,ELSET=BLOCK 1,10,2,1,10,80,10 *ELSET,ELSET=OUTE 1,23,45,67,89,100 *ELSET,ELSET=QA_TEST_ELEM 1, 67, 100 *SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE .02, *MATERIAL,NAME=CORE *PERMEABILITY,SPECIFIC=10000. 3.7E-4, *SORPTION -100000.,.04 -10000.,.05 -4500.,.1 -3500.,.18 -2000.,.45 -1000.,.91 0.,1. *SORPTION,TYPE=EXSORPTION -100000.,.09 -10000.,.1 -8000.,.11 -6000.,.18 -4500.,.33 -3000.,.79 -2000.,.91 0.,1. *elastic 100000.0, 0.0 *POROUS BULK MODULI ,100000000. *GEL .0005,.0015,1.E8,500. *INITIAL CONDITIONS,TYPE=SATURATION ALLN,.05 *NSET,NSET=PORN,GENERATE 1,21,2 81,101,2 161,181,2 241,261,2 321,341,2 401,421,2 481,501,2 561,581,2 641,661,2 721,741,2 801,821,2 *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. 1.,600.,,100.0 *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. 1.,600.,,100.0 *BOUNDARY ALLN,PINNED 1,8,,-10000. *END STEP *STEP,INC=100,AMPLITUDE=RAMP *SOILS,CONSOLIDATION,UTOL=10000. 1.,800.,,100.0 *BOUNDARY ALLN,PINNED 1,8,,-100. *END STEP