Next Article in Journal
Blockchain-Federated and Deep-Learning-Based Ensembling of Capsule Network with Incremental Extreme Learning Machines for Classification of COVID-19 Using CT Scans
Next Article in Special Issue
Gait Event Prediction Using Surface Electromyography in Parkinsonian Patients
Previous Article in Journal
Advances in Animal Models for Studying Bone Fracture Healing
Previous Article in Special Issue
Simultaneous Estimation of the Vertical Stiffness in the Knee and Hip for Healthy Human Subjects during Walking
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Bioengineering
Volume 10
Issue 2
10.3390/bioengineering10020202
bioengineering-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

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

by
Tanja Lerchl
1,2,*,
Kati Nispel
1,2,
Thomas Baum
2,
Jannis Bodden
2,
Veit Senner
1 and
Jan S. Kirschke
2
1
Sport Equipment and Sport Materials, School of Engineering and Design, Technical University of Munich, 85748 Garching, Germany
2
Department of Diagnostic and Interventional Neuroradiology, School of Medicine, Klinikum Rechts der Isar, Technical University of Munich, 81675 Munich, Germany
*
Author to whom correspondence should be addressed.
Bioengineering 2023, 10(2), 202; https://doi.org/10.3390/bioengineering10020202
Submission received: 29 December 2022 / Revised: 30 January 2023 / Accepted: 31 January 2023 / Published: 3 February 2023
(This article belongs to the Special Issue Biomechanics-Based Motion Analysis)
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

Chronic back pain is one of the major health issues worldwide. Though general risk factors such as occupation, obesity or anthropometric parameters could be identified in the past years [1], the specification of individual biomechanical indicators for the prediction of symptoms and chronicity is challenging, as it requires an in-depth knowledge of spinal kinematics and resulting loads. Even though experimental methods are essential to help build this knowledge, they come with limitations. In vitro studies can help understand segment mechanics but are not applicable when it comes to the investigation of complex in vivo biomechanics of the whole torso [2]. The invasive character of the in vivo measurement of these parameters via intradiscal pressure sensors [3,4] or instrumented vertebral implants [5,6] makes these methods unsuitable for clinical analysis. Computational, biomechanical models can provide a valuable alternative when it comes to the estimation of spinal loads. There are two approaches for the numerical analysis of spinal loading. While finite element models (FEM) hold the potential to investigate internal stress states in flexible bodies and their underlying or resulting deformation, multibody models (multibody system, MBS) can help analyze mechanical loads on the musculoskeletal system at a holistic level. Breaking the system down to its essential mechanical components, classic MBS models incorporate rigid bodies connected by joints and, depending on the respective research question, force elements representing flexible structures such as intervertebral discs (IVD), ligaments, cartilage, and other connective tissue. This way, MBS models represent a valuable tool to increase a profound understanding of healthy and pathological biomechanics. Gould et al. published a review on FEM and MBS models of the healthy and scoliotic spine in 2021 [7]. Focusing on the latter one, the authors state that their review provides solely a brief overview on MBS models of the healthy spine and refer the reader to the review on MBS modeling of the cervical spine by Alizadeh et al. [8] and the review by Dreischarf et al. on in vivo studies and computational models, published in 2016 [9].
The wide range of applications, improved technical capabilities, and increasing knowledge of spinal biomechanics, which answer old questions and raise new ones, mean that the demand for high-quality MBS models is not abating. As a consequence, the number of published models is increasing every year providing new opportunities and insight.
In recent years, models have been introduced that extend the classic notion of a multibody or musculoskeletal models. These models incorporate flexible bodies such as beam elements into rigid body models and thus soften the boundary between FEM and MBS models [10,11]. However, within the scope of this work, we want to review the developments in the field of multibody models of the healthy thoracolumbar spine, focusing on classical rigid body models. Hereby, we shed light on common modeling methods and applications, as well as identify and discuss related limitations and challenges in state-of-the-art spine modeling.

2. Methods

To generate a list of potentially relevant publications, a systematic search was carried out in PubMed and Scopus in November and December 2022. The search included the keywords “spine AND model AND ((multi AND body) OR musculoskeletal)”. Excluding results prior to 2013 left 1288 publications on PubMed and 1304 on Scopus. However, relevant citations in the articles were also included, if they were published before 2013. Subsequently, duplicates were removed by identical PubMedIDs and titles. Remaining articles were then filtered by title and abstract and the full text eventually analyzed. Publications were excluded if they featured at least one of the following topics:
  • Finite element modeling;
  • Models of the cervical spine;
  • Models without muscle incorporation;
  • Models of the scoliotic spine;
  • Models of the nonhuman spine;
  • Studies with a medical scope other than biomechanics.
Inclusion criteria were set to
  • Musculoskeletal models;
  • Multibody models;
  • Models of the thoracolumbar spine;
  • Models of the healthy spine.
We analyzed the remaining studies systematically according to the represented modeling methods and applications and identified existing limitations and challenges.

3. Multibody Modeling of the Healthy Spine

After filtering a total of 2592 articles, 81 articles remained, which were included in this review. Focusing on extensive musculoskeletal models of the thoracolumbar spine, we discuss models with reduced complexity, such as abstracted models [12,13,14,15,16], skeletal models neglecting muscular effects [17,18] or models of the lumbar spine [19,20,21,22,23,24,25,26,27,28,29] only in passing.
Overall, our literature review revealed that a large proportion of published studies was based on a few original models [30,31,32,33]. Due to the accessibility of these models via the commercially available software AnyBody (AnyBody Technology A/S, Aalborg, Denmark) [30,33] or the open-source software OpenSim [31,32,34], numerous studies can be found that used, modified, and extended these models, beyond the boundaries of the respective research groups as well [35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58]. Apart from these widely reused models, further original models can be found in the literature using alternative software [59,60,61,62,63,64]. Table 1 provides an overview of the original models found and subsequent studies associated with them.

3.1. General Model Setup and Kinematics

In the past two decades, simplified models of the whole torso with a detailed lumbar spine were developed to investigate lumbar loads [30,31,59,61,69]. One of the first generic models for lumbar load estimation was introduced by de Zee et al. in 2007 [30], which comprised seven rigid bodies for the pelvis including the sacrum, five lumbar vertebrae, and one lumped segment representing the thoracic spine including the rib cage and cervical spine. The model anatomy was based on publications by Hansen et al. [78] and Bodguk et al. [79]. De Zee defined intervertebral joints as spherical joints with their respective center of rotation (COR) located in the intersection of the instantaneous axis of rotation and the midsagittal plane according to Pearcy and Bodguk [80]. A total of 154 actuators representing muscle fascicles for the erector spinae (ES), rectus abdominis (RA), internal obliques (IO), external obliques (EO), psoas major (PM), quadratus lumborum (QL), and multifidus (MF) were implemented in the model either as a straight line between insertion and origin or, in order to mimic more realistic lines of action, redirected using so-called via points or wrapping surfaces [30].
Inspired by de Zee’s model, Christophy et al. published a generic multibody model of the lumbar spine in 2012 [31], incorporating a more detailed muscle architecture regarding the latissimus dorsi (LD) and the MF muscle. Using the open-source software OpenSim [34], the model has been widely used and extended in the past years [31,37,40,41,48,49,50,51,52,81,82]. In recent years, other models with simplified thorax have been published [59,61,69].
Favier et al. published a full-body model with a detailed lumbar spine in 2021 [69]. The model was created in OpenSim and included in total 20 rigid bodies including the head–neck, three-segment thoracic and cervical spine (spherical joints in T7-T8 and C7-T1), five lumbar vertebrae, pelvis with sacrum, as well as upper and lower extremities. The model incorporated a total of 538 muscle actuators for the lower limbs and lumbar spine [69].
Lerchl et al. introduced a pipeline for the semiautomated generation of individualized MBS models with a detailed lumbar spine created in the commercial multibody modeling software Simpack (Dassault Systèmes, France) in 2022 [59]. Based on CT data, the models included individual vertebrae T1-L5 with a fused thoracic part and rib cage and spherical lumbar intervertebral joints, and generic segments for the head–neck, pelvis, sacrum, and simplified arms. A total number of 103 actuators representing the muscles of the lower back were incorporated [59].
Research devoted to the loading of the thoracic spine is less common and therefore, only few models incorporating a detailed thoracic spine and rib cage can be found in the literature [32,33,72]. As opposed to musculoskeletal models with a rigid thorax, these models allow a comprehensive analysis of spinal loading for load cases involving thoracic movement. Based on the generic model of the lumbar spine by de Zee et al. [30], Ignasiak et al. introduced a musculoskeletal model of the thoracolumbar spine with a detailed articulated rib cage [33]. Ignasiak et al. extended the model by individual rigid bodies of 12 vertebrae, 10 pairs of ribs, and a sternum. Intervertebral thoracic joints were defined as six-DOF joints and lumbar joints, originally modeled as spherical joints [30], were also modified, respectively. Costovertebral (CV) and costotransverse (CT) joints were defined as revolute joints with the rotation axis in the frontal direction and all joints between the ribs and the sternum were modeled with six DOFs, except the first pair, which were modeled as spherical joints. The model was validated against in vivo data and used in follow-up studies [33,39,67,68].
A comprehensive model of the upper body including 60 segments (vertebrae, ribs, skull, sternum, hyoid, thyrohyoid, clavicles, scapulas, humeri, sacrum, and pelvis) created in AnyBody was published by Bayoglu et al. in 2019 [72].
Based on the lumbar spine model of Christophy et al. [31], Bruno et al. developed and validated a fully articulated model of the thoracolumbar spine in OPENSIM including individual vertebrae, ribs, and sternum [32]. Like Christophy’s model, the thoracolumbar model of Bruno et al. has been widely used and adapted since its publication [32,43,54,56,57,58,83,84].
In biomechanical MBS modeling, intersegmental connections are usually implemented as joints with defined DOFs, which can either be defined directly in the joint or are implemented as constraints, limiting the joint’s effective degrees of freedom to its relevant components. It is common practice to model intervertebral joints as spherical joints allowing rotation around three spatial axes [31,62]. Few models exist, that defined intervertebral joints with six DOFs, additionally accounting for translational motion [33,37,41,61]. The centers of rotation are located either in the geometrical center of the IVD [33,59,62] or in the instantaneous axis of rotation according to Pearcy and Bodguk [30,31,69,80]. CV joints are modeled as pin joints rotating around the vector between the costovertebral and costotransverse joints [32,33,72] or spherical joints [64] and CS joints as six DOFs [33,64,72]. Depending on the simulation approach (Section 3.4), kinematic data have been most commonly assigned to the respective DOFs according to findings from our own experimental studies or the literature (Section 3.3). This way, model kinematics are usually described using relative minimum coordinates. However, for inverse kinematic approaches, absolute coordinates are assigned to the end link of the kinematic chain. Providing stable boundary conditions for the mechanical analyses, the models are usually connected to the inertial frame of reference and therefore leaving the head–neck complex as the end link of the open kinematic chain. Upper-body weight is either combined and included in the center of mass of the lumped thoracic body [61], distributed according to the literature [85,86] or derived from patient-specific CT or MRI data and distributed levelwise along the thoracolumbar spine [59,62].

