A Finite Element Model of the Lower Limb for Simulating Automotive Impacts

Modeling and simulation credibility assessments of whole-body finite element computational models for use in NASA extravehicular activity applications

Computational finite element (FE) models are used in suited astronaut injury risk assessments; however, these models' verification, validation, and credibility (VV&C) procedures for simulating injuries in altered gravity environments are limited. Our study conducts VV&C assessments of THUMS and Elemance whole-body FE models for predicting suited astronaut injury biomechanics using eight credibility factors, as per NASA-STD-7009A. Credibility factor ordinal scores are assigned by reviewing existing documentation describing VV&C practices, and credibility sufficiency thresholds are assigned based on input from subject matter experts. Our results show the FE models are credible for suited astronaut injury investigation in specific ranges of kinematic and kinetic conditions correlating to highway and contact sports events. Nevertheless, these models are deficient when applied outside these ranges. Several credibility elevation strategies are prescribed to improve models' credibility for the NASA-centric application domain.

Mapping the knowledge of traffic collision Reconstruction: A scientometric analysis in CiteSpace, VOSviewer, and SciMAT

Traffic collisions are incidents with high fatality rate which generate billions of US dollars of loss worldwide each year. Post-collision scene reconstruction, which involves knowledge of multiple disciplines, is an important approach to restore the traffic collision and infer the cause of it. This paper uses software CiteSpace, VOSviewer, and SciMAT to conduct a visualization study of knowledge mapping on the literature of traffic collision scene reconstruction from 2001 to 2021 based on the Web of Science database. Knowledge mapping is a cutting-edge research method in scientometric, which has been widely applied in medicine and informatics. Compared with traditional literature review, knowledge mapping with visual techniques identifies hot keywords and key literature in the field more scientifically, and displays them in schematic diagrams intuitively which allows to further predict potential hotspots. A total of 803 original papers are retrieved to analyze and discuss the evolution of the field in the past 20 years, from macro to micro, in term of background information, popular themes, and knowledge structure. Results indicate the number of publications in this field is limited, and collaborations among authors and among institutions are insufficient. In the meantime, mappings imply the top three hot themes being scene reconstruction, computer technology, and injuries. The introduction of AI related technologies, such as neural networks and genetic algorithms, into collision reconstruction would be a potential research direction.

Cortical bone continuum damage mechanics constitutive model with stress triaxiality criterion to predict fracture initiation and pattern

A primary objective of finite element human body models (HBMs) is to predict response and injury risk in impact scenarios, including cortical bone fracture initiation, fracture pattern, and the potential to simulate post-fracture injury to underlying soft tissues. Current HBMs have been challenged to predict the onset of failure and bone fracture patterns owing to the use of simplified failure criteria. In the present study, a continuum damage mechanics (CDM) model, incorporating observed mechanical response (orthotropy, asymmetry, damage), was coupled to a novel phenomenological effective strain fracture criterion based on stress triaxiality and investigated to predict cortical bone response under different modes of loading. Three loading cases were assessed: a coupon level notched shear test, whole bone femur three-point bending, and whole bone femur axial torsion. The proposed material model and fracture criterion were able to predict both the fracture initiation and location, and the fracture pattern for whole bone and specimen level tests, within the variability of the reported experiments. There was a dependence of fracture threshold on finite element mesh size, where higher mesh density produced similar but more refined fracture patterns compared to coarser meshes. Importantly, the model was functional, accurate, and numerically stable even for relatively coarse mesh sizes used in contemporary HBMs. The proposed model and novel fracture criterion enable prediction of fracture initiation and resulting fracture pattern in cortical bone such that post-fracture response can be investigated in HBMs.

Age-Related Study and Collision Response of Material Properties of Long Bones in Chinese Pedestrian Lower Limbs

