*Heading
 Surface Ablation Test: 3D Perf-A
** Job name: perf_3D-A Model name: Model-1
**
** PARTS
**
*Part, name=Formation
*End Part
*PREPRINT,MODEL=YES,HISTORY=YES,ECHO=YES
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Formation-1, part=Formation
*Node
*Element, type=C3D8P
*Element, type=M3D4
** Region: (Rock:Picked)
*Elset, elset=_I1, internal
** Section: Rock
*Solid Section, elset=_I1, material=Rock
1.,
** Region: (Casing:Picked)
*Elset, elset=_I2, internal
** Section: Casing
*Solid Section, elset=_I2, material=Rock
1.,
** Region: (Cement:Picked)
*Elset, elset=_I3, internal
** Section: Cement
*Solid Section, elset=_I3, material=Rock
1.,
** Region: (Cement:Picked)
*Elset, elset=_I4, internal
** Section: Cement
*Solid Section, elset=_I4, material=Rock
1.,
** Region: (Casing ID Face:Picked)
*Elset, elset=_I5, internal, generate
** Section: Casing ID Face
*Membrane Section, elset=_I5, material="Low Stifness"
0.001,
** Region: (Perforation Face:Picked)
*Elset, elset=_I6, internal, generate
** Section: Perforation Face
*Membrane Section, elset=_I6, material="Low Stifness"
0.001,
** Region: (Well Bore Face 2:Picked)
*Elset, elset=_I7, internal, generate
** Section: Well Bore Face 2
*Membrane Section, elset=_I7, material="Low Stifness"
0.001,
*End Instance
**
*Elset, Elset="Well Bore Skin", instance=Formation-1, gen
*Nset, nset="Well Bore", instance=Formation-1
*Elset, elset="Well Bore", instance=Formation-1
*Nset, nset=Cement, instance=Formation-1
*Elset, elset=Cement, instance=Formation-1
*Nset, nset=Casing, instance=Formation-1
*Elset, elset=Casing, instance=Formation-1
*Nset, nset=Perf_vol_Rock, instance=Formation-1
*Elset, elset=Perf_vol_Rock, instance=Formation-1
*Nset, nset="Well Bore Face", instance=Formation-1
*Elset, elset="Well Bore Face", instance=Formation-1
*Nset, nset=Face-1, instance=Formation-1
*Elset, elset=Face-1, instance=Formation-1
*Nset, nset=Face-2, instance=Formation-1
*Elset, elset=Face-2, instance=Formation-1
*Nset, nset=Face-radial, instance=Formation-1
*Elset, elset=Face-radial, instance=Formation-1
*Nset, nset=Face-Bottom, instance=Formation-1
*Elset, elset=Face-Bottom, instance=Formation-1
*Nset, nset="Well Bore Vol-Face", instance=Formation-1
*Elset, elset="Well Bore Vol-Face", instance=Formation-1
*Nset, nset=perf_vol_Cement-Casing, instance=Formation-1
*Elset, elset=perf_vol_Cement-Casing, instance=Formation-1
*Nset, nset="Plot Group A", instance=Formation-1
*Elset, elset="Plot Group A", instance=Formation-1
*Nset, nset=Perf_vol_Nodes, instance=Formation-1
*Elset, elset=Perf_vol_Nodes, instance=Formation-1
*Nset, nset=Face_Cement-Casing-Perf, instance=Formation-1
*Elset, elset=Face_Cement-Casing-Perf, instance=Formation-1
*Nset, nset=Face_Casing_w-Perf, instance=Formation-1
*Elset, elset=Face_Casing_w-Perf, instance=Formation-1
*Nset, nset=Face_Casing_ID, instance=Formation-1
*Elset, elset=Face_Casing_ID, instance=Formation-1
*Nset, nset=Face_Rock_Perf, instance=Formation-1
*Elset, elset=Face_Rock_Perf, instance=Formation-1
*Nset, nset=Face-Top, instance=Formation-1
*Elset, elset=Face-Top, instance=Formation-1
*Nset, nset="All Rock", instance=Formation-1, generate
*Elset, elset="All Rock", instance=Formation-1, generate
*Nset, nset="Adaptive Zone 1", instance=Formation-1
*Elset, elset="Adaptive Zone 1", instance=Formation-1
*Nset, nset="Center Line", instance=Formation-1
*Elset, elset="Center Line", instance=Formation-1
*Elset, elset=_Top_face_S2, internal, instance=Formation-1
*Elset, elset=_Top_face_S4, internal, instance=Formation-1
*Elset, elset=_Top_face_S3, internal, instance=Formation-1
*Surface, type=ELEMENT, name=Top_face
_Top_face_S2, S2
_Top_face_S4, S4
_Top_face_S3, S3
*Elset, elset=_Face-2_S2, internal, instance=Formation-1
*Elset, elset=_Face-2_S6, internal, instance=Formation-1
*Elset, elset=_Face-2_S4, internal, instance=Formation-1
*Elset, elset=_Face-2_S3, internal, instance=Formation-1
*Surface, type=ELEMENT, name=Face-2
_Face-2_S2, S2
_Face-2_S6, S6
_Face-2_S4, S4
_Face-2_S3, S3
*Elset, elset=_Face-1_S4, internal, instance=Formation-1
*Elset, elset=_Face-1_S3, internal, instance=Formation-1
*Elset, elset=_Face-1_S6, internal, instance=Formation-1
*Elset, elset=_Face-1_S2, internal, instance=Formation-1
*Elset, elset=_Face-1_S1, internal, instance=Formation-1, generate
*Surface, type=ELEMENT, name=Face-1
_Face-1_S4, S4
_Face-1_S3, S3
_Face-1_S6, S6
_Face-1_S2, S2
_Face-1_S1, S1
*Elset, elset=_Face-Bottom_S4, internal, instance=Formation-1
*Elset, elset=_Face-Bottom_S3, internal, instance=Formation-1
*Elset, elset=_Face-Bottom_S2, internal, instance=Formation-1
*Elset, elset=_Face-Bottom_S1, internal, instance=Formation-1, generate
*Elset, elset=_Face-Bottom_S6, internal, instance=Formation-1, generate
*Surface, type=ELEMENT, name=Face-Bottom
_Face-Bottom_S4, S4
_Face-Bottom_S3, S3
_Face-Bottom_S2, S2
_Face-Bottom_S6, S6
_Face-Bottom_S1, S1
**
**
*NSET, NSET=NODE_boundary
Face-1, Face-2, Face-radial, "Well Bore Face"
**
*Transform, nset=Node_boundary, type=C
**
**
*Nset, nset=outPOre, instance=Formation-1
*Nset, Nset=LAG,   instance=Formation-1
**
*End Assembly
************************************************************************
** MATERIALS
**  Units: inch, lbf, psi
**
**
*Material, name=Rock
*Drucker Prager
*Drucker Prager Hardening
**2490 , 0.2
*Elastic
1.32e6, 0.22
*Permeability, specific=0.67
150     ,0.271597
180     ,0.285892
225     ,0.300939
275     ,0.316778
325     ,0.33345
400     ,0.351
440     ,0.35802
490     ,0.36518
**
*Material, name="Low Stifness"
*Elastic
 0.001, 0.2
