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

复杂边界条件下近场水下爆炸对固支方板的毁伤效应

邓硕 赖志超 秦健 孟祥尧 迟卉 黄瑞源

马彦, 袁福平, 武晓雷. 双相高强钢FeNiAlC的动态剪切行为及微结构机理[J]. 爆炸与冲击, 2021, 41(1): 011404. doi: 10.11883/bzycj-2020-0224
引用本文: 邓硕, 赖志超, 秦健, 孟祥尧, 迟卉, 黄瑞源. 复杂边界条件下近场水下爆炸对固支方板的毁伤效应[J]. 爆炸与冲击, 2023, 43(11): 112204. doi: 10.11883/bzycj-2023-0164
MA Yan, YUAN Fuping, WU Xiaolei. Dynamic shear behaviors and microstructural deformation mechanisms in FeNiAlC dual-phase high strength alloy[J]. Explosion And Shock Waves, 2021, 41(1): 011404. doi: 10.11883/bzycj-2020-0224
Citation: DENG Shuo, LAI Zhichao, QIN Jian, MENG Xiangyao, CHI Hui, HUANG Ruiyuan. Damage effects of clamped square plates by near-field underwater explosion with complex boundary conditions[J]. Explosion And Shock Waves, 2023, 43(11): 112204. doi: 10.11883/bzycj-2023-0164

复杂边界条件下近场水下爆炸对固支方板的毁伤效应

doi: 10.11883/bzycj-2023-0164
基金项目: 国家自然科学基金(12172178)
详细信息
    作者简介:

    邓 硕(1999- ),男,硕士研究生,shuodeng1999@163.com

    通讯作者:

    黄瑞源(1984- ),男,博士,研究员,huangruiyuan1984@163.com

  • 中图分类号: O383.1

