Lower Extremity Injury Risk Curve Development for a Human Body Model in the Underbody Blast Environment

Computational human body models (HBMs) provide the ability to explore numerous candidate injury metrics ranging from local strain based criteria to global combined criteria such as the Tibia Index. Despite these efforts, there have been relatively few studies that focus on determining predicted injury risk from HBMs based on observed postmortem human subjects (PMHS) injury data. Additionally, HBMs provide an opportunity to construct risk curves using measures that are difficult or impossible to obtain experimentally. The Global Human Body Models Consortium (GHBMC) M50-O v 6.0 lower extremity was simulated in 181 different loading conditions based on previous PMHS tests in the underbody blast (UBB) environment and 43 different biomechanical metrics were output. The Brier Metric Score were used to determine the most appropriate metric for injury risk curve development. Using survival analysis, three different injury risk curves (IRC) were developed: "any injury," "calcaneus injury," and "tibia injury." For each injury risk curve, the top three metrics selected using the Brier Metric Score were tested for significant covariates including boot use and posture. The best performing metric for the "any injury," "calcaneus injury" and "tibia injury" cases were calcaneus strain, calcaneus force, and lower tibia force, respectively. For the six different injury risk curves where covariates were considered, the presence of the boot was found to be a significant covariate reducing injury risk in five out of six cases. Posture was significant for only one curve. The injury risk curves developed from this study can serve as a baseline for model injury prediction, personal protective equipment (PPE) evaluation, and can aid in larger scale testing and experimental protocols.

Modular incorporation of deformable spine and 3D neck musculature into a simplified human body finite element model

Spinal injuries are a concern for automotive applications, requiring large parametric studies to understand spinal injury mechanisms under complex loading conditions. Finite element computational human body models (e.g. Global Human Body Models Consortium (GHBMC) models) can be used to identify spinal injury mechanisms. However, the existing GHBMC detailed models (with high computational time) or GHBMC simplified models (lacking vertebral fracture prediction capabilities) are not ideal for studying spinal injury mechanisms in large parametric studies. To overcome these limitations, a modular 50th percentile male simplified occupant model combining advantages of both the simplified and detailed models, M50-OS + DeformSpine, was developed by incorporating the deformable spine and 3D neck musculature from the detailed GHBMC model M50-O (v6.0) into the simplified GHBMC model M50-OS (v2.3). This new modular model was validated against post-mortem human subject test data in four rigid hub impactor tests and two frontal impact sled tests. The M50-OS + DeformSpine model showed good agreement with experimental test data with an average CORrelation and Analysis (CORA) score of 0.82 for the hub impact tests and 0.75 for the sled impact tests. CORA scores were statistically similar overall between the M50-OS + DeformSpine (0.79 ± 0.11), M50-OS (0.79 ± 0.11), and M50-O (0.82 ± 0.11) models (p > 0.05). This new model is computationally 6 times faster than the detailed M50-O model, with added spinal injury prediction capabilities over the simplified M50-OS model.

A simulation-based study for optimizing proportional-integral-derivative controller gains for different control strategies of an active upper extremity model using experimental data

This study investigates the effect of PID controller gains, reaction time, and initial muscle activation values on active human model behavior while comparing three different control strategies. The controller gains and reaction delays were optimized using published experimental data focused on the upper extremity. The data describes the reaction of five male subjects in four tests based on two muscle states (relaxed and tensed) and two states of awareness (open and closed eye). The study used a finite element model of the left arm isolated from the Global Human Body Models Consortium (GHBMC) average male simplified occupant model for simulating biomechanical simulations. Major skeletal muscles of the arm were modeled as 1D beam elements and assigned a Hill-type muscle material. Angular position control, muscle length control, and a combination of both were used as a control strategy. The optimization process was limited to 4 variables; three Proportional-Integral-Derivative (PID) controller gains and one reaction delay time. The study assumed the relaxed and tensed condition require distinct sets of controller gains and initial activation and that the closed-eye simulations can be achieved by increasing the reaction delay parameter. A post-hoc linear combination of angle and muscle length control was used to arrive at the final combined control strategy. The premise was supported by variation in the controller gains depending on muscle state and an increase in reaction delay based on awareness. The CORA scores for open-eye relaxed, closed-eye relaxed, open-eye tensed, and closed-eye tensed was 0.95, 0.90, 0.95, and 0.77, respectively using the combined control strategy.

Contrast enhanced computed tomography of small ruminants: Caprine and ovine

The use of small ruminants, mainly sheep and goats, is increasing in biomedical research. Small ruminants are a desirable animal model due to their human-like anatomy and physiology. However, the large variability between studies and lack of baseline data on these animals creates a barrier to further research. This knowledge gap includes a lack of computed tomography (CT) scans for healthy subjects. Full body, contrast enhanced CT scans of caprine and ovine subjects were acquired for subsequent modeling studies. Scans were acquired from an ovine specimen (male, Khatadin, 30-35 kg) and caprine specimen (female, Nubian 30-35 kg). Scans were acquired with and without contrast. Contrast enhanced scans utilized 1.7 mL/kg of contrast administered at 2 mL/s and scans were acquired 20 seconds, 80 seconds, and 5 minutes post-contrast. Scans were taken at 100 kV and 400 mA. Each scan was reconstructed using a bone window and a soft tissue window. Sixteen full body image data sets are presented (2 specimens by 4 contrast levels by 2 reconstruction windows) and are available for download through the form located at: https://redcap.link/COScanData. Scans showed that the post-contrast timing and scan reconstruction method affected structural visualization. The data are intended for further biomedical research on ruminants related to computational model development, device prototyping, comparative diagnostics, intervention planning, and other forms of translational research.

Preliminary validation of the GHBMC average male occupant models and 70YO aged model in far-side impact

The objective of the current study was to perform a preliminary validation of the Global Human Body Models Consortium (GHBMC) average male occupant models, simplified (M50-OS) and detailed (M50-O) and the 70YO aged model in Far-side impacts and compare the head kinematics against the PMHS responses published by Petit et al. (2019). The buck used to simulate the far-side impacts comprised a seat, headrest, center console plate, leg support plate, and footrest plate with rigid material properties. The three occupant models were gravity settled onto the rigid seat and belted with a 3-point seatbelt. Positioning details of the PMHS were followed in the model setup process. A deceleration pulse with ΔV of 8 m/s was applied. The far-side crash simulations were performed with and without the addition of a plexiglass cover around the setup similar to the experimental setup. The head kinematics were extracted from the models for comparison against the PMHS data. Peak head displacements in Y and Z axes from the three models were compared to the PMHS data in addition to the head rotation along X axes. The peak head displacement values for the M50-OS, M50-O, and M50-O 70YO aged models are 594.10 mm, 568.44 mm, and 567.90 mm along Y and 325.21 mm, 402.66 mm, and 375.92 mm respectively along Z when the plexiglass cover is included in the test. The peak head rotation values for the M50-OS, M50-O, and M50-O 70YO aged models are 95.64°, 122.15°, and 129.08° respectively when the plexiglass cover is included in the test. The three occupant models capture the general trend of the PMHS data. The detailed occupant models have higher head rotation compared to the simplified model because of the deformable structure of the spine and intervertebral discs modeled. These three occupant models can be used for further parametric studies in this condition to study the influence of restraint parameters.

Assessment of finite element human body and ATD models in estimating injury risk in far-side impacts using field-based injury risk

The objective of this study was to assess the ability of finite element human body models (FEHBMs) and Anthropometric Test Device (ATD) models to estimate occupant injury risk by comparing it with field-based injury risk in far-side impacts. The study used the Global Human Body Models Consortium midsize male (M50-OS+B) and small female (F05-OS+B) simplified occupant models with a modular detailed brain, and the ES-2Re and SID-IIs ATD models in the simulated far-side crashes. A design of experiments (DOE) with a total of 252 simulations was conducted by varying lateral ΔV (10-50kph; 5kph increments), the principal direction of force (PDOF 50°, 60°, 65°, 70°, 75°, 80°, 90°), and occupant models. Models were gravity-settled and belted into a simplified vehicle model (SVM) modified for far-side impact simulations. Acceleration pulses and vehicle intrusion profiles used for the DOE were generated by impacting a 2012 Camry vehicle model with a mobile deformable barrier model across the 7 PDOFs and 9 lateral ΔV's in the DOE for a total of 63 additional simulations. Injury risks were estimated for the head, chest, lower extremity, pelvis (AIS 2+; AIS 3+), and abdomen (AIS 3+) using logistic regression models. Combined AIS 3+ injury risk for each occupant was calculated using AIS 3+ injury risk estimations for the head, chest, abdomen, and lower extremities. The injury risk calculated using computational models was compared with field-based injury risk derived from NASS-CDS by calculating their correlation coefficient. The field-based injury risk was calculated using risk curves that were created based on real-world crash data in a previous study (Hostetler et al., 2020). Occupant age (40 years), seatbelt use (belted occupant), collision deformation classification, lateral ΔV, and PDOF of the crash event were used in these curves to estimate field injury risk. Large differences in the kinematics were observed between HBM and ATD models. ATD models tended to overestimate risk in almost every case whereas HBMs yielded better risk estimates overall. Chest and lower extremity risks were the least correlated with field injury risk estimates. The overall risk of AIS 3+ injury risk was the strongest comparison to the field data-based risk curves. The HBMs were still not able to capture all the variance but future studies can be carried out that are focused on investigating their shortfalls and improving them to estimate injury risk closer to field injury risk in far-side crashes.

