Development and validation of a discretised multi-body spine model in LifeMOD for biodynamic behaviour simulation

Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges

Numerical models of the musculoskeletal system as investigative tools are an integral part of biomechanical and clinical research. While finite element modeling is primarily suitable for the examination of deformation states and internal stresses in flexible bodies, multibody modeling is based on the assumption of rigid bodies, that are connected via joints and flexible elements. This simplification allows the consideration of biomechanical systems from a holistic perspective and thus takes into account multiple influencing factors of mechanical loads. Being the source of major health issues worldwide, the human spine is subject to a variety of studies using these models to investigate and understand healthy and pathological biomechanics of the upper body. In this review, we summarize the current state-of-the-art literature on multibody models of the thoracolumbar spine and identify limitations and challenges related to current modeling approaches.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Coordination Between Trunk Muscles, Thoracolumbar Fascia, and Intra-Abdominal Pressure Toward Static Spine Stability

STUDY DESIGN

Numerical in-silico human spine stability finite element analysis.

OBJECTIVE

The purpose of this study was to investigate the contribution of major torso tissues toward static spine stability, mainly the thoracolumbar fascia (TLF), abdominal wall with its intra-abdominal pressure (IAP), and spinal muscles inclusive of their intramuscular pressure.

SUMMARY OF BACKGROUND DATA

Given the numerous redundancies involved in the spine, current methodologies for assessing static spinal stability are limited to specific tissues and could lead to inconclusive results. A three-dimensional finite element model of the spine, with structured analysis of major torso tissues, allows for objective investigation of static spine stability.

METHODS

A novel previously fully validated spine model was employed. Major torso tissues, mainly the muscles, TLF, and IAP were individually, and in combinations, activated under a 350N external spine perturbation. The stability contribution exerted by these tissues, or their ability to restore the spine to the unperturbed position, was assessed in different case-scenarios.

RESULTS

Individual activations recorded significantly different stability contributions, with the highest being the TLF at 75%. Combined or synergistic activations showed an increase of up to 93% stability contribution when all tissues were simultaneously activated with a corresponding decrease in the tensile load exerted by the tissues themselves.

CONCLUSION

This investigation demonstrated torso tissues exhibiting different roles toward static spine stability. The TLF appeared able to dissipate and absorb excessive loads, the muscles acted as antagonistic to external perturbations, and the IAP played a role limiting movement. Furthermore, the different combinations explored suggested an optimized engagement and coordination between different tissues to achieve a specific task, while minimizing individual work.Level of Evidence: N/A.

Musculoskeletal modeling and humanoid control of robots based on human gait data

The emergence of exoskeleton rehabilitation training has brought good news to patients with limb dysfunction. Rehabilitation robots are used to assist patients with limb rehabilitation training and play an essential role in promoting the patient's sports function with limb disease restoring to daily life. In order to improve the rehabilitation treatment, various studies based on human dynamics and motion mechanisms are still being conducted to create more effective rehabilitation training. In this paper, considering the human biological musculoskeletal dynamics model, a humanoid control of robots based on human gait data collected from normal human gait movements with OpenSim is investigated. First, the establishment of the musculoskeletal model in OpenSim, inverse kinematics, and inverse dynamics are introduced. Second, accurate human-like motion analysis on the three-dimensional motion data obtained in these processes is discussed. Finally, a classic PD control method combined with the characteristics of the human motion mechanism is proposed. The method takes the angle values calculated by the inverse kinematics of the musculoskeletal model as a benchmark, then uses MATLAB to verify the simulation of the lower extremity exoskeleton robot. The simulation results show that the flexibility and followability of the method improves the safety and effectiveness of the lower limb rehabilitation exoskeleton robot for rehabilitation training. The value of this paper is also to provide theoretical and data support for the anthropomorphic control of the rehabilitation exoskeleton robot in the future.

Development and validation of a timely and representative finite element human spine model for biomechanical simulations

