Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle

Many investigators have performed studies on specific defect situations or determined the contribution on isolated structures. Investigating the contribution of functional structures requires obtaining the kinematic response directly on spinal segments. The purpose of this study was to quantify the function of anatomical components on lumbar segments for different loading magnitudes. Eight spinal segments (L4-5) with a median age of 52 years (ranging from 38 to 59 years) and a low degree of disc degeneration were utilized for the in vitro testing. Specimens were mounted in a custom-built spine tester and loaded with pure moments (1-10 N m) to move within three anatomical planes at a loading rate of 1.0 degrees /s. Anatomy was successively reduced by: ligaments, facet capsules, joints and nucleus. Data were evaluated for range of motion, neutral zone and lordosis angle. Transection of posterior ligaments predominantly increased specimen flexion for all bending moments applied. Supraspinous ligament also indicated to resist in extension slightly, whereas the facet capsules did not. Facet joints contributed to axial rotation, but not in lateral bending. The anterior longitudinal ligament was found to slightly resist in axial rotation, but strongly in extension. Nucleotomy caused largest increase of all movements. The unloaded posture of the specimens changed after ligament dissection, indicating ligament pretension. The region of lumbar spine is interesting for finite element (FE) simulation due to the high evidence of disc degeneration and injuries. This study may help to understand the function of specific anatomical structures and assists in FE model calibration. We suggest to start a calibration procedure for such models with the smallest functional structure (annulus) and to cumulatively add further structures.     

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.     

The effect of follower load on the intersegmental coupled motion characteristics of the human thoracic spine: An in vitro study using entire rib cage specimens

The mechanical coupling behaviour of the thoracic spine is still not fully understood. For the validation of numerical models of the thoracic spine, however, the coupled motions within the single spinal segments are of importance to achieve high model accuracy. In the present study, eight fresh frozen human thoracic spinal specimens (C7-L1, mean age 54 ± 6 years) including the intact rib cage were loaded with pure bending moments of 5 Nm in flexion/extension (FE), lateral bending (LB), and axial rotation (AR) with and without a follower load of 400 N. During loading, the relative motions of each vertebra were monitored. Follower load decreased the overall ROM (T1-T12) significantly (p < 0.01) in all primary motion directions (extension: −46%, left LB: −72%, right LB: −72%, left AR: −26%, right AR: −26%) except flexion (−36%). Substantial coupled motion was found in lateral bending with ipsilateral axial rotation, which increased after a follower load was applied, leading to a dominant axial rotation during primary lateral bending, while all other coupled motions in the different motion directions were reduced under follower load. On the monosegmental level, the follower load especially reduced the ROM of the upper thoracic spine from T1-T2 to T4-T5 in all motion directions and the ROM of the lower thoracic spine from T9-T10 to T11-T12 in primary lateral bending. The facet joints, intervertebral disc morphologies, and the sagittal curvature presumably affect the thoracic spinal coupled motions depending on axial compressive preloading. Using these results, the validation of numerical models can be performed more accurately.     

How does the geometry affect the internal biomechanics of a lumbar spine bi-segment finite element model? Consequences on the validation process

Numerical modelling can provide a thorough understanding of the mechanical influence on the spinal tissues and may offer explanations to mechanically linked pathologies. Such objective might be achieved only if the models are carefully validated. Sensitivity study must be performed in order to evaluate the influence of the approximations inherent to modelling. In this study, a new geometrically accurate L3-L5 lumbar spine bi-segmental finite-element model was acquired by modifying a previously existing model. The effect of changes in bone geometry, ligament fibres distribution, nucleus position and disc height was investigated in flexion and extension by comparison of the results obtained from the model before and after the geometrical update. Additional calculations were performed in axial rotation and lateral bending in order to compare the computed ranges of motion (ROM) with experimental results. It was found that the geometrical parameters affected the stress distribution and strain energy in the zygapophysial joints, the ligaments, and the intervertebral disc, changing qualitatively and quantitatively their relative role in resisting the imposed loads. The predicted ROM were generally in good agreement with the experimental results, independently of the geometrical changes. Hence, although the model update affected its internal biomechanics, no conclusions could be drawn from the experimental data about the validation of a particular geometry. Hence the validation of the lumbar spine model should be based on the relative role of its structural components and not only on its global mobility.     

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.     

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.     

