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.     

Experimental analysis of the lower cervical spine in flexion with a focus on facet tracking

Cervical traumas are among the most common events leading to serious spinal cord injuries. While models are often used to better understand injury mechanisms, experimental data for their validation remain sparse, particularly regarding articular facets. The aim of this study was to assess the behavior of cervical FSUs under quasi-static flexion with a specific focus on facet tracking. 9 cadaveric cervical FSUs were imaged and loaded under a 10 Nm flexion moment, exerted incrementally, while biplanar X-rays were acquired at each load increment. The relative vertebral and facet rotations and displacements were assessed using radio-opaque markers implanted in each vertebra and CT-based reconstructions registered on the radiographs. The only failures obtained were due to specimen preparation, indicating a failure moment of cervical FSUs greater than 10 Nm in quasistatic flexion. Facet motions displayed a consistent anterior sliding and a variable pattern regarding their normal displacement. The present study offers insight on the behavior of cervical FSUs under quasi-static flexion beyond physiological thresholds with accurate facet tracking. The data provided should prove useful to further understand injury mechanisms and validate models.     

Moment-rotation relationships of the ligamentous occipito-atlanto-axial complex

The relationships between applied pure moments at the occiput (C0) and the resulting rotations at the atlanto-occipital (C0-C1) and atlanto-axial (C1-C2) joints are quantified. In axial twist, with a moment of 0.3 Nm, a mean rotation of about 2.5 degrees and 23.3 degrees was observed at C0-C1 and C1-C2 units respectively. Both the atlas and axis contributed to produce lateral bending motion. The ratio between extension and flexion rotations at C0-C1 was 2.5:1. Lateral bending and axial rotations were strongly coupled to each other. The occipito-atlanto-axial complex exhibited a large 'neutral zone' compared to lower cervical spine segments. The likely clinical significance of these findings are discussed.     

An experimental method for measuring force on the spinal facet joint: description and application of the method.

A technique is described for measuring load magnitude and resultant load contact location in the facet joint in response to applied loads and moments, and the technique applied to the canine lumbar spine motion segment. Due to the cantilever beam geometry of the cranial articular process, facet joint loads result in surface strains on the lateral aspect of the cranial articular process. Strains were quantified by four strain gages cemented to the bony surface of the process. Strain measured at any one gage depended on the loading site on the articular surface of the caudal facet and on the magnitude of the facet load. Determination of facet loads during in vitro motion segment testing required calibration of the strains to known loads of various magnitudes applied to multiple sites on the caudal facet. The technique is described in detail, including placement of the strain gages. There is good repeatability of strains to applied facet loads and the strains appear independent of load distribution area. Error in the technique depends on the location of the applied facet loads, but is only significant in nonphysiologic locations. The technique was validated by two independent methods in axial torsion. Application of the technique to five in vitro canine L2-3 motion segments testing resulted in facet loads (in newtons, N) of 74+ / -23 N (mean + / -STD) in 2 newton-meter, Nm, extension, to unloaded in flexion. Lateral bending resulted in loads in the right facet of 40+ / -32 N for 1 Nm right lateral bending and 54+ / -29 N for 1 Nm left lateral bending. 4 Nm Torsion with and without 100 N axial compression resulted in facet loads of 92+ / -27 N and 69+ / -19 N, respectively. The technique is applicable to dynamic and in vivo studies.     

Comparative strengths and structural properties of the upper and lower cervical spine in flexion and extension

The purpose of this study is to test the hypothesis that the upper cervical spine is weaker than the lower cervical spine in pure flexion and extension bending, which may explain the propensity for upper cervical spine injuries in airbag deployments. An additional objective is to evaluate the relative strength and flexibility of the upper and lower cervical spine in an effort to better understand injury mechanisms, and to provide quantitative data on bending responses and failure modes. Pure moment flexibility and failure testing was conducted on 52 female spinal segments in a pure-moment test frame. The average moment at failure for the O-C2 segments was 23.7+/-3.4Nm for flexion and 43.3+/-9.3Nm for extension. The ligamentous upper cervical spine was significantly stronger in extension than in flexion (p=0.001). The upper cervical spine was significantly stronger than the lower cervical spine in extension. The relatively high strength of the upper cervical spine in tension and in extension is paradoxical given the large number of upper cervical spine injuries in out-of-position airbag deployments. This discrepancy is most likely due to load sharing by the active musculature.     

