Image driven subject-specific finite element models of spinal biomechanics

Differences in lumbar spine intradiscal pressure between standing and sitting postures: a comprehensive literature review

Background

Musculoskeletal disorders (MSDs), especially in the lumbar spine, are a leading concern in occupational health. Work activities associated with excessive exposure are a source of risk for MSDs. The optimal design of workplaces requires changes in both sitting and standing postures. In order to secure such a design scientifically proved quantitative data are needed that would allow for the assessment of differences in spine load due to body posture and/or exerted force. Intradiscal pressure (IP) measurement in the lumbar spine is the most direct method of estimating spinal loads. Hence, this study aims at the quantitative evaluation of differences in lumbar spine load due to body posture and exerted forces, based on IP reported in publications obtained from a comprehensive review of the available literature.

Methodology

In order to collect data from studies measuring IP in the lumbar spine, three databases were searched. Studies with IP for living adults, measured in various sitting and standing postures, where one of these was standing upright, were included in the analysis. For data to be comparable between studies, the IP for each position was referenced to upright standing. Where different studies presented IP for the same postures, those relative IPs (rIP) were merged. Then, an analysis of the respective outcomes was conducted to find the possible relationship of IPs dependent on a specific posture.

Results

A preliminary analysis of the reviewed papers returned nine items fulfilling the inclusion and exclusion criteria. After merging relative IPs from different studies, rIP for 27 sitting and 26 standing postures was yielded. Some of the data were useful for deriving mathematical equations expressing rIP as a function of back flexion angle and exerted force in the form of a second degree polynomial equation for the standing and sitting positions. The equations showed that for the standing posture, the increase in IP with increasing back flexion angle is steeper when applying an external force than when maintaining body position only. In a sitting position with the back flexed at 20°, adding 10 kg to each hand increases the IP by about 50%. According to the equations developed, for back flexion angles less than 20°, the IP is greater in sitting than in standing. When the angle is greater than 20°, the IP in the sitting position is less than in the standing position at the same angle of back flexion.

Conclusions

Analysis of the data from the reviewed papers showed that: sitting without support increases IP by about 30% in relation to upright standing; a polynomial of the second degree defines changes in IP as a function of back flexion for for both postures. There are differences in the pattern of changes in IP with a back flexion angle between sitting and standing postures, as back flexion in standing increases IP more than in sitting.

Biomechanical response of lumbar intervertebral disc in daily sitting postures: a poroelastic finite element analysis

This study aims to establish and validate a poroelastic L4-L5 finite element model to evaluate the effect of different sitting postures and their durations on the mechanical responses of the disc. During the sustained loading conditions, the height loss, fluid loss and von-Mises stress gradually increased, but the intradiscal pressure decreased. The varying rates of aforementioned parameters were more significant at the initial loading stage and less so at the end. The predicted values in the flexed sitting posture were significantly greater than other postures. The extended sitting posture caused an obvious von-Mises stress concentration in the posterior region of the inter-lamellar matrix. From the biomechanical perspective, prolonged sitting may pose a high risk of lumbar disc degeneration, and therefore adjusting the posture properly in the early stage of sitting time may be useful to mitigate that. Additionally, upright sitting is a safer posture, while flexed sitting posture is more harmful.

Novel force-displacement control passive finite element models of the spine to simulate intact and pathological conditions; comparisons with traditional passive and detailed musculoskeletal models

Passive finite element (FE) models of the spine are commonly used to simulate intact and various pre- and postoperative pathological conditions. Being devoid of muscles, these traditional models are driven by simplistic loading scenarios, e.g., a constant moment and compressive follower load (FL) that do not properly mimic the complex in vivo loading condition under muscle exertions. We aim to develop novel passive FE models that are driven by more realistic yet simple loading scenarios, i.e., in vivo vertebral rotations and pathological-condition dependent FLs (estimated based on detailed musculoskeletal finite element (MS-FE) models). In these novel force-displacement control FE models, unlike the traditional passive FE models, FLs vary not only at different spine segments (T12-S1) but between intact, pre- and postoperative conditions. Intact, preoperative degenerated, and postoperative fused conditions at the L4-L5 segment for five static in vivo activities in upright and flexed postures were simulated by the traditional passive FE, novel force-displacement control FE, and gold-standard detailed MS-FE spine models. Our findings indicate that, when compared to the MS-FE models, the force-displacement control passive FE models could accurately predict the magnitude of disc compression force, intradiscal pressure, annulus maximal von Mises stress, and vector sum of all ligament forces at adjacent segments (L3-L4 and L5-S1) but failed to predict disc shear and facet joint forces. In this regard, the force-displacement control passive FE models were much more accurate than the traditional passive FE models. Clinical recommendations made based on traditional passive FE models should, therefore, be interpreted with caution.

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.

