Lateral Impact Validation of a Geometrically Accurate Full Body Finite Element Model for Blunt Injury Prediction

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.

Development and Validation of a Whole Human Body Finite Element Model with Detailed Lumbar Spine

OBJECTIVE

Investigations showed that low back pain of occupational drivers might be closely related to the whole-body vibration. Restricted by ethical concerns, the finite element method had become a viable alternative to invasive human experiments. Many mechanical behaviors of the human spine inside of the human body were unclear; therefore, a human whole-body finite element model might be required to better understand the lumbar behavior under whole-body vibration.

METHODS

In this study, a human whole-body finite element model with a detailed lumbar spine segment was developed. Several validations were performed to ensure the correctness of this model.

RESULTS

The results of anthropometry and geometry validation, static validation, and dynamic validation were presented in this study. The validation results showed that the whole human body model was reasonable and valid by comparing with published data.

CONCLUSIONS

The model developed in this study could reflect the biomechanical response of the human lumbar spine under vibration and could be used in further vibration analysis and offer proposals for protecting human body under whole-body vibration environment.

Longitudinal Tibia Stress Fracture Risk During High-Volume Training: a Multi-Scale Modeling Pipeline Incorporating Bone Remodeling

Tibia stress fractures are prevalent during high-intensity training, yet a mechanistic model linking longitudinal training intensity, bone health, and long-term injury risk has yet to be demonstrated. The objective of this study was to develop and validate a multi-scale model of gross and tissue level loading on the tibia including bone remodeling on a timescale of week. Peak tensile tibial strain (3517 µstrain) during 4 m/s running was below injury thresholds, and the peak anteromedial tibial strain (1248 µstrain) was 0.17 standard deviations away from the mean of reported literature values. An initial study isolated the effects of cortical density and stiffness on tibial strain during a simulated eight week training period. Tibial strains and cortical microcracking correlated to initial cortical modulus, with all simulations presenting peak anteromedial tensile strains (1047-1600 µstrain) near day 11. Average cortical densities decreased by 7-8 percent of their nominal value by day 11, but the overall density change was <2% by the end of the simulated training period, in line with reported results. This study demonstrates the benefits of multi-scale models for investigating stress fracture risk and indicates that peak tibial strain, and thus injury risk, may increase early in a high intensity training program. Future studies could optimize training volume and recovery time to reduce injury risk during the most vulnerable training periods.

Effects of Loading Conditions on the Pelvic Fracture Biomechanism and Discrimination of Forensic Injury Manners of Impact and Run-Over Using the Finite Element Pelvic Model

This study aimed to systematically simulate the responses of pelvic fracture under impact and run-over to clarify the effects of boundary and loading conditions on the pelvic fracture mechanism and provide complementary quantitative evidence for forensic practice. Based on the THUMS finite element model, we have validated the simulation performance of the model by a real postmortem human pelvis side impact experiment. A total of 54 simulations with two injury manners (impact and run-over), seven loading directions (0°, 30°, 60°, 90°, 270°, 300°, 330°), and six loading velocities (10, 20, 30, 40, 50, and 60 km/h) were conducted. Criteria of effective strain, Von-Mises stress, contact force, and self-designed normalized eccentricity were used to evaluate the biomechanism of pelvic fracture. Based on our simulation results, it’s challenging to distinguish impact from run-over only rely on certain characteristic fractures. Loads on the front and back were less likely to cause pelvic fractures. In the 30°, 60°, 300° load directions, the overall deformation caused a “diagonal” pelvic fracture. The higher is the velocity (kinetic energy), the more severe is the pelvic fracture. The contact force will predict the risk of fracture. In addition, our self-designed eccentricity will distinguish the injury manner of impact and run-over under the 90° loads. The “biomechanical fingerprints” based on logistic regression of all biomechanical variables have an AUC of 0.941 in discriminating the injury manners. Our study may provide simulation evidence and new methods for the forensic community to improve the forensic identification ability of injury manners.

Subject-specific rib finite element models with material data derived from coupon tests under bending loading