Investigation of impact loading rate effects on the ligamentous cervical spinal load-partitioning using finite element model of functional spinal unit C2–C3

The cervical spine functions as a complex mechanism that responds to sudden loading in a unique manner, due to intricate structural features and kinematics. The spinal load-sharing under pure compression and sagittal flexion/extension at two different impact rates were compared using a bio-fidelic finite element (FE) model of the ligamentous cervical functional spinal unit (FSU) C2–C3. This model was developed using a comprehensive and realistic geometry of spinal components and material laws that include strain rate dependency, bone fracture, and ligament failure. The range of motion, contact pressure in facet joints, failure forces in ligaments were compared to experimental findings. The model demonstrated that resistance of spinal components to impact load is dependent on loading rate and direction. For the loads applied, stress increased with loading rate in all spinal components, and was concentrated in the outer intervertebral disc (IVD), regions of ligaments to bone attachment, and in the cancellous bone of the facet joints. The highest stress in ligaments was found in capsular ligament (CL) in all cases. Intradiscal pressure (IDP) in the nucleus was affected by loading rate change. It increased under compression/flexion but decreased under extension. Contact pressure in the facet joints showed less variation under compression, but increased significantly under flexion/extension particularly under extension. Cancellous bone of the facet joints region was the only component fractured and fracture occurred under extension at both rates. The cervical ligaments were the primary load-bearing component followed by the IVD, endplates and cancellous bone; however, the latter was the most vulnerable to extension as it fractured at low energy impact.     

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.     

Flexion and extension structural properties and strengths for male cervical spine segments

New vehicle safety standards are designed to limit the amount of neck tension and extension seen by out-of-position motor vehicle occupants during airbag deployments. The criteria used to assess airbag injury risk are currently based on volunteer data and animal studies due to a lack of bending tolerance data for the adult cervical spine. This study provides quantitative data on the flexion-extension bending properties and strength on the male cervical spine, and tests the hypothesis that the male is stronger than the female in pure bending. An additional objective is to determine if there are significant differences in stiffness and strength between the male upper and lower cervical spine. Pure-moment flexibility and failure testing was conducted on 41 male spinal segments (O-C2, C4-C5, C6-C7) in a pure-moment test frame and the results were compared with a previous study of females. Failures were conducted at approximately 90 N-m/s. In extension, the male upper cervical spine (O-C2) fails at a moment of 49.5 (s.d. 17.6)N-m and at an angle of 42.4 degrees (s.d. 8.0 degrees). In flexion, the mean moment at failure is 39.0 (s.d. 6.3 degrees) N-m and an angle of 58.7 degrees (s.d. 5.1 degrees). The difference in strength between flexion and extension is not statistically significant. The difference in the angles is statistically significant. The upper cervical spine was significantly stronger than the lower cervical spine in both flexion and extension. The male upper cervical spine was significantly stiffer than the female and significantly stronger than the female in flexion. Odontoid fractures were the most common injury produced in extension, suggesting a tensile mechanism due to tensile loads in the odontoid ligamentous complex.     

Morphological aspects and biomechanical properties of the vertebroaxial joint (C2-C3).

The C2-C3 intervertebral joint must be regarded as a transitional area situated between the upper cervical spine where most rotation of the neck and little flexion and extension occur and the lower cervical spaces where chiefly motion in the sagittal plane and also somewhat rotation take place. Under normal circumstances the range of flexion-extension reaches 11 degrees, slighter than below (19.5 degrees at C5-C6); on the opposite, the range of rotation attains 7 degrees; less than above but much more than below (0 degrees at C5-C6). The motion in the sagittal and coronal planes is relatively poor because of the location of Penning's motor-axis of C2 which runs far from the vertebral body and the lowness of the intervertebral disc. However, the rotation of C2 with respect to C3 is fair by the peculiar inclination of the articular facets which slope sagittally but also coronally and trace a sphere whereupon C2 may move around its motor-centre in any plane. When C2-C3 is surgically fixed by bone graft, the lack of motion is completed by a "compensatory movement" in the upper cervical spaces and especially at the atlantooccipital joint for flexion-extension. In the same way, C2-C3 may improve its mobility especially in the sagittal plane when the inferior partner is blocked by surgical or arthritic fusion.     

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.     

