An omni-directional model of injury risk in planar crashes with application for autonomous vehicles.

Objective: Automated driving systems (ADS) are actively being deployed within the driving fleet. ADS are designed to safely navigate roadways, which entails an expectation of encountering varying degrees of potential conflict with other road users. The ADS design and evaluation process benefits from estimating injury severity probabilities for collisions that may occur. Current regression models in the literature are typically bespoke analyses involving targeted principal directions of force (PDOFs) and occupant positions. It is preferable to rely on injury severity models derived from a single source to provide a continuous function of risk for all planar collisions, while also accounting for specific vehicle and occupant characteristics. The novel feature of the proposed models is continuous, parametric injury risk surfaces that encompass the full spectrum of available United States field data.Methods: We used years 2001-2015 of the National Automotive Sampling System, Crashworthiness Data System (NASS-CDS) and years 2017-2019 of the Crash Investigation Sampling System (CISS) to estimate injury risk at the maximum abbreviated injury scale (MAIS) 3 and higher (3+) and 5 and higher (5+) levels for all adult occupants traveling in 2002 or newer passenger vehicles which were less than 10 years old at the time of the crash. The models account for occupant, vehicle, and crash characteristics. Interactions with vulnerable road users (e.g., pedestrian, bicyclist) were not considered.Results: We present statistical models suitable to predict injury in all non-rollover crashes at the maximum MAIS3+ and 5+ levels, and show that these models can be comparable to similar single scenario (e.g., frontal) crash models. We discuss challenges with imputing missing field data, and discuss handling of covariates that may not be known at the time of the crash.Conclusions: Collision severity assessment is a vital component of the ADS design process. We developed a novel injury risk function that can assess occupant injury risks across the spectrum of foreseeable planar collisions. These models can provide insight on potential outcomes of counterfactual simulations, injury risk and crashworthiness considerations for human driven vehicles, and provide an evaluation tool that can be applied in ADS safety impact evaluation.

Comparison of Laboratory and On-Field Performance of American Football Helmets.

The relationship between laboratory and on-field performance of football helmets was assessed for 31 football helmet models selected from those worn by players in the 2015-2019 National Football League (NFL) seasons. Linear impactor tests were conducted with helmets placed on an instrumented Hybrid III head and neck assembly mounted on a sliding table. Based on impacts to each helmet at six impact locations and three velocities, a helmet performance score (HPS) was calculated using a linear combination of the head injury criterion (HIC) and the diffuse axonal multi-axis general evaluation (DAMAGE). To determine the on-field performance of helmets, helmet model usage, player participation, and incident concussion data were collected from the five NFL seasons and used to calculate helmet model-specific concussion rates. Comparison of laboratory HPS to the helmet model-specific concussion rates on a per play basis showed a positive correlation (r2 = 0.61, p < 0.001) between laboratory and on-field performance of helmet models, indicating that helmets which exhibited reduced impact severity in the laboratory tests were also generally associated with lower concussion rates on-field. Further analysis showed that NFL-prohibited helmet models exhibited a significantly higher odds of concussion (OR 1.24; 95% CI 1.04-1.47; p = 0.017) relative to other helmet models.

Position-Specific Circumstances of Concussions in the NFL: Toward the Development of Position-Specific Helmets.

Consideration of position-specific features of the NFL concussion environment could enable improved risk mitigation through the design of position-specific helmets to improve self-protection as well as protection for the other player with whom the contact occurs. The purpose of this paper is to quantify position-specific features of scenarios resulting in concussions to NFL players, and the players they contact, by reviewing all game footage (broadcast and non-broadcast) over 4 seasons. Position-specific features were documented for 647 concussions in which a primary exposure could be visualized, including impact source, helmet impact location, activity, and the other player with whom the contact occurred. Findings include the over-representation of helmet-to-ground impacts to the rear of the quarterback's helmet, the high frequency of impacts to the side (upper) location of both concussed players and the players they contacted regardless of position, and distinct differences in the circumstances of concussions to cornerbacks and safeties. The study shows that some features of concussion scenarios are common to all positions, but several position-specific features exist and can inform the design of position-specific helmets for NFL players.

