留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

饱和冲量及其等效方法在舱室内爆炸中的应用

朱凌 田岚仁 李德聪 王佳月 张霄

朱凌, 田岚仁, 李德聪, 等. 饱和冲量及其等效方法在舱室内爆炸中的应用[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01876
引用本文: 朱凌, 田岚仁, 李德聪, 等. 饱和冲量及其等效方法在舱室内爆炸中的应用[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01876
ZHU L, TIAN L R, LI D C, et al. Saturated impulse and application of saturation equivalent method in cabin explosion[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01876
Citation: ZHU L, TIAN L R, LI D C, et al. Saturated impulse and application of saturation equivalent method in cabin explosion[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01876

饱和冲量及其等效方法在舱室内爆炸中的应用

doi: 10.19693/j.issn.1673-3185.01876
基金项目: 国家自然科学基金资助项目(51579199)
详细信息
    作者简介:

    朱凌,男,1962年生,博士,教授。研究方向:船海结构安全可靠性与完整性,舰船抗爆抗冲击。E-mail:zl79111@hotmail.com

    田岚仁,男,1994年生,博士生。研究方向:舰船抗爆抗冲击。E-mail:libratian@126.com

    李德聪,男,1979年生,博士,高级工程师。研究方向:舰船抗爆抗冲击。E-mail:simon_ldc@sina.com

    通讯作者:

    朱凌

  • 中图分类号: U661.4;O383.3

Saturated impulse and application of saturation equivalent method in cabin explosion

  • 摘要:   目的  对于受到爆炸脉冲载荷冲击作用的船体结构,基于饱和冲量现象的相关研究表明,仅根据最大载荷幅值和脉冲总冲量来设计船体结构是不合理的,需探究工程应用中的饱和冲量现象。  方法  首先,总结饱和冲量概念的提出及研究发展;然后,以舱室内爆炸为典型算例,分析内爆炸载荷的曲线特性及结构响应特征;最后,基于饱和等效方法将复杂的内爆炸载荷等效为矩形脉冲载荷,采用理论及数值方法对等效载荷进行计算。  结果  结果表明:在舱室内爆炸准静态超压情况下普遍存在饱和冲量现象,实际工程应用中爆炸载荷会对结构造成较大的塑性变形,通常超过10倍板厚;而运用基于饱和冲量的等效方法分析,所得结果与数值仿真结果的误差小于10%。  结论  运用此方法可更准确地得出结构塑性动力响应结果,在结构抗冲击设计优化时,还可减少繁琐的复杂非线性数值计算,使设计更高效。
  • 图  1  简化后的舱室结构计算模型

    Figure  1.  Simplified calculation model of cabin structure

    图  2  舱室内爆炸有限元模型

    Figure  2.  Finite element model of cabin explosion

    图  3  横舱壁冲击波压力时历曲线(20 kg TNT)

    Figure  3.  Pressure time histories of blast wave on the transverse bulkhead (20 kg TNT)

    图  4  压力峰值和准静态超压随炸药量的变化曲线

    Figure  4.  Variation of pressure peak values and quasi-static overpressure with explosive mass

    图  5  不同炸药量下冲量时历曲线

    Figure  5.  Time histories of impulse with different explosive mass

    图  6  不同炸药量下板中心挠度时历曲线

    Figure  6.  Deflection time histories of plate center with different explosive mass

    图  7  饱和时间随炸药量的变化曲线

    Figure  7.  Variation of saturated time with explosive mass

    图  8  饱和冲量与板的中心最大挠度随炸药量的变化曲线

    Figure  8.  Variation of saturated impulse and maximum deflection of plate center with explosive mass

    图  9  饱和等效方法示意图

    Figure  9.  Schematic diagram of saturated equivalent method

    图  10  矩形板有限元计算模型

    Figure  10.  Finite element model of rectangular plate

    图  11  采用不同方法计算的矩形板最大挠度

    Figure  11.  Maximum deflection of rectangular plate calculated by different methods

    表  1  Cowper-Symonds强度模型参数

    Table  1.   Parameters of Cowper-Symonds strength model

    参数数值
    初始屈服应力${\sigma _0}$/MPa 235
    硬化模量${E_{\rm{h}}}$/MPa 250
    硬化指数 1
    常数D/s−1 40.4
    常数q 5
    下载: 导出CSV

    表  2  TNT炸药JWL状态方程参数

    Table  2.   Parameters of TNT explosives JWL state equation

    参数数值参数数值
    A1/kPa3.737 7×108w0.35
    A2/kPa3.747 1×106C-J爆速DCJ/(m·s−1)6.93e×103
    R14.15E0/(kJ·m−3)6.0×106
    R20.9C-J爆压PCJ/kPa2.1×107
    下载: 导出CSV

    表  3  等效后的矩形脉冲载荷参数

    Table  3.   Parameters of equivalent rectangular pulse loading

    m/kgP0/kPatsat/msIsat'/(kPa·ms)2tmean/msPe/kPa
    10 2 013.16 7.469 0 4 478.844 4 6.688 7 669.62
    15 3 116.80 7.207 4 6 666.832 4 6.576 2 1 013.78
    20 3 822.60 7.029 4 8 779.069 3 6.301 1 1 393.26
    25 4 745.55 6.926 9 10 769.058 9 6.130 9 1 756.51
    下载: 导出CSV
  • 侯海量, 朱锡, 梅志远. 舱内爆炸载荷及舱室板架结构的失效模式分析[J]. 爆炸与冲击, 2007, 27(2): 151–158. doi:  10.3321/j.issn:1001-1455.2007.02.010

    HOU H L, ZHU X, MEI Z Y. Study on the blast load and failure mode of ship structure subject to internal explosion[J]. Explosion and Shock Waves, 2007, 27(2): 151–158 (in Chinese). doi:  10.3321/j.issn:1001-1455.2007.02.010
    李德聪, 段宏, 吴国民, 等. 船内爆炸载荷特性及对舰船结构毁伤研究综述[J]. 中国舰船研究, 2018, 13(1): 7–16. doi:  10.3969/j.issn.1673-3185.2018.01.002

    LI D C, DUAN H, WU G M, et al. Advances in the research of warship structural damage due to inner explosion[J]. Chinese Journal of Ship Research, 2018, 13(1): 7–16 (in Chinese). doi:  10.3969/j.issn.1673-3185.2018.01.002
    姚熊亮, 屈子悦, 姜子飞, 等. 舰船舱内爆炸载荷特征与板架毁伤规律分析[J]. 中国舰船研究, 2018, 13(3): 140–148.

    YAO X L, QU Z Y, JIANG Z F, et al. Analysis on characteristics of blast loading and stiffened plate damage due to internal blast in ship[J]. Chinese Journal of Ship Research, 2018, 13(3): 140–148 (in Chinese).
    MARTIN J B, SYMONDS P S. Mode approximations for impulsively loaded rigid plastic structures[C]. Journal of the Engineering Mechanics Division. 1966, 92(5): 43-66.
    LEE L S S. Mode responses of dynamically loaded structures[J]. Journal of Applied Mechanics, 1972, 39(4): 904–910. doi:  10.1115/1.3422889
    JONES N. A theoretical study of the dynamic plastic behavior of beams and plates with finite-deflections[J]. International Journal of Solids and Structures, 1971, 7(8): 1007–1029. doi:  10.1016/0020-7683(71)90078-3
    PERRONE N, BHADRA P. A simplified method to account for plastic rate sensitivity with large deformations[J]. Journal of Applied Mechanics, 1979, 46(4): 811–816. doi:  10.1115/1.3424659
    SYMONDS P S. The optimal mode in the mode approximation technique[J]. Mechanics Research Communications, 1980, 7(1): 1–6. doi:  10.1016/0093-6413(80)90015-4
    YMONDS P S, WIERZBICKI T. Membrane mode solutions for impulsively loaded circular plates[J]. Journal of Applied Mechanics, 1979, 46(1): 58–64. doi:  10.1115/1.3424528
    BAKER W E. Approximate techniques for plastic deformation of structures under impulsive loading, III[J]. Shock & Vibration Digest, 1982, 14(11): 3–11.
    NURICK G N, MARTIN J B. Deformation of thin plates subjected to impulsive loading—A review: part I: theoretical considerations[J]. International Journal of Impact Engineering, 1989, 8(2): 159–170. doi:  10.1016/0734-743X(89)90014-6
    JONES N. Dynamic inelastic response of strain rate sensitive ductile plates due to large impact, dynamic pressure and explosive loadings[J]. International Journal of Impact Engineering, 2014, 74: 3–15. doi:  10.1016/j.ijimpeng.2013.05.003
    黄震球. 固支加筋方板的大挠度塑性动力响应[J]. 固体力学学报, 1995, 16(1): 7–12.

    HUANG Z Q. Large deflection dynamic plastic response of clamped square plates with stiffeners[J]. Acta Mechanica Solida Sinica, 1995, 16(1): 7–12 (in Chinese).
    颜丰, 刘敬喜. 爆炸载荷下固支矩形板的大挠度塑性动力响应[J]. 中国舰船研究, 2013, 8(1): 47–53.

    YAN F, LIU J X. The large deflection dynamic plastic response of rectangular plates subjected to blast load[J]. Chinese Journal of Ship Research, 2013, 8(1): 47–53 (in Chinese).
    吴有生, 彭兴宁, 赵本立. 爆炸载荷作用下舰船板架的变形与破损[J]. 中国造船, 1996(2): 86–92.

    WU Y S, PENG X N, ZHAO B L. Plastic deformation and damage of naval panels subjected to explosion loading[J]. Shipbuildig of China, 1996(2): 86–92 (in Chinese).
    牟金磊, 朱锡, 张振华, 等. 爆炸冲击作用下加筋板结构变形研究[J]. 海军工程大学学报, 2007, 19(6): 12–16. doi:  10.3969/j.issn.1009-3486.2007.06.003

    MU J L, ZHU X, ZHANG Z H, et al. A study on deformation of blast-loaded stiffened plates[J]. Journal of Naval University of Engineering, 2007, 19(6): 12–16 (in Chinese). doi:  10.3969/j.issn.1009-3486.2007.06.003
    ZHAO Y P, YU T X, FANG J. Large dynamic plastic deflection of a simply supported beam subjected to rectangular pressure pulse[J]. Archive of Applied Mechanics, 1994, 64(3): 223–232.
    ZHAO Y P, YU T X, FANG J. Saturation impulses for dynamically loaded structures with finite-deflections[J]. Structural Engineering and Mechanics, 1995, 3(6): 583–592. doi:  10.12989/sem.1995.3.6.583
    ZHU L. Dynamic inelastic response of ship plates in collision[D]. Glasgow: University of Glasgow, 1990.
    ZHU L, YU T X. Saturated impulse for pulse-loaded elastic-plastic square plates[J]. International Journal of Solids and Structures, 1997, 34(14): 1709–1718. doi:  10.1016/S0020-7683(96)00111-4
    席丰, 杨嘉陵. 强脉冲载荷作用下弹-塑性薄圆板的大挠度动力响应[J]. 爆炸与冲击, 2000, 20(4): 379–384. doi:  10.3321/j.issn:1001-1455.2000.04.016

    XI F, YANG J L. Dynamic response analysis of elastic-plastic thin circular plates under impulse loading with considenration of large deflection[J]. Explosion and Shock Waves, 2000, 20(4): 379–384 (in Chinese). doi:  10.3321/j.issn:1001-1455.2000.04.016
    席丰, 张云. 脉冲载荷作用下钢梁动力响应及反常行为的应变率效应[J]. 爆炸与冲击, 2012, 32(1): 34–42. doi:  10.3969/j.issn.1001-1455.2012.01.006

    XI F, ZHANG Y. The effects of strain rate on the dynamic response and abnormal behavior of steel beams under pulse loading[J]. Explosion and Shock Waves, 2012, 32(1): 34–42 (in Chinese). doi:  10.3969/j.issn.1001-1455.2012.01.006
    ZHU L, HE X, YU T X, et al. Scaling effect on saturated impulse for square plates under rectangular pulse loading[C]//Proceedings of the ASME 2016 35th International Conference on Ocean, Offshore and Arctic Engineering. Busan, South Korea: ASME, 2016.
    BAI X Y, ZHU L, YU T X. Saturated impulse for pulse-loaded rectangular plates with various boundary conditions[J]. Thin-Walled Structures, 2017, 119: 166–177. doi:  10.1016/j.tws.2017.03.030
    ZHU L, HE X, CHEN F L, et al. Effects of the strain rate sensitivity and strain hardening on the saturated impulse of plates[J]. Latin American Journal of Solids and Structures, 2017, 14(7): 1273–1292. doi:  10.1590/1679-78253664
    ZHU L, BAI X Y, YU T X. The saturated impulse of fully clamped square plates subjected to linearly decaying pressure pulse[J]. International Journal of Impact Engineering, 2017, 110: 198–207. doi:  10.1016/j.ijimpeng.2016.12.012
    BAI X Y, ZHU L, YU T X. Saturated impulse for fully clamped square plates under blast loading[J]. International Journal of Mechanical Sciences, 2018, 146-147: 417–431. doi:  10.1016/j.ijmecsci.2017.08.047
    TIAN L R, CHEN F L, ZHU L, et al. Saturated analysis of pulse-loaded beams based on membrane factor method[J]. International Journal of Impact Engineering, 2019, 131: 17–26. doi:  10.1016/j.ijimpeng.2019.04.021
    YOUNGDAHL C K. Correlation parameters for eliminating the effect of pulse shape on dynamic plastic deformation[J]. Journal of Applied Mechanics, 1970, 37(3): 744–752. doi:  10.1115/1.3408605
    陈攀, 刘志忠. 舱室内爆载荷作用下夹层板优化设计研究[J]. 舰船科学技术, 2016, 38(11): 14–20.

    CHEN P, LIU Z Z. Optimal design of anti-shock capability of sandwich panel under explosion loading inside closed cabin[J]. Ship Science and Technology, 2016, 38(11): 14–20 (in Chinese).
    鄢顺伟, 杜茂华, 王伟力, 等. 战斗部内爆对舰艇舱室的毁伤效应仿真[J]. 海军航空工程学院学报, 2013, 28(2): 181–188. doi:  10.3969/j.issn.1673-1522.2013.02.020

    YAN S W, DU M H, WANG W L, et al. Damage effect simulation of warhead inner explosion in warship cabin[J]. Journal of Naval Aeronautical and Astronautical, 2013, 28(2): 181–188 (in Chinese). doi:  10.3969/j.issn.1673-1522.2013.02.020
    JONES N. Structural impact[M]. Cambridge: Cambridge University Press, 1989.
    ZHU L, FAULKNER D, ATKINS A G. The impact of rectangular plates made from strain-rate sensitive materials[J]. International Journal of Impact Engineering, 1994, 15(3): 245–255. doi:  10.1016/S0734-743X(05)80016-8
  • [1] 王喜梦, 刘均, 陈长海, 程远胜, 张攀2020: 近距空爆载荷下钢板/聚脲复合结构动响应特性仿真 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01833
    [2] 蔡伟, 朱凌2020: 冰碰载荷下船体板弹塑性动力响应分析 , 中国舰船研究, 15, 155-161. doi: 10.19693/j.issn.1673-3185.01804
    [3] 陆潇杨, 刘志全, 高迪驹, 褚建新2020: 带漂角和输入饱和的水面船舶航向控制 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01775
    [4] 焦立启, 张权, 李茂, 张磊, 张春辉2020: 典型舱内爆炸载荷对加筋板的毁伤特性 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01841
    [5] 苏智伟, 黄修长, 吴静波, 华宏星2019: 含负刚度动力吸振的混合隔振系统振动冲击响应特性分析 , 中国舰船研究, 14, 59-65. doi: 10.19693/j.issn.1673-3185.01194
    [6] 赵鹏铎, 黄松, 尹建平, 张磊, 黄阳洋, 徐豫新2019: 泄爆薄板厚度对舱室破坏及舱内载荷影响的数值仿真 , 中国舰船研究, 14, 58-65,82. doi: 10.19693/j.issn.1673-3185.01295
    [7] 金勇, 邝俊信, 田相玉, 劳坤胜, 欧阳武, 刘正林2019: 饱和含液穿孔板条水润滑艉轴承减振性能研究 , 中国舰船研究, 14, 58-63. doi: 10.19693/j.issn.1673-3185.01402
    [8] 吴广明, 许亚东, 李正国, 王志凯, 姚熊亮2019: 舱室复合材料结构的抗爆性能数值仿真 , 中国舰船研究, 14, 14-21. doi: 10.19693/j.issn.1673-3185.01222
    [9] 李春鹏, 张攀, 刘均, 程远胜2018: 空爆载荷下功能梯度泡沫铝夹层板动响应数值仿真 , 中国舰船研究, 13, 77-84. doi: 10.19693/j.issn.1673-3185.01172
    [10] 徐维铮, 吴卫国2018: 多舱室内气云爆炸爆轰产物运动二维数值模拟 , 中国舰船研究, 13, 92-99. doi: 10.19693/j.issn.1673-3185.01186
    [11] 吴林杰, 朱锡, 侯海量, 陈长海2016: 冲击载荷下舰船夹芯舱室结构动态响应及吸能特性 , 中国舰船研究, 11, 70-76,96. doi: 10.3969/j.issn.1673-3185.2016.06.011
    [12] 葛亮, 田正东, 袁利毫2015: 气垫船碰撞冲击动力学响应三维数值模拟 , 中国舰船研究, 10, 57-62. doi: 10.3969/j.issn.1673-3185.2015.03.010
    [13] 颜丰, 刘敬喜2013: 爆炸载荷下固支矩形板的大挠度塑性动力响应 , 中国舰船研究, 8, 47-53. doi: 10.3969/j.issn.1673-3185.2013.01.008
    [14] 郑杰, 谢伟, 李慧敏2011: 三维流体动力方法计算穿浪双体船的船体响应 , 中国舰船研究, 6, 29-34. doi: 10.3969/j.issn.1673-3185.2011.02.006
    [15] 李伦, 许建, 黄国兵2011: 外压作用下带初挠度加筋圆柱壳极限载荷及变形研究 , 中国舰船研究, 6, 13-16. doi: 10.3969/j.issn.1673-3185.2011.06.003
    [16] 王佳颖, 张世联, 徐敏2010: 接触爆炸下舰船强力甲板塑性动态响应特性研究 , 中国舰船研究, 5, 10-14,21. doi: 10.3969/j.issn.1673-3185.2010.05.003
    [17] 陈晓洪, 张光辉, 刘辉2010: 潜艇舱室内的火灾模拟及消防决策研究 , 中国舰船研究, 5, 34-38. doi: 10.3969/j.issn.1673-3185.2010.01.008
    [18] 刘敬喜, 刘尧, 李威2010: 爆炸载荷作用下单向加筋板的塑性动力响应分析 , 中国舰船研究, 5, 6-9. doi: 10.3969/j.issn.1673-3185.2010.05.002
    [19] 孔祥韶, 吴卫国, 李晓彬, 徐双喜, 黄涛2009: 舰船舱室内部爆炸的数值模拟研究 , 中国舰船研究, 4, 7-11. doi: 10.3969/j.issn.1673-3185.2009.04.002
    [20] 彭英, 杨平2006: 面内流固冲击载荷下加筋板的动力响应 , 中国舰船研究, 1, 16-20. doi: 10.3969/j.issn.1673-3185.2006.04.004
  • 加载中
图(11) / 表 (3)
计量
  • 文章访问数:  35
  • HTML全文浏览量:  4
  • PDF下载量:  1142
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-05
  • 修回日期:  2020-03-22
  • 网络出版日期:  2020-12-10

饱和冲量及其等效方法在舱室内爆炸中的应用

doi: 10.19693/j.issn.1673-3185.01876
    基金项目:  国家自然科学基金资助项目(51579199)
    作者简介:

    朱凌,男,1962年生,博士,教授。研究方向:船海结构安全可靠性与完整性,舰船抗爆抗冲击。E-mail:zl79111@hotmail.com

    田岚仁,男,1994年生,博士生。研究方向:舰船抗爆抗冲击。E-mail:libratian@126.com

    李德聪,男,1979年生,博士,高级工程师。研究方向:舰船抗爆抗冲击。E-mail:simon_ldc@sina.com

    通讯作者: 朱凌
  • 中图分类号: U661.4;O383.3

摘要:   目的  对于受到爆炸脉冲载荷冲击作用的船体结构,基于饱和冲量现象的相关研究表明,仅根据最大载荷幅值和脉冲总冲量来设计船体结构是不合理的,需探究工程应用中的饱和冲量现象。  方法  首先,总结饱和冲量概念的提出及研究发展;然后,以舱室内爆炸为典型算例,分析内爆炸载荷的曲线特性及结构响应特征;最后,基于饱和等效方法将复杂的内爆炸载荷等效为矩形脉冲载荷,采用理论及数值方法对等效载荷进行计算。  结果  结果表明:在舱室内爆炸准静态超压情况下普遍存在饱和冲量现象,实际工程应用中爆炸载荷会对结构造成较大的塑性变形,通常超过10倍板厚;而运用基于饱和冲量的等效方法分析,所得结果与数值仿真结果的误差小于10%。  结论  运用此方法可更准确地得出结构塑性动力响应结果,在结构抗冲击设计优化时,还可减少繁琐的复杂非线性数值计算,使设计更高效。

English Abstract

朱凌, 田岚仁, 李德聪, 等. 饱和冲量及其等效方法在舱室内爆炸中的应用[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01876
引用本文: 朱凌, 田岚仁, 李德聪, 等. 饱和冲量及其等效方法在舱室内爆炸中的应用[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01876
ZHU L, TIAN L R, LI D C, et al. Saturated impulse and application of saturation equivalent method in cabin explosion[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01876
Citation: ZHU L, TIAN L R, LI D C, et al. Saturated impulse and application of saturation equivalent method in cabin explosion[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01876
    • 船舶与海洋工程结构物在运行过程中会受到脉冲载荷的冲击作用,例如爆炸、砰击等。在结构响应研究和实际结构设计中,通常将爆炸载荷简化为具有初速度的冲击载荷或者脉冲载荷。而随着各种反舰武器性能的提高,其能够穿透船体外板在舱室内部发生爆炸,使得爆炸载荷对船体结构的威胁日益增加,因此研究约束空间(有限空间)内的抗爆、抗冲击影响越来越重要。

      然而,相比敞开环境下的爆炸,约束空间的内爆炸载荷远比自由场的空爆载荷复杂,原因在于有限空间内的爆炸冲击波的壁面反射会造成多峰效应,经反射波叠加,其爆炸峰值会高于敞开环境下的峰值,作用时间也显著变长,且通常会存在准静态超压的现象[1-3]。迄今,针对舱室内爆炸冲击波的多峰效应及其存在的准静态超压等现象的研究成果颇丰,但尚未完全探明内爆炸冲击载荷作用下结构响应的力学机理。

      20世纪50年代以来,许多学者研究了脉冲载荷作用下的结构动态塑性响应特性,所采用的大部分理论方法也都是基于理想的刚塑性(R-PP)假设,这是因为在忽略弹性变形时可显著简化理论模型,从而得到塑性变形的解析表达式。目前,在大挠度下的结构动态响应刚塑性理论求解方面应用得较广泛的是模态近似分析法[4],运用该方法避免了考虑瞬态响应阶段的移行铰模型在求解中塑性铰两侧物理量不连续性和非线性的问题,使理论分析变得更简单。之后,国外许多学者又进一步发展了模态近似分析法[5-12]。在国内,黄震球[13]和颜丰等[14]基于动量和动量矩守恒发展了加筋板结构塑性变形的解析方法,所得结果经与试验结果对比,吻合较好。吴有生等[15]和牟金磊等[16]采用能量法考虑大变形时的应变关系及中面膜力的影响,求解了船体板结构塑性大变形。

      上述研究爆炸冲击载荷作用下的结构响应特性针对的脉冲载荷持续时间较短,即结构变形将在载荷结束之前结束。然而,如果载荷持续时间较长,结构动力响应中将出现饱和冲量现象,从而影响采用理论方法计算所得结果的合理性和准确性。因此,本文将系统地总结饱和冲量概念的提出及研究发展成果,并以舱室内爆炸为典型算例,分析内爆炸载荷曲线的特性和结构响应特征,基于饱和等效方法将复杂的内爆炸载荷等效为矩形脉冲载荷,采用理论及数值方法计算等效载荷,给出相关的设计曲线和计算公式。通过研究舱室内爆炸作用下的饱和冲量现象及饱和等效方法,以更准确地得到结构塑性动力响应结果,并预测结构变形的终止时间与载荷间的关系。运用上述方法还可减少频繁的复杂非线性数值计算,更有效地开展船体结构的抗冲击设计优化。

    • 20世纪90年代,Zhao等[17]定义了刚塑性结构在中等强度脉冲载荷作用下的一种特征现象—“饱和冲量”,其反映的是板在塑性动力响应中的一个特性,即板受到强烈的横向压力脉冲载荷作用时会产生较大的变形,而大变形诱发的膜力对板的承载能力起到了增强作用。若板受到足够长时间的矩形压力脉冲作用,就只有脉冲载荷达到饱和时间前的冲量对最大及最终变形有影响,之后的加载脉冲不会进一步增加变形量。此后,Zhao等[18]又将饱和冲量概念拓展到基于理想的刚塑性模型的简支圆板、简支和固支方板以及圆柱壳中,利用模态分析法得到了封闭形式的解析解。然而,“饱和冲量”概念只对应于理想的刚塑性结构最大变形情况,未考虑弹性影响。鉴此,Zhu[19]提出一种运用有限差分方法来分析固支方板的弹塑性响应数值程序,基于该程序,Zhu等[20]提出了分别对应于“最大挠度”和“最终挠度”的“饱和冲量”,从而完善了基于刚塑性和弹塑性模型的饱和冲量现象的定义。

      针对基于刚塑性−弹塑性饱和冲量现象的问题,席丰等[21]利用最小加速度原理,通过建立数值方法,分析了脉冲载荷作用下简支圆板的动力响应,指出在高载荷范围内也存在脉冲载荷作用下的“饱和冲量”现象。此后,席丰等[22]又分析了脉冲载荷作用下的钢梁动力响应及反常行为的应变率效应,指出在脉冲载荷作用下发生塑性变形的钢梁总是存在“饱和冲量”现象,且发生时与载荷强度及载荷作用时间相关。

      近年来,武汉理工大学朱凌教授的研究团队针对饱和冲量问题开展了进一步研究,并更系统地予以了分析,例如:尺度效应对方板饱和冲量的影响[23];长宽比和边界条件对矩形板饱和冲量的影响[24];材料应变率敏感性和应变硬化对饱和冲量的影响[25];不同脉冲载荷作用下的方板饱和冲量[26-27];考虑移行铰的方板饱和冲量[27];同时考虑瞬态响应阶段和准确屈服面的梁饱和冲量[28]。此外,还基于饱和冲量的研究,提出了物理意义更明确且计算结果更准确的脉冲载荷等效方法[26-27],改进了Youngdahl[29]于70年代提出的半经验等效方法。然而,上述研究尚未在饱和冲量的应用层面开展更深入的探讨。

    • 结构响应的理论分析通常是在假定载荷已知的基础上开展的,不同于敞开环境下的爆炸,约束空间内爆炸载荷曲线通常很难被写成某种较为通用的函数形式。因此,可以针对需要分析的模型,利用有限元软件计算出爆炸所产生的脉冲载荷。此外,过去几十年以来,学者们开展了大量的爆炸载荷实验和数值研究,若针对的是某一系列结构,可通过前人的研究成果确定较为合理的载荷曲线形式,并直接对结构进行动力响应的理论分析。本文采用AUTODYN有限元软件计算舱室内爆炸载荷及结构动态响应,选取图1所示舰船右舷典型舱室的结构计算模型[1],模型的尺寸为5 000 mm×3 000 mm×2 500 mm。

      图  1  简化后的舱室结构计算模型

      Figure 1.  Simplified calculation model of cabin structure

      为简化研究和定量分析,仅采用无加筋板进行原理性探讨,并参考文献[30-31]将舱室内爆炸简化为炸药在舱室内中心爆炸。如图2所示,可变形壁面(横舱壁)为3 000 mm×2 500 mm×14 mm(2L(长)×2B(宽)×H(厚))的矩形板,其他5个壁面为刚性壁面,设定舱室为完全封闭空间。矩形板的4个端部为固支,空气域网格划分为110×60×50,网格数330 000,采用高阶单物质Euler-FCT求解器。矩形板网格划分为120×100,网格数12 000,板与空气采用完全耦合实现相互作用。本文仿真计算选取4组TNT球状炸药。

      图  2  舱室内爆炸有限元模型

      Figure 2.  Finite element model of cabin explosion

    • 根据文献[1]中的材料参数,选取Q235钢作为矩形板材料,密度为7 800 kg/m3,杨氏模量E=210 GPa,泊松比υ=0.3。采用双线性弹塑性本构模型,材料的应变率效应由Cowper-Symonds模型描述,相关参数见表1

      表 1  Cowper-Symonds强度模型参数

      Table 1.  Parameters of Cowper-Symonds strength model

      参数数值
      初始屈服应力${\sigma _0}$/MPa 235
      硬化模量${E_{\rm{h}}}$/MPa 250
      硬化指数 1
      常数D/s−1 40.4
      常数q 5

      自由空气采用的是γ律状态方程描述。TNT炸药爆轰产物的压力则由JWL状态方程描述:

      $$P = {A_1}\left(1 - \frac{w}{{{R_1}V}}\right){{\rm{e}}^{( - {R_1}V)}} + {A_2}\left(1 - \frac{w}{{{R_2}V}}\right){{\rm{e}}^{( - {R_2}V)}} + \frac{{w{E_0}}}{V}$$ (1)

      式中:P为爆轰产物的压力;V=ρ0/ρ,为相对体积,其中ρ0为炸药初始密度,ρ0=1 630 kg/m3ρ为炸药爆轰产物的密度;E0为单位体积炸药的初始内能;A1,A2,R1,R2,w为与炸药性质相关的材料常数;e为自然数。其他参数见表2

      表 2  TNT炸药JWL状态方程参数

      Table 2.  Parameters of TNT explosives JWL state equation

      参数数值参数数值
      A1/kPa3.737 7×108w0.35
      A2/kPa3.747 1×106C-J爆速DCJ/(m·s−1)6.93e×103
      R14.15E0/(kJ·m−3)6.0×106
      R20.9C-J爆压PCJ/kPa2.1×107
    • 图3所示为采用20 kg炸药在内爆炸作用下横舱壁的冲击波压力时历曲线。由图可见:壁面反射作用使冲击波压力含有多个峰值;因舱室模型未设置泄爆口,内爆炸后产生的准静态超压Ps始终作用在结构上,板上不同位置的压力峰值有所差别,但准静态压力值基本一致。由于受到复杂的反射波影响,约束空间内爆炸下的最大压力峰值Pm尚无可广泛使用的经验公式。对于舱室内爆炸准静态超压,文献[2]总结了4种经验公式。图4为本文有限元计算得到的准静态超压值(带圆形黑色实线)与这4种经验公式的对比。图中红色叉点在横坐标上的时间即对应于炸药量下的饱和时间tsat。由图可见:本文计算结果与Moir和Carlson经验公式得到的结果较为接近。同时,还可发现最大压力峰值(带矩形黑色实线)和准静态超压随着炸药量的增加都大致呈线性上升的趋势。

      图  3  横舱壁冲击波压力时历曲线(20 kg TNT)

      Figure 3.  Pressure time histories of blast wave on the transverse bulkhead (20 kg TNT)

      图  4  压力峰值和准静态超压随炸药量的变化曲线

      Figure 4.  Variation of pressure peak values and quasi-static overpressure with explosive mass

      图5所示为3组不同炸药量下的爆炸载荷冲量I随时间的变化曲线。由图可见,在反射波作用下,冲量经过初始时间段的波动后大致呈线性上升的趋势,这主要是后期准静态超压所致。图6所示为3组不同炸药量下板中心挠度W0随时间的变化曲线。结合图3可发现,板在饱和时间tsat(蓝色点划线所示)后仍受到大于板的静态塑性极限载荷的压力作用,冲量也始终处于上升的趋势,但板的变形已达到最大值且未再增加。这说明只有早期压力脉冲载荷所形成的冲量对板的实际变形有效,后期的冲击压力并未引起板的变形值进一步增加,此即为“饱和冲量”现象,对应时间tsat为饱和时间,对应冲量Isat为饱和冲量。

      图  5  不同炸药量下冲量时历曲线

      Figure 5.  Time histories of impulse with different explosive mass

      图  6  不同炸药量下板中心挠度时历曲线

      Figure 6.  Deflection time histories of plate center with different explosive mass

      上述分析表明,舱室内爆炸下结构存在“饱和冲量”现象,若预测由此导致的结构大变形应基于“饱和冲量”的相关方法。但是,在弹性效应的影响下,由永久挠度来确定tsat的精确值较难,故可选择最大挠度进行分析。需要指出的是,基于最大挠度分析的tsat值通常偏小,但随着载荷幅值的增大,即结构塑性变形的增加,2种挠度确定的tsat的差值会逐渐缩小。

      图7图8分别给出了饱和时间tsat、饱和冲量Isat和板中心点最大挠度$ {W}_{0}^{\mathrm{s}\mathrm{a}\mathrm{t}} $随炸药量m的变化关系。由图可知:饱和时间随炸药量的增加逐渐减少,但减少的趋势逐渐变缓,这与刚塑性的理论分析中饱和时间不随压力峰值变化有所出入,主要是由弹性效应及载荷曲线发生变化所致。此外,最大挠度和饱和冲量则随炸药量的增加而大致呈线性上升的趋势。尽管本文计算算例中的载荷曲线较复杂,但两者的趋势与文献中采用的矩形脉冲[24]、线性衰减脉冲[26]或者线性上升指数衰减脉冲[27]的图线趋势一致。

      图  7  饱和时间随炸药量的变化曲线

      Figure 7.  Variation of saturated time with explosive mass

      图  8  饱和冲量与板的中心最大挠度随炸药量的变化曲线

      Figure 8.  Variation of saturated impulse and maximum deflection of plate center with explosive mass

    • 在工程实际中,脉冲形状通常十分复杂,例如舱室内爆炸的压力时历曲线具有多个峰值且形状不规则,这对于求解脉冲载荷作用下的结构响应和结构设计造成了极大困难。而工程中有很多结构需优化设计和反复计算,若每次都采用非线性有限元计算,计算资源消耗太大。因此,基于有限元计算或实验测量等方法得到载荷时历曲线,可将其引入到载荷等效方法中,结合理论分析快速计算结构响应,以减少繁琐的复杂非线性数值计算。文献[29]提出的消除脉冲载荷形状影响的等效方法,可将一个任意形状的脉冲由有效冲量Ie和有效载荷Pe这2个量来表征:

      $${I_{\rm{e}}} = \int_{{t_{\rm{y}}}}^{{t_{\rm{f}}}} {P(t){\rm{d}}t} $$ (2)
      $${t_{{\rm{mean}}}} = \frac{{\int_{{t_{\rm{y}}}}^{{t_{\rm{f}}}} {(t - {t_{\rm{y}}})P(t){\rm{d}}t} }}{{{I_{\rm{e}}}}}$$ (3)
      $${P_{\rm{e}}} = \frac{{{I_{\rm{e}}}}}{{2{t_{{\rm{mean}}}}}}$$ (4)

      式中:P(t)为脉冲载荷;tytf分别为塑性变形开始和结束时刻;2tmean为等效后的矩形脉冲长度。通过式(3)找出该有效脉冲的重心,作为等效矩形脉冲的中心,再由式(4)确定等效矩形脉冲的幅值,从而对原脉冲载荷曲线进行“掐头去尾”后得到有效冲量Ie

      在实际应用中,通常很难确定tf,文献[29]提出采用经验公式(5)进行“去尾”处理。

      $$({t_{\rm{f}}} - {t_{\rm{y}}}){P_{\rm{y}}} = \int_{{t_{\rm{y}}}}^{{t_{\rm{f}}}} {P(t){\rm{d}}t} $$ (5)

      式中,Py为结构的塑性极限载荷。

      在过去几十年内,Youngdahl这种半经验等效方法被广泛用来简化处理实际的工程问题。然而,该等效方法并未给出理论依据。近几年,文献[24-25]在饱和冲量研究的基础上,提出了物理意义更明确,且计算结果更准确的解析脉冲载荷等效方法,从而改进了Youngdahl半经验等效方法。改进后的等效方法由饱和时间tsat代替式(5)中的tf,划分出脉冲中的有效部分,得到饱和冲量${I^{{\rm{sat}}}}$,然后利用文献[29]提出的“重心等效”方法获得等效的矩形脉冲载荷。

      在本文算例中,采用饱和等效方法对实际舱室的内爆炸载荷进行了等效,原理如图9所示。图中,对于固支矩形板,静态塑性极限载荷由式(6)计算[32]

      图  9  饱和等效方法示意图

      Figure 9.  Schematic diagram of saturated equivalent method

      $${P_{\rm{y}}} = \frac{{12{M_0}}}{{{{(\sqrt {3 + {\beta ^2}} - \beta )}^2}{B^2}}}$$ (6)

      将塑性极限弯矩${M_0} = {{{\sigma _0}{H^2}}/4}$和宽长比$\beta {\rm{ = }} B/L$代入,可以计算得到本文模型的静态塑性极限载荷为Py=74.6 kPa。

      图9可见,载荷曲线初始段上升得很快。为简化分析,ty可以取为0,计算结果的误差将随着载荷幅值(炸药量)的增大而减小。Youngdahl等效方法建议采用式(5)近似地确定塑性变形结束时刻tf,但从本文算例可见,在存在准静态超压的情况下,整条载荷曲线的值长时间处于静态塑性极限载荷Py之上,使得式(5)无法成立。这说明Youngdahl等效方法不适用于舱室内爆炸载荷的等效计算。表3给出了等效后的矩形脉冲载荷参数。需要注意的是,基于等效后的载荷采用理论方法或者有限元法计算响应时,应考虑是否有大气压的影响。表3中,P0为实际压力峰值,Isat'为减去大气压影响后的饱和冲量。

      表 3  等效后的矩形脉冲载荷参数

      Table 3.  Parameters of equivalent rectangular pulse loading

      m/kgP0/kPatsat/msIsat'/(kPa·ms)2tmean/msPe/kPa
      10 2 013.16 7.469 0 4 478.844 4 6.688 7 669.62
      15 3 116.80 7.207 4 6 666.832 4 6.576 2 1 013.78
      20 3 822.60 7.029 4 8 779.069 3 6.301 1 1 393.26
      25 4 745.55 6.926 9 10 769.058 9 6.130 9 1 756.51
    • 得到矩形等效载荷后,可非常方便地对结构响应进行理论分析及求解。白雪玉等[24]针对矩形脉冲载荷作用下矩形板的饱和冲量现象进行刚塑性理论分析,给出了板的无量纲饱和挠度及饱和时间的计算公式:

      $${\left( {\frac{{{W_0}}}{H}} \right)^{{\rm{sat}}}} = \frac{{3 - \xi }}{{{\xi ^2} - 3\xi + 3}}\left[ {\left( {3 - 2\xi } \right)\lambda - \frac{{1 + \eta }}{2}} \right]$$ (7)
      $${t^{{\rm{sat}}}} = \frac{{\rm{\pi }}}{{{a_3}}},a_3^{} = \sqrt {\frac{{{P_{{\rm{y}}0}}}}{{\mu H}}\frac{{2({\xi ^2} - 3\xi + 3)}}{{(2 - \xi )(3 - 2\xi )}}} $$ (8)

      式中:(W0/H)sat为无量纲饱和挠度;$\xi = \; \beta \sqrt {{\beta ^2} + 3} - {\beta ^2}$,为绞线位置参数;$\lambda {\rm{ = }}{P_{\rm{e}}}/{P_{{\rm{y}}0}}$,为无量纲载荷幅值,其中基准载荷${P_{{\rm{y}}0}}{\rm{ = }}12{M_0}/{B^2}$$\mu $为单位面积质量;$\eta $为边界条件参数,若是简支,$\eta = 0$,若是固支,$\eta = 1$

      式(7)和式(8)适用于2tmeantsat的情况,但等效后的矩形脉冲长度2tmean<tsat时,即在矩形脉冲载荷作用下不发生饱和冲量现象,此时,板的变形过程将分为2相:第1相中载荷一直作用在结构上,第2相中载荷消失而结构继续变形直至停止。此时,可见文献[32]给出的固支矩形板的最终挠度(W0/Hf结果:

      $$\begin{split} & \qquad\qquad\qquad\qquad\qquad{\left( {\frac{{{W_0}}}{H}} \right)^{\rm{f}}} = \\&\frac{{(3 \!-\! \xi )(\sqrt {1\! \!+\!\! 2({{(3\! -\! 2\xi )}^2}{\lambda ^2} \!\!-\!\! (3 \!- \!2\xi )\lambda )(1\! \!-\! \!\cos (2{a_3}{t_{{\rm{mean}}}}))} \! - \!\!1)}}{{2({\xi ^2} - 3\xi + 3)}} \end{split}$$ (9)

      然而,在爆炸载荷作用下,由于船体钢结构响应通常受应变率的影响很大,所以为得到更准确的结构大变形的解,分析时计及应变率效应具有工程实际意义。文献[33]假设应变率服从Cowper-Symonds模型,从对刚体撞击下的矩形板不同区域的应变率进行计算,给出了矩形板的整体应变率加强因子。为量化应变率效应对饱和冲量现象的影响,文献[25]针对脉冲加载的方板进一步开展了研究。这些方法也可扩展到脉冲加载的矩形板研究中。

      对于2tmeantsat的情况,假设应变率服从Cowper-Symonds模型,空间平均应变率为

      $${\dot \varepsilon _1} = {W_0}{\dot W_0}/{(\xi L)^2},{\dot \varepsilon _2} = {W_0}{\dot W_0}/{B^2}$$ (10)

      式中,${\dot \varepsilon _1}$${\dot \varepsilon _2}$分别为文献[33]图2中区域I和区域II的应变率;${\dot W_0}$为板中心点变形速度。由此,可求得矩形板的整体应变率强化因子n:

      $$\begin{split} & n = \left\{ {1 + {{\left( {\frac{{3\sqrt 3 }}{{32}}\frac{{{{(W_0^{{\rm{sat}}})}^{\rm{2}}}{a_3}}}{{{{(\xi L)}^{\rm{2}}}D}}} \right)}^{\frac{1}{q}}}} \right\}\frac{\xi }{2} +\\&\;\;\; \left\{ {1 + {{\left( {\frac{{3\sqrt 3 }}{{32}}\frac{{{{(W_0^{{\rm{sat}}})}^{\rm{2}}}{a_3}}}{{{B^2}D}}} \right)}^{\frac{1}{q}}}} \right\}\frac{{2{\rm{ - }}\xi }}{2} \end{split}$$ (11)

      对于2tmean<tsat的情况,可得矩形板的整体应变率强化因子为

      $$\begin{split} & n = \left\{ {1 + {{\left( {\frac{4}{{3\sqrt 2 }}\frac{{W_0^{{\rm{sat}}}{P_{\rm{e}}}{t_{{\rm{mean}}}}}}{{{{(\xi L)}^{\rm{2}}}D\mu }}} \right)}^{\frac{1}{q}}}} \right\}\frac{\xi }{2} +\\&\;\;\; \left\{ {1 + {{\left( {\frac{4}{{3\sqrt 2 }}\frac{{W_0^{{\rm{sat}}}{P_{\rm{e}}}{t_{{\rm{mean}}}}}}{{{B^2}D\mu }}} \right)}^{\frac{1}{q}}}} \right\}\frac{{2{\rm{ - }}\xi }}{2} \end{split}$$ (12)

      则考虑了应变率效应的无量纲饱和挠度为

      $${\left( {\frac{{{W_0}}}{H}} \right)^{{\rm{sat'}}}} = \frac{{3 - \xi }}{{{\xi ^2} - 3\xi + 3}}\left[ {\left( {3 - 2\xi } \right)\frac{\lambda }{n} - 1} \right]$$ (13)

      因此,可得等效后的矩形脉冲载荷作用下船体结构的最大挠度理论结果。由于文献[20]前期研究中指出饱和时间与板厚无关,所以可将等效载荷代入到上述理论公式中进行结构的板厚设计,从而可减少复杂的非线性数值计算。

    • 鉴于AUTODYN软件中不便于在板结构上准确施加均布载荷,本文采用Abaqus/Explicit软件对板结构施加等效后的矩形脉冲载荷进行仿真计算。计算时,采用与AUTODYN中相同的模型,但无需设置空气域。矩形板网格取60×50,网格类型为S4R四节点减缩积分壳单元,如图10所示。矩形板材料设置为AUTODYN软件计算中相同的材料参数,计及应变率效应。

      图  10  矩形板有限元计算模型

      Figure 10.  Finite element model of rectangular plate

    • 图11所示为采用不同方法计算得到的简化矩形板舱室横舱壁的最大挠度值Wsat随炸药量的变化曲线。图中,3种计算方法所得结果分别为:

      图  11  采用不同方法计算的矩形板最大挠度

      Figure 11.  Maximum deflection of rectangular plate calculated by different methods

      1) 采用AUTODYN模拟实际炸药爆炸过程得到的结构变形(实际载荷有限元计算值);

      2) 对板结构施加等效后的矩形脉冲载荷,采用Abaqus计算得到的结构变形(等效载荷有限元计算值);

      3) 对板结构施加等效后的矩形脉冲载荷,采用刚塑性理论方法计算得到的结构变形(等效载荷理论计算值)。

      图11可知,当对板结构施加等效后的矩形脉冲载荷时,采用有限元方法和理论方法计算的结果与实际结果十分接近,尤其是将等效载荷输入到有限元中计算时,与实际模拟结果的偏差小于10%;对于刚塑性理论方法,炸药量较小时,误差较大,而当变形接近或超过10倍板厚时,理论方法的预测结果较好。

      理论计算值与实际结果之间产生偏差的主要原因是:等效脉冲与实际脉冲载荷形状差异的影响;理论模型中未考虑弹性效应;分析过程采用的是方形近似屈服面等。

      然而,对于工程中的实际应用,理论公式在初步设计时优势更明显。

    • 在实际工程应用中,通常敞开环境下近距离爆炸载荷作用的时间很短,但对于约束空间爆炸而言,例如舰船舱室内的爆炸,因存在准静态超压,载荷作用持续时间会较长。此时,研究饱和冲量现象能更准确地给出结构塑性动力响应结果,以及预测结构变形的停止时间与载荷间的关系,从而更有效地指导船体结构的抗冲击设计。

      本文介绍了“饱和冲量”概念提出后近几年的主要研究进展,并以舱室内爆炸为典型算例,计算了有限空间内爆炸载荷下的船体结构响应,基于饱和等效方法将复杂的舱室内爆炸载荷等效为矩形脉冲载荷,运用理论和数值方法对等效载荷进行了计算,进一步验证了饱和等效方法在工程应用中的有效性。通过分析,得到以下结论:

      1) 炸药量(或载荷峰值)对舱室内爆炸下的结构饱和时间具有一定的影响,即随着载荷的增大,饱和时间将逐渐减小;饱和挠度和饱和冲量随炸药量的变化大致呈线性上升的趋势,这为工程应用提供了很大的便利。

      2) 对于舱室内爆炸载荷,在采用饱和等效法等效时,为简化分析,也可将塑性变形的开始时间ty取为0,计算结果的误差随载荷幅值(炸药量)的增大而减小;

      3) 将载荷等效为矩形脉冲载荷后,使用理论公式可快速计算,但应计及材料的应变率效应,否则理论计算的结果会偏大;

      4) Youngdahl等效方法在舱室内爆炸的计算中并不适用,而饱和等效方法提供了一种便利的解决途径。

参考文献 (33)

目录

    /

    返回文章
    返回