Realistic loading conditions for upper body bending

Different modes of load applications are used to simulate flexion and extension of the upper body. It is not clear which loading modes deliver realistic results and allow the comparison of different studies.In a numerical study, a validated finite element model of the lumbar spine, ranging from the vertebra L1 to the disc L5–S1 was employed. Each of six different loading modes was studied for simulating flexion and extension, including pure moments, an eccentric axial force, using a wedged fixture, and applying upper body weight plus follower load plus muscle forces. Intersegmental rotations, intradiscal pressures and facet joint contact forces were calculated. Where possible, results were compared to data measured in vivo.The results of the loading modes studied show a large variance for some values. Outcome measures such as flexion angle and intradiscal pressure differed at a segment by up to 44% and 88%, respectively, related to their maximum values. Intradiscal pressure is mainly determined by the magnitude of the applied compressive force. For flexion maximum contact forces between 0 and 69 N are predicted in each facet joint for different loading modes. For both flexion and extension, applying upper body weight plus follower load plus muscle forces as well as a follower load together with a bending moment delivers results which agreed well with in vivo data from the literature.Choosing an adequate loading mode is important in spine biomechanics when realistic results are required for intersegmental rotations, intradiscal pressure and facet joint contact forces. Only then will results of different studies be comparable.     

Spinal constraint modulates head instantaneous center of rotation and dictates head angular motion

The head is kinematically constrained to the torso through the spine and thus, the spine dictates the amount of output head angular motion expected from an input impact. Here, we investigate the spinal kinematic constraint by analyzing the head instantaneous center of rotation (HICOR) with respect to the torso in head/neck sagittal extension and coronal lateral flexion during mild loads applied to 10 subjects. We found the mean HICOR location was near the C5-C6 intervertebral joint in sagittal extension, and T2-T3 intervertebral joint in coronal lateral flexion. Using the impulse-momentum relationship normalized by subject mass and neck length, we developed a non-dimensional analytical ratio between output angular velocity and input linear impulse as a function of HICOR location. The ratio was 0.65 and 0.50 in sagittal extension and coronal lateral flexion respectively, implying 30% greater angular velocities in sagittal extension given an equivalent impulse. Scaling to subject physiology also predicts larger required impulses given greater subject mass and neck length to achieve equivalent angular velocities, which was observed experimentally. Furthermore, the HICOR has greater motion in sagittal extension than coronal lateral flexion, suggesting the head and spine can be represented with a single inverted pendulum in coronal lateral flexion, but requires a more complex representation in sagittal extension. The upper cervical spine has substantial compliance in sagittal extension, and may be responsible for the complex motion and greater extension angular velocities. In analyzing the HICOR, we can gain intuition regarding the neck’s role in dictating head motion during external loading.     

Effect of lumbar fasciae on the stability of the lower lumbar spine

The biomechanical effect of tensioning the lumbar fasciae (LF) on the stability of the spine during sagittal plane motion was analysed using a validated finite element model of the normal lumbosacral spine (L4-S1). To apply the tension in the LF along the direction of the fibres, a local coordinate was allocated using dummy rigid beam elements that originated from the spinous process. Up to 10 Nm of flexion and 7.5 Nm of extension moment was applied with and without 20 N of lateral tension in the LF. A follower load of 400 N was additionally applied along the curvature of the spine. To identify how the magnitude of LF tension related to the stability of the spine, the tensioning on the fasciae was increased up to 40 N with an interval of 10 N under 7.5 Nm of flexion/extension moment. A fascial tension of 20 N produced a 59% decrease in angular motion at 2.5 Nm of flexion moment while there was a 12.3% decrease at 10 Nm in the L5-S1 segment. Its decrement was 53 and 9.6% at 2.5 Nm and 10 Nm, respectively, in the L4-L5 segment. Anterior translation was reduced by 12.1 and 39.0% at the L4-L5 and L5-S1 segments under 10 Nm of flexion moment, respectively. The flexion stiffness shows an almost linear increment with the increase in fascial tension. The results of this study showed that the effect of the LF on the stability of the spine is significant.     

Development of a model based method for investigating facet articulation

