JCFpm

Asked by Mahdeyeh on 2019-08-10

Hi guys

I want to model crack and micro joints in rock and I am completely amateur in YADE.
this is my code:

rad = 1e-1 # unit is meter
c=6.54e6 # unit is Pa
Q=55.12 # unit is degree should convert to radian
v=0.20
T=12.21e6 # unit is Pa
E=24.36e9 # unit is Pa
param = {'cohesion':c, 'frictionAngle':Q*2*math.pi/360, 'poisson':v, 'tensileStrength':T, 'young':E}
Mat=O.materials.append(JCFpmMat(density=2.89e3,jointNormalStiffness=2.5e6,jointShearStiffness=1e6,jointCohesion=1e6,jointDilationAngle=0.3,jointFrictionAngle=0.5,**param))

O.bodies.append([
   utils.sphere(center=(0,0,0),radius=rad,material=Mat),
   utils.sphere((0,1.5*rad,0),radius=rad,material=Mat)
])
O.bodies.append(utils.wall(position=(0,0,-1),axis=2,sense=+1,material=Mat))

O.forces.setPermF(0,Vector3(10,20,3))
O.forces.setPermF(1,Vector3(10,20,3))
O.dt = 1e-6

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
   InteractionLoop(
     [Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
     [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys()],
     [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM()]
     ),
   NewtonIntegrator(damping=0.1,gravity=(0,0,-9.81))
]
O.saveTmp()

1- How I can add " JCFpmState " with its arguments to body?? (https://www.yade-dem.org/doc//yade.wrapper.html#yade.wrapper.JCFpmState)

2- How I can add a joint with the following properties to first sphere?
dip=84 degree
dip direction=280 degree
persistence=0.8 m
aperture=1 mm

3- How I can change arguments of JCFpmPhys ?? (https://www.yade-dem.org/doc//yade.wrapper.html#yade.wrapper.JCFpmPhys)

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
2019-08-10
Last reply:
2019-08-13
Luc Scholt├Ęs (luc) said : #1

Hi,

Welcome to YADE.

I don't see the point for changing JCFpmState or the JCFpmPhys arguments...

What do you want to simulate exactly (process wise)?

There are an examples in the "examples/jointedCohesiveFrictionalPM/" folder [1] which should help.

Please have a look at the different files and let us know if anything is unclear.

Luc

[1] https://gitlab.com/yade-dev/trunk/tree/master/examples/jointedCohesiveFrictionalPM

Mahdeyeh (mahdiye.sky) said : #2

Thank you Luc for your response.

My 3 question were just for more learning and not related to the model I said.

Can you help with this problem?

Provide an answer of your own, or ask Mahdeyeh for more information if necessary.

To post a message you must log in.