Rib fractures are common thoracic injuries in motor vehicle crashes. Several human finite element (FE) human models have been created to numerically assess thoracic injury risks. However, the accurate prediction of rib biomechanical response has shown to be challenging due to human variation and modeling approaches. The main objective of this study was to better understand the role of modeling approaches on the biomechanical response of human ribs in anterior-posterior bending. Since the development of subject specific rib models is a time-consuming process, the second objective of this study was to develop an accurate morphing approach to quickly generate high quality subject specific rib meshes. The exterior geometries and cortical-trabecular boundaries of five human 6th-level ribs were extracted from CT-images. One rib mesh was developed in a parametric fashion and the other four ribs were developed with an in-house morphing algorithm. The morphing algorithm automatically defined landmarks on both the periosteal and endosteal boundaries of the cortical layer, which were used to morph the template nodes to target geometries. Three different cortical bone material models were defined based on the stress-strain data obtained from subject-specific tensile coupon tests for each rib. Full rib anterior-posterior bending tests were simulated based on data recorded in testing. The results showed similar trends to test data with some sensitivity relative to the material modeling approach. Additionally, the FE models were substantially more resistant to failure, highlighting the need for better techniques to model rib fracture. Overall, the results of this work can be used to improve the biofidelity of human rib finite element models.

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.

Biomechanical properties of abdominal organs under tension with special reference to increasing strain rate

Currently, abdominal finite element models overlook the organs such as gallbladder, bladder, and intestines, which instead are modeled as a simple bag that is not included in the analysis. Further characterization of the material properties is required for researchers to include these organs into models. This study characterized the mechanical properties of human and porcine gallbladder, bladder, and intestines using uniaxial tension loading from the rates of 25%/s to 500%/s. Small differences were observed between human and porcine gallbladder elastic modulus, failure stress, and failure strain. Strain rate was determined to be a significant factor for predicting porcine gallbladder elastic modulus and failure stress which were found to be 9.03 MPa and 1.83 MPa at 500%/s. Human bladder was observed to be slightly stiffer with a slightly lower failure stress than porcine specimens. Both hosts, however, demonstrated a strain rate dependency with the elastic modulus and failure stress increasing as the rate increased with the highest elastic modulus (2.16 MPa) and failure stress (0.65 MPa) occurring at 500%/s. Both human and porcine intestines were observed to be affected by the strain rate. Failure stress was found to be 1.6 MPa and 1.42 MPa at 500%/s for the human and porcine intestines respectively. For all properties found to be strain rate dependent, a numerical model was created to quantify the impact. These results will enable researchers to create more detailed finite element models that include the gallbladder, bladder, and intestines.

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.

Comparing FE human body model rib geometry to population data

Finite element human body models (HBMs) are used to assess injury risk in a variety of impact scenarios. The ribs are a key structural component within the chest, so their accuracy within HBMs is vitally important for modeling human biomechanics. We assessed the geometric correspondence between the ribs defined within five widely used HBMs and measures drawn from population-wide studies of rib geometry, focusing on (1) rib global shape, (2) rib cross-sectional size and shape, and (3) rib cortical bone thickness. A parametric global shape model fitted to all HBM ribs was compared to expected rib parameters calculated for each HBM's subject demographic using population reference data. The GHBMC M50 and THUMS M50 male HBMs showed 24% and 50% of their fitted rib shape parameters (6 parameters per each 12 ribs) falling outside 1SD from population expected values, respectively. For female models the GHBMC F05, THUMS F05, and VIVA F50 models had 21%, 26%, and 19% of their rib shape parameters falling outside 1SD, respectively. Cross-sectional areas and inertial moments obtained along the HBM ribs were compared to average ± 1SD corridors for male and female ribs drawn from reference population data. The GHBMC M50, THUMS M50, and VIVA F50 model ribs were all larger in overall cross-sectional area than their targeted average population values by 0.9SDs (average across the rib's full length), 1.7SDs, and 1.3SDs, respectfully. When considering cortical bone cross-sectional area, the THUMS and VIVA models-which each define a constant bone thickness value across the entire rib-overestimated bone content on average by 1.1SDs and 1.2SDs, respectively. HBMs have traditionally performed poorly when predicting rib fracture onset or fracture site, and in all HBMs in this study the rib regions with the most extreme cortical bone thickness and cross-sectional area discrepancies (compared to average reference data) corresponded to regions toward the sternal end of the ribs where rib fractures most frequently occur. Results from this study highlight geometrical components of current HBM ribs that differ from the rib geometry that would be expected from within those models' target demographics, and help researchers prioritize improvements to their biofidelity.

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.