In forensic examination cases, lower limb injuries are common, and pedestrians of different ages suffer different injuries when they are hit by vehicles, especially the injuries to the long bones of the lower limbs. Aging remains a challenging issue for the material properties and injury biomechanical properties of pedestrian lower limb long bones. We analyzed the regression relationship between the age of 50 Chinese pedestrians and the material properties of the lower limb long bones (femur, tibia). We compared them with previous studies to propose a regression model suitable for Chinese human long bone material properties. Through the established Human Active Lower Limb (HALL) model that conforms to the Chinese human anatomy, seven pedestrians’ (20/30/40/50/60/70/80 years old (YO)) lower limbs were parameterized to assign long bone material properties. In the finite element analysis, the Hall model was side-impacted by a family car (FCR) at speeds of 30/40/50/60/70 km/h, respectively. The results showed that an increase in age was negatively correlated with a decrease in the material properties of each long bone. Moreover, with an increase in age, the tolerance limit of long bones gradually decreases, but there will be a limit, and there is no obvious positive correlation with age. During a standing side impact, the stress change in the femur was significantly smaller than that of the tibia, and the stress of the femur and tibia decreased with age. Age is a more significant influencing factor for lower limb injuries. Older pedestrians have a higher risk of lower limb injuries. Forensic experts should pay attention to the critical factor of age when encountering lower limb traffic accident injuries in forensic identification work.

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.

A Review of Finite Element Models of Ligaments in the Foot and Considerations for Practical Application

PURPOSE

Finite element (FE) modeling has been used as a research tool for investigating underlying ligaments biomechanics and orthopedic applications. However, FE models of the ligament in the foot have been developed with various configurations, mainly due to their complex 3D geometry, material properties, and boundary conditions. Therefore, the purpose of this review was to summarize the current state of finite element modeling approaches that have been used in the ?eld of ligament biomechanics, to discuss their applicability to foot ligament modeling in a practical setting, and also to acknowledge current limitations and challenges.

METHODS

A comprehensive literature search was performed. Each article was analyzed in terms of the methods used for: (a) ligament geometry, (b) material property, (c) boundary and loading condition related to its application, and (d) model verification and validation.

RESULTS

Of the reviewed studies, 80% of the studies used simplified representations of ligament geometry, the non-linear mechanical behavior of ligaments was taken into account in only 19.2% of the studies, 33% of included studies did not include any kind of validation of the FE model.

CONCLUSION

Further refinement in the functional modeling of ligaments, the micro-structure level characteristics, nonlinearity, and time-dependent response, may be warranted to ensure the predictive ability of the models.

Preliminary study on the mechanisms of ankle injuries under falling and impact conditions based on the THUMS model

Ankle injuries are common in forensic practice, which are mainly caused by falling and traffic accidents. Determining the mechanisms and manners of ankle injuries is a critical and challenging problem for forensic experts. The identification of the injury mechanism is still experience-based and strongly subjective. There also lacks systematic research in current practice. In our study, based on the widely used Total Human Model of Safety 4.0 (THUMS 4.0), we utilized the finite element (FE) method to simulate ankle injuries caused by falls from different heights (5 m, 10 m and 20 m) with different landing postures (natural posture, inversion, eversion, plantar-flexion and dorsi-flexion) and injuries caused by impacts from different directions (anterior-posterior, lateral-medial and posterior-anterior) with different speeds (10 m/s, 15 m/s and 20 m/s) at different sites (ankle and lower, middle and upper sections of leg). We compared the injury morphology and analyzed the mechanisms of ankle injuries. The results showed that falling causes a specific compression fracture of the distal tibia, while fractures of the tibia and fibula diaphysis and ligament injuries caused by falling from a lower height or inversion, planter flexion or dorsiflexion at a large angle are not distinguishable from the similar injury patterns caused by impact on the middle and upper segments of the leg. No obvious compression fracture of the tibia distal was caused by the impacts, whereas ligament injuries and avulsion fractures of the medial or lateral condyle and fractures of the diaphysis of the tibia and fibula were observed. Systematic studies will be helpful in reconstructing the ankle injury processes and analyzing the mechanisms in forensic practice, providing a deeper understanding of ankle injury mechanisms for forensic experts.

Development and Application of Digital Human Models in the Field of Vehicle Collisions: A Review

In the human-vehicle-road system of collisions, the human is the most important factor, and digital human models (DHMs) are developed with the aim of preventing or at least reducing human injury. Because most of the relevant literature is focused mainly on collisions in traffic accidents (TAs), only some of the literature reviewed in this paper involves research results on other aspects of collisions. In this review, based on the background of DHMs and the application of DHMs regarding human injury biomechanics in collisions field, research results regarding the development of DHMs are described, the methods for verifying such models are introduced, and the application of the research results is discussed based on the aspect of human injury biomechanics. From the research literature, the development and validation of DHMs and their application in human injury biomechanics are summarized, and future research trends are proposed and discussed.