Development and preliminary validation of computationally efficient and detailed 50th percentile female human body models

OBJECTIVE

No vehicle testing standard (physical or computational) employs a mid-sized female human surrogate, despite discrepancies related to injury outcomes for female occupants amongst all vehicle users. We detail the design and preliminary validation of 50th percentile female (F50) computational human body models (HBMs) based on Global Human Body Models Consortium (GHBMC) models.

METHOD

Data for the target geometry was collected as part of the initial generation of GHBMC models. Imaging, surface data, and 15 anthropomorphic measures from a living female subject (60.8 kg and 1.61 m) served as the baseline for model development. Due to the role rib cage geometry plays in biomechanical loading, rib cage morphology from secondary retrospective data was leveraged to identify an average female rib cage based on gross anatomical features. A female rib cage was selected from an existing dataset closest to the mean depth, height, and width of the set, considering only those aged 20 - 50 years. The selected subject among this secondary set also exhibited a 7th rib angle and sternum angle within 5% of the mean measurements, and within the range of previously reported studies. The GHBMC 5th percentile, small female detailed (high biofidelity) and simplified (computationally efficient) models were morphed to match the F50 subject body surface, selected bones, and mean rib cage using established thin plate spline techniques. The models were validated vs. previously published literature studies with an emphasis on rib cage response. Model data was compared to 47 channels of experimental data across four biomechanical hub simulations, two sled test simulations (one of which included all female PMHS), and two robustness simulations to test stability. Model results were mass scaled to the average of the reported corridors. Objective evaluation was conducted using CORA. IRB approval was obtained for all prospective and retrospective data collected or used. The target rib cage was selected from retrospective image data used in prior studies (n = 339 chest CT scans).

RESULTS

The morphed HBMs closely matched the target geometry. The detailed and simplified models had masses and element counts of 61.2 kg and 61.8 kg, and 2.8 million and 0.3 million, respectively. The mass difference is due to a coarser mesh in the simplified model. The simplified model ran 23 times faster than the detailed model on the same hardware. Each model exhibited stability in robustness tests, and the average CORA scores were 0.80 and 0.72 in the detailed and simplified models, respectively. The models performed well in frontal impacts against PMHS corridors after mass scaling.

CONCLUSIONS

Numerous recent studies underscore poorer injury outcomes for female vehicle occupants compared to males. While such outcomes are multifactorial, the average female models introduced in this work offer a novel tool within a widely used family of HBMs to reduce the outcome gap in terms of injury for all drivers. HBMs can be deployed in safety studies or in future regulatory requirements faster and more economically than a resized or newly designed ATDs aimed at the same target population.

Development and Validation of an Active Muscle Simplified Finite Element Human Body Model in a Standing Posture

Active muscles play an important role in postural stabilization, and muscle-induced joint stiffening can alter the kinematic response of the human body, particularly that of the lower extremities, under dynamic loading conditions. There are few full-body human body finite element models with active muscles in a standing posture. Thus, the objective of this study was to develop and validate the M50-PS+Active model, an average-male simplified human body model in a standing posture with active musculature. The M50-PS+Active model was developed by incorporating 116 skeletal muscles, as one-dimensional beam elements with a Hill-type material model and closed-loop Proportional Integral Derivative (PID) controller muscle activation strategy, into the Global Human Body Models Consortium (GHBMC) simplified pedestrian model M50-PS. The M50-PS+Active model was first validated in a gravity standing test, showing the effectiveness of the active muscles in maintaining a standing posture under gravitational loading. The knee kinematics of the model were compared against volunteer kinematics in unsuited and suited step-down tests from NASA's active response gravity offload system (ARGOS) laboratory. The M50-PS+Active model showed good biofidelity with volunteer kinematics with an overall CORA score of 0.80, as compared to 0.64 (fair) in the passive M50-PS model. The M50-PS+Active model will serve as a useful tool to study the biomechanics of the human body in vehicle-pedestrian accidents, public transportation braking, and space missions piloted in a standing posture.

Micro-CT Imaging and Mechanical Properties of Ovine Ribs

The use of ovine animal models in the study of injury biomechanics and modeling is increasing, due to their favorable size and other physiological characteristics. Along with this increase, there has also been increased interest in the development of in silico ovine models for computational studies to compliment physical experiments. However, there remains a gap in the literature characterizing the morphological and mechanical characteristics of ovine ribs. The objective of this study therefore is to report anatomical and mechanical properties of the ovine ribs using microtomography (micro-CT) and two types of mechanical testing (quasi-static bending and dynamic tension). Using microtomography, young ovine rib samples obtained from a local abattoir were cut into approximately fourteen 38 mm sections and scanned. From these scans, the cortical bone thickness and cross-sectional area were measured, and the moment of inertia was calculated to enhance the mechanical testing data. Based on a standard least squares statistical model, the cortical bone thickness varied depending on the region of the cross-section and the position along the length of the rib (p < 0.05), whereas the cross-sectional area remained consistent (p > 0.05). Quasi-static three-point bend testing was completed on ovine rib samples, and the resulting force-displacement data was analyzed to obtain the stiffness (44.67 ± 17.65 N/mm), maximum load (170.54 ± 48.28 N) and displacement at maximum load (7.19 ± 2.75 mm), yield load (167.81 ± 48.12 N) and displacement at yield (6.10 ± 2.25 mm), and the failure load (110.90 ± 39.30 N) and displacement at failure (18.43 ± 2.10 mm). The resulting properties were not significantly affected by the rib (p > 0.05), but by the animal they originated from (p < 0.05). For the dynamic testing, samples were cut into coupons and tested in tension with an average strain rate of 18.9 strain/sec. The resulting dynamic testing properties of elastic modulus (5.16 ± 2.03 GPa), failure stress (63.29 ± 14.02 MPa), and failure strain (0.0201 ± 0.0052) did not vary based on loading rate (p > 0.05).

Sensitivity Analysis for Multidirectional Spaceflight Loading and Muscle Deconditioning on Astronaut Response

A sensitivity analysis for loading conditions and muscle deconditioning on astronaut response for spaceflight transient accelerations was carried out using a mid-size male human body model with active musculature. The model was validated in spaceflight-relevant 2.5-15 g loading magnitudes in seven volunteer tests, showing good biofidelity (CORA: 0.69). Sensitivity analysis was carried out in simulations varying pulse magnitude (5, 10, and 15 g), rise time (32.5 and 120 ms), and direction (10 directions: frontal, rear, vertical, lateral, and their combination) along with muscle size change (± 15% change) and responsiveness (pre-braced, relaxed, vs. delayed response) changes across 600 simulations. Injury metrics were most sensitive to the loading direction (50%, partial-R²) and least sensitive to muscle size changes (0.2%). The pulse magnitude also had significant effect on the injury metrics (16%), whereas muscle responsiveness (3%) and pulse rise time (2%) had only slight effects. Frontal and upward loading directions were the worst for neck, spine, and lower extremity injury metrics, whereas rear and downward directions were the worst for head injury metrics. Higher magnitude pulses and pre-bracing also increased the injury risk.

Quantifying Cardiothoracic Variation with Posture and Respiration to Inform Cardiac Device Design

PURPOSE

With extravascular implantable cardioverter defibrillator leads placed beneath the sternum, it is important to quantify heart motion relative to the rib cage with postural changes and respiration.

METHODS

MRI scans from five males and five females were collected in upright and supine postures at end inspiration [n = 10 each]. Left and right decubitus [n = 8 each] and prone [n = 5] MRIs at end inspiration and supine MRIs at end expiration [n = 5] were collected on a subset. Four cardiothoracic measurements, six cardiac measurements, and six cardiac landmarks were collected to measure changes across different postures and stages of respiration.

RESULTS

The relative location of the LV apex to the nearest intercostal space was significantly different between the supine and decubitus postures (average ± SD difference: - 15.7 ± 11.4 mm; p < 0.05). The heart centroid to xipho-sternal junction distance was 9.7 ± 7.9 mm greater in the supine posture when compared to the upright posture (p < 0.05). Cardiac landmark motion in the lateral direction was largest due to postural movement (range 23-50 mm) from the left decubitus to the right decubitus posture, and less influenced by respiration (5-17 mm). Caudal-cranial displacement was generally larger due to upright posture (13-23 mm caudal) and inspiration (7-20 mm cranial).

CONCLUSIONS

This study demonstrates that the location of the heart with respect to the rib cage varies with posture and respiration. The gravitational effects of postural shifts on the heart position are roughly 2-3 times larger than the effects of normal respiration.

Application of survival analysis to model proliferation likelihood of Escherichia coli biofilm following laser-induced hyperthermia treatment

