Influence of bone microstructure on the mechanical properties of skull cortical bone - A combined experimental and computational approach.

The strength and compliance of the dense cortical layers of the human skull have been examined since the beginning of the 20th century with the wide range in the observed mechanical properties attributed to natural biological variance. Since this variance may be explained by the difference in structural arrangement of bone tissue, micro-computed tomography (µCT) was used in conjunction with mechanical testing to study the relationship between the microstructure of human skull cortical coupons and their mechanical response. Ninety-seven bone samples were machined from the cortical tables of the calvaria of ten fresh post mortem human surrogates and tested in dynamic tension until failure. A linear response between stress and strain was observed until close to failure, which occurred at 0.6% strain on average. The effective modulus of elasticity for the coupons was 12.01 ± 3.28GPa. Porosity of the test specimens, determined from µCT, could explain only 51% of the variation of their effective elastic modulus. Finite element (FE) models of the tested specimens built from µCT images indicated that modeling the microstructural arrangement of the bone, in addition to the porosity, led to a marginal improvement of the coefficient of determination to 54%. Modulus for skull cortical bone for an element size of 50µm was estimated to be 19GPa at an average. Unlike the load bearing bones of the body, almost half of the variance in the mechanical properties of cortical bone from the skull may be attributed to differences at the sub-osteon (< 50µm) level. ANOVA tests indicated that effective failure stress and strain varied significantly between the frontal and parietal bones, while the bone phase modulus was different for the superior and inferior aspects of the calvarium. The micro FE models did not indicate any anisotropy attributable to the pores observable under µCT.

Lumbar spine endplate fractures: Biomechanical evaluation and clinical considerations through experimental induction of injury.

Lumbar endplate fractures were investigated in different experimental scenarios, however the biomechanical effect of segmental alignment was not outlined. The objectives of this study were to quantify effects of spinal orientation on lumbar spine injuries during single-cycle compressive loads and understand lumbar spine endplate injury tolerance. Twenty lumbar motion segments were compressed to failure. Two methods were used in the preparation of the lumbar motion segments. Group 1 (n = 7) preparation maintained pre-test sagittal lordosis, whereas Group 2 (n = 13) specimens had a free-rotational end condition for the cranial vertebra, allowing sagittal rotation of the cranial vertebra to create parallel endplates. Five Group 1 specimens experienced posterior vertebral body fracture prior to endplate fracture, whereas two sustained endplate fracture only. Group 2 specimens sustained isolated endplate fractures. Group 2 fractures occurred at approximately 41% of the axial force required for Group 1 fracture (p < 0.05). Imaging and specimen dissection indicate endplate injury consistently took place within the confines of the endplate boundaries, away from the vertebral periphery. These findings indicate that spinal alignment during compressive loading influences the resulting injury pattern. This investigation identified the specific mechanical conditions under which an endplate breach will take place. Development of endplate injuries has significant clinical implication as previous research identified internal disc disruption (IDD) and degenerative disc disease (DDD) as long-term consequences of the axial load-shift that occurs following a breach of the endplate. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1084-1091, 2016.

Characterization of Lumbar Spine Annular Disruption in PMHS Using MRI, Cryomicrotomy and Histology Techniques.

Internal intervertebral disc disruption is involved in the onset of a wide range of spinal dysfunction, ultimately affecting not only the disc itself but the surrounding osseous and neural structures as well. The ability of disc to withstand and effectively distribute axial load is dependent upon whether peripherally located annular fibers provide the support necessary to contain and corral the pressure sensitive nucleus. Any alteration in the structures immediate to the nucleus jeopardize this ability. While annular tears and fissures have been thoroughly investigated, one form of internal disc disruption is less well-understood. A network of elastin cross-bridges provides resistance to delamination of the collagenous sheets that comprise the annulus. The current investigation utilized a Nitrogen gas-induced pressure mechanism to disrupt elastin cross links that exist between annular lamellae. Twenty five cadaveric lumbar spine motion segments (mean age: 52±12 yr.) were subjected to the annular disruption protocol. Damage to the annulus was assessed using MRI, cryomicrotome and histological staining procedures. MRI images were compared to cryomicrotome images to determine the ability of standard clinical MRI scans to determine annular damage. In many cases MRI was moderately revealing in terms of damage. Future studies will quantify biomechanical consequences of these low level annular disruptions relative to segmental stability.

Variation of bone layer thicknesses and trabecular volume fraction in the adult male human calvarium.

