Finite element model of the human neck during omni-directional impacts. Part II: relation between cervical curvature and risk of injury

A detailed 3D FE model of the human neck was used to assess a possible relationship between risk of injury and cervical spine curvature for various impacts. A FE model was previously developed, representing the head and neck of a 50th percentile human with a normal lordotic curvature. The model behaviour was omni-directionally validated for various impacts using published results. For the present study, the model was deformed in order to obtain a straight and a kyphotic curvature, and for each geometry, rear-end, frontal, lateral and oblique impact were simulated. Although results showed similar kinematic patterns, significant differences were found in the distribution and peak values of ligament elongations, forces and moments along the cervical spine for the three configurations. It was concluded that the variability observed on the curvature of the human cervical spine may have a significant influence both on the behaviour and on the risk of injury of the neck during impact.     

A female head–neck model for rear impact simulations

Several mathematical cervical models of the 50th percentile male have been developed and used for impact biomechanics research. However, for the 50th percentile female no similar modelling efforts have been made, despite females being subject to a higher risk of soft tissue neck injuries. This is a limitation for the development of automotive protective systems addressing Whiplash Associated Disorders (WADs), most commonly caused in rear impacts, as the risk for females sustaining WAD symptoms is double that of males.In this study, a finite element head and neck model of a 50th percentile female was validated in rear impacts. A previously validated ligamentous cervical spine model was complemented with a rigid body head, soft tissues and muscles. In both physiological flexion-extension motions and simulated rear impacts, the kinematic response at segment level was comparable to that of human subjects. Evaluation of ligament stress levels in simulations with varied initial cervical curvature revealed that if an individual assumes a more lordotic posture than the neutral, a higher risk of WAD might occur in rear impact.The female head and neck model, together with a kinematical whole body model which is under development, addresses a need for tools for assessment of automotive protection systems for the group which is at the highest risk to sustain WAD.     

Analysis of control strategies for VIVA OpenHBM with active reflexive neck muscles

Modeling muscle activity in the neck muscles of a finite element (FE) human body model can be based on two biological reflex systems. One approach is to approximate the Vestibulocollic reflex (VCR) function, which maintains the head orientation relative to a fixed reference in space. The second system tries to maintain the head posture relative to the torso, similar to the Cervicocolic reflex (CCR). Strategies to combine these two neck muscle controller approaches in a single head-neck FE model were tested, optimized, and compared to rear-impact volunteer data. The first approach, Combined-Control, assumed that both controllers simultaneously controlled all neck muscle activations. In the second approach, Distributed-Control, one controller was used to regulate activation of the superficial muscles while a different controller acted on deep neck muscles. The results showed that any muscle controller that combined the two approaches was less effective than only using one of VCR- or CCR-based systems on its own. A passive model had the best objective rating for cervical spine kinematics, but the addition of a single active controller provided the best response for both head and cervical spine kinematics. The present study demonstrates the difficulty in completely capturing representative head and cervical spine responses to rear-impact loading and identified a controller capturing the VCR reflex as the best candidate to investigate whiplash injury mechanisms through FE modeling.

     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

Investigation of whiplash injuries in the upper cervical spine using a detailed neck model

Whiplash injuries continue to have significant societal cost; however, the mechanism and location of whiplash injury is still under investigation. Recently, the upper cervical spine ligaments, particularly the alar ligament, have been identified as a potential whiplash injury location. In this study, a detailed and validated explicit finite element model of a 50th percentile male cervical spine in a seated posture was used to investigate upper cervical spine response and the potential for whiplash injury resulting from vehicle crash scenarios. This model was previously validated at the segment and whole spine levels for both kinematics and soft tissue strains in frontal and rear impact scenarios. The model predicted increasing upper cervical spine ligament strain with increasing impact severity. Considering all upper cervical spine ligaments, the distractions in the apical and alar ligaments were the largest relative to their failure strains, in agreement with the clinical findings. The model predicted the potential for injury to the apical ligament for 15.2 g frontal or 11.7 g rear impacts, and to the alar ligament for a 20.7 g frontal or 14.4 g rear impact based on the ligament distractions. Future studies should consider the effect of initial occupant position on ligament distraction.     

The influence of surface padding properties on head and neck injury risk

A validated computational head-neck model was used to understand the mechanical relationships between surface padding characteristics and injury risk during impacts near the head vertex. The study demonstrated that injury risk can be decreased by maximizing the energy-dissipating ability of the pad, choosing a pad stiffness that maximizes pad deformation without bottoming out, maximizing pad thickness, and minimizing surface friction. That increasing pad thickness protected the head without increasing neck loads suggests that the increased cervical spine injury incidence previously observed in cadaveric impacts to padded surfaces relative to lubricated rigid surfaces was due to increased surface friction rather than pocketing of the head in the pad.     

