水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性_张轶凡.pdf
《水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性_张轶凡.pdf》由会员分享,可在线阅读,更多相关《水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性_张轶凡.pdf(15页珍藏版)》请在咨信网上搜索。
1、第 44 卷第 2 期2 0 2 3 年 2 月兵工学报ACTA AMAMENTAIIVol 44 No 2Feb2023DOI:10 12382/bgxb 2021 0598水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性张轶凡1,刘亮涛1,王金相1,李恒2,唐奎1(1 南京理工大学 瞬态物理国家重点实验室,江苏 南京 210094;2 海军研究院,北京 100161)摘要:为研究近场爆炸冲击波及气泡载荷对缩比后的鱼雷典型圆柱壳结构的毁伤特性,并探讨装药距离、装药方位等相关参数对圆柱壳结构变形特征及毁伤特性的影响规律,采用任意拉格朗日欧拉算法对近场爆炸冲击波及气泡载荷对圆柱壳结构的毁伤特性
2、进行了数值模拟分析。将壁面附近爆炸气泡演化过程的仿真结果和试验结果进行了对比,验证了数值模拟方法的有效性。利用该数值分析方法对多组不同装药距离和装药方位下的爆炸过程进行了研究,深入分析了冲击波、气泡脉动载荷及射流载荷等不同形式载荷对圆柱壳结构的毁伤作用机理。研究结果表明:冲击波对圆柱壳结构的毁伤效果受装药距离的影响较为明显,装药距离的增大会急剧削弱冲击波在圆柱壳上造成的破坏,而圆柱壳的方位改变对冲击波的毁伤作用影响较小;水射流载荷对圆柱壳结构产生的毁伤受方位因素的影响较为明显,当药包位于圆柱壳下方时圆柱壳迎爆面在气泡脉动及射流载荷联合作用下产生的塑性应变最大。关键词:水下爆炸;圆柱壳;应力应变
3、;爆距;气泡脉动中图分类号:U661.44;TB125;O359+.1文献标志码:A文章编号:1000-1093(2023)02-0345-15收稿日期:2021-09-03基金项目:国家自然青年科学基金项目(12002168);江苏省自然青年科学基金项目(BK20200469);总装备部装备预先研究基金项目(6142604190405)Damage Characteristics of Underwater Explosion Shock Waveand Bubble Load on Typical Cylindrical Shell StructureZHANG Yifan1,LIU Li
4、angtao1,WANG Jinxiang1,LI Heng2,TANG Kui1(1 National Key Laboratory of Transient Physics,Nanjing University of Science and Technology,Nanjing 210094,Jiangsu,China;2 Naval Academy of Armament,Beijing 100161,China)Abstract:To investigate the damage characteristics of the near-field explosion shock wav
5、e load andbubble load on the typical cylindrical shell structure of a scaled-down torpedo,the ALE method wasconducted and the effects of explosion distance and charge position on the deformation and damagecharacteristics of the cylindrical shell structure were explored The validity of the numerical
6、simulationmethod was verified by comparing it with the experimental results of the evolution process of the explosionbubble near the wall The interaction between the explosion bubble and the cylindrical shell was studiedunder different explosion distances and charge positions,and the damage mechanis
7、m of different forms ofload such the shock wave,bubble pulsation,and water jet on the cylindrical shell structure was analyzedin detail The results showed that the damage effect of the shock wave on the cylindrical shell structure is兵工学报第 44 卷strongly influenced by explosion distance;the increase of
8、 explosion distance sharply weakens the damagecaused by the shock wave on the cylindrical shell,while the charge position of the cylindrical shell has aweaker effect on the damage of the shock wave;the damage to the cylindrical shell is significantlyinfluenced by charge position;when the charge is a
9、rranged below the cylindrical shell,the shellgenerated the maximum plastic strain under the combined effects of bubble pulsation and water jetKeywords:underwater explosion;cylindrical shell;stress and strain;explosion distance;bubble pulsa-tion0引言随着国际形势的日益严峻和地区冲突的不断升级,海上作战成为一种重要形式,鱼雷是目前海军水下作战的重要武器装备
10、,同时水下爆炸也成为水中作战的主要毁伤方式之一。水中结构物在遭受水下爆炸作用时,冲击波及气泡载荷是对结构产生毁伤效果的主要因素。对于舰船等水下结构物,鱼雷在其近场爆炸造成的危害十分巨大,爆炸载荷会对结构物的整体强度造成巨大破坏。目前针对水下爆炸对舰船单位的毁伤研究,国内外学者已经做了大量的工作 1 5,得出了很多具有指导性的结论。但与此同时,鱼雷这一海军武器本身也作为一种水下结构物存在,在航行过程中同样会受到反鱼雷武器的影响和破坏。因此,研究水下爆炸冲击波和气泡载荷对鱼雷典型圆柱壳结构的毁伤特性,具有重要的军事意义。姚熊亮等6 对冲击波与圆柱壳结构相互作用过程中产生的压力场进行仿真和分析,给出
11、了在不同板厚与爆距下圆柱壳周围的压力场分布情况。刘东等7 基于内嵌的水下爆炸载荷计算方法和声固耦合方法,对不同壳体厚度及肋骨间距下的环肋圆柱壳在非接触爆炸中的变形响应进行了研究。仝志永等8 采用声固耦合算法分析了填沙薄壁圆柱壳结构水下爆炸作用下的弹塑性形变和损伤的破坏形式。Wu 等9 针对破片和爆炸载荷联合作用下金属圆柱壳结构的动态响应以及破坏效应进行了实验和数值模拟分析。沈晓乐等10 采用应变测量的方法 对 加 肋 圆 柱 壳 模 型 进 行 了 模 态 分 析。周游等11 利用流固耦合法分析了不同装药设置条件下圆柱壳的爆炸变形模式特征,分析了爆源条件对圆 柱 壳 动 力 响 应 及 变 形
12、 模 式 的 影 响 规 律。刘云龙等12 基于不可压缩势流理论,运用边界元方法研究了攻角对圆柱壳附近水下爆炸气泡动态特性的影响,并将双渐进法与有限元法相耦合,研究了圆柱壳结构受水下爆炸载荷作用的冲击损伤问题和主要损伤模式13。纪冲等14 对充液及内空圆柱壳在爆炸载荷下动力屈曲响应特性及不同工况下圆柱壳的变形破坏模式进行了研究。汪浩等15 研究了近距离非接触爆炸作用下的气泡脉动载荷以及气泡溃破高速射流载荷对内加筋圆柱壳结构的毁伤机理,以及结构参数、药包位置对结果的影响。综上所述,目前已有的研究工作对爆炸距离及方位相关参数影响规律的研究较少,仍有大量工作需要完成。为研究鱼雷中部位置在爆炸作用下的
13、响应情况,本文将鱼雷头部和尾部的典型特征去除,简化为圆柱壳结构进行分析。采用 LS-Dyna 软件对其在近场爆炸冲击作用下的毁伤特性进行研究,详细探讨爆距、药包方位等相关物理量对圆柱壳结构变形特征与毁伤特性的影响。1水下爆炸毁伤算法设置及材料模型缩比后的鱼雷典型圆柱壳结构长度为 1.1 m,半径 0.081 m,壳体厚度 0.003 m;在圆柱壳上均匀分布有 3 条加强肋,加强肋间距为 0.5 m,宽度为0.03 m,厚度 0.006 m。1.1水下爆炸圆柱壳结构模型算法设置在 LS-DYNA 软件中建立圆柱壳结构水下爆炸计算模型,采用任意拉格朗日欧拉(ALE)算法处理圆柱壳结构与水介质、空气
14、之间的相互作用,其中空气和水介质采用欧拉算法,圆柱壳结构采用拉格朗日算法。为保证计算的稳定性,对圆柱壳结构进行网格划分时尽可能满足对网格形状的优化;同时为保证流固耦合的准确性,对水域中部区域的网格进行加密,以保证此区域内的欧拉单元尺寸与圆柱壳中部的拉格朗日单元尺寸基本一致,加密区域网格大小为 0.01 m 0.01 m。各部分网格划分如图 1 所示。计算模型示意图如图 2 所示,数值计算采用国际单位制。其中欧拉域大小为4 m 4 m 6 m,空气层厚度为 0.3 m,对空气域上表面设置无反射边界条件,将空气域、水域的侧表面以及水域的底面设置为固壁边界条件,TNT 炸药被放置在计算域中心位643
15、第 2 期水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性图 1各部分网格划分Fig 1Meshing of various parts置,爆炸深度为 2.7 m。图 2数值模拟计算模型示意图Fig 2Schematic diagram of numerical simulation model最终得到各部分网格数量如表 1 所示。表 1各部分网格数量Table 1The number of grids in each part空气域水域圆柱壳总计71 8245 045 6364 5085 121 9681.2材料模型圆柱壳的本构模型选用 Johnson-Cook 模型,失效模型选用 Joh
16、nson-Cook 失效模型。Johnson-Cook模型是一种基于实验得到的、能够反映应变率强化效应和温升软化效应的材料模型16,其应力表达式为=(A+Bnp)(1+Cln*eq()1 T TrefTm T)ref(1)式中:为 von Mises 等效流变应力;A 代表材料屈服极限;B 为加工硬化模量;p为塑性应变;C 为应变速率常数;n 为硬化系数;*eq=p/0,p为等效塑性应变率,0为应变率参考值;T 为材料温度;Tref为温度参考值;Tm为材料融化温度。采用单元最大失效应变控制单元强度,等效塑性应变累计损伤度控制材料的破坏断裂,当单元的有效累计塑性应变达到失效应变时,发生失效侵蚀,
17、单元损伤系数表达式为D=pf(2)式中:D 为损伤系数,在 0 1 之间变化,初始状态下D=0,当 D=1 时材料失效;p为单个时间步的塑性应变增量;f为当前时间步对应的破坏应变,是当前应力、应变、应变率及温度的函数,f=D1+D2exp(D3*)(1+D4ln*)(1+D5T*)(3)D1 D5为材料参数;*=p/eff为应力状态参数,p 为正压力,eff为等效应力;T*为温度状态参数。模型中圆柱壳采用铝合金材料17 18,如表 2 所示。水介质和空气的材料本构模型采用常见的空物质模型,水介质采用 Gruneisen 状态方程19,其表达式为pw=wc2ww(1+1 0)2w22w1 (S1
18、1)w S22ww+1 S(33ww)+122+(0+w)Ew(4)式中:pw为介质中的压力;w为材料密度;cw为材料声速;w=1/Vw 1,Vw为相对体积;0为材料Gruneisen 系数;为 0的修正系数;S1、S2、S3为Gruneisen 状态方程材料常数;Ew为单位体积初始内能。水的材料参数20 如表 3 所示。空气采用 Liner Polynomial 多项式状态方程19,由式(5)定义:pa=C0+C1a+C22a+C33a+(C4a+C5a+C62a)Ea(5)式中:pa为介质中的压力;a=1/Va1,Va为相对体743兵工学报第 44 卷积;Ea为单位体积初始内能;Ci(i=
19、0 6)为状态方程参数。空气材料参数20 如表 4 所示。表 2铝合金材料参数Table 2Material parameters of aluminum alloyAl/(kg m3)G0/PaTm/Kcv/(J Kg1 K1)A/PaB/PanCmD1D2D3D4D52 7802.86 10100.347758752.4 1084.26 1080.340.01510.130.131.5 0.0110注:Al为铝合金的密度,G0为剪切模量,cv为比热容,m 为热软化强化参数。表 3水的材料参数Table 3Material parameters of waterw/(kg m3)cw/(m
20、s1)S1S2S30Ew/Pa1 0001 4801.9210.09600.3502.895 105表 4空气材料参数Table 4Material parameters of aira/(kg m3)C0C1C2C3C4C5C6Ea/Pa1.2900000.40.402.5 105爆轰产物的压力采用 JWL 状态方程19 来描述,其一般形式为pt=A(t1 1)Ve 1V+B(t1 2)Ve 2V+EtV(6)式中:pt为爆轰产物的压力;V 为爆轰产物和炸药的相对体积,为无量纲量;Et为单位体积内能;At、Bt、1、2、为待拟合参数。TNT 炸药的 JWL 参数20 如表 5 所示。表 5T
21、NT 材料参数Table 5Material parameters of TNTAt/PaBt/Pa12Et/(J kg1)3.712 10113.23 1094.150.950.33.8 1062算法有效性验证通过与文献 21 中的实验进行对比,对 ALE 数值模拟方法进行验证。采用4 g 的 PETN 炸药在2 m的正方体水箱中进行爆炸试验,炸药安置在厚度0.02 m 水平钢板下方 0.19 m 处。在数值模拟中采用 TNT 炸药对实验中的 PETN 炸药进行等效替代,PETN 和 TNT 的单位体积爆轰能量分别为9 109J/kg和7 109J/kg,通过计算得到等能量情况下 4 gPE
22、TN 炸药对应的 TNT 炸药质量为 5.2 g。不同时刻下气泡形态的实验和仿真结果的对比如图 3 所示。爆炸初期,气泡在内外压差的作用下呈球形膨胀,膨胀到一定大小后,逐渐受到上部平板的影响,气泡上部形态逐渐变得平坦,这一现象随着气泡的继续膨胀变得更加明显,通过图 3(d)可以直观地观察到这一现象;气泡膨胀到最大后在内外压差的作用下开始收缩,如图 3(e)所示;由于气泡上部钢板的存在,收缩阶段气泡上部的收缩速度明显小于气泡下部,表现为水平钢板对气泡的吸引作用,使得气泡在收缩阶段呈现为近似圆锥的形态,如图 3(h)所示;气泡在收缩的最后阶段形成了一个贯穿气泡中部的水射流,对水平钢板造成了射流抨击
23、。通过与实验中气泡形态演变规律的对照,发现数值模拟方法得出的气泡形态结果与实验结果能够吻合;同时,实验中爆炸气泡的第 1 次脉动周期为45.062 ms,通过数值模拟得到的气泡周期为51.5 ms,在气泡周期这一结果上存在一定误差,究其原因可能是实验中的炸药未完全起爆。根据 TNT 炸药无限水域爆炸气泡周期经验公式22:Tb=2.11W1/3(Z+10.1)5/6(7)式中:Tb为爆炸气泡周期;W 为装药量;Z 为爆炸水深,可以计算出无限水域中 5.2 g TNT 炸药在 0.8 m水深下的爆炸气泡周期为50.3 ms,而此数值模拟算法下同工况气泡周期为50.9 ms,数值模拟结果与理论计算值
24、的误差为 1.19%;由此表明此数值模拟方法的计算精度较高,能够满足研究的需求。3计算结果分析3.1不同爆距下冲击波和气泡载荷对圆柱壳的毁伤特性分析在保持 TNT 装药量、爆炸水深、边界条件不变的情况下,通过改变 TNT 药包与圆柱壳的相对距离,研究不同爆距下冲击波和气泡载荷对圆柱壳的毁伤特性。3.1.1工况设置设置 TNT 药量为 27.0 g,爆炸深度 2.7 m,药包到圆柱壳的距离从 0.05 m 到 0.30 m 进行设置,得到 6 组不同工况,如表 6 所示。843第 2 期水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性图 3不同时刻下气泡形态的实验与数值模拟结果对比Fig 3C
25、omparison of experimental and numerical results of bubble shape表 6计算工况Table 6Simulation conditions工况编号参数药量/g爆炸水深/m药包相对圆柱壳位置圆柱壳相对药包距离/m127.02.7下方0.05227.02.7下方0.10327.02.7下方0.15427.02.7下方0.20527.02.7下方0.25627.02.7下方0.303.1.2水下爆炸气泡与圆柱壳结构的响应过程分析取工况 1 进行过程分析,当药包位于圆柱壳下方 0.05 m 的位置时,气泡形态随时间的变化情况及圆柱壳在爆炸气泡作
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 水下 爆炸 冲击波 气泡 载荷 典型 圆柱 结构 毁伤 特性 张轶凡
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。