Creep for CpmMat?

Asked by behzad

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()

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

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jan Stránský
Solved:
Last query:
Last reply:
Revision history for this message
Best Jan Stránský (honzik) said :
#1

Hi Behzad,
sorry for late answer. Yes, you are right, the trunk version of Cpm model
does not have creep.
cheers
Jan

2014-06-27 16:47 GMT+02:00 behzad <email address hidden>:

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

Revision history for this message
behzad (behzad-majidi) said :
#2

Hi Jan,

no problem. Thanks for the answer.

Revision history for this message
behzad (behzad-majidi) said :
#3

Thanks Jan Stránský, that solved my question.