An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.     

Biomechanical comparison between fusion of two vertebrae and implantation of an artificial intervertebral disc

Surgical treatments for lower back pain can be distributed into two main groups: fusion (arthrodesis) and disc replacement (arthroplasty). The objective of this study was to compare, under severe loading conditions, the biomechanics of the lumbar spine treated either by fusion or total disc replacement (TDR). A three-dimensional model of a two-level ligamentous lumbar segment was created and simulated through static analyses with the finite-element method (FEM) software ABAQUS. The model was validated by comparing mobility, pressure on the facets, force in the ligaments, maximum stresses, disc bulge, and endplate deflection with measured data given in the literature. The FEM analysis predicted that the mobility of the model after arthrodesis on the upper level was reduced in all rotational degrees of freedom by an average of approximately 44%, relative to healthy normal discs. Conversely, the mobility of the model after TDR on the upper level was increased in all rotational degrees of freedom by an average of approximately 52%. The level implanted with the artificial disc showed excessive ligament tensions (greater than 500 N), high facet pressures (greater than 3 MPa), and a high risk of instability. The mobility and the stresses in the level adjacent to the arthroplasty were also increased. In conclusion, the model for an implanted movable artificial disc illustrated complications common to spinal arthroplasty and showed greater risk of instability and further degeneration than predicted for the fused model. This modeling technique provides an accurate means for assessing potential biomechanical risks and can be used to improve the design of future artificial intervertebral discs.     

High pressures and asymmetrical stresses in the scoliotic disc in the absence of muscle loading

Background

Loads acting on scoliotic spines are thought to be asymmetric and involved in progression of the scoliotic deformity; abnormal loading patterns lead to changes in bone and disc cell activity and hence to vertebral body and disc wedging. At present however there are no direct measurements of intradiscal stresses or pressures in scoliotic spines. The aim of this study was to obtain quantitative measurements of the intradiscal stress environment in scoliotic intervertebral discs and to determine if loads acting across the scoliotic spine are asymmetric. We performed in vivo measurements of stresses across the intervertebral disc in patients with scoliosis, both parallel (termed horizontal) and perpendicular (termed vertical) to the end plate, using a side mounted pressure transducer (stress profilometry)

Methods

Stress profilometry was used to measure horizontal and vertical stresses at 5 mm intervals across 25 intervertebral discs of 7 scoliotic patients during anterior reconstructive surgery. A state of hydrostatic pressure was defined by identical horizontal and vertical stresses for at least two consecutive readings. Results were compared with similar stress profiles measured during surgery across 10 discs of 4 spines with no lateral curvature and with data from the literature.

Results

Profiles across scoliotic discs were very different from those of normal, young, healthy discs of equivalent age previously presented in the literature. Hydrostatic pressure regions were only seen in 14/25 discs, extended only over a short distance. Non-scoliotic discs of equivalent age would be expected to show large centrally placed hydrostatic nuclear regions in all discs. Mean pressures were significantly greater (0.25 MPa) than those measured in other anaesthetised patients (<0.07 MPa). A stress peak was seen in the concave annulus in 13/25 discs. Stresses in the concave annulus were greater than in the convex annulus indicating asymmetric loading in these anaesthetised, recumbent patients.

Conclusion

Intradiscal pressures and stresses in scoliotic discs are abnormal, asymmetrical and high in magnitude even in the absence of significant applied muscle loading. The origin of these abnormal stresses is unclear.     

Biomechanical analysis of lumbar interbody fusion cages with various lordotic angles: a finite element study

Inappropriate lordotic angle of lumbar fusion cage could be associated with cage damage or subsidence. The biomechanical influence of cage lordotic angle on lumbar spine has not been fully investigated. Four surgical finite element models were constructed by inserting cages with various lordotic angles at L3-L4 disc space. The four motion modes were simulated. The range of motion (ROM) decreased with increased lordotic angle of cage in flexion, extension, and rotation, whereas it was not substantially changed in bending. The maximum stress in cage decreased with increased lordotic angle of cage in all motion modes. The maximum stress in endplate at surgical level increased with increased lordotic angle of cage in flexion and rotation, whereas it was not substantially changed in extension and bending. The facet joint force (FJF) was much smaller than that for the intact conditions in extension, bending, and rotation, while it was not substantially changed in flexion. In conclusion, the ROM, stresses in the cage and endplate at surgical level are sensitive to the lordotic angle of cage. The increased cage lordotic angle may provide better stability and reduce the risk of cage damage, whereas it may increase the risk of subsidence in flexion and rotation.     

Sitting versus standing: Does the intradiscal pressure cause disc degeneration or low back pain?

