A high performance crashworthiness simulation system based on GPU

Crashworthiness simulation system is one of the key computer-aided engineering (CAE) tools for the automobile industry and implies two potential conflicting requirements: accuracy and efficiency. A parallel crashworthiness simulation system based on graphics processing unit (GPU) architecture and the explicit finite element (FE) method is developed in this work. Implementation details with compute unified device architecture (CUDA) are considered. The entire parallel simulation system involves a parallel hierarchy-territory contact-searching algorithm (HITA) and a parallel penalty contact force calculation algorithm. Three basic GPU-based parallel strategies are suggested to meet the natural parallelism of the explicit FE algorithm. Two free GPU-based numerical calculation libraries, cuBLAS and Thrust, are introduced to decrease the difficulty of programming. Furthermore, a mixed array and a thread map to element strategy are proposed to improve the performance of the test pairs searching. The outer loop of the nested loop through the mixed array is unrolled to realize parallel searching. An efficient storage strategy based on data sorting is presented to realize data transfer between different hierarchies with coalesced access during the contact pairs searching. A thread map to element pattern is implemented to calculate the penetrations and the penetration forces; a double float atomic operation is used to scatter contact forces. The simulation results of the three different models based on the Intel Core i7-930 and the NVIDIA GeForce GTX 580 demonstrate the precision and efficiency of this developed parallel crashworthiness simulation system.     

Homogenization with uncertainty in Poisson ratio for polymers with rubber particles

The main aim of this work is dual computer analysis of probabilistic coefficients for the homogenized tensor of the polymer filled with the rubber particles having randomized Poisson ratios of both constituents. The major issue is to verify an influence of a randomness in rubber Poisson ratio close to the compressibility limit on the uncertainty of the effective tensor probabilistic characteristics. Probabilistic analysis presented here is carried out using mainly the stochastic perturbation technique provided by the common application of the traditional FEM commercial code ABAQUS and the symbolic computations package MAPLE. This FEM-based technique employs polynomial response function of the optimum order recovered from the weighted least squares method and following a set of deterministic solutions obtained for various values of the randomized input parameter. Optimization procedure is released entirely into a symbolic environment, where maximization of the correlation factor together with minimization of the fitting variance and approximation error are applied. Homogenization technique consists in equating of deformation energies for the real composite and the artificial one characterized by the effective elasticity tensor with uncertainty.     

Design and Fabrication of Concentration‐Gradient Generators with Two and Three Inlets in Microfluidic Chips

A simple and low‐cost method for designing and fabricating concentration‐gradient generators with two and three inlets is proposed which can generate different concentration gradients at varying flow velocities. The microchannel structure was designed in S‐shape and left‐right symmetry. The concentration‐gradient generator was simulated based on the finite element method. The microchannels were processed on a computer numerical control (CNC) engraving and milling machine on poly(methylmethacrylate) substrate, and then two concentration‐gradient generators were fabricated by hot bonding technology. The results of experiment and simulation were compared to prove the feasibility of the method. Flow velocity was an important factor for generating different concentration gradients. The concentration‐gradient profiles of the generators with two and three inlets present approximately linear and quadratic curves.     

DP980+Z烘烤硬化值与退火温度的关系

摘要：为了研究退火温度对镀锌DP980+Z烘烤硬化值的影响，退火温度控制在760～820℃之间，系统分析退火温度对烘烤硬化值的影响。通过准静态拉伸试验机测量烘烤硬化值及抗拉强度，采用lepara试剂对组织中的马氏体进行着色，利用金相显微镜及图像处理软件测量马氏体的体积分数；采用扫描电镜观察DP980+Z的双相组织特点，并且将组织图片通过CAD转化成有限元图进行网格划分，建立代表性体积单元（RVE），通过有限元分析铁素体、马氏体强度对烘烤硬化值的影响。在同样的变形量情况下，DP980+Z的原始屈服强度越高，烘烤硬化值越高。     

Passive target detection and tracking from electromagnetic field measurements

