• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

激光选区熔化增材制造GP1不锈钢动态拉伸力学响应与层裂破坏

史同亚 刘东升 陈伟 谢普初 汪小锋 王永刚

王振宁, 尹建平, 伊建亚, 李旭东. 柱形装药近地动爆冲击波周向传播规律研究[J]. 爆炸与冲击, 2023, 43(6): 063201. doi: 10.11883/bzycj-2022-0313
引用本文: 史同亚, 刘东升, 陈伟, 谢普初, 汪小锋, 王永刚. 激光选区熔化增材制造GP1不锈钢动态拉伸力学响应与层裂破坏[J]. 爆炸与冲击, 2019, 39(7): 073101. doi: 10.11883/bzycj-2019-0015
WANG Zhenning, YIN Jianping, YI Jianya, LI Xudong. A study on the circumferential propagation law of the shock waves produced by the near ground dynamic explosion of cylindrical charge[J]. Explosion And Shock Waves, 2023, 43(6): 063201. doi: 10.11883/bzycj-2022-0313
Citation: SHI Tongya, LIU Dongsheng, CHEN Wei, XIE Puchu, WANG Xiaofeng, WANG Yonggang. Dynamic tensile behavior and spall fracture of GP1 stainless steel processed by selective laser melting[J]. Explosion And Shock Waves, 2019, 39(7): 073101. doi: 10.11883/bzycj-2019-0015

激光选区熔化增材制造GP1不锈钢动态拉伸力学响应与层裂破坏

doi: 10.11883/bzycj-2019-0015
基金项目: 科学挑战专题(TZ2018001)
详细信息
    作者简介:

    史同亚(1992- ),男,硕士,819330522@qq.com

    通讯作者:

    王永刚(1976- ),男,博士,教授,wangyonggang@nbu.edu.cn

  • 中图分类号: O347.3

