Dynamic stiffness and damping of human intervertebral disc using axial oscillatory displacement under a free mass system

Anatomy-aware computed tomography-to-ultrasound spine registration

BACKGROUND

Ultrasound (US) has demonstrated to be an effective guidance technique for lumbar spine injections, enabling precise needle placement without exposing the surgeon or the patient to ionizing radiation. However, noise and acoustic shadowing artifacts make US data interpretation challenging. To mitigate these problems, many authors suggested using computed tomography (CT)-to-US registration to align the spine in pre-operative CT to intra-operative US data, thus providing localization of spinal landmarks.

PURPOSE

In this paper, we propose a deep learning (DL) pipeline for CT-to-US registration and address the problem of a need for annotated medical data for network training. Firstly, we design a data generation method to generate paired CT-US data where the spine is deformed in a physically consistent manner. Secondly, we train a point cloud (PC) registration network using anatomy-aware losses to enforce anatomically consistent predictions.

METHODS

Our proposed pipeline relies on training the network on realistic generated data. In our data generation method, we model the properties of the joints and disks between vertebrae based on biomechanical measurements in previous studies. We simulate the supine and prone position deformation by applying forces on the spine models. We choose the spine models from 35 patients in VerSe dataset. Each spine is deformed 10 times to create a noise-free data with ground-truth segmentation at hand. In our experiments, we use one-leave-out cross-validation strategy to measure the performance and the stability of the proposed method. For each experiment, we choose generated PCs from three spines as the test set. From the remaining, data from 3 spines act as the validation set and we use the rest of the data for training the algorithm. To train our network, we introduce anatomy-aware losses and constraints on the movement to match the physics of the spine, namely, rigidity loss and bio-mechanical loss. We define rigidity loss based on the fact that each vertebra can only transform rigidly while the disks and the surrounding tissue are deformable. Second, by using bio-mechanical loss we stop the network from inferring extreme movements by penalizing the force needed to get to a certain pose.

RESULTS

To validate the effectiveness of our fully automated data generation pipeline, we qualitatively assess the fidelity of the generated data. This assessment involves verifying the realism of the spinal deformation and subsequently confirming the plausibility of the simulated ultrasound images. Next, we demonstrate that the introduction of the anatomy-aware losses brings us closer to state-of-the-art (SOTA) and yields a reduction of 0.25 mm in terms of target registration error (TRE) compared to using only mean squared error (MSE) loss on the generated dataset. Furthermore, by using the proposed losses, the rigidity loss in inference decreases which shows that the inferred deformation respects the rigidity of the vertebrae and only introduces deformations in the soft tissue area to compensate the difference to the target PC. We also show that our results are close to the SOTA for the simulated US dataset with TRE of 3.89 mm and 3.63 mm for the proposed method and SOTA respectively. In addition, we show that our method is more robust against errors in the initialization in comparison to SOTA and significantly achieves better results (TRE of 4.88 mm compared to 5.66 mm) in this experiment.

CONCLUSIONS

In conclusion, we present a pipeline for spine CT-to-US registration and explore the potential benefits of utilizing anatomy-aware losses to enhance registration results. Additionally, we propose a fully automatic method to synthesize paired CT-US data with physically consistent deformations, which offers the opportunity to generate extensive datasets for network training. The generated dataset and the source code for data generation and registration pipeline can be accessed via https://github.com/mfazampour/medphys_ct_us_registration.

Variability of intervertebral joint stiffness between specimens and spine levels

Introduction: Musculoskeletal multibody models of the spine can be used to investigate the biomechanical behaviour of the spine. In this context, a correct characterisation of the passive mechanical properties of the intervertebral joint is crucial. The intervertebral joint stiffness, in particular, is typically derived from the literature, and the differences between individuals and spine levels are often disregarded. Methods: This study tested if an optimisation method of personalising the intervertebral joint stiffnesses was able to capture expected stiffness variation between specimens and between spine levels and if the variation between spine levels could be accurately captured using a generic scaling ratio. Multibody models of six T12 to sacrum spine specimens were created from computed tomography data. For each specimen, two models were created: one with uniform stiffnesses across spine levels, and one accounting for level dependency. Three loading conditions were simulated. The initial stiffness values were optimised to minimize the kinematic error. Results: There was a range of optimised stiffnesses across the specimens and the models with level dependent stiffnesses were less accurate than the models without. Using an optimised stiffness substantially reduced prediction errors. Discussion: The optimisation captured the expected variation between specimens, and the prediction errors demonstrated the importance of accounting for level dependency. The inaccuracy of the predicted kinematics for the level-dependent models indicated that a generic scaling ratio is not a suitable method to account for the level dependency. The variation in the optimised stiffnesses for the different loading conditions indicates personalised stiffnesses should also be considered load-specific.

Biodynamic responses of adolescent idiopathic scoliosis exposed to vibration

