Linking continuous and discrete intervertebral disc models through homogenisation

A new method to design energy-conserving surrogate models for the coupled, nonlinear responses of intervertebral discs

The aim of this study was to design physics-preserving and precise surrogate models of the nonlinear elastic behaviour of an intervertebral disc (IVD). Based on artificial force-displacement data sets from detailed finite element (FE) disc models, we used greedy kernel and polynomial approximations of second, third and fourth order to train surrogate models for the scalar force-torque -potential. Doing so, the resulting models of the elastic IVD responses ensured the conservation of mechanical energy through their structure. At the same time, they were capable of predicting disc forces in a physiological range of motion and for the coupling of all six degrees of freedom of an intervertebral joint. The performance of all surrogate models for a subject-specific L4 | 5 disc geometry was evaluated both on training and test data obtained from uncoupled (one-dimensional), weakly coupled (two-dimensional), and random movement trajectories in the entire six-dimensional (6d) physiological displacement range, as well as on synthetic kinematic data. We observed highest precisions for the kernel surrogate followed by the fourth-order polynomial model. Both clearly outperformed the second-order polynomial model which is equivalent to the commonly used stiffness matrix in neuro-musculoskeletal simulations. Hence, the proposed model architectures have the potential to improve the accuracy and, therewith, validity of load predictions in neuro-musculoskeletal spine models.

Recent Advances in Coupled MBS and FEM Models of the Spine—A Review

How back pain is related to intervertebral disc degeneration, spinal loading or sports-related overuse remains an unanswered question of biomechanics. Coupled MBS and FEM simulations can provide a holistic view of the spine by considering both the overall kinematics and kinetics of the spine and the inner stress distribution of flexible components. We reviewed studies that included MBS and FEM co-simulations of the spine. Thereby, we classified the studies into unidirectional and bidirectional co-simulation, according to their data exchange methods. Several studies have demonstrated that using unidirectional co-simulation models provides useful insights into spinal biomechanics, although synchronizing the two distinct models remains a key challenge, often requiring extensive manual intervention. The use of a bidirectional co-simulation features an iterative, automated process with a constant data exchange between integrated subsystems. It reduces manual corrections of vertebra positions or reaction forces and enables detailed modeling of dynamic load cases. Bidirectional co-simulations are thus a promising new research approach for improved spine modeling, as a main challenge in spinal biomechanics is the nonlinear deformation of the intervertebral discs. Future studies will likely include the automated implementation of patient-specific bidirectional co-simulation models using hyper- or poroelastic intervertebral disc FEM models and muscle forces examined by an optimization algorithm in MBS. Applications range from clinical diagnosis to biomechanical analysis of overload situations in sports and injury prediction.

Effects of geometric individualisation of a human spine model on load sharing: neuro-musculoskeletal simulation reveals significant differences in ligament and muscle contribution

In spine research, two possibilities to generate models exist: generic (population-based) models representing the average human and subject-specific representations of individuals. Despite the increasing interest in subject specificity, individualisation of spine models remains challenging. Neuro-musculoskeletal (NMS) models enable the analysis and prediction of dynamic motions by incorporating active muscles attaching to bones that are connected using articulating joints under the assumption of rigid body dynamics. In this study, we used forward-dynamic simulations to compare a generic NMS multibody model of the thoracolumbar spine including fully articulated vertebrae, detailed musculature, passive ligaments and linear intervertebral disc (IVD) models with an individualised model to assess the contribution of individual biological structures. Individualisation was achieved by integrating skeletal geometry from computed tomography and custom-selected muscle and ligament paths. Both models underwent a gravitational settling process and a forward flexion-to-extension movement. The model-specific load distribution in an equilibrated upright position and local stiffness in the L4/5 functional spinal unit (FSU) is compared. Load sharing between occurring internal forces generated by individual biological structures and their contribution to the FSU stiffness was computed. The main finding of our simulations is an apparent shift in load sharing with individualisation from an equally distributed element contribution of IVD, ligaments and muscles in the generic spine model to a predominant muscle contribution in the individualised model depending on the analysed spine level.

