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.     

Occurrence of selected pharmaceuticals in river water in Japan and assessment of their environmental risk

The existence of pharmaceuticals in the water environment is thought to be a potential problem for aquatic organisms. In this study, we conducted a nationwide survey to clarify the occurrence of 24 selected pharmaceuticals in major Japanese rivers and evaluated their environmental risk to aquatic organisms. We found a total of 22 substances in river waters at concentrations from several nanograms per liter to several micrograms per liter. We found the highest, which was 2.4 μg/L of caffeine, followed by 1.5 μg/L of crotamiton and 1.4 μg/L of sulpiride. We conducted an environmental risk assessment of the 22 pharmaceuticals detected in river water, for which predicted no-effect concentration (PNEC) values for crustacea and algae had been obtained. The measured environmental concentration/PNEC values of four substances, caffeine, carbamazepine, clarithromycin, and ketoprofen, exceeded 0.1 with the maximum value of 9.0 for clarithromycin. As clarithromycin exhibits a high environmental risk to aquatic organisms, particular attention is required.     

Integrated biogeochemical modelling of nitrogen load from anthropogenic and natural sources in Japan

This study proposed an integrated biogeochemical modelling of nitrogen loads from anthropogenic and natural sources in Japan. Firstly, the nitrogen load (NL) from different sources such as crop, livestock, industrial plant, urban and rural resident was calculated by using datasets of fertilizer utilization, population distribution, land use map, and social census. Then, the nitrate leaching from soil layers in farmland, grassland and natural conditions was calculated by using a terrestrial nitrogen cycle model (TNCM). The Total Runoff Integrating Pathways (TRIP) was used to transport nitrogen from natural and anthropogenic sources through river channels, as well as collect and route nitrogen to the river mouths. The forcing meteorological and hydrological data is a 30-year (1976–2005) dataset for Japan which were obtained by the land surface model, Minimal Advanced Treatments of Surface Interaction and Runoff (MATSIRO). For the model validation, we collected total nitrogen (TN) concentration data from 59 rivers in Japan. As a comparison result, calculated TN concentration values were in good agreement with the observed ones, which shows the reliability of the proposed model. Finally, the TN loads from point and nonpoint sources were summarized and evaluated for 59 river basins in Japan. The proposed modelling framework can be used as a tool for quantitative evaluation of nitrogen load in terrestrial ecosystems at a national scale. The connection to land use and climate data provides a possibility to use this model for analysis of climate change and land use change impacts on hydrology and water quality.     

Estimating monthly total nitrogen concentration in streams by using artificial neural network

Artificial Neural Network (ANN) is a flexible and popular tool for predicting the non-linear behavior in the environmental system. Here, the feed-forward ANN model was used to investigate the relationship among the land use, fertilizer, and hydrometerological conditions in 59 river basins over Japan and then applied to estimate the monthly river total nitrogen concentration (TNC). It was shown by the sensitivity analysis, that precipitation, temperature, river discharge, forest area and urban area have high relationships with TNC. The ANN structure having eight inputs and one hidden layer with seven nodes gives the best estimate of TNC. The 1:1 scatter plots of predicted versus measured TNC were closely aligned and provided coefficients of errors of 0.98 and 0.93 for ANNs calibration and validation, respectively. From the results obtained, the ANN model gave satisfactory predictions of stream TNC and appears to be a useful tool for prediction of TNC in Japanese streams. It indicates that the ANN model was able to provide accurate estimates of nitrogen concentration in streams. Its application to such environmental data will encourage further studies on prediction of stream TNC in ungauged rivers and provide a useful tool for water resource and environment managers to obtain a quick preliminary assessment of TNC variations.     

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.     

Polybrominated diphenyl ethers and persistent organochlorines in Japanese human adipose tissues

The present study determined concentrations of polybrominated diphenyl ethers (PBDEs) and persistent organochlorines (OCs) in Japanese human adipose tissues collected during 2003–2004. Concentrations of PBDEs in adipose tissues were 1–2 orders of magnitude lower than those of OCs. However, observed PBDE congener levels in this study were relatively higher than those in Japanese human adipose tissues collected during 2000 reported previously, while OC levels were comparable to those in specimens collected during 1999 reported by our group. In addition, no age-dependent accumulation of PBDEs was observed, while OC levels except chlordane compounds increased with age. These results indicate recent human exposure to PBDEs in Japan. Among PBDE congeners accumulated in Japanese adipose tissues, BDE-153 was dominant, but this trend was different from those in human milk (BDE-47) and blood (BDE-209) reported previously in Japan, implying the congener-specific kinetics in human bodies. The significant positive correlations between PBDEs and OCs were observed in Japanese adipose tissues, indicating the similar exposure route of these contaminants for Japanese citizens, probably via fish intake.     

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.