Patients with adolescent idiopathic scoliosis suffer severe health issues. The unclear dynamic biomechanical characteristics of scoliosis were needed to be explored to improve the prevention and treatment in clinics. Validated 3D finite element (FE) models of thoracolumbosacral spine (T1-S1) both with and without scoliosis were developed from computed tomography (CT) images. Modal and harmonic analyses were performed to investigate the biomechanical responses of the spinal models to vibration. Resonant frequencies of the scoliotic model were lower than those of the model without scoliosis. Peak amplitudes occurred at vibrational frequencies close to the modal resonant frequencies, which caused the deformed thoracic segment in scoliosis suffered the maximum amplitude. The stresses on vertebrae and intervertebral discs in the scoliotic model derived from vibrations were significantly larger than those in the non-scoliosis model, and heterogeneously concentrated on the scoliotic thoracic segment. In conclusion, the scoliotic spine in the patients with Lenke 1BN scoliosis is more prone to injuries than the non-scoliotic spine while vibrating. Scoliotic thoracic segments in patients with Lenke 1BN scoliosis were the more vulnerable and sensitive component of the T1-S1 spine to vibration than lumbar spines. This study suggested that vibration would impair the scoliotic spines, and patients with Lenke 1BN scoliosis should avoid exposure to vibration, especially the low-frequency vibration.

Region- and degeneration dependent stiffness distribution in intervertebral discs derived by shear wave elastography

Information on the local stiffness characteristics of the intervertebral disc (IVD) is crucial for the understanding of its structure-function properties in health and disease and may improve numerical modeling. Previous studies have attempted to map local tissue stiffness by sectioning the disc and performing mechanical testing on these discrete tissue units, which is technically challenging and may bias the results. Shear wave elastography (SWE) represents a nondestructive alternative that can provide spatially continuous elasticity estimates. We investigated the feasibility of SWE for human intervertebral disc elasticity mapping in a laboratory setting. To this end, global spinal segment mechanical behavior was determined in 6 loading directions and served as ground truth data for the validation of the approach. Subsequently, the cranial spinal vertebra was removed and shear wave elastographic scans of the IVD were acquired. SWE-measurements were reconstructed into three-dimensional elastographic maps, discretized into distinct IVD regions and correlated with global segment mechanical parameters. SWE-derived Young's modulus estimates were compared among different regions and as a function of their state of degeneration. We found annulus shear wave speed to be moderately correlated with segment mechanical behavior irrespective of the loading direction whereas shear wave speed in the nucleus pulposus showed a very weak association (mean (SD) absolute Pearson correlation coefficients: 0.51 (0.14) and 0.17 (0.12), respectively). Young's modulus mapping of the intervertebral disc revealed stiffness to be highest in the ventral annulus with a stiffness decrease both circumferentially towards the dorsal aspect as well as towards the center of the disc. SWE hence provides a valid alternative to disc sectioning and piecewise mechanical testing.

Finite element investigation on the dynamic mechanical properties of low-frequency vibrations on human L2-L3 spinal motion segments with different degrees of degeneration

Exposure to low-frequency vibration is harmful to human lumbar health. However, the dynamic mechanical properties of lumbar spines with varying degrees of degeneration during time-dependent vibration remain incompletely understood. In this study, four poroelastic finite element models of human L2-L3 spinal motion segments, including the non-degeneration and the mild, moderate, and serious degeneration, were established. One-hour low-frequency vibrations with different frequencies were applied. Then, the dynamic mechanical properties of different degenerated lumbar models under the same vibration and the same lumbar model under vibrations at different frequencies were investigated. The results indicated and implied that the negative influences of 1-h vibration on the dynamic mechanical properties of the non-degenerated and mildly degenerated models were similar, but became obvious for the moderately and seriously degenerated models with time. Therefore, the damage caused by low-frequency vibration on the degenerated spinal motion segments was more serious compared with that on the healthy one. Meanwhile, the dynamic mechanical properties of the same lumbar model under vibrations at different frequencies expressed the negligible differences when the vibration frequency was not close to the lumbar natural frequency. Thus, the effects of the 1-h vibrations at different frequencies on one spinal motion segment were similar. Vibration frequency sensitivity analysis on the dynamic characteristics of human L2-L3 spinal motion segments with different degrees of degeneration.

Differential response to vibration of three forms of scoliosis during axial cyclic loading: a finite element study

Background

Scoliosis is a serious disease that can affect all segments of society. Few studies have investigated the response to vibration of differing sinusoidal axial cyclic loading frequencies for different forms of scoliosis in the lumbar spine.

Methods

In this study, four finite element models, comprising a healthy spine, Lenke-A, Lenke-B and Lenke-C scoliosis of the lumbar S1-L1 region were developed. Modal analysis extracted resonant frequencies of the FE models with an upper body mass of 40 kg and 400 N preload. A transient dynamic analysis was performed to obtain the response to vibration of models under a sinusoidal axial loading of ± 40N at frequencies of 3, 5, 7, 9, 11 and 13 Hz using an upper body mass of 40 kg and 400 N preload.

Results

The first-order resonant frequencies of healthy, Lenke-A, Lenke-B and Lenke-C spines were 9.2, 3.9, 4.6 and 5.7 Hz, respectively. A Lenke-A lumbar spine was more likely to deform at a lower vibration frequency and Lenke-C deformed more easily at a higher vibration frequency. Furthermore, the vibration amplitude in the Y-direction (left-right) was greatest and least in the Z-direction (top-bottom). The frequency of cyclic loading closest to the resonant frequency resulted in a maximum value of peak-to-peak vibrational displacement. Furthermore, the vibrational amplitudes in patients with scoliosis were larger than they were in healthy subjects. In addition, axial displacement of the vertebrae in the healthy spine changed steadily whereas fluctuations in the scoliotic vertebrae in scoliosis patients were greater than that of other vertebrae.

Conclusions

