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

回转体高速倾斜入水的流场特性及结构响应

高英杰 孙铁志 张桂勇 尤天庆 殷志宏 宗智

章浪, 赵丰鹏, 张钰忠, 邓勇军, 李继承. 钨纤维增强金属玻璃复合材料的长杆弹斜侵彻/穿甲性能[J]. 爆炸与冲击, 2025, 45(3): 033302. doi: 10.11883/bzycj-2024-0158
引用本文: 高英杰, 孙铁志, 张桂勇, 尤天庆, 殷志宏, 宗智. 回转体高速倾斜入水的流场特性及结构响应[J]. 爆炸与冲击, 2020, 40(12): 123301. doi: 10.11883/bzycj-2020-0014
ZHANG Lang, ZHAO Fengpeng, ZHANG Yuzhong, DENG Yongjun, LI Jicheng. Ballistic performance of tungsten fiber-reinforced metallic glass composite in the long rod oblique penetration/perforation[J]. Explosion And Shock Waves, 2025, 45(3): 033302. doi: 10.11883/bzycj-2024-0158
Citation: GAO Yingjie, SUN Tiezhi, ZHANG Guiyong, YOU Tianqing, YIN Zhihong, ZONG Zhi. Flow characteristics and structure response of high-speed oblique water-entry for a revolution body[J]. Explosion And Shock Waves, 2020, 40(12): 123301. doi: 10.11883/bzycj-2020-0014

回转体高速倾斜入水的流场特性及结构响应

doi: 10.11883/bzycj-2020-0014
基金项目: 国家自然科学基金(51639003,51709042);中央高校基本科研业务费专项(DUT2017TB05);中国博士后科学基金(2018M631791,2019T120211);工信部高技术船舶科研项目(2017-614);辽宁省自然科学基金(20180550619);辽宁省“兴辽英才计划”(XLYC1908027);海洋工程国家重点实验室开放基金(1803)
详细信息
    作者简介:

    高英杰(1994- ),男,硕士研究生,yingjie_gao@foxmail.com

    通讯作者:

    孙铁志(1986- ),男,博士,副教授,suntiezhi@dlut.edu.cn

  • 中图分类号: O352

