Mechanisms of changes in the human spinal column in response to static and dynamic axial mechanical loading

The purpose of this study was to analyze the characteristics of the reaction of the human spinal column to axial static and dynamic loading. Segments of human column from dorsal vertebra XI to lumber vertebra III were exposed to axial static (20 mm/min) and dynamic (200 and 500 mm/min) loading. The measured variables included the load value, whole segment deformation, deformation of the anterior surfaces of intervertebral disk of Th_XI–Th_XII and dorsal vertebra XII, and acoustic emission signals (indicator of spongy bone microdestruction). It was found that the rate of increase in body deformation is lower than that in the intervertebral disk deformation and that the central parts of the spinal end plate compress greater than the peripheral parts. This difference was larger during static loading than that during dynamic loading. To produce the same deformation of a spinal segment due to dynamic loading as that due to static one, it is necessary to overcome a stronger resistance of a greater number of trabecular bones than that due to static loading. Therefore, it can be concluded that, first, the value of dynamic loading causing the same segment compression should be higher than the value of static loading, and the main practical significance here is that the dynamic strength of the column is markedly higher than the static one; second, spinal hardness during impact is higher than that under the static condition; and, third, the same degree of deformation should result in a larger volume of microdestructions in the case of dynamic loading than that in the case of static loading, which is evidenced by a significant difference in the number of AE signals accumulated prior to fracture. The number of AE signals is 444.2 ± 308.2 and 85.0 ± 36.6 in the case of dynamic and static loading, respectively (p < 0.05 according to Student’s t test).     

Contribution of vertebral [corrected] bodies, endplates, and intervertebral discs to the compression creep of spinal motion segments

Spinal segments show non-linear behavior under axial compression. It is unclear to what extent this behavior is attributable to the different components of the segment. In this study, we quantified the separate contributions of vertebral bodies and intervertebral discs to creep of a segment. Secondly, we investigated the contribution of bone and osteochondral endplate (endplates including cartilage) to the deformation of the vertebral body. From eight porcine spines a motion segment, a disc and a vertebral body were dissected and subjected to mechanical testing. In an additional test, cylindrical samples, machined from the lowest thoracic vertebrae of 11 porcine spines, were used to compare the deformation of vertebral bone and endplate. All specimens were subjected to three loading cycles, each comprising a loading phase (2.0 MPa, 15 min) and a recovery phase (0.001 MPa, 30 min). All specimens displayed substantial time-dependent height changes. Average creep was the largest in motion segments and smallest in vertebral bodies. Bone samples with endplates displayed substantially more creep than samples without. In the early phase, behavior of the vertebra was similar to that of the disc. Visco-elastic deformation of the endplate therefore appeared dominant. In the late creep phase, behavior of the segment was similar to that of isolated discs, suggesting that in this phase the disc dominated creep behavior, possibly by fluid flow from the nucleus. We conclude that creep deformation of vertebral bodies contributes substantially to creep of motion segments and that within a vertebral body endplates play a major role.     

Can periods of static loading be used to enhance the resistance of the spine to cumulative compression?

Results of in vitro studies conducted on isolated bone specimens have indicated a higher tolerance to static load than exists when exposed to cyclic loading, when controlled for creep rate. If this difference in load tolerance exists, it may be exploited to extend the life of vertebral bone exposed to repetitive compression, and potentially alter the development of spinal injury. However, little work has been conducted on functional spinal units to determine if bone displays this characteristic within an intact joint. Additionally, static loading may result in load redistribution within the intervertebral disc forcing more of the compressive load towards the periphery of the endplate away from the nucleus. In order to examine these potential mechanisms, 218 osteoligamentous porcine functional spinal units were assigned to one of 15 loading scenarios. This involved one of three normalized peak load magnitudes (50%, 70% and 90% of estimated compressive tolerance) and one of five normalized static load applications (0%, 50%, 100%, 200% and 1000% of the total dynamic work duration). Load magnitude significantly altered the resistance to cumulative compression with decreased peak magnitudes corresponding to both increased cumulative load tolerance and increased height loss. Static load periods did not alter the resistance of the spinal unit to cumulative compression or impact the number of cycles tolerated to failure. The insertion of static load periods impacted the total survival time to failure, but only for the 1000% static load group, an exposure unlikely to occur for most in vivo exposures. The insertion of static load periods decreased the amount of height loss during testing which may play a protective role by allowing load redistribution within the vertebral bone and intervertebral disc.     