A Novel Approach to Scaling Age-, Sex-, and Body Size-Dependent Thoracic Responses using Structural Properties of Human Ribs

Thoracic injuries are frequently observed in motor vehicle crashes, and rib fractures are the most common of those injuries. Thoracic response targets have previously been developed from data obtained from post-mortem human subject (PMHS) tests in frontal loading conditions, most commonly of mid-size males. Traditional scaling methods are employed to identify differences in thoracic response for various demographic groups, but it is often unknown if these applications are appropriate, especially considering the limited number of tested PMHS from which those scaling factors originate. Therefore, the objective of this study was to establish a new scaling approach for generating age-, sex-, and body size- dependent thoracic responses utilizing structural properties of human ribs from direct testing of various demographics. One-hundred forty-seven human ribs (140 adult; 7 pediatric) from 132 individuals (76 male; 52 female; 4 pediatric) ranging in age from 6 to 99 years were included in this study. Ribs were tested at 2 m/s to failure in a frontal impact scenario. Force and displacement for individual ribs were used to develop new scaling factors, with a traditional mid-size biomechanical target as a baseline response. This novel use of a large, varied dataset of dynamic whole rib responses offers vast possibilities to utilize existing biomechanical data in creative ways to reduce thoracic injuries in diverse vehicle occupants.

Detailed subject-specific FE rib modeling for fracture prediction

Objective: The current state of the art human body models (HBMs) underpredict the number of fractured ribs. Also, it has not been shown that the models can predict the fracture locations. Efforts have been made to create subject specific rib models for fracture prediction, with mixed results. The aim of this study is to evaluate if subject-specific finite element (FE) rib models, based on state-of-the-art clinical CT data combined with subject-specific material data, can predict rib stiffness and fracture location in anterior-posterior rib bending.Method: High resolution clinical CT data was used to generate detailed subject-specific geometry for twelve FE models of the sixth rib. The cortical bone periosteal and endosteal surfaces were estimated based on a previously calibrated cortical bone mapping algorithm. The cortical and the trabecular bone were modeled using a hexa-block algorithm. The isotropic material model for the cortical bone in each rib model was assigned subject-specific material data based on tension coupon tests. Two different modeling strategies were used for the trabecular bone.The capability of the FE model to predict fracture location was carried out by modeling physical dynamic anterior-posterior rib bending tests. The rib model predictions were directly compared to the results from the tests. The predicted force-displacement time history, strain measurements at four locations, and rotation of the rib ends were compared to the results from the physical tests by means of CORA analysis. Rib fracture location in the FE model was estimated as the position for the element with the highest first principle strain at the time corresponding to rib fracture in the physical test.Results: Seven out of the twelve rib models predicted the fracture locations (at least for one of the trabecular modeling strategies) and had a force-displacement CORA score above 0.65. The other five rib models, had either a poor force-displacement CORA response or a poor fracture location prediction. It was observed that the stress-strain response for the coupon test for these five ribs showed significantly lower Young's modulus, yield stress, and elongation at fracture compared to the other seven ribs.Conclusion: This study indicates that rib fracture location can be predicted for subject specific rib models based on high resolution CT, when loaded in anterior-posterior bending, as long as the rib's cortical cortex is of sufficient thickness and has limited porosity. This study provides guidelines for further enhancements of rib modeling for fracture location prediction with HBMs.

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.

Development and Validation of Dummies and Human Models Used in Crash Test

The crash test dummy, an important tool for car crash safety tests, is of great significance to explore the injury biomechanics of the occupants and improve the safety performance of the vehicle. The article mainly consists of four parts: brief introduction of injury mechanism, early experiments for obtaining biomechanical response (animal tests, cadaver tests, and volunteer tests), and development and validation of mechanical dummies and computational models. This study finds that the current crash test dummies are generally designed based on European and American, so they have limitations on the damage prediction of other regions. Further research in the crash test dummy needs the participation of various countries in order to develop a crash test dummy that meets the national conditions of each country. Simultaneously, it is necessary to develop dummies of vulnerable groups, such as the elderly dummy and obese people dummy.

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.

