Frontal crash simulations using parametric human models representing a diverse population

Abstract

Objective: Analyses of crash data have shown that older, obese, and/or female occupants have a higher risk of injury in frontal crashes compared to the rest of the population. The objective of this study was to use parametric finite element (FE) human models to assess the increased injury risks and identify safety concerns for these vulnerable populations.

Methods: We sampled 100 occupants based on age, sex, stature, and body mass index (BMI) to span a wide range of the U.S. adult population. The target anatomical geometry for each of the 100 models was predicted by the statistical geometry models for the rib cage, pelvis, femur, tibia, and external body surface developed previously. A regional landmark-based mesh morphing method was used to morph the Global Human Body Models Consortium (GHBMC) M50-OS model into the target geometries. The morphed human models were then positioned in a validated generic vehicle driver compartment model using a statistical driving posture model. Frontal crash simulations based on U.S. New Car Assessment Program (U.S. NCAP) were conducted. Body region injury risks were calculated based on the risk curves used in the US NCAP, except that scaling was used for the neck, chest, and knee–thigh–hip injury risk curves based on the sizes of the bony structures in the corresponding body regions. Age effects were also considered for predicting chest injury risk.

Results: The simulations demonstrated that driver stature and body shape affect occupant interactions with the restraints and consequently affect occupant kinematics and injury risks in severe frontal crashes. U-shaped relations between occupant stature/weight and head injury risk were observed. Chest injury risk was strongly affected by age and sex, with older female occupants having the highest risk. A strong correlation was also observed between BMI and knee–thigh–hip injury risk, whereas none of the occupant parameters meaningfully affected neck injury risks.

Conclusions: This study is the first to use a large set of diverse FE human models to investigate the combined effects of age, sex, stature, and BMI on injury risks in frontal crashes. The study demonstrated that parametric human models can effectively predict the injury trends for the population and may now be used to optimize restraint systems for people who are not similar in size and shape to the available anthropomorphic test devices (ATDs). New restraints that adapt to occupant age, sex, stature, and body shape may improve crash safety for all occupants.     

Development and Validation of the Total HUman Model for Safety (THUMS) Toward Further Understanding of Occupant Injury Mechanisms in Precrash and During Crash

Objective: Active safety devices such as automatic emergency brake (AEB) and precrash seat belt have the potential to accomplish further reduction in the number of the fatalities due to automotive accidents. However, their effectiveness should be investigated by more accurate estimations of their interaction with human bodies. Computational human body models are suitable for investigation, especially considering muscular tone effects on occupant motions and injury outcomes. However, the conventional modeling approaches such as multibody models and detailed finite element (FE) models have advantages and disadvantages in computational costs and injury predictions considering muscular tone effects. The objective of this study is to develop and validate a human body FE model with whole body muscles, which can be used for the detailed investigation of interaction between human bodies and vehicular structures including some safety devices precrash and during a crash with relatively low computational costs.

Methods: In this study, we developed a human body FE model called THUMS (Total HUman Model for Safety) with a body size of 50th percentile adult male (AM50) and a sitting posture. The model has anatomical structures of bones, ligaments, muscles, brain, and internal organs. The total number of elements is 281,260, which would realize relatively low computational costs. Deformable material models were assigned to all body parts. The muscle–tendon complexes were modeled by truss elements with Hill-type muscle material and seat belt elements with tension-only material. The THUMS was validated against 35 series of cadaver or volunteer test data on frontal, lateral, and rear impacts. Model validations for 15 series of cadaver test data associated with frontal impacts are presented in this article. The THUMS with a vehicle sled model was applied to investigate effects of muscle activations on occupant kinematics and injury outcomes in specific frontal impact situations with AEB.

Results and Conclusions: In the validations using 5 series of cadaver test data, force–time curves predicted by the THUMS were quantitatively evaluated using correlation and analysis (CORA), which showed good or acceptable agreement with cadaver test data in most cases. The investigation of muscular effects showed that muscle activation levels and timing had significant effects on occupant kinematics and injury outcomes. Although further studies on accident injury reconstruction are needed, the THUMS has the potential for predictions of occupant kinematics and injury outcomes considering muscular tone effects with relatively low computational costs.     

Chest injuries of elderly postmortem human surrogates (PMHSs) under seat belt and airbag loading in frontal sled impacts: Comparison to matching THOR tests

Abstract

Objective: The goal of the study was to develop experimental chest loading conditions that would cause up to Abbreviated Injury Scale (AIS) 2 chest injuries in elderly occupants in moderate-speed frontal crashes. The new set of experimental data was also intended to be used in the benchmark of existing thoracic injury criteria in lower-speed collision conditions.

Methods: Six male elderly (age >63) postmortem human subjects (PMHS) were exposed to a 35?km/h (nominal) frontal sled impact. The test fixture consisted of a rigid seat, rigid footrest, and cable seat back. Two restraint conditions (A and B) were compared. Occupants were restrained by a force-limited (2.5?kN [A] and 2?kN [B]) seat belt and a preinflated (16?kPa [A] and 11?kPa [B]; airbag). Condition B also incorporated increased seat friction. Matching sled tests were carried out with the THOR-M dummy. Infra-red telescoping rod for the assessment of chest compression (IRTRACC) readings were used to compute chest injury risk. PMHSs were exposed to a posttest injury assessment. Tests were carried out in 2 stages, using the outcome of the first one combined with a parametric study using the THUMS model to adjust the test conditions in the second. All procedures were approved by the relevant ethics board.