Intuitive assessment of modeled lumbar spinal motion by clustering and visualization of finite helical axes

A variety of medical imaging procedures, cadaver experiments, and computer models have been utilized to capture, depict, and understand the motion of the human lumbar spine. Particular interest lies in assessing the relative movement between two adjacent vertebrae, which can be represented by a temporal evolution of finite helical axes (FHA). Mathematically, this FHA evolution constitutes a seven-dimensional quantity: one dimension for the time, two for the (normalized) direction vector, another two for the (unique) position vector, as well as one for each the angle of rotation around and the amount of translation along the axis. Predominantly in the literature, however, movements are assumed to take place in certain physiological planes on which FHA are projected. The resulting three-dimensional quantity - the so-called centrode - is easily presentable but leaves out substantial pieces of available data. Here, we investigate and assess several possibilities to visualize subsets of FHA data of increasing dimensionality. Finally, we utilize an agglomerative hierarchical clustering algorithm and propose a novel visualization technique, namely the quiver principal axis plot (QPAP), to depict the entirety of information inherent to hundreds or thousands of FHA. The QPAP method is applied to flexion-extension, lateral bending, and axial rotation movements of a lumbar spine within both a reduced model as well as a complex upper body system.

Embodiment of intra-abdominal pressure in a flexible multibody model of the trunk and the spinal unloading effects during static lifting tasks

The role of intra-abdominal pressure (IAP) in spinal load reduction has remained controversial, partly because previous musculoskeletal models did not introduce the pressure generating mechanism. In this study, an integrated computational methodology is proposed to combine the IAP change with core muscle activations. An ideal gas relationship was introduced to calculate pressure distribution within the abdominal cavity. Additionally, based on flexible multibody dynamics, a muscle membrane element was developed by incorporating the muscular fiber deformation, inter-fiber stiffness, and volume constancy. This element was then utilized in discretizing the diaphragm and transversus abdominis, forming an IAP-muscle coupling system of the abdominal cavity. Based on this methodology, a forward dynamic simulation of spinal flexion was presented to examine the unloading effect of abdominal breathing. The results confirm that core muscle contraction during the abdominal breathing cycle can substantially reduce the forces of spinal compression together with trunk extensor muscles, and this effect is more pronounced when the IAP increase is produced by contraction of the transversus abdominis. This unloading effect still holds even with the co-activation of other abdominal muscles, providing a potential choice when designing trunk movements during weight-lifting tasks.

Estimation of Trunk Muscle Forces Using a Bio-Inspired Control Strategy Implemented in a Neuro-Osteo-Ligamentous Finite Element Model of the Lumbar Spine

Low back pain (LBP), the leading cause of disability worldwide, remains one of the most common and challenging problems in occupational musculoskeletal disorders. The effective assessment of LBP injury risk, and the design of appropriate treatment modalities and rehabilitation protocols, require accurate estimation of the mechanical spinal loads during different activities. This study aimed to: (1) develop a novel 2D beam-column finite element control-based model of the lumbar spine and compare its predictions for muscle forces and spinal loads to those resulting from a geometrically matched equilibrium-based model; (2) test, using the foregoing control-based finite element model, the validity of the follower load (FL) concept suggested in the geometrically matched model; and (3) investigate the effect of change in the magnitude of the external load on trunk muscle activation patterns. A simple 2D continuous beam-column model of the human lumbar spine, incorporating five pairs of Hill's muscle models, was developed in the frontal plane. Bio-inspired fuzzy neuro-controllers were used to maintain a laterally bent posture under five different external loading conditions. Muscle forces were assigned based on minimizing the kinematic error between target and actual postures, while imposing a penalty on muscular activation levels. As compared to the geometrically matched model, our control-based model predicted similar patterns for muscle forces, but at considerably lower values. Moreover, irrespective of the external loading conditions, a near (<3°) optimal FL on the spine was generated by the control-based predicted muscle forces. The variation of the muscle forces with the magnitude of the external load within the simulated range at the L1 level was found linear. This work presents a novel methodology, based on a bio-inspired control strategy, that can be used to estimate trunk muscle forces for various clinical and occupational applications toward shedding light on the ever-elusive LBP etiology.

