*Heading
** Job name: ale_cpe4p_washout_2 Model name: Well Bore Wash Out
*Preprint, echo=NO, model=NO, history=NO, contact=NO
*Part, name=formation-current
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=formation-current-1, part=formation-current
*Node
*Element, type=CPE4P
*Nset, nset=_PickedSet2, internal, generate
*Elset, elset=_PickedSet2, internal, generate
** Section: Rock
*Solid Section, elset=_PickedSet2, material=Rock
1.,
*End Instance
**
*Nset, nset="Formation Face", instance=formation-current-1
*Elset, elset="Formation Face", instance=formation-current-1, generate
*Nset, nset="All Rock", instance=formation-current-1, generate
*Elset, elset="All Rock", instance=formation-current-1, generate
*Nset, nset="Far Field U_p", instance=formation-current-1
*Elset, elset="Far Field U_p", instance=formation-current-1
*Nset, nset=ale_domain, instance=formation-current-1
*Elset, elset=ale_domain, instance=formation-current-1, generate
*Nset, nset=xsymm, instance=formation-current-1
*Elset, elset=xsymm, instance=formation-current-1
*Nset, nset=ysymm, instance=formation-current-1
*Elset, elset=ysymm, instance=formation-current-1
*Elset, elset=_formationFace_S3, internal, instance=formation-current-1, generate
*Surface, type=ELEMENT, name=formationFace
_formationFace_S3, S3
*Elset, elset=_xPressure_S2, internal, instance=formation-current-1, generate
*Elset, elset=_xPressure_S1, internal, instance=formation-current-1, generate
*Surface, type=ELEMENT, name=xPressure
_xPressure_S2, S2
_xPressure_S1, S1
*Elset, elset=_yPressure_S2, internal, instance=formation-current-1, generate
*Elset, elset=_yPressure_S1, internal, instance=formation-current-1, generate
*Surface, type=ELEMENT, name=yPressure
_yPressure_S2, S2
_yPressure_S1, S1
*Nset, nset="All Rock_PP_", internal, instance=formation-current-1, generate
*Nset, nset="Far Field U_p_PP_", internal, instance=formation-current-1
*Nset, nset="Formation Face_PP_", internal, instance=formation-current-1
*End Assembly
**
** MATERIALS
**
*Material, name=Rock
*Drucker Prager
*Drucker Prager Hardening
*Elastic
1000., 0.22
*Permeability, specific=0.0076
 2.34833e-08, 0.04
 1.17341e-07, 0.07
 3.36324e-07,  0.1
 7.45578e-07, 0.13
 1.42563e-06, 0.16
  2.4774e-06, 0.19
 4.02852e-06, 0.22
 6.24154e-06, 0.25
**
** BOUNDARY CONDITIONS
**
** Name: Dir-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
xsymm, XSYMM
** Name: Dir-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
ysymm, YSYMM
** Name: Initial Pore Pressure Type: Pore pressure
*Boundary
"All Rock_PP_", 8, 8
*map solution
** ----------------------------------------------------------------
**
** STEP: Equilibrate
**
*Step, name=Equilibrate, nlgeom=YES
*Geostatic
**
** BOUNDARY CONDITIONS
**
** Name: Dir-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
xsymm, XSYMM
** Name: Dir-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
ysymm, YSYMM
** Name: Far Field Pore Type: Pore pressure
*Boundary, op=NEW
"Far Field U_p_PP_", 8, 8
** Name: Initial Pore Pressure Type: Pore pressure
*Boundary, op=NEW
** Name: underbalance drilling Type: Pore pressure
*Boundary, op=NEW
"Formation Face_PP_", 8, 8, -3.4475
*Adaptive Mesh Controls, name=Ada-1, geometric enhancement=YES
**
** LOADS
**
** Name: Load-1   Type: Pressure
*Dsload
formationFace, P, 6.895
** Name: xPressure   Type: Pressure
*Dsload
xPressure, P, 24.7
** Name: yPressure   Type: Pressure
*Dsload
yPressure, P, 19.5
**
** CONTROLS
**
*Controls, reset
*Controls, parameters=field, field=displacement
, 1., , , , , ,
*Controls, parameters=field, field=pore fluid pressure
, 1., , 0.01, , , ,
**
** OUTPUT REQUESTS
*Restart, write, frequency=0
**
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, freq=9999
*Node Output
CF, POR, RF, U
*Element Output, directions=YES, ELSET=ALE_DOMAIN
LE, PEEQ, S, SAT, VOIDR
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: Continue underbalanced
**
*Step, name="Continue underbalanced", nlgeom=YES, inc=10000
Continue underbalanced
*Soils, consolidation, end=PERIOD, utol=10.
*Adaptive Mesh, elset=ale_domain, controls=Ada-1, frequency=1, mesh sweeps=5, op=NEW
**
** ADAPTIVE MESH CONSTRAINTS
**
** Name: Ada-Cons-1 Type: Velocity/Angular velocity
*Adaptive Mesh Constraint, user, type=VELOCITY
"Formation Face",1,,2.e-3
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, freq = 9999
*Node Output
CF, POR, RF, U
*Element Output, directions=YES, ELSET=ALE_DOMAIN
FLVEL, LE, PEEQ, S, SAT, VOIDR
**
** HISTORY OUTPUT: H-Output-2
**
*Output, history
*Element Output, elset=ale_domain
VOLC,
**
** HISTORY OUTPUT: H-Output-1
**
*print,adaptive mesh=yes
*End Step