*HEADING
 EXAMPLE - AXISYMMETRIC SIMULATION OF OIL WELL     
 UNITS: F = lbs, L = in., T = days
*REST, W, FREQ=10
*NODE
  1,     0.0, 0.0, 0.0
  9,  3180.0, 0.0, 0.0
 11,  9216.0, 0.0, 0.0
 13, 21504.0, 0.0, 0.0
 15, 46080.0, 0.0, 0.0
 17, 95232.0, 0.0, 0.0
*NGEN
  1,  9, 1
  9, 11, 1
 11, 13, 1
 13, 15, 1
 15, 17, 1
*NSET,NSET=CORNODES, GENERATE
  1, 17, 2
*NSET,NSET=MIDNODES, GENERATE
  2, 16, 2
*NSET,NSET=BASE
 CORNODES,
 MIDNODES,
*NCOPY, CHANGE NUMBER=100, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,7200.

*NCOPY, CHANGE NUMBER=200, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,14400.

*NCOPY, CHANGE NUMBER=300, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,18000.

*NCOPY, CHANGE NUMBER=400, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,21600.

*NCOPY, CHANGE NUMBER=500, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,25200.

*NCOPY, CHANGE NUMBER=600, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,28800.

*NCOPY, CHANGE NUMBER=700, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,31200.

*NCOPY, CHANGE NUMBER=800, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,33600.

*NCOPY, CHANGE NUMBER=900, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,34800.

*NCOPY, CHANGE NUMBER=1000, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,36000.

*NCOPY, CHANGE NUMBER=1100, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,38400.

*NCOPY, CHANGE NUMBER=1200, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,40800.

*NCOPY, CHANGE NUMBER=1300, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,42000.

*NCOPY, CHANGE NUMBER=1400, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,43200.

*NCOPY, CHANGE NUMBER=1500, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,43800.

*NCOPY, CHANGE NUMBER=1600, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,44400.

*NCOPY, CHANGE NUMBER=1700, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,46800.

*NCOPY, CHANGE NUMBER=1800, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,49200.

*NCOPY, CHANGE NUMBER=1900, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,52500.

*NCOPY, CHANGE NUMBER=2000, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,55800.

*NCOPY, CHANGE NUMBER=2100, OLD SET=CORNODES, NEW SET=MID, SHIFT
 0,56700.

*NCOPY, CHANGE NUMBER=2200, OLD SET=CORNODES, NEW SET=CORNER, SHIFT
 0,57600.

*NCOPY, CHANGE NUMBER=200, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,14400.

*NCOPY, CHANGE NUMBER=400, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,21600.

*NCOPY, CHANGE NUMBER=600, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,28800.

*NCOPY, CHANGE NUMBER=800, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,33600.

*NCOPY, CHANGE NUMBER=1000, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,36000.

*NCOPY, CHANGE NUMBER=1200, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,40800.

*NCOPY, CHANGE NUMBER=1400, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,43200.

*NCOPY, CHANGE NUMBER=1600, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,44400.

*NCOPY, CHANGE NUMBER=1800, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,49200.

*NCOPY, CHANGE NUMBER=2000, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,55800.

*NCOPY, CHANGE NUMBER=2200, OLD SET=MIDNODES, NEW SET=MID, SHIFT
 0,57600.

*NSET,NSET=AXIS,GENERATE
 1, 2201, 100
*NSET, NSET=NALL
 MID,
 MIDNODES,
 CORNER,
 CORNODES,
*NSET,NSET=FREESURF,GENERATE
 2201, 2215, 2
**
** REGION WHERE STEAM IS INJECTED
**
*NSET, NSET=TARGET
 601,602,603,604,605,606,607,608,609
 701,702,703,704,705,706,707,708,709
 801,802,803,804,805,806,807,808,809
 901,902,903,904,905,906,907,908,909
 1001,1002,1003,1004,1005,1006,1007,1008,1009
 1101,1102,1103,1104,1105,1106,1107,1108,1109
 1201,1202,1203,1204,1205,1206,1207,1208,1209
 1301,1302,1303,1304,1305,1306,1307,1308,1309
 1401,1402,1403,1404,1405,1406,1407,1408,1409
 1501,1502,1503,1504,1505,1506,1507,1508,1509
 1601,1602,1603,1604,1605,1606,1607,1608,1609
**
** ELEMENTS GENERATION
**
*ELEMENT,TYPE=CAX8RP
 1, 1,3,203,201,2,103,202,101
*ELGEN, ELSET=ESOIL
 1, 7,2,1,11,200,100
*ELEMENT, TYPE=CINAX5R, ELSET=FAR
    8,  215,   15,   17,  217,  115
  108,  415,  215,  217,  417,  315
  208,  615,  415,  417,  617,  515
  308,  815,  615,  617,  817,  715
  408, 1015,  815,  817, 1017,  915
  508, 1215, 1015, 1017, 1217, 1115
  608, 1415, 1215, 1217, 1417, 1315
  708, 1615, 1415, 1417, 1617, 1515
  808, 1815, 1615, 1617, 1817, 1715
  908, 2015, 1815, 1817, 2017, 1915
 1008, 2215, 2015, 2017, 2217, 2115
