A New Continuum-Based Thick Shell Finite Element for Soft Biological Tissues in Dynamics: Part 2 - Anisotropic Hyperelasticity and Incompressibility Aspects
For biomedical engineers, it provides a reliable and efficient tool for simulating large 3D deformations of soft tissues without penalty factors or Lagrange multipliers.
This paper extends a thick continuum-based shell finite element to model anisotropic hyperelastic and incompressible soft biological tissues in dynamics, achieving excellent agreement with analytical and experimental data for rubber-like materials and tissues.
In a companion article (Part 1), we presented the development of a thick continuum-based (CB) shell finite element (FE) based on Mindlin-Reissner theory. We verified the accuracy, efficiency and locking insensitivity of the element in modeling large 3D deformations, using linear elastic material properties. In the present article, we developed and implemented the kinetics description, within the updated Lagrangian (UL) formulation, of anisotropic incompressible hyperelastic constitutive relations that enable the CB shell FE to accurately model very large 3D strains and deformations. Specifically, we developed the measures of deformation in the lamina coordinate system, presented three techniques to model nonlinear hyperelastic strains, and enabled the direct enforcement of incompressibility and of the zero normal stress condition without using a penalty factor or a Lagrange multiplier. Moving towards the application of the present work to the biomedical realm, we performed multiple experiments concerning mechanical behavior of rubber-like materials and soft biological tissues in different geometries and loading conditions. Excellent agreements between the present FE results and the analytical and/or experimental data proved the CB shell FE combined with the present constitutive techniques to be a highly reliable and efficient tool for modeling, analyzing, and predicting mechanical behavior of soft biological tissues.