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.
_______________________________________________
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