Video Analysis of Reported Concussion Events in the National Football League During the 2015-2016 and 2016-2017 Seasons.

BACKGROUND: Concussions in American football remain a high priority of sports injury prevention programs. Detailed video review provides important information on causation, the outcomes of rule changes, and guidance on future injury prevention strategies. PURPOSE: Documentation of concussions sustained in National Football League games played during the 2015-2016 and 2016-2017 seasons, including consideration of video views unavailable to the public. STUDY DESIGN: Descriptive epidemiology study. METHODS: All reported concussions were reviewed with all available video footage. Standardized terminology and associated definitions were developed to describe and categorize the details of each concussion. RESULTS: Cornerbacks sustained the most concussions, followed by wide receivers, then linebackers and offensive linemen. Half (50%) of concussions occurred during a passing play, 28% during a rushing play, and 21% on a punt or kickoff. Tackling was found to be the most common activity of concussed players, with the side of the helmet the most common helmet impact location. The distribution of helmet impact source-the object that contacted the concussed player's helmet-differed from studies of earlier seasons, with a higher proportion of helmet-to-body impacts (particularly shoulder) and helmet-to-ground impacts and with a lower proportion of helmet-to-helmet impacts. Helmet-to-ground concussive impacts were notable for the high prevalence of impacts to the back of the helmet and their frequency during passing plays. CONCLUSION: Concussion causation scenarios in the National Football League have changed over time. CLINICAL RELEVANCE: The results of this study suggest the need for expanded evaluation of concussion countermeasures beyond solely helmet-to-helmet test systems, including consideration of impacts with the ground and with the body of the opposing player. It also suggests the possibility of position-specific countermeasures as part of an ongoing effort to improve safety.

Kinematics and kinetics of vigorous head shaking.

Previous studies on neck muscle strength and motion have assumed or imposed varying constraints on the heads and bodies of the subjects. In this study, we asked 20 subjects to vigorously shake their heads 5-10 times in a completely unconstrained manner. The kinematics and kinetics of the head and neck were measured from video analysis and instrumentation mounted inside the mouth. Subjects shook their heads at self-selected tempos ranging from 1.9-4.7 Hz over a 20-91° range of motion. The motion of each subject's head could be approximated by a fixed center of rotation that was typically located in the midcervical spine, but varied widely among subjects. Significant differences between men and women were observed. Peak head accelerations were low (4.3 ± 1.1 g and 250 ± 103 rad/s2 for men, 3.0 ± 0.9 g and 182 ± 58 rad/s2 for women) and estimated peak generated neck moments at C7/T1 were comparable to values reported in isometric neck strength studies (47 ± 14 N·m in extension and 22 ± 9 N·m in flexion for men, 25 ± 8 N·m in extension and 9 ± 7 N·m in flexion for women).

Factors affecting ejection risk in rollover crashes.

Ejection greatly increases the risk of injury and fatality in a rollover crash. The purpose of this study was to determine the crash, vehicle, and occupant characteristics that affect the risk of ejection in rollovers. Information from real world rollover crashes occurring from 2000 - 2010 was obtained from the National Automotive Sampling System (NASS) in order to analyze the effect of the following parameters on ejection risk: seatbelt use, rollover severity, vehicle type, seating position, roof crush, side curtain airbag deployment, glazing type, and occupant age, gender, and size. Seatbelt use was found to reduce the risk of partial ejection and virtually eliminate the risk of complete ejection. For belted occupants, the risk of partial ejection risk was significantly increased in rollover crashes involving more roof inversions, light trucks and vans (LTVs), and larger occupants. For unbelted occupants, the risk of complete ejection was significantly increased in rollover crashes involving more roof inversions, LTVs, far side occupants, and higher levels of roof crush. Roof crush was not a significant predictor of ejection after normalizing for rollover severity. Curtain airbag deployment was associated with reduced rates of partial and complete ejection, but the effect was not statistically significant, perhaps due to the small sample size (n = 89 raw cases with curtain deployments). A much greater proportion of occupants who were ejected in spite of curtain airbag deployment passed through the sunroof and other portals as opposed to the adjacent side window compared to occupants who were ejected in rollovers without a curtain airbag deployment. The primary factors that reduce ejection risk in rollover crashes are, in generally decreasing order of importance: seatbelt use, fewer roof inversions, passenger car body type, curtain airbag deployment, near side seating position, and small occupant size.

