In this new version, the user can define any number of fixation sets and name them as he/she wants. It is thus possible to define much more complex senarios than before.
The previous commands should be gathered in a single list of fixations:
name: a user-defined string, mainly used to name the group of nodes,
nodes: as for fixations, it can be a point defined by its 3 coordinates between square brackets ([x, y, z]), or a list of points defined in the same manner (e.g. [[x1, y1, z1], [x2, y2, z2], ...]) or a mesh file containing one or more vertices.
values: the force vector applied to each node.
A new model named models/dolicorhynchops/dolicorhynchops_10k_extloads.py has been created to test this new feature.
️ New load distribution models on surface meshes
In addition to the "Ad Hoc Uniform Traction Model" ('U'), the "Tangential-Traction Model" (T) and the "Tangential-Plus-Normal-Traction Model" (T+N), the user can now use 2 new load distributions:
Directional-Traction Model (D)
The forces are distributed on the nodes of a given surface (file key). They are all pointing in the focal_pt direction seen as a vector. They are scaled so that the total resulting force has the magnitude prescribed by the user (with the key force)
Example:
p['muscles'] = [ { 'method': 'D', 'file': f'{path}/Lmuscle.stl', # surface mesh 'force': 100., # total force 'focalpt': [-1, 1, 1], # <= used as a direction } ]
Note: the keyword focal_pt can be misleading here. The vector given as argument is a direction. However we have kept the string focal_pt as a keyword to have the same keywords for all the methods.
:dart:Normal-Traction Model (N)
In this new method, the nodal forces are applied in the perpendicular direction to each facet of the surface mesh. The focalpt keyword is ignored. Forces are scaled to have the prescribed resulting magnitude (the sum of all the forces is force).
Boman Romain authored