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

平头锥型回转体高速入水结构强度数值分析

黄志刚 孙铁志 杨碧野 张桂勇 宗智

引用本文:
Citation:

平头锥型回转体高速入水结构强度数值分析

    作者简介: 黄志刚(1992- ),男,硕士研究生,1059283140@qq.com;
    通讯作者: 张桂勇, gyzhang@dlut.edu.cn
  • 中图分类号: O352

Numerical analysis on structural strength of a cone-shaped flatted revolution body during high-speed water-entry

    Corresponding author: ZHANG Guiyong, gyzhang@dlut.edu.cn ;
  • CLC number: O352

  • 摘要: 为探究回转体在高速入水过程中的结构强度,基于非线性有限元LS-DYNA软件中流固耦合任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian, ALE)方法,分析了不同壁厚的回转体以100 m/s的初速度入水过程中的冲击力特性和结构强度。结果表明:数值计算得到的入水冲击压强峰值和速度衰减曲线与相应的理论值吻合较好,从而验证了数值方法的有效性;入水冲击载荷峰值出现在结构入水瞬间,结构入水后冲击载荷急剧变小且微小震荡;回转体的结构形式对其在高速入水过程中的结构强度有重要影响,尤其回转体头部厚度影响回转体结构强度,当回转体头部厚度为8 mm、后体壁厚大于2.5 mm时,可以保证回转体强度要求。
  • 图 1  计算域

    Figure 1.  Computational domain

    图 2  回转体几何模型

    Figure 2.  The geometrical model for a revolution body

    图 3  回转体头部中心区域压强曲线

    Figure 3.  Pressure intensity curve in the central head region of the revolution body

    图 4  回转体速度衰减曲线

    Figure 4.  Velocity attenuation of the revolution body over time

    图 5  流-构耦合力曲线

    Figure 5.  Fluid-structure interaction force varying with time

    图 6  t=1.1 ms时回转体的应力分布

    Figure 6.  Stress distribution in the revolution body at t=1.1 ms

    图 7  t=3.0 ms时回转体的应力分布

    Figure 7.  Stress distribution in the revolution body at t=3.0 ms

    图 8  回转体头部边缘单元应变曲线

    Figure 8.  Strain-time curve of the edge element for the head of the revolution body

    图 9  回转体头部边缘单元等效应力曲线

    Figure 9.  Equivalent stress-time curve of the edge element for the head of the revolution body

    图 10  t=0.89 ms时刻回转体的应力分布

    Figure 10.  Stress distribution in the revolution body at t=0.89 ms

    图 11  t=2.00 ms回转体应力图

    Figure 11.  Stress distribution in the revolution body at t=2.00 ms

    图 12  回转体头部边缘及中心单元应变随时间的变化

    Figure 12.  Strain-time curves of the edge and central elements for the head of the revolution body

    图 13  回转体头部边缘单元等效应力曲线

    Figure 13.  Equivalent stress-time curve of the edge element for the head of the revolution body

    图 14  回转体头部中心单元等效应力曲线

    Figure 14.  Equivalent stress-time curve of the central element for the head of the revolution body

  • [1] 王永虎, 石秀华. 入水冲击问题研究的现状与进展 [J]. 爆炸与冲击, 2008, 28(3): 276–282. DOI: 10.11883/1001-1455(2008)03-0276-07
    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.11883/1001-1455(2008)03-0276-07
    [2] WORTHINGTON A M, COLE R S. Impact with a liquid surface studied by the aid of instantaneous photography: paper II [J]. Philosophical Transactions of the Royal Society of London, 1900, 194: 175–199. DOI: 10.1098/rsta.1897.0005.
    [3] VON KARMAN T. The impact on seaplane floats during landing: 20140804-120956389 [R]. 1929.
    [4] WAGNER V H. Phenomena associated with impacts and sliding on liquid surfaces [J]. ZAMM: Journal of Applied Mathematics and Mechanics, 1932, 12(4): 193–215. doi: 10.1002/(ISSN)1521-4001
    [5] MAYO W L. Hydrodynamic impact of a system with a single elastic mode: I: theory and generalized solution with an application to an elastic airframe: NACA-report-1074 [R]. Technical Report Archive & Image Library, 1952.
    [6] MAY A. Review of water-entry theory and data [J]. Journal of Hydronautics, 1970, 4(4): 140–142. DOI: 10.2514/3.62851.
    [7] ANGHILERI M, SPIZZICA A. Experimental validation of finite element models for water impacts [C] // Proceedings of the Second International Crash Users Seminar. Cranfield, UK, 1995.
    [8] FALTINSEN O M, ZHAO R. Water entry of ship sections and axisymmetric bodies: AGARD report, 827 [R]. Neuilly-sur-Seine, WaterKiev, Ukraine: Advisory Group for Aerospace Research and Development, 1997: 24−1−24−11.
    [9] KOROBKIN A A, WU G X. Impact on a floating circular cylinder [J]. Proceedings: Mathematical, Physical and Engineering Sciences, 2000, 456(2002): 2489–2514. DOI: 10.1098/rspa.2000.0622.
    [10] DONGUY B, PESEUX B, GORNET L, et al. Three-dimensional hydro-elastic water entry: preliminary results [C] // The Eleventh International Offshore and Polar Engineering Conference. Stavanger, Norway: International Society of Offshore and Polar Engineers, 2001: 324−330.
    [11] 郑金伟, 宗智. 三维刚体椭圆头结构高速倾斜入水冲击模拟 [J]. 船海工程, 2012, 41(3): 7–9 doi: 10.3963/j.issn.1671-7953.2012.03.003
    ZHENG Jinwei, ZONG Zhi. 3-dimensional numerical simulation of rigid elliptic structure inclined water-entry at high speed [J]. Ship and Ocean Engineering, 2012, 41(3): 7–9 doi: 10.3963/j.issn.1671-7953.2012.03.003
    [12] 施红辉, 高见卓也, 伊藤基之. 钝体入水时的水下声场的测量 [J]. 流体力学实验与测量, 2001, 15(2): 78–84 doi: 10.3969/j.issn.1672-9897.2001.02.011
    SHI Honghui, TAKAMI Takuya, ITOH Motoyuki. Measurement of the underwater acoustic field in water entry of blunt body [J]. Experiments and Measurements in Fluid Mechanics, 2001, 15(2): 78–84 doi: 10.3969/j.issn.1672-9897.2001.02.011
    [13] SHI Honghui, KUME Makoto. Underwater acoustics and cavitating flow of water entry [J]. Acta Mechanica Sinica, 2004, 20(4): 374–382. DOI: 10.1007/bf02489375.
    [14] 王云, 袁绪龙, 吕策. 弹体高速入水弯曲弹道实验研究 [J]. 兵工学报, 2014, 35(12): 1998–2002 doi: 10.3969/j.issn.1000-1093.2014.12.010
    WANG Yun, YUAN Xulong, LYU Ce. Experimental research on curved trajectory of high-speed water-entry missile [J]. Acta Armamentarii, 2014, 35(12): 1998–2002 doi: 10.3969/j.issn.1000-1093.2014.12.010
    [15] 潘光, 杨悝. 空投鱼雷入水载荷 [J]. 爆炸与冲击, 2014, 34(5): 521–526. DOI: 10.11883/1001-1455(2014)05-0521-06
    PAN Guang, YANG Kui. Impact force encountered by water-entry airborne torpedo [J]. Explosion and Shock Waves, 2014, 34(5): 521–526. DOI: 10.11883/1001-1455(2014)05-0521-06
    [16] 黄凯, 乐述文. 不同头型弹体模型入水现象的实验研究 [J]. 物理实验, 2016, 36(5): 13–18 doi: 10.3969/j.issn.1005-4642.2016.05.003
    HUANG Kai, YUE Shuwen. Experimental research on the behavior of water-entry of different head shape projectile models [J]. Physics Experimentation, 2016, 36(5): 13–18 doi: 10.3969/j.issn.1005-4642.2016.05.003
    [17] 李佳川, 魏英杰, 王聪, 等. 不同扰动角速度高速射弹入水弹道特性 [J]. 哈尔滨工业大学学报, 2017, 49(4): 131–136
    LI Jiachuan, WEI Yingjie, WANG Cong, et al. Water entry trajectory characteristics of high-speed projectile with various turbulent angular velocity [J]. Journal of Harbin Institute of Technology, 2017, 49(4): 131–136
    [18] 张伟, 郭子涛, 肖新科, 等. 弹体高速入水特性实验研究 [J]. 爆炸与冲击, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06
    ZHANG Wei, GUO Zitao, XIAO Xinke, 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
    [19] 郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [J]. 爆炸与冲击, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07
    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. DOI: 10.11883/1001-1455(2017)04-0727-07
    [20] 马庆鹏, 何春涛, 王聪, 等. 球体垂直入水空泡实验研究 [J]. 爆炸与冲击, 2014, 34(2): 174–180. DOI: 10.11883/1001-1455(2014)02-0174-07
    MA Qingpeng, HE Chuntao, WANG Cong, et al. Experimental investigation on vertical water-entry of sphere [J]. Explosion and Shock Waves, 2014, 34(2): 174–180. DOI: 10.11883/1001-1455(2014)02-0174-07
    [21] 王珂, 王自力, 王志东, 等. 弹性回转体入水砰击载荷预报 [J]. 船海工程, 2011, 40(5): 20–22 doi: 10.3963/j.issn.1671-7953.2011.05.006
    WANG Ke, WANG Zili, WANG Zhidong, et al. Prediction of the slamming pressure on a 3-D elastic axisymmetric structure [J]. Ship and Ocean Engineering, 2011, 40(5): 20–22 doi: 10.3963/j.issn.1671-7953.2011.05.006
  • [1] 王健赵庆彬陶钢吴军基 . 火箭橇水刹车高速入水冲击数值模拟. 爆炸与冲击, 2010, 30(6): 628-632. doi: 10.11883/1001-1455(2010)06-0628-05
    [2] 王永虎石秀华 . 入水冲击问题研究的现状与进展. 爆炸与冲击, 2008, 28(3): 276-282. doi: 10.11883/1001-1455(2008)03-0276-07
    [3] 潘光杨悝 . 空投鱼雷入水载荷. 爆炸与冲击, 2014, 34(5): 521-526. doi: 10.11883/1001-1455(2014)05-0521-06
    [4] 郭君杨文山姚熊亮张阿漫任少飞 . 基于场分离的水下爆炸流固耦合计算方法. 爆炸与冲击, 2011, 31(3): 295-299. doi: 10.11883/1001-1455(2011)03-0295-05
    [5] 马庆鹏何春涛王聪魏英杰路中磊孙健 . 球体垂直入水空泡实验研究. 爆炸与冲击, 2014, 34(2): 174-180. doi: 10.11883/1001-1455(2014)02-0174-07
    [6] 陈洋吴亮曾国伟周俊汝 . 带环形密闭气囊弹体入水冲击过程的数值分析. 爆炸与冲击, 2018, 38(5): 1155-1164. doi: 10.11883/bzycj-2017-0387
    [7] 陈晓坤李海涛王秋红金永飞邓军张嬿妮 . 瓦斯爆炸载荷作用下矿用柱壳救生舱抗爆性分析及结构优化. 爆炸与冲击, 2018, 38(1): 124-132. doi: 10.11883/bzycj-2016-0248
    [8] 廖华林李根生 . 孔隙流体耦合效应对射流冲击应力分布的影响分析. 爆炸与冲击, 2006, 26(1): 84-90. doi: 10.11883/1001-1455(2006)01-0084-07
    [9] 刘云龙汪玉张阿漫 . 基于二阶双渐近法的双层圆柱壳在水下爆炸作用下的鞭状运动. 爆炸与冲击, 2014, 34(6): 691-700. doi: 10.11883/1001-1455(2014)06-0691-10
    [10] 周杰陶钢潘保青张洪伟 . 爆炸冲击波对人体胸部创伤机理的有限元方法研究. 爆炸与冲击, 2013, 33(3): 315-321. doi: 10.11883/1001-1455(2013)03-0315-06
    [11] 周杰陶钢王健 . 爆炸冲击波对肺损伤的数值模拟. 爆炸与冲击, 2012, 32(4): 418-422. doi: 10.11883/1001-1455(2012)04-0418-05
    [12] 曹源金先龙李政 . 冲击载荷下柔性储液罐动态响应数值模拟及规律分析. 爆炸与冲击, 2011, 31(5): 469-474. doi: 10.11883/1001-1455(2011)05-0469-06
    [13] 陈洋吴亮许锋鲁帅 . 爆破开挖振动下既有大型储油罐的动力响应. 爆炸与冲击, 2018, 38(6): 1394-1403. doi: 10.11883/bzycj-2017-0128
    [14] 康德严平 . 基于LS-DYNA的高速破片水中运动特性流固耦合数值模拟. 爆炸与冲击, 2014, 34(5): 534-538. doi: 10.11883/1001-1455(2014)05-0534-05
    [15] 张伟郭子涛肖新科王聪 . 弹体高速入水特性实验研究. 爆炸与冲击, 2011, 31(6): 579-584. doi: 10.11883/1001-1455(2011)06-0579-06
    [16] 高全臣陆华王东何广骥 . 多孔隙流固耦合砂岩的冲击损伤效应. 爆炸与冲击, 2012, 32(6): 629-634. doi: 10.11883/1001-1455(2012)06-0629-06
    [17] 郭攀武文华刘君吴志刚 . 激波作用下含缺陷固体火箭装药的流固耦合数值模拟. 爆炸与冲击, 2014, 34(1): 93-98. doi: 10.11883/1001-1455(2014)01-0093-06
    [18] 浦锡锋王仲琦白春华周刚 . 用于爆炸流场与结构间相互作用分析的Euler-Lagrange耦合模拟技术. 爆炸与冲击, 2011, 31(1): 6-10. doi: 10.11883/1001-1455(2011)01-0006-05
    [19] 陈先富 . 弹丸入水空穴的试验研究. 爆炸与冲击, 1985, 5(4): 70-73.
    [20] 顾建农张志宏范武杰 . 旋转弹丸入水侵彻规律. 爆炸与冲击, 2005, 25(4): 341-349. doi: 10.11883/1001-1455(2005)04-0341-09
  • 加载中