Damage effects of clamped square plates by near-field underwater explosion with complex boundary conditions

  • 摘要: 近场水下爆炸会产生复杂的载荷模式,而复杂的边界条件使结构在近场水下爆炸作用下的毁伤形态更加难以预测。因此,采用耦合的欧拉-拉格朗日算法探究了水下爆炸气泡在多边界耦合作用下(自由面、弹塑性板、泥沙边界)的演化过程及其对固支方板的毁伤效应。首先,开展了2.5 g TNT在不同尺寸(板边长为0.46、0.92和1.61倍最大气泡直径)固支方板底部10 cm起爆的水下爆炸试验,验证了有限元方法的准确性。然后,结合试验和有限元结果分析了不同边界条件下板的毁伤机理。最后,通过系列数值模拟发现:随着板尺寸和爆距的增大,气泡会出现溃散、下射流和上射流3种不同的演化方式;随着板尺寸的增大,爆距对板中心最终变形的影响减小;泥沙边界能减缓气泡收缩,使气泡从中部塌陷形成方向相反的对射流,降低固支方板的位移和应变,对于气泡提前溃散的工况,泥沙边界基本无影响。
  • 高应变速率载荷下,金属材料往往因局部化形变而失效破坏[1]。在军事工业领域,子弹的穿甲侵彻与材料的局域化绝热剪切损伤紧密相关。对实心穿甲弹,弹头的绝热剪切有利于在侵彻过程中自锐化而增大穿刺深度[2]。在航空航天领域,飞行器的安全威胁主要来自空间碎片的高速撞击,材料的绝热剪切直接导致飞行器构件的灾难性破坏。所以,金属材料的动态剪切行为是材料在特殊使用环境下必须考虑的关键问题[1-3]

    高应变速率载荷诱发材料的局部化形变,有两个基本特征:加载速度快变形过程接近绝热状态;变形区域小易形成剪切带。由于绝热剪切带(adiabatic shear band, ASB)的形成和扩展在约105 s内完成且宽度一般在102 μm内,目前仍缺乏有效的实验技术实现原位直接观测[2]。现阶段,主要手段是通过表征材料变形失效后的微观结构来分析和推演ASB形成和扩展过程。Meyers等[4]在低层错能304L不锈钢动态剪切实验中,发现ASB内形成晶粒尺寸100~200 nm的等轴晶,通过形变织构分析,得出这些等轴晶是通过变形过程中初始晶粒的旋转、滑移和孪生机制形成的。Yang等[5]则通过应变冻结实验,系统地研究1Cr18Ni9Ti奥氏体型不锈钢形成ASB过程中微结构的演化,证实ASB内晶粒的动态旋转再结晶机制。Li等[6]在高锰钢的子弹侵彻试验中,发现ASB内产生明显的纳米晶区域和非晶区域,表明在ASB内由于极高的温升产生了金属熔化和快速冷却。Wei等[7-8]则分别研究了粗晶(coarse grained, CG)、超细晶(ultrafine grained, UFG)、纳米晶(nanocrystalline, NC)结构铁在高应变率下的压缩变形,结果表明,减小晶粒尺寸会导致材料形成ASB的概率提高,同时定量分析了应变硬化率和率敏感性对形成ASB的影响。

    目前,关于金属材料ASB失效的基本准则是由Zener等[9]提出的材料应变硬化与热软化相互竞争的机制,即当材料丧失应变硬化能力时发生ASB失效。Batra等[10]提出ASB形成判据(最大应力准则:ASB在流变应力达到最大值时形成并快速扩展,导致流变应力快速跌落至峰值的80%),但该准则往往仅适用于均质材料。在梯度材料[11]和多层材料[12]中,ASB形成于材料中硬化能力较弱的区域,但由于材料的非均匀性,当局部形成ASB后材料整体依然具有硬化能力,当ASB贯穿整个样品时材料才会丧失承载能力。在这些异构材料[11-12]中,局部微区的力学属性差异能够延缓ASB的形成和传播,从而延缓材料失效的发生。由此可见,材料的微观结构对形成ASB失效非常重要。本文中,以FeNiAlC双相钢为例,讨论不同微结构材料的剪切变形行为,进一步理解微结构对材料形成ASB的影响以及温升和热软化对材料发生局部化变形的影响。

    FeNiAlC合金由纯度为99.99%的铁、镍、铝金属单质进行熔炼而成。首先,在真空电弧炉中经过5次反复熔炼形成20 mm的铸锭,测得铸锭中铁以外化学成分的质量分数为w(Ni)=24.86%, w(Al)=5.8%, w(C)=0.38%,记为Fe-24.86Ni-5.8Al-0.38C。接着,将铸锭加热至1 100 ℃后热锻成10 mm×10 mm×50 mm的钢板,在950 ℃下进行热轧,终轧厚度为6 mm,冷却至室温。然后,将热轧板在真空马弗炉中进行均匀化固熔,固熔温度为1 150 ℃,固熔24 h后进行快速淬火处理。最后,将固熔后的金属板进行大应变冷轧处理,压下量(截面积减少量)分别为39%、53%(记为CR39、CR53),未进行冷轧处理的材料为粗晶态(记为CG)。

    动态冲击样品采用强迫剪切构形试样(帽形试样),强迫剪切区宽度为200 μm,剪切方向垂直于板材轧制方向,帽形试样的厚度即为冷轧后板材的厚度,如图1所示。实验时,将试样紧密放置于用马氏体时效钢制成的限位环中,限位环在冲击过程中变形量极小,保证样品主要受剪应力的作用。由有限元模拟结果[13],剪切区的高度(2 mm)远大于剪切区的宽度(0.2 mm),压应力相对剪应力是小量,因此帽形样品的变形过程可以近似等效为纯剪切变形。在13 mm分离式霍普金森压杆(split Hopkinson pressure bar, SHPB)装置上进行冲击实验,子弹速度约20 m/s,加载应变率范围为10 000~30 000 s−1,每个状态的样品重复5次实验,重复性相对较好,实验误差小。通过杆上的应变片和动态应变仪测量,得到入射杆和透射杆中的反射波应变εR和透射波应变εT,进行空载标定应变片敏感系数。利用两波法和一维弹性理论,计算样品在一维应力加载下的剪切应力τs、剪切应变γs、应变率˙γs与和位移us

    图  1  动态剪切实验装置和样品
    Figure  1.  Dynamic shear experimental device and its sample
    {τs=E(A/As)εTus=2c0t0εRdtγs=Us/ts˙γs=2c0tsεR
    (1)

    式中:EAc0分别为霍普金森杆的杨氏模量(210 GPa)、横截面积(132.7 mm2)和纵波波速(5 050 m/s),Asts分别为帽形样品的剪切区面积和宽度[14]

    实验前后样品的微观结构采用Zeiss Gemini 300扫描电镜(scanning electron microscope, SEM)和JEOL 2100F透射电镜(transmission electron microscope, TEM)进行表征。使用电子背散射衍射(electron backscatter diffraction, EBSD)技术获取材料的晶体学取向及位错密度等,选用Oxford Symmetry探头进行扫描,扫描步长为50~1 000 nm,使用CHANNEL 5软件进行数据后处理。观测前的样品准备,详见文献[15]。

    冲击实验前,对CG、CR53样品进行微结构表征。图2为CG样品晶粒的反极图(inverse pole figure, IPF)和相图(phase figure, PF)。固熔后样品由奥氏体相(austenite, γ)和马氏体相(martensite, α)构成,奥氏体平均晶粒尺寸较大,约为11.5 μm,马氏体相析出在奥氏体晶界处,平均晶粒尺寸约为3.55 μm,马氏体相约占15%,属于一种双相合金。由图2(a)可见,固熔后晶粒取向随机分布,无明显织构。图3为CR53样品的IPF和PF。轧制过程中发生少量马氏体相变,马氏体相比例略有增加,约占18%。奥氏体晶粒被拉长,统计晶粒的长轴和短轴比约为2.1,据晶粒被拉长计算等效应变约为0.49[16],而冷轧53%换算成等效应变约为0.50。这说明冷轧过程中,变形主要发生在奥氏体相中,马氏体相几乎不承担变形。由图3(a)可见,与粗晶态相比,冷轧后出现明显的〈111〉丝织构[17]。为了研究样品冷轧前后的位错形貌和位错密度,对CG、CR53样品进行TEM观察。图4(a)为CG样品奥氏体相的形貌和[011]晶向衍射谱。固熔后样品内部有少量的平形位错,长度较短,位错密度极低。图4(b)为CR53样品的晶粒衍衬像、α相[001]晶向的衍射谱和γ相[011]晶向的孪晶衍射谱,样品内部位错密度明显增大,位错大多缠结在晶界和相界附近,奥氏体晶粒内部产生少量变形孪晶。CG和CR53样品中晶粒形貌、尺寸和位错密度的差异,影响两种微结构样品的动态剪切行为,这将在以下作详细介绍。

    图  2  实验前CG样品的微观结构
    Figure  2.  Microstructures of CG sample before experiment
    图  3  实验前CR53样品的微观结构
    Figure  3.  Microstructures of CR53 sample before experiment
    图  4  实验前CG和CR53样品的TEM
    Figure  4.  TEM observations of CG and CR53 sample before experiment

    用动态剪切实验所得的数据,通过式(1)获得的剪切应力与变形位移、剪切应变的关系曲线如图5所示。圆点为名义动态剪切屈服强度,CG样品约为480 MPa,CR39、CR53样品分别为910、1 250 MPa,冷轧量越大,样品的剪切屈服强度越高。方点为样品的最大剪切应力,根据最大剪应力准则[10],此时样品剪切区内部形成ASB,此后剪切应力快速下降,样品发生失稳破坏。在形成剪切带前,变形主要集中在宽200 μm的剪切区内,一般认为这个过程中剪切区内的变形是均匀的[14]图5(b)中,曲线所包围的面积可以近似为材料的动态剪切韧性。

    图  5  不同微结构的动态剪切性能
    Figure  5.  Dynamic shear properties of various microstructures

    与其他合金一样,利用帽形试样获得FeNiAlC合金的动态剪切性能,并进行比较。图6为FeNiAlC合金和其他常用合金(301、316不锈钢等)的均匀剪切应变和剪切韧性与动态剪切强度的关系[13-14, 18-24]。在同等动态剪切强度情况下,FeNiAlC合金剪切韧性优于其他合金,拥有较优异的动态剪切性能。

    图  6  不同金属的动态剪切性能
    Figure  6.  Dynamic shear properties of various metals

    为了研究变形过程中微结构的变化,采用SEM中的EBSD技术对变形后的样品进行表征。由图7(a)(b)可见,CG样品中剪切区的中心有一条识别率较低的窄带(两条虚线之间)是剪切失稳后形成的ASB。ASB两侧约100 μm区域是主要的剪切变形区,剪切变形区内产生大量小角晶界,在奥氏体和马氏体的相界附近有大量的小晶粒产生。在等轴的奥氏体晶粒内形成大量的变形孪晶(deformation twin, DT),宽度为2~3 μm。由图7(b)可见,马氏体相所占比例约14%,与剪切变形前相比几乎没有变化,说明变形过程中没有发生马氏体相变。在图7(c)中,红线为该点与相邻点之间的取向差,蓝线为该点与初始点的取向差。每跨过一个(孪晶)界面时,取向角改变60°,这是典型的FCC晶体Σ3孪晶取向关系[25]。CG样品剪切变形后孪晶密度较高,发生孪生变形的区域占总体积的60%以上,说明在此应力状态下,奥氏体晶粒大部分发生孪生变形。

    图  7  实验后CG样品的微观结构
    Figure  7.  Microstructures of CG sample after experiment

    图8(a)(b)可见,CR53样品中靠近ASB区域产生大量的小晶粒,与变形前相比晶粒明显细化,远离ASB区域产生少量的孪晶。由图8(b)可见,α相所占比例约18%,与变形前相近,说明冷轧态样品也几乎没有发生马氏体相变。由图8(c)可见,当跨过孪晶界面时,CR53样品界面两侧取向角相差60°,同CG样品结果类似。

    图  8  实验后CR53样品的微观结构
    Figure  8.  Microstructures of CR53 sample after experiment

    KAM (kernel average misorientation)即mka,为晶粒内一点与周围相邻点的取向差的平均值,是一个与几何必需位错(geometrically necessary dislocation, GND)密度正相关的物理量[26]图910分别为CG和CR53样品剪切变形前后局部晶粒的KAM及分布。分布图中,统计分布时取向差步长为0.25°。对比CG样品剪切变形前后局部晶粒的KAM,可以得到两个有价值的结论。(1)变形前,由于经过高温长时间的固熔处理,晶粒内部取向差极小,说明此时晶粒内部位错密度极低,仅在晶界周围小于1 μm内存在明显的取向差分布(0.5°~1.1°)。这可能是因固熔后进行快速淬火处理,在界面处会产生残余应变,引起GND密度升高。(2)与变形前相比,剪切变形后奥氏体晶粒承担主要变形,因此在奥氏体晶粒内部出现了高密度的GND,KAM在0.8°~2.0°,尤其在孪晶界周围,KAM约为2°。说明变形产生大量的位错钉扎在孪晶界附近,位错与孪晶的相互作用为CG样品在剪切变形过程中提供了较大的应变硬化能力。对比CG态样品变形前后γ相KAM分布可见,由于样品中α相含量较少,硬度较大[15],几乎不参与变形,变形前后KAM几乎没有变化,所以这里只进行γ相KAM的比较。变形前KAM平均值约为0.96°,变形后为1.26°,根据KAM与GND密度的关系[26],可以估计变形后GND增加约30%。这也说明,CG样品均匀剪应变之所以能达到5,晶粒内部有充足的可动位错参与变形。

    图  9  实验前CG和CR53样品的KAM及分布
    Figure  9.  KAM values and distributions of CG and CR53 sample before experiment
    图  10  实验后CG和CR53样品的KAM及分布
    Figure  10.  KAM values and distributions of CG and CR53 sample after experiment

    对比CR53样品剪切变形前后局部晶粒的KAM可见,由于经过大应变冷轧过程(等效应变约0.50),样品表现出很高的KAM,尤其是在晶界附近KAM为1.8°~3.0°,远高于固熔后的CG样品。不同于CG样品,CR53样品变形后晶粒内部的KAM明显小于变形前的KAM。同样,对比CR53样品变形前后γ相KAM分布可见,变形前分布范围为0.8°~3.1°,变形后减小为0.8°~2.3°,变形前平均值为1.73°,变形后减小至1.33°,根据KAM与GND密度的线性关系,GND密度减小约23%。通常,变形会使位错密度增大,但由于CR53样品变形前的位错密度已经很高,且动态实验应变率高,变形时间短,剪切区内发生绝热剪切,有较大的温升,会使局部晶粒发生部分回复。同时,该样品均匀剪切应变较小,约为1,所以剪切变形后,回复的作用大于变形的效果,KAM不增反降,GND密度降低,致使CR样品应变硬化能力降低,过早的发生塑性失稳(形成ASB)。

    在动态加载过程中,由于加载速度极高,变形时间在100 μs内,热传导来不及进行,所以整个剪切过程可近似成绝热过程[1-3]。同时,认为在均匀剪切阶段,剪切区内的应变是均匀的。计算剪切过程中剪切区温升的塑性功转化温升公式为[14]

    ΔT=ηρcVγ0τdγp
    (2)

    式中:η为总的非弹性功转化为热量的系数(Taylor-Quinney系数),本文中取0.8;ρcV为材料的密度和比热容,τγp为剪切应力和剪切塑性应变。对FeNiAlC合金,测得ρ为9.8 g/cm3cV取540 J/(kg·K)[27]图11为用实验所测的剪切应力和应变、据式(2)计算的温升曲线。CR样品的温升斜率明显高于CG样品,这是因为CR样品由于冷轧过程中引入大量的可动位错及位错胞等亚结构,在塑性变形初期剪切流变应力大,导致温升的斜率高于CG样品。当剪应变达到0.5时,CR样品的温升约为100 K,CG样品只有50 K,而当剪应变达到最大均匀剪应变时,CR53样品最终温升达到190 K,而CG样品最终温升高达720 K。这表明,绝热温升可能不是ASB形成的主导或者唯一因素[27-28]

    图  11  最大应力点前剪切区塑性功引起的温度升高
    Figure  11.  Temperature rise due to plastic dissipation work in shear zone before maximum stress point

    一般认为,ASB是材料受到冲击载荷作用时的一种失效模式[1-2]。当材料受到高速率剪切力时,通常认为材料的失效是塑性应变硬化和绝热温升所致的热软化相互竞争的结果[9]。变形初期,温升较小,材料的应变硬化占据主导,发生均匀变形,剪切应力不断提高。当剪切区温升达到一定值时,材料的应变硬化不足以抵抗温升带来的热软化,材料发生局部化变形,剪切应力骤然下降,形成ASB。图78中未能解析的带状部分就是变形形成的ASB。因为带内发生了比带外高1~2个数量级的变形,所以当前的扫描步长不足以解析其微结构。对于ASB内的微结构采用40 nm的步长重新扫描,如图12所示。CG样品ASB内晶粒不同于剪切区的晶粒发生孪生变形或局部晶粒细化,而是形成了大小几乎相等的等轴晶。晶粒尺寸集中在200~500 nm,远远小于初始的晶粒尺寸(13 μm)。由此可得,在ASB形成时带内组织发生了动态再结晶(dynamic recrystallization, DRX),与大多数金属材料类似[14],在大应变和高温的作用下,ASB内形成取向随机、大小相等的再结晶晶粒。同CG样品相似,CR53样品同样发生了DRX。CG和CR53两种样品ASB内虽然应变极大,但均不发生相变,只是在原有马氏体相的基础上发生晶粒细化,没有形成新的马氏体相。

    图  12  绝热剪切带内微观结构
    Figure  12.  Microstructures within ASBs

    根据热塑模型理论[1],材料发生ASB失效时,剪切带的宽度可以通过材料所受最大剪应力和应变率等进行估算。Dodd等[1]提出剪切带宽度的计算公式:

    δ=2(λTτ˙γ)1/2
    (3)

    式中:τ˙γ分别为形成剪切带时的剪应力和剪应变率,T为试样温度,λ为热传导系数(不锈钢约为15 W/(m·K)[28]。基于前面的实验和计算结果,CG样品剪切带的理论宽度为29.7 μm(T=1 018 K,τ=1 100 MPa,˙γ=63 000 s−1),CR53样品的为14.6 μm (T=468 K,τ=1 340 MPa,˙γ=98 000 s−1)。图78中的剪切带宽度并不是均匀的,受到局部微结构和受力状态的影响较大。通过划线法统计剪切带的宽度,结果如图13所示。CG样品剪切带的宽度范围为12~19 μm,平均值为15.8 μm,而CR53样品剪切带宽度范围为9~15 μm,平均值为12.3 μm。对比来看,CR53样品的理论计算ASB宽度更接近于实验所测真实宽度,而CG样品的实际剪切带宽度范围小于理论计算值。这可能是根据塑性功计算样品温度带来的影响。热塑模型的基础是理想的绝热条件,而实际实验过程中:CR53样品从加载到ASB失效过程较短,更接近绝热条件;CG样品具有较好的应变硬化能力,因此从加载至产生ASB失效经历时间更长,考虑到热耗散因素,实际的剪切带宽度小于理论值。

    图  13  绝热剪切带宽度分布
    Figure  13.  ASBs width distributions

    在FeNiAlC合金准静态拉伸变形的研究[15]中,CG样品的主要变形方式为位错滑移和少量孪生变形,硬化能力较弱。而在本文中,由于应变率提高约108倍,面心立方结构一般具有正的率敏感系数,动态条件下会显著提升流变应力水平,极大地促进孪晶的形成[29-30],变形方式以孪晶与位错交互作用为主导,因此剪切韧性得到充分提升。在准静态拉伸下,CR样品变形方式以位错滑移和马氏体相变为主,相变诱导塑性(transformation induced plasticity,TRIP)效应显著,所以保持了较高的均匀塑性。但在动态剪切加载下,由于应变速率的提高,局部温升导致材料层错能升高,马氏体相变被抑制[31],因此应变硬化能力不足导致过早地发生局部失稳。同时,CR样品冷轧后,由于存在明显的织构不易于发生孪生变形[17],所以在动态剪切时,屈服后硬化能力不足,较早地发生失稳。

    ASB的萌生、绝热温升和应力骤降三者的因果关系,一直是材料绝热剪切变形研究的热点问题之一[3, 27]。在热塑模型理论中,认为绝热温升是形成ASB失效的主要原因,最大应力准则也是建立在温升导致的ASB形成后产生应力的骤然下降上。但是,由于动态冲击的变形特点(高应变率和极短时间),ASB范围极窄,形成速度极快与载荷下降时间相差无几,所以一直没有明确的实验证据证实三者的先后关系。Guo等[27-28]试图用实验证实材料剪切带失效与温升的关系,结果表明:在Ti合金中ASB的萌生先于塑性变形产生明显温升,而材料的最大应力也出现在ASB萌生前。该结果与传统的热塑模型所给的预测结论相悖。本文中关于绝热剪切带热塑模型分析计算的结果与实验结果对比显示,当剪切变形时间较长时(材料的硬化能力较强),变形过程与热塑模型的理论前提(理想绝热条件)有所差异,材料的热耗散和微观结构变形不均匀等原因会影响ASB的形成和扩展,使热塑模型的结论与真实实验结果相差较大。同时,本文中也证实了绝热温升不是双相钢材料形成ASB的主导因素或唯一因素,还应考虑共同影响材料ASB失效过程的局部化应变和材料微观结构演化等。

    对FeNiAlC合金进行帽形冲击实验,研究其动态剪切变形行为和微结构机理。实验主要对均匀化处理的粗晶结构样品和冷轧后的层状结构样品进行强迫剪切,研究不同结构样品均匀变形过程和发生ASB失效前后的微观结构变化,利用热塑模型对ASB宽度和塑性功的热转化温升进行分析,主要结论如下。

    (1)FeNiAlC双相钢拥有较优异的动态剪切性能,剪切强度为400~1 300 MPa,剪切应变为1.5~5.1。同等剪切强度的剪切韧性均高于传统的TWIP钢、304、316不锈钢等合金材料。

    (2)孪生与位错主导剪切变形过程中,CG样品在变形过程中发生大量的孪生变形,位错密度显著提高,位错与孪晶界的交互作用使材料获得良好的硬化能力,提高了剪切韧性。CR样品因在剪切变形前冷轧引入大量的位错,变形初期位错密度极高,屈服强度达1.3 GPa,同时因冷轧织构的存在抑制孪生变形;变形过程中,绝热温升使位错发生部分回复,材料的应变硬化能力不及CG样品,屈服后在较小的均匀应变下形成绝热剪切带,发生失稳。动态剪切下,材料的马氏体相变被抑制。

    (3)绝热剪切带内发生动态再结晶,形成晶粒取向随机的超细晶晶粒,同时因剪切带内极高的温度,带内不发生马氏体相变。

    (4)热塑模型在应变较小的情况能够与实验结果高度吻合;绝热温升不是ASB形成的主导因素或唯一因素,ASB的形成可能由局部化变形、材料的微观结构演化等影响。

  • 图  1  耦合欧拉-拉格朗日方法的求解过程

    Figure  1.  The solution procedure of the coupled Eulerian-Lagrangian method

    图  2  场地布置及试件设计

    Figure  2.  Site layout and specimen design

    图  3  有限元模型

    Figure  3.  Finite element models

    图  4  不同网格尺寸下距药包中心不同距离处的冲击波峰值压力对比

    Figure  4.  Comparison of peak shock wave pressures at different distances from the charge center under different element sizes

    图  5  TNT在20 cm×20 cm板底起爆引起的气泡演化

    Figure  5.  Evolution images of bubbles caused by the explosion of a TNT charge at the bottom of the 20 cm×20 cm plate

    图  6  TNT在40 cm×40 cm板底起爆引起的气泡演化

    Figure  6.  Evolution images of bubbles caused by explosion of the TNT charge at the bottom of the 40 cm×40 cm plate

    图  7  TNT在70 cm×70 cm板底起爆引起的气泡演化

    Figure  7.  Evolution images of bubbles caused by the explosion of a TNT charge at the bottom of the 70 cm×70 cm plate

    图  8  气泡半径-时程曲线

    Figure  8.  The time-history curve of equivalent bubble radius

    图  9  20 cm×20 cm板的毁伤模式

    Figure  9.  Damage mode of the 20 cm×20 cm plate

    图  10  40 cm×40 cm板的毁伤模式

    Figure  10.  Damage mode of the 40 cm×40 cm plate

    图  11  70 cm×70 cm板的毁伤模式

    Figure  11.  Damage mode of the 70 cm×70 cm plate

    图  12  参数示意图

    Figure  12.  Schematic diagram of parameters

    图  13  下射流的形成过程

    Figure  13.  The formation process of the downward jet

    图  14  水射流的不同形态

    Figure  14.  Different forms of water jets

    图  15  位移时程曲线

    Figure  15.  Displacement-time curves

    图  16  固支方板的最终毁伤特征曲线

    Figure  16.  Final damage characteristic curves of clamped square plates

    图  17  气泡的演化过程(μ=0.46, da=10 cm)

    Figure  17.  The evolution process of bubbles (μ=0.46, da=10 cm)

    图  18  气泡的演化过程(μ=1.38, da=1 cm)

    Figure  18.  The evolution process of bubbles (μ=1.38, da=1 cm)

    图  19  吃水深度对板动态响应的影响

    Figure  19.  Effects of draught depth on the dynamic responses of the plates

    图  20  泥沙边界对气泡演化的影响(μ=0.46)

    Figure  20.  Effects of sediment boundaries on bubble evolution (μ=0.46)

    图  21  泥沙边界对气泡演化的影响(μ=1.38)

    Figure  21.  Effects of sediment boundaries on bubble evolution (μ=1.38)

    图  22  泥沙边界对射流形成的影响

    Figure  22.  Effects of sediment boundaries on jet formation

    图  23  泥沙边界对板动态响应的影响

    Figure  23.  Effects of sediment boundaries on dynamic responses of plates

    表  1  试验工况

    Table  1.   Test conditions

    工况药包质量/g爆距/cmL/cml/cm厚度/mm
    12.51070602
    22.51040302
    32.51020102
    下载: 导出CSV

    表  2  数值模拟工况设置

    Table  2.   Numerical simulation condition settings

    工况 方板尺寸 μ 方板爆距/cm γ
    1 20 cm×20 cm 0.46 10 0.46
    2 15 0.69
    3 20 0.92
    4 25 1.15
    5 30 1.38
    6 40 cm×40 cm 0.92 10 0.46
    7 15 0.69
    8 20 0.92
    9 25 1.15
    10 30 1.38
    11 50 cm×50 cm 1.15 10 0.46
    12 15 0.69
    13 20 0.92
    14 25 1.15
    15 30 1.38
    16 60 cm×60 cm 1.38 10 0.46
    17 15 0.69
    18 20 0.92
    19 25 1.15
    20 30 1.38
    21 70 cm×70 cm 1.61 10 0.46
    22 15 0.69
    23 20 0.92
    24 25 1.15
    25 30 1.38
    下载: 导出CSV

    表  3  吃水深度工况设置

    Table  3.   Settings of draught depth

    工况 方板尺寸 μ da/cm
    1 20 cm×20 cm 0.46 0
    2 1
    3 3
    4 5
    5 10
    6 60 cm×60 cm 1.38 0
    7 1
    8 3
    9 5
    10 10
    下载: 导出CSV

    表  4  泥沙边界工况设置

    Table  4.   Sediment boundary condition settings

    工况 方板尺寸 μ 药包距泥沙边界距离/cm λ
    1 20 cm×20 cm 0.46 10 0.46
    2 20 0.92
    3 30 1.38
    4 40 1.84
    5 60 cm×60 cm 1.38 10 0.46
    6 20 0.92
    7 30 1.38
    8 40 1.84
    下载: 导出CSV
  • [1] 秦业志, 王莹, 王志凯, 等. 小当量柱型装药水下近场爆炸固支单层方形钢板毁伤特性研究 [J]. 振动与冲击, 2021, 40(7): 29–36. DOI: 10.13465/j.cnki.jvs.2021.07.004.

    QIN Y Z, WANG Y, WANG Z K, et al. Damage characteristics of fixed single-layer square steel plate under near-field underwater explosion of small equivalent column charge [J]. Journal of Vibration and Shock, 2021, 40(7): 29–36. DOI: 10.13465/j.cnki.jvs.2021.07.004.
    [2] 张弛, 刘凯, 李海涛, 等. 水下爆炸下典型舰船结构整体损伤模式表征方法及图谱研究 [J]. 爆炸与冲击, 2022, 42(6): 065101. DOI: 10.11883/bzycj-2021-0200.

    ZHANG C, LIU K, LI H T, et al. Study on the characterization method and mode map of overall damage of typical warship structures subjected to underwater explosions [J]. Explosion and Shock Waves, 2022, 42(6): 065101. DOI: 10.11883/bzycj-2021-0200.
    [3] LIU L T, YAO X L, ZHANG A M, et al. Research on the estimate formulas for underwater explosion bubble jet parameters [J]. Ocean Engineering, 2018, 164: 563–576. DOI: 10.1016/j.oceaneng.2018.06.070.
    [4] RAMAJEYATHILAGAM K, VENDHAN C P, RAO V B. Non-linear transient dynamic response of rectangular plates under shock loading [J]. International Journal of Impact Engineering, 2000, 24(10): 999–1015. DOI: 10.1016/S0734-743X(00)00018-X.
    [5] 代利辉, 吴成, 安丰江. 水下爆炸载荷下固支方板的动态毁伤模式 [J]. 兵工学报, 2020, 41(S2): 111–119. DOI: 10.3969/j.issn.1000-1093.2020.S2.015.

    DAI L H, WU C, AN F J. Dynamic damage mode of clamped square plates subjected to underwater explosive loading [J]. Acta Armamentarii, 2020, 41(S2): 111–119. DOI: 10.3969/j.issn.1000-1093.2020.S2.015.
    [6] 汪俊, 孟利平, 伍星星, 等. 水面浮体结构底部水下爆炸射流试验研究 [J]. 船舶力学, 2022, 26(9): 13. DOI: 10.3969/j.issn.1007-7294.2022.09.014.

    WANG J, MENG L P, WU X X, et al. Experimental investigation on water-jets resulting from bubble collapse of underwater explosion under surface floating structures [J]. Journal of Ship Mechanics, 2022, 26(9): 13. DOI: 10.3969/j.issn.1007-7294.2022.09.014.
    [7] LI H T, ZHENG X Y, ZHANG C, et al. Sagging damage characteristics of hull girder with trapezoidal cross-section subjected to near-field underwater explosion [J]. Defence Technology, 2021, 21: 1–13. DOI: 10.1016/j.dt.2021.10.004.
    [8] 赖志超, 邓硕, 秦健, 等. 不同类型炸药近场水下爆炸下固支方板动态响应研究 [J/OL]. 工程力学[2023-05-04]. http://kns.cnki.net/kcms/detail/11.2595.o3.20221226.1340.003.html.

    LAI Z C, DENG S, QIN J, et al. Study on dynamic response of clamped square plates under near-field underwater explosion with different explosives [J/OL]. Engineering Mechanics[2023-05-04]. http://kns.cnki.net/kcms/detail/11.2595.o3.20221226.1340.003.html.
    [9] GAN N, LIU L T, YAO X L, et al. Experimental and numerical investigation on the dynamic response of a simplified open floating slender structure subjected to underwater explosion bubble [J]. Ocean Engineering, 2021, 219: 108308. DOI: 10.1016/j.oceaneng.2020.108308.
    [10] ZHANG A M, YAO X L, LI J. The interaction of an underwater explosion bubble and an elastic-plastic structure [J]. Applied Ocean Research, 2008, 30(3): 159–171. DOI: 10.1016/j.apor.2008.11.003.
    [11] 王诗平, 孙士丽, 张阿漫, 等. 冲击波和气泡作用下舰船结构动态响应的数值模拟 [J]. 爆炸与冲击, 2011, 31(4): 367–372. DOI: 10.11883/1001-1455(2011)04-0367-06.

    WANG S P, SUN S L, ZHANG A M, et al. Numerical simulation of dynamic response of warship structures subjected to underwater explosion shock waves and bubbles [J]. Explosion and Shock Waves, 2011, 31(4): 367–372. DOI: 10.11883/1001-1455(2011)04-0367-06.
    [12] 文彦博, 胡亮亮, 秦健, 等. 近场水下爆炸气泡脉动及水射流的实验与数值模拟研究 [J]. 爆炸与冲击, 2022, 42(5): 053203. DOI: 10.11883/bzycj-2021-0206.

    WEN Y B, HU L L, Q J, et al. Experimental study and numerical simulation on bubble pulsation and water jet in near-field underwater explosion [J]. Explosion and Shock Waves, 2022, 42(5): 053203. DOI: 10.11883/bzycj-2021-0206.
    [13] 王树山, 李梅, 马峰. 爆炸气泡与自由水面相互作用动力学研究 [J]. 物理学报, 2014, 63(19): 194703. DOI: 10.7498/aps.63.194703.

    WANG S S, LI M, MA F. Dynamics of the interaction between explosion bubble and free surface [J]. Acta Physica Sinica, 2014, 63(19): 194703. DOI: 10.7498/aps.63.194703.
    [14] LIU N N, CUI P, REN S F, et al. Study on the interactions between two identical oscillation bubbles and a free surface in a tank [J]. Physics of Fluids, 2017, 29(5): 052104. DOI: 10.1063/1.4984080.
    [15] HUNG C F, HWANGFU J J. Experimental study of the behaviour of mini-charge underwater explosion bubbles near different boundaries [J]. Journal of Fluid Mechanics, 2010, 651: 55–80. DOI: 10.1017/S0022112009993776.
    [16] ZHANG A M, YAO X L, FENG L H. The dynamic behavior of a gas bubble near a wall [J]. Ocean Engineering, 2009, 36(3/4): 295–305. DOI: 10.1016/j.oceaneng.2008.12.006.
    [17] 张之凡, 谢宇杰, 王成, 等. 近自由面水下爆炸气泡与破损结构耦合作用机理研究 [J]. 北京理工大学学报, 2022, 42(9): 909–917. DOI: 10.15918/j.tbit1001-0645.2022.103.

    ZHANG Z F, XIE Y J, WANG C, et al. Coupling mechanism between damaged structure and underwater explosion bubble near free surface [J]. Transactions of Beijing institute of Technology, 2022, 42(9): 909–917. DOI: 10.15918/j.tbit1001-0645.2022.103.
    [18] 贺铭, 张阿漫, 刘云龙. 近场水下爆炸气泡与双层破口结构的相互作用 [J]. 爆炸与冲击, 2020, 40(11): 111402. DOI: 10.11883/bzycj-2020-0110.

    HE M, ZHANG A M, LIU Y L. Interaction of the underwater explosion bubbles and nearby double-layer structures with circular holes [J]. Explosion and Shock Waves, 2020, 40(11): 111402. DOI: 10.11883/bzycj-2020-0110.
    [19] 金辉, 张庆明, 高春生, 等. 不同边界条件水下爆炸气泡脉动对比的试验研究 [J]. 兵工学报, 2009, 30(S2): 213–217. DOI: CNKI:SUN:BIGO.0.2009-S2-045.

    JIN H, ZHANG Q M, GAO C S, et al. Comparison experimental study of underwater explosion bubble pulse among the different boundaries [J]. Acta Armamentarii, 2009, 30(S2): 213–217. DOI: CNKI:SUN:BIGO.0.2009-S2-045.
    [20] 金辉, 李兵, 权琳, 等. 不同边界条件下炸药水中爆炸的能量输出结构 [J]. 爆炸与冲击, 2013, 33(3): 325–330. DOI: 10.11883/1001-1455(2013)03-0325-05.

    JIN H, LI B, QUAN L, et al. Configuration of explosive energy output in different underwater boundary conditions [J]. Explosion and Shock Waves, 2013, 33(3): 325–330. DOI: 10.11883/1001-1455(2013)03-0325-05.
    [21] LINDAU O, LAUTERBORN W. Cinematographic observation of the collapse and rebound of a laser-produced cavitation bubble near a wall [J]. Journal of Fluid Mechanics, 2003, 479: 327–348. DOI: 10.1017/S0022112002003695.
    [22] JAYAPRAKASH A, HSIAO C, CHAHINE G. Numerical and experimental study of the interaction of a spark-generated bubble and a vertical wall [J]. Massachusetts Institute of Technology, 2012, 134(3): 381–382. DOI: 10.1115/1.4005688.
    [23] MA X, HUANG B, ZHAO X, ET AL. Comparisons of spark-charge bubble dynamics near the elastic and rigid boundaries [J]. Ultrasonics Sonochemistry, 2018, 43: 80–90. DOI: 10.1016/j.ultsonch.2018.01.005.
    [24] ZHANG A M, CUI P, WANG Y. Experiments on bubble dynamics between a free surface and a rigid wall [J]. Experiments in Fluids, 2013, 54: 1602. DOI: 10.1007/s00348-013-1602-7.
    [25] HUANG G H, ZHANG M D, MA X J, et al. Dynamic behavior of a single bubble between the free surface and rigid wall [J]. Ultrasonics Sonochemistry, 2020, 67: 105147. DOI: 10.1016/j.ultsonch.2020.105147.
    [26] 陈志鹏. 气泡与复杂边界耦合机理研究[D]. 江苏镇江: 江苏科技大学, 2019.

    CHEN Z P. A study on the coupling effect of interaction between bubbles and complex boundaries [D]. Zhenjiang, Jiangsu, China: Jiangsu University of Science and Technology, 2019.
    [27] TIAN Z L, LIU Y L, ZHANG A M, et al. Jet development and impact load of underwater explosion bubble on solid wall [J]. Applied Ocean Research, 2020, 95: 102013. DOI: 10.1016/j.apor.2019.102013.
    [28] 张桂夫, 朱雨建, 杨基明. 水下爆炸冲击凹陷液面诱导射流研究 [J]. 爆炸与冲击, 2018, 38(2): 241–249. DOI: 10.11883/bzycj-2016-0238.

    ZHANG G F, ZHU Y J, YANG J M. A study on jet flow induced by underwater explosion at a pit-interface [J]. Explosion and Shock Waves, 2018, 38(2): 241–249. DOI: 10.11883/bzycj-2016-0238.
    [29] XU L Y, WANG S P, LIU Y L, et al. Numerical simulation on the whole process of an underwater explosion between a deformable seabed and a free surface [J]. Ocean Engineering, 2020, 219: 108311. DOI: 10.1016/j.oceaneng.2020.108311.
    [30] LEE E, FINGER M, COLLINS W. JWL equation of state coefficients for high explosives [R]//Office of Scientific and Technical Information Technical Reports, 1973. DOI: 10.2172/4479737.
    [31] 李晓杰, 张程娇, 王小红, 等. 水的状态方程对水下爆炸影响的研究 [J]. 工程力学, 2014, 31(8): 46–52. DOI: 10.6052/j.issn.1000-4750.2013.03.0180.

    LI X J, ZHANG C J, WANG X H, et al. Numerical study on the effect of equations of state of water on underwater explosions [J]. Engineering Mechanics, 2014, 31(8): 46–52. DOI: 10.6052/j.issn.1000-4750.2013.03.0180.
    [32] YANG X Q, YANG H, GARDNER L, et al. A continuous dynamic constitutive model for normal-and high-strength structural steels [J]. Journal of Constructional Steel Research, 2022, 192: 107254. DOI: 10.1016/j.jcsr.2022.107254.
    [33] 孙远翔, 田俊宏, 张之凡, 等. 含铝炸药近场水下爆炸冲击波的实验及数值模拟 [J]. 振动与冲击, 2020, 39(14): 171–178, 193. DOI: 10.13465/j.cnki.jvs.2020.14.025.

    SUN Y X, TIAN J H, ZHANG Z F, et al. Experiment and numerical simulation study on the near-field underwater explosion of aluminized explosive [J]. Journal of Vibration and Shock, 2020, 39(14): 171–178, 193. DOI: 10.13465/j.cnki.jvs.2020.14.025.
    [34] COLE R H. Underwater explosion [M]. New Jersey: Princeton University Press, 1948.
    [35] HU J, CHEN Z Y, ZHANG X D, et al. Underwater explosion in centrifuge. part I: validation and calibration of scaling laws [J]. Science China Technological Sciences, 2017, 60(11): 1638–1657. DOI: 10.1007/s11431-017-9083-0.
  • 期刊类型引用(1)

    1. 陈钰帆,黄锡堤,赖志超,秦健,孟祥尧,文彦博,黄瑞源. 板结构水下爆炸下的空化特性. 火炸药学报. 2024(08): 714-729 . 百度学术

    其他类型引用(1)

  • 加载中
图(23) / 表(4)
计量
  • 文章访问数:  261
  • HTML全文浏览量:  65
  • PDF下载量:  194
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-05-05
  • 修回日期:  2023-09-12
  • 网络出版日期:  2023-09-13
  • 刊出日期:  2023-11-17

目录

/

返回文章
返回