Results: Restraint condition A resulted in an unexpected high number of rib fractures (fx; 10, 14, 15 fx). Under condition B, the adjustment of the relative airbag/occupant position combined with a lower airbag pressure and lower seat belt load limit resulted in a reduced pelvic excursion (85 vs. 110?mm), increased torso pitch and a substantially lower number of rib fractures (1, 0, 4 fx) as intended.

Conclusions: The predicted risk of rib fractures provided by the THOR dummy using the C_max and PC Score injury criteria were lower than the actual injuries observed in the PMHS tests (especially in restraint condition A). However, the THOR dummy was capable of discriminating between the 2 restraint scenarios. Similar results were obtained in the parametric study with the THUMS model.     

THOR dummy chest deflection response in oblique and lateral far-side sled tests

Abstract

Objective: The focus of this study is side impact. Though occupant injury assessment and protection in nearside impacts has received considerable attention and safety standards have been promulgated, field studies show that a majority of far-side occupant injuries are focused on the head and thorax. The 50th percentile male Test Device for Human Occupant Restraint (THOR) has been used in oblique and lateral far-side impact sled tests, and regional body accelerations and forces and moments recorded by load cells have been previously reported. The aim of this study is to evaluate the chestband-based deflection responses from these tests.

Methods: The 3-point belt–restrained 50th percentile male THOR dummy was seated upright in a buck consisting of a rigid flat seat, simulated center console, dashboard, far-side side door structure, and armrest. It was designed to conduct pure lateral and oblique impacts. The center console, dashboard, simulated door structure, and armrest were covered with energy-absorbing materials. A center-mounted airbag was mounted to the right side of the seat. Two 59-gage chestbands were routed on the circumference of the thorax, with the upper and lower chestbands at the level of the third and sixth ribs, respectively, following the rib geometry. Oblique and pure lateral far-side impact tests with and without airbags were conducted at 8.3 m/s. Maximum chest deflections were computed by processing temporal contours using custom software and 3 methods: Procedures paralleling human cadaver studies, using the actual anchor point location and actual alignment of the InfraRed Telescoping Rods for the Assessment of Chest Compression (IR-TRACC) in the dummy on each aspect—that is, right or left,—and using the same anchor location of the internal sensor but determining the location of the peak chest deflection on the contour confined to the aspect of the sensor; these were termed the SD, ID, and TD metrics, respectively.

Results: All deformation contours at the upper and lower thorax levels and associated peak deflections are given for all tests. Briefly, the ID metrics were the lowest in magnitude for both pure lateral and oblique modes, regardless of the presence or absence of an airbag. This was followed by the TD metric, and the SD metric produced the greatest deflections.

Conclusion: The chestbands provide a unique opportunity to compute peak deflections that parallel current IR-TRACC-type deflections and allow computation of peak deflections independent of the initial point of attachment to the rib. The differing locations of the peak deflection vectors along the rib contours for different test conditions suggest that a priori attachment is less effective. Further, varying magnitudes of the differences between ID and TD metrics underscore the difficulty in extrapolating ID outputs under different conditions: Pure lateral versus oblique, airbag presence, and thoracic levels. Deflection measurements should, therefore, not be limited to an instrument that can only track from a fixed point. For improved predictions, these results suggest the need to investigate alternative techniques, such as optical methods to improve chest deflection measurements for far-side occupant injury assessment and mitigation.     

Overview of a combined computational–experimental evaluation for the assessment of panoramic sunroof impact characteristics for ejection mitigation

Abstract

Objective: To meet increasing customer demand, many vehicle manufacturers are now offering a panoramic sunroof option in their vehicle lineup. Currently, there is no regulatory or consumer test aimed at assessing the potential for ejection mitigation of roof glazing, which leaves manufacturers to develop internal performance standards to guide designs. The goal of this study was to characterize the variety of occupant-to-roof impacts involving unbelted occupants in rollover crashes to determine the ranges of possible effective masses and impact velocities. This information can be used to define occupant retention requirements and performance criteria for roof glazing in occupant ejection protection.

Methods: This study combined computational (MADYMO and LS-Dyna) simulations of occupant kinematics in rollover crashes with laboratory rollover crash tests using the dynamic rollover test system (DRoTS) and linked them through controlled anthropomorphic test device (ATD)-to-roof (“drop”) impact tests. The DRoTS and the ATD drop tests were performed to explore impact scenarios and estimate dummy-to-roof impact impulses. Next, 13 sets of vehicle kinematics and deformation data were extracted from a combination of vehicle dynamics and finite element model simulations that reconstructed variations of rollover crash cases from the field data. Then occupant kinematics data were extracted from a full-factorial sensitivity study that used MADYMO simulations to investigate how changes in anthropometry and seating position would affect occupant–roof impacts across all 13 cases. Finite element (FE) simulations of ATD and Global Human Body Models Consortium (GHBMC) human body model (HBM) roof impacts were performed to investigate the most severe cases from the MADYMO simulations to generate a distribution of head-to-roof impact energies.

