External Landmark, Body Surface, and Volume Data of a Mid-Sized Male in Seated and Standing Postures

A Review of Image-Based Simulation Applications in High-Value Manufacturing

Image-Based Simulation (IBSim) is the process by which a digital representation of a real geometry is generated from image data for the purpose of performing a simulation with greater accuracy than with idealised Computer Aided Design (CAD) based simulations. Whilst IBSim originates in the biomedical field, the wider adoption of imaging for non-destructive testing and evaluation (NDT/NDE) within the High-Value Manufacturing (HVM) sector has allowed wider use of IBSim in recent years. IBSim is invaluable in scenarios where there exists a non-negligible variation between the 'as designed' and 'as manufactured' state of parts. It has also been used for characterisation of geometries too complex to accurately draw with CAD. IBSim simulations are unique to the geometry being imaged, therefore it is possible to perform part-specific virtual testing within batches of manufactured parts. This novel review presents the applications of IBSim within HVM, whereby HVM is the value provided by a manufactured part (or conversely the potential cost should the part fail) rather than the actual cost of manufacturing the part itself. Examples include fibre and aggregate composite materials, additive manufacturing, foams, and interface bonding such as welding. This review is divided into the following sections: Material Characterisation; Characterisation of Manufacturing Techniques; Impact of Deviations from Idealised Design Geometry on Product Design and Performance; Customisation and Personalisation of Products; IBSim in Biomimicry. Finally, conclusions are drawn, and observations made on future trends based on the current state of the literature.

Body mass index influence on lap belt position and abdominal injury in frontal motor vehicle crashes

OBJECTIVE

As obesity rates climb, it is important to study its effects on motor vehicle safety due to differences in restraint interaction and biomechanics. Previous studies have shown that an abdominal seatbelt sign (referred hereafter as seatbelt sign) sustained from motor vehicle crashes (MVCs) is associated with abdominal trauma when located above the anterior superior iliac spine (ASIS). This study investigates whether placement of the lap belt causing a seatbelt sign is associated with abdominal organ injury in occupants with increased body mass index (BMI). We hypothesized that higher BMI would be associated with a higher incidence of superior placement of the lap belt to the ASIS level, and a higher incidence of abdominal organ injury.

METHODS

A retrospective data analysis was performed using 230 cases that met inclusion criteria (belted occupant in a frontal collision that sustained at least one abdominal injury) from the Crash Injury Research and Engineering Network (CIREN) database. Computed tomography (CT) scans were rendered to visualize fat stranding to determine the presence of a seatbelt sign. 146 positive seatbelt signs were visualized. ASIS level was measured by adjusting the transverse slice of the CT to the visualized ASIS level, which was used to determine seatbelt sign location as superior, on, or inferior to the ASIS.

RESULTS

Obese occupants had a significantly higher incidence of superior belt placement (52%) vs on-ASIS placement (24%) compared to their normal (27% vs 67%) BMI counterparts (p < 0.001). Notable trends included obese occupants with superior placement having less abdominal organ injury incidence than those with on-ASIS belt placement (42% superior placement vs 55% on-ASIS). In non-obese occupants, there was a higher incidence of abdominal organ injury with superior lap belt placement compared to on-ASIS placement counterparts (Normal BMI: 62% vs 41%, Overweight: 57% vs 43%).

CONCLUSIONS

In CIREN occupants with abdominal injury, those with obesity are more prone to positioning the lap belt superior to the ASIS, though the impact on abdominal injury incidence remains a key point for continued exploration into how occupant BMI affects crash safety and belt design.

The relationship of body mass index, belt placement, and abdominopelvic injuries in motor vehicle crashes: A Crash Injury Research and Engineering Network (CIREN) study

OBJECTIVE

Obesity has important implications for motor vehicle safety due to altered crash injury responses from increased mass and improper seatbelt placement. Abdominal seatbelt signs (ASBS) above the anterior superior iliac spine (ASIS) in motor vehicle crashes (MVCs) often correlate with abdominopelvic trauma. We investigated the relationship of body mass index (BMI), lap belt placement, and the incidence of abdominopelvic injury using computed tomography (CT) evaluation for subcutaneous ASBS mark and its location relative to the ASIS.

METHODS

A retrospective analysis of 235 Crash Injury Research and Engineering Network (CIREN) cases and their associated abdominal injuries was conducted. CT Scans were analyzed to visualize fat stranding. 150 positive ASBS were found and their ASBS mark location was classified as superior, on, or inferior to the ASIS.