Different forms of scoliosis may have different vibrational characteristics, the scoliotic vertebrae being the weak link in scoliosis under loading condition of whole body vibration. Scoliosis was more sensitive to this form of vibration. Where the frequency of axial cyclic vibrational loading of the lumbar spine was closer to its resonant frequency, the vibrational amplitude was larger. These results suggest that vibration will exacerbate the degree of scoliosis and so such patients should reduce their exposure to vibration. Clinical treatment should pay attention to the scoliotic vertebrae and reduce their vibration. These findings may assist in the clinical prevention and treatment of scoliosis.

Biomechanics of Lumbar Motion-Segments in Dynamic Compression

Recent epidemiology studies have reported increase in lumbar spine injuries in frontal crashes. Whole human body finite element models (FEHBM) are frequently used to delineate mechanisms of such injuries. However, the accuracy of these models in mimicking the response of human spine relies on the characterization data of the spine model. The current study set out to generate characterization data that can be input to FEHBM lumbar spine, to obtain biofidelic responses from the models. Twenty-five lumbar functional spinal units were tested under compressive loading. A hydraulic testing machine was used to load the superior ends of the specimens. A 75N load was placed on the superior PMMA to remove the laxity in the joint and mimic the physiological load. There were three loading sequences, namely, preconditioning, 0.5 m/s (non-injurious) and 1.0 m/s (failure). Forces and displacements were collected using six-axis load cell and VICON targets. In addition, acoustic signals were collected to identify the times of failures. Finally, response corridors were generated for the two speeds. To demonstrate the corridors, GHBMC FE model was simulated in frontal impact condition with the default and updated lumbar stiffness. Bi-linear trend was observed in the force versus displacement plots. In the 0.5 m/s tests, mean toe- and linear-region stiffnesses were 0.96±0.37 and 2.44±0.92 kN/mm. In 1.0 m/s tests, the toe and linear-region stiffnesses were 1.13±0.56 and 4.6±2.5 kN/mm. Lumbar joints demonstrated 2.5 times higher stiffness in the linear-region when the loading rate was increased by 0.5 m/s.

The Use of Body Worn Sensors for Detecting the Vibrations Acting on the Lower Back in Alpine Ski Racing

