← Back to team overview

yade-users team mailing list archive

[Question #682755]: JCFpm

 

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

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)



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