RESULTS

Obese occupants had a higher incidence rate of belt placement superior to the ASIS, and occupants with normal BMI had a higher incidence of proper belt placement (p < 0.05). Trends of interest developed, notably that non-obese occupants with superior belt placement had increased incidence of internal abdominopelvic organ injury compared to those with proper belt placement (Normal BMI: 53.3% superior vs 39.4% On-ASIS, Overweight: 47.8% superior vs 34.7% On-ASIS).

CONCLUSIONS

Utilizing CT scans to confirm ASBS and lap belt placement relative to the ASIS, superior belt placement above the ASIS was associated with elevated BMI and a trend of increasing incidence of internal abdominopelvic organ injury.

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.

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.

Objective Evaluation of Whole Body Kinematics in a Simulated, Restrained Frontal Impact

The use of human body models as an additional data point in the evaluation of human-machine interaction requires quantitative validation. In this study a validation of the Global Human Body Models Consortium (GHBMC) average male occupant model (M50-O v. 4.5) in a restrained frontal sled test environment is presented. For vehicle passengers, frontal crash remains the most common mode, and the most common source of fatalities. A total of 55-time history traces of reaction loads and kinematics from the model were evaluated against corresponding PMHS data (n = 5). Further, the model's sensitivity to the belt path was studied by replicating two documented PMHS cases with prominent lateral and medial belt paths respectively. Results were quantitatively evaluated using open source CORA software. A tradeoff was observed; better correlation scores were achieved on gross measures (e.g. reaction loads), whereas better corridor scores were achieved on localized measures (rib deflections), indicating that subject specificity may dominate the comparison at localized anatomical regions. On an overall basis, the CORA scores were 0.68, 0.66 and 0.60 for force, body kinematics and chest wall kinematics. Belt force responses received the highest grouped CORA score of 0.85. Head and sternum kinematics earning a 0.8 and 0.7 score respectively. The model demonstrated high sensitivity to belt path, resulting in a 20-point increase in CORA score when the belt was routed closer to analogous location of data collection. The human model demonstrated overall reasonable biofidelity and sensitivity to countermeasures in frontal crash kinematics.

Quantification of Skeletal and Soft Tissue Contributions to Thoracic Response in a Dynamic Frontal Loading Scenario

Thoracic injuries continue to be a major health concern in motor vehicle crashes. Previous thoracic research has focused on 50^th percentile males and utilized scaling techniques to apply results to different demographics. Individual rib testing offers the advantage of capturing demographic differences; however, understanding of rib properties in the context of the intact thorax is lacking. Therefore, the objective of this study was to obtain the data necessary to develop a transfer function between individual rib and thoracic response. A series of non-injurious frontal impacts were conducted on six PMHS, creating a loading environment commensurate to previously published individual rib testing. Each PMHS was tested in four tissue states: intact, intact with upper limbs removed, denuded, and eviscerated. Following eviscerated thoracic testing, eight individual mid-level ribs from each PMHS were removed and loaded to failure. A simplified model in which ribs of each thorax are treated as parallel springs was utilized to evaluate the ability of individual rib response data to predict each subject's eviscerated thoracic response. On average across subjects, denuded thoraces retained 89% and eviscerated thoraces retained 46% of intact force. Similarly, denuded thoraces retained 70% and eviscerated thoraces retained 30% of intact stiffness. The rib model did not adequately predict eviscerated thoracic response but provided a better understanding of the influence of connective tissue on a rib's behavior with-in the thorax. Results of this study could be used in conjunction with the database of individual rib test results to improve thoracic response targets and help assess biofidelity of current anthropomorphic test devices.

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.

New Reference PMHS Tests to Assess Whole-Body Pedestrian Impact Using a Simplified Generic Vehicle Front-End

