yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #07336
[Question #226370]: get python-access to energy information of Hertz-Mindlin contact model
New question #226370 on Yade:
https://answers.launchpad.net/yade/+question/226370
Hi,
I want to investigate energy dissipation by damping and friction with Hertz-Mindlin model. I modified examples/simple-scene/simple-scene-energy-tracking.py a little bit:
...
O.trackEnergy=True
O.engines=[
...
[Ip2_FrictMat_FrictMat_MindlinPhys(betan=0.9,betas=0.9,label='ContactModel')],
[Law2_ScGeom_MindlinPhys_Mindlin(neverErase=False,label='Mindlin')]#,calcEnergy=True)]
),
...
]
...
Mindlin.calcEnergy=True
def myAddPlotData():
normal_Work = Mindlin.normElastEnergy()
shear_Work = Mindlin.shearEnergy()
norm_damp_dissip = Mindlin.normDampDissip()
shear_damp_dissip = Mindlin.shearDampDissip()
fric_dissip = Mindlin.frictionDissipation()
...
Now I am facing a TypeError: 'float' object is not callable
Looking at the source code [1] gives me the hint, that these informations are not defined as usual ...
((OpenMPAccumulator<Real>,frictionDissipation,,Attr::noSave,"Energy dissipation due to sliding"))
((OpenMPAccumulator<Real>,shearEnergy,,Attr::noSave,"Shear elastic potential energy"))
((OpenMPAccumulator<Real>,normDampDissip,,Attr::noSave,"Energy dissipated by normal damping"))
((OpenMPAccumulator<Real>,shearDampDissip,,Attr::noSave,"Energy dissipated by tangential damping"))
[1] https://github.com/yade/trunk/blob/master/pkg/dem/HertzMindlin.hpp#L146
Regards,
christian
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.