Neutral zone and range of motion in the spine are greater with stepwise loading than with a continuous loading protocol. An in vitro porcine investigation

Biomechanical testing of the spine has traditionally been performed to help understand the normal function of the spine as well as to evaluate the effects of injury and surgical procedures on spinal behaviour. The overall objective of this investigation was to compare traditional stepwise loading with the recently introduced continuous loading protocol, determining the effect of loading protocol on the mechanical behaviour of the spine. For all tests, a custom spine testing machine was used to apply pure moments of flexion extension, axial rotation, and lateral bending to a maximum of 2 Nm, using six porcine cervical spine specimens (C2-C4). Motions of C2 with respect to C4 were measured with an optoelectronic camera system. Motion parameters calculated were range of motion (ROM), neutral zone (NZ), and the ratio of NZ and ROM. The continuous loading protocol had smaller values for all motion parameters in each loading direction (p<0.05). ROM for the continuous test ranged between 88% and 93% of that of stepwise for the three loading directions. The continuous protocol NZ was 56-75% of that of the stepwise test. The findings of the study demonstrate that the two loading protocols provide differing spinal behaviours.     

Lubricated squeezing flow: a useful method for measuring the viscoelastic properties of soft tissues

Conducting experiments on very soft biological tissues can be difficult. Traditionally, unconfined compression and shear have been used. Here, an improved method of compression testing, lubricated squeezing flow is described. This gives a uniform compression along the squeezing axis and almost uniform equi-biaxial elongation at right angles to the squeezing axis, with minimal shear deformation due to the constant lubrication of the sample surfaces during testing. Sample results for porcine liver obtained using this method are described here.     

An MRI-based method to align the compressive loading axis for human cadaveric knees

There is a need to align the mechanical axis of the tibia with the axis of loading for studies involving tibiofemoral compression to interpret results and to ensure repeatability of loading within and among specimens. Therefore, the objectives of this study were (1) to develop a magnetic resonance imaging (MRI)-based alignment method for use with apparatuses applying tibiofemoral joint compression, (2) to demonstrate the usefulness of the method by aligning cadaveric knees in an apparatus that could apply tibiofemoral joint compression, and (3) to quantify the error associated with the alignment method. A four degree-of-freedom adjustable device was constructed to allow determination and alignment of the mechanical axis of the tibia of cadaveric knee joints with the axis of loading of an apparatus applying tibiofemoral joint compression. MRI was used to determine the locations of bony landmarks in three dimensions defining the mechanical axis of the tibia relative to an initial orientation of the four degree-of-freedom device. Adjustment values of the device were then computed and applied to the device to align the mechanical axis of the tibia with the axis of a compressive loading apparatus. To demonstrate the usefulness of the method, four cadaveric knees were aligned in the compressive loading apparatus. The vectors describing the mechanical axis of the tibia and the loading axis of the apparatus before and after adjustment of the four degree-of-freedom device were computed for each cadaveric knee. After adjustment of the four degree-of-freedom device, the mechanical axis of the tibia was collinear with the loading axis of the apparatus for each cadaveric knee. The errors in the adjustment values introduced by inaccuracies in the MR images were quantified using the Monte Carlo technique. The precisions in the translational and rotational adjustments were 1.20 mm and 0.90 deg respectively. The MR-based alignment method will allow consistent interpretation of results obtained during tibiofemoral compressive studies conducted using the apparatus described in this paper by providing a well-defined loading axis. The alignment method can also be adapted for use with other apparatuses applying tibiofemoral compression.     

