含随机填充孔圆形蜂窝结构的面内冲击性能*

何 强,马大为,张震东

(南京理工大学机械工程学院,江苏 南京210094)

含随机填充孔圆形蜂窝结构的面内冲击性能*

何 强,马大为,张震东

(南京理工大学机械工程学院,江苏 南京210094)

研究多孔材料细观结构与宏观力学性能之间的关系,建立具有固定相对密度的含随机固体填充孔的圆形蜂窝结构模型。在此模型的基础上具体讨论了不同孔洞填充比和冲击速度对圆形蜂窝结构变形模式、动态冲击平台应力以及能量吸收性能的影响。研究结果表明:填充孔在蜂窝变形过程中有局部牵制作用,蜂窝材料变形模式仍为准静态模式、过渡模式、动态模式;当变形模式为过渡模式或动态模式时,结构的平台应力与速度的平方成线性关系,存在明显的速度效应;高速冲击下,含固体填充孔的蜂窝结构单位质量吸收的能量高于规则蜂窝结构。研究结果可为蜂窝材料的研究和设计提供参考。

固体力学;变形模式;能量吸收;平台应力;填充孔;圆形蜂窝

蜂窝及多孔泡沫材料作为具有重量轻、比强度和比刚度高 、稳定性好 、隔热隔声性好等一系列优点的复合材料,已 经广泛 应用在 航空、航天 、建筑 、包 装 等 很 多 领 域[1-2]。 多 孔 材 料 的 力 学 性 能 除 了 与密度和基体材料特性有关外,还与胞元的排列方式、孔壁的缺失等细观结构密切相关。因此,研究多孔材料细观结构与宏观力学性能之间的关系,对于蜂窝及多孔材料的研究和设计都具有很重要的意义。

在蜂窝及多孔材料的实际生产过程中,由于孔壁的坍塌或发泡不完全导致材料中存在实体堆积的现象并不少见,开 孔泡沫 金属制 备中也 存在颗 粒堆 积 不 理 想 或 未 完 全 消 除 导 致 材 料 密 实 的 情 况[3-4],但是目前关于这方面缺陷对多孔材料力学性能影响的研究还很有限。O.Prakash等[5]对蜂窝材料的实验研究表明,部分孔的实体填充导致蜂窝局部强化,其弹性模量和应变强化也因此提高,同时结构的密实化应变减小。C.Chen等[6]用有限元方法研究了二维蜂窝结构 准静态单向和平面 静 水压加载下实 体 填充孔对结构弹性模量和屈服应力的影响,发现实体填充孔使蜂窝结构弹性模量略有提高,但对单向屈服强度和平面静 水屈服 强度影 响不大 。I.Jeon等[3]则通 过实验 研究发 现因泡 沫垮塌 而出现 实体堆 积的 闭孔泡沫铝弹性模量低于相同密度的无垮塌缺陷泡沫铝,材料的塑性垮塌应力则基本不受影响。H.Nakamoto等[7-8]研究了具有线性排列填充孔的六边形 蜂 窝 结 构 的 面 内 冲 击 性 能,并 将 其 得 出 的 结 论 与 实验数据进行了对比。

目前关于含固体填充孔蜂窝结构的研究主要集中在其静态及准静态力学性能方面,本文中利用有限元分析软件 Abaqus对冲击载荷下具有随机固体填充孔的圆形蜂窝结构变形过程进行数值模拟,并从蜂窝结构的变形模式、动态冲击平台应力及能量吸收能力3个方面研究随机固体填充孔对蜂窝结构动态力学性能的影响,研究结果可为蜂窝材料的研究和设计提供参考。

1 有限元模型

面内冲击载荷作用下的规则圆形蜂窝材料计算模型如图1所示。为了消除尺寸效应,本文所选的规则蜂窝材料结构与 Sun De-qiang等[9]提出的多 层 规 则 排布圆 形 蜂 窝模 型 相 似。 试 件由 具 有 相同 半径和壁厚的圆形蜂窝胞元方形堆列而成。其中单个胞元的半径R=3 mm,壁厚t=0.15mm,试件的尺寸为L1×L2=96 mm×96 mm。基体材料为金属铝,采用理想弹塑性模型,泊松比ν=0.3,杨氏模量E=69 GPa,密度ρs=2 700 kg/m3,则圆形蜂窝材 料 结 构 的 密 度ρ0=π/2(t/R)ρs=212 kg/m3。采用显示 动 力 学 软 件 Abaqus/Explicit对 蜂 窝材料进行动力学特性分析,计算中选用S4R单元(4节点减缩积分壳单元)。为了收敛和计算精度的需要,沿壳厚度方向定义5个高斯积分点。整个压缩过程中所有可能接触的表面均定义为自由接触表面,且无摩擦。采用与文献[9]相同的边界条件,即当顶端刚性板沿y方向冲击蜂窝试件时,底端刚性板固定,结构所有面外位移均被限制以保证其处于平面应变状态。上下刚性板与蜂窝试件的外表面均视为光滑,两者接触无摩擦。模型的面外(沿z方向)厚度b =10 mm。