Eighty percent of bacterial infections associated with living tissue and medical devices are linked to drug-resistant biofilms, leading to lengthy and costly recoveries. Laser-induced hyperthermia can disrupt cell proliferation within biofilms and increase susceptibility to antibiotics. However, there can be bacterial survival differences dependent upon laser irradiation times, and prolonged time at elevated temperature can damage healthy tissue. The objective of this study was to use survival analysis to model the impact of temperature increases on reducing viable biofilm bacteria. In vitro biofilms of Escherichia coli were grown on silicone discs or silicone doped with photothermal poly(3,4-ethylenedioxythiophene) hydrate (PEDOT) nanotubes, and subjected to laser-induced hyperthermia, using a 3 W continuous wave laser at 800 nm for varying times. The number of colony forming units per milliliter (CFU/mL) and maximum temperature were measured after each trial. Survival analysis was employed to estimate bacterial cell proliferation post-treatment to provide a quantitative framework for future studies evaluating photothermal inactivation of bacterial biofilms. The results demonstrate the first application of survival analysis for predicting the likelihood of bacterial cell proliferation based on temperature.

Effect of Active Muscles on Astronaut Kinematics and Injury Risk for Piloted Lunar Landing and Launch While Standing

While astronauts may pilot future lunar landers in a standing posture, the response of the human body under lunar launch and landing-related dynamic loading conditions is not well understood. It is important to consider the effects of active muscles under these loading conditions as muscles stabilize posture while standing. In the present study, astronaut response for a piloted lunar mission in a standing posture was simulated using an active human body model (HBM) with a closed-loop joint-angle based proportional integral derivative controller muscle activation strategy and compared with a passive HBM to understand the effects of active muscles on astronaut body kinematics and injury risk. While head, neck, and lumbar spine injury risk were relatively unaffected by active muscles, the lower extremity injury risk and the head and arm kinematics were significantly changed. Active muscle prevented knee-buckling and spinal slouching and lowered tibia injury risk in the active vs. passive model (revised tibia index: 0.02-0.40 vs. 0.01-0.58; acceptable tolerance: 0.43). Head displacement was higher in the active vs. passive model (11.6 vs. 9.0 cm forward, 6.3 vs. 7.0 cm backward, 7.9 vs. 7.3 cm downward, 3.7 vs. 2.4 cm lateral). Lower arm movement was seen with the active vs. passive model (23 vs. 35 cm backward, 12 vs. 20 cm downward). Overall simulations suggest that the passive model may overpredict injury risk in astronauts for spaceflight loading conditions, which can be improved using the model with active musculature.

Response of small female and midsize male models with active musculature in pre-crash maneuvers and low-speed impacts

OBJECTIVE

The objectives of this study were to evaluate computationally efficient small female (54.1 kg, 149.9 cm) and midsize male (78.4 kg, 174.9 cm) models with active muscles using volunteer sled test data in a frontal-oblique loading direction and check their response in crash mitigating maneuvers using field test data.

METHODS

The Global Human Body Models Consortium small female (F05-OS+Active) and midsize male (M50-OS+Active) simplified occupant models with active musculature were used in this study. The data from a total of 48 previously published sled test experiments were used to simulate a total of 16 simulations. The experimental study recorded occupant responses of six small female and six midsize male volunteers (n = 12 total) in two muscle conditions (relaxed and braced) at two acceleration pulses representing pre-crash braking (1.0 g) and a low-speed impact (2.5 g). Each model's kinematics and reaction forces were compared with experimental data. Along with sled test simulations, both of these models were simulated in abrupt braking, lane change, and turn and brake events using literature data. A total of 36 field test simulations were carried out. A CORA analysis was carried out using reaction load and displacement time-history data for sled test simulations and head CG displacement time-history was used for field test simulations.

RESULTS

The occupant peak forward and lateral excursion results of both active models reasonably matched the volunteer data in the low-speed sled test simulations for both pulse severities. The differences between the active and control models were statistically significant (p-value < 0.05) based on the results of Wilcoxon signed-rank tests using peak forward and lateral excursion data. The average CORA scores calculated for the sled test (sled test: M50-OS+Active= 0.543, male control= 0.471, F05-OS+Active= 0.621, female control= 0.505) and field test (M50-OS+Active= 0.836, male control= 0.466, F05-OS+Active= 0.832, female control= 0.787) simulations were higher for active models than control.

CONCLUSIONS

The responses of the F05-OS+Active and M50-OS+Active models were better than control models based on overall CORA scores calculated using both sled and field tests. The results highlight their ability to predict occupant kinematics in crash-mitigating maneuvers and low-speed impacts in the frontal, lateral and frontal-oblique directions.

Effects of Standing, Upright Seated, vs. Reclined Seated Postures on Astronaut Injury Biomechanics for Lunar Landings

Astronauts may pilot a future lunar lander in a standing or upright/reclined seated posture. This study compared kinematics and injury risk for the upright/reclined (30°; 60°) seated vs. standing postures for lunar launch/landing using human body modeling across 30 simulations. While head metrics for standing and upright seated postures were comparable to 30 cm height jumps, those of reclined postures were closer to 60 cm height jumps. Head linear acceleration for 60° reclined posture in the 5 g/10 ms pulse exceeded NASA's tolerance (10.1 g; tolerance: 10 g). Lower extremity metrics exceeding NASA's tolerance in the standing posture (revised tibia index: 0.36-0.53; tolerance: 0.43) were lowered in seated postures (0.00-0.04). Head displacement was higher in standing vs. seated (9.0 cm vs. 2.4 cm forward, 7.0 cm vs. 1.3 cm backward, 2.1 cm vs. 1.2 cm upward, 7.3 cm vs. 0.8 cm downward, 2.4 cm vs. 3.2 cm lateral). Higher arm movement was seen with seated vs. standing (40 cm vs. 25 cm forward, 60 cm vs. 15 cm upward, 30 cm vs. 20 cm downward). Pulse-nature contributed more than 40% to the injury metrics for seated postures compared to 80% in the standing posture. Seat recline angle contributed about 22% to the injury metrics in the seated posture. This study established a computational methodology to simulate the different postures of an astronaut for lunar landings and generated baseline injury risk and body kinematics data.

Lower Extremity Validation of a Human Body Model for High Rate Axial Loading in the Underbody Blast Environment

While the use of Human Body Models (HBMs) in the underbody blast (UBB) environment has increased and shown positive results, the potential of these models has not been fully explored. Obtaining accurate kinematic and kinetic response are necessary to better understand the injury mechanisms for military safety applications. The objective of this study was to validate the Global Human Body Models Consortium (GHBMC) M50 lower extremity using a combined objective rating scheme in vertical and horizontal high-rate axial loading. The model's lower extremity biomechanical response was compared to Post Mortem Human Subjects (PMHS) subjects for vertically and horizontally-applied high rate axial loading. Two distinct experimental setups were used for model validation, comprising a total of 33 distinct end points for validation. A combined Correlation and Analysis (CORA) score that incorporates CORA, time-to-peak (TTP) and peak magnitude of the experimental signals and ISO TS 18571 was used to evaluate the model response. For the horizontal impacts, the combined CORA scores were 0.80, 0.84, and 0.81 for compression, force, and strain respectively. For the vertical impacts combined CORA scores for the knee Z force, compression and heel Z displacement ranged from 0.70-0.81, 0.87-0.91, and 0.82-0.99 respectively. The GHBMC lower extremity model showed good agreement with PMHS experimental data in the horizontal and vertical loading environment in 33 unique tests. The accuracy is demonstrated by using the ISO TS 18571 standard and a combined CORA score that takes into consideration the peak and time to peak of the signal. The results of this study show that GHBMC v 6.0 HBM lower extremity can be used for kinetic and kinematic predictions in the UBB environment.

Implementation and calibration of active small female and average male human body models using low-speed frontal sled tests

OBJECTIVE

The objective of this study was to implement active muscles in a computationally efficient small female finite element model (54.1 kg, 149.9 cm) suitable for predicting occupant response during precrash braking and low-speed frontal sled tests. We further calibrate and compare its results against an average male model (78.4 kg, 174.9 cm) using the same developmental approach.

METHODS

The active female model (F05-OS + Active) was developed by adding active skeletal muscle elements (n = 232) to the Global Human Body Models Consortium (GHBMC) 5th percentile female simplified occupant model (F05-OS v2.3). The muscle properties and physiological cross-sectional area (PCSA) for each muscle were taken from the M50-OS + Active v2.3 model but PCSAs were mass scaled to a 5th percentile female. A total of 8 simulations were conducted; 2 acceleration pulses (1.0 g and 2.5 g), 2 models (F05-OS + Active and M50-OS + Active), and 2 muscle states (activation and control; e.g., no activation). Each model's kinematics and reaction forces were compared with experimental data. Occupant responses of 6 5th percentile female and 6 50th percentile male volunteers (n = 12 total) were used. The data depict occupant response in precrash braking and low-speed frontal sled tests in a rigid test buck. All procedures were reviewed and approved by the Virginia Tech institutional review board. Each volunteer was in a relaxed state before the applied acceleration.