An apparatus for tensile and bending tests of perinatal, neonatal, pediatric and adult cadaver osteoligamentous cervical spines

Investigations of biomechanical properties of pediatric cadaver cervical spines subjected to tensile or bending modes of loading are generally limited by a lack of available tissue and limiting sample sizes, both per age and across age ranges. It is therefore important to develop fixation techniques capable of testing individual cadavers in multiple modes of loading to obtain more biomechanical data per subject. In this study, an experimental apparatus and fixation methodology was developed to accommodate cadaver osteoligamentous head-neck complexes from around birth (perinatal) to full maturation (adult) [cervical length: 2.5-12.5 cm; head breadth: 6-15 cm; head length: 6-19 cm] and sequentially test the whole cervical spine in tension, the upper cervical spine in bending and the upper cervical spine in tension. The experimental apparatus and the fixation methodology provided a rigid casting of the head during testing and did not compromise the skull. Further testing of the intact skull and sub-cranial material was made available due to the design of the apparatus and fixation techniques utilized during spinal testing. The stiffness of the experimental apparatus and fixation technique are reported to better characterize the cervical spine stiffness data obtained from the apparatus. The apparatus and fixation technique stiffness was 1986 N/mm. This experimental system provides a stiff and consistent platform for biomechanical testing across a broad age range and under multiple modes of loading.     

A continuous pure moment loading apparatus for biomechanical testing of multi-segment spine specimens

An apparatus is described that enables the application of continuous pure moment loads to multi-segment spine specimens. This loading apparatus allows continuous cycling of the spine between specified flexion and extension (or right and left lateral bending) maximum load endpoints. Using a six-degree-of-freedom load cell and three-dimensional optoelectronic stereophotogrammetry, characteristic displacement versus load hysteresis curves can be generated and analyzed for different spinal constructs of interest. Unlike quasi-static loading, the use of continuous loading permits the analysis of the spine's behaviour within the neutral zone. This information is of particular clinical significance given that the instability of a spinal segment is related to its flexibility within the neutral zone. Representative curves for the porcine lumbar spine in flexion-extension and lateral bending are presented to illustrate the capabilities of this system.     

Mechanism of fractures of adjacent and augmented vertebrae following simulated vertebroplasty

Percutaneous vertebroplasty (VP) is a minimally invasive procedure that is used to treat osteoporosis-induced vertebral compression fractures (OVCFs). Frequently observed complications are fractures of adjacent and augmented vertebrae. In the present work, mechanisms for these fractures are presented. Fresh 4-level osteoporotic thoracic motion segments were tested. Both ends of the specimen were mounted. The lower level of the free vertebra was compressively fractured and followed by an injection of a 3.5 mL of a PMMA bone cement. Three steps of fatigue loading (5 Hz for 5 h) were incrementally and vertically applied on the specimens from 650 N to 950 N to 1150 N. Specimens of intact, compressively fractured, cement augmented and post-fatigued loading were radiographed for the measurement of deformations of the vertebra, the canal, and the foramen. At the end of fatigue loading, the vertebrae were sliced for micro morphologic analysis. The largest height loss after fatigue loading was at the posterior region of the augmented vertebra. In the augmented vertebra, fissures were found along the bone-cement interface. These fissures split the cement and the trabeculae and propagated into the vertebrae and the endplates. The compactness ratio of the trabeculae region of the adjacent cranial vertebra was higher than that for intact and adjacent caudal ones. We attribute the fracture of the augmented vertebra, following simulated VP, to the initiation of fissures along the cement-bone interface, which, in turn, may be due to uneven deformation of the vertebra. Fracture of the adjacent cranial vertebra is attributed to collapse of its trabeculae.     

3-dimensional registration of micro-movements of vertebral implants by a specially adapted measuring system using magnetic transducers]