**
**
** Porosity 26%
*Initial Conditions, type=Ratio
"All Rock" ,  0.3513514
**
*Initial Conditions, type=Pore
"All Rock" ,  5500.00
*Initial Conditions, type=Stress
"All Rock" ,  -5000., -5000., -7500.0, 0.0, 0.0, 0.0
** ----------------------------------------------------------------
**
** STEP: geostatic
**
*Step, name=geostatic, nlgeom
Establish geostatic equilibrium
*Geostatic
*Boundary
Face-1, 1, 2, 0.00
Face-2, 1, 2, 0.00
Face-Radial, 1, 1, 0.00
Face-Bottom, 3, 3, 0.00
"All Rock", 8, 8, 5500.
"Center Line", 1, 1, 0.000
**
*DSLOAD
TOP_Face , P, 13000.00
**
** OUTPUT REQUESTS
**
***Restart, write, overlay, frequency=1
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*ELEMENT OUTPUT
PE, PEEQ, S
*NODE OUTPUT
U,POR
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=999999
*End Step
***********************************************************
** STEP 2: Drilling Step
**
**
** Elset _I5 is casing ID
**       _I6 is perforation face
**       _I7 is well bore face
**
***********************************************************
*Step, name="Drilling", inc=6, nlgeom=yes
*SOILS, utol=300.0
***
*** Borehole and Perf removal
***
*Model Change, remove
"Well Bore", Perf_vol_Rock, Formation-1._I5
***
**
*Boundary
"Well Bore Face", 1, 1
*************************************************
*** LOADS
****
*** Borehole pressure applied to borehole face
*Dload
Formation-1._I6, P, -5800.0
**
***********************************************************
** OUTPUT REQUESTS
***********************************************************
***Restart, write, frequency=1, overlay
*Output, field
*Element Output, elset="all rock"
PEEQ, PE,  S
*NODE OUTPUT
U,POR
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=999999
*End Step
**
********************************
**  Step 3: Changing Boundary conditions
**
********************************
*Step, name="DD_1000", inc=1, nlgeom=yes
*SOILS, utol=10.0
***********************************************************
*** BOUNDARY CONDITIONS
***********************************************************
*Boundary, op=NEW
Face-1, 2, 2, 0.00
Face-2, 2, 2, 0.00
Face-Radial, 1, 1, 0.00
"Well Bore Face", 1, 1, 0.00
Face-Bottom, 3, 3, 0.00
Face-Radial, 8, 8, 5500.
Face_Rock_PERF, 8, 8, 5500.0
***********************************************************
***********************************************************
** OUTPUT REQUESTS
***********************************************************
***Restart, write, frequency=1, overlay
*Output, field
*Element Output, elset="all rock"
PEEQ, PE, S
*NODE OUTPUT
U, POR
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=1, nset=outPOre
POR,
*End Step
************************************************************
**  Step 4: Go to the start of Material Removal
**
*************************************************************
*Step, name="DD_1500", inc=100, nlgeom=yes
*SOILS, utol=10.0
***********************************************************
*** BOUNDARY CONDITIONS
***********************************************************
*Boundary
Face_Rock_PERF, 8, 8, 4960.0
***********************************************************
*** LOADS
***********************************************************
*** Borehole pressure applied to borehole face
*Dload
Formation-1._I6, P, -4960.0
***********************************************************
** OUTPUT REQUESTS
***********************************************************
***Restart, write, frequency=1, overlay
*Output, field
*Element Output, elset="all rock"
PEEQ, PE, S
*NODE OUTPUT
U,POR
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=1, nset=outPOre
POR,
*El Print, Elset="Adaptive Zone 1"
VOLC
*Adaptive Mesh, Elset="Adaptive Zone 1", Freq=100, Mesh=20
*Adaptive Mesh Constraint, Constraint Type=Lagrangian
LAG
*End Step
***
************************************************************
**  Step 5: Erode material at "Constant" Draw Down
**          Pressure using "ADAPTIVE MESHING"
*************************************************************
*Step, name="Erode", inc=10000, nlgeom=yes
*SOILS, CONSOLIDATION, UTOL=100
**THE TIME PERIOD FOR THIS STEP HAS BEEN REDUCED FROM 4.5 DAYS
**TO 0.6 DAYS TO SHORTEN THE RUNNING TIME. TO OBTAIN THE RESULTS
**IN THE EXAMPLE PROBLE MANUAL, USE THE DATA LINE BELOW TO
**INCREASE THE TIME PERIOD TO 4.5 DAYS
**0.3, 4.5, 1e-4, 0.3
***********************************************************
*** BOUNDARY CONDITIONS
***********************************************************
*Boundary
Face_Rock_PERF, 8, 8, 4960.0
***********************************************************
*** LOADS
***********************************************************
*** Borehole pressure applied to borehole face
*Dload
Formation-1._I6, P,  -4960.0
***********************************************************
** OUTPUT REQUESTS
***********************************************************
***Restart, write, frequency=1, overlay
*Output, field
*Element Output, elset="all rock"
PEEQ, PE, FLVEL, IVOL, S
*NODE OUTPUT
U,POR
*Output, history
*Element Output, Elset="Adaptive Zone 1"
VOLC
*El Print, freq=999999
*Node Print, freq=1, nset=outPOre
POR,
*El Print, Elset="Adaptive Zone 1"
VOLC
*Adaptive Mesh, Elset="Adaptive Zone 1", Freq=1, Mesh=40
*Adaptive Mesh Constraint, Constraint Type=Lagrangian
LAG
*Adaptive Mesh Constraint, Type=Velocity, User
Face_Rock_Perf,1,1,1.0
*End Step