在规则蜂窝中引入随机固体填充孔,假定结构的相对密度保持不变,定义孔洞填充比为:

图1 规则圆形蜂窝材料计算模型Fig.1 Calculation model of regular circular honeycombs

式中:m为蜂窝结构中的随机固体填充数,N 为蜂窝结构中的胞元总数。容易得到孔洞填充比为α时未填充胞孔的壁厚t满足下面的方程,即:

式 中 :s0=πR2为 单 个 胞 孔 的 面 积 ,ρ为 蜂 窝 结 构 相 对 密 度 ,即 :ρ=ρ0/ρs=π/2(t/R)。

为简化模型,固体填充孔的细观结构与未填充胞孔相同,将固体填充孔质量均匀分布于胞壁上,可以算得相应的固体填充孔胞元胞壁材料密度为:

式 中 :s=πR2-π(R-t)2,为 胞 壁 的 面 积 。 不 同 孔 洞 填 充 比α的 蜂 窝 结 构 与 未 填 充 胞 孔 壁 厚t及 胞 壁 材料密度ρe的对应关系如表1所示。

本文中所研究的二维蜂窝结构的相对密度为0.1,其刚度和强度远低于固体填充孔。因此,固体填充孔可通过改变胞壁材料杨氏模量和屈服强度为规则蜂窝结构胞壁材料相应值的1 000倍引入,与 C.Chen等[6]的有限元模 型 相同。本文中对同一孔洞填充比计算了3组不同随机样本的含随机固体填充孔圆形蜂窝结构。

表1 含随机固体填充孔蜂窝结构参数Table 1 Parameters of honeycombs with randomly arranged solid inclusions

2 计算结果与分析

2.1 变形模式

限于篇幅,图2仅给出了α=0,0.04时,不同冲击速度下圆形蜂窝结构的变形模式图。当冲击速度为3 m/s时,蜂窝结构冲击端和固定端的应力基本相同,变形模式均为准静态模式,在初始压缩变形阶段,整个蜂窝试件的变形是均匀的,各个胞元均具有横向扩张。紧接着,变形不再是均匀的,部分胞元开始交替平行地发生坍塌,很多胞元被4个椭圆坍塌胞元围住,结构总体变形呈现“分布崩塌带”模式,结构中的填充孔对崩塌带的分布有微弱影响,所有固体填充孔必然会被坍塌的椭圆形胞元包围。随着相对压缩量的增大,胞元逐层压缩直至蜂窝材料密实化。

图2 圆形蜂窝材料变形轮廓图Fig.2 Deformation profile of circular honeycomb

随 着 冲 击 速度的提高 (v=40 m/s),规则蜂窝结构变形由准静态模式转变为过渡模式。变形初始阶段,惯性效应使得蜂窝结构在靠近冲击端处出现变形局部化,结构冲击端应力处于应力应变曲线的压实应变侧,因此高于准静态模式下的压缩应力值。随着相对压缩量的增大,规则圆形蜂窝结构中部的胞元变形要早于蜂窝结构两边的胞元,使得蜂窝结构整体变形出现“V”形变形带,而固体填充孔在变形过程中具有一定的牵制效应[5],使含填充孔蜂窝结构靠近冲击端的变形局部化现象更加明显,且变形带向结构两侧扩展,呈现出倒置“W”形。

图3 圆形蜂窝结构的变形模式Fig.3 Deformation mode of the circular honeycomb

冲击速度进一步增加为100 m/s时,含填充孔的圆形蜂窝结构 的变形过程与规则蜂窝材料的 变 形模式类似,只产生一个呈“I”形状的局部变形带,对应着从冲击端到固定端的逐层压溃变形模式。与含随 机 固 体 填 充 孔 六 边 形 蜂 窝 结 构 的 动 态 压 缩 变 形 模 式 不 同[7],圆 形 蜂 窝 结 构 密 实 化 时,所 有 的 胞 元(除填充孔外)均被压实,并没有出现渗漏效应(当蜂窝结构压缩至密实化应变时,部分被固体填充孔包围的未完全坍塌的胞元将不再发生变形)。