Quantifying the force transmission through the pelvic joints during total hip arthroplasty: A pilot cadaveric study

BACKGROUND

Total hip arthroplasty is one of the most successful and cost effective procedures in orthopedics. The purpose of this study is to investigate force transmission through the sacroiliac joint as a possible source of post-operative pain after total hip arthroplasty through the following three questions: Does the ipsilateral sacroiliac joint, contralateral sacroiliac joint, or pubic symphysis experience more force during placement? Does the larger mallet used to seat the implant generate a higher force? Does the specimen's bone density or BMI alter force transmission?

METHODS

A solid design acetabular component was impacted into five human cadaver pelves with intact soft tissues. The pressure at both sacroiliac joints and the pubic symphysis was measured during cup placement. This same procedure was replicated using an existing pelvis finite element model to use for comparison.

FINDINGS

The location of the peak force for each hammer strike was found to be specimen specific. The finite model results indicated the ipsilateral sacroiliac joint had the highest pressure and strain followed by the pubic symphysis over the course of the full simulation. The heft of the mallet and bone mineral density did not predict force values or locations. The largest median force was generated in extremely obese specimens.

INTERPRETATION

Contrary to previous ideas, it is highly unlikely that forces experienced at the pelvic joints are large enough to contribute post-operative pain during impaction of an acetabular component. These results indicate more force is conveyed to the pubic symphysis compared to the sacroiliac joints.

True compression of pelvic fractures under lateral impact

PURPOSE

To promote the understanding of pelvic fracture mechanism and make more accurate evaluation of maximal deformity at the moment of fracture, kinematic response of pelvis to lateral impact and the difference between peak and final displacement were investigated.

METHODS

A total of three human cadaver pelves were seated uprightly on a sled test table, explored to horizontal lateral impact by a 22.1-kg impactor at a speed of 5.2, 4.0, and 4.8 m/s. Kinematic data of pelvic osseous interesting points (POIP) were measured by the motion capture system. Trajectories of POIP, duration of impact, and deflection of pelvis were calculated as well as rotational movement of pelvis was evaluated. After impact, autopsy and CT scan were made to validate the motion capture data.

RESULTS

The peak deflection of pelvis under lateral impact was 31.9, 30.1, and 18.5%, while final deflection was 19.6, 13.8, and 13.8%. The final deflection was only 61.5, 45.9, and 74.46% of the peak deflection.

CONCLUSIONS

In clinical practice, pelvic fracture displacement tends to be underestimated. The peak compression can be 1.3-2.2 times of final compression appearing on images in hospital. Clinicians shall give adequate estimation of displacement and related injuries.

Development of Subject-Specific Proximal Femur Finite Element Models Of Older Adults with Obesity to Evaluate the Effects of Weight Loss on Bone Strength

Study background

Recommendation of intentional weight loss in older adults remains controversial, due in part to the loss of bone mineral density (BMD) known to accompany weight loss. While finite element (FE) models have been used to assess bone strength, these methods have not been used to study the effects of weight loss. The purpose of this study is to develop subject-specific FE models of the proximal femur and study the effect of intentional weight loss on bone strength.

Methods

Computed tomography (CT) scans of the proximal femur of 25 overweight and obese (mean BMI=29.7 ± 4.0 kg/m²), older adults (mean age=65.6 ± 4.1 years) undergoing an 18-month intentional weight loss intervention were obtained at baseline and post-intervention. Measures of volumetric BMD (vBMD) and variable cortical thickness were derived from each subject CT scan and directly mapped to baseline and post-intervention models. Subject-specific FE models were developed using morphing techniques. Bone strength was estimated through simulation of a single-limb stance and sideways fall configuration.

Results

After weight loss intervention, there were significant decreases from baseline to 18 months in vBMD (total hip: -0.024 ± 0.013 g/cm³; femoral neck: -0.012 ± 0.014 g/cm³), cortical thickness (total hip: -0.044 ± 0.032 mm; femoral neck: -0.026 ± 0.039 mm), and estimated strength (stance: -0.15 ± 0.12 kN; fall: -0.04 ± 0.06 kN). Adjusting for baseline bone measures, body mass, and gender, correlations were found between weight change and change in total hip and femoral neck cortical thickness (all p<0.05). For every 1 kilogram of body mass lost cortical thickness in the total hip and femoral neck decreased by 0.003 mm and 0.004 mm, respectively. No significant correlations were present for the vBMD or strength data.

