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.

Ejection injury to the spine in small aviators: sled tests of manikins vs. post mortem specimens.

INTRODUCTION: This study presents the results of seven aerospace manikin and three post mortem human surrogate (PMHS) horizontal deceleration sled tests. The objective of this study was to establish a body of baseline data that examines the ability of small (fifth percentile) manikins to predict whole-body kinematics associated with aircraft ejection, and whether currently available head and neck injury criteria are applicable in these situations. METHODS: Subjects were exposed to a short-duration local z-axis sled pulse while horizontally seated and restrained in an ejection seat. Test subjects included instrumented fifth percentile female and male manikins, and two small (163.8 cm, 48.3 kg; 143.5 cm, 48.6 kg) female and one small (166.2 cm, 54.3 kg) male PMHS. RESULTS: The anterior (local x-axis) translations of the PMHS heads were less than those observed in the manikin tests, but the local z-axis translations of the PMHS heads were greater than those of the manikins. Z-axis translations of the manikins' T1 were generally similar to those of the PMHS T1, but the anterior x-axis translations of T1 were greater in the PMHS. The neck injury criterion (Nij) tended to under-predict observed injury (primarily ruptures of the posterior ligaments at C4-5, T2-3), and the Beam Criterion (BC) tended to over-predict observed injury for small occupants. The USN/USAF neck injury criteria (NIC) performed best in predicting the observed injuries. DISCUSSION: Present manikin designs do not predict the kinematics of PMHS in ejection tests. Further refinement of existing injury criteria is required to accurately predict location and severity of ejection-induced injuries.

Physical effects of ejection on the head-neck complex: demonstration of a cadaver model.

Vertebral fracture is the most common severe injury during high-speed pilot ejection. However, the loading paradigm experienced by pilots may also lead to soft-tissue spinal injuries that are more difficult to quantify and can lead to long-term deficits. This manuscript describes a new experimental protocol to simulate the effects of pilot ejection on the tissues of the head-neck complex. The model permits precise control of head-neck complex initial positioning, detailed analysis of head and spinal kinematics and upper and lower neck loads, and the ability to thoroughly investigate and identify soft-tissue injuries through upper and lower neck injury criteria, radiography, manual palpation, and cryomicrotomy. For the current test, peak acceleration of +14.8 Gz was similar to actual ejection events and duration of the acceleration pulse was approximately 100 ms. The specimen was oriented in flexion prior to initiation of inferior-to-superiorly directed acceleration. Subfailure ligamentum flavum injuries were sustained at the C4-C5 and C5-C6 cervical spinal levels and identified by increased segmental motions during the simulated ejection, increased laxity following testing, and cryomicrotomy. Upper and lower neck injury criteria did not predict these soft-tissue injuries. This experimental model can be used for detailed analysis of the effects of gender, head-neck orientation, helmet instrumentation, and acceleration pulse characteristics on cervical spine injury potential during pilot ejection events.

Gender dependent cervical spine anatomical differences in size-matched volunteers - biomed 2009.

The objective was to examine significant differences in the bony structure of cervical spine vertebrae based on gender and spinal level that may influence injury risk in women following automotive rear impact. Male and female subjects were recruited for a separate study and data from two subsets were selected for inclusion in this study. Subjects were size-matched based on sitting height (17 males, 11 females) and head circumference (9 males, 18 females). Axial CT scans were obtained of the cervical spine from the C1 through C6. Bony boundaries of cervical vertebrae were defined using image-analysis software and biomechanically-relevant dimensions were derived at spinal levels C2 through C6. Six of seven vertebral dimensions were significantly dependent upon gender and spinal level in both subgroups. Male vertebrae had larger dimensions for each metric. Depth dimensions were greatest at caudal and cranial extents, whereas width dimensions were smallest at C2 and increased caudally. Greater linear and areal dimensions in size-matched male subjects indicates a more stable cervical spinal column that may be more capable of resisting inertial loading of the head-neck complex during automotive rear impacts. Although the explanation for greater injury susceptibility in females is likely multi-factorial, including differences in spinal material properties, soft tissue tolerance thresholds, occupant-seatback orientation, and neck muscle size/orientations, the present study has identified significant differences in cervical spine anatomical dimensions that may contribute to greater rates of whiplash injury in that population.