Human shoulder response to side impacts: a finite element study

This study aimed at developing a shoulder finite element (FE) model able to simulate the dynamic behaviour and to predict injuries in case of side impacts. This model is an updated version of the initial Human Model for Safety (HUMOS) FE model of the human body. Simulations performed with the model have been compared to experimental results of side impact tests conducted previously at INRETS. The shoulder model response under lateral impact appears to be in good agreement with experimental data such as impact force and shoulder deflections for different impact speeds and impact directions. These results seem promising for future applications such as shoulder injury prediction in simulated car crashes.     

An automated method to morph finite element whole-body human models with a wide range of stature and body shape for both men and women

Field data analyses have shown that small female, obese, and/or older occupants are at increased risks of death and serious injury in motor-vehicle crashes compared with mid-size young men. The current adult finite element (FE) human models represent occupants in the same three body sizes (large male, mid-size male, and small female) as those for the contemporary adult crash dummies. Further, the time needed to develop an FE human model using the traditional method is measured in months or even years. In the current study, an improved regional mesh morphing method based on landmark-based radial basis function (RBF) interpolation was developed to rapidly morph a mid-size male FE human model into different geometry targets. A total of 100 human models with a wide range of human attributes were generated. A pendulum chest impact condition was applied to each model as an initial assessment of the resulting variability in response. The morphed models demonstrated mesh quality similar to the baseline model. The peak impact forces and chest deflections in the chest pendulum impacts varied substantially with different models, supportive of consideration of population variation in evaluating the occupant injury risks. The method developed in this study will enable future safety design optimizations targeting at various vulnerable populations that cannot be considered with the current models.     

Development and validation of a CO-C7 FE complex for biomechanical study

In this study, the digitized geometrical data of the embalmed skull and vertebrae (C0-C7) of a 68-year old male cadaver were processed to develop a comprehensive, geometrically accurate, nonlinear C0-C7 FE model. The biomechanical response of human neck under physiological static loadings, near vertex drop impact and rear-end impact (whiplash) conditions were investigated and compared with published experimental results. Under static loading conditions, the predicted moment-rotation relationships of each motion segment under moments in midsagittal plane and horizontal plane agreed well with experimental data. In addition, the respective predicted head impact force history and the S-shaped kinematics responses of head-neck complex under near-vertex drop impact and rear-end conditions were close to those observed in reported experiments. Although the predicted responses of the head-neck complex under any specific condition cannot perfectly match the experimental observations, the model reasonably reflected the rotation distributions among the motion segments under static moments and basic responses of head and neck under dynamic loadings. The current model may offer potentials to effectively reflect the behavior of human cervical spine suitable for further biomechanics and traumatic studies.     

The effect of axial load on the sagittal plane curvature of the upright human spine in vivo

Determining the effect of load carriage on the human spine in vivo is important for determining spinal forces and establishing potential mechanisms of back injury. Previous studies have suggested that the natural curvature of the spine straightens under load, but are based on modelling and external measurements from the surface of the back. In the current study, an upright positional MRI scanner was used to acquire sagittal images of the lumbar and lower thoracic spine of 24 subjects. The subjects were imaged in standing whilst supporting 0, 8 and 16kg of load which was applied axially across the shoulders using an apron. An active shape model of the vertebral bodies from T10 to S1 was created and used to characterise the effect of load. The results from the shape model showed that the behaviour of the average-shaped spine was to straighten slightly. However, the shape model also showed that the effect of load exhibited systematic variation between individuals. Those who had a smaller than average curvature before loading straightened under load, whereas those who had a greater than average curvature before loading showed an increase in curvature under load. The variation in behaviour of differently shaped spines may have further implications for the effects of load in lifting manoeuvres and in understanding the aetiology of back pain.     

Developmental biomechanics of the cervical spine: Tension and compression