研究不同填充比的蜂窝结构在不同冲击速度下的变形模式,结果表明随机填充孔圆形蜂窝材料的变形模式包括准静态模式、过渡模式及动态模式,即固体填充孔并没有改变圆形蜂窝结构y方向动态冲击时的变形模式。图2则给出了这几种典型变形模式的变形轮廓图。

图3则给出了不同填充比下蜂窝结构变形模式转变时的冲击速度值。填充孔的引入对临界速度有一定的影响,根据图3中的数据,得到含填充孔蜂窝结构由准静态变形模式向过渡模式转变的临界速度v1(c)和 由 过 渡 模 式 向 动 态 模 式 转 变 的 临 界 速 度v2(c)的 近 似 公 式:

2.2 冲击端平台应力的速度效应

作为表征多胞材料能量吸收能力大小的重要参数,平台应力σp定义为冲击端刚性板压缩应力从第1个 应力峰 值到试 件密实 化前最 大应变 所对应 的压缩 反力的 平均名 义应力[10],即定 义平台 应力为 :

式中:名义应力σ为刚性板作用在试件上的压缩反力F与初始横截面积A(L2×b)的比值,名义应变ε定 义 为 试 件 在y 方 向 的 压 缩 量 与 初 始 长 度L1之 比 ,εd为 蜂 窝 结 构 的 密 实 化 应 变 值 ,εy为 屈 服 应 变 ,本 文中取εy=0.02[9]。

分析可得,当冲击速度大于由准静 态 变 形 模 式向 过 渡 模 式 转 变的 临 界 速 度v(c)1时,平 台 应 力σp与v2呈线性关系,并且可以由以下公式拟合:

式中:a1和a2为拟合系数。

图4(a)给出了含填充孔的圆形蜂窝结构冲击端平台应力与速度参量v2之间的关系,图4(b)所示为局部放大图。从局部放大图中也能够看出,当冲击速度小于临界速度v(c)1,即蜂窝结构的变形仍处于准静态模式时,冲击 端平台 应力曲 线 表 现 出 明 显 的 非 线 性。S.R.Reid等[11]基 于 一 维 冲 击 波 理 论 建 立 了均 匀多孔材料动态冲击平台应力σp与冲击速度之间的关系,并指出当冲击速度大于一定值时,规则蜂窝结构的平台应力存在速度效应,平台应力值与速度的平方成线性关系。各蜂窝结构的拟合参数值如表2所示,本文中的拟合结果表明,当蜂窝结构中含有随机固体填充孔时这样的速度效应同样存在,但此时拟合参数与填充比的大小有关。从表中可知,随着填充比的增大,拟合参数a2也有所增大,这说明填充比的增加使得圆形蜂窝结构冲击端平台应力的速度效应更加明显。

图4 冲击速度对平台应力的影响Fig.4 Effect of impact velocity on the plateau stress of honeycombs

表2 平台应力拟合关系式中的系数Table 2 Coefficients in fitting equation of plateau stress

2.3 固体填充孔对冲击端平台应力的影响

图5给 出 了 低 速 冲 击 时(v=3 m/s)4 种 蜂 窝 结 构 冲 击 端 的 名 义 应 力 应 变 曲 线 。研 究 的 蜂 窝 结 构 模型相对密度保持不变,随着孔洞填充比的增大,变形胞元孔的壁厚明显变薄,圆形蜂窝结构冲击端的应力平台降低。此时蜂窝结构均处于准静态变形模式,规则蜂窝冲击端压缩应力曲线较为平缓,含固体填充孔的圆形蜂窝结构因为填充孔附近的崩塌变形以及变形局部化区域的宏观密度突然提高,使得蜂窝结构的应力明显变大,冲击端应力应变曲线出现了应力尖峰值,当固体填充孔与冲击面的速度一致时,应力幅值回落。另外,随着孔洞填充比的增大,蜂窝结构的密实化应变值变小,由于未出现渗漏现象,其密实化应变值减小的幅度小于含随机固体填充孔的六边形蜂窝结构[7]。

定义相对平台应力为含固体填充孔圆形蜂窝结构冲击端平台应力σp*与规则圆形蜂窝结构冲击端平台 应 力σp(0)之比。图6所示为蜂窝结构在y方向上受冲击时冲击端相对平台应力随冲击速度的变化关系。低速冲击下,随着孔洞填充比的增大,变形胞孔壁厚的减小使得蜂窝结构冲击端的承载能力显著下 降,如 图 中 所 示 ,冲 击 速 度 为3 m/s时 ,α=0.06 的 蜂 窝 结 构 平 台 应 力 约 为 相 同 相 对 密 度 的 规 则 蜂 窝结构平台应力的50%。

