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.     

Effect of Aging on Brain Injury Prediction in Rotational Head Trauma—A Parameter Study with a Rat Finite Element Model

Objective: The aim of this study was to investigate the possible effects of age-related intracranial changes on the potential outcome of diffuse axonal injuries and acute subdural hematoma under rotational head loading.

Methods: A simulation-based parametric study was conducted using an updated and validated finite element model of a rat head. The validation included a comparison of predicted brain cortex sliding with respect to the skull. Further, model material properties were modified to account for aging; predicted tissue strains were compared with experimental data in which groups of rats in 2 different lifecycle stages, young adult and mature adult, were subjected to rotational trauma. For the parameter study, 2 age-dependent factors—brain volume and region-specific brain material properties—were implemented into the model. The models young adult and old age were subjected to several injurious and subinjurious sagittal plane rotational acceleration levels.

Results: Sequential analysis of the simulated trauma progression indicates that an increase in acute subdural hematoma injury risk indicator occurs at an early stage of the trauma, whereas an increase in diffuse axonal injury risk indicators occurs at a later stage. Tissue stiffening from young adult to mature adult rats produced an increase in strain-based thresholds accompanied by a wider spread of strain distribution toward the rear part of the brain, consistent with rotational trauma experiments with young adult and mature adult rats. Young adult to old age brain tissue softening and brain atrophy resulted in an increase in diffuse axonal injuries and acute subdural hematoma injury risk indicators, respectively.

Conclusions: The findings presented in this study suggest that age-specific injury thresholds should be developed to enable the development of superior restraint systems for the elderly. The findings also motivate other further studies on age-dependency of head trauma.     

三维点云基础上排土场垮塌影响范围研究

高台阶排土场下游常分布有居民区、厂房等,一旦垮塌,后果严重,因此完善的防治与应急准备工作十分关键,而排土场垮塌影响范围的预测是这些工作的重要前提依据。为提升高台阶排土场垮塌影响范围的预测准确性,需要获取准确的地形边界基础数据。通过利用三维激光扫描仪对新建排土场及其下游地形精细扫描,对获取的三维点云数据去噪、拼接、建模等处理,生成精度更高的数字高程模型与三维实体模型。在此基础上,结合新建高台阶排土场地质勘查与设计资料,建立颗粒离散元数值模型,利用PFC3D计算高台阶排土场垮塌影响范围,并与经验公式法计算结果对比分析,结果表明排土场颗粒离散元数值模型精度符合要求,可快速、准确地计算影响范围,从而为相关部门提供科学决策依据。     

Application of CFD on the sensitivity analyses of some parameters of the modified Hartmann tube

The aim of this work is to determine the influence of operating parameters such as the dispersion pressure, the ignition delay and height on the dust flammability. A Computational Fluid Dynamics (CFD) simulation, based on an Euler–Lagrange approach, was developed with Ansys Fluent™ and validated experimentally. Such analysis will facilitate the choice of the most conservative conditions for a flammability test. This paper is focused on a case study performed on wheat starch with the modified Hartmann tube. The dispersion process of the powder was studied with granulometric analyses performed in situ and high speed videos. Tests were performed with injections at gas pressure ranging from 3 to 6 bars and the evolution of the particle size distribution (PSD) was recorded at different ignition heights (5, 10 and 15 cm over the dispersion nozzle). The observations highlighted the presence of agglomeration/deagglomeration processes and dust segregation. Besides, a CFD simulation analysis was aimed at evaluating the impact of a set of parameters on the PSD and the local turbulence, which are closely linked to some flammability parameters. For this computational analysis, the CFD simulation was coupled with a collision treatment based on a Discrete Element Method (DEM) in order to consider the cohesive behavior of the combustible dust. Thus the results suggest performing the injection of the gases at approximately 5 bars for the flammability tests of wheat starch in order to obtain the finest PSD at a given ignition height. It is also shown that the finest PSD are obtained at 5 cm over the dispersion nozzle. However, the local instabilities and turbulence levels are so high during the first stages of the dispersion that the flame growth can be disturbed for short ignition delays. Moreover, the stabilization of the bulk of the dust cloud requires longer periods of time when the ignition sources are located at 15 cm. As a result, the recommended height to perform a flammability test is 10 cm in this case. Finally, this study proposes some tools that might improve the procedure of dust flammability testing.     

Composition and source apportionment of dust fall around a natural lake