Results: From the multiparameter design of experiment and experimental study, kinematics and energy output were extracted and analyzed. Based on dummy-to-roof impact force and dummy-to-roof impact velocity, the most severe rollover scenarios were identified. In the DRoTS experiments followed by the drop tests, the range of identified impact velocities was between 2 and 5.8 m/s. However, computational simulations of the rollover crashes showed higher impact velocities and similar effective masses. The largest dummy-to-roof impact velocity was 11 m/s.

Conclusions: This study combined computational and experimental analyses to determine a range of possible unbelted occupant-to-roof impact energies. These results can be used to determine design parameters for an impactor for the assessment of the risk of roof glazing ejection for unbelted occupants in rollover crashes.     

Experimental investigation of the effect of occupant characteristics on contemporary seat belt payout behavior in frontal impacts

Objective: The goal of this study was to investigate the influence of the occupant characteristics on seat belt force vs. payout behavior based on experiment data from different configurations in frontal impacts.

Methods: The data set reviewed consists of 58 frontal sled tests using several anthropomorphic test devices (ATDs) and postmortem human subjects (PMHS), restrained by different belt systems (standard belt, SB; force-limiting belt, FLB) at 2 impact severities (48 and 29 km/h). The seat belt behavior was characterized in terms of the shoulder belt force vs. belt payout behavior. A univariate linear regression was used to assess the factor significance of the occupant body mass or stature on the peak tension force and gross belt payout.

Results: With the SB, the seat belt behavior obtained by the ATDs exhibited similar force slopes regardless of the occupant size and impact severities, whereas those obtained by the PMHS were varied. Under the 48 km/h impact, the peak tension force and gross belt payout obtained by ATDs was highly correlated to the occupant stature (P =.03, P =.02) and body mass (P =.05, P =.04), though no statistical difference with the stature or body mass were noticed for the PMHS (peak force: P =.09, P =.42; gross payout: P =.40, P =.48). With the FLB under the 48 km/h impact, highly linear relationships were noticed between the occupant body mass and the peak tension force (R² = 0.9782) and between the gross payout and stature (R² = 0.9232) regardless of the occupant types.

Conclusions: The analysis indicated that the PMHS characteristics showed a significant influence on the belt response, whereas the belt response obtained with the ATDs was more reproducible. The potential cause included the occupant anthropometry, body mass distribution, and relative motion among body segments specific to the population variance. This study provided a primary data source to understand the biomechanical interaction of the occupant with the restraint system. Further research is necessary to consider these effects in the computational studies and optimized design of the restraint system in a more realistic manner.     

The effect of precrash velocity reduction on occupant response using a human body finite element model

Objective: The objective of this study is to use a validated finite element model of the human body and a certified model of an anthropomorphic test dummy (ATD) to evaluate the effect of simulated precrash braking on driver kinematics, restraint loads, body loads, and computed injury criteria in 4 commonly injured body regions.

Methods: The Global Human Body Models Consortium (GHBMC) 50th percentile male occupant (M50-O) and the Humanetics Hybrid III 50th percentile models were gravity settled in the driver position of a generic interior equipped with an advanced 3-point belt and driver airbag. Fifteen simulations per model (30 total) were conducted, including 4 scenarios at 3 severity levels: median, severe, and the U.S. New Car Assessment Program (U.S.-NCAP) and 3 extra per model with high-intensity braking. The 4 scenarios were no precollision system (no PCS), forward collision warning (FCW), FCW with prebraking assist (FCW+PBA), and FCW and PBA with autonomous precrash braking (FCW + PBA + PB). The baseline ΔV was 17, 34, and 56.4 kph for median, severe, and U.S.-NCAP scenarios, respectively, and were based on crash reconstructions from NASS/CDS. Pulses were then developed based on the assumed precrash systems equipped. Restraint properties and the generic pulse used were based on literature.

Results: In median crash severity cases, little to no risk (<10% risk for Abbreviated injury Scale [AIS] 3+) was found for all injury measures for both models. In the severe set of cases, little to no risk for AIS 3+ injury was also found for all injury measures. In NCAP cases, highest risk was typically found with No PCS and lowest with FCW + PBA + PB. In the higher intensity braking cases (1.0–1.4 g), head injury criterion (HIC), brain injury criterion (BrIC), and chest deflection injury measures increased with increased braking intensity. All other measures for these cases tended to decrease. The ATD also predicted and trended similar to the human body models predictions for both the median, severe, and NCAP cases. Forward excursion for both models decreased across median, severe, and NCAP cases and diverged from each other in cases above 1.0 g of braking intensity.

Conclusions: The addition of precrash systems simulated through reduced precrash speeds caused reductions in some injury criteria, whereas others (chest deflection, HIC, and BrIC) increased due to a modified occupant position. The human model and ATD models trended similarly in nearly all cases with greater risk indicated in the human model. These results suggest the need for integrated safety systems that have restraints that optimize the occupant's position during precrash braking and prior to impact.     

