yade-users team mailing list archive

[Jan Stránský <question248562@xxxxxxxxxxxxxxxxxxxxx>]

Question #248562 on Yade changed:
https://answers.launchpad.net/yade/+question/248562

    Status: Open => Answered

Jan Stránský proposed the following answer:
Bi Behzad,
have a look at [1], sigmaT has no effect on normal direction.
Cheers
Jan

[1] http://beta.arcig.cz/~eudoxos/smilauer2010-phd-thesis.pdf
Dne 12.5.2014 17:02 "behzad" <question248562@xxxxxxxxxxxxxxxxxxxxx>
napsal(a):

> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.