Study on the Long Bone Failure Behaviors Under the Indenter Rigid-Contact by Experiment Analysis and Subject-Specific Simulation

The bending fracture behaviors of long bone have gained great attention due to the high bending fracture risk during sports events, traffic accidents, and falling incidents, etc. For evaluating bone bending behaviors, most of the previous studies used an indenter in three point bending experiments while the effect of its rigidity was never considered. In this work, using the porcine long bones, the three point bending tests were conducted to explore the bone fracture behaviors under a rigid indenter. In addition to collecting the force applied, the bone fracture dynamic process was recorded by high-speed photography, and the fracture surface profile in mesoscale was observed by the scanning electron microscope (SEM). Based on CT scanning of long bones, the cross section properties of test specimens were calculated by a homemade matlab script for correlating with their failure strengths. Also, a subject-specific finite element (FE) model was developed to identify the outcomes induced by a rigid indenter on simulation. Findings led to conclusions as follows: (1) The tension fracture came with fracture path deflection, which was caused by the bone indentation induced mesoscale crack-opening. Due to this damage before the whole bone fracture, a bone fracture moment correction was established to compensate experimental data. (2) The plastic indentation caused the force fluctuation as suggested by correlation analysis. (3) The bone failure moment correlated with the inertial moment of the bone cross section at the fracture location higher than the traditional cross section area. (4) In the subject-specific simulation, the indentation caused compression fracture under a much lower failure force. Removing the element erosion on the indenter-contacted area only during the validation was verified as a good option to solve this issue.

Biomechanical evaluation of the stability of extra-articular distal radius fractures fixed with volar locking plates according to the length of the distal locking screw

Surgeons usually used short screws to avoid extensor tendon problems during volar locking plate fixation in distal radius fracture. However, the stability according to the length of distal locking screws have not been fully understood. We investigated this issue through finite element analysis and compression test using synthetic radius. Our results demonstrated that the bi-cortical full-length fixation does not contribute to the stiffness increase in the axial compression direction, and a reduction in length of up to more than 50% length can still provide similar stability to full-length screws. Our data can support that surgeon should undersize the distal screw.

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.

Monte carlo method for estimating whole-body injury metrics from pedestrian impact simulation results

The goal of the current study was to develop a method to estimate whole-body injury metrics (WBIMs), which measure the overall impact of injuries, using stochastic injury prediction results from a computational human surrogate. First, hospitalized pedestrian data was queried to identify injuries sustained by pedestrians and their frequencies. Second, with consideration for an understanding of injury mechanisms and the capability of the computational human surrogate, the whole-body was divided into 17 body regions. Then, an injury pattern database was constructed for each body region for various maximum abbreviated injury scale (MAIS) levels. Third, a two-step Monte Carlo sampling process was employed to generate N virtual pedestrians with an assigned list of injuries in AIS codes. Then, the expected values of WBIMs such as injury severity score (ISS), probability of death, whole-body functional capacity index (WBFCI), and lost years of life (LYL), were estimated. Lastly, the proposed method was verified using injury information from the inpatient pedestrian database. Also, the proposed method was applied to pedestrian impact simulations with various impact speeds to estimate the probability of death with respect to the impact speed. The probability of death from the proposed method was compared with those from epidemiological studies. The proposed method accurately estimated WBIMs such as ISS and WBFCI using either for a given distribution of injury risk or MAIS levels. The predicted probability of death with respect to the impact speed showed a good correlation with those from the epidemiological study. These results imply that if we have a human surrogate that can predict the risk of injury accurately, we can accurately estimate WBIMs using the proposed method. The proposed method can simplify a vehicle design optimization process by transforming the multi-objective optimization problem into the single-objective one. Lastly, the proposed method can be applied to other human surrogates such as occupant models.

An examination of the performance of damaged energy-absorbing end terminals

OBJECTIVES