An energy-absorbing sliding seat for reducing neck injury risks in rear impact—analysis for prototype built

Objective: This study investigated overall performance of an energy-absorbing sliding seat concept for whiplash neck injury prevention. The sliding seat allows its seat pan to slide backward for some distance under certain restraint force to absorb crash energy in rear impacts.

Methods: A numerical model that consisted of vehicle interior, seat, seat belt, and BioRID II dummy was built in MADYMO to evaluate whiplash neck injury in rear impact. A parametric study of the effects of sliding seat parameters, including position and cushion stiffness of head restraint, seatback cushion stiffness, recliner characteristics, and especially sliding energy-absorbing (EA) restraint force, on neck injury criteria was conducted in order to compare the effectiveness of the sliding seat concept with that of other existing anti-whiplash mechanisms. Optimal sliding seat design configurations in rear crashes of different severities were obtained. A sliding seat prototype with bending of a steel strip as an EA mechanism was fabricated and tested in a sled test environment to validate the concept. The performance of the sliding seat under frontal and rollover impacts was checked to make sure the sliding mechanism did not result in any negative effects.

Results: The protective effect of the sliding seat with EA restraint force is comparable to that of head restraint–based and recliner stiffness–based anti-whiplash mechanisms. EA restraint force levels of 3 kN in rear impacts of low and medium severities and 6 kN in impacts of high severity were obtained from optimization. In frontal collision and rollover, compared to the nonsliding seat, the sliding seat does not result in any negative effects on occupant protection. The sled test results of the sliding seat prototype have shown the effectiveness of the concept for reducing neck injury risks.

Conclusion: As a countermeasure, the sliding seat with appropriate restraint forces can significantly reduce whiplash neck injury risk in rear impacts of low, medium, and high severities with no negative effects on other crash load cases.     

Small female rib cage fracture in frontal sled tests

Objectives: The 2 objectives of this study are to (1) examine the rib and sternal fractures sustained by small stature elderly females in simulated frontal crashes and (2) determine how the findings are characterized by prior knowledge and field data.

Methods: A test series was conducted to evaluate the response of 5 elderly (average age 76 years) female postmortem human subjects (PMHS), similar in mass and size to a 5th percentile female, in 30 km/h frontal sled tests. The subjects were restrained on a rigid planar seat by bilateral rigid knee bolsters, pelvic blocks, and a custom force-limited 3-point shoulder and lap belt. Posttest subject injury assessment included identifying rib cage fractures by means of a radiologist read of a posttest computed tomography (CT) and an autopsy. The data from a motion capture camera system were processed to provide chest deflection, defined as the movement of the sternum relative to the spine at the level of T8.