This paper presents a novel approach to the localization of moving targets in a complex environment based on the measurement of the perturbations induced by the target presence on an independently‐generated time‐varying electromagnetic field. Field perturbations are measured via a set of sensors deployed over the domain of interest and used to detect and track a possible target by resorting to a particle Bernoulli filter (PBF). To comply with real‐time operation, the PBF works along with an artificial neural network (ANN) model of the environment trained offline via finite elements (FEs). The performance of the proposed algorithm is assessed via simulation experiments.     

Pull out of FRP reinforcement from masonry pillars: Experimental and numerical results

In this paper, debonding phenomena between carbon fiber reinforced polymer (CFRP) strips and masonry support were investigated on the basis of single-lap shear tests, considering different dimensions of the bond length. To capture the post-peak response of the CFRP–masonry joint, the slip between the support and the reinforcement strip was controlled using a clip gauge positioned at the end of the reinforcement. The tests were simulated by means of a finite element model able to capture the post-peak snap-back behavior due to the failure process. The numerical model is based on zero-thickness interface elements and on a proper non-linear cohesive law. The comparison between experimental and numerical results was performed in terms of overall response, measured by both the machine stroke and the clip gauge positioned at the free end of the reinforcement. The cases of effective bond length greater and lesser than the minimum anchorage length, suggested by the CNR Italian recommendation, were considered.     

A comparison of transformation methods for evaluating two‐dimensional weakly singular integrals

Accurate numerical evaluation of integrals arising in the boundary element method is fundamental to achieving useful results via this solution technique. In this paper, a number of techniques are considered to evaluate the weakly singular integrals which arise in the solution of Laplace's equation in three dimensions and Poisson's equation in two dimensions. Both are two‐dimensional weakly singular integrals and are evaluated using (in a product fashion) methods which have recently been used for evaluating one‐dimensional weakly singular integrals arising in the boundary element method. The methods used are based on various polynomial transformations of conventional Gaussian quadrature points where the transformation polynomial has zero Jacobian at the singular point. Methods which split the region of integration into sub‐regions are considered as well as non‐splitting methods. In particular, the newly introduced and highly accurate generalized composite subtraction of singularity and non‐linear transformation approach (GSSNT) is applied to various two‐dimensional weakly singular integrals. A study of the different methods reveals complex relationships between transformation orders, position of the singular point, integration kernel and basis function. It is concluded that the GSSNT method gives the best overall results for the two‐dimensional weakly singular integrals studied. Copyright © 2002 John Wiley & Sons, Ltd.     

水力射孔对套管强度的影响研究

水力射孔技术是一种新型的完井方式，深穿透水力射孔技术辅助定向水力压裂可以实现油气层改造和油气井增产。但是，水力射孔套管会在孔眼处造成应力集中，使得局部应力过大而削弱了套管和水泥环的强度，可能会影响套管的使用寿命。因此，文章利用有限元软件建立了射孔套管的三维有限元模型，分析了水力射孔后套管管体的应力分布，重点研究了水力射孔参数-孔密、孔径、孔深、射孔方位角等对套管强度的影响规律。分析结果表明，沿着最大水平地应力方向深穿透水力射孔，选择合适的孔密、孔径可将射孔对套管的损害降至最小。研究结果为提高水力射孔完井设计水平提供了一定的参数依据。     

燃气管网动态仿真的研究及应用

研究和了解燃气在管道中的流动，是进行燃气管网系统优化和提高管网输运系统安全性的前提和基础。以流体力学三大守恒方程为基础，建立了等温和非等温条件下的燃气管网稳态和动态仿真理论模型。以有限差分法为基础，得到了上述模型的求解方法，并加入了耗散项，提高了求解方法的稳定性。通过增加初始点和延长出口点的方法，简化了模型的求解。通过在北京六环天然气管网中的初步应用，分别得到了与实际值综合相对误差为1.21％和2.62％的压力和流量仿真值，从而验证了模型及求解方法的可靠性。     

大型弧形闸门静力特性有限元分析

以某水电站溢洪道弧形闸门为例,建立了大型弧形闸门的有限元分析模型,应用ANSYS软件对其在各种工况下的支铰反力、闸门应力、位移挠度以及屈曲失稳等静力特性进行了计算分析,为完善闸门结构的设计提供了依据,对同类闸门结构的设计计算有参考意义。