Validation of a booted finite element model of the WIAMan ATD lower limb in component and whole-body vertical loading impacts with an assessment of the boot influence model on response

Objective: A novel anthropomorphic test device (ATD) representative of the 50th percentile male soldier is being developed to predict injuries to a vehicle occupant during an underbody blast (UBB). The main objective of this study was to develop and validate a finite element (FE) model of the ATD lower limb outfitted with a military combat boot and to insert the validated lower limb into a model of the full ATD and simulate vertical loading experiments.

Methods: A Belleville desert combat boot model was assigned contacts and material properties based on previous experiments. The boot model was fit to a previously developed model of the barefoot ATD. Validation was performed through 6 matched pair component tests conducted on the Vertically Accelerated Loads Transfer System (VALTS). The load transfer capabilities of the FE model were assessed along with the force-mitigating properties of the boot. The booted lower limb subassembly was then incorporated into a whole-body model of the ATD. Two whole-body VALTS experiments were simulated to evaluate lower limb performance in the whole body.

Results: The lower limb model accurately predicted axial loads measured at heel, tibia, and knee load cells during matched pair component tests. Forces in booted simulations were compared to unbooted simulations and an amount of mitigation similar to that of experiments was observed. In a whole-body loading environment, the model kinematics match those recorded in experiments. The shape and magnitude of experimental force–time curves were accurately predicted by the model. Correlation between the experiments and simulations was backed up by high objective rating scores for all experiments.

Conclusion: The booted lower limb model is accurate in its ability to articulate and transfer loads similar to the physical dummy in simulated underbody loading experiments. The performance of the model leads to the recommendation to use it appropriately as an alternative to costly ATD experiments.     

现代生物力学在煤矿安全科技中的应用研究

笔者根据国内外现代生物力学的发展情况 ,结合我国煤矿安全的现状 ,探讨现代生物力学在煤矿安全方面的应用 ,并提出了现代生物力学在安全科技两个方面的应用 :①研究人体结构在冲击事故的破坏极限 ,对安全帽的防护性能做简略的生物力学分析 ;②骨力———电性质的研究对骨折愈合与康复的意义。笔者提出随着现代科学技术的迅猛发展 ,煤矿安全方面的研究能否从现代生物力学在人体安全保护、骨折治疗与康复、职业病防治机理等方面的应用去研究的新思想和新方法 ,使煤矿安全的研究达到一个新的水平。     

Development of a Computationally Efficient Full Human Body Finite Element Model

Introduction: A simplified and computationally efficient human body finite element model is presented. The model complements the Global Human Body Models Consortium (GHBMC) detailed 50th percentile occupant (M50-O) by providing kinematic and kinetic data with a significantly reduced run time using the same body habitus.

Methods: The simplified occupant model (M50-OS) was developed using the same source geometry as the M50-O. Though some meshed components were preserved, the total element count was reduced by remeshing, homogenizing, or in some cases omitting structures that are explicitly contained in the M50-O. Bones are included as rigid bodies, with the exception of the ribs, which are deformable but were remeshed to a coarser element density than the M50-O. Material models for all deformable components were drawn from the biomechanics literature. Kinematic joints were implemented at major articulations (shoulder, elbow, wrist, hip, knee, and ankle) with moment vs. angle relationships from the literature included for the knee and ankle. The brain of the detailed model was inserted within the skull of the simplified model, and kinematics and strain patterns are compared.

Results: The M50-OS model has 11 contacts and 354,000 elements; in contrast, the M50-O model has 447 contacts and 2.2 million elements. The model can be repositioned without requiring simulation. Thirteen validation and robustness simulations were completed. This included denuded rib compression at 7 discrete sites, 5 rigid body impacts, and one sled simulation. Denuded tests showed a good match to the experimental data of force vs. deflection slopes. The frontal rigid chest impact simulation produced a peak force and deflection within the corridor of 4.63 kN and 31.2%, respectively. Similar results vs. experimental data (peak forces of 5.19 and 8.71 kN) were found for an abdominal bar impact and lateral sled test, respectively. A lateral plate impact at 12 m/s exhibited a peak of roughly 20 kN (due to stiff foam used around the shoulder) but a more biofidelic response immediately afterward, plateauing at 9 kN at 12 ms. Results from a frontal sled simulation showed that reaction forces and kinematic trends matched experimental results well. The robustness test demonstrated that peak femur loads were nearly identical to the M50-O model. Use of the detailed model brain within the simplified model demonstrated a paradigm for using the M50-OS to leverage aspects of the M50-O. Strain patterns for the 2 models showed consistent patterns but greater strains in the detailed model, with deviations thought to be the result of slightly different kinematics between models. The M50-OS with the deformable skull and brain exhibited a run time 4.75 faster than the M50-O on the same hardware.