Dynamic tensile behavior and spall fracture of GP1 stainless steel processed by selective laser melting

  • 摘要: 采用选择性激光熔化增材制造技术,制备了GP1不锈钢单轴拉伸板条试样和层裂圆片试样,并对材料微观结构进行了表征。借助Zwick-HTM5020 高速拉伸试验机,并结合数字图像相关性全场应变测量技术,开展了增材制造GP1不锈钢材料的轴向拉伸力学性能实验研究,得到了不同应变率下材料的拉伸应力-应变曲线,结果显示:(1) GP1不锈钢流动应力具有比较显著的应变强化效应;(2)通过回收试样的电子背散射衍射表征,发现GP1不锈钢在拉伸变形过程中会发生奥氏体与马氏体之间的相变;(3) GP1不锈钢的屈服应力随着应变率呈幂指数增大,断裂应变在中低应变率下保持不变,但在高应变率下则显著减小。采用一级轻气炮实验装置和激光干涉粒子速度测量技术,开展了增材制造GP1不锈钢的层裂实验,发现GP1不锈钢的层裂强度随着飞片撞击速度增大而减小。单轴拉伸试样断口和层裂试样断口的显微分析结果表明:随着应变率增大,单轴拉伸断裂模式和断裂机理都发生了转变;层裂损伤易成核于激光熔池边界线的交汇处,断口韧窝形貌明显区别于单向拉伸断口。
  • 在常规炮弹、火箭弹、导弹等弹种中,战斗部的装药形状多为圆柱形或近似为圆柱形,装药抵近地面时的姿态和速度可能对产生的近地爆炸冲击波的传播造成不同的影响。研究柱形装药在多种条件下的近地爆炸冲击波传播特征、压力参数对于战斗部威力的评估及冲击波防护具有重要意义。在装药的动爆研究方面,有较全面的动爆冲击波传播理论,并首先在球形装药的动爆研究中得到了理论推导。陈龙明等[1]研究发现,球形装药的超压、比冲量、正压作用时间以90°为分界线增强或减弱。聂源等[2]得到了球形装药的动态爆炸冲击波超压场计算模型表达式。蒋海燕等[3]研究了球形装药动爆速度影响冲击波场空间和位置参数的规律,其中关于柱形装药的动爆特征研究较少,且普遍为二维平面的数值模拟研究。

    在柱形装药的爆炸冲击波研究中,郗洪柱等[4]基于镜像法研究了柱装药近地空爆冲击波峰值超压空间转换模型,并进行了实验验证;郑监等[5]对不同圆锥角和不同柱形装药质量下锥形激波管内的冲击波传播过程进行了模拟。在近地爆炸研究中,有较多关于近地爆炸冲击波的反射理论,廖真等[6]对不同装药类型和形状产生的三波点轨迹进行了分析,并结合实验得到了不同装药峰值压力的差别;陈鑫等[7]研究了地面介质对近地马赫波超压的影响,并与刚性地面进行了对比。Wang等[8]基于ANSYS对近地爆炸进行了数值模拟,分析了不同地面对冲击波的反射作用;Lu等[9]研究了近地爆炸对地面建筑物结构响应的过程分析,并分析了土壤中冲击波的作用;Deng等[10]通过研究带壳装药的近地爆炸,得到了马赫波的传播与温度场变化的对比;Li等[11]通过AME(adaptive mesh enlargement)方法对网格进行了自适应调整,从而可以对大规模流体场进行计算;You等[12]研究了一种改进的有限测点爆炸超压数据插值重建方法,从而可以得到地面附近爆炸的峰值压力分布;任朋飞等[13]对炸药在刚性地面爆炸进行了实验与数值模拟研究,验证了ALE(arbitrary Lagrange-Euler)方法算法的正确;赵蓓蕾等[14]利用LS-DYNA分析了地面冲击波传播规律,得到了马赫杆、地面超压的衰减规律;成凤生等[15]比较了柱形与球形装药超压变化,研究了波阵面的传播规律;段晓瑜等[16]开展了不同材质炸药的近地反射超压的实验研究;杜红棉等[17]拟合出了近地爆炸冲击波马赫杆的三维三波点马赫波轨迹,验证了数值模拟结果的可靠性;姚成宝等[18]使用网格自适应技术,对炸药的空中强爆炸在不同爆炸高度下的冲击波地面反射过程进行了数值模拟;吴赛等[19]得出刚性地面的反射作用会显著增强爆炸效应,可对结构或构件造成更大的危害,提高了对工程中的反射效应最低炸高理论值。

    综合之前的研究成果可以看出,目前关于柱形装药的冲击波传播理论研究较少,且没有考虑实际情况下柱形装药战斗部的动态参数特征。因此,本文中,基于AUTODYN-3D软件对不同末弹道参数的柱形装药近地爆炸进行数值模拟,研究动爆条件下的柱形装药各参数对近地冲击波传播的规律,并对峰值压力、反射波产生、马赫波等冲击波参数和传播过程进行分析。

    以常规榴弹战斗部装药量作为参考(依据口径不同,装药量一般在1.75~5.00 kg),装药药量选取4.00 kg,柱形装药近地爆炸示意图如图1所示。炸药选择TNT炸药,起爆方式为中心点起爆,材料参数为AUTODYN材料库中的默认参数,TNT与空气的材料模型均采用Euler网格方法计算,TNT的状态方程参数见表1,表中:A、B、R1、R2ω为状态方程独立参数,E0、V分别为炸药比内能和比容。为避免网格尺寸影响计算精度,将0.859 kg的装药在网格尺寸分别为0.5、1.0、1.5、2.0 cm进行空爆数值模拟。将数值模拟结果与Henrych和Brode空爆公式[20]进行对比,随着网格精度的提高,近地压力值与Henrych公式计算结果逐渐逼近。结合廖真等[6]在近地爆炸数值模拟中的网格尺寸分析,验证了1.0 cm的网格尺寸与郭炜等[21]实验中的三波点高度具有较好的一致性。同时在三维模型中要兼顾空气域尺寸和计算速度,故选取空气域网格尺寸为1.0 cm。

    图  1  柱形装药近地动爆示意图
    Figure  1.  Schematic diagram of near ground dynamic explosion of cylindrical charge
    表  1  TNT炸药状态方程参数
    Table  1.  Parameters of equation of state for TNT explosive
    A/GPaB/GPaR1R2E0/GPaωV
    3.73773.7474.150.90.060.31
    下载: 导出CSV 
    | 显示表格

    建模分为2个方向:第1个方向为柱形装药前后方冲击波传播方向,以x轴为基准方向,在x轴的距原点5 m处放置柱形炸药;第2个方向为炸药侧方,以y轴为基准方向,在原点位置放置炸药。2个方向的建模网格数量分别为2 000 000和3 600 000。x方向与y方向的空气域宽度均可将炸药完全覆盖,前后方向建模尺寸分别为10.00、0.10和2.00 m;侧方向建模尺寸分别为5.00、0.36和2.00 m。空气的材料参数如表2所示,材料模型见表3。表中:ρ为空气密度,γ为绝热指数,e为空气的初始比内能。

    表  2  空气状态方程参数
    Table  2.  Parameters of equation of state for air
    ρ/(kg·m−3)γe/(MJ·kg−1)
    1.2251.4206.8
    下载: 导出CSV 
    | 显示表格
    表  3  材料模型
    Table  3.  Material model
    材料状态方程强度模型失效模型
    空气理想气体
    TNTJWL
    下载: 导出CSV 
    | 显示表格

    TNT炸药和空气分别采用JWL状态方程和理想气体状态方程。

    p=A(1ωR1V)eR1V+B(1ωR2V)eR2V+ωE0V
    (1)

    式中:p为爆轰产物压力。

    在空气域前后面、侧面和上面设置物质流出边界,剩余两面为自由反射边界,模拟对称方向和刚性地面。设置以装药为中心的辐射型测点不能满足不同姿态近地冲击波的传播特征,所以在相同高度上设置等距的测点,并沿z方向间隔相同高度设置,测点设置如图2所示。

    图  2  柱形装药近地爆炸建模及测点布置
    Figure  2.  Modeling and measuring point diagram of near ground explosion of cylindrical charge

    将数值计算模型与任朋飞等[13]所做的刚性地面爆炸冲击场实验进行对比,在混凝土路面上放置90 mm厚的钢板来模拟刚性地面,并对长径比为1∶3的0.2、0.8和2.0 kg等3种装药量的柱形装药进行了实验。在距爆心位置处的1.12和1.95 m设置了测点,建模如图3所示。测点压力的计算值与实验值的对比如表4所示。

    图  3  实验验证示意图
    Figure  3.  Schematic diagram of experimental verification
    表  4  测点压力的计算值与实验值对比
    Table  4.  Comparison between calculated values and experimental values on measuring point pressure
    TNT质量/kg测点距离/m压力/MPa误差/%
    实验值计算值
    2.01.121.231.4114.56
    2.01.950.690.724.30
    0.81.120.800.9012.86
    0.81.950.480.6024.17
    0.21.120.300.326.67
    0.21.950.090.1455.56
    下载: 导出CSV 
    | 显示表格

    实际情况下的刚性地面反射强度弱于数值模拟中的绝对反射,装药质量小导致网格精度变差,故数值模拟值与实际值存在误差。可以看出,随着装药质量的增加,误差不断减小,2.0 kg装药的平均误差为9.43%,所以对4.0 kg装药的数值模拟符合实际的刚性地面反射情况。

    1.3.1   弹前激波对近地爆炸的影响

    飞行状态下的柱形装药在前端面产生压缩波,随着运动速度的不断提高,在弹前产生的波逐渐从弱的压缩波变为强的压缩波,该冲击波是由于弹前压缩波的聚集、反射与叠加[22]形成的。在数值模拟中,先让弹丸在空气中飞行1.0 ms以上,使得装药在空气中产生较稳定的激波及压力场,然后对其进行起爆。在数值模拟中对比了静爆、动爆、静爆加激波和动爆加激波4种起爆类型,图45给出了长径比为2∶1的柱形装药在600和1 200 m/s速度下的稳定激波压力云图和矢量图。

    图  4  不同飞行速度下柱形装药的压力云图
    Figure  4.  Pressure nephogram of cylindrical charge at different flight speeds
    图  5  不同飞行速度下柱形装药的压力矢量图
    Figure  5.  Pressure vector diagrams of cylindrical charge at different flight speeds

    柱形装药在距离地面1.2 m处起爆,并在刚性界面反射传播形成马赫波,观测上述4种类型产生的爆炸压力场及峰值压力,4.0 ms时1 200 m/s速度下4种起爆类型产生的马赫波如图6所示。

    图  6  不同起爆条件下柱形装药的压力云图
    Figure  6.  Pressure nephograms of cylindrical charge under different initiation conditions

    由于爆炸产生的入射压力较大,是激波最大压力的200倍以上,在与地面接触产生的反射波是其压力的20倍,所以激波对马赫波的传播过程、马赫波的峰值压力影响有限,仅对起爆点上方入射波的传播有一定干扰。

    表5中可以看出,激波和动爆的联合作用与动爆的峰值压力差别很小,平均差0.12%,并且对马赫杆高度的主要影响为装药动爆速度的施加,因此在数值模拟中忽略弹前激波的影响。

    表  5  不同起爆条件下柱形装药的峰值压力
    Table  5.  Peak pressure of cylindrical charge under different initiation conditions
    初始速度/(m·s−1)2 ms时的压力/MPa4 ms时的压力/MPa
    激波加动爆动爆静爆激波加动爆动爆静爆
    6001.041.040.980.400.400.39
    1 2000.980.981.000.390.390.40
    下载: 导出CSV 
    | 显示表格
    1.3.2   地面沙土对近地爆炸的影响

    实际战场环境地面情况多样,沙土等软质地面不能视作刚体。AUTODYN中沙土选择Linear材料模型和von Mises强度模型[23],密度为1.6 g/cm3。计算中沙土采用Lagrange模型,数值模拟采用流固耦合计算方法,在0.2、0.4、0.6、0.8和1.0 m等5种近地高度下与刚性地面进行对比。

    刚性和沙土地面的对比如图7所示,可以看出沙土介质对马赫波的影响有限。

    图  7  不同地面条件下柱形装药的压力云图
    Figure  7.  Pressure nephograms of cylindrical charge under different ground conditions

    表6中看出,沙土对近地爆炸冲击波的峰值压力影响较为有限,0.5 ms时,峰值压力平均差1.76%,马赫杆高度平均差6.48%;1 ms时,峰值压力平均差1.17%,马赫杆高度平均差4.79%,说明数值模拟中沙土对近地冲击波的影响较小。

    表  6  不同地面条件下柱形装药的峰值压力和马赫杆高度
    Table  6.  Peak pressure and Mach stem height of cylindrical charge under different ground conditions
    地面类型起爆高度/m峰值压力/MPa马赫杆高度/mm
    0.5 ms1.0 ms0.5 ms1.0 ms
    刚性地面0.24.212.146501 380
    沙土地面0.24.122.066101 340
    刚性地面1.07.953.23 20 90
    沙土地面1.08.853.27 20 90
    下载: 导出CSV 
    | 显示表格

    近地爆炸时冲击波到达地面并反射形成反射波,反射波由于空气的压缩加热其速度大于正常产生的入射波,将逐渐追上入射波;之后反射波与入射波叠加形成马赫波;马赫波在水平方向传播的同时其高度不断扩大,反射波及马赫波形成的示意如图8所示。

    图  8  近地爆炸马赫波的产生过程示意图
    Figure  8.  Schematic diagram of Mach wave generation process of near ground explosion

    图9为炸高1.0 m,装药质量为4 kg,长径比为2 : 1的柱形装药在倾斜角为0°的马赫波产生过程。

    图  9  不同时刻下柱形装药的马赫波压力云图
    Figure  9.  Mach wave pressure nephogram of cylindrical charge at different times

    以0.2 m为间隔设置了0.2~1.8 m范围内的炸高,柱形装药长径比为2∶1的近地爆炸数值模拟对照。起爆高度会直接影响冲击波到达地面的时间,导致入射波到达地面的峰值压力随起爆高度增加而减少。国防工程设计草案[19]中,是否考虑刚性地面的反射增强效应的表达式为:

    HW1/30.35m/kg1/3
    (2)

    式中:H为爆心距地面的距离,m;W为TNT炸药当量,kg。

    按照式(2)计算得到高度为0.56 m,实际上反射效应在此高度仍然较大,因此对柱形装药在不同起爆高度下的近地峰值压力进行对比,如图10所示。从图10可以看出,峰值压力随爆心距和起爆高度整体呈现下降趋势,在2.0 m爆心距后的峰值压力差距明显减小。当爆心距为0.5 m时,0.4 m爆高的峰值压力是1.8 m爆高的3.6倍。峰值压力的不同造成了马赫杆的高度差距,2.0 ms时0.4 m爆高的马赫杆高度与1.6 m相比相差1.93 m。图11为爆高为1.0 m时水平方向上不同高度测点柱形装药的峰值压力-时间曲线,从0 m高度峰值压力曲线中可以看到,3.5 m测点处存在明显的2个波峰,并向后传递;在1.0 m高度的峰值压力有一定程度的衰减,而在2.0 m高度处测点的峰值压力几乎与自由空气中爆炸相近。

    图  10  不同爆心高度下柱形装药的峰值压力的对比
    Figure  10.  Comparison of peak pressures of the cylindrical charge under different heights of the explosion center
    图  11  不同高度测点柱形装药的压力曲线
    Figure  11.  Pressure curves of cylindrical charge at measuring points of different heights

    马赫波的产生与装药的比例炸高相关,依据极限入射角与比例炸高倒数的关系曲线来确定马赫波反射区的范围。当装药质量大并且炸心高度较低时,比例炸高倒数较大。从图12中看出,当比例炸高倒数大于0.8时,入射角将趋近于40°,所以马赫波在水平方向的发生位置与炸高近似。

    图  12  入射角与比例炸高的关系
    Figure  12.  Relationship between incident angle and proportional explosion height

    近地面冲击波斜反射的传播情况较复杂,将入射波质点速度u1分解为平行于地面的分量和垂直于地面的分量,垂直于地面的分量形成的反射波在被压缩和加热的空气中传播从而速度增大。反射波在入射角θ1较小时,水平方向的速度分量较小,随着θ1的增加,反射波在水平方向的速度分量增大。从某个θ1开始,反射波质点速度u2的垂直分量将大于u1的垂直分量,2个波阵面的交点A将离开地面。斜反射过程的示意图如图13所示。

    图  13  近地爆炸冲击波斜反射示意图
    Figure  13.  Schematic diagram of oblique reflection of shock wave produced by near ground explosion

    柱形装药径向冲击波比轴向冲击波传播速度快,在装药倾斜角度不大时,轴向入射波将率先到达地面并进行反射。随着倾斜角度的增加,轴向入射波与地面接触的位置不断向前移动,使得斜反射与地面的夹角不断减小,水平方向的分量增加而反射波的峰值压力不断减弱。所以装药前侧的马赫杆高度在开始大于后侧的马赫杆,之后被逐渐追上并超越。以15°和75°静爆条件下马赫波传播过程为例可以看出这一趋势,如图14所示。

    图  14  倾斜角度不同时柱形装药的马赫波压力云图
    Figure  14.  Mach wave pressure nephograms of cylindrical charges with different inclined angles

    在大于45°附近的倾斜装药,不会再产生两个波同时抵达地面的情况,此时径向入射波将率先在后侧产生反射波。装药后侧产生的马赫杆高度将始终大于前侧的马赫杆高度,图15所示为2种不同的入射波抵近地面的状态。

    图  15  倾斜角度不同时柱形装药的近地状态压力云图
    Figure  15.  Near ground pressure nephograms of cylindrical charges with different inclined angles

    随着倾斜角度的增加,装药两侧的峰值压力较为接近,传播速度及马赫杆高度随着角度的增加而逐渐增大,在3.5 ms时,90°比0°的马赫杆高出77.1%。通过倾斜角相同时装药前方、后方和侧方3个方向的传播过程对比,可以看出在同一时刻冲击波的传播差别,图16所示为45°倾斜角、4.0 ms时刻3个方向的冲击波传播过程。

    图  16  倾斜角为45°时不同方向柱形装药的压力云图
    Figure  16.  Pressure nephograms of cylindrical charge in different directions when the inclination angle is 45°

    图17给出了4.0 m处随倾斜角度变化3个方向的峰值压力变化。在装药前方的峰值压力随着倾斜角度的增加不断增大,并在75°时达到最大值;而在装药的侧方和后方的峰值压力随着倾斜角度的增加整体呈现递减的趋势,分别在60°及90°取得最小值。

    图  17  柱形装药的峰值压力随倾斜角度的变化
    Figure  17.  Peak pressure of cylindrical charge varied with inclined angle

    表7为倾斜角度对马赫杆高度的影响,其中3.0 ms时前方马赫杆高度最大增长幅度为24.8%,后方马赫杆高度最大增长幅度为88.1%。

    表  7  不同倾斜角度下柱形装药在前、后方的马赫杆高度
    Table  7.  Mach stem heights of cylindrical chargeswith different inclined angles at the front and rear
    角度/(°) 马赫杆高度/mm
    1.0 ms2.0 ms3.0 ms
    前方后方前方后方前方后方
    0 90 90480 4801 0101 010
    303101707401 0401 1801 850
    601303907501 1801 2601 900
    90130130610 6101 2501 250
    下载: 导出CSV 
    | 显示表格

    由于柱形装药形成的冲击波在径向的传播速度比轴向的快,随着倾斜角度的增大,冲击波在各方向呈现出不同的波形和峰值压力。图18为不同倾斜角度下柱形装药在3个方向的峰值压力变化。

    图  18  不同方向柱形装药的峰值压力随倾斜角度的变化
    Figure  18.  Variation of peak pressure of cylindrical charge with inclined angle in different directions

    在实际作战中战斗部一般为动态爆炸,使整个冲击波场向速度方向产生移动,影响冲击波流场的演化过程,造成装药沿速度方向前端压力升,后端压力降。以倾斜角为30°,爆心高度为1.0 m的2∶1长径比装药进行研究,柱形装药的初速度依次设置为0、300、600、900、1 200和1 500 m/s。

    图19为0和300 m/s的冲击波场生成过程对比,动爆会使整个冲击波场向速度方向移动,增大近地侧入射波的峰值压力而减小远地侧的入射波压力。动爆会加速装药前方的轴向冲击波,驱使其率先抵达地面,使得冲击波能量向前侧聚集,从而前侧马赫波峰值压力及马赫杆高度追上后侧。

    图  19  不同动爆速度下柱形装药的冲击波压力云图
    Figure  19.  Shock wave pressure nephogram of cylindrical charge under different dynamic explosion velocities

    当装药特性、测点位置等条件都相同时,在压力上升区域,装药速度越大,超压峰值的上升幅度越大,在压力下降区域,装药速度越大,超压峰值的下降幅度也越大。图20为4.0 m处不同方向的峰值压力随速度变化关系,从图20可以看出,随着动爆速度的提升,前方峰值压力显著增大。比较3.0 ms时不同速度的马赫杆高度发现,前方马赫杆高度最大增长幅度为156.2%,后方马赫杆高度最大幅度为76.2%,侧方马赫杆高度最大幅度为388.2%。

    图  20  柱形装药的峰值压力随动爆速度变化
    Figure  20.  Diagram of peak pressure of cylindrical charge varying with dynamic explosion velocity

    由于动爆速度指向地面,侧方冲击波受此影响使得压力区域快速下移并与地面接触,速度越快区域与地面的夹角越大,使得速度快的侧向马赫波高度将迅速升高。300和1 500 m/s的侧方冲击波传播过程对比见图21

    图  21  不同动爆速度下柱形装药在侧方的马赫波压力云图
    Figure  21.  Mach wave pressure nephogram of cylindrical charge on the side under different dynamic explosion velocities

    对比不同动爆速度在3个方向的峰值压力变化曲线,如图22所示。随着动爆速度的提升,在装药的前方及侧方的不同距离峰值压力都呈现增加的趋势,而装药后方趋势较不明显,且在距爆心2.0 m距离外的峰值压力无较大变化。

    图  22  不同方向柱形装药的峰值压力随动爆速度的变化
    Figure  22.  Peak pressure varying with dynamic explosion velocity for cylindrical charge in different directions

    装药长径比会对自由空气场中冲击波波形产生影响,由于柱形装药的入射波特点为径向波传播快、轴向波传播慢,随着高度的增加其波的差异会减弱。对1∶1、1.5∶1、2 ∶1、2.5∶1的4种长径比柱形装药在1.0 m高度、0 m/s、30°倾斜角条件下对装药前部、后部、侧部3个方向进行分析,图23为2.0 ms时刻不同长径比装药在前后方向的传播对比。

    图  23  不同长径比柱形装药在前、后方的马赫波压力云图
    Figure  23.  Mach wave pressure nephogram of cylindrical charges with different aspect ratios at the front and rear

    随着柱形装药长径比的增加,装药前侧的马赫杆高度逐渐降低,而装药后侧的马赫杆高度逐渐增加。前侧马赫杆高度差距在3.0 ms时可达1倍,后侧差距则不到0.5倍。图24为4.0 ms时刻装药侧方的冲击波传播过程,随着装药长径比的增加,侧方的马赫杆高度逐渐降低,3.0 ms时装药比例1∶1的马赫杆平均高度是2.5∶1高度的1.43倍。

    图  24  不同长径比柱形装药在侧方的马赫波压力云图
    Figure  24.  Mach wave pressure nephogram of cylindrical charges with different aspect ratios on the side

    图25所示为不同比例装药在前、后、侧3个方向的峰值超压随距爆心距离的变化关系,不同的比例装药在各个方向呈现出不同的变化趋势,在前侧的变化最显著。从表8不同长径比装药的马赫杆高度对比可以看出,随着柱形装药长径比的增加,在前方和侧方的马赫杆高度均不断下降,而后方的马赫杆高度不断增加。

    图  25  不同方向柱形装药的峰值压力随长径比变化
    Figure  25.  Peak pressure varying with length-to-diameter ratio for cylindrical charge in different directions
    表  8  不同长径比柱形装药的马赫杆高度
    Table  8.  Mach stem height of cylindrical charge with different aspect ratio
    长径比 马赫杆高度/mm
    1.0 ms2.0 ms3.0 ms
    前方后方侧方前方后方侧方前方后方侧方
    1∶1480150901 250 8705201 9701 5201 060
    1.5∶138017060 940 9904101 4801 740 890
    2∶131017050 7401 0503501 1801 850 790
    2.5∶126017040 6101 080310 9601 920 740
    下载: 导出CSV 
    | 显示表格

    对爆心高度、倾斜角度、动爆速度、装药长径比进行16组正交优化的数值模拟,对各组选取距爆心2.0、3.0、4.0 m处测点,对装药前后方向的峰值压力进行分析,并测量在柱形装药前方、后方马赫杆高度。装药前方、后方距爆心水平距离4.0 m处的近地峰值超压的各因素影响如图26~27所示。

    图  26  正交优化下对装药前方近地峰值压力的各因素影响
    Figure  26.  Influence of various factors of near ground peak pressure in front of the charge under orthogonal optimization
    图  27  正交优化下对装药后方近地峰值压力的各因素影响
    Figure  27.  Influence of various factors of near ground peak pressure at the rear of the charge under orthogonal optimization

    图26中动爆速度的峰值超压极差最大,平均极差为1.21 MPa。动爆速度越高,前侧的峰值压力越大,其次影响因素为装药长径比和爆心高度。图27中各因素对压力的影响趋势基本相同,其中倾斜角度的极差最大,对峰值压力的影响占据主导。其余因素较倾斜角度的变化幅度较小,动爆速度及装药长径比对后侧峰值压力基本呈现负影响的趋势,爆心高度对后侧峰值压力的影响不明显。

    由于建模尺寸的影响,对1.0 ms时刻装药前方马赫杆高度、后方马赫杆高度进行极差分析,结果如图28~29所示。

    图  28  正交优化下对装药前方马赫杆高度的各因素影响
    Figure  28.  Influence of various factors on Mach stem height in front of the charge under orthogonal optimization
    图  29  正交优化下对装药后方马赫杆高度的各因素影响
    Figure  29.  Influence of various factors on Mach stem height at the rear of the charge under orthogonal optimization

    从装药前方和装药后方马赫杆高度的极差可以看出,爆心越低,则马赫杆越高。倾斜角度是影响马赫杆高度的次要因素,且在30°附近取得最大值。动爆速度对前后高度呈现相反的趋势,这点与前面的结论相同。

    通过分析柱形装药近地爆炸中的4个参数:爆心高度、倾斜角度、动爆速度和装药长径比,得到了各参数对近地爆炸中各指标的影响规律。分析前、后、侧3个方向的近地冲击波传播规律及峰值压力和马赫杆高度,对各个方向的压力变化和马赫波传播进行了量化对比。综合4个变量的正交优化结果,各因素对不同的近地爆炸指标有较强的规律性,在柱形装药的近地爆炸研究中得到以下主要结论。

    (1) 柱形装药静爆时,随装药倾斜角的增大,前方的峰值压力先减小后增大,后方呈现下降,而侧方呈上升趋势。马赫杆高度在装药前后的比值先增大后减小,侧方的马赫杆逐渐增高。倾斜角度与装药长径比是影响静爆冲击波前后侧传播差异的主要因素,3 ms时对前、后、侧3个方向的马赫杆高度的最大提升幅度依次为29.7%、88.1%和96.4%。

    (2) 柱形装药动爆的趋势较一致,即动爆速度越高,装药前方和侧方的峰值压力及马赫杆高度越大,但对后方的影响较小。3.0 ms时,前方马赫杆高度最大增长幅度为156.2%,后方马赫杆高度最大幅度为76.2%,侧方马赫杆高度最大提升为388.2%。

    (3) 动爆速度对装药前方马赫波的增益最大,平均极差为1.21 MPa;倾斜角度对装药后方的马赫波的增益最大,平均极差为1.38 MPa;爆心高度则对前后侧的马赫杆高度影响最大,平均极差为1 051.25 mm。

  • 图  1  GP1不锈钢粉末扫描电子显微镜形貌

    Figure  1.  Scanning electron microscope morphology of GP1 stainless steel powder

    图  2  粒径分布图

    Figure  2.  Powder size distribution for GP1 stainless steel powder

    图  3  试样设计尺寸、成型后试样照片以及成型工艺过程示意图

    Figure  3.  Sample design dimension, sample photograph after moulding and schematic diagram of moulding process

    图  4  激光选区熔化法制备的GP1不锈钢光学显微结构

    Figure  4.  Optical images of GP1 stainless steel processed by selective laser melting

    图  5  激光选区熔化法制备的GP1不锈钢显微结构的电子背散射衍射表征

    Figure  5.  Electron backscattered diffraction images of GP1 stainless steel processed by selective laser melting

    图  6  HTM-5020高速拉伸材料试验机、试样夹具及试样散斑图

    Figure  6.  HTM-5020 high-speed machine, sample fixture, and speckle image of specimen

    图  7  应变片粘贴在GP1不锈钢试样的位置

    Figure  7.  Position of strain gauge pasted on GP1 stainless steel specimen

    图  8  层裂实验装置示意图

    Figure  8.  Schematics of the setup for spallation experiment

    图  9  闭环控制时力与速度时程曲线

    Figure  9.  Force and velocity profiles under closed loading control mode

    图  10  不同应变率下的加载速度时程曲线

    Figure  10.  Velocity profiles at different strain rates

    图  11  不同应变率下的力时程曲线

    Figure  11.  Force profiles at different strain rates

    图  12  不同时刻的GP1不锈钢试样标距段应变分布云图(应变率为10−2 s−1)

    Figure  12.  Strain distributions of GP1 stainless sample at different times (strain rate is 10−2 s−1)

    图  13  不同应变率下GP1不锈钢的真实应力-应变曲线

    Figure  13.  True stress-true strain curves of GP1 stainless steel at different strain rates

    图  14  GP1不锈钢中奥氏体和马氏体体积分数电子背散射衍射表征(黄色代表奥氏体,红色代表马氏体)

    Figure  14.  Electron backscattered diffraction characterizations of volume fraction for austenite and martensite in GP1 stainless steel (yellow represents austenite, red represents martensite)

    图  15  GP1不锈钢屈服应力随着应变率的变化曲线

    Figure  15.  Yield stress variation of GP1 stainless steel with strain rate

    图  16  不同初始速度下GP1不锈钢的自由面速度剖面

    Figure  16.  Free-surface velocity profiles of GP1 stainless steel at different initial velocities

    图  17  不同应变率下GP1不锈钢的宏观断口以及微观形貌

    Figure  17.  Macro-fracture and micro-morphology of GP1 stainless steel at different strain rates

    图  18  不同初始速度下的层裂剖面

    Figure  18.  Spall profiles at different initial velocities

    图  19  初始层裂的微观金相

    Figure  19.  Microscopic metallographic phase of initial spallation

    图  20  GP1不锈钢层裂断口微观形貌

    Figure  20.  Micrographs of ductile fractures in GP1 stainless steel spallation

    表  1  一维应变平板撞击实验结果

    Table  1.   Results of one-dimensional strain plane impact test

    实验编号初始速度/(m·s−1)Δu/(m·s−1)σHEL/GPaσy/GPaσs/GPa
    1250176.231.781.013.91
    2270174.971.791.023.88
    3350169.991.851.053.76
    下载: 导出CSV
  • [1] SAMES W J, LIST F A, PANALA S, et al. The metallurgy and processing science of metal additive manufacturing [J]. International Materials Reviews, 2016, 61(5): 1–46. DOI: 10.1080/09506608.2015.1116649.
    [2] ZHAI Y, LADOS D A, LAGOY J L. Additive manufacturing: making imagination the major limitation [J]. Journal of metals, 2014, 66(5): 808–816. DOI: 10.1007/s11837-014-0886-2.
    [3] YADOLLAHI A, SHAMSAEI N. Additive manufacturing of fatigue resistant materials: challenges and opportunities [J]. International Journal of Fatigue, 2017, 98(1): 14–31. DOI: 10.1016/j.ijfatigue.2017.01.001.
    [4] 王沛, 黄正华, 戚文军, 等. 基于SLM技术的3D打印工艺参数对316不锈钢组织缺陷的影响 [J]. 机械制造文摘: 焊接分册, 2016, 1(2): 2–7.

    WANG Pei, HUANG Zhenghua, QI Wenjun, et al. Effects of 3D printing process parameters based on SLM technology on structural defects of 316 stainless steel [J]. Mechanical Manufacturing Abstracts: Welding Brochures, 2016, 1(2): 2–7.
    [5] 吕豪, 杨志斌, 王鑫, 等. 激光增材制造GH4099合金热处理后的显微组织及拉伸性能 [J]. 中国激光, 2018, 45(10): 3–9. DOI: 10.3788/cjl.201845.1002003.

    LÜ Hao, YANG Zhibin, WANG Xin, et al. Microstructure and tensile properties of GH4099 alloy fabricated by laser additive manufacturing after heat treatment [J]. Chinese Journal of Lasers, 2018, 45(10): 3–9. DOI: 10.3788/cjl.201845.1002003.
    [6] 尹燕, 刘鹏宇. 路超., et al 选区激光熔化成形316L不锈钢微观组织及拉伸性能分析 [J]. 焊接学报, 2018, 39(8): 77–81. DOI: 10.12073/j.hjxb.2018390205.

    YIN Yan, LIU Pengyu, LU Chao, et al. Microstructure and tensile properties of selective laser melting forming 316L stainless steel [J]. Transactions of the China Welding Institution, 2018, 39(8): 77–81. DOI: 10.12073/j.hjxb.2018390205.
    [7] 王志会, 王华明, 刘栋. 激光增材制造AF1410超高强度钢组织与力学性能研究 [J]. 中国激光, 2016, 43(4): 59–65. DOI: 10.3788/CJL201643.0403001.

    WANG Zhihui, WANG Huaming, LIU Dong. Microstructure and mechanical properties of AF1410 ultra-high strength steel using laser additive manufacture technique [J]. Chinese Journal of Lasers, 2016, 43(4): 59–65. DOI: 10.3788/CJL201643.0403001.
    [8] YADOLLAHI A, SHAMSAEI N, THOMPSONS M, et al. Effects of building orientation and heat treatment on fatigue behavior of selective laser melted 17-4 PH stainless steel [J]. International Journal of Fatigue, 2017, 94(11): 218–235. DOI: 10.1016/j.ijfatigue.2016.03.014.
    [9] SURYAWANSHI J. Mechanical behavior of selective laser melted 316L stainless steels [J]. Materials Science and Engineering: A, 2017, 696(7): 113–121. DOI: 10.1016/j.msea.2017.04.058.
    [10] WANG Y M, VOISIN T, MCKEOWN J T, et al. Additively manufactured hierarchical stainless steels with high strength and ductility [J]. Nature Materials, 2017, 17(1): 63–71. DOI: 10.1038/nmat5021.
    [11] YU S, HEBERT R J, MARK A. Effect of heat treatments on microstructural evolution of additively manufactured and wrought 17-4PH stainless steel [J]. Materials and Design, 2018, 156(10): 429–440. DOI: 10.1016/j.matdes.2018.07.015.
    [12] GRAY G T, LIVESCU V, RIGG P A, et al. Structure/property (constitutive and spallation response) of additive manufactured 316L stainless steel [J]. Acta Materialia, 2017, 138(10): 140–149. DOI: 10.1016/j.actamat.2017.07.045.
    [13] SONG B, NISHIDA E, SANBORN B, et al. Compressive and tensile stress-strain responses of additively manufactured (AM) 304L stainless steel at high strain rates [J]. Journal of Dynamic Behavior of Materials, 2017, 3(3): 412–425. DOI: 10.1007/s40870-017-0122-6.
    [14] BRANDON M W, BRAHMANNADA P, ANDELLE K, et al. High strain rate compressive deformation behavior of an additively manufactured stainless steel [J]. Additive Manufacturing, 2018, 9(24): 432–439. DOI: 1010.1016/j.addma.2018.09.16.
    [15] 丁利, 李怀学, 王玉岱, 等. 热处理对激光选区熔化成形316不锈钢组织与拉伸性能的影响 [J]. 中国激光, 2015, 42(4): 187–193. DOI: 10.3788/CJL201542.0406003.

    DING Li, LI Huaixue, WANG Yudai, et al. Heat treatment on microstructure and tensile strength of 316 stainless steel by selective laser melting [J]. Chinese Journal of Lasers, 2015, 42(4): 187–193. DOI: 10.3788/CJL201542.0406003.
    [16] WOOD P, SCHLEY C A, WILLIAMS M A, et al. A method to calibrate a specimen with strain gauges to measure force over the full-force range in high rate testing [C] // SCHLEY C A. DYMAT 2009: 9th International Conference on the Mechanical and Physical Behaviour of Materials under Dynamic Loading. Belgium: Experimental Techniques?, 2009: 265−273. DOI: 10.1051/dymat/2009036.
    [17] 申海艇, 蒋招绣, 王贝壳, 等. 基于超高速相机的数字图像相关性全场应变分析在SHTB实验中的应用 [J]. 爆炸与冲击, 2017, 37(1): 15–20. DOI: 10.11883/1001-1455(2017)01-0015-06.

    SHEN Haiting, JIANG Zhaoxiu, WANG Beike, et al. Full field strain measurement in split Hopkinson tension bar experiments by using ultra-high-speed camera with digital image correlation [J]. Explosion and Shock Waves, 2017, 37(1): 15–20. DOI: 10.11883/1001-1455(2017)01-0015-06.
    [18] PIERRON F, SUTTON M A, TIWARI V. Ultra high speed DIC and virtual fields method analysis of a three point bending impact test on an aluminum bar [J]. Experimental Mechanics, 2011, 51(4): 537–563. DOI: 10.1007/s11340-010-9402-y.
    [19] 王楠, 李恩普, 汤忠斌, 等. 二维数字图像相关方法的拉伸实验误差分析 [J]. 光学仪器, 2012, 34(3): 5–12. DOI: 10.3969/j.issn.1005-5630.2012.03.002.

    WANG Nan, LI Enpu, TANG Zhongbin, et al. An investigation of the experimental error of 2-D DIC method applied to tensile strain measurement [J]. Optical Instruments, 2012, 34(3): 5–12. DOI: 10.3969/j.issn.1005-5630.2012.03.002.
    [20] CHEVRIER P, KLEPACZKO J R. Spall fracture: mechanical and micro-structural aspects [J]. Engineering Fracture Mechanics, 1999, 63(3): 273–294. DOI: 10.1016/S0013-7944(99)00022-3.
    [21] 张万甲, 曾元金. 不锈钢(00Cr18Ni9)动态累积损伤研究 [J]. 爆炸与冲击, 1999, 19(4): 309–314. doi: 10.3321/j.issn:1001-1455.1999.04.004

    ZHANG Wanjia, ZENG Yuanjin. Study on the dynamic accumulation-damage for the stainless steel (00Cr18Ni9) [J]. Explosion and Shock Waves, 1999, 19(4): 309–314. doi: 10.3321/j.issn:1001-1455.1999.04.004
    [22] WENG J, TAN H, WANG X, et al. Optical-fiber interferometer for velocity measurements with picosecond resolution [J]. Applied Physics Letters, 2006, 89(11): 111101-0. DOI: 10.1063/1.2335948.
    [23] CLAUSEN B, BROWN D W, CARPENTER J S, et al. Deformation behavior of additively manufactured GP1 stainless steel [J]. Materials Science and Engineering: A, 2017, 696(4): 331–340. DOI: 10.1016/j.msea.2017.04.081.
    [24] 刘超, 王磊, 刘杨. 应变速率对Q&P钢拉伸变形行为的影响 [J]. 特钢技术, 2012, 18(3): 18–22. DOI: 10.3969/j.issn.1674-0971.2012.03.007.

    LIU Chao, WANG Lei, LIU Yang. Effect of strain rates on tensile deformation behavior of quenching and partitioning steel [J]. Special Steel Technology, 2012, 18(3): 18–22. DOI: 10.3969/j.issn.1674-0971.2012.03.007.
    [25] COWPER G R, SYMONDS P S. Strain hardening and strain rate effects in the impact loading of cantilever beams [R] // Applied Mathematics Report. Brown University, 1958.
    [26] 蒋招绣, 辛铭之, 王永刚. 高强铝合金的动态拉伸断裂行为实验研究 [J]. 固体力学学报, 2014, 35(6): 552–558. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2014.06.007.

    JIANG Zhaoxiu, XIN Mingzhi, WANG Yonggang. Experimental study on dynamic tensile fracture of aluminum alloy [J]. Chinese Journal of Solid Mechanics, 2014, 35(6): 552–558. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2014.06.007.
    [27] 彭辉, 李平, 裴晓阳, 等. 平面冲击下铜的拉伸应变率相关特性研究 [J]. 物理学报, 2014, 63(19): 281–287. DOI: 10.7498/aps.63.196202.

    PENG Hui, LI Ping, PEI Xiaoyang, et al. Rate-dependent characteristics of copper under plate impact [J]. Acta Physica Sinica, 2014, 63(19): 281–287. DOI: 10.7498/aps.63.196202.
    [28] ANTOUN T, SEAMAN L, CURRAN D R, et al. Spall fracture [M]. New York, USA: Springer, 2003: 90−92. DOI: 10.1007/b97226.
  • 期刊类型引用(1)

    1. 王明涛,程月华,吴昊. 柱形装药空中运动爆炸冲击波荷载计算模型. 爆炸与冲击. 2024(07): 107-125 . 本站查看

    其他类型引用(2)

  • 加载中
图(20) / 表(1)
计量
  • 文章访问数:  5856
  • HTML全文浏览量:  1931
  • PDF下载量:  74
  • 被引次数: 3
出版历程
  • 收稿日期:  2019-01-16
  • 修回日期:  2019-03-25
  • 网络出版日期:  2019-06-25
  • 刊出日期:  2019-07-01

目录

/

返回文章
返回