3.2. Passive (Visco)elastic Components

Various approaches have been taken regarding the modeling of viscoelastic structures that passively stabilize the spine, such as IVD, spinal ligaments, or the (cartilage) tissue of the thorax. The modeling approach can vary both in the level of detail and in the mechanical characteristics considered. Thus, some models neglect the effects of these components entirely [30,31,32,62], whereas others combine them partially or completely into one single stabilizing element per joint [60,69,72], or even integrate individual components explicitly [59,61,64,77]. The majority of approaches simplify the mechanical properties of connective tissue to linear elastic force elements, which produce corresponding forces and moments exclusively depending on their deformation. In multibody models, such material behavior is described via spring elements with constant stiffness for the corresponding DOFs. Only a few models incorporate the nonlinear mechanical behavior of biological passive structures [87]. However, modeling these components as purely elastic does not account for viscous effects that influence the mechanical response as a function of the deformation rate, also known as damping behavior. A detailed nonlinear viscoelastic modeling of IVDs and spinal ligaments, such as the anterior and posterior longitudinal ligament, the flavum ligament, as well as the interspinal and supraspinal ligament, can be found in only a few models [59,61]. The respective parameters are usually taken from in vitro studies available in the literature [88,89,90,91,92].
To examine thoracic loads, models require an appropriate force transmission from the rib cage to the thoracic spine in addition to intervertebral passive structures. In this context, costosternal (CS), costotransverse (CT) and costovertebral (CV) articulations are a central issue. Commonly, these connections are constrained and modeled as linear elastic elements according to the resulting DOFs. Stiffness parameters are usually taken from in vitro studies or adapted from previously published in silico studies. Bruno et al. included point-to-point actuators, which were placed between the ends of the ribs and the sternum (ribs 1–7) or between the ends of adjacent ribs (ribs 8–10) to represent forces transmitted by costal cartilage. As a result of a sensitivity analysis, forces generated by the actuators were set to 1000 N allowing the costal cartilage to provide a high supporting force to the end of the ribs [32].
Mechanical properties are usually incorporated either directly from mechanical testing, such as ligament tensile tests [88,93] or by simulating in vitro protocols, such as stepwise reduction studies, where individual connective structures are gradually removed from functional biological units, such as the FSU or the rib cage, while measuring the mechanical properties of the units after every resection [89,94,95]. However, due to the high level of intra- and interindividual variability regarding the mechanical characteristics of biological materials, the resulting parameters usually come with high standard deviations [96].

3.3. Scaling and Individualization

