• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

横向爆炸载荷下薄壁圆管的动态响应

于博丽 冯根柱 李世强 刘志芳

引用本文:
Citation:

横向爆炸载荷下薄壁圆管的动态响应

    作者简介: 于博丽(1994- ),女,硕士研究生,964855831@qq.com;
    通讯作者: 刘志芳, liuzhifang@tyut.edu.cn
  • 中图分类号: O347.3

Dynamic response of thin-wall circular tubes under transverse blast loading

    Corresponding author: LIU Zhifang, liuzhifang@tyut.edu.cn
  • CLC number: O347.3

  • 摘要: 采用实验研究、理论分析和有限元模拟相结合的方法,研究了横向爆炸载荷作用下薄壁圆管的动态响应。利用弹道冲击摆锤系统,对圆管在爆炸载荷下的动力响应进行了实验研究,分析了薄壁圆管的变形模式;基于地基梁模型,建立了横向爆炸载荷作用下圆管跨中挠度的理论模型,并进行了无量纲化;通过有限元模拟,分析了圆管的几何参数对其变形模式和跨中挠度的影响,并与理论结果进行了对比。研究结果表明:随着TNT药量增加圆管的变形区域和跨中挠度增大;圆管的长径比、厚度及爆炸载荷参数对圆管的变形模式有较大影响;理论预测、有限元模拟结果与实验结果吻合较好。
  • 图 1  冲击摆锤系统

    Figure 1.  Ballistic pendulum system

    图 2  圆管残余变形模式

    Figure 2.  Residual deformation modes of circular tubes

    图 3  爆炸载荷下的圆管几何模型

    Figure 3.  Geometry of circular tube under blast loading

    图 4  爆炸载荷下圆管刚塑性地基梁

    Figure 4.  Rigid-plastic beam-on-foundation of the circular tube under impact loading

    图 5  In=0.35时无量纲挠度-无量纲量$\eta $

    Figure 5.  Non-dimensional deflection andnon-dimensional quantity $\eta $ at In =0.35

    图 6  In=0.35时无量纲挠度-无量纲量$\lambda $

    Figure 6.  Non-dimensional deflection andnon-dimensional quantity $\lambda $ at In=0.35

    图 7  实验与数值模拟变形模式对比

    Figure 7.  Comparison of numerical deformation modes with experimental result

    图 8  圆管跨中挠度的对比

    Figure 8.  Comparison of mid-span deflection of the circular tube

    图 9  不同In下无量纲挠度的数值模拟与理论预测的对比

    Figure 9.  Comparison of non-dimensional deflection between numerical results and theoretical results at different values of In

    表 1  试件几何参数和实验冲量与理论冲量对比

    Table 1.  Geometrical parameters and comparison of the experimental and theoretical impulses

    试件D/mh/mw/mmW/g$\bar H$/(m·kg−1/3)ΔpΦ/(kg·cm−2)ΔpΦr/(kg·cm−2)p0/MPaIE/(N·s)IT/(N·s)[(IT-IE)/IE]/%
    890.9 5.5200.5522.7148.214.8 9.8 9.2−6.12
    890.829.1350.4634.0236.023.614.914.7−1.34
    760.826.8350.4634.0236.023.613.612.6−7.35
    760.729.5350.4634.0236.023.613.412.6−5.97
    下载: 导出CSV

    表 2  横向爆炸载荷下圆管跨中挠度

    Table 2.  Mid-span deflection circular tube under transverse blast loading

    D/mmh/mmp0/MPaw/mm
    TheoryFEAExperiment
    760.723.631.6834.2129.5
    760.823.624.1928.8426.8
    760.923.619.0622.86
    890.714.810.7214.19
    890.814.8 8.19 9.80
    890.914.8 6.46 7.355.5
    890.723.627.2538.34
    890.823.620.8329.2429.1
    890.923.616.4323.32
    下载: 导出CSV
  • [1] PAYTON R G. Transient interaction of an acoustic wave with a circular cylindrical elastic shell [J]. Journal of the Acoustical Society of America, 1959, 32(6): 722–729. DOI: 10.1121/1.1908197.
    [2] MENKES S B, OPAT H J. Broken beams [J]. Experimental Mechanics, 1973, 13(11): 480–486. DOI: 10.1007/BF02322734.
    [3] 杨桂通. 塑性动力学 [M]. 北京: 清华大学出版社, 1984.
    [4] WEGENER R B, MARTIN J B. Predictions of permanent deformation of impulsively loaded simply supported square tube steel beams [J]. International Journal of Mechanical Sciences, 1985, 27(1): 55–69. DOI: 10.1016/0020-7403(85)90066-9.
    [5] WIERZBICKI T. Damage assessment of cylinders due to impact and explosive loading [J]. International Journal of Impact Engineering, 1993, 13(2): 215–241. DOI: 10.1016/0734-743X(93)90094-N.
    [6] HOOFATT M, WIERZBICKI T. Damage of plastic cylinders under localized pressure loading [J]. International Journal of Mechanical Sciences, 1991, 33(12): 999–1016. DOI: 10.1016/0020-7403(91)90055-8.
    [7] YU T X, STRONGE W J. Large deflection of a rigid-plastic beam-on-foundation from impact [J]. International Journal of Impact Engineering, 1990, 9(1): 115–126. DOI: 10.1016/0734-743X(90)90025-Q.
    [8] YUEN S C K, NURICK G N, BRINCKMANN H B, et al. Response of cylindrical shells to lateral blast load [J]. International Journal of Protective Structures, 2013, 4(3): 209–230. DOI: 10.1260/2041-4196.4.3.209.
    [9] JAMA H H, NURICK G N, BAMBACH M R, et al. Steel square hollow sections subjected to transverse blast loads [J]. Thin-Walled Structures, 2012, 53: 109–122. DOI: 10.1016/j.tws.2012.01.007.
    [10] BAMBACH M R. Behavior and design of aluminum hollow sections subjected to transverse blast loads [J]. Steel Construction, 2008, 46(12): 1370–1381. DOI: 10.1016/j.tws.2008.03.010.
    [11] KARAGIOZOVA D, YU T X, LU G. Transverse blast loading of hollow beams with square cross-sections [J]. Thin-Walled Structures, 2013, 62(1): 169–178. DOI: 10.1016/j.tws.2012.09.004.
    [12] KARAGIOZOVA D, YU T X, LU G, et al. Response of a circular metallic hollow beam to an impulsive loading [J]. Thin-Walled Structures, 2014, 80: 80–90. DOI: 10.1016/j.tws.2014.02.021.
    [13] 李世强, 李鑫, 吴桂英, 等. 梯度蜂窝夹芯板在爆炸荷载作用下的动力响应 [J]. 爆炸与冲击, 2016, 36(3): 333–339. DOI: 10.11883/1001-1455(2016)03-0333-07.
    LI Shiqiang, LI Xin, WU Guiying, et al. Dynamic response of functionally graded honeycomb sandwich plates under blast loading [J]. Explosion and Shock Waves, 2016, 36(3): 333–339. DOI: 10.11883/1001-1455(2016)03-0333-07.
    [14] Cole R H, Weller R. Underwater explosions [J]. Physics Today, 1948, 1(6): 35–35. DOI: 10.1063/1.3066176.
    [15] KARAGIOZOVA D, NURICK G N, LANGDON G S. Behavior of sandwich panels subject to intense air blasts. Part 2: Numerical simulation [J]. Composite Structures, 2009, 91(4): 442–450. DOI: 10.1016/j.compstruct.2009.04.010.
    [16] HENRYCH J. The dynamics of explosion and its use [M]. Elsevier Scientific Publishing Company, 1979: 562. DOI: 10.1115/1.3153619.
    [17] MARTIN J B, SYMONDS P S. Mode approximations for impulsively loaded rigid plastic structures [J]. Journal of the Engineering Mechanics Division, 1966, 92(5): 43–66. DOI: 10.21236/ad0621580.
    [18] WALTERS R M, JONES N. An approximate theoretical study of the dynamic plastic behavior of shells [J]. International Journal of Non-Linear Mechanics, 1972, 7(3): 255–273. DOI: 10.1016/0020-7462(72)90049-2.
    [19] 余同希, 华云龙. 结构塑性动力学引论 [M]. 合肥: 中国科技大学出版社, 1994: 88−89.
  • [1] 俞鑫炉董新龙潘顺吉 . 不同爆炸载荷下TA2钛合金圆管膨胀破坏过程. 爆炸与冲击, 2018, 38(1): 148-154. doi: 10.11883/bzycj-2017-0014
    [2] 张世文刘仓理于锦泉 . 微缺陷对圆管膨胀断裂的影响. 爆炸与冲击, 2008, 28(4): 316-323. doi: 10.11883/1001-1455(2008)04-0316-08
    [3] 李猛深李杰李宏施存程张宁 . 爆炸荷载下钢筋混凝土梁的变形和破坏. 爆炸与冲击, 2015, 35(2): 177-183. doi: 10.11883/1001-1455(2015)02-0177-07
    [4] 范志庚陈常青万强 . 泡沫铝率相关性能的有限元模拟. 爆炸与冲击, 2014, 34(6): 742-747. doi: 10.11883/1001-1455(2014)06-0742-06
    [5] 谢凡张涛陈继恩刘土光 . 应力三轴度的有限元计算修正. 爆炸与冲击, 2012, 32(1): 8-14. doi: 10.11883/1001-1455(2012)01-0008-07
    [6] 张晓天贾光辉黄海 . 基于节点分离Lagrange有限元方法的超高速碰撞碎片云数值模拟. 爆炸与冲击, 2010, 30(5): 499-504. doi: 10.11883/1001-1455(2010)05-0499-06
    [7] 郑波王安稳 . 直杆碰撞刚性壁弹塑性动力后屈曲有限元分析. 爆炸与冲击, 2007, 27(2): 126-130. doi: 10.11883/1001-1455(2007)02-0126-05
    [8] 王吉王肖钧卞梁 . 光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用. 爆炸与冲击, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
    [9] 颜峰姜福兴 . 爆炸冲击载荷作用下岩石的损伤实验. 爆炸与冲击, 2009, 29(3): 275-280. doi: 10.11883/1001-1455(2009)03-0275-06
    [10] 冯慧平刘洪兵左兴惠浪浪 . 地下坑道对其临界震塌爆距处钻地武器爆炸荷载的动力响应. 爆炸与冲击, 2014, 34(5): 539-546. doi: 10.11883/1001-1455(2014)05-0539-08
    [11] 张旭红王志华赵隆茂 . 爆炸载荷作用下铝蜂窝夹芯板的动力响应. 爆炸与冲击, 2009, 29(4): 356-360. doi: 10.11883/1001-1455(2009)04-0356-05
    [12] 侯海量朱锡古美邦 . 爆炸载荷作用下加筋板的失效模式分析及结构优化设计. 爆炸与冲击, 2007, 27(1): 26-33. doi: 10.11883/1001-1455(2007)01-0026-08
    [13] 侯海量朱锡梅志远 . 舱内爆炸载荷及舱室板架结构的失效模式分析. 爆炸与冲击, 2007, 27(2): 151-158. doi: 10.11883/1001-1455(2007)02-0151-08
    [14] 李忠盛吴护林陈韵如易同斌李立潘复生 . 内爆炸载荷作用下7A55铝合金的动态性能及断裂行为. 爆炸与冲击, 2012, 32(2): 190-195. doi: 10.11883/1001-1455(2012)02-0190-06
    [15] 陈鹏宇侯海量吴林杰朱锡 . 水下舷侧多层防护隔舱接触爆炸毁伤载荷特性分析. 爆炸与冲击, 2017, 37(2): 283-290. doi: 10.11883/1001-1455(2017)02-0283-08
    [16] 苏兴亚敬霖赵隆茂 . 爆炸载荷下分层梯度泡沫铝夹芯板的失效模式与抗冲击性能. 爆炸与冲击, 2019, 39(6): 063103-1-063103-8. doi: 10.11883/bzycj-2018-0198
    [17] 刘文祥张德志钟方平李焰 . 球形和等长径圆柱装药在爆炸近区内 载荷差异的实验研究. 爆炸与冲击, 2013, 33(3): 330-336. doi: 10.11883/1001-1455(2013)03-0330-07
    [18] 徐海斌钟方平杨军张德志秦学军刘峻岭史国凯梁志刚沈兆武 . 水及盛水容器对近距离爆炸载荷影响的实验研究. 爆炸与冲击, 2016, 36(4): 525-531. doi: 10.11883/1001-1455(2016)04-0525-07
    [19] 史博张晓颖李阔李志强 . 爆炸载荷下中空夹层玻璃的动力响应影响因素. 爆炸与冲击, 2018, 38(1): 119-123. doi: 10.11883/bzycj-2017-0018
    [20] 李强李鹏辉赵君官茹站勇 . 冲击载荷下自增强身管残余应力变化规律. 爆炸与冲击, 2011, 31(6): 635-640. doi: 10.11883/1001-1455(2011)06-0635-06
  • 加载中