This study explored the use of body worn sensors to evaluate the vibrations that act on the human body in alpine ski racing from a general and a back overuse injury prevention perspective. In the course of a biomechanical field experiment, six male European Cup-level athletes each performed two runs on a typical giant slalom (GS) and slalom (SL) course, resulting in a total of 192 analyzed turns. Three-dimensional accelerations were measured by six inertial measurement units placed on the right and left shanks, right and left thighs, sacrum, and sternum. Based on these data, power spectral density (PSD; i.e., the signal's power distribution over frequency) was determined for all segments analyzed. Additionally, as a measure expressing the severity of vibration exposure, root-mean-square (RMS) acceleration acting on the lower back was calculated based on the inertial acceleration along the sacrum's longitudinal axis. In both GS and SL skiing, the PSD values of the vibrations acting at the shank were found to be largest for frequencies below 30 Hz. While being transmitted through the body, these vibrations were successively attenuated by the knee and hip joint. At the lower back (i.e., sacrum sensor), PSD values were especially pronounced for frequencies between 4 and 10 Hz, whereas a corresponding comparison between GS and SL revealed higher PSD values and larger RMS values for GS. Because vibrations in this particular range (i.e., 4 to 10 Hz) include the spine's resonant frequency and are known to increase the risk of structural deteriorations/abnormalities of the spine, they may be considered potential components of mechanisms leading to overuse injuries of the back in alpine ski racing. Accordingly, any measure to control and/or reduce such skiing-related vibrations to a minimum should be recognized and applied. In this connection, wearable sensor technologies might help to better monitor and manage the overall back overuse-relevant vibration exposure of athletes in regular training and or competition settings in the near future.

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.

Design Requirements for Annulus Fibrosus Repair: Review of Forces, Displacements, and Material Properties of the Intervertebral Disk and a Summary of Candidate Hydrogels for Repair

There is currently a lack of clinically available solutions to restore functionality to the intervertebral disk (IVD) following herniation injury to the annulus fibrosus (AF). Microdiscectomy is a commonly performed surgical procedure to alleviate pain caused by herniation; however, AF defects remain and can lead to accelerated degeneration and painful conditions. Currently available AF closure techniques do not restore mechanical functionality or promote tissue regeneration, and have risk of reherniation. This review determined quantitative design requirements for AF repair materials and summarized currently available hydrogels capable of meeting these design requirements by using a series of systematic PubMed database searches to yield 1500+ papers that were screened and analyzed for relevance to human lumbar in vivo measurements, motion segment behaviors, and tissue level properties. We propose a testing paradigm involving screening tests as well as more involved in situ and in vivo validation tests to efficiently identify promising biomaterials for AF repair. We suggest that successful materials must have high adhesion strength (∼0.2 MPa), match as many AF material properties as possible (e.g., approximately 1 MPa, 0. 3 MPa, and 30 MPa for compressive, shear, and tensile moduli, respectively), and have high tensile failure strain (∼65%) to advance to in situ and in vivo validation tests. While many biomaterials exist for AF repair, few undergo extensive mechanical characterization. A few hydrogels show promise for AF repair since they can match at least one material property of the AF while also adhering to AF tissue and are capable of easy implantation during surgical procedures to warrant additional optimization and validation.

Viscoelastic properties of a spinal posterior dynamic stabilisation device

The purpose of this study was to quantify the frequency dependent viscoelastic properties of two types of spinal posterior dynamic stabilisation devices. In air at 37°C, the viscoelastic properties of six BDyn 1 level, six BDyn 2 level posterior dynamic stabilisation devices (S14 Implants, Pessac, France) and its elastomeric components (polycarbonate urethane and silicone) were measured using Dynamic Mechanical Analysis. The viscoelastic properties were measured over the frequency range 0.01-30Hz. The BDyn devices and its components were viscoelastic throughout the frequency range tested. The mean storage stiffness and mean loss stiffness of the BDyn 1 level device, BDyn 2 level device, silicone component and polycarbonate urethane component all presented a logarithmic relationship with respect to frequency. The storage stiffness of the BDyn 1 level device ranged from 95.56N/mm to 119.29N/mm, while the BDyn 2 level storage stiffness ranged from 39.41N/mm to 42.82N/mm. BDyn 1 level device and BDyn 2 level device loss stiffness ranged from 10.72N/mm to 23.42N/mm and 4.26N/mm to 9.57N/mm, respectively. No resonant frequencies were recorded for the devices or its components. The elastic property of BDyn 1 level device is influenced by the PCU and silicone components, in the physiological frequency range. The viscoelastic properties calculated in this study may be compared to spinal devices and spinal structures.

Mechanical Impedance and Its Relations to Motor Control, Limb Dynamics, and Motion Biomechanics

The concept of mechanical impedance represents the interactive relationship between deformation kinematics and the resulting dynamics in human joints or limbs. A major component of impedance, stiffness, is defined as the ratio between the force change to the displacement change and is strongly related to muscle activation. The set of impedance components, including effective mass, inertia, damping, and stiffness, is important in determining the performance of the many tasks assigned to the limbs and in counteracting undesired effects of applied loads and disturbances. Specifically for the upper limb, impedance enables controlling manual tasks and reaching motions. In the lower limb, impedance is responsible for the transmission and attenuation of impact forces in tasks of repulsive loadings. This review presents an updated account of the works on mechanical impedance and its relations with motor control, limb dynamics, and motion biomechanics. Basic questions related to the linearity and nonlinearity of impedance and to the factors that affect mechanical impedance are treated with relevance to upper and lower limb functions, joint performance, trunk stability, and seating under dynamic conditions. Methods for the derivation of mechanical impedance, both those for within the system and material-structural approaches, are reviewed. For system approaches, special attention is given to methods aimed at revealing the correct and sufficient degree of nonlinearity of impedance. This is particularly relevant in the design of spring-based artificial legs and robotic arms. Finally, due to the intricate relation between impedance and muscle activity, methods for the explicit expression of impedance of contractile tissue are reviewed.

The effect of creep on human lumbar intervertebral disk impact mechanics

The intervertebral disk (IVD) is a highly hydrated tissue, with interstitial fluid making up 80% of the wet weight of the nucleus pulposus (NP), and 70% of the annulus fibrosus (AF). It has often been modeled as a biphasic material, consisting of both a solid and fluid phase. The inherent porosity and osmotic potential of the disk causes an efflux of fluid while under constant load, which leads to a continuous displacement phenomenon known as creep. IVD compressive stiffness increases and NP pressure decreases as a result of creep displacement. Though the effects of creep on disk mechanics have been studied extensively, it has been limited to nonimpact loading conditions. The goal of this study is to better understand the influence of creep and fluid loss on IVD impact mechanics. Twenty-four human lumbar disk samples were divided into six groups according to the length of time they underwent creep (tcreep = 0, 3, 6, 9, 12, 15 h) under a constant compressive load of 400 N. At the end of tcreep, each disk was subjected to a sequence of impact loads of varying durations (timp = 80, 160, 320, 400, 600, 800, 1000 ms). Energy dissipation (ΔE), stiffness in the toe (ktoe) and linear (klin) regions, and neutral zone (NZ) were measured. Analyzing correlations with tcreep, there was a positive correlation with ΔE and NZ, along with a negative correlation with ktoe. There was no strong correlation between tcreep and klin. The data suggest that the IVD mechanical response to impact loading conditions is altered by fluid content and may result in a disk that exhibits less clinical stability and transfers more load to the AF. This could have implications for risk of diskogenic pain as a function of time of day or tissue hydration.

Experimental testing on free vibration behaviour for silicone rubbers proposed within lumbar disc prosthesis

This research was focused on the damping capacity study of two types of silicone rubbers proposed as layers within total lumbar disc prostheses of ball-and-socket model. In order to investigate the damping capacity, the two silicone rubber types mainly differing by the molecular mass of polymeric matrix and the filler content, as was emphasized by scanning electron microscopy and differential scanning calorimetry, were subjected to free vibration testing. Using an adapted experimental installation, three kinds of damping testing were realised: tests without samples and tests with three samples of each type of silicone rubber (69 ShA and 99 ShA). The free vibration tests were performed at a frequency of about 6 Hz using a weight of 11.8 kg. The relative damping coefficient was determined by measuring of two successive amplitudes on the vibrogram and calculating of the logarithmic decrement. The test results with silicone rubber samples showed a relative damping coefficient of 0.058 and respectively 0.077, whilst test results without samples showed a relative damping coefficient of 0.042. These silicone rubbers were found to have acceptable damping properties to be used as layers placed inside the prosthetic components.

Delivering mesenchymal stem cells in collagen microsphere carriers to rabbit degenerative disc: reduced risk of osteophyte formation

Mesenchymal stem cells (MSCs) have the potential to treat early intervertebral disc (IVD) degeneration. However, during intradiscal injection, the vast majority of cells leaked out even in the presence of hydrogel carrier. Recent evidence suggests that annulus puncture is associated with cell leakage and contributes to osteophyte formation, an undesirable side effect. This suggests the significance of developing appropriate carriers for intradiscal delivery of MSCs. We previously developed a collagen microencapsulation platform, which entraps MSCs in a solid microsphere consisting of collagen nanofiber meshwork. These solid yet porous microspheres support MSC attachment, survival, proliferation, migration, differentiation, and matrix remodeling. Here we hypothesize that intradiscal injection of MSCs in collagen microspheres will outperform that of MSCs in saline in terms of better functional outcomes and reduced side effects. Specifically, we induced disc degeneration in rabbits and then intradiscally injected autologous MSCs, either packaged within collagen microspheres or directly suspended in saline, into different disc levels. Functional outcomes including hydration index and disc height were monitored regularly until 6 months. Upon sacrifice, the involved discs were harvested for histological, biochemical, and biomechanical evaluations. MSCs in collagen microspheres showed advantage over MSCs in saline in better maintaining the dynamic mechanical behavior but similar performance in hydration and disc height maintenance and matrix composition. More importantly, upon examination of gross appearance, radiograph, and histology of IVD, delivering MSCs in collagen microspheres significantly reduced the risk of osteophyte formation as compared to that in saline. This work demonstrates the significance of using cell carriers during intradiscal injection of MSCs in treating disc degeneration.

A Comparison of the Human Lumbar Intervertebral Disc Mechanical Response to Normal and Impact Loading Conditions

Thirty-four percent of U.S. Navy high speed craft (HSC) personnel suffer from lower back injury and low back pain, compared with 15 to 20% of the general population. Many of these injuries are specifically related to the intervertebral disc, including discogenic pain and accelerated disc degeneration. Numerous studies have characterized the mechanical behavior of the disc under normal physiological loads, while several have also analyzed dynamic loading conditions. However, the effect of impact loads on the lumbar disc—and their contribution to the high incidence of low back pain among HSC personnel—is still not well understood. An ex vivo study on human lumbar anterior column units was performed in order to investigate disc biomechanical response to impact loading conditions. Samples were subjected to a sequence of impact events of varying duration (Δt = 80, 160, 320, 400, 600, 800, and 1000 ms) and the level of displacement (0.2, 0.5, and 0.8 mm), stiffness k, and energy dissipation ΔE were measured. Impacts of Δt = 80 ms saw an 18–21% rise in k and a 3–7% drop in ΔE compared to the 1000 ms baseline, signaling an abrupt change in disc mechanics. The altered disc mechanical response during impact likely causes more load to be transferred directly to the endplates, vertebral bodies, and surrounding soft tissues and can help begin to explain the high incidence of low back pain among HSC operators and other individuals who typically experience similar loading environments. The determination of a “safety range” for impacts could result in a refinement of design criteria for shock mitigating systems on high-speed craft, thus addressing the low back injury problem among HSC personnel.

Dynamic biomechanical examination of the lumbar spine with implanted total spinal segment replacement (TSSR) utilizing a pendulum testing system

Background

Biomechanical investigations of spinal motion preserving implants help in the understanding of their in vivo behavior. In this study, we hypothesized that the lumbar spine with implanted total spinal segment replacement (TSSR) would exhibit decreased dynamic stiffness and more rapid energy absorption compared to native functional spinal units under simulated physiologic motion when tested with the pendulum system.

Methods

Five unembalmed, frozen human lumbar functional spinal units were tested on the pendulum system with axial compressive loads of 181 N, 282 N, 385 N, and 488 N before and after Flexuspine total spinal segment replacement implantation. Testing in flexion, extension, and lateral bending began by rotating the pendulum to 5°; resulting in unconstrained oscillatory motion. The number of rotations to equilibrium was recorded and bending stiffness (N-m/°) was calculated and compared for each testing mode.

Results

The total spinal segment replacement reached equilibrium with significantly fewer cycles to equilibrium compared to the intact functional spinal unit at all loads in flexion (p<0.011), and at loads of 385 N and 488 N in lateral bending (p<0.020). Mean bending stiffness in flexion, extension, and lateral bending increased with increasing load for both the intact functional spinal unit and total spinal segment replacement constructs (p<0.001), with no significant differences in stiffness between the intact functional spinal unit and total spinal segment replacement in any of the test modes (p>0.18).

Conclusions

Lumbar functional spinal units with implanted total spinal segment replacement were found to have similar dynamic bending stiffness, but absorbed energy at a more rapid rate than intact functional spinal units during cyclic loading with an unconstrained pendulum system. Although the effects on clinical performance of motion preserving devices is not fully known, these results provide further insight into the biomechanical behavior of this device under approximated physiologic loading conditions.

Damping properties of the nucleus pulposus

BACKGROUND

The nucleus pulposus is extremely deformable and it is not uncommon to observe strain amplitudes as large as 12.5% in physiological loading conditions. It has been shown that the nucleus pulposus contributes to the damping properties of the intervertebral disc. The quantification of the damping properties of the nucleus pulposus under physiological large deformations is then a key aspect for its mechanical characterization and for the design of nucleus replacement devices.

METHODS

A specific mechanical device has been developed to encapsulate nucleus pulposus tissues into a deformable and permeable device, while quantifying its water content. The specific damping capacity was defined by dividing the energy loss by the work input. With this device and definition, the specific damping capacity of the bovine coccygeal nucleus pulposus was quantified in large compressive deformations (12.5%) and for frequencies ranging between 10(-2) and 10(1)Hz.

FINDINGS

It is found that the specific damping capacity of the nucleus pulposus of the bovine coccygeal ranged between 18 and 36%. The lowest values of specific damping capacity are found for frequencies corresponding to the dynamics of loads in all day activities such as walking (0.1 to 1Hz).

INTERPRETATION

The nucleus pulposus contributes to dissipate energy under physiological large deformations. However, it seems that the nucleus pulposus is designed so that damping is minimal for frequencies corresponding to moderate daily activities.

Dynamic biomechanical examination of the lumbar spine with implanted total disc replacement using a pendulum testing system

STUDY DESIGN

Biomechanical cadaver investigation.

OBJECTIVE

To examine dynamic bending stiffness and energy absorption of the lumbar spine with and without implanted total disc replacement (TDR) under simulated physiological motion.

SUMMARY OF BACKGROUND DATA

The pendulum testing system is capable of applying physiological compressive loads without constraining motion of functional spinal units (FSUs). The number of cycles to equilibrium observed under pendulum testing is a measure of the energy absorbed by the FSU.

METHODS

Five unembalmed, frozen human lumbar FSUs were tested on the pendulum system with axial compressive loads of 181 N, 282 N, 385 N, and 488 N before and after Synthes ProDisc-L TDR implantation. Testing in flexion, extension, and lateral bending began by rotating the pendulum to 5º resulting in unconstrained oscillatory motion. The number of rotations to equilibrium was recorded and bending stiffness (N·m/º) was calculated and compared for each testing mode.

RESULTS

In flexion/extension, the TDR constructs reached equilibrium with significantly (P < 0.05) fewer cycles than the intact FSU with compressive loads of 282 N, 385 N, and 488 N. Mean dynamic bending stiffness in flexion, extension, and lateral bending increased significantly with increasing load for both the intact FSU and TDR constructs (P < 0.001). In flexion, with increasing compressive loading from 181 N to 488 N, the bending stiffness of the intact FSUs increased from 4.0 N·m/º to 5.5 N·m/º, compared with 2.1 N·m/º to 3.6 N·m/º after TDR implantation. At each compressive load, the intact FSU was significantly stiffer than the TDR (P < 0.05).

CONCLUSION

Lumbar FSUs with implanted TDR were found to be less stiff, but absorbed more energy during cyclic loading with an unconstrained pendulum system. Although the effects on clinical performance of motion-preserving devices are not fully known, these results provide further insight into the biomechanical behavior of these devices under approximated physiological loading conditions.

Transmission of force in the lumbosacral spine during backward falls

STUDY DESIGN

Mathematical model, combined with and verified using human subject data.

OBJECTIVE

(1) To develop and verify a lumped-parameter mathematical model for prediction of spine forces during backward falls; (2) to use this model to evaluate the effect of floor stiffness on spine forces during falls; and (3) to compare predicted impact forces with forces previously measured to fracture the spine.

SUMMARY OF BACKGROUND DATA

Vertebral fractures are the most common osteoporotic fractures and commonly result from falls from standing height. Compliant flooring reduces the force at the ground during a backward fall from standing; however, the effect on spine forces is unknown.

METHODS

A 6-df model of the body was developed and verified using data from 10 human subjects falling from standing onto 3 types of compliant floors (soft: 59 kN/m, medium: 67 kN/m, and firm: 95 kN/m). The simulated ground forces were compared with those measured experimentally. The model was also used to assess the effect of floor stiffness on spine forces at various intervertebral levels.

RESULTS

There was less than 14% difference between model predictions and experimentally measured peak ground reaction forces, when averaged over all floor conditions. When compared with the rigid floor, average peak spine force attenuations of 46%, 43%, and 41% were achieved with the soft, medium, and firm floors, respectively (3.7, 3.9, 4.1 kN vs. 6.9 kN at L4/L5). Spine forces were lower than those at the ground and decreased cranially (4.9, 3.9, 3.7, 3.5 kN at the ground, L5/S1, L4/L5, and L3/L4, respectively, for the soft floor).

CONCLUSION

Lowering the floor stiffness (from 400 to 59 kN/m) can attenuate peak lumbosacral spine forces in a backward fall onto the buttocks from standing by 46% (average peak from 6.9 to 3.7 kN at L4/L5) to values closer to the average tolerance of the spine to fracture (3.4 kN).

Investigation of the intervertebral disc and fused joint dynamics through experimental modal analysis and the receptance coupling method

Identification of intervertebral disc (IVD) dynamics is important in understanding the spine mechanism and behavior. This paper experimentally identifies the dynamics of the bovine caudal IVD using experimental modal analysis and the inverse receptance coupling method. Experimental modal analysis was performed on free-free mounted bovine caudal vertebrae joined by an IVD and a fused IVD joint. Shear, rotational, and axial dynamics of the joints are identified by curve fitting of the frequency response functions, and identifying the damping ratio, stiffness, and modal frequency in each axis. The identified dynamics are compared with the IVD joints with and without fusion. Results provide important insight into IVD dynamics and fused IVD dynamics. This method can be extended to identify human IVD joint dynamics.

Is vibration truly an injurious stimulus in the human spine?

Epidemiological data at one time was taken to suggest that chronic vibrations--for example operating vehicles with low-quality seats--contributed to intervertebral disc degeneration and lower back pain. More recent discussions, based in part upon extended twin studies, have cast doubt upon this interpretation, and question how much of the vibration is actually transmitted to the spine during loading. This review summarizes our recent survey of the current state of knowledge. In particular, we note that current studies are lacking a detailed factorial exploration of frequency, amplitude, and duration; this may be the primary cause for inconclusive and/or contradictory studies. It is our conclusion that vibrations are still an important consideration in discogenic back pain, and further controlled studies are warranted to definitively examine the underlying hypothesis: that chronic vibration can influence IVD cell biology and tissue mechanics.

Development and validation of a predictive bone fracture risk model for astronauts

There are still many unknowns in the physiological response of human beings to space, but compelling evidence indicates that accelerated bone loss will be a consequence of long-duration spaceflight. Lacking phenomenological data on fracture risk in space, we have developed a predictive tool based on biomechanical and bone loading models at any gravitational level of interest. The tool is a statistical model that forecasts fracture risk, bounds the associated uncertainties, and performs sensitivity analysis. In this paper, we focused on events that represent severe consequences for an exploration mission, specifically that of spinal fracture resulting from a routine task (lifting a heavy object up to 60 kg), or a spinal, femoral or wrist fracture due to an accidental fall or an intentional jump from 1 to 2 m. We validated the biomechanical and bone fracture models against terrestrial studies of ground reaction forces, skeletal loading, fracture risk, and fracture incidence. Finally, we predicted fracture risk associated with reference missions to the moon and Mars that represented crew activities on the surface. Fracture was much more likely on Mars due to compromised bone integrity. No statistically significant gender-dependent differences emerged. Wrist fracture was the most likely type of fracture, followed by spinal and hip fracture.

Linear time delay methods and stability analyses of the human spine. Effects of neuromuscular reflex response

Linear stability methods were applied to a biomechanical model of the human musculoskeletal spine to investigate effects of reflex gain and reflex delay on stability. Equations of motion represented a dynamic 18 degrees-of-freedom rigid-body model with time-delayed reflexes. Optimal muscle activation levels were identified by minimizing metabolic power with the constraints of equilibrium and stability with zero reflex time delay. Muscle activation levels and associated muscle forces were used to find the delay margin, i.e., the maximum reflex delay for which the system was stable. Results demonstrated that stiffness due to antagonistic co-contraction necessary for stability declined with increased proportional reflex gain. Reflex delay limited the maximum acceptable proportional reflex gain, i.e., long reflex delay required smaller maximum reflex gain to avoid instability. As differential reflex gain increased, there was a small increase in acceptable reflex delay. However, differential reflex gain with values near intrinsic damping caused the delay margin to approach zero. Forward-dynamic simulations of the fully nonlinear time-delayed system verified the linear results. The linear methods accurately found the delay margin below which the nonlinear system was asymptotically stable. These methods may aid future investigations in the role of reflexes in musculoskeletal stability.

Rest cannot always recover the dynamic properties of fatigue-loaded intervertebral disc

STUDY DESIGN

An in vitro biomechanical study using adolescent porcine discs.

OBJECTIVE

To find the effect of fatigue loading and rest on the dynamic properties of healthy intervertebral disc.

SUMMARY OF BACKGROUND DATA

The fatigue loading is a risk factor for low back pain. The disc dynamic properties describe the shock attenuation capability of disc. Knowledge of effect of fatigue loading and rest on the disc dynamic properties can be beneficial for the study of fatigue loading induced spinal disorder.

METHODS

Specimens were divided into short-term (0.5 hours at 5 Hz, n = 9) and long-term (2 hours at 5 Hz, n = 9) fatigue loading groups. The specimen was applied with fatigue loading, followed by a 12-hours rest, and then applied with the second fatigue loading again. The impulse loading was applied at 0, 10, 20, and 30 minutes during short-term group and at 0, 0.5, 1, and 2 hours during long-term group using a drop-tower apparatus. The stiffness, damping coefficient, and ratio were calculated using impulse loading information. Dynamic properties between first and second fatigue loading were compared.

RESULTS

The stiffness increased, but the damping coefficient and ratio decreased with fatigue loadings. A 12-hours rest can fully restore the dynamic properties of fatigue loaded disc to original status. However, the degradation of dynamic properties during the subsequent fatigue loading was faster than the previous fatigue loading if the disc experienced a long-term fatigue loading.

CONCLUSION

A 12-hours rest can largely restore disc dynamic properties during fatigue loading; hence fully recover disc dynamic properties. A long-term fatigue loading squeezes the disc fluid and injures the disc integrity. The 12-hours rest can recover the fluid loss but not the injury of disc. Hence, the disc degradation during the second fatigue loading was faster than the first fatigue loading.

Subject-specific models of the hindfoot reveal a relationship between morphology and passive mechanical properties

The morphology of the bones, articular surfaces and ligaments and the passive mechanical characteristics of the ankle complex were reported to vary greatly among individuals. The goal of this study was to test the hypothesis that the variations observed in the passive mechanical properties of the healthy ankle complex are strongly influenced by morphological variations. To evaluate this hypothesis six numerical models of the ankle joint complex were developed from morphological data obtained from MRI of six cadaveric lower limbs, and from average reported data on the mechanical properties of ligaments and articular cartilage. The passive mechanical behavior of each model, under a variety of loading conditions, was found to closely match the experimental data obtained from each corresponding specimen. Since all models used identical material properties and were subjected to identical loads and boundary conditions, it was concluded that the observed variations in passive mechanical characteristics were due to variations in morphology, thus confirming the hypothesis. In addition, the average and large variations in passive mechanical behavior observed between the models were similar to those observed experimentally between cadaveric specimens. The results suggest that individualized subject-specific treatment procedures for ankle complex disorders are potentially superior to a one-size-fits-all approach.

Vertebral endplate trauma induces disc cell apoptosis and promotes organ degeneration in vitro

There is a major controversy whether spinal trauma with vertebral endplate fractures can result in post-traumatic disc degeneration. Intervertebral discs, which are adjacent to burst endplates, are frequently removed and an intercorporal spondylodesis is performed. In any case, the biological effects within the discs following endplate fractures are poorly elucidated to date. The aim of our investigations was therefore to establish a novel disc/endplate trauma culture model to reproducibly induce endplate fractures and investigate concurrent disc changes in vitro. This model is based on a full-organ disc/endplate culture system, which has been validated by the authors before. Intervertebral disc/endplate specimens were isolated from Burgundy rabbits and cultured in standard media (DMEM/F12, 10%FCS). Burst endplate fractures were induced in half of the specimens with a custom-made fracture device and subsequently cultured for 9 days. The biological effects such as necrotic or apoptotic cell death and the expression of pro-apoptotic genes and other genes involved in organ degeneration, e.g. matrix metalloproteinases (MMPs) were analyzed. Cell damage was assessed by quantification of the lactate dehydrogenase (LDH) activity in the supernatant. The expression of genes involved in the cellular apoptotic pathway (caspase 3) and the pro-apoptotic proteins FasL and TNF-alpha were monitored. The results demonstrate that LDH levels increased significantly post trauma compared to the control and remained elevated for 3 days. Furthermore, a constant up-regulation of the caspase 3 gene in both disc compartments was present. The pro-apoptotic proteins FasL and TNF-alpha were up regulated predominantly in the nucleus whereas the MMP-1 and -13 transcripts (collagenases) were increased in both disc structures. From this study we can conclude that endplate burst fractures result in both necrotic and apoptotic cell death in nucleus and annulus tissue. Moreover, FasL and TNF-alpha expression by nucleus cells may lead to continued apoptosis induced by Fas- and TNF-alpha receptor bearing cells. In addition TNF-alpha over-expression has potentially deleterious effects on disc metabolism such as over-expression of matrix proteinases. Taken together, the short term biological response of the disc following endplate fracture exhibits characteristics, which may initiate the degeneration of the organ.

The dynamic flexion/extension properties of the lumbar spine in vitro using a novel pendulum system

The biomechanical properties of the ligamentous cadaver spine have been previously examined using a variety of experimental testing protocols. Ongoing technical challenges in the biomechanical testing of the spine include the application of physiologic compressive loads and the application of dynamic bending moments while allowing unconstrained three-dimensional motion. The purpose of this study was to report the development of a novel pendulum apparatus that addressed these challenges and to determine the effects of various axial compressive loads on the dynamic biomechanical properties of the lumbar functional spinal unit (FSU). Lumbar FSUs were tested in flexion and extension under five axial compressive loads chosen to represent physiologic loading conditions. After an initial rotation, the FSUs behaved as a dynamic, underdamped vibrating elastic system. Bending stiffness and coefficient of damping increased significantly as the compressive pendulum load increased. The apparatus described herein is a relatively simple approach to determining the dynamic bending properties of the FSU, and potentially disc arthroplasty devices. It is capable of applying physiologic compressive loads at dynamic rates without constraining the kinematics of the joints, crucial requirements for testing FSUs in vitro.

Role of reflex gain and reflex delay in spinal stability--a dynamic simulation

The goal of this study was to investigate the role of reflex and reflex time delay in muscle recruitment and spinal stability. A dynamic biomechanical model of the musculoskeletal spine with reflex response was implemented to investigate the relationship between reflex gain, co-contraction, and stability in the spine. The first aim of the study was to investigate how reflex gain affected co-contraction predicted in the model. It was found that reflexes allowed the model to stabilize with less antagonistic co-contraction and hence lower metabolic power than when limited to intrinsic stiffness alone. In fact, without reflexes there was no feasible recruitment pattern that could maintain spinal stability when the torso was loaded with 200N external load. Reflex delay is manifest in the paraspinal muscles and represents the time from a perturbation to the onset of reflex activation. The second aim of the study was to investigate the relationship between reflex delay and the maximum tolerable reflex gain. The maximum acceptable upper bound on reflex gain decreased logarithmically with reflex delay. Thus, increased reflex delay and reduced reflex gain requires greater antagonistic co-contraction to maintain spinal stability. Results of this study may help understanding of how patients with retarded reflex delay utilize reflex for stability, and may explain why some patients preferentially recruit more intrinsic stiffness than healthy subjects.