Quantitative MRI to Characterize the Nucleus Pulposus Morphological and Biomechanical Variation According to Sagittal Bending Load and Radial Fissure, an ex vivo Ovine Specimen Proof-of-Concept Study

Background and context: Low back pain is a dramatic burden worldwide. Discography studies have shown that 39% of chronic low back pain patients suffer from discogenic pain due to a radial fissure of intervertebral disc. This can have major implications in clinical therapeutic choices. The use of discography is restricted because of its invasiveness and interest in it remains low as it represents a static condition of the disc morphology. Magnetic Resonance Imaging (MRI) appears to be less invasive but does not describe the biomechanical dynamic behavior of the fissure. Purpose: We aimed to seek a quantitative MRI protocol combined with ex vivo sagittal loading to analyze the morphological and biomechanical changes of the intervertebral disc structure and stress distribution. Study design: Proof of concept. Methods: We designed a proof-of-concept ovine study including 3 different 3.0 T-MRI sequences (T₂-weighted, T₁ and T₂ mapping). We analyzed 3 different mechanical states (neutral, flexion and extension) on a fresh ovine spine specimen to characterize an intervertebral disc before and after puncturing the anterior part of the annulus fibrosus. We used a mark tracking method to calculate the bending angles and the axial displacements of the discal structures. In parallel, we created a finite element model to calculate the variation of the axial stress and the maximal intensity shear stress, extrapolated from our experimental boundary conditions. Results: Thanks to an original combination of specific nuclear relaxation time quantifications (T₁, T₂) of the discal tissue, we characterized the nucleus movement/deformation into the fissure according to the synchronous mechanical load. This revealed a link between disc abnormality and spine segment range of motion capability. Our finite element model highlighted significant variations within the stress distribution between intact and damaged disc. Conclusion: Quantitative MRI appears to provide a new opportunity to characterize intra-discal structural morphology, lesions and stress changes under the influence of mechanical load. This preliminary work could have substantial implications for non-invasive disc exploration and could help to validate novel therapies for disc treatment.

Spine biomechanical testing methodologies: The controversy of consensus vs scientific evidence

Biomechanical testing methodologies for the spine have developed over the past 50 years. During that time, there have been several paradigm shifts with respect to techniques. These techniques evolved by incorporating state-of-the-art engineering principles, in vivo measurements, anatomical structure-function relationships, and the scientific method. Multiple parametric studies have focused on the effects that the experimental technique has on outcomes. As a result, testing methodologies have evolved, but there are no standard testing protocols, which makes the comparison of findings between experiments difficult and conclusions about in vivo performance challenging. In 2019, the international spine research community was surveyed to determine the consensus on spine biomechanical testing and if the consensus opinion was consistent with the scientific evidence. More than 80 responses to the survey were received. The findings of this survey confirmed that while some methods have been commonly adopted, not all are consistent with the scientific evidence. This review summarizes the scientific literature, the current consensus, and the authors' recommendations on best practices based on the compendium of available evidence.

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

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

Quantification of lumbar spine load transfer is important for understanding low back pain, especially among persons with a lower limb amputation. Computational modeling provides a helpful solution for obtaining estimates of in vivo loads. A multiscale model was constructed by combining musculoskeletal and finite element (FE) models of the lumbar spine to determine tissue loading during daily activities. Three-dimensional kinematic and ground reaction force data were collected from participants with ([Formula: see text]) and without ([Formula: see text]) a unilateral transtibial amputation (TTA) during 5 sit-to-stand trials. We estimated tissue-level load transfer from the multiscale model by controlling the FE model with intervertebral kinematics and muscle forces predicted by the musculoskeletal model. Annulus fibrosis stress, intradiscal pressure (IDP), and facet contact forces were calculated using the FE model. Differences in whole-body kinematics, muscle forces, and tissue-level loads were found between participant groups. Notably, participants with TTA had greater axial rotation toward their intact limb ([Formula: see text]), greater abdominal muscle activity ([Formula: see text]), and greater overall tissue loading throughout sit-to-stand ([Formula: see text]) compared to able-bodied participants. Both normalized (to upright standing) and absolute estimates of L4-L5 IDP were close to in vivo values reported in the literature. The multiscale model can be used to estimate the distribution of loads within different lumbar spine tissue structures and can be adapted for use with different activities, populations, and spinal geometries.

