Implementation of controlling strategy in a biomechanical lower limb model with active muscles for coupling multibody dynamics and finite element analysis

Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion

The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.

Development of a three-dimensional muscle-driven lower limb model developed using an improved CFD-FE method

Analysis of the musculoskeletal movements (gait analysis) is needed in many scenarios. The in vivo method has some difficulties. For example, recruiting human subjects for the gait analysis is challenging due to many issues. In addition, when plenty of subjects are required, the follow-up experiments take a long period and the dropout of subjects always occurs. An efficient and reliable in silico simulation platform for gait analysis has been desired for a long time. Therefore, a technique using three-dimensional (3D) muscle modeling to drive the 3D musculoskeletal model was developed and the application of the technique in the simulation of lower limb movements was demonstrated. A finite element model of the lower limb with anatomically high fidelity was developed from the MRI data, where the main muscles, the bones, the subcutaneous tissues, and the skin were reconstructed. To simulate the active behavior of 3D muscles, an active, fiber-reinforced hyperelastic muscle model was developed using the user-defined material (VUMAT) model. Two typical movements, that is, hip abduction and knee lifting, were simulated by activating the responsible muscles. The results show that it is reasonable to use the improved CFD-FE method proposed in the present study to simulate the active contraction of the muscle, and it is feasible to simulate the movements by activating the relevant muscles. The results from the present technique closely match the physiological scenario and thus the technique developed has a great potential to be used in the in silico human simulation platform for many purposes.

Biomechanical effects of typical lower limb movements of Chen-style Tai Chi on knee joint

The load and stress distribution on cartilage and meniscus of the knee joint in typical lower limb movements of Chen-style Tai Chi (TC) and deep squat (DS) were analyzed using finite element (FE) analysis. The loadings for this analysis consisted of muscle forces and ground reaction force (GRF), which were calculated through the inverse dynamic approach based on kinematics and force plate measurements obtained from motion capture experiments. Thirteen experienced practitioners performed four typical TC movements, namely, single whip (SW), brush knee and twist step (BKTS), stretch down (SD), and part the wild horse's mane (PWHM), which exhibit lower posture and greater lower limb force compared to other TC styles. The results indicated that TC required greater lower limb muscle strength than DS, resulting in greater knee joint forces. The stress on the medial cartilage in SW and BKTS fell within a range conductive to maintaining the balance between anabolism and catabolism of cartilage matrix. This was due to the fact that SW and BKTS reduce the medial to total tibiofemoral contact force ratios through knee abduction, which may effectively alleviate mild medial knee osteoarthritis (KOA). However, the greater medial contact force ratios observed in SD and PWHM resulted in great contact stresses that may aggravate the pain of patients with KOA. To mitigate these effects, practitioners should consider elevating their postures appropriately to reduce knee flexion angles, especially during the single-leg support phase. This adjustment can decrease the required muscle strength, load and stress on the knee joint.

Development of a finite element full spine model with active muscles for quantitatively analyzing sarcopenia effects on lumbar load

BACKGROUND AND OBJECTIVE

The musculoskeletal imbalance caused by disease is one of the most critical factors leading to spinal injuries, like sarcopenia. However, the effects of musculoskeletal imbalances on the spine are difficult to quantitatively investigate. Thus, a complete finite element spinal model was established to analyze the effects of musculoskeletal imbalance, especially concerning sarcopenia.

METHODS

A finite element spinal model with active muscles surrounding the vertebrae was established and validated from anatomic verification to the whole spine model in dynamic loading at multiple levels. It was then coupled with the previously developed neuromuscular model to quantitatively analyze the effects of erector spinae (ES) and multifidus (MF) sarcopenia on spinal tissues. The severity of the sarcopenia was classified into three levels by changing the physiological cross-sectional area (PCSA) of ES and MF, which were mild (60% PCSA of ES and MF), moderate (48% PCSA of ES and MF), and severe (36% PCSA of ES and MF).

RESULTS

The stress and strain levels of most lumbar tissues in the sarcopenia models were more significant than those of the normal model during spinal extension movement. The sarcopenia caused load concentration in several specific regions. The stress level of the L4-L5 intervertebral disc and L1 vertebra significantly increased with the severity of sarcopenia and showed relatively larger values than other segments. From the normal model to a severe sarcopenia model, the stress value of the L4-L5 intervertebral disc and L1 vertebra increased by 128% and 113%, respectively. The strain level of L5-S1 also inclined significantly with the severity of sarcopenia, and the relatively larger capsule strain values occurred at lower back segments from L3 to S1.

CONCLUSIONS

In summary, the validated spinal coupling model can be used for spinal injury risk analysis caused by musculoskeletal imbalance. The results suggested that sarcopenia can primarily lead to high injury risk of the L4-L5 intervertebral disc, L1 vertebrae, and L3-S1 joint capsule regarding significant stress or strain variance.

Effect of hip flexion angle on lower limb injuries of occupants in autonomous vehicle crashes