?A complementary field data investigation involved querying the NASS-CDS database over the years 1997–2012. The targeted cases involved belted front seat small female passenger vehicle occupants over 40 years old who were injured in 25 to 35 km/h delta-V frontal crashes (11 to 1 o'clock).

Results: Peak upper shoulder belt tension averaged 1,970 N (SD = 140 N) in the sled tests. For all subjects, the peak x-axis deflection was recorded at the sternum with an average of ?44.5 mm or 25% of chest depth. The thoracic injury severity based on the number and distribution of rib fractures yielded 4 subjects coded as Abbreviated Injury Scale (AIS) 3 (serious) and one as AIS 5 (critical). The NASS-CDS field data investigation of small females identified 205 occupants who met the search criteria. Rib fractures were reported for 2.7% of the female occupants.

Conclusions: The small elderly test subjects sustained a higher number of rib cage fractures than expected in what was intended to be a minimally injurious frontal crash test condition. Neither field studies nor prior laboratory frontal sled tests conducted with 50th percentile male PMHS predicted the injury severity observed. Although this was a limited study, the results justify further exploration of the risk of rib cage injury for small elderly female occupants.     

Prediction of probability of fatality due to brain injury in traffic accidents

Abstract

Objective: Fatal brain injuries result from physiological changes in brain tissues, subsequent to primary damage caused by head impact. Although efforts have been made in past studies to estimate the probability of brain injury, none of them involved prediction of such physiological changes. The goal of this study was to evaluate the fatality prediction capability of a novel approach that predicts an increase in intracranial pressure (ICP) due to primary head injury to estimate the fatality rate using clinical data that correlate ICP with fatality rate.

Methods: A total of 12 sets of head acceleration time histories were used to represent no, severe, and fatal brain injury. They were obtained from the literature presenting head kinematics data in noninjurious volunteer sled tests or from accident reconstruction for severe and fatal injury cases. These were first applied to a Global Human Body Models Consortium (GHBMC) head–brain model to predict nodal displacement time histories of the brain, which were then fed into FEBio to predict ICP. A Weibull distribution was applied to the data for the relationship between fatality rate and ICP obtained from a clinical paper to estimate fatality rate from ICP (procedure A). Fatality rate was also estimated by applying the temporal and spatial maximum value of maximum principal strain (MPS_max) obtained from the GHBMC simulation to an injury probability function for MPS_max (procedure B). Estimated fatality rates were compared between the 2 procedures.

Results: Both procedures estimated higher average fatality rate for higher injury severity. The average fatality rate for procedure A without ischemia representation and procedure B was 72.4 and 51.0% for the fatal injury group and 8.2 and 21.7% for the severe injury group, respectively, showing that procedure A provides more distinct classification between fatal and nonfatal brain injury. It was also found that representation of ischemia in procedure A provides results sensitive to injury severity and impact conditions, requiring further validation of the initial estimate for the relationship between brain compression and ischemic cell death.

Conclusions: Prediction of the probability of fatality by means of a combination of simulations of the primary brain deformation and subsequent ICP increase was found to be more distinct compared to the prediction of primary injury alone combined with the injury probability function from a past study in the select 12 head impact cases.     

Computational investigation of the effects of knee airbag design on the interaction with occupant lower extremity in frontal and oblique impacts

Objective: The lower extremity of the occupant represents the most frequently injured body region in motor vehicle crashes. Knee airbags (KABs) have been implemented as a potential countermeasure to reduce lower extremity injuries. Despite the increasing prevalence of KABs in vehicles, the biomechanical interaction of the human lower extremity with the KAB has not been well characterized. This study uses computational models of the human body and KABs to explore how KAB design may influence the impact response of the occupant's lower extremities.

Methods: The analysis was conducted using a 50th percentile male occupant human body model with deployed KABs in a simplified vehicle interior. The 2 common KAB design types, bottom-deploy KAB (BKAB) and rear-deploy KAB (RKAB), were both included. A state-of-the-art airbag modeling technique, the corpuscular particle method, was adopted to represent the deployment dynamics of the unfolding airbags. Validation of the environment model was performed based on previously reported test results. The kinematic responses of the occupant lower extremities were compared under both KAB designs, 2 seating configurations (in-position and out-of-position), and 3 loading conditions (static, frontal, and oblique impacts). A linear statistical model was used to assess factor significance considering the impact responses of the occupant lower extremities.

Results: The presence of a KAB had a significant influence on the lower extremity kinematics compared to no KAB (P <.05) by providing early restraint and distributing contact force on the legs during airbag deployment. For in-position occupants, the KAB generally tended to decrease tibia loadings. The RKAB led to greater lateral motion of the legs compared to the BKAB, resulting in higher lateral displacement at the knee joint and abduction angle change (51.2 ± 21.7 mm and 15° ± 6.0°) over the dynamic loading conditions. Change in the seating position led to a significant difference in occupant kinematic and kinetic parameters (P <.05). For the out-of-position (forward-seated) occupant, the earlier contact between the lower extremity and the deploying KAB resulted in 28.4° ± 5.8° greater abduction, regardless of crash scenarios. Both KAB types reduced the axial force in the femur relative to no KAB. Overall, the out-of-position occupant sustained a raised axial force and bending moment of the tibia by 0.8 ± 0.2 kN and 21.1 ± 8.7 Nm regardless of restraint use.

Conclusions: The current study provided a preliminary computational examination on KAB designs based on a limited set of configurations in an idealized vehicle interior. Results suggested that the BKAB tended to provide more coverage and less leg abduction compared to the RKAB in oblique impact and/or the selected out-of-position scenario. An out-of-position occupant was associated with larger abduction and lower extremity loads over all occupant configurations. Further investigations are recommended to obtain a full understanding of the KAB performance in a more realistic vehicle environment.     

Assessment of Bilateral Thoracic Loading on the Near-Side Occupant Due to Occupant-to-Occupant Interaction in Vehicle Crash Tests

Objective: This study aims, by means of the WorldSID 50th percentile male, to evaluate thoracic loading and injury risk to the near-side occupant due to occupant-to-occupant interaction in combination with loading from an intruding structure.

Method: Nine vehicle crash tests were performed with a 50th percentile WorldSID male dummy in the near-side (adjacent to the intruding structure) seat and a THOR or ES2 dummy in the far-side (opposite the intruding structure) seat. The near-side seated WorldSID was equipped with 6 + 6 IR-Traccs (LH and RH) in the thorax/abdomen enabling measurement of bilateral deflection. To differentiate deflection caused by the intrusion, and the deflection caused by the neighboring occupant, time history curves were analyzed. The crash tests were performed with different modern vehicles, equipped with thorax side airbags and inflatable curtains, ranging from a compact car to a large sedan, and in different loading conditions such as car-to-car, barrier, and pole tests. Lateral delta V based on vehicle tunnel acceleration and maximum residual intrusion at occupant position were used as a measurement of crash severity to compare injury measurements.

Result: In the 9 vehicle crash tests, thoracic loading, induced by the intruding structure as well as from the far-side occupant, varied due to the size and structural performance of the car as well as the severity of the crash. Peak deflection on the thoracic outboard side occurred during the first 50 ms of the event. Between 70 to 150 ms loading induced by the neighboring occupant occurred and resulted in an inboard-side peak deflection and viscous criterion. In the tests where the target vehicle lateral delta V was below 30 km/h and intrusion less than 200 mm, deflections were low on both the outboard (20–40 mm) and inboard side (10–15 mm). At higher crash severities, delta V 35 km/h and above as well as intrusions larger than 350 mm, the inboard deflections (caused by interaction to the far-side occupant) were of the same magnitude or even higher (30–70 mm) than the outboard deflections (30–50 mm).

Conclusion: A WorldSID 50th percentile male equipped with bilateral IR-Traccs can detect loading to the thorax from a neighboring occupant making injury risk assessment feasible for this type of loading. At crash severities resulting in a delta V above 35 km/h and intrusions larger than 350 mm, both the inboard deflection and VC resulted in high risks of Abbreviated Injury Scale (AIS) 3+ injury, especially for a senior occupant.     

Driver Injury Risk Variability in Finite Element Reconstructions of Crash Injury Research and Engineering Network (CIREN) Frontal Motor Vehicle Crashes

Objective: A 3-phase real-world motor vehicle crash (MVC) reconstruction method was developed to analyze injury variability as a function of precrash occupant position for 2 full-frontal Crash Injury Research and Engineering Network (CIREN) cases.

Method: Phase I: A finite element (FE) simplified vehicle model (SVM) was developed and tuned to mimic the frontal crash characteristics of the CIREN case vehicle (Camry or Cobalt) using frontal New Car Assessment Program (NCAP) crash test data. Phase II: The Toyota HUman Model for Safety (THUMS) v4.01 was positioned in 120 precrash configurations per case within the SVM. Five occupant positioning variables were varied using a Latin hypercube design of experiments: seat track position, seat back angle, D-ring height, steering column angle, and steering column telescoping position. An additional baseline simulation was performed that aimed to match the precrash occupant position documented in CIREN for each case. Phase III: FE simulations were then performed using kinematic boundary conditions from each vehicle's event data recorder (EDR). HIC15, combined thoracic index (CTI), femur forces, and strain-based injury metrics in the lung and lumbar vertebrae were evaluated to predict injury.

Results: Tuning the SVM to specific vehicle models resulted in close matches between simulated and test injury metric data, allowing the tuned SVM to be used in each case reconstruction with EDR-derived boundary conditions. Simulations with the most rearward seats and reclined seat backs had the greatest HIC15, head injury risk, CTI, and chest injury risk. Calculated injury risks for the head, chest, and femur closely correlated to the CIREN occupant injury patterns. CTI in the Camry case yielded a 54% probability of Abbreviated Injury Scale (AIS) 2+ chest injury in the baseline case simulation and ranged from 34 to 88% (mean = 61%) risk in the least and most dangerous occupant positions. The greater than 50% probability was consistent with the case occupant's AIS 2 hemomediastinum. Stress-based metrics were used to predict injury to the lower leg of the Camry case occupant. The regional-level injury metrics evaluated for the Cobalt case occupant indicated a low risk of injury; however, strain-based injury metrics better predicted pulmonary contusion. Approximately 49% of the Cobalt occupant's left lung was contused, though the baseline simulation predicted 40.5% of the lung to be injured.

Conclusions: A method to compute injury metrics and risks as functions of precrash occupant position was developed and applied to 2 CIREN MVC FE reconstructions. The reconstruction process allows for quantification of the sensitivity and uncertainty of the injury risk predictions based on occupant position to further understand important factors that lead to more severe MVC injuries.     

Influence of Seat Properties on Occupant Dynamics in Severe Rear Crashes

Seat performance in retaining an occupant, transferring energy, and controlling neck responses is often questioned after severe rear crashes when fatal or disabling injury occur. It is argued that a stiffer seat would have improved occupant kinematics. However, there are many factors in occupant interactions with the seat. This study evaluates four different seat types in 26 and 32 mph (42 and 51 km/h), rear crash delta Vs. Two seats were yielding with k = 20 kN/m occupant load per displacement. One represented a 1970s yielding seat with j = 3.4°/kN frame rotation per occupant load, and 3 kN maximum load (660 Nm moment), and the other a high retention seat phased into production since 1997 with j = 1.4°/kN, and 10 kN maximum load (2200 Nm). Two seats were stiff with k = 40 kN/m. One represented a 1990s foreign benchmark with j = 1.8°/kN and a 7.7 kN maximum load (1700 Nm), and the other an all belts to seat (ABTS) with j = 1.0°/kN and 20 kN maximum load (4400 Nm). The crash was a constant acceleration of 11.8 g, or 14.5 g for 100 ms. Occupant interactions with the seat were modeled using a torso mass, flexible neck and head mass. By analysis of the equations of motion, the initial change in seatback angle (Δθ) is proportional to jk(y ? x), the product jk and the differential motion between the vehicle (seat cushion) and occupant. The transition from 1970s–80s yielding seats to stronger seats of the 1990s involved an increase in k stiffness; however, the jk property did not change as frame structures became stronger. The yielding seats of the 1970s had jk = 68°/m, while the stiff foreign benchmark seat had jk = 72°/m. The foreign benchmark rotated about the same as the 1970s seat up to 50 ms in the severe rear crashes. While it was substantially stronger, it produced higher loads on the occupant, and the higher loads increased seatback rotations and neck responses. The ABTS seat had the lowest rotations but also caused high neck responses because of the greater loads on the torso. Neck displacement (d) is initially proportional to (k/m_T) ∫∫ y, seat stiffness times the second integral of vehicle displacement divided by torso mass. As seat stiffness increases, head-torso acceleration, velocity, and neck displacement increase. This study shows that the jk seat property determines the initial seatback rotation in rear crashes. If a stronger seat has a higher stiffness, it rotates at higher loads on the occupant, reducing the overall benefit of the stronger frame, while increasing neck responses related to whiplash or neck extension prior to subsequent impacts. The aim of seat designs should be to reduce jk, provide pocketing of the pelvis, and give head-neck support for the best protection in severe rear crashes. For low-speed crashes, a low k is important to reduce early neck responses related to whiplash.     

Injuries in Full-Scale Vehicle Side Impact Moving Deformable Barrier and Pole Tests Using Postmortem Human Subjects

Objective: To conduct near-side moving deformable barrier (MDB) and pole tests with postmortem human subjects (PMHS) in full-scale modern vehicles, document and score injuries, and examine the potential for angled chest loading in these tests to serve as a data set for dummy biofidelity evaluations and computational modeling.

Methods: Two PMHS (outboard left front and rear seat occupants) for MDB and one PMHS (outboard left front seat occupant) for pole tests were used. Both tests used sedan-type vehicles from same manufacturer with side airbags. Pretest x-ray and computed tomography (CT) images were obtained. Three-point belt-restrained surrogates were positioned in respective outboard seats. Accelerometers were secured to T1, T6, and T12 spines; sternum and pelvis; seat tracks; floor; center of gravity; and MDB. Load cells were used on the pole. Biomechanical data were gathered at 20 kHz. Outboard and inboard high-speed cameras were used for kinematics. X-rays and CT images were taken and autopsy was done following the test. The Abbreviated Injury Scale (AIS) 2005 scoring scheme was used to score injuries.

Results: MDB test: male (front seat) and female (rear seat) PMHS occupant demographics: 52 and 57 years, 177 and 166 cm stature, 78 and 65 kg total body mass. Demographics of the PMHS occupant in the pole test: male, 26 years, 179 cm stature, and 84 kg total body mass. Front seat PMHS in MDB test: 6 near-side rib fractures (AIS = 3): 160–265 mm vertically from suprasternal notch and 40–80 mm circumferentially from center of sternum. Left rear seat PMHS responded with multiple bilateral rib fractures: 9 on the near side and 5 on the contralateral side (AIS = 3). One rib fractured twice. On the near and contralateral sides, fractures were 30–210 and 20–105 mm vertically from the suprasternal notch and 90–200 and 55–135 mm circumferentially from the center of sternum. A fracture of the left intertrochanteric crest occurred (AIS = 3). Pole test PMHS had one near-side third rib fracture. Thoracic accelerations of the 2 occupants were different in the MDB test. Though both occupants sustained positive and negative x-accelerations to the sternum, peak magnitudes and relative changes were greater for the rear than the front seat occupant. Magnitudes of the thoracic and sternum accelerations were lower in the pole test.

Conclusions: This is the first study to use PMHS occupants in MDB and pole tests in the same recent model year vehicles with side airbag and head curtain restraints. Injuries to the unilateral thorax for the front seat PMHS in contrast to the bilateral thorax and hip for the rear seat occupant in the MDB test indicate the effects of impact on the seating location and restraint system. Posterolateral locations of fractures to the front seat PMHS are attributed to constrained kinematics of occupant interaction with torso side airbag restraint system. Angled loading to the rear seat occupant from coupled sagittal and coronal accelerations of the sternum representing anterior thorax loading contributed to bilateral fractures. Inward bending initiated by the distal femur complex resulting in adduction of ipsilateral lower extremity resulted in intertrochanteric fracture to the rear seat occupant. These results serve as a data set for evaluating the biofidelity of the WorldSID and federalized side impact dummies and assist in validating human body computational models, which are increasingly used in crashworthiness studies.     

Is sledding safe? Beyond acute,traumatic injuries: Estimating the spine loads sustained by children while sledding

ObjectiveCompare acceleration level between two sleds used for downhill winter sliding. The hypotheses include (a) acceleration levels will differ significantly between sled types, (b) acceleration levels will be unsafe, (c) acceleration levels will not differ significantly between genders.DesignCross sectional study, within subject design. A cross correlation was used to verify similar terrain and path of travel between the two sleds. The trials with the highest correlation, between sleds, were then used in a t-test to analyze the differences in the mean, maximum, minimum, peak accelerations, average and maximum peak levels between the two sleds and genders. A one-way ANOVA evaluated the relationship among gender and all variables.SettingIn a randomized order, participants slid down the hill three times with one sled, followed by three times with the other. Acceleration was captured with a PBC piezoelectric accelerometer connected to a Larson Davis Human Vibration Meter 100.ParticipantsFive males and seven females, age 7–14 years, were recruited for this study. Participants were a sample of convenience.ResultsThe measured accelerations exceeded safe levels with values ranging from 1.07 m/s² to 1330 m/s² (mean = 241.97 m/s²). Two-tailed t-test showed no significant difference in acceleration between sleds. A one-way ANOVA determined that there is no significant difference between gender or sled for any measured variables.ConclusionAcceleration magnitudes did not differ significantly between sled types or between genders. Peak acceleration levels were determined to be unsafe, in particular with regard to the forces transmitted to the spine.     

Influence of seat properties on occupant dynamics in severe rear crashes

Seat performance in retaining an occupant, transferring energy, and controlling neck responses is often questioned after severe rear crashes when fatal or disabling injury occur. It is argued that a stiffer seat would have improved occupant kinematics. However, there are many factors in occupant interactions with the seat. This study evaluates four different seat types in 26 and 32 mph (42 and 51 km/h), rear crash delta Vs. Two seats were yielding with k = 20 kN/m occupant load per displacement. One represented a 1970s yielding seat with j = 3.4 degrees /kN frame rotation per occupant load, and 3 kN maximum load (660 Nm moment), and the other a high retention seat phased into production since 1997 with j = 1.4 degrees /kN, and 10 kN maximum load (2200 Nm). Two seats were stiff with k = 40 kN/m. One represented a 1990s foreign benchmark with j = 1.8 degrees /kN and a 7.7 kN maximum load (1700 Nm), and the other an all belts to seat (ABTS) with j = 1.0 degrees /kN and 20 kN maximum load (4400 Nm). The crash was a constant acceleration of 11.8 g, or 14.5 g for 100 ms. Occupant interactions with the seat were modeled using a torso mass, flexible neck and head mass. By analysis of the equations of motion, the initial change in seatback angle (Deltatheta) is proportional to jk(y - x), the product jk and the differential motion between the vehicle (seat cushion) and occupant. The transition from 1970s-80s yielding seats to stronger seats of the 1990s involved an increase in k stiffness; however, the jk property did not change as frame structures became stronger. The yielding seats of the 1970s had jk = 68 degrees /m, while the stiff foreign benchmark seat had jk = 72 degrees /m. The foreign benchmark rotated about the same as the 1970s seat up to 50 ms in the severe rear crashes. While it was substantially stronger, it produced higher loads on the occupant, and the higher loads increased seatback rotations and neck responses. The ABTS seat had the lowest rotations but also caused high neck responses because of the greater loads on the torso. Neck displacement (d) is initially proportional to (k/m(T)) integral integral y, seat stiffness times the second integral of vehicle displacement divided by torso mass. As seat stiffness increases, head-torso acceleration, velocity, and neck displacement increase. This study shows that the jk seat property determines the initial seatback rotation in rear crashes. If a stronger seat has a higher stiffness, it rotates at higher loads on the occupant, reducing the overall benefit of the stronger frame, while increasing neck responses related to whiplash or neck extension prior to subsequent impacts. The aim of seat designs should be to reduce jk, provide pocketing of the pelvis, and give head-neck support for the best protection in severe rear crashes. For low-speed crashes, a low k is important to reduce early neck responses related to whiplash.     

大客车正面碰撞乘客约束系统仿真分析

为探究约束系统在全承载客车正面碰撞事故中对乘客损伤的影响,利用有限元分析软件LS_DYNA建立某大客车正面碰撞仿真模型,并开展整车50 km/h正面100％重叠碰撞固定刚性壁障试验;从车身变形、加速度曲线和乘员损伤等3方面验证仿真模型;基于已验证的仿真模型,开展不同座椅间距、车厢位置及安全带类型的乘员运动响应和损伤等综...     

Relationship between frontal car-to-car test result and vehicle crash compatibility evaluation in mobile progressive deformable barrier test

Abstract

Objective: In 2020, the world’s first crash compatibility rating test will be introduced in the European mobile progressive deformable barrier (MPDB) test. In this research, the quantitative change in partner protection performance of large vehicles in car-to-car (C2C) impacts was studied if these large vehicles were designed in future based on MPDB tests addressing crash compatibility ratings.

Methods: Representative vehicles of the European fleet were selected and a Computer Aided Engineering (CAE) parameter study was conducted. In particular, by changing an indicator of structural interaction performance (SD; i.e., the degree of uniformity of barrier deformation)/mass/stiffness of large vehicles systematically in a step-by-step approach, the compatibility evaluation results of large vehicles in MPDB and the occupant injury score of small vehicles in C2C impacts were compared. The CAE result was evaluated compared to that of C2C physical impact tests.

Results: The CAE parameter study showed that in the C2C impact condition, the effects on occupant injury in a small vehicle due to changes in the large vehicle were as follows: (1) SD change: The effect was minor except for small overlap condition. (2) Mass and stiffness change: The effect was relatively major.

On the other hand, compatibility evaluation in the MPDB showed a tendency to overestimate the effect of SD change in comparison with the above-mentioned C2C impact condition.

In addition, physical impact tests showed that, based on SD evaluation, the large vehicle with a relatively inferior compatibility rating compared to those with superior compatibility ratings showed a contradicting trend of better compatibility performance in the C2C test.

Conclusions: The currently proposed compatibility evaluation method of the MPDB test showed some tendency to overestimate the effect of SD change and resulted in quantitatively inconsistent outcomes regarding occupant injury in the partner car in C2C impact conditions.