The aim of this study was to develop a biomechanical experimental set-up to quantify motion of ventrally inserted spinal implants at the implant/bone interface. The model we used was the vertebral column of the calf. Lumbar vertebrae L2 to L4 were "instrumented" with a screw-rod system. The adjacent vertebrae L1 and L5 were connected to a servohydraulic testing machine and axial compression applied. Shortening of the specimen and three-dimensional movement of the most cranial implant relative to the bone was recorded using 3 electromagnetic transducers. 100,000 cycles of axial loading varying between -0.5 kN and -1 kN were applied. Static shortening of the specimen of 8.5 mm and an elastic movement of 180 to 280 microns were measured. The greatest amplitude of single-plane motion was recorded in the sagittal plane in both static and elastic modes. Implant motion within each cycle was recorded accurately as load-displacement curves within a range of 1.35 to 30 microns. With this test set-up, primary stability of different spinal implant systems can be compared. The use of electromagnetic transducers permits three-dimensional implant motion analysis even when only a mono-axial servohydraulic testing machine is available.     

Simulation of the behaviour of the L1 vertebra for different material properties and loading conditions

Three-dimensional finite element models of the thoracolumbar junction (T12–L2) and isolated L1 vertebra were developed to investigate the role of material properties and loading conditions on vertebral stresses and strains to predict fracture risk. The geometry of the vertebrae was obtained from computed tomography images. The isolated vertebra model included an L1 vertebra loaded through polymethylmethacrylate plates located at the top and bottom of the vertebra, and the segment model included T12 to L2 vertebrae and seven ligaments, fibrous intervertebral discs and facet joints. Each model was examined with both homogeneous and spatially varying bone tissue properties. Stresses and strains were compared for uniform compression and flexion. Including material heterogeneity remarkably reduced the stiffness of the isolated L1 vertebra and increased the magnitudes of the minimum principal strains and stresses in the mid-transverse section. The stress and strain distributions further changed when physiological loading was applied to the L1 vertebra. In the segment models, including heterogeneous material properties increased the magnitude of the minimum principal strain by 158% in the centre of the mid-transverse section. Overall, the inclusion of heterogeneity and physiological loading increased the magnitude of the strains up to 346% in flexion and 273% in compression.     

Modelling and simulation of porcine liver tissue indentation using finite element method and uniaxial stress–strain data

We hypothesize that both compression and elongation stress–strain data should be considered for modeling and simulation of soft tissue indentation. Uniaxial stress–strain data were obtained from in vitro loading experiments of porcine liver tissue. An axisymmetric finite element model was used to simulate liver tissue indentation with tissue material represented by hyperelastic models. The material parameters were derived from uniaxial stress–strain data of compressions, elongations, and combined compression and elongation of porcine liver samples. in vitro indentation tests were used to validate the finite element simulation. Stress–strain data from the simulation with material parameters derived from the combined compression and elongation data match the experimental data best. This is due to its better ability in modeling 3D deformation since the behavior of biological soft tissue under indentation is affected by both its compressive and tensile characteristics. The combined logarithmic and polynomial model is somewhat better than the 5-constant Mooney–Rivlin model as the constitutive model for this indentation simulation.     

Radiography used to measure internal spinal cord deformation in an in vivo rat model

Little is known about the internal mechanics of the in vivo spinal cord during injury. The objective of this study was to develop a method of tracking internal and surface deformation of in vivo rat spinal cord during compression using radiography. Since neural tissue is radio-translucent, radio-opaque markers were injected into the spinal cord.Two tantalum beads (260 µm) were injected into the cord (dorsal and ventral) at C5 of nine anesthetized rats. Four beads were glued to the lateral surface of the cord, caudal and cranial to the injection site. A compression plate was displaced 0.5 mm, 2 mm, and 3 mm into the spinal cord and lateral X-ray images were taken before, during, and after each compression for measuring bead displacements. Potential bead migration was monitored for by comparing displacements of the internal and glued surface beads.Dorsal beads moved significantly more than ventral beads with a range in averages of 0.57–0.71 mm and 0.31–0.35 mm respectively. Bead displacements during 0.5 mm compressions were significantly lower than 2 mm and 3 mm compressions. There was no statistically significant migration of the internal beads.The results indicate the merit of this technique for measuring in vivo spinal cord deformation. The pattern of bead displacements illustrates the complex internal and surface deformations of the spinal cord during transverse compression. This information is needed for validating physical and finite element spinal cord surrogates and to define relationships between loading parameters, internal cord deformation, and biological and functional outcomes.     