Conclusion

The developed subject-specific FE models could be used to better understand the effects of intentional weight loss on bone health.

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.

Finite Element Based Pelvic Injury Metric Creation and Validation in Lateral Impact for a Human Body Model

Pelvic fractures are serious injuries resulting in high mortality and morbidity. The objective of this study is to develop and validate local pelvic anatomical, cross-section-based injury risk metrics for a finite element (FE) model of the human body. Cross-sectional instrumentation was implemented in the pelvic region of the Global Human Body Models Consortium (GHBMC M50-O) 50th percentile detailed male FE model (v4.3). In total, 25 lateral impact FE simulations were performed using input data from cadaveric lateral impact tests performed by Bouquet et al. The experimental force-time data was scaled using five normalization techniques, which were evaluated using log rank, Wilcoxon rank sum, and correlation and analysis (CORA) testing. Survival analyses with Weibull distribution were performed on the experimental peak force (scaled and unscaled) and the simulation test data to generate injury risk curves (IRCs) for total pelvic injury. Additionally, IRCs were developed for regional injury using cross-sectional forces from the simulation results and injuries documented in the experimental autopsies. These regional IRCs were also evaluated using the receiver operator characteristic (ROC) curve analysis. Based on the results of the all the evaluation methods, the Equal Stress Equal Velocity (ESEV) and ESEV using effective mass (ESEV-EM) scaling techniques performed best. The simulation IRC shows slight under prediction of injury in comparison to these scaled experimental data curves. However, this difference was determined to not be statistically significant. Additionally, the ROC curve analysis showed moderate predictive power for all regional IRCs.

Development of the CAVEMAN Human Body Model: Validation of Lower Extremity Sub-Injurious Response to Vertical Accelerative Loading

Improving injury prediction accuracy and fidelity for mounted Warfighters has become an area of focus for the U.S. military in response to improvised explosive device (IED) use in both Iraq and Afghanistan. Although the Hybrid III anthropomorphic test device (ATD) has historically been used for crew injury analysis, it is only capable of predicting a few select skeletal injuries. The Computational Anthropomorphic Virtual Experiment Man (CAVEMAN) human body model is being developed to expand the injury analysis capability to both skeletal and soft tissues. The CAVEMAN model is built upon the Zygote 50^th percentile male human CAD model and uses a finite element modeling approach developed for high performance computing (HPC). The lower extremity subset of the CAVEMAN human body model presented herein includes: 28 bones, 26 muscles, 40 ligaments, fascia, cartilage and skin. Sensitivity studies have been conducted with the CAVEMAN lower extremity model to determine the structures critical for load transmission through the leg in the underbody blast (UBB) environment. An evaluation of the CAVEMAN lower extremity biofidelity was also carried out using 14 unique data sets derived by the Warrior Injury Assessment Manikin (WIAMan) program cadaveric lower leg testing. Extension of the CAVEMAN lower extremity model into anatomical tissue failure will provide additional injury prediction capabilities, beyond what is currently achievable using ATDs, to improve occupant survivability analyses within military vehicles.

DEVELOPMENT AND BIOFIDELITY EVALUATION OF AN OCCUPANT BIOMECHANICAL MODEL OF A CHINESE 50TH PERCENTILE MALE FOR SIDE IMPACT

