*Heading
  Hydralic-driven fracture simulation using XFEM
*Preprint, echo=NO, model=yes, history=yes, contact=NO
**
** units kN and KPa and m
*parameter
 ultI   = 0.3200e+3
 ultII  = 1.0e+10
 Emod   = 8.5E7
 GIc    = 28.
 GIIc   = 28.
 eta    = 2.284
 thick  = 0.01325e-1
 width  = 1.0
** fluid viscosity:  approx 1 centerpoise
 amu    = 1.e-6
** flow rate applied at crack mouth:
 qbound = 5.0e-6
** crack face (filter cake) permeability
 c_u    = 5.879e-10
 c_l    = 5.879e-10
**
perm_matrix = 2.418e-8
**  Perm 2.5 mDarcy ; spec grav 9.8 kN/m^3; mu 1 cp
**
*Node, nset=nall
*Element, type=CPE4P,elset=solid_eall
**
*elset,elset=bc-e1
*Nset, nset=_PickedSet13
*Elset, elset=_PickedSet13
*Nset, nset=_PickedSet14
*Elset, elset=_PickedSet14
**
*enrichment,elset=solid_eall,interaction=lower,type=propagation crack, name=abc
*Solid Section, elset=solid_eall, material=Material-1, CONTROLS=CTR
1.,
**SECTION CONTROLS,NAME=CTR,INITIAL GAP OPENING=0.002
**
** MATERIALS
**
*Material, name=Material-1
*Elastic
 1.294e+06, 0.25
** Modulus of 1.29 GPa or 1.29e6 kPa   Compressibility ~ 5e-6 1/psi
*Permeability, specific=9.8
<perm_matrix>,0.
*GAP FLOW
<amu>
*FLUID LEAKOFF
<c_u>, <c_l>
*DAMAGE INITIATION,CRITERION=MAXPS,tol=0.2
<ultI>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=<eta>
<GIc>, <GIIc>
*DAMAGE STABILIZATION
**
*surface interaction,name=lower
1.0
** specific gravity of water in kN/m^3
**
*INITIAL CONDITIONS, TYPE=RATIO
 nall,0.333333
** Initial Porosity is 25%
*INITIAL CONDITIONS, TYPE=PORE
 nall, 793.00e+0
**
*INITIAL CONDITIONS, TYPE=STRESS
solid_eall, -12.547e+3, -8.547e+3, -15.00e+3
**
*initial condition,type=enrichment
1235,1,abc, 1.0, 0.0
1235,2,abc, 1.0,-1.0
1235,3,abc,-1.0,-1.0
1235,4,abc,-1.0, 0.0
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet13, 1, 1
_PickedSet13, 2, 2
** Name: BC-2 Type: Displacement/Rotation
*Boundary
_PickedSet14, 1, 1
**
*Boundary,phantom=node
214,1,1
215,1,1
**
*Amplitude, name=volumerate
0.0,0.0, 63.0,-200.0
**
**
*restart,write,freq=5
** STEP: Step-1
**
*Step, name=Step-1,unsymm=YES
*geostatic
**
** Name: BC-2 Type: Pore pressure
*Boundary
_PickedSet13, 8, 8, 793.00e+0
**
** OUTPUT REQUESTS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*node output
U,POR,PHILSM
*element output
S
*End Step
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom,inc=1000, unsymm=YES
Injection ~0.02m^3 over 140 second
*Soils, consolidation, end=PERIOD, utol=100000.00
** time in seconds
**2.5, 140., 0.001, 5.0,
**
*cflow,amplitude=volumerate,phantom=edge
214,215,<qbound>
**
** Name: BC-2 Type: Pore pressure
*Boundary
_PickedSet13, 8, 8, 793.00e+0
**Controls, PARAM=FIELD, FIELD=displacement
**,1.0
**Controls, PARAM=FIELD, FIELD=pore
**10.0, 0.5
**
** OUTPUT REQUESTS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field,freq=5
*node output
U,
**POR,
PHILSM
*element output
S
**PFOPENXFEM
**LEAKVRTXFEM
**LEAKVRBXFEM
**ALEAKVRTXFEM
**ALEAKVRBXFEM
**PORPRES
**GFVRXFEM
**
*Output, history,freq=20
*element output,elset=bc-e1
PFOPENXFEM
LEAKVRTXFEM
LEAKVRBXFEM
ALEAKVRTXFEM
ALEAKVRBXFEM
PORPRES
GFVRXFEM
PFOPENXFEMCOMP
PORPRESCOMP
CRDCUTXFEM
*End Step