Rodin/Solid/Integrators/InternalForce.h file

Internal force vector integrator for hyperelastic formulations.

Assembles the nonlinear residual vector (internal force) contribution:

\[ R_{\text{int}}(\mathbf{v}) = \int_\Omega \mathbf{P}(\mathbf{u}) : \nabla \mathbf{v} \, dX \]

where $ \mathbf{P} $ is the first Piola-Kirchhoff stress and $ \mathbf{v} $ is the test function.

The integrator is generic: it obtains the finite element basis from the FE space (not hardcoded to P1), supports arbitrary quadrature rules, and builds a ConstitutivePoint (composed over Geometry::Point) at each quadrature point for constitutive evaluation.

Namespaces

namespace Rodin
The Rodin library for finite element methods and shape optimization.
namespace Rodin::Solid
Hyperelastic solid mechanics module for large-deformation problems.

Classes

template<class LawDerived, class FES>
class Rodin::Solid::InternalForce
Linear form integrator for the internal force vector in hyperelastic problems.