Studies of lumbar intradiscal pressure (IDP) in standing and upright sitting have mostly reported higher pressures in sitting. It was assumed clinically that flexion of the lumbar spine in sitting relative to standing, caused higher IDP, disc degeneration or rupture, and low back pain. IDP indicates axial compressive load upon a non-degenerate disc, but provides little or no indication of shear, axial rotation or bending. This review is presented in two main parts. First, in vivo IDP data in standing and upright sitting for non-degenerate discs are comprehensively reviewed. As methodology, results and interpretations varied between IDP studies, in vivo studies measuring spinal shrinkage and spinal internal-fixator loads to infer axial compressive load to the discs are also reviewed. When data are considered together, it is clear that IDP is often similar in standing and sitting. Secondly, clinical assumptions related to IDP in sitting are considered in light of basic and epidemiologic studies. Current studies indicate that IDP in sitting is unlikely to pose a threat to non-degenerate discs, and sitting is no worse than standing for disc degeneration or low back pain incidence. If sitting is a greater threat for development of low back pain than standing, the mechanism is unlikely to be raised IDP.     

Load-displacement properties of lower cervical spine motion segments

The load-displacement behavior of 35 fresh adult cervical spine motion segments was measured in compression, shear, flexion, extension, lateral bending and axial torsion tests. Motion segments were tested both intact and with posterior elements removed. Applied forces ranged to 73.6 N in compression and to 39 N in shear, while applied moments ranged to 2.16 Nm. For each mode of loading, principal and coupled motions were measured and stiffnesses were calculated. The effect of disc degeneration on motion segment stiffnesses and the moments required for motion segment failure were also measured. In compression, the stiffnesses of the cervical motion segments were similar to those of thoracic and lumbar motion segments. In other modes of loading, cervical stiffnesses were considerably smaller than thoracic or lumbar stiffnesses. Removal of the posterior elements decreased cervical motion segment stiffnesses by as much as 50%. Degenerated cervical discs were less stiff in compression and stiffer in shear than less degenerated discs, but in bending or axial torsion, no statistically significant differences were evident. Bending moments causing failure were an order of magnitude lower than those for lumbar segments.     

The influence of different magnitudes and methods of applying preload on fusion and disc replacement constructs in the lumbar spine: a finite element analysis

In a finite element (FE) analysis of the lumbar spine, different preload application methods that are used in biomechanical studies may yield diverging results. To investigate how the biomechanical behaviour of a spinal implant is affected by the method of applying the preload, hybrid-controlled FE analysis was used to evaluate the biomechanical behaviour of the lumbar spine under different preload application methods. The FE models of anterior lumbar interbody fusion (ALIF) and artificial disc replacement (ADR) were tested under three different loading conditions: a 150 N pressure preload (PP) and 150 and 400 N follower loads (FLs). This study analysed the resulting range of motion (ROM), facet contact force (FCF), inlay contact pressure (ICP) and stress distribution of adjacent discs. The FE results indicated that the ROM of both surgical constructs was related to the preload application method and magnitude; differences in the ROM were within 7% for the ALIF model and 32% for the ADR model. Following the application of the FL and after increasing the FL magnitude, the FCF of the ADR model gradually increased, reaching 45% at the implanted level in torsion. The maximum ICP gradually decreased by 34.1% in torsion and 28.4% in lateral bending. This study concluded that the preload magnitude and application method affect the biomechanical behaviour of the lumbar spine. For the ADR, remarkable alteration was observed while increasing the FL magnitude, particularly in the ROM, FCF and ICP. However, for the ALIF, PP and FL methods had no remarkable alteration in terms of ROM and adjacent disc stress.     

Biomechanics of the C5-C6 Spinal Unit Before and After Placement of a disc prosthesis

The study consists of a biomechanical comparison between the intact C5–C6 spinal segment and the same segment implanted with the Bryan^TM artificial disc prosthesis (Medtronic Ltd., Memphis, TN, USA), by the use of the finite element (FE) method. Our target is the prediction of the influence of prosthesis placement on the resulting mechanics of the C5–C6 spine unit. A FE model of the intact C5–C6 segment was built, employing realistic models of the vertebrae, disc and ligaments. Simulations were conducted imposing a compression preload combined to a flexion/extension moment, a pure lateral bending moment and a pure torsion moment, and the calculated results were compared to data from literature. The model was then modified to include the Bryan^TM cervical disc prosthesis, and the simulations were repeated. The location of the instantaneous center of rotation (ICR) of C5 with respect to C6 throughout flexion/extension was calculated in both models. In general, the moment–rotation curves obtained from the disc prosthesis-implanted model were comparable to the curves obtained from the intact model, except for a slightly greater stiffness induced by the artificial disc. The position of the calculated ICRs was rather stable throughout flexion-extension and was generally confined to a small area, qualitatively matching the corresponding physiological region, in both models. These results imply that the Bryan^TM disc prosthesis allows to correctly reproduce a physiological flexion/extension at the implanted level. The results of this study have quantified aspects that may assist in optimizing cervical disc replacement primarily from a biomechanical point of view.     