*ELSET, ELSET=L1, GENERATE
 1,7,1
*ELSET, ELSET=U1, GENERATE
 101,107,1
*ELSET, ELSET=D7, GENERATE
 201,207,1
*ELSET, ELSET=D6, GENERATE
 301,307,1
*ELSET, ELSET=D5, GENERATE
 401,407,1
*ELSET, ELSET=D4, GENERATE
 501,507,1
*ELSET, ELSET=D3, GENERATE
 601,607,1
*ELSET, ELSET=D2, GENERATE
 701,707,1
*ELSET, ELSET=D1, GENERATE
 801,807,1
*ELSET, ELSET=T1, GENERATE
 901,907,1
*ELSET, ELSET=S1, GENERATE
 1001,1007,1
**
** SOIL STRATA L1 (3600 ft - 4800 ft)
** ----------------------------------
*MATERIAL,NAME=L1
*ELASTIC
 360000.,0.29
*DRUCKER PRAGER, SHEAR CRITERION=LINEAR
 38.,1.,0.
*DRUCKER PRAGER HARDENING
 286.,0.
 460.,0.0037
 358.,0.04
*PERMEABILITY,SPECIFIC=3.61e-2
 15.2, 
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA U1 (3000 ft - 3600 ft)
** ----------------------------------
*MATERIAL,NAME=U1
*ELASTIC
 68000.,0.22
*DRUCKER PRAGER
 38.,1.,0.
*DRUCKER PRAGER HARDENING
 286.,0.
 460.,0.0037
 358.,0.04
*PERMEABILITY,SPECIFIC=3.61e-2
 15.2, 
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D7 (2400 ft - 3000 ft)
** ----------------------------------
*MATERIAL,NAME=D7
*ELASTIC
 162600.,.17
*CAP PLASTICITY
 247.,58.5,0.23,0.02,0.,1.
*CAP HARDENING
  500.,.0
  600.,.006
 1100.,.012
 9800.,.03
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-4, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D6 (2000 ft - 2400 ft)
** -----------------------------------
*MATERIAL,NAME=D6
*ELASTIC
 112400.,.17
*CAP PLASTICITY
 247.,50.2,0.23,0.02,0.,1.
*CAP HARDENING
  400.,.0
  600.,.005
 1100.,.02
 9000.,.05
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-4, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D5 (1800 ft - 2000 ft)
** ----------------------------------
*MATERIAL,NAME=D5
*ELASTIC
 71600.,.17
*CAP PLASTICITY
 200.,40.4,0.30,0.02,0.,1.
*CAP HARDENING
  400.,.0
  700.,.02
 2000.,.04
 9000.,.08
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-4, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D4 (1400 ft - 1800 ft)
** ----------------------------------
*MATERIAL,NAME=D4
*ELASTIC
 59600.,.2
*CAP PLASTICITY
 174.,40.1,0.30,0.02,0.,1.
*CAP HARDENING
  200.,.0
  730.,.02
 1000.,.04
 9000.,.1
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-4, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D3 (1200 ft - 1400 ft)
** ----------------------------------
*MATERIAL,NAME=D3
*ELASTIC
 79200.,.19
*CAP PLASTICITY
 200.,42.0,0.34,0.02,0.,1.
*CAP HARDENING
  200.,.0
  500.,.02
 2000.,.04
 9000.,.06
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-4, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA D2 (1100 ft - 1200 ft)
** ----------------------------------
*MATERIAL,NAME=D2
*ELASTIC
 63000.,.17
*CAP PLASTICITY
 200.,39.4,0.33,0.02,0.,1.
*CAP HARDENING
  200.,.0
  600.,.02
 1000.,.04
 8000.,.1
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-5, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
** 
** SOIL STRATA D1 (700 ft - 1100 ft)
** ---------------------------------
*MATERIAL,NAME=D1
*ELASTIC
 47700.,.17
*CAP PLASTICITY
 200.,36.9,0.33,0.02,0.,1.
*CAP HARDENING
 400.,.0
 600.,.02
 800.,.05
 900.,.09
*CAP CREEP, LAW=SINGHM, MECHANISM=CONSOLIDATION, TIME=TOTAL
 2.2e-7, 2.1e-2, 1.0, 1.0, 50.0
 3.5e-5, 2.1e-2, 1.0, 1.0, 212.0
*PERMEABILITY,SPECIFIC=3.07e-2
 16.8,,0
 32.0,,250
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA T1 (150 ft - 700 ft)
** --------------------------------
*MATERIAL,NAME=T1
*ELASTIC
 300000.,0.25
*DRUCKER PRAGER
 36.,1.,0.
*DRUCKER PRAGER HARDENING
 70.,0.
 90.,.058
 70.,.116
*PERMEABILITY,SPECIFIC=3.61e-2
 40, 
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
**
** SOIL STRATA S1 (0 ft - 150 ft)
** ------------------------------
*MATERIAL,NAME=S1
*ELASTIC
 18000.,0.3
*DRUCKER PRAGER
 42.,1.,0.