Guardrail end terminals are designed to gradually decelerate vehicles during impact and protect vehicle occupants from severe injuries. It has been observed that some in-service end terminals are damaged, and it is unclear if their safety performance is still acceptable. The objectives of this study were to examine the conditions of in-service end terminals, and to evaluate the performance of damaged relative to undamaged end terminals in simulated impacts.

METHODS

Common damage patterns of guardrail end terminals were investigated by using post-crash pictures collected from the National Automotive Sampling System-Crashworthiness Data System (NASS-CDS). Conditions of in-service end terminals mounted along roads in portions of six U.S. states were examined by using a sample from the second Strategic Highway Research Program-Roadway Information Database (SHRP2-RID). Finite Element (FE) models of two minorly and three severely damaged ET-Plus systems, a commonly used energy-absorbing guardrail end terminal along U.S. roads, were developed. To evaluate the performance of the damaged ET-Plus systems, we performed impact simulations with vehicle-to-damaged ET-Plus models according to the National Cooperative Highway Research Program (NCHRP) 350, test conditions 3-30.

RESULTS

Of the 1000 in-service end terminal cases we investigated, 73% were undamaged, 18% had minor damage, and 8% had major damage. Increases in the average vehicle deceleration rates, maximum vehicle yaw angles, and vehicle local deformations were observed in simulated impacts with damaged ET-Plus end terminals relative to impacts with undamaged ET-Plus end terminals. For one damaged ET-Plus, a secondary collision was observed. Overall, we found that the damaged end terminals usually increased collision severity when compared with undamaged end terminals.

CONCLUSIONS

The findings of this study point out the need for in-service performance evaluations and proper maintenance and repair practices of end terminals. The simulation models developed in this study could be further employed to investigate device performance in crash situations that are physically impractical to test and investigate the effects of site characteristics on device performance. The simulation models could also supplement crash tests to certify new hardware designs.

The Mechanical Response and Tolerance of the Anteriorly-Tilted Human Pelvis Under Vertical Loading

Military vehicle underbody blast (UBB) is the cause of many serious injuries in theatre today; however, the effects of these chaotic events on the human body are not well understood. The purpose of this research was to replicate UBB loading conditions on the human pelvis and investigate the resulting response in a controlled laboratory setting. In addition to better understanding the response of the human pelvis to high rate vertical loading, this test series also aimed to identify high rate injury thresholds. Twenty-seven post mortem human surrogate (PMHS) component pelvis tests were completed using the University of Virginia's (UVa) simulated blast rig under a range of loading conditions and postures. Of those tests, 17 were in the anteriorly-tilted posture and used to construct the human injury probability curve. Average seat pan (rigid) accelerations for this subset of tests ranged from 300 to 2400 g over 2 to 3 ms of positive phase duration. Post-test computed tomography (CT) scans and necropsies were performed to determine injuries and revealed a frequent occurrence of anterior and posterior injuries, resulting in unstable pelvis ring fractures. The resulting Human Injury Probability Curve (HIPC) yielded mean forces of 5529, 8516, and 12431 N as measured by mass compensated seat platen loadcells applied through the rigid seat to the bilateral ischium are associated with a 10, 25, and 50% risk for unstable pelvic ring sacrum fractures in an anteriorly-tilted pelvis (28° from vertical), respectively.

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.

Implementation of controlling strategy in a biomechanical lower limb model with active muscles for coupling multibody dynamics and finite element analysis

Computational biomechanics for human body modeling has generally been categorized into two separated domains: finite element analysis and multibody dynamics. Combining the advantages of both domains is necessary when tissue stress and physical body motion are both of interest. However, the method for this topic is still in exploration. The aim of this study is to implement unique controlling strategies in finite element model for simultaneously simulating musculoskeletal body dynamics and in vivo stress inside human tissues. A finite element lower limb model with 3D active muscles was selected for the implementation of controlling strategies, which was further validated against in-vivo human motion experiments. A unique feedback control strategy that couples together a basic Proportion-Integration-Differentiation (PID) controller and generic active signals from Computed Muscle Control (CMC) method of the musculoskeletal model or normalized EMG singles was proposed and applied in the present model. The results show that the new proposed controlling strategy show a good correlation with experimental test data of the normal gait considering joint kinematics, while stress distribution of local lower limb tissue can be also detected in real-time with lower limb motion. In summary, the present work is the first step for the application of active controlling strategy in the finite element model for concurrent simulation of both body dynamics and tissue stress. In the future, the present method can be further developed to apply it in various fields for human biomechanical analysis to monitor local stress and strain distribution by simultaneously simulating human locomotion.

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.

