*Heading
  mini-frac simulation based on SPE 63239 & SPE 56135
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** units kN and KPa and m
*parameter
 ultI   = 0.3200e+3
 ultII  = 1.0e+10
** ultII  = 0.3200e+3
 Emod   = 8.5E7
** GIc    = 0.33e2
** GIIc   = 0.80e2
 GIc    = 28.
 GIIc   = 28.
** GIIc   = 0.00001
 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    = 0.00e+0
** c_l    = 0.00e+0
 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
*ngen,nset=inter
*nset,nset=nall
 nall,inter
*nset,nset=mouth
*nset,nset=mcod
*ncopy,shift,old set=nall,new set=nall_new,change number=10000
*ngen,nset=inter_new
*nset,nset=inter
 inter,inter_new
*nset,nset=mouth_new
*nset,nset=mouth
 mouth,mouth_new
*nset,nset=mcod_new
*nset,nset=mcod
 mcod,mcod_new
*nset,nset=nall
 nall,nall_new
**
*Element, type=C3D8RP,elset=solid_eall
*ELEMENT,TYPE=COH3D8P,ELSET=inter
*COHESIVE SECTION,ELSET=inter,MATERIAL=MAT1,RESPONSE=TRACTION SEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO
<thick>
*SECTION CONTROLS,NAME=VISCO,VISCOSITY=0.01
*MATERIAL,NAME=MAT1
*ELASTIC,TYPE=TRACTION
<Emod>,<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=<eta>
<GIc>, <GIIc>, <GIIc>
*GAP FLOW
<amu>
*FLUID LEAKOFF
<c_u>, <c_l>
**
*Solid Section, elset=solid_eall, material=Material-1
**
** 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.
** specific gravity of water in kN/m^3
**
*Nset, nset=_PickedSet90
*Nset, nset=_PickedSet91
**
*INITIAL CONDITIONS, TYPE=RATIO
 nall,0.333333
** Initial Porosity is 25%
*INITIAL CONDITIONS, TYPE=PORE
 nall, 793.00e+0
*INITIAL CONDITIONS,TYPE=INITIAL GAP
 2022
**
*INITIAL CONDITIONS, TYPE=STRESS
solid_eall, -12.547e+3, -8.547e+3, -15.00e+3
**
*elset,elset=one
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet90, 1, 1
_PickedSet90, 2, 2
nall, 3, 3
** Name: BC-2 Type: Displacement/Rotation
*Boundary
_PickedSet91, 1, 1
**
*Amplitude, name=volumerate
0.0,0.0, 63.0,-100.0
**0.0,0.0, 63.0,-18.0
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES, inc=1000, unsymm=YES
Injection ~0.1085m^3 over 140 second
*Soils, consolidation, end=PERIOD, utol=28000.00
** time in seconds
**
*Controls, PARAM=FIELD, FIELD=PORE
, , 1.0e-3
**
*cflow,amplitude=volumerate
2119,,<qbound>
12119,,<qbound>
**
** Name: BC-2 Type: Pore pressure
*Boundary
_PickedSet90, 8, 8, 793.00e+0
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=20, overlay
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
**Output, field
*Node Output
U
**POR,U
*element output
LE
**S,LE
GFVR
**SDEG
DMICRT
PFOPEN
**LEAKVRT
**LEAKVRB
ALEAKVRT
ALEAKVRB
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*node output,nset=mouth
por
*node output,nset=mcod
u2
*element output,elset=one
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
*el print
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
*node print,nset=inter
por
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom=YES, inc=1000, unsymm=YES
Falloff after injection 900 seconds
*Soils, consolidation, end=PERIOD, utol=28000.00
** time in seconds
**
*cflow,op=new
**
** Name: BC-2 Type: Pore pressure
*Boundary
_PickedSet90, 8, 8, 793.00e+0
**
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet90, 1, 1
_PickedSet90, 2, 2
nall,3
** Name: BC-2 Type: Displacement/Rotation
*Boundary
_PickedSet91, 1, 1
** OUTPUT REQUESTS
**
*Restart, write, frequency=20, overlay
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
POR,U
*element output
S,NFORC,  LE
GFVR
**SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*node output,nset=mouth
por
*node output,nset=mcod
u2
*element output,elset=one
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
*el print
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
*End Step