New question #694718 on Yade:
https://answers.launchpad.net/yade/+question/694718
Dear Yae Users,
I am trying to understand and use an implementation of a combined
Mohr-Coulomb–Weibull (MCW) failure criterion [1].
I have created an example based on the "splitter.py" example prepared by
François Kneib. My example consists of a cube compressed between two walls. I
created a polyhedraMat, which is splittable and has defined all the properties
of Mohr-Coulomb–Weibull failure criterion. I set all the strengths to 1,
probability P=0.63, and high Weibull modulus m=1000 (just in case). I expect to
get strength close to 1 Pa. Instead of that, the strength is around 1.6 e5. It
is independent of any parameter dedicated for MCW criterion. On the other hand,
it depends on "strength" material property (it is not included in the example,
but the default value is 100).
Yade version is 2018.02b
The major issue is:
- Do I have to activate the MCW criterion somehow (or deactivate the previous
one [2])?
And there are two minor issues that I do not understand:
- Do I have to set "Wei_S0". From what I understood from [1] it is some general
value that could be used if the tension and compression strengths were the
same. But if I declare both strengths (sigmaCD, sigmaCZ), I expect that Wei_S0
is calculated based on them. Am I right?
- Since the default "strength" is 100, why do I get 1.6e5 in my "experiment"?
The scale effect should even decrease strength. In the experiment, the contact
force between the wall and the cube is divided by crosssection area (1m2) to
calculate stress.
Cheers,
Karol
[1] Gladkyy, A., & Kuna, M. (2017). DEM simulation of polyhedral particle
cracking using a combined Mohr-Coulomb–Weibull failure criterion. Granular
Matter, 19(3), 41.
[2] Eliáš, J. (2014). Simulation of railway ballast using crushable polyhedral
particles. Powder Technology, 264, 458-465.
---------------------------------------
from yade import polyhedra_utils
from yade import plot
# FUNCTION DENINITIONS
def myAddData():
force = O.forces.f(1)[1]# cube wall interaction
sig = force/(d**2)
plot.addData(t = O.time, force = force, sig = sig)
d = 1 # cube length
#cube material
gravel1 = PolyhedraMat()
gravel1.IsSplitable = True
gravel1.Wei_P = 0.63
gravel1.Wei_m = 1000
gravel1.Wei_S0 = 1
gravel1.Wei_V0 = d**3
gravel1.sigmaCD = 1
gravel1.sigmaCZ = 1
#wall material
steel = PolyhedraMat()
steel.young = 1e10
steel.frictionAngle = 0
# create cube
p1 = polyhedra_utils.polyhedra(gravel1, v=
[(0,0,0),(0,d,0),(d,d,0),(d,0,0),(0,0,d),(0,d,d),(d,d,d),(d,0,d)])
w1 = utils.wall(+d*1.001, axis=1, sense=-1, material=steel)
w2 = utils.wall(-0.001*d, axis=1, sense=+1, material=steel)
v = 5e-4
w1.state.vel = (0,-v,0)
O.bodies.append((p1,w1,w2))
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb()]),
InteractionLoop(
[Ig2_Wall_Polyhedra_PolyhedraGeom(),
Ig2_Polyhedra_Polyhedra_PolyhedraGeom()],
[Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()],
[Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()]
),
NewtonIntegrator(),
PolyhedraSplitter(iterPeriod=100),
PyRunner(command='myAddData()',iterPeriod=2000),
]
O.dt= 0.1*polyhedra_utils.PWaveTimeStep()
plot.plots={'t':'sig'}
plot.plot(subPlots =False)
try:
from yade import qt
qt.Controller()
v = qt.View()
v.ortho = True
except:
pass
O.run(50000),
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