RESULTS

The occupant peak forward excursion results of both active models reasonably match the volunteer data for both pulse severities. The differences between active and control models were found to be significant by Wilcoxon signed-rank test (p < .05). The reaction loads of the active and control models lie within the experimental corridors.

CONCLUSIONS

To the authors' knowledge, this study is the first to concurrently calibrate and compare equivalently developed computational models of females and males in precrash and low-speed impacts. The modeling approach is capable of capturing the varied kinematics observed in the relaxed condition, which may be an important factor in studies focused on the effects of low-g vehicle dynamics on the occupant position. Finally, the computationally efficient modeling approach is imperative given the long duration (>500 ms) of the events simulated.

Age targeted human body models indicate increased thoracic injury risk with aging

OBJECTIVE

The objective of this study is to generate age targeted versions of the male and female Global Human Body Models Consortium (GHBMC) occupant human body models (HBMs), to validate each in frontal impacts, and to assess rib fracture probability of each.

METHODS

Six age targeted models were developed based on the GHBMC average male and small female occupant models (M50-O v6.0 and F05-O v6.0, respectively). All age targeted models were modified to represent population means for height, weight, shape, and relevant material properties. The thin plate spline method was used to morph models, and material properties were modified using available literature. Validation focused on chest response. Models were evaluated in a rigid body frontal chest impact at 6.7 m/s. Furthermore, the male and female age targeted models were evaluated against published data from 40 km/hr and 30 km/hr frontal sled tests respectively.

RESULTS

Chest deflections and landmark kinematics reasonably matched the respective corridors in the M50-O and F05-O aged models. Regional probability of rib fracture was assessed using probabilistic methods based on cortex strain. Increasing rib fracture with age was observed in both impacts for both sexes. For the rigid chest impact, the M50-O 70YO resulted in 10 ribs exceeding 50% probability of fracture whereas the younger ages reported 4 to 6 ribs exceeding the same probability. In the same simulation, the F05-O 70YO resulted in 8 regions exceeding 50% probability of rib fracture as opposed to 3 and 0 such regions at the youngest ages. Sled simulation demonstrated similar trends. The 70YO age adjusted models best aligned with the reported extent of fractures from the referenced PMHS studies, which tend to be composed of subjects of advanced age.

CONCLUSIONS

Age targeted HBMs demonstrated increased fracture probability with age when subjected to equivalent impacts. Gross model kinematics approximate PMHS data but showed little difference between targeted age models. The findings indicate that while gross kinematics are unaffected by age-targeting models, such models can capture trends of increased thoracic injury risk observed in experimental and field studies, and further suggest their potential use to target interventions for vulnerable driving populations, such as older adults.

Development of a computational fluid dynamic model to investigate the hemodynamic impact of REBOA

Background: Resuscitative endovascular balloon occlusion of the aorta (REBOA) is a lifesaving intervention for major truncal hemorrhage. Balloon-tipped arterial catheters are inserted via the femoral artery to create a temporary occlusion of the aorta, which minimizes the rate of internal bleeding until definitive surgery can be conducted. There is growing concern over the resultant hypoperfusion and potential damage to tissues and organs downstream of REBOA. To better understand the acute hemodynamic changes imposed by REBOA, we developed a three-dimensional computational fluid dynamic (CFD) model under normal, hemorrhage, and aortic occlusion conditions. The goal was to characterize the acute hemodynamic changes and identify regions within the aortic vascular tree susceptible to abnormal flow and shear stress. Methods: Hemodynamic data from established porcine hemorrhage models were used to build a CFD model. Swine underwent 20% controlled hemorrhage and were randomized to receive a full or partial aortic occlusion. Using CT scans, we generated a pig-specific aortic geometry and imposed physiologically relevant inlet flow and outlet pressure boundary conditions to match in vivo data. By assuming non-Newtonian fluid properties, pressure, velocity, and shear stresses were quantified over a cardiac cycle. Results: We observed a significant rise in blood pressure (∼147 mmHg) proximal to REBOA, which resulted in increased flow and shear stress within the ascending aorta. Specifically, we observed high levels of shear stress within the subclavian arteries (22.75 Pa). Alternatively, at the site of full REBOA, wall shear stress was low (0.04 ± 9.07E-4 Pa), but flow oscillations were high (oscillatory shear index of 0.31). Comparatively, partial REBOA elevated shear levels to 84.14 ± 19.50 Pa and reduced flow oscillations. Our numerical simulations were congruent within 5% of averaged porcine experimental data over a cardiac cycle. Conclusion: This CFD model is the first to our knowledge to quantify the acute hemodynamic changes imposed by REBOA. We identified areas of low shear stress near the site of occlusion and high shear stress in the subclavian arteries. Future studies are needed to determine the optimal design parameters of endovascular hemorrhage control devices that can minimize flow perturbations and areas of high shear.

Development and implementation of a time- and computationally-efficient methodology for reconstructing real-world crashes using finite element modeling to improve crash injury research investigations

Eleven Crash Injury Research and Engineering Network (CIREN) frontal crashes were reconstructed using a novel, time-efficient methodology involving a simplified vehicle model. Kinematic accuracy was assessed using novel kinematic scores between 0-1 and chest injury was assessed using literature-defined injury metric time histories. The average kinematic score across all simulations was 0.87, indicating good kinematic accuracy. Time histories for chest compression, rib strain, shoulder belt force, and steering column force discerned the most causative components of chest injury in all cases. Abbreviated Injury Scale (AIS) 2+ and AIS 3+ chest injury risk functions using belt force identified chest injury with 81.8% success.

Effect of body size and enhanced helmet systems on risk for motorsport drivers

OBJECTIVE

Computational modeling has been shown to be a useful tool for simulating representative motorsport impacts and analyzing data for relative injury risk assessment. Previous studies have used computational modeling to analyze the probability of injury in specific regions of a 50th percentile male driver. However, NASCAR drivers can represent a large range in terms of size and female drivers are becoming increasingly more common in the sport. Additionally, motorsport helmets can be outfitted with external attachments, or enhanced helmet systems (EHS), whose effect is unknown relative to head and neck kinematics. The current study expands on this previous work by incorporating the F05-OS and M95-OS into the motorsport environment in order to determine correlations between metrics and factors such as PDOF, resultant ΔV occupant size, and EHS.

METHODS

A multi-step computational process was used to integrate the Global Human Body Models Consortium family of simplified occupant models into a motorsport environment. This family included the 5th percentile female (F05-OS), 50th percentile male (M50-OS), and 95th percentile male (M95-OS), which provide a representative range for the size and sex of drivers seen in NASCAR's racing series'. A series of 45 representative impacts, developed from real-world crash data, and set of observed on-track severe impacts were conducted with these models. These impacts were run in triplicate for three helmet configurations: bare helmet, helmet with visor, helmet with visor and camera. This resulted in 450 total simulations. A paired t-test was initially performed as an exploratory analysis to study the effect of helmet configuration on 10 head and neck injury metrics. A mixed-effects model with unstructured covariance matrix was then utilized to correlate the effect between five independent variables (resultant ΔV, body size, helmet configuration, impact quadrant, and steering wheel position) and a selection of 25 metrics. All simulations were conducted in LS-Dyna R. 9.1.

RESULTS

Risk estimates from the M50-OS with bare helmet were used as reference values to determine the effect of body size and helmet configuration. The paired t-test found significance for helmet configuration in select head-neck metrics, but the relative increase in these metrics was low and not likely to increase injury risk. The mixed-effects model analyzed statistical relationships across multiple types of variables. Within the mixed-effects model, no significance was found between helmet configuration and metrics. The greatest effect was found from resultant ΔV, body size, and impact quadrant.

CONCLUSIONS

Overall, smaller drivers showed statistically significant reductions in injury metrics, while larger drivers showed statistically significant increases. Lateral impacts showed the greatest effect on neck metrics and, on average, showed decreases for head metrics related to linear acceleration and increases for head metrics related to angular velocity. HBM parametric studies such as this may provide an avenue to assist injury detection for motorsport incidents, improve triage effectiveness, and assist in the development of safety standards.

Deployment and testing of an automated medical equipment service communication and documentation system at a rural hospital in Kenya

Background

The Medical Equipment Network Documentation System (MENDS) provides a simple communication network for equipment servicing from failure to restoration. It is a text messaging-based platform, designed to use existing technologies in place in low- and middle-income settings. The system gathers and relays information about equipment service requests and reports and automatically saves them to an online database.

Methods

MENDS was deployed at a high volume, rural, charity medical facility in Kijabe, Kenya for a 3-mo pilot test.

Results

The results show MENDS more than tripled documentation and enhanced ease and speed of communication.

Conclusions

Comprehensive data provided by MENDS created more accurate measures of equipment performance, which can be used to decrease the time that equipment is out of service and improve the efficiency of repairs, equipment quality and procurement.

A detailed finite element model of a mid-sized male for the investigation of traffic pedestrian accidents