A Novel Modelling Methodology Which Predicts the Structural Behaviour of Vertebral Bodies under Axial Impact Loading: A Finite Element and DIC Study

Cervical spine injuries (CSIs) arising from collisions are uncommon in contact sports, such as rugby union, but their consequences can be devastating. Several FE modelling approaches are available in the literature, but a fully calibrated and validated FE modelling framework for cervical spines under compressive dynamic-impact loading is still lacking and material properties are not adequately calibrated for such events. This study aimed to develop and validate a methodology for specimen-specific FE modelling of vertebral bodies under impact loading. Thirty-five (n = 35) individual vertebral bodies (VBs) were dissected from porcine spine segments, potted in bone cement and μCT scanned. A speckle pattern was applied to the anterior faces of the bones to allow digital image correlation (DIC), which monitored the surface displacements. Twenty-seven (n = 27) VBs were quasi-statically compressively tested to a load up to 10 kN from the cranial side. Specimen-specific FE models were developed for fourteen (n = 14) of the samples in this group. The material properties were optimised based on the experimental load-displacement data using a specimen-specific factor (kGSstatic) to calibrate a density to Young’s modulus relationship. The average calibration factor arising from this group was calculated (K¯GSstatic) and applied to a control group of thirteen (n = 13) samples. The resulting VB stiffnesses was compared to experimental findings. The final eight (n = 8) VBs were subjected to an impact load applied via a falling mass of 7.4kg at a velocity of 3.1ms⁻¹. Surface displacements and strains were acquired from the anterior VB surface via DIC, and the impact load was monitored with two load cells. Specimen-specific FE models were created for this dynamic group and material properties were assigned again based on the density–Young’s modulus relationship previously validated for static experiments, supplemented with an additional factor (KGSdynamic). The optimised conversion factor for quasi-static loading, K¯GSstatic, had an average of 0.033. Using this factor, the validation models presented an average numerical stiffness value 3.72% greater than the experimental one. From the dynamic loading experiments, the value for KGSdynamic was found to be 0.14, 4.2 times greater than K¯GSstatic. The average numerical stiffness was 2.3% greater than in the experiments. Almost all models presented similar stiffness variations and regions of maximum displacement to those observed via DIC. The developed FE modelling methodology allowed the creation of models which predicted both static and dynamic behaviour of VBs. Deformation patterns on the VB surfaces were acquired from the FE models and compared to DIC data, achieving high agreement. This methodology is now validated to be fully applied to create whole cervical spine models to simulate axial impact scenarios replicating rugby collision events.

Modelling tri-cortical pedicle screw fixation in thoracic vertebrae under osteoporotic condition: A finite element analysis based on computed tomography

BACKGROUND AND OBJECTIVE

The technique of tri-cortical pedicle screw (TCPS) has been used to improve the anchoring strength in the sacral vertebrae. However, no studies have reported their application in the thoracic vertebrae. Our research is aimed to assess the stability and strength of the TCPS in thoracic vertebrae under osteoporotic condition by three dimensional (3D) finite element method on the basis of medical image reconstruction using computed tomography (CT), and verifying its effectiveness in clinical application.

MATERIALS AND METHODS

The 3D finite element models were constructed using Mimcs to transfer two dimensional CT images into 3D models by marching cubes algorithm of six-partition. Six physiological activities were simulated in 3D finite element models. Compared with the strength and stability of the uni-cortical pedicle screw (UCPS) and bi-cortical pedicle screw (BCPS), the effectiveness of TCPS was assessed. The stress distribution and maximum stress were measured to evaluate the strength. The maximum displacement and the range of motion were analysed to assessed the stability.

EXPERIMENTAL RESULTS

Four geometrically accurate and nonlinear T7-T9 finite element models were constructed successfully by 3D finite element method based on the CT images. Three kinds of internal fixation methods in the osteoporotic thoracic vertebral body can improved the maximum stress, decrease the maximum displacement and range of motion in six physiological activities. The range of motion and maximum displacement of TCPS decreased more significantly than that of UCPS and BCPS. The maximum von Mises stress of TCPS was minimum and UCPS was maximum under the condition of extension, right lateral bending, left rotation and right rotation.

