*heading
*preprint, model=yes, echo=yes, parsubstitution=yes, parvalues=yes
*parameter
 numElem   = 10
 length    = 1.
 width     = 0.1
 height    = 0.1

 x0        = 0.
 y0        = 0.
 z0        = 0.
 x1        = width
 y1        = 0.
 z1        = 0.

 xL0       = 0.
 yL0       = length
 zL0       = 0.
 xL1       = width
 yL1       = length
 zL1       = 0.

 numNodeL  = numElem+1
 numNodeT  = 6*numNodeL
 nf1       = 1
 nf2       = nf1 + numNodeL
 nf3       = nf2 + numNodeL
 nf4       = nf3 + numNodeL
 nf5       = nf4 + numNodeL
 nf6       = nf5 + numNodeL
 nb1       = numNodeL
 nb2       = nb1 + numNodeL
 nb3       = nb2 + numNodeL
 nb4       = nb3 + numNodeL
 nb5       = nb4 + numNodeL
 nb6       = nb5 + numNodeL

 ne1       = 1
 ne2       = nf2
 ne3       = ne2+1
 ne4       = ne1+1
 ne5       = nf3
 ne6       = nf4
 ne7       = ne6+1
 ne8       = ne5+1
 ne9       = nf5
 ne10      = nf6
 ne11      = ne10+1
 ne12      = ne9+1
*preprint, model=yes
*node
 <nf1>,  <x0>, <y0>, <z0>
 <nb1>, <xL0>,<yL0>,<zL0>
 <nf2>,  <x1>, <y1>, <z1>
 <nb2>, <xL1>,<yL1>,<zL1>
 <nf3>,  <x0>, <y0>, <z0>
 <nb3>, <xL0>,<yL0>,<zL0>
 <nf4>,  <x1>, <y1>, <z1>
 <nb4>, <xL1>,<yL1>,<zL1>
 <nf5>,  <x0>, <y0>, <z0>
 <nb5>, <xL0>,<yL0>,<zL0>
 <nf6>,  <x1>, <y1>, <z1>
 <nb6>, <xL1>,<yL1>,<zL1>
*ngen, nset=nEdge1
 <nf1>,<nb1>
*ngen, nset=nEdge2
 <nf2>,<nb2>
*ngen, nset=nEdge3
 <nf3>,<nb3>
*ngen, nset=nEdge4
 <nf4>,<nb4>
*ngen, nset=nEdge5
 <nf5>,<nb5>
*ngen, nset=nEdge6
 <nf6>,<nb6>
*nset, nset=nAll
 nEdge1,nEdge2,nEdge3,nEdge4,nEdge5,nEdge6
*nset, nset=nBot
 nEdge1,nEdge2
*nset, nset=nTop
 nEdge3,nEdge4
*nset, nset=nMid
 nEdge5, nEdge6
*nset, nset=corFront
 <nf1>,<nf2>,<nf3>,<nf4>
*nset, nset=corBack
 <nb1>,<nb2>,<nb3>,<nb4>
*nset, nset=midFront
 <nf5>,<nf6>
*nset, nset=midBack
 <nb5>,<nb6>
*element, type=cod3d8p
 1, <ne1>,<ne2>,<ne3>,<ne4>,<ne5>,<ne6>,<ne7>,<ne8>,<ne9>,<ne10>,<ne11>,<ne12>
*elgen, elset=block
 1,<numElem>,1,1,1
*cohesive section, elset=block, material=core, response=traction separation, thickness=specified
 1e-1,1.0
*material,name=core
*density
 1
*density,pore fluid
 1
*elastic,type=traction
 50,50
*damage initiation,criterion=quads
 0.32, 0.32
*damage evolution,type=energy,mixed mode behavior=bk,power=2.284
 0.012, 0.012
*permeability, specific=10
 2.418e-8
***gap flow,type=HERSCHEL
**1e-4,1.0,6
*gap flow,type=BINGHAM
1e-4,6.0
***gap flow,type=NEWTONIAN
** 1.e-4
***fluid leakoff,user
*fluid leakoff
0,0
*initial conditions,type=pore pressure
 nAll, 0.0
*initial conditions,type=ratio
 nAll, 5.
*initial conditions,type=initial gap
 block
*amplitude, name=ampNT11
 0, 1, 3, 1
*amplitude, name=ampCFLOW
 0.0, 0.0, 0.07,0.0,0.1, 1.0
*amplitude, name=ampDISP
 0.0, 0, 0.06, 1
*step,name=step-1, nlgeom=yes, inc=10000, unsymm=yes,amplitude=step
*soils, consolidation, end=period, utol=28000.
 0.01,0.1,1e-10,1e-2
*boundary
 nAll,     1,2,0.
 nTop,     8,8,0.
 nBot,     8,8,0.
 midBack,  8,8,0.
*boundary, amp=ampDISP
 nBot, 3,3,0.
 nTop, 3,3,0.1
*cflow,amp=ampCFLOW
 midFront,, -3.75e-3
***
***
*output, field
*node output
 por,u,v,cff,rvf,rvt,fldvel
*element output
 gfvr
 por
 fldvel
 fldpgrad
 fluvr
 pfopen
 s,nforc, le
 sdeg
 dmicrt
 sat
 leakvrt
 leakvrb
 aleakvrt
 aleakvrb
*node print, nset=nEdge6
 por
*el print, elset=block
 1,3
 gfvr
*el print, elset=block
 1,3
 leakvrt
*end step
*step,name=step-2, nlgeom=yes, inc=10000, unsymm=yes,amplitude=step
*soils, consolidation, end=period, utol=28000.
 0.1,1.,1e-10
*cflow
 midFront,, -3.75e-3
*end step