枪弹穿甲过程仿真的有限元接触模型与分区并行计算误差特性

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

清华大学辅仁网/2017-07-07

吕振华, 刘赛

清华大学汽车工程系, 北京 100084

Finite element contact modeling method and error characteristics of partitioned parallel computations for bullet penetration simulations

LÜ Zhenhua, LIU Sai

Department of Automotive Engineering, Tsinghua University, Beijing 100084, China

摘要:

输出: BibTeX | EndNote (RIS)

摘要为了提高小口径普通钢芯弹和穿甲燃烧弹侵彻均质钢靶板的仿真分析精度，分析了接触力罚函数算法的接触刚度参数和接触摩擦系数对弹道极限速度计算值的影响，分析了通过不同的并行计算系统得到的弹道极限速度计算值的差异及其统计规律。研究结果表明：多部件间的接触力罚函数的接触刚度参数分别设置与优化方法可以提高弹道极限速度值的计算精度；接触摩擦系数的不同设置值对穿甲燃烧弹侵彻较厚靶板的弹道极限速度计算值有显著影响，需探索多部件间的接触摩擦系数的准确定义方法；由于分区并行计算过程中多个区域的数据交换及参与计算的顺序的随机变化会导致计算结果的不一致，并行计算的弹道极限速度值存在离散性，需在同一入射速度下进行多次重复计算，并基于统计分析确定可靠的弹道极限速度值。

关键词 ：枪弹穿甲,有限元仿真分析,接触力罚函数算法,接触摩擦系数,弹道极限速度,分区并行计算误差

Abstract：Improved ballistic simulations were developed for homogeneous steel sheets against small-caliber bullets. The simulations analyzed the influences of the contact-force penalty function parameters and contact friction coefficients on the contact interface state and ballistic limit velocity. The inconsistencies among computed ballistic limit velocities predicted by different partitioned parallel systems were investigated statistically. The interface stiffnesses within the contact-force penalty function should be defined between different contact pairs rather than assigning an equal interface stiffness to all contact regions to improve the accuracy of the ballistic limit velocity. Different friction coefficients between contact pairs can seriously affect the computed ballistic limit velocity for thicker steel targets against armor piercing bullets. Further study is needed to determine the correct friction coefficients for multiple contact pairs. The exchange of regional data and the use of this data in different sequences make the partitioned parallel computations of repeated computation tasks inconsistent, which affects the ballistic limit velocities predicted by the parallel computations. Therefore, repeated calculations using the same model for the same impact velocity are needed to obtain reliable ballistic limit velocity statistics.

Key words：bullet penetration finite element analysis contact-force penalty method contact friction coefficient ballistic limit velocity partitioned parallel computation error

收稿日期: 2016-07-24 出版日期: 2017-05-20

ZTFLH:

O385

引用本文:

吕振华, 刘赛. 枪弹穿甲过程仿真的有限元接触模型与分区并行计算误差特性[J]. 清华大学学报（自然科学版）, 2017, 57(5): 483-490.
LÜ Zhenhua, LIU Sai. Finite element contact modeling method and error characteristics of partitioned parallel computations for bullet penetration simulations. Journal of Tsinghua University(Science and Technology), 2017, 57(5): 483-490.

链接本文:

http://jst.tsinghuajournals.com/CN/10.16511/j.cnki.qhdxxb.2017.22.025或 http://jst.tsinghuajournals.com/CN/Y2017/V57/I5/483

图表:

图1 接触力罚函数算法示意图

图2 小口径枪弹侵彻均质钢靶板的模型剖视图

图3 接触刚度系数对最大接触穿入功和最大穿入量的影响

表1 接触刚度系数的优化值

图4 接触摩擦系数对弹道极限速度计算值的影响

表2 摩擦系数值为0.05时接触刚度系数的优化值

表3 不同摩擦系数下的侵彻结果

图5 模型1的2核并行计算区域分割

表4 20次重复的2核并行计算得到的穿甲状态(模型1)

表5 不同数值实验次数下的穿透频率(模型1)

