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

带环形密闭气囊弹体入水冲击过程的数值分析

陈洋 吴亮 曾国伟 周俊汝

引用本文:
Citation:

带环形密闭气囊弹体入水冲击过程的数值分析

    作者简介: 陈洋(1994-), 男, 硕士研究生;
    通讯作者: 吴亮, wuliangwust@sina.com
  • 基金项目: 国家自然科学基金项目 51479147
    湖北省自然科学基金项目 2014CFB822
    国家自然科学基金项目 11602178
    国家自然科学基金项目 51004079

  • 中图分类号: O39;V244

Numerical analysis of the water entry process of a projectile with a circular airbag

    Corresponding author: WU Liang, wuliangwust@sina.com ;
  • CLC number: O39;V244

  • 摘要: 针对带环形密闭气囊弹体入水冲击问题,基于LS-DYNA,运用控制体积法模拟环形密闭气囊,结合流固耦合算法,模拟了某弹体及附带环形密闭气囊入水过程。将入水过程分为弹体砰水、气囊着水、入水减速、水中悬停、缓慢上浮、上浮出水、水面漂浮7个主要阶段,对比分析了垂直与倾斜入水过程中不同阶段弹体和气囊的姿态变化、减速特性及入水深度等特征的异同。从气囊内压变化、流体对气囊的作用合力、气囊内压与入水速度的关系等方面研究了流体与气囊的相互作用,发现入水过程中气囊内压的变化主要受入水深度、运动速度、连接绳拉力等因素影响。通过计算不同初始内压条件下弹体的入水深度、减速时间及连接绳的拉力峰值,发现囊压越高,入水深度越小,减速时间越短,但是相应连接绳对弹体外壳的拉力峰值越大。因此,在进行入水回收气囊参数设计时,需要综合考虑缓冲效果、减速效果及气囊安全性等因素。
  • 图 1  气囊与弹体模型

    Figure 1.  Airbag and projectile model

    图 2  计算模型示意

    Figure 2.  Schematic diagram of calculation model

    图 3  有限元模型

    Figure 3.  Finite element model

    图 4  流体压力静平衡状态

    Figure 4.  Static equilibrium state of fluid pressure

    图 5  垂直入射全过程

    Figure 5.  Processes of vertical incidence

    图 6  弹体头部节点的竖向速度时程曲线

    Figure 6.  Vertical velocity history curve of projectile's head node

    图 7  弹体头部节点的竖向位移时程曲线

    Figure 7.  Vertical displacement history curve of projectile's head node

    图 8  斜入射全过程

    Figure 8.  Process of oblique incidence

    图 9  弹体头部节点的竖向和水平速度时程曲线

    Figure 9.  Vertical and horizontal velocity history curves of projectile's head node

    图 10  弹体头部节点的竖向和水平位移时程曲线

    Figure 10.  Vertical and horizontal displacement history curves of projectile's head node

    图 11  弹体角速度时程曲线

    Figure 11.  Angular velocity history curve of projectile

    图 12  弹体角位移时程曲线

    Figure 12.  Angular displacement curve of projectile

    图 13  气囊内压变化曲线

    Figure 13.  Internal pressure change curve of airbag

    图 14  流体对气囊的作用力合力时程曲线

    Figure 14.  Time history curve of fluid force on airbag

    图 15  囊压峰值与初始速度的关系曲线

    Figure 15.  Relationship between peak value of airbag internal pressure and initial velocity

    图 16  入水深度与初始囊压的关系曲线

    Figure 16.  Relationship between water entry depth and initial airbag internal pressure

    图 17  减速时间与初始囊压的关系曲线

    Figure 17.  Relationship between deceleration time and initial airbag internal pressure

    图 18  连接绳拉力峰值与初始囊压的关系曲线

    Figure 18.  Relationship between peak tension of connecting rope and initial airbag internal pressure

    表 1  气囊和连接绳的材料参数

    Table 1.  Material parameters of the airbag and corresponding rope

    材料 ρ/(kg·m-3) μ E/GPa
    气囊 875 0.2 0.557
    连接绳 840 0.2 21.9
    下载: 导出CSV

    表 2  水和空气参数

    Table 2.  Parameters of water and air

    材料 ρ/(kg·m-3) pc/Pa ν/(mPa·s) C/(m·s-1) S1 S2 S3 γ0
    998 -10 000 870 1 480 2.56 -1.99 0.227 0.5
    空气 1.185 -10 0.018 4 340 0 0 0 1.4
    下载: 导出CSV
  • [1] 陈帅, 李斌, 温金鹏, 等.软着陆气囊缓冲特性与参数设置的理论研究[J].振动与冲击, 2009, 28(4):25-28. doi: 10.3969/j.issn.1000-3835.2009.04.006
    CHEN Shuai, LI Bin, WEN Jinpeng, et al. Cushioning characteristic and parameter design of a soft landing airbag[J]. Journal of Vibration and Shock, 2009, 28(4):25-28. doi: 10.3969/j.issn.1000-3835.2009.04.006
    [2] 温金鹏, 李斌, 谭德伟, 等.考虑织布弹性的软着陆气囊缓冲特性研究[J].振动与冲击, 2010, 29(2):79-83. doi: 10.3969/j.issn.1000-3835.2010.02.018
    WEN Jinpeng, LI Bin, TAN Dewei, et al. Cushioning characteristics of a soft landing airbag with elastic fabric[J]. Journal of Vibration and Shock, 2010, 29(2):79-83. doi: 10.3969/j.issn.1000-3835.2010.02.018
    [3] 卫剑征, 谭惠丰, 万志敏, 等.缓冲气囊展开与缓冲着陆过程的仿真分析[J].航天返回与遥感, 2010, 31(5):1-8. doi: 10.3969/j.issn.1009-8518.2010.05.001
    WEI Jianzheng, TAN Huifeng, WAN Zhimin, et al. Simulation for airbag deployment and landing process of inflatable landing vechiles[J]. Spacecraft Recovery & Remote Sensing, 2010, 31(5):1-8. doi: 10.3969/j.issn.1009-8518.2010.05.001
    [4] 邓春燕, 裴锦华.全向式着陆缓冲气囊的折叠建模与充气过程仿真[J].南京航空航天大学学报, 2009, 41(增刊1):97-102.
    DENG Chunyan, PEI Jinhua. Modeling and process simulating about folding and inflating of omni-directional-type airbag[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2009, 41(Suppl 1):97-102.
    [5] 蔡文, 李斌, 温金鹏, 等.无人机软着陆气囊缓冲特性研究[J].兵工学报, 2014, 35(11):1867-1875. doi: 10.3969/j.issn.1000-1093.2014.11.019
    CAI Wen, LI Bin, WEN Jinpeng, et al. Research on cushioning characteristics of UAV soft landing airbags[J]. Acta Armamentarii, 2014, 35(11):1867-1875. doi: 10.3969/j.issn.1000-1093.2014.11.019
    [6] 邵志建, 裴锦华.某无人机横向圆柱排气式气囊着陆装置缓冲过程研究[J].航天返回与遥感, 2016, 37(2):26-33.
    SHAO Zhijian, PEI Jinghua. Simulation of bi-cylindrical airbag cushioning system for pilotless aircraft[J]. Spacecraft Recovery & Remote Sensing, 2008, 20(5):1325-1327.
    [7] 温金鹏, 李斌, 杨智春.缓冲气囊冲击减缓研究进展[J].宇航学报, 2010, 31(11):2438-2447. doi: 10.3873/j.issn.1000-1328.2010.11.002
    WEN Jinpeng, LI Bin, YANG Zhichun. Progress of study on impact attenuation capability of airbag cushion system[J]. Journal of Astronautics, 2010, 31(11):2438-2447. doi: 10.3873/j.issn.1000-1328.2010.11.002
    [8] 秦洪德, 赵林岳, 申静.入水冲击问题综述[J].哈尔滨工业大学学报, 2011(增刊1):152-157.
    QIN Hongde, ZHAO Linyue, SHEN Jing. Review of water entry problem[J]. Journal of Harbin Institute of Technology, 2011(Suppl 1):152-157.
    [9] WORTHINGTON A M, COLE R S. Impact with a liquid surface studied with aid of instantaneous photography[J]. Philosophical Transactions of the Royal Society of London, 1900, 194:175-199. doi: 10.1098/rsta.1900.0016
    [10] 王永虎, 石秀华.入水冲击问题研究的现状与进展[J].爆炸与冲击, 2008, 28(3):276-282. doi: 10.3321/j.issn:1001-1455.2008.03.014
    WANG Yonghu, SHI Xiuhua. Review on research and development of water-entry impact problem[J]. Explosion and Shock Waves, 2008, 28(3):276-282. doi: 10.3321/j.issn:1001-1455.2008.03.014
    [11] 卢炽华, 何友声.二维弹性结构入水冲击过程中的流固耦合效应[J].力学学报, 2000, 32(2):129-140. doi: 10.3321/j.issn:0459-1879.2000.02.001
    LU Chihua, HE Yousheng. Coupled analysis of nonlinear interaction between fluid and structure during impact[J]. Acta Mechanica Sinica, 2000, 32(2):129-140. doi: 10.3321/j.issn:0459-1879.2000.02.001
    [12] 李飞, 孙凌玉, 张广越, 等.圆柱壳结构入水过程的流固耦合仿真与试验[J].北京航空航天大学学报, 2007, 33(9):1117-1120. doi: 10.3969/j.issn.1001-5965.2007.09.027
    LI Fei, SUN Lingyu, ZHANG Guangyue, et al. Simulation and experiment of cylinder shell structure dropping into water based on fluid structure interaction[J]. Journal of Beijing University of Aeronautics and Astronautics, 2007, 33(9):1117-1120. doi: 10.3969/j.issn.1001-5965.2007.09.027
    [13] 施红辉, 胡青青, 陈波, 等.钝体倾斜和垂直冲击入水时引起的超空泡流动特性实验研究[J].爆炸与冲击, 2015, 35(5):617-624.
    SHI Honghui, HU Qingqing, CHEN Bo, et al. Experimental study of supercavitating flows induced by oblique and vertical water entry of blunt bodies[J]. Explosion and Shock Waves, 2015, 35(5):617-624.
    [14] 程涵.气囊工作过程仿真研究[D].南京: 南京航空航天大学, 2009.
    CHENG Han. Numerical simulation research on airbag working process[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2009.
    [15] 李裕春, 时党勇, 赵远.ANSYS 11.0/LS-DYNA基础理论与工程实践[M].北京:中国水利水电出版社, 2008.
    [16] 蒋克强, 张遵鸥, 张培成.基于ANSYS/LS-DYNA的示位标入水冲击仿真分析[J].电子机械工程, 2012, 28(2):15-17, 21. doi: 10.3969/j.issn.1008-5300.2012.02.005
    JIANG Keqiang, ZHANG Zun'ou, ZHANG Peicheng. Water-entry impact simulation analysis of a position indicating radio beacon based on ANSYS/LS-DYNA[J]. Electro-Mechanical Engineering, 2012, 28(2):15-17, 21. doi: 10.3969/j.issn.1008-5300.2012.02.005
    [17] 刘树红, 吴玉林, 左志钢.应用流体力学[M].北京:清华大学出版社, 2012.
    [18] 郭子涛, 张伟, 郭钊, 等.截卵形弹水平入水的速度衰减及空泡扩展特性[J].爆炸与冲击, 2017, 37(4):727-733.
    GUO Zitao, ZHANG Wei, GUO Zhao, 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.
    [19] 何春涛, 王聪, 魏英杰, 等.圆柱体垂直入水空泡形态试验[J].北京航空航天大学学报, 2012, 38(11):1542-1546.
    HE Chuntao, WANG Cong, WEI Yingjie, et al. Vertical water entry cavity of cylinder body[J]. Journal of Beijing University of Aeronautics and Astronautics, 2012, 38(11):1542-1546.
    [20] 邓春燕, 裴锦华.全向式气囊着陆装置缓冲过程的仿真研究[J].中国空间科学技术, 2010, 30(1):78-83.
    DENG Chunyan, PEI Jinhua. Simulation about buffer process of omni-directional-type airbag landing device[J]. Chinese Space Science and Technology, 2010, 30(1):78-83.
    [21] 李名琦.应急气囊着水冲击特性的试验研究与数值分析[D].南京: 南京航空航天大学, 2008.
    [22] 戴华杰, 胡振东, 咸奎成, 等.火星探测器气囊缓冲系统着陆过程仿真[J].力学季刊, 2010, 31(4):555-561.
    DAI Huajie, HU Zhendong, XIAN Kuicheng, et al. Simulation analysis for landing process of a Mars detector with airbag buffer system[J]. Chinese Quarterly of Mechanics, 2010, 31(4):555-561.
  • [1] 宋歌龙源钟明寿王敏吴建宇 . 离心机和密闭加压容器中水下爆炸相似关系. 爆炸与冲击, 2019, 39(2): 024102-1-024102-9. doi: 10.11883/bzycj-2017-0321
    [2] 徐浩铭顾文彬胡亚锋王振雄陈江海 . 串联EFP装药隔爆结构和延时起爆控制. 爆炸与冲击, 2014, 34(6): 723-729. doi: 10.11883/1001-1455(2014)06-0723-07
    [3] 龚曙光饶刚伍贤洪 . 基于无网格Galerkin法的穿甲侵彻数值模拟. 爆炸与冲击, 2011, 31(6): 658-663. doi: 10.11883/1001-1455(2011)06-0658-06
    [4] 李利莎杜建国张洪海谢清粮 . 爆炸冲击震动对砖墙破坏作用的数值模拟. 爆炸与冲击, 2015, 35(4): 459-466. doi: 10.11883/1001-1455(2015)04-0459-08
    [5] 侯日立周平彭建祥 . 冲击波作用下LY12铝合金结构毁伤的数值模拟. 爆炸与冲击, 2012, 32(5): 470-474. doi: 10.11883/1001-1455(2012)05-0470-05
    [6] 李勇肖伟程远胜刘均张攀 . 冲击波与破片对波纹杂交夹层板的联合毁伤数值研究. 爆炸与冲击, 2018, 38(2): 279-288. doi: 10.11883/bzycj-2016-0224
    [7] 廖华林李根生 . 孔隙流体耦合效应对射流冲击应力分布的影响分析. 爆炸与冲击, 2006, 26(1): 84-90. doi: 10.11883/1001-1455(2006)01-0084-07
    [8] 刘云龙汪玉张阿漫 . 基于二阶双渐近法的双层圆柱壳在水下爆炸作用下的鞭状运动. 爆炸与冲击, 2014, 34(6): 691-700. doi: 10.11883/1001-1455(2014)06-0691-10
    [9] 王永虎魏照宇 . 楔形体入水时域解的复边界元数值分析. 爆炸与冲击, 2012, 32(1): 55-60. doi: 10.11883/1001-1455(2012)01-0055-06
    [10] 王健赵庆彬陶钢吴军基 . 火箭橇水刹车高速入水冲击数值模拟. 爆炸与冲击, 2010, 30(6): 628-632. doi: 10.11883/1001-1455(2010)06-0628-05
    [11] 郭君杨文山姚熊亮张阿漫任少飞 . 基于场分离的水下爆炸流固耦合计算方法. 爆炸与冲击, 2011, 31(3): 295-299. doi: 10.11883/1001-1455(2011)03-0295-05
    [12] 黄志刚孙铁志杨碧野张桂勇宗智 . 平头锥型回转体高速入水结构强度数值分析. 爆炸与冲击, 2019, 39(4): 043201-1-043201-8. doi: 10.11883/bzycj-2017-0330
    [13] 周杰陶钢王健 . 爆炸冲击波对肺损伤的数值模拟. 爆炸与冲击, 2012, 32(4): 418-422. doi: 10.11883/1001-1455(2012)04-0418-05
    [14] 曹源金先龙李政 . 冲击载荷下柔性储液罐动态响应数值模拟及规律分析. 爆炸与冲击, 2011, 31(5): 469-474. doi: 10.11883/1001-1455(2011)05-0469-06
    [15] 周杰陶钢潘保青张洪伟 . 爆炸冲击波对人体胸部创伤机理的有限元方法研究. 爆炸与冲击, 2013, 33(3): 315-321. doi: 10.11883/1001-1455(2013)03-0315-06
    [16] 吴鹤翔刘颖 . 梯度变化对密度梯度蜂窝材料力学性能的影响. 爆炸与冲击, 2013, 33(2): 163-168. doi: 10.11883/1001-1455(2013)02-0163-06
    [17] 田玉滨李朝张春巍 . 爆炸荷载作用下配筋砌体结构的动力响应. 爆炸与冲击, 2012, 32(6): 658-662. doi: 10.11883/1001-1455(2012)06-0658-05
    [18] 赖鸣冯顺山黄广炎边江楠 . 不同加筋结构在水中接触爆炸下的破损规律. 爆炸与冲击, 2012, 32(6): 599-604. doi: 10.11883/1001-1455(2012)06-0599-05
    [19] 陈少辉李治源雷彬吕庆敖 . 大气/钢靶板交界处对平行入射聚能射流干扰的数值模拟. 爆炸与冲击, 2011, 31(6): 630-634. doi: 10.11883/1001-1455(2011)06-0630-05
    [20] 蒋宝权李玉龙刘元镛虞青俊 . SiC增强颗粒分布规律对梯度装甲板抗侵彻过程的影响. 爆炸与冲击, 2005, 25(6): 493-498. doi: 10.11883/1001-1455(2005)06-0493-06
  • 加载中
图(18)表(2)
计量
  • 文章访问数:  199
  • HTML全文浏览量:  47
  • PDF下载量:  70
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-10-26
  • 录用日期:  2018-01-02
  • 刊出日期:  2018-09-25

带环形密闭气囊弹体入水冲击过程的数值分析

    作者简介:陈洋(1994-), 男, 硕士研究生
    通讯作者: 吴亮, wuliangwust@sina.com
  • 1. 武汉科技大学理学院工程力学系, 湖北 武汉 430065
  • 2. 中铁港航-武汉科技大学爆破技术研究中心, 湖北 武汉 430065
基金项目:  国家自然科学基金项目 51479147湖北省自然科学基金项目 2014CFB822国家自然科学基金项目 11602178国家自然科学基金项目 51004079

摘要: 针对带环形密闭气囊弹体入水冲击问题,基于LS-DYNA,运用控制体积法模拟环形密闭气囊,结合流固耦合算法,模拟了某弹体及附带环形密闭气囊入水过程。将入水过程分为弹体砰水、气囊着水、入水减速、水中悬停、缓慢上浮、上浮出水、水面漂浮7个主要阶段,对比分析了垂直与倾斜入水过程中不同阶段弹体和气囊的姿态变化、减速特性及入水深度等特征的异同。从气囊内压变化、流体对气囊的作用合力、气囊内压与入水速度的关系等方面研究了流体与气囊的相互作用,发现入水过程中气囊内压的变化主要受入水深度、运动速度、连接绳拉力等因素影响。通过计算不同初始内压条件下弹体的入水深度、减速时间及连接绳的拉力峰值,发现囊压越高,入水深度越小,减速时间越短,但是相应连接绳对弹体外壳的拉力峰值越大。因此,在进行入水回收气囊参数设计时,需要综合考虑缓冲效果、减速效果及气囊安全性等因素。

English Abstract

  • 缓冲气囊在现代航空航天软着陆技术、重装空投、无人机回收等领域中有着广泛的应用,这些领域的研究成果是缓冲气囊参数设计的重要依据。陈帅等[1]基于热力学和刚体运动学方程,建立了软着陆气囊缓冲过程的理论模型,并结合数值方法分析得出了气囊参数设计的理论依据;温金鹏等[2]在考虑气囊织布弹性势能的基础上,结合能量守恒与热力学基本方程建立了缓冲气囊的简化理论模型,研究发现织物弹性模量、初始充气压力、排气口大小等因素对气囊的缓冲效果均有显著的影响;卫剑征等[3]运用显式有限元软件模拟了气囊展开和缓冲过程,验证了数值模拟方法的准确性。由此可知,关于各类缓冲气囊软着陆过程中冲击特性的研究一直是研究热点问题[4-7]

    入水冲击问题是固、液、气三者相互耦合作用的一种复杂的物理过程[8]。1900年,Worthington[9]利用闪光摄影技术对小球落入不同液体时的飞溅和空泡现象进行了实验观测。最初,研究者主要通过实验研究入水冲击过程,了解入水冲击荷载的规律,随后又逐渐提出了小斜升角模型的近似平板理论、自相似解法等。20世纪中期,各国在航空航天及军事领域开展了大量有关入水冲击的研究,该问题日益受到研究者的关注。

    随着数值方法和计算机技术的发展,复杂三维几何结构入水问题也逐渐得到有效解决[10]。近年来,我国学者也对各类结构冲击入水问题中进行了大量的研究工作[11-13]。虽然有关各类结构入水冲击问题的研究成果已经比较多,但是目前针对气囊入水冲击问题的研究仍有欠缺。由于气囊结构的特殊性,气囊入水冲击过程比气囊着陆缓冲过程或者一般结构的入水冲击过程都要复杂的多。例如,在水面对某飞行器进行无动力回收时,在到达水面之前,缓冲气囊迅速充气弹出,飞行器和气囊接触水面缓冲减速,入水之后在浮力的作用下气囊带着飞行器上浮,最终稳定漂浮在水面。以上过程中涉及到囊内气体、气囊壁、囊外液体三者之间复杂的相互作用,对该问题的研究在航空航天、救生、船舶等领域均有重要的实际意义。有鉴于此,本文中基于LS-DYNA,运用控制体积(control volume, CV)法模拟环形密闭气囊,结合流固耦合算法,模拟某飞行器回收入水姿态及减速上浮过程,研究气囊的入水冲击特性并分析囊内气压、气囊壁、囊外液体压力三者之间的相互作用机理,以期为入水回收气囊参数的优化设计提供依据。

    • 基于理想气体均压模型,CV法认为气囊内压由理想气体状态方程决定,囊内各处压力相等[14]。CV法不拘泥于气囊内部流场细节,在气囊壁围成的气囊体积内,通过质量流量和温度两个与时间相关的参数描述气流变化。当不需要考虑气囊充气过程外形和流场变化时,CV法是一种简单高效的方法。例如本文中假设,环形密闭气囊内部为理想气体,且其热容恒定;气囊入水过程绝热,内部温度与压强均匀。根据以上假设,气囊的控制方程可以写为:

      $ \left\{ \begin{gathered} pV = nRT \hfill \\ p = \left( {k-1} \right)U/V \hfill \\ U = n{C_V}T \hfill \\ \end{gathered} \right. $

      式中:p为气囊内压;V为气囊体积;n为囊内气体物质的量;R为理想气体常数;T为温度;k=Cp/CV为气体绝热指数,其中CpCV分别为定压和定容热容量;U为气囊内气体的内能。

    • 采用LS-DYNA进行显式动力分析主要有Lagrange、Euler和ALE 3种算法[15]。其中Lagrange方法主要应用于固体的结构分析,这种方法描述的网格与结构是一体的,单元节点即为物质点,网格的变化与结构变形完全一致。Lagrange方法主要优点是能够精确描述结构边界的运动,但是在处理大变形问题时,会由于单元的严重畸变而导致计算终止。Euler方法以空间坐标为基础,空间网格与物质相分离,单元节点为空间点,物质可以在网格之间流动。这种方法的优点在于网格大小和空间位置不变,计算中具有恒定的计算精度,但难以准确描述物质边界,多用于流体分析中。ALE兼具以上两种方法的优点,首先通过Lagrange方法处理物质边界的运动变形,然后执行ALE时步计算,过程如下:(1)保持变形后的物体边界条件,重新划分内部网格,网格的拓扑关系保持不变;(2)将变形网格中的单元变量(密度、能量、应力张量等)和节点速度矢量输运到重分后的新网格中。目前,LS-DYNA主要是通过多物质ALE算法实现流体与结构的耦合,具体参考文献[15],在此不加赘述。

    • 某实验弹上安装有回收气囊,在水面对其进行无动力回收,当弹体接近水面时,气囊迅速充气弹出,保证弹体落水后能够稳定并长时间漂浮,避免落水后击中水下装备,保障实验弹可靠回收。展开状态下的气囊与弹体模型如图 1所示。弹体呈圆柱状,环形气囊展开后体积大约是2.8 m3,气囊壁厚大约是0.5 mm,气囊与弹体之间通过4根连接绳相连,连接绳直径为1.5 cm。模拟计算中将弹体设为刚体,其质量为2 523 kg,当以环形气囊中心为原点,以弹体轴线为z轴时,弹体相对于xyz轴的转动惯量JxcJycJzc分别为3 011、3 011、225 kg·m2。气囊壁和连接绳的材料参数见表 1,其中ρ为密度,μ为泊松比,E为弹性模量。

      图  1  气囊与弹体模型

      Figure 1.  Airbag and projectile model

      材料 ρ/(kg·m-3) μ E/GPa
      气囊 875 0.2 0.557
      连接绳 840 0.2 21.9

      表 1  气囊和连接绳的材料参数

      Table 1.  Material parameters of the airbag and corresponding rope

      气囊展开状态下弹体及气囊入水冲击过程,涉及固、液、气三态耦合,是一类非常复杂的数值计算问题。LS-DYNA作为世界上著名的显式非线性动力分析软件,在求解各类非线性结构高速碰撞、爆炸和金属成型等非线性动力冲击问题有极大的优势,同时具有求解热传导、流体及流固耦合问题的能力,能够有效模拟真实物理世界的各种复杂问题,在汽车安全设计、武器系统设计、金属成形、跌落仿真等工程领域得到了广泛应用[15]。用LS-DYNA模拟气囊入水冲击问题是非常合适的。

    • 图 2所示,弹体和气囊在重力作用下,从水域上方以一定初速度坠入水中,计算中涉及到弹体、气囊、连接绳、水和空气的相互作用。运用LS-DYNA进行显式动力学分析时,流体采用ALE算法,固体采用Lagrange算法,固体和流体之间的相互作用通过流固耦合关键字定义,气囊与弹体之间通过自动单面接触关键字定义接触,连接绳与气囊之间通过生成节点刚性体的方式连接,与弹体之间采用共节点方式连接。有限元建模时,兼顾计算效率和精度,液体和空气采用solid 164单元划分成六面体网格,网格尺寸为0.2 m,气囊材料简化为各向同性的线弹性无弯矩薄膜材料,采用四边形薄膜单元划分网格,网格尺寸为0.1 m。连接绳采用柔性索单元划分网格,网格尺寸为0.1 m。弹体简化成刚性体,采用shell单元划分弹体外壳网格,并定义弹体质量特性。为模拟无限水域环境,取10 m×10 m×10 m水域范围,约束底部竖向自由度,其他水域边界及空气边界均设为无反射边界条件。最终生成的有限元模型如图 3所示,一共划分了166 220个单元。

      图  2  计算模型示意

      Figure 2.  Schematic diagram of calculation model

      图  3  有限元模型

      Figure 3.  Finite element model

    • 水和空气均采用MAT_NULL材料模型,通过Grüneisen状态方程描述,其压力为:

      $ p = \frac{{{\rho _0}{C^2}\mu \left[{1 + \left( {1-{\gamma _0}/2} \right)\mu-a{\mu ^2}/2} \right]}}{{{{\left[{1-\left( {{S_1}-1} \right)\mu-{S_2}\frac{{{\mu ^2}}}{{\mu + 1}} - {S_3}\frac{{{\mu ^3}}}{{{{\left( {\mu + 1} \right)}^2}}}} \right]}^2}}} + \left( {{\gamma _0} + a\mu } \right)E $

      式中:p为压力;C为声速;S1S2S3为系数;γ0E分别是Grüneisen常数和初始内能,均取为默认值零;a未来一阶体积修正系数;μ= ρ/ρ0-1,ρρ0分别为当前密度和初始密度。水和空气参数设置如表 2所示[16],其中pc为截断压力,一般假定为负值;ν为材料动力黏度。

      材料 ρ/(kg·m-3) pc/Pa ν/(mPa·s) C/(m·s-1) S1 S2 S3 γ0
      998 -10 000 870 1 480 2.56 -1.99 0.227 0.5
      空气 1.185 -10 0.018 4 340 0 0 0 1.4

      表 2  水和空气参数

      Table 2.  Parameters of water and air

    • 气囊入水过程中,随着深度增加,静水压力呈线性增加。为准确模拟自然环境中在重力作用下水域静压强环境(本文中水域压强、气囊内压等均以相对压强进行分析,即超出大气压的超压值),在弹体和气囊入水之前,对水域压力进行静平衡处理。在LS-DYNA中用LOAD_BODY关键字模拟重力,结合INITIAL_ALE_HYDROSTATIC关键字对静水压强进行初始化。如图 4所示,流体压强在0.05 s以内即达到平衡状态,平衡状态符合静水压强公式[17]

      $ {p_{\text{w}}} = \rho gh $

      图  4  流体压力静平衡状态

      Figure 4.  Static equilibrium state of fluid pressure

    • 气囊和弹体的入水过程主要分为弹体砰水、气囊着水、入水减速、水中悬停、缓慢上浮、上浮出水、水面漂浮7个阶段。气囊的缓冲作用主要体现在气囊着水阶段;入水减速阶段主要受流体黏滞阻力、浮力和重力影响;当弹体到达最低点附近,在一段时间内将处于悬停状态(此阶段速度非常小,可以认为是静止状态);随后,浮力作用使弹体和气囊开始缓慢上浮,并最终稳定漂浮在水面。各阶段的特征受初始条件影响较大,取气囊初始囊压为50 kPa,弹体初始速度为50 m/s,数值模拟弹体从距离水面3 m处垂直入射或从距离水面垂直高度为2.14 m处斜入射(弹体轴线与水面夹角为45°)的入水过程。

    • 图 5所示,在0.06 s时刻,弹体头部接触水面,此时为弹体砰水阶段。0.09 s时,气囊接触水面,在冲击挤压及水压作用下,气囊迅速压缩,起到减速与缓冲作用。由于入水速度较快,弹体和气囊周围形成较大范围的空泡现象[13, 18-19]。入水之后,在浮力和流体黏滞阻力作用下,气囊和弹体持续减速,直到0.81 s,速度减为零,气囊和弹体在水中基本处于悬停状态;此时周围的空泡逐渐闭合,水面溅起较高的水花。从2.81 s开始,气囊和弹体缓慢加速上浮,上浮最大速度小于2 m/s,最终在6.72 s时浮出并漂浮在水面,整个过程弹体始终保持直立姿态。图 5中弹体周围空泡的生成、发展和闭合过程及形态与文献[18]和文献[19]中弹体垂直入水的实验照片基本一致,但图 5中弹体带有气囊,形成的空泡范围更大。图 6图 7分别为入水过程中弹体头部节点的速度和位移时程曲线,反映出各阶段的速度和位移变化特征。由图 6~7可知:弹体从入水到浮出水面历时约为6.7 s,最大入水深度为6.47 m(垂直入射时,零时刻弹体距离水面高度为3 m)。

      图  5  垂直入射全过程

      Figure 5.  Processes of vertical incidence

      图  6  弹体头部节点的竖向速度时程曲线

      Figure 6.  Vertical velocity history curve of projectile's head node

      图  7  弹体头部节点的竖向位移时程曲线

      Figure 7.  Vertical displacement history curve of projectile's head node

    • 图 5对应,斜入射时弹体和气囊入水的全过程如图 8所示。0.06 s弹体砰水;0.09 s气囊着水;0.14 s气囊完全入水;随后一直到0.56 s为入水减速阶段;0.56 ~1.25 s为水中悬停阶段,此时弹体头部速度极小,主要是尾部在运动,弹体由倾斜转变成直立姿态;1.25 s时刻弹体和气囊开始上浮;2.25 s完全浮出水面。

      图  8  斜入射全过程

      Figure 8.  Process of oblique incidence

      斜入射时弹体在水中的姿态与垂直入射有较大差异,以弹体头部节点为基点,图 9图 10分别是斜入水过程中弹体的平动速度和平动位移时程曲线,图 11图 12分别是弹体的角速度和角位移(顺时针为正,逆时针为负)时程曲线。入水后弹体角速度迅速增大到正向峰值,弹体轴线与水平方向的夹角逐渐减小,这是由于入水减速后弹体头部和尾部速度不一致导致的:头部所受阻力较大,并且受到气囊和连接绳的拉力,所以竖向速度衰减较快,而空泡主要靠近尾部,所以尾部所受阻力较小。这样,入水之后弹体头部竖向速度将小于尾部,弹体将发生顺时针旋转。直到到达最低点,气囊开始上浮,在连接绳拉力作用下,弹体角速度逐渐由正转为负,弹体调整到直立状态,并上浮出水,上浮速度小于1 m/s,最终弹体和气囊一起在水面漂浮。从入水到浮出水面历时约2.3 s,最大入水深度为3.67 m(斜入射时,零时刻弹体距离水面高度垂直高度为2.14 m)。

      图  9  弹体头部节点的竖向和水平速度时程曲线

      Figure 9.  Vertical and horizontal velocity history curves of projectile's head node

      图  10  弹体头部节点的竖向和水平位移时程曲线

      Figure 10.  Vertical and horizontal displacement history curves of projectile's head node

      图  11  弹体角速度时程曲线

      Figure 11.  Angular velocity history curve of projectile

      图  12  弹体角位移时程曲线

      Figure 12.  Angular displacement curve of projectile

    • 在着陆过程中,气囊着陆缓冲系统会发生多次弹跳,受摩擦力、气囊结构阻尼和内流阻尼的影响,弹跳速度会逐步降低,着陆动能在与地面碰撞过程中逐步衰减[20]。因此,着陆过程中气囊内压会出现多次波动峰值,并且峰值大小逐步衰减。气囊入水过程与着陆过程有明显差异:(1)气囊着水阶段是气囊与流体碰撞的流固耦合冲击问题,该过程中流体本身对碰撞有一定的缓冲作用,和气囊与地面碰撞不同,不会出现弹跳现象;(2)流体的黏滞阻力与速度v呈正比,随着弹体和气囊在水中运动速度的减小,阻力逐渐减小;(3)随着入水深度h的增加,流体对气囊的压力逐渐增大;(4)弹体和气囊之间连接绳拉力T也是影响囊压的重要因素。因此,在入水后到上浮过程中,气囊内压可用以下函数表示:

      $ p\left( t \right) = {p_0} + p\left( {h, v, {T_{{\text{rope}}}}, \rho, \mu } \right) $

      式中:p0为气囊初始内压,hvTrope均是与时间相关的函数,分别表示入水深度、运动速度、连接绳拉力,ρμ分别为流体密度和黏滞系数。

      图 13为气囊内压的变化曲线。初始囊压为50 kPa,弹体以50 m/s的初始速度从水面上方入射,分为垂直入射和45°斜入射两种情况。在气囊着水阶段,气囊与水面猛烈碰撞,气囊急剧变形,囊压迅速升高到峰值,囊压变化曲线与文献[21]一致。随后,由于入水速度较大,气囊和弹体排水形成空泡,水对气囊的压力减小,气囊内压迅速减小。随着空泡的闭合,水猛烈撞击气囊,气囊内压再次迅速上升到峰值,然后随着运动速度的衰减而缓慢衰减。可以看到,垂直入射时,在1.5~3.0 s时间段内,气囊内压基本稳定在约68 kPa,此时间段对应于气囊和弹体到达最低点,处于悬停阶段。此后,随着气囊上浮,静水压力减小,气囊内压逐渐减小。无论是斜入射还是垂直入射,浮出水面后,漂浮状态下,气囊的内压都稳定在约52.5 kPa。

      图  13  气囊内压变化曲线

      Figure 13.  Internal pressure change curve of airbag

      图 14为流体对气囊的作用力合力的时程曲线,垂直入射时,合力主要沿z方向,另外两个方向基本为零;而斜入射主时合力主要沿z方向和x方向,且z方向合力略大于x方向。结合弹体的速度变化进行分析可以发现:流体合力峰值出现的区间正是气囊着水和水中减速阶段,并且与囊压峰值出现的时间段一致;当初始速度衰减到零之后,流体对气囊的作用合力主要是向上的浮力。

      图  14  流体对气囊的作用力合力时程曲线

      Figure 14.  Time history curve of fluid force on airbag

      由以上分析可知,入射速度对气囊压力峰值的大小有很大影响。设气囊初始内压为50 kPa,弹体以不同初速度垂直入射到水中,在此条件下计算得到气囊内压峰值与初始速度之间的关系曲线,如图 15所示。可以看出:初始速度越大,囊压峰值越高,即气囊受水的冲击力越大。

      图  15  囊压峰值与初始速度的关系曲线

      Figure 15.  Relationship between peak value of airbag internal pressure and initial velocity

    • 气囊内压是进行缓冲气囊优化设计的重要参数。在软着陆气囊设计中,囊压太低容易引起气囊变形过大,从而导致“硬着陆”而损坏结构;气压太高又会导致缓冲加速度上升,缓冲效果不佳[22]。在气囊入水问题中,随着入水深度的增加,气囊内压在流体压力作用下会升高。根据CV法原理,密闭气囊体积与内压呈反比,因此,若气囊初始囊压过低,入水压缩后气囊体积偏小,所产生的浮力和阻力不够,则会导致入水过深无法短时间内浮出水面,甚至沉入水底。

      图 16为弹体以50 m/s的速度垂直入水时入水深度与初始囊压之间的关系曲线。可以看出,入水深度随着囊压的降低而增大。充气完成后,气囊体积约为2.8 m3,弹体本身的排水体积约为1.9 m3,根据浮力计算公式[17],所能产生的浮力大于弹体及结构重力。但是当气囊囊压过低时,随着入水深度的增大,气囊体积持续被压缩,最终所能产生的浮力将不足以抵消重力,气囊将无法上浮甚至沉入水底。在压力为20~30 kPa区间内按照二分法选取囊压进行计算,可以确定气囊的临界压力大约为23.75 kPa,低于此压力,弹体将沉没。

      图  16  入水深度与初始囊压的关系曲线

      Figure 16.  Relationship between water entry depth and initial airbag internal pressure

      图 17为初始囊压与弹体速度减到零所耗时间的关系曲线。可以看出:初始囊压越高,减速时间越短,减速效果越好。但是囊压的设计值并不是越高越好。一方面,囊压过高可能导致回收过程中气囊破裂,影响整个弹体回收过程的安全性与可靠性;另一方面,气囊与弹体之间连接绳的拉力峰值随着囊压的升高而增大,即囊压越高,连接绳对弹体的冲击力越大,这样极有可能导致回收过程中弹体外壳被拉坏,影响回收效果,如图 18所示。因此,在进行入水回收气囊参数设计时,有必要根据实际情况综合考虑缓冲效果、减速效果以及气囊安全性等因素。

      图  17  减速时间与初始囊压的关系曲线

      Figure 17.  Relationship between deceleration time and initial airbag internal pressure

      图  18  连接绳拉力峰值与初始囊压的关系曲线

      Figure 18.  Relationship between peak tension of connecting rope and initial airbag internal pressure

    • 基于LS-DYNA,运用CV法模拟环形密闭气囊,结合流固耦合算法,模拟了某弹体及附带环形密闭气囊在展开状态下的入水过程,分析了多种工况下气囊和弹体入水姿态,探讨了影响气囊和弹体入水回收效果的多种因素,可得到以下结论。

      (1) 气囊和弹体的入水过程分为弹体砰水、气囊着水、入水减速、水中悬停、缓慢上浮、上浮出水、水面漂浮7个阶段。弹体垂直入水后能始终保持直立姿态;斜入水初期,弹体轴线与水平方向夹角会有所减小,但随着运动速度的衰减,弹体在悬停阶段开始逐渐调整姿态,并在上浮阶段中逐渐转变成直立姿态,浮出水面后弹体能保持直立姿态漂浮在水面。

      (2) 气囊入水过程中囊压的变化主要受入水深度、运动速度、连接绳拉力等因素影响。囊压峰值出现在气囊着水阶段,此时气囊与水面猛烈碰撞,水对气囊的冲击力较大。在上浮阶段随着入水深度减小,气囊内压也逐渐减小。另外,入水冲击过程中气囊内压峰值随着入水速度的增大而升高。

      (3) 弹体入水深度随着气囊初始内压的降低而增大。初始内压过低可能导致气囊减速时间过长,或气囊入水后产生的浮力不足以抵消弹体重力而沉入水底。气囊内压越高,入水冲击过程中连接绳对弹体的拉力峰值越大。当囊压过高时,连接绳可能会拉坏弹体外壳,从而影响回收效果。因此,在进行入水回收气囊参数设计时,需要综合考虑缓冲效果、减速效果及气囊安全性等因素。

参考文献 (22)

目录

    /

    返回文章
    返回