This study aims to provide a set of reference post-mortem human subject tests which can be used, with easily reproducible test conditions, for developing and/or validating pedestrian dummies and computational human body models against a road vehicle. An adjustable generic buck was first developed to represent vehicle front-ends. It was composed of four components: two steel cylindrical tubes screwed on rigid supports in V-form represent the bumper and spoiler respectively, a quarter of a steel cylindrical tube represents the bonnet leading edge, and a steel plate represents the bonnet. These components were positioned differently to represent three types of vehicle profile: a sedan, a SUV and a van. Eleven post-mortem human subjects were then impacted laterally in a mid-gait stance by the bucks at 40 km/h: three tests with the sedan, five with the SUV, and three with the van. Kinematics of the subjects were recorded via high speed videos, impact forces between the subjects and the bucks were measured via load cells behind each tube, femur and tibia deformation and fractures were monitored via gauges on these bones. Based on these tests, biofidelity corridors were established in terms of: 1) displacement time history and trajectory of the head, shoulder, T1, T4, T12, sacrum, knee and ankle, 2) impact forces between the subjects and the buck. Injury outcome was established for each PMHS via autopsy. Simplicity of its geometry and use of standard steel tubes and plates for the buck will make it easy to perform future, new post-mortem human subject tests in the same conditions, or to assess dummies or computational human body models using these reference tests.

A Finite Element Model of a Midsize Male for Simulating Pedestrian Accidents

Pedestrians represent one of the most vulnerable road users and comprise nearly 22% the road crash-related fatalities in the world. Therefore, protection of pedestrians in car-to-pedestrian collisions (CPC) has recently generated increased attention with regulations involving three subsystem tests. The development of a finite element (FE) pedestrian model could provide a complementary component that characterizes the whole-body response of vehicle–pedestrian interactions and assesses the pedestrian injuries. The main goal of this study was to develop and to validate a simplified full body FE model corresponding to a 50th male pedestrian in standing posture (M50-PS). The FE model mesh and defined material properties are based on a 50th percentile male occupant model. The lower limb-pelvis and lumbar spine regions of the human model were validated against the postmortem human surrogate (PMHS) test data recorded in four-point lateral knee bending tests, pelvic\abdomen\shoulder\thoracic impact tests, and lumbar spine bending tests. Then, a pedestrian-to-vehicle impact simulation was performed using the whole pedestrian model, and the results were compared to corresponding PMHS tests. Overall, the simulation results showed that lower leg response is mostly within the boundaries of PMHS corridors. In addition, the model shows the capability to predict the most common lower extremity injuries observed in pedestrian accidents. Generally, the validated pedestrian model may be used by safety researchers in the design of front ends of new vehicles in order to increase pedestrian protection.

The Role of Fluid Dynamics in Distributing Ankle Stresses in Anatomic and Injured States

BACKGROUND

In 1976, Ramsey and Hamilton published a landmark cadaveric study demonstrating a dramatic 42% decrease in tibiotalar contact area with only 1 mm of lateral talar shift. An increase in maximum principal stress of at least 72% is predicted based on these findings though the delayed development of arthritis in minimally misaligned ankles does not appear to be commensurate with the results found in dry cadaveric models. We hypothesized that synovial fluid could be a previously unrecognized factor that contributes significantly to stress distribution in the tibiotalar joint in anatomic and injured states.

METHODS

As it is not possible to directly measure contact stresses with and without fluid in a cadaveric model, finite element analysis (FEA) was employed for this study. FEA is a modeling technique used to calculate stresses in complex geometric structures by dividing them into small, simple components called elements. Four test configurations were investigated using a finite element model (FEM): baseline ankle alignment, 1 mm laterally translated talus and fibula, and the previous 2 bone orientations with fluid added. The FEM selected for this study was the Global Human Body Models Consortium-owned GHBMC model, M50 version 4.2, a model of an average-sized male (distributed by Elemance, LLC, Winston-Salem, NC). The ankle was loaded at the proximal tibia with a distributed load equal to the GHBMC body weight, and the maximum principal stress was computed.

RESULTS

All numerical simulations were stable and completed with no errors. In the baseline anatomic configuration, the addition of fluid between the tibia, fibula, and talus reduced the maximum principal stress computed in the distal tibia at maximum load from 31.3 N/mm² to 11.5 N/mm². Following 1 mm lateral translation of the talus and fibula, there was a modest 30% increase in the maximum stress in fluid cases. Qualitatively, translation created less high stress locations on the tibial plafond when fluid was incorporated into the model.

CONCLUSIONS

The findings in this study demonstrate a meaningful role for synovial fluid in distributing stresses within the ankle that has not been considered in historical dry cadaveric studies. The increase in maximum stress predicted by simulation of an ankle with fluid was less than half that projected by cadaveric data, indicating a protective effect of fluid in the injured state. The trends demonstrated by these simulations suggest that bony alignment and fluid in the ankle joint change loading patterns on the tibia and should be accounted for in future experiments.

CLINICAL RELEVANCE

