yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #10793
[Question #261082]: Save/compute average strain 0.5*(∇u+∇u^t)
New question #261082 on Yade:
https://answers.launchpad.net/yade/+question/261082
Dear all,
The TesselationWrapper is able to compute deformation from given chronological states/positions of a packing. Though it writes particle-based strain tensors to in vtk format and returns the average strain (sym_grad_u), I didn't find a way to save the strain in Matrix3 format so that I can compute the deviatoric part.
An alternative way is to triangulate the packing (Regular/Delaunay) following Bagi's method such that the products would be similar to what TesselationWrapper produces. The Delaunay triangulation can be easily imported from scipy [1], but the regular triangulation [2] looks more preferable because it is also used in Yade. Unfortunately I don't know how to import it to Yade.
Is it possible to save the average strain returned from the TesselationWrapper? If not, could anyone tell me how to import the regular triangulation package to Yade without touching the source code?
[1] http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.spatial.Delaunay.html
[2] http://doc.cgal.org/latest/Triangulation_3/
Cheers
Alex
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.