Structural Changes of Active Skeletal Muscles: Modelling, Validation and Numerical Experiments

Hadi Rahemi,Nilima Nigam,James M. Wakeling
DOI: https://doi.org/10.48550/arXiv.1504.00413
2015-04-02
Abstract:The purpose of this study was to report numerical validation of a 3D finite element model of contracting muscle. The model was based on continuum theory for fibre-reinforced composite materials. Here we simulated contractions for an idealized medial gastrocnemius muscle in man, using the model. Simulations were performed to test the force-length relation of the whole muscle, to evaluate the changes in internal fascicle geometry during contractions, and to assess the importance of material formulations for the aponeurosis and tendon. The simulation results were compared to previously published experimental values. The force-length profile for the whole muscle showed a realistic profile. As the muscle contracted the fascicles curved into S-shaped trajectories and curled around 3D paths, both of which matched previous experimental findings. As the fascicles shortened they increased in their cross-sectional area, but this increase was asymmetric with the smaller increase occurring within the fascicle-plane: the Poisson's ratio in this plane matched that previously shown from ultrasound imaging. The distribution of strains in the aponeurosis and tendon was shown to be a function of their material properties. This study demonstrated that the model could replicate realistic patterns of whole muscle-force, and changes to the internal muscle geometry, and so will be useful for testing mechanisms that affect the structural changes within contracting muscle.
Quantitative Methods,Tissues and Organs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to verify the validity of a three - dimensional finite element model (3D FEM) in simulating contracting muscles, especially for simulating the structural changes of the human medial gastrocnemius muscle. Specifically, the researchers hope to test and verify the following points through this model: 1. **Force - length relationship**: Evaluate the force output of the entire muscle at different lengths. 2. **Internal muscle fiber geometric changes**: Assess the geometric changes of internal muscle fibers (fascicle) during muscle contraction, including their curved paths and curling patterns. 3. **Importance of material properties**: Evaluate the influence of the material properties of the aponeurosis and tendon on the structural changes of the muscle. To achieve these goals, the researchers conducted a variety of numerical experiments and compared the simulation results with published experimental data. In this way, they hope to prove that the model can realistically reproduce the overall force output and internal structural changes during muscle contraction, thereby providing a reliable tool for future research. ### Summary of main problems - **Verify the validity of the model**: By comparing with experimental data, ensure that the model can accurately reflect the behavior of the muscle during contraction. - **Understand the relationship between muscle structure and function**: Explore how the internal structure of the muscle (such as muscle fiber arrangement, bending, etc.) affects its function (such as force output). - **Evaluate the influence of material properties**: Study the influence of different material properties of the aponeurosis and tendon on the muscle contraction behavior. ### Key formulas The material property formulas involved in the paper are as follows: - **Stress - stretch relationship of tendon along muscle fiber**: \[ \sigma_{\text{Tend}}(\lambda) = \begin{cases} 10^4\times1.904\times(\lambda^{68.8}- 1), & \text{if }1\leq\lambda\leq1.07 \\ 10^4\times1.904\times(6758\times(\lambda - 1.07)+ 104.1), & \text{if }\lambda > 1.07 \end{cases} \] - **Strain energy of tendon matrix material**: \[ \Psi_{\text{Tend}} = 10^4\times2.857\times(I_1 - 3) \] - **Stress - stretch relationship of aponeurosis along muscle fiber**: \[ \sigma_{\text{Apo}}(\lambda) = \begin{cases} 10^6\times3.053\times(\lambda^{124.6}- 1), & \text{if }1\leq\lambda\leq1.025 \\ 10^6\times3.053\times(17375\times(\lambda - 1.025)+ 20.7), & \text{if }\lambda > 1.025 \end{cases} \] - **Strain energy of aponeurosis matrix material**: \[ \Psi_{\text{Apo}} = 10^4\times57.84\times e^{579.6\times(I_1 - 3)} \] These formulas are used to describe the mechanical behaviors of tendons and aponeuroses and are the key mathematical expressions in the model.