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

铆接油箱的水锤毁伤效应

张景飞 贾豪博 任柯融 卿华 郭攀 杜晓伟 陈荣 卢芳云

张景飞, 贾豪博, 任柯融, 卿华, 郭攀, 杜晓伟, 陈荣, 卢芳云. 铆接油箱的水锤毁伤效应[J]. 爆炸与冲击, 2023, 43(7): 073301. doi: 10.11883/bzycj-2022-0275
引用本文: 张景飞, 贾豪博, 任柯融, 卿华, 郭攀, 杜晓伟, 陈荣, 卢芳云. 铆接油箱的水锤毁伤效应[J]. 爆炸与冲击, 2023, 43(7): 073301. doi: 10.11883/bzycj-2022-0275
ZHANG Jingfei, JIA Haobo, REN Kerong, QING Hua, GUO Pan, DU Xiaowei, CHEN Rong, LU Fangyun. Damage of hydrodynamic ram effect to riveted fuel tanks[J]. Explosion And Shock Waves, 2023, 43(7): 073301. doi: 10.11883/bzycj-2022-0275
Citation: ZHANG Jingfei, JIA Haobo, REN Kerong, QING Hua, GUO Pan, DU Xiaowei, CHEN Rong, LU Fangyun. Damage of hydrodynamic ram effect to riveted fuel tanks[J]. Explosion And Shock Waves, 2023, 43(7): 073301. doi: 10.11883/bzycj-2022-0275

铆接油箱的水锤毁伤效应

doi: 10.11883/bzycj-2022-0275
基金项目: 湖南省自然科学基金(2022JJ10058)
详细信息
    作者简介:

    张景飞(1974- ) 男,硕士生导师,副教授,zjf723@163.com

    通讯作者:

    任柯融(1993- ),男,博士研究生,讲师,renkerong@nudt.edu.cn

  • 中图分类号: O385;V221.91