This study aims to determine the influence of the hip flexion angle on the injury trends of lower limbs. An impact model was established using a hybrid human body model and an accurate vehicle model. Simulations were performed in two boundary environments of 25 and 40% overlap impacts under different hip flexion angles. The analysis of the dynamic responses indicated that the hip flexion angle significantly affected the injury trends. The maximum femur index of different overlaps was all found at the minimum hip angle, except for the left femur at 25% overlap rate. Meanwhile, the maximum acetabular stress was all found at the minimum hip angle (approximately 0.09-0.20 GPa). This study provides mechanistic insights into the lower limb injuries associated with complex human postures.

Modeling of active skeletal muscles: a 3D continuum approach incorporating multiple muscle interactions

Skeletal muscles have a highly organized hierarchical structure, whose main function is to generate forces for movement and stability. To understand the complex heterogeneous behaviors of muscles, computational modeling has advanced as a non-invasive approach to evaluate relevant mechanical quantities. Aiming to improve musculoskeletal predictions, this paper presents a framework for modeling 3D deformable muscles that includes continuum constitutive representation, parametric determination, model validation, fiber distribution estimation, and integration of multiple muscles into a system level for joint motion simulation. The passive and active muscle properties were modeled based on the strain energy approach with Hill-type hyperelastic constitutive laws. A parametric study was conducted to validate the model using experimental datasets of passive and active rabbit leg muscles. The active muscle model with calibrated material parameters was then implemented to simulate knee bending during a squat with multiple quadriceps muscles. A computational fluid dynamics (CFD) fiber simulation approach was utilized to estimate the fiber arrangements for each muscle, and a cohesive contact approach was applied to simulate the interactions among muscles. The single muscle simulation results showed that both passive and active muscle elongation responses matched the range of the testing data. The dynamic simulation of knee flexion and extension showed the predictive capability of the model for estimating the active quadriceps responses, which indicates that the presented modeling pipeline is effective and stable for simulating multiple muscle configurations. This work provided an effective framework of a 3D continuum muscle model for complex muscle behavior simulation, which will facilitate additional computational and experimental studies of skeletal muscle mechanics. This study will offer valuable insight into the future development of multiscale neuromuscular models and applications of these models to a wide variety of relevant areas such as biomechanics and clinical research.

Formulation and exploration of novel, intramuscular pressure based, muscle activation strategies in a spine model

Optimization models are often devised to assess spinal stability via estimating individual muscle forces. However, neglecting muscles' fluidic behavior remains an approximation due to the role of muscle pressure in force transmission. The purpose of this study was to leverage a validated Finite Element (FE) model of the spine, inclusive of Intra-Muscular Pressure (IMP), to explore muscle activation strategies towards maintaining equilibrium spinal stability. Three conventional strategies governing minimizing muscle effort, minimizing IVD compressive forces, and maintaining stability at all costs were first investigated to explore model's validity. Thereafter, two novel IMP-based strategies were devised and explored, specifically minimizing and maximizing IMP. The model was previously shown valid in light of in vivo and in silico observations with an average discrepancy of 6%. This being the case, the conventional strategies dictated efficacy in muscular activations whilst maintaining an equilibrium stable position, as quantified in the present paper, with a difference of 9.8% from documented data. In addition, the explored novel IMP-based strategies suggested the presence of a threshold individual muscles IMP, approximately 272 mmHg for the longissimus muscle for example, beyond which muscles potentially start to share radial loads with surrounding tissues, whilst limiting the contraction of the underlying muscles. In conclusion, this study theoretically supports the possibility of activation strategies based on muscular pressure, which the developed, verified, and validated FE spine model was leveraged to investigate. The explored novel IMP-based strategies may have significance in informing clinical applications such as motion analysis and functional electrical stimulation of muscles.

Quantitative cervical spine injury responses in whiplash loading with a numerical method of natural neural reflex consideration

BACKGROUND AND OBJECTIVE

Neural reflex is hypothesized as a regulating step in spine stabilizing system. However, neural reflex control is still in its infancy to consider in the previous finite element analysis of head-neck system for various applications. The purpose of this study is to investigate the influences of neural reflex control on neck biomechanical responses, then provide a new way to achieve an accurate biomechanical analysis for head-neck system with a finite element model.

METHODS

A new FE head-neck model with detailed active muscles and spinal cord modeling was established and globally validated at multi-levels. Then, it was coupled with our previously developed neuromuscular head-neck model to analyze the effects of vestibular and proprioceptive reflexes on biomechanical responses of head-neck system in a typical spinal injury loading condition (whiplash). The obtained effects were further analyzed by comparing a review of epidemiologic data on cervical spine injury situations.

RESULT

