yade-users team mailing list archive

[behzad <question248562@xxxxxxxxxxxxxxxxxxxxx>]

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.