CONCLUSIONS

Effectively, TCPS can provide better stability and strength than that of UCPS and BCPS techniques in the osteoporotic thoracic vertebrae. In practice, the technique of TCPS can be applied in the osteoporotic thoracic vertebral body to enhance the griping strength of the screws and reduce the risk of pedicle screw loosening. However, further cadaver experiments and more biomechanical analysis are necessary to confirmed our findings.

Subject-specific loads on the lumbar spine in detailed finite element models scaled geometrically and kinematic-driven by radiography images

Traditional load-control musculoskeletal and finite element (FE) models of the spine fail to accurately predict in vivo intervertebral joint loads due mainly to the simplifications and assumptions when estimating redundant trunk muscle forces. An alternative powerful protocol that bypasses the calculation of muscle forces is to drive the detailed FE models by image-based in vivo displacements. Development of subject-specific models, however, both involves the risk of extensive radiation exposures while imaging in supine and upright postures and is time consuming in terms of the reconstruction of the vertebrae, discs, ligaments, and facets geometries. This study therefore aimed to introduce a remedy for the development of subject-specific FE models by scaling the geometry of an existing detailed FE model of the T12-S1 lumbar spine. Five subject-specific scaled models were driven by their own radiography image-based displacements in order to predict joint loads, ligament forces, facet joint forces, and disc fiber strains during relaxed upright as well as moderate flexion and extension tasks. The predicted intradiscal pressures were found in adequate agreement with in vivo data for upright, flexion, and extension tasks. There were however large intersubject variations in the estimated joint loads and facet forces.

Presentation of an Approach on Determination of the Natural Frequency of Human Lumbar Spine Using Dynamic Finite Element Analysis

Occurring resonance may negatively affect the health of the human lumbar spine. Hence, vibration generated in working and living environments should be optimized to avoid resonance when identifying the natural frequency of the human lumbar spine. The range of the natural frequency of the human lumbar spine has been investigated, but its specific numerical value has not been determined yet. This study aimed at presenting an approach based on resonance for predicting the specific numerical value of the natural frequency of the human lumbar spine. The changes in the numerical fluctuation amplitudes and the cycles of lumbar mechanical parameters during resonance are greater than those during nonresonant vibration. Given that the range of the natural frequency has been identified, vibrations at different excitation frequencies within this range can be applied in a human lumbar finite element model for dynamic finite element analysis. When the excitation frequency is close to the natural frequency, resonance occurs, causing great changes in the numerical fluctuation amplitudes and the cycles of lumbar mechanical parameters. Therefore, the natural frequency of the lumbar finite element model could be back-calculated. Results showed that the natural frequency of the established model was 3.5 Hz. Meanwhile, the closer the excitation frequency was to the natural frequency, the greater the changes in the numerical fluctuation amplitudes and cycles in the parameters would be. This study presented an approach for predicting the specific numerical value of the natural frequency of the human lumbar spine. Identifying the natural frequency assists in finding preventive measures for lumbar injury caused by vibration and in designing the vibration source in working and living environments to avoid approximating to the natural frequency of the human lumbar spine.

Data-driven Classification of the 3D Spinal Curve in Adolescent Idiopathic Scoliosis with an Applications in Surgical Outcome Prediction

Adolescent idiopathic scoliosis (AIS) is a three-dimensional (3D) deformity of the spinal column. For progressive deformities in AIS, the spinal fusion surgery aims to correct and stabilize the deformity; however, common surgical planning approaches based on the 2D X-rays and subjective surgical decision-making have been challenged by poor clinical outcomes. As the suboptimal surgical outcomes can significantly impact the cost, risk of revision surgery, and long-term rehabilitation of adolescent patients, objective patient-specific models that predict the outcome of different treatment scenarios are in high demand. 3D classification of the spinal curvature and identifying the key surgical parameters influencing the outcomes are required for such models. Here, we show that K-means clustering of the isotropically scaled 3D spinal curves provides an effective, data-driven method for classification of patients. We further propose, and evaluate in 67 right thoracic AIS patients, that by knowing the patients’ pre-operative and early post-operation clusters and the vertebral levels which were instrumented during the surgery, the two-year outcome cluster can be determined. This framework, once applied to a larger heterogeneous patient dataset, can further isolate the key surgeon-modifiable parameters and eventually lead to a patient-specific predictive model based on a limited number of factors determinable prior to surgery.