Finite element modeling has played a significant role in the study of human body biomechanical responses and injury mechanisms during vehicle impacts. However, there are very few reports on similar studies conducted in China for the Chinese population. In this study, a high-precision human body finite element model of the Chinese 50th percentile male was developed. The anatomical structures and mechanical characteristics of real human body were replicated as precise as possible. In order to analyze the model’s biofidelity in side-impact injury prediction, a global technical standard, ISO/TR 9790, was used that specifically assesses the lateral impact biofidelity of anthropomorphic test devices (ATDs) and computational models. A series of model simulations, focusing on different body parts, were carried out against the tests outlined in ISO/TR 9790. Then, the biofidelity ratings of the full human body model and different body parts were evaluated using the ISO/TR 9790 rating method. In a 0–10 rating scale, the resulting rating for the full human body model developed is 8.57, which means a good biofidelity. As to different body parts, the biofidelity ratings of the head and shoulder are excellent, while those of the neck, thorax, abdomen and pelvis are good. The resulting ratings indicate that the human body model developed in this study is capable of investigating the side-impact responses of and injuries to occupants’ different body parts. In addition, the rating of the model was compared with those of the other human body finite element models and several side-impact dummy models. This allows us to assess the robustness of our model and to identify necessary improvements.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

An investigation of human body model morphing for the assessment of abdomen responses to impact against a population of test subjects

OBJECTIVE

Human body models have the potential to better describe the human anatomy and variability than dummies. However, data sets available to verify the human response to impact are typically limited in numbers, and they are not size or gender specific. The objective of this study was to investigate the use of model morphing methodologies within that context.

METHODS

In this study, a simple human model scaling methodology was developed to morph two detailed human models (Global Human Body Model Consortium models 50th male, M50, and 5th female, F05) to the dimensions of post mortem human surrogates (PMHS) used in published literature. The methodology was then successfully applied to 52 PMHS tested in 14 impact conditions loading the abdomen. The corresponding 104 simulations were compared to the responses of the PMHS and to the responses of the baseline models without scaling (28 simulations). The responses were analysed using the CORA method and peak values.

RESULTS

The results suggest that model scaling leads to an improvement of the predicted force and deflection but has more marginal effects on the predicted abdominal compressions. M50 and F05 models scaled to the same PMHS were also found to have similar external responses, but large differences were found between the two sets of models for the strain energy densities in the liver and the spleen for mid-abdomen impact simulations. These differences, which were attributed to the anatomical differences in the abdomen of the baseline models, highlight the importance of the selection of the impact condition for simulation studies, especially if the organ location is not known in the test.

CONCLUSIONS

While the methodology could be further improved, it shows the feasibility of using model scaling methodologies to compare human models of different sizes and to evaluate scaling approaches within the context of human model validation.

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.

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.

Validation of Shoulder Response of Human Body Finite-Element Model (GHBMC) Under Whole Body Lateral Impact Condition

In previous shoulder impact studies, the 50th-percentile male GHBMC human body finite-element model was shown to have good biofidelity regarding impact force, but under-predicted shoulder deflection by 80% compared to those observed in the experiment. The goal of this study was to validate the response of the GHBMC M50 model by focusing on three-dimensional shoulder kinematics under a whole-body lateral impact condition. Five modifications, focused on material properties and modeling techniques, were introduced into the model and a supplementary sensitivity analysis was done to determine the influence of each modification to the biomechanical response of the body. The modified model predicted substantially improved shoulder response and peak shoulder deflection within 10% of the observed experimental data, and showed good correlation in the scapula kinematics on sagittal and transverse planes. The improvement in the biofidelity of the shoulder region was mainly due to the modifications of material properties of muscle, the acromioclavicular joint, and the attachment region between the pectoralis major and ribs. Predictions of rib fracture and chest deflection were also improved because of these modifications.

Construction and evaluation of thoracic injury risk curves for a finite element human body model in frontal car crashes

There is a need to improve the protection to the thorax of occupants in frontal car crashes. Finite element human body models are a more detailed representation of humans than anthropomorphic test devices (ATDs). On the other hand, there is no clear consensus on the injury criteria and the thresholds to use with finite element human body models to predict rib fractures. The objective of this study was to establish a set of injury risk curves to predict rib fractures using a modified Total HUman Model for Safety (THUMS). Injury criteria at the global, structural and material levels were computed with a modified THUMS in matched Post Mortem Human Subjects (PMHSs) tests. Finally, the quality of each injury risk curve was determined. For the included PMHS tests and the modified THUMS, DcTHOR and shear stress were the criteria at the global and material levels that reached an acceptable quality. The injury risk curves at the structural level did not reach an acceptable quality.

Comparison of Kriging and Moving Least Square Methods to Change the Geometry of Human Body Models

