层合多孔圆柱壳的轴向冲击吸能特性

贺璞 邓庆田 李新波

赵亮, 李明哲, 吴楠楠, 王金龙, 梁晓波, 谷洲之, 李怀勇. 缠绕离散式大腔体超高压模具的设计计算及数值模拟[J]. 高压物理学报, 2025, 39(3): 033301. doi: 10.11858/gywlxb.20240851
引用本文: 贺璞, 邓庆田, 李新波. 层合多孔圆柱壳的轴向冲击吸能特性[J]. 高压物理学报, 2022, 36(4): 044203. doi: 10.11858/gywlxb.20210909
ZHAO Liang, LI Mingzhe, WU Nannan, WANG Jinlong, LIANG Xiaobo, GU Zhouzhi, LI Huaiyong. Calculation and Numerical Simulation of Winding Discreted Large Cavity of Ultra-High Pressure Die[J]. Chinese Journal of High Pressure Physics, 2025, 39(3): 033301. doi: 10.11858/gywlxb.20240851
Citation: HE Pu, DENG Qingtian, LI Xinbo. Energy Absorption Characteristics of Laminated Cellular Cylindrical Shell under Axial Impact[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 044203. doi: 10.11858/gywlxb.20210909

层合多孔圆柱壳的轴向冲击吸能特性

doi: 10.11858/gywlxb.20210909
基金项目: 中央高校基本科研业务费专项资金(310812163504)
详细信息
    作者简介:

    贺 璞(1996-),男,硕士研究生,主要从事多孔结构承载与吸能研究.E-mail:hepu96@126.com

    通讯作者:

    邓庆田(1980-),男,博士,副教授,主要从事多孔材料与结构力学性能研究.E-mail:dengqt@chd.edu.cn

  • 中图分类号: O347.1

Energy Absorption Characteristics of Laminated Cellular Cylindrical Shell under Axial Impact

  • 摘要: 为提高圆柱壳结构在轴向冲击载荷下的力学响应,设计了多层层合多孔圆柱壳结构,通过改变层合方式以及孔单元数目研究其力学行为。通过准静态压缩实验验证有限元分析的有效性,使用有限元模拟方法研究了结构在轴向冲击载荷下的变形模式和吸能特性。从层合多孔圆柱壳的压缩力-位移曲线得到总吸能、峰值压缩力和平均压缩力等指标。对比不同模型之间的吸能和压缩力特性,结果表明:改变层合方式对圆柱壳吸能效果影响较显著,增加孔单元数目对多孔圆柱壳吸能效果影响较小。对比正、负泊松比两种多孔圆柱壳,同等质量内凹圆柱壳总吸能较六边圆柱壳平均提高了17%。

     

  • 全浸水带间隙发射作为一种新的水下发射方式[1],通过间隙燃气在膛内贴壁运动,卷吸回流后逐渐汇聚成弹前气幕,排出弹前水柱,将射弹在膛口的发射环境由水介质转化为气体介质。当气体射流流出枪口后,膛口处射流迅速膨胀成球形气体空腔,射弹穿过气体空腔与水介质接触,产生强烈的冲击载荷。水下高速射弹的弹体主要由硬铝合金尾杆和钨合金头部组成,两者镶嵌连接,连接强度有限,入水瞬间弹头会承受巨大冲击,因此射弹的入水冲击载荷成为水下射弹设计中的一个重要问题。

    Karman[2]最先开始对入水冲击现象进行研究,采用动量定理并引入附加质量的概念,推导出入水冲击载荷的计算公式。Wagner[3]将Karman的方法理论化,提出了近似平板理论及自相似解法,得出了冲击压力在结构沾湿面的分布情况,使理论分析更加符合实际情况,为后来学者的理论研究奠定了基础。在国内,秦洪德等[4]、王永虎等[5]对入水冲击问题的现状和进展进行了详细的分析。宋保维等[6]基于不可压缩的非定常势流理论,建立了空投水雷入水冲击计算的数学模型。卢炽华等[7]利用不同浸深的附加质量,对刚性细长体斜姿态落水冲击进行建模,得出其入水角很小,会使弹体处在最危险的状态。王永虎等[8-9]先后对刚性尖拱体垂直姿态高速入水和斜入水的冲击理论进行了建模和仿真。魏卓慧等[10]建立了刚性截锥弹体垂直入水冲击载荷的数学模型,并对其进行了数值计算。陈诚等[11]对超空泡航行器倾斜入水冲击载荷进行了试验研究,得出了峰值时刻的阻力系数。朱珠等[12]利用商业软件FLUENT,建立了柱体回转体高速入水冲击的数值模拟模型,得到了速度对入水冲击载荷的影响规律。然而以上研究主要针对由空中入水的冲击载荷分析,对于水下入水冲击问题研究较少。本工作在此前提条件下,计算分析全水下发射高速射弹入水的冲击载荷,这对于水下发射武器研究具有一定的现实意义。

    本研究拟建立锥形弹体水平及斜入水的冲击载荷理论模型,模型中考虑弹体重力、弹体浮力、附加质量、弹头锥角及入水攻角的影响,对不同头部结构参数的锥形弹体以不同入水速度入水的冲击载荷进行计算,分析入水速度、弹头锥角和入水攻角对冲击载荷的影响。研究结果对于弹体入水冲击载荷的预测及全水下发射方式发射的射弹头部结构设计具有参考价值。

    假设射弹为刚体,不考虑射弹入水时空泡的影响,根据动量定理,射弹高速入水冲击时的动量方程为

    Mv0=(M+m)v+Mgt+Fbt+Fdt
    (1)

    式中:M为射弹的质量,m为射弹的附加质量,Fd为射弹入水时受到的阻力,Fb为射弹所受的浮力。

    对式(1)等号两边进行微分,得到射弹入水冲击时的动力学方程

    (M+m)d2hdt2+dmdtdhdt+12ρCdA(dhdt)2=0
    (2)

    式中:h为射弹侵入的距离,A为阻力面积,Cd为阻力系数,ρ为水的密度。

    射弹沿 x 轴水平入水时,由于入水冲击过程的瞬时性,入水初期其运动方向基本保持不变,流体动力主要作用于射弹轴线方向,在射弹轴线方向上重力和浮力对射弹影响很小,基本可忽略不计。

    射弹入水冲击过程中,射弹浸没在水中的体积会排挤液面流体产生隆起现象,如图1所示,有效液面决定了自由水平液面的抬高程度,这取决于射弹锥头的外形和入水角等初始状态。沾湿因子定义为有效液面与实际液面的比值。

    图  1  射弹水平入水冲击示意图
    Figure  1.  Schematic diagram of horizontal impact of the projectile

    利用轴长体假设和Tayler关于在不同浸深时附加质量的表达式[13],参考垂直入水相关文献[10, 14],求出锥形射弹入水的附加质量为

    m=43ρ(Cwh)3tan3β
    (3)
    A=π(Cwh)2tan2β
    (4)
    dmdt=dmdhdhdt=4ρtan3β(Cwh)2Cwdhdt
    (5)

    式中:Cw为沾湿因子,β为射弹头部的半锥角,如图1所示。

    将式(3)~式(5)代入式(2),得到总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2=0
    (6)

    采用MATLAB软件,利用龙格-库塔方法进行求解,可以得到弹体入水时的冲击载荷。

    图2所示,采用全新的发射方式时,燃气排出并在膛口形成球形气体空腔。以地面为坐标系,枪口斜向上(即y轴正方向)、斜向下(即y轴负方向)发射时,根据圆切线定理,穿过球形气腔仍可看作垂直于液面入水,然而射弹入水在有攻角的情况下,轴向上会受到重力和浮力的作用分力影响,攻角正负值相反时,重力与浮力在射弹轴向上的作用分力方向也相反。

    图  2  射弹攻角示意图
    Figure  2.  Angle of attack of the projectile

    其他条件与水平入水时保持不变,攻角为正时,根据动量方程,得到动力学方程

    (M+m)d2hdt2+dmdtdhdt+12ρCdA(dhdt)2+MgsinαFbsinα=0
    (7)

    式中:α为射弹攻角,射弹斜向上发射时为正值,斜向下时为负值。

    射弹所受浮力为

    Fb=13πρg(Cwh)3tan2β
    (8)

    当攻角为正值时,得到的总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2+Mgsinα13πρg(Cwh)3tan2βsinα=0
    (9)

    同理,攻角为负值时的总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2Mgsinα+13πρg(Cwh)3tan2βsinα=0
    (10)

    对式(9)、式(10)进行求解,可以得到不同攻角下的入水冲击载荷。

    由于水下射弹的质量较轻,入水速度较大,水下超空泡射弹入水平均速度约为600 m/s,因此入水速度对射弹入水冲击载荷的影响很大。本计算中,设射弹质量为0.14 kg,射弹头部半锥角为6°,计算得到不同速度时射弹的入水冲击载荷曲线,如图3所示,其中用射弹加速度反映入水冲击载荷。可以看出:射弹的入水冲击载荷先增大后减小,载荷变化主要发生在射弹入水前1 ms,最后渐渐趋于稳定;入水速度越大,冲击载荷峰值越大,入水后达到峰值的时间越短。此外,计算了不同速度下的入水冲击载荷峰值,如图4所示。入水速度在400~700 m/s范围时,射弹入水冲击载荷峰值一般为103g量级(g为重力加速度)。从冲击载荷峰值与速度的关系可以得出入水冲击载荷的峰值与速度基本呈线性关系。

    图  3  不同速度时的入水冲击载荷
    Figure  3.  Driving impact loads at different velocities
    图  4  不同速度时的冲击载荷峰值
    Figure  4.  Peak impact loads at different velocities

    设射弹质量为0.14 kg,通过计算获得了不同锥角的锥形射弹以600 m/s入水时的冲击载荷曲线,如图5所示。可以看出:射弹的入水冲击载荷先增大后减小,最后趋于稳定;半锥角β越大,冲击载荷峰值越大,并且入水后达到峰值的时间越短。改变半锥角,计算出不同半锥角情况下射弹入水冲击载荷峰值,如图6所示。从冲击载荷峰值与半锥角的关系可以得出入水冲击载荷峰值与半锥角基本呈线性关系。

    图  5  不同锥角时的入水冲击载荷
    Figure  5.  Driving impact loads at different cone angles
    图  6  不同锥角时的冲击载荷峰值
    Figure  6.  Peak impact loads at different cone angles

    当射弹质量为0.14 kg,入水速度为600 m/s,攻角α分别取45°、−45 °和0°(即水平入水)时,计算得到的射弹入水冲击载荷曲线如图7所示。可见,3条曲线基本重叠,差值在1g量级,相比于速度和半锥角对冲击载荷的影响,攻角对冲击载荷的影响几乎可以忽略不计。

    图  7  不同攻角时的入水冲击载荷
    Figure  7.  Inlet impact loads at different angles of attack

    此外,射弹在高速入水状态下,从液面进入水中的实际深度与理想深度存在一定的偏差,为了表达该误差,引入沾湿因子,沾湿因子的取值对冲击载荷的结果也存在影响。沾湿因子不同时射弹的入水冲击载荷如图8所示。从图8可以看出:沾湿因子越大,入水冲击载荷峰值越大;但沾湿因子对入水冲击载荷的影响较小,当沾湿因子变化值为0.1时,入水冲击载荷的变化在10%以内。

    图  8  沾湿因子不同时的入水冲击载荷
    Figure  8.  Inlet impact loads with different wet factors

    (1)射弹头部锥角相同时,入水冲击载荷峰值与速度呈正线性相关,入水速度越大,冲击载荷达到峰值的时间越短;射弹入水速度相同时,入水冲击载荷的大小与锥角呈正线性相关,锥角越大,冲击载荷达到峰值的时间越短。

    (2)射弹锥角和入水速度相同、入水攻角不同时,入水冲击载荷曲线与水平入水曲线基本重合,说明射弹重力和浮力在轴向上对入水冲击载荷的影响相对入水阻力几乎可以忽略不计。当射弹带有攻角入水后,重力和浮力更多的是对射弹产生径向力矩影响。

    (3)沾湿因子越大,入水冲击载荷峰值越大,但对入水冲击载荷影响较小,当沾湿因子变化为0.1时,入水冲击载荷的变化在10%以内。沾湿因子作为一个变量,其大小反过来也取决于入水冲击载荷,入水冲击载荷越大,沾湿因子也越大。

    (4)理论模型借鉴了高速弹体垂直入水的理论模型,数值模拟计算结果与已报道的高速射弹垂直入水冲击载荷数值计算和仿真结果高度一致,验证了本数学建模和数值计算的准确性。

    (5)探讨了全水下发射高速射弹入水瞬间的冲击载荷,高速射弹入水后形成超空泡,冲击载荷迅速减小。本工作对射弹未形成超空泡的情况进行了模型推导和数值模拟计算,对形成超空泡之前的理论研究具有重要意义。

  • 图  不同层合方式的圆柱壳模型

    Figure  1.  Cylindrical shell models withdifferent laminated ways

    图  不同孔单元数目的圆柱壳模型

    Figure  2.  Cylindrical shell models withdifferent cells number

    图  PLA层合多孔圆柱壳

    Figure  3.  Laminated cellular cylindrical shell of PLA material

    图  ABAQUS有限元模型

    Figure  4.  Finite element model in ABAQUS

    图  PLA材料的压缩实验与应力-应变曲线

    Figure  5.  Compression experiment and stress-strain curve of PLA material

    图  实验与数值模拟得到的力-位移曲线对比

    Figure  6.  Comparisons of compressive force-displacement curves between experiment and simulation

    图  实验与数值模拟得到的变形模式对比

    Figure  7.  Comparisons of deformation mode between experiment and simulation

    图  3种尺寸网格的有限元仿真模型

    Figure  8.  Finite element simulation models of three kinds of meshes

    图  网格敏感性验证

    Figure  9.  Validation of mesh sensitivity

    图  10  伪应变能与内能之比

    Figure  10.  Ratio of artificial strain energy to internal energy

    图  11  六边多孔圆柱壳在不同冲击速度下的变形模式

    Figure  11.  Deformation mode of hexagonal cellular cylindrical shell at different impact velocities

    图  12  内凹多孔圆柱壳变形模式

    Figure  12.  Deformation mode of re-entrant cellular cylindrical shell

    图  13  不同层合方式多孔圆柱壳的力-位移曲线

    Figure  13.  Compressive force-displacement curves of cellular cylindrical shell with different laminated ways

    图  14  不同孔数多孔圆柱壳的力-位移曲线

    Figure  14.  Compressive force-displacement curves of cellular cylindrical shell with different cells number

    表  1  层合多孔圆柱壳的几何尺寸

    Table  1.   Geometric dimensions of laminated cellular cylindrical shell

    Number of axial cellNumber of circumferential cellThickness of each layer/mmHeight/mmDiameter/mm
    8242/210678.4
    12 362/210678.4
    16 482/210678.4
    8241/2/110678.4
    8241/1/1/110678.4
    下载: 导出CSV

    表  2  不同层合方式下的冲击吸能指标

    Table  2.   Impact energy absorption parameters of different laminated way

    Cellular unitsLayer numberEt/JFp/kNFa/kN
    Hexagon21 206.029.912.2
    31 032.143.911.0
    41 414.332.614.5
    Re-entrant21 416.730.813.7
    31 264.648.611.8
    41 673.532.617.0
    下载: 导出CSV

    表  3  不同孔数多孔圆柱壳的冲击吸能指标

    Table  3.   Impact energy absorption parameters of of cellular cylindrical shell with different cells number

    Cellular unitsCell numberEt/JFp/kNFa/kN
    Hexagon8×241 206.029.912.2
    12×361 169.725.811.2
    16×481 266.429.411.9
    Re-entrant8×241 416.730.813.7
    12×361 383.428.913.4
    16×481 394.039.613.3
    下载: 导出CSV
  • [1] 解江, 张雪晗, 苏璇, 等. 铺层顺序对复合材料薄壁圆管轴向压溃吸能特性的影响研究 [J]. 工程力学, 2018, 35(6): 231–239.

    XIE J, ZHANG X H, SU X, et al. Influence of layer sequence on energy absorption characteristics of thin-walled composite circular tubes under axial compression [J]. Engineering Mechanics, 2018, 35(6): 231–239.
    [2] 李金柱, 齐凯丽, 袁婷, 等. 冲击载荷下金属薄壁圆柱壳动态响应特性研究 [J]. 北京理工大学学报, 2020, 40(4): 362–368.

    LI J Z, QI K L, YUAN T, et al. Dynamic response of thin-walled cylindrical shells suffering impact loading [J]. Transactions of Beijing Institute of Technology, 2020, 40(4): 362–368.
    [3] 何强, 王勇辉, 史肖娜, 等. 引入Sierpinski层级特性的新型薄壁多胞管轴向冲击吸能特性 [J]. 爆炸与冲击, 2020, 40(12): 123101.

    HE Q, WANG Y H, SHI X N, et al. Energy absorption of new thin-walled, multi-cellular, tubular structures with Sierpinski hierarchical characteristics under axial impact [J]. Explosion and Shock Waves, 2020, 40(12): 123101.
    [4] 张雨, 李应刚, 沈云龙, 等. 蜂窝金属夹芯板重复冲击动态响应研究 [J]. 振动与冲击, 2021, 40(4): 255–260.

    ZHANG Y, LI Y G, SHEN Y L, et al. Dynamic responses of honeycomb sandwich panels under repeated impacts [J]. Journal of Vibration and Shock, 2021, 40(4): 255–260.
    [5] 虞科炯, 徐峰祥, 华林. 正弦曲边负泊松比蜂窝结构面内冲击性能研究 [J]. 振动与冲击, 2021, 40(13): 51–59.

    YU K J, XU F X, HUA L. In plane impact performance of honeycomb structure with sinusoidal curved edge and negative Poisson’s ratio [J]. Journal of Vibration and Shock, 2021, 40(13): 51–59.
    [6] 王璠, 何一帆, 宋毅, 等. 引发方式、铺层对纤维增强复合材料圆柱壳吸能特性影响的冲击试验研究 [J]. 振动工程学报, 2013, 26(1): 33–40. doi: 10.3969/j.issn.1004-4523.2013.01.006

    WANG F, HE Y F, SONG Y, et al. The impact experimental study of evoking type and Layup type on absorption energy of fiber reinforced composite cylindrical shell [J]. Journal of Vibration Engineering, 2013, 26(1): 33–40. doi: 10.3969/j.issn.1004-4523.2013.01.006
    [7] 罗国强, 费细欢, 喻寅, 等. 孔洞排布对PMMA多孔材料冲击响应行为的影响 [J]. 高压物理学报, 2020, 34(5): 054202. doi: 10.11858/gywlxb.20200542

    LUO G Q, FEI X H, YU Y, et al. Effect of voids arrangement on behavior of PMMA cellular materials on impact loading [J]. Chinese Journal of High Pressure Physics, 2020, 34(5): 054202. doi: 10.11858/gywlxb.20200542
    [8] BICH D H, DUNG D V, NAM V H, et al. Nonlinear static and dynamic buckling analysis of imperfect eccentrically stiffened functionally graded circular cylindrical thin shells under axial compression [J]. International Journal of Mechanical Sciences, 2013, 74: 190–200. doi: 10.1016/j.ijmecsci.2013.06.002
    [9] ZHANG X, ZHANG H. Energy absorption of multi-cell stub columns under axial compression [J]. Thin-Walled Structures, 2013, 68: 156–163. doi: 10.1016/j.tws.2013.03.014
    [10] 杨欣, 范晓文, 许述财, 等. 仿虾螯结构薄壁管设计及耐撞性分析 [J]. 爆炸与冲击, 2020, 40(4): 043301. doi: 10.11883/bzycj-2019-0280

    YANG X, FAN X W, XU S C, et al. Design and crashworthiness analysis of thin-walled tubes based on a shrimp chela structure [J]. Explosion and Shock Waves, 2020, 40(4): 043301. doi: 10.11883/bzycj-2019-0280
    [11] 刘亚军, 何玉龙, 刘姗姗, 等. 正多边形基多胞薄壁管的吸能特性 [J]. 爆炸与冲击, 2020, 40(7): 071404. doi: 10.11883/bzycj-2019-0423

    LIU Y J, HE Y L, LIU S S, et al. Energy absorption capacity of regular polygon-based multi-cell tubes [J]. Explosion and Shock Waves, 2020, 40(7): 071404. doi: 10.11883/bzycj-2019-0423
    [12] 张天辉, 邓健强, 刘志芳, 等. 弧形折纸模式薄壁结构的压缩变形与能量吸收 [J]. 爆炸与冲击, 2020, 40(7): 071405. doi: 10.11883/bzycj-2019-0355

    ZHANG T H, DENG J Q, LIU Z F, et al. Compression deformation and energy absorption of thin-walled structures with arc-shaped origami patterns [J]. Explosion and Shock Waves, 2020, 40(7): 071405. doi: 10.11883/bzycj-2019-0355
    [13] 杨智春, 邓庆田. 负泊松比材料与结构的力学性能研究及应用 [J]. 力学进展, 2011, 41(3): 335–350. doi: 10.6052/1000-0992-2011-3-lxjzJ2011-013

    YANG Z C, DENG Q T. Mechanical property and application of materials and structures with negative Poisson’s ratio [J]. Advances in Mechanics, 2011, 41(3): 335–350. doi: 10.6052/1000-0992-2011-3-lxjzJ2011-013
    [14] 张越, 李世强, 王志华, 等. 二阶层级自相似四边形蜂窝动力压缩行为数值模拟 [J]. 高压物理学报, 2017, 31(4): 358–363. doi: 10.11858/gywlxb.2017.04.002

    ZHANG Y, LI S Q, WANG Z H, et al. Dynamic crushing response of seif-similar second order hierarchical square honeycombs [J]. Chinese Journal of High Pressure Physics, 2017, 31(4): 358–363. doi: 10.11858/gywlxb.2017.04.002
    [15] 王根伟, 刘冕, 宋辉, 等. 冲击载荷下径向密度排布对泡沫金属力学性能影响的研究 [J]. 爆炸与冲击, 2020, 40(7): 071401. doi: 10.11883/bzycj-2019-0403

    WANG G W, LIU M, SONG H, et al. Influence of radial density arrangement on mechanical properties of metal foam under impact loading [J]. Explosion and Shock Waves, 2020, 40(7): 071401. doi: 10.11883/bzycj-2019-0403
    [16] 王仁, 韩铭宝, 黄筑平, 等. 受轴向冲击的圆柱壳塑性动力屈曲实验研究 [J]. 力学学报, 1983, 19(5): 509–514.

    WANG R, HAN M B, HUANG Z P, et al. Dynamic plastic buckling of cylindrical shell under axial impact: an experimental study [J]. Chinese Journal of Theoretical and Applied Mechanics, 1983, 19(5): 509–514.
  • 加载中
图(14) / 表(3)
计量
  • 文章访问数:  325
  • HTML全文浏览量:  165
  • PDF下载量:  34
出版历程
  • 收稿日期:  2021-11-23
  • 修回日期:  2021-12-11
  • 录用日期:  2022-04-07
  • 网络出版日期:  2022-06-23
  • 刊出日期:  2022-07-28

目录

/

返回文章
返回