Are coupled rotations in the lumbar spine largely due to the osseo-ligamentous anatomy?--a modeling study

Prior studies have found that primary rotations in the lumbar spine are accompanied by coupled out-of-plane rotations. However, it is not clear whether these accompanying rotations are primarily due to passive (discs, ligaments and facet joints) or active (muscles) spinal anatomy. The aim of this study was to use a finite element (FE) model of the lumbar spine to predict three-dimensional coupled rotations between the lumbar vertebrae, due to passive spinal structures alone. The FE model was subjected to physiologically observed whole lumbar spine rotations about in vivo centres of rotation. Model predictions were validated by comparison of intra-discal pressures and primary rotations with in vivo measurements and these showed close agreement. Predicted coupled rotations matched in vivo measurements for all primary motions except lateral bending. We suggest that coupled rotations accompanying primary motions in the sagittal (flexion/extension) and transverse (axial rotation) planes are primarily due to passive spinal structures. For lateral bending the muscles most likely play a key role in the coupled rotation of the spine.     

Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments—A finite element model study

The current study investigated mechanical predictors for the development of adjacent disc degeneration. A 3-D finite element model of a lumbar spine was modified to simulate two grades of degeneration at the L4–L5 disc. Degeneration was modeled by changes in geometry and material properties. All models were subjected to follower preloads of 800 N and moment loads in the three principal directions of motion using a hybrid protocol. Degeneration caused changes in the loading and motion patterns of the segments above and below the degenerated disc. At the level (L3–L4) above the degenerated disc, the motion increased due to moderate degeneration by 21% under lateral bending, 26% under axial rotation and 28% under flexion/extension. At the level (L5-S1) below the degenerated disc, motion increased only during lateral bending by 20% due to moderate degeneration. Both the L3–L4 and L5-S1 segment showed a monotonic increase in both the maximum von Mises stress and shear stress in the annulus as degeneration progressed for all loading directions, expect extension at L3–L4. The most significant increase in stress was observed at the L5-S1 level during axial rotation with nearly a ten-fold increase in the maximum shear stress and 103% increase in the maximum von Mises stress. The L5-S1 segment also showed a progressive increase in facet contact force for all loading directions with degeneration. Nucleus pressure did not increase significantly for any loading direction at either the caudal or cephalic adjacent segment. Results suggest that single-level degeneration can increase the risk for injury at the adjacent levels.     

Direct measurement of intervertebral disc maximum shear strain in six degrees of freedom: motions that place disc tissue at risk of injury

Human intervertebral disc specimens were tested to determine the regions of largest maximum shear strain (MSS) experienced by disc tissues in each of three principal displacements and three rotations, and to identify the physiological rotations and displacements that may place the disc at greatest risk for large tissue strains and injury. Tearing of disc annulus may be initiated by large interlamellar shear strains. Nine human lumbar discs were tagged with radiographic markers on the endplates, disc periphery and with a grid of wires in the mid-transverse plane and subjected to each of the six principal displacements and rotations. Stereo-radiographs were taken in each position and digitized for reconstruction of the three-dimensional position of each marker. Maximum tissue shear strains were calculated from relative marker displacements and normalized by the input displacement or rotation. Lateral shear, compression, and lateral bending were the motions that produced the mean (95% confidence interval) largest mean MSS of 9.6 (0.7)%/mm, 9.0 (0.5)%/mm, and 5.8 (1.6)%/ degrees , respectively, and which occurred in the posterior, posterolateral and lateral peripheral regions of the disc. After taking into account the reported maximum physiological range of motion for each degree of freedom, motions producing the highest physiological MSS were lateral bending (57.8 (16.2)%) and flexion (38.3 (3.3)%), followed by lateral shear (14.4 (1.1)%) and compression (12.6 (0.7)%).     

Smart garment for trunk posture monitoring: A preliminary study

Background

Hueter-Volkmann's law regarding growth modulation suggests that increased pressure on the end plate of bone retards the growth (Hueter) and conversely, reduced pressure accelerates the growth (Volkmann). Literature described the same principle in Rat-tail model. Human spine and its deformity i.e. scoliosis has also same kind of pattern during the growth period which causes wedging in disc or vertebral body.