Finite Element Human Body Models (HBM) have become powerful tools to study the response to impact. However, they are typically only developed for a limited number of sizes and ages. Various approaches driven by control points have been reported in the literature for the non-linear scaling of these HBM into models with different geometrical characteristics. The purpose of this study is to compare the performances of commonly used control points based interpolation methods in different usage scenarios. Performance metrics include the respect of target, the mesh quality and the runability. For this study, the Kriging and Moving Least square interpolation approaches were compared in three test cases. The first two cases correspond to changes of anthropometric dimensions of (1) a child model (from 6 to 1.5 years old) and (2) the GHBMC M50 model (Global Human Body Models Consortium, from 50th to 5th percentile female). For the third case, the GHBMC M50 ribcage was scaled to match the rib cage geometry derived from a CT-scan. In the first two test cases, all tested methods provided similar shapes with acceptable results in terms of time needed for the deformation (a few minutes at most), overall respect of the targets, element quality distribution and time step for explicit simulation. The personalization of rib cage proved to be much more challenging. None of the methods tested provided fully satisfactory results at the level of the rib trajectory and section. There were corrugated local deformations unless using a smooth regression through relaxation. Overall, the results highlight the importance of the target definition over the interpolation method.

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.

Development of a computationally efficient full human body finite element model

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Age- and sex-specific thorax finite element model development and simulation

OBJECTIVE

The shape, size, bone density, and cortical thickness of the thoracic skeleton vary significantly with age and sex, which can affect the injury tolerance, especially in at-risk populations such as the elderly. Computational modeling has emerged as a powerful and versatile tool to assess injury risk. However, current computational models only represent certain ages and sexes in the population. The purpose of this study was to morph an existing finite element (FE) model of the thorax to depict thorax morphology for males and females of ages 30 and 70 years old (YO) and to investigate the effect on injury risk.

METHODS

Age- and sex-specific FE models were developed using thin-plate spline interpolation. In order to execute the thin-plate spline interpolation, homologous landmarks on the reference, target, and FE model are required. An image segmentation and registration algorithm was used to collect homologous rib and sternum landmark data from males and females aged 0-100 years. The Generalized Procrustes Analysis was applied to the homologous landmark data to quantify age- and sex-specific isolated shape changes in the thorax. The Global Human Body Models Consortium (GHBMC) 50th percentile male occupant model was morphed to create age- and sex-specific thoracic shape change models (scaled to a 50th percentile male size). To evaluate the thoracic response, 2 loading cases (frontal hub impact and lateral impact) were simulated to assess the importance of geometric and material property changes with age and sex.

RESULTS

Due to the geometric and material property changes with age and sex, there were observed differences in the response of the thorax in both the frontal and lateral impacts. Material property changes alone had little to no effect on the maximum thoracic force or the maximum percent compression. With age, the thorax becomes stiffer due to superior rotation of the ribs, which can result in increased bone strain that can increase the risk of fracture. For the 70-YO models, the simulations predicted a higher number of rib fractures in comparison to the 30-YO models. The male models experienced more superior rotation of the ribs in comparison to the female models, which resulted in a higher number of rib fractures for the males.

CONCLUSION

In this study, age- and sex-specific thoracic models were developed and the biomechanical response was studied using frontal and lateral impact simulations. The development of these age- and sex-specific FE models of the thorax will lead to an improved understanding of the complex relationship between thoracic geometry, age, sex, and injury risk.

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

OBJECTIVE

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

METHOD

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

RESULTS

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

CONCLUSIONS

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

Biomechanical model of the thorax under blast loading: a three dimensional numerical study

Injury mechanisms due to high speed dynamic loads, such as blasts, are not well understood. These research fields are widely investigated in the literature, both at the experimental and numerical levels, and try to answer questions about the safety and efficiency of protection devices or biomechanical traumas. At a numerical level, the development of powerful mathematical models tends to study tolerance limits and injury mechanisms in order to avoid experimental tests which cannot be easily conducted. In a military framework, developing a fighter/soldier numerical model can help to the understanding of many traumas which are specific to soldier injuries, like mines, ballistic impacts or blast traumas. The aim of this study is to investigate the consequences of violent loads in terms of human body response, submitting a developed and validated three-dimensional thorax finite element (FE) model to blast loadings. Specific formulations of FE methods are used to simulate this loading, and its consequence on the biomechanical model. Mechanical parameters such as pressure in the air field and also in internal organs are observed, and these values are compared to the experimental data in the literature. This study gives encouraging results and allows going further in soldier trauma investigations.

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.

