金属蜂窝夹芯结构抗水下冲击性能

李汶蔚 黄威

引用本文:
Citation:

金属蜂窝夹芯结构抗水下冲击性能

    作者简介: 李汶蔚(1985-),男,硕士,副研究员,主要从事爆轰物理研究. Email:wenweili@gmail.com;
    通讯作者: 黄威, weihuang@hust.edu.cn
  • 中图分类号: O347

Impulsive Resistance of Metallic Honeycomb Sandwich Structures Subjected to Underwater Impulsive Loading

    Corresponding author: HUANG Wei, weihuang@hust.edu.cn
  • CLC number: O347

  • 摘要: 为揭示高强度水下爆炸冲击载荷作用下金属夹芯结构的抗冲击性能,在实验室开展小尺寸水下爆炸加载技术对金属蜂窝夹芯结构性能影响的实验研究。基于实验结果,开展了全尺寸数值模拟金属蜂窝夹芯结构在水下冲击载荷作用下的动态响应和抗冲击性能研究。结果表明,数值模拟、实验和理论模型计算的结果具有良好的一致性。由于蜂窝芯材相对密度对夹芯结构能量耗散方式和载荷传递机制的影响,结构动态响应、失效模式以及抗冲击性能随着冲击强度的变化表现出较为明显的不同。通过抗冲击参数分析,建立了反映金属蜂窝夹芯结构抗冲击性能的结构横向变形、固支反力、透射脉冲和塑性能耗随冲击强度和芯材相对密度变化的结构-载荷-性能量化关系。
  • 图 1  高强度水下冲击加载实验装置示意图(a)和数值分析模型(b)

    Figure 1.  Intensive underwater explosive simulator (a) and the numerical model (b)

    图 2  蜂窝夹芯结构准静态压缩下的应力-应变关系

    Figure 2.  Stress-strain relationship compression of honeycomb sandwich structure under quasi-static compression

    图 3  实验和数值模拟的失效模式对比

    Figure 3.  Comparison of the deformation modes obtained from the simulation and experiment

    图 4  蜂窝夹芯结构背板中点变形时程曲线(a)和最大变形与冲击强度的关系(b)

    Figure 4.  Mid-point deflection history of honeycomb sandwish structure (a) and relationship between impulsive intensities and maximum deflection (b)

    图 5  蜂窝夹芯结构在${10^3}{\overline I _{\rm{t}}}$ = 3.58作用下的等效应变分布和塑性铰运动

    Figure 5.  Dynamic deformation and propagation of plastic hinges of the honeycomb sandwich,${10^3}{\overline I _{\rm{t}}} $ = 3.58

    图 6  蜂窝夹芯结构前后面板中点的速度响应过程

    Figure 6.  Velocity time history at the mid-point of front and back face-sheet of honeycomb sandwich structure

    图 7  相同冲击强度下不同芯材相对密度的蜂窝夹芯结构等效应力分布

    Figure 7.  Equivalent strain distribution of honeycomb sandwich under the same impact loading

    图 8  蜂窝夹芯结构横向变形与冲击强度及芯材相对密度的关系

    Figure 8.  The relationship of transverse deflections for honeycomb sandwich structure load intensity and relative core density

    图 9  夹芯结构支座的反力时程曲线

    Figure 9.  Reaction forces time histories for sandwich panels

    图 10  夹芯结构透射脉冲强度与冲击强度及几何特性的关系

    Figure 10.  Relationship of transmitted impulses to incident impulse intensity and geometries

    图 11  蜂窝夹芯结构各部分的塑性能耗时程曲线(a)及塑性能耗与冲击强度和芯材相对密度关系(b)

    Figure 11.  Plastic dissipation histories at different places of honeycomb sandwich structure (a) and the relationship of plastic dissipations to incident impulse intensity and relative core density (b)

    表 1  面板和芯材材料力学性能参数

    Table 1.  Mechanical parameters of aluminum materials

    MaterialsYoung’s modulus/GPaDensity/(kg·m−3)Parameters
    A/MPaB/MPaCn
    5A06 aluminium alloy74.02 780167.0443.70.0200.44
    3003 aluminium alloy74.22 700 85.2170.00.0380.44
    下载: 导出CSV

    表 2  Mie-Grüneisen状态方程参数

    Table 2.  Parameters for the Mie-Grüneisen equation of state

    Density/(kg·m−3)Sound speed in water/(m·s−1)γ
    1 0001 1060.05
    下载: 导出CSV
  • [1] WADLEY H, DHARMASENA K, CHEN Y, et al. Compressive response of multilayered pyramidal lattices during underwater shock loading [J]. International Journal of Impact Engineering, 2008, 35(9): 1102–1114. doi: 10.1016/j.ijimpeng.2007.06.009
    [2] FLECK N A, DESHPANDE V S. The resistance of clamped sandwich beams to shock loading [J]. Journal of Applied Mechanics, 2004, 71(3): 386–401. doi: 10.1115/1.1629109
    [3] DESHPANDE V S, FLECK N A. One-dimensional response of sandwich plates to underwater shock loading [J]. Journal of the Mechanics and Physics of Solids, 2005, 53(11): 2347–2383. doi: 10.1016/j.jmps.2005.06.006
    [4] MCSHANE G J, DESHPANDE V S, FLECK N A. The underwater blast resistance of metallic sandwich beams with prismatic lattice cores [J]. Journal of Applied Mechanics, 2007, 74(2): 352–364. doi: 10.1115/1.2198549
    [5] HUANG W, ZHANG W, YE N, et al. Dynamic response and failure of PVC foam core metallic sandwich subjected to underwater impulsive loading [J]. Composites Part B: Engineering, 2016, 97: 226–238. doi: 10.1016/j.compositesb.2016.05.015
    [6] HUANG W, ZHANG W, HUANG X L, et al. Dynamic response of aluminum corrugated sandwich subjected to underwater impulsive loading: experiment and numerical modeling [J]. International Journal of Impact Engineering, 2017, 109: 78–91. doi: 10.1016/j.ijimpeng.2017.06.002
    [7] HUANG W, ZHANG W, CHEN T, et al. Dynamic response of circular composite laminates subjected to underwater impulsive loading [J]. Composites Part A: Applied Science and Manufacturing, 2018, 109: 63–74. doi: 10.1016/j.compositesa.2018.02.043
    [8] AVACHAT S, ZHOU M. High-speed digital imaging and computational modeling of dynamic failure in composite structures subjected to underwater impulsive loads [J]. International Journal of Impact Engineering, 2015, 77: 147–165. doi: 10.1016/j.ijimpeng.2014.11.008
    [9] TAYLOR G. The pressure and impulse of submarine explosion waves on plates [M]//BATCHELOR G K. The Scientific Papers of Sir Geoffrey Ingram Taylor. Cambridge: Cambridge University Press, 1963: 287–303.
    [10] LEE S, BARTHELAT F, HUTCHINSON J W, et al. Dynamic failure of metallic pyramidal truss core materials-experiments and modeling [J]. International Journal of Plasticity, 2006, 22(11): 2118–2145. doi: 10.1016/j.ijplas.2006.02.006
    [11] MCSHANE G J, RADFORD D D, DESHPANDE V S, et al. The response of clamped sandwich plates with lattice cores subjected to shock loading [J]. European Journal of Mechanics A Solids, 2006, 25(2): 215–229. doi: 10.1016/j.euromechsol.2005.08.001
    [12] RADFORD D D, MCSHANE G J, DESHPANDE V S, et al. The response of clamped sandwich plates with metallic foam cores to simulated blast loading [J]. International Journal of Solids and Structures, 2006, 43(7/8): 2243–2259. doi: 10.1016/j.ijsolstr.2005.07.006
    [13] QIN Q H, WANG T J, ZHAO S Z. Large deflections of metallic sandwich and monolithic beams under locally impulsive loading [J]. International Journal of Mechanical Sciences, 2009, 51(11/12): 752–773.
    [14] CUI X D, ZHAO L M, WANG Z H, et al. Dynamic response of metallic lattice sandwich structures to impulsive loading [J]. International Journal of Impact Engineering, 2012, 43: 1–5. doi: 10.1016/j.ijimpeng.2011.11.004
    [15] HUANG W, ZHANG W, LI D C, et al. Analytical model of the dynamic response of clamped metallic sandwich beam subjected to underwater impulsive loading [J]. Marine Structures, 2019, 63: 333–350. doi: 10.1016/j.marstruc.2018.08.008
    [16] 任鹏, 张伟, 刘建华, 等. 高强度水下爆炸等效冲击波加载特性研究 [J]. 兵工学报, 2015, 36(4): 716–722. doi: 10.3969/j.issn.1000-1093.2015.04.021
    REN P, ZHANG W, LIU J H, et al. Characteristics of high strength underwater explosion equivalent shock loading [J]. Acta Armamentarii, 2015, 36(4): 716–722. doi: 10.3969/j.issn.1000-1093.2015.04.021
    [17] 任鹏, 张伟, 黄威, 等. 非药式水下爆炸冲击波加载装置研究 [J]. 爆炸与冲击, 2014, 34(3): 334–339. doi: 10.11883/1001-1455(2014)03-0334-06
    REN P, ZHANG W, HUANG W, et al. Research on non-explosive underwater shock loading device [J]. Explosion and Shock Waves, 2014, 34(3): 334–339. doi: 10.11883/1001-1455(2014)03-0334-06
    [18] HUANG W, ZHANG W, LI D C, et al. Dynamic failure of honeycomb-core sandwich structures subjected to underwater impulsive loads [J]. European Journal of Mechanics A Solids, 2016, 60: 39–51. doi: 10.1016/j.euromechsol.2016.06.006
    [19] DESHPANDE V S, HEAVER A, FLECK N A. An underwater shock simulator [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2006, 462(2067): 1021–1041. doi: 10.1098/rspa.2005.1604
    [20] ESPINOSA H D, LEE S, MOLDOVAN N. A novel fluid structure interaction experiment to investigate deformation of structural elements subjected to impulsive loading [J]. Experimental Mechanics, 2006, 46(6): 805–824. doi: 10.1007/s11340-006-0296-7
    [21] HUANG W, ZHANG W, REN P, et al. An experimental investigation of water-filled tank subjected to horizontal high speed impact [J]. Experimental Mechanics, 2015, 55(6): 1123–1138. doi: 10.1007/s11340-015-0012-6
    [22] TILBROOK M T, DESHPANDE V S, FLECK N A. Underwater blast loading of sandwich beams: regimes of behaviour [J]. International Journal of Solids and Structures, 2009, 46(17): 3209–3221. doi: 10.1016/j.ijsolstr.2009.04.012
  • [1] 丁铁张晓晴姚小虎 . 筋条形状对复合材料加筋壁板低速冲击动态响应的影响. 高压物理学报, 2017, 31(6): 769-777. doi: 10.11858/gywlxb.2017.06.012
    [2] 贾宪振胡毅亭董明荣许学忠刘家骢 . 深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究. 高压物理学报, 2008, 22(2): 208-214 . doi: 10.11858/gywlxb.2008.02.016
    [3] 李平李大红宁建国经福谦 . 冲击载荷下Al2O3陶瓷的动态响应. 高压物理学报, 2002, 16(1): 22-28 . doi: 10.11858/gywlxb.2002.01.004
    [4] 刘姗姗刘亚军张英杰李志强 . 碳纤维-泡沫铝夹芯板低速冲击响应. 高压物理学报, 2020, 34(3): 034202-1-034202-10. doi: 10.11858/gywlxb.20190872
    [5] 李焰张向荣谭红梅刘晓玲裴明敬 . 国产PVDF压电薄膜的冲击加载及卸载响应研究. 高压物理学报, 2004, 18(3): 261-266 . doi: 10.11858/gywlxb.2004.03.011
    [6] 段卓平关智勇 . 冲击波中横向、纵向应力计响应初步研究. 高压物理学报, 2002, 16(4): 265-270 . doi: 10.11858/gywlxb.2002.04.005
    [7] 朱锡李海涛牟金磊黄晓明 . 水下近距爆炸作用下船体梁的动态响应特性. 高压物理学报, 2010, 24(5): 343-350 . doi: 10.11858/gywlxb.2010.05.005
    [8] 宿华祥易伟建 . 钢筋混凝土墙抗冲击性能的数值模拟分析. 高压物理学报, 2020, 34(1): 014201-1-014201-14. doi: 10.11858/gywlxb.20190772
    [9] 李汶蔚梅杰黄威 . 碳纤维增强复合材料层合板的抗冲击性能. 高压物理学报, 2020, 34(2): 024101-1-024101-8. doi: 10.11858/gywlxb.20190822
    [10] 郑松林 . 晶体塑性有限元在材料动态响应研究中的应用进展. 高压物理学报, 2019, 33(3): 030108-1-030108-21. doi: 10.11858/gywlxb.20190725
    [11] 王光勇曹安生余锐林加剑 . 顶爆和拱腰侧爆同时作用下锚固洞室的动态响应. 高压物理学报, 2020, 34(2): 025202-1-025202-8. doi: 10.11858/gywlxb.20190812
    [12] 徐明利张若棋张光莹 . 用元胞法研究层合材料中的应力波. 高压物理学报, 2000, 14(3): 182-188 . doi: 10.11858/gywlxb.2000.03.005
    [13] 任会兰郭婷婷宁建国 . 聚能射流对氧化铝陶瓷靶的侵彻特性研究. 高压物理学报, 2011, 25(6): 526-532. doi: 10.11858/gywlxb.2011.06.008
    [14] 陈鹏侯秀慧张凯 . 面内冲击荷载下半凹角蜂窝的抗冲击特性. 高压物理学报, 2019, 33(6): 064104-1-064104-13. doi: 10.11858/gywlxb.20190759
    [15] 段卓平郁锐张连生黄风雷 . 陶瓷材料在压剪联合冲击加载下动态响应的实验研究. 高压物理学报, 2007, 21(4): 337-341 . doi: 10.11858/gywlxb.2007.04.001
    [16] 吴成Gupta Y M . 冲击载荷下非晶态碳材料的动态响应模型研究. 高压物理学报, 2002, 16(1): 34-41 . doi: 10.11858/gywlxb.2002.01.006
    [17] 谷卓伟金孝刚张清福孙悦 . 材料预加热冲击压缩实验技术及高温下不锈钢的动态响应. 高压物理学报, 1998, 12(3): 190-198 . doi: 10.11858/gywlxb.1998.03.005
    [18] 冯根柱于博丽李世强刘志芳 . 多层级夹芯结构的变形与能量吸收. 高压物理学报, 2019, 33(5): 055902-1-055902-10. doi: 10.11858/gywlxb.20180707
    [19] 朱志武宁建国刘煦 . 冲击载荷下土的动态力学性能研究. 高压物理学报, 2011, 25(5): 444-450 . doi: 10.11858/gywlxb.2011.05.010
    [20] 张新春刘颖 . 缺陷对金属蜂窝材料面内冲击性能的影响. 高压物理学报, 2012, 26(6): 645-652. doi: 10.11858/gywlxb.2012.06.008
  • 加载中
图(11)表(2)
计量
  • 文章访问数:  771
  • 阅读全文浏览量:  718
  • PDF下载量:  5
出版历程
  • 收稿日期:  2019-06-10
  • 录用日期:  2019-06-27
  • 网络出版日期:  2020-04-25
  • 刊出日期:  2020-06-01

金属蜂窝夹芯结构抗水下冲击性能

    作者简介:李汶蔚(1985-),男,硕士,副研究员,主要从事爆轰物理研究. Email:wenweili@gmail.com
    通讯作者: 黄威, weihuang@hust.edu.cn
  • 1. 中国原子能科学研究院,北京 102413
  • 2. 华中科技大学船舶与海洋工程学院,湖北 武汉 430074

摘要: 为揭示高强度水下爆炸冲击载荷作用下金属夹芯结构的抗冲击性能,在实验室开展小尺寸水下爆炸加载技术对金属蜂窝夹芯结构性能影响的实验研究。基于实验结果,开展了全尺寸数值模拟金属蜂窝夹芯结构在水下冲击载荷作用下的动态响应和抗冲击性能研究。结果表明,数值模拟、实验和理论模型计算的结果具有良好的一致性。由于蜂窝芯材相对密度对夹芯结构能量耗散方式和载荷传递机制的影响,结构动态响应、失效模式以及抗冲击性能随着冲击强度的变化表现出较为明显的不同。通过抗冲击参数分析,建立了反映金属蜂窝夹芯结构抗冲击性能的结构横向变形、固支反力、透射脉冲和塑性能耗随冲击强度和芯材相对密度变化的结构-载荷-性能量化关系。

English Abstract

  • 在水下爆炸冲击载荷作用下,舰船局部防护结构的响应是一个高度非线性动力学过程,涉及流固耦合、断裂力学、塑性动力学等多个学科。舰船材料创新、优化船体结构及关键部位的抗冲击性能是目前船体防护性能研究的主要方向[1]。与单层板相比,夹芯结构具有轻质、比刚度大、比强度高、能量吸收性能和隔声隔热性能优越等特点[2],在冲击载荷作用下,其缓冲性和抗穿透性也具有明显优势[3-6]。相较于其他复合材料结构[7-8],拓扑芯材金属夹芯结构在准静态、抗侵彻以及爆炸性能研究等方面已取得了丰富的研究成果。近年来,金属夹芯结构也因其缓冲吸能优势广泛应用于航空、航海领域。

    为了进一步开展水下防爆结构研究,Fleck等[2]在Taylor[9]一维水下冲击波理论的基础上,设计了模拟远场水下爆炸的高强度水下冲击加载装置,通过该装置可产生与水下爆炸相似的加载脉冲信号。基于该加载装置,水下冲击载荷作用下金属夹芯结构的动态响应和失效机理研究获得了重大进展。Fleck等[2]将冲击载荷作用下的夹芯结构响应分为流固耦合阶段、芯材压缩阶段和结构响应阶段3个阶段,并预测了夹芯结构的变形模型。后续,研究人员通过一系列实验、数值模拟和理论分析对该结构动态响应模型进行修正[3, 10-12],通过对3个阶段的耦合和解耦,逐步深入揭示该结构的响应机理。国内学者Qin等[13]和Cui等[14]对空气爆炸下芯材强度对解耦响应的影响进行了研究;Huang等[15-17]采用不同水下爆炸模拟装置对多种芯材结构开展理论和实验研究。理论研究主要是分析芯材结构的变形机制,并且均基于解耦分析获得;实验研究则由于测试手段所限,很难获得完整的结构响应参数。因此,通过数值模拟手段进一步揭示水下冲击载荷作用下芯材结构的抗冲击特性十分必要。

    本研究在已完成的实验研究基础上,采用全尺寸三维数值模拟方法开展了金属蜂窝夹芯结构在高强度水下爆炸模拟装置加载下的动态响应和抗冲击特性研究。通过改变蜂窝夹芯结构的相对密度,研究了不同冲击载荷和芯材相对密度影响下夹芯结构的抗冲击特性,并与相应的实验和理论模型进行对比。

    • 通过数值模拟方法对文献[18]的实验结果进行深入分析,利用ABAQUS/EXPLICIT建立了如图1所示的基于拉格朗日方法的三维数值分析模型,模型所有尺寸与实验保持一致。高强度水下冲击加载模拟实验装置主要由图1(a)所示的一级轻气炮系统和圆桶状激波管组成,通过轻气炮发射高速飞片与激波管前置活塞平面撞击,在激波管中形成一系列平面冲击波,冲击波由激波管前端传播到尾端。这种可控性强、强度高的水下冲击波使得在实验室开展夹芯结构的抗水下冲击性能研究更为方便[19-20]。该加载装置的有效性已经被大量学者证实[8, 16, 19-20]

      图  1  高强度水下冲击加载实验装置示意图(a)和数值分析模型(b)

      Figure 1.  Intensive underwater explosive simulator (a) and the numerical model (b)

      该模型采用的飞片和前置活塞厚度分别为8.1 mm和12.0 mm,质量分别为0.22 kg和0.31 kg。实验过程中,飞片的初始速度vf为20~220 m/s,通过压力传感器测得水下脉冲峰值强度为10~300 MPa。靶板的有效加载面是一个直径(2L)为 66 mm的圆。金属蜂窝夹芯结构的蜂窝芯材为六边形蜂窝芯子,壁厚为0.01 mm,横向厚度为10 mm,材料为3003H18铝合金;前后面板厚度相同,均为0.5 mm,材料为5A06铝合金。考虑应变率效应的两种铝合金材料基于Johnson-Cook模型的材料性能参数见表1。加载应变率为7.41 × 10−4 s−1,横向屈服强度为1.92 MPa,压实应变为0.54。设计了4种不同相对密度的蜂窝夹芯结构进行数值模拟研究,通过将蜂窝芯子的单胞边长a分别设置为4、6和8 mm实现不同的结构相对密度,蜂窝单胞的其余尺寸一致。其中,通过多组重复实验得到的典型蜂窝夹芯结构(a = 4 mm)在准静态压缩下的应力-应变关系如图2所示。

      MaterialsYoung’s modulus/GPaDensity/(kg·m−3)Parameters
      A/MPaB/MPaCn
      5A06 aluminium alloy74.02 780167.0443.70.0200.44
      3003 aluminium alloy74.22 700 85.2170.00.0380.44

      表 1  面板和芯材材料力学性能参数

      Table 1.  Mechanical parameters of aluminum materials

      图  2  蜂窝夹芯结构准静态压缩下的应力-应变关系

      Figure 2.  Stress-strain relationship compression of honeycomb sandwich structure under quasi-static compression

      为了简化分析过程,忽略压力波传播过程中水容器及活塞可能与水发生的流固耦合现象,采用无摩擦的接触算法。在流体和结构接触处,由于压力降低在流固界面处易形成空化效应。为了模拟空化现象,水单元采用C3D8R单元。在水下冲击作用下表现为线弹性,并设置其对应的拉伸及剪切强度为零。为了不影响冲击波在水下的传播,水的横向单元尺寸必须受到严格控制。设置沿波为传播方向上水的单元尺寸为0.2 mm,径向尺寸为0.4 mm。将飞片、活塞和水容器均设置为刚体。在单层5A06铝合金靶板抗冲击实验中,靶板单元同样采用C3D8R单元。在受冲击区域环向单元尺寸为0.4 mm,沿厚度方向为3个单元。为了与实验一致,靶板背面用完全固定边界条件的刚性金属圆环对靶板进行固定。

      金属蜂窝结构面板和芯材均采用Johnson-Cook本构模型,模型参数见表1。二者的失效模型均采用延性失效模型,假设累计等效塑性应变达到临界值时材料开始失效。水的模型采用线性Hugoniot形式的Mie-Grüneisen状态方程,其参数见表2,其中水下声速由实验结果测得[21]$\gamma $为Grüneisen常数。

      Density/(kg·m−3)Sound speed in water/(m·s−1)γ
      1 0001 1060.05

      表 2  Mie-Grüneisen状态方程参数

      Table 2.  Parameters for the Mie-Grüneisen equation of state

    • 为了验证数值分析模型的有效性和准确性,对文献[18]中蜂窝单胞边长为4 mm的蜂窝夹芯结构实验结果进行对比。为了与实验结果保持一致,本研究采用无量纲冲量表示加载强度

      $ \overline I = \frac{{{I_{\rm{t}}}}}{{L\sqrt {{\rho _{\rm{f}}}{\sigma _{\rm{f}}}} }} $

      式中:It为透射脉冲强度,L为靶板半跨长度,ρfσf分别为面板的密度和强度。进行无量纲变形

      $\overline w = w/L$

      $\overline t = \frac{t}{{L\sqrt {{\rho _{\rm{f}}}{\sigma _{\rm{f}}}} }}$

      式中:w为横向变形大小,t为结构响应时间。

      图3展示了蜂窝夹芯结构在受到初始速度为48.0 m/s的飞片撞击时,产生无量纲冲量103${\overline I _{\rm{t}}}$ = 2.68时,实验和数值模拟得到的失效模式对比。由图3可知,数值模拟结果与实验结果具有较好的一致性,包括失效发生的位置和分布,以及涉及单胞屈曲的失效模式。随着冲击强度发生变化,芯材的压缩也随之变化,伴随结构响应的变化。

      图  3  实验和数值模拟的失效模式对比

      Figure 3.  Comparison of the deformation modes obtained from the simulation and experiment

      图4(a)为蜂窝结构在不同冲击强度下背板中点的响应时程曲线。随着冲击强度增加,中点响应速度明显增加,蜂窝结构发生更大的塑性变形。尽管数值模拟的响应速度比实验结果更大,但是二者的中点最大变形相近。在数值模拟分析中,假设水下压力为零即发生空化,结构发生横向变形导致空化出现,这两个阶段的解耦现象较为清晰,空化出现的位置均位于流固界面处。在实验加载过程中,由于流固耦合效应,流固界面处的空化演化更为复杂,空化出现的位置和时间均与数值模拟存在差异。这种简化的数值模拟分析是响应时间不同的主要原因。背板中点变形随冲击强度变化的总体趋势与实验结果有较好的一致性,并且与Fleck等[2]采用的外包法理论分析结果较为吻合,如图4(b)所示。因此,该数值模型能够较为准确地模拟蜂窝夹芯结构在水下冲击载荷作用下的动态响应。

      图  4  蜂窝夹芯结构背板中点变形时程曲线(a)和最大变形与冲击强度的关系(b)

      Figure 4.  Mid-point deflection history of honeycomb sandwish structure (a) and relationship between impulsive intensities and maximum deflection (b)

    • 图5为数值模拟得到的胞元边长为4 mm的蜂窝夹芯结构在水下冲击载荷作用下横向变形的动态过程以及对应前面板上塑性铰的运动。对比实验和模拟得到的蜂窝结构失效模式发现:冲击强度较小时不足以驱动前面板的塑性铰运动至终点,最大变形以平台值的形式显现在面板的中心区域,如图3所示的实验结果,对应图5中0.40 ms时的作用结果。芯材的中心区域整体横向移动,两侧芯材胞元由边界向中心发生胞元渐进形式的屈曲。增加冲击强度,前面板以初始速度继续运动,在0.44 ms时两侧塑性铰汇聚于中点并继续横向运动。动态塑性铰链的运动及失效模式与实验结果有较好的一致性[18]

      图  5  蜂窝夹芯结构在${10^3}{\overline I _{\rm{t}}}$ = 3.58作用下的等效应变分布和塑性铰运动

      Figure 5.  Dynamic deformation and propagation of plastic hinges of the honeycomb sandwich,${10^3}{\overline I _{\rm{t}}} $ = 3.58

      Tilbrook等[22]通过夹芯结构前后面板的速度响应过程反映不同芯材强度金属夹芯结构在受到冲击载荷作用时的变形机制,并将其按芯材的强度区域将金属夹芯结构的变形机制分为4类。图6所示为夹芯结构在${10^3}{\overline I _{\rm{t}}}$ = 4.28时前后面板中点的速度响应过程。在冲击的初始时刻,前面板获得速度v0,并在芯材和背板的作用下开始减速,与此同时背板在芯材作用下开始加速;在时刻teq二者获得共同速度;之后,二者同时开始减速运动直至速度为零。由此判定含两种面板结构的芯材为高强度芯材。结合Tilbrook等[22]的结果与本研究的数值模拟结果可以看出,在不高于图6所示的冲击强度时,前、后面板在芯材未发生完全压实时已经获得共同速度。

      图  6  蜂窝夹芯结构前后面板中点的速度响应过程

      Figure 6.  Velocity time history at the mid-point of front and back face-sheet of honeycomb sandwich structure

    • 选用面板相同、蜂窝单胞边长不同的夹芯结构,对比分析芯材相对密度对结构抗冲击性能的影响,单胞边长分别为4、6和8 mm。采用无量纲质量${\overline m _{\rm{c}}} = {m_{\rm{c}}}/m$${m_{\rm{c}}}$为芯材质量,m为结构总体质量)。图7为不同相对质量芯材的蜂窝夹芯结构受到相同冲击载荷作用时的等效塑性应变分布。当芯材的相对质量较低时,其失效机理主要为面板的拉伸和蜂窝单胞壁的屈曲;而当${m_{\rm{c}}}$ = 0.21时,蜂窝芯材的失效除了单胞壁的屈曲外还有纵向的拉伸失效。这种纵向拉伸由外侧向中心逐渐增大,当冲击强度不断增加,芯材首先在中心位置发生局部压实和完全压实的失效模式。芯材相对密度变化呈现出的不同的芯材失效模式说明了芯材密度对结构响应的重要作用,进一步可表现为其抗冲击性能的差异。

      图  7  相同冲击强度下不同芯材相对密度的蜂窝夹芯结构等效应力分布

      Figure 7.  Equivalent strain distribution of honeycomb sandwich under the same impact loading

    • 在实验结果的基础上,利用建立的有效数值模型,通过改变响应夹芯结构芯材的相对密度可以获取更多的有效参数以评估结构的抗冲击性能,进而通过得到的趋势获取该抗冲击参数与冲击强度及相关几何特性间的量化关系,最终通过优化参数获得性能更好的结构。对于受到爆炸载荷加载的结构,可主要通过背板的中点横向变形$\overline w $、结构固支端反力FR、透射脉冲强度${\overline I _{\rm{T}}}$以及塑性能耗${\overline E _{\rm{P}}}$等几个方面对金属夹芯结构在水下冲击载荷作用下的结构响应和能量耗散机理进行对比分析,从而评估其抗冲击性能。结构-载荷-性能的关系[8]可直观反映结构的抗水下冲击性能,即抗冲击参数Z ($\overline w $FR${\overline I _{\rm{T}}}$${\overline E _{\rm{P}}}$)、相对质量x (${\overline m _{\rm{c}}}$)和脉冲加载强度y (${\overline I _{\rm{t}}}$)之间的量化关系

      $Z = K{x^m}{y^n}$

      式中:Kmn分别为拟合参数。

      通过数值仿真分析,得到图8所示的不同相对密度蜂窝夹芯结构在冲击强度发生变化下的无量纲变形。与图4相似,随着冲击强度增加,结构变形程度增大。在相同冲击强度下,中点变形随着芯材相对密度的降低而增加。通过式(4)对结果拟合,得到数值模拟后的横向变形和冲击强度,结构参数间的量化关系

      图  8  蜂窝夹芯结构横向变形与冲击强度及芯材相对密度的关系

      Figure 8.  The relationship of transverse deflections for honeycomb sandwich structure load intensity and relative core density

      $\overline w = 0.02{ {{{\overline m }_{\rm{c}}}} ^{ - 0.36}}{\left( {{{10}^3}{{\overline I }_{\rm{t}}}} \right)^{0.90}}$

      利用反力计算得到夹芯结构在不同强度水下冲击载荷作用下透射到固支端的脉冲强度。如图9所示,结构固支端反力时程曲线的变化趋势与冲击强度和芯材相对密度无关,幅值随冲击强度的增加而增大,随芯材相对密度的增加而降低,在单一冲击强度或芯材相对密度影响下结构有明显的反力差异。因此,夹芯结构受到水下冲击载荷作用后,作用于固支端的反力受到芯材相对密度和冲击强度的影响较大,在评估结构的抗冲击过程中,需要将相连结构形成的反力作为重要因素进行考虑。与相同单位面积质量金属单层板的反力对比,二者呈现明显不同的趋势,单层板的反力峰值大且持续时间长。

      图  9  夹芯结构支座的反力时程曲线

      Figure 9.  Reaction forces time histories for sandwich panels

      图10为反力积分得到蜂窝夹芯结构的透射脉冲强度、冲击强度和芯材相对密度的关系,其变化趋势与反力趋势相同。芯材相对密度最小时无量纲透射脉冲最大,冲击强度最高时透射脉冲最小。对计算结果进行拟合,得到数值模拟下的横向变形和冲击强度、结构参数间的量化关系

      图  10  夹芯结构透射脉冲强度与冲击强度及几何特性的关系

      Figure 10.  Relationship of transmitted impulses to incident impulse intensity and geometries

      ${\overline I _{\rm{T}}} = 0.05{ {{{\overline m }_{\rm{c}}}} ^{ - 2.2}}{\left( {{{10}^3}{{\overline I }_{\rm{t}}}} \right)^{ - 0.3}}$

      分析夹芯结构面板和芯材在结构受到冲击载荷作用下发生变形和失效过程中塑性能的耗散比例,对于分析芯材和面板对结构的防护性能有着十分重要的作用。对夹芯结构而言,塑性能的无量纲形式为

      ${\overline E _{\rm{P}}} = \frac{{{E_{\rm{P}}}}}{{A{\sigma _{\rm{f}}}\left( {m_{\rm{z}}/{\rho _{\rm{f}}}} \right)}}$

      式中:Ep为直接获取的结构消耗的塑性能,A为加载区域面积,$m_{\rm{z}}$为结构整体的单位面积质量。根据数值模拟结果,结合不同结构组成部分的塑性能量耗散时程曲线,从式(7)可得到相应的量化关系。

      图11所示为冲击载荷作用下蜂窝夹芯结构各部分以及总体塑性能耗时程曲线。由图11可知,前面板和芯材的初始响应时间基本同步,背板的响应时间明显大于二者。前面板的初始运动使得芯材整体迅速发生压缩和屈曲,面板发生拉伸和弯曲。在经历相同初始增长速率后,芯材的塑性能耗增长速率明显高于前面板的增长速率,进而使得在较低冲击强度下,芯材的最终塑性能耗高于前面板的能耗。但是,随着冲击强度增加,芯材发生压实或者局部压实失效后,前面板的能耗增长速率又超过了芯材,并最终发生最大的塑性变形。对于拓扑结构和几何尺寸不同的夹芯结构,其前面板获得最大塑性变形的临界冲击强度各不相同。由模拟结果可知,芯材相对密度降低导致该临界冲击强度降低。由于前面板和芯材消耗了塑性能,背板发生的塑性变形始终最小。这种特点使得总体塑性能耗不像前述几种抗冲击参数随冲击强度和相对密度呈现单调变化,而是表现为:在相同冲击载荷作用下,不同芯材密度引起的单位塑性能耗变化微小;随着冲击强度增加,由于其耗能机理不同,相同结构的塑性能耗出现降低,但整体塑性能耗仍然随冲击强度增强而增加。拟合得到塑性能耗、冲击强度和相对密度的量化关系

      图  11  蜂窝夹芯结构各部分的塑性能耗时程曲线(a)及塑性能耗与冲击强度和芯材相对密度关系(b)

      Figure 11.  Plastic dissipation histories at different places of honeycomb sandwich structure (a) and the relationship of plastic dissipations to incident impulse intensity and relative core density (b)

      ${\overline E _{\rm{P}}} = 0.02{ {{{\overline m }_{\rm{c}}}} ^{0.18}}{\left( {{{10}^3}\overline I _{\rm{t}}} \right)^{2.24}}$

    • 对比实验与数值模拟结果和经典的金属夹芯结构动态响应理论,三者有较好的一致性。在结构发生完全失效前,主要由面板拉伸和芯材压缩引起了一系列失效模式。芯材的屈曲从边界到中心呈现渐进型,随着冲击强度增加,芯材首先在中心位置发生局部压实和完全压实的失效模式。获得了冲击强度、芯材相对密度对水下冲击载荷下蜂窝夹芯结构横向变形、反力、透射脉冲以及塑性变形能耗几种抗水下冲击参数的结构、载荷、性能的影响规律与量化关系。

参考文献 (22)

目录

    /

    返回文章
    返回