Comparison of risk factors for cervical spine, head, serious, and fatal injury in rollover crashes.

Previous epidemiological studies of rollover crashes have focused primarily on serious and fatal injuries in general, while rollover crash testing has focused almost exclusively on cervical spine injury. The purpose of this study was to examine and compare the risk factors for cervical spine, head, serious, and fatal injury in real world rollover crashes. Rollover crashes from 1995-2008 in the National Automotive Sampling System-Crashworthiness Data System (NASS-CDS) were investigated. A large data set of 6015 raw cases (2.5 million weighted) was generated. Nonparametric univariate analyses, univariate logistic regression, and multivariate logistic regression were conducted. Complete or partial ejection, a lack of seatbelt use, a greater number of roof inversions, and older occupant age significantly increased the risk of all types of injuries studied (p<0.05). Far side seating position increased the risk of fatal, head, and cervical spine injury (p<0.05), but not serious injury in general. Higher BMI was associated with an increased risk of fatal, serious, and cervical spine injury (p<0.05), but not head injury. Greater roof crush was associated with a higher rate of fatal and cervical spine injury (p<0.05). Vehicle type, occupant height, and occupant gender had inconsistent and generally non-significant effects on injury. This study demonstrates both common and unique risk factors for different types of injuries in rollover crashes.

Head and neck loading in everyday and vigorous activities.

The purpose of this study was to document head and neck loading in a group of ordinary people engaged in non-injurious everyday and more vigorous physical activities. Twenty (20) volunteers that were representative of the general population were subjected to seven test scenarios: a soccer ball impact to the forehead, a self-imposed hand strike to the forehead, vigorous head shaking, plopping down in a chair, jumping off a step, a seated drop onto the buttocks, and a vertical drop while seated supine in a chair. Some scenarios involved prescribed and well-controlled stimuli, while others allowed the volunteers to perform common activities at a self-selected level of intensity. Head accelerations up to 31 g and 2888 rad/s(2) and neck loads up to 268 N in posterior shear, 526 N in compression, and 36 Nm in extension were recorded. Most head and neck injury criteria predicted a low risk of injury in all activities. However, rotational head accelerations and Neck Injury Criterion (NIC) values were much higher than some proposed tolerance limits in a large number of tests, all of which were non-injurious. The data from this study help us to establish an envelope of head and neck loading that is commonly encountered and presents a minimal risk of injury.

Effect of delta-V errors in NASS on frontal crash risk calculations.

The most important factor in predicting the risk of injury or death in a frontal crash is the crash severity, which is expressed as the velocity change, or delta-V, experienced by the vehicle during the crash. The National Automotive Sampling System (NASS) is the largest database in the world linking injury outcomes with delta-Vs, which are obtained from field reconstructions. The accuracy of these reconstructions was assessed by analyzing 228 NASS cases involving single event frontal crashes in which the vehicle's frontal delta-V was also measured directly by an onboard event data recorder (EDR). Compared to the EDR measurements, the delta-V values in NASS averaged 19% lower with a standard deviation of 8.6 kph. The effect of this error on injury and fatality risk calculations was investigated using NASS data from 1997 - 2006 for frontal crashes with a known delta-V. Injury and fatality risk functions were calculated by curve fitting the distributions of the delta-V values associated with injury and fatality incidence normalized by the fitted crash exposure distribution. Individual delta-V values were linearly scaled to correct for the bias error, and the delta-V distributions were corrected for scatter error using a numerical deconvolution technique. Correcting for delta-V bias error shifted the calculated risk curves to the right and correcting for delta-V scatter error shifted the curves back to the left, but to a lesser extent. The effects of occupant age, gender, and belt use on injury and fatality risk were substantial.

