Bending

Definition

Name	Mesh>Bending
Type	Bending
Defined in	Edges

Use normals that belong to the two triangles that share an edge and calculate bending energy as [SN88]:

\[U = \frac{1}{2} \kappa \vert \hat{\mathbf{n}}_k- \hat{\mathbf{n}}_l \vert^2\]

where \(( \kappa )\) is the bending rigidity and \(( \hat{\mathbf{n}}_i, \hat{\mathbf{n}}_j )\) are the normal unit vectors of the triangles meeting at the edge.

Python calling

evolver.add_force("Mesh>Bending", {"kappa": {"1": str(kappa_value)}})