yade-users team mailing list archive

[Anton Gladky <question251712@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Anton Gladky proposed the following answer:
Hi,

2014-07-16 5:11 GMT+02:00 Hongyang Alex Cheng
<question251712@xxxxxxxxxxxxxxxxxxxxx>:
> 1. With VTKRecorder it is possible to export normal stress and shear stress of a particle but what are the physical meanings of normStress and shearStress and how are these stresses computed? I checked the source code and found the following sentences might give some hints, but I failed to fully understand them.
>
> In VTKRecorder.cpp: Line 454 and 455
>
>                         const Vector3r& stress = bodyStates[b->getId()].normStress;
>                         const Vector3r& shear = bodyStates[b->getId()].shearStress;

Implementation is here [1]. Feel free to check it and let us know
whether it is correctly implemented. I have doubts, crossSection
is defined as Pi*R^2, should it not be the volume?

We should probably return those to values from bodyStressTensors().

[1] https://github.com/yade/trunk/blob/master/pkg/dem/Shop_02.cpp#L176


Anton

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.