yade-users team mailing list archive

[drasszkusz <question695966@xxxxxxxxxxxxxxxxxxxxx>]

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

drasszkusz posted a new comment:
Thank you Jan, I updated Yade and it works as dict now :) Here is an MWE for those who will have the same problem in the future:
###
def vtkexporter():
   vtk.exportFacets(ids=facetid, what=dict(vel='b.state.vel',pos='b.state.pos',dspl='b.state.displ()',coordnumber='len(b.intrs())',forceLen='O.forces.f(b.id).norm()',forceVec='O.forces.f(b.id)',stress='O.forces.f(b.id).norm()/b.shape.area'))
###

About lenForce, stress, stressLen, etc , I'm still figuring them out.
According to [5], stressLen is the sum of normal and shear stresses and
the norm of the summed value. So yes, Len is the magnitude. I don't know
if stress='O.forces.f(b.id).norm()/b.shape.area' cause it is being
retrieved from another file's algorithm with the following commands
in[5]:

##
const Vector3r& stress = bodyStates[b->getId()].normStress;
const Vector3r& shear  = bodyStates[b->getId()].shearStress;
##

I'm still searching that how are they calculated.

Cheers,
D

[5] https://gitlab.com/yade-
dev/trunk/-/blob/master/pkg/dem/VTKRecorder.cpp

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.