Flow characteristics and structure response of high-speed oblique water-entry for a revolution body

  • 摘要: 回转体高速入水过程涉及液体和固体的耦合作用,是一个复杂的非线性、非定常过程。为研究回转体高速入水的结构动响应及流场演变规律,本文中基于STAR-CCM+和ABAQUS平台,建立了回转体高速入水的双向流固耦合数值模型,开展了不同入水速度的回转体高速倾斜入水流固耦合数值计算。结果表明:数值计算的入水速度、位移曲线和空泡形态与实验结果良好吻合,验证了流固耦合方法的有效性;回转体倾斜高速入水的载荷先集中在触水部分边缘处,后集中于回转体底部中心处;流固耦合方法的入水冲击载荷峰值小于刚体的,弹性回转体的载荷曲线产生明显波动;撞水阶段,回转体空泡呈现不对称形态,随着入水加深,空泡不对称性变弱;入水速度60 m/s下,空泡发生表面闭合,回转体入水初速度越快,空泡表面闭合越晚;冲击载荷与入水速度有关,入水速度越大,峰值出现越早,震荡越明显,速度超过100 m/s时,回转体产生塑性形变。
  • 动能穿甲弹的穿甲能力取决于材料强度、质量以及撞击速度等多种因素。在现代穿甲技术中,钨合金和贫铀合金是2种常用的弹芯材料,然而它们在应用过程中均存在不足:钨合金弹体在穿甲过程中因材料具有良好的塑性,弹头易形成类似“蘑菇头”的膨胀结构,这种结构会增大穿甲阻力,进而削弱穿甲效果[1];贫铀合金弹体则常发生“自锐”现象[2],即弹头前端保持尖锐,但贫铀材料伴随的放射性限制了其应用[3]。鉴于上述瓶颈,20余年来,军事领域一直致力于探索对人和环境无放射性危害且具有“自锐”特性的新型动能穿甲弹芯材料,其中钨纤维增强金属玻璃复合材料具有巨大潜力[4]。该复合材料弹体得益于其金属玻璃(metallic glass,MG)基体的高剪切敏感性[5-6],在侵彻/穿甲过程中弹头易呈现出尖锐构型,即“自锐”行为,其侵彻/穿甲能力得到明显提升[7-21]

    一般而言,侵彻/穿甲试验成本高昂,且试验难以捕捉撞击过程中弹靶的变形和破坏细节,也难以测量高速冲击条件下弹体的应力-应变情况。随着计算机技术和材料动态本构模型的发展,有限元模拟已成为侵彻/穿甲试验的良好补充,有助于更深入的机理研究[15-16, 22-24]。Li等[25-26]引入材料微结构自由体积参量,以准确描述金属玻璃基体特有的高强度和高剪切敏感性,进一步发展了热力耦合本构模型,并根据钨纤维增强金属玻璃复合材料的细观结构建立了弹体有限元模型,基于相关穿甲试验,开展了复合材料长杆弹在正侵彻条件下侵彻/穿甲钢靶的有限元模拟,分析了复合材料弹体的“自锐”机理,并系统地探讨了多个关键因素对复合材料弹体“自锐”行为的影响。

    在实际工程应用中,动能穿甲弹受发射条件和装甲结构设计等因素的影响,对目标的打击多为斜撞击而非理想的正撞击。在斜撞击过程中,弹体受到靶板的非对称作用力,其弹道特征可能出现偏转甚至跳弹现象[27-28]。然而,对于钨纤维增强金属玻璃复合材料长杆弹斜撞击靶板过程中,弹体在非对称作用力下的变形和破坏特征及其“自锐”特性将可能与正撞击情形截然不同,试验[9]也显示,相较于正撞击,斜撞击条件下复合材料弹体的侵彻/穿甲性能降低,有必要开展进一步研究。

    基于Li等[15-16, 22-23]、陈建良等[24]提出的有限元模拟方法,结合相关侵彻/穿甲试验,针对钨纤维增强金属玻璃复合材料和钨合金材料2种材质的长杆弹开展斜侵彻/穿甲钢靶的三维有限元模拟分析,研究2种材质弹体在斜侵彻过程中的弹靶变形和破坏特征,并对比复合材料弹体在正/斜侵彻条件下的“自锐”特征和弹道行为,分析不同撞击速度下撞击倾角对复合材料长杆弹“自锐”行为、侵彻/穿甲性能及相应弹道特征等的影响。

    参考相关穿甲试验[9],采用弹道炮发射8 mm×88 mm的复合材料弹体和钨合金弹体,正/斜侵彻50 mm厚的30CrMnMo钢靶。复合材料弹体的形貌如图1所示,弹体头部为半球形,弹体中增强钨纤维体积约占80%,直径约为0.3 mm。弹体斜侵彻/穿甲靶板的初始状态如图2所示,其中:D为弹体直径;L为弹体长度;v0为弹体的初始速度;靶板倾角θ为靶体表面法线方向与弹体初始速度方向的夹角;H为靶板厚度(试验中,H=50 mm);S为斜侵彻/穿甲时的靶板有效厚度,S=H/cosθ

    图  1  穿甲试验中复合材料弹体的初始形貌[9]
    Figure  1.  Initial composite long rod in penetration test[9]
    图  2  弹靶结构与撞击姿态示意图
    Figure  2.  Schematic diagram of the projectile and target structures and the impact attitude

    对于复合材料弹体正侵彻/穿甲靶板,二维数值模拟能够描述弹靶变形和破坏的主要特征[15-16, 22-24]。针对斜撞击条件,本文中,将开展三维模拟分析。考虑到弹靶结构和撞击条件的平面对称特征,为减小计算规模,数值模型将采用1/2弹靶结构。复合材料弹体的三维有限元模型如图3(a)所示,其中:MG matrix为金属玻璃基体,WF为钨纤维。弹体和靶板网格均采用八节点六面体单元进行划分,网格尺寸约为75 μm,网格之间以共节点的方式进行固接,弹身横截面的网格划分如图3(b)所示。Zhang等[29]的研究表明,钨纤维与金属玻璃基体之间的界面结合良好,扫描电子显微镜(scanning electron microscope,SEM)直观地显示了实际复合材料的横截面(图3(c))。由图3(a)~(c)可知,所构建的有限元模型能够准确描述复合材料的细观结构特征,验证了模型的准确性。

    图  3  复合材料长杆弹的三维有限元模型
    Figure  3.  3D finite element geometrical model of a composite long rod

    对钨合金弹体进行数值模拟时,将复合材料弹体模型中的基体材料替换为钨合金材料即可。

    1.2.1   金属玻璃材料

    Li等[25-26]推导的三维修正热力耦合本构模型综合考虑了多种因素对金属玻璃变形和破坏的影响,如应变率、自由体积浓度、温度、应力以及静水应力等,从而能够准确描述金属玻璃材料在不同应变率和初始温度下的力学行为特性。试验弹体中,金属玻璃基体的本构模型参数如表1所示,其中:˙ε为应变率,Λc为压缩条件下的静水应力敏感因子,Λt为拉伸条件下的静水应力敏感因子。

    表  1  锆基金属玻璃的修正热力耦合模型参数[25-26]
    Table  1.  Parameters in the modified coupled thermo-mechanical constitutive model for Zr-based metallic glass[25-26]
    参量 符号 单位 数值
    弹性模量 E GPa 96
    泊松比 ν 0.36
    密度 ρ kg/m3 6125
    熔化温度 Tm K 993
    玻璃转变温度 Tg K 625
    初始温度 T0 K 300
    比定容热容 cV J/(kg·K) 400
    临界体积 v* m3 2.0×10−29
    平均原子体积 Ω m3 2.5×10−29
    原子振动频率 f s−1 1×1013
    临界破坏自由体积浓度 ξc 0.065
    初始自由体积浓度 ξ0 0.05
    运动激活能 ΔGm eV ΔGm(˙ε)
    几何因子 α 0.05
    所需跃迁次数 nD 3
    静水应力敏感因子 Λ 0.05(Λc
    0.35(Λt
    下载: 导出CSV 
    | 显示表格
    1.2.2   金属材料

    30CrMnMo钢靶、钨合金弹体以及复合材料弹体中的增强钨纤维等均为常规的晶态合金材料。为了描述高应变率下合金材料的力学行为,综合Johnson-Cook本构模型[30]和累积损伤失效模型[31]来表征金属材料的力学性能,采用Grüneisen状态方程[32]来描述穿甲过程中金属材料的压力状态。钨合金和30CrMnMo钢的Johnson-Cook模型参数如表2所示,其中:˙ε0为参考应变率,Tr为参考温度,A为参考应变率条件下的屈服强度,B为应变硬化系数,C为应变率敏感系数,m为温度敏感系数,n为应变硬化指数,D1D5为材料损伤参数,c0为体积声速,S1vs-vp曲线斜率的系数(vs为冲击波速度,vp为粒子速度),γ0为Grüneisen系数,a为对γ0的一阶体积修正。其他材料模型参数可参考文献[15-16, 22-23]。

    表  2  金属材料的Johnson-Cook模型参数
    Table  2.  Johnson-Cook model parameters of metallic materials
    材料 ρ/(kg·m−3) ν E/GPa ˙ε0/s−1 Tr/K Tm/K cV/(J·kg−1·K−1)
    95W钨合金 17900 0.28 410 1 300 1752 134
    30CrMnMo钢 7850 0.29 200 1 300 1793 477
    材料 A/MPa B/MPa C m n D1 D2
    95W钨合金 1650 450 0.016 1.00 0.12 3.00 0
    30CrMnMo钢 1200 310 0.014 1.03 0.26 3.20 0
    材料 D3 D4 D5 c0/(m·s−1) S1 γ0 a
    95W钨合金 0 0 0 3850 1.44 1.58 0
    30CrMnMo钢 0 0 0 4578 1.38 1.67 0.47
    下载: 导出CSV 
    | 显示表格

    Li等[15]对文献[9]中的弹体正侵彻进行了模拟分析,本文中,对文献[9]中的2种斜侵彻试验进行数值模拟,试验和模拟结果如表3所示,其中:带“*”号数据是从文献中的图直接测量得到。

    表  3  侵彻试验[9]以及相应的数值模拟结果
    Table  3.  Penetration test[9] and the corresponding simulation results
    弹材 撞击速度/
    (m·s−1)
    倾角/
    (°)
    侵彻深度/mm 剩余弹长/mm
    试验[9] 模拟 试验[9] 模拟
    钨合金 852.9 0 23.4 22.2[15] 14.1* 13.2[15]
    1076.2 0 背面鼓包 47.7[15] 14.1[15]
    1235.1 50 背面鼓包 49.7 12.8
    复合材料 857.5 0 36.4 32.3[15] 18.8* 17.1[15]
    1066.3 0 穿透 穿透[15] 28.0* 26.0[15]
    1263.9 50 穿透 穿透 14.5
    下载: 导出CSV 
    | 显示表格

    倾角θ=50°、撞击速度v0=1235.1 m/s时,钨合金弹斜侵彻靶板的最终弹靶变形和破坏形貌如图4所示,其中:图4(a)~(c)分别为靶板纵截面、正面和背面的数值模拟结果,图4(d)~(e)分别为实际靶板正面和背面的最终形貌。可以看出,数值模拟结果中,靶板前端的椭圆形开孔和花瓣形撕裂形貌、靶板后端的隆起鼓包形貌等均与试验结果较为相符,侵彻深度和剩余弹长的模拟值与试验值也较一致(表3)。由于钨合金弹体具有良好的塑性,在侵彻过程中其头部发生镦粗变形,导致靶板内形成的弹孔直径相较于原始弹体直径明显增大,特别是在靶板的开坑区域,由于倾斜撞击作用,靶板的开孔口径显著增大,约为弹径的3倍,且靶孔呈现显著的非对称特征,下端孔口相对光滑而上端孔口发生显著的花瓣形撕裂(图4(a)~(b))。另外,受靶板两侧的挤压和摩擦作用,弹体发生质量侵蚀,即弹体碎片自弹头外边缘剥落并向后端流出,碎片的外流也表现出非对称特征,主要沿靶孔上表面流出(图4(a))。由于弹头镦粗变形导致较大的侵彻阻力,钨合金弹体未能穿透靶板,而是在靶板背面形成较大鼓包,其侵彻能力与钨合金弹在θ=0°(正侵彻)、v0=1076.2 m/s条件下的侵彻能力相近(表3);弹体侵蚀也较严重,剩余弹体长度仅为12.8 mm。特别地,由于斜侵彻(θ=50°),弹体在侵彻后期临近穿透靶板时,靶板背面施加于弹体的侵彻阻力并不对称,导致弹道轨迹偏转,同时弹体自身发生明显弯曲,弹头的镦粗变形也呈现出非对称特征(图4(a))。

    图  4  钨合金弹斜侵彻30CrMnMo钢靶的弹靶变形和破坏形貌(θ=50°,v0=1235.1 m/s)
    Figure  4.  Deformation and failure morphologies of projectile and target materials after tungsten alloy long rod penetrating a 30CrMnMo steel target (θ=50°, v0=1235.1 m/s)

    θ=50°、v0=1263.9 m/s时,复合材料弹体斜侵彻钢靶的最终弹靶变形和破坏形貌如图5所示,结合表3中的侵彻深度和剩余弹长可知,数值模拟与试验结果符合较好。与图4中的钨合金弹斜侵彻相比,弹体在侵彻初期的镦粗变形程度减弱,进而导致靶板前端的开孔口径减小(图5(a)),上端孔口的花瓣形撕裂程度也减弱;在随后的侵彻过程中,复合材料弹体头部发生显著的剪切“自锐”行为,弹头表现为尖头构型,并导致靶板孔径减小,基本与弹径相当;在侵彻后期临近穿透靶板时,受靶板背面的非对称侵彻阻力作用,弹道轨迹偏转,弹头的锐化尖头构型也呈现出非对称特征(图5(a)),但弹体自身未弯曲。由于复合材料弹体的“自锐”行为,其受到的靶板阻力小于钨合金弹工况,因此,弹体侵彻性能显著提高,在相近的撞击速度下直接穿透靶板,同时弹体侵蚀程度也有所减小,剩余长度约为14.5 mm,比钨合金弹工况下的剩余弹长(12.8 mm)长。

    图  5  复合材料弹斜侵彻30CrMnMo钢靶的弹靶变形和破坏形貌(θ=50°,v0=1263.9 m/s)
    Figure  5.  Deformation and failure morphologies of projectile and target materials after a composite long rod penetrating a 30CrMnMo steel target (θ=50°, v0=1263.9 m/s)

    钨合金弹和复合材料弹在不同撞击倾角和速度下的残余弹体形貌如图67所示。可以看出,在正侵彻条件下,2种材质弹体的后端部分均较好地保持了初始形态,未出现明显变形;对于弹体头部,钨合金弹头镦粗为“蘑菇头”形状,临近弹头的弹身直径也明显增大(图6(a)~(b)和(d)),而复合材料弹头在“边缘层”[11, 14]区域锐化为尖头结构,且尖头结构后端区域的弹身与初始形貌基本相同(图7(a)~(b)和(d))。在斜侵彻条件下,因为侵彻后期的弹道轨迹偏转,2种弹体的后端区域均与靶板发生碰撞和摩擦作用,进而产生变形和破坏,钨合金弹体甚至发生了一定程度的弯曲(图6(c));对于弹体头部,钨合金弹头也镦粗为“蘑菇头”形状(图6(c)),复合材料弹头锐化为尖头构型(图7(c)),且二者均呈现出显著的非对称特征,与正侵彻工况明显不同,结合图4(a)和图5(a)可知,钨合金弹和复合材料弹头部的变形主要发生在与靶孔下表面接触的区域,即弹头的非对称变形和破坏形貌主要源于靶板作用的非对称特征。

    图  6  钨合金弹在不同撞击倾角和速度下的残余弹体形貌
    Figure  6.  Residual morphology of tungsten alloy projectiles after the impact at different oblique angles and initial velocities
    图  7  复合材料弹在不同撞击倾角和速度下的残余弹体形貌
    Figure  7.  Residual morphology of composite projectiles after the impact at different oblique angles and initial velocities

    本节进一步分析斜侵彻/穿甲条件下复合材料弹体的变形和破坏特征,详细讨论其独特的“自锐”行为及作用机制,并与钨合金弹进行对比。

    钨合金弹(θ=50°,v0=1235.1 m/s)斜侵彻50 mm厚的30CrMnMo钢靶时,弹靶内部的塑性应变发展历程如图8所示。可以看出,在斜侵彻的开坑阶段,由于钨合金具有较高的塑性,并且弹体受到靶板施加的非对称侵彻阻力,因此,弹头材料主要向阻力较小的上侧发生侵蚀和流动,初始的半球形弹头发生非对称镦粗变形,其上侧区域的镦粗程度显著大于下侧区域,导致靶板前端的靶孔直径远大于弹径,也呈现出显著的非对称特征(图8(a))。随后,弹头完全侵入靶板,此时弹体受到的靶板阻力较为对称,弹头的镦粗逐渐呈现出相对对称的“蘑菇头”形状,但形貌的对称程度较正侵彻/穿甲工况(图6(b))有所降低,该侵彻过程相对稳定,受靶板两侧的挤压作用,弹头材料从镦粗的“蘑菇头”两侧侵蚀剥落并向弹体后端流失,靶板孔径较初期开坑阶段明显减小(图8(b))。弹体临近靶板后端时,由于弹体前端的靶体材料厚度显著降低,使得靶板的变形模式转变为薄靶的整体响应,靶板背面开始出现鼓包变形;类似地,由于靶板倾斜放置,靶板后端未侵蚀区域的厚度不均匀,后端的鼓包变形不对称,靶体材料施加到弹体的侵彻阻力也不对称,弹道轨迹向侵彻阻力较小的下侧发生偏转,同时弹体发生一定程度的弯曲变形(图8(c)~(d)),弹头的镦粗变形再次呈现出非对称特征(图6(c))。相应地,弹体弯曲变形和弹道偏转进一步增大了弹体的侵彻阻力,且靶板倾斜放置增加了有效厚度,弹体最终未能穿透靶板,弹体侵蚀也更为严重,剩余弹长进一步减小(表3)。

    图  8  钨合金弹斜侵彻靶板时弹靶的塑性应变发展历程(θ=50°,v0=1235.1 m/s)
    Figure  8.  Development of effective plastic strain in the projectile and target materials during the penetration of tungsten alloy long rod (θ=50°, v0=1235.1 m/s)

    复合材料弹(θ=50°,v0=1263.9 m/s)斜侵彻50 mm厚的30CrMnMo钢靶时弹靶内部的塑性应变发展历程如图9所示。在斜侵彻开坑阶段,同样受到靶板的非对称侵彻阻力,弹头材料也主要向上侧侵蚀和流动,弹头发生非对称的镦粗变形,导致靶板前端形成相对较大的孔径(图9(a)),但弹头的镦粗程度和靶板的孔径均明显小于钨合金弹工况(图8(a))。主要原因在于,高速冲击下,复合材料弹体中的金属玻璃基体具有极高的强度和剪切敏感性[15-16, 22-23, 28-29],在破坏之前可较好地约束钨纤维的变形,而在损伤破坏之后,其局域化剪切带或剪切裂纹对钨纤维产生强烈冲击,导致弹头前端形成一个“边缘层”[15-16],“边缘层”内的钨纤维受靶板的侧向挤压而迅速发生剪切破坏并从弹体剥落,因此,弹头的镦粗程度相对较弱。在随后的弹头完全侵入靶板的过程中,靶板的阻力逐渐对称,弹头的非对称变形和破坏特征减弱,弹头逐渐“自锐”为尖头结构(图9(b)),但形貌的对称程度较正侵彻/穿甲工况(图7(b))有所降低,相应地,侵彻过程中靶板的孔径逐渐与弹径相当。在弹体临近靶板后端时,靶板背面也发生不对称的鼓包变形,导致弹道轨迹向靶板下侧偏转,同时弹体发生弯曲变形,且弹头的“自锐”变形再次呈现非对称特征(图9(c)~(d)),由于弹体内金属玻璃基体的高强度和高硬度,弹体的弯曲程度较钨合金弹工况有所降低(图8(c))。由于弹头锐化导致侵彻阻力较低,弹体的侵彻能力提升,最终穿透靶板(图5(a)),弹体侵蚀程度严重,剩余弹长明显减小(表3),弹体最终的尖头形貌呈现出显著的非对称特征(图7(c))。

    图  9  复合材料弹斜侵彻靶板时弹靶的塑性应变发展历程(θ=50°,v0=1263.9 m/s)
    Figure  9.  Development of effective plastic strain in the projectile and target materials during the penetration of composite long rod (θ=50°, v0=1263.9 m/s)

    图10为2种弹体斜侵彻后靶孔最终的轮廓形貌。可以看出:对于2种弹体,靶板的弹道形貌均表现出靶板前端面区域的开坑和临近后端面区域的弹道偏转特性;在整个侵彻过程中,钨合金弹形成的靶板孔径均大于复合材料弹形成的靶板孔径,前者约为后者的1.4倍,钨合金弹体的侵彻深度小于复合材料弹体的侵彻深度。相关差异也主要源于钨合金弹体的头部镦粗变形和复合材料弹体的“自锐”行为。

    图  10  不同弹体斜侵彻/穿甲后靶板孔洞的轮廓形貌(θ=50°)
    Figure  10.  Contour morphologies of penetrating hole in the target after the impact of different long rods (θ=50°)

    在弹体侵彻/穿甲过程中,弹体速度是靶板对弹体产生阻力的一个重要参量。图11给出了2种弹体在不同穿甲条件下的速度变化曲线。由图11表3可知:撞击倾角一致时,侵彻前期,2种弹体的速度变化几乎一致;但在侵彻后期,钨合金弹体的速度下降明显快于复合材料弹体,这主要归因于复合材料弹的头部展现出的“自锐”特性,使得弹体在侵彻时所受阻力显著小于弹头呈镦粗状态的钨合金弹体。另外,相较于正侵彻/穿甲,斜侵彻/穿甲工况下弹体速度下降更快,这主要源于在斜侵彻/穿甲过程中弹头的非均匀变形和弹道弯曲特征。

    图  11  不同弹体在不同侵彻/穿甲条件下的速度变化曲线
    Figure  11.  Variation of velocity of different projectiles under different long rods and impact conditions

    综上所述,长杆弹斜撞击时,在开坑阶段,受到靶板的非对称阻力,弹头材料向阻力较小的一侧侵蚀和流动,进而导致弹头发生非对称的镦粗变形;在随后相对稳定的侵彻过程中,弹头的非对称变形和破坏特征逐渐减弱,且复合材料弹头开始发生“自锐”行为;在侵彻后期,弹道轨迹向靶板阻力较小的一侧偏转,同时弹体发生弯曲变形,弹头结构的非对称特征再次增强。相较于正撞击,斜侵彻/穿甲工况下,由于受力的非对称性,复合材料弹头的“自锐”特征减弱,易发生非对称变形和破坏,弹体侵蚀和弹道偏转明显加剧。在相同斜撞击工况下,复合材料弹的头部镦粗、弯曲变形和弹道偏转等均小于钨合金弹,其侵彻/穿甲性能仍高于钨合金弹。

    如上所述,复合材料弹在斜侵彻/穿甲过程中将受到非对称的靶板阻力,进而导致与正侵彻/穿甲情形显著不同的弹头“自锐”特征和弹道特性。不同的撞击倾角下,弹体在开坑阶段和后期穿靶阶段受到靶板上下两侧所施加的阻力差异将有所不同,这可能使得弹靶的变形和破坏特征以及弹道偏转产生显著差别,最终导致不同的侵彻/穿甲性能。为了分析撞击倾角的影响,选取v0=857.5 m/s为典型中低速撞击速度,v0=1263.9 m/s为典型高速撞击速度,模拟撞击倾角θ在0°~70°范围内复合材料弹体撞击钢靶的全过程。

    θ=30°时,2种典型速度下复合材料弹体撞击钢靶的弹靶变形和破坏历程如图12所示。由图12(a)~(c)可知,高速侵彻时,因基体材料的屈服强度较高,基体材料在破坏之前对钨纤维形成较强约束,使得在开坑阶段,弹头区域材料的侵蚀较小,基体材料破坏后其内部剪切带或剪切裂纹快速切断临近钨纤维,使其向后端流失,因此,弹头的镦粗变形并不明显(图12(a))。在后续侵彻过程中,随着靶板上下两侧的挤压作用逐渐对称,弹头的镦粗变形程度降低,结构的非对称性程度减弱(图12(b))。在后期弹体临近靶板后端时,由于靶板的倾斜放置,靶板后端的隆起鼓包呈非对称性,靶板两侧对弹头的挤压作用也不对称,使得弹头上侧表现为剪切破坏,而下侧钨纤维难以被切断,弹头又呈现出一定程度的非对称镦粗变形(图12(c))。此外,在高速侵彻过程中,侵彻初期弹体头部承受非对称阻力作用的时间较短,并且弹头镦粗程度较低,导致靶板阻力较小,因此,弹道的偏转程度较低,总体来说相对笔直(图12(c))。

    图  12  复合材料弹体在不同速度下撞击钢靶的弹靶变形和破坏历程(θ=30°)
    Figure  12.  Deformation and failure process of composite long rod projectiles and targets during the penetration at different impact velocities (θ=30°)

    图12(d)~(f)所示,中低速侵彻时,金属玻璃基体对钨纤维的前期约束作用和后期冲击剪切作用减弱,使得弹头明显镦粗,同时弹头承受非对称阻力作用的时间较长,弹头的镦粗变形呈现出显著的非对称特征,且弹道也开始发生偏转(图12(d))。在后续侵彻过程中,弹头的镦粗程度和非对称特征减弱,但弹道持续偏转(图12(e))。在后期,弹头发生“自锐”行为,锐化的尖头构型呈现非对称特征,弹道继续偏转,弹体整体侵蚀程度严重(图12(f))。相较于高速侵彻,中低速侵彻时,复合材料弹体的侵彻能力显著降低。

    θ取50°和70°、v01263.9和857.5 m/s时,复合材料弹体撞击钢靶的最终弹靶变形和破坏形貌如图13所示。由图13(b)可知,θ=50°、v0=857.5 m/s时,弹道严重偏转,弹体严重侵蚀破坏,未形成有效侵彻。θ=70°工况下,弹体的侵彻性能进一步衰减,v0=1263.9 m/s时,弹道严重偏转(图13(c)),v0=857.5 m/s时,弹体直接跳飞(图13(d))。

    图  13  复合材料弹体在不同倾角和速度条件下撞击钢靶的最终弹靶变形和破坏形貌
    Figure  13.  Final deformation and failure morphologies of composite long rod projectiles and targets after the impact at different oblique angles and initial velocities

    图14给出了复合材料弹体以不同撞击速度和撞击倾角侵彻/穿甲钢靶后靶板孔洞的轮廓形貌,相应的剩余弹体形貌和弹体速度变化历程曲线如图1516所示。θ=70°时,弹体未能形成有效侵彻甚至发生跳飞,因此,不列出相关结果。由图14(a)可知,v0=1263.9 m/s时,随着撞击倾角的增大,靶板前端的开孔口径逐渐增大,孔洞的非对称特征也更加明显,表明随着撞击倾角的增大,弹头材料加剧向靶孔上侧侵蚀和流动,弹头的镦粗变形程度增强;在后续侵彻过程中,不同撞击倾角工况下靶板的孔径均逐渐减小,并最终与弹径一致,弹道相对笔直,表明弹体均发生了“自锐”行为;在侵彻后期弹体接近靶板后端面时,弹道轨迹向下偏转逐渐加剧,后端靶孔的非对称特征更加明显。v0=857.5 m/s(图14(b))时,靶板前端的开坑孔径随着倾角的增加而不断增大;在后续侵彻过程中,不同撞击倾角工况下的靶孔均逐渐变小,最终的靶孔底端呈现出尖锐的构型特征,构型的非对称特征随着倾角的增加而增强,弹道偏转也更为显著;当倾角达到50°以上时,弹体已无法形成有效侵彻。

    图  14  复合材料弹体以不同速度和倾角侵彻/穿甲钢靶后靶板孔洞的轮廓形貌
    Figure  14.  Contour morphologies of target plate holes after composite long rod impact at different initial velocities and oblique angles
    图  15  复合材料弹体以不同速度和倾角侵彻/穿甲钢靶后残余弹体的形貌
    Figure  15.  Residual morphology of projectiles after the penetration/perforation of composite long rods at different initial velocities and oblique angles
    图  16  复合材料弹体以不同速度和倾角侵彻/穿甲钢靶时弹体的速度变化曲线
    Figure  16.  Variation of velocity of composite long rods at different initial velocities and oblique angles

    图15(a)可以看出,在v0=1263.9 m/s工况下,当θ<50°时,在侵彻后期临近穿透靶板时,弹体均受薄靶响应特征影响而发生一定的镦粗变形,且弹头的镦粗程度和构型的非对称特征均随着撞击倾角的增大而增强。当θ=50°时,受较大的靶板倾斜角度和弹体后期较低的穿透靶板速度影响,弹头又呈现出非对称的尖锐头型特征;此外,倾角增加使得靶板的有效厚度增大,随着倾角的增加,弹体侵蚀更为严重,剩余长度逐渐减小,且减小幅度更为明显。由图15(b)可知,v0=857.5 m/s工况下,随着撞击倾角的增大,弹体下侧面的侵蚀和弹体的弯曲变形程度加剧,剩余弹体头形的非对称特征也增强,但剩余弹长基本相同,主要是因为弹道偏转加剧,弹体的侵彻路程差距较小,弹体的侵蚀也基本相同。另外,对比图15(a)和图15(b)可知,在θ<50°范围内,相同撞击倾角工况下,撞击速度越低,弹体的侵蚀程度越严重,与雷波[9]的试验结果相同,这也源于金属玻璃基体强度和剪切敏感性的应变率相关性[15-16, 22-23, 25-26];当θ=50°时,由于v0=857.5 m/s时弹体未能形成有效侵彻,剩余弹长反而相对较长。

    图16(a)可知,高速侵彻时,随着撞击倾角的增大,弹体速度的下降幅度增大,但变化特征基本一致。这主要是因为,在靶板倾斜放置产生非对称阻力的条件下,弹头仍能较好地发生“自锐”行为(图15(a)),且弹道相对笔直(图14(a))。由图16(b)可知,中低速侵彻时,相较于正侵彻,θ<40°范围内弹体速度下降更快,这是受弹头“自锐”能力降低、弹道偏转(图14(b))以及弹身与靶板孔壁相互作用(图15(b))等因素的综合影响;当θ≥40°时,随着撞击倾角的增大,弹体速度减速变小,这是因为弹道偏转严重,弹体主要在靶板的前端面附近穿行,受到的靶板总体阻力变小,此时弹体已难以形成有效侵彻。

    综上所述,在侵彻初期开坑阶段,不同撞击速度下,随着撞击倾角的增大,由于弹体受力的非对称性,复合材料弹体头部的“自锐”行为更难发生;在侵彻后期弹体临近穿透靶板时,高速侵彻工况下,随着倾角的增大,弹道偏转程度以及弹头变形和破坏的非对称特征加剧,弹体侵蚀更严重。相应地,弹体的侵彻/穿甲性能随倾角增加而逐渐降低,弹体可形成有效撞击的速度阈值也逐渐提高,当倾角达到50°后,弹体在900 m/s速度范围内难以形成有效侵彻。总体来说,撞击倾角对复合材料弹体的侵彻/穿甲主要起到消极影响,实际工程应用中应尽量减小撞击倾角。

    基于细观有限元模拟,研究了钨纤维增强金属玻璃复合材料长杆弹对50 mm厚30CrMnMo钢靶的斜侵彻/穿甲性能,并与钨合金长杆弹进行对比分析,讨论了撞击倾角对复合材料弹体的“自锐”特性、穿甲性能以及相应弹道行为的影响,得到以下主要结论。

    (1) 复合材料长杆弹对钢靶的斜侵彻/穿甲行为与正撞击行为完全不同。在弹体斜撞击开坑阶段,受靶板的非对称阻力作用,弹头发生非对称的镦粗变形;在后续侵彻过程中,弹头的非对称变形和破坏特征减弱,并发生“自锐”行为;在侵彻后期,非对称的靶板阻力使弹体弯曲,弹头结构又呈现非对称特征,弹道轨迹偏转。相较于正撞击,斜撞击时,复合材料弹体的“自锐”性能降低,弹体侵蚀和弹道偏转程度加剧,弹体速度加速下降,弹体的侵彻/穿甲性能降低,但仍高于钨合金弹的侵彻/穿甲性能。

    (2) 随着撞击倾角的增大,在弹体初期开坑阶段弹头的“自锐”行为更难发生,在侵彻后期弹道的偏转以及弹头的非对称变形和破坏特征加剧,弹体侵蚀更严重,弹体侵彻/穿甲性能降低, 甚至难以形成有效侵彻。

    总体来说,撞击倾角不利于复合材料弹体的侵彻/穿甲,在实际工程应用中,应尽量减小撞击倾角。

  • 图  1  双向耦合求解过程

    Figure  1.  Computing process of fluid-structure interaction method

    图  2  计算域

    Figure  2.  Computational domain

    图  3  流体网格模型

    Figure  3.  Mesh of fluid domain

    图  4  速度衰减曲线

    Figure  4.  Velocity attenuation curves

    图  5  侵蚀位移曲线

    Figure  5.  Erosion displacement curves

    图  6  不同时刻的空泡形态

    Figure  6.  Cavity features at different times

    图  7  回转体轴向、底部径向的单元分布

    Figure  7.  Element distribution of revolution body

    图  8  回转体的应力沿轴向分布曲线

    Figure  8.  Equivalent stress curves of elements along axis

    图  9  不同时刻回转体底部应力分布

    Figure  9.  Stress distributions of revolution body bottom at different times

    图  10  回转体底部的应力沿径向分布曲线

    Figure  10.  Equivalent stress curves of bottom elements along radius

    图  11  回转体底部的应变沿径向分布曲线

    Figure  11.  Strain curves of bottom elements along radius

    图  12  不同时刻的空泡演变

    Figure  12.  Cavity features at different times

    图  13  不同时刻的应力分布

    Figure  13.  Stress distributions at different times

    图  14  速度60 m/s倾斜入水的空泡演变

    Figure  14.  Cavity features for oblique water-entry at v0=60 m/s

    图  15  入水速度80 m/s时回转体底部中心点压力曲线

    Figure  15.  Pressure curve of bottom central point at v0=80 m/s

    图  16  速度80 m/s入水时速度曲线

    Figure  16.  Velocity curves at v0=80 m/s

    图  17  归一化的速度衰减曲线

    Figure  17.  Normalized velocity attenuation curves

    图  18  回转体中心单元的有效应力曲线

    Figure  18.  Equivalent stress curves of central elements

    图  19  空泡形态对比

    Figure  19.  Comparisons of cavity features

    表  1  不同入水速度的压力峰值

    Table  1.   Peak pressures at different initial water entry velocities

    入水速度/(m·s−1)压力峰值/MPa
    CFDFSI
    60 9.93 9.42
    8016.0315.59
    10027.9225.12
    下载: 导出CSV
  • [1] 王永虎, 石秀华. 入水冲击问题研究的现状与进展 [J]. 爆炸与冲击, 2008, 28(3): 276–282. DOI: 10.11883/1001-1455(2008)03-0276-07.

    WANG Y H, SHI X H. Review on research and development of water-entry impact problem [J]. Explosion and Shock Waves, 2008, 28(3): 276–282. DOI: 10.11883/1001-1455(2008)03-0276-07.
    [2] Von KARMAN T H. The impact on seaplane floats during landing [R]. Washington DC: National Advisory Committee on Aeronautics, 1929: 309−313. DOI: 10.1115/1.4023571.
    [3] WAGNER H. Über stoß-und gleitvorgänge an der oberfläche von flüssigkeiten [J]. Zeitschrift für Angewandte Mathematik und Mechanik, 1932, 12(4): 193–215. DOI: 10.1002/zamm.19320120402.
    [4] GAVRILENKO V N, KUBENKO V D. Plane problem of rigid body penetration into a compressible fluid [J]. Soviet Applied Mechanics, 1985, 21(4): 345–352. DOI: 10.1007/BF00886581.
    [5] WORTHINGTON A M, COLE R S. Impact with a liquid surface studied by the aid of instantaneous photography. Paper Ⅱ [J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1900, 194: 175–199. DOI: 10.1098/rsta.1900.0016.
    [6] MAY A. Effect of surface condition of a sphere on its water-entry cavity [J]. Journal of Applied Physics, 1951, 22(10): 1219–1222. DOI: 10.1063/1.1699831.
    [7] MAY A, WOODHULL J C. The virtual mass of a sphere entering water vertically [J]. Journal of Applied Physics, 1950, 21(12): 1285–1289. DOI: 10.1063/1.1699592.
    [8] ABELSON H I. Pressure measurements in the water-entry cavity [J]. Journal of Fluid Mechanics, 1970, 44(1): 129–144. DOI: 10.1017/S0022112070001738.
    [9] PANAHI R. Simulation of water-entry and water-exit problems using a moving mesh algorithm [J]. Journal of Theoretical and Applied Mechanics, 2012, 42(2): 79–92.
    [10] ERFANIAN M R, ANBARSOOZ M, RAHIMI N, et al. Numerical and experimental investigation of a three dimensional spherical-nose projectile water entry problem [J]. Ocean Engineering, 2015, 104: 397–404. DOI: 10.1016/j.oceaneng.2015.05.024.
    [11] NEAVES M D, EDWARDS J R. All-speed time-accurate underwater projectile calculations using a preconditioning algorithm [J]. Journal of Fluids Engineering, 2006, 128(2): 284–296. DOI: 10.1115/1.2169816.
    [12] 孙辉, 卢炽华, 何友声. 二维楔形体冲击入水时的流固耦合响应的实验研究 [J]. 水动力学研究与进展, 2003, 18A(1): 104–109. DOI: 10.3969/j.issn.1000-4874.2003.01.018.

    SUN H, LU C H, HE Y S. Experimental research on the fluid-structure interaction in water entry of 2D elastic wedge [J]. Journal of Hydrodynamics, 2003, 18A(1): 104–109. DOI: 10.3969/j.issn.1000-4874.2003.01.018.
    [13] 张伟, 郭子涛, 肖新科, 等. 弹体高速入水特性实验研究 [J]. 爆炸与冲击, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06.

    ZHANG W, GUO Z T, XIAO X K, et al. Experimental investigations on behaviors of projectile high-speed water entry [J]. Explosion and Shock Waves, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06.
    [14] 郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [J]. 爆炸与冲击, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07.

    GUO Z T, ZHANG W, GUO Z, et al. Characteristics of velocity attenuation and cavity expansion induced by horizontal water-entry of truncated-ogive nosed projectiles [J]. Explosion and Shock Waves, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07.
    [15] 陈诚, 袁绪龙, 党建军, 等. 超空泡航行器20°角倾斜入水冲击载荷特性试验研究 [J]. 兵工学报, 2018, 39(6): 1159–1164. DOI: 10.3969/j.issn.1000-1093.2018.06.016.

    CHEN C, YUAN X L, DANG J J, et al. Experimental investigation into impact load during oblique water-entry of a supercavitating vehicle at 20° [J]. Acta Armamentarii, 2018, 39(6): 1159–1164. DOI: 10.3969/j.issn.1000-1093.2018.06.016.
    [16] YAN G X, PAN G, SHI Y, et al. Experimental and numerical investigation of water impact on air-launched AUVs [J]. Ocean Engineering, 2018, 167: 156–168. DOI: 10.1016/j.oceaneng.2018.08.044.
    [17] 钱铖铖, 余春华, 穆青, 等. 发射速度和发射角度对射弹高速入水流动的影响 [J]. 兵器装备工程学报, 2019, 40(7): 35–39; 50. DOI: 10.11809/bqzbgcxb2019.07.008.

    QIAN C C, YU C H, MU Q, et al. Numerical research of effects of launch speed and launch angle on water entry of high-speed projectile [J]. Journal of Ordnance Equipment Engineering, 2019, 40(7): 35–39; 50. DOI: 10.11809/bqzbgcxb2019.07.008.
    [18] 张佳悦, 李达钦, 吴钦, 等. 航行体回收垂直入水空泡流场及水动力特性研究 [J]. 力学学报, 2019, 51(3): 803–812. DOI: 10.6052/0459-1879-18-364.

    ZHANG J Y, LI D Q, WU Q, et al. Numerical investigation on cavity structures and hyrodynamics of the vehicle during vertical water-entry [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(3): 803–812. DOI: 10.6052/0459-1879-18-364.
    [19] 何春涛, 王聪, 闵景新, 等. 回转体匀速垂直入水早期空泡数值模拟研究 [J]. 工程力学, 2012, 29(4): 237–243.

    HE C T, WANG C, MIN J X, et al. Numerical simulation of early air-cavity of cylinder cone with vertical water-entry [J]. Engineering Mechanics, 2012, 29(4): 237–243.
    [20] 侯昭, 孙铁志, 张桂勇, 等. 回转体倾斜入水空泡试验及六自由度数值计算研究 [J]. 宇航总体技术, 2017, 1(4): 38–45.

    HOU Z, SUN T Z, ZHANG G Y, et al. Experimental investigation and 6-DOF simulation of oblique water-entry cavity of revolution body [J]. Astronautical Systems Engineering Technology, 2017, 1(4): 38–45.
    [21] 黄志刚, 孙铁志, 杨碧野, 等. 平头锥型回转体高速入水结构强度数值分析 [J]. 爆炸与冲击, 2019, 39(4): 043201. DOI: 10.11883/bzycj-2017-0330.

    HUANG Z G, SUN T Z, YANG B Y, et al. Numerical analysis on structural strength of a cone-shaped flatted revolution body during high-speed water-entry [J]. Explosion and Shock Waves, 2019, 39(4): 043201. DOI: 10.11883/bzycj-2017-0330.
    [22] 陈宇, 曹红松, 王智军, 等. 一种新型防空MEFP的设计与仿真 [J]. 现代防御技术, 2013, 41(2): 202–206; 211. DOI: 10.3969/j.issn.1009-086x.2013.02.037.

    CHEN Y, CAO H S, WANG Z J, et al. New type of air defense MEFP [J]. Modern Defense Technology, 2013, 41(2): 202–206; 211. DOI: 10.3969/j.issn.1009-086x.2013.02.037.
    [23] AGMELL M, AHADI A, STÅHL J E. A numerical and experimental investigation of the deformation zones and the corresponding cutting forces in orthogonal cutting [J]. Advanced Materials Research, 2011, 223: 152–161. DOI: 10.4028/www.scientific.net/AMR.223.152.
    [24] CHEN T, HUANG W, ZHANG W, et al. Experimental investigation on trajectory stability of high-speed water entry projectiles [J]. Ocean Engineering, 2019, 175: 16–24. DOI: 10.1016/j.oceaneng.2019.02.021.
  • 加载中
图(19) / 表(1)
计量
  • 文章访问数:  1442
  • HTML全文浏览量:  771
  • PDF下载量:  88
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-13
  • 修回日期:  2020-05-27
  • 刊出日期:  2020-12-05

目录

/

返回文章
返回