Reported investigations of facet articulation in the human spine have often been conducted through the insertion of pressure sensitive film into the joint space, which requires incision of the facet capsule and may alter the characteristics of interaction between the facet surfaces. Load transmission through the facet has also been measured using strain gauges bonded to the articular processes. While this method allows for preservation of the facet capsule, it requires extensive instrumentation of the spine, as well as strain-gauge calibration, and is highly sensitive to placement and location of the strain gauges. The inherently invasive nature of these techniques makes it difficult to translate them into medical practice. A method has been developed to investigate facet articulation through the application of test kinematics to a specimen-specific rigid-body model of each vertebra within a lumbar spine segment. Rigid-body models of each vertebral body were developed from CT scans of each specimen. The distances between nearest-neighboring points on each facet surface were calculated for specific time frames of each specimen's flexion/extension test. A metric describing the proportion of each facet surface within a distance (2 mm) from the neighboring surface, the contact area ratio (CAR), was calculated at each of these time frames. A statistically significant difference (p<0.037) was found in the CAR between the time frames corresponding to full flexion and full extension in every level of the lumbar spine (L1-L5) using the data obtained from the seven specimens evaluated in this study. The finding that the contact area of the facet is greater in extension than flexion corresponds to other findings in the literature, as well as the generally accepted role of the facets in extension. Thus, a biomechanical method with a sufficiently sensitive metric is presented as a means to evaluate differences in facet articulation between intact and treated or between healthy and pathologic spines.     

Variations of stiffness and strength along the human cervical spine

The load-displacement response and strength of the mid (C2-C5) and lower (C5-T1) cervical regions were determined for combinations of sagittal loads, in vitro. In unpaired t-test comparisons, the mid cervical region was significantly stiffer in compression and extension than the lower region. In tests to failure, failure in six out of seven mid cervical specimens resulted from flexion alone, while combined compression-flexion was required to fail five of the eight lower cervical specimens. Post-test dissections revealed no regional differences in the pattern of failure. In addition to sagittal tests, the load-displacement responses of three-vertebrae cervical specimens were measured with the upper body axially rotated with respect to the lower body. The effect of this pre-torsion was to diminish the zone of low slope near zero load for axial, shear, and flexion motion. Three of the four axially rotated specimens failed in flexion without added compression. These controlled load-displacement measurements of cervical spine specimens describe for the first time the continuous flexion-compression response up to failure, and suggest that consideration of the biomechanics of three apparently distinct mobile regions of the cervical spine (C1-C2, C2-C5, C5-T1) may facilitate the interpretation of hazardous conditions and the diagnosis of injury. These data also provide basic information for the in vitro investigation of passive cervical spine protection such as helmets and head-rests, suggesting that the head should be kept in a non-rotated position to reduce risk of injury.     

Compressive force magnitude and intervertebral joint flexion/extension angle influence shear failure force magnitude in the porcine cervical spine

Despite the findings that peak anterior shear load is highly correlated with low-back pain reporting, very little research has been conducted to determine how vertebral shear injury potential is influenced. The current study quantified the combined effects of vertebral joint compression and flexion/extension postural deviation from neutral on ultimate shear failure. Ninety-six porcine cervical specimens (48C3-C4, 48C5-C6) were tested. Each specimen was randomly assigned to one of twelve combinations of compressive force (15%, 30%, 45%, or 60% of predicted compressive failure force) and flexion/extension postural deviation (extended, neutral, or flexed). Vertebral joint shear failure was induced by applying posterior shear displacement of the caudal vertebra at a constant rate of 0.15 mm/s. Throughout shear failure tests, vertebral joint kinematics were measured using an optoelectronic camera and a series of infrared light emitting diodes while shear force was measured from load cells rigidly interfaced in series with linear actuators that applied the shear displacement. Measurements of shear stiffness, ultimate force, displacement, and energy stored were made from the force-displacement data. Compressive force and postural deviation demonstrated main effects without a statistically significant interaction for any of the measurements. Shear failure force increased by 11.1% for each 15% increment in compressive force (p<0.05). Postural deviation from neutral impacted ultimate shear failure force by a 12.8% increase with extension (p<0.05) and a 13.2% decrease with flexion (p<0.05). Displacement at ultimate failure was not significantly altered by either compressive force or postural deviation. These results demonstrate that shear failure force may be governed by changes in facet articulation, either by postural deviation or by reducing vertebral joint height through compression that alter the moment arm length between the center of facet contact pressure and the pars interarticularis location. However, objective evidence of this alteration currently does not exist. Both compression and flexion/extension postural deviation should be equally considered while assessing shear injury potential.     