Effect of Exercise Modality During Weight Loss on Bone Health in Older Adults With Obesity and Cardiovascular Disease or Metabolic Syndrome: A Randomized Controlled Trial

The objective of this study was to determine the ability of either aerobic or resistance training to counter weight-loss-associated bone loss in older adults. There were 187 older adults (67 years, 70% women, 64% white) with obesity (BMI = 34.5 ± 3.7 kg/m² ) and cardiovascular disease and/or metabolic syndrome who were randomized to participate in an 18-month, community-based trial, with a follow-up assessment at 30 months. Intervention arms included: weight loss alone (WL; 7% to 10% baseline weight), WL plus aerobic training (WL + AT), and WL plus resistance training (WL + RT), as well as DXA-acquired total hip, femoral neck, and lumbar spine areal bone mineral density (aBMD), and trabecular bone score (TBS). Biomarkers of bone turnover (procollagen type 1 N-terminal propeptide, C-terminal telopeptide of type 1 collagen) were measured at baseline, 6, 18, and 30 (aBMD and TBS only) months. CT-acquired hip and spine volumetric BMD (vBMD), cortical thickness, and bone strength were measured in a subset at baseline (n = 55) and 18 months. Total hip aBMD was reduced by 2% in all groups at 18 months, with a primary analysis showing no significant treatment effects for any DXA, biomarker, or CT outcome. After adjustment for WL and follow-up at 30 months, secondary analyses revealed that total hip [-0.018 (-0.023 to -0.012) g/cm² versus -0.025 (-0.031 to -0.019) g/cm² ; p = 0.05] and femoral neck [-0.01 (-0.009 to 0.008) g/cm² versus -0.011 (-0.020 to -0.002) g/cm² ; p = 0.06] aBMD estimates were modestly attenuated in the WL + RT group compared with the WL group. Additionally, lumbar spine aBMD was increased in the WL [0.015 (0.007 to 0.024) g/cm² ] and the WL + RT [0.009 (0.000 to 0.017) g/cm² ] groups compared with the WL + AT [-0.003 (-0.012 to 0.005)g/cm² ] group; both p ≤ 0.01. Community-based exercise does not prevent bone loss during active WL in older adults; however, adding RT may help minimize long-term hip bone loss. © 2018 American Society for Bone and Mineral Research.

A finite element model of an anthropomorphic test device lower limb to assess risk of injuries during vertical accelerative loading

Improvised explosive devices (IEDs) were used extensively to target occupants of military vehicles during the conflicts in Iraq and Afghanistan (2003-2011). War fighters exposed to an IED attack were highly susceptible to lower limb injuries. To appropriately assess vehicle safety and make informed improvements to vehicle design, a novel Anthropomorphic Test Device (ATD), called the Warrior Injury Assessment Manikin (WIAMan), was designed for vertical loading. The main objective of this study was to develop and validate a Finite Element (FE) model of the WIAMan lower limb (WIAMan-LL). Appropriate materials and contacts were applied to realistically model the physical dummy. Validation of the model was conducted based on experiments performed on two different test rigs designed to simulate the vertical loading experienced during an under-vehicle explosion. Additionally, a preliminary evaluation of the WIAMan and Hybrid-III test devices was performed by comparing force responses to post-mortem human surrogate (PMHS) corridors. The knee axial force recorded by the WIAMan-LL when struck on the plantar surface of the foot (2 m/s) fell mostly within the PMHS corridor, but the corresponding data predicted by the Hybrid-III was almost 60% higher. Overall, good agreements were observed between the WIAMan-LL FE predictions and experiments at various pre-impact speeds ranging from 2 m/s up to 5.8 m/s. Results of the FE model were backed by mean objective rating scores of 0.67-0.76 which support its accuracy relative to the physical lower limb dummy. The observations and objective rating scores show the model is validated within the experimental loading conditions. These results indicate the model can be used in numerical studies related to possible dummy design improvements once additional PMHS data is available. The numerical lower limb is currently incorporated into a whole body model that will be used to evaluate the vehicle design for underbody blast protection.