The influence of body mass index on thoracic injuries in frontal impacts.

For this study, a comprehensive analysis was performed to assess the influence of body mass index on thoracic injury potential. The data for this study were obtained from the National Automotive Sampling System-Crashworthiness Data System (NASS-CDS) database for years 1993-2005. Obese occupants had a 26 and 33% higher risk of AIS > or = 2 and AIS > or = 3 thoracic injury when compared to lean occupants. The increased risk of AIS > or = 3 injury due to obesity was slightly higher for older occupants, but the influence of age was greater than that of obesity. The increase in injury potential was higher for unbelted obese occupants than unbelted. Non-parametric and parametric risk curves were developed to estimate the risk of thoracic injury based on occupant BMI, belt use and delta-V. Overall, increase in thoracic injury risk due to obesity is more prominent in males and older occupants and for occupants sustaining AIS > or = 3 thoracic injuries.

Forearm fracture bending risk functin for the 50th percentile male.

The increase in upper extremity injuries in automobile collisions, because of the widespread implantation of airbags, has lead to a better understanding of forearm injury criteria. Risk functions for upper extremity injury that can be used in instrumented upper extremities would be useful. This paper presents a risk function for forearm injury for the 50th percentile male based on bending fracture moment data gathered from previous studies. The data was scaled using two scaling factors, one for orientation and one for mass, and the Weibull survival analysis model was then used to develop the risk function. It was determined that a 25% risk of injury corresponds to an 82 Nm bending load, a 50% risk of injury corresponds to a 100 Nm bending load, and a 75% risk of injury corresponds to a 117 Nm bending load. It is believed the risk function can be used with an instrumented upper extremity during vehicle testing.

Relationship between linear and rotational head acceleration in various activities.

Head injury is typically predicted using linear or rotational acceleration-based injury criteria. In many cases, the linear components of head acceleration can be determined more easily than the rotational components. Peak rotational head acceleration (apeak) can be calculated from the peak linear head acceleration (apeak) by assuming a value for the effective radius of rotation (r) of the head (apeak = apeak / r). Empirical values for the effective radius of rotation were calculated using linear and angular head acceleration data from 20 human volunteers subjected to a wide variety of test scenarios, including a soccer ball impact to the forehead, a voluntary hand strike to the forehead, voluntary shaking of the head, plopping down in a chair, and a vertical drop while seated supine in a chair. In addition, values for the effective radius of rotation of the head were calculated for American football, boxing, and frontal and rear end automotive impacts using data from the literature. The range of values for the effective radius of rotation of the head for each activity was characterized statistically by calculating median, middle 50%, and middle 95% values from the cumulative distribution. Median values for the effective radius of rotation of the head ranged from 84 mm for boxing to 376 mm for voluntary hand strikes to the forehead. It is important to note that the concept of an effective radius of rotation of the head is simply a convenient method for expressing an empirical relationship between peak linear and peak rotational head acceleration, and does not represent an accurate model of the kinematics of the head. In most of the activities studied, the head kinematics were complex, and the center of gravity of the head did not rotate at a constant radius about a fixed point.

Humerus fracture bending risk function for the 50th percentile male.

The increase in upper extremity injuries in automobile collisions, because of the widespread implantation of airbags, has lead to an increased focus in humerus injury criteria. Risk functions for upper extremity injury that can be used in instrumented upper extremities would be useful. This paper presents a risk function for humerus injury for the 50th percentile male based on bending fracture moment data gathered from previous studies. The data was scaled using two scaling factors, one for mass and one for rate, and the Weibull survival analysis model was then used to develop the risk function. It was determined that a 25% risk of injury corresponds to a 214 Nm bending load, a 50% risk of injury corresponds to a 257 Nm bending load, and a 75% risk of injury corresponds to a 296 Nm bending load. It is believed the risk function can be used with an instrumented upper extremity during vehicle testing.