Numerous spine Finite Element (FE) models have been developed to assess spinal tolerances, spinal loadings and low back pain-related issues. However, justified simplifications, in terms of tissue decomposition and inclusion, for such a complex system may overlook crucial information. Thus, the purpose of this research was to develop and validate a comprehensive and representative spine FE model inclusive of an accurate representation of all major torso elements. A comprehensive model comprised of 273 tissues was developed via a novel FE meshing method to enhance computational feasibility. A comprehensive set of indirect validation tests were carried out to validate every aspect of the model. Under an increasing angular displacement of 24°-41°, the lumbar spine recorded an increasing moment from 5.5 to 9.3 Nm with an increase in IVD pressures from 0.41 to 0.66 MPa. Under forward flexion, vertical vertebral displacements simulated a 6% and 13% maximum discrepancy for intra-abdominal and intramuscular pressure results, all closely resembling previously documented in silico measured values. The developed state-of-the-art model includes most physiological tissues known to contribute to spinal loadings. Given the simulation's accuracy, confirmed by its validation tests, the developed model may serve as a reliable spinal assessment tool.

Finite element analysis of knee and ankle joint during gait based on motion analysis

Contact pressures in the articular cartilage during gait affect injuries and the degenerative arthritis of knee and ankle joints. However, only contact forces at the knee and ankle joints during gait can be estimated by using a rigid body dynamic model. The contact pressure distribution can be obtained quantitatively for a static posture by using finite element (FE) analysis in most cases. The purpose of this study is to develop a new method to obtain the contact pressure distribution at the knee and ankle joints during gait by integrating FE analysis with rigid body dynamic analysis. In this method, a reference FE model of the lower extremity is constructed first and is then transformed to each stance phase of the gait obtained from dynamic analysis by using homogeneous transformation. The muscle forces and ground reaction force (GRF) during gait obtained from the dynamic analysis were used as loading conditions for FE analysis. Finally, the contact pressure distribution at the tibia plateau cartilage and talus cartilage were estimated at the 1st peak, mid-stance, and the 2nd peak at the same time. The present method can provide the contact pressure distribution at the knee and ankle joints over the entire gait.

Application of MR-derived cross-sectional guideline of cervical spine muscles to validate neck surface electromyography placement

The importance of surface-EMG placement for development and interpretation of EMG-assisted biomechanical models is well established. Since MR has become a reliable noninvasive cervical spine musculoskeletal diagnostic tool, this investigation attempted to illustrate the anatomical relationships of individual cervical spine muscles with their paired surface-EMG electrodes. The secondary purpose of this investigation was to provide an MR cross-sectional pictorial and descriptive guideline of the cervical spine musculature. MR scans were performed on a healthy adult male subject from skull to manubrium of the sternum. Prior to scanning, MR safe markers were placed over neck muscles following surface EMG placement recommendations. Twenty-three neck muscles were traced manually in each of 267 scan slices. 3-D models of the neck musculoskeletal structure were constructed to aid with understanding the complex anatomy of the region as well as to identify correct EMG electrode locations and to identify muscles' curved lines-of-action. 3D models of the MR-safe markers were constructed relative to the target muscles. Based on the findings of this study, muscle palpation and bony landmarks can be used to effectively identify appropriate surface EMG electrode locations to record upper trapezius, middle trapezius, semispinalis capitis, splenius capitis, levator scapulae, scalenus, sternocleidomastoid and hyoid muscles activities.

Prediction of Cervical Spinal Joint Loading and Secondary Motion Using a Musculoskeletal Multibody Dynamics Model Via Force-Dependent Kinematics Approach

STUDY DESIGN

A cervical spine biomechanical investigation using multibody dynamics.

OBJECTIVE

To develop a comprehensive cervical spine multibody dynamics model incorporated with the force-dependent kinematics (FDK) approach, and to study the influence of soft tissue deformation on the joint loading prediction.

SUMMARY OF BACKGROUND DATA

Musculoskeletal multibody dynamics models have been widely used to analyze joint loading. Current cervical spine musculoskeletal models, however, neglect the joint internal motion caused by soft tissue deformation. A novel FDK approach is introduced, which can predict joint internal motion and spinal joint loading simultaneously.

METHODS

