Finite element simulations are widely used to study the non-linear mechanical behavior of various biomaterials. Constructing an anisotropic strain energy function within the framework of hyperelasticity is an effective approach to describe the material constitutive behavior. The constructed strain energy function needs to satisfy several mathematical requirements to ensure the framework indifference, the material (near) incompressibility, and the material stability. While the framework indifference, or objectivity, can be naturally satisfied by the invariant-based constitutive formulation, how to enforce the material incompressibility and stability requires detailed discussions and careful treatments. The scope of this paper is to examine in detail the impacts of the mathematical requirements on the constitutive formulation for non-linear anisotropic biomaterials. Particularly, theoretical analyses and numerical simulations are combined to investigate the influences of the material incompressibility as a mathematical constraint and various convexity conditions on the invariant-based constitutive modeling. Through a constructed boundary value problem, analytical solutions are derived via the energy minimization and further used to quantify the influences of the material anisotropic and isotropic components on the material responses. The impact of the volumetric–deviatoric split is demonstrated separately for the strictly incompressible and nearly incompressible materials. In order to ensure the material stability, two commonly used convexity conditions, including the strong ellipticity condition and the polyconvexity condition, are discussed in detail. Several numerical examples are provided to demonstrate their impacts on the material stability under different loading conditions. These discussions are particularly relevant to model biomaterials that exhibit non-linear and anisotropic behaviors under complex loading conditions.
Read full abstract