A numerical model of passive and active behavior of skeletal muscles

doi:10.1016/S0045-7825(97)00162-X

Computer Methods in Applied Mechanics and Engineering

Volume 151, Issues 3–4, 20 January 1998, Pages 419-433

https://doi.org/10.1016/S0045-7825(97)00162-X Get rights and content

Abstract

In this paper a numerical model is developed for the simulation of the 3-D deformations of skeletal muscles. Muscles are distinct from other materials, particularly from other biological tissues, because of their capability of active contraction. The constitutive relation adopted in this work is a generalization of the model proposed by Humphrey and Yin for the passive behavior of cardiac muscle. This generalization is done in such a way as to make Humphrey's model compatible with Zajac's 1-D model of the passive and active behavior of skeletal muscles. The finite element method is used to obtain numerical results for some simple tests.

References (21)

D.T. Davy et al.
A dynamic optimization technique for predicting muscle forces in the swing phase of gate
J. Biomech.
(1987)
M.G. Pandy et al.
An optimal control model for maximum-height human jumping
J. Biomech.
(1990)
J.C. Simo
On a fully three-dimensional finite-strain viscoelastic damage model: Formulation and computational aspects
Comput. Methods Appl. Mech. Engrg.
(1987)
J.C. Simo et al.
Quasi-incompressible finite elasticity in principal stretches: Continuum basis and numerical algorithms
Comput. Methods Appl. Mech. Engrg.
(1991)
J.A. Weiss et al.
Finite element implementation of incompressible transversely isotropic hyperelasticity
Comput. Methods Appl. Mech. Engrg.
(1996)
D.T. Chen et al.
Pump it up: Computer animation of a biomechanically based model of muscle using the finite element method
Computer Graphics
(1992)
R. Crowningsfield et al.
A physiologically based criterion of muscle force prediction in locomotion
J. Biomech.
(1981)
G. Engel
Optimization techniques in the analysis of the shoulder and arm mechanism
Y.C. Fung
Biomechanics: Mechanical Properties of Living Tissues
(1981)
A.W.J. Gielen et al.
A finite element formulation of muscle contraction

There are more references available in the full text version of this article.

Cited by (175)