Determination of elastomeric foam parameters for simulations of complex loading

BACKGROUND: Finite element (FE) analysis has shown promise for the evaluation of elastomeric foam personal protection devices. Although appropriate representation of foam materials is necessary in order to obtain realistic simulation results, material definitions used in the literature vary widely and often fail to account for the multi-mode loading experienced by these devices. This study aims to provide a library of elastomeric foam material parameters that can be used in FE simulations of complex loading scenarios. METHOD OF APPROACH: Twelve foam materials used in footwear were tested in uni-axial compression, simple shear and volumetric compression. For each material, parameters for a common compressible hyperelastic material model used in FE analysis were determined using: (a) compression; (b) compression and shear data; and (c) data from all three tests. RESULTS: Material parameters and Drucker stability limits for the best fits are provided with their associated errors. The material model was able to reproduce deformation modes for which data was provided during parameter determination but was unable to predict behavior in other deformation modes. CONCLUSIONS: Simulation results were found to be highly dependent on the extent of the test data used to determine the parameters in the material definition. This finding calls into question the many published results of simulations of complex loading that use foam material parameters obtained from a single mode of testing. The library of foam parameters developed here presents associated errors in three deformation modes that should provide for a more informed selection of material parameters.     

Effects of freeze–thaw and micro-computed tomography irradiation on structure–property relations of porcine trabecular bone

We study the effects of freeze–thaw and irradiation on structure–property relations of trabecular bone. We measure the porosity, apparent density, mineral content, trabecular orientation, trabecular thickness, fractal dimension, surface area, and connectivity of trabecular bone using micro-computed tomography (micro-CT) and relate them to Young?s modulus and ultimate strength measured by uniaxial compression testing. The analysis is done on six-month porcine trabecular bone from femoral heads. The effects of freeze–thaw are studied by using bones from three different groups: fresh bone and bones frozen for one and five years. We find that the porosity and apparent density have most dominant influence on the elastic modulus and strength of fresh bone. Also, five years of freezing lowers both Young?s modulus and ultimate strength of trabecular bone. Additionally, the effects of radiation are investigated by comparing Young?s modulus before and after micro-CT exposure. We find that the micro-CT irradiation has a negligible effect on the Young?s modulus of trabecular bone. These findings provide insights on the effects of tissue preservation and imaging on properties of trabecular bone.     

Effect of sulfated glycosaminoglycan digestion on the transverse permeability of medial collateral ligament

Dermatan and chondroitin sulfate glycosaminoglycans (GAGs) comprise over 90% of the GAG content in ligament. Studies of their mechanical contribution to soft tissues have reported conflicting results. Measuring the transient compressive response and biphasic material parameters of the tissue may elucidate the contributions of GAGs to the viscoelastic response to deformation. The hypotheses of the current study were that digestion of sulfated GAGs would decrease compressive stress and aggregate modulus while increasing the permeability of porcine medial collateral ligament (MCL). Confined compression stress relaxation experiments were carried out on porcine MCL and tissue treated with chondroitinase ABC (ChABC). Results were fit to a biphasic constitutive model to derive permeability and aggregate modulus. Bovine articular cartilage was used as a benchmark tissue to verify that the apparatus provided reliable results. GAG digestion removed up to 88% of sulfated GAGs from the ligament. Removal of sulfated GAGs increased the permeability of porcine MCL nearly 6-fold versus control tissues. Peak stress decreased significantly. Bovine articular cartilage exhibited the typical reduction of GAG content and resultant decreases in stress and modulus and increases in permeability with ChABC digestion. Given the relatively small amount of GAG in ligament (<1% of tissue dry weight) and the significant change in peak stress and permeability upon removal of GAGs, sulfated GAGs may play a significant role in maintaining the apposition of collagen fibrils in the transverse direction, thus supporting dynamic compressive loads experienced by the ligament during complex joint motion.     

