
Citation: | DU Ning, ZHANG Xianfeng, XIONG Wei, DING Li, WANG Jipeng, LIU Chuang. Experimental Study on the Penetration of Steel Fragments with Different Hardness into Q235A Steel Plate[J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 055102. doi: 10.11858/gywlxb.20180631 |
光滑粒子流体动力学(Smoothed Particle Hydrodynamics, SPH)方法中的搜索算法较耗时,即每个时间步都要对领域粒子进行搜索,粒子越多,耗时情况越突出,与有限元法相比,SPH方法的计算效率要低得多。为了解决SPH方法计算效率低的问题,Johnson等[1-2]和Attaway等[3]将有限元与SPH方法相结合,提出了SPH-FEM耦合的算法,即:在小变形区域使用有限元法,大变形区域仍使用SPH方法。该方法不仅提高了计算效率,而且适应性较强。
目前,采用SPH方法对爆炸焊接进行数值模拟的相关报道较少,而且多采用二维SPH方法。Tanaka[4]采用SPH方法对爆炸焊接的斜碰撞过程进行了数值模拟,成功地模拟出射流、波形和涡旋,波长的模拟结果相对实验结果偏大。李晓杰等[5]采用SPH方法及热塑性流体力学模型对爆炸复合板的斜碰撞过程中出现的界面波进行了数值模拟,模拟结果与张登霞等[6-7]实验结果的一致性较好。刘江等[8]利用AUTODYN软件中的SPH方法模拟了爆炸复合的斜碰撞,结合模拟中有效塑性变形、温度及剪切应力呈现的变化规律发现,爆炸复合的结合机理集塑性变形、熔化和扩散为一体。本研究将采用三维SPH方法对双面爆炸焊接过程进行模拟,将其结果与实验及理论结果进行对比,分析SPH-FEM耦合方法对爆炸焊接模拟的有效性。
以前期45钢/Q235钢双面爆炸焊接实验[9]为基础,考虑到计算效率,利用LS-DYNA建立如图 1及图 2所示的两组双面爆炸焊接SPH-FEM耦合的三维真实计算模型,选用的炸药为乳化炸药(玻璃微球的质量分数为5%),计算模型中基板和复板的材料、尺寸、间隙(δ)及药厚如表 1所示。起爆方式为点起爆。
Calculationmodel | Flyer plate | Base plate | Gap δ/mm |
Size of explosive/(mm×mm×mm) | |||
Material | Size/(mm×mm×mm) | Material | Size/(mm×mm×mm) | ||||
Ⅰ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×10 | |
Ⅱ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×5 |
基、复板采用3D Solid 164实体单元,单元边长为0.1 cm;炸药划分为光滑粒子,粒子的大小Δr取为0.1 cm。考虑到模型的对称性,为了提高计算效率,采用1/2模型进行计算。单位制为cm-g-μs。
数值计算中乳化炸药采用高能燃烧模型[10-11]及JWL状态方程[12]。JWL状态方程表达式为
p=AJWL(1−ωR1v)e−R1v+BJWL(1−ωR2v)e−R2v+ωE0v |
(1) |
式中:AJWL、BJWL、R1、R2和ω为材料参数;p为爆轰产物压力,GPa;E0为初始比内能,kJ/cm3;v为爆轰气体产物的相对比容,为无量纲量。炸药的相关参数见表 2,其中:ρ为密度,D为炸药爆速。
数值计算中,基、复板均采用Mie-Grüneisen状态方程[14]和Johnson-Cook材料模型[15]。Johnson-Cook材料模型的形式如下
σ=(A+Bεnp)(1+Cln˙ε∗p)(1−T∗m) |
(2) |
式中:εp为有效塑性应变;˙ε∗p=˙εp/˙ε0p为有效塑性应变率,其中˙ε0p为参考应变率;A、B、C、m及n为与材料相关的常数;无量纲温度T*表示为T*=(T-Tr)/(Tm-Tr),其中Tr为室温, Tm为熔点。45钢选用与Q235钢相同的Johnson-Cook材料模型参数,具体参数如表 3所示。
图 3所示是爆炸焊接结束时复板的竖向位移云图。由图 3可看出,复板的位移大致均为6 mm,表明基、复板已完全复合。为了更直观地观察复板单元位移的变化情况,在复板上选择3个特征单元(431 806、437 359、444 788),输出其位移-时间曲线,如图 4所示。由图 4可看出,特征单元的竖向位移均略大于间隙(6 mm),这是由于在爆炸载荷作用下复板有一定程度的减薄率所致。
图 5所示是一对分别取自基板与复板结合界面处的特征单元(基板单元:798 751;复板单元:416 251),特征单元的选取与前期实验[9]中金相试样的取样位置一致。
图 6所示是这对特征单元的速度-时间曲线。可以看出,基板在碰撞前有一个正的速度峰;该现象的产生如文献[17]所述,是由于爆轰产物不断堆积以及前碰撞点在待复合区产生的振动能所致。复板上所取单元的最大碰撞速度为897 m/s。
图 7所示是在复板结合界面处所选取的3个特征单元(410 476、416 251、420 976)。图 8所示是这3个特征单元的速度-时间曲线。
由图 8可以看出,随着距起爆端距离的增加,复板的碰撞速度增大。由文献[17]的结论可知,该现象是由于基板与复板的碰撞在金属板的待复合区产生了强烈振动引起的。
图 9所示是在结合界面处选取的3个特征单元(415 576、418 051、419 776),单元415 576取在复板中心处,与前期实验[7]中取样做金相观察的位置一致。图 10所示是3个特征单元的压力历程。
由图 10可以看出,随着距起爆端距离的增加,复板的碰撞压力增大。由文献[17]的结论可知,该现象是爆轰产物不断堆积以及前碰撞点在金属板待复合区振动能不断增加的共同作用结果。
图 11所示是爆炸焊接结束时复板的竖向位移云图。由图 11可看出,复板的位移大致均为6 mm,表明基、复板已完全复合。为了更加直观地观察复板单元位移的变化情况,在复板上选择3个特征单元(432 182、438 034、443 960),输出其位移-时间曲线,如图 12所示。由图 12可看出,特征单元的竖向位移均略大于6 mm,但较10 mm药厚下的竖向位移小。这是由于5 mm药厚下的爆炸载荷作用比10 mm药厚下小,导致5 mm药厚下的复板减薄率比10 mm药厚下低。
图 13所示是一对分别取自基板与复板结合界面处的特征单元(基板单元:799 201;复板单元:416 701),特征单元的选取与前期实验[9]中金相试样的取样位置一致。
图 14所示是这对特征单元的速度-时间曲线, 可以看出,基板在碰撞前也有一个正的速度峰。复板上所取单元的最大碰撞速度为565 m/s。
图 15所示是在复板结合界面处所选取的3个特征单元(411 976、417 001、423 826)。图 16所示是这3个特征单元的速度-时间曲线。由图 16可以看出,随着距起爆端距离的增加,复板的碰撞速度增大。
图 17所示是在结合界面处选取的3个特征单元(416 326、418 801、422 776),单元416 326取在复板中心处,与前期实验[7]中取样做金相观察的位置一致。图 18所示是这3个特征单元的压力历程。由图 18可以看出,随着距起爆端距离的增加,复板的碰撞压力增大。
由图 6可以看出,10 mm药厚下复板的最大碰撞速度为897 m/s。由图 14可以看出,5 mm药厚下复板的最大碰撞速度为565 m/s。利用前期工作[18]中提到的3种理论公式(Gurney公式、Aziz公式、Deribas公式)计算了复板的碰撞速度,如表 4、表 5所示,并与数值模拟结果进行了比较。由表 4和表 5可以看出:Gurney公式和Aziz公式的计算结果均存在较大的偏差;而由Deribas公式计算的两组结果与数值模拟结果较接近,误差均未超过5%,且与前期实验结果较吻合,证明了SPH-FEM耦合算法的可靠性。
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.75 | 1 089 | 897 | -21.0 |
Aziz | 0.75 | 711 | 897 | 20.0 |
Deribas | 0.75 | 853 | 897 | 4.9 |
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.45 | 863 | 565 | -52.7 |
Aziz | 0.45 | 480 | 565 | 15.0 |
Deribas | 0.45 | 576 | 565 | -1.9 |
由图 10可以看出,10 mm药厚下复板单元415 576处的碰撞压力为17.08 GPa。由图 18可以看出,5 mm药厚下复板单元416 326处的碰撞压力为11.25 GPa。
Ezra等提出的碰撞压力的计算公式为[16]
p=ρ1vs,1vp1+ρ1vs,1ρ2vs, 2 |
(3) |
式中:vs, 1、vs, 2分别表示复板、基板的声速,m·s-1;ρ1、ρ2分别表示复板、基板的密度,g·cm-3;vp表示复板的碰撞速度,m·s-1。
结合表 4和表 5中3种理论公式计算得到的碰撞速度,通过(3)式可得到复板的碰撞压力,表 6及表 7为其理论计算值与数值模拟结果的比较。可见:Gurney公式和Aziz公式的计算结果均存在较大的偏差;而由Deribas公式计算的两组结果与数值模拟结果较接近,误差均未超过5%,说明Deribas公式和SPH-FEM耦合方法对双面爆炸焊接具有较好的指导意义。
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 22.08 | 17.08 | -29.3 |
Aziz | 14.42 | 17.08 | 15.6 |
Deribas | 17.30 | 17.08 | -1.3 |
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 17.50 | 11.25 | -55.6 |
Aziz | 9.73 | 11.25 | 13.5 |
Deribas | 11.68 | 11.25 | -3.8 |
利用LS-DYNA软件和SPH-FEM耦合方法对前期双面爆炸焊接实验进行了三维数值模拟,并将模拟结果与实验及理论计算结果进行了对比,得到如下结论。
(1) 10 mm药厚和5 mm药厚下复板位移均略大于间隙值6 mm,这是由于爆轰载荷作用下复板有一定的减薄率所致。
(2) 10 mm药厚下,复板中部的最大碰撞速度为897 m/s,碰撞压力为17.08 GPa;5 mm药厚下,复板中部的最大碰撞速度为565 m/s,碰撞压力为11.25 GPa。通过与3种理论公式(Gurney公式、Aziz公式、Deribas公式)计算得到的碰撞速度进行比较发现,数值模拟结果与Deribas公式的计算结果较接近,误差较小,且与实验结果较吻合,证明了SPH-FEM耦合方法用于双面爆炸复合模拟的有效性,同时Deribas公式和SPH-FEM耦合方法对双面爆炸复合具有较好的指导意义。
(3) 10 mm药厚和5 mm药厚下复板的碰撞速度及碰撞压力均随着距起爆端距离的增加而增大,该现象是由于爆轰产物的不断堆积和前碰撞点在金属板待复合区振动能的不断增加共同作用的结果。
[1] |
金山, 汤铁钢, 孙学林, 等. 不同热处理条件下45钢柱壳的动态性能 [J]. 爆炸与冲击, 2006, 26(5): 423–428. doi: 10.3321/j.issn:1001-1455.2006.05.006
JIN S, TANG T G, SUN X L, et al. Dynamic characteristics of 45 steel cylinder shell by different heat treatment conditions [J]. Explosion and Shock Waves, 2006, 26(5): 423–428. doi: 10.3321/j.issn:1001-1455.2006.05.006
|
[2] |
DEMIR T, ÜBEYLI M, YILDIRIM R O. Effect of hardness on the ballistic impact behavior of high-strength steels against 7.62-mm armor piercing projectiles [J]. Journal of Materials Engineering and Performance, 2009, 18(2): 145–153. doi: 10.1007/s11665-008-9288-3
|
[3] |
DIKSHIT S N, KUTUMBARAO V V, SUNDARARAJAN G. The influence of plate hardness on the ballistic penetration of thick steel plates [J]. International Journal of Impact Engineering, 1995, 16(2): 293–320. doi: 10.1016/0734-743X(94)00041-T
|
[4] |
吴广, 冯顺山, 董永香, 等. 不同硬度弹丸对中厚钢板作用的试验研究 [J]. 兵工学报, 2010, 31(1): 260–263.
WU G, FENG S S, DONG Y X, et al. Experimental study on penetration effects of projectiles with different hardnesses on target [J]. Acta Armamentarii, 2010, 31(1): 260–263.
|
[5] |
陈小伟, 张方举, 梁斌, 等. A3钢钝头弹侵彻45钢板破坏模式的试验研究 [J]. 爆炸与冲击, 2006, 26(3): 199–207. doi: 10.3321/j.issn:1001-1455.2006.03.002
CHEN X W, ZHANG F J, LIANG B, et al. Three modes of penetration mechanics of A3 steel cylindrical projectiles impact onto 45 steel plates [J]. Explosion and Shock Waves, 2006, 26(3): 199–207. doi: 10.3321/j.issn:1001-1455.2006.03.002
|
[6] |
徐豫新, 任杰, 王树山. 钨球正侵彻下低碳钢板的极限贯穿厚度研究 [J]. 北京理工大学学报, 2017, 37(6): 551–556.
XU Y X, REN J, WANG S S. Research on perforation limit thickness of low carbon steel plates impacted normally by tungsten spheres [J]. Transaction of Beijing Institute of Technology, 2017, 37(6): 551–556.
|
[7] |
任杰, 徐豫新, 王树山. 超高强度平头圆柱形弹体对低碳合金钢板的高速撞击实验 [J]. 爆炸与冲击, 2017, 37(4): 629–636. doi: 10.11883/1001-1455(2017)04-0629-08
REN J, XU Y X, WANG S S. High velocity impact of low-carbon alloy steel plates by ultra-high strength blunt projectiles [J]. Explosion and Shock Waves, 2017, 37(4): 629–636. doi: 10.11883/1001-1455(2017)04-0629-08
|
[8] |
赵丽俊, 焦志刚, 李晓婕, 等. 预制破片侵彻均质装甲钢的极限穿透速度 [J]. 爆炸与冲击, 2018, 38(1): 183–190. doi: 10.11883/bzycj-2016-0116
ZHAO L J, JIAO Z G, LI X J, et al. Critical penetration velocity of prefabricated fragment in penetrating homogeneous armor steel plate [J]. Explosion and Shock Waves, 2018, 38(1): 183–190. doi: 10.11883/bzycj-2016-0116
|
[9] |
周楠, 王金相, 张亚宁, 等. 球形破片侵彻下钢/铝复合钢板的失效模式与吸能机理 [J]. 爆炸与冲击, 2018, 38(1): 66–75. doi: 10.11883/bzycj-2016-0131
ZHOU N, WANG J X, ZHANG Y N, et al. Failure mode and energy absorption mechanism of steel/aluminum composite plates impacted by spherical fragment [J]. Explosion and Shock Waves, 2018, 38(1): 66–75. doi: 10.11883/bzycj-2016-0131
|
[10] |
BACKMAN M E, GOLDSMITH W. The mechanics of penetration of projectiles into targets [J]. International Journal of Engineering Science, 1978, 16(1): 1–99. doi: 10.1016/0020-7225(78)90002-2
|
[11] |
陈继恩. 基于应力三轴度的材料失效研究 [D]. 武汉: 华中科技大学, 2009: 23–34.
CHEN J E. Research of material failure basic on stress triaxiality [D]. Wuhan: Huazhong University of Science and Technology of China, 2009: 23–34.
|
[12] |
谈庆明. 量纲分析 [M]. 合肥: 中国科学技术大学出版社, 2005: 9–18.
TAN Q M. Dimensional analysis [M]. Hefei: University of Science and Technology of China Press, 2005: 9–18.
|
[13] |
毛亮, 王华, 姜春兰, 等. 钨合金球形破片侵彻陶瓷/DFRP复合靶的弹道极限速度 [J]. 振动与冲击, 2015, 34(13): 1–5.
MAO L, WANG H, JIANG C L, et al. Ballistic limit velocity of tungsten alloy spherical fragment penetrating ceramic/DFRP composite target plates [J]. Journal of Vibration and Shock, 2015, 34(13): 1–5.
|
[14] |
BANERJEE B. The mechanical threshold stress model for various tempers of AISI 4340 steel [J]. International Journal of Solids and Structures, 2005, 44(3): 834–859.
|
[1] | WANG Xiao-Wei, DONG Lian-Ke, LONG Qi-Wei. The Computer Simulation of Fractal Dimension-Ality with Perimeter-Area Relation[J]. Chinese Journal of High Pressure Physics, 1991, 5(2): 124-129 . doi: 10.11858/gywlxb.1991.02.007 |
[2] | SU Fang, CHE Rong-Zheng, XU Wei, YAN Fei-Na. Measurement and Analysis of Strain and Calibration of Pressure in Piston-Cylinder Type Apparatus with 4.5 GPa Hydrostatic Pressure[J]. Chinese Journal of High Pressure Physics, 1990, 4(1): 63-71 . doi: 10.11858/gywlxb.1990.01.010 |
[3] | PENG Jian, ZHAO Yi-Jun, ZHANG Zheng-Wen, YUAN Jian-Min. Calculation for the Thermal Stress Damage in the Thin Plate of Al Irradiated by CW Laser Beam[J]. Chinese Journal of High Pressure Physics, 1990, 4(2): 130-136 . doi: 10.11858/gywlxb.1990.02.009 |
[4] | YUAN Yong-Hua, LIU Chang-Ling, HAN Li-Shi, GUI Yuan-Zhen, LI Qi-Min. Measurement of the Vaporization Pressure on Aluminum Target Irradiated by Laser Beam[J]. Chinese Journal of High Pressure Physics, 1990, 4(2): 114-117 . doi: 10.11858/gywlxb.1990.02.006 |
[5] | ZHOU Zhi-Kui, CHEN Qiu-Hua. Numerical Simulation of the Performance of Electromagnetic Railguns[J]. Chinese Journal of High Pressure Physics, 1989, 3(4): 308-314 . doi: 10.11858/gywlxb.1989.04.008 |
[6] | DONG Yu-Bin, SU Lin-Xiang, CHEN Da-Nian, JING Fu-Qian, HAN Jun-Wan, FENG Jia-Bo. Numerical Simulation on the Spallation of a Steel Cylindrical Shell Imploded under Slipping Detonation[J]. Chinese Journal of High Pressure Physics, 1989, 3(1): 1-10 . doi: 10.11858/gywlxb.1989.01.001 |
[7] | WANG Yi-Feng, SU Wen-Hui, QIAN Zheng-Nan, MA Xian-Feng, YAN Xue-Wei. A Study on the High Temperature-High Pressure Stability and the X-Ray Phase Analysis for Compounds R2Fe4/3W2/3O7 with Pyrochlore Structure[J]. Chinese Journal of High Pressure Physics, 1988, 2(4): 296-304 . doi: 10.11858/gywlxb.1988.04.002 |
[8] | SHEN Zhu-Tong, DAI Shou-Yu, CHEN Li-Quan, CHI Yong-Qian. Some Investigations on the High Pressure Synthesis of Niobium Trisulfide and XPS Analysis[J]. Chinese Journal of High Pressure Physics, 1988, 2(1): 10-16 . doi: 10.11858/gywlxb.1988.01.002 |
[9] | DONG Yu-Bin, ZAHNG Wan-Jia, JING Fu-Qian, HAN Jun-Wan, CHEN Da-Nian, SU Lin-Xiang, FENG Jia-Bo. Numerical Analysis for Dynamic Damage Processes and LY-12 Aluminum Spallations[J]. Chinese Journal of High Pressure Physics, 1988, 2(4): 305-312 . doi: 10.11858/gywlxb.1988.04.003 |
[10] | CHEN Dong-Quan, XIE Guo-Qiang. Molecular Dynamics Simulation of Polymorphous Transitions[J]. Chinese Journal of High Pressure Physics, 1987, 1(1): 50-57 . doi: 10.11858/gywlxb.1987.01.007 |
Calculationmodel | Flyer plate | Base plate | Gap δ/mm |
Size of explosive/(mm×mm×mm) | |||
Material | Size/(mm×mm×mm) | Material | Size/(mm×mm×mm) | ||||
Ⅰ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×10 | |
Ⅱ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×5 |
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.75 | 1 089 | 897 | -21.0 |
Aziz | 0.75 | 711 | 897 | 20.0 |
Deribas | 0.75 | 853 | 897 | 4.9 |
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.45 | 863 | 565 | -52.7 |
Aziz | 0.45 | 480 | 565 | 15.0 |
Deribas | 0.45 | 576 | 565 | -1.9 |
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 22.08 | 17.08 | -29.3 |
Aziz | 14.42 | 17.08 | 15.6 |
Deribas | 17.30 | 17.08 | -1.3 |
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 17.50 | 11.25 | -55.6 |
Aziz | 9.73 | 11.25 | 13.5 |
Deribas | 11.68 | 11.25 | -3.8 |
Calculationmodel | Flyer plate | Base plate | Gap δ/mm |
Size of explosive/(mm×mm×mm) | |||
Material | Size/(mm×mm×mm) | Material | Size/(mm×mm×mm) | ||||
Ⅰ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×10 | |
Ⅱ | 45 steel | 300×150×2 | Q235 | 300×150×16 | 6 | 300×150×5 |
ρ/(g·cm-3) | D/(m·s-1) | AJWL/GPa | BJWL/GPa | R1 | R2 | ω | E0/(kJ·cm-3) |
1.12 | 4 510 | 326.42 | 5.808 9 | 5.80 | 1.56 | 0.57 | 3.323 |
ρ/(g·cm-3) | G/GPa | A/GPa | B/GPa | C | n | m | Tm/K | Tr/K |
7.83 | 77 | 0.792 | 0.51 | 0.014 | 0.26 | 1.03 | 1 793 | 294 |
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.75 | 1 089 | 897 | -21.0 |
Aziz | 0.75 | 711 | 897 | 20.0 |
Deribas | 0.75 | 853 | 897 | 4.9 |
Theoreticalformula | Massfraction | Collision velocity/(m·s-1) | Error/% | |
Theoretical calculation[18] | Simulation | |||
Gurney | 0.45 | 863 | 565 | -52.7 |
Aziz | 0.45 | 480 | 565 | 15.0 |
Deribas | 0.45 | 576 | 565 | -1.9 |
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 22.08 | 17.08 | -29.3 |
Aziz | 14.42 | 17.08 | 15.6 |
Deribas | 17.30 | 17.08 | -1.3 |
Theoreticalformula | Collision pressure/GPa | Error/% | |
Calculation | Simulation | ||
Gurney | 17.50 | 11.25 | -55.6 |
Aziz | 9.73 | 11.25 | 13.5 |
Deribas | 11.68 | 11.25 | -3.8 |