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