*DRUCKER PRAGER HARDENING
 11.,0.
 12.,.058
 11.,.116
*PERMEABILITY,SPECIFIC=3.61e-2
 40., 
*EXPANSION
 0.32E-05, 
*EXPANSION, POREFLUID 
 0.32E-06, 
** FAR FIELD MATERIAL DATA
** -----------------------
*MATERIAL,NAME=FAR
*ELASTIC
 47700.,.17
**
*SOLID SECTION, ELSET=L1,  MATERIAL=L1
*SOLID SECTION, ELSET=U1,  MATERIAL=U1
*SOLID SECTION, ELSET=D7,  MATERIAL=D7
*SOLID SECTION, ELSET=D6,  MATERIAL=D6
*SOLID SECTION, ELSET=D5,  MATERIAL=D5
*SOLID SECTION, ELSET=D4,  MATERIAL=D4
*SOLID SECTION, ELSET=D3,  MATERIAL=D3
*SOLID SECTION, ELSET=D2,  MATERIAL=D2
*SOLID SECTION, ELSET=D1,  MATERIAL=D1
*SOLID SECTION, ELSET=T1,  MATERIAL=T1
*SOLID SECTION, ELSET=S1,  MATERIAL=S1
*SOLID SECTION, ELSET=FAR, MATERIAL=FAR
**
*BOUNDARY
 FREESURF, 8,,0.
**
*INITIAL CONDITIONS, TYPE=RATIO
 NALL, 1.5
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
 NALL, 50
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
  L1, -2151.2,     0.,-1615.7, 14400., .85 , .85
  U1, -1615.7, 14400.,-1346.4, 21600., .85 , .85
  D7, -1346.4, 21600.,-1077.1, 28800., .85 , .85
  D6, -1077.1, 28800., -897.6, 33600., .85 , .85
  D5,  -897.6, 33600., -807.8, 36000., .85 , .85
  D4,  -807.8, 36000., -628.3, 40800., .85 , .85
  D3,  -628.3, 40800., -538.6, 43200., .85 , .85
  D2,  -538.6, 43200., -493.7, 44400., .85 , .85
  D1,  -493.7, 44400., -314.2, 49200., .85 , .85
  T1,  -314.2, 49200.,  -67.3, 55800., .85 , .85
  S1,   -67.3, 55800. ,   0.0, 57600., .85 , .85
  FAR,-2151.2,     0.,    0.0, 57600., .85 , .85
**
** 
**  GEOSTATIC STEP TO BALANCE GRAVITY LOADING
**
*STEP, UNSYMM=YES, NLGEOM
*GEOSTATIC
*MONITOR,NODE=2201,DOF=2
*BOUNDARY
 BASE, 2
 AXIS, 1
*DLOAD
 S1,  BZ, -0.0374
 T1,  BZ, -0.0374
 D1,  BZ, -0.0374
 D2,  BZ, -0.0374
 D3,  BZ, -0.0374
 D4,  BZ, -0.0374
 D5,  BZ, -0.0374
 D6,  BZ, -0.0374
 D7,  BZ, -0.0374
 U1,  BZ, -0.0374
 L1,  BZ, -0.0374
*END STEP
**
** CONSOLIDATION ANALYSIS (1 MONTH - BALANCE CREEP)
**
*STEP,NLGEOM,INC=100,UNSYMM=YES
*SOILS,CONSOLIDATION,CETOL=0.0001,UTOL=10
 1.,31.,
*EL FILE, FREQ=0, ELSET=ESOIL
 S,E,PE,POR,VOIDR,CE
*END STEP
**
** THERMAL LOADING - INJECTING STEAM
**
*STEP, NLGEOM, INC=100, UNSYMM=YES
*SOILS, CONSOLIDATION, UTOL=25, NO CREEP
 1.,1.
*TEMPERATURE,OP=MOD
 TARGET, 212.
*EL FILE, FREQ=0, ELSET=ESOIL
 S,E,PE,POR,VOIDR,CE
*END STEP
** 
** INSTALL PUMP
**
*STEP, NLGEOM, INC=100, UNSYMM=YES, AMP=RAMP
*SOILS, CONSOLIDATION, UTOL=25
 1.,1.
*BOUNDARY, OP=MOD
 1001,8,8,-170
 1201,8,8,-170
*EL FILE, FREQ=0, ELSET=ESOIL
 S,E,PE,POR,VOIDR,CE
*END STEP
**
** PUMPING OIL OVER FIVE YEAR PERIOD
**
*STEP,NLGEOM,INC=200,UNSYMM=YES
*SOILS,CONSOLIDATION,CETOL=0.0001,UTOL=20
 1.,1825.,
*CONTROLS, PARAMETER=FIELD, FIELD=DISPLACEMENT
 ,1.
*CONTROLS, PARAMETER=FIELD, FIELD=PORE FLUID PRESSURE
 ,1.
*EL FILE,FREQ=200, ELSET=ESOIL
 S,E,PE,CE,POR,VOIDR
*EL PRINT,FREQ=0, ELSET=ESOIL
 S,
 CE,
 PE,
*NODE PRINT,FREQ=0
 U,RVF,RVT
*END STEP