随着冲击速度的增大,蜂窝结构变形向冲击端附近移动,惯性效应引起的变形局部化使得结构平台应 力 增 加 。 对 于α=0.04,0.06的 蜂 窝 结 构,填 充 比 的 增 大 使 得 胞 元 壁 厚 变 薄,一 定 程 度 上 会 导 致 结 构压缩应力的下降,但是固体填充孔的牵制效应使得动态冲击的惯性效应影响增强,导致蜂窝结构冲击端应力应变曲线中出现应力尖峰,提高了蜂窝结构的吸能能力,这一增强效应很明显,因而在冲击速度小于40 m/s时,随 着 冲 击 速 度 的 提 高,相 对 平 台 应 力 会 显 著 增 强 。 当 冲 击 速 度 达 到40 m/s时,含 固 体 填充孔蜂窝结构的平台应力与规则蜂窝结构平台应力基本相同。继续加大蜂窝结构的冲击加载速度(v> 80 m/s),惯性效应的影响起主要作用,胞元的主要变形由低速时的交替平行胞元坍塌变为由接触 点 处开始向内凹陷崩塌的逐层压缩变形,且蜂窝结构的变形开始慢慢地集中在冲击端附近,这就使得规则蜂窝结构的平台应力也明显增强,从而削弱了相对平台应力的增 长,当速度达到100 m/s时候,相对 平 台应力甚至有下降的趋势。

图5 低速冲击下圆形蜂窝上材料冲击端动态响应曲线Fig.5 Dynamic response curve of circular honeycomb at the crushing end

图6 固体填充孔对y方向相对平台应力的影响Fig.6 Effect of solid inclusions on the relative plateau stress in y direction

2.4 随机填充孔对能量吸收特性的影响

作为一种高效的阻尼材料,多孔胞元材料已被广泛地应用于各种能量吸收结构中。评价多胞材料与结构能量吸收能力大小的重要指标就是单位质量吸收的能量大小[12]:

图7 不同冲击速度下蜂窝结构单位质量吸收能量Fig.7 Energy absorption of honeycombs structures at different impact velocities

由图7可见,低速冲击下(v=3 m/s)因为各蜂窝结构均处于准静态压缩模式,单位质量能量吸收曲线较为平滑,固体填充孔并没有提高蜂窝结构的能量吸收能力,相反,胞元孔壁厚度的减小降低了冲击端的应力值,导 致含填 充孔蜂窝 结 构 的 吸 能 能 力 下 降,当ε=0.8 时,α=0.06 的 蜂 窝 结构单位 质 量 吸收能量大小仅为规则圆形蜂窝结构的37.5%。当冲击速度为40 m/s时,蜂窝结构的单位质量吸收能量与低速冲击时相比均有所提高,固体填充孔的存在提高了蜂窝结构的吸能能力,减小了与规则蜂窝结构单位质量吸收能量大小之间的差距,当ε=0.8 时,α=0.06 的 蜂 窝 结 构 单 位 质 量 吸 收 能 量 大 小 为 规 则圆形蜂窝结构的72%。当冲击速度为80 m/s时,蜂窝结构都处于动态冲击模式,惯性效应起到了主要的作用,固体填充孔的牵制效应使得蜂窝结构的变形局部化更加明显,尽管含随机固体填充孔的蜂窝结构胞壁较薄,但是其冲击端应力值仍大于规则圆形蜂窝结构,单位质量吸能能力大大提升。由图7(b) ~(c)也能看出,含随机固体填充孔蜂窝结构冲击端应力应变曲线中出现的应力峰值使得单位质量吸收能量曲线体现了局部的凹凸性。

3 结 论

本文中讨论了随机固体填充孔对圆形蜂窝结构材料的面内冲击性能的影响。计算结果表明:

(1)含固体填充孔的圆形蜂窝结构与相同相对密度的规则圆形蜂窝结构具有相同的变形模式,即准静态模式、过渡模式及动态模式,但填充比对变形模式的临界速度有影响。

(2)当含随机固体填充孔的圆形蜂窝结构变形处于准静态模式时,冲击端平台应力曲线表现出明显的非线性,冲击速度大于由准静态变形 模 式 向过渡模式转 变 的 临 界速度v(c)1时,平台 应 力σp与v2呈 线性关系,存在明显的速度效应。

