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的原始屈服强度越高，烘烤硬化值越高。     

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.     

构皮滩水电站泄洪洞优化设计

构皮滩水电站为满足2008年12月工程提前发电的需要,结合现有泄洪洞地质情况和施工情况对泄洪洞布置进行了专题研究,结果表明利用降低进口高程为550 m的泄洪洞和坝身放空底孔联合参与后期导流,可缓解2007年导流隧洞下闸封堵后大坝在2008年汛期施工的压力,为大坝施工增加了2个月左右的工期。     

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.     

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

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

Fatigue crack propagation of multiple coplanar cracks with the coupled extended finite element/fast marching method

A numerical technique for modeling fatigue crack propagation of multiple coplanar cracks is presented. The proposed method couples the extended finite element method (X-FEM) [Int. J. Numer. Meth. Engng. 48 (11) (2000) 1549] to the fast marching method (FMM) [Level Set Methods & Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science, Cambridge University Press, Cambridge, UK, 1999]. The entire crack geometry, including one or more cracks, is represented by a single signed distance (level set) function. Merging of distinct cracks is handled naturally by the FMM with no collision detection or mesh reconstruction required. The FMM in conjunction with the Paris crack growth law is used to advance the crack front. In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity [Comput. Meth. Appl. Mech. Engng. 139 (1996) 289]. This enables the domain to be modeled by a single fixed finite element mesh with no explicit meshing of the crack surfaces. In an earlier study [Engng. Fract. Mech. 70 (1) (2003) 29], the methodology, algorithm, and implementation for three-dimensional crack propagation of single cracks was introduced. In this paper, simulations for multiple planar cracks are presented, with crack merging and fatigue growth carried out without any user-intervention or remeshing.     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.