Methods

This cross sectional study in 150 patients of adolescent idiopathic scoliosis was done to evaluate vertebral body and disc wedging in scoliosis and to compare the extent of differential wedging of body and disc, in thoracic and lumbar area. We measured wedging of vertebral bodies and discs, along with two adjacent vertebrae and disc, above and below the apex and evaluated them according to severity of curve (curve < 30° and curve > 30°) to find the relationship of vertebral body or disc wedging with scoliosis in thoracic and lumbar spine. We also compared the wedging and rotations of vertebrae.

Results

In both thoracic and lumbar curves, we found that greater the degree of scoliosis, greater the wedging in both disc and body and the degree of wedging was more at apex supporting the theory of growth retardation in stress concentration area. However, the degree of wedging in vertebral body is more than the disc in thoracic spine while the wedging was more in disc than body in lumbar spine. On comparing the wedging with the rotation, we did not find any significant relationship suggesting that it has no relation with rotation.

Conclusion

From our study, we can conclude that wedging in disc and body are increasing with progression on scoliosis and maximum at apex; however there is differential wedging of body and disc, in thoracic and lumbar area, that is vertebral body wedging is more profound in thoracic area while disc wedging is more profound in lumbar area which possibly form 'vicious cycle' by asymmetric loading to spine for the progression of curve.     

Morphometric analysis of vertebrae and intervertebral discs as a basis of disc replacement

In this study measurements were obtained from 359 vertebrae and 215 intervertebral discs in an attempt to classify discs by their size. At the cervical and thoracic levels, this attempt was unproductive because of extensive variations. In the lumbar spine, discs were allocated to six size-matched groups and to two height-matched groups. The breadth of marginal rims were also measured with a view to provide surgeons operating on the spine with precise data on disc morphology to facilitate disc replacement.     

Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites

Abstract

The current paper aims at assessing the sensitivity of muscle and intervertebral disc force computations against potential errors in modeling muscle attachment sites. We perturbed each attachment location in a complete and coherent musculoskeletal model of the human spine and quantified the changes in muscle and disc forces during standing upright, flexion, lateral bending, and axial rotation of the trunk. Although the majority of the muscles caused minor changes (less than 5%) in the disc forces, certain muscle groups, for example, quadratus lumborum, altered the shear and compressive forces as high as 353% and 17%, respectively. Furthermore, percent changes were higher in the shear forces than in the compressive forces. Our analyses identified certain muscles in the rib cage (intercostales interni and intercostales externi) and lumbar spine (quadratus lumborum and longissimus thoracis) as being more influential for computing muscle and disc forces. Furthermore, the disc forces at the L4/L5 joint were the most sensitive against muscle attachment sites, followed by T6/T7 and T12/L1 joints. Presented findings suggest that modeling muscle attachment sites based on solely anatomical illustrations might lead to erroneous evaluation of internal forces and promote using anatomical datasets where these locations were accurately measured. When developing a personalized model of the spine, certain care should also be paid especially for the muscles indicated in this work.     

Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression,flexion and axial rotation