A Review of Pediatric Lower Extremity Data for Pedestrian Numerical Modeling: Injury Epidemiology, Anatomy, Anthropometry, Structural, and Mechanical Properties

Pedestrian injuries are the fourth leading cause of unintentional injury-related death among children aged 1 to 19. The lower extremity represents the most frequently injured body region in car-to-pedestrian accidents. The goal of this study was to perform a systematic review of the data related to pedestrian lower extremity injuries, anatomy, anthropometry, structural, and mechanical properties, which can be used in the development of new pediatric computational models. The study began with a review of epidemiologic data related to pediatric pedestrian accidents. Anatomy of the child lower extremity and age-related anthropometry data were presented as well. Then, both the mechanical and structural properties of the lower extremity main components (e.g., bones, cartilages, knee ligaments, muscles, tendons, and growth plates) available in literature were summarized. The study concluded with a brief description of current child pedestrian models, which included a discussion about their limitations. We believe that data included in this review study can help in improving the biofidelity of current child models and support the development and validation of new child models used by safety researchers for protection of pediatric population.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

A methodology for the generation and non-destructive characterisation of transverse fractures in long bones

Long bone fractures are common and although treatments are highly effective in most cases, it is challenging to achieve successful repair for groups such as open and periprosthetic fractures. Previous biomechanical studies of fracture repair, including computer and experimental models, have simplified the fracture with a flat geometry or a gap, and there is a need for a more accurate fracture representation to mimic the situation in-vivo. The aims of this study were to develop a methodology for generating repeatable transverse fractures in long bones in-vitro and to characterise the fracture surface using non-invasive computer tomography (CT) methods. Ten porcine femora were fractured in a custom-built rig under high-rate loading conditions to generate consistent transverse fractures (angle to femoral axis < 30 degrees). The bones were imaged using high resolution peripheral quantitative CT (HR-pQCT). A method was developed to extract the roughness and form profiles of the fracture surface from the image data using custom code and Guassian filters. The method was tested and validated using artificially generated waveforms. The results revealed that the smoothing algorithm used in the script was robust but the optimum kernel size has to be considered.

•

A new method to generate and quantify long bone fracture surface geometry is proposed.

•

Reliable transverse fractures were generated and imaged using microCT.

•

A non-destructive method to characterise the surface from the images was developed

•

Validation of the method using benchmark waveforms showed the method to be robust

•

By tuning the filtering conditions, specific characteristics of the surface can be extracted.

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.

Development and Validation of an Age-Specific Lower Extremity Finite Element Model for Simulating Pedestrian Accidents

The objective of the present study is to develop an age-specific lower extremity finite element model for pedestrian accident simulation. Finite element (FE) models have been used as a versatile tool to simulate and understand the pedestrian injury mechanisms and assess injury risk during crashes. However, current computational models only represent certain ages in the population, the age spectrum of the pedestrian victims is very large, and the geometry of anatomical structures and material property of the lower extremities changes with age for adults, which could affect the injury tolerance, especially in at-risk populations such as the elderly. The effects of age on the material mechanical property of bone and soft tissues of the lower extremities as well as the geometry of the long bone were studied. Then an existing 50th percentile male pedestrian lower extremity model was rebuilt to depict lower extremity morphology for 30- to 70-year-old (YO) individuals. A series of PMHS tests were simulated to validate the biofidelity and stability of the created age-specific models and evaluate the lower extremity response. The development of age-specific lower extremity models will lead to an improved understanding of the pedestrian lower extremity injury mechanisms and injury risk prediction for the whole population in vehicle-pedestrian collision accidents.

Skin mechanical properties and modeling: A review

The mechanical properties of the skin are important for various applications. Numerous tests have been conducted to characterize the mechanical behavior of this tissue, and this article presents a review on different experimental methods used. A discussion on the general mechanical behavior of the skin, including nonlinearity, viscoelasticity, anisotropy, loading history dependency, failure properties, and aging effects, is presented. Finally, commonly used constitutive models for simulating the mechanical response of skin are discussed in the context of representing the empirically observed behavior.

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.

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 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.

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

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

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

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

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.

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.