Synovial fluid may play a protective role in ankle injuries, thus delaying the onset of arthritis. Reactive joint effusions may also function to additionally redistribute stresses with higher volumes of viscous fluid.

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.

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.

Thoracoabdominal organ volumes for small women

OBJECTIVE

Thoracoabdominal injuries commonly occur as a result of motor vehicle crashes. In order to design occupant protection systems that reduce risk of injury, researchers are using a variety of tools, including computational human body models. Though research has been conducted to provide morphological and volumetric data for the thoracoabdominal cavity of the average male, there is currently an interest in developing models for a wider range of occupants. One particular cohort of interest is the small female by stature and weight because of their use in restraint system development. Geometric data on thoracoabdominal organs are needed to construct accurate representations of female occupants. This study aimed to gather information on organ volumes from clinical medical imaging studies of small females.

METHODS

Anonymized clinical computed tomography (CT) and magnetic resonance images were used to segment organs relevant to crash-induced injuries: namely, the liver, spleen, left kidney, right kidney, pancreas, gallbladder, lungs, and heart. Segmentations were conducted using semi-automatic techniques. Additionally, diametric measurements of the vena cava, aorta, trachea, and colon were obtained from the medical images at discrete locations using linear measurement tools.

RESULTS

A total of 14 adult scans were selected with stature and weight ranges of 145.0 to 162.6 cm and 43.7 to 65.5 kg, respectively. The following are the average thoracoabdominal organ volumes: liver (1,224.5 ± 220.7 mL), spleen (151.6 ± 42.1 mL), left kidney (123.7 ± 20.1 mL), right kidney (115.4 ± 20.9 mL), heart (417.8 ± 36.6 mL), pancreas (54.1 ± 11.8 mL), and gallbladder (20.6 ± 13.4 mL). The average diameters were 19.7 ± 3.2 mm and 17.7 ± 5.1 mm for the vena cava and aorta, respectively. The colon had an average diameter of 37.9 ± 7.1 mm.

CONCLUSION

Data characterizing the small female are important to validate the geometries used in computational models, including models derived from scaling techniques and those developed using subject-specific medical imaging. The goal of this study was to use a sample of subjects anthropometrically representative of small females to evaluate the average volume for organs commonly injured in motor vehicle crashes. Based on these data, the right and left lungs were strongly correlated with stature and the heart was strongly correlated with weight. Ultimately, these measurements will be useful for the validation of computational models of the small female.

Application of Radial Basis Function Methods in the Development of a 95th Percentile Male Seated FEA Model

Human body finite element models (FEMs) are a valuable tool in the study of injury biomechanics. However, the traditional model development process can be time-consuming. Scaling and morphing an existing FEM is an attractive alternative for generating morphologically distinct models for further study. The objective of this work is to use a radial basis function to morph the Global Human Body Models Consortium (GHBMC) average male model (M50) to the body habitus of a 95th percentile male (M95) and to perform validation tests on the resulting model. The GHBMC M50 model (v. 4.3) was created using anthropometric and imaging data from a living subject representing a 50th percentile male. A similar dataset was collected from a 95th percentile male (22,067 total images) and was used in the morphing process. Homologous landmarks on the reference (M50) and target (M95) geometries, with the existing FE node locations (M50 model), were inputs to the morphing algorithm. The radial basis function was applied to morph the FE model. The model represented a mass of 103.3 kg and contained 2.2 million elements with 1.3 million nodes. Simulations of the M95 in seven loading scenarios were presented ranging from a chest pendulum impact to a lateral sled test. The morphed model matched anthropometric data to within a rootmean square difference of 4.4% while maintaining element quality commensurate to the M50 model and matching other anatomical ranges and targets. The simulation validation data matched experimental data well in most cases.

Investigation of the Mass Distribution of a Detailed Seated Male Finite Element Model

Accurate mass distribution in computational human body models is essential for proper kinematic and kinetic simulations. The purpose of this study was to investigate the mass distribution of a 50th percentile male (M50) full body finite element model (FEM) in the seated position. The FEM was partitioned into 10 segments, using segment planes constructed from bony landmarks per the methods described in previous research studies. Body segment masses and centers of gravity (CGs) of the FEM were compared with values found from these studies, which unlike the present work assumed homogeneous body density. Segment masses compared well to literature while CGs showed an average deviation of 6.0% to 7.0% when normalized by regional characteristic lengths. The discrete mass distribution of the FEM appears to affect the mass and CGs of some segments, particularly those with low-density soft tissues. The locations of the segment CGs are provided in local coordinate systems, thus facilitating comparison with other full body FEMs and human surrogates. The model provides insights into the effects of inhomogeneous mass on the location of body segment CGs.

