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.     

面积坐标法构造含转角自由度的四结点膜元

以四边形面积坐标作为工具,构造了两个含转角自由度的广义协调四边形单元AQ4和lAQ4。它们通过强式分片检验,与同类单元相比,具有很高的计算精度,能消除梯形闭锁现象,有很强的抗网格畸变的能力。     

径向水平井弯曲转向机构影响因素仿真研究

针对径向水平井钻进中钻杆弯曲转向前进运动困难的问题,通过建立弯曲转向机构仿真有限元模型,对其主要影响因素:滑道轨迹曲率、间隙、工作压力、滚轮形状与位置和钻杆壁厚等用ANSYS软件进行仿真研究。研究表明:影响截面变形的主要因素是滑道曲率的改变和滑道摩擦因数,钻井失败时钻杆所受阻力的增大不是因为升高工作压力导致截面变形过大与滑道干涉所致,钻杆经校直段,钻头中轴线与校直中心线存在一定的角度,即钻杆的前进轨迹就会偏离理论值,或与地层干涉,或可能再次发生弯曲变形;另一主要原因是转向器滑道工作恶化。增加钻杆壁厚,降低工作压力和滑道摩擦因数,设计合理弯曲转向机构是解决问题的关键。     

Development of geometrically exact new shell elements based on general curvilinear co‐ordinates

In the present study first‐order shear deformable shell finite elements based on general curvilinear co‐ordinates are proposed. For the development of the present shell elements, a partial mixed variational functional with independently assumed strains is provided in order to avoid the severe locking troubles known as transverse shear and membrane lockings. Bubble functions are included in the shape function of displacement to improve the performance of the developed element. The proposed assumed strain four‐ and nine‐node elements based on the general tensor shell theory provide an efficient linkage framework for shell surface modelling and finite element analysis. In the several benchmark problems, the present shell elements with exact geometric representations demonstrate their performance compared to previously reported results. Copyright © 2002 John Wiley & Sons, Ltd.     

在经典极限下量子斯特林制冷机的优化

经济利润率是评价一个实际热力装置的主要指标之一。将有限时间热力学，非平衡量子统计理论和yong经济学相结合，导出了量子斯特林制冷机的最大利润率以及对应的性能界限，其结果与实际斯特林制冷机的优化设计和模型评估提供了一个最佳的预选方案。     

压型钢板组合楼板承载力的有限元分析

压型钢板组合楼板是 2 0世纪发展起来的一种新型结构形式。它具有显著的优势 ,如轻质、高强 ,施工快捷等 ,将在未来的建筑领域中大有作为。组合楼板的承载力主要取决于钢板 -混凝土界面处的抗剪强度。建立了一个二维有限元模型来研究单向组合板的性能。并将计算结果与试验数据进行对比 ,验证了模型的有效性     

HP700圆锥破碎机动锥的有限元分析

对HP700圆锥破碎机的主要部件动锥建立了有限元分析模型，用ANSYS软件对动锥进行了有限元强度和刚度分析，结果输出了动锥分布及变形图，对圆锥破碎机的结够设计提供了依据。     

Adhesive element modelling and weighted static shape control of composite plates with piezoelectric actuator patches

A novel finite element model is presented for static and dynamic analysis of composite plates integrated with a laminated piezoelectric layer, a host laminated composite plate and an adhesive layer between them. A new adhesive element is developed which includes both peel and shear effects in the adhesive layer based on first‐order shear deformation plate theory. The thin adhesive layer between the piezoelectric layer and the host plate is modelled by assuming that it carries constant shear and peel strains throughout its thickness. In addition, a weighted static shape control scheme for finding the optimal voltage distribution for static shape control is given. By selecting different weighting matrices, a variety of items such as displacements, slopes, curvatures, strains and even generalized forces, can be included in finding the optimal actuating voltage for static shape control. The present model is validated by comparing with those results available in the literature. The numerical results show that the weighted linear least method can give a satisfactory voltage distribution to best match the desired shape. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct stiffness analysis of a composite beam with partial interaction

The use of the conventional semi-analytical stiffness method in finite element analysis, in which interpolation polynomials are used to develop the stiffness relationships, leads to problems of curvature locking when beam-type elements are developed for composite members with partial interaction between the materials of which it is comprised. The curvature locking phenomenon that occurs for composite steel–concrete members is quite well reported, and the general approach to minimizing the undesirable ramifications of curvature locking has been to use higher-order polynomials with increasing numbers of internal nodes. This paper presents an alternate formulation based on a direct stiffness approach rather than starting from pre-defined interpolation polynomials, and which does not possess the undesirable locking characteristics. The formulation is based on a more general approach for a bi-material composite flexural member, whose constituent materials are joined by elastic shear connection so as to provide partial interaction. The stiffness relationships are derived, and these are applied to a simply supported and a continuous steel–concrete composite beam to demonstrate the efficacy of the method, and in particular its ability to model accurately both very flexible and very stiff shear connection that causes difficulties when implemented in competitive semi-analytical algorithms. Copyright © 2004 John Wiley & Sons, Ltd.