Summary

Based on the temporally and spatially discretized weak form, we've explored methods to compute the mass matrix, deformation gradient, and elasticity force under the linear finite element setting, all of which align with our implementation in Section Inversion-Free Elasticity.

With linear finite elements, the world space coordinates are approximated as a piecewise linear function of . This approximation, , is a linear function inside each triangle and is -continuous at the edges. By using two parameters, and , to represent points on each triangle, we can identify the linear shape functions that interpolate the displacements at the triangle vertices and derive the deformation gradient . The mass matrix entries and elasticity terms can then be computed via integration with respect to and .