yade-users team mailing list archive

[behzad <question250860@xxxxxxxxxxxxxxxxxxxxx>]

New question #250860 on Yade:
https://answers.launchpad.net/yade/+question/250860

Hi,

The material model of CpmMat cannot have creep, right? 
This is not supposed to be viscoelastic?


==========================

from yade import utils, plot
from yade import pack, qt


bMat=O.materials.append(CpmMat(young=12e5,frictionAngle=.8,poisson=.2,density=3800,sigmaT=0.0,relDuctility=2,neverDamage=False, epsCrackOnset=5e-2,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=False,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(omegaThreshold=0.5),Law2_ScGeom_FrictPhys_CundallStrack()],
),
NewtonIntegrator(damping=0.7,gravity=[0,0,-10]),
PyRunner(command='AutoData()',iterPeriod=1000)
]



O.forces.addF(1,(0,0,2.21e5),permanent=True)
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.