越野汽车驾驶室底板对爆炸冲击的响应研究
石秉良1,2 王显会3 盘朝奉1 张明3 张元正2 212013;2. 210028;3. 210094) (1. 江苏大学,镇江 陆军汽车试验场,南京 南京理工大学,南京
- ALE)【摘要】针对地雷等简易爆炸装置在车辆底部非接触爆炸问题,引入任意拉格朗日 欧拉( 算法,对某型越野汽车驾驶室底板结构的响应进行模拟。在整车环境中仿真分析了炸药当量、材料厚度、材料参数等引起的驾驶室底板结构响应的变化。研究结果表明:随材料厚度的减小和炸药当量的增加,驾驶室底板的应力和加速度响应增大;响应目标函数对底板材料的弹性模量和泊松比的灵敏度较高。同时,找出了响应敏感点,可为驾驶室底板结构抗爆炸优化设计提供参考。
1 前言
越野汽车由于具备优越的机动能力,在各国军队中得以广泛应用,随着战争形态的发展变化,其受到来自底部的爆炸袭击威胁日显突出。地雷等简易爆炸装置在驾驶室底部爆炸引起的损伤,是非对称战争中车辆及其乘员受到的主要损伤模式之一[1]。车辆应对爆炸冲击的能力越来越受到人们重视,驾驶室底板是车辆受到爆炸冲击时对乘员形成保护的第一道防线,对爆炸冲击波作用下越野汽车驾驶室底板结构的响应进行分析与研究是开展车辆底部防护能力研究的基础[2]。
对爆炸冲击作用下车辆驾驶室底板结构的塑性动力响应进行仿真分析,能获得较好的虚拟试验结果,可 1953为车辆驾驶室底部防护结构设计提供参考依据[3]。
Bruce
年, 等[4]模拟了一维气相不稳定径向和线形流,数
J. O. Hallquist 1976值仿真技术诞生; 博士[5]在 年主持并
LS-DYNA
完成了 软件的开发,较好地解决了非弹性材料结构在高速碰撞和爆炸冲击载荷下的塑性大变形动
K. Willianms力学响应的三维求解问题。 等[6]通过数字
7.5 kg C- 4
仿真对预埋的 地雷替代品在底部加装防护
M113 MTVL
组件的 底部爆炸时的作用效果进行了系统
LS- DYNA 9的仿真分析研究,利用 软件将车体分为 部
Westine分划分有限元单元格,并用 经验公式计算节点
Craig Barker Westine
速度。 等[7]同样利用 爆炸冲击经验公式对车辆底部受到爆炸冲击时车辆底板的速度响应进行了分析,为进一步研究直接与底板接触的乘员下肢
主题词:越野汽车 驾驶室底板 爆炸冲击 结构响应U463.83 A中图分类号: 文献标识码:
10.19620/j.cnki.1000-3703.20170280 DOI:
运动提供了理论依据。此方面的研究工作,近年来国内才得以重视,采用
- Arbitary Lagrange Euler,ALE)任意拉格朗日 欧拉( 算
LS-DYNA
法,利用 软件对爆炸冲击作用下典型越野汽车驾驶室底板结构响应的影响因素和不同当量炸药爆炸作用下的响应进行分析研究,可对车辆底部结构的防护能力进行预测和优化分析。
2 仿真计算理论模型 2.1 简化系统几何模型
车辆底板结构通常较为复杂,本文将车辆底部板结构简化为四边约束的标准靶板,地雷等简易爆炸装置的爆炸产物对车辆底部板壳结构的作用简化为爆炸冲击-
波对靶板的作用,作用系统简化为炸药 靶板系统几何1
模型[8],如图 所示。 - 500 mm炸药 靶板系统几何模型中,靶板为边长 、
2 mm Q235
厚 的方形结构,靶板材料选择 号钢, =7.8 g/cm3 E= 210 GPa,
密度ρ ,弹性模量 切线模量Et= 466 MPa, v= 0.30, σs= 235 MPa;泊松比 屈服应力 炸药
TNT 700 g, 1∶1,为裸装圆柱形 药柱,药量为 长径比为采用中心引爆的方式。靶板采用四边约束的方式固定
470 mm
在距炸药中心点 的位置,中心正对炸药。
-
炸药 靶板系统几何模型中,炸药爆轰采用高能燃
JWL[烧模型,产物状态方程采用 9-10]方程:
式中, A、B、R1、R2、ω为材料参数; p为高能炸药材料内压力; V为相对体积; E0为初始比内能。炸药的材料参数
1
和状态方程参数如表 所示。
2.2 ALE算法炸药-靶板系统基本模型
在仿真计算中,利用法向约束的方法在模型的对称
1/4 2面中设置对称约束,建立 有限元模型,如图 所示。
ALE
算法需要模拟爆炸冲击波,计算时将流固耦合方法用于冲击波对靶板作用的响应[11- 14],为了提高计算效率,模型中仅建了局部的空气域,并在空气网格表面设 置无反射边界条件以模拟无限空气域。炸药物质与空气介质为欧拉单元,靶板为拉格朗日单元,且炸药与空气网格的节点融合,空气单元覆盖靶板,炸药的初始引
2 -爆点设在炸药中心位置处。如图 所示,炸药 靶板爆
309 796,炸仿真有限元模型网格单元数为 节点数为
317 542
。
2.3 地雷-空气-车辆系统仿真模型
为较准确地分析爆炸冲击作用下某型越野汽车驾驶室底板的响应,应在整车环境中进行研究,本文建立
- - 3
了地雷 空气 车辆系统仿真模型,如图 所示。
CAD
建模时将整车三维建模软件中的 模型导入有Hypermesh
限元前处理软件 中。车身结构件模型中包含了大量倒角与内、外饰安装孔,由于显式计算中网格尺寸直接影响了计算的时间步长,为了避免一些不重要的细微几何信息造成网格划分困难,导致计算无法收CAD
敛,在进行网格划分之前,对 模型进行一定的简化
4与几何清理,清理后驾驶室底板模型如图 所示。
3 驾驶室底板响应研究
炸药在车辆底部爆炸时,驾驶室底板是车辆底部直接对乘员形成保护的关键部件,乘员座椅通常直接安装在底板上,脚部也通常直接放置在底板上。爆炸冲击的巨大能量通过驾驶室底板、座椅等间接作用于车内乘员,从而造成伤害。不同炸药爆炸后,驾驶室底板的冯
von Mises Stress) 5米斯应力( 云图和位移云图如图 所示。
5 326 g
从图 中可以看出,在规格为 地雷替代品的爆炸冲击波作用下,驾驶室底板已经出现了较多的塑性变形区域。从应力云图可知:驾驶室底板在结构的边缘、
450 MPa接合处、加强筋布置处均出现了超过 的应力区域;随着当量的增加,驾驶室底板的变形增加明显,在规
755 g
格为 地雷替代品的爆炸冲击波作用下,驾驶室底部中心变形最大,底板在爆炸冲击波作用下最大变形挠
90 mm
度接近 。但驾驶室没有出现结构贯穿的损伤。为进一步分析爆炸冲击对驾驶室底板的影响,在驾
5 6驶室底板上设置了 个测量点,如图 所示。仿真过程
1235
中采集位置 、、、的加速度时间历程和应力时间历程
755 g 2进行分析。 炸药爆炸时反应较为剧烈的位置 的
7, 8加速度时间历程见图 各测量点应力时间历程见图 。 4
通过对驾驶室底板的 个测量点的加速度和应力时间历程进行分析可知,地雷爆炸冲击时驾驶室底板中1)
心(位置 具有最大的应力峰值,此处容易产生结构变2) 70
形;驾驶室中间后端(位置 处,出现了较大的加速度响应,且该位置在爆炸冲击作用下加速度具有明显跃升趋势,此处容易给车内乘员相关部位直接或间接造成较大的冲击。
4 爆炸冲击对驾驶室底板影响分析
根据初步分析,对车辆驾驶室底板受到爆炸冲击后的响应影响较大的因素主要有爆炸物当量、驾驶室底板厚度和底板材料。本文主要从不同炸药当量、不同材料厚度方面分析爆炸冲击对驾驶室底板影响,同时对驾驶室底板响应特性对不同材料特性参数的灵敏
480 948)度进行分析。在驾驶室底板模型上选取A( 、
483 538) 487 878) 498 628) 4
B( 、C( 、D( 等 个测量点,如
9
图 所示。仿真分析中,重点分析驾驶室底板整体应力情况和不同测点的变形情况。
4.1 爆炸物当量
为了分析爆炸物当量对驾驶室底板响应的影响,
5 mm 552 g研究中选择了材料厚度为 的钢板,分别以 、
998 g 1 500 g 2 029 g 4
、 、 等 种不同当量的炸药进行仿
10真。不同当量炸药爆炸时驾驶室底板应力情况如图
4 11
所示;个测量点变形情况如图 所示。 10 998 g
从图 可以看出,当炸药当量超过 时,驾驶室底板边缘处已经出现了局部撕裂,驾驶室中后部也出
5 mm现了贯穿性损伤,这说明采用 钢板的该型驾驶室底板,在不采取任何底部防护措施的条件下,将无法抵
998 g
御超过 当量炸药的爆炸冲击。
4 4
通过分析 种不同炸药当量条件下 个测量点的位移变化情况进行比较可知:不同当量炸药爆炸时,驾驶室底板关键测量点的变形时间历程曲线的形态相似,随着地雷当量的增加,测量点位移幅度逐渐变大,
9 ms例如测量点A在不同炸药条件下均是在 时刻达到最大变形。
4
对 个测量点在不同炸药条件下的变形最大幅
12 552 g 998 g
值进行统计分析,结果如图 所示。在 、、
1 500 g
当量炸药的爆炸作用下,测量点变形幅值变化
1 500 g趋势基本相同(斜率近似)。但炸药当量从 增加
2 029 g
到 时,测量点的位移幅值出现了两种情况:测量点C的位移变化斜率突然增加,说明该点的变形幅值变化对炸药当量较为敏感,在车辆底部防爆炸冲击结构优化设计中,可考虑对此处的结构进行适当加强;测量点 A和B变形变化斜率随着炸药当量的增加而减小,说明测量点A和B的位移幅值变化对炸药当量变化敏感度较小,在车辆底部防爆炸冲击结构优化设计中,可考虑适当减小此处材料厚度。
4.2 材料厚度
为了研究驾驶室底板材料厚度对其底部防护能力
998 g
的影响,以 当量炸药爆炸冲击作为仿真输入,分析
1 mm 3 mm 5 mm 7 mm材料厚度分别为 、 、 、 时驾驶室底板
4
的应力与变形情况。 种不同厚度时驾驶室底板的冯
13
米斯应力云图如图 所示。从仿真分析结果可知,随着底板材料厚度不断增加,底板产生的最大应力不断减
4小,但图中这一现象并不明显,这是因为在 种不同厚度的爆炸分析中底板边缘的约束部分均出现了较大的
998 g应力集中,而驾驶室底板中心位置在 当量炸药的爆炸冲击下均没有出现贯穿损伤,也没有出现超过强度极限的位置。 1 mm 3 mm 7 mm材料厚度分别为 、 、 时驾驶室底板
14 4变形测量点的位移时间历程曲线如图 所示。 个测
998 g
量点在 当量炸药的爆炸冲击作用下,变形结果有较大区别。其中,测量点B由于更靠近爆炸中心,变形最大,随着底板材料厚度不断增加,测量点的变形减小,节点位移历程曲线变化趋势符合预期。
4
将不同材料厚度情况下 个测量点的位移峰值
15 1 mm进行统计,如图 所示,对于测量点B、C、D,在
3 mm
与 的不同料厚情况下,其位移变化衰减明显。其中,测量点C的位移变化对料厚最为敏感,但不同材料厚度条件下的测量点A处的位移变化不明显,这说明在后期优化设计中,在测量点BCD、、 对应的位置应适当增加材料厚度,才能有效提升驾驶室底部防护能力,而对于测量点A对应的位置可适当考虑相应的轻量化措施。
4.3 材料特性参数
材料特性参数影响着材料的各项性能,同样对材料的抗爆炸冲击能力有着较大的影响,但各种参数影响程度不同。研究材料不同特性参数变化对驾驶室底板对
998 g炸药爆炸冲击波响应的影响时,选择当量为 的炸药,将驾驶室底板材料的弹性模量E、切线模量Et、泊松比μ及屈服强度σ作为设计变量,以驾驶室底板的变形响应作为目标函数,进行参数灵敏度分析。驾驶室底板
2材料参数初始值及变化范围如表 所示。
由于整车爆炸仿真试验计算量大、耗时较长,仿真试验设计时应尽量缩减试验中水平值,从而减少试验设计中爆炸仿真的计算成本。本文使用多变量二水平筛
Plackett- Burman
选的 试验建立参数筛选样本空间,各
1.5变量的水平差值不能过大(高水平为低水平的 倍以
Plackett-Burman内),不考虑因子之间的交互效应,通过
Analysis of Variance,ANOVA)采样结合方差分析( 能高效地对多因子空间的结构变量进行一阶灵敏度分析,从而有效地判断各个变量对响应函数的贡献,缩减变量的
4
样本空间[15]。对车辆驾驶室底板的 个特性参数进行变
16 16形响应的全局灵敏度分析,结果如图 所示。由图
4
可知,个变量中泊松比μ对响应目标函数结果影响最大,其次是弹性模量E。
5 结束语
-
采用任意拉格朗日 欧拉算法,对某型车辆驾驶室底板抗爆炸冲击能力进行了仿真分析与研究,对炸药当量、底板材料厚度对驾驶室底板抗爆炸冲击能力影响规律进行了研究,分析了材料特性参数对材料变形目标函数的灵敏度,主要结论有:
a.
随着炸药当量的增加和材料厚度的减小,该车辆
988 g驾驶室底板的破坏程度均会加大,当炸药当量为 时,
5 mm该车辆驾驶室底板即便采用 的钢板也会被击穿;
b.
在材料厚度不变的情况下,测量点AB、的位移幅值变化对炸药当量变化敏感度较小,在车辆底部防爆炸冲击结构优化设计中,可考虑适当减小此处材料厚度,同时应适当增加CD、点处厚度。
c.
在炸药当量不变的情况下,测量点BCD、、 对应的位置应适当增加材料厚度,才能有效提升驾驶室底部防护能力,而对于测量点A对应的位置可适当考虑相应的轻量化措施。
d.
弹性模量E、切线模量Et、泊松比μ及屈服强度
4 σ等 个材料特性参数中, Eμ、对变形响应的目标函数灵敏度较高,车辆驾驶室底板设计时应合理选择材料。
参考文献
[1] Grujicic M, Bell W C. A Computational Analysis of Survivability of a Pick- up Truck Subjected to Mine Detonation Loads[J]. Multidiscipline Modeling in Materials & Structures, 2011, 7(4): 386-423.
[2] , , , .
石秉良 王显会 张云 等 军用车辆底部防护发展与研究
[J]. , 2016, 37(10):1902-1914.
综述 兵工学报
[3] Müller M, Dierkes U, Hampel J. Blast Protection in Military Land Vehicle Programmes: Approach, Methodology and Testing[C]// Structures under Shock and Impact. 2012: 247257.
[4] Bruce G H, Peaceman D W, Rachford H H, et al. Calculation of Unsteady- State Gas Flow Through Porous Media[J]. Journal of Petroleum Technology, 1953, 5(3): 7992.
[5] , , . LS- DYNA HyperWorks胡远志 曾必强 谢书港 基于 和
[M]. , 2011:的汽车安全仿真分析 北京:清华大学出版社
11-17.
[6] Willianms K, Poon K. A Numerical Analysis of the Effect of Surrogate Anti- tank Mine Blasts on the M133,DREV TM2000-007[R]. Canada: National Defense, 2004: 1-19.
[7] Barker C, Howle D. Predicting Floor Velocity of a Vehicle Subjected to an Underbody Blast Event[C]// Craig B, et al. Eightyth Shock & Vibration Symposium, Shock and Vibration Information Analysis Center (SAVIAC), San Diego, USA, 2009: 719-729.
[8] , , , . SPH
石秉良 周孔亢 张云 等 基于 算法的驾驶室底部结
[J]. , 2016,构对爆炸冲击波响应数值仿真 机械工程学报
52(16):132-139.
[9] . SPH [D].
陈娟 水下爆炸冲击载荷的 算法研究 哈尔滨:哈
, 2010.
尔滨工程大学
[10] , , .
徐定海 王善 杨世全 板壳结构接触爆炸数值仿真分析
[J]. , 2006, 27(1):53-56.
哈尔滨工程大学学报
[11] Genevieve T. Finite Element Simulation Using SPH Particles as Loading on Typical Light Armoured Vehicles [C]. 10th International LS-DYNA Users Conference.
[12] Clough R W. Early History of the Finite Element Method from the View Point of a Pioneer[J]. International Journal for Numerical Methods in Engineering, 2004, 60(60): 283287.
[13] . ANSYS11.0LS- DYNA [M].李裕春 基础理论与工程实践
, 2008:453-458.北京:水利水电出版社
[14] , .
王芳 冯顺山 爆炸冲击波作用下靶板的塑性大变形响
[J]. , 2003, 13(3):59-61.
应研究 中国安全学报
[15] Chen F, Zhang Q, Fei S, et al. Optimization of Ultrasonic Circulating Extraction of Samara Oil from Acer Saccharum Using Combination of Plackett- Burman Design and BoxBehnken Design[J]. Ultrasonics Sonochemistry, 2017, 35 (Pt A):161-175.
(责任编辑 斛畔) 2017 12 18修改稿收到日期为 年 月 日。