图(9)表(2)
计量
  • 文章访问数:  56
  • HTML全文浏览量:  170
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-08
  • 录用日期:  2018-12-27
  • 网络出版日期:  2019-09-25
  • 刊出日期:  2019-10-01

横向爆炸载荷下薄壁圆管的动态响应

    作者简介:于博丽(1994- ),女,硕士研究生,964855831@qq.com
    通讯作者: 刘志芳, liuzhifang@tyut.edu.cn
  • 1. 太原理工大学应用力学研究所,山西 太原 030024
  • 2. 太原理工大学力学学院材料强度与结构冲击山西省重点实验室,山西 太原 030024

摘要: 采用实验研究、理论分析和有限元模拟相结合的方法,研究了横向爆炸载荷作用下薄壁圆管的动态响应。利用弹道冲击摆锤系统,对圆管在爆炸载荷下的动力响应进行了实验研究,分析了薄壁圆管的变形模式;基于地基梁模型,建立了横向爆炸载荷作用下圆管跨中挠度的理论模型,并进行了无量纲化;通过有限元模拟,分析了圆管的几何参数对其变形模式和跨中挠度的影响,并与理论结果进行了对比。研究结果表明:随着TNT药量增加圆管的变形区域和跨中挠度增大;圆管的长径比、厚度及爆炸载荷参数对圆管的变形模式有较大影响;理论预测、有限元模拟结果与实验结果吻合较好。