A finite element-based machine learning framework to predict the mechanical behavior of the pelvic floor muscles during childbirth
2024, Expert Systems with Applications
The medical community has been focusing on gaining a deeper understanding of birth trauma, which affects millions of women worldwide. Maternal lesions can be challenging to diagnose and expensive to examine. To better comprehend the mechanism of injuries occurring in the pelvic floor muscles (PFM), biomechanical simulations can be a valuable tool. However, utilizing the finite element method (FEM) to conduct simulations can be a time-consuming process. To overcome this issue, the present study aims to develop a machine learning (ML) framework to predict stresses on the PFM during childbirth by training ML algorithms on FEM simulation data. To generate the dataset for the ML algorithm’s training, childbirth simulations were performed using different material properties to characterize the PFM. Four ML algorithms were employed, namely Random Forest (RF), Extreme Gradient Boosting (XGBT), Support Vector Regression (SVR), and Artificial Neural Networks (ANN), considering two scenarios: (1) stress prediction for the maximum stretch level of the muscle, and (2) for multiple levels of fetal descent. Results showed that the ANN performed best in the former, with a mean absolute error (MAE) of 0.191 MPa. In the latter, XGBT provided lower errors for 20 and 35 mm of fetal descent, with MAE values of 0.002 and 0.028 MPa, respectively. Nevertheless, the ANN yielded better predictions for 50 and 65 mm, with MAE values of 0.214 and 0.187 MPa, respectively. The present work represents the first attempt to use FEM-based ML algorithms with childbirth simulations to obtain near-real-time predictions in routine clinical procedures.
A biomechanical perspective on perineal injuries during childbirth
2024, Computer Methods and Programs in Biomedicine
Childbirth trauma is a major health concern that affects millions of women worldwide. Severe degrees of perineal trauma, designated as obstetric anal sphincter injuries (OASIS), and levator ani muscle (LAM) injuries are associated with long-term morbidity. While significant research has been conducted on LAM avulsions, less attention has been given to perineal trauma and OASIS, which affect up to 90% and 11% of vaginal deliveries, respectively. Despite being widely discussed, childbirth trauma remains unpredictable. This work aims to enhance the modeling of the maternal musculature during childbirth, with a particular focus on understanding the mechanisms underlying the often overlooked perineal injuries.
A geometrical model of the pelvic floor muscles (PFM) and perineum (including the perineal body, ischiocavernosus, bulbospongiosus, superficial and deep transverse perineal muscles) was created. The muscles were characterized by a transversely isotropic visco-hyperelastic constitutive model. Two simulations of vaginal delivery were conducted with the fetus in the vertex presentation and occipito-anterior position, with and without the perineum.
The simulation that considered the perineum exhibited higher stresses over an extended area of the PFM, which suggests that including additional structures can impact the obtained results. The maximum stretch of the urogenital hiatus was 2.94 and the maximum stress was 23.86 kPa. The perineal body reached a maximum stretch of 1.95, which was more pronounced near the urogenital hiatus, where perineal tears may occur. The external anal sphincter's transverse diameter decreased by 51% and the maximum principal stresses were observed in the area close to the perineal body, where OASIS can occur.
The present study emphasizes the importance of including most structures involved in vaginal delivery in its biomechanical analysis and represents another step further in the understanding of perineal injuries and OASIS. The superior region of the perineal body and its connection to the urogenital hiatus and anal sphincter have been identified as the most critical regions, highly susceptible to injury.
Model of calf muscle tear during a simulated eccentric contraction, comparison between ex-vivo experiments and discrete element model
2023, Journal of the Mechanical Behavior of Biomedical Materials
The tearing of the muscle-tendon complex (MTC) is one of the common sports-related injuries. A better understanding of the mechanisms of rupture and its location could help clinicians improve the way they manage the rehabilitation period of patients. A new numerical approach using the discrete element method (DEM) may be an appropriate approach, as it considers the architecture and the complex behavior of the MTC. The aims of this study were therefore: first, to model and investigate the mechanical elongation response of the MTC until rupture with muscular activation. Secondly, to compare results with experimental data, ex vivo tensile tests until rupture were done on human cadavers {triceps surae muscle + Achilles tendon}. Force/displacement curves and patterns of rupture were analyzed. A numerical model of the MTC was completed in DEM. In both numerical and experimental data, rupture appeared at the myotendinous junction (MTJ). Moreover, force/displacement curves and global rupture strain were in agreement between both studies. The order of magnitude of rupture force was close between numerical (858 N for passive rupture and 996 N–1032 N for rupture with muscular activation) and experimental tests (622 N ± 273 N) as for the displacement of the beginning of rupture (numerical: 28–29 mm, experimental: 31.9 mm ± 3.6 mm). These differences could be explained by choices of DEM model and mechanical properties of MTC's components or their rupture strain values. Here we show that he MTC was broken by fibers' delamination at the distal MTJ and by tendon disinsertion at the proximal MTJ in agreement with experimental data and literature.
Numerical investigation of the biomechanical effects of orthodontic intermaxillary elastics on the temporomandibular joint
2023, Journal of the Mechanical Behavior of Biomedical Materials
Temporomandibular joint disorder (TMD) often coincides with malocclusion, and in some cases TMDs are reported after orthodontic treatment. Intermaxillary elastics (also known as orthodontic elastics, OE) are a common way to apply force during orthodontic treatment, and they might cause mechanical effects on the temporomandibular joint (TMJ), thereby lead to joint remodeling. It is still a controversial topic whether the adapted remodeling of the TMJ or of the alveolar bone is the main cause for the alteration of occlusion after treatment with OEs. It was the aim of this study to analyze whether variations of OEs would develop harmful effects on the healthy TMJ. A TMJ model with a masticatory system based on Hill-type muscle actuators was established. Mouth opening and closure with and without OEs were simulated, and maximum principal stresses in the disc and condylar cartilage as well as the displacement of the mandible were analyzed. We found no considerably difference in the mandibular movement without and with symmetrical OEs during mouth opening and closing. At full mouth opening, stresses in the disc and condylar cartilage of some models with OEs were much smaller than without OEs, but we did not find consistency in the results from the left and right sides of the same model (e.g. the lowest compressive stress on the left side of disc from the model with Class II OEs is much smaller than without OEs, −17.3 MPa compared with −28.2, while on the right side, there was no obvious difference). Hence, we could not conclude that OEs would develop deleterious effects on the healthy TMJ.
A comparative study of 85 hyperelastic constitutive models for both unfilled rubber and highly filled rubber nanocomposite material
2022, Nano Materials Science
Nonlinear finite element analysis is widely used for structural optimization of the design and the reliability analysis of complex elastomeric components. However, high-precision numerical results cannot be achieved without reliable strain energy functions (SEFs) of the rubber or rubber nanocomposite material. Although hyperelastic constitutive models have been studied for nearly 80 years, selecting one that accurately describes rubber's mechanical response is still a challenge. This work reviews 85 isotropic SEFs based on both the phenomenological theory and the micromechanical network theory proposed from the 1940s to 2019. A fitting algorithm which can realize the automatic fitting optimization and determination of the parameters of all SEFs reviewed is developed. The ability of each SEF to reproduce the experimental data of both the unfilled and highly filled rubber nanocomposite is quantitatively assessed based on a new proposed evaluation index. The top 30 SEFs for the unfilled rubber and the top 14 SEFs for the highly filled rubber nanocomposite are presented in the ranking lists. Finally, some suggestions on how to select an appropriate hyperelastic constitutive model are given, and the perspective on the future progress of constitutive models is summarized.
Effect of mesh anchoring technique in uterine prolapse repair surgery: A finite element analysis
2021, Journal of Biomechanics
The female pelvic cavity involves muscles, ligaments, endopelvic fasciae and multiple organs where different pathologies may occur, namely the pelvic organ prolapse (POP). The synthetic implants are used for the reconstructive surgery of POP, but severe complications associated with their use have been reported, mainly related to their mechanical properties (e.g., implant stiffness) and microstructure.
In this study, we mimicked a transvaginal reconstructive surgery to repair the apical ligaments (uterosacral ligaments (USLs) and cardinal ligaments (CLs)), by modeling, their impairment (90% and 50%) and/or total rupture. The implants to reinforce/replace these ligaments were built based on literature specifications and their mechanical properties were obtained through uniaxial tensile tests. The main aim of this study was to simulate the effect of mesh anchoring technique (simple stich and continuous stitch), and compare the displacement magnitude of the pelvic tissues, during Valsalva maneuver.
The absence/presence of the synthetic implant was simulated when total rupture of the CLs and USLs occurs, causing a variation of the vaginal displacement (9% for the CLs and 27% for the USLs). Additionally, the simulations showed that there was a variation of the supero-inferior displacement of the vaginal wall between different anchoring techniques (simple stich and continuous stitch) being approximately of 10% for the simulation USLs and CLs implant.
The computational simulation was able to mimic the biomechanical behavior of the USLs and CLs, in response to different anchoring techniques, which can be help improving the outcomes of the prolapse surgery.

View all citing articles on Scopus

View full text

A numerical model of passive and active behavior of skeletal muscles

Abstract

J. Biomech.

J. Biomech.

Comput. Methods Appl. Mech. Engrg.

Comput. Methods Appl. Mech. Engrg.

Comput. Methods Appl. Mech. Engrg.

Pump it up: Computer animation of a biomechanically based model of muscle using the finite element method

Computer Graphics

A physiologically based criterion of muscle force prediction in locomotion

J. Biomech.

Optimization techniques in the analysis of the shoulder and arm mechanism

Biomechanics: Mechanical Properties of Living Tissues

A finite element formulation of muscle contraction