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)%).     

Effect of a pedicle-screw-based motion preservation system on lumbar spine biomechanics: A probabilistic finite element study with subsequent sensitivity analysis

Pedicle-screw-based motion preservation systems are often used to support a slightly degenerated disc. Such implants are intended to reduce intradiscal pressure and facet joints forces, while having a minimal effect on the motion patterns.In a probabilistic finite element study with subsequent sensitivity analysis, the effects of 10 input parameters, such as elastic modulus and diameter of the elastic rod, distraction of the segment, level of bridged segments, etc. on the output parameters intervertebral rotations, intradiscal pressures, and facet joint forces were determined. A validated finite element model of the lumbar spine was employed. Probabilistic studies were performed for seven loading cases: upright standing, flexion, extension, left and right lateral bending and left and right axial rotation.The simulations show that intervertebral rotation angles, intradiscal pressures and facet joint forces are in most cases reduced by a motion preservation system. The influence on intradiscal pressure is small, except in extension. For many input parameter combinations, the values for intervertebral rotations and facet joint forces are very low, which indicates that the implant is too stiff in these cases. The output parameters are affected most by the following input parameters: loading case, elastic modulus and diameter of the elastic rod, distraction of the segment, and angular rigidity of the connection between screws and rod.The designated functions of a motion preservation system can best be achieved when the longitudinal rod has a low stiffness, and when the connection between rod and pedicle screws is rigid.     

Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4–L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.     

Shape and orientation of articular facets of cervical vertebrae (C3-C7) in dogs denoting axial rotational ability: an osteological study

Macerated cervical vertebrae of 139 adult and 39 juvenile dogs were assessed for morphological evidence of axial rotation, which is denoted by the shape and orientation of the articular facets. The potential effect of the presence of caudal costal foveae at C7 on facet geometry was also investigated. Four variations of caudal facet shape were evident, namely, plane, concave, convex and sigmoid. The earliest onset of change in shape from plane to curved was noticed in juvenile dogs who were older than 8 weeks. The incidence of curved facets denoting axial rotation ability was observed to increase on descending down the cervical spine. Concave caudal facets were more frequently present in males and in large dogs (P < 0.001), compared to Dachshunds and small breeds. The degree of concavity was not related to age but was associated with the transverse distance between the most medial aspects of the caudal and corresponding cranial facets (P < 0.05). In large breeds, the presence of concave caudal facets at C7 was inversely related to the presence of caudal costal foveae (P < 0.01). The effect of axial rotation is discussed in context with clinical literature suggesting that axial rotation might attribute to the development of the wobbler syndrome.     

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.     

Quantitative evaluation of facet deflection,stiffness, strain and failure load during simulated cervical spine trauma

Traumatic cervical facet dislocation (CFD) is often associated with devastating spinal cord injury. Facet fractures commonly occur during CFD, yet quantitative measures of facet deflection, strain, stiffness and failure load have not been reported. The aim of this study was to determine the mechanical response of the subaxial cervical facets when loaded in directions thought to be associated with traumatic bilateral CFD – anterior shear and flexion. Thirty-one functional spinal units (6 × C2/3, C3/4, C4/5, and C6/7, 7 × C5/6) were dissected from fourteen human cadaver cervical spines (mean donor age 69 years, range 48–92; eight male). Loading was applied to the inferior facets of the inferior vertebra to simulate the in vivo inter-facet loading experienced during supraphysiologic anterior shear and flexion motion. Specimens were subjected to three cycles of sub-failure loading (10–100 N, 1 mm/s) in each direction, before being failed in a randomly assigned direction (10 mm/s). Facet deflection, surface strains, stiffness, and failure load were measured. Linear mixed-effects models (α = 0.05; random effect of cadaver) accounted for variations in specimen geometry and bone density. Specimen-specific parameters were significantly associated with most outcome measures. Facet stiffness and failure load were significantly greater in the simulated flexion loading direction, and deflection and surface strains were higher in anterior shear at the non-destructive analysis point (47 N applied load). The sub-failure strains and stiffness responses differed between the upper and lower subaxial cervical regions. Failure occurred through the facet tip during anterior shear loading, while failure through the pedicles was most common in flexion.     

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 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.     

Tensile cervical facet capsule ligament mechanics: failure and subfailure responses in the rat

