Study on Deformation of Perforated Plates under Blast Loading
-
摘要: 利用非线性动力学软件AUTODYN,对爆炸冲击波作用下无孔及预制孔靶板的塑性变形进行了数值模拟,获得了不同孔数、孔径情况下靶板中心点挠度的变化规律;设计了模型实验,其结果与数值结果符合较好;并以爆炸冲击波作用下无孔靶板中心点挠度计算公式和数值计算结果为基础,通过拟合给出了预制孔靶板中心点挠度与孔密度、孔径之间的函数关系,该关系可为有孔平板目标的毁伤评估及破片、冲击波对平板目标联合毁伤研究提供参考。Abstract: Deformation of perforated plates under impulsive loading of blast wave was simulated by AUTODYN software and the deflection relationship of plates with different hole numbers and diameters was obtained.Experiments were designed and the consistency between model computation and experiment results was validated.A general perforated plate deflection formula was proposed based on deflection formula of non-hole thin plate and the computation results, which can provide reference for damage estimation of perforated plate and damage study on effects of fragment and blast wave.
-
Key words:
- blast wave /
- perforated plate /
- deflection /
- simulation
-
1. 引言
无论在军事还是民用中, 爆炸冲击波载荷作用下平板结构的动态响应一直是重要的工程研究课题。在对结构进行安全性防护设计时, 通常将目标等效为梁或者板等简单结构, 结构的塑性变形是判定结构破坏、失效的重要参考指标。近年来国内外学者对四边固定的平板在爆炸冲击波作用下的变形问题进行了大量研究, 获得了平板塑性变形相关参数的多种计算模型[1-4]。
Nurick和Martin以经典的冲击波场参数计算方法和冲击动力学理论为基础, 对方形平板塑性变形进行了理论研究和分析, 提出了冲击波作用下方形平板中心点挠度的计算公式[4]。其后, Chung Kim Yuen等人[5]对爆炸载荷下低碳钢方板的塑性动力学响应问题进行了大量的实验研究, 实验装药量范围从100到26 288 kg, 获得了不同炸药高度情况下的冲击波参数和方形靶板中心点挠度值, 实验结果与Nurick和Martin的理论计算结果符合得较好。
北京理工大学的王芳等人以能量守恒为基础, 采用理论分析和实验相结合的方法, 建立了冲击波载荷作用下方形靶板中心点挠度的半经验计算模型[1], 可较准确地预测载荷作用后靶板中心点挠度值, 从而对目标的毁伤程度进行评估。该半经验公式为
D0=0.0756(Weρ)1332reai+π3h√2σsρt (1) 式中:We为装药质量, ρ为装药密度, a为方形靶板半边长, i+为比冲量, h为靶板厚度, σs为靶板屈服强度, ρt为靶板密度, 反射系数re可由下式得到
re={2+6Δp1Δp1+7p0Δp1⩽0.8MPa50.8MPa⩽Δp1⩽5MPa (2) 式中:Δp1为入射的冲击波超压值, p0为标准大气压。
与无孔靶板相比, 有孔靶板的动态响应过程更为复杂。孔一方面降低了靶板的结构强度, 另一方面也减少了靶板上冲击波的承载面积。冲击波作用时孔附近存在应力集中, 靶板更容易变形。响应过程涉及爆轰产物的泄露, 冲击波的反射和叠加及预制孔靶板强度等效等问题, 通常具有几何非线性和物理非线性, 难以通过理论分析、工程计算手段对结果做精确的预测。
本研究采用数值计算与实验相结合的方法, 对装药爆炸后预制孔靶板的变形过程进行研究, 得到预制孔靶板中心点挠度与孔数、孔径之间的关系, 并在无孔靶板中心点挠度计算公式的基础上, 通过拟合给出预制孔靶板中心点挠度的计算公式。该公式可为有孔平板目标的毁伤评估、防爆安全距离的计算等提供依据, 也可作为破片、冲击波对平板目标联合毁伤研究的参考。
2. 数值计算
2.1 物理模型
物理模型如图 1所示, 尺寸为400 mm×400 mm×2 mm的Q235钢靶板四边固定, 1 kg长径比为1:1的柱形TNT装药底部距离靶板1 m, 装药单端点起爆。为方便研究, 靶板上预制圆孔均匀分布, 选取靶板上的孔数为4、9、16、25, 孔径取4、6、8、10、12 mm, 组合可得20种不同的工况。
2.2 数值模拟计算模型
2.2.1 网格模型及计算方法
利用AUTODYN-3D对问题进行建模和求解, 考虑到对称性, 4、16孔靶板可只考虑模型中1/4区域, 其余工况下则对全部区域进行计算[6]。空气域采用Euler算法, 除对称面外, 各边界设置空气流出边界条件, 实现Euler场边界能量的流出。为了提高计算效率, 利用AUTODYN-2D首先进行无限域中柱装药冲击波场计算, 然后通过FCT映射方法, 将冲击波场信息重新映射到空气域。这种方法既减少了从装药起爆到冲击波传播至靶板附近的计算量和时间, 也减小了Euler场的尺寸。为了准确计算孔附近的应力集中, 在孔附近20 mm范围内将网格细化尺寸为1 mm, 在z方向靶板附近25 mm范围内将网格细化尺寸为0.7 mm; 其他位置则采用渐变网格, 网格模型如图 2所示。
靶板采用Lagrange算法, 通过设置固定边界条件, 对四边节点各方向运动进行约束, 以描述靶板四周固定的特点。靶板采用均匀网格, 网格尺寸为2 mm。数值计算模型如图 3所示。
2.2.2 材料模型
空气采用理想气体状态方程(Equation of State, EOS)描述, 有
pg=(γ−1)ρgeg+pshift (3) 式中:pg为压力; eg为气体单位体积内能; V为相对体积; γ为多方指数, pshift为压力偏移量。对于空气模型:γ=1.4, ρg为1.225 kg/m3, eg在mm/mg/ms单位制下取2.068×105。
TNT炸药采用Jones-Wilkins-Lee状态方程描述
pexp=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωEexpV (4) 式中:pexp为爆轰压力, V为相对体积, Eexp为炸药单位体积内能, 其余为材料参数。对于TNT, 材料直接从AUTODYN材料库中选取。
靶板采用Q235钢材料, 状态方程选用Shock模型
p=pH+Γρ(e−eH) (5) 式中:Γρ为Grüneisen常数, 压力pH和能量eH计算方法为
pH=ρ0c20μ(μ+1)/[1−(s−1)μ]2 (6) eH=μpH/[2(μ+1)ρ0] (7) 式中:μ=ρ/ρ0-1, ρ0为材料初始密度, ρ为材料密度, c0为材料声速。
质点速度方程为
us=C1+S1uc+S2u2c (8) 式中:uc为冲击波速度, C1、S1和S2为实验常数。
靶板强度模型选用Bilinear Hardening模型, 其应力方程如下
σ=σs+Esε (9) Es=EEt/(E−Et) (10) 式中:σ为材料应力, ε为有效塑性应变, σs为屈服极限, E为杨氏模量, Et为切线模量。Q235材料参数及Bilinear Hardening模型参数见表 1
3. 模型实验
实验选择在开阔平整场地中进行。现场布置如图 4所示, 与数值计算中的参数保持一致。靶板由上压板及螺钉固定在靶架上; 悬挂支架用于悬挂炸药, 保证炸药位于靶板中心正上方, 药柱底面距靶板中心点高度H为1 m。
炸药选取长径比约1:1的圆柱形TNT, 质量1 kg。靶板材料为2 mm厚Q235钢。靶板有无孔和预制孔两种, 预制孔靶板上的孔均匀分布, 孔数可取4、9、16, 孔径可取4、6和8 mm, 共实验9种不同工况的预制孔靶板。
4. 结果与分析
实验后靶板状态如图 5所示, 靶板在爆炸冲击波作用下沿着载荷作用方向出现显著的塑性变形。靶板中心点挠度最大, 随着与中心点间距的增加, 变形区挠度逐渐减小。如图 6所示, 将靶板沿中轴线切开后, 放置于水平测量台, 对靶板中心点的挠度进行测量。
数值计算与实验得到的不同工况下靶板中心点的挠度值及相对误差见表 2。由表 2可知, 数值计算的最大误差为7.6%, 与实验结果符合较好。但数值计算结果相比实验结果均偏小, 这是因为数值计算中靶板边界为完全约束, 而实验中只能采用有限数目螺钉约束。在受到爆炸载荷冲击时, 靶板螺钉孔处变形较大(见图 6), 不能完全模拟理想状态, 因此实验得到的中心点的挠度值也较大。
表 2 部分实验数据与计算结果Table 2. Results of experiment and simulationNo. Hole number Hole diameters/
(mm)Deflection/(mm) Relative error/
(%)Experiment Simulation 1 0 72.1 68.3 5.6 2 4 4 74.2 70.2 5.7 3 4 6 76.2 70.8 7.6 4 9 4 74.3 72.4 2.6 5 9 6 78.6 73.7 6.6 6 16 4 76.2 75.3 1.2 7 16 6 81.1 78.3 3.6 图 7为数值计算所得无孔靶板和∅10 mm、16孔的预制孔靶板在爆炸冲击波作用0.5 ms时的应力云图。由图 7可见, 预制孔靶板的应力云图更复杂, 孔附近存在应力集中区域, 靶板上高应力区域的面积也明显增加。
数值计算得到的一些典型靶板中心点变形量随时间变化曲线见图 8。在孔数不同的情况下, 靶板中心点的挠度值随孔径的变化曲线见图 9。由图 9可知, 在孔密度较低时(4孔、9孔), 与无孔靶板相比, 预制孔靶板中心点的挠度有所增加, 但增加幅度并不明显; 孔密度较高时, 靶板中心点的挠度值有明显增加。靶板中心点的挠度增益随靶板孔直径的变化也有类似的规律。
对以上结果进行分析, 引入系数Kf来表征预制孔靶板相比于无孔靶板的变形增益, 此系数为靶板孔密度nf(个/米2)、孔直径df(mm)和的靶板厚度h0(mm)的函数, 即
Kf=Df/D0=f(nf,df/h0)=f(x1,x2) (11) 式中:Df为靶板中心点的挠度值, D0为无孔靶板中心点的挠度值, 自变量x1为nf, x2为df/h0。考虑到孔密度和孔径相互无影响, 因此假设x1和x2为独立参量, 利用数值计算结果拟合得到的Kf的二项式表达式为
Kf={1x1=0 或 x2=0f(x1,x2)=(c1+c2x1+c3x21)(d1+d2x2+d3x22)0<x1<160,0<x2<6 (12) 式中:c1=0.148 8, c2=3.07×10-4, c3=-8.569 4, d1=6.120 5, d2=0.235 65, d3=-0.015 06。(12)式与爆炸冲击波载荷下无孔靶板中心点挠度(1)式相结合, 可用于计算无孔以及不同孔密度、孔径情况下靶板中心点挠度。由于(12)式是通过对有限的数值计算结果进行拟合得到的, 因此只在一定范围内适用。而超出应用边界时, (12)式的可靠性则还需对更多的工况进行数值计算和实验验证。
5. 结论
利用AUTODYN软件对爆炸冲击波作用下预制孔靶板的动态响应过程进行了数值计算和实验验证, 得到了不同孔数、孔径工况下靶板中心点的挠度值; 并以爆炸冲击波作用下无孔靶板中心点挠度公式为基础, 利用数值计算结果, 通过拟合建立了可计算不同预制孔密度、孔径靶板中心点挠度的计算公式。该公式可用于预测爆炸冲击波载荷下有孔靶板中心点的挠度及毁伤评估, 也可为破片对平板目标先穿孔、再由冲击波作用的联合毁伤模型的研究提供参考, 具有一定的工程应用价值。
-
Material Material parameter EOS or strength model ρ/(g/cm3) σs/(MPa) E/(GPa) Poisson ratio S1 Γρ C1/(m/s) Et/(MPa) Q235 steel 7.83 235 206 0.3 1.49 2.17 4 569 500 表 2 部分实验数据与计算结果
Table 2. Results of experiment and simulation
No. Hole number Hole diameters/
(mm)Deflection/(mm) Relative error/
(%)Experiment Simulation 1 0 72.1 68.3 5.6 2 4 4 74.2 70.2 5.7 3 4 6 76.2 70.8 7.6 4 9 4 74.3 72.4 2.6 5 9 6 78.6 73.7 6.6 6 16 4 76.2 75.3 1.2 7 16 6 81.1 78.3 3.6 -
[1] 王芳, 冯顺山, 俞为民.爆炸冲击波作用下靶板的塑性大变形响应研究[J].中国安全科学学报, 2003, 13(3): 59-61. http://d.wanfangdata.com.cn/periodical/zgaqkxxb200303016Wang F, Feng S S, Yu W M. Study on large plastic deformation response of target plate under explosive blast wave[J]. China Safety Science Journal, 2003, 13(3): 59-61. (in Chinese) http://d.wanfangdata.com.cn/periodical/zgaqkxxb200303016 [2] Jones N, Uran T O, Tekin S A. The dynamic plastic behavior of fully clamped rectangular plates[J]. Int J Solids Structures, 1970(6): 1499-1512. [3] 何建, 肖玉凤, 陈振勇, 等.空爆载荷作用下固支矩形钢板的塑性极限变形[J].哈尔滨工业大学学报, 2007, 39(2): 310-312. http://www.cqvip.com/Main/Detail.aspx?id=23902325He J, Xiao Y F, Chen Z Y, et al. Plastic limited deformation analysis of the clamped steel plate subjected to air non-contact explosions[J]. Journal of Harbin Institute of Technology, 2007, 39(2): 310-312. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=23902325 [4] Nurick G N, Martin J B. Deformation of thin plates subjected to impulsive loading[J]. Int J Impact Eng, 1989, 8(2): 159-169. doi: 10.1016/0734-743X(89)90014-6 [5] Chung Kim Yuen S, Nurick G N, Versterb W, et al. Deformation of mild steel plates subjected to large-scale explosions[J]. Int J Impact Eng, 2008, 35(8): 684-703. doi: 10.1016/j.ijimpeng.2008.02.001 [6] 石磊, 杜修力, 樊鑫.爆炸冲击波数值计算网格划分方法研究[J].北京工业大学学报, 2010, 36(11): 1465-1470.Shi L, Du X L, Fan X. A study on the mesh generation method for numerical simulation of blast wave[J]. Journal of Beijing Polytechnic University, 2010, 36(11): 1465-1470. (in Chinese) [7] Century Dynamics Inc. AUTODYN Theory Manual R4.3[M]. Concord, California, 2005. [8] 陈燕, 王振生, 米双山.爆炸冲击波对靶板的损伤仿真与评估研究[J].军械工程学院学报, 2009, 21(1): 10-14. http://d.wanfangdata.com.cn/Periodical/jxgcxyxb200901003Chen Y, Wang Z S, Mi S S. Simulation and assessment of target plates damage by shock wave of certain explosion[J]. Journal o f Ordnance Engineering College, 2009, 21(1): 10-14. (in Chinese) http://d.wanfangdata.com.cn/Periodical/jxgcxyxb200901003 [9] 张世臣, 米双山.爆炸冲击波损伤靶板问题中的流固耦合算法[J].科学技术与工程, 2008, 8(7): 1655-1660. http://d.wanfangdata.com.cn/Periodical/kxjsygc200807003Zhang S C, Mi S S. Researching of fluid-solid coupling arithmetic in target damnifying by explosive shock wave[J]. Science Technology and Engineering, 2008, 8(7): 1655-1660. (in Chinese) http://d.wanfangdata.com.cn/Periodical/kxjsygc200807003 -