Please wait a minute...
 首页  期刊介绍 期刊订阅 联系我们 横山亮次奖 百年刊庆
 
最新录用  |  预出版  |  当期目录  |  过刊浏览  |  阅读排行  |  下载排行  |  引用排行  |  横山亮次奖  |  百年刊庆
清华大学学报(自然科学版)  2017, Vol. 57 Issue (5): 483-490    DOI: 10.16511/j.cnki.qhdxxb.2017.22.025
  汽车工程 本期目录 | 过刊浏览 | 高级检索 |
枪弹穿甲过程仿真的有限元接触模型与分区并行计算误差特性
吕振华, 刘赛
清华大学 汽车工程系, 北京 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
全文: PDF(1364 KB)  
输出: BibTeX | EndNote (RIS)      
摘要 为了提高小口径普通钢芯弹和穿甲燃烧弹侵彻均质钢靶板的仿真分析精度,分析了接触力罚函数算法的接触刚度参数和接触摩擦系数对弹道极限速度计算值的影响,分析了通过不同的并行计算系统得到的弹道极限速度计算值的差异及其统计规律。研究结果表明:多部件间的接触力罚函数的接触刚度参数分别设置与优化方法可以提高弹道极限速度值的计算精度;接触摩擦系数的不同设置值对穿甲燃烧弹侵彻较厚靶板的弹道极限速度计算值有显著影响,需探索多部件间的接触摩擦系数的准确定义方法;由于分区并行计算过程中多个区域的数据交换及参与计算的顺序的随机变化会导致计算结果的不一致,并行计算的弹道极限速度值存在离散性,需在同一入射速度下进行多次重复计算,并基于统计分析确定可靠的弹道极限速度值。
服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
关键词 枪弹穿甲有限元仿真分析接触力罚函数算法接触摩擦系数弹道极限速度分区并行计算误差    
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 wordsbullet penetration    finite element analysis    contact-force penalty method    contact friction coefficient    ballistic limit velocity    partitioned parallel computation error
收稿日期: 2016-07-24      出版日期: 2017-05-15
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!
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
版权所有 © 《清华大学学报(自然科学版)》编辑部
本系统由北京玛格泰克科技发展有限公司设计开发 技术支持:support@magtech.com.cn