yade-users team mailing list archive

[Hongyang Alex Cheng <question251712@xxxxxxxxxxxxxxxxxxxxx>]

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

Hello,

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;

2. The function bodyStressTensors() calculates per-particle micro-stress which is a 3 by 3 matrix. I think it would be better to calculate mean principal stress p and deviatoric stress q from it, which will need a nonlinear equations solver. I tried to import scipy.optimize [1] but it led to core dumped. Does eigen library has some modules to solve nonlinear equations? Are there any software that is importable to Yade will work for this?

I am using git version 2014-06-29.git-de4c01a.

[1] http://docs.scipy.org/doc/scipy/reference/optimize.html#module-scipy.optimize

Thank you,

Alex

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