On 10/24/2017 06:08 PM, Kathrin Sobe wrote:
I was looking for a postprocessing method to calculate the displacements for points on the surface of the mesh that are not necessarily node points. Precisely I need a interpolation function that gives me the displacement of these points. The input would be the already calculated linear elastic problem from my earlier example with the node displacements as result.
In the sfepy documentation I found the probing described here: http://sfepy.org/doc-devel/primer.html#sec-primer
Is this the right way to solve the interpolation? Do you recommend something different?
Yes, the best is to use FieldVariable.evaluate_at(), which can evaluate a FE field variable in physical points coordinates (i.e. not reference element coordinates) - check tests/test_projections.py or tests/test_mesh_interp.py. Unfortunately, there is no example showing this feature, just the tests.
Thank you and regards, Kathrin _______________________________________________ SfePy mailing list firstname.lastname@example.org https://mail.python.org/mm3/mailman3/lists/sfepy.python.org/