(3)低速冲击下,含固体填充孔的蜂窝结构平台应力随孔洞填充比的增大而显著降低,导致含填充孔蜂窝结构的吸能能力下降,冲 击速度 为3 m/s时,α=0.06 的蜂窝 结构平 台 应 力 约 为 相 同 相 对 密 度 的规则蜂窝结构平台应力的50%,单位质量吸收能量大小仅为规则圆形蜂窝结构的37.5%。随着冲击速度的提高,固体填充孔的牵制效应使得动态冲击的惯性效应影响增强,冲击端应力应变曲线出现明显的应力尖峰,提高了蜂窝结构的单位质量吸能能力。

本文中对含随机固体填充孔的圆形蜂窝结构的面内冲击性能的研究可为有缺陷蜂窝材料的研究提供参考,但关于随机固体填充孔与蜂窝结构动态响应特性更深层次的探讨还有待进一步展开。

[1]Gibson L J,Ashby M F.Cellular solids:structure and properties[M].2nd ed,Cambridge:Cambridge University Press,1997:93-147.

[2]卢 天 健,何 德 坪,陈 常 青,等.超 轻 多 孔 金 属 材 料 的 多 功 能 特 性 及 应 用[J].力 学 进 展 ,2006,36(4):517-535. Lu Tian-jian,He De-ping,Chen Chang-qing,et al.The multi-functionality of ultra-light porous metals and their applications[J].Advances in Mechanics,2006,36(4):517-535.

[3]Jeon I,Asahina T.The effect of structural defects on the compressive behavior of closed-cell A1 foam[J].Acta Materialia,2005,53:3415-3423.

[4]Kepets M,Lu T J,Dowling A P.Modeling of the role of defects in sintered FeCr Al Y foams[J].Acta Mechanica Sinica,2007,23(5):511-529.

[5]Prakash O,Bichebois P,Brechet Y,et al.A note on the deformation behaviour of two-dimensional model cellular structures[J].Philosophical Magazine A:Physics of Condensed Matter Structure Defects and Mechanical Properties,1996,73:739-751.

[6]Chen C,Lu T J,Fleck N A.Effect of inclusions and holes on the stiffness and strength of honeycombs[J].International Journal of Mechanical Sciences,2001,43(2):487-504.

[7]Nakamoto H,Adachi T,Araki W.In-plane impact behavior of honeycomb structures randomly filled with rigid inclusions[J].International Journal of Impact Engineering,2009,36(1):73-80.

[8]Nakamoto H,Adachi T,Araki W.In-plane impact behavior of honeycomb structures filled with linearly arranged inclusions[J].International Journal of Impact Engineering,2009,36(8):1019-1026.

[9]Sun De-qiang,Zhang Wei-hong.In-plane crushing and energy absorption performance of multi-layer regularly arranged circular honeycombs[J].Composite Structures,2013,96:726-735.

[10]Tan P J,Reid S R,Harrigan J J,et al.Dynamic compressive strength properties of aluminum foams:PartⅠ: Experimental data and observations[J].Journal of the Mechanics and Physics of Solids,2005,53(10):2174-2205.

[11]Reid S R,Peng C.Dynamic uniaxial crushing of wood[J].International Journal of Impact Engineering,1997,19 (5/6):531-570.

[12]Kooistra G W,Deshpande V S,Wadley H N G.Compressive behavior of age hardenable tetrahedral lattice truss structures made from aluminum[J].Acta Materialia,2004,52(14):4229-4237.

In-plane impact behavior of circular honeycomb structures randomly filled with rigid inclusions

He Qiang,Ma Da-wei,Zhang Zhen-dong
(School of Mechanical Engineering,Nanjing University of Science and Technology, Nanjing 210094,Jiangsu,China)

The model of circular honeycomb structures randomly filled with rigid inclusions which keeps the relative density as a constant is developed.And then the effects of impact velocity and packing ratio on the deformation modes,dynamic plateau stress and energy absorption capacities are discussed in detail.Research results show that the rigid inclusions have pinning effect in the process of deformation and the deformation modes can still be classified as quasi-static mode,transitional mode and dynamic mode.The plateau stress is proportional to the square of the impact velocity when the honeycombs are deformed at transitional mode or dynamic mode,which shows obvious speed effect. The energy absorption capacities of circular honeycombs are higher than these of the regular honeycombs at high-velocity impact.These results can provide valuable suggestions in the study and design of the functionally gradient honeycombs.

solid mechanics;deformation modes;energy absorption;plateau stress;inclusions;circular honeycomb

O347国标学科代码:13015

:A

10.11883/1001-1455-(2015)03-0401-08

(责任编辑 王易难)

2013-09-17;

2014-02-21

何 强(1989— ),男,博士研 究生,18260098162@163.com。