Clinical, epidemiological, and biomechanical studies suggest the involvement of the cervical facet joint in neck pain. Mechanical studies have suggested the facet capsular ligament to be at risk for subfailure tensile injury during whiplash kinematics of the neck. Ligament mechanical properties can be altered by subfailure injury and such loading can induce cellular damage. However, at present, there is no clear understanding of the physiologic context of subfailure facet capsular ligament injury and mechanical implications for whiplash-related pain. Therefore, this study aimed to define a relationship between mechanical properties at failure and a subfailure condition associated with pain for tension in the rat cervical facet capsular ligament. Tensile failure studies of the C6/C7 rat cervical facet capsular ligament were performed using a customized vertebral distraction device. Force and displacement at failure were measured and stiffness and energy to failure were calculated. Vertebral motions and ligament deformations were tracked and maximum principal strains and their directions were calculated. Mean tensile force at failure (2.96 +/- 0.69 N) was significantly greater (p < 0.005) than force at subfailure (1.17 +/- 0.48 N). Mean ligament stiffness to failure was 0.75 +/- 0.27 N/mm. Maximum principal strain at failure (41.3 +/- 20.0%) was significantly higher (p = 0.003) than the corresponding subfailure value (23.1 +/- 9.3%). This study determined that failure and a subfailure painful condition were significantly different in ligament mechanics and findings provide preliminary insight into the relationship between mechanics and pain physiology for this ligament. Together with existing studies, these findings offer additional considerations for defining mechanical thresholds for painful injuries.     

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.     

Effect of annular incision type on the change in biomechanical properties in a herniated lumbar intervertebral disc

The technique used to incise the disc during discectomy may play a role in the subsequent healing and change in biomechanical stiffness of the disc. Several techniques of lumbar disc annulotomy have been described in clinical reports. The purpose of this paper was to study the influence of annulotomy technique on motion segment stiffness using a finite element model. Four incision methods (square, circular, cross, and slit) were compared. The analyses showed that each of the annular incisions produced increase in motions under axial moment loadings with circular incision producing the largest change in the corresponding rotational motion. Under shear loading mode, cross and slit-type annular incisions produced slightly larger changes in the principal motions of the disc than square and circular incisions. All other incision types considered in the current study produced negligibly small increase in motion under rest of the loading conditions. In addition to annulotomy, when nucleotomy was also included in the analyses, once again cross and slit incisions produced larger change in motion under shear loading mode as compared to the other two incision types. A comparison between the four types of annular incisions showed that cross incision produced an increase in motion larger than those produced by the other three incisions under flexion/extension and lateral moment loading and both shear force loadings. Circular incision produced the largest increase in motion under axial moment load in comparison to those produced by square, cross, and slit incisions. Sagittal plane symmetry was influenced by the incision injury to the motion segment leading to coupled motions as well as increased facet loads. From the study it can be concluded that the increase inflexibility of the disc due to annulotomy depends on the type of annulotomy and the annulotomy also produce asymmetrical deformations leading to increased facet loading.     

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.     

In vitro analysis of kinematics and elastostatics of the human rib cage during thoracic spinal movement for the validation of numerical models

Neither kinematic nor stiffness properties of the rib cage during thoracic spinal motion were investigated in previous studies, while being essential for the accurate validation of numerical models of the whole thorax. The aim of this in vitro study therefore was to quantify the kinematics and elastostatics of the human rib cage under defined boundary conditions. Eight fresh frozen human thoracic spine specimens (C7-L1, median age 55 years, ranging from 40 to 60 years) including entire rib cages were loaded quasi-statically in flexion/extension, lateral bending, and axial rotation using pure moments of 5 Nm. Relative motions of ribs, thoracic vertebrae, and sternal structures as well as strains on the ribs were measured using optical motion tracking of 150 reflective markers per specimen, while specimens were loaded displacement-controlled with a constant rate of 1°/s for 3.5 cycles. The third full cycle was used to determine relative angles and strains at full loading of the spine for all motion directions. Largest relative angles were found in the main loading directions with only small motions at the mid-thoracic levels. Highest strains of the intercostal spaces were detected in the anterior section of the lowest fourth of the rib cage, showing compressions and elongations of more than 10% in all spinal motion planes. Elastostatic rib deformation was generally less than 1%. Rib-sternum relative motions exhibited complex motion patterns, overall showing relative angles below 2°. The results indicate that rib cage structures are not macroscopically deformed during spinal motion, but exhibit characteristic reproducible kinematics patterns.     

A videofluoroscopy-based tracking algorithm for quantifying the time course of human intervertebral displacements