Validation of simulated chestband data in frontal and lateral loading using a human body finite element model

OBJECTIVE

Finite element (FE) computer models are an emerging tool to examine the thoracic response of the human body in the simulated environment. In this study, a recently developed human body model, the Global Human Body Models Consortium (GHBMC) mid-sized male, was used to examine chestband contour deformations in a frontal and lateral impact. The objective of this study was 2-fold. First, a methodology for extracting and analyzing virtual chestband data from a full-body FE model is presented. Then, this method is applied to virtual chestband data from 2 simulations to evaluate the model's performance against experimental data.

METHODS

One frontal and one lateral impact case were simulated using the FE model, which was preprogrammed with upper, middle, and lower chestbands. Maximum compression was determined using established techniques. Furthermore, a quadrant-based analysis technique for the results was introduced that enabled regional comparisons between the model and the experimental data in the anterior, posterior, right, and left sections of the chestband.

RESULTS

For the frontal case at 13.3 m/s, the model predicted a peak compression of 13.6 and 12.9 percent for the upper and middle chestbands. For the lateral case at 6.7 m/s, the model predicted peak compression of the upper, middle, and lower chestbands of 27.9, 26.0, and 20.4 percent. Regional analysis showed average differences at maximum deformation between the model and experiments ranging from 0.9 percent (posterior) to 6.3 percent (anterior) in the frontal case and 2.3 percent (posterior) to 10.8 percent (anterior) in the lateral case. The greatest difference between model and experimental findings was found in the anterior quadrant.

CONCLUSIONS

Though this work was focused on techniques to extract and analyze chestband data from FE models, the comparative results provide further validation of the model used in this study. The results suggest the importance of evaluating comparisons between virtual and experimental chestband data on a regional basis. These data also provide the potential to correlate chestband deformations to the loading of underlying thoraco-abdominal structures. Supplemental materials are available for this article. Go to the publisher's online edition of Traffic Injury Prevention to view the supplemental file.

Development of brain injury criteria (BrIC)

Rotational motion of the head as a mechanism for brain injury was proposed back in the 1940s. Since then a multitude of research studies by various institutions were conducted to confirm/reject this hypothesis. Most of the studies were conducted on animals and concluded that rotational kinematics experienced by the animal's head may cause axonal deformations large enough to induce their functional deficit. Other studies utilized physical and mathematical models of human and animal heads to derive brain injury criteria based on deformation/pressure histories computed from their models. This study differs from the previous research in the following ways: first, it uses two different detailed mathematical models of human head (SIMon and GHBMC), each validated against various human brain response datasets; then establishes physical (strain and stress based) injury criteria for various types of brain injury based on scaled animal injury data; and finally, uses Anthropomorphic Test Devices (ATDs) (Hybrid III 50th Male, Hybrid III 5th Female, THOR 50th Male, ES-2re, SID-IIs, WorldSID 50th Male, and WorldSID 5th Female) test data (NCAP, pendulum, and frontal offset tests) to establish a kinematically based brain injury criterion (BrIC) for all ATDs. Similar procedures were applied to college football data where thousands of head impacts were recorded using a six degrees of freedom (6 DOF) instrumented helmet system. Since animal injury data used in derivation of BrIC were predominantly for diffuse axonal injury (DAI) type, which is currently an AIS 4+ injury, cumulative strain damage measure (CSDM) and maximum principal strain (MPS) were used to derive risk curves for AIS 4+ anatomic brain injuries. The AIS 1+, 2+, 3+, and 5+ risk curves for CSDM and MPS were then computed using the ratios between corresponding risk curves for head injury criterion (HIC) at a 50% risk. The risk curves for BrIC were then obtained from CSDM and MPS risk curves using the linear relationship between CSDM - BrIC and MPS - BrIC respectively. AIS 3+, 4+ and 5+ field risk of anatomic brain injuries was also estimated using the National Automotive Sampling System - Crashworthiness Data System (NASS-CDS) database for crash conditions similar to the frontal NCAP and side impact conditions that the ATDs were tested in. This was done to assess the risk curve ratios derived from HIC risk curves. The results of the study indicated that: (1) the two available human head models - SIMon and GHBMC - were found to be highly correlated when CSDMs and max principal strains were compared; (2) BrIC correlates best to both - CSDM and MPS, and rotational velocity (not rotational acceleration) is the mechanism for brain injuries; and (3) the critical values for angular velocity are directionally dependent, and are independent of the ATD used for measuring them. The newly developed brain injury criterion is a complement to the existing HIC, which is based on translational accelerations. Together, the two criteria may be able to capture most brain injuries and skull fractures occurring in automotive or any other impact environment. One of the main limitations for any brain injury criterion, including BrIC, is the lack of human injury data to validate the criteria against, although some approximation for AIS 2+ injury is given based on the angular velocities calculated at 50% probability of concussion in college football players instrumented with 5 DOF helmet system. Despite the limitations, a new kinematic rotational brain injury criterion - BrIC - may offer a way to capture brain injuries in situations when using translational accelerations based HIC alone may not be sufficient.