Epidemiological data and clinical indicia reveal devastating consequences associated with pediatric neck injuries. Unfortunately, neither injury prevention nor clinical management strategies will be able to effectively reduce these injuries or their effects on children, without an understanding of the cervical spine developmental biomechanics. Thus, we investigated the relationship between spinal development and the functional (stiffness) and failure biomechanical characteristics of the cervical spine in a baboon model. A correlation study design was used to define the relationships between spinal tissue maturation and spinal biomechanics in both tension and compression. Eighteen baboon cervical spine specimens distributed across the developmental spectrum (1–26 human equivalent years) were dissected into osteoligamentous functional spinal units. Using a servo-hydraulic MTS, these specimens (Oc–C2, C3–C4, C5–C6, C7–T1) were non-destructively tested in tension and compression and then displaced to failure in tension while measuring the six-axes of loads and displacements. The functions describing the developmental biomechanical response of the cervical spine for stiffness and normalized stiffness exhibited a significant direct relationship in both tension and compression loading. Similarly, the tensile failure load and normalized failure load demonstrated significant maturational increases. Further, differences in biomechanical response were observed between the spinal levels examined and all levels exhibited clinically relevant failure patterns. These data support our understanding of the child cervical spine from a developmental biomechanics perspective and facilitate the development of injury prevention or management schema for the mitigation of child spine injuries and their deleterious effects.     

A biomechanical evaluation of whiplash using a multi-body dynamic model

This paper presents a biomechanical evaluation of whiplash injury potential during the initial extension motion of the head in a rear-end collision. A four-segment dynamic model is developed in the sagittal plane for the analysis. The model response is validated using the existing experimental data and is shown to simulate the "S-shape" kinematics of the cervical spine and the resulting dynamics observed in human and cadaver experiments. The model is then used to evaluate the effects of parameters such as collision severity, head/headrest separation, and the initial head orientation in the sagittal plane on the "S-shape" kinematics of the cervical spine and the resulting neck loads. It is shown, for example, that the cervical spine forms an "S-shape" for a range of change in speeds and that at lower and higher speeds changes the spine does not form the "S-shape." Furthermore, it is shown that the "S-shape" formation also depends on the head to headrest separation distance.     

Inertial properties and loading rates affect buckling modes and injury mechanisms in the cervical spine

Cervical spine injuries continue to be a costly societal problem. Future advancements in injury prevention depend on improved physical and computational models which, in turn, are predicated on a better understanding of the responses of the neck during dynamic loading. Previous studies have shown that the tolerance of the neck is dependent on its initial position and its buckling behavior. This study uses a computational model to examine the mechanical factors influencing buckling behavior during impact to the neck. It was hypothesized that the inertial properties of the cervical spine influence the dynamics during compressive axial loading. The hypothesis was tested by performing parametric analyses of vertebral mass, mass moments of inertia, motion segment stiffness, and loading rate. Increases in vertebral mass resulted in increasingly complex kinematics and larger peak loads and impulses. Similar results were observed for increases in stiffness. Faster loading rates were associated with higher peak loads and higher-order buckling modes. The results demonstrate that mass has a great deal of influence on the buckling behavior of the neck, particularly with respect to the expression of higher-order modes. Injury types and mechanisms may be substantially altered by loading rate because inertial effects may influence whether the cervical spine fails in a compressive mode, or a bending mode.     

An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.     

Flexion and extension structural properties and strengths for male cervical spine segments

New vehicle safety standards are designed to limit the amount of neck tension and extension seen by out-of-position motor vehicle occupants during airbag deployments. The criteria used to assess airbag injury risk are currently based on volunteer data and animal studies due to a lack of bending tolerance data for the adult cervical spine. This study provides quantitative data on the flexion-extension bending properties and strength on the male cervical spine, and tests the hypothesis that the male is stronger than the female in pure bending. An additional objective is to determine if there are significant differences in stiffness and strength between the male upper and lower cervical spine. Pure-moment flexibility and failure testing was conducted on 41 male spinal segments (O-C2, C4-C5, C6-C7) in a pure-moment test frame and the results were compared with a previous study of females. Failures were conducted at approximately 90 N-m/s. In extension, the male upper cervical spine (O-C2) fails at a moment of 49.5 (s.d. 17.6)N-m and at an angle of 42.4 degrees (s.d. 8.0 degrees). In flexion, the mean moment at failure is 39.0 (s.d. 6.3 degrees) N-m and an angle of 58.7 degrees (s.d. 5.1 degrees). The difference in strength between flexion and extension is not statistically significant. The difference in the angles is statistically significant. The upper cervical spine was significantly stronger than the lower cervical spine in both flexion and extension. The male upper cervical spine was significantly stiffer than the female and significantly stronger than the female in flexion. Odontoid fractures were the most common injury produced in extension, suggesting a tensile mechanism due to tensile loads in the odontoid ligamentous complex.     

Quantitative comparison of ligament formulation and pre-strain in finite element analysis of the human lumbar spine