A statistical geometrical description of the human liver for probabilistic occupant models

Realistic numerical assessments of liver injury risk for the entire occupant population require incorporating inter-subject variations into numerical models. Statistical shape models of the abdominal organs have been shown to be useful tools for the investigation of the organ variations and could be applied to the development of statistical computational models. The main objective of this study was to establish a standard procedure to quantify the shape variations of a human liver in a seated posture, and construct three-dimensional (3D) statistical shape boundary models. Statistical shape analysis was applied to construct shape models of 15 adult human livers. Principal component analysis (PCA) was then utilized to obtain the modes of variation, the mean model, and a set of statistical boundary shape models, which were constructed using the q-hyper-ellipsoid approach. The first five modes of a human liver accounted for the major anatomical variations. The modes were highly correlated to the height, thickness, width, and curvature of the liver, and the concavity of the right lobe. The mean model and the principal components were utilized to construct four boundary models of human liver. The statistical boundary model approach presented in this study could be used to develop probabilistic finite element (FE) models. In the future, the probabilistic liver models could be used in FE simulations to better understand the variability in biomechanical responses and abdominal injuries under impact loading.

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.

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.

Development, calibration, and validation of a head-neck complex of THOR mod kit finite element model

INTRODUCTION/OBJECTIVE

In an effort to continually improve upon the design of the test device for human occupant restraint (THOR) dummy, a series of modifications have recently been applied. The first objective of this study was to update the THOR head-neck finite element (FE) model to the specifications of the latest dummy modifications. The second objective was to develop and apply a new optimization-based methodology to calibrate the FE head-neck model based on experimental test data. The calibrated head-neck model was validated against both frontal and lateral impact test data. Finally, the sensitivities of the model, in terms of head and neck injury criteria, to pretest positioning conditions were evaluated in a frontal crash test simulation.

METHODS

The updated parts of the head-neck THOR FE model were remeshed from CAD geometries of the modified parts. In addition, further model modifications were made to improve the effectiveness of the model (e.g., model stability). A novel calibration methodology, which incorporates the CORA (CORelation and Analysis) rating system with an optimization algorithm implemented in Isight software, was developed to improve both kinematic and kinetic responses of the model in various THOR dummy certification and biomechanical response tests. A parametric study was performed to evaluate head and neck injury criteria values in the calibrated head-neck model during a 40 km/h frontal crash test with respect to variation in the THOR model upper body and belt pretest position.

RESULTS

Material parameter optimization was shown to greatly improve the updated model response by increasing the average rating score from 0.794 ± 0.073 to 0.964 ± 0.019. The calibrated neck showed the biggest improvement in the pendulum flexion simulation from 0.681 in the original model up to 0.96 in the calibrated model. The fully calibrated model proved to be effective at predicting dummy response in frontal and lateral loading conditions during the validation phase (0.942 average score). Upper body position was shown to have a greater effect on head-neck response than belt position. The pretest positioning variation resulted in a 10 percent maximum change in HIC36 values and 14 percent maximum change in NIJ values.

CONCLUSION

The optimization-based calibration methodology was effective as it markedly improved model performance. The calibrated head-neck model demonstrated application in a crash safety analysis, showing slight head-neck injury sensitivity to pretest positioning in a frontal crash impact scenario.

A numerical investigation on the variation in hip injury tolerance with occupant posture during frontal collisions

OBJECTIVE

More than half of occupant lower extremity (LEX) injuries during automotive frontal crashes are in the knee-thigh-hip (KTH) complex. The objective of this study is to develop a detailed and biofidelic finite element (FE) occupant LEX model that may improve current understanding of mechanisms and thresholds of KTH injuries.

METHODS

Firstly, the pelvis, thigh-knee-hip, and foot models developed in our previous studies were connected into an occupant lower limb model. Further validations, including posterior cruciate ligament (PCL) stretching, thigh lateral loading, KT, and KTH impact loading were then performed to verify the injury predictability of the model under complex frontal and lateral loading corresponding to automotive impacts. Finally, a sensitivity study was performed with the whole lower limb model to investigate the effect of the hip joint angle to acetabulum injury tolerance in frontal impacts.

RESULTS

