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.     

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.     

The association between prenatal exposure to polycyclic aromatic hydrocarbons and childhood intelligence: a systematic review of observational studies

Environmental Science and Pollution Research - Exposure to polycyclic aromatic hydrocarbons (PAHs) during pregnancy has been associated with many adverse child health. However, the evidence on such...     

Simulation and mathematical programming decision-making support for smallholder farming

Many mathematical programs have been developed over the past 50 years to aid agricultural experts and other farming decision-makers. The application of these mathematical programs has seen limited success because their development has focused on mathematical theory as opposed to the requirements needed for application. This paper describes the development of two mathematical programs that were designed to integrate with a visualization simulation that aids a nontraditional group of agricultural decision-makers: illiterate Sri Lankan subsistence farmers. The simulation was designed to help these illiterate farmers make business decisions about their crop selection choices which, in turn, will help them develop their business plans required for obtaining bank micro-loans. This paper’s focus is on the use of linear programming as a potential tool to demonstrate the benefits of crop diversification and rotation to the farmer based on various available crop types. It also highlights the issues using such an approach.     

Stochastic growth reduces population fluctuations in Daphnia-algal systems

Deterministic, size-structured models are widely used to describe consumer-resource interactions. Such models typically ignore potentially large random variability in juvenile development rates. We present simple representations of this variability and show five approaches to calculating the model parameters for Daphnia pulex interacting with its algal food. Using our parameterized models of growth variability, we investigate the robustness of a recently proposed stabilizing mechanism for Daphnia populations. Growth rate variability increases the range of enrichments over which small-amplitude cycles or quasi-cycles occur, thus increasing the plausibility that the underlying mechanism contributes to the prevalence of small-amplitude cycles in the field and in experiments. More generally, our approach allows us to relate commonly available information on variance of development times to population stability.     

The Effect of Estrogens,River Water,and Treated Wastewater on Vitellogenin Induction of Japanese Medaka1

Abstract: The objective of this study is to evaluate the effect of estrogens on fish endocrine disruption in river water and treated wastewater. Endogenous estrogen estrone (E1), 17β‐estradiol (E2), river water, and treated wastewater were used for exposure tests on male Japanese medaka. Vitellogenin induction in male medaka was regarded as the endpoint of endocrine disruption. The effective concentrations of E1 and E2 on vitellogenin induction in medaka were evaluated by breeding medaka for 14 days in tanks with various concentrations of E1 or E2. Vitellogenin induction increased with elapse of time during exposure, with higher estrogen concentrations causing greater vitellogenin induction. According to the test results, the lowest observed effect concentrations (LOECs) of E1 and E2 were estimated to be 31.6 and 5.0 ng/l, respectively; and the E2 equivalent (EEQ) LOEC value of E1 was 5.2 ng‐E2/l, derived by multiplying the relative potency of E1 by that of E2 as estimated in this study. In continuous exposure tests using river water or treated wastewater, the estrogenic activity, an index of total estrogenic potential measured by yeast screen assay and expressed as EEQ, varied widely during the exposure tests, and significant vitellogenin induction was observed after several days of high levels of estrogenic activity in water. Vitellogenin concentration tended to significantly increase if EEQ exceeded the level of 5 ng‐E2/l. The threshold value was substantially the same as the results for E1 and E2 exposure tests of medaka. Consequently, EEQ over 5 ng‐E2/l was revealed to have the potential to cause endocrine disruption of male medaka. As estrogenic activity exceeding 5 ng‐E2/l was observed in some rivers in Japan, total estrogenic potential needs to be further decreased at wastewater treatment plants to prevent fish endocrine disruption.     

Fate estimation of estrogenic substances in an urban river by flux calculation.

This study investigated the fate of estrogenic substances in an urban river receiving discharge from wastewater treatment plants (WTPs) by flux calculation, focusing on the middle reaches of the Tama River in Tokyo, which is one of the most urbanized rivers in Japan. The level of estrogenic activity flux was almost negligibly small at the upstream station. The level was considerably raised after inflows from the WTPs and then the level declined in the lower reaches of the river. When contributions of estrogenic substances to estrogenic activity were estimated, estrone (E1) was the primary contributor to the total estrogenic activity in all the sampling stations, followed by estradiol (E2). The contribution of nonylphenol to estrogenic activity was small. The E1 and E2 accounted for approximately 90% or more of estrogenic activity in the Tama River. As for the total fluxes of the estrogenic substances in the study area in the Tama River, the proportion of flux associated with WTP discharge was approximately 100% of the total fluxes, and the effects of the tributaries flowing into the river were almost negligible. When the reduction ratios of estrogenic activity were calculated by the flux, the ratios were found to increase toward the lower reaches of the river. Similar changes were observed for E1. Meanwhile, the change of reduction ratios for E2 was different from that observed for estrogenic activity.     

Synergistic effect on the performance of ash-based bricks with glass wastes and granite tailings along with strength prediction by adopting machine learning approach

Environmental Science and Pollution Research - The study proposes a novel and sustainable method to appropriately utilize wastes from granite as well as glass industries in brick manufacturing. An...