The human calvarium is a sandwich structure with two dense layers of cortical bone separated by porous cancellous bone. The variation of the three dimensional geometry, including the layer thicknesses and the volume fraction of the cancellous layer across the population, is unavailable in the current literature. This information is of particular importance to mathematical models of the human head used to simulate mechanical response. Although the target geometry for these models is the median geometry of the population, the best attempt so far has been the scaling of a unique geometry based on a few median anthropometric measurements of the head. However, this method does not represent the median geometry. This paper reports the average three dimensional geometry of the calvarium from X-ray computed tomography (CT) imaging and layer thickness and trabecular volume fraction from micro CT (µCT) imaging of ten adult male post-mortem human surrogates (PMHS). Skull bone samples have been obtained and µCT imaging was done at a resolution of 30 µm. Monte Carlo simulation was done to estimate the variance in these measurements due to the uncertainty in image segmentation. The layer thickness data has been averaged over areas of 5mm(2). The outer cortical layer was found to be significantly (p < 0.01; Student's t test) thicker than the inner layer (median of thickness ratio 1.68). Although there was significant location to location difference in all the layer thicknesses and volume fraction measurements, there was no trend. Average distribution and the variance of these metrics on the calvarium have been shown. The findings have been reported as colormaps on a 2D projection of the cranial vault.

Rate-dependent fracture characteristics of lumbar vertebral bodies.

Experimental testing incorporating lumbar columns and isolated components is essential to advance the understanding of injury tolerance and for the development of safety enhancements. This study incorporated a whole column axial acceleration model and an isolated vertebral body model to quantify compression rates during realistic loading and compressive tolerance of vertebrae. Eight lumbar columns and 53 vertebral bodies from 23 PMHS were used. Three-factor ANOVA was used to determine significant differences (p<0.05) in physiologic and failure biomechanics based on compression rate, spinal level, and gender. Results demonstrated a significant increase in ultimate force (i.e., fracture) from lower to higher compression rates. Ultimate stress also increased with compression rate. Displacement and strain to failure were consistent at both compression rates. Differences in ultimate mechanics between vertebral bodies obtained from males and females demonstrated non-significant trends, with female vertebral bodies having lower ultimate force that would be associated with decreased injury tolerance. This was likely a result of smaller vertebrae in that population. Combined with existing literature, results presented in this manuscript contribute to the understanding of lumbar spine tolerance during axial loading events that occur in both military and civilian environments with regard to effects of compression rate and gender.

Effects of acceleration level on lumbar spine injuries in military populations.

BACKGROUND CONTEXT: Clinical studies have indicated that thoracolumbar trauma occurs in the civilian population at its junction. In contrast, injury patterns in military populations indicate a shift to the inferior vertebral levels of the lumbar spine. Controlled studies offering an explanation for such migrations and the associated clinical biomechanics are sparse in literature. PURPOSE: The goals of this study were to investigate the potential roles of acceleration loading on the production of injuries and their stability characteristics using a human cadaver model for applications to high-speed aircraft ejection and helicopter crashes. STUDY DESIGN: Biomechanical laboratory study using unembalmed human cadaver lumbar spinal columns. METHODS: Thoracolumbar columns from post-mortem human surrogates were procured, x-rays taken, intervertebral joints and bony components evaluated for degeneration, and fixed using polymethylmethacrylate. The inferior end was attached to a platform via a load cell and uniaxial accelerometer. The superior end was attached to the upper metal platform via a semi-circular cylinder. The pre-flexed specimen was preloaded to simulate torso mass. The ends of the platform were connected to the vertical post of a custom-designed drop tower. The specimen was dropped inducing acceleration loading to the column. Axial force and acceleration data were gathered at high sampling rates, filtered, and peak accelerations and inertia-compensated axial forces were obtained during the loading phase. Computed tomography images were used to identify and classify injuries using the three-column concept (stable vs. unstable trauma). RESULTS: The mean age, total body mass, and stature of the five healthy degeneration-free specimens were 42 years, 73 kg, and 167 cm. The first two specimens subjected to peak accelerations of approximately 200 m/sec(2) were classified as belonging to high-speed aircraft ejection-type and the other three specimens subjected to greater amplitudes (347-549 m/sec(2)) were classified as belonging to helicopter crash-type loadings. Peak axial forces for all specimens ranged from 4.8 to 7.2 kN. Ejection-type loaded specimens sustained single-level injuries to the L1 vertebra; one injury was stable and the other was unstable. Helicopter crash-type loaded specimens sustained injuries at inferior levels, including bilateral facet dislocation at L4-L5 and L2-L4 compression fractures, and all specimens were considered unstable at least at one spinal level. CONCLUSIONS: These findings suggest that the severity of spinal injuries increase with increasing acceleration levels and, more importantly, injuries shift inferiorly from the thoracolumbar junction to lower lumbar levels. Acknowledging that the geometry and load carrying capacity of vertebral bodies increase in the lower lumbar spine, involvement of inferior levels in trauma sparing the superior segments at greater acceleration inputs agree with military literature of caudal shift in injured levels. The present study offers an experimental explanation for the clinically observed caudal migration of spinal trauma in military populations as applied to high-speed aircraft ejection catapult and helicopter crashes.

Implementation and validation of probabilistic models of the anterior longitudinal ligament and posterior longitudinal ligament of the cervical spine.

The objective of this investigation was to develop probabilistic finite element (FE) models of the anterior longitudinal ligament (ALL) and posterior longitudinal ligament (PLL) of the cervical spine that incorporate the natural variability of biological specimens. In addition to the model development, a rigorous validation methodology was developed to quantify model performance. Experimental data for the geometry and dynamic properties of the ALL and PLL were used to create probabilistic FE models capable of predicting not only the mean dynamic relaxation response but also the observed experimental variation of that response. The probabilistic FE model uses a quasilinear viscoelastic material constitutive model to capture the time-dependent behaviour of the ligaments. The probabilistic analysis approach yields a statistical distribution for the model-predicted response at each time point rather than a single deterministic quantity (e.g. ligament force) and that response can be statistically compared to experimental data for validation. A quantitative metric that compares the cumulative distribution functions of the experimental data and model response is computed for both the ALL and PLL throughout the time histories and is used to quantify model performance.

A new PMHS model for lumbar spine injuries during vertical acceleration.

Ejection from military aircraft exerts substantial loads on the lumbar spine. Fractures remain common, although the overall survivability of the event has considerably increased over recent decades. The present study was performed to develop and validate a biomechanically accurate experimental model for the high vertical acceleration loading to the lumbar spine that occurs during the catapult phase of aircraft ejection. The model consisted of a vertical drop tower with two horizontal platforms attached to a monorail using low friction linear bearings. A total of four human cadaveric spine specimens (T12-L5) were tested. Each lumbar column was attached to the lower platform through a load cell. Weights were added to the upper platform to match the thorax, head-neck, and upper extremity mass of a 50th percentile male. Both platforms were raised to the drop height and released in unison. Deceleration characteristics of the lower platform were modulated by foam at the bottom of the drop tower. The upper platform applied compressive inertial loads to the top of the specimen during deceleration. All specimens demonstrated complex bending during ejection simulations, with the pattern dependent upon the anterior-posterior location of load application. The model demonstrated adequate inter-specimen kinematic repeatability on a spinal level-by-level basis under different subfailure loading scenarios. One specimen was then exposed to additional tests of increasing acceleration to induce identifiable injury and validate the model as an injury-producing system. Multiple noncontiguous vertebral fractures were obtained at an acceleration of 21 g with 488 g/s rate of onset. This clinically relevant trauma consisted of burst fracture at L1 and wedge fracture at L4. Compression of the vertebral body approached 60% during the failure test, with -6,106 N axial force and 168 Nm flexion moment. Future applications of this model include developing a better understanding of the vertebral injury mechanism during pilot ejection and developing tolerance limits for injuries sustained under a variety of different vertical acceleration scenarios.

Determination of normative neck muscle morphometry using upright MRI with comparison to supine data.

BACKGROUND: Neck muscles are important in the static and dynamic stability of the head-neck complex. Deep neck muscles act to maintain upright posture and superficial muscles are responsible for gross movements. Previous studies have quantified neck muscle geometry using traditional supine magnetic resonance imaging (MRI). However, supine orientation removes the vertical load on the cervical spine from the head-neck complex and changes the relative orientation of the spine and neck muscles. Therefore, the purpose of this study was to demonstrate the feasibility of upright MRI to obtain neck muscle morphometric data on a spinal level-by-level basis for subjects in upright seated positions. METHODS: Upright MRI scans were obtained of the neck region for six younger male volunteers in neutral and flexed positions. Planar images were oriented parallel to the intervertebral disc space at each level. Cross-sectional area (CSA) and orientation of neck muscles were quantified at four spinal levels. RESULTS: Area and position of all four muscles were significantly dependent upon spinal level. Average CSA of the sternocleidomastoid, longus colli, levator scapulae, and trapezius muscles in neutral position were 512, 113, 281, and 174 mm2. Head-neck position significantly affected area and position of the sternocleidomastoid and position of posterior neck muscles. DISCUSSION: Comparison of neck muscle areas from the present study to a previous study incorporating supine MRI demonstrated differing trends between anterior and posterior neck muscles that may be attributable to upright orientation of volunteers and planar image orientation in the present study. Differences between supine and upright MRI identified in the present study may warrant incorporation of this technique in future spinal imaging studies.

Kinematic response of the spine during simulated aircraft ejections.

INTRODUCTION: Military aviators are susceptible to spinal injuries during high-speed ejection scenarios. These injuries commonly arise as a result of strains induced by extreme flexion or compression of the spinal column. This study characterizes the vertebral motion of two postmortem human surrogates (PMHS) during a simulated catapult phase of ejection on a horizontal decelerator sled. METHODS: During testing, the PMHS were restrained supinely to a mock ejection seat and subjected to a horizontal deceleration profile directed along the local z-axis. Two midsized males (175.3 cm, 77.1 kg; 185.4 cm, 72.6 kg) were tested. High-rate motion capture equipment was used to measure the three-dimensional displacement of the head, vertebrae, and pelvis during the ejection event. RESULTS: The two PMHS showed generally similar kinematic motion. Head injury criterion (HIC) results were well below injury threshold levels for both specimens. The specimens both showed compression of the spine, with a reduction in length of 23.9 mm and 45.7 mm. Post-test autopsies revealed fractures in the C5, T1, and L1 vertebrae. DISCUSSION: This paper provides an analysis of spinal motion during an aircraft ejection.The injuries observed in the test subjects were consistent with those seen in epidemiological studies. Future studies should examine the effects of gender, muscle tensing, out-of-position (of head from neutral position) occupants, and external forces (e.g., windblast) on spinal kinematics during aircraft ejection.