• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

负泊松比蜂窝材料抗爆炸特性及优化设计研究

孙晓旺 陶晓晓 王显会 李进军 王利辉

孙晓旺, 陶晓晓, 王显会, 李进军, 王利辉. 负泊松比蜂窝材料抗爆炸特性及优化设计研究[J]. 爆炸与冲击, 2020, 40(9): 095101. doi: 10.11883/bzycj-2020-0011
引用本文: 孙晓旺, 陶晓晓, 王显会, 李进军, 王利辉. 负泊松比蜂窝材料抗爆炸特性及优化设计研究[J]. 爆炸与冲击, 2020, 40(9): 095101. doi: 10.11883/bzycj-2020-0011
SUN Xiaowang, TAO Xiaoxiao, WANG Xianhui, LI Jinjun, WANG Lihui. Research on explosion-proof characteristics and optimization design of negative Poisson’s ratio honeycomb material[J]. Explosion And Shock Waves, 2020, 40(9): 095101. doi: 10.11883/bzycj-2020-0011
Citation: SUN Xiaowang, TAO Xiaoxiao, WANG Xianhui, LI Jinjun, WANG Lihui. Research on explosion-proof characteristics and optimization design of negative Poisson’s ratio honeycomb material[J]. Explosion And Shock Waves, 2020, 40(9): 095101. doi: 10.11883/bzycj-2020-0011

负泊松比蜂窝材料抗爆炸特性及优化设计研究

doi: 10.11883/bzycj-2020-0011
基金项目: 国家自然科学基金(11802140)
详细信息
    作者简介:

    孙晓旺(1987- ),男,博士,讲师,xwsun@njust.edu.cn

    通讯作者:

    王显会(1968- ),男,博士,教授,13770669850@139.com

  • 中图分类号: O389

