yade-users team mailing list archive

[Liu Changdong <question691703@xxxxxxxxxxxxxxxxxxxxx>]

New question #691703 on Yade:
https://answers.launchpad.net/yade/+question/691703

Hello Everyone：
I am a new novice in Yade. Now i want to model hydraulic fracture with DFNFlowEngine and JCFpmMat. The model 1 is a cube, and I will inject water in the center of the cube. But I don't see any cracks in the model,  just the movement of the particles. Then 
 i try the model 2, the model 2 a cube with a joint in the center. And when i inject water in the center, yade reports the following error like: Vh==NULL!! id=0 Point=3.07081 4.72667 8.26223 rad=0.533742
1 : Vh==NULL!! id=1 Point=2.93626 7.27644 5.58063 rad=0.533623
2 : Vh==NULL!! id=2 Point=8.7494 2.35821 3.8193 rad=0.533504
 My question is based on the following:
1) How can cracks caused by water pressure be seen in the cube
2)How to correctly inject water into joints
The script of model 1 I wrote is as follows:

==================
from yade import pack,plot,utils
import math

intR=0.1
DENS=2500
YOUNG=1800
FRICT=7
ALPHA=0.1
TENS=100
COH=1000
iterper=100

O.materials.append(JCFpmMat(type=1,density=DENS,young=YOUNG,poisson=ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH,label='mat1'))

mn,mx=Vector3(0,0,0),Vector3(1,1,.1)
walls=aabbWalls([mn,mx],thickness=0,material='mat1')
wallIds=O.bodies.append(walls)

sp=pack.randomDensePack(pack.inAlignedBox(mn,mx),spheresInCell=2000,radius=0.01,memoizeDb='/tmp/hftest.sqlite',returnSpherePack=True)
sp.toSimulation(color=(0,0,250),material='mat1')

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb()]),
    InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR)],
     [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
     [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=True,label='interactionLaw', recordCracks=True, Key='identifier',neverErase=True),Law2_ScGeom_FrictPhys_CundallStrack()]
    ),
    DFNFlowEngine(dead=1,label="flow"),
    GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5, defaultDt=utils.PWaveTimeStep()),
    NewtonIntegrator(damping=0.1),
]

flow.dead=0
flow.viscosity=10
flow.defTolerance=0.3
flow.meshUpdateInterval=1
flow.updateTriangulation=True
flow.breakControlledRemesh=True
flow.bndCondIsPressure=[0,0,0,0,0,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
flow.imposePressure(Vector3(0.5,0.5,0.05),30000000)
flow.useSolver=3
O.dt=0.1e-5
O.dynDt=False

=============

Can anyone help me with this? Or any tips or tricks.

Thank you

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.