The results showed that the active model (AM) with neural reflex control obviously presented both rational head-neck kinematics and tissue injury risk referring to the previous experimental and epidemiologic studies, when compared with the passive model (PM) without it. Tissue load concentration locations as well as stress/strain levels were both changed due to the muscle activation forces caused by neural reflex control during the whole loading process. For the bony structures, the AM showed a peak stress level accounting for only about 25% of the PM. For the discs, the stress concentrated location was transferred from C2-C6 in the PM to C4-C6 in the AM. For the spinal cord, the strain concentrated locations were transferred from C1 segment to around C4 segment when the effects of neural reflex control were implemented, while the gray matter and white matter peak strains were reduced to 1/3 and 1/2 of the PM, respectively. All these were well correlated with epidemiological studies on clinical cervical spine injuries.

CONCLUSION

In summary, the present work demonstrated necessity of considering neural reflex in FE analysis of a head-neck system as well as our model biofidelity. Overall results also verified the previous hypothesis and further quantitatively indicated that the muscle activation caused by neural reflex is providing a protection for the neck in impact loading by decreasing the strain level and changing the possible injury to lower spinal cord level to reduce injury severity.

Analysis of Foot-Ankle-Leg Injuries in Various Under-Foot Impact Loading Environments with a Human Active Lower Limb (HALL) Model

Lower limb injuries caused by under-foot impacts often appear in sport landing, automobile collision, and anti-vehicular landmine blasts. The purpose of the present study was to evaluate a foot-ankle-leg model of the Human Active Lower Limb (HALL) model, and used it to investigate lower leg injury responses in different under-foot loading environments to provide a theoretical basis for the design of physical dummies adapted to multiple loading conditions. The model was first validated in allowable rotation loading conditions, like dorsiflexion, inversion/eversion, and external rotation. Then, its sensitivity to loading rates and initial postures was further verified through experimental data concerning both biomechanical stiffness and injury locations. Finally, the model was used to investigate the biomechanical responses of the foot-ankle-leg region in different under-foot loading conditions covering the loading rate from sport landing to blast impact. The results showed that from -15° plantarflexion to 30° dorsiflexion, the neutral posture always showed the largest tolerance, and more than 1.5 times tolerance gap was achieved between neutral posture and dorsiflexion 30°. Under-foot impacts from 2 m/s to 14 m/s, the peak tibia force increased at least 1.9 times in all postures. Thus, we consider that it is necessary to include initial posture and loading rate factors in the definition of the foot-ankle-leg injury tolerance for under-foot impact loading.

A Novel Neuromuscular Head-Neck Model and Its Application on Impact Analysis

OBJECTIVE

Neck muscle activation plays an important role in maintaining posture and preventing trauma injuries of the head-neck system, levels of which are primarily controlled by the neural system. Thus, the present study aims to establish and validate a neuromuscular head-neck model as well as to investigate the effects of realistic neural reflex control on head-neck behaviors during impact loading.

METHODS

The neuromuscular head-neck model was first established based on a musculoskeletal model by including neural reflex control of the vestibular system and proprioceptors. Then, a series of human posture control experiments was implemented and used to validate the model concerning both joint kinematics of the cervical spine and neck muscle activations. Finally, frontal impact experiments of varying loading severities were simulated with the newly established model and compared with an original model to investigate the influences of the implanted neural reflex controllers on head-neck kinematic responses.

RESULTS

The simulation results using the present neuromuscular model showed good correlations with in-vivo experimental data while the original model even cannot reach a correct balance status. Furthermore, the vestibular reflex is noted to dominate the muscle activation in less severe impact loadings while both vestibular and proprioceptive controllers have a lot of effect in higher impact loading severity cases.

CONCLUSIONS

In summary, a novel neuromuscular head-model was established and its application demonstrated the significance of the neural reflex control in predicting in vivo head-neck responses and preventing related injury risk due to impact loading.

In vitro compressive properties of skeletal muscles and inverse finite element analysis: Comparison of human versus animals

Virtual finite element human body models have been widely used in biomedical engineering, traffic safety injury analysis, etc. Soft tissue modeling like skeletal muscle accounts for a large portion of a human body model establishment, and its modeling method is not enough explored. The present study aims to investigate the compressive properties of skeletal muscles due to different species, loading rates and fiber orientations, in order to obtain available parameters of specific material laws as references for building or improving the human body model concerning both modeling accuracy and computational cost. A series of compressive experiments of skeletal muscles were implemented for human gastrocnemius muscle, bovine and porcine hind leg muscle. To avoid long-time preservation effects, all experimental tests were carried out in 24 h after that the samples were harvested. Considering computational cost and generally used in the previous human body models, one-order hyperelastic Ogden model and three-term simplified viscoelastic quasi-linear viscoelastic (QLV) were selected for numerical analysis. Inverse finite element analysis was employed to obtain corresponding material parameters. With good fitting records, the simulation results presented available material parameters for human body model establishment, and also indicated significant differences of muscle compressive properties due to species, loading rates and fiber orientations. When considering one-order Ogden law, it is worthy of noting that the inversed material parameters of the porcine muscles are similar to those of the human gastrocnemius regardless of fiber orientations. In conclusion, the obtained material parameters in the present study can be references for global human body and body segment modeling.