Load-bearing role of facets in a lumbar segment under sagittal plane loadings

In the present work, the load-bearing role of the facet joints in a lumbar I2-3 segment is quantitatively determined by means of a three dimensional nonlinear finite element program. The analysis accounts for both material and geometric nonlinearities and treats the facet articulation as a nonlinear moving contact problem. The disc nucleus is considered as an inviscid incompressible fluid and the annulus as a composite of collagenous fibres embedded in a matrix of ground substance. The spinal ligaments are modelled as a collection of nonlinear axial elements. The loadings consist of axial compression and sagittal plane shears and bending moments, acting alone or combined. The results show that in pure compression, the external axial force is transmitted primarily by the intervertebral disc. The facet joints carry only a small percentage of the force. However, the facet joints carry large forces in extension, whereas in small flexion they carry none. Addition of compression tends to increase these contact forces in extension while it has no effect on them in flexion. In extension, the forces on the facet joints are transmitted by both the articular surfaces and the capsular ligaments. Although in small flexion the facets carry no load, large contact forces are predicted to develop as the segment is flexed beyond 7-8 degrees. These forces are of the same magnitude as those computed under large extension rotation and are oriented nearly in the horizontal plane with negligible component in the axial direction. The horizontal components of the contact forces generated during articulation are often larger than the axial components which directly resist the applied compressive force. The axial components of the contact forces, therefore, grossly underestimate the total forces acting on the facets. The transfer of forces from one facet to the adjacent one occurs through distinct areas in flexion and in extension loadings. That is, on the superior articular surface, the contact area shifts from the upper tip in large flexion to the lower margin in extension. On the inferior articular surface, the contact area shifts from the upper and central regions in large flexion to the lower tip in extension.     

A technique for quantifying the bending moment acting on the lumbar spine in vivo

The bending properties of cadaveric lumbar spines were measured and used to convert in vivo measurements of lumbar flexion into bending moments ('stresses'). Forty-two lumbar motion segments were subjected to complex physiological loading and graphs were obtained of bending moment vs flexion angle. Variability was reduced by expressing both variables as a percentage of their values at the elastic limit. Data were averaged for each lumbar level, and a composite bending curve was compiled for the lumbar spine, L1-S1. A linear relationship was established between lumbar flexion measured in vitro and in vivo. This enabled values of 'per cent lumbar flexion' measured in vivo to be converted into 'per cent maximum bending moment' with a maximum likely error of about +/- 8%, which is equivalent to about +/- 5 Nm at L5-S1 for an average person. The technique was applied to 28 subjects, using dynamic measurements of lumbar flexion obtained with the '3-Space Isotrack' system. The bending moment at L5-S1 was 12 Nm on average when picking a pen up off the floor. Highly significant increases in bending moment were observed when heavier and bulkier objects were lifted.     

Experimental flexion/extension data corridors for validation of finite element models of the young, normal cervical spine

Finite element (FE) modeling is an important tool for studying the cervical spine in normal, injured and diseased conditions. To understand the role of mechanical changes on the spine as it goes from a normal to a diseased or injured state, experimental studies are needed to establish the external response of young, normal cervical spinal segments compared to injured or degenerated cervical spinal segments under physiologic loading. It is important to differentiate injured or degenerated specimens from young, normal specimens to provide accurate experimental results necessary for the validation of FE models. This study used seven young, normal fresh adult cadaver cervical spine segments C2-T1 ranging in age from 20 to 51 years. Prior to testing, the spines were graded in three ways: specimen quality, facet degeneration and disc degeneration. Spine segments were tested in flexion/extension, and the range of loads applied to the specimens was 0.33, 0.5, 1.0, 1.5 and 2.0 Nm. These loads resulted in rotations in the direction of loading as the primary response to loading. In general, results for young, normal specimens showed greater flexibility in flexion and less flexibility in extension than results previously reported in the literature. The flexion/extension curves are asymmetric with a greater magnitude in flexion than in extension. These experimental results will be used to validate FE models of young, normal cervical spines.     

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.