图(14)
计量
  • 文章访问数:  147
  • HTML全文浏览量:  414
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-09-08
  • 录用日期:  2018-01-10
  • 网络出版日期:  2019-03-25
  • 刊出日期:  2019-04-01

平头锥型回转体高速入水结构强度数值分析

    作者简介:黄志刚(1992- ),男,硕士研究生,1059283140@qq.com
    通讯作者: 张桂勇, gyzhang@dlut.edu.cn
  • 1. 大连理工大学船舶工程学院辽宁省深海浮动结构工程实验室,辽宁 大连 116024
  • 2. 大连理工大学工业装备结构分析国家重点实验室,辽宁 大连 116024
  • 3. 高新船舶与深海开发装备协同创新中心,上海 200240

摘要: 为探究回转体在高速入水过程中的结构强度,基于非线性有限元LS-DYNA软件中流固耦合任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian, ALE)方法,分析了不同壁厚的回转体以100 m/s的初速度入水过程中的冲击力特性和结构强度。结果表明:数值计算得到的入水冲击压强峰值和速度衰减曲线与相应的理论值吻合较好,从而验证了数值方法的有效性;入水冲击载荷峰值出现在结构入水瞬间,结构入水后冲击载荷急剧变小且微小震荡;回转体的结构形式对其在高速入水过程中的结构强度有重要影响,尤其回转体头部厚度影响回转体结构强度,当回转体头部厚度为8 mm、后体壁厚大于2.5 mm时,可以保证回转体强度要求。