Determination of elastomeric foam parameters for simulations of complex loading

Background: Finite element (FE) analysis has shown promise for the evaluation of elastomeric foam personal protection devices. Although appropriate representation of foam materials is necessary in order to obtain realistic simulation results, material definitions used in the literature vary widely and often fail to account for the multi-mode loading experienced by these devices. This study aims to provide a library of elastomeric foam material parameters that can be used in FE simulations of complex loading scenarios.

Method of Approach: Twelve foam materials used in footwear were tested in uni-axial compression, simple shear and volumetric compression. For each material, parameters for a common compressible hyperelastic material model used in FE analysis were determined using: (a) compression; (b) compression and shear data; and (c) data from all three tests.

Results: Material parameters and Drucker stability limits for the best fits are provided with their associated errors. The material model was able to reproduce deformation modes for which data was provided during parameter determination but was unable to predict behavior in other deformation modes.

Conclusions: Simulation results were found to be highly dependent on the extent of the test data used to determine the parameters in the material definition. This finding calls into question the many published results of simulations of complex loading that use foam material parameters obtained from a single mode of testing. The library of foam parameters developed here presents associated errors in three deformation modes that should provide for a more informed selection of material parameters.     

Prior storage conditions and loading rate affect the in vitro fracture response of spinal segments under impact loading

Traumatic injuries of the spine are mostly the consequence of rapid overload e.g. impact loading. In vitro investigations on this topic usually encompass biomechanical testing using frozen/thawed specimens and employ quasi-static loading conditions. It is generally accepted that a freezing/thawing cycle does not alter mechanical properties for slow loading rates. However, this has never been investigated for high impact velocities. In order to assess the effects of freezing/thawing and the influence of different impact velocities, we loaded 27 fresh and 15 frozen/thawed cadaveric rabbit spinal segments (intervertebral disc with one third of the adjacent vertebrae) with different impact energies and velocities using a custom-made, dropped-weight loading device. Endplate fractures were assessed by micro-CT scans. Specimen dimensions (disk, bone, and total height) and vertebrae bone density (BV/TV) were compared pre- and post-trauma. Energy absorption by spinal segments was quantified by measuring the initial ball rebound. We found that freezing/thawing increased endplate fracture frequency and decreased the energy absorption of the segments. Higher impact velocities increased the energy absorption, while higher impact energy increased both energy absorption and fracture frequency. Two conclusions are drawn: first, under impact loading, freezing alters permanently the biomechanical response, and second, for different impact velocities, different fracture initiation mechanisms apply. Therefore, quasi-static loading of frozen/thawed spinal segments is not a valid model for traumatic endplate injuries. However, caution should be exercised in extrapolating these findings to human vertebrae until tests on larger vertebrae are performed.     

A novel experimental procedure based on pure shear testing of dermatome-cut samples applied to porcine skin

This paper communicates a novel and robust method for the mechanical testing of thin layers of soft biological tissues with particular application to porcine skin. The key features include the use of a surgical dermatome and the highly defined deformation kinematics achieved by pure shear testing. Thin specimens of accurate thickness were prepared using a dermatome and were subjected to different quasi-static and dynamic loading protocols. Although simple in its experimental realisation, pure shear testing provides a number of advantages over other classic uni- and biaxial testing procedures. The preparation of thin specimens of porcine dermis, the mechanical tests as well as first representative results are described and discussed in detail. The results indicate a pronounced anisotropy between the directions along and across the cleavage lines and a strain rate-dependent response.