图6 2核并行计算的穿透概率分布(模型1)

图7 多核并行计算的穿透概率分布

表6 弹道极限速度计算值概率分布为5%~95%的速度区间

表7 采用不同计算系统求得的弹道极限速度统计值

图8 模型1的2核并行计算的非对称区域分割

图9 非对称区域分割的2核并行计算的穿透概率分布

参考文献:

[1]	Hallquist J O, Goudreau G L, Benson D J. Sliding interfaces with contact-impact in large-scale Lagrangian computations [J]. Computer Methods in Applied Mechanics and Engineering, 1985, 51(1-3): 107-137.
[2]	Peric D, Owen D R J. Computational model for 3-D contact problems with friction based on the penalty method [J]. International Journal for Numerical Methods in Engineering, 1992, 35(6): 1289-1309.
[3]	BLørvik T, Olovsson L, Dey S, et al. Normal and oblique impact of small arms bullets on AA6082-T4 aluminium protective plates [J]. International Journal of Impact Engineering, 2011, 38(7): 577-589.
[4]	Holmen J K, Johnsen J, Jupp S, et al. Effects of heat treatment on the ballistic properties of AA6070 aluminium alloy [J]. International Journal of Impact Engineering, 2013, 57: 119-133.
[5]	BLørvik T, Dey S, Olovsson L. Penetration of granular materials by small-arms bullets [J]. International Journal of Impact Engineering, 2015, 75: 123-139.
[6]	Jankowiak T, Rusinek A, Wood P. A numerical analysis of the dynamic behaviour of sheet steel perforated by a conical projectile under ballistic conditions [J]. Finite Elements in Analysis and Design, 2013, 65: 39-49.
[7]	Kpenyigba K M, Jankowiak T, Rusinek A, et al. Influence of projectile shape on dynamic behavior of steel sheet subjected to impact and perforation [J]. Thin-Walled Structures, 2013, 65: 93-104.
[8]	Arias A, Rodríguez-Martínez J A, Rusinek A. Numerical simulations of impact behaviour of thin steel plates subjected to cylindrical, conical and hemispherical non-deformable projectiles [J]. Engineering Fracture Mechanics, 2008, 75(6): 1635-1656.
[9]	Rusinek A, Rodríguez-Martínez J A, Arias A, et al. Influence of conical projectile diameter on perpendicular impact of thin steel plate [J]. Engineering Fracture Mechanics, 2008, 75(10): 2946-2967.
[10]	Farhat C. A simple and efficient automatic FEM domain decomposer [J]. Computers and Structures, 1988, 28(5): 579-602.
[11]	Pothen A, Simon H D, Liou K P. Partitioning sparse matrices with eigenvectors of graphs [J]. SIAM Journal on Matrix Analysis and Applications, 1990, 11(3): 430-452.
[12]	Hallquist J O. LS-DYNA Theory Manual[Z/OL]. [2016-06-21]. http://www.lstc.com/pdf/ls-dyna_theory_manual_2006.pdf.
[13]	Zukas J A, Hohler V, Jameson R L, et al. High Velocity Impact Dynamics [M]. New York: Wiley, 1990.
[14]	BLørvik T, Dey S, Clausen A H. Perforation resistance of five different high-strength steel plates subjected to small-arms projectiles [J]. International Journal of Impact Engineering, 2009, 36(7): 948-964.
[15]	Ravid M, Bodner S R. Dynamic perforation of viscoplastic plates by rigid projectiles [J]. International Journal of Engineering Science, 1983, 21(6): 577-591.
[16]	Ozawa K, Miyazaki M. A summation algorithm with error correction for parallel computers [J]. Systems and Computers in Japan, 1993, 24(7): 62-68.
[17]	GJB 4300—2002. 军用防弹衣安全技术性能要求 [S]. 北京: 中国人民解放军总后勤部, 2002.GJB 4300—2002. Requirements of Safety Technical Performance for Military Body Armor [S]. Beijing: General Logistics Department of CPLA, 2002. (in Chinese)

No related articles found!