The pedestrian is one of the most vulnerable road users and comprises approximately 23% of the road crash-related fatalities in the world. To protect pedestrians during Car-to-Pedestrian Collisions (CPC), subsystem impact tests are used in regulations. These tests provide insight but cannot characterize the complex vehicle-pedestrian interaction. The main purpose of this study was to develop and validate a detailed pedestrian Finite Element (FE) model corresponding to a 50th percentile male to predict CPC induced injuries. The model geometry was reconstructed using a multi-modality protocol from medical images and exterior scan data corresponding to a mid-sized male volunteer. To investigate injury response, this model included internal organs, muscles and vessels. The lower extremity, shoulder and upper body of the model were validated against Post Mortem Human Surrogate (PMHS) test data in valgus bending, and lateral/anterior-lateral blunt impacts, respectively. The whole-body pedestrian model was validated in CPC simulations using a mid-sized sedan and simplified generic vehicles bucks and previously unpublished PMHS coronal knee angle data. In the component validations, the responses of the FE model were mostly within PMHS test corridors and in whole body validations the kinematic and injury responses predicted by the model showed similar trends to PMHS test data. Overall, the detailed model showed higher biofidelity, especially in the upper body regions, compared to a previously reported simplified pedestrian model, which recommends using it in future pedestrian automotive safety research.

An Improved Method for Developing Injury Risk Curves Using the Brier Metric Score

Many injury metrics are routinely proposed from measured or derived quantities from biomechanical experiments using post mortem human subjects (PMHS). The existing literature did not provide guidance on deciding between parameters collected in an experiment that would be best to use for the development of human injury probability curves (HIPC). The objective of this study was to use the Brier Metric Score (BMS) to identify the most appropriate metric from an experiment that predicts injury outcomes. The Brier Metric Score assesses how well a metric predicts the outcome for a censored data point (a lower BMS is better). Survival analysis was then conducted with the selected metric and the best distribution was selected using Akaike information criterion (AIC). Confidence intervals (CIs) and the normalized confidence interval width (NCIS) were calculated for the injury probability curve. The testing and validation of the methods described were performed using biomechanics data in the open literature. The methods for the HIPC development procedure detailed herein have been rigorously tested and used in the generation of WIAMan HIPCs and Injury Assessment Reference Curves (IARCs) for the WIAMan ATD, but can also be used in other ATD or PMHS injury risk curve development.

Injury risk curves in far-side lateral motor vehicle crashes by AIS level, body region and injury code

OBJECTIVE

The objective of this study was to develop injury risk curves as a function of change in vehicle velocity for occupants in far-side lateral motor vehicle crashes (MVCs) by AIS level, body region, and specific AIS codes that commonly occur in this crash mode.

METHODS

The National Automotive Sampling System-Crashworthiness Data System (NASS-CDS) years 2000-2015 database was queried, resulting in 4,495 non-weighted far-side crashes. For each case, occupant age, sex, and the following metadata were collected: vehicle model year, vehicle body type, lateral delta-v, normalized PDOF, multiple impacts, belt use, seat position, object contacted, striking vehicle body type, maximum crush extent and side airbag deployment. Multivariable logistic regression was used to develop risk curves for AIS 2+ through 5+ injuries, AIS 2+ injuries by body region (head, thorax, lower extremity), and for each of the 10 most frequent far-side AIS 2+ injuries. Significant covariates were determined by backwards elimination (p < 0.05). The full dataset and a subsampled dataset of only cases with side airbag deployment were used to develop risk curves.

RESULTS

