Quantifying the centre of rotation pattern in a multi-body model of the lumbar spine

Identification of a lumped-parameter model of the intervertebral joint from experimental data

Through predictive simulations, multibody models can aid the treatment of spinal pathologies by identifying optimal surgical procedures. Critical to achieving accurate predictions is the definition of the intervertebral joint. The joint pose is often defined by virtual palpation. Intervertebral joint stiffnesses are either derived from literature, or specimen-specific stiffnesses are calculated with optimisation methods. This study tested the feasibility of an optimisation method for determining the specimen-specific stiffnesses and investigated the influence of the assigned joint pose on the subject-specific estimated stiffness. Furthermore, the influence of the joint pose and the stiffness on the accuracy of the predicted motion was investigated. A computed tomography based model of a lumbar spine segment was created. Joints were defined from virtually palpated landmarks sampled with a Latin Hypercube technique from a possible Cartesian space. An optimisation method was used to determine specimen-specific stiffnesses for 500 models. A two-factor analysis was performed by running forward dynamic simulations for ten different stiffnesses for each successfully optimised model. The optimisations calculated a large range of stiffnesses, indicating the optimised specimen-specific stiffnesses were highly sensitive to the assigned joint pose and related uncertainties. A limited number of combinations of optimised joint stiffnesses and joint poses could accurately predict the kinematics. The two-factor analysis indicated that, for the ranges explored, the joint pose definition was more important than the stiffness. To obtain kinematic prediction errors below 1 mm and 1° and suitable specimen-specific stiffnesses the precision of virtually palpated landmarks for joint definition should be better than 2.9 mm.

Ultrasound-guided disc pain induction test for diagnosis of discogenic lumbar pain: a cross-sectional study

BACKGROUND

Several methods can be used to diagnose discogenic pain, but only discoblock can diagnose discogenic pain definitively. This study aimed to examine the usefulness of an ultrasound-guided disc pain induction test for a simple and accurate diagnosis of the culprit lesion.

METHODS

We included 41 patients with lumbar pain in whom pain was induced by an ultrasound-guided disc pain induction test. All patients had confirmed pain at L1/2 to L5/S1 based on an ultrasound-guided disc pain induction test and underwent X-ray photography and magnetic resonance imaging. Seventeen patients who required injection due to severe pain underwent discoblock procedures for discs with the most intense pain, and visual analogue scale (VAS) scores were obtained before and after the procedure for these patients. We analysed the association between painful discs and radiological findings.

RESULTS

Pain induction was noted in a total of 65 discs, and the pain was induced in 23 patients in only one disc. All patients had disc degeneration of Pfirrmann classification grade 1 or higher, with more significant disc degeneration in painful discs than in painless discs. There was no significant relationship between the presence or absence of pain and Modic type. The average VAS measurements improved significantly from 9.5 (pre-procedure) to 2.5 (post-procedure). These results suggest that the most painful discs were the causes of discogenic lumbar pain.

CONCLUSIONS

Our ultrasound-guided disc pain induction test may help diagnose disc degeneration and identify culprit lesions, even when multiple discs exhibit findings of degeneration.

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

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

Computational lumbar spine models: A literature review

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

Muscle-driven and torque-driven centrodes during modeled flexion of individual lumbar spines are disparate

Lumbar spine biomechanics during the forward-bending of the upper body (flexion) are well investigated by both in vivo and in vitro experiments. In both cases, the experimentally observed relative motion of vertebral bodies can be used to calculate the instantaneous center of rotation (ICR). The timely evolution of the ICR, the centrode, is widely utilized for validating computer models and is thought to serve as a criterion for distinguishing healthy and degenerative motion patterns. While in vivo motion can be induced by physiological active structures (muscles), in vitro spinal segments have to be driven by external torque-applying equipment such as spine testers. It is implicitly assumed that muscle-driven and torque-driven centrodes are similar. Here, however, we show that centrodes qualitatively depend on the impetus. Distinction is achieved by introducing confidence regions (ellipses) that comprise centrodes of seven individual multi-body simulation models, performing flexion with and without preload. Muscle-driven centrodes were generally directed superior–anterior and tail-shaped, while torque-driven centrodes were located in a comparably narrow region close to the center of mass of the caudal vertebrae. We thus argue that centrodes resulting from different experimental conditions ought to be compared with caution. Finally, the applicability of our method regarding the analysis of clinical syndromes and the assessment of surgical methods is discussed.

ICR in human cadaveric specimens: An essential parameter to consider in a new lumbar disc prosthesis design

STUDY DESIGN

Biomechanical study in cadaveric specimens.

BACKGROUND

The commercially available lumbar disc prostheses do not reproduce the intact disc's Instantaneous centre of Rotation (ICR), thus inducing an overload on adjacent anatomical structures, promoting secondary degeneration.

AIM

To examine biomechanical testing of cadaveric lumbar spine specimens in order to evaluate and define the ICR of intact lumbar discs.

MATERIAL AND METHODS

Twelve cold preserved fresh human cadaveric lumbosacral spine specimens were subjected to computerized tomography (CT), magnetic resonance imaging (MRI) and biomechanical testing. Kinematic studies were performed to analyse range of movements in order to determine ICR.

RESULTS

Flexoextension and lateral bending tests showed a positive linear correlation between the angle rotated and the displacement of the ICR in different axes.

DISCUSSION

ICR has not been taken into account in any of the available literature regarding lumbar disc prosthesis. Considering our results, neither the actual ball-and-socket nor the withdrawn elastomeric nucleus models fit the biomechanics of the lumbar spine, which could at least in part explain the failure rates of the implants in terms of postoperative failed back syndrome (low back pain). It is reasonable to consider then that an implant should also adapt the equations of the movement of the intact ICR of the joint to the post-surgical ICR.

CONCLUSIONS

This is the first cadaveric study on the ICR of the human lumbar spine. We have shown that it is feasible to calculate and consider this parameter in order to design future prosthesis with improved clinical and biomechanical characteristics.