Novel model to analyze the effect of a large compressive follower pre-load on range of motions in a lumbar spine

A 3-D finite element model (FEM) of the lumbar spine (L1-S1) was used to determine the effect of a large compressive follower pre-load on range of motions (ROM) in all three planes. The follower load modeled in the FEM produced minimal vertebral rotations in all the three planes. The model was validated by comparing the disc compression at all levels in the lumbar spine with the corresponding results obtained by compressing 10 cadevaric lumbar spines (L1-S1) using the follower load technique described by Patwardhan et al. [1999. A follower load increases the load-carrying capacity of the lumbar spine in compression. Spine 24(10), 1003-1009]. Further validation of the model was performed by comparing the lateral bending and torsion response without pre-load and the flexion-extension response without pre-load and with an 800 N follower pre-load with those obtained using cadaver lumbar spines. Following validation, the FEM was subjected to bending moments in all three planes with and without compressive follower pre-loads of up to 1200 N. Disc compression values and the flexion-extension range of motion under 800 N follower pre-load predicted by the FEM compared well with in vitro results. The current model showed that compressive follower pre-load decreased total as well as segmental ROM in flexion-extension by up to 18%, lateral bending by up to 42%, and torsion by up to 26%.     

The role of the facet capsular ligament in providing spinal stability

Abstract

Low back pain (LBP) is the most common type of pain in America, and spinal instability is a primary cause. The facet capsular ligament (FCL) encloses the articulating joints of the spine and is of particular interest due to its high innervation – as instability ensues, high stretch values likely are a cause of this pain. Therefore, this work investigated the FCL's role in providing stability to the lumbar spine. A previously validated finite element model of the L4-L5 spinal motion segment was used to simulate pure moment bending in multiple planes. FCL failure was simulated and the following outcome measures were calculated: helical axes of motion, range of motion (ROM), bending stiffness, facet joint space, and FCL stretch. ROM increased, bending stiffness decreased, and altered helical axis patterns were observed with the removal of the FCL. Additionally, a large increase in FCL stretch was measured with diminished FCL mechanical competency, providing support that the FCL plays an important role in spinal stability.     

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.     

Quantitative comparison of ligament formulation and pre-strain in finite element analysis of the human lumbar spine

Data has been published that quantifies the nonlinear, anisotropic material behaviour and pre-strain behaviour of the anterior longitudinal, supraspinous (SSL), and interspinous ligaments of the human lumbar spine. Additionally, data has been published on localized material properties of the SSL. These results have been incrementally incorporated into a previously validated finite element model of the human lumbar spine. Results suggest that the effects of increased ligament model fidelity on bone strain energy were moderate and the effects on disc pressure were slight, and do not justify a change in modelling strategy for most clinical applications. There were significant effects on the ligament stresses of the ligaments that were directly modified, suggesting that these phenomena should be included in FE models where ligament stresses are the desired metric.     

Analysis of simulated single ligament transection on the mechanical behaviour of a lumbar functional spinal unit.

The influence of the different lumbar spinal ligaments on intersegmental rotation is not fully understood. In order to explore this effect, a finite element model of the functional spinal unit L3/L4 was loaded with pure moments in the three main anatomic planes. The two extremes--minimum and maximum--ligament stiffness values reported in the literature were applied. After virtual transection of each of the spinal ligaments in turn, the intersegmental rotation and forces in the remaining ligaments were calculated. On flexion, the highest force was found for the posterior longitudinal ligament; on extension and lateral bending for the anterior longitudinal ligament; and on axial rotation for the facet capsular ligament. The strongest influence on intersegmental rotation is exerted by the interspinous ligament on flexion, by the anterior longitudinal ligament on extension and lateral bending, and by the facet capsular ligaments on axial rotation. Ligament stiffness has a strong influence on intersegmental rotation and forces in the ligaments, so that finite element models of spinal segments must be validated by experimental data. This study should help to elucidate the role of the various ligaments.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

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.     