For AIS 2+ through 5+ injury, greater delta-V was associated with greater injury risk (OR's: 2.48-3.66 per 11.9 kph increase) and belt use was associated with lower risk (OR's: 0.04-0.36 compared to unbelted). Multiple impacts were significant predictors of increased AIS 3+, 4+ and 5+ injury risk (OR's: 2.56, 2.27 and 2.83 compared to single impact). For AIS 2+ body region injuries, lateral delta-V and maximum CDC extent were positively associated with increased head, thorax and lower extremity injury risk while belt use was associated with lower risk. Increased lateral delta-v, unbelted status, and greater maximum CDC extent frequently increased injury risk for the most common far-side injuries. Side airbag deployment was not a significant covariate for the injury risk models.

CONCLUSIONS

The resulting risk models expand upon previous literature gaps to provide a more comprehensive view of contributors to injury risk for occupants in far-side MVCs. This study yields risk curves based on the latest available NASS-CDS data.

Simulation-based assessment of injury risk for an average male motorsport driver

OBJECTIVE

While well-protected through a variety of safety countermeasures, motorsports drivers can be exposed to a large variety of crash modes and severities. Computational human body models (HBMs) are currently used to assess occupant safety for the general driving public in production vehicles. The purpose of this study was to incorporate a HBM into a motorsport environment using a simulation-based approach and provide quantitative data on relative risk for on-track motorsport crashes.

METHODS

Unlike a traditional automotive seat, the NASCAR driver environment is driver-customized and form-fitting. A multi-step process was developed to integrate the Global Human Body Models Consortium (GHBMC) 50th percentile male simplified occupant into a representative motorsport environment which includes a donned helmet, a 7-point safety belt system, head and neck restraint (HNR), poured-foam seat, steering wheel, and leg enclosure. A series of 45 representative impacts, developed from real-world crash data, of varying severity (10 kph ≤ ΔV ≤ 100 kph) and impact direction (∼290° ≤ PDOF ≤ 20°) were conducted with the GHBMC 50th percentile male simplified occupant (M50-OS v2.2). Kinematic and kinetic data, and a variety of injury criteria, were output from each of the simulations and used to calculate AIS 1+, 2+, and 3+ injury risk. All simulations were conducted in LS-Dyna R. 9.1.

RESULTS

Injury risk of the occupant using the previously mentioned injury criteria was calculated for the head, neck, thorax, and lower extremity, and the probability of injury for each region was plotted. Of the simulated impacts, five had a maximum AIS 1+ injury risk >20%, six had a maximum AIS 2+ injury risk >10%, and no cases had a maximum AIS 3+ injury >1%. Overall, injury risk estimates were reasonable compared to on-track data reported from Patalak et al. (2020).

CONCLUSIONS

Beyond injury risk, the study is the first of its kind to provide mechanical loading values likely experienced during motorsports crash incidents with crash pulses developed from real-world data. Given the severity of the crash pulses, the simulated environments reinforce the need for the robust safety environment implemented by NASCAR.

Comparison of Neck Injury Criteria Values Across Human Body Models of Varying Complexity

Due to the severity and frequency of cervical spine injuries, the neck injury criterion (Nij) was developed to provide a quantitative relationship between forces and moments of the upper neck with accompanied injury risk. The present study was undertaken to evaluate differences in calculated Nij for the Global Human Body Model Consortium's detailed and simplified average 50th percentile male models. The simplified model is a computationally light version of the more detailed model and therefore it is of interest to achieve similar Nij values between the two models. These models were compared in 15 match paired conditions of rigid head impact and a mixture of seven full body rigid hub and sled pulses, for 44 total simulations. Collectively, Nij values of the simplified model were found to exhibit a second-degree polynomial fit, allowing for a conversion to the prediction of the detailed model. Correlates were also derived for impact and inertial loading cases individually, for which the latter may be the subject of future work. The differences in Nij may be attributed to a variety of modeling approach differences related to neck muscles (attachment location and morphometric implementation), localization of head mass within the M50-OS, head geometry, and intervertebral joint space properties. With a primary focus on configurations in the anterior-posterior direction, there is a potential limitation in extensibility to lateral loading cases. In response to the relatively low Nij values exhibited, future work should evaluate the appropriateness of the established critical intercepts of Nij for computational human body models.

Lumbar Spine Response of Computational Finite Element Models in Multidirectional Spaceflight Landing Conditions

The goals of this study are to compare the lumbar spine response variance between the hybrid III, test device for human occupant restraint (THOR), and global human body models consortium simplified 50th percentile (GHBMC M50-OS) finite element models and evaluate the sensitivity of lumbar spine injury metrics to multidirectional acceleration pulses for spaceflight landing conditions. The hybrid III, THOR, and GHBMC models were positioned in a baseline posture within a generic seat with side guards and a five-point restraint system. Thirteen boundary conditions, which were categorized as loading condition variables and environmental variables, were included in the parametric study using a Latin hypercube design of experiments. Each of the three models underwent 455 simulations for a total of 1365 simulations. The hybrid III and THOR models exhibited similar lumbar compression forces. The average lumbar compression force was 45% higher for hybrid III (2.2 ± 1.5 kN) and 51% higher for THOR (2.0 ± 1.6 kN) compared to GHBMC (1.3 ± 0.9 kN). Compared to hybrid III, THOR sustained an average 64% higher lumbar flexion moment and an average 436% higher lumbar extension moment. The GHBMC model sustained much lower bending moments compared to hybrid III and THOR. Regressions revealed that lumbar spine responses were more sensitive to loading condition variables than environmental variables across all models. This study quantified the intermodel lumbar spine response variations and sensitivity between hybrid III, THOR, and GHBMC. Results improve the understanding of lumbar spine response in spaceflight landings.

Head injury metric response in finite element ATDs and a human body model in multidirectional loading regimes

Objective: The objective was to quantify head injury metric sensitivity of the 50th percentile male Hybrid III, THOR, and Global Human Body Models Consortium simplified occupant (GHBMC M50-OS) to changes in loading conditions in loading regimes that may be experienced by occupants of spaceflight vehicles or highly autonomous vehicles (HAVs) with nontraditional seating configurations.Methods: A Latin hypercube (LHD) design of experiments (DOE) was employed to develop boundary conditions for 455 unique acceleration profiles. Three previously validated finite element (FE) models of the Hybrid III anthropomorphic test device (ATD), THOR ATD, and GHBMC M50-OS were positioned in an upright 90°-90°-90° seat and with a 5-point belt. Acceleration pulses were applied to each of the three occupants in the ± X, +Y, and ± Z directions, with peak resultant acceleration magnitudes ranging from 5 to 20 G and times to peak ranging from 32.5 to 120.8 ms with duration 250 ms, resulting in 1,248 simulations. Head injury metrics included peak linear head acceleration, peak rotational head acceleration, head injury criteria (HIC15), and brain injury criteria (BrIC). Injury metrics were regressed against boundary condition parameters using 2nd order multiple polynomial regression, and compared between occupants using matched pairs Wilcoxon signed rank analysis.Results: Across the 416 matched-simulations that reached normal termination with all three models, HIC15 values ranged from 1.0-396.5 (Hybrid III), 1.2-327.9 (THOR), and 0.6-585.6 (GHBMC). BrIC ranged from 0.03-0.95 (Hybrid III), 0.03-1.21 (THOR), and 0.04-0.84 (GHBMC). Wilcoxon signed rank analysis demonstrated significant pairwise differences between each of the three occupant models for head injury metrics. For HIC15, the largest divergence between GHBMC and the ATDs was observed in simulations with components of combined underbody and rear impact loading. The three models performed most similarly with respect to BrIC output when loaded in a frontal direction. Both the GHBMC and the Hybrid III produced lower values of BrIC than the THOR on average, with the differences most pronounced in rear impact loading.Conclusion: In conclusion, observed differences between the occupant models' head injury metric output were quantified. Loading direction had a large effect on metric outcome and metric comparability across models, with frontal and rear impacts with low vertical acceleration tending to be the most similar. One explanation for these differences could be the differences in neck stiffness between the models that allowed more rotation in the GHBMC and THOR. Care should be taken when using ATDs as human volunteer surrogates in these low energy events.

Validation of a simplified human body model in relaxed and braced conditions in low-speed frontal sled tests

Objective: The goal of this study was to implement active musculature into the Global Human Body Models Consortium (GHBMC) average male simplified occupant model (M50-OS v2) and validate its performance in low-speed frontal crash scenarios.Methods: Volunteer and postmortem human subjects (PMHS) data from low-speed frontal sled tests by Beeman et al., including 2.5 and 5.0 g acceleration pulses, were used to simulate events in LS-DYNA. All muscles were modeled as 1D beam elements and assigned a Hill-type muscle material. From the output of proportional-integral-derivative (PID) controllers, the activation level for each muscle was calculated using a sigmoid function, representing the firing rate of motor neurons. The PID controller attempts to preserve the initial posture of the model. Percentage muscle contribution for all skeletal muscles was precalculated using the M50-OS with active muscles (M50-OS + Active). The M50-OS + Active employs varying levels of neural delays to represent volunteer relaxed and braced conditions, taken from literature. Braced condition experiments were simulated using elevated joint angle set values for the PID controller. The M50-OS + Active model was used to simulate 2 muscle conditions (relaxed and braced) at 2 pulse severities (2.5 and 5.0 g). A control set of simulations was conducted to compare the effect of adding active muscle. Ten whole-body simulations were conducted.Results: The results from volunteer simulations showed a strong dependence of reaction loads and kinematics on muscle activation. Compared to baseline, M50-OS, at 5.0 g acceleration, 33.3% and 7.6% decreases were observed in the overall head kinematics of the M50-OS + Active for the braced and relaxed conditions, respectively. Regarding the anterior direction, similar reductions in overall kinematics were observed for both volunteer test conditions. In comparison to control simulations in which no active muscle was implemented, objective evaluation scores increased markedly at both speeds for the braced condition. Little to no gain was found in the relaxed condition.Conclusions: The results justify the need for use of an active human body model for predicting low-speed frontal kinematics, particularly in the braced condition. Head kinematics were reduced when using active modeling for all simulations (braced and relaxed).

Validated thoracic vertebrae and costovertebral joints increase biofidelity of a human body model in hub impacts

Objective: A recent emphasis on nontraditional seating and omnidirectional impact directions has motivated the need for deformable representation of the thoracic spine (T-spine) in human body models. The goal of this study was to develop and validate a deformable T-spine for the Global Human Body Models Consortium (GHBMC) M50-O (average male occupant) human model and to demonstrate improved biofidelity.Methods: Eleven functional spinal units (FSUs) were developed with deformable vertebrae (cortical and trabecular), spinal and costovertebral ligaments, and intervertebral discs. Material properties for all parts were obtained from the literature.FSUs were subjected to quasistatic loads per Panjabi et al. (1976) in 6 degrees of freedom. Stiffness values were calculated for each moment (Nm/°) and translational force (N/µm). Updated costovertebral (CV) joints of ribs 2, 6, and 10 were subjected to moments along 3 axes per Duprey et al. (2010). The response was optimized by maximum force and laxity in the ligaments. In both cases, updated models were compared to the baseline approach, which employed rigid bodies and joint-like behavior. The deformable T-spine and CV joints were integrated into the full M50-O model Ver. 5.0β and 2 full-body cases were run: (1) a rear pendulum impact per Forman et al. (2015) at speeds up to 5.5 m/s. and (2) a lateral shoulder impact per Koh (2005) at 4.5 m/s. Quantitative evaluation protocols were used to evaluate the time history response vs. experimental data, with an average correlation and analysis (CORA) score of 0.76.Results: All FSU responses showed reduced stiffness vs. baseline. Tension, extension, torsion, and lateral bending became more compliant than experimental data. Like the experimental results, no trend was observed for joint response by level. CV joints showed good biofidelity. The response at ribs 2, 6, and 10 generally followed the experimental data.Conclusions: Deformable T-spine and CV joint validation has not been previously published and yielded high biofidelity in rear impact and notable improvement in lateral impact at the full body level. Future work will focus on localized T-spine injury criteria made possible by the introduction of this fully deformable representation of the anatomy.

Evaluation of finite element human body models for use in a standardized protocol for pedestrian safety assessment

Objective: Finite element human body models (HBMs) must be certified for use within the EuroNCAP pedestrian safety assessment protocol. We demonstrate that the Global Human Body Model Consortium (GHBMC) simplified pedestrian series of HBMs meet all criteria set forth in Technical Bulletin (TB) 024 (v 1.1 Jan. 2019) for model certification. We further explore variation in head contact time (HIT) and location by HBM size and impact speed across 48 full body impact simulations.Methods: The EuroNCAP Pedestrian Protocol (v. 8.5, Oct. 2018) assesses the overall safety of adult and child pedestrians by outlining a variety of physical tests and finite element simulations using HBMs. These tests are designed to assess the efficacy of vehicle safety technology such as active bonnets. The 50th percentile male simplified pedestrian model (M50-PS, H:175 cm, W:74.5 kg), six-year-old (6YO-PS, H:117 cm, W:23.4 kg), 5th percentile female (F05-PS, H:150 cm, W:50.7 kg), and 95th percentile male (M95-PS, H:190 cm, W:102 kg) were simulated through the suite of cases totaling 48 simulations (12 each). The process gathers three kinematic trajectories and contact force data from designated anatomical locations. The impacting vehicles include a family car (FCR), multi-purpose vehicle (MPV), roadster (RDS), and sports utility vehicle (SUV), all provided by TU Graz, Vehicle Safety Institute as part of the Coherent Project, each simulated at 30 kph, 40 kph, and 50 kph. Each simulation underwent a 23-point pre-simulation check and post-simulation model response comparison. The posture of all HBMs met criteria consisting of 15 measures. All simulations were conducted in LS-Dyna R. 7.1.2.Results and Conclusions: All simulations normal terminated. For each of the simulations, sagittal plane coordinate histories of the center of the head, 12th thoracic vertebrae, and center of acetabulum were compared with standard corridors and did not exceed the tolerance of 50 mm deviation. Head contact time was also compared with the reference values and did not exceed the tolerance interval of +3.5% and -7%. Comparison of contact forces was required for monitoring purposes only. The head contact time of the models for each simulation was recorded and compared by model size, impact speed, and vehicle geometry. Head contact times varied by roughly 3-fold, were lowest for the child model, and showed the greatest sensitivity for the tallest stature model (M95-PS). As stated in the certification process, other body sizes within a model family qualify for certification if the 50th percentile male model passes, provided that model sizes meet the required posture.

Modeling Human Volunteers in Multidirectional, Uni-axial Sled Tests Using a Finite Element Human Body Model

A goal of the Human Research Program at National Aeronautics and Space Administration (NASA) is to analyze and mitigate the risk of occupant injury due to dynamic loads. Experimental tests of human subjects and biofidelic anthropomorphic test devices provide valuable kinematic and kinetic data related to injury risk exposure. However, these experiments are expensive and time consuming compared to computational simulations of similar impact events. This study aimed to simulate human volunteer biodynamic response to unidirectional accelerative loading. Data from seven experimental studies involving 212 volunteer tests performed at the Air Force Research Laboratory were used to reconstruct 13 unique loading conditions across four different loading directions using finite element human body model (HBM) simulations. Acceleration pulses and boundary conditions from the experimental tests were applied to the Global Human Body Models Consortium (GHBMC) simplified 50th percentile male occupant (M50-OS) using the LS-Dyna finite element solver. Head acceleration, chest acceleration, and seat belt force traces were compared between the experimental and matched simulation signals using correlation and analysis (CORA) software and averaged into a comprehensive response score ranging from 0 to 1 with 1 representing a perfect match. The mean comprehensive response scores were 0.689 ± 0.018 (mean ± 1 standard deviation) in two frontal simulations, 0.683 ± 0.060 in four rear simulations, 0.676 ± 0.043 in five lateral simulations, and 0.774 ± 0.013 in two vertical simulations. The CORA scores for head and chest accelerations in these simulations exceeded mean scores reported in the original development and validation of the GHBMC M50-OS model. Collectively, the CORA scores indicated that the HBM in these boundary conditions closely replicated the kinematics of the human volunteers across all loading directions.

Applying dynamic contrast enhanced MSOT imaging to intratumoral pharmacokinetic modeling

Examining the dynamics of an agent in the tumor microenvironment can offer critical insights to the influx rate and accumulation of the agent. Intratumoral kinetic characterization in the in vivo setting can further elicudate distribution patterns and tumor microenvironment. Dynamic contrast-enhanced Multispectral Optoacoustic Tomographic imaging (DCE-MSOT) acquires serial MSOT images with the administration of an exogenous contrast agent over time. We tracked the dynamics of a tumor-targeted contrast agent, HypoxiSense 680 (HS680), in breast xenograft mouse models using MSOT. Arterial input function (AIF) approach with MSOT imaging allowed for tracking HS680 dynamics within the mouse. The optoacoustic signal for HS680 was quantified using the ROI function in the ViewMSOT software. A two-compartment pharmacokinetics (PK) model constructed in MATLAB to fit rate parameters. The contrast influx (k_in) and outflux (k_out) rate constants predicted are k_in = 1.96 × 10^-2 s^-1 and k_out = 9.5 × 10^-3 s^-1 (R = 0.9945).

Similitude assessment method for comparing PMHS response data from impact loading across multiple test devices

Biological tissue testing is inherently susceptible to the wide range of variability specimen to specimen. A primary resource for encapsulating this range of variability is the biofidelity response corridor or BRC. In the field of injury biomechanics, BRCs are often used for development and validation of both physical, such as anthropomorphic test devices, and computational models. For the purpose of generating corridors, post-mortem human surrogates were tested across a range of loading conditions relevant to under-body blast events. To sufficiently cover the wide range of input conditions, a relatively small number of tests were performed across a large spread of conditions. The high volume of required testing called for leveraging the capabilities of multiple impact test facilities, all with slight variations in test devices. A method for assessing similitude of responses between test devices was created as a metric for inclusion of a response in the resulting BRC. The goal of this method was to supply a statistically sound, objective method to assess the similitude of an individual response against a set of responses to ensure that the BRC created from the set was affected primarily by biological variability, not anomalies or differences stemming from test devices.

A preliminary study of human model head and neck response to frontal loading in nontraditional occupant seating configurations

OBJECTIVE

Computational human body models (HBMs) are nominally omnidirectional surrogates given their structural basis in human anatomy. As a result, such models are well suited for studies related to occupant safety in anticipated highly automated vehicles (HAVs). We utilize a well-validated HBM to study the head and neck kinematics in simulations of nontraditional occupant seating configurations.

METHODS

The GHBMC M50-O v. 4.4 HBM was gravity settled into a generic seat buck and situated in a seated posture. The model was simulated in angular increments of 15 degrees clockwise from forward facing to rear facing. A pulse of 17.0 kph (NASS median) was used in each to simulate a frontal impact for each of the 13 seating configurations. Belt anchor points were rotated with the seat; the airbag was appropriately powered based on delta-V, and was not used in rear-facing orientations. Neck forces and moments were calculated.

RESULTS

The 30-degree oblique case was found to result in the maximum neck load and sagittal moment, and thus Neck Injury Criteria (NIJ). Neck loads were minimized in the rear facing condition. The moments and loads, however, were greatest in the lateral seating configuration for these frontal crash simulations.

CONCLUSIONS

In a recent policy statement on HAVs, the NHTSA indicated that vehicle manufacturers will be expected to provide countermeasures that will fully protect occupants given any planned seating or interior configurations. Furthermore, the agency indicated that virtual tests using human models could be used to demonstrate such efficacy. While the results presented are only appropriate for comparison within this study, they do indicate that human models provide reasonable biomechanical data for nontraditional occupant seating arrangements.

Failed rib region prediction in a human body model during crash events with precrash braking

OBJECTIVE

The objective of this study is 2-fold. We used a validated human body finite element model to study the predicted chest injury (focusing on rib fracture as a function of element strain) based on varying levels of simulated precrash braking. Furthermore, we compare deterministic and probabilistic methods of rib injury prediction in the computational model.

METHODS

The Global Human Body Models Consortium (GHBMC) M50-O model was gravity settled in the driver position of a generic interior equipped with an advanced 3-point belt and airbag. Twelve cases were investigated with permutations for failure, precrash braking system, and crash severity. The severities used were median (17 kph), severe (34 kph), and New Car Assessment Program (NCAP; 56.4 kph). Cases with failure enabled removed rib cortical bone elements once 1.8% effective plastic strain was exceeded. Alternatively, a probabilistic framework found in the literature was used to predict rib failure. Both the probabilistic and deterministic methods take into consideration location (anterior, lateral, and posterior). The deterministic method is based on a rubric that defines failed rib regions dependent on a threshold for contiguous failed elements. The probabilistic method depends on age-based strain and failure functions.

RESULTS

Kinematics between both methods were similar (peak max deviation: ΔX_head = 17 mm; ΔZ_head = 4 mm; ΔX_thorax = 5 mm; ΔZ_thorax = 1 mm). Seat belt forces at the time of probabilistic failed region initiation were lower than those at deterministic failed region initiation. The probabilistic method for rib fracture predicted more failed regions in the rib (an analog for fracture) than the deterministic method in all but 1 case where they were equal. The failed region patterns between models are similar; however, there are differences that arise due to stress reduced from element elimination that cause probabilistic failed regions to continue to rise after no deterministic failed region would be predicted.

CONCLUSIONS

Both the probabilistic and deterministic methods indicate similar trends with regards to the effect of precrash braking; however, there are tradeoffs. The deterministic failed region method is more spatially sensitive to failure and is more sensitive to belt loads. The probabilistic failed region method allows for increased capability in postprocessing with respect to age. The probabilistic failed region method predicted more failed regions than the deterministic failed region method due to force distribution differences.

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.

Modular use of human body models of varying levels of complexity: Validation of head kinematics

OBJECTIVE

The significant computational resources required to execute detailed human body finite-element models has motivated the development of faster running, simplified models (e.g., GHBMC M50-OS). Previous studies have demonstrated the ability to modularly incorporate the validated GHBMC M50-O brain model into the simplified model (GHBMC M50-OS+B), which allows for localized analysis of the brain in a fraction of the computation time required for the detailed model. The objective of this study is to validate the head and neck kinematics of the GHBMC M50-O and M50-OS (detailed and simplified versions of the same model) against human volunteer test data in frontal and lateral loading. Furthermore, the effect of modular insertion of the detailed brain model into the M50-OS is quantified.

METHODS

Data from the Navy Biodynamics Laboratory (NBDL) human volunteer studies, including a 15g frontal, 8g frontal, and 7g lateral impact, were reconstructed and simulated using LS-DYNA. A five-point restraint system was used for all simulations, and initial positions of the models were matched with volunteer data using settling and positioning techniques. Both the frontal and lateral simulations were run with the M50-O, M50-OS, and M50-OS+B with active musculature for a total of nine runs.

RESULTS

Normalized run times for the various models used in this study were 8.4 min/ms for the M50-O, 0.26 min/ms for the M50-OS, and 0.97 min/ms for the M50-OS+B, a 32- and 9-fold reduction in run time, respectively. Corridors were reanalyzed for head and T1 kinematics from the NBDL studies. Qualitative evaluation of head rotational accelerations and linear resultant acceleration, as well as linear resultant T1 acceleration, showed reasonable results between all models and the experimental data. Objective evaluation of the results for head center of gravity (CG) accelerations was completed via ISO TS 18571, and indicated scores of 0.673 (M50-O), 0.638 (M50-OS), and 0.656 (M50-OS+B) for the 15g frontal impact. Scores at lower g levels yielded similar results, 0.667 (M50-O), 0.675 (M50-OS), and 0.710 (M50-OS+B) for the 8g frontal impact. The 7g lateral simulations also compared fairly with an average ISO score of 0.565 for the M50-O, 0.634 for the M50-OS, and 0.606 for the M50-OS+B. The three HBMs experienced similar head and neck motion in the frontal simulations, but the M50-O predicted significantly greater head rotation in the lateral simulation.

CONCLUSION

The greatest departure from the detailed occupant models were noted in lateral flexion, potentially indicating the need for further study. Precise modeling of the belt system however was limited by available data. A sensitivity study of these parameters in the frontal condition showed that belt slack and muscle activation have a modest effect on the ISO score. The reduction in computation time of the M50-OS+B reduces the burden of high computational requirements when handling detailed HBMs. Future work will focus on harmonizing the lateral head response of the models and studying localized injury criteria within the brain from the M50-O and M50-OS+B.

Development and Validation of a Brain Phantom for Therapeutic Cooling Devices

Tissue cooling has been proven as a viable therapy for multiple conditions and injuries and has been applied to the brain to treat epilepsy and concussions, leading to improved long-term outcomes. To facilitate the study of temperature reduction as a function of various cooling methods, a thermal brain phantom was developed and analyzed. The phantom is composed of a potassium-neutralized, superabsorbent copolymer hydrogel. The phantom was tested in a series of cooling trials using a cooling block and 37 deg water representing nondirectional blood flow ranging up to 6 gph, a physiologically representative range based on the prototype volume. Results were compared against a validated finite difference (FD) model. Two sets of parameters were used in the FD model: one set to represent the phantom itself and a second set to represent brain parenchyma. The model was then used to calculate steady-state cooling at a depth of 5 mm for all flow rates, for both the phantom and a model of the brain. This effort was undertaken to (1) validate the FD model against the phantom results and (2) evaluate how similar the thermal response of the phantom is to that of a perfused brain. The FD phantom model showed good agreement with the empirical phantom results. Furthermore, the empirical phantom agreed with the predicted brain response within 3.5% at physiological flow, suggesting a biofidelic thermal response. The phantom will be used as a platform for future studies of thermally mediated therapies applied to the cerebral cortex.

A finite element model of a six-year-old child for simulating pedestrian accidents

Child pedestrian protection deserves more attention in vehicle safety design since they are the most vulnerable road users who face the highest mortality rate. Pediatric Finite Element (FE) models could be used to simulate and understand the pedestrian injury mechanisms during crashes in order to mitigate them. Thus, the objective of the study was to develop a computationally efficient (simplified) six-year-old (6YO-PS) pedestrian FE model and validate it based on the latest published pediatric data. The 6YO-PS FE model was developed by morphing the existing GHBMC adult pedestrian model. Retrospective scan data were used to locally adjust the geometry as needed for accuracy. Component test simulations focused only the lower extremities and pelvis, which are the first body regions impacted during pedestrian accidents. Three-point bending test simulations were performed on the femur and tibia with adult material properties and then updated using child material properties. Pelvis impact and knee bending tests were also simulated. Finally, a series of pediatric Car-to-Pedestrian Collision (CPC) were simulated with pre-impact velocities ranging from 20km/h up to 60km/h. The bone models assigned pediatric material properties showed lower stiffness and a good match in terms of fracture force to the test data (less than 6% error). The pelvis impact force predicted by the child model showed a similar trend with test data. The whole pedestrian model was stable during CPC simulations and predicted common pedestrian injuries. Overall, the 6YO-PS FE model developed in this study showed good biofidelity at component level (lower extremity and pelvis) and stability in CPC simulations. While more validations would improve it, the current model could be used to investigate the lower limb injury mechanisms and in the prediction of the impact parameters as specified in regulatory testing protocols.

Development and Full Body Validation of a 5th Percentile Female Finite Element Model

To mitigate the societal impact of vehicle crash, researchers are using a variety of tools, including finite element models (FEMs). As part of the Global Human Body Models Consortium (GHBMC) project, comprehensive medical image and anthropometrical data of the 5th percentile female (F05) were acquired for the explicit purpose of FEM development. The F05-O (occupant) FEM model consists of 981 parts, 2.6 million elements, 1.4 million nodes, and has a mass of 51.1 kg. The model was compared to experimental data in 10 validation cases ranging from localized rigid hub impacts to full body sled cases. In order to make direct comparisons to experimental data, which represent the mass of an average male, the model was compared to experimental corridors using two methods: 1) post-hoc scaling the outputs from the baseline F05-O model and 2) geometrically morphing the model to the body habitus of the average male to allow direct comparisons. This second step required running the morphed full body model in all 10 simulations for a total of 20 full body simulations presented. Overall, geometrically morphing the model was found to more closely match the target data with an average ISO score for the rigid impacts of 0.76 compared to 0.67 for the scaled responses. Based on these data, the morphed model was then used for model validation in the vehicle sled cases. Overall, the morphed model attained an average weighted score of 0.69 for the two sled impacts. Hard tissue injuries were also assessed and the baseline F05-O model was found to predict a greater occurrence of pelvic fractures compared to the GHBMC average male model, but predicted fewer rib fractures.