Predicting zygoma fractures from baseball impact.

The purpose of this study is to develop injury risk functions that predict zygoma fracture based on baseball type and impact velocity. Zygoma fracture strength data from published experiments were mapped with the force exerted by a baseball on the orbit as a function of ball velocity. Using a normal distribution, zygoma fracture risk functions were developed. Experimental evaluation of these risk functions was performed using six human cadaver tests and two baseballs of different stiffness values. High speed video measured the baseball impact velocity. Post test analysis of the cadaver skulls was performed using CT imaging including three-dimensional reconstruction as well as autopsy. The developed injury risk functions accurately identify the risk of zygoma fracture as a result of baseball impact. The experimental results validated the zygoma risk functions at the lower and upper levels. The injuries observed in the post test analysis included fractures of the zygomatic arch, frontal process and the maxilla, zygoma suture, with combinations of these creating comminuted, tripod fractures of the zygoma. Tests with a softer baseball did result in injury but these had fewer resulting zygoma bone fragments and occurred at velocities 50% higher than the major league ball.

Regional variation in the structural response and geometrical properties of human ribs.

By incorporating material and geometrical properties into a model of the human thorax one can develop an injury criterion that is a function of stress and strain of the material and not a function of the global response of the thorax. Previous research on the mechanical properties of ribs has focused on a limited set of specific ribs. For this study a total of 52 rib specimens were removed from four cadaver subjects. Variation in peak moment by thoracic region was significant (p < 0.01) with average values of 2, 2.9 and 3.9 N-m for the anterior, lateral and posterior regions respectively. Two geometrical properties, radius of gyration and distance from the neutral axis, showed significant variation by region (p < 0.0001) as well as by rib level (p = < 0.01, 0.05). The results of this study can be used to update current models of the human thorax to account for the variation in strength and geometrical properties throughout the rib cage. Accounting for the variation in rib properties by region will improve injury predictive measures and, therefore, the ability to design systems to prevent thoracic injury.

Lateral and posterior dynamic bending of the mid-shaft femur: fracture risk curves for the adult population.

The purpose of this study was to develop injury risk functions for dynamic bending of the human femur in the lateral-to-medial and posterior-to-anterior loading directions. A total of 45 experiments were performed on human cadaver femurs using a dynamic three-point drop test setup. An impactor of 9.8 kg was dropped from 2.2 m for an impact velocity of 5 m/s. Five-axis load cells measured the impactor and support loads, while an in situ strain gage measured the failure strain and subsequent strain rate. All 45 tests resulted in mid-shaft femur fractures with comminuted wedge and oblique fractures as the most common fracture patterns. In the lateral-to-medial bending tests the reaction loads were 4180 +/- 764 N, and the impactor loads were 4780 +/- 792 N. In the posterior-to-anterior bending tests the reaction loads were 3780 +/- 930 N, and the impactor loads were 4310 +/- 1040 N. The difference between the sum of the reaction forces and the applied load is due to inertial effects. The reaction loads were used to estimate the mid-shaft bending moments at failure since there was insufficient data to include the inertial effects in the calculations. The resulting moments are conservative estimates (lower bounds) of the mid-shaft bending moments at failure and are appropriate for use in the assessment of knee restraints and pedestrian impacts with ATD measurements. Regression analysis was used to identify significant parameters, and parametric survival analysis was used to estimate risk functions. Femur cross-sectional area, area moment of inertia (I), maximum distance to the neutral axis (c), I/c, occupant gender, and occupant mass are shown to be significant predictors of fracture tolerance, while no significant difference is shown for loading direction, bone mineral density, leg aspect and age. Risk functions are presented for femur cross-sectional area and I/c as they offer the highest correlation to peak bending moment. The risk function that utilizes the most highly correlated (R2 = 0.82) and significant (p = 0.0001) variable, cross-sectional area, predicts a 50 percent risk of femur fracture of 240 Nm, 395 Nm, and 562 Nm for equivalent cross-sectional area of the 5(th) percentile female, 50(th) percentile male, and 95(th) percentile male respectively.

