Dear developers of SfePy,
I am relatively new to SfePy. So this question may be very trivial. So how do I compute the linear traction force on surface given a set of nodal pressure values. I am aware that we have a term called "dw_surface_ltr" but the example shows the 'material' are computed on quad nodes. Thanks!
I think I found the key, I can call cmesh.get_incident to get the conn of a specific region of facets and I assume that the quad nodes are constructed based on the order of this conn, right?
Hi,
On 11/29/2017 06:13 AM, benechiao@gmail.com wrote:
I think I found the key, I can call cmesh.get_incident to get the conn of a specific region of facets and I assume that the quad nodes are constructed based on the order of this conn, right?
The surface connectivity is directly available in the SurfaceData of the field
sd = virtual.field.surface_data[self.region.name]
Note that the connectivity from CMesh is not the same as the field connectivity
It should be also possible to construct a surface field for the region, define a corresponding variable, set its DOFs to the nodal pressure values, and then set the material parameter of dw_surface_ltr by a function, where you just evaluate the variable in the required coordinates. Let me know if you need more advice (send the example script).
Otherwise for your purpose, maybe you can use dw_point_load term?
r.
Robert,
Thanks for your advice. I will test several things out based on the information you gave.