An Objective Evaluation of Mass Scaling Techniques Utilizing Computational Human Body Finite Element Models

Biofidelity response corridors developed from post-mortem human subjects are commonly used in the design and validation of anthropomorphic test devices and computational human body models (HBMs). Typically, corridors are derived from a diverse pool of biomechanical data and later normalized to a target body habitus. The objective of this study was to use morphed computational HBMs to compare the ability of various scaling techniques to scale response data from a reference to a target anthropometry. HBMs are ideally suited for this type of study since they uphold the assumptions of equal density and modulus that are implicit in scaling method development. In total, six scaling procedures were evaluated, four from the literature (equal-stress equal-velocity, ESEV, and three variations of impulse momentum) and two which are introduced in the paper (ESEV using a ratio of effective masses, ESEV-EffMass, and a kinetic energy approach). In total, 24 simulations were performed, representing both pendulum and full body impacts for three representative HBMs. These simulations were quantitatively compared using the International Organization for Standardization (ISO) ISO-TS18571 standard. Based on these results, ESEV-EffMass achieved the highest overall similarity score (indicating that it is most proficient at scaling a reference response to a target). Additionally, ESEV was found to perform poorly for two degree-of-freedom (DOF) systems. However, the results also indicated that no single technique was clearly the most appropriate for all scenarios.