English Abstract

  • 回转体从空中高速跨越空气和水交界面时,受到水介质对其作用的巨大冲击载荷,从而引起回转体发生形变和破坏,这就需要设计的回转体强度满足结构安全要求。同时由于回转体入水过程涉及气-液-固三相耦合作用[1],在高速条件下这种流固耦合作用变得更复杂,从而大大增加了研究高速回转体入水问题的难度。因此,开展回转体高速入水过程结构强度研究对回转体结构设计具有重要意义。

    入水问题研究长期以来一直受到广泛关注。最早Worthington等[2]利用闪光摄影技术,观测到了小球落水时的液体飞溅和空泡现象。在理论研究上,1929年,Von Karman[3]最早提出使用附加质量法计算了入水冲击载荷,并采用动量守恒定律推导了相关的计算公式。随后Wagner[4]在此基础上运用伯努利方程,提出了平板理论,此理论在研究结构入水过程受力和流体动能方面得到了广泛应用。Mayo[5]考虑了波浪对附加质量、压力分布、入水冲击载荷的影响。May[6]研究了导弹在垂直和倾斜入水过程中的入水空泡形成、发展和溃灭的过程,同时研究了入水过程的流场变化情况。近些年来,随着数值计算能力的提高,有限元将流体和结构耦合解决了大型三维几何入水问题。Anghileri等[7]利用有限元分析了刚性球在垂直入水过程中所受的载荷。Faltinsen等[8]建立了入水问题的数值和理论模型,并通过相关的入水实验验证其模型的准确性。Korobkin等[9]运用边界元和有限差分法计算悬浮体入水冲击过程中引起的液流非定常问题。Donguy等[10]对二维、三维、刚性和弹性体的入水耦合过程进行了详细分析,其采用的是一种变化的有限元耦合方法,通过耦合矩阵处理流体和结构之间的作用。

    郑金伟等[11]利用LS-DYNA程序计算了三维刚体结构倾斜入水过程的流体压力和冲击载荷。施红辉等[12]、Shi等[13]通过钝体入水实验测量了水下声场的变化,获得了激波的空间能量分布以及传播速度变化,同时其还研究了钝体入水自由液面变形问题和超空泡现象。王云等[14]开展了模型实验,利用高速摄像机拍摄了回转体入水过程的空泡形态演变,分析了头型、入水角度、入水速度和水下弹道的影响。潘光等[15]利用MSC. DYTRAN软件计算了入水过程中鱼雷所受到的冲击压力以及压力分布情况。黄凯等[16]通过实验研究了回转体头型对反弹水花的影响,发现回转体顶角越大,空泡越大。李佳川等[17]基于入水弹道学和超空泡理论,研究了不同扰动角速度下的回转体入水轨迹、速度、俯仰角的变化。张伟等[18] 、郭子涛等[19]开展了水平回转体入水实验,研究了回转体入水过程的弹道稳定性和空泡拓展特性,计算了入水阻力系数,并将实验中回转体的速度衰减曲线与理论解进行对比,两者吻合良好。马庆鹏等[20]通过球体入水实验,分析了球体入水过程中的入水空泡演变过程,得到了不同的入水速度以及表面黏湿状态对入水空泡流场的影响。

    目前针对入水问题已开展了大量的研究,并取得了丰硕的成果,然而在研究回转体速度达到100 m/s及以上的高速入水问题时,多以入水空泡演化以及回转体入水稳定性作为研究重点。由于高速入水问题是冲击作用时间极短的非线性问题,对高速回转体入水过程结构强度方面的研究仍显不足。本文中采用数值计算的方法,分析入水过程中的冲击载荷特点,并开展回转体在不同壳体厚度情况下高速入水过程结构强度研究,以期获得的研究成果可为高速回转体模型结构设计提供参考,同时丰富入水过程机理内涵。

    • 基于LS-DYNA通用非线性有限元程序,使用ALE(arbitrary Lagrangian-Eulerian)方法模拟回转体入水,计算入水冲击载荷、平均压力以及速度等的变化,着重考虑此过程中回转体的形变情况。

      在计算中,首先进行Lagrange时步计算,单元网格随材料流动开始变形,然后进行ALE时步计算:

      (1) 保持变形后物体边界条件不变,内部单元进行重分网格,称为光滑步。

      (2) 变形后的网格中的单元变量(密度、能量、应力张量等)和节点速度矢量输运到重新划分的新网格中,称为对流步。

      在LS-DYNA中,流体介质的压力由状态方程进行描述,本文中研究三维回转体-水-空气相互耦合求解过程。对于空气和水,采用LS-DYNA中的NULL材料,对空气介质压力,使用多项式状态方程描述,在LS-DYNA中通过*EOS_LINEAR_POLYNOMIAL关键字施加,状态方程压力公式用下式表示:

      ${p_{\rm{a}}} = {C_0} + {C_1}{\mu _{\rm{a}}} + {C_2}\mu _{\rm{a}}^2 + {C_3} + ({C_4} + {C_5}{\mu _{\rm{a}}} + {C_6}\mu _{\rm{a}}^2){E_{\rm{a0}}}$

      ${\mu _{\rm{a}}} = \frac{1}{{{V_{\rm{a}}}}} - 1$

      式中:pa为气体压力,${C_0}$${C_6}$为自定义常数,${V_{\rm{a}}}$为相对体积,${\mu _{\rm{a}}}$为体积变化率。对理想气体,${C_0}$${C_{\rm{3}}}$${C_6}$均为0,${C_{\rm{4}}}$=${C_{\rm{5}}}$=${\gamma _{\rm{a}}} - 1$, 其中${\gamma _{\rm{a}}}$为单位热值率,初始气体内能${E_{{\rm{a0}}}}$=253.3 kJ

      对于水介质,其压力使用Grüneisen状态方程描述,在LS-DYNA中所使用的关键字为*EOS_Grüneisen,水压力方程表示为:

      ${p_{\rm{w}}} = \frac{{{\rho _{{\rm{w0}}}}c_{\rm{w}}^2{\mu _{\rm{w}}}\left[ {1 + \left( {1 - \displaystyle\frac{{{\gamma _0}}}{2}} \right){\mu _{\rm{w}}} - \displaystyle\frac{\alpha }{2}\mu _{\rm{w}}^2} \right]}}{{{{\left[ {1 - \left( {{S_1} - 1} \right){\mu _{\rm{w}}} - {S_2}\displaystyle\frac{{\mu _{\rm{w}}^2}}{{{\mu _{\rm{w}}} + 1}} - {S_3}\displaystyle\frac{{\mu _{\rm{w}}^3}}{{{{\left( {{\mu _{\rm{w}}} + 1} \right)}^2}}}} \right]}^2}}} + \left( {{\gamma _0} + \alpha {\mu _{\rm{w}}}} \right){E_{{\rm{w0}}}}$

      式中:pw为水压力,${c_{\rm{w}}}$为声音在水中的传播速度,${\mu _{\rm{w}}}$为体积变化率,$\alpha $为对Grüneisen系数${\gamma _0}$的一阶体积修正,${S_1}$${S_2}$${S_3}$分别为${u_{\rm{s}}}$-${u_{\rm{p}}}$曲线斜率无量纲系数,${u_{\rm{s}}}$为冲击波速度,${u_{\rm{p}}}$为流体质点的速度,${E_{{\rm{w0}}}}$为材料初始内能。模拟中水状态方程参数: ${c_{\rm{w}}}$=1 647 m/s,${S_1}$=1.921,${S_2}$=−0.096, S3=0,${\gamma _0}$=0.35,${E_{{\rm{w0}}}}$=289.5 kJ,相对初始体积${V_{{\rm{w0}}}}$=1.0。

    • 建立计算域模型如图1所示,计算域包括空气域和水域。空气域尺寸为1.2 m×1.2 m×1.0 m,水域尺寸为1.2 m×1.2 m×1.5 m。水域和空气域宽度为回转体最大直径的10倍。在模拟中采用了2种回转体,其结构形式分别为:(1)全回转体等厚度;(b)回转体头部部分厚度为8 mm,后体区域赋予一种厚度。回转体模型整体为平头锥形结构,最大直径为120 mm,头部直径为100 mm,回转体长度为0.508 m,其首端倾斜斜面与回转体头部平面所呈角度为104.7°,回转体外形如图2所示。

      图  1  计算域

      Figure 1.  Computational domain

      图  2  回转体几何模型

      Figure 2.  The geometrical model for a revolution body

      在模型单元类型的选择上,弹壳体采用4节点SHELL单元,水和空气采用单元类型为LS-DYNA中模拟ALE物质的SOLID_ALE六面体单元。在入水的初始时刻保证回转体在空气中垂直水面向下,头部贴近水面位置。模拟中应用了ALE多物质算法,空气和水为2种ALE物质。结构和流体通过*CONSTRAINED_LAGRANGE_IN_SOLID进行耦合,同时使用*INITIAL_FRACTION_GEOMETRY关键字将回转体的内部物质设置为空气,这样符合实际情况,模拟更准确。

      由于高速入水过程时间极短,重力对整个入水过程的影响很小,所以在模拟中忽略重力的作用。

      在数值计算过程中,回转体为弹塑性材料,外壳采用45钢材,四边形网格大小为0.005 m。水和空气六面体网格大小为0.02 m,计算模拟中总网格量为520 000。45钢材的密度为7 850 kg/m3,杨氏模量为210 GPa,泊松比为0.3,屈服应力为355 MPa,塑性失效应变为0.35。

    • 为了验证数值方法的准确性,选取第1种工况下板厚为5 mm的回转体进行入水计算。在入水的瞬间回转体头部会受到巨大的冲击载荷,其中心区域有最大的压强峰值。通过数值模拟得到的回转体中心区域在入水抨击过程中的压强曲线,如图3所示,其抨击压强的最大值为189 MPa。参考王珂等[21]对弹性回转体入水抨击载荷峰值的预报,弹性回转体入水抨击压强峰值与回转体的厚度、入水速度以及材料的弹性模量有关,抨击过程峰值压强$p$可以表示为:

      图  3  回转体头部中心区域压强曲线

      Figure 3.  Pressure intensity curve in the central head region of the revolution body

      $ p = {\rho _{\rm{w}}}{k_{d}}{k_{{E}}}{k_{{v}}} $

      式中:${k_{{d}}} = 0.004\;23{d^2} + 0.003\;758d + 0.42$${k_{{E}}} = \ln(E/{\rm{Pa}}) + 15$${k_{{v}}} = {v^2_0} + 5$${\rho _{\rm{w}}}$为水的密度;${k_{{d}}}$为厚度相关系数,d为回转体厚度(m);${k_{{E}}}$为弹性模量相关系数,E为材料弹性模量;${k_{{v}}}$为速度相关系数,v0为回转体入水的初速度(m/s)。由式(4)可以求得在5 mm板厚下的回转体入水冲击载荷压强峰值为173 MPa,和本文求得的数值解比较接近,但数值模拟结果略大于公式求得的抨击压强峰值。其原因是在数值模拟中未考虑空气垫的作用,即高速时的空气压缩效应,这也在一定程度上解释了数值模拟所得压强数值偏大的现象,空气垫效应是数值模拟过程中较难考虑的一个因素,在数值模拟方面还有很大的研究空间,亟待将来的研究。

      除将入水时回转体头部抨击压力的峰值与拟合公式相对比外,还进行了回转体入水过程中速度衰减的对比验证,为了更好地对比速度衰减过程,在此数值模拟中延长了入水过程的求解时间,增大水域深度至3.5 m,计算时间为0.11 s。参考文献[18],根据牛顿第二定律,忽略重力效应。

      假定回转体在入水过程中空泡内外压差保持不变,同时空化数为非定值,有:

      $\sigma = \frac{{\Delta p}}{{{\rho _{\rm{w}}}v_{\rm{p}}^2/2}} = \frac{{{\rho _{\rm{w}}}gh + {C_{\rm{a}}}{\rho _{\rm{a}}}v_0^2/2}}{{{\rho _{\rm{w}}}v_{\rm{p}}^2/2}} \approx \frac{{{C_{\rm{a}}}{\rho _{\rm{a}}}}}{{{\rho _{\rm{w}}}}}\frac{{v_0^2}}{{v_{\rm{p}}^2}} = {\sigma _0}\frac{{v_0^2}}{{v_{\rm{p}}^2}}$

      式中:${\Delta }p$为回转体入水空泡内外压差,${\rho _{\rm{a}}}$${\rho _{\rm{w}}}$分别为空气和水的密度,$v_0^{}$为回转体入水初速度,${v_{\rm{p}}}$为回转体在水中的速度;${C_{\rm{a}}}$为气流压力降因数,${C_{\rm{a}}} = 5 \text{~} 15$${\sigma _0}$为初始空化数,${\sigma _0} = 0.006 \text{~} 0.018$

      柱形平头回转体阻力因数和空化数关系式:

      ${C_{\rm{d}}} = {C_{\rm{p}}}(1 + \sigma )$

      式中:${C_{\rm{p}}} = 0.82 \text{~} 0.83$

      忽略入水过程中重力的影响,回转体有以下运动方程:

      ${m_{\rm{p}}}\frac{{{\rm{d}}{v_{\rm{p}}}}}{{{\rm{d}}t}} = - \frac{1}{2}{\rho _{\rm{w}}}{A_0}{C_{\rm{d}}}v_{\rm{p}}^2$

      式中:${m_{\rm{p}}}$为回转体的质量,${A_0}$为回转体头部面积。

      根据式(5)~(6)得到速度衰减和时间的关系式:

      ${v_{\rm{p}}} = \sqrt {{\sigma _0}} {v_0}\tan \left[ {\arctan \left( {\frac{1}{{\sqrt {{\sigma _0}} }}} \right) - k\sqrt {{\sigma _0}} {v_0}t} \right]$

      式中:衰减系数$k = {\rho _{\rm{w}}}{A_0}{C_{\rm{p}}}/\left( {2{m_{\rm{p}}}} \right)$。依据式(8)计算回转体入水速度衰减曲线并与数值模拟结果进行对比,如图4所示。从图4可以看出,数值模拟结果与理论速度衰减曲线吻合较好,从而进一步验证了所建立的数值方法的有效性。

      图  4  回转体速度衰减曲线

      Figure 4.  Velocity attenuation of the revolution body over time

      在入水数值验证工作的同时,图5给出了当前回转体结构形式下入水过程中回转体头部表面作用力曲线图,可以进一步分析回转体入水过程受力特征。从图5可以看出,回转体头部受到的作用力在入水瞬间急剧上升,达到峰值后迅速下降,最后以小幅度震荡趋势变化。可见在入水瞬间产生了巨大冲击力,所以在设计回转体模型时要重点考虑入水瞬间的结构强度。下面将对2种结构形式的回转体进行强度的分析计算。

      图  5  流-构耦合力曲线

      Figure 5.  Fluid-structure interaction force varying with time

    • 在模拟中,首先进行回转体等厚度工况分析,分别计算回转体厚度为2、3、4和5 mm的情况。计算中发现:只有回转体厚度为5 mm时,该种结构的回转体是安全的;此种结构的回转体发生破坏时,破坏部位均发生在回转体头部和锥形斜面连接处。对回转体厚度为2 mm的计算结果进行分析,其中1.1 ms和3.0 ms时刻的应力云图分别如图6图7所示。由图67可以清晰地看出,破坏部位为头部和锥形斜面连接处,回转体头部中心处未发生破坏现象。提取回转体头部和锥形斜面连接处的塑性应变和单元应力曲线,分别如图8图9所示。由这2条曲线可以得到:在回转体入水过程中,等效应力达到了材料的屈服应力355 MPa,塑性应变也达到了材料的塑性失效应变0.35,因此回转体出现破坏现象。可见,回转体头部和锥形斜面连接处应予以加强,来承受入水瞬间高数量级的冲击载荷。同时因为入水过程中回转体头部为主要的冲击载荷作用面,所以对回转体头部进行加强是十分必要的。

      图  6  在t=1.1 ms时回转体的应力分布

      Figure 6.  Stress distribution in the revolution body at t=1.1 ms

      图  7  在t=3.0 ms时回转体的应力分布

      Figure 7.  Stress distribution in the revolution body at t=3.0 ms

      图  8  回转体头部边缘单元应变曲线

      Figure 8.  Strain-time curve of the edge element for the head of the revolution body

      图  9  回转体头部边缘单元等效应力曲线

      Figure 9.  Equivalent stress-time curve of the edge element for the head of the revolution body

    • 由2.2节中计算结果可知,回转体头部和锥形斜面连接处强度脆弱,因此对回转体头部进行加强,考虑到工程安全系数要求,将头部赋予8 mm的厚度,回转体其余板厚度分别赋予2.0、2.5、3.0、3.5和4.0 mm。在初速度相同的条件下分别进行计算,由计算结果可知:在头部厚度为8 mm、后体区域厚度为2.0 mm时,回转体无法满足强度要求;在后体区域厚度为2.5 mm及以上时,结构强度满足要求。选取后体区域厚度为2.5 mm的情况加以分析。回转体入水过程中0.89 ms和2.00 ms时刻的应力云图分别如图10图11所示,同时图12给出了回转体头部边缘和头部中心单元塑性应变曲线,图13图14分别给出了回转体头部和头部中心单元的等效应力变化曲线。由应力、应变变化曲线可知,在回转体入水过程中,无论其边缘还是中心单元的等效应力很快达到屈服应力355 MPa,随后等效应力下降并表现为波动变化。回转体头部边缘的塑性应变由0迅速到达0.14,并最终趋近0.16,回转体头部中心的塑性应变有同样的变化规律并最终达到0.25。这个塑性应变小于材料的塑性失效应变0.35,因此,材料虽然进入塑性阶段,但还未达到破坏,此时回转体结构满足强度要求。

      图  10  t=0.89 ms时刻回转体的应力分布

      Figure 10.  Stress distribution in the revolution body at t=0.89 ms

      图  11  t=2.00 ms回转体应力图

      Figure 11.  Stress distribution in the revolution body at t=2.00 ms

      图  12  回转体头部边缘及中心单元应变随时间的变化

      Figure 12.  Strain-time curves of the edge and central elements for the head of the revolution body

      图  13  回转体头部边缘单元等效应力曲线

      Figure 13.  Equivalent stress-time curve of the edge element for the head of the revolution body

      图  14  回转体头部中心单元等效应力曲线

      Figure 14.  Equivalent stress-time curve of the central element for the head of the revolution body

    • 采用有限元方法对平头锥头回转体在高速入水过程中的结构强度进行了数值分析,得到的主要结论如下:

      (1) 在回转体入水的瞬间,其头部受到高速冲击载荷,且作用时间极短,在研究回转体高速入水时要重点考虑入水瞬间的作用力。

      (2) 在回转体等厚度结构形式下,等效应力和塑性应变均较大,在入水初速度为100 m/s、回转体厚度小于5.0 mm时,回转体应力达到了材料的屈服应力355 MPa,其塑性应变也达到失效应变0.35,此种形式下回转体厚度低于5.0 mm时无法保证回转体的强度要求。

      (3) 当回转体入水初速度为100 m/s、其头部厚度为8 mm、其后体壁厚为2.5 mm及以上时,结构未发生破坏,此时回转体可以满足强度要求,但材料同样进入了塑性阶段,建议对实际结构还应该考虑一定的安全系数。

参考文献 (21)

目录

    /

    返回文章
    返回