search

GCMLab / GCMLab-FEM

General purpose finite element program
0 stars 0 forks source link

Create a function to calculate nodal forces from a traction along an arbitrary edge #45

Open b3gee opened 3 years ago

b3gee commented 3 years ago

Right now we are manually calculating and prescribing nodal forces in the config files. This is easy enough to implement for first-order elements, but becomes more of a pain for higher order elements.

A function which accepts the mesh, traction, and the equation of a line, and returns the vector of nodal forces and corresponding nodes/DoFs would be a nice quality of life enhancement, and make debugging easier.

b3gee commented 3 years ago

Some code to this effect exists in PatchTestC.