Finite element comparison of human and Hybrid III responses in a frontal impact

The improvement of finite element (FE) Human Body Models (HBMs) has made them valuable tools for investigating restraint interactions compared to anthropomorphic test devices (ATDs). The objective of this study was to evaluate the effect of various combinations of safety restraint systems on the sensitivity of thoracic injury criteria using matched ATD and Human Body Model (HBM) simulations at two crash severities. A total of seven (7) variables were investigated: 3-point belt with two (2) load limits, frontal airbag, knee bolster airbag, a buckle pretensioner, and two (2) delta-v's - 40kph and 50kph. Twenty four (24) simulations were conducted for the Hybrid III ATD FE model and repeated with a validated HBM for 48 total simulations. Metrics tested in these conditions included sternum deflection, chest acceleration, chest excursion, Viscous Criteria (V*C) criteria, pelvis acceleration, pelvis excursion, and femur forces. Additionally, chest band deflection and rib strain distribution were measured in the HBM for additional restraint condition discrimination. The addition of a frontal airbag had the largest effect on the occupant chest metrics with an increase in chest compression and acceleration but a decrease in excursion. While the THUMS and Hybrid III occupants demonstrated the same trend in the chest compression measurements, there were conflicting results in the V*C, acceleration, and displacement metrics. Similarly, the knee bolster airbag had the largest effect on the pelvis with a decrease in acceleration and excursion. With a knee bolster airbag the simulated occupants gave conflicting results, the THUMS had a decrease in femur force and the ATD had an increase. Preferential use of dummies or HBM's is not debated; however, this study highlights the ability of HBM metrics to capture additional chest response metrics.

Head and neck response of a finite element anthropomorphic test device and human body model during a simulated rotary-wing aircraft impact

A finite element (FE) simulation environment has been developed to investigate aviator head and neck response during a simulated rotary-wing aircraft impact using both an FE anthropomorphic test device (ATD) and an FE human body model. The head and neck response of the ATD simulation was successfully validated against an experimental sled test. The majority of the head and neck transducer time histories received a CORrelation and analysis (CORA) rating of 0.7 or higher, indicating good overall correlation. The human body model simulation produced a more biofidelic head and neck response than the ATD experimental test and simulation, including change in neck curvature. While only the upper and lower neck loading can be measured in the ATD, the shear force, axial force, and bending moment were reported for each level of the cervical spine in the human body model using a novel technique involving cross sections. This loading distribution provides further insight into the biomechanical response of the neck during a rotary-wing aircraft impact.