Conclusions: The simplified GHBMC model is intended to complement rather than replace the detailed M50-O model. It exhibited, on average, a 35-fold reduction in run time for a set of rigid impacts. The model can be used in a modular fashion with the M50-O and more broadly can be used as a platform for parametric studies or studies focused on specific body regions.     

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.     

The Contribution of Pre-impact Posture on Restrained Occupant Finite Element Model Response in Frontal Impact

Objective: The objective of this study was to discuss the influence of the pre-impact posture to the response of a finite element human body model (HBM) in frontal impacts.

Methods: This study uses previously published cadaveric tests (PMHS), which measured six realistic pre-impact postures. Seven postured models were created from the THUMS occupant model (v4.0): one matching the standard UMTRI driving posture as it was the target posture in the experiments, and six matching the measured pre-impact postures. The same measurements as those obtained during the cadaveric tests were calculated from the simulations, and biofidelity metrics based on signals correlation (CORA) were established to compare the response of the seven models to the experiments.

Results: The HBM responses showed good agreement with the PMHS responses for the reaction forces (CORA = 0.80 ± 0.05) and the kinematics of the lower part of the torso but only fair correlation was found with the head, the upper spine, rib strains (CORA= 0.50 ± 0.05) and chest deflections (CORA = 0.67 ± 0.08). All models sustained rib fractures, sternal fracture and clavicle fracture. The average number of rib fractures for all the models was 5.3 ± 1.0, lower than in the experiments (10.8 ± 9.0).

Variation in pre-impact posture greatly altered the time histories of the reaction forces, deflections and the rib strains, mainly in terms of time delay, but no definite improvement in HBM response or injury prediction was observed. By modifying only the posture of the HBM, the variability in the impact response was found to be equivalent to that observed in the experiments. The postured HBM sustained from 4 to 8 rib fractures, confirming that the pre-impact posture influenced the injury outcome predicted by the simulation.

Conclusions: This study tries to answer an important question: what is the effect of occupant posture on kinematics and kinetics. Significant differences in kinematics observed between HBM and PMHS suggesting more coupling between the pelvis and the spine for the models which makes the model response very sensitive to any variation in the spine posture. Consequently, the findings observed for the HBM cannot be extended to PMHS. Besides, pre-impact posture should be carefully quantified during experiments and the evaluation of HBM should take into account the variation in the predicted impact response due to the variation in the model posture.     

Biomechanics of Brain and Spinal-Cord Injury: Analysis of Neuropathologic and Neurophysiology Experiments

This study analyzes 46 brain and 48 spinal-cord impact experiments. The velocity of brain impact was 2.0-10.0 m/s and displacement, 0.75-5.0 mm (5.3-33% compression) using a controlled pneumatic impact. The velocity of spinal-cord impact was 1.5-6.0 m/s and displacement, 1.25-3.25 mm (25-65% compression). Brain injury varied from cortical contusion, diffuse axonal injury (DAI), to fatalities, and spinal-cord injury from temporary to complete loss of somatosensory-evoked potentials. Logist functions were determined for each injury severity and various biomechanical parameters, VC, C, V, and combinations. Brain and spinal-cord injury is most strongly correlated to VC, the viscous response. The goodness-of-fit was x² = 22.1, R-0.84 and p< 0.0000 for fatal brain injury, x² = 27.5, R = 0.96 and p< 0.0000 for cortical contusion, and x² = 17.7, R = 0.49 and p < 0.0001 for partial recovery of spinal-cord conduction. Neural tissue is viscoelastic, with a rate-dependent tolerance related to energy absorption. VC is a measure of energy absorption by impact deformation and is predictive of neural contusion, DAI, long-duration coma, spinal-cord dysfunction, and death. Tolerances for various severities of neural injury are presented. At the tissue level, VC is the product of strain and strain-rate, ε dε/dt. The research shows that strain is not a sufficient parameter of neural injury risk, and that the product of strain and strain-rate is a key biomechanical parameter for brain and spinal-cord injury.     

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.     

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.     

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.     

Load Acceleration and Footstep Strategies in Asymmetrical Lifting and Lowering

Accelerated execution effects for lifting and lowering a 12-kg box using two footstep strategies associated with experienced workers were studied. Eight healthy male participants performed a normal and an accelerated execution of a lifting task and a lowering task, using a minimal feet displacement strategy (oblique-step) and a strategy with a step (crossed-step). It was hypothesized that the accelerated executions, as compared to the normal executions, would have a different effect on L5/S1 resultant moment, body posture, and other kinematic variables. A tridimensional dynamic rigid body model was used to compute L5/S1 resultant moments. Results showed that the accelerated condition did not reduce body asymmetry of posture, but it reduced the length of the path of the global center of gravity and the duration of the supporting phase of the box, and it did not significantly affect L5/S1 maximal resultant moments for lifting but increased them for lowering. These results indicate that the net work production for accelerated strategies might be smaller, which may represent an economy of energy. Furthermore, the results showed that the use of an accelerated strategy for lowering should be avoided.