Research on explosion-proof characteristics and optimization design of negative Poisson’s ratio honeycomb material

  • 摘要: 为了深入研究车辆底部防护组件爆炸冲击下的结构响应,提高防护型车辆的抗爆炸冲击性能,建立了某车辆底部防护组件在爆炸冲击下的有限元模型,并进行爆炸冲击台架试验验证了有限元模拟的可靠性;将内凹六边形负泊松比蜂窝材料作为防护组件的夹芯部分,分析负泊松比蜂窝材料在爆炸冲击下的变形模式,并对比了同等质量的其他3种防护组件的抗爆炸冲击性能。结果表明,含有负泊松比蜂窝夹芯的防护组件具有更优的抗爆性能。建立了以内凹六边形负泊松比蜂窝胞元尺寸参数为设计变量的多目标优化问题的数学模型,采用多目标遗传算法获得胞元几何参数的最优方案,有效降低了防护组件基板的最大挠度和最大动能。
  • 当特种车辆面对地雷、简易爆炸装置等威胁时,车辆底部防护性能的优劣将直接影响车内乘坐人员的安全性[1]。三明治防雷组件是一种有效的抗爆炸冲击结构,其由两块平板及中间的夹芯层组成。前面板能够将爆炸冲击载荷分配在可压溃的芯层上,吸收大量的能量减缓冲击波[2]。Li等[3]通过有限元模拟研究了空爆下蜂窝铝夹芯板结构的动态响应,发现夹芯层是吸收爆炸冲击能量的主要部分。

    负泊松比材料是一种新型的多胞材料,在受到轴向压缩时,材料在垂直于外力方向会产生收缩现象,在受到轴向拉伸时,材料在垂直于外力方向会产生膨胀现象[4],如图1所示。负泊松比蜂窝是一种特殊的多孔材料,由于其具有优异的剪切模量、断裂韧性等独特的性能,在许多领域有着广泛的应用[5]。Zhou等[6]对双箭头负泊松比结构的非线性几何特性进行了充分的研究,并采用HAM方法成功地实现了对负泊松比结构的优化设计;杨德庆等[7]研究了星形负泊松比超材料防护结构的抗爆抗冲击性能,并发现相比于普通防护结构,负泊松比蜂窝夹芯防护结构具有更加优良的水下抗爆性能;裴连政[8]通过数值模拟与试验相结合的方法,研究了含有负泊松比蜂窝的夹芯板在空中爆炸载荷下的动态响应和吸能特性,发现负泊松比夹芯结构具有良好的爆炸防护性能;Lan等[9]通过数值模拟比较了泡沫铝芯、六角形蜂窝芯及负泊松比蜂窝芯圆柱夹芯板在爆炸环境下的动态响应,研究发现含负泊松比蜂窝芯的圆柱面板具有更好的抗爆性能;Jin等[10]研究了爆炸冲击下具有功能梯度的负泊松比蜂窝芯层的结构响应,并与未分级的蜂窝芯层和规则排列的蜂窝芯层进行了对比,结果显示,分级蜂窝芯层及交叉排列蜂窝芯层能够显著提高夹芯层结构的抗爆性能。研究人员前期做了一些负泊松比材料的理论和应用分析,但还很少将负泊松比蜂窝材料应用到车辆底部防护组件进行相关的爆炸分析研究。

    图  1  负泊松比材料受载变形示意图
    Figure  1.  Schematic deformation of negative Poisson’s ratio material under load

    本文中通过有限元模拟及试验介绍某车辆底部防护组件的爆炸冲击台架。然后,基于负泊松比蜂窝材料优异的抗冲击性能及吸能效果,将内凹六边形负泊松比蜂窝材料作为三明治防护组件的夹芯层部分,分析其在爆炸冲击下的结构响应及吸能特性,并与其他3种防护组件进行对比。此外,基于遗传优化算法,建立以负泊松比胞元尺寸参数为设计变量的多目标优化模型进行优化设计。

    利用有限元软件建立某车辆底部防护组件的爆炸冲击台架模型,如图2所示,该模型中包括土壤、空气、炸药、台架、防护组件以及配重。台架整体尺寸为1 500 mm×1 860 mm×760 mm,其支撑结构采用Q235钢焊接构成,台架上方配重8 t。防护组件布置于配重支撑梁与台架支撑脚之间,防护组件最低点离地面高度为330 mm。2 kg柱形炸药埋于土下,炸点位于防护组件中心正下方,按照标准,炸药上表面距离土壤表面100 mm。

    图  2  爆炸冲击台架有限元模型
    Figure  2.  A finite element model for the explosive impact bench
    1.1.1   有限元模型

    防护组件主要由10 mm厚面板、8 mm厚背板及8 mm厚基板组成,材料分别为np500钢、np500钢和960E钢,具体材料参数如表1所示。背板与面板间距105 mm,背板与基板紧紧贴合,如图3所示。在面板与背板中部布置一根工字支撑梁,工字梁尺寸为1 000 mm×80 mm×100 mm,材料为KS700钢。台架各部分结构采用二维面网格进行模拟,配重以实体单元的形式进行模拟,网格单元基本尺寸为10 mm。防护组件结构之间采用*CONTACT_AUTOMATIC_SINGLE_SURFACE来定义接触,摩擦因数设定为0.2。空气与土壤均采用六面体实体网格进行模拟,采用*CONSRTAINED_LAGRANGE_IN_SOLID关键字来定义防护组件与空气和土壤域的流固耦合。用基板来模拟战术车辆中的底甲板,若底甲板受到较大冲击,必将传递到车身,从而影响车内乘员的安全性。因此,对该试验台架主要考察爆炸冲击后防护组件的破坏情况及基板的变形情况。

    表  1  防护组件各部分材料参数
    Table  1.  Material parameters for each part of the protective component
    材料密度ρ/(kg·m−3)杨氏模量E/GPa屈服强度σy/MPa泊松比µ抗拉强度σt/MPa
    np500钢7.8 × 10321013820.31757
    960E钢7.8 × 103210 9860.31150
    KS700钢7.8 × 103210 7000.3 752
    下载: 导出CSV 
    | 显示表格
    图  3  防护组件示意图
    Figure  3.  Schematic diagram of protection component

    基于多物质任意拉格朗日欧拉流固耦合算法(multi-material arbitrary Lagrangian-Eulerian fluid-structure interaction method,ALE-FSI)能够模拟炸药从引爆到对目标结构产生冲击的整个过程,其计算精度最高,是国内公认最成熟的爆炸数值算法[11]。因此,基于ALE-FSI爆炸冲击数值算法,利用LS-DYNA软件,对上述爆炸冲击台架进行数值模拟。

    1.1.2   爆炸后防护组件的结构响应

    数值模拟中,防护组件在受到爆炸冲击后,防护组件的面板、背板和基板都未出现结构损坏,3块钢板的结构响应如表2所示。爆炸产生后,面板首先受到爆炸冲击载荷作用,面板的最大挠度、最大加速度、最大动能和内能都远大于背板的。爆炸冲击波产生的能量经面板和背板吸收后作用在基板上,基板的最大挠度达到91.68 mm,图4为基板挠度云图,最大加速度为448.52g,爆炸过程中基板的最大动能高达16.81 kJ。

    表  2  防护组件的结构响应与能量
    Table  2.  Structural response and energy of protective component
    结构最大挠度d/mm最大加速度a/g最大动能Ek/kJ最大内能Ei/kJ
    面板132.821018.15123.2869.39
    背板 89.38 470.87 16.3531.72
    基板 91.68 448.52 16.8135.14
    下载: 导出CSV 
    | 显示表格
    图  4  基板挠度云图
    Figure  4.  Cloud diagram of substrate deflection

    为验证数值模型的准确性,布置了如图5所示的爆炸冲击试验台架。试验中台架结构、材料、配置及炸药布置等均与数值计算中保持一致。在台架基板上方中部安装应变梳,测量基板在爆炸过程中的动态变形,如图6所示。

    图  5  爆炸冲击试验台架
    Figure  5.  Explosive impact test bench
    图  6  应变梳布置
    Figure  6.  Strain comb arrangement

    试验后,台架整体未出现明显的结构损坏与焊缝失效等现象,计算结果与试验结果比较吻合。测量基板上方安装的应变梳的最大挠度约为89 mm,如图7所示,可以作为爆炸过程中基板的最大动态变形。计算中,基板中心处产生了较大弹性变形,最大挠度为91.96 mm,数值计算与试验相比误差为3.11%,误差在可接受范围内。综上所述,对比试验现象,计算精度能够满足要求,通过ALE-FSI算法来模拟爆炸冲击下防护组件的结构响应是合理的。

    图  7  试验后应变梳
    Figure  7.  Strain comb after test

    上述防护组件中,虽然面板与背板吸收了爆炸冲击产生的大量能量,但基板挠度仍然较大,基板的最大加速度高达448.52g,最大动能超过16 kJ,若将此防护组件布置到车辆底部、车身侧围和地板等关键位置必将受到较高冲击,从而增大车内乘员受伤的风险。

    为了提高防护组件的防护能力,在不改变其整体结构的情况下,在背板与面板之间添加负泊松比蜂窝夹芯材料,并在两侧增加夹芯固定件,如图8所示,内凹六边形负泊松比蜂窝材料局部结构示意图如图9所示。内凹六边形蜂窝与夹芯固定件的基体材料都为H14铝,材料参数如表3所示。计算中,防护组件各部分之间的接触类型及模型的边界条件保持不变。

    表  3  H14铝材料参数
    Table  3.  H14 aluminum material parameters
    材料密度ρ/(kg·m−3)杨氏模量E/GPa屈服强度度σy/MPa泊松比µ抗拉强度σt/MPa
    H14铝2.7×103701880.3271
    下载: 导出CSV 
    | 显示表格
    图  8  负泊松比蜂窝夹芯防护组件
    Figure  8.  Negative Poisson’s ratio honeycomb sandwich protection component
    图  9  负泊松比蜂窝夹芯材料局部结构
    Figure  9.  Local structure of honeycomb sandwich material with negative Poisson’s ratio

    内凹六边形负泊松比蜂窝材料胞元结构几何参数如图10所示。其中,水平胞壁的长度为L1,弯曲胞壁的长度为L2,水平胞壁与弯曲胞壁之间的夹角为胞壁夹角θ,弯曲胞壁与水平胞壁的壁厚相等,称为胞元壁厚tc

    图  10  胞元几何参数示意图
    Figure  10.  Schematic diagram of cell geometric parameters
    L12L2cosθ
    (1)

    图10中所示几何参数必须满足式(1)时,才能使胞元结构满足内凹特性且具有一定的可变行程[12]

    研究中,在防护组件背板与面板之间布置的内凹六边形蜂窝材料胞元结构几何参数为:水平胞壁长度L1=22.5mm,弯曲胞壁长度L2=17.32mm,胞壁厚度tc=0.32mm,胞元夹角θ=60;蜂窝在垂直于胞元方向上的延伸长度为1000 mm,在z向上布置胞元3层,蜂窝结构总高度为90 mm,芯层整体从中间位置被支撑梁分成两部分,总重量为15.03 kg,如图9所示。

    内凹六边形负泊松比蜂窝防护组件在爆炸冲击下的典型变形模式如图11所示,图12为负泊松比蜂窝芯层中心部分的变形图,结合图13所示的防护组件各部分结构的动能时程曲线,来更好地理解负泊松比蜂窝及整个防护组件在爆炸冲击下的结构响应。爆炸产生的球形冲击波在0.4 ms时首先作用在面板上,面板上动能也随之开始增大;在0.1 ms之后蜂窝芯层底部中心开始被压缩,蜂窝芯层动能开始增大;在爆炸发生后的1.5 ms,蜂窝芯层中心区域继续被压缩,由于蜂窝材料的负泊松比效应,中心断开处出现向内收缩的趋势,如图12所示,背板、基板开始产生变形,动能开始增大;随着爆炸冲击波继续扩散,在爆炸发生后的2.0 ms~3.5 ms,面板部分区域出现回弹,其动能开始减小,远离爆炸中心的蜂窝芯层边缘处也开始被压缩,两侧的蜂窝结构开始向内收缩,负泊松比蜂窝致密到中心加载区域,蜂窝芯层的动能也达到峰值,背板、面板的变形及动能继续增大;在爆炸发生后的4.0 ms,芯层边缘处的材料快速向内收缩,而蜂窝在载荷冲击方向被压缩量已达峰值且基本维持不变,蜂窝整体的动能也开始减小,与此同时,背板与基板的变形持续增大,传递到背板与面板的动能继续增加;直到爆炸发生后5.0 ms,蜂窝芯层在载荷冲击方向和垂直载荷方向的压缩量基本维持不变,至此蜂窝芯层完全被压实,在横向和纵向上皆无可变形空间,此时背板与基板的动能也达到峰值。

    图  11  爆炸冲击下负泊松比蜂窝防护组件典型变形模式
    Figure  11.  Typical deformation mode of Poisson’s ratio honeycomb protection component under explosion impact
    图  12  负泊松比蜂窝芯层中心变形图
    Figure  12.  Deformation diagram of negative Poisson’s ratio honeycomb core center
    图  13  防护组件各部分结构动能时程曲线
    Figure  13.  Time history curve of the kinetic energy of each part of the protective component

    含有负泊松比蜂窝夹芯层的防护组件在2 kg柱形炸药爆炸冲击下的结构响应与能量的变化,如表4所示。爆炸发生后,内凹六边形负泊松比蜂窝被压缩,材料拉入局部加载区域,更加适应动态载荷,与此同时夹芯层的内能持续增大,最大超过40 kJ,如图14所示。在爆炸冲击吸能结构设计中,比吸能(specific energy absorption,SEA)意为结构单位质量吸收的能量[13],是衡量抗爆性的重要指标,可用下式表示:

    表  4  防护组件的结构响应与能量
    Table  4.  Structural response and energy of protective components
    结构最大挠度d/mm最大加速度a/g最大动能Ek/kJ最大内能Ei/kJ比吸能η/(kJ·kg−1)
    面板130.4 1821.8048.2520.7310.095
    夹芯层 208.39 2.9046.69 3.107
    背板78.35 350.67 7.9810.99 0.063
    基板77.96 340.56 8.0411.48 0.078
    下载: 导出CSV 
    | 显示表格
    图  14  防护组件各部分结构内能时程曲线
    Figure  14.  Time history curve of the internal energy of each part of the protective components
    η=EtotalMtotal
    (2)

    式中:Etotal为结构整体吸收的总能量,Mtotal为结构的总质量。

    比吸能越大,说明该结构的吸能效率越高。由表4可知,内凹六边形负泊松比蜂窝夹芯层的吸能效率远大于防护组件其他结构的吸能效率。

    通过数值模拟提取的原始方案防护组件与含有负泊松比蜂窝夹芯材料的防护组件基板中心挠度及基板动能时程曲线如图1516所示,相比于原始防护组件,经过负泊松比蜂窝夹芯层防护组件吸能的基板,其最大挠度降低了13.72 mm,最大动能降低了52.17%。

    图  15  两种防护组件基板中心挠度时程曲线
    Figure  15.  Time history curve of the center deflection of the two kinds of protective component’s substrates
    图  16  两种防护组件基板动能时程曲线
    Figure  16.  Time history curve of the kinetic energy of the two kinds of protective component’s substrates

    通过对含有负泊松比蜂窝材料夹芯层的防护组件进行数值分析,研究发现,蜂窝夹芯层的吸能效率远高于防护组件中其他结构的吸能效率,且能够有效降低基板的挠度及动能。若将此类防护组件布置于车辆底部,车辆在面临爆炸威胁时,该防护组件能够起到很好的吸能效果,从而降低车身结构受到的冲击,进而降低车内乘员受到的损伤。

    为了验证内凹六边形负泊松比蜂窝材料在抗爆炸冲击台架中防护效果的优越性,在防护组件总质量相同的前提下进行了以下3组2 kg炸药台架爆炸的数值模拟。T1:面板厚度增大至10.69 mm;T2:背板厚度增大至8.69 mm;T3:将内凹六边形负泊松比蜂窝芯层替换为等质量的正六边形蜂窝芯层。

    其中,图17所示为正六边形蜂窝芯层部分结构,其整体结构设置为3层,整体尺寸与内凹六边形负泊松比蜂窝材料基本保持一致,基体材料为H14铝。正六边形蜂窝胞元结构几何尺寸如图18所示,正六边形边长L1=22.5mm,胞元夹角θ=120,胞元壁厚tc=0.67mm

    图  17  正六边形蜂窝芯层部分结构
    Figure  17.  Partial structure of regular hexagonal honeycomb core layer
    图  18  正六边形蜂窝胞元结构示意图
    Figure  18.  Schematic diagram of a regular hexagonal honeycomb cell structure

    图1920分别为数值模拟提取的4种防护组件基板中心的挠度时程曲线和基板动能时程曲线。其中,将内凹六边形负泊松比蜂窝芯层材料的质量附在面板和背板上的防护组件在受到爆炸冲击后,基板的最大挠度分别为88.42、87.28 mm,传递到基板的最大动能分别为12.44、12.48 kJ。正六边形蜂窝芯层防护组件在受到爆炸冲击后,基板的最大挠度为80.83 mm,基板的最大动能为9.38 kJ。在上述3种防护组件中,正六边形蜂窝芯层的防护组件的防护效果优于其他两种不加芯层的防护组件,但在同样工况同等质量情况下,与内凹六边形负泊松比蜂窝结构的防护组件相比,其基板的最大挠度及最大动能仍然明显更高。

    图  19  4种防护组件基板中心挠度时程曲线
    Figure  19.  Time history curves of center deflection of four kinds of protective component’s substrates
    图  20  4种防护组件基板动能时程曲线
    Figure  20.  Time history curves of kinetic energy of four kinds of protective component’s substrates

    韩会龙等[14]通过数值分析负泊松比蜂窝材料的动力学响应及能量吸收特性,发现改变胞元微结构可以提升负泊松比蜂窝的能量吸收能力。针对前文中的防护组件爆炸冲击台架,通过对内凹六边形负泊松比结构的胞元几何参数进行优化设计,提高负泊松比蜂窝结构的整体吸能效果,进而降低基板的最大挠度及最大动能。在确定设计变量、变量范围及优化目标的基础上,采用D-optimal试验设计对变量进行采样,然后根据计算得到的样本采用Kinging法构建代理模型,通过非支配排序遗传算法(NSGA-Ⅱ)来求解优化目标的Pareto解集,得到内凹六边形负泊松比蜂窝结构多目标优化的最优解。

    负泊松比胞元的水平胞壁的长度L1、弯曲胞壁的长度L2、胞元壁厚tc、胞壁夹角θ为设计变量,其中,L1L2tc为连续变量,θ为离散变量,且变量之间应满足L12L2cosθ。内凹六边形负泊松比蜂窝结构的优化设计目的是降低基板的最大挠度及最大动能,同时使蜂窝夹芯层的总质量限制在一定范围内。因此,根据设计目标、约束条件、设计变量确定了多目标优化数学模型如下式:

    minF(X)={Fd(x),FE(x)}s.t.{MM20.25mmL124.75mm15.58mmL219.09mm0.28mmtc0.36mmθ=40,45,50,55,60L12L2cosθ
    (3)

    式中:Fd(x)为基板的z向最大变形量,FE(x)为基板的最大动能,M为蜂窝夹芯的总质量,M为蜂窝夹芯的约束质量,设定为20 kg。在优化设计中,使得各目标函数最小,从而达到防护组件性能最好的要求。

    试验设计(是结构优化设计中重要的统计方法,主要用于选择合适的试验参数、确定最佳的参数组合、分析设计参数与响应结果之间的关系以及构建相应的代理模型来代替大计算量的实际模型进行求解等。本文中采用D-optimal试验设计方法来选取试验样点,为得到更加可靠的试验设计,要让选取的试验点可使模型的渐进协方差矩阵的行列式最小[15]。优化设计的设计变量有4个,采用D-optimal采样方法获得28组设计变量的组合形式及计算结果如表5所示。

    表  5  基于D-optimal采样的试验设计及结果
    Table  5.  Experimental design and results based on D-optimal sampling
    序号L1/mmL2/mmtc/mmθ/(°)M/kgd/mmEk/kJ
    121.1516.970.345518.6177.668.10
    222.9516.280.316013.0076.767.78
    2721.3316.560.305016.8178.838.59
    2823.1315.870.345514.7577.768.16
    下载: 导出CSV 
    | 显示表格

    Kringing代理模型是一种估计方差最小的无偏代理模型。该方法可以描述防护组件爆炸仿真这样的高度非线性动态过程[16]。利用有限元软件将每组采样点计算得出响应值,采用Kringing法构造代理模型后,一般采用决定系数R2来评估代理模型的拟合精度,其中决定系数表达式为:

    R2=ni=1(~yi¯yi)2ni=1(yi¯yi)2
    (4)

    式中:n为样本点数量,yi为第i个目标的实际响应值,~yi为预测目标的响应值,¯yi为平均的实际响应值。

    R2值越接近于1,代理模型精度越高[17]。由表6可知,代理模型各目标响应的决定系数R2都大于0.95,因此代理模型符合精度要求。

    表  6  目标响应的决定系数
    Table  6.  Decision coefficients of target response
    MFd(x)FE(x)
    R20.9630.9870.959
    下载: 导出CSV 
    | 显示表格

    在构建代理模型后,本文利用非支配排序遗传算法(NSGA-Ⅱ)来求解优化目标的Pareto解集。NSGA-Ⅱ算法运行速度快、计算复杂度低、易于实现,已成为多目标优化算法的基准之一,广泛应用于科学研究及工程实践领域[18]。在求解多目标优化问题时,通常不能获得唯一的最优解,得到的多个无法简单进行比较的解集即为帕累托(Pareto)最优解。在运用NSGA-Ⅱ算法求解时设置Pareto解集数为1 000个,遗传代数为50代,每代精英数量占样本空间的10%,遗传变异率为0.01。最终得到帕累托解集858个,帕累托前沿如图21所示。在车辆底部爆炸防护中,最重要的是降低车内乘员损伤,起决定作用的是传递到车身的冲击能量,因此本文优化结果更加注重降低基板的动能,从而筛选出第50代第156组为内凹六边形负泊松比蜂窝夹芯优化的最优解,如表7所示。最优解对应的水平胞壁长度为20.43 mm,弯曲胞壁长度为18.11 mm,胞壁厚度为0.32 mm,胞壁夹角为50°。

    表  7  第50代Pareto解集(部分)及变量参数
    Table  7.  The 50th generation Pareto solution set (part) and variable parameters
    解集序号L1/mmL2/mmtc/mmθ/(°)M/kgd/mmEk/kJ
    124.3515.780.336012.4174.567.99
    222.3615.780.325012.4776.647.84
    15524.4415.790.315512.5375.837.70
    15620.4318.110.325018.1876.236.52
    85721.2717.850.326017.8976.166.96
    85820.7817.440.345016.2376.176.99
    下载: 导出CSV 
    | 显示表格
    图  21  帕累托前沿
    Figure  21.  Pareto front

    为了验证内凹六边形负泊松比蜂窝夹芯材料多目标优化结果的准确性,根据3.4节优化结果得到的设计变量参数,建立爆炸冲击台架有限元模型,求解后与优化结果对比。基板最大挠度的优化结果和数值模拟结果分别为76.23和74.58 mm,相对误差为2.16%;基板最大动能分别为6.52和6.41 kJ,相对误差为1.69%,误差皆小于5%,因此可以认为优化结果具有较高的准确性。如图2223所示,优化后新方案的基板最大挠度为74.58 mm,相比于优化前降低了4.34%;基板的最大动能为6.41 kJ,相比于优化前降低了20.27%。

    图  22  优化前后基板挠度时程曲线
    Figure  22.  Time history curves of substrate deflection before and after optimization
    图  23  优化前后基板动能时程曲线
    Figure  23.  Time history curves of substrate kinetic energy before and after optimization

    本文基于车辆底部防护组件爆炸冲击台架,提出了一种将内凹六边形负泊松比蜂窝材料作为其夹芯层的防护组件,分析了防护组件在2 kg柱形炸药爆炸冲击下的结构响应;通过NSGA-Ⅱ遗传算法求解基于负泊松比胞元尺寸参数的多目标优化问题,得到基板最大挠度与基板最大动能的Pareto前沿并选出最优解。研究表明:(1)相比于不含夹芯结构的防护组件,内凹六边形负泊松比蜂窝防护组件能够有效降低基板的最大挠度及动能,且负泊松比蜂窝夹芯在整个防护组件中的比吸能远大于面板和背板。(2)同等质量的前提下,相比于正六边形蜂窝夹芯等其他3种防护组件,内凹六边形负泊松比蜂窝夹芯防护组件的防护效果更优。(3)在满足优化要求的前提下,成功实现了对负泊松比夹芯层结构的多目标优化,优化后基板的最大挠度降低了4.34%,基板的最大动能降低了20.27%。

  • 图  1  负泊松比材料受载变形示意图

    Figure  1.  Schematic deformation of negative Poisson’s ratio material under load

    图  2  爆炸冲击台架有限元模型

    Figure  2.  A finite element model for the explosive impact bench

    图  3  防护组件示意图

    Figure  3.  Schematic diagram of protection component

    图  4  基板挠度云图

    Figure  4.  Cloud diagram of substrate deflection

    图  5  爆炸冲击试验台架

    Figure  5.  Explosive impact test bench

    图  6  应变梳布置

    Figure  6.  Strain comb arrangement

    图  7  试验后应变梳

    Figure  7.  Strain comb after test

    图  8  负泊松比蜂窝夹芯防护组件

    Figure  8.  Negative Poisson’s ratio honeycomb sandwich protection component

    图  9  负泊松比蜂窝夹芯材料局部结构

    Figure  9.  Local structure of honeycomb sandwich material with negative Poisson’s ratio

    图  10  胞元几何参数示意图

    Figure  10.  Schematic diagram of cell geometric parameters

    图  11  爆炸冲击下负泊松比蜂窝防护组件典型变形模式

    Figure  11.  Typical deformation mode of Poisson’s ratio honeycomb protection component under explosion impact

    图  12  负泊松比蜂窝芯层中心变形图

    Figure  12.  Deformation diagram of negative Poisson’s ratio honeycomb core center

    图  13  防护组件各部分结构动能时程曲线

    Figure  13.  Time history curve of the kinetic energy of each part of the protective component

    图  14  防护组件各部分结构内能时程曲线

    Figure  14.  Time history curve of the internal energy of each part of the protective components

    图  15  两种防护组件基板中心挠度时程曲线

    Figure  15.  Time history curve of the center deflection of the two kinds of protective component’s substrates

    图  16  两种防护组件基板动能时程曲线

    Figure  16.  Time history curve of the kinetic energy of the two kinds of protective component’s substrates

    图  17  正六边形蜂窝芯层部分结构

    Figure  17.  Partial structure of regular hexagonal honeycomb core layer

    图  18  正六边形蜂窝胞元结构示意图

    Figure  18.  Schematic diagram of a regular hexagonal honeycomb cell structure

    图  19  4种防护组件基板中心挠度时程曲线

    Figure  19.  Time history curves of center deflection of four kinds of protective component’s substrates

    图  20  4种防护组件基板动能时程曲线

    Figure  20.  Time history curves of kinetic energy of four kinds of protective component’s substrates

    图  21  帕累托前沿

    Figure  21.  Pareto front

    图  22  优化前后基板挠度时程曲线

    Figure  22.  Time history curves of substrate deflection before and after optimization

    图  23  优化前后基板动能时程曲线

    Figure  23.  Time history curves of substrate kinetic energy before and after optimization

    表  1  防护组件各部分材料参数

    Table  1.   Material parameters for each part of the protective component

    材料密度ρ/(kg·m−3)杨氏模量E/GPa屈服强度σy/MPa泊松比µ抗拉强度σt/MPa
    np500钢7.8 × 10321013820.31757
    960E钢7.8 × 103210 9860.31150
    KS700钢7.8 × 103210 7000.3 752
    下载: 导出CSV

    表  2  防护组件的结构响应与能量

    Table  2.   Structural response and energy of protective component

    结构最大挠度d/mm最大加速度a/g最大动能Ek/kJ最大内能Ei/kJ
    面板132.821018.15123.2869.39
    背板 89.38 470.87 16.3531.72
    基板 91.68 448.52 16.8135.14
    下载: 导出CSV

    表  3  H14铝材料参数

    Table  3.   H14 aluminum material parameters

    材料密度ρ/(kg·m−3)杨氏模量E/GPa屈服强度度σy/MPa泊松比µ抗拉强度σt/MPa
    H14铝2.7×103701880.3271
    下载: 导出CSV

    表  4  防护组件的结构响应与能量

    Table  4.   Structural response and energy of protective components

    结构最大挠度d/mm最大加速度a/g最大动能Ek/kJ最大内能Ei/kJ比吸能η/(kJ·kg−1)
    面板130.4 1821.8048.2520.7310.095
    夹芯层 208.39 2.9046.69 3.107
    背板78.35 350.67 7.9810.99 0.063
    基板77.96 340.56 8.0411.48 0.078
    下载: 导出CSV

    表  5  基于D-optimal采样的试验设计及结果

    Table  5.   Experimental design and results based on D-optimal sampling

    序号L1/mmL2/mmtc/mmθ/(°)M/kgd/mmEk/kJ
    121.1516.970.345518.6177.668.10
    222.9516.280.316013.0076.767.78
    2721.3316.560.305016.8178.838.59
    2823.1315.870.345514.7577.768.16
    下载: 导出CSV

    表  6  目标响应的决定系数

    Table  6.   Decision coefficients of target response

    MFd(x)FE(x)
    R20.9630.9870.959
    下载: 导出CSV

    表  7  第50代Pareto解集(部分)及变量参数

    Table  7.   The 50th generation Pareto solution set (part) and variable parameters

    解集序号L1/mmL2/mmtc/mmθ/(°)M/kgd/mmEk/kJ
    124.3515.780.336012.4174.567.99
    222.3615.780.325012.4776.647.84
    15524.4415.790.315512.5375.837.70
    15620.4318.110.325018.1876.236.52
    85721.2717.850.326017.8976.166.96
    85820.7817.440.345016.2376.176.99
    下载: 导出CSV
  • [1] 李补莲, 原树兴. 现代军用车辆的防护措施 [J]. 国外坦克, 2011(5): 41–43. DOI: CNKI:SUN:GWTK.0.2011-05-019.

    LI B L, YUAN S X. Protective measures of modern military vehicles [J]. Foreign Tanks, 2011(5): 41–43. DOI: CNKI:SUN:GWTK.0.2011-05-019.
    [2] IMBALZANO G, LINFORTH S, NGO T D, et al. Blast resistance of auxetic and honeycomb sandwich panels: comparisons and parametric designs [J]. Composite Structures, 2018, 183: 242–261. DOI: 10.1016/j.compstruct.2017.03.018.
    [3] LI X, ZHANG P, WANG Z, et al. Dynamic behavior of aluminum honeycomb sandwich panels under air blast: experiment and numerical analysis [J]. Composite Structures, 2014, 108: 1001–1008. DOI: 10.1016/j.compstruct.2013.10.034.
    [4] 周冠. 新型负泊松比结构关键技术研究及其在车身设计中的应用[D]. 长沙: 湖南大学, 2015: 6−7.
    [5] QIAO J, CHEN C Q. Analyses on the in-plane impact resistance of auxetic double arrowhead honeycombs [J]. Journal of Applied Mechanics, 2015, 82(5): 051007. DOI: 10.1115/1.4030007.
    [6] ZHOU G, MA Z D, GU J, et al. Design optimization of a NPR structure based on HAM optimization method [J]. Structural and Multidisciplinary Optimization, 2016, 53(3): 635–643. DOI: 10.1007/s00158-015-1341-x.
    [7] 杨德庆, 吴秉鸿, 张相闻. 星型负泊松比超材料防护结构抗爆抗冲击性能研究 [J]. 爆炸与冲击, 2019, 39(6): 124–135. DOI: 10.11883/bjycj-2018-0060.

    YANG D Q, WU B H, ZHANG X W. Explosion and impact resistance of star negative Poisson’s ratio metamaterial protective structure [J]. Explosion and Shock Waves, 2019, 39(6): 124–135. DOI: 10.11883/bjycj-2018-0060.
    [8] 裴连政. 负泊松比夹芯板抗爆性能实验与仿真研究[D]. 大连: 大连理工大学, 2016: 11−34.
    [9] LAN X K, FENG S S, HUANG Q, et al. A comparative study of blast resistance of cylindrical sandwich panels with aluminum foam and auxetic honeycomb cores [J]. Aerospace Science and Technology, 2019, 87: 37–47. DOI: 10.1016/j.ast.2019.01.031.
    [10] JIN X C, WANG Z H, NING J G, et al. Dynamic response of sandwich structures with graded auxetic honeycomb cores under blast loading [J]. Composites Part B: Engineering, 2016, 106: 206–217. DOI: 10.1016/j.compositesb.2016.09.037.
    [11] 魏然. 爆炸冲击下车身结构防护机理及多学科优化研究[D]. 南京: 南京理工大学, 2019: 3−4.
    [12] 王陶. 负泊松比结构力学特性研究及其在商用车耐撞性优化设计中的应用[D]. 南京: 南京理工大学, 2018: 21−22.
    [13] ZHOU G, MA Z D, LI G, et al. Design optimization of a novel NPR crash box based on multi-objective genetic algorithm [J]. Structural and Multidisciplinary Optimization, 2016, 54(3): 673–684. DOI: 10.1007/s00158-016-1452-z.
    [14] 韩会龙, 张新春, 王鹏. 负泊松比蜂窝材料的动力学响应及能量吸收特性 [J]. 爆炸与冲击, 2019, 39(1): 47–57. DOI: 10.11883/bzycj-2017-0281.

    HAN H L, ZHANG X C, WANG P. Dynamic response and energy absorption characteristics of negative Poisson’s ratio honeycomb materials [J]. Explosion and Shock Waves, 2019, 39(1): 47–57. DOI: 10.11883/bzycj-2017-0281.
    [15] 田永军, 段国林, 夏晓光, 等. 响应面模型与混合优化算法相结合的锯片参数优化设计 [J]. 中国机械工程, 2016, 27(22): 3025–3031. DOI: 10.3969/j.issn.1004-132X.2016.22.008.

    TIAN Y J, DUAN G L, XIA X G, et al. Optimization design of saw blade parameters based on response surface model and hybrid optimization algorithm [J]. China Mechanical Engineering, 2016, 27(22): 3025–3031. DOI: 10.3969/j.issn.1004-132X.2016.22.008.
    [16] SIMPSON T W, MAUERY T M, KORTE J J, et al. Kriging models for global approximation in simulation-based multidisciplinary design optimization [J]. AIAA Journal, 2001, 39(12): 2233–2241. DOI: 10.2514/2.1234.
    [17] 任朋飞. 抗冲击型车身柔性底部结构防护技术研究[D]. 南京: 南京理工大学, 2018: 54−55.
    [18] 詹长书, 王清. 基于改进遗传算法电动汽车变速器参数设计与优化 [J]. 重庆理工大学学报(自然科学版), 2020, 34(2): 1–5. DOI: 10.3969/j.issn.1674-8425(z).2020.02.001.

    ZHAN C S, WANG Q. Design and optimization of transmission parameters of electric vehicle based on improved genetic algotithm [J]. Journal of Chongqing Institute of Technology (Natural Science), 2020, 34(2): 1–5. DOI: 10.3969/j.issn.1674-8425(z).2020.02.001.
  • 期刊类型引用(14)

    1. 沈鲁豫,温垚珂,董方栋,覃彬,徐浩然. TPU材料负泊松比结构缓冲性能研究. 兵器装备工程学报. 2024(02): 291-299 . 百度学术
    2. 贺锋,梁一鸣,秦海,李季,高超,刘三丰. 负泊松比结构的工程应用现状及展望. 防护工程. 2024(01): 68-78 . 百度学术
    3. 蒋舟顺,徐峰祥,邹震,周谦谋. 爆炸载荷下正弦曲边三维负泊松比夹芯板的动态响应和吸能特性. 爆炸与冲击. 2024(02): 3-20 . 本站查看
    4. 张旭,王子沱,李伟. 基于静态凝聚法的蜂窝芯轴向变胞元非均匀排布优化. 应用力学学报. 2024(03): 557-566 . 百度学术
    5. 贺锋,梁一鸣,李季,高超,慕孟,白家琦,祁少博. 异型陶瓷负泊松比复合结构的防弹防爆性能数值模拟研究. 安全与环境学报. 2024(08): 2919-2928 . 百度学术
    6. 孔祥清,李若男,常雅慧,付莹. 泡沫填充负泊松比蜂窝夹层结构的抗爆性能数值模拟. 兵工学报. 2024(09): 3091-3104 . 百度学术
    7. 时圣波,王韧之,唐佳宾,甘云丹,袁建飞,陈勇. 复合点阵结构强爆炸冲击载荷下的损伤机理与动态响应特性. 爆炸与冲击. 2023(06): 39-53 . 本站查看
    8. 崔飞蝶,董焱章,李青峰. 多型负泊松比材料-结构协同设计. 湖北汽车工业学院学报. 2023(04): 64-69+75 . 百度学术
    9. 赵雨薇,张宏伟,孙晓旺,王显会. 爆炸载荷下负泊松比下肢保护装置参数分析. 兵器装备工程学报. 2022(03): 60-66+106 . 百度学术
    10. 孙魁远,孙晓旺,张宏伟,王显会,彭兵,张绍彦. 厚度梯度型负泊松比蜂窝抗爆炸特性及优化. 兵器装备工程学报. 2022(04): 190-197 . 百度学术
    11. 赵颖,王凯锋,施劲余,桑叶,李云伍,陈宇,殷举伞. 微元胞缺失下三星型微结构面内动态性能分析. 振动与冲击. 2022(10): 115-123+185 . 百度学术
    12. 魏子涵,赵振宇,叶帆,裴轶群,王昕,张钱城,卢天健. 金属蜂窝夹层结构抗水下爆炸特性. 爆炸与冲击. 2021(08): 61-77 . 本站查看
    13. 沙康康,王显会,郭峤,彭兵,孙晓旺. 某型扫雷辊在爆炸环境下的仿真与试验研究. 科学技术与工程. 2021(25): 10576-10581 . 百度学术
    14. 乔及森,李明,苗红丽. 串联梯度蜂窝结构的面内力学性能. 塑性工程学报. 2021(11): 115-123 . 百度学术

    其他类型引用(25)

  • 加载中
图(23) / 表(7)
计量
  • 文章访问数:  3471
  • HTML全文浏览量:  1704
  • PDF下载量:  125
  • 被引次数: 39
出版历程
  • 收稿日期:  2020-01-13
  • 修回日期:  2020-06-05
  • 刊出日期:  2020-09-01

目录

/

返回文章
返回