A comprehensive cervical spine musculoskeletal model with the posterior facet joints and essential ligaments was developed. To quantify the influence of soft tissue structures on joint loading prediction, four different models with different features were created. These newly developed models were validated, under flexion-extension movement. The predicted intervertebral disc loads (from C3-C4 to C5-C6) were compared with the published cadaveric experimental results. Moreover, the predicted facet joint forces, ligament forces, and anterior-posterior translations of instantaneous centers of rotation were also studied.

RESULTS

The obtained intervertebral disc loads were varied among different models. Model 3 provided the closest prediction of joint loading to the experimental results. Moreover, the facet joint and ligament forces were in similar range of magnitude as literature findings. The predicted instantaneous centers of rotation translational changes were in accordance with the in vivo kinematics observation.

CONCLUSION

In the present study, a validated cervical spine musculoskeletal model was developed, using multibody dynamics and FDK approach. It can simulate the function of musculature and consider joint internal motion, and thus provides more reliable joint loading prediction. This newly developed cervical model can be used as an efficient tool to study the biomechanical behaviors of human cervical spine, and to understand the fundamental pathologies of spinal pains.

LEVEL OF EVIDENCE

N /A.

RIGID MUSCULOSKELETAL MODELS OF THE HUMAN BODY SYSTEMS: A REVIEW

Advancing the knowledge of the biomechanics of the human body is essential to improve the clinical decision-makings of musculoskeletal disorders in the framework of in silico medicine. An impressive number of research projects focused on the development of rigid-body musculoskeletal models have been conducted over the world thanks to the new research directives. However, the application of these models in clinical practices remains a challenging issue. The objective of this review paper was to present the most current rigid-body musculoskeletal models of the human body systems and to analyze their trends and weaknesses for clinical applications. Then, recommendations were proposed for future researches toward fully clinical decision support. A systematic review process was performed. Well-selected studies related to the most current rigid-body 3D musculoskeletal models for each body system component (jaw, cervical spine, upper limbs, lumbar spine, and lower limbs) were summarized and explored. Trends in rigid musculoskeletal modeling are highlighted as personalization, new imaging techniques for specific joint kinematics, and computational efficiency. Weaknesses are highlighted as modeling assumptions, use of generic model, lack of modeling consensus, model validation, and parameter and model uncertainties. Future directions related to joint and muscle modeling, neuro-musculoskeletal modeling, model validation, data and model uncertainty quantification are recommended.

Incorporating Six Degree-of-Freedom Intervertebral Joint Stiffness in a Lumbar Spine Musculoskeletal Model-Method and Performance in Flexed Postures

Intervertebral translations and rotations are likely dependent on intervertebral stiffness properties. The objective of this study was to incorporate realistic intervertebral stiffnesses in a musculoskeletal model of the lumbar spine using a novel force-dependent kinematics approach, and examine the effects on vertebral compressive loading and intervertebral motions. Predicted vertebral loading and intervertebral motions were compared to previously reported in vivo measurements. Intervertebral joint reaction forces and motions were strongly affected by flexion stiffness, as well as force-motion coupling of the intervertebral stiffness. Better understanding of intervertebral stiffness and force-motion coupling could improve musculoskeletal modeling, implant design, and surgical planning.

A Cervico-Thoraco-Lumbar Multibody Dynamic Model for the Estimation of Joint Loads and Muscle Forces

Computational musculoskeletal models have been developed to predict mechanical joint loads on the human spine, such as the forces and moments applied to vertebral and facet joints and the forces that act on ligaments and muscles because of difficulties in the direct measurement of joint loads. However, many whole-spine models lack certain elements. For example, the detailed facet joints in the cervical region or the whole spine region may not be implemented. In this study, a detailed cervico-thoraco-lumbar multibody musculoskeletal model with all major ligaments, separated structures of facet contact and intervertebral disk joints, and the rib cage was developed. The model was validated by comparing the intersegmental rotations, ligament tensile forces, facet joint contact forces, compressive and shear forces on disks, and muscle forces were to those reported in previous experimental and computational studies both by region (cervical, thoracic, or lumbar regions) and for the whole model. The comparisons demonstrated that our whole spine model is consistent with in vitro and in vivo experimental studies and with computational studies. The model developed in this study can be used in further studies to better understand spine structures and injury mechanisms of spinal disorders.