yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #20303
[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.