The aim of this study was to determine the source apportionment of dust fall around Lake Chini, Malaysia. Samples were collected monthly between December 2012 and March2013 at seven sampling stations located around Lake Chini. The samples were filtered to separate the dissolved and undissolved solids. The ionic compositions(NO-3, SO2-4, Cl-and NH+4) were determined using ion chromatography(IC) while major elements(K, Na, Ca and Mg) and trace metals(Zn, Fe, Al, Ni, Mn, Cr, Pb and Cd) were determined using inductively coupled plasma mass spectrometry(ICP-MS). The results showed that the average concentration of total solids around Lake Chini was 93.49 ± 16.16 mg/(m2·day). SO2-4, Na and Zn dominated the dissolved portion of the dust fall. The enrichment factors(EF) revealed that the source of the trace metals and major elements in the rain water was anthropogenic, except for Fe. Hierarchical agglomerative cluster analysis(HACA) classified the seven monitoring stations and 16 variables into five groups and three groups respectively. A coupled receptor model, principal component analysis multiple linear regression(PCA-MLR), revealed that the sources of dust fall in Lake Chini were dominated by agricultural and biomass burning(42%),followed by the earth's crust(28%), sea spray(16%) and a mixture of soil dust and vehicle emissions(14%).     

超高速瞬态测试系统软硬件构架设计研究

目的研究采样频率为2~10 MHz的64通道超高速同步瞬态测试系统的设计技术,实现两类典型超高速瞬态测试系统的硬件架构设计与软件架构设计。方法一类采用PXI-Express高带宽总线和高速RAID磁盘阵列架构构建持续流盘存储的连续高速测试系统,另一类是采用大容量板载数据缓存和PXI总线事后下载传输数据的架构构造高速测试系统。在高性能测试软件设计方面,主要应用生产者/消费者结构与有限状态机相结合的软件架构进行高性能测试系统软件设计。结果目前64通道下基于持续流盘架构的测试系统受数据记录的速度限制系统最高采样频率仅达2.5 MHz,而基于板载缓存数据与PXI总线事后下载数据架构的测试系统最高采样频率可达10 MHz,测试时长可达5 s。结论当前两类架构的测试系统均可满足超高速瞬态测试需求,设计时需根据需求的最高采样频率决定使用的架构形式。     

深埋水平岩层隧道开挖稳定性分析及控制*

针对深埋高地应力水平岩层掌子面开挖稳定性及支护结构失效问题,以大峡谷隧道为工程背景,通过现场测试、室内试验、数值模拟等方法,探究深埋高地应力水平岩层失稳机理及控制措施。研究结果表明:坚硬岩体被节理面切割后,在高地应力作用下容易发生挤压破碎,破碎岩体遇水发生软化,导致掌子面发生大范围塌方,初支和超前支护失效;隧道开挖后岩层发生不均匀沉降,浅部岩层最先发生弯折破坏,层内块体错动滑移,继而向上方岩层发展,并伴随层间分离和层内裂隙发育,最终形成宏观破裂面;提出的台阶法、2 m开挖进尺、砼喷层、双层小导管、提高初支强度的整体优化控制措施,可有效提高现场支护效果。     

焙烧炉筒体结构强度有限元分析

对焙烧炉筒体结构进行有限元建模仿真,考虑热应力、压力和重力对筒体结构受力的影响,计算了焙烧炉关键部位的蠕变疲劳寿命,同时探究进气温度、保温层及对流系数等参数的影响。结果表明:热应力是影响结构强度的关键因素;焙烧炉筒体在许用蠕变变形量为5%的情况下,其使用寿命超过设计寿命;在筒体加热过渡段铺设保温层,降低与外界空气对流和提高进气温度,有利于降低筒体应力。     

车载接地装置的冲击接地阻抗研究

为解决车载装置在不同路面的接地问题,我们研制了可迅速安装和撤收的便携式接地装置,并选取一种组合在不同岩土介质地面上进行了实验,测得各自的冲击接地阻抗。然后,采用时域有限差分法对实验装置进行了数值分析,得出了该接地装置冲击接地阻抗的时域特性。与实验结果进行对比的结果表明,两者具有较好的一致性。     

关于区域可持续发展的系统分析

本文从更高、更新、更综合的高度和角度,对区域社会可持续发展系统进行了单元要素分析、结构层次分析。提出了三单元要素、三结构层次、三种发展类型的设想。这种分析尽管还是初步的、定性的、不很深刻的,但是为研究可持续发展系统提供了一种新的思路。