Effects of cervical arthrodesis and arthroplasty on neck response during a simulated frontal automobile collision

BACKGROUND CONTEXT

Whereas arthrodesis is the most common surgical intervention for the treatment of symptomatic cervical degenerative disc disease, arthroplasty has become increasingly more popular over the past decade. Although literature exists comparing the effects of anterior cervical discectomy and fusion and cervical total disc replacement (CTDR) on neck kinematics and loading, the vast majority of these studies apply only quasi-static, noninjurious loading conditions to a segment of the cervical spine.

PURPOSE

The objective of this study was to investigate the effects of arthrodesis and arthroplasty on biomechanical neck response during a simulated frontal automobile collision with air bag deployment.

STUDY DESIGN

This study used a full-body, 50th percentile seated male finite element (FE) model to evaluate neck response during a dynamic impact event. The cervical spine was modified to simulate either an arthrodesis or arthroplasty procedure at C5-C6.

METHODS

Five simulations of a belted driver, subjected to a 13.3 m/s ΔV frontal impact with air bag deployment, were run in LS-DYNA with the Global Human Body Models Consortium full-body FE model. The first simulation used the original model, with no modifications to the neck, whereas the remaining four were modified to represent either interbody arthrodesis or arthroplasty of C5-C6. Cross-sectional forces and moments at the C5 and C6 cervical levels of the neck, along with interbody and facet forces between C5 and C6, were reported.

RESULTS

Adjacent-level, cross-sectional neck loading was maintained in all simulations without exceeding any established injury thresholds. Interbody compression was greatest for the CTDRs, and interbody tension occurred only in the fused and nonmodified spines. Some interbody separation occurred between the superior and inferior components of the CTDRs during flexion-induced tension of the cervical spine, increasing the facet loads.

CONCLUSIONS

This study evaluated the effects of C5-C6 cervical arthrodesis and arthroplasty on neck response during a simulated frontal automobile impact. Although cervical arthrodesis and arthroplasty at C5-C6 did not appear to significantly alter the adjacent-level, cross-sectional neck responses during a simulated frontal automobile impact, key differences were noted in the interbody and facet loading.

Deflection Corridors of Abdomen and Thorax in Oblique Side Impacts Using Equal Stress Equal Velocity Approach: Comparison With Other Normalization Methods

The first objective of the study was to determine the thorax and abdomen deflection time corridors using the equal stress equal velocity approach from oblique side impact sled tests with postmortem human surrogates fitted with chestbands. The second purpose of the study was to generate deflection time corridors using impulse momentum methods and determine which of these methods best suits the data. An anthropometry-specific load wall was used. Individual surrogate responses were normalized to standard midsize male anthropometry. Corridors from the equal stress equal velocity approach were very similar to those from impulse momentum methods, thus either method can be used for this data. Present mean and plus/minus one standard deviation abdomen and thorax deflection time corridors can be used to evaluate dummies and validate complex human body finite element models.

Development of structural and material clavicle response corridors under axial compression and three point bending loading for clavicle finite element model validation

Clavicle injuries were frequently observed in automotive side and frontal crashes. Finite element (FE) models have been developed to understand the injury mechanism, although no clavicle loading response corridors yet exist in the literature to ensure the model response biofidelity. Moreover, the typically developed structural level (e.g., force-deflection) response corridors were shown to be insufficient for verifying the injury prediction capacity of FE model, which usually is based on strain related injury criteria. Therefore, the purpose of this study is to develop both the structural (force vs deflection) and material level (strain vs force) clavicle response corridors for validating FE models for injury risk modeling. 20 Clavicles were loaded to failure under loading conditions representative of side and frontal crashes respectively, half of which in axial compression, and the other half in three point bending. Both structural and material response corridors were developed for each loading condition. FE model that can accurately predict structural response and strain level provides a more useful tool in injury risk modeling and prediction. The corridor development method in this study could also be extended to develop corridors for other components of the human body.