Damage of hydrodynamic ram effect to riveted fuel tanks

  • 摘要: 为了研究高速侵彻体撞击飞机油箱等充液容器所产生的水锤效应对容器结构产生的灾难性破坏,以铆接油箱为对象,通过开展弹道射击实验,结合数字图像相关测试技术,获取了铆接油箱在射弹冲击作用下的箱体变形、破孔直径等数据;并建立流-固耦合有限元模型,分析射弹入射速度与射弹动能损失、箱体变形、流体动压、铆钉失效之间的关系。结果表明:有限元模拟结果与实验结果基本吻合,数值模型可用于描述油箱在水锤作用下的动力学行为;射弹动能损失、箱体变形量、液体压力峰值与射弹入射速度呈正比例关系;当射弹入射速度达到1400 m/s之后,油箱后壁板开始出现裂纹,并呈花瓣式破孔损伤;当射弹入射速度达到1600 m/s时,铆钉开始发生断裂。
  • 燃油系统是飞机易损部件中暴露面积最大的系统,其安全性对于飞机作战生存力至关重要。当射弹以极高的速度冲击充满液体的容器时,射弹携带的动能将转化为液体的冲击波能量,并最终以液体压力波的形式传递到容器壁面,从而对容器结构造成巨大破坏,这种流-固耦合现象称为水锤效应[1]。据统计,在美军开展的“沙漠风暴”行动中,有75%的飞机战斗损伤事件都与油箱受损有关,且在这些飞机战斗受损事件中水锤效应导致的战损比例最大[2]

    为了研究水锤效应对油箱的毁伤机制,自20世纪70年代开始就有学者开展了大量的实验。Ball[3] 首先使用不同能量等级的射弹对充液箱体进行了水锤效应的实验研究,测量了射弹侵彻油箱过程中的液体压力和入射板应变,发现影响射弹能量衰减的主要因素是射弹的质量和入射角度。Disimile等[4]使用不同材料的射弹冲击模拟油箱,记录了油箱中不同位置的压力数据,分析了射弹冲击油箱时流体压力的空间分布规律。Varas等[5]使用质量为8 g的钢球以不同的速度、不同的充液比进行了射弹撞击充水铝管的实验,研究了射弹速度和油箱充液比对水锤效应的影响。Ren等[6]采用不同形状的射弹以不同的速度撞击油箱,建立了射弹动能、形状与容器壁板变形的关系。李营等[7]设计了水间隔靶板在球形弹体撞击下动态响应的实验装置,对比了不同靶板在弹体侵彻作用下的变形特点。陈安然等[8]等通过开展高速破片侵彻充液容器的实验,得到了水锤效应引起的液体喷溅特性及其影响因素。近些年,随着计算机技术的快速发展,有限元方法成为了一种实用高效的数值分析方法,并逐渐应用于各个科学领域。李亚智等[9]通过MSC. DYTRAN软件对空箱和充液油箱在高速弹丸冲击下的动态响应进行了数值模拟,得到了2种箱体壁板的应力、应变时间历程。Varas等[10]使用了ALE算法模拟了钢球弹丸冲击充水铝管所产生的水锤效应,并与实验结果进行对比,验证了ALE算法可以较为真实地模拟水锤效应的各个阶段。Mansoori等[11]模拟了射弹冲击充水容器和不充水容器的整个过程,分析了容器壁面的失效形式。蓝肖颖等[12]建立了双射弹冲击充水油箱的有限元模型,研究了液压水锤的叠加效应。韩璐等[13-14]对飞机油箱水锤效应的影响因素和影响程度进行了综合分析。

    在上述研究中,所关注的研究对象主要集中于通过螺栓固定的模拟等效油箱,而对具有飞机结构特征油箱关注较少。真实的飞机油箱结构大多是通过铆接进行连接固定的。陈钢[15]研究发现,在常规金属结构中,构件之间的连接接头是最容易被破坏的部位。飞机结构中的铆钉本来设计为以剪切力的形式进行载荷传递;然而当机翼油箱壁板(或机翼蒙皮)被液压水锤垂直作用时,会发生弯曲变形,从而导致连接这些壁板的铆钉将处于纵向拉应力和横向剪切应力同时作用的复杂应力状态下。如果铆钉无法承受这种复杂应力加载,将发生断裂失效,最终导致油箱壁板(或机翼蒙皮)脱落[16]。为了对液压水锤载荷下的油箱接头强度进行定量评估,莱特-帕特森空军基地开发出了水锤模拟器[17],该设备可以对飞机结构接头进行可控的液压水锤加载,从而获取了接头的动态破坏极限,验证了以承受剪力为主的铆接接头在水锤载荷的正向加载下并不具备理想的动态强度,但是他们并没有考虑水锤对铆接油箱的整体毁伤效应。

    为了探究铆接油箱在水锤作用下的动态响应机制,本文开展对铆接油箱的弹道射击实验,结合数字图像相关测试技术,测量铆接油箱在射弹冲击作用下的箱体变形、破孔直径等数据;并进行流固耦合的有限元数值模拟,分析射弹入射速度与射弹动能损失、箱体变形、流体动压、铆钉失效之间的规律;结合实验和模拟结果,探索铆接油箱在速度为780~1600 m/s的钨合金射弹冲击下的动态力学响应行为及水锤对铆接油箱的毁伤机理,建立射弹入射参量-水锤参量-铆接油箱毁伤参量的映射关系。

    采用弹道射击实验对铆接油箱进行水锤加载,实验现场整体布置如图1所示。实验装置主要包括弹道枪、测速靶、铆接油箱以及高速摄像机。

    图  1  实验场景布置
    Figure  1.  Experimental layout

    实验中使用的射弹质量为69.7 g,长度为51.6 mm,直径为10.8 mm,形状为尖头形,材质为钨合金。射弹由一门12.7 mm口径的弹道枪发射,设计着靶速度为780 m/s,射击瞄准点为油箱前壁板的几何中心位置。采用断靶法测试射弹的入射速度,将2张测速靶布置于油箱正前方约10 m处。

    铆接油箱为长方体结构,尺寸为440 mm×190 mm×416 mm,由6块2024铝合金板铆接而成,内部填充95%容量(约23.4 L)的航空煤油。前后壁板厚度为2 mm,油箱其余壁板厚度为3 mm。油箱上的铆钉分为2种,一种是用来连接油箱上下壁板和侧壁板的304L不锈钢铆钉,另一种是连接前后壁板和其余壁板的5052铝合金铆钉,其直径均为4 mm。

    实验中,使用1台高速摄像机记录整个撞击过程,拍摄频率为300 000 Hz,位于油箱前方45°方向10 m处。同时,为了记录油箱后壁板的变形情况,将2台高速摄像机设置在油箱侧后方3 m处,拍摄频率为50 000 Hz,对拍摄的结果进行数字图像相关(digital image correlation, DIC)技术处理[18]以获取后壁板的应变场历史。

    为了更好地模拟射弹高速冲击下油箱箱体及航空煤油的动力学响应过程,采用ALE流-固耦合算法进行仿真计算[5]图2为所建立的与实验相同尺寸的有限元模型,模型包括油箱箱体、射弹、煤油、空气和铆钉5部分,其中流体域包括煤油和空气,使用多物质Euler单元描述;固体域包括油箱箱体、射弹和铆钉,使用Lagrange单元描述,模型中的所有单元均采用六面体单元划分。为了提高计算精度,并控制计算成本,本文采用1/2平面对称模型进行计算,并将射弹运动路径周围的单元进行加密划分,未加密区域单元尺寸约为4 mm,加密区域单元尺寸约为1 mm,整个模型共包含1 602 566个单元。空气域四周设为无反射边界条件,油箱箱体之间以及油箱箱体和铆钉之间使用关键字*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE定义为自动面-面接触,射弹和油箱前后壁板之间使用关键字*CONTACT_ERODING_SURFACE_TO_SURFACE定义为侵蚀面-面接触。在建立有限元模型时,将流体域的空气和煤油单元之间进行共节点处理,利用关键字*CONSTRAINED_LAGRANGE_IN_SOLID实现Lagrange单元和Euler单元之间的相互作用,进而实现流-固耦合计算。

    图  2  有限元模型
    Figure  2.  Finite element model

    在材料模型方面,射弹和箱体分别为钨合金和2024铝合金,使用 Johnson-Cook本构模型、Johnson-Cook失效模型和Grüneisen状态方程描述,具体取值详见文献[19-20]。铆钉材料为304L不锈钢和5052铝合金,使用塑性随动(plastic kinematic, PK)本构模型和等效塑形应变(equivalent plastic strain, EQPS)失效准则描述,具体取值详见文献[21-22]。在EQPS失效准则中,当单元的EQPS值达到临界值时,就会被删除,不参与下一步的计算。煤油和空气分别采用Grüneisen状态方程和先行多项式状态方程描述,具体取值详见文献[23]。

    Johnson-Cook本构模型[20]广泛的应用于描述材料的大变形行为,其屈服应力表达式为:

    σ=(ABˉεnp)(1+Cjln˙ε)(1Tmh)
    (1)

    式中:A为材料初始屈服应力;B为材料应变硬化模量;n为硬化指数;Cj为材料应变率参数;m为材料热软化参数;ˉεp为等效塑形应变,˙ε=˙ˉεp/˙ε0为应变率˙ε0=1 s−1的塑性比,˙ˉεp为等效塑性应变率;相对温度Th=(TTr)/(TmTr)Tr为室温,Tm为熔点。

    在Johnson-Cook失效模型中,材料的失效应变εf表达式[20]为:

    εf=[D1+D2exp(D3σhσe)](1+D4ln˙ε)[1+D5θθ0θmθ0]
    (2)

    式中:σeσh为等效应力和静水压力;当损伤参数Σ(Δˉεp/εf)=1时材料发生断裂,Σ代表对断裂前的变形过程求和;Δˉεp为等效塑性应变增量;D1D5为常数。

    Grüneisen 状态方程[24]用于定义物质在压缩和膨胀状态下的压力变化,其表达式为:

    p1=ρ0c02u[1 + (1γ0/2)uau2/2][1(S11)uS2u2/(u+1)S3u2/(u+1)2]2+(γ0+au)e
    (3)
    p2=ρ0c20u+(γ0+au)e
    (4)

    式中:p1p2分别为压缩材料和膨胀材料的压力;c0为声速体积;S1、S2、S3为材料参数;γ0为 Grüneisen系数,aγ0的一阶体积修正系数;u=ρ/ρ0−1,ρ为某一静水压力下材料的瞬时密度,ρ0为材料的初始密度,e为单位体积的初始内能。

    在PK本构模型[21]中,塑性流动应力与塑形应变、应变率具有如下关系:

    σ=(σ0+βEpεeffp)[1+(˙ε/Cpk)1/P]
    (5)

    式中:σ0为初始屈服应力,˙ε为应变率,εeffp为有效塑形应变,β为硬化参数,Ep为塑性硬化模量,CpkP为应变率参数。

    图3为实验测试得到的不同时刻油箱后壁板z轴方向的位移云图,同时给出了有限元计算结果进行对比。起始时刻为射弹击中油箱前的时刻。从图中可以看出随着时间的推移,油箱后壁板的变形将逐渐增大。

    图  3  不同时刻实验与模拟中的油箱后壁板变形云图对比
    Figure  3.  Comparison of deformation contours of the back wall of the fuel tank in experiment and simulation at different time

    图4给出了实验和模拟得到的后壁板变形结果对比。其中,实验和模拟结果的相关系数使用Pearson公式[25]计算,表达式如下:

    图  4  不同时刻油箱后壁板z轴变形
    Figure  4.  z-axis deformation of the back wall of the fuel tank at different times
    R=cov(X,Y)σXσY=E(XY)E(X)E(Y)σXσY
    (6)

    式中:cov(X,Y)为两组数据XY的协方差; σXσY分别为XY的标准差;E(X)、E(X)、E(XY)分别为X、Y、XY的期望值。|R|越接近1,表示两组数据的相关性越强。

    图5所示,测量实验和仿真结果中后壁板的破孔尺寸分别为15和14 mm,相对误差为6.67%。图6给出油箱后壁板弹孔上方20 mm处壁板单元的速度时程曲线,从图可以看出,该点在射弹入射油箱后0.28 ms左右达到了速度峰值,其中实验和仿真所得速度峰值分别为51.81和54.84 m/s,误差为5.81%。因此,本文所建的有限元模型计算精度较好,能够较好地描述实验结果。

    图  5  后壁板破孔实验结果和模拟结果的对比
    Figure  5.  Comparison of the experiment result and the simulation result of the broken hole in the back wall
    图  6  后壁板弹孔上方20 mm的速度
    Figure  6.  Velocity of units at 20 mm above the projectile hole in the back wall
    3.2.1   射弹运动过程及动能衰减分析

    为了进一步探究铆接油箱在具有不同初始动能射弹作用下的毁伤规律,开展了射弹以780、1000、1200、1400、1600 m/s等5个速度工况垂直入射油箱的有限元计算。当射弹击中油箱时会在油箱箱体和煤油的作用下做减速运动,并将一部分动能转化为自身的变形能、煤油的动能和内能以及油箱箱体的动能和变形能。仿真计算所得的射弹加速度随时间变化曲线如图7所示,从图中可以看出射弹在侵彻油箱的过程中主要经历3个阶段:

    图  7  射弹加速度时程曲线
    Figure  7.  Projectile acceleration history curves

    (1) 射弹击中油箱前壁板时加速度急剧增大,速度快速下降;根据动量方程可知射弹撞击油箱壁时所产生的冲击力与射弹的速度成正比,因此当射弹入射速度为1600 m/s时,加速度最大,达到了106 m/s2

    (2) 射弹在煤油中运动阶段,加速度保持稳定,速度稳步下降;射弹在液体中的运动阻力表达式[23]如下

    f=mdvdt=12cdρlv2S
    (7)

    式中:f为射弹在液体中运动所受阻力;cd为阻力系数,与液体雷诺数有关,在数值计算中通常假设液体不可压缩,因此cd为常数,取值0.275[14]ρl为液体密度;v为射弹速度;S为射弹横截面积。由式(7)可知,同一种形状和材质的射弹在煤油中运动时所受阻力与射弹速度的平方成正比,因此射弹入射速度为1600 m/s时在煤油中运动的加速度相较1400 m/s入射时提高了33.8%;

    (3) 在射弹穿出油箱后壁板前,流体内的冲击波提前到达后壁板,并在后壁板处产生预应力。

    图8给出了油箱前、后壁板对射弹施加的侵彻阻力以及射弹到达后壁板前冲击波对后壁板中心与弹道线交汇处单元施加的预应力。从图中可以看出油箱前壁板对射弹施加的侵彻阻力明显大于后壁板,当射弹入射速度从780 m/s提升到1600 m/s,油箱前壁板对射弹施加的侵彻阻力从16.75 kN增加到56.8 kN, 而油箱后壁板对射弹造成的侵彻阻力则维持在20 kN左右,同时煤油内的冲击波对油箱后壁板造成的预应力也从441 MPa提高到537 MPa。可见,射弹的速度越大,在煤油中内形成的冲击波越强,对油箱后壁板产生的预拉应力也就越大;而预拉应力越大使得后壁板更容易被侵彻[26],因此导致后壁板的对不同速度射弹的侵彻阻力都恒为一个定值。

    图  8  油箱前后壁板对射弹的阻力及冲击波对后壁板单元的预应力
    Figure  8.  Projectile resistance from the front and back walls of the tank and prestress of the back wall elements by the shock wave

    图9给出射弹在不同入射速度下的动能损失变化关系。从图7图9中可以看出,当射弹入射速度为1600、1400、1200、1000 m/s时,在煤油内运动的加速度相较与780 m/s分别提高了405.9%、259.2%、160.9%、66.3%;出射时间分别提前了0.125、0.109、0.084、0.049 ms;动能损失增加了1341、901、504、210 J。因此,射弹入射速度越高,在煤油中运动所受的阻力越大,贯穿油箱所用的时间越短,动能损失越多。

    图  9  射弹入射速度对动能损失的影响
    Figure  9.  Influence of projectile incident velocity on kinetic energy loss
    3.2.2   油箱液体压力分析

    图10是在射弹在780~1600 m/s的速度下煤油内的压力云图,可以看出当射弹撞击煤油时,会在煤油内产生初始冲击波。其中,射弹头部周围形成的冲击波强度最大,并且随着冲击波传递距离的增加,强度不断衰减。

    图  10  不同射弹入射速度下煤油的压力分布演化
    Figure  10.  Pressure contours of kerosene at different projectile impact velocities

    图11给出了射弹入射速度、初始冲击波强度与持续时间的变化关系,可以看出射弹入射速度越大,在煤油中形成的初始冲击波强度也越大,同时随着射弹入射速度的增加,穿透油箱所用的时间更短,因此初始冲击波持续的时间也越短。

    图  11  射弹入射速度对初始冲击波压力和持续时间的影响
    Figure  11.  Infuence of impact velocity of projectile on pressure of initial shock wave and its duration

    为对比不同射弹入射速度下油箱侧壁面所受载荷大小,选取油箱右壁板中心处液体单元观察其压力变化。图12为射弹以不同速度撞击油箱时油箱右壁板中心处煤油压力变化曲线。由于射弹撞击油箱后产生的冲击波波速与射弹穿透入射壁板所携带的能量呈正比例关系[27],因此射弹入射速度越快,煤油压力峰值出现的时间越早,且峰值压力越大。如图13所示,射弹入射速度为1600 m/s时,峰值压力出现的时间相较780 m/s时提前了0.0075 ms,而峰值压力大小提高了392.29%。同时可以发现侧向冲击波由于在传播过程中出现了衰减,导致其压力值迅速下降。例如,射弹速度为1600 m/s时,煤油中形成的初始冲击波强度高达550 MPa,但是当冲击波传递到油箱右壁板时峰值压力衰减至仅7 MPa,因此油箱侧壁板的变形较小。

    图  12  壁面测点处的煤油压力时程曲线
    Figure  12.  Time-history curves of kerosene pressure at measuring point on the side wall
    图  13  油箱前后壁板挠度变化
    Figure  13.  Deflection evolutions of the front and back walls
    3.2.3   油箱后壁板变形分析

    当射弹高速侵彻油箱时,煤油吸收射弹动能后会转化为流体动压载荷,进而对油箱造成严重的变形与破坏。图13给出了不同射弹着速下油箱前后壁板挠度变化的时程曲线图。对比图13(a)和图13(b)可以看出:油箱后壁板的变形程度明显大于油箱前壁板,如当射弹入射速度为1600 m/s时,前壁板最大挠度(t=0.5 ms)为13.44 mm,后壁板最大挠度(t=0.5 ms)为27.64 mm,因此后壁板的毁伤程度甚于前壁板。

    图14图15分别为射弹在不同入射速度下油箱后壁板同一时刻的等效应力云图和破孔尺寸、裂纹尺寸、翘曲尺寸变化情况。从图中可以看出,后壁板的损伤主要有以下特点:(1) 当射弹入射速度小于1400 m/s时,后壁板的损伤模式为冲塞失效、翘曲变形;当射弹入射速度达到1400 m/s时,后壁板开始出现裂纹;并且随着射弹入射速度的增大,油箱后壁板裂纹逐渐向四周扩展;当射弹入射速度达到1600 m/s时,裂纹半径扩展到14.7 mm;(2) 后壁板的应力分布呈环状间隔分布,在破孔和裂纹处出现应力集中,并且射弹入射速度越大,后壁板的应力就越大;(3) 随着射弹入射速度的增大,后壁板的变形程度有着明显的增加,当射弹入射速度为1600、1400、1200、1000 m/s时,后壁板最大翘曲尺寸是射弹入射速度为780 m/s的1.13、1.45、1.81、2.3倍。可见射弹入射速度的增加增强了水锤效应对油箱的破坏作用。

    图  14  后壁板等效应力云图
    Figure  14.  Equivalent stress contours of the back wall
    图  15  不同射弹入射速度下后壁板的损伤模式演化
    Figure  15.  Evolution of the damage patterns of the back wall at different projectile impact velocities
    3.2.4   铆钉变形分析

    对于铆接油箱来说,射弹在侵彻油箱的过程中,油箱箱体在射弹和液压水锤的作用下发生变形,同时还将带动着铆钉发生变形。因此射弹除了将自身的动能传递给煤油和油箱箱体外,还会将一部分能量转化为铆钉的变形能。当这一部分能量超过铆钉材料所能承受的极限时,铆钉就会发生断裂,从而造成油箱壁板脱离。本次实验及数值模拟中射弹的入射位置为油箱正中心,当射弹入射速度为1600 m/s时,油箱后壁板弹孔下方的铆钉最先发生断裂,之后相邻两侧的铆钉相继发生断裂 ,其余的铆钉发生了变形。图16图17分别给出油箱后壁板下侧中间位置铆钉发生最大变形时的等效应力云图和不同方向变形情况。可以看出:随着射弹入射速度的增大,铆钉的塑形变形也在增大,且铆钉的变形程度与射弹入射速度呈正比例关系;当射弹入射速度达到1600 m/s时,铆钉中间部位单元达到了等效塑性应变的极限,因此发生断裂,此时单元的轴向拉伸应变为0.275,横向剪切应变为0.368。如图18所示,当射弹速度到达1600 m/s时,铆钉发生断裂失效,油箱后壁板有脱落风险。

    图  16  铆钉发生最大变形时的等效应力云图
    Figure  16.  Equivalent stress of the rivets contours at maximum deformation
    图  17  铆钉不同方向变形
    Figure  17.  Rivet deflection at different directions
    图  18  不同射弹入射速度下铆钉变形对比
    Figure  18.  Comparison of rivet deformation under different projectile impact velocities

    图19给出了射弹入射速度为1600 m/s时,煤油、油箱前壁板、油箱后壁板以及铝合金铆钉的总能量(动能+变形能)时程曲线,图20给出了射弹入射速度为1600 m/s时,油箱后壁板和铝合金铆钉的动能、变形能(由于铆钉的动能太小,忽略不计)时程曲线。图中射弹入射时刻为0时刻。可以发现煤油、油箱以及铆钉的能量变化存在以下特征:(1) 在射弹入射后至0.0995 ms时,煤油的总能量急剧增加,同时前壁板能量也不断增加,说明在射弹入射初期,射弹的动能主要是转换为煤油的总能量以及前壁板的总能量;(2) 从0.0995至0.1080 ms时,煤油总能量陡峭下跌,同时后壁板总能量(主要是动能)迅速增加,此时水锤已经运动抵达后壁板,引起了后壁板的运动与变形;(3) 从0.1080至0.6310 ms时刻,射弹从后板穿出,由于失去了来自于射弹的能量输入,后壁板的总能量增长率明显放缓,并且在0.6310 ms时刻后,后壁板的总能量开始减少,同时可以发现后壁板总能量明显高于前壁板总能量,即射弹动能绝大多数是被后壁板所吸收;(4) 在射弹于0.1080 ms时从后壁板射出之后,后板的动能经过一小段时间的上升后就开始下降,变形能开始上升,同时铆钉的变形能也开始有了明显的上升,此时油箱后壁板在受到液压水锤的作用下发生塑形变形,并将自身的动能转移为变形能和铆钉的变形能;(5) 在0.715 ms时,铆钉的变形能达到了34.5 J并且开始小幅度的下降,此时铆钉发生了断裂。

    图  19  油箱各构件总能量时程曲线
    Figure  19.  Total energy curves for various parts of the tank
    图  20  壁板和铆钉能量变化时程曲线
    Figure  20.  Energy curves of the back wall and the rivet

    本文对铆接油箱进行了弹道射击实验与有限元数值仿真,并通过DIC测试结果对比验证了数值模型的有效性。更进一步,通过有限元模拟研究了射弹在780~1600 m/s速度冲击下铆接油箱的毁伤规律。主要结论如下:

    (1) 射弹在侵彻油箱前壁板时速度急剧下降,之后在煤油内运动时的加速度趋于平缓;煤油内的冲击波相较于射弹将会提前到达油箱后壁板对其产生一个预应力,当射弹入射速度从780 m/s提升到1600 m/s时,冲击波对油箱后壁板造成的预拉应力也从441 MPa提高到537 MPa,导致虽然射弹入射速度逐渐提升,但是油箱后壁板对射弹的侵彻阻力一直维持在20 kN左右;

    (2) 射弹高速侵彻油箱会在煤油中产生初始冲击波,射弹入射速度从780 m/s提升至1600 m/s,煤油内的初始冲击波强度从150 MPa提升为550 MPa;初始冲击波以球面波的形式以入射点为中心向外扩展侧向冲击波衰减迅速,在射弹入射速度为1600 m/s的情况下,煤油中形成的初始冲击波幅值高达550 MPa,到达右壁板时冲击波幅值仅为7 MPa,因此油箱侧壁板的变形较小;

    (4) 在射弹侵彻与液压水锤的共同作用下,油箱后壁板会产生严重大变形,当射弹速度低于1400 m/s时,后壁板损伤模式以冲塞失效、翘曲变形为主;当射弹速度高于1400 m/s时,后壁板损伤模式以花瓣型开裂、翘曲大变形为主;对于铆接油箱而言,由于油箱后壁板并不是固支边界,其变形将会带动边界变形,即铆钉变形;当射弹入射速度达到1600 m/s,铆钉发生断裂,油箱后壁板将出现脱落,有煤油大量泄漏的风险。

  • 图  1  实验场景布置

    Figure  1.  Experimental layout

    图  2  有限元模型

    Figure  2.  Finite element model

    图  3  不同时刻实验与模拟中的油箱后壁板变形云图对比

    Figure  3.  Comparison of deformation contours of the back wall of the fuel tank in experiment and simulation at different time

    图  4  不同时刻油箱后壁板z轴变形

    Figure  4.  z-axis deformation of the back wall of the fuel tank at different times

    图  5  后壁板破孔实验结果和模拟结果的对比

    Figure  5.  Comparison of the experiment result and the simulation result of the broken hole in the back wall

    图  6  后壁板弹孔上方20 mm的速度

    Figure  6.  Velocity of units at 20 mm above the projectile hole in the back wall

    图  7  射弹加速度时程曲线

    Figure  7.  Projectile acceleration history curves

    图  8  油箱前后壁板对射弹的阻力及冲击波对后壁板单元的预应力

    Figure  8.  Projectile resistance from the front and back walls of the tank and prestress of the back wall elements by the shock wave

    图  9  射弹入射速度对动能损失的影响

    Figure  9.  Influence of projectile incident velocity on kinetic energy loss

    图  10  不同射弹入射速度下煤油的压力分布演化

    Figure  10.  Pressure contours of kerosene at different projectile impact velocities

    图  11  射弹入射速度对初始冲击波压力和持续时间的影响

    Figure  11.  Infuence of impact velocity of projectile on pressure of initial shock wave and its duration

    图  12  壁面测点处的煤油压力时程曲线

    Figure  12.  Time-history curves of kerosene pressure at measuring point on the side wall

    图  13  油箱前后壁板挠度变化

    Figure  13.  Deflection evolutions of the front and back walls

    图  14  后壁板等效应力云图

    Figure  14.  Equivalent stress contours of the back wall

    图  15  不同射弹入射速度下后壁板的损伤模式演化

    Figure  15.  Evolution of the damage patterns of the back wall at different projectile impact velocities

    图  16  铆钉发生最大变形时的等效应力云图

    Figure  16.  Equivalent stress of the rivets contours at maximum deformation

    图  17  铆钉不同方向变形

    Figure  17.  Rivet deflection at different directions

    图  18  不同射弹入射速度下铆钉变形对比

    Figure  18.  Comparison of rivet deformation under different projectile impact velocities

    图  19  油箱各构件总能量时程曲线

    Figure  19.  Total energy curves for various parts of the tank

    图  20  壁板和铆钉能量变化时程曲线

    Figure  20.  Energy curves of the back wall and the rivet

  • [1] 纪杨子燚, 李向东, 周兰伟, 等. 高速侵彻体撞击充液容器形成的液压水锤效应研究进展 [J]. 振动与冲击, 2019, 38(19): 242–252. DOI: 10.13465/j.cnki.jvs.2019.19.036.

    JI Y Z Y, LI X D, ZHOU L W, et al. Review of study on hydrodynamic ram effect generated due to high-velocity penetrator impacting fluid-filled container [J]. Journal of Vibration and Shock, 2019, 38(19): 242–252. DOI: 10.13465/j.cnki.jvs.2019.19.036.
    [2] ADDESSIO F L, SCHRAAD M W, LEWIS M W. Physics-based damage predictions for simulating testing and evaluation (T and E) experiments: LA-UR-99-484 [R]. New Mexico: Los Alamos National Laboratory, 1999.
    [3] BALL R E. Structural response of fluid containing tanks to penetrating projectiles (Hydraulic Ram): a comparison of experimental and analytical results: NPS-57BP76051 [R]. Monterey: Naval Postgraduate School, 1976.
    [4] DISIMILE P J, SWANSON L A, TOY N. The hydrodynamic ram pressure generated by spherical projectiles [J]. International Journal of Impact Engineering, 2009, 36(6): 821–829. DOI: 10.1016/j.ijimpeng.2008.12.009.
    [5] VARAS D, ZAERA R, LÓPEZ-PUENTE J. Numerical modelling of the hydrodynamic ram phenomenon [J]. International Journal of Impact Engineering, 2009, 36(3): 363–374. DOI: 10.1016/j.ijimpeng.2008.07.020.
    [6] REN P, ZHOU J Q, TIAN A L, et al. Experimental investigation on dynamic failure of water-filled vessel subjected to projectile impact [J]. International Journal of Impact Engineering, 2018, 117: 153–163. DOI: 10.1016/j.ijimpeng.2018.03.009.
    [7] 李营, 张玮, 杜志鹏, 等. 球形弹体打击作用下宽距水间隔铝板的动态响应特性 [J]. 振动与冲击, 2018, 37(1): 106–110. DOI: 10.13465/j.cnki.jvs.2018.01.017.

    LI Y, ZHANG W, DU Z P, et al. Dynamic responses of wide interval water-spacing aluminum plates under sphere projectile impact [J]. Journal of Vibration and Shock, 2018, 37(1): 106–110. DOI: 10.13465/j.cnki.jvs.2018.01.017.
    [8] 陈安然, 李向东, 周兰伟, 等. 液压水锤效应引起液体喷溅特性及其影响因素试验研究 [J]. 国防科技大学学报, 2021, 43(5): 144–152. DOI: 10.11887/j.cn.202105017.

    CHEN A R, LI X D, ZHOU L W, et al. Experimental study on the characteristics and influencing factors of liquid spurt caused by hydrodynamic ram [J]. Journal of National University of Defense Technology, 2021, 43(5): 144–152. DOI: 10.11887/j.cn.202105017.
    [9] 李亚智, 陈钢. 充液箱体受弹丸撞击下动态响应的数值模拟 [J]. 机械强度, 2007, 29(1): 143–147. DOI: 10.3321/j.issn:1001-9669.2007.01.029.

    LI Y Z, CHEN G. Numerical simulation of liquid-filled tank response to projectile impact [J]. Journal of Mechanical Strength, 2007, 29(1): 143–147. DOI: 10.3321/j.issn:1001-9669.2007.01.029.
    [10] VARAS D, ZAERA R, LÓPEZ-PUENTE J. Numerical modelling of partially filled aircraft fuel tanks submitted to hydrodynamic ram [J]. Aerospace Science and Technology, 2012, 16(1): 19–28. DOI: 10.1016/j.ast.2011.02.003.
    [11] MANSOORI H, ZAREI H. FSI simulation of hydrodynamic ram event using LS-Dyna software [J]. Thin-Walled Structures, 2019, 134: 310–318. DOI: 10.1016/j.tws.2018.10.002.
    [12] 蓝肖颖, 李向东, 周兰伟, 等. 双破片撞击充液容器时液体内压力分布研究 [J]. 振动与冲击, 2019, 38(19): 191–197. DOI: 10.13465/j.cnki.jvs.2019.19.029.

    LAN X Y, LI X D, ZHOU L W, et al. Pressure distribution inside liquid during a liquid-filled vessel impacted by double-fragment [J]. Journal of Vibration and Shock, 2019, 38(19): 191–197. DOI: 10.13465/j.cnki.jvs.2019.19.029.
    [13] 韩璐, 韩庆, 杨爽. 飞机油箱水锤效应影响因素及其影响程度研究 [J]. 航空工程进展, 2018, 9(4): 489–500. DOI: 10.16615/j.cnki.1674-8190.2018.04.005.

    HAN L, HAN Q, YANG S. Simulation analysis of hydrodynamic ram factors and effects in aircraft fuel tank [J]. Advances in Aeronautical Science and Engineering, 2018, 9(4): 489–500. DOI: 10.16615/j.cnki.1674-8190.2018.04.005.
    [14] 韩璐, 韩庆, 杨爽. 多破片高速冲击下飞机油箱水锤效应数值模拟 [J]. 爆炸与冲击, 2018, 38(3): 473–484. DOI: 10.11883/bzycj-2017-0230.

    HAN L, HAN Q, YANG S. Simulation analysis of hydrodynamic ram in an aircraft fuel tank subjected to high-velocity multi-fragment impact [J]. Explosion and Shock Waves, 2018, 38(3): 473–484. DOI: 10.11883/bzycj-2017-0230.
    [15] 陈钢. 高速弹丸冲击下油箱动态响应的数值模拟 [D]. 西安: 西北工业大学, 2005: 3–20. DOI: 10.7666/d.y843974.

    CHEN G. Numerical simulation of fuel tankresponse to projectile impact [D]. Xi’an: Northwestern Polytechnical University, 2005: 3–20. DOI: 10.7666/d.y843974.
    [16] BALL R E. The fundamentals of aircraft combat survivability: analysis and design [M]. 2nd ed. Reston: AIAA Education, 2003: 667–668. DOI: 10.2514/4.862519.
    [17] HULL B T, SEDALOR T, MIFSUD T. Utilization of hydrodynamic ram simulator to determine the dynamic strength thresholds of structural joints [C]//AIAA Scitech 2019 Forum. San Diego, California: American Institute of Aeronautics and Astronautics, 2019. DOI: 10.2514/6.2019-0524.
    [18] 崔新男, 汪旭光, 王尹军, 等. 爆炸加载下混凝土表面的裂纹扩展 [J]. 爆炸与冲击, 2020, 40(5): 052203. DOI: 10.11883/bzycj-2019-0364.

    CUI X N, WANG X G, WANG Y J, et al. External crack propagation of concrete surface under explosive loading [J]. Explosion and Shock Waves, 2020, 40(5): 052203. DOI: 10.11883/bzycj-2019-0364.
    [19] LIDÉN E, HELTE A. Fracture mechanics of long rod projectiles subjected to oblique moving plates [C]// Proceedings of the 26th International Symposium on Ballistics, 2011: 1736–1747.
    [20] BUYUK M, KURTARAN H, MARZOUGUI D, et al. Automated design of threats and shields under hypervelocity impacts by using successive optimization methodology [J]. International Journal of Impact Engineering, 2008, 35(12): 1449–1458. DOI: 10.1016/j.ijimpeng.2008.07.057.
    [21] NAYAK S K, SINGH A K, BELEGUNDU A D, et al. Process for design optimization of honeycomb core sandwich panels for blast load mitigation [J]. Structural and Multidisciplinary Optimization, 2013, 47(5): 749–763. DOI: 10.1007/s00158-012-0845-x.
    [22] 余海燕, 王友. 5052铝合金冲压成形过程中韧性断裂的仿真研究 [J]. 中国有色金属学报, 2015, 25(11): 2975–2981. DOI: 10.19476/j.ysxb.1004.0609.2015.11.003.

    YU H Y, WANG Y. Bulging simulation of ductile fracture of 5052 aluminum alloy [J]. The Chinese Journal of Nonferrous Metals, 2015, 25(11): 2975–2981. DOI: 10.19476/j.ysxb.1004.0609.2015.11.003.
    [23] 马丽英, 李向东, 周兰伟, 等. 高速破片撞击充不同介质液体容器的数值计算及试验研究 [J]. 振动与冲击, 2018, 37(24): 115–122. DOI: 10.13465/j.cnki.jvs.2018.24.018.

    MA L Y, LI X D, ZHOU L W, et al. Numerical simulation and experimental study on high-speed fragment impact filling different liquid containers [J]. Journal of Vibration and Shock, 2018, 37(24): 115–122. DOI: 10.13465/j.cnki.jvs.2018.24.018.
    [24] MEDINA S F, HERNANDEZ C A. General expression of the Zener-Hollomon parameter as a function of the chemical composition of low alloy and microalloyed steels [J]. Acta Materialia, 1996, 44(1): 137–148. DOI: 10.1016/1359-6454(95)00151-0.
    [25] 纪德洋, 金锋, 冬雷, 等. 基于皮尔逊相关系数的光伏电站数据修复 [J]. 中国电机工程学报, 2022, 42(4): 1514–1522. DOI: 10.13334/j.0258-8013.

    JI D Y, JIN F, DONG L, et al. Data repairing of photovoltaic power plant based on pearson correlation coefficient [J]. Proceedings of the CSEE, 2022, 42(4): 1514–1522. DOI: 10.13334/j.0258-8013.
    [26] 蓝肖颖. 双破片作用下液压水锤叠加效应研究 [D]. 南京: 南京理工大学, 2019: 54–56. DOI: 10.27241/d.cnki.gnjgu.2019.000834.
    [27] CHOU P C, CHEN S. Hypervelocity impact of bumper-protected fuel tanks [J]. Journal of Spacecraft and Rockets, 1970, 7(12): 1412–1418. DOI: 10.2514/3.30183.
  • 加载中
图(20)
计量
  • 文章访问数:  363
  • HTML全文浏览量:  121
  • PDF下载量:  77
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-06-27
  • 修回日期:  2023-05-05
  • 网络出版日期:  2023-06-05
  • 刊出日期:  2023-07-05

目录

/

返回文章
返回