A statistical factorial analysis approach was conducted on a poroelastic finite element model of a lumbar intervertebral disc to analyse the influence of six material parameters (permeabilities of annulus, nucleus, trabecular vertebral bone, cartilage endplate and Young's moduli of annulus and nucleus) on the displacement, fluid pore pressure and velocity fields. Three different loading modes were investigated: compression, flexion and axial rotation. Parameters were varied considering low and high levels in agreement with values found in the literature for both healthy and degenerated lumbar discs. Results indicated that annulus stiffness and cartilage endplate permeability have a strong effect on the overall fluid- and solid-phase responses in all loading conditions studied. Nucleus stiffness showed its main relevance in compression while annulus permeability influenced mainly the annular pressure field. This study confirms the permeability's central role in biphasic modelling and highlights for the lumbar disc which experiments of material property characterization should be performed. Moreover, such sensitivity study gives important guidelines in poroelastic material modelling and finite element disc validation.     

Ultimate strength of the lumbar spine in flexion--an in vitro study

The ultimate strength in flexion of 16 lumbar functional spinal units (FSU) was determined. The specimens were exposed to a combined static load of bending and shearing in the sagittal plane until overt rupture occurred (simulated flexion-distraction injuries). The biomechanical response of the FSU was measured with a force and moment platform. Mechanical displacement gauges were used to measure vertical displacements (flexion angulation) of the specimens. Photographs were taken after each loading step for determination of horizontal displacements and the centre of rotation. The lumbar FSU could resist a combination of bending moment and shear force of 156 Nm and 620 N respectively, before complete disruption occurred. The tension force acting on the posterior structures was 2.8 kN. The flexion angulation just before failure was 20 degrees and the anterior horizontal displacement between the upper and lower vertebrae was 9 mm. The centre of rotation was located in the posterior part of the lower vertebral body. The bone mineral content in the vertebrae appeared to be a good predictor of ultimate strength of the lumbar FSU. Knowledge of the biomechanical response of the lumbar spine under different static traumatic loads is a first step to better understand the injury mechanisms of the spine in traffic accidents.     

Effects of lumbo-pelvic rhythm on trunk muscle forces and disc loads during forward flexion: A combined musculoskeletal and finite element simulation study

Previous in-vivo studies suggest that the ratio of total lumbar rotation over pelvic rotation (lumbo-pelvic rhythm) during trunk sagittal movement is essential to evaluate spinal loads and discriminate between low back pain and asymptomatic population. Similarly, there is also evidence that the lumbo-pelvic rhythm is key for evaluation of realistic muscle and joint reaction forces and moments predicted by various computational musculoskeletal models. This study investigated the effects of three lumbo-pelvic rhythms defined based on in-vivo measurements on the spinal response during moderate forward flexion (60°) using a combined approach of musculoskeletal modeling of the upper body and finite element model of the lumbosacral spine. The muscle forces and joint loads predicted by the musculoskeletal model, together with the gravitational forces, were applied to the finite element model to compute the disc force and moment, intradiscal pressure, annular fibers strain, and load-sharing. The results revealed that a rhythm with high pelvic rotation and low lumbar flexion involves more global muscles and increases the role of the disc in resisting spinal loads, while its counterpart, with low pelvic rotation, recruits more local muscles and engages the ligaments to lower the disc loads. On the other hand, a normal rhythm that has balanced pelvic and lumbar rotations yields almost equal disc and ligament load-sharing and results in more balanced synergy between global and local muscles. The lumbo-pelvic rhythm has less effect on the intradiscal pressure and annular fibers strain. This work demonstrated that the spinal response during forward flexion is highly dependent on the lumbo-pelvic rhythm. It is therefore, essential to adapt this parameter instead of using the default values in musculoskeletal models for accurate prediction of muscle forces and joint reaction forces and moments. The findings provided by this work are expected to improve knowledge of spinal response during forward flexion, and are clinically relevant towards low back pain treatment and disc injury prevention.     

Prediction of biomechanical parameters in the lumbar spine during static sagittal plane lifting.

A combined approach involving optimization and the finite element technique was used to predict biomechanical parameters in the lumbar spine during static lifting in the sagittal plane. Forces in muscle fascicles of the lumbar region were first predicted using an optimization-based force model including the entire lumbar spine. These muscle forces as well as the distributed upper body weight and the lifted load were then applied to a three-dimensional finite element model of the thoracolumbar spine and rib cage to predict deformation, the intradiskal pressure, strains, stresses, and load transfer paths in the spine. The predicted intradiskal pressures in the L3-4 disk at the most deviated from the in vivo measurements by 8.2 percent for the four lifting cases analyzed. The lumbosacral joint flexed, while the other lumbar joints extended for all of the four lifting cases studied (rotation of a joint is the relative rotation between its two vertebral bodies). High stresses were predicted in the posterolateral regions of the endplates and at the junctions of the pedicles and vertebral bodies. High interlaminar shear stresses were found in the posterolateral regions of the lumbar disks. While the facet joints of the upper two lumbar segments did not transmit any load, the facet joints of the lower two lumbar segments experienced significant loads. The ligaments of all lumbar motion segments except the lumbosacral junction provided only marginal moments. The limitations of the current model and possible improvements are discussed.     

Biomechanical implications of lumbar spinal ligament transection

Many lumbar spine surgeries either intentionally or inadvertently damage or transect spinal ligaments. The purpose of this work was to quantify the previously unknown biomechanical consequences of isolated spinal ligament transection on the remaining spinal ligaments (stress transfer), vertebrae (bone remodelling stimulus) and intervertebral discs (disc pressure) of the lumbar spine. A finite element model of the full lumbar spine was developed and validated against experimental data and tested in the primary modes of spinal motion in the intact condition. Once a ligament was removed, stress increased in the remaining spinal ligaments and changes occurred in vertebral strain energy, but disc pressure remained similar. All major biomechanical changes occurred at the same spinal level as the transected ligament, with minor changes at adjacent levels. This work demonstrates that iatrogenic damage to spinal ligaments disturbs the load sharing within the spinal ligament network and may induce significant clinically relevant changes in the spinal motion segment.