The motions of individual intervertebral joints can affect spine motion, injury risk, deterioration, pain, treatment strategies, and clinical outcomes. Since standard kinematic methods do not provide precise time-course details about individual vertebrae and intervertebral motions, information that could be useful for scientific advancement and clinical assessment, we developed an iterative template matching algorithm to obtain this data from videofluoroscopy images. To assess the bias of our approach, vertebrae in an intact porcine spine were tracked and compared to the motions of high-contrast markers. To estimate precision under clinical conditions, motions of three human cervical spines were tracked independently ten times and vertebral and intervertebral motions associated with individual trials were compared to corresponding averages. Both tests produced errors in intervertebral angular and shear displacements no greater than 0.4° and 0.055 mm, respectively. When applied to two patient cases, aberrant intervertebral motions in the cervical spine were typically found to correlate with patient-specific anatomical features such as disc height loss and osteophytes. The case studies suggest that intervertebral kinematic time-course data could have value in clinical assessments, lead to broader understanding of how specific anatomical features influence joint motions, and in due course inform clinical treatments.     

Comparing non-invasive scapular tracking methods across elevation angles,planes of elevation and humeral axial rotations

Altered scapular motions premeditate shoulder impingement and other musculoskeletal disorders. Divergent experimental conditions in previous research precludes rigorous comparisons of non-invasive scapular tracking techniques. This study evaluated scapular orientation measurement methods across an expanded range of humeral postures. Scapular medial/lateral rotation, anterior/posterior tilt and protraction/retraction was measured using an acromion marker cluster (AMC), a scapular locator, and a reference stylus. Motion was captured using reflective markers on the upper body, as well as on the AMC, locator and stylus. A combination of 5 arm elevation angles, 3 arm elevation planes and 3 arm axial rotations was examined. Measurement method interacted with elevation angle and plane of elevation for all three scapular orientation directions (p < 0.01). Method of measurement interacted with axial rotation in anterior/posterior tilt and protraction/retraction (p < 0.01). The AMC had strong agreement with the reference stylus than the locator for the majority of humeral elevations, planes and axial rotations. The AMC underestimated lateral rotation, with the largest difference of ∼2° at 0° elevation. Both the locator and AMC overestimated posterior tilt at high arm elevation by up to 7.4°. Misestimations from using the locator could be enough to potentially obscure meaningful differences in scapular rotations.     

Uncertainty analysis of material properties and morphology parameters in numerical models regarding the motion of lumbar vertebral segments

Abstract

The kinematics of a spinal motion segment is determined by the material properties of the soft-tissue and the morphology. The material properties can vary within subjects and between vertebral levels, leading to a wide possible range of motion of a spinal segment independently on its morphology. The goal of this numerical study was to identify the most influential material parameters concerning the kinematics of a spinal motion segment and their plausible ranges. Then, a method was tested to deduce the material properties automatically, based on a given ROM and morphology. A fully parametric finite element model of the morphology and material properties of a lumbar spinal motion segment was developed. The impact of uncertainty of twelve spinal material parameters, as well as the size of the gap between the articular surfaces of the facet joints was examined. The simulation results were compared to our own in vitro data. The flexibility of a lumbar segment was especially influenced by the properties of the anterior annulus region, the facet gap size and the interspinous ligament. The high degree of uncertainty in the material properties and facet gap size published in the literature can lead to a wide scatter in the motion of a spinal segment, with a range of 6°-17° in the intact condition in flexion/extension, from 5°-22° in lateral bending and from 3°-14° in axial rotation. Statistical analysis of the variability might help to estimate the sensitivity and total uncertainty propagated through biomechanical simulations, affecting the reliability of the predictions.     

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.     

Measured and modelled static and dynamic axial trunk torsion during twisting in males and females

Study of the mechanics of trunk twisting is special interest given epidemiological evidence linking occupational twisting to increased incidence of low back pain. An anatomically detailed, three-dimensional model of the trunk (rib cage, pelvis, five lumbar vertebrae and 50 muscles), was used to predict maximum axial trunk torque. Predicted axial torques were compared with measured torques. Thirty-one (10 male and 21 female) subjects performed maximum effort isometric twisting exertions, at 0° of twist and ±30° of twist together with dynamic exertions, at 30° s⁻¹ and 60° s⁻¹. Females were able to generate approximately two-thirds of the torque of males (males, 97 Nm; females 60 Nm, isometric at 0°). When the trunk was prerotated to 30°, subjects were able to generate greater torque when the effort was toward the 0° position (approximately 105 Nm by males and 68 Nm by females). Experimental data indicated that velocity of rotation and amount of twist are important modulators of axial torque. Changes in muscle length were demonstrated to be minimal from model output as most muscle length changes during a twist from 0° to 30°, measured between the pelvis and the shoulder harness, were less than 1%, although some portions of the abdominal obliques underwent a length excursion of 5%. The small changes in the individual muscle force components that contribute to twist, i.e. the muscle unit vector about the axial twist axis and its moment arm that change as a function of twisted position, do not entirely account for the measured differences in torque, suggesting that additional mechanisms influence axial torque generation.