A model-based approach for estimation of changes in lumbar segmental kinematics associated with alterations in trunk muscle forces

The kinematics information from imaging, if combined with optimization-based biomechanical models, may provide a unique platform for personalized assessment of trunk muscle forces (TMFs). Such a method, however, is feasible only if differences in lumbar spine kinematics due to differences in TMFs can be captured by the current imaging techniques. A finite element model of the spine within an optimization procedure was used to estimate segmental kinematics of lumbar spine associated with five different sets of TMFs. Each set of TMFs was associated with a hypothetical trunk neuromuscular strategy that optimized one aspect of lower back biomechanics. For each set of TMFs, the segmental kinematics of lumbar spine was estimated for a single static trunk flexed posture involving, respectively, 40° and 10° of thoracic and pelvic rotations. Minimum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 0° to 0.7° and 0 mm to 0.04 mm, respectively. Maximum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 2.4° to 7.6° and 0.11 mm to 0.39 mm, respectively. The differences in kinematics of lumbar segments between each combination of two sets of TMFs in 97% of cases for angular deformation and 55% of cases for translational deformation were within the reported accuracy of current imaging techniques. Therefore, it might be possible to use image-based kinematics of lumbar segments along with computational modeling for personalized assessment of TMFs.

The effect of screw tunnels on the biomechanical stability of vertebral body after pedicle screws removal: a finite element analysis

PURPOSE

Posterior reduction and pedicle screw fixation is a widely used procedure for thoracic and lumbar vertebrae fractures. Usually, the pedicle screws would be removed after the fracture healing and screw tunnels would be left. The aim of this study is to evaluate the effect of screw tunnels on the biomechanical stability of the lumbar vertebral body after pedicle screws removal by finite element analysis (FEA).

METHODS

First, the CT values of the screw tunnels wall in the fractured vertebral bodies were measured in patients whose pedicle screws were removed, and they were then compared with the values of vertebral cortical bone. Second, an adult patient was included and the CT images of the lumbar spine were harvested. Three dimensional finite element models of the L1 vertebra with unilateral or bilateral screw tunnels were created based on the CT images. Different compressive loads were vertically acted on the models. The maximum loads which the models sustained and the distribution of the force in the different parts of the models were recorded and compared with each other.

RESULTS

The CT values of the tunnels wall and vertebral cortical bone were 387.126±62.342 and 399.204±53.612, which were not statistically different (P=0.149). The models of three dimensional tetrahedral mesh finite element of normal lumbar 1 vertebra were established with good geometric similarity and realistic appearance. After given the compressive loads, the cortical bone was the first one to reach its ultimate stress. The maximum loads which the bilateral screw tunnels model, unilateral screw tunnel model, and normal vertebral model can sustain were 3.97 Mpa, 3.83 Mpa, and 3.78 Mpa, respectively. For the diameter of the screw tunnels, the model with a diameter of 6.5 mm could sustain the largest load. In addition, the stress distributing on the outside of the cortical bone gradually decreased as the thickness of the tunnel wall increased.

CONCLUSIONS

Based on the FEA, pedicle screw tunnels would not decrease the biomechanical stability and strength of the vertebral body. A large diameter of screw tunnel and thick tunnel wall were helpful for the biomechanical stability of the vertebral body.

Impact of material and morphological parameters on the mechanical response of the lumbar spine - A finite element sensitivity study

Finite element models are frequently used to study lumbar spinal biomechanics. Deterministic models are used to reflect a certain configuration, including the means of geometrical and material properties, while probabilistic models account for the inherent variability in the population. Because model parameters are generally uncertain, their predictive power is frequently questioned. In the present study, we determined the sensitivities of spinal forces and motions to material parameters of intervertebral discs, vertebrae, and ligaments and to lumbar morphology. We performed 1200 model simulations using a generic model of the human lumbar spine loaded under pure moments. Coefficients of determination and of variation were determined for all parameter and response combinations. Material properties of the vertebrae displayed the least impact on results, whereas those of the discs and morphology impacted most. The most affected results were the axial compression forces in the vertebral body and in several ligaments during flexion and the facet-joint forces during extension. Intervertebral rotations were considerably affected only when several parameters were varied simultaneously. Results can be used to decide which model parameters require careful consideration in deterministic models and which parameters might be omitted in probabilistic studies. Findings allow quantitative estimation of a model׳s precision.