English Abstract

  • 薄壁构件广泛应用于建筑、航空航天、采矿、军事等领域。由于现代技术工程的发展安全要求,此类结构在冲击载荷作用下的动态响应也引起了关注。发生爆炸时工程结构的安全可靠性已逐渐成为人们关注的重要问题,但由于结构响应复杂、随机变量多,使得爆炸载荷作用下工程结构的理论分析有一定难度。工程结构中圆管结构受各种冲击载荷的作用,人们对它的变形机理和能量吸收进行了研究。

    Payton[1]基于爆炸冲击波为平面波的假定,研究了无限长圆柱壳在瞬间声波冲击下的弹性动力响应问题。Menkes等[2]研究了横向冲击载荷作用下的梁,结果表明,随着冲量的增大,结构的破坏模式发生变化,且变化规律一般是从塑性失稳过渡为剪切破坏。结构的弹塑性动力响应问题已经得到比较广泛而深入的研究[3]。Wegener等[4]研究了冲击荷载作用下方形截面空心简支钢质梁的力学行为,采用数值分析导出了半经验解析解来确定梁的变形模式。Wierzbicki[5]提出了基于地基梁模型预测冲击和爆炸载荷作用下圆柱体破坏的分析方法。Hoofatt等[6]提出一种将圆柱壳的二维边值问题转化为等效的非线性塑性地基上塑性弦的一维问题的通用方法。对于矩形压力载荷和高斯压力载荷两种压力载荷,得到了壳体轴向变形和速度场的封闭解及各载荷情况下的应变分布和最大应变,并简要预测了壳体断裂的可能性。Yu等[7]利用膜力因子法,分析了质量块作用下刚塑性地基梁的动态变形。Yuen等[8]通过一系列实验和数值模拟,研究了外部局部爆炸载荷作用下圆柱壳的响应,解析解低估了永久中点挠度,与数值模拟结果具有良好的相关性。Jama等[9]通过实验研究了横向爆炸载荷下钢质空心方形截面梁的动态响应,并进行了理论分析;Bambach[10]在横向爆炸载荷作用下对铝质梁进行实验研究,提出了一种半经验的解析解,给出了一种基于最终永久变形或失效载荷的通用设计方法。Karagiozova等[11]研究了横向爆炸载荷作用下空心正方形截面梁的变形与能量吸收,理论分析结果与已有文献中的实验结果取得了较好的一致性;Karagiozova等[12]采用刚塑性分析法,研究了冲击载荷对空心圆形梁的动态响应的两相变形模型,将爆炸冲击近似为脉冲载荷,假设在第一变形阶段同时发生局部和整体变形,第二阶段只发生整体弯曲变形。

    综上所述,爆炸载荷下薄壁管动态响应的研究主要集中于理论分析与数值模拟,理论分析主要是通过刚塑性梁模型或模态分析法来预测横向爆炸载荷下梁中点的挠度。本文中,采用实验研究、理论分析和有限元模拟相结合的方法,研究横向爆炸载荷下圆管的动态响应。利用弹道冲击摆锤系统,对圆管在横向爆炸载荷下的动态响应进行实验研究;基于理想塑性地基上刚塑性梁模型,采用模态分析法研究圆管变形模式,推导横向爆炸载荷作用下圆管跨中挠度的解析解;采用有限元模拟,分析脉冲加载下圆管的变形模式,研究不同的脉冲与圆管的直径、厚度对其跨中挠度的影响。

    • 利用弹道冲击摆锤系统[13]对横向爆炸载荷下的圆管进行了实验研究,如图1所示。试件所受的冲量通过测定摆锤系统摆动的位移来确定。圆管两端固定在摆锤的夹具中,距离圆管跨中H=150 mm的位置放置TNT炸药。爆炸后测量圆管跨中挠度,圆管变形后的侧视图、径向截面图和轴向截面图如图2所示。试件几何参数、TNT药量$W$表1。试件材料参数分别为:密度$\rho $=7 830 kg/m3,弹性模量E=193 GPa,泊松比$\mu $=0.25,屈服应力${\sigma _s}$=205 MPa,在${\varepsilon _{{\rm{uts}}}}$=0.4时,${\sigma _{{\rm{uts}}}}$=520 MPa。圆管试件壁厚分别为0.7、0.8、0.9 mm,外径为89、76 mm,长度为280 mm。

      图  1  冲击摆锤系统

      Figure 1.  Ballistic pendulum system

      图  2  圆管残余变形模式

      Figure 2.  Residual deformation modes of circular tubes

      试件D/mh/mw/mmW/g$\bar H$/(m·kg−1/3)ΔpΦ/(kg·cm−2)ΔpΦr/(kg·cm−2)p0/MPaIE/(N·s)IT/(N·s)[(IT-IE)/IE]/%
      890.9 5.5200.5522.7148.214.8 9.8 9.2−6.12
      890.829.1350.4634.0236.023.614.914.7−1.34
      760.826.8350.4634.0236.023.613.612.6−7.35
      760.729.5350.4634.0236.023.613.412.6−5.97

      表 1  试件几何参数和实验冲量与理论冲量对比

      Table 1.  Geometrical parameters and comparison of the experimental and theoretical impulses

      由激光位移传感器Micro-Epsilon LD1607-200测得圆管跨中位移,可计算系统冲量:

      ${I_{\rm{E}}} = M{x_1}\frac{{2{\text{π}}}}{T}{{\rm{e}}^{\textstyle\frac{1}{4}\ln \frac{{{x_1}}}{{{x_2}}}}}$

      式中:摆锤的总质量M=151.3 kg,周期T=3.14 s,${x_1}$${x_2}$分别为T1=T/4和T2=3T/4时的位移。摆线的长度为2.6 m。

      为了消除变量的量纲效应,定义一个无量纲冲量:

      ${I_{\rm{n}}} = \frac{{2{I_{\rm{E}}}}}{{{\text{π}}DLh\sqrt {\rho {\sigma _{\rm{s}}}} }}$

      式中:L、D、h分别为圆管的长度、外径与厚度。

    • 实验结果分为两组:(1)定量结果,包括圆管上的冲量IE、跨中的永久挠度w;(2)试件的变形破坏模式。圆管在不同爆炸加载条件下的实验结果见表1

      圆管的失效有3种不同的失效形式。第1类是局部塑性变形模式(见图2试件Ⅰ变形);第2类是菱形形状的大的塑性变形(见图2试件Ⅱ、Ⅲ、Ⅳ变形);第3类是圆管大的非弹性变形[8]。当TNT放置距离一定时,前面区域的破坏区和跨中永久挠度随炸药量的增加而增加。从图2可以知,在圆管上半部分顶部有一个扁平的区域并存在着面内屈曲,而圆管下半部分仍然保持原始形状。对称面上半部的挠度沿轴线逐渐减小,圆管下半部分没有明显的变形。

    • 对于爆炸载荷作用下的圆管结构,假设由爆炸引起的脉冲的近似方程为[14]

      $p\left( {x,\theta ,t} \right) = {p_0}{{\rm{e}}^{ - \textstyle{ \frac{t}{\tau }}}}f\left( x \right)g\left( \theta \right)$

      式中:${p_0}$为爆炸载荷峰值,衰减常数τ=50 µs[15]$f\left( x \right)$$g\left( x \right)$分别为轴向和圆周方向的初始压力分布函数。取$f\left( x \right) = {\cos ^2}\left( {{\text{π}}x/L} \right)$[5]$g\left( \theta \right) = \sin \theta $[12],则爆炸冲击载荷为:

      $p\left( {x,\theta ,t} \right) = {p_0}{{\rm{e}}^{ - \textstyle{\frac{t}{\tau }}}}{\cos ^2}({\text{π}}x/L)\sin \theta $

      式中:$0 {\text{<} }x {\text{<}} {L/2}$$0 {\text{<}}\theta {\text{<}} {{\text{π}} /2}$。圆管受爆炸载荷作用如图3(a)所示,周向压力分布如图3(b)所示。

      图  3  爆炸载荷下的圆管几何模型

      Figure 3.  Geometry of circular tube under blast loading

      圆管上施加的总脉冲表示为:

      ${I_{\rm{T}}} = 4\int_0^{L/2} {\int_0^{{\text{π}}/2} {\int_0^\infty {p(x,\theta ,t)\frac{D}{2}} } {\rm{d}}x{\rm{d}}\theta {\rm{d}}t} = \frac{1}{2}{p_0}DL\tau $

      根据大量爆炸实验[16],爆炸产生的超压$\Delta {p_\Phi }$(单位为kg/cm2)为比例距离$\bar H$(单位为m/kg1/3)的经验函数:

      $\begin{aligned} \Delta {p_\Phi } = \frac{{6.193\;8}}{{\bar H}} - \frac{{0.326\;2}}{{{{\bar H}^2}}} + \frac{{2.132\;4}}{{{{\bar H}^3}}}\quad\quad\quad0.3{\text{≤}}\bar H{\text{≤}}1.0,\quad \bar H = H/\sqrt[3]{W},\quad \Delta {p_\Phi } = {p_\Phi } - {p_{\rm{a}}} \end{aligned}$

      式中:H为炸药距离,${p_\Phi }$为冲击波阵面压力,${p_{\rm{a}}}$为大气压。反射爆炸超压为:

      $\Delta {p_{\Phi {\rm{r}}}} = 2\Delta {p_\Phi } + \frac{{6\Delta p_\Phi ^2}}{{\Delta {p_\Phi } + 7.2}} = \frac{{8\Delta p_\Phi ^2 + 14\Delta {p_\Phi }}}{{\Delta {p_\Phi } + 7.2}}\quad\quad\quad\Delta {p_\Phi } {\text{≤}} 40$

      实验结果与理论分析结果的对比见表1,其中冲量${I_{\rm{E}}}$为实验结果,压力峰值${p_0} = \Delta {p_{\Phi {\rm{r}}}}= \Delta {p_{\Phi {\rm{r}}}} \times $0.1 MPa,超压$\Delta {p_\Phi }$和比例距离$\bar H$由式(6)计算得到,总冲量${I_{\rm{T}}}$$\Delta {p_{\Phi {\rm{r}}}}$分别由式(5)、(7)给出。由表1可知,超压$\Delta {p_\Phi }$的计算结果均小于40 kg/cm2,式(7)适用。当$\Delta {p_\Phi } \text{>}40\;{\rm{kg}}/{{\rm{cm}}^2}$时,不在本文的研究范围内,式(7)也不再适用。

    • 圆管受爆炸载荷作用时,在载荷作用区域及其附近周围的壳面将产生大的塑性变形和转动,精确的分析相当复杂。通常利用薄壁圆管受局部冲击作用下的响应主要是轴向拉伸及环向变形的特点,将它类比为地基梁模型来计算。对于有限长薄壁圆管在横向爆炸载荷下的动态响应,采用理想塑性地基上刚塑性梁的模型进行分析,如图4所示。

      图  4  爆炸载荷下圆管刚塑性地基梁

      Figure 4.  Rigid-plastic beam-on-foundation of the circular tube under impact loading

      根据方程(4)给出的爆炸载荷形状,假设初始速度场${v^{(0)}}(x,\theta )$为:

      ${v^{(0)}}\left( {x,\theta } \right) = {v_0}{f^*}\left( x \right){g^*}\left( \theta \right) = {v_0}{\cos ^2}\left(\frac{{{\text{π}}x}}{L}\right)\sin \theta \quad\quad\quad 0 \text{≤} x \text{≤} \frac{L}{2},\quad\;\;0 \text{≤} \theta \text{≤} \frac{{\text{π}}}{2}$

      式中:${v_0}$为跨中的初始速度,${f^*}\left( x \right)$${g^*}\left( x \right)$分别为轴向和圆周方向的初始速度场函数,并且根据爆炸载荷方程取为${f^*}\left( x \right) = {\cos ^2}({\text{π}}x/L)$${g^*}\left( \theta \right) = \sin \theta $。作用在圆管上的横向动量$\varOmega $为:

      $\varOmega = 4\int_0^{L/2} {\int_0^{{\text{π}}/2} {{m_{\rm{A}}}\frac{D}{2}{v_0}{{\cos }^2}\left( {\frac{{{\text{π}}x}}{L}} \right)\sin \theta } {\rm{d}}x{\rm{d}}\theta } $

      式中:${m_{\rm{A}}} = \rho h$为圆管上半部分的面积质量。根据${I_{\rm{T}}} = \varOmega $,可给出跨中的等效初始速度:

      ${v_0} = \frac{{{p_0}\tau }}{{{m_{\rm{A}}}}}$

    • 假设地基梁符合如下屈服准则:

      $\left| {\frac{M}{{{M_{\rm{p}}}}}} \right| + {\left( {\frac{N}{{{N_{\rm{p}}}}}} \right)^2} = 1$

      式中:${M_{\rm{p}}} = {\sigma _{\rm{s}}}{D^2}h$为圆管截面的塑性极限弯矩,${N_{\rm{p}}} = {\sigma _{\rm{s}}}{\text{π}} hD$为圆管的塑性极限轴力。

      由于圆管结构和承受载荷的对称性,取圆管的1/2进行分析。基于爆炸载荷作用下圆管实验研究观察到的变形模式(见图2),假设圆管变形的形函数$\varphi _i^*(x) = \cos ({\text{π}}x/L)$,则圆管跨中模态解[17]的初始速度$\dot w_*^{(0)}$为:

      $\dot w_*^{(0)} = \frac{{\displaystyle\int_V {\rho {v_0}{f^*}\varphi _i^*{\rm{d}}V} }}{{\displaystyle\int_V {\rho \varphi _i^*\varphi _i^*{\rm{d}}V} }} = \frac{{2\displaystyle\int_0^{L/2} {{m_{\rm{L}}}{v_0}{{\cos }^2}\left( {\dfrac{{{\text{π}}x}}{L}} \right)\cos \left( {\dfrac{{{\text{π}}x}}{L}} \right){\rm{d}}x} }}{{2\displaystyle\int_0^{L/2} {{m_{\rm{L}}}{{\cos }^2}\left( {\dfrac{{{\text{π}}x}}{L}} \right){\rm{d}}x} }} = \frac{{8{v_0}}}{{3{\text{π}}}}$

      式中:${m_{\rm{L}}} = {\text{π}}Dh\rho /2$为圆管上半部分的长度质量。相应的减速度${a_*}$为:

      ${a_*} = \frac{{\displaystyle\int_V {Q_j^*{{\dot q}_j}{\rm{d}}V} }}{{\displaystyle\int_V {{m_{\rm{L}}}{{\dot \omega }_*}\varphi _i^*\varphi _i^*{\rm{d}}V} }}$

      式中:广义应力$Q_j^*$为塑性极限弯矩,$Q_1^* = M$,轴力$Q_2^* = N$,地基反力$Q_3^* = q$$\dot q_j^*$为广义应变率,$\dot q_1^*$对应塑性铰的转动角速度为$\dot \theta = {\dot w_*}(t)/\left( {L/2} \right)$。结合${\dot q_j} = {\dot w_*}(t)k_j^*(x)$${\dot q_3} = {\dot w_*}(t)\cos \left( {{\text{π}}x/L} \right)$,可得:

      $k_1^*(x) = {2/L}$

      $k_3^*(x) = \cos \left( {\frac{{{\text{π}}x}}{L}} \right)$

      根据塑性理论中关联流动法则,广义应力与应义应变率具有相互正交的关系,由Drucker公设,得出广义应变率${\dot q_1}$和轴向应变率${\dot q_2}$之间的比值:

      $\frac{{{{\dot q}_1}}}{{{{\dot q}_2}}} = \frac{{N_{\rm{p}}^2}}{{2N{M_{\rm{p}}}}}$

      因此,$k_2^*(x)$可以表示为:

      $k_2^*(x) = \frac{{4N{M_{\rm{p}}}}}{{LN_{\rm{p}}^2}}$

      圆管在爆炸载荷作用下发生了显著的塑性变形,弹性变形的影响可忽略,因此在分析过程中将材料简化为理想刚塑性材料。采用理想塑性地基上刚塑性梁的屈服条件[18],假设$M = 0$$N = {N_{\rm{p}}}$$q = {{16{M_0}}/D}$,其中${M_0} = {{{\sigma _{\rm{s}}}{h^2}}/4}$为凹形环模型单位长度的塑性极限弯矩。

      将式(14)、(15)和(17)和屈服准则代入式(13),可得:

      ${a_*} = \frac{{32{M_{\rm{p}}}/L + 64{M_0}L/\left( {{\rm{\pi }}D} \right)}}{{{m_{\rm{L}}}L}}$

      对于理想刚塑性结构,弹性效应、应变强化和应变率效应可以忽略,广义应力$Q_j^*$和广义应变率$k_j^*$都与时间无关。结构的总响应时间为$t_{\rm{f}}^* = \dot w_*^{(0)}/{a_*}$,因此跨中的最大挠度为:

      $w_*^{\rm{f}} = \frac{{\rm{1}}}{{\rm{2}}}\dot w_*^{(0)}t_{\rm{f}}^* = \frac{{I_{\rm{T}}^2}}{{\dfrac{9}{2}\rho {\sigma _{\rm{s}}}{\text{π}} {D^3}{h^2} + \dfrac{9}{4}\rho {\sigma _{\rm{s}}}{L^2}{h^3}}}$

      式(19)中,用实验冲量IE代替假设冲量IT,可得到实验冲量产生的近似挠度$\tilde w_*^{\rm{f}}$,引入无量纲量$\lambda = {L/D}$$\eta = {h/D}$,将In代入可得跨中无量纲挠度为:

      $\bar w_*^{\rm{f}} = \frac{{\tilde w_*^{\rm{f}}}}{h} = \frac{{I_{\rm{E}}^2}}{{\dfrac{9}{2}\rho {\sigma _{\rm{s}}}{\text{π}}{D^3}{h^3} + \dfrac{9}{4}\rho {\sigma _{\rm{s}}}{L^2}{h^4}}} = \frac{{{{\text{π}}^2}{\lambda ^2}}}{{18{\text{π}}\eta + 9{\lambda ^2}{\eta ^2}}}I_{\rm{n}}^2$

      式(20)表明,跨中挠度取决于$\lambda $$\eta $In图5给出了当In=0.35、$\lambda $为3.15和3.68时的对比图,跨中挠度随$\eta $的增大而减小。图6给出了In=0.35、$\eta $为0.008和0.011时的对比图,跨中挠度随$\lambda $的增大而增大。

      图  5  In=0.35时无量纲挠度-无量纲量$\eta $

      Figure 5.  Non-dimensional deflection andnon-dimensional quantity $\eta $ at In =0.35

      图  6  In=0.35时无量纲挠度-无量纲量$\lambda $

      Figure 6.  Non-dimensional deflection andnon-dimensional quantity $\lambda $ at In=0.35

    • 采用有限元软件ABAQUS对不同壁厚的圆管在横向爆炸载荷作用下进行了数值分析,研究了圆管的变形模式和动态响应。圆管两端固定,上半部分承受爆炸载荷,载荷以式(4)的形式施加在圆管上。圆管的材料参数选取与实验研究相同的参数,材料模型选取Cowper-Symonds本构关系,其中应变率敏感率$D = 100\;{{\rm s}^{ - 1}}$$q = 10$[19]。圆管采用S4R壳单元,并对有限元模型进行了网格敏感性验证,最终网格尺寸选为2 mm。

    • 图7给出了试件Ⅲ的数值模拟最终变形模式与实验最终变形模式。由图7可知,结构整体变形模式与有限元结果吻合较好,但由于实验中炸药产生的载荷非常复杂,同时伴随有高温高压环境,有限元结果中并未观察到与实验类似的面内屈曲变形。图8给出了4组不同壁厚的试件在不同脉冲荷载下跨中挠度的理论预测、实验测试和数值模拟结果。由图8可知,理论预测、数值模拟结果与实验测试吻合较好,表明所建立理论分析模型是合理的。在此基础上,研究了不同脉冲载荷下圆管的几何参数对圆管跨中挠度的影响,详细的参数与模拟结果见表2。由表2可知,随着圆管厚度及脉冲荷载的增大跨中挠度增大,随着圆管外径的增大跨中挠度减小,因此脉冲载荷和几何参数对圆管变形和跨中挠度有较大影响。由表2可以看出,理论结果与数值模拟结果、实验结果之间存在一定的误差:原因可能是由于实验中材料发生较大的塑性变形时产生了一定的应变硬化;并且理论与数值模拟中并未考虑装药量、形状以及装药距离与实际的超压峰值及分布规律之间的非线性关系,这也是造成误差的原因。

      图  7  实验与数值模拟变形模式对比

      Figure 7.  Comparison of numerical deformation modes with experimental result

      图  8  圆管跨中挠度的对比

      Figure 8.  Comparison of mid-span deflection of the circular tube

      D/mmh/mmp0/MPaw/mm
      TheoryFEAExperiment
      760.723.631.6834.2129.5
      760.823.624.1928.8426.8
      760.923.619.0622.86
      890.714.810.7214.19
      890.814.8 8.19 9.80
      890.914.8 6.46 7.355.5
      890.723.627.2538.34
      890.823.620.8329.2429.1
      890.923.616.4323.32

      表 2  横向爆炸载荷下圆管跨中挠度

      Table 2.  Mid-span deflection circular tube under transverse blast loading

      根据无量纲分析理论,当$\lambda $=3,$\eta $=0.01时,取D=100 mm,L=300 mm,h=1 mm和D=50 mm,L=150 mm,h=0.5 mm两种尺寸的试件进行数值模拟。通过无量纲冲量In计算p0,由式(2)可得:

      ${I_{\rm{E}}} = {I_{\rm n}}{\text{π}}DLh\sqrt {\rho {\sigma _{\rm{s}}}} /2$

      IEIT近似相等,代入式(5)可得:

      ${p_0} = \frac{{{I_{\rm n}}{\text{π}}h\sqrt {\rho {\sigma _{\rm{s}}}} }}{\tau }$

      通过数值模拟计算,对结果进行无量纲化(${W_{{\rm{FEA}}}} = {w_{{\rm{FEA}}}}/h$),结果如图9所示。

      图  9  不同In下无量纲挠度的数值模拟与理论预测的对比

      Figure 9.  Comparison of non-dimensional deflection between numerical results and theoretical results at different values of In

      图9可知,当冲量较小时数值模拟结果与理论结果吻合较好。当冲量较大时,结构会产生较大的挠度和塑性变形,挠度越大,膜力做功所占比例越大。但是由于在理论计算中,为了工程计算的简易性,采用了圆管截面梁的全塑性极限弯矩(${M_{\rm{p}}} = {\sigma _{\rm{s}}}{D^2}h$),较高估计了减速度${a_*}$,这导致了在较大冲量下理论计算结果小于数值模拟计算结果。

    • 采用实验研究、理论分析和有限元模拟相结合的方法,分析了在横向爆炸载荷作用下圆管的动态响应,对比分析了圆管的几何参数与圆管跨中挠度之间的关系,并通过实验研究、数值模拟验证了理论分析模型的合理性。结果表明:随着TNT药量增加,圆管的变形区域和跨中挠度增大;基于完全塑性地基上刚塑性梁模型,利用模态分析法给出了横向爆炸载荷下圆管中点挠度的解析解,圆管的长径比、厚度及爆炸载荷参数是影响圆管的变形模式的主要因素,通过无量纲分析得到了无量纲长度$\lambda $、壁厚$\eta $以及无量纲冲量In对跨中挠度的影响规律,随着无量纲长度$\lambda $及无量纲冲量In的增大,跨中挠度增大,随着无量纲量壁厚$\eta $增大,跨中挠度减小;理论预测、数值模拟结果与实验研究结果吻合较好。

参考文献 (19)

目录

    /

    返回文章
    返回