Spinal loading is highly dependent on a variety of subject-specific characteristics, such as spinal alignment, anthropometry, body weight distribution, or kinematics. While finite element models exist that account for individual characteristics [97,98,99,100,101,102,103,104], multibody models are predominantly generic in nature. In the past years, an increasing number of studies have been published, putting an emphasis on the individualization of the models [54,55,58,59,62].
A wide range of MBS models are based on measurements available in respective databases, e.g., in the OpenSim database (https://simtk.org/projects/osimdatabase, accessed on 27 December 2022). To gain reliable insights for the examined load cases, it is important to match the subject characteristics to the investigated kinematics as congruently as possible. It is common scientific practice to use available data based on measurements of bone geometries derived from imaging data or cadaver studies of individuals and scale and adapt the relevant parameters to the desired anthropometry depending on the characteristics of the studied target group. The need to make use of various sources in this regard makes it essential to be clear about the underlying data sets, in order to draw meaningful conclusions from simulation results. Thus, segment masses and body weight distribution and simplified kinematics are usually taken from the literature [85,86,105]. Some studies include experimental data collection of kinematics to scale the existing model appropriately [45,51,83] and include muscle activity from electromyography (EMG) measurement to drive the model [52]. This usually does not incorporate individual bone geometries, muscle morphology, or the mechanical properties of viscoelastic components.
However, the neglect or only limited consideration of interindividual variation makes these models poorly suitable for a detailed subject-specific analysis. Models based on coherent datasets regarding bone geometry, anthropometry, and muscle architecture, and kinematics are rare in the literature. Bayoglu et al. built a model based on extensive measurements of one cadaver, incorporating general kinematic data from the literature [72,73,74].
Dao et al. published a patient-specific model based on CT and MRI data [20] of the lumbar spine. Bruno et al. used their generic model [32] for the investigation of the impact of the integration of subject-specific properties [42]. Therefore, they incorporated CT-based measurements of trunk anatomy, such as spinal alignment and muscle morphology, indicating the relevance of considering these factors [42]. Based on this publication, Banks et al. investigated lumbar load in a patient-specific MBS model using CT data and marker-based motion capturing to combine individual musculoskeletal geometry and coherent kinematics [58]. However, the individualization of those models usually involves a time-consuming, manual, or semiautomated process which requires expert knowledge. To the best of our knowledge, only two publications can be found that deal with the topic of automating the individualization of MBS models [59,62].
Fasser at al. used annotated bi-planar radiography images (EOS imaging, Paris, France) for the automated generation of semi-subject-specific MBS models of the torso. The models included individual size and the alignment of bony structures as well as an individual body mass distribution. In the process, 112 and 109 points were marked in the frontal and sagittal plane, respectively, and converted into 3D coordinates. The body mass distribution was determined using the individual body contour of the imaging data. Individual bone geometries, muscle morphology, and passive elements were not included in the model [62].
Based on the use of artificial neural networks (ANN), Lerchl et al. introduced a pipeline for the automated segmentation of vertebrae [106] and soft tissue of the torso, as well as the generation of the points of interest defining muscles and ligaments’ attachment points and the location and orientation of intervertebral joints. All data were derived from CT imaging and the model generation required minimal manual interaction, making it suitable for the analysis of large patient cohorts. However, the individual characteristics of the muscles and connective tissue could not yet be integrated in the process [59].

3.4. Muscle Force Estimation

A mechanical analysis with multibody systems can follow two approaches, which define the necessary input data. Forward dynamic simulations (FD) require kinetic data to drive the model to generate specific kinematics. This usually means that muscle forces are applied directly or indirectly to the model to produce a desired motion. This is contrasted with the idea of inverse dynamic simulations(ID), which use kinematic data as input to calculate the required kinetic data. Thus, joint kinematics during as specific movement is imposed to the model and necessary joint moments and therefore, associated muscle forces are calculated. However, having more control variables, namely, muscle fascicles, than DOFs, the human musculoskeletal system is redundant. This leads to an infinite number of solutions for each load case. In order to determine the most suitable solutions, a mathematical optimization is a commonly used method. Numerous algorithms are available to find the optimal solution. Hereby, depending on the chosen algorithm, control variables, namely, muscle activation, excitation, or forces are varied in a deterministic or stochastic way until some given optimality criteria and constraints are met. Most commonly, a combination of inverse dynamics and static optimization (SO) is used [30,32,45], sometimes including inverse kinematics (IK) to determine individual joint kinematics [62,63,69]. The inverse dynamic simulation provides joint moments necessary to generate the simulated movement. Subsequently, the static optimization solves the redundancy problem for each time frame sequentially under the consideration of meeting equilibrium conditions.
In MBS models of the spine, muscles of interest are usually modeled as multiple fascicles, which comprehensively consider the respective lines of action (Section 3.1). Individual fascicles are modeled either as simple force actuators or, more complex, as Hill type muscles [107]. The classic muscle model according to Hill comprises serial and parallel elastic elements, representing passive elastic properties of the muscle–tendon complex as well as a contractile element representing the active component, namely, the function of myofilaments. This element can include muscle-specific characteristics, such as the force–length and force–velocity relationship as well as activation dynamics. Depending on how far these dynamics are taken into account, the muscle excitation, activation, or force can drive the model and therefore represent control variables for optimization routines. Detailed definitions of muscle-specific dynamics can be found in the literature [108,109].

4. Applications of MBS Models

MBS models can be used to address a wide range of questions. There are numerous publications devoted to the evaluation of methods in numerical modeling, including sensitivity analyss or validation studies. Furthermore, validated models can help to gain valuable insights into biomechanically or clinically relevant load cases. However, depending on the investigated load case and subject collective, model extensions, and modifications are usually necessary. Table 2 provides an overview of the most relevant studies using existing models to address specific research questions.

4.1. Studies with Methodological Focus

Various publications can be found in the literature evaluating and validating new approaches in MBS modeling [19,30,31,32,45,63,64,69]. For the purpose of validating these approaches, it is common scientific practice to compare simulation results with existing results from in vivo or in vitro measurements. Of note, those comparisons are mainly relative, as few in vivo measurements are available and exact boundary conditions are hard to control. Frequently used in vivo studies to validate results on spinal loading from simulation are intradiscal pressure measurements [4,114]. Estimated muscle forces are usually compared to EMG measurements from one’s own experimental studies [48] or the literature [59].
Apart from evaluating the validity of the modeling approach, the simulation results of generated MBS models can be used to validate novel methods in data processing regarding the derivation of both relevant modeling data from imaging [19,20,21] and kinematic data motion capture [54]. Due to the usually extensive effort connected to the processing of individual data, recent publications have focused on the automation of the process [59,62].
Simplifications are an integral part of any model and have to be taken into consideration when it comes to the interpretation of the results. To understand and evaluate their influence, MBS models have been used to systematically investigate common assumptions, such as the reduction of complex mechanics of the functional spine unit (FSU) [37,115]. Further, the sensitivity of the model accuracy to assumed positions of intervertebral centers of rotation [23,36] or muscle insertions [75] have been analyzed. Rockenfeller et al. investigated the effect of muscle- or torque-driven centrodes using an MBS model of the lumbar spine.
Furthermore, a systematic model-based analysis can help standardize clinical procedures, such as the classification of spinal shapes [116] or to define boundary conditions for experimental protocols [24].

4.2. Studies with Biomechanical or Clinical Focus

Validated models are used to comprehensively investigate biomechanical and clinical aspects of a wide range from routine scenarios to nonphysiological, or even traumatic events.
The relevance of low-dynamic everyday or work-related activities for the general population, as well as their experimental accessibility, make these scenarios among the most studied in biomechanical simulations. Therefore, numerous models exist that deal with the mechanical effects of lifting [12,13,25,46,76,77,82], everyday activities such as walking, flexion, extension, or lateral bending [15,43,69] or work-related situations such as high-frequency axial loading [17,18]. In this context, different lifting techniques were evaluated [50,51,83,117]. Accident situations were investigated by Wei et al. [16] for snowboarding and for frontal impact by Valdano et al. [14]. Incorporating noncritical higher dynamics, Raabe et al. combined a generic model of the lumbar spine [31] with a model of the lower limbs [111] to analyze the biomechanics of jogging [40]. Studies investigating specific kinematic boundary conditions usually involve an experimental setup to collect kinematic data in a healthy adult population [46,47,52,58,83]. Comparably few studies target more vulnerable populations, such as amputees [48,53] or children [27,56], who used validated models of adults and scales them according to the literature to match the average anthropometric data of children.
Regarding the influence of healthy anatomical and anthropometric and anatomical characteristics, biomechanical modeling have been used to determine the effect of spinal alignment [28,43,47], to gain insight into load sharing of passive structures of the FSU [22], the effect of ligament stiffness [65] or muscle strengthening [118].
Furthermore, MBS models can help to understand and treat pathological developing or surgically induced pathological biomechanics. Kuai et al. analyzed the impact of disc herniation on the kinetics of the spine and lower extremities during everyday activities [44].
Surgical interventions always represent a major intervention in the natural biomechanics of the musculoskeletal system. Thus, several studies on the effects of spinal fusion can be found in the literature [29,71,119]. The resulting kinematic effects of spinal fusion were investigated by Ignasiak et al., who proposed a method for the prediction of a full-body sagittal alignment including reciprocal changes as a reaction to spinal fusion [68].

5. Limitations and Challenges

It is in the nature of numerical models that they come with limitations. One of the great challenges is to keep the balance between necessary accuracy and reasonable complexity. This requires not only in-depth knowledge of the object to be modeled but also the corresponding data from experimental studies and the appropriate technical solutions for implementation. During our literature research, we were able to identify several core limitations that could be found in a wide range of MBS models of the spine and the related challenges when it came to addressing these limitations.

5.1. Database

Any model can only be as good as its input data. In the context of biomechanical models, this comprises bony geometry, anthropometry, muscle architecture, the mechanical parameters of viscoelastic components and kinematic data. Due to the necessary measurements to determine these parameters, it is currently not possible to build models based on fully consistent datasets. While anthropometric and kinematic data can be determined via noninvasive measures in biomechanics labs, such as marker-based motion capturing, the derivation of bony geometries, muscle architecture, and a detailed distribution of soft tissue usually need medical imaging or are performed in cadaver studies. However, the mechanical properties of viscoelastic components such as ligaments or the IVD can currently only be determined with the help of in vitro studies, which require the isolation of the structure of interest to mount them in respective testing machines. Consequently, these measurements are also usually performed with specimens from cadaver studies and highly dependent on the experimental conditions.
In the past years, more studies including widely individualized models were published [55,59,62]. However, even these models can only offer a limited customization.
In order to obtain consistent data sets for biomechanical models, alternative, noninvasive methods must be developed to determine these parameters in large subject cohorts. Here, the combination of experimental studies, multimodal imaging, and ANNs could be a possible solution to increase the level of model individualization beyond its anthropometric and skeletal characteristics. Thus, the individual mechanical condition of functional components can be evaluated partly on the basis of imaging data. For instance, according to the Pfirrmann scale, a potential degradation of the IVD can be determined via the height and signal intensity from MRI data [120]. Correlating this degradation with the mechanical alteration of IVD [121], this can be used to consider the individual mechanical state of connective tissue, when it is implemented in respective models. Training ANNs with these data will provide large, more diverse datasets for individualized multibody models.
Furthermore, invasive experimental studies on spinal loading for model validation are rare and are not widely feasible due to ethical reasons. Accordingly, even consistently constructed models cannot ultimately be validated against data pertaining to the individual in question. Additionally, the high level of variability in mechanical properties of biological materials as mentioned in Section 3.2, and therefore, the integration of parameters with high standard deviations inevitably leads to models containing inaccuracies. Depending on the complexity of the model, these inaccuracies can accumulate and further blur the generated results. It is necessary to be aware of existing inconsistencies and imprecision when interpreting simulation results in order not to draw incorrect conclusions.

5.2. Joint Definition

Intervertebral connections are a complex combination of the IVD, ligaments, facet joints, and articulated capsules. Depending on the applied load, this leads to complicated kinematics in which the instantaneous center of rotation migrates in the course of the motion [122]. However, in the vast majority of spine models, intervertebral joints are simplified to spherical joints allowing three rotational DOFs around a fixed center of rotation. The sensitivity of this assumption has been subject to several in silico studies [23,113,123], indicating that the effect of this assumption on the calculated muscle forces and spinal loading should not be neglected. Detailed modeling requires six degrees of freedom and the consideration of appropriate stabilizing structures, the validity of which depends primarily on the definition of their mechanical parameters (Section 5.1). There are some models to be found in the literature incorporating such detailed representation of intervertebral connection [22], mainly focusing on load sharing in passive structures.
Larger data sets could also help to better understand intervertebral dynamics in order to develop corresponding valid modeling approaches. As already mentioned in Section 5.1, the combination of imaging, machine learning for process automation, and in vitro studies can contribute to progress.

5.3. Intra-Abdominal Pressure

The stabilizing influence of intra-abdominal pressure (IAP) on the spine has been widely studied [124,125]. However, only a few MBS models consider its effects [38,60,63,70,77]. In consequence, spinal loads in lifting tasks or the inclination of the upper body are assumed to be overestimated in the MBS modeling of the spine. Arshad et al. observed a decrease of up to 514 N in lumbar compression force and 279 N in global muscle force due to the inclusion of intra-abdominal pressure [38]. These results indicated that it was necessary to consider the effects of IAP to obtain reliable quantitative results on spinal loads.

5.4. Muscle Modeling and Muscle Force Estimation

A valid representation of relevant muscles is crucial to gain meaningful findings on the biomechanics of the spine. Most of the models contain a detailed muscle architecture consisting of multiple fascicles spanning between origin and insertion according to the literature. Deploying modeling components, that are usually defined as point-to-point force elements, can lead to nonphysiological lever arms depending on the imposed movement. De Zee’s model used so-called via points to redirect the lines of action of the modeled long muscle fascicles along the rib cage and thus create more realistic lines of action compared to simple straight lines [30]. However, this approach came with an increased computational cost, making it only conditionally suited for a systematic analysis of large participant cohorts.
Another aspect that has to be critically discussed is the applied muscle model. While simple force actuators are considered sufficient for a static investigation, high-dynamic load situations require the consideration of activation and contraction dynamics. This requires an in-depth knowledge of the characteristics of individual muscle morphology such as optimal fiber length, physiological cross-sectional area (PCSA), or pennation angle. Again, the need for subject-specific solutions is evident, as muscle morphology is highly dependent on the individual.
The vast majority of currently published models use a combination of inverse dynamics and static optimization for muscle force calculation. This approach provides a sufficient accuracy in static and quasi-static simulations but is dependent on the defined cost function, constraints, and used algorithm. Most commonly used are criteria for minimum fatigue [126], or the sum of squared muscle strength [127] or activation [34], and the maximum muscle stress is defined as the upper-bound constraint, which is usually set to 100 MPa [32,49,59] to guarantee that equilibrium conditions are met reliably. However, this value does not correspond to a physiological value [49]. Furthermore, SO neglects cocontraction, which incorporates the activation of the antagonist in addition to the agonist stabilizing the respective joint and therefore increasing muscle activation. This is in contradiction to the idea of static optimization, which aims at minimizing the defined cost function (e.g., muscle activation) [128]. In high-dynamic load cases, where the role of cocontraction is more evident, this leads to an underestimation of spinal loading.
One way to address this problem is to use dynamic optimization (DO). In contrast to static optimization, the entire time history of the motion under investigation is taken into account [128]. Integrating the respective criteria in the optimization objective, stabilizing effects such as cocontraction can come into play [25]. However, this method comes with a massive increase of computational cost [129]. Another possibility would be to train models with the help of artificial intelligence. However, such training requires large quantities of data, which is not possible due to the still widely manual and therefore time-consuming process of modeling [128]. Anderson et al. compared both approaches for the simulation of normal gait in 2001, stating that both provided equivalent results for low-dynamic simulations [129]. A similar comparison was made by Morrow et al. for wheelchair propulsion, noticing significant differences in estimated muscle activations [130]. Keeping in mind that wheelchair propulsion comprises higher dynamics than normal gait, these findings indicate that the validity of the chosen approach was largely dependent on the investigated load case.

6. Conclusions

Multibody models are a powerful tool to gain insight into the healthy and pathological musculoskeletal system. They can promote a general understanding of the patho-biomechanics of a large set of medical impairments and might even be able to support diagnostics and therapy planning in the future. Although simplifications and assumptions are an integral part of any model, it is essential to look closely at the implications of these assumptions, potential interactions, and possible solutions. Modern technology holds the potential to provide some of these solutions. Thus, artificial intelligence and state-of-the-art medical imaging can provide the necessary extensive data basis to systematically investigate critical parameters to derive appropriate solutions. These technical approaches coupled with a distinct awareness of existing limitations will lead us towards a growing, more profound understanding of musculoskeletal mechanics.

Author Contributions

Conceptualization, T.L. and J.S.K.; writing—original draft preparation, T.L.; writing—review and editing, K.N., T.B., J.B., V.S. and J.S.K.; supervision, J.S.K.; project administration, J.S.K.; funding acquisition, J.S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program. Grant no.: 101045128—iBack-epic—ERC-2021-COG.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

J.S.K. is a co-founder of Bonescreen GmbH. All other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MBSMultibody system
FEMFinite element method
DOFDegree of freedom
FSUFunctional spine unit
IAPIntra-abdominal pressure
EMGElectromyography
CORCenter of rotation
IVDIntervertebral disc
IVIntervertebral
CSCostosternal
CVCostovertebral
CTCostotransversal
FDForward dynamic
IDInverse dynamic
IKInverse kinematic
SOStatic optimization
DODynamic optimization
ANNArtificial neural network
ALEArbitrary Langrangian–Eulerian

References

  1. Murtezani, A.; Ibraimi, Z.; Sllamniku, S.; Osmani, T.; Sherifi, S. Prevalence and risk factors for low back pain in industrial workers. Folia Med. 2011, 53, 68–74. [Google Scholar] [CrossRef]
  2. Fu, L.; Ma, J.; Lu, B.; Jia, H.; Zhao, J.; Kuang, M.; Feng, R.; Xu, L.; Bai, H.; Sun, L.; et al. Biomechanical effect of interspinous process distraction height after lumbar fixation surgery: An in vitro model. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2017, 231, 663–672. [Google Scholar] [CrossRef]
  3. Sato, K.; Kikuchi, S.; Yonezawa, T. In vivo intradiscal pressure measurement in healthy individuals and in patients with ongoing back problems. Spine 1999, 24, 2468. [Google Scholar] [CrossRef]
  4. Wilke, H.J.; Neef, P.; Hinz, B.; Seidel, H.; Claes, L. Intradiscal pressure together with anthropometric data—A data set for the validation of models. Clin. Biomech. 2001, 16, S111–S126. [Google Scholar] [CrossRef] [PubMed]
  5. Dreischarf, M.; Rohlmann, A.; Graichen, F.; Bergmann, G.; Schmidt, H. In vivo loads on a vertebral body replacement during different lifting techniques. J. Biomech. 2016, 49, 890–895. [Google Scholar] [CrossRef] [PubMed]
  6. Rohlmann, A.; Graichen, F.; Kayser, R.; Bender, A.; Bergmann, G. Loads on a Telemeterized Vertebral Body Replacement Measured in Two Patients. Spine 2008, 33, 1170–1179. [Google Scholar] [CrossRef]
  7. Gould, S.L.; Cristofolini, L.; Davico, G.; Viceconti, M. Computational Modelling of the Scoliotic Spine: A Literature Review. Int. J. Numer. Methods Biomed. Eng. 2021, 37, e3503. [Google Scholar] [CrossRef] [PubMed]
  8. Alizadeh, M.; Knapik, G.G.; Mageswaran, P.; Mendel, E.; Bourekas, E.; Marras, W.S. Biomechanical musculoskeletal models of the cervical spine: A systematic literature review. Clin. Biomech. 2020, 71, 115–124. [Google Scholar] [CrossRef]
  9. Dreischarf, M.; Shirazi-Adl, A.; Arjmand, N.; Rohlmann, A.; Schmidt, H. Estimation of loads on human lumbar spine: A review of in vivo and computational model studies. J. Biomech. 2016, 49, 833–845. [Google Scholar] [CrossRef]
  10. Heidari, E.; Arjmand, N.; Kahrizi, S. Comparisons of lumbar spine loads and kinematics in healthy and non-specific low back pain individuals during unstable lifting activities. J. Biomech. 2022, 144, 111344. [Google Scholar] [CrossRef] [PubMed]
  11. Khoddam-Khorasani, P.; Arjmand, N.; Shirazi-Adl, A. Effect of changes in the lumbar posture in lifting on trunk muscle and spinal loads: A combined in vivo, musculoskeletal, and finite element model study. J. Biomech. 2020, 104, 109728. [Google Scholar] [CrossRef] [PubMed]
  12. Breloff, S.P.; Chou, L.S. Three-dimensional multi-segmented spine joint reaction forces during common workplace physical demands/activities of daily living. Biomed. Eng.-Appl. Basis Commun. 2017, 29, 1750025. [Google Scholar] [CrossRef]
  13. Zaman, R.; Xiang, Y.; Cruz, J.; Yang, J. Three-dimensional asymmetric maximum weight lifting prediction considering dynamic joint strength. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2021, 235, 437–446. [Google Scholar] [CrossRef] [PubMed]
  14. Valdano, M.; Asensio-Gil, J.M.; Jiménez-Octavio, J.R.; Cabello-Reyes, M.; Vasserot-Tolmos, R.; López-Valdés, F.J. Parametric Analysis of The Effect of CRS Seatback Angle in Dummy Measurements in Frontal Impacts. In Proceedings of the IRCOBI Conference 2022, Porto, Portugal, 14–16 September 2022; Volume 2022, pp. 519–531. [Google Scholar]
  15. Panero, E.; Digo, E.; Ferrarese, V.; Dimanico, U.; Gastaldi, L. Multi-segments kinematic model of the human spine during gait. In Proceedings of the 2021 IEEE International Symposium on Medical Measurements and Applications (MeMeA), Lausanne, Switzerland, 23–25 June 2021. [Google Scholar] [CrossRef]
  16. Wei, W.; Evin, M.; Bailly, N.; Arnoux, P.J. Biomechanical evaluation of Back injuries during typical snowboarding backward falls. Scand. J. Med. Sci. Sport. 2022, 1–11. [Google Scholar] [CrossRef]
  17. Valentini, P.P.; Pennestrì, E. An improved three-dimensional multibody model of the human spine for vibrational investigations. Multibody Syst. Dyn. 2016, 36, 363–375. [Google Scholar] [CrossRef]
  18. Low, L.; Newell, N.; Masouros, S. A Multibody Model of the Spine for Injury Prediction in High-Rate Vertical Loading. In Proceedings of the IRCOBI Conference 2022, Porto, Portugal, 14–16 September 2022. [Google Scholar]
  19. Dao, T.T.; Pouletaut, P.; Charleux, F.; Lazáry, Á.; Eltes, P.; Varga, P.P.; Tho, M.C.H.B. Estimation of patient specific lumbar spine muscle forces using multi-physical musculoskeletal model and dynamic MRI. In Knowledge and Systems Engineering; Springer International: Basel, Switzerland, 2014; pp. 411–422. [Google Scholar]
  20. Dao, T.T.; Pouletaut, P.; Charleux, F.; Lazáry, Á.; Eltes, P.; Varga, P.P.; Tho, M.C.H.B. Multimodal medical imaging (CT and dynamic MRI) data and computer-graphics multi-physical model for the estimation of patient specific lumbar spine muscle forces. Data Knowl. Eng. 2015, 96, 3–18. [Google Scholar] [CrossRef]
  21. Dao, T.T.; Pouletaut, P.; Lazáry, Á.; Tho, M.C.H.B. Multimodal Medical Imaging Fusion for Patient Specific Musculoskeletal Modeling of the Lumbar Spine System in Functional Posture. J. Med. Biol. Eng. 2017, 37, 739–749. [Google Scholar] [CrossRef]
  22. Abouhossein, A.; Weisse, B.; Ferguson, S.J. A multibody modelling approach to determine load sharing between passive elements of the lumbar spine. Comput. Methods Biomech. Biomed. Eng. 2011, 14, 527–537. [Google Scholar] [CrossRef]
  23. Abouhossein, A.; Weisse, B.; Ferguson, S.J. Quantifying the centre of rotation pattern in a multi-body model of the lumbar spine. Comput. Methods Biomech. Biomed. Eng. 2013, 16, 1362–1373. [Google Scholar] [CrossRef] [PubMed]
  24. Borrelli, S.; Putame, G.; Pascoletti, G.; Terzini, M.; Zanetti, E.M. In Silico Meta-Analysis of Boundary Conditions for Experimental Tests on the Lumbar Spine. Ann. Biomed. Eng. 2022, 50, 1243–1254. [Google Scholar] [CrossRef]
  25. Ghiasi, M.S.; Arjmand, N.; Boroushaki, M.; Farahmand, F. Investigation of trunk muscle activities during lifting using a multi-objective optimization-based model and intelligent optimization algorithms. Med Biol. Eng. Comput. 2016, 54, 431–440. [Google Scholar] [CrossRef] [PubMed]
  26. Bauer, S.; Hausen, U.; Gruber, K. Effects of individual spine curvatures—A comparative study with the help of computer modelling. Biomed. Tech. Biomed. Eng. 2012, 57 (Suppl. 1), 132–135. [Google Scholar] [CrossRef] [PubMed]
  27. Bauer, S.; Wasserhess, C.; Paulus, D. Quantification of loads on the lumbar spine of children with different body weight—A comparative study with the help of computer modelling. Biomed. Tech. 2014, 59, S913–S916. [Google Scholar]
  28. Müller, A.; Rockenfeller, R.; Damm, N.; Kosterhon, M.; Kantelhardt, S.R.; Aiyangar, A.K.; Gruber, K. Load Distribution in the Lumbar Spine During Modeled Compression Depends on Lordosis. Front. Bioeng. Biotechnol. 2021, 9, 661258. [Google Scholar] [CrossRef] [PubMed]
  29. Kantelhardt, S.; Hausen, U.; Kosterhon, M.; Amr, A.; Gruber, K.; Giese, A. Computer simulation and image guidance for individualised dynamic spinal stabilization. Int. J. Comput. Assist. Radiol. Surg. 2015, 10, 1325–1332. [Google Scholar] [CrossRef]
  30. de Zee, M.; Hansen, L.; Wong, C.; Rasmussen, J.; Simonsen, E.B. A generic detailed rigid-body lumbar spine model. J. Biomech. 2007, 40, 1219–1227. [Google Scholar] [CrossRef]
  31. Christophy, M.; Faruk Senan, N.A.; Lotz, J.C.; O’Reilly, O.M. A musculoskeletal model for the lumbar spine. Biomech. Model. Mechanobiol. 2012, 11, 19–34. [Google Scholar] [CrossRef]
  32. Bruno, A.G.; Bouxsein, M.L.; Anderson, D.E. Development and Validation of a Musculoskeletal Model of the Fully Articulated Thoracolumbar Spine and Rib Cage. J. Biomech. Eng. 2015, 137, 081003. [Google Scholar] [CrossRef]
  33. Ignasiak, D.; Dendorfer, S.; Ferguson, S.J. Thoracolumbar spine model with articulated ribcage for the prediction of dynamic spinal loading. J. Biomech. 2016, 49, 959–966. [Google Scholar] [CrossRef]
  34. Delp, S.L.; Anderson, F.C.; Arnold, A.S.; Loan, P.; Habib, A.; John, C.T.; Guendelman, E.; Thelen, D.G. OpenSim: Open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Bio-Med. Eng. 2007, 54, 1940–1950. [Google Scholar] [CrossRef] [Green Version]
  35. Han, K.S.; Zander, T.; Taylor, W.R.; Rohlmann, A. An enhanced and validated generic thoraco-lumbar spine model for prediction of muscle forces. Med. Eng. Phys. 2012, 34, 709–716. [Google Scholar] [CrossRef] [PubMed]
  36. Han, K.S.; Kim, K.; Park, W.M.; Lim, D.S.; Kim, Y.H. Effect of centers of rotation on spinal loads and muscle forces in total disk replacement of lumbar spine. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2013, 227, 543–550. [Google Scholar] [CrossRef]
  37. Meng, X.; Bruno, A.G.; Cheng, B.; Wang, W.; Bouxsein, M.L.; Anderson, D.E. Incorporating Six Degree-of-Freedom Intervertebral Joint Stiffness in a Lumbar Spine Musculoskeletal Model-Method and Performance in Flexed Postures. J. Biomech. Eng. 2015, 137, 101008. [Google Scholar] [CrossRef] [PubMed]
  38. Arshad, R.; Zander, T.; Dreischarf, M.; Schmidt, H. Influence of lumbar spine rhythms and intra-abdominal pressure on spinal loads and trunk muscle forces during upper body inclination. Med. Eng. Phys. 2016, 38, 333–338. [Google Scholar] [CrossRef] [PubMed]
  39. Ignasiak, D.; Ferguson, S.J.; Arjmand, N. A rigid thorax assumption affects model loading predictions at the upper but not lower lumbar levels. J. Biomech. 2016, 49, 3074–3078. [Google Scholar] [CrossRef]
  40. Raabe, M.E.; Chaudhari, A.M. An investigation of jogging biomechanics using the full-body lumbar spine model: Model development and validation. J. Biomech. 2016, 49, 1238–1243. [Google Scholar] [CrossRef]
  41. Senteler, M.; Weisse, B.; Rothenfluh, D.A.; Snedeker, J.G. Intervertebral reaction force prediction using an enhanced assembly of OpenSim models. Comput. Methods Biomech. Biomed. Eng. 2016, 19, 538–548. [Google Scholar] [CrossRef]
  42. Bruno, A.G.; Mokhtarzadeh, H.; Allaire, B.T.; Velie, K.R.; De Paolis Kaluza, M.C.; Anderson, D.E.; Bouxsein, M.L. Incorporation of CT-based measurements of trunk anatomy into subject-specific musculoskeletal models of the spine influences vertebral loading predictions. J. Orthop. Res. Off. Publ. Orthop. Res. Soc. 2017, 35, 2164–2173. [Google Scholar] [CrossRef]
  43. Bruno, A.G.; Burkhart, K.; Allaire, B.; Anderson, D.E.; Bouxsein, M.L. Spinal loading patterns from biomechanical modeling explain the high incidence of vertebral fractures in the thoracolumbar region. J. Bone Miner. Res. 2017, 32, 1282–1290. [Google Scholar] [CrossRef]
  44. Kuai, S.; Zhou, W.; Liao, Z.; Ji, R.; Guo, D.; Zhang, R.; Liu, W. Influences of lumbar disc herniation on the kinematics in multi-segmental spine, pelvis, and lower extremities during five activities of daily living. BMC Musculoskelet. Disord. 2017, 18, 1–13. [Google Scholar] [CrossRef] [Green Version]
  45. Bassani, T.; Stucovitz, E.; Qian, Z.; Briguglio, M.; Galbusera, F. Validation of the AnyBody full body musculoskeletal model in computing lumbar spine loads at L4L5 level. J. Biomech. 2017, 58, 89–96. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, X.Y.; Kim, H.K.; Zhang, Y. Development of an enhanced musculoskeletal model for simulating lumbar spine loading during manual lifting tasks. Lect. Notes Comput. Sci. (Incl. Subser. Lect. Notes Artif. Intell. Lect. Notes Bioinform.) 2017, 10286, 229–237. [Google Scholar] [CrossRef]
  47. Bassani, T.; Casaroli, G.; Galbusera, F. Dependence of lumbar loads on spinopelvic sagittal alignment: An evaluation based on musculoskeletal modeling. PLoS ONE 2019, 14, e0207997. [Google Scholar] [CrossRef] [PubMed]
  48. Actis, J.A.; Honegger, J.D.; Gates, D.H.; Petrella, A.J.; Nolasco, L.A.; Silverman, A.K. Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine. J. Biomech. 2018, 68, 107–114. [Google Scholar] [CrossRef] [PubMed]
  49. Beaucage-Gauvreau, E.; Robertson, W.S.P.; Brandon, S.C.E.; Fraser, R.; Freeman, B.J.C.; Graham, R.B.; Thewlis, D.; Jones, C.F. Validation of an OpenSim full-body model with detailed lumbar spine for estimating lower lumbar spine loads during symmetric and asymmetric lifting tasks. Comput. Methods Biomech. Biomed. Eng. 2019, 22, 451–464. [Google Scholar] [CrossRef]
  50. Beaucage-Gauvreau, E.; Brandon, S.C.; Robertson, W.S.; Fraser, R.; Freeman, B.J.; Graham, R.B.; Thewlis, D.; Jones, C.F. A braced arm-to-thigh (BATT) lifting technique reduces lumbar spine loads in healthy and low back pain participants. J. Biomech. 2020, 100, 109584. [Google Scholar] [CrossRef]
  51. Beaucage-Gauvreau, E.; Brandon, S.C.; Robertson, W.S.; Fraser, R.; Freeman, B.J.; Graham, R.B.; Thewlis, D.; Jones, C.F. Lumbar spine loads are reduced for activities of daily living when using a braced arm-to-thigh technique. Eur. Spine J. 2021, 30, 1035–1042. [Google Scholar] [CrossRef]
  52. Molinaro, D.D.; King, A.S.; Young, A.J. Biomechanical analysis of common solid waste collection throwing techniques using OpenSim and an EMG-assisted solver. J. Biomech. 2020, 104, 109704. [Google Scholar] [CrossRef]
  53. Honegger, J.D.; Actis, J.A.; Gates, D.H.; Silverman, A.K.; Munson, A.H.; Petrella, A.J. Development of a multiscale model of the human lumbar spine for investigation of tissue loads in people with and without a transtibial amputation during sit-to-stand. Biomech. Model. Mechanobiol. 2021, 20, 339–358. [Google Scholar] [CrossRef]
  54. Burkhart, K.; Grindle, D.; Bouxsein, M.L.; Anderson, D.E. Between-session reliability of subject-specific musculoskeletal models of the spine derived from optoelectronic motion capture data. J. Biomech. 2020, 112, 110044. [Google Scholar] [CrossRef]
  55. Overbergh, T.; Severijns, P.; Beaucage-Gauvreau, E.; Jonkers, I.; Moke, L.; Scheys, L. Development and validation of a modeling workflow for the generation of image-based, subject-specific thoracolumbar models of spinal deformity. J. Biomech. 2020, 110, 109946. [Google Scholar] [CrossRef] [PubMed]
  56. Schmid, S.; Burkhart, K.A.; Allaire, B.T.; Grindle, D.; Anderson, D.E. Musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6–18 years. J. Biomech. 2020, 102, 109305. [Google Scholar] [CrossRef]
  57. Schmid, S.; Connolly, L.; Moschini, G.; Meier, M.L.; Senteler, M. Skin marker-based subject-specific spinal alignment modeling: A feasibility study. J. Biomech. 2022, 137, 111102. [Google Scholar] [CrossRef]
  58. Banks, J.J.; Alemi, M.M.; Allaire, B.T.; Lynch, A.C.; Bouxsein, M.L.; Anderson, D.E. Using static postures to estimate spinal loading during dynamic lifts with participant-specific thoracolumbar musculoskeletal models. Appl. Ergon. 2023, 106, 103869. [Google Scholar] [CrossRef] [PubMed]
  59. Lerchl, T.; El Husseini, M.; Bayat, A.; Sekuboyina, A.; Hermann, L.; Nispel, K.; Baum, T.; Löffler, M.T.; Senner, V.; Kirschke, J.S. Validation of a Patient-Specific Musculoskeletal Model for Lumbar Load Estimation Generated by an Automated Pipeline From Whole Body CT. Front. Bioeng. Biotechnol. 2022, 10, 862804. [Google Scholar] [CrossRef]
  60. Malakoutian, M.; Street, J.; Wilke, H.J.; Stavness, I.; Fels, S.; Oxland, T. A musculoskeletal model of the lumbar spine using ArtiSynth–development and validation. Comput. Methods Biomech. Biomed. Eng. Imaging Vis. 2018, 6, 483–490. [Google Scholar] [CrossRef]
  61. Rupp, T.K.; Ehlers, W.; Karajan, N.; Günther, M.; Schmitt, S. A forward dynamics simulation of human lumbar spine flexion predicting the load sharing of intervertebral discs, ligaments, and muscles. Biomech. Model. Mechanobiol. 2015, 14, 1081–1105. [Google Scholar] [CrossRef]
  62. Fasser, M.R.; Jokeit, M.; Kalthoff, M.; Gomez Romero, D.A.; Trache, T.; Snedeker, J.G.; Farshad, M.; Widmer, J. Subject-Specific Alignment and Mass Distribution in Musculoskeletal Models of the Lumbar Spine. Front. Bioeng. Biotechnol. 2021, 9, 745. [Google Scholar] [CrossRef]
  63. Huynh, K.; Gibson, I.; Jagdish, B.; Lu, W. Development and validation of a discretised multi-body spine model in LifeMOD for biodynamic behaviour simulation. Comput. Methods Biomech. Biomed. Eng. 2015, 18, 175–184. [Google Scholar] [CrossRef]
  64. Khurelbaatar, T.; Kim, K.; Kim, Y.H. A cervico-thoraco-lumbar multibody dynamic model for the estimation of joint loads and muscle forces. J. Biomech. Eng. 2015, 137, 111001. [Google Scholar] [CrossRef]
  65. Putzer, M.; Auer, S.; Malpica, W.; Suess, F.; Dendorfer, S. A numerical study to determine the effect of ligament stiffness on kinematics of the lumbar spine during flexion. BMC Musculoskelet. Disord. 2016, 17, 1–7. [Google Scholar] [CrossRef] [PubMed]
  66. Andersen, M.S.; Damsgaard, M.; Rasmussen, J. Force-dependent kinematics: A new analysis method for non-conforming joints. In Proceedings of the XIII International Symposium on Computer Simulation in Biomechanics, Leuven, Belgium, 30 June–2 July 2011. [Google Scholar]
  67. Ignasiak, D.; Valenzuela, W.; Reyes, M.; Ferguson, S.J. The effect of muscle ageing and sarcopenia on spinal segmental loads. Eur. Spine J. 2018, 27, 2650–2659. [Google Scholar] [CrossRef]
  68. Ignasiak, D. A novel method for prediction of postoperative global sagittal alignment based on full-body musculoskeletal modeling and posture optimization. J. Biomech. 2020, 102, 109324. [Google Scholar] [CrossRef]
  69. Favier, C.D.; Finnegan, M.E.; Quest, R.A.; Honeyfield, L.; McGregor, A.H.; Phillips, A.T.M. An open-source musculoskeletal model of the lumbar spine and lower limbs: A validation for movements of the lumbar spine. Comput. Methods Biomech. Biomed. Eng. 2021, 24, 1310–1325. [Google Scholar] [CrossRef] [PubMed]
  70. Malakoutian, M.; Sanchez, C.A.; Brown, S.H.; Street, J.; Fels, S.; Oxland, T.R. Biomechanical properties of paraspinal muscles influence spinal loading—A musculoskeletal simulation study. Front. Bioeng. Biotechnol. 2022, 10, 852201. [Google Scholar] [CrossRef] [PubMed]
  71. Fasser, M.R.; Gerber, G.; Passaplan, C.; Cornaz, F.; Snedeker, J.G.; Farshad, M.; Widmer, J. Computational model predicts risk of spinal screw loosening in patients. Eur. Spine J. 2022, 31, 2639–2649. [Google Scholar] [CrossRef]
  72. Bayoglu, R.; Galibarov, P.E.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente Spine Model: A thorough investigation of the spinal loads in a complete and coherent musculoskeletal model of the human spine. Med. Eng. Phys. 2019, 68, 35–45. [Google Scholar] [CrossRef]
  73. Bayoglu, R.; Geeraedts, L.; Groenen, K.H.J.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the lumbar region of the human spine. J. Biomech. 2017, 53, 111–119. [Google Scholar] [CrossRef]
  74. Bayoglu, R.; Geeraedts, L.; Groenen, K.H.J.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the thoracic and cervical regions of the human spine. J. Biomech. 2017, 58, 52–63. [Google Scholar] [CrossRef]
  75. Bayoglu, R.; Guldeniz, O.; Verdonschot, N.; Koopman, B.; Homminga, J. Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites. Comput. Methods Biomech. Biomed. Eng. 2019, 22, 1135–1143. [Google Scholar] [CrossRef] [Green Version]
  76. Huang, M.; Hajizadeh, K.; Gibson, I.; Lee, T. Analysis of compressive load on intervertebral joint in standing and sitting postures. Technol. Health Care 2016, 24, 215–223. [Google Scholar] [CrossRef] [PubMed]
  77. Guo, J.; Guo, W.; Ren, G. Embodiment of intra-abdominal pressure in a flexible multibody model of the trunk and the spinal unloading effects during static lifting tasks. Biomech. Model. Mechanobiol. 2021, 20, 1599–1626. [Google Scholar] [CrossRef] [PubMed]
  78. Hansen, L.; de Zee, M.; Rasmussen, J.; Andersen, T.B.; Wong, C.; Simonsen, E.B. Anatomy and biomechanics of the back muscles in the lumbar spine with reference to biomechanical modeling. Spine 2006, 31, 1888–1899. [Google Scholar] [CrossRef] [PubMed]
  79. Bogduk, N. Clinical Anatomy of the Lumbar Spine and Sacrum; Elsevier Health Sciences: Amsterdam, The Netherlands, 1997. [Google Scholar]
  80. Pearcy, M.J.; Bogduk, N. Instantaneous axes of rotation of the lumbar intervertebral joints. Spine 1988, 13, 1033–1041. [Google Scholar] [CrossRef]
  81. Byrne, R.M.; Aiyangar, A.K.; Zhang, X. Sensitivity of musculoskeletal model-based lumbar spinal loading estimates to type of kinematic input and passive stiffness properties. J. Biomech. 2020, 102, 109659. [Google Scholar] [CrossRef]
  82. Kim, H.K.; Zhang, Y. Estimation of lumbar spinal loading and trunk muscle forces during asymmetric lifting tasks: Application of whole-body musculoskeletal modelling in OpenSim. Ergonomics 2017, 60, 563–576. [Google Scholar] [CrossRef]
  83. von Arx, M.; Liechti, M.; Connolly, L.; Bangerter, C.; Meier, M.L.; Schmid, S. From Stoop to Squat: A comprehensive analysis of lumbar loading among different lifting styles. Front. Bioeng. Biotechnol. 2021, 9, 769117. [Google Scholar] [CrossRef]
  84. Wang, W.; Wang, D.; De Groote, F.; Scheys, L.; Jonkers, I. Implementation of physiological functional spinal units in a rigid-body model of the thoracolumbar spine. J. Biomech. 2020, 98, 109437. [Google Scholar] [CrossRef]
  85. Pearsall, D.J.; Reid, J.G.; Livingston, L.A. Segmental inertial parameters of the human trunk as determined from computed tomography. Ann. Biomed. Eng. 1996, 24, 198–210. [Google Scholar] [CrossRef]
  86. Winter, D.A. Biomechanics and Motor Control of Human Movement; Wiley: Weilheim, Germany, 2009. [Google Scholar]
  87. Fung, Y.C. Biomechanics: Mechanical Properties of Living Tissues; Springer Science & Business Media: Luxemburg, 2013. [Google Scholar]
  88. Pintar, F.A.; Yoganandan, N.; Myers, T.; Elhagediab, A.; Sances, A., Jr. Biomechanical properties of human lumbar spine ligaments. J. Biomech. 1992, 25, 1351–1356. [Google Scholar] [CrossRef]
  89. Heuer, F.; Schmidt, H.; Klezl, Z.; Claes, L.; Wilke, H.J. Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle. J. Biomech. 2007, 40, 271–280. [Google Scholar] [CrossRef] [PubMed]
  90. Ashton-Miller, J.A.; Schultz, A.B. Biomechanics of the human spine. Basic Orthop. Biomech. 1997, 2, 353–385. [Google Scholar]
  91. Panjabi, M.M.; Brand, R., Jr.; White, A., 3rd. Mechanical properties of the human thoracic spine as shown by three-dimensional load-displacement curves. JBJS 1976, 58, 642–652. [Google Scholar] [CrossRef]
  92. White, A.A. Clinical Biomechanics of the Spine; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2022. [Google Scholar]
  93. Myklebust, J.B.; Pintar, F.; Yoganandan, N.; Cusick, J.F.; Maiman, D.; Myers, T.J.; Sances, A., Jr. Tensile strength of spinal ligaments. Spine 1988, 13, 526–531. [Google Scholar] [CrossRef]
  94. Liebsch, C.; Graf, N.; Appelt, K.; Wilke, H.J. The rib cage stabilizes the human thoracic spine: An in vitro study using stepwise reduction of rib cage structures. PLoS ONE 2017, 12, e0178733. [Google Scholar] [CrossRef]
  95. Wilke, H.J.; Grundler, S.; Ottardi, C.; Mathew, C.E.; Schlager, B.; Liebsch, C. In vitro analysis of thoracic spinal motion segment flexibility during stepwise reduction of all functional structures. Eur. Spine J. 2020, 29, 179–185. [Google Scholar] [CrossRef]
  96. Cook, D.; Julias, M.; Nauman, E. Biological variability in biomechanical engineering research: Significance and meta-analysis of current modeling practices. J. Biomech. 2014, 47, 1241–1250. [Google Scholar] [CrossRef]
  97. Akhavanfar, M.H.; Kazemi, H.; Eskandari, A.H.; Arjmand, N. Obesity and spinal loads; a combined MR imaging and subject-specific modeling investigation. J. Biomech. 2018, 70, 102–112. [Google Scholar] [CrossRef]
  98. El Ouaaid, Z.; Shirazi-Adl, A.; Plamondon, A. Effects of variation in external pulling force magnitude, elevation, and orientation on trunk muscle forces, spinal loads and stability. J. Biomech. 2016, 49, 946–952. [Google Scholar] [CrossRef]
  99. Eskandari, A.H.; Arjmand, N.; Shirazi-Adl, A.; Farahmand, F. Hypersensitivity of trunk biomechanical model predictions to errors in image-based kinematics when using fully displacement-control techniques. J. Biomech. 2019, 84, 161–171. [Google Scholar] [CrossRef]
  100. Ghezelbash, F.; Shirazi-Adl, A.; Arjmand, N.; El-Ouaaid, Z.; Plamondon, A. Subject-specific biomechanics of trunk: Musculoskeletal scaling, internal loads and intradiscal pressure estimation. Biomech. Model. Mechanobiol. 2016, 15, 1699–1712. [Google Scholar] [CrossRef] [PubMed]
  101. Little, J.P.; Adam, C.J. Geometric sensitivity of patient-specific finite element models of the spine to variability in user-selected anatomical landmarks. Comput. Methods Biomech. Biomed. Eng. 2015, 18, 676–688. [Google Scholar] [CrossRef] [PubMed]
  102. Naserkhaki, S.; Jaremko, J.L.; El-Rich, M. Effects of inter-individual lumbar spine geometry variation on load-sharing: Geometrically personalized Finite Element study. J. Biomech. 2016, 49, 2909–2917. [Google Scholar] [CrossRef] [PubMed]
  103. Périé, D.; Sales De Gauzy, J.; Hobatho, M.C. Biomechanical evaluation of Cheneau-Toulouse-Munster brace in the treatment of scoliosis using optimisation approach and finite element method. Med. Biol. Eng. Comput. 2002, 40, 296–301. [Google Scholar] [CrossRef] [PubMed]
  104. Vergari, C.; Courtois, I.; Ebermeyer, E.; Bouloussa, H.; Vialle, R.; Skalli, W. Experimental validation of a patient-specific model of orthotic action in adolescent idiopathic scoliosis. Eur. Spine J. 2016, 25, 3049–3055. [Google Scholar] [CrossRef]
  105. Wong, K.W.N.; Luk, K.D.K.; Leong, J.C.Y.; Wong, S.F.; Wong, K.K.Y. Continuous Dynamic Spinal Motion Analysis. Spine 2006, 31, 414–419. [Google Scholar] [CrossRef]
  106. Sekuboyina, A.; Husseini, M.E.; Bayat, A.; Löffler, M.; Liebl, H.; Li, H.; Tetteh, G.; Kukačka, J.; Payer, C.; Štern, D. VerSe: A vertebrae labelling and segmentation benchmark for multi-detector CT images. arXiv 2020, arXiv:2001.09193. [Google Scholar] [CrossRef]
  107. Hill, A.V. The heat of shortening and the dynamic constants of muscle. Proc. R. Soc. Lond. Ser. B-Biol. Sci. 1938, 126, 136–195. [Google Scholar]
  108. Thelen, D.G. Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults. J. Biomech. Eng. 2003, 125, 70–77. [Google Scholar] [CrossRef]
  109. Millard, M.; Uchida, T.; Seth, A.; Delp, S.L. Flexing computational muscle: Modeling and simulation of musculotendon dynamics. J. Biomech. Eng. 2013, 135, 021005. [Google Scholar] [CrossRef]
  110. Delp, S.L.; Loan, J.P.; Hoy, M.G.; Zajac, F.E.; Topp, E.L.; Rosen, J.M. An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans. Biomed. Eng. 1990, 37, 757–767. [Google Scholar] [CrossRef]
  111. Hamner, S.R.; Seth, A.; Delp, S.L. Muscle contributions to propulsion and support during running. J. Biomech. 2010, 43, 2709–2716. [Google Scholar] [CrossRef] [Green Version]
  112. Anderson, F.C.; Pandy, M.G. A dynamic optimization solution for vertical jumping in three dimensions. Comput. Methods Biomech. Biomed. Eng. 1999, 2, 201–231. [Google Scholar] [CrossRef] [PubMed]
  113. Senteler, M.; Aiyangar, A.; Weisse, B.; Farshad, M.; Snedeker, J.G. Sensitivity of intervertebral joint forces to center of rotation location and trends along its migration path. J. Biomech. 2018, 70, 140–148. [Google Scholar] [CrossRef] [PubMed]
  114. Takahashi, I.; Kikuchi, S.i.; Sato, K.; Sato, N. Mechanical load of the lumbar spine during forward bending motion of the trunk-a biomechanical study. Spine 2006, 31, 18–23. [Google Scholar] [CrossRef] [PubMed]
  115. Wang, Q.D.; Guo, L.X. Biomechanical role of osteoporosis in the vibration characteristics of human spine after lumbar interbody fusion. Int. J. Numer. Methods Biomed. Eng. 2020, 36, e3402. [Google Scholar] [CrossRef]
  116. Rockenfeller, R.; Müller, A. Augmenting the Cobb angle: Three-dimensional analysis of whole spine shapes using Bézier curves. Comput. Methods Programs Biomed. 2022, 225, 107075. [Google Scholar] [CrossRef]
  117. Kim, J.W.; Eom, G.M.; Kwon, Y.R. Analysis of maximum joint moment during infant lifting-up motion. Technol. Health Care 2022, 30, S441–S450. [Google Scholar] [CrossRef]
  118. Nowakowska-Lipiec, K.; Michnik, R.; Linek, P.; Myśliwiec, A.; Jochymczyk-Woźniak, K.; Gzik, M. A numerical study to determine the effect of strengthening and weakening of the transversus abdominis muscle on lumbar spine loads. Comput. Methods Biomech. Biomed. Eng. 2020, 23, 1287–1296. [Google Scholar] [CrossRef]
  119. Bauer, S.; Paulus, D. Analysis of the biomechanical effects of spinal fusion to adjacent vertebral segments of the lumbar spine using multi body simulation. Int. J. Simul. Syst. Sci. Technol. 2014, 15, 1–7. [Google Scholar] [CrossRef]
  120. Pfirrmann, C.W.; Metzdorf, A.; Zanetti, M.; Hodler, J.; Boos, N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine 2001, 26, 1873–1878. [Google Scholar] [CrossRef]
  121. Foltz, M.H.; Kage, C.C.; Johnson, C.P.; Ellingson, A.M. Noninvasive assessment of biochemical and mechanical properties of lumbar discs through quantitative magnetic resonance imaging in asymptomatic volunteers. J. Biomech. Eng. 2017, 139, 111002. [Google Scholar] [CrossRef]
  122. Bogduk, N.; Macintosh, J.E.; Pearcy, M.J. A universal model of the lumbar back muscles in the upright position. Spine 1992, 17, 897–913. [Google Scholar] [CrossRef] [PubMed]
  123. Aiyangar, A.; Zheng, L.; Anderst, W.; Zhang, X. Instantaneous centers of rotation for lumbar segmental extension in vivo. J. Biomech. 2017, 52, 113–121. [Google Scholar] [CrossRef]
  124. Daggfeldt, K.; Thorstensson, A. The mechanics of back-extensor torque production about the lumbar spine. J. Biomech. 2003, 36, 815–825. [Google Scholar] [CrossRef] [PubMed]
  125. Hodges, P.W.; Cresswell, A.G.; Daggfeldt, K.; Thorstensson, A. In vivo measurement of the effect of intra-abdominal pressure on the human spine. J. Biomech. 2001, 34, 347–353. [Google Scholar] [CrossRef] [PubMed]
  126. Rasmussen, J.; Damsgaard, M.; Voigt, M. Muscle recruitment by the min/max criterion—A comparative numerical study. J. Biomech. 2001, 34, 409–415. [Google Scholar] [CrossRef]
  127. Crowninshield, R.D.; Brand, R.A. A physiologically based criterion of muscle force prediction in locomotion. J. Biomech. 1981, 14, 793–801. [Google Scholar] [CrossRef]
  128. Ezati, M.; Ghannadi, B.; McPhee, J. A review of simulation methods for human movement dynamics with emphasis on gait. Multibody Syst. Dyn. 2019, 47, 265–292. [Google Scholar] [CrossRef]
  129. Anderson, F.C.; Pandy, M.G. Static and dynamic optimization solutions for gait are practically equivalent. J. Biomech. 2001, 34, 153–161. [Google Scholar] [CrossRef]
  130. Morrow, M.M.; Rankin, J.W.; Neptune, R.R.; Kaufman, K.R. A comparison of static and dynamic optimization muscle force predictions during wheelchair propulsion. J. Biomech. 2014, 47, 3459–3465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Table
Table 1. Overview of original models of the musculoskeletal thoracolumbar spine and related modeling methods. Semi-individualized models are those that contain both individualized and generic musculoskeletal components. Joint definitions include potentially assigned constraints.
Table 1. Overview of original models of the musculoskeletal thoracolumbar spine and related modeling methods. Semi-individualized models are those that contain both individualized and generic musculoskeletal components. Joint definitions include potentially assigned constraints.
ReferenceIncluded SegmentsJoint DefinitionGeneric/Indiv.Passive Force ElementsMuscle Model and Force EstimationSoftwareRelated Studies
de Zee et al. [30]Pelvis, sacrum, L1-L5, thorax3 rot. DOFs (IV)Generic-Act., ID, SOAnyBody[33,35,36,39,44,45,47,65]
Christophy et al. [31]Pelvis, sacrumL1-L5, thorax3 rot. DOFs (IV)Generic-Hill typeOpenSim[37,40,41,46,48,49,50,51,52,53]
Bruno et al. [32]Pelvis, sacrumT1-L5, ribs, sternum, upper limbs, head–neck3 rot. DOFs (IV)1 rot. DOFs (CV)Generic-Hill type, ID, SOOpenSim[38,42,43,54,55,56,57,58]
Ignasiak et al. [33]Pelvis, sacrum T1-L5, ribs, sternum head–neck6 rot. DOFs (IV)1 rot. DOFs (CV/CT)3 rot. DOFs (CS I)6 rot. DOFs (CS II-X)GenericCS, CT, CV, IV joint (lin.)Act., ID, FSK [66], SOAnyBody[39,67,68]
Lerchl et al. [59]Pelvis, sacrum, L1-L5, thorax, upper limbs, head–neck3 rot. DOFs (IV)Semi-indiv.Lig. (nonlin.)IVD (nonlin.)Actuators, ID, SOSimpack-
Favier et al. [69]Lower limbspelvis, sacrum, L1-L5, thorax (3 segments), upper limbs, head–neck3 rot. DOFs (IV)Semi-indiv.Joint (lin.)Hill type, IK, ID, SOOpenSim-
Malakoutian et al. [60]Pelvis, sacrum, L1-L5, thorax, humeri6 DOFs (IV)genericJoint, IAPHill type, FD-assisted SOAriSynth[70]
Rupp et al. [61]Pelvis, sacrum, L1-L5, thorax6 DOFs (IV)GenericLig. (nonlin.)IVD (nonlin.)Hill type, FDIn-house-
Fasser et al. [62]Pelvis, sacrum, L1-L5, thorax3 rot. DOFs (IV)Semi-indiv.-Hill type, IK, ID, SOMatlab[71]
Bayoglu et al. [72]Pelvis, sacrum, C1-L5, ribs, sternum, skull (3 segments), shoulder (3 Segments)3 rot. DOFs (IV)6 DOFs (CS)1 DOF (CV/CT)Individ.Joint (lin.)Act., ID, SOAnyBody[73,74,75]
Huynh et al. [63]Full-body, C1-L53 rot. DOFs (IV)GenericLig. (lin.)IVD (lin.), IAPIK, ID, SOLifeMOD[76]
Khurelbaatar et al. [64]Pelvis, sacrum, C1-L5, ribs, sternum, upper limbs, head6 DOFs (IV/CS), 3 rot. DOFs (CV)Semi-indiv. (bones)Lig. (nonlin.), IVD (nonlin.), CS cartilage (lin.), facet jointsAct., ID, SORECURDYN-
Guo et al. [77]Pelvis, sacrum, C1-L5, ribs, sternum, upper limbs, head6 DOFs (IV)GenericLig. (nonlin.), IVD (lin.), facet joints, IAPHill type, ALE, FDOpenSim-
The definition of the abbreviations can be found at the end of this article.
Table
Table 2. Overview of representative studies using available original models to address methodological or biomechanical research questions.
Table 2. Overview of representative studies using available original models to address methodological or biomechanical research questions.
StudyFocusModificationsOriginal Model
Actis et al. [48]Methodological Validation for flexion, extension, lateral bending, axial rotation for participants with and without transtibial amputationmodel extension by lower body [110], muscle strength [32], and body mass distribution [86] inclusion of experimental protocol for EMG and kinematic data collection[31]
Arshad et al. [38]Biomechanical Influence of spinal rhythm and IAP on lumbar loads during trunk inclinationAdapted spinal rhythm, inclusion of ligaments, IVD, and IAP[30]
Arx et al. [83]Biomechanical Lumbar loading during different lifting stylesIntegration of measured kinematic data[32]
Banks et al. [58]Biomechanical Comparison of static and dynamic vertebral loading during lifting patient-specific models in an older study populationCT-based individualization and integration of patient-specific kinematic data[32]
Bassani et al. [45]Methodological Model validation for various loading tasks via spinopelvic rhythm and IDP according to [4]Integration of kinematic data[30]
Bassani et al. [47]Biomechanical Effect of spinopelvic sagittal alignment on lumbar loadsVariation of spinal alignment based on four parameters[30].
Bayoglu et al. [75]Methodological Sensitivity of muscle and IV disc force computations to variations in muscle attachment sitesVariation of the location of muscle insertion[72]
Raabe et al. [40]Biomechanical Jogging biomechanicsCombination with full-body model by [111][31]
Beaucage-Gauvreau et al. [49,50,51]Biomechanical Effects of lifting techniques on lumbar loadsAdjust all spinal joints with 3 DOFs and inclusion of kinematic data from motion capturing during lifting[31,40]
Burkhart et al. [54]Methodological Reliability of optoelectronic motion capturing for subject-specific spine model generationCombination with model of lower limbs [110][32]
Malakoutian et al. [70]Methodological Effect of muscle parameters on spinal loadingVariation of biomechanical parameters of paraspinal muscles[60]
Senteler et al. [41]Methodological Joint reaction forces for flexion and liftingCombination with models of upper limbs and neck, IV joints set to 6 DOFs, added passive lin. joint stiffness[31]
Meng et al. [37]Methodological Force-motion coupling in 6-DOF joint6 DOFs (IV), added 6-DOF stiffness[31]
Molinaro et al. [52]Biomechanical Effects of throwing technique solid waste collection occupation on lumbar loadsIncorporation of collected kinematics and EMG data, EMG-assisted muscle force estimation and SO[49]
Schmid et al. [56]Methodological Validation of a thoracolumbar model for children and adolescentsCombination with model of the lower limbs [112], scaling to anthropometry of children and adolescents[32]
Schmid et al. [57]Methodological Feasibility of a skin-marker based method for spinal alignment modelingReduction of muscle architecture, implementation of skin-marker derived alignment[56]
Wang et al. [84]Methodological Implementation of a physiological FSUAdaption of IV joints to represent passive properties of a physiological FSU[32]
Overbergh et al. [55]Methodological Workflow for generation of an image-based (CT), subject-specific thoracolumbar model of spinal deformityAddition of kinematic coupling constraints, personalization of bone geometries, alignment, IV joint definitions and kinematics[32]
Han et al. [36]Methodological Effect of centers of rotation on spinal loads and muscle forces in total disc replacement of lumbar spineLigaments and facet joints added, altering location of CoR[30]
Zhu et al. [46]Biomechanical Effects of lifting techniques on lumbar loadsCombining with models of upper and lower limbs, 6-DOF IV joint, integration of a customized marker set[31]
Kuai et al. [44]Biomechanical Influence of disc herniation on kinematics of the spine and lower limbsIntegration of kinematic data from patients with lumbar disc herniation[30]
Senteler et al. [113]Methodological Sensitivity of intervertebral joint forces to CoR locationAltering location of CoR[41]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lerchl, T.; Nispel, K.; Baum, T.; Bodden, J.; Senner, V.; Kirschke, J.S. Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges. Bioengineering 2023, 10, 202. https://doi.org/10.3390/bioengineering10020202

AMA Style

Lerchl T, Nispel K, Baum T, Bodden J, Senner V, Kirschke JS. Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges. Bioengineering. 2023; 10(2):202. https://doi.org/10.3390/bioengineering10020202

Chicago/Turabian Style

Lerchl, Tanja, Kati Nispel, Thomas Baum, Jannis Bodden, Veit Senner, and Jan S. Kirschke. 2023. "Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges" Bioengineering 10, no. 2: 202. https://doi.org/10.3390/bioengineering10020202

APA Style

Lerchl, T., Nispel, K., Baum, T., Bodden, J., Senner, V., & Kirschke, J. S. (2023). Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges. Bioengineering, 10(2), 202. https://doi.org/10.3390/bioengineering10020202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop