Effect of loading rate and hydration on the mechanical properties of the disc

An Extended Neck Position is Likely to Produce Cervical Spine Injuries Through Buckling in Accidental Head-First Impacts During Rugby Tackling

Catastrophic cervical spine injuries in rugby often occur during tackling. The underlying mechanisms leading to these injuries remain unclear, with neck hyperflexion and buckling both proposed as the causative factor in the injury prevention literature. The aim of this study was to evaluate the effect of pre-impact head-neck posture on intervertebral neck loads and motions during a head-first rugby tackle. Using a validated, subject-specific musculoskeletal model of a rugby player, and computer simulations driven by in vivo and in vitro data, we examined the dynamic response of the cervical spine under such impact conditions. The simulations demonstrated that the initial head-neck sagittal-plane posture affected intervertebral loads and kinematics, with an extended neck resulting in buckling and supraphysiologic intervertebral shear and flexion loads and motions, typical of bilateral facet dislocation injuries. In contrast, an initially flexed neck increased axial compression forces and flexion angles without exceeding intervertebral physiological limits. These findings provide objective evidence that can inform injury prevention strategies or rugby law changes to improve the safety of the game of rugby.

Geometric determinants of the mechanical behavior of image-based finite element models of the intervertebral disc

The intervertebral disc is an important structure for load transfer through the spine. Its injury and degeneration have been linked to pain and spinal fractures. Disc injury and spine fractures are associated with high stresses; however, these stresses cannot be measured, necessitating the use of finite element (FE) models. These models should include the disc's complex structure, as changes in disc geometry have been linked to altered mechanical behavior. However, image-based models using disc-specific structures have yet to be established. This study describes a multiphasic FE modeling approach for noninvasive estimates of subject-specific intervertebral disc mechanical behavior based on medical imaging. The models (n = 22) were used to study the influence of disc geometry on the predicted global mechanical response (moments and forces), internal local disc stresses, and tractions at the interface between the disc and the bone. Disc geometry was found to have a strong influence on the predicted moments and forces on the disc (R2 = 0.69-0.93), while assumptions regarding the side curvature (bulge) of the disc had only a minor effect. Strong variability in the predicted internal disc stresses and tractions was observed between the models (mean absolute differences of 5.1%-27.7%). Disc height had a systematic influence on the internal disc stresses and tractions at the disc-to-bone interface. The influence of disc geometry on mechanics highlights the importance of disc-specific modeling to estimate disc injury risk, loading on the adjacent vertebral bodies, and the mechanical environment present in disc tissues.

Effect of muscle activation on dynamic responses of neck of pilot during emergency ejection: a finite element study

To determine the effect of muscle activation on the dynamic responses of the neck of a pilot during simulated emergency ejections. A complete finite element model of the pilot's head and neck was developed and dynamically validated. Three muscle activation curves were designed to simulate different activation times and levels of muscles during pilot ejection: A is the unconscious activation curve of the neck muscles, B is the pre-activation curve, and C is the continuous activation curve. The acceleration-time curves obtained during ejection were applied to the model, and the influence of the muscles on the dynamic responses of the neck was investigated by analyzing both angles of rotation of the neck segments and disc stresses. Muscle pre-activation reduced fluctuations in the angle of rotation in each phase of the neck. Continuous muscle activation caused a 20% increase in the angle of rotation compared to pre-activation. Moreover, it resulted in a 35% increase in the load on the intervertebral disc. The maximum stress on the disc occurred in the C4-C5 phase. Continuous muscle activation increased both the axial load on the neck and the posterior extension angle of rotation of the neck. Muscle pre-activation during emergency ejection has a protective effect on the neck. However, continuous muscle activation increases the axial load and rotation angle of the neck. A complete finite element model of the pilot's head and neck was established and three neck muscle activation curves were designed to investigate the effects of muscle activation time and level on the dynamic response of the pilot's neck during ejection. This increased insights into the protection mechanism of neck muscles on the axial impact injury of the pilot's head and neck.

Importance of Matrix Cues on Intervertebral Disc Development, Degeneration, and Regeneration

Back pain is one of the leading causes of disability worldwide and is frequently caused by degeneration of the intervertebral discs. The discs’ development, homeostasis, and degeneration are driven by a complex series of biochemical and physical extracellular matrix cues produced by and transmitted to native cells. Thus, understanding the roles of different cues is essential for designing effective cellular and regenerative therapies. Omics technologies have helped identify many new matrix cues; however, comparatively few matrix molecules have thus far been incorporated into tissue engineered models. These include collagen type I and type II, laminins, glycosaminoglycans, and their biomimetic analogues. Modern biofabrication techniques, such as 3D bioprinting, are also enabling the spatial patterning of matrix molecules and growth factors to direct regional effects. These techniques should now be applied to biochemically, physically, and structurally relevant disc models incorporating disc and stem cells to investigate the drivers of healthy cell phenotype and differentiation. Such research will inform the development of efficacious regenerative therapies and improved clinical outcomes.

In Vitro Studies for Investigating Creep of Intervertebral Discs under Axial Compression: A Review of Testing Environment and Results

Creep responses of intervertebral discs (IVDs) are essential for spinal biomechanics clarification. Yet, there still lacks a well-recognized investigation protocol for this phenomenon. Current work aims at providing researchers with an overview of the in vitro creep tests reported by previous studies, specifically specimen species, testing environment, loading regimes and major results, based on which a preliminary consensus that may guide future creep studies is proposed. Specimens used in creep studies can be simplified as a “bone–disc–bone” structure where three mathematical models can be adopted for describing IVDs’ responses. The preload of 10–50 N for 30 min or three cycles followed by 4 h-creep under constant compression is recommended for ex vivo simulation of physiological condition of long-time sitting or lying. It is worth noticing that species of specimens, environment temperature and humidity all have influences on biomechanical behaviors, and thus are summarized and compared through the literature review. All factors should be carefully set according to a guideline before tests are conducted to urge comparable results across studies. To this end, this review also provides a guideline, as mentioned before, and specific steps that might facilitate the community of biomechanics to obtain more repeatable and comparable results from both natural specimens and novel biomaterials.

An Artificial PVA-BC Composite That Mimics the Biomechanical Properties and Structure of a Natural Intervertebral Disc

Disc herniation is one of the most ubiquitous healthcare problems in modern cities—severe patients eventually require surgical intervention. However, the existing operations—spinal fusion and artificial disc replacement—alter the biomechanics of the spine, leaving much room for improvement. The appropriateness of polyvinyl alcohol (PVA) for biomedical applications has been recognised due to its high water content, excellent biocompatibility, and versatile mechanical properties. In this study, a newly-designed PVA–bacterial cellulose (PVA-BC) composite was assembled to mimic both the biomechanics and annular structure of natural intervertebral discs (IVDs). PVA-BC composites of various concentrations were fabricated and tested under unconfined compression and compressive creep in order to acquire the values of the normalised compressive stiffness and whole normalised deformation. The normalised compressive stiffness increased considerably with an increasing PVA concentration, spanning from 1.82 (±0.18) to 3.50 (±0.14) MPa, and the whole normalised deformation decreased from 0.25 to 0.13. Formulations of 40% PVA provided the most accurate mimicry of natural human IVDs in normalised whole deformation, and demonstrated higher dimensional stability. The biocompatible results further confirmed that the materials had excellent biocompatibility. The novel bionic structure and formulations of the PVA-BC materials mimicked the biomechanics and structure of natural IVDs, and ensured dimensional stability under prolonged compression, reducing the risk of impingement on the surrounding tissue. The PVA-BC composite is a promising material for third-generation artificial IVDs with integrated construction.

Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.

Multiscale Validation of Multiple Human Body Model Functional Spinal Units

A validation comparing five human body model (HBM) lumbar spines is carried out across two load cases, with the objective to use and apply HBMs in high strain rate applications such as car occupant simulation. The first load case consists of an individual intervertebral disc (IVD) loaded in compression at a strain rate of 1/s by a material testing machine. The second load case is a lumbar functional spine unit (FSU) loaded in compression using a drop tower setup, producing strain rates of up to 48/s. The IVD simulations were found to have a better agreement with the experiments than the FSU simulations, and the ranking of which HBMs matched best to the experiment differed by load case. These observations suggest the need for more hierarchical validations of the lumbar spine for increasing the utility of HBMs in high strain rate loading scenarios.

Measurement of Three-Dimensional Internal Dynamic Strains in the Intervertebral Disc of the Lumbar Spine With Mechanical Loading and Golden-Angle Radial Sparse Parallel-Magnetic Resonance Imaging

BACKGROUND

Noninvasive measurement of internal dynamic strain can be potentially useful to characterize spine intervertebral disc (IVD) in the setting of injury or degenerative disease.

PURPOSE

To develop and demonstrate a noninvasive technique to quantify three-dimensional (3D) internal dynamic strains in the IVD using a combination of static mechanical loading of the IVD using a magnetic resonance imaging (MRI)-compatible ergometer.

STUDY TYPE

Prospective.

SUBJECTS

Silicone gel phantom studies were conducted to assess strain variation with load and repeatability. Mechanical testing was done on the phantoms to confirm MR results. Eight healthy human volunteers (four men and four woman, age = 29 ± 5 years) underwent MRI using a rest, static loading, and recovery paradigm. Repeatability tests were conducted in three subjects.

FIELD STRENGTH/SEQUENCE

MRI (3 T) with 3D continuous golden-angle radial sparse parallel (GRASP) and compressed sensing (CS) reconstruction.

ASSESSMENT

CS reconstruction of the images, motion deformation, and Lagrangian strain maps were calculated for five IVD segments from L1/L2 to L5/S1.

STATISTICAL TESTS

Ranges of displacement and strain in each subject and the resulting mean and standard deviation were calculated. Student t-tests were used to calculate changes in strain from loading to recovery. The correlation coefficient (CC) in the repeatability study was calculated.

RESULTS

The most compressive strain experienced by the IVD segments under loaded conditions was in the L4/L5 segment (-7.5 ± 2.9%). The change in minimum strain from load to recovery was the most for the L4/L5 segment (-7.5% to -5.0%, P = 0.026) and the least for the L1/L2 segment (-4.4% to -3.9%, P = 0.51). In vivo repeatability in three subjects shows strong correlation between scans in subjects done 6 months apart, with CCs equal to 0.86, 0.94, and 0.94 along principal directions.

DATA CONCLUSION

This study shows the feasibility of using static mechanical loading with continuous GRASP-MRI acquisition with CS reconstruction to measure 3D internal dynamic strains in the spine IVD.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

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

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

In vivo intervertebral disc deformation: intratissue strain patterns within adjacent discs during flexion-extension

The biomechanical function of the intervertebral disc (IVD) is a critical indicator of tissue health and pathology. The mechanical responses (displacements, strain) of the IVD to physiologic movement can be spatially complex and depend on tissue architecture, consisting of distinct compositional regions and integrity; however, IVD biomechanics are predominately uncharacterized in vivo. Here, we measured voxel-level displacement and strain patterns in adjacent IVDs in vivo by coupling magnetic resonance imaging (MRI) with cyclic motion of the cervical spine. Across adjacent disc segments, cervical flexion-extension of 10° resulted in first principal and maximum shear strains approaching 10%. Intratissue spatial analysis of the cervical IVDs, not possible with conventional techniques, revealed elevated maximum shear strains located in the posterior disc (nucleus pulposus) regions. IVD structure, based on relaxometric patterns of T₂ and T_1ρ images, did not correlate spatially with functional metrics of strain. Our approach enables a comprehensive IVD biomechanical analysis of voxel-level, intratissue strain patterns in adjacent discs in vivo, which are largely independent of MRI relaxometry. The spatial mapping of IVD biomechanics in vivo provides a functional assessment of adjacent IVDs in subjects, and provides foundational biomarkers for elastography, differentiation of disease state, and evaluation of treatment efficacy.

The Effect of the Loading Rate on the Full-Field Strain Distribution on the Surface on the Intervertebral Discs

Contrasting results are reported when the spine is tested at different strain rates. Tissue specimens from the ligaments or the intervertebral discs (IVD, including annulus fibrosus and nucleus pulposus) exhibit higher stiffness and lower dissipation at high strain rates. Counterintuitively, when spine segments are tested at high rates, the hysteresis area and loop width increase. It is unclear how the load is shared between the different structures at different loading rates. The hypotheses of this study were: (i) As the IVD stiffens at higher loading rates, the strain distribution around the disc would be different depending on the loading rate; (ii) Preconditioning attenuates the strain-rate dependency of the IVD, thus making differences in strain distribution smaller at the different rates. Six segments of three vertebrae (L4-L6) were extracted from porcine spines and tested in presso-flexion at different loading rates (reaching full load in 0.67, 6.7, and 67 s). The full-field strain maps were measured using digital image correlation on the surface of the IVDs from lateral. The posterior-to-anterior trends of the strain were computed in detail for each IVD, and compared between loading rates. The values and the direction of principal strain on the surface of the IVDs, vertebrae, and endplates remained unchanged at different rates. In the transition zone between IVD and vertebra, only slight differences due to the loading rate appeared but with no statistical significance. These findings will allow better understanding of the rate-dependent behavior and failure of the IVD.

Effect of Strain Rates on Failure of Mechanical Properties of Lumbar Intervertebral Disc Under Flexion

Objective

To evaluate the strain‐rate‐dependent viscoelastic properties of the intervertebral disc by in vitro experiments.

Method

The biomechanical experiments were conducted from September 2019 to December 2019. The lumbar spines of sheep were purchased within 4–6 hours from the local slaughterhouse, and the intervertebral disc samples were divided into three groups. In rupture group, the samples were used to test the mechanical behavior of the intervertebral disc rupture at different strain rates. In fatigue injury group, the samples were used to test the mechanical behavior of fatigue injury on the intervertebral disc under different strain rates. In internal displacement group, the samples were used to test the internal displacement distribution of the intervertebral disc at different strain rates by applying an optimized digital image correlation (DIC) technique.

Results

Both the yielding and cracking phenomenon occurs at fast and medium loading rates, while only the yielding phenomenon occurs at a slow loading rate. The yield stress, compressive strength, and elastic modulus all increase with the increase of the strain rate, while the yield strain decreases with the increase of the strain rate. The logarithm of the elastic modulus in the intervertebral disc is approximately linear with the logarithm of the strain rate under different strain rates. Both before and after fatigue loading, the stiffness in the loading and unloading curves of the intervertebral disc is inconsistent, forming a hysteresis loop, which is caused by the viscoelastic effect. The strain rate has no significant effect on the internal displacement distribution of the intervertebral disc. Based on the experimental data, the constitutive relationship of the intervertebral disc at different strain rates is obtained. The fitting curves are well coupled with the experimental data, while the fitting parameters are approximately linear with the logarithm of the strain rate.

Conclusions

These experiments indicate that the strain rate has a significant effect on the mechanical behavior of the intervertebral disc rupture and fatigue injury, while the constitutive equation can predict the rate‐dependent mechanical behavior of lumbar intervertebral disc under flexion very well. These results have important theoretical guiding significance for preventing lumbar disc herniation in daily life.

Occlusion of the lumbar spine canal during high-rate axial compression

BACKGROUND CONTEXT

While burst fracture is a well-known cause of spinal canal occlusion with dynamic, axial spinal compression, it is unclear how such loading mechanisms might cause occlusion without fracture.

PURPOSE

To determine how spinal canal occlusion during dynamic compression of the lumbar spine is differentially caused by fracture or mechanisms without fracture and to examine the influence of spinal level on occlusion.

STUDY DESIGN

A cadaveric biomechanical study.

METHODS

Twenty sets of three-vertebrae specimens from all spinal levels between T12 and S1 were subjected to dynamic compression using a hydraulic loading apparatus up to a peak velocity between 0.1 and 0.9 m/s. The presence of canal occlusion was measured optically with a high-speed camera. This was repeated with incremental increases of 4% compressive strain until a vertebral fracture was detected using acoustic emission measurements and computed tomographic imaging.

RESULTS

For axial compression without fracture, the peak occlusion (O_max) was 29.9±10.0%, which was deduced to be the result of posterior bulging of the intervertebral disc into the spinal canal. O_max correlated significantly with lumbar spinal level (p<.001), the compressive displacement (p<.001) and the cross-sectional area of the vertebra (p=.031).

CONCLUSIONS

Spinal canal occlusion observed without vertebral fracture involves intervertebral disc bulging. The lower lumbar spine tended to be more severely occluded than more proximal levels.

CLINICAL SIGNIFICANCE

Clinically, intermittent canal occlusion from disc bulging during dynamic compression may not show any radiographic features. The lower lumbar spine should be a focus of injury prevention intervention in cases of high-rate axial compression.

Influence of testing environment and loading rate on intervertebral disc compressive mechanics: An assessment of repeatability at three different laboratories

In vitro mechanical testing of intervertebral discs is crucial for basic science and pre-clinical testing. Generally, these tests aim to replicate in vivo conditions, but simplifications are necessary in specimen preparation and mechanical testing due to complexities in both structure and the loading conditions required to replicate in vivo conditions. There has been a growing interest in developing a consensus of testing protocols within the spine community to improve comparison of results between studies. The objective of this study was to perform axial compression experiments on bovine bone-disc-bone specimens at three institutions. No differences were observed between testing environment being air, with PBS soaked gauze, or a PBS bath (P > .206). A 100-fold increase in loading rate resulted in a small (2%) but significant increase in compressive mechanics (P < .017). A 7% difference in compressive stiffness between Labs B and C was eliminated when values were adjusted for test system compliance. Specimens tested at Lab A, however, were found to be stiffer than specimens from Lab B and C. Even after normalizing for disc geometry and adjusting for system compliance, an ∼35% difference was observed between UK based labs (B and C) and the USA based lab (A). Large differences in specimen stiffness may be due to genetic differences between breeds or in agricultural feed and use of growth hormones; highlighting significant challenges in comparing mechanics data across studies. This research provides a standardized test protocol for the comparison of spinal specimens and provides steps towards understanding how location and test set-up may affect biomechanical results.

Effect of Changing Heart Rate on the Ocular Pulse and Dynamic Biomechanical Behavior of the Optic Nerve Head

Purpose

To study the effect of changing heart rate on the ocular pulse and the dynamic biomechanical behavior of the optic nerve head (ONH) using a comprehensive mathematical model.

Methods

In a finite element model of a healthy eye, a biphasic choroid consisted of a solid phase with connective tissues and a fluid phase with blood, and the lamina cribrosa (LC) was viscoelastic as characterized by a stress-relaxation test. We applied arterial pressures at 18 ocular entry sites (posterior ciliary arteries), and venous pressures at four exit sites (vortex veins). In the model, the heart rate was varied from 60 to 120 bpm (increment: 20 bpm). We assessed the ocular pulse amplitude (OPA), pulse volume, ONH deformations, and the dynamic modulus of the LC at different heart rates.

Results

With an increasing heart rate, the OPA decreased by 0.04 mm Hg for every 10 bpm increase in heart rate. The ocular pulse volume decreased linearly by 0.13 µL for every 10 bpm increase in heart rate. The storage modulus and the loss modulus of the LC increased by 0.014 and 0.04 MPa, respectively, for every 10 bpm increase in heart rate.

Conclusions

In our model, the OPA, pulse volume, and ONH deformations decreased with an increasing heart rate, whereas the LC became stiffer. The effects of blood pressure/heart rate changes on ONH stiffening may be of interest for glaucoma pathology.

Osmo-inelastic response of the intervertebral disc annulus fibrosus tissue

The aim of this article is to provide some insights on the osmo-inelastic response under stretching of annulus fibrosus of the intervertebral disc. Circumferentially oriented specimens of square cross section, extracted from different regions of bovine cervical discs (ventral-lateral and dorsal-lateral), are tested under different strain-rates and saline concentrations within normal range of strains. An accurate optical strain measuring technique, based upon digital image correlation, is used in order to determine the full-field displacements in the lamellae and fibers planes of the layered soft tissue. Annulus stress-stretch relationships are measured along with full-field transversal strains in the two planes. The mechanical response is found hysteretic, rate-dependent and osmolarity-dependent with a Poisson's ratio higher than 0.5 in the fibers plane and negative (auxeticity) in the lamellae plane. While the stiffness presents a regional-dependency due to variations in collagen fibers content/orientation, the strain-rate sensitivity of the response is found independent on the region. A significant osmotic effect is found on both the auxetic response in the lamellae plane and the stiffness rate-sensitivity. These local experimental observations will result in more accurate chemo-mechanical modeling of the disc annulus and a clearer multi-scale understanding of the disc intervertebral function.

Finite element models can reproduce the effect of nucleotomy on the multi-axial compliance of human intervertebral discs

Finite element (FE) models can unravel the link between intervertebral disc (IVD) degeneration and its mechanical behaviour. Nucleotomy may provide the data required for model verification. Three human IVDs were scanned with MRI and tested in multiple loading scenarios, prior and post nucleotomy. The resulting data was used to generate, calibrate, and verify the FE models. Nucleotomy increased the experimental range of motion by 26%, a result reproduced by the FE simulation within a 5% error. This work demonstrates the ability of FE models to reproduce the mechanical compliance of human IVDs prior and post nucleotomy.

How pre-strain affects the chemo-torsional response of the intervertebral disc

BACKGROUND

The role of the axial pre-strain on the torsional response of the intervertebral disc remains largely undefined. Moreover, the chemo-mechanical interactions in disc tissues are still unclear and corresponding data are rare in the literature. The paper deals with an in-vitro study of the pre-strain effect on the chemical sensitivity of the disc torsional response.

METHODS

Fifteen non-frozen 'motion segments' (two vertebrae and the intervening soft tissues) were extracted from the cervical spines of mature sheep. The motion segments were loaded in torsion at various saline concentrations and axial pre-strain levels in order to modulate the intradiscal pressure. After preconditioning with successive low-strain compressions at a magnitude of 0.1 mm (10 cycles at 0.05 mm/s), the motion segment was subjected to a cyclic torsion until a twisting level of 2 deg. at 0.05 deg./s while a constant axial pre-strain (in compression or in tension) is maintained, the saline concentration of the surrounding fluid bath being changed from hypo-osmotic condition to hyper-osmotic condition.

FINDINGS

Analysis of variance shows that the saline concentration influences the torsional response only when the motion segments are pre-compressed (p < .001) with significant differences between hypo-osmotic condition and hyper-osmotic condition.

INTERPRETATION

The combination of a compressive pre-strain with twisting amplifies the nucleus hydrostatic pressure on the annulus and the annulus collagen fibers tensions. The proteoglycans density increases with the compressive pre-strain and leads to higher chemical imbalances, which would explain the increase in chemical sensitivity of the disc torsional response.

The strain distribution in the lumbar anterior longitudinal ligament is affected by the loading condition and bony features: An in vitro full-field analysis

The role of the ligaments is fundamental in determining the spine biomechanics in physiological and pathological conditions. The anterior longitudinal ligament (ALL) is fundamental in constraining motions especially in the sagittal plane. The ALL also confines the intervertebral discs, preventing herniation. The specific contribution of the ALL has indirectly been investigated in the past as a part of whole spine segments where the structural flexibility was measured. The mechanical properties of isolated ALL have been measured as well. The strain distribution in the ALL has never been measured under pseudo-physiological conditions, as part of multi-vertebra spine segments. This would help elucidate the biomechanical function of the ALL. The aim of this study was to investigate in depth the biomechanical function of the ALL in front of the lumbar vertebrae and of the intervertebral disc. Five lumbar cadaveric spine specimens were subjected to different loading scenarios (flexion-extension, lateral bending, axial torsion) using a state-of-the-art spine tester. The full-field strain distribution on the anterior surface was measured using digital image correlation (DIC) adapted and validated for application to spine segments. The measured strain maps were highly inhomogeneous: the ALL was generally more strained in front of the discs than in front of the vertebrae, with some locally higher strains both imputable to ligament fibers and related to local bony defects. The strain distributions were significantly different among the loading configurations, but also between opposite directions of loading (flexion vs. extension, right vs. left lateral bending, clockwise vs. counterclockwise torsion). This study allowed for the first time to assess the biomechanical behaviour of the anterior longitudinal ligament for the different loading of the spine. We were able to identify both the average trends, and the local effects related to osteophytes, a key feature indicative of spine degeneration.

Material properties of human lumbar intervertebral discs across strain rates

BACKGROUND CONTEXT

The use of finite element (FE) methods to study the biomechanics of the intervertebral disc (IVD) has increased over recent decades due to their ability to quantify internal stresses and strains throughout the tissue. Their accuracy is dependent upon realistic, strain-rate dependent material properties, which are challenging to acquire.

PURPOSE

The aim of this study was to use the inverse FE technique to characterize the material properties of human lumbar IVDs across strain rates.

STUDY DESIGN

A human cadaveric experimental study coupled with an inverse finite element study.

METHODS

To predict the structural response of the IVD accurately, the material response of the constituent structures was required. Therefore, compressive experiments were conducted on 16 lumbar IVDs (39±19 years) to obtain the structural response. An FE model of each of these experiments was developed and then run through an inverse FE algorithm to obtain subject-specific constituent material properties, such that the structural response was accurate.

RESULTS

Experimentally, a log-linear relationship between IVD stiffness and strain rate was observed. The material properties obtained through the subject-specific inverse FE optimization of the annulus fibrosus (AF) fiber and AF fiber ground matrix allowed a good match between the experimental and FE response. This resulted in a Young modulus of AF fibers (-MPa) to strain rate (ε˙, /s) relationship of YMAF=31.5ln(ε˙)+435.5, and the C₁₀ parameter of the Neo-Hookean material model of the AF ground matrix was found to be strain-rate independent with an average value of 0.68 MPa.

CONCLUSIONS

These material properties can be used to improve the accuracy, and therefore predictive ability of FE models of the spine that are used in a wide range of research areas and clinical applications.

CLINICAL SIGNIFICANCE

Finite element models can be used for many applications including investigating low back pain, spinal deformities, injury biomechanics, implant design, design of protective systems, and degenerative disc disease. The accurate material properties obtained in this study will improve the predictive ability, and therefore clinical significance of these models.

Synchrotron tomography of intervertebral disc deformation quantified by digital volume correlation reveals microstructural influence on strain patterns

The intervertebral disc (IVD) has a complex and multiscale extracellular matrix structure which provides unique mechanical properties to withstand physiological loading. Low back pain has been linked to degeneration of the disc but reparative treatments are not currently available. Characterising the disc's 3D microstructure and its response in a physiologically relevant loading environment is required to improve understanding of degeneration and to develop new reparative treatments. In this study, techniques for imaging the native IVD, measuring internal deformation and mapping volumetric strain were applied to an in situ compressed ex vivo rat lumbar spine segment. Synchrotron X-ray micro-tomography (synchrotron CT) was used to resolve IVD structures at microscale resolution. These image data enabled 3D quantification of collagen bundle orientation and measurement of local displacement in the annulus fibrosus between sequential scans using digital volume correlation (DVC). The volumetric strain mapped from synchrotron CT provided a detailed insight into the micromechanics of native IVD tissue. The DVC findings showed that there was no slipping at lamella boundaries, and local strain patterns were of a similar distribution to the previously reported elastic network with some heterogeneous areas and maximum strain direction aligned with bundle orientation, suggesting bundle stretching and sliding. This method has the potential to bridge the gap between measures of macro-mechanical properties and the local 3D micro-mechanical environment experienced by cells. This is the first evaluation of strain at the micro scale level in the intact IVD and provides a quantitative framework for future IVD degeneration mechanics studies and testing of tissue engineered IVD replacements. STATEMENT OF SIGNIFICANCE: Synchrotron in-line phase contrast X-ray tomography provided the first visualisation of native intact intervertebral disc microstructural deformation in 3D. For two annulus fibrosus volumes of interest, collagen bundle orientation was quantified and local displacement mapped as strain. Direct evidence of microstructural influence on strain patterns could be seen such as no slipping at lamellae boundaries and maximum strain direction aligned with collagen bundle orientation. Although disc elastic structures were not directly observed, the strain patterns had a similar distribution to the previously reported elastic network. This study presents technical advances and is a basis for future X-ray microscopy, structural quantification and digital volume correlation strain analysis of soft tissue.

Musculoskeletal modelling of the human cervical spine for the investigation of injury mechanisms during axial impacts

Head collisions in sport can result in catastrophic injuries to the cervical spine. Musculoskeletal modelling can help analyse the relationship between motion, external forces and internal loads that lead to injury. However, impact specific musculoskeletal models are lacking as current viscoelastic values used to describe cervical spine joint dynamics have been obtained from unrepresentative quasi-static or static experiments. The aim of this study was to develop and validate a cervical spine musculoskeletal model for use in axial impacts. Cervical spine specimens (C2-C6) were tested under measured sub-catastrophic loads and the resulting 3D motion of the vertebrae was measured. Specimen specific musculoskeletal models were then created and used to estimate the axial and shear viscoelastic (stiffness and damping) properties of the joints through an optimisation algorithm that minimised tracking errors between measured and simulated kinematics. A five-fold cross validation and a Monte Carlo sensitivity analysis were conducted to assess the performance of the newly estimated parameters. The impact-specific parameters were integrated in a population specific musculoskeletal model and used to assess cervical spine loads measured from Rugby union impacts compared to available models. Results of the optimisation showed a larger increase of axial joint stiffness compared to axial damping and shear viscoelastic parameters for all models. The sensitivity analysis revealed that lower values of axial stiffness and shear damping reduced the models performance considerably compared to other degrees of freedom. The impact-specific parameters integrated in the population specific model estimated more appropriate joint displacements for axial head impacts compared to available models and are therefore more suited for injury mechanism analysis.

Crystallinity dependency of the time-dependent mechanical response of polyethylene: application in total disc replacement

Degeneration of the intervertebral disc (IVD) is a leading source of chronic low back pain or neck pain, and represents the main cause of long-term disability worldwide. In the aim to relieve pain, total disc replacement (TDR) is a valuable surgical treatment option, but the expected benefit strongly depends on the prosthesis itself. The present contribution is focused on the synthetic mimic of the native IVD in the aim to optimally restore its functional anatomy and biomechanics, and especially its time-dependency. Semi-crystalline polyethylene (PE) materials covering a wide spectrum of the crystallinity are used to propose new designs of TDR. The influence of the crystallinity on various features of the time-dependent mechanical response of the PE materials is reported over a large strain range by means of dynamic mechanical thermo-analysis and video-controlled tensile mechanical tests. The connection of the stiffness and the yield strength with the microstructure is reported in the aim to propose a model predicting the crystallinity dependency of the response variation with the frequency. New designs of TDR are proposed and implemented into an accurate computational model of a cervical spine segment in order to simulate the biomechanical response under physiological conditions. Predicted in-silico motions are found in excellent agreement with experimental data extracted from published in-vitro studies under compression and different neck movements, namely, rotation, flexion/extension and lateral bending. The simulation results are also criticized by analyzing the local stresses and the predicted biomechanical responses provided by the different prosthetic solutions in terms of time-dependency manifested by the hysteretic behavior under a cyclic movement and the frequency effect.

Biomechanical test protocols to detect minor injury effects in intervertebral discs

Intervertebral discs (IVDs) maintain flexibility of the spine and bear mechanical load. Annulus fibrosus (AF) defects are associated with IVD degeneration and herniation which disrupt biomechanical function and can cause pain. AF puncture injuries can induce IVD degeneration but are needed to inject therapies. Identifying small AF defects with biomechanical testing can be difficult because IVDs have a complex, composite structure and nonlinear biomechanical properties that are dependent on AF fiber tension. It remains unclear how choice of biomechanical testing protocols affect the sensitivity of biomechanical properties to AF injuries. This study determined whether axial preload or magnitude of cyclic axial or torsional testing affected the ability to detect minor AF defects in rat caudal motion segments using ex vivo biomechanical testing. Intact and injured motion segments were subjected to a repeated measures study design with multiple biomechanical testing protocols that varied axial tension-compression force amplitude (±1.6 N, ±8.0 N, ±16.0 N), axial preload (-1.6 N, -8.0 N, -16.0 N, corresponding to -0.1 MPa, -0.5 MPa, and -1.0 MPa, respectively), and torsional rotation angle (±10°, ±15°, and ±20°). Biomechanical properties obtained from the lowest force testing conditions for axial tension-compression (±1.6 N), axial preload (-1.6 N), and angular rotation (±10°) exhibited the largest differences in biomechanical properties between intact and injured conditions. Biomechanical properties determined under low axial force or torsion amplitudes involve less AF fiber tension and were most sensitive to injury. Low force testing protocols are recommended for detecting minor structural AF defects and may enable more precise assessments of IVD injuries, healing or repair.

Osmo-inelastic response of the intervertebral disc

The intervertebral disc exhibits a complex inelastic response characterized by relaxation, hysteresis during cyclic loading and rate dependency. All these inelastic phenomena depend on osmotic interactions between disc tissues and their surrounding chemical environment. Coupling between osmotic and inelastic effects is not fully understood, so this article aimed to study the influence of chemical conditions on the inelastic behaviour of the intervertebral disc in response to different modes of loading. A total of 18 non-frozen ‘motion segments’ (two vertebrae and the intervening soft tissues) were dissected from the cervical spines of mature sheep. The motion segments were loaded in tension, compression and torsion at various loading rates and saline concentrations. Analysis of variance showed that saline concentration significantly influenced inelastic effects in tension and especially in compression (p < 0.05), but not in torsion. Opposite effects were seen in tension and compression. An interpretation of the underlying osmo-inelastic mechanisms is proposed in which two sources of inelastic effects are identified, that is, extracellular matrix rearrangements and fluid exchange created by osmosis.

The Influence of Concordant Complex Posture and Loading Rate on Motion Segment Failure: A Mechanical and Microstructural Investigation

STUDY DESIGN

Microstructural investigation of compression-induced herniation of a lumbar disc held in a concordant complex posture.

OBJECTIVE

To explore the significance of loading rate in a highly asymmetric concordant posture, comparing the mechanisms of failure to an earlier study using a nonconcordant complex posture.

SUMMARY OF BACKGROUND DATA

A recent study with a nonconcordant complex posture (turning in the opposite direction to that which the load is applied) demonstrated the vulnerability of the disc to loading that is borne by one set of oblique-counter oblique fiber sets in the alternating lamellae of the annulus, and aggravated by an elevated loading rate. Given the strain rate-dependent properties of the disc it might be expected that the outcome differs if the posture is reversed.

METHODS

Forty-one motion segments from ovine 16 spines were split into two cohorts; adopting the previously employed low rate (40 mm/min) and surprise rate (400 mm/min) of loading. Both groups of damaged discs were then analyzed microstructurally.

RESULTS

With the lower rate loading the concordant posture significantly reduced the load required to cause disc failure than earlier described for nonconcordant posture (6.9 vs. 8.4 kN), with more direct tears and alternate lamella damage extending to the anterior disc. Contrary to this result, with a surprise rate, the load at failure was significantly increased with the concordant posture (8.08 vs. 6.96 kN), although remaining significantly less than that from a simple flexed posture (9.6 kN). Analysis of the damage modes and postures suggest facet engagement plays a significant role.

CONCLUSION

This study confirms that adding shear to the posture lowers the load at failure, and causes alternate lamella rupture. Load at failure in a complex posture is not determined by loading rate alone. Rather, the strain rate-dependent properties of the disc influence which elements of the system are brought into play.

LEVEL OF EVIDENCE

N/A.

Dynamic shear analysis: a novel method to determine mechanical integrity of normal and torn human acetabular labra: Implications for prediction of outcome of repair

Objectives

The aim of this study was to investigate the structural integrity of torn and non-torn human acetabular labral tissue.

Methods

A total of 47 human labral specimens were obtained from a biobank. These included 22 torn specimens and 25 control specimens from patients undergoing total hip arthroplasty with macroscopically normal labra. The specimens underwent dynamic shear analysis using a rheometer to measure storage modulus, as an indicator of structural integrity.

Results

There was a significant difference in the storage modulus between torn (mean modulus = 2144.08 Pa) and non-torn (3178.1 Pa) labra (p = 0.0001).

Conclusion

The acetabular labrum of young patients with a tear has significantly reduced structural integrity compared with a non-torn labrum in older patients with end-stage osteoarthritis. This study contributes to the understanding of the biomechanics of labral tears, and the observation of reduced structural integrity in torn labra may explain why some repairs fail. Our data demonstrate that labral tears probably have a relatively narrow phenotype, presenting a basis for further investigations that will provide quantifiable data to support their classification and a means to develop a standardized surgical technique for their repair. This study also demonstrates the value of novel biomechanical testing methods in investigating pathological tissues of orthopaedic interest.Cite this article: A. K. Woods, J. Broomfield, P. Monk, F. Vollrath, S. Glyn-Jones. Dynamic shear analysis: a novel method to determine mechanical integrity of normal and torn human acetabular labra: Implications for prediction of outcome of repair. Bone Joint Res 2018;7:440-446. DOI: 10.1302/2046-3758.77.BJR-2017-0282.R2.

Promise, progress, and problems in whole disc tissue engineering

Intervertebral disc degeneration is frequently implicated as a cause of back and neck pain, which are pervasive musculoskeletal complaints in modern society. For the treatment of end stage disc degeneration, replacement of the disc with a viable, tissue-engineered construct that mimics native disc structure and function is a promising alternative to fusion or mechanical arthroplasty techniques. Substantial progress has been made in the field of whole disc tissue engineering over the past decade, with a variety of innovative designs characterized both in vitro and in vivo in animal models. However, significant barriers to clinical translation remain, including construct size, cell source, culture technique, and the identification of appropriate animal models for preclinical evaluation. Here we review the clinical need for disc tissue engineering, the current state of the field, and the outstanding challenges that will need to be addressed by future work in this area.

Does an Annular Puncture Influence the Herniation Path?: An In Vitro Mechanical and Structural Investigation

STUDY DESIGN

A study of mechanically induced herniation in punctured ovine discs followed by structural analysis.

OBJECTIVE

To investigate whether an annular puncture influences the path that herniation takes by providing direct passage for nucleus through the annulus and therefore whether it increases the risk of acute herniation from overload at the site of damage independent of any longer-term degeneration.

SUMMARY OF BACKGROUND DATA

Ten years after treatment with discography both degenerative changes and frequency of herniation have been shown to increase compared to untreated discs. Although the effect of an annular puncture over time has been widely investigated the question of whether it increases the risk of acute herniation has not been resolved.

METHODS

The posterolateral annuli of healthy ovine lumbar discs were punctured with either a 25-gauge (n = 8) or a larger 18-gauge (n = 8) needle and then compressed in a flexed posture of 10° until initial indications of failure. The entire volume of the disc was visually assessed for structural damage by obtaining progressive, full transverse cross-sections of its entire height thus exposing all regions of the disc.

RESULTS

There was no association between the 25-gauge puncture and disc disruption and herniation. In contrast, nuclear material was observed to migrate through the 18-gauge needle puncture. Disruption of the lateral inner annulus was observed in 12 out of the 16 discs tested.

CONCLUSION

The risk of acute herniation through the puncture site is dependent on the needle diameter used. Under the conditions employed the lateral inner annulus remains the site most vulnerable to disruption independent of the presence of a posterolateral puncture.

LEVEL OF EVIDENCE

N /A.

Stiffening of the nucleus pulposus upon axial loading of the intervertebral disc: An experimental in situ study

Mechanical loading is inherently related to the function and degeneration of the intervertebral disc. We present a series of experiments aimed at measuring the effect of a loading/unloading cycle of the intervertebral disc on the mechanical properties of the nucleus pulposus. The study relies on our new minimally invasive microindenter, which allows us to quantify the storage and loss moduli of the nucleus pulposus by inserting an optomechanical probe in an intact (resected) intervertebral disk through the annulus fibrosis via a small needle. Our results indicate that, under the influence of compressive loading, the nucleus pulposus exhibits a more solid-like behavior.

MRI evaluation of the effects of extension exercises on the disc fluid content and location of the centroid of the fluid distribution

BACKGROUND

McKenzie prone press-up exercises have been hypothesised to reduce intradiscal pressure, allowing fluid to be reabsorbed into the disc, which could improve the internal stability and local chemical milieu of the disc, potentially reducing symptoms.

OBJECTIVE

To investigate the immediate effects of prone press-up exercises on lumbar disc fluid content and movement.

DESIGN

Quantification of MRI changes before and after a single exercise session.

METHODS

The mid-sagittal T₂-weighted MR images of 22 volunteers with low back pain were obtained before and immediately after performing press-up exercises. The whole disc and nucleus regions of the L4-5 and L5-S1 discs were then segmented, and their mean signal intensity (MSI) and signal intensity weighted centroid (SIWC) were computed to estimate disc fluid content and displacement.

RESULTS

There were no significant differences between the MSI and the vertical position of the SIWC of the whole disc before and after extension at either disc level (effect size [ES]: -0.23 to 0.09). There was a significant anterior displacement (0.1 ± 5.4 mm) of the location of the SIWC of the disc after extension exercise at L4-5 (ES: 0.22), but not at L5-S1 (ES: 0.00) or at either level for the nucleus region (ES: -0.06; 0.16).

CONCLUSION

Little evidence was found supporting the hypothesis that press-up exercises affect disc fluid content and distribution. Novel parameters reflecting fluid distribution detected similar or larger effects of the extension than MSI. If such exercises are effective in reducing symptoms, it is likely through other mechanisms than by changing fluid content or distribution.

Compression-rate-dependent nonlinear mechanics of normal and impaired porcine knee joints

Background

The knee joint performs mechanical functions with various loading and unloading processes. Past studies have focused on the kinematics and elastic response of the joint with less understanding of the rate-dependent load response associated with viscoelastic and poromechanical behaviors.

Methods

Forty-five fresh porcine knee joints were used in the present study to determine the loading-rate-dependent force-compression relationship, creep and relaxation of normal, dehydrated and meniscectomized joints.

Results

The mechanical tests of all normal intact joints showed similar strong compression-rate-dependent behavior: for a given compression-magnitude up to 1.2 mm, the reaction force varied 6 times over compression rates. While the static response was essentially linear, the nonlinear behavior was boosted with the increased compression rate to approach the asymptote or limit at approximately 2 mm/s. On the other hand, the joint stiffness varied approximately 3 times over different joints, when accounting for the maturity and breed of the animals. Both a loss of joint hydration and a total meniscectomy greatly compromised the load support in the joint, resulting in a reduction of load support as much as 60% from the corresponding intact joint. However, the former only weakened the transient load support, but the latter also greatly weakened the equilibrium load support. A total meniscectomy did not diminish the compression-rate-dependence of the joint though.

Conclusions

These findings are consistent with the fluid-pressurization loading mechanism, which may have a significant implication in the joint mechanical function and cartilage mechanobiology.

Development of an experimental model of burst fracture with damage characterization of the vertebral bodies under dynamic conditions

BACKGROUND

Burst fractures represent a significant proportion of fractures of the thoracolumbar junction. The recent advent of minimally invasive techniques has revolutionized the surgical treatment of this type of fracture. However mechanical behaviour and primary stability offered by these solutions have to be proved from experimental validation tests on cadaveric specimens. Therefore, the aim of this study was to develop an original and reproducible model of burst fracture under dynamic impact.

METHODS

Experimental tests were performed on 24 cadaveric spine segments (T11-L3). A system of dynamic loading was developed using a modified Charpy pendulum. The mechanical response of the segments (strain measurement on vertebrae and discs) was obtained during the impact by using an optical method with a high-speed camera. The production of burst fracture was validated by an analysis of the segments by X-ray tomography.

FINDINGS

Burst fracture was systematically produced on L1 for each specimen. Strain analysis during impact highlighted the large deformation of L1 due to the fracture and small strains in adjacent vertebrae. The mean reduction of the vertebral body of L1 assessed for all the specimens was around 15%. No damage was observed in adjacent discs or vertebrae.

INTERPRETATION

With this new, reliable and replicable procedure for production and biomechanical analysis of burst fractures, comparison of different types of stabilization systems can be envisaged. The loading system was designed so as to be able to produce loads leading to other types of fractures and to provide data to validate finite element modelling.

Yield Strength Testing in Human Cadaver Nasal Septal Cartilage and L-Strut Constructs

Importance

To our knowledge, yield strength testing in human nasal septal cartilage has not been reported to date. An understanding of the basic mechanics of the nasal septum may help surgeons decide how much of an L-strut to preserve and how much grafting is needed.

Objectives

To determine the factors correlated with yield strength of the cartilaginous nasal septum and to explore the association between L-strut width and thickness in determining yield strength.

Design, Setting, and Participants

In an anatomy laboratory, yield strength of rectangular pieces of fresh cadaver nasal septal cartilage was measured, and regression was performed to identify the factors correlated with yield strength. To measure yield strength in L-shaped models, 4 bonded paper L-struts models were constructed for every possible combination of the width and thickness, for a total of 240 models. Mathematical modeling using the resultant data with trend lines and surface fitting was performed to quantify the associations among L-strut width, thickness, and yield strength. The study dates were November 1, 2015, to April 1, 2016.

Main Outcomes and Measures

The factors correlated with nasal cartilage yield strength and the associations among L-strut width, thickness, and yield strength in L-shaped models.

Results

Among 95 cartilage pieces from 12 human cadavers (mean [SD] age, 67.7 [12.6] years) and 240 constructed L-strut models, L-strut thickness was the only factor correlated with nasal septal cartilage yield strength (coefficient for thickness, 5.54; 95% CI, 4.08-7.00; P < .001), with an adjusted R2 correlation coefficient of 0.37. The mean (SD) yield strength R2 varied with L-strut thickness exponentially (0.93 [0.06]) for set widths, and it varied with L-strut width linearly (0.82 [0.11]) or logarithmically (0.85 [0.17]) for set thicknesses. A 3-dimensional surface model of yield strength with L-strut width and thickness as variables was created using a 2-dimensional gaussian function (adjusted R2 = 0.94). Estimated yield strengths were generated from the model to allow determination of the desired yield strength with different permutations of L-strut width and thickness.

Conclusions and Relevance

In this study of human cadaver nasal septal cartilage, L-strut thickness was significantly associated with yield strength. In a bonded paper L-strut model, L-strut thickness had a more important role in determining yield strength than L-strut width. Surgeons should consider the thickness of potential L-struts when determining the amount of cartilaginous septum to harvest and graft.

Level of Evidence

NA.

Biomechanics of the human intervertebral disc: A review of testing techniques and results

Many experimental testing techniques have been adopted in order to provide an understanding of the biomechanics of the human intervertebral disc (IVD). The aim of this review article is to amalgamate results from these studies to provide readers with an overview of the studies conducted and their contribution to our current understanding of the biomechanics and function of the IVD. The overview is presented in a way that should prove useful to experimentalists and computational modellers. Mechanical properties of whole IVDs can be assessed conveniently by testing 'motion segments' comprising two vertebrae and the intervening IVD and ligaments. Neural arches should be removed if load-sharing between them and the disc is of no interest, and specimens containing more than two vertebrae are required to study 'adjacent level' effects. Mechanisms of injury (including endplate fracture and disc herniation) have been studied by applying complex loading at physiologically-relevant loading rates, whereas mechanical evaluations of surgical prostheses require slower application of standardised loading protocols. Results can be strongly influenced by the testing environment, preconditioning, loading rate, specimen age and degeneration, and spinal level. Component tissues of the disc (anulus fibrosus, nucleus pulposus, and cartilage endplates) have been studied to determine their material properties, but only the anulus has been thoroughly evaluated. Animal discs can be used as a model of human discs where uniform non-degenerate specimens are required, although differences in scale, age, and anatomy can lead to problems in interpretation.

Fatigue responses of the human cervical spine intervertebral discs

Numerous studies have been conducted since more than fifty years to understand the behavior of the human lumbar spine under fatigue loading. Applications have been largely driven by low back pain and human body vibration problems. The human neck also sustains fatigue loading in certain type of civilian occupational and military operational activities, and research is very limited in this area. Being a visco-elastic structure, it is important to determine the stress-relaxation properties of the human cervical spine intervertebral discs to enable accurate simulations of these structures in stress-analysis models. While finite element models have the ability to incorporate viscoelastic material definitions, data specific to the cervical spine are limited. The present study was conducted to determine these properties and understand the responses of the human lower cervical spine discs under large number of cyclic loads in the axial compression mode. Eight disc segments consisting of the adjacent vertebral bodies along with the longitudinal ligaments were subjected to compression, followed by 10,000 cycles of loading at 2 or 4Hz frequency by limiting the axial load to approximately 150 N, and subsequent to resting period, subjected to compression to extract the stress-relaxation properties using the quasi-linear viscoelastic (QLV) material model. The coefficients of the model and disc displacements as a function of cycles and loading frequency are presented. The disc responses demonstrated a plateauing effect after the first 2000 to 4000 cycles, which were highly nonlinear. The paper compares these responses with the "work hardening" phenomenon proposed in clinical literature for the lumbar spine to explain the fatigue behavior of the discs. The quantitative results in terms of QLV coefficients can serve as inputs to complex finite element models of the cervical spine to delineate the local and internal load-sharing responses of the disc segment.

Restoring Segmental Biomechanics Through Nucleus Augmentation: An In Vitro Study

STUDY DESIGN

In vitro biomechanical laboratory study.

OBJECTIVES

The purpose of this study is to evaluate a mechanical treatment to create a degenerative motion segment and the ability of nucleus augmentation to restore biomechanics.

SUMMARY OF BACKGROUND

In cases with an intact annulus fibrosus, the replacement or augmentation of the nucleus pulposus alone may provide a less invasive option to restore normal biomechanics and disk height when compared with spinal fusion or total disk replacement. Laboratory testing allows these changes to be fully characterized. However, without preexisting pathology, nucleus augmentation therapies are difficult to evaluate in vitro.

METHODS

The present study evaluated pure moment bending and compressive biomechanics in 3 states (n=6): (1) intact, (2) after creep loading and nucleus disruption to induce degenerative biomechanical changes, and (3) after nucleus augmentation through an injectable polymer (DiscCell).

RESULTS

Neutral zone and ROM were increased in all modes of bending after the degenerative treatment. The most sensitive mode of bending was lateral bending, with intact ROM (20.0±2.9 degrees) increased to 22.3±2.6 degrees after degenerative treatment and reduced to 18.4±1.6 degrees after injection of the polymer. All bending ROM and NZ changes induced by the degenerative treatment were reversed by nucleus augmentation.

CONCLUSIONS

This material was shown to be effective at altering motion segment biomechanics and restoring disk height during time zero tests. This technique may provide a model to examine the time zero performance of a nucleus augmentation device/material.

Material properties of bovine intervertebral discs across strain rates

The intervertebral disc (IVD) is a complex structure responsible for distributing compressive loading to adjacent vertebrae and allowing the vertebral column to bend and twist. To study the mechanical behaviour of individual components of the IVD, it is common for specimens to be dissected away from their surrounding tissues for mechanical testing. However, disrupting the continuity of the IVD to obtain material properties of each component separately may result in erroneous values. In this study, an inverse finite element (FE) modelling optimisation algorithm has been used to obtain material properties of the IVD across strain rates, therefore bypassing the need to harvest individual samples of each component. Uniaxial compression was applied to ten fresh-frozen bovine intervertebral discs at strain rates of 10^-3-1/s. The experimental data were fed into the inverse FE optimisation algorithm and each experiment was simulated using the subject specific FE model of the respective specimen. A sensitivity analysis revealed that the IVD's response was most dependent upon the Young's modulus (YM) of the fibre bundles and therefore this was chosen to be the parameter to optimise. Based on the obtained YM values for each test corresponding to a different strain rate (ε̇), the following relationship was derived:YM=35.5lnε̇+527.5. These properties can be used in finite element models of the IVD that aim to simulate spinal biomechanics across loading rates.

The effect of six degree of freedom loading sequence on the in-vitro compressive properties of human lumbar spine segments

The complex, direction-dependent, poro-viscoelastic properties of the intervertebral disc (disc) suggest that investigations of the six degree of freedom (6DOF) behaviour may be susceptible to inter-test variation in mechanical response if the disc does not return to initial conditions between loading directions. No studies have quantified the effects of sequential multi-directional loading on the consistency of the compressive response of the disc throughout a 6DOF testing protocol. Therefore, the objective of this study was to determine the effect of 6DOF loading on the compressive properties (stiffness and phase angle) of human discs, as evaluated by a reference compression test performed after each single DOF test. Fourteen intact human functional spinal units (FSU) were tested in each of ±6DOFs (shear directions followed by bending and compression) across four orders of magnitude loading frequencies (0.001-1Hz), followed by reference compression tests while subjected to physiological preload, hydration, and body temperature conditions in a hexapod robot. Repeated measures ANOVA revealed significant within-subjects effects between the reference compression tests for modulus (p<0.001), stiffness (p<0.001), and phase angle (p=0.008). Significant post-hoc pairwise comparisons were initially seen between the control and other reference compression tests for stiffness and modulus after the shear DOFs, however, no significant differences were present after the final reference compression test compared to control. More pronounced effects were seen for stiffness in comparison to modulus and phase angle. These effects may be due to three potentials factors, which include the sequence of testing, the cohort of degenerative specimens, and/or cumulative creep due to the constant application of a follower load. While the sequence of test directions was chosen to minimise the biphasic effect, there may be other sequences, which could result in minimal changes in compressive properties.

Review of the fluid flow within intervertebral discs - How could in vitro measurements replicate in vivo?

By maintaining a balance between external mechanical loads and internal osmotic pressure, fluid content of intervertebral discs constantly alters causing fluctuations in disc hydration, height, diameter and pressure that govern disc temporal response. This paper reviews and discusses the relevant findings of earlier studies on the disc fluid flow with the aim to understand and remedy discrepancies between in vivo and in vitro observations. New results of finite element model studies are also exploited in order to help identify the likely causes for such differences and underlying mechanisms observed in vitro. In vivo measurements of changes in spinal height and disc fluid content/pressure via stadiometry, magnetic resonance imaging and intradiscal pressure measurements have been carried out. They have demonstrated that the disc volume, fluid content, height and nucleus pressure alter depending to a large extent on prior-current external load conditions. Although the diurnal loading lasts on average nearly twice longer than the subsequent resting (16 vs. 8h), the disc completely recovers its height and volume during the latter period through fluid inflow. In view of much longer periods required to recover disc height and pressure in vitro in ovine, porcine, caprine, bovine and rat discs, concerns have been raised on the fluid inflow through the endplates that might be hampered by clogged blood vessels post mortem. Analyses of discrepancies in the flow-dependent recoveries in vivo and in vitro highlight an excessive fluid content in the latter as a likely cause. To replicate in vivo conditions as closely as possible in vitro, preparation and preconditioning of specimens and/or pressure and osmolarity of the culture media in which specimens are immersed should hence be designed in a manner as to diminish disc hydration level and/or fluid transport.

The Effect of Local Hydration Environment on the Mechanical Properties and Unloaded Temporal Changes of Isolated Porcine Annular Samples

Preventing dehydration during in vitro testing of isolated layers of annulus fibrosus tissue may require different test conditions than functional spine units. The purpose of the study was twofold: (A) to quantify changes in mass and thickness of multilayer annulus samples in four hydration environments over 120 min; and (B) to quantify cycle-varying biaxial tensile properties of annulus samples in the four environments. The environments included a saline bath, air, relative humidity control, and misting combined with controlled humidity. The loading protocol implemented 24 cycles of biaxial tensile loading to 20% strain at a rate of 2%/s with 3-, 8-, and 13-min of intermittent rest. Specimen mass increased an average (standard deviation) 72% (11) when immersed for 120 min (p < 0.0001). The air condition and the combined mist and relative humidity conditions reduced mass by 45% (15) and 25% (23), respectively, after 120 min (p < 0.0014). Stress at 16% stretch in the air condition was higher at cycle 18 (18 min of exposure) and cycle 24 (33 min of exposure) compared to all other environments in both the axial and circumferential directions (p < 0.0460). There was no significant change in mass or thickness over time in the relative humidity condition and the change in circumferential stress at 16% stretch between cycles 6 and 24 was a maximum of 0.099 MPa and not statistically significant. Implementation of a controlled relative humidity environment is recommended to maintain hydration of isolated annulus layers during cyclic tensile testing.

Failure of the human lumbar motion-segments resulting from anterior shear fatigue loading

An in-vitro experiment was designed to investigate the mode of failure following shear fatigue loading of lumbar motion-segments. Human male lumbar motion-segments (age 32-42 years, n=6) were immersed in Ringer solution at 37°C and repeatedly loaded, using a modified materials testing machine. Fatigue loading consisted of a sinusoidal shear load from 0 N to 1,500 N (750 N±750 N) applied to the upper vertebra of the motion-segment, at a frequency of 5 Hz. During fatigue experiments, several failure events were observed in the dynamic creep curves. Post-test x-ray, CT and dissection revealed that all specimens had delamination of the intervertebral disc. Anterior shear fatigue predominantly resulted in fracture of the apophyseal processes of the upper vertebrae (n=4). Exposure to the anterior shear fatigue loading caused motion-segment instability and resulted in vertebral slip corresponding to grade I and 'mild' grade II spondylolisthesis, as observed clinically.

Effect of degeneration on the six degree of freedom mechanical properties of human lumbar spine segments

While the effects of disc degeneration on compression and rotation motions have been studied, there is no data for shear loading. Clinical research has shown that those with low back pain (a potential consequence of degeneration) experience a 75% greater lateral shear force than those without it. Therefore, the aim was to compare the effect of degeneration on spine segment stiffness and phase angle in each of six degree of freedom (6DOF) loading directions. Fourteen intact functional spinal units (FSU) were dissected from human lumbar spines (mean (SD) age 76.2 (11) years, Thompson grades 3 (N = 5, mild), 4 (N = 6, moderate), 5 (N = 3, severe)). Each FSU was tested in ±6DOFs while subjected to a physiological preload, hydration, and temperature (37°C) conditions in a hexapod robot. A one-way ANOVA between degenerated groups was performed on stiffness and phase angle for each DOF. Significant differences in stiffness were found between mild and moderate degenerative groups in lateral shear (p = 0.001), and axial rotation (p = 0.001), where moderate degeneration had decreased stiffness. For phase angle, significant differences were seen in anterior shear (p = 0.017), and axial rotation (p = 0.026), where phase angle for mild degeneration was less than moderate. Trends of stiffness and phase angle changes between degenerative groups were similar within each DOF. Clinically, the identification of the DOFs that are most affected by degeneration could be used in rehabilitation to improve supplemental stabilization of core muscle groups. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1399-1409, 2016.

The Impact of Posture on the Mechanical Properties of a Functional Spinal Unit During Cyclic Compressive Loading

To assess how posture affects the transmission of mechanical energy up the spinal column during vibration, 18 porcine functional spinal units (FSUs) were exposed to a sinusoidal force (1500 ± 1200 N) at 5 Hz for 120 min in either a flexed, extended, or neutral posture. Force and FSU height were measured continuously throughout the collection. From these data, specimen height loss, dynamic stiffness, hysteresis, and parameters from a standard linear solid (SLS) model were determined and analyzed for differences between postures. Posture had an influence on all of these parameters. In extension, the FSU had higher dynamic stiffness values than when neutral or flexed (p < 0.0001). In flexion, the FSU had higher hysteresis than both an extended or neutral posture (p < 0.0001). Height loss was greatest in a flexed posture and smallest in an extended posture (p < 0.0001). In extension, the series spring element in the SLS model had a stiffness value higher than both flexed and neutral posture conditions, whereas the stiffness in the parallel spring was the same between extension and neutral (p < 0.01), both higher than in flexion. Viscosity coefficients were highest in extension compared to both flexed and neutral (p < 0.01). Based on these results, it was determined that posture had a significant influence in determining the mechanical properties of the spine when exposed to cyclic compressive loading.