Comparison of organ location, morphology, and rib coverage of a midsized male in the supine and seated positions

The location and morphology of abdominal organs due to postural changes have implications in the prediction of trauma via computational models. The purpose of this study is to use data from a multimodality image set to devise a method for examining changes in organ location, morphology, and rib coverage from the supine to seated postures. Medical images of a male volunteer (78.6 ± 0.77 kg, 175 cm) in three modalities (Computed Tomography, Magnetic Resonance Imaging (MRI), and Upright MRI) were used. Through image segmentation and registration, an analysis between organs in each posture was conducted. For the organs analyzed (liver, spleen, and kidneys), location was found to vary between postures. Increases in rib coverage from the supine to seated posture were observed for the liver, with a 9.6% increase in a lateral projection and a 4.6% increase in a frontal projection. Rib coverage area was found to increase 11.7% for the spleen. Morphological changes in the organs were also observed. The liver expanded 7.8% cranially and compressed 3.4% and 5.2% in the anterior-posterior and medial-lateral directions, respectively. Similar trends were observed in the spleen and kidneys. These findings indicate that the posture of the subject has implications in computational human body model development.

An evaluation of objective rating methods for full-body finite element model comparison to PMHS tests

OBJECTIVE

Objective evaluation methods of time history signals are used to quantify how well simulated human body responses match experimental data. As the use of simulations grows in the field of biomechanics, there is a need to establish standard approaches for comparisons. There are 2 aims of this study. The first is to apply 3 objective evaluation methods found in the literature to a set of data from a human body finite element model. The second is to compare the results of each method, examining how they are correlated to each other and the relative strengths and weaknesses of the algorithms.

METHODS

In this study, the methods proposed by Sprague and Geers (magnitude and phase error, SGM and SGP), Rhule et al. (cumulative standard deviation, CSD), and Gehre et al. (CORrelation and Analysis, or CORA, size, phase, shape, corridor) were compared. A 40 kph frontal sled test presented by Shaw et al. was simulated using the Global Human Body Models Consortium midsized male full-body finite element model (v. 3.5). Mean and standard deviation experimental data (n = 5) from Shaw et al. were used as the benchmark. Simulated data were output from the model at the appropriate anatomical locations for kinematic comparison. Force data were output at the seat belts, seat pan, knee, and foot restraints.

RESULTS

Objective comparisons from 53 time history data channels were compared to the experimental results. To compare the different methods, all objective comparison metrics were cross-plotted and linear regressions were calculated. The following ratings were found to be statistically significantly correlated (P < .01): SGM and CORrelation and Analysis (CORA) size, R (2) = 0.73; SGP and CORA shape, R (2) = 0.82; and CSD and CORA's corridor factor, R (2) = 0.59. Relative strengths of the correlated ratings were then investigated. For example, though correlated to CORA size, SGM carries a sign to indicate whether the simulated response is greater than or less than the benchmark signal. A further analysis of the advantages and drawbacks of each method is discussed.

CONCLUSIONS

The results demonstrate that a single metric is insufficient to provide a complete assessment of how well the simulated results match the experiments. The CORA method provided the most comprehensive evaluation of the signal. Regardless of the method selected, one primary recommendation of this work is that for any comparison, the results should be reported to provide separate assessments of a signal's match to experimental variance, magnitude, phase, and shape. Future work planned includes implementing any forthcoming International Organization for Standardization standards for objective evaluations. Supplemental materials are available for this article. Go to the publisher's online edition of Traffic Injury Prevention to view the supplemental file.