New question #248562 on Yade:
https://answers.launchpad.net/yade/+question/248562
Hi,
I have the following script for two contacting spheres with CpmMat.
I have assigned a strain rate in tension to see at which force the bond between
the spheres break.
The code is plotting the interaction force vs. the displacement of the upper
sphere.
The problem is that the graph and the force at which the bond breaks is not a
function of sigmaT.
In this case the obtained graph always has a peak at around 7.5e3 even if we
change sigmaT. Changing epsCrackOnset is changing the strain at which the
maximum force happens.
What's the point I'm missing here?
=====================
O.reset()
from yade import utils, plot
from yade import pack, qt
bMat=O.materials.append(CpmMat(young=24e6,frictionAngle=.8,poisson=.2,density=3800,sigmaT=1e5,relDuctility=30,epsCrackOnset=1e-4,isoPrestress=0))
id_Mat2=O.materials.append(FrictMat(young=1e10,poisson=0.3,density=1000,frictionAngle=1))
Rigid=O.materials[id_Mat2]
s1=utils.sphere([0.0,0.0,1.0],1.0,fixed=True,material=bMat)
s2=utils.sphere([0.0,0.0,3.0],1.0,fixed=True,material=bMat)
O.bodies.append(s1)
O.bodies.append(s2)
box_1 = O.bodies.append(box((0,0,0),(1,1,.001),fixed=True,material=Rigid))
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_CpmMat_CpmMat_CpmPhys(), Ip2_FrictMat_CpmMat_FrictPhys()],
[Law2_ScGeom_CpmPhys_Cpm(),Law2_ScGeom_FrictPhys_CundallStrack()],
),
NewtonIntegrator(damping=0.7,gravity=[0,0,0]),
PyRunner(command='AutoData()',iterPeriod=1000)
]
O.bodies[1].state.vel=(0,0,0.000001)
qt.View()
def AutoData():
i = O.interactions[0,1]
plot.addData(
d=O.bodies[1].state.pos[2]-3.0,
F=i.phys.normalForce.norm())
plot.plots={'d':('F')}
plot.plot()
====================
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