The effect of three-dimensional geometrical changes during adolescent growth on the biomechanics of a spinal motion segment

During adolescent growth, vertebrae and intervertebral discs undergo various geometrical changes. Although such changes in geometry are well known, their effects on spinal stiffness remains poorly understood. However, this understanding is essential in the treatment of spinal abnormalities during growth, such as scoliosis.A finite element model of an L3–L4 motion segment was developed, validated and applied to study the quantitative effects of changing geometry during adolescent growth on spinal stiffness in flexion, extension, lateral bending and axial rotation. Height, width and depth of the vertebrae and intervertebral disc were varied, as were the width of the transverse processes, the length of the spinous process, the size of the nucleus, facet joint areas and ligament size. These variations were based on average growth data for girls, as reported in literature.Overall, adolescent growth increases the stiffness with 36% (lateral bending and extension) to 44% (flexion). Two thirds of this increase occurs between 10 and 14 years of age and the last third between 14 years of age and maturity.Although the height is the largest geometrical change during adolescent growth, its effect on the biomechanics is small. The depth increase of the disc and vertebrae significantly affects the stiffness in all directions, while the width increase mainly affects the lateral bending stiffness. Hence, when analysing the biomechanics of the growing adolescent spine (for instance in scoliosis research), the inclusion of depth and width changes, in addition to the usually implemented height change, is essential.     

Effects of the rib cage on thoracic spine flexibility.

Besides protecting the internal organs of the thorax, the rib cage is the site of numerous muscle attachments. It also decreases the overall flexibility of the thoracic spine. This study developed finite element (FE) models of the thoracic spine with and without the rib cage, and the effects of the rib cage on thoracic spine flexibility were determined. The numerical models were validated by comparing the maximum rotation of the models for several loading cases with experimental data in the literature. After adapting the material properties for the discs and ligaments, the calculated maximum rotations differed from the measured median values by less than 1 degrees without the rib cage and by less than 2.5 degrees with it. The rib cage decreased the mean flexibility of the thoracic spine by 23% to 47%, depending on the loading plane. Assuming the ribs to be rigid beams required a corresponding reduction of ligament stiffnesses in order to achieve the same agreement of the maximum rotations with the measured median values. Interconnecting the FE thoracic spine model plus rib cage with the existing detailed FE lumbar spine model improves the simulation of force directions of muscles attached to the rib cage or thoracolumbar spine. In addition, such a model is suitable for determining the effects of lumbar spine implants on spinal balance.     

Finite Element Analysis of a New Pedicle Screw-Plate System for Minimally Invasive Transforaminal Lumbar Interbody Fusion

Purpose

Minimally invasive transforaminal lumbar interbody fusion (MI-TLIF) is increasingly popular for the surgical treatment of degenerative lumbar disc diseases. The constructs intended for segmental stability are varied in MI-TLIF. We adopted finite element (FE) analysis to compare the stability after different construct fixations using interbody cage with posterior pedicle screw-rod or pedicle screw-plate instrumentation system.

Methods

A L3–S1 FE model was modified to simulate decompression and fusion at L4–L5 segment. Fixation modes included unilateral plate (UP), unilateral rod (UR), bilateral plate (BP), bilateral rod (BR) and UP+UR fixation. The inferior surface of the S1 vertebra remained immobilized throughout the load simulation, and a bending moment of 7.5 Nm with 400N pre-load was applied on the L3 vertebra to recreate flexion, extension, lateral bending, and axial rotation. Range of motion (ROM) and Von Mises stress were evaluated for intact and instrumentation models in all loading planes.

Results