The response surface method-genetic algorithm for identification of the lumbar intervertebral disc material parameters

Long-term compressive load on the lumbar intervertebral disc (IVD) might lead to lumbar IVD herniation. Exploring the material parameters of normal and degenerative enucleated IVDs is the basis for studying their mechanical behavior. According to the inverse analysis principle of the parameter estimation, an optimization method was proposed to identify the parameters of the porous material of the lumbar IVD based on finite element inverse analysis. The poroelastic finite element models were established in line with the compression creep experiment. The material parameters were combined by Box-Behnken design (BBD), and the response surface (RS) models were constructed using a quadratic polynomial with cross terms and optimized by genetic algorithm (GA). The results showed that the simulation result of the best material parameter combination had a good agreement with the experiment. Compared with the normal lumbar IVD, the elastic modulus and permeability decreased, and Poisson's ratio increased for the enucleated disk, resulting in a significant difference in mechanical properties. The algorithm used in this study can reduce the parameter identification error compared with only the RS method, and decrease the number of finite element simulations compared with only the GA.

Lumbar spinal ligament characteristics extracted from stepwise reduction experiments allow for preciser modeling than literature data

Lumbar ligaments play a key role in stabilizing the spine, particularly assisting muscles at wide-range movements. Hence, valid ligament force–strain data are required to generate physiological model predictions. These data have been obtained by experiments on single ligaments or functional units throughout the literature. However, contrary to detailed spine geometries, gained, for instance, from CT data, ligament characteristics are often inattentively transferred to multi-body system (MBS) or finite element models. In this paper, we use an elaborated MBS model of the lumbar spine to demonstrate how individualized ligament characteristics can be obtained by reversely reenacting stepwise reduction experiments, where the range of motion (ROM) was measured. We additionally validated the extracted characteristics with physiological experiments on intradiscal pressure (IDP). Our results on a total of in each case 160 ROM and 49 IDP simulations indicated superiority of our procedure (seven and eight outliers) toward the incorporation of classical literature data (on average 71 and 31 outliers).

Musculoskeletal modelling of the human cervical spine for the investigation of injury mechanisms during axial impacts

Head collisions in sport can result in catastrophic injuries to the cervical spine. Musculoskeletal modelling can help analyse the relationship between motion, external forces and internal loads that lead to injury. However, impact specific musculoskeletal models are lacking as current viscoelastic values used to describe cervical spine joint dynamics have been obtained from unrepresentative quasi-static or static experiments. The aim of this study was to develop and validate a cervical spine musculoskeletal model for use in axial impacts. Cervical spine specimens (C2-C6) were tested under measured sub-catastrophic loads and the resulting 3D motion of the vertebrae was measured. Specimen specific musculoskeletal models were then created and used to estimate the axial and shear viscoelastic (stiffness and damping) properties of the joints through an optimisation algorithm that minimised tracking errors between measured and simulated kinematics. A five-fold cross validation and a Monte Carlo sensitivity analysis were conducted to assess the performance of the newly estimated parameters. The impact-specific parameters were integrated in a population specific musculoskeletal model and used to assess cervical spine loads measured from Rugby union impacts compared to available models. Results of the optimisation showed a larger increase of axial joint stiffness compared to axial damping and shear viscoelastic parameters for all models. The sensitivity analysis revealed that lower values of axial stiffness and shear damping reduced the models performance considerably compared to other degrees of freedom. The impact-specific parameters integrated in the population specific model estimated more appropriate joint displacements for axial head impacts compared to available models and are therefore more suited for injury mechanism analysis.

A multi-body model for comparative study of cervical traction simulation - development, improvement and validation

A computer simulation model was developed to study the dynamic behavior of the cervical spine during cervical traction therapy in inclined and sitting traction positions. The model improved upon an old model with additional components to represent the behavior of the intervertebral discs and the posterior ligaments. The simulation result of the new model was compared against the cervical traction data from a radiographic experiment in both positions. The simulation results of the old model and new model were compared to illustrate the improvement. Using the new model, we compared the timing response of cervical traction in the inclined and sitting positions.