The whole lower limb model proved to be stable under severe impacts along the knee, foot, and lateral components. In addition, the biomechanical and injury responses predicted by the model correlated well with the corresponding test data. An increase in hip joint extension angle from -30 to +20° relative to neutral posture showed an increase of 19 to 58 percent hip injury tolerance.

CONCLUSIONS

The stability and biofidelity response of the pelvis-lower limb (PLEX) model indicates its potential application in future frontal and lateral impact FE simulations.

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.

Statistical modeling of human liver incorporating the variations in shape, size, and material properties

The liver is one of the most frequently injured abdominal organs during motor vehicle crashes. Realistic numerical assessments of liver injury risk for the entire occupant population require incorporating inter-subject variations into numerical models. The main objective of this study was to quantify the shape variations of human liver in a seated posture and the statistical distributions of its material properties. Statistical shape analysis was applied to construct shape models of the livers of 15 adult human subjects, recorded in a typical seated (occupant) posture. The principal component analysis was then utilized to obtain the modes of variation, the mean model, and 95% statistical boundary shape models. In addition, a total of 52 tensile tests were performed on the parenchyma of three fresh human livers at four loading rates (0.01, 0.1, 1, and 10 s^-1) to characterize the rate-dependent and failure properties of the human liver. A FE-based optimization approach was employed to identify the material parameters of an Ogden material model for each specimen. The mean material parameters were then determined for each loading rate from the characteristic averages of the stress-strain curves, and a stochastic optimization approach was utilized to determine the standard deviations of the material parameters. Results showed that the first five modes of the human liver shape models account for more than 60% of the overall anatomical variations. The distributions of the material parameters combined with the mean and statistical boundary shape models could be used to develop probabilistic finite element (FE) models, which may help to better understand the variability in biomechanical responses and injuries to the abdominal organs under impact loading.

Biomechanical and Injury Response of Human Foot and Ankle Under Complex Loading

Ankle and subtalar joint injuries of vehicle front seat occupants are frequently recorded during frontal and offset vehicle crashes. A few injury criteria for foot and ankle were proposed in the past; however, they addressed only certain injury mechanisms or impact loadings. The main goal of this study was to investigate numerically the tolerance of foot and ankle under complex loading which may appear during automotive crashes. A previously developed and preliminarily validated foot and leg finite element (FE) model of a 50th percentile male was employed in this study. The model was further validated against postmortem human subjects (PMHS) data in various loading conditions that generates the bony fractures and ligament failures in ankle and subtalar regions observed in traffic accidents. Then, the foot and leg model were subjected to complex loading simulated as combinations of axial, dorsiflexion, and inversion loadings. An injury surface was fitted through the points corresponding to the parameters recorded at the time of failure in the FE simulations. The compelling injury predictions of the injury surface in two crash simulations may recommend its application for interpreting the test data recorded by anthropometric test devices (ATD) during crash tests. It is believed that the methodology presented in this study may be appropriate for the development of injury criteria under complex loadings corresponding to other body regions as well.

Statistical shape analysis of clavicular cortical bone with applications to the development of mean and boundary shape models

During car collisions, the shoulder belt exposes the occupant's clavicle to large loading conditions which often leads to a bone fracture. To better understand the geometric variability of clavicular cortical bone which may influence its injury tolerance, twenty human clavicles were evaluated using statistical shape analysis. The interior and exterior clavicular cortical bone surfaces were reconstructed from CT-scan images. Registration between one selected template and the remaining 19 clavicle models was conducted to remove translation and rotation differences. The correspondences of landmarks between the models were then established using coordinates and surface normals. Three registration methods were compared: the LM-ICP method; the global method; and the SHREC method. The LM-ICP registration method showed better performance than the global and SHREC registration methods, in terms of compactness, generalization, and specificity. The first four principal components obtained by using the LM-ICP registration method account for 61% and 67% of the overall anatomical variation for the exterior and interior cortical bone shapes, respectively. The length was found to be the most significant variation mode of the human clavicle. The mean and two boundary shape models were created using the four most significant principal components to investigate the size and shape variation of clavicular cortical bone. In the future, boundary shape models could be used to develop probabilistic finite element models which may help to better understand the variability in biomechanical responses and injuries to the clavicle.