K. Garikipati,S. Göktepe,C. Miehe
Abstract:Continuum strain energy functions are developed for soft biological tissues that possess long fibrillar components. The treatment is based on the model of an elastica, which is our fine scale model, and is homogenized in a simple fashion to obtain a continuum strain energy function. Notably, we avoid solving the full fourth-order, nonlinear, partial differential equation for the elastica by resorting to other assumptions, kinematic and energetic, on the response of the individual, elastica-like fibrils.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a new strain - energy density function for soft biological tissues with slender fibrous components. Specifically, the author aims to describe the mechanical behavior of these tissues through an elastomer - based model approach without relying on the entropic - elastic model. The following is a detailed analysis of this problem:
### 1. Research Background
Currently, the strain - energy density functions used to describe soft biological tissues (such as tendons, ligaments, and muscles) are mainly derived from rubber elasticity and polymer elasticity theories, or are specific mathematical forms designed to reproduce the locking behavior observed in experiments. However, for many soft biological tissues (such as tendons), the entropic - elastic model may not be applicable. For example, the collagen molecules in tendons have a triple - helical structure with a diameter of about 1.5 nm and a fully extended length of about 300 nm. These molecules form fibers with a diameter of about 300 nm and a length of up to several hundred micrometers, and are further arranged into fiber bundles that can run through the entire length of the tendon (several centimeters). This hierarchical structure has extensive cross - linking, which limits the degrees of freedom of molecular thermal fluctuations.
### 2. Research Objectives
In order to better describe the mechanical behavior of these tissues, the author proposes a new method based on the elastomer model. Specifically, they hope to improve the existing strain - energy density functions through the following points:
- **Avoid using complex partial differential equations**: Instead of directly solving the fourth - order nonlinear partial differential equations that describe the deformation of elastomers, simplify the problem by introducing other assumptions (kinematic and energy assumptions).
- **Based on internal energy rather than entropy effects**: Considering the influence of cross - linking on molecular configurations, adopt a model based on internal energy instead of a model relying on entropy effects.
- **Combine the micro - and macro - scales**: Homogenize the elastomer model at the micro - scale to obtain the strain - energy density function at the macro - scale, thereby describing the overall mechanical response of soft tissues.
### 3. Method Overview
The author proposes a general and powerful procedure to develop the strain - energy density function of soft tissues based on the elastomer model. The specific steps include:
- **Define the basic problem**: Consider the parametric representation of the elastomer curve Γ in the reference configuration and its position vector in the deformed configuration.
- **Decompose the strain energy**: Assume that the strain energy can be decomposed into two parts, bending and stretching, that is, \( W(\kappa, \lambda)=K(\kappa)+U(\lambda) \), where \( K(\kappa) \) and \( U(\lambda) \) represent the bending and stretching energies respectively.
- **Introduce simplifying assumptions**: In order to avoid solving complex partial differential equations, the author introduces some kinematic and energy assumptions, such as the inextensibility assumption and the planar incompressibility assumption.
- **Solve the force - displacement relationship**: Through the above assumptions, derive the force - displacement responses of elastomers under different conditions and compare them with experimental data.
### 4. Conclusions
Through this method, the author hopes to develop a strain - energy density function that is more in line with the actual mechanical behavior of soft biological tissues. This model can not only explain the characteristic soft - tissue responses observed in experiments, but also provide a theoretical basis for subsequent biomechanical research.
### Formula Summary
The key formulas involved in the paper are as follows:
- Bending energy: \[ K(\kappa)=L\int_{0}^{L}\frac{1}{2}B(\kappa - \kappa_0)^2dS \]
- Stretching energy: \[ U(\lambda)=L\int_{0}^{L}\frac{1}{2}EA(\lambda^2 - 1)^2dS \]
- Free - energy functional: \[ G[u]=L\int_{0}^{L}\left(\frac{1}{2}B\left(\left\|\frac{d^{2}x}{ds^{2}}\right\|-\kappa_0(S)\right)^2+\frac{1}{2}EA\left(\left\|\frac{dx}{dS}\right\|^2 - 1\right)^2 - q\cdot u\right)dS \]
These formulas are used to describe the deformation behavior of elastomers and their energy changes, and are the basis for constructing the strain - energy density function.