All reconstructive conditions displayed decreased motion at L4–L5. The pedicle screw-plate system offered equal ROM to pedicle screw-rod system in unilateral or bilateral fixation modes respectively. Pedicle screw stresses for plate system were 2.2 times greater than those for rod system in left lateral bending under unilateral fixation. Stresses for plate were 3.1 times greater than those for rod in right axial rotation under bilateral fixation. Stresses on intervertebral graft for plate system were similar to rod system in unilateral and bilateral fixation modes respectively. Increased ROM and posterior instrumentation stresses were observed in all loading modes with unilateral fixation compared with bilateral fixation in both systems.

Conclusions

Transforaminal lumbar interbody fusion augmentation with pedicle screw-plate system fixation increases fusion construct stability equally to the pedicle screw-rod system. Increased posterior instrumentation stresses are observed in all loading modes with plate fixation, and bilateral fixation could reduce stress concentration.     

Biomechanical contributions of upper cervical ligamentous structures in Type II odontoid fractures

Fractures of the odontoid present frequently in spinal trauma, and Type II odontoid fractures, occurring at the junction of the odontoid process and C2 vertebrae, represent the bulk of all traumatic odontoid fractures. It is currently unclear what soft-tissue stabilizers contribute to upper cervical motion in the setting of a Type II odontoid fracture, and evaluation of how concomitant injury contributes to cervical stability may inform surgical decision-making as well as allow for the creation of future, accurate, biomechanical models of the upper cervical spine. The objective of the current study was to determine the contribution of soft-tissue stabilizers in the upper cervical spine following a Type II odontoid fracture. Eight cadaveric C0-C2 specimens were evaluated using a robotic testing system with motion tracking. The unilateral facet capsule (UFC) and anterior longitudinal ligament (ALL) were serially resected to determine their biomechanical role following odontoid fracture. Range of motion (ROM) and moment at the end of intact specimen replay were the primary outcomes. We determined that fracture of the odontoid significantly increases motion and decreases resistance to intact motion for flexion–extension (FE), axial rotation (AR), and lateral bending (LB). Injury to the UFC increased AR by 3.2° and FE by 3.2°. ALL resection did not significantly increase ROM or decrease end-point moment. The UFC was determined to contribute to 19% of intact flexion resistance and 24% of intact AR resistance. Overall, we determined that Type II fracture of the odontoid is a significant biomechanical destabilizer and that concurrent injury to the UFC further increases upper cervical ROM and decreases resistance to motion in a cadaveric model of traumatic Type II odontoid fractures.     

Between-session reliability of opto-electronic motion capture in measuring sagittal posture and 3-D ranges of motion of the thoracolumbar spine

This study evaluated between-session reliability of opto-electronic motion capture to measure trunk posture and three-dimensional ranges of motion (ROM). Nineteen healthy participants aged 24–74 years underwent spine curvature, pelvic tilt and trunk ROM measurements on two separate occasions. Rigid four-marker clusters were attached to the skin overlying seven spinous processes, plus single markers on pelvis landmarks. Rigid body rotations of spine marker clusters were calculated to determine neutral posture and ROM in flexion, extension, total lateral bending (left-right) and total axial rotation (left-right). Segmental spine ROM values were in line with previous reports using opto-electronic motion capture. Intraclass correlation coefficients (ICC) and standard error of measurement (SEM) were calculated as measures of between-session reliability and measurement error, respectively. Retroreflective markers showed fair to excellent between-session reliability to measure thoracic kyphosis, lumbar lordosis, and pelvic tilt (ICC = 0.82, 0.63, and 0.54, respectively). Thoracic and lumbar segments showed highest reliabilities in total axial rotation (ICC = 0.78) and flexion-extension (ICC = 0.77–0.79) ROM, respectively. Pelvic segment showed highest ICC values in flexion (ICC = 0.78) and total axial rotation (ICC = 0.81) trials. Furthermore, it was estimated that four or fewer repeated trials would provide good reliability for key ROM outcomes, including lumbar flexion, thoracic and lumbar lateral bending, and thoracic axial rotation. This demonstration of reliability is a necessary precursor to quantifying spine kinematics in clinical studies, including assessing changes due to clinical treatment or disease progression.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion-extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion–extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.