Data has been published that quantifies the nonlinear, anisotropic material behaviour and pre-strain behaviour of the anterior longitudinal, supraspinous (SSL), and interspinous ligaments of the human lumbar spine. Additionally, data has been published on localized material properties of the SSL. These results have been incrementally incorporated into a previously validated finite element model of the human lumbar spine. Results suggest that the effects of increased ligament model fidelity on bone strain energy were moderate and the effects on disc pressure were slight, and do not justify a change in modelling strategy for most clinical applications. There were significant effects on the ligament stresses of the ligaments that were directly modified, suggesting that these phenomena should be included in FE models where ligament stresses are the desired metric.     

A three-dimensional model of the human cervical spine for impact simulation

A three-dimensional analytical model of the cervical spine is described. The cervical vertebrae and the head are modeled as rigid bodies which are interconnected by deformable elements representing the intervertebral disks, facet joints, ligaments and muscles. A special pentahedral continuum element for representing the articular facets is described which effectively maintains stability of the cervical spine in both lateral and frontal plane accelerations, which is very difficult with multi-spring models of the facets. A simplified representation is used for the spine and body below the level of T1. The neck musculature is modeled by over 100 muscle elements representing 22 major muscle groups in the neck. The model has been validated for frontal and sideways impact accelerations by simulating published experimental data. Results are also presented to show the effects of the stretch reflex response on the dynamics of the head and neck under moderate acceleration.     

Brain tissue analysis of impacts to American football helmets

Concussion in American football is a prevalent concern. Research has been conducted examining frequencies, location, and thresholds for concussion from impacts. Little work has been done examining how impact location may affect risk of concussive injury. The purpose of this research was to examine how impact site on the helmet and type of impact, affects the risk of concussive injury as quantified using finite element modelling of the human head and brain. A linear impactor was used to impact a helmeted Hybrid III headform in several locations and using centric and non-centric impact vectors. The resulting dynamic response was used as input for the Wayne State Brain Injury Model to determine the risk of concussive injury by utilizing maximum principal strain as the predictive variable. The results demonstrated that impacts that occur primarily to the side of the head resulted in higher magnitudes of strain in the grey and white matter, as well as the brain stem. Finally, commonly worn American football helmets were used in this research and significant risk of injury was incurred for all impacts. These results suggest that improvements in American football helmets are warranted, in particular for impacts to the side of the helmet.     

The neutral posture of the cervical spine is not unique in human subjects

Cervical spine injuries often happen in dynamic environments (e.g., sports and motor vehicle crashes) where individuals may be moving their head and neck immediately prior to impact. This motion may reposition the cervical vertebrae in a way that is dissimilar to the upright resting posture that is often used as the initial position in cadaveric studies of catastrophic neck injury. Therefore our aim was to compare the “neutral” cervical alignment measured using fluoroscopy of 11 human subjects while resting in a neutral posture and as their neck passed through neutral during the four combinations of active flexion and extension movements in both an upright and inverted posture. Muscle activation patterns were also measured unilaterally using surface and indwelling electromyography in 8 muscles and then compared between the different conditions. Overall, the head posture, cervical spine alignment and muscle activation levels were significantly different while moving compared to resting upright. Compared to the resting upright condition, average head postures were 6–13° more extended, average vertebral angles varied from 11° more extended to 10° more flexed, and average muscle activation levels varied from unchanged to 10% MVC more active, although the exact differences varied with both direction of motion and orientation. These findings are important for ex vivo testing where the head and neck are statically positioned prior to impact – often in an upright neutral posture with negligible muscle forces – and suggest that current cadaveric head-first impact tests may not reflect many dynamic injury environments.     

A new acceleration apparatus for the study of whiplash with human cadaveric cervical spine specimens

The biomechanics of whiplash is often studied using cadaveric cervical spine specimens. One of the most important points in this kind of study is to create realistic loading conditions. The aim of the present project therefore was to develop an acceleration apparatus, which allows the study of whiplash with human cadaveric cervical spine specimens under as realistic loading conditions as possible. The new acceleration apparatus mainly consisted of a sled, a pneumatic acceleration unit and a railtrack and offered several unique features to create more realistic loading conditions. Among these features, the possibility to simulate the passive movements of the trunk is of capital importance. In this new apparatus, first, the general feasibility of whiplash experiments was studied, second, the reproducibility of the impacts was quantified and third, the effect of simulated movements of the trunk on accelerations and loads was examined. In the new acceleration apparatus various types of collisions could reproducibly be simulated. Simulated passive movements of the trunk strongly influenced the loading pattern of the neck. Without pivoting a steep increase of all loading parameters could be observed. This increase was less pronounced if pivoting was allowed. In conclusion, biomechanical aspects of whiplash could reproducibly be examined in the new acceleration apparatus. Due to its significant effects on the loading of the neck, pivoting of the trunk should always be taken into account in future experiments on the biomechanics of whiplash in which isolated cervical spine specimens are used.