Defining regional variation in the material properties of human rib cortical bone and its effect on fracture prediction.

This paper presents the results of dynamic material tests and computational modeling that elucidate the effects of regional rib mechanical properties on thoracic fracture patterns. First, a total of 80 experiments were performed using small cortical bone samples from 23 separate locations on the rib cages of four cadavers (2 male, 2 female). Each specimen was subjected to dynamic three-point bending resulting in an average strain rate of 5 +/- 1.5 strain/s. Test coupon modeling was used to verify the test setup. Regional variation was defined by location as anterior, lateral, or posterior as well as by rib level 1 through 12. The specimen stiffness and ultimate stress and strain were analyzed by location and rib level. Second, these material properties were incorporated into a human body computational model. The rib cage was partitioned into anterior, lateral, and posterior segments and the material properties were varied by location using an elastic-plastic material model. A total of 12 simulations with a rigid impactor were performed including 2 separate material assumptions, original and modified rib properties for regional variations, 3 separate impactor velocities, and 2 directions, anterior and lateral. The data from the material tests for all subjects indicate a statistically significant increase in the average stiffness and average ultimate stress for the cortical bone specimens located in the lateral (11.9 GPa modulus, 153.5 MPa ultimate stress) portion of the ribs versus the anterior (7.51 GPa, 116.7 MPa) and posterior (10.7 GPa, 127.7 MPa) rib locations. In addition, the stiffness, ultimate stress, and ultimate strain for all subjects are significantly different by rib level with each variable generally increasing with increasing rib number. The results from the computational modeling for both frontal and lateral impacts illustrate that the location and number of rib fractures are altered by the inclusion of rib material properties that vary by region.

A nonlinear finite element model of the eye with experimental validation for the prediction of globe rupture.

Over 2.4 million eye injuries occur each year in the US, with over 30,000 patients left blind as a result of the trauma. The majority of these injuries occur in automobile crashes, military operations and sporting activities. This paper presents a nonlinear finite element model of the eye and the results of 22 experiments using human eyes to validate for globe rupture injury prediction. The model of the human eye consists of the cornea, sclera, lens, ciliary body, zonules, aqueous humor and vitreous body. Lagrangian membrane elements are used for the cornea and sclera, Lagrangian bricks for the lens, ciliary, and zonules, and Eulerian brick elements comprise the aqueous and vitreous. Nonlinear, isotropic material properties of the sclera and cornea were gathered from uniaxial tensile strip tests performed up to rupture. Dynamic modeling was performed using LS-Dyna. Experimental validation tests consisted of 22 tests using three scenarios: impacts from foam particles, BB's, and baseballs onto fresh eyes used within 24 hours postmortem. The energies of the projectiles were chosen so as to provide both globe rupture and no rupture tests. Displacements of the eye were recorded using high speed color video at 7100 frames per second. The matched simulations predicted rupture of the eye when rupture was seen in the BB and baseball tests, and closely predicted displacements of the eye for the foam tests. Globe rupture has previously been shown to occur at peak stresses of 9.4 MPa using the material properties included in the model. Because of dynamic effects and improvements in boundary conditions resulting from a more realistic modeling of the fluid in the anterior and posterior chambers, the stresses can be much higher than those previously predicted, with the globe remaining intact. The model is empirically verified to predict globe rupture for stresses in the corneoscleral shell exceeding 23 MPa, and local dynamic pressures exceeding 2.1 MPa. The model can be used as a predictive aid to reduce the burden of eye injury, and can serve as a validated model to predict globe rupture.