In-plane impact behavior of circular honeycomb structures randomly filled with rigid inclusions
-
摘要: 研究多孔材料细观结构与宏观力学性能之间的关系, 建立具有固定相对密度的含随机固体填充孔的圆形蜂窝结构模型。在此模型的基础上具体讨论了不同孔洞填充比和冲击速度对圆形蜂窝结构变形模式、动态冲击平台应力以及能量吸收性能的影响。研究结果表明:填充孔在蜂窝变形过程中有局部牵制作用, 蜂窝材料变形模式仍为准静态模式、过渡模式、动态模式; 当变形模式为过渡模式或动态模式时, 结构的平台应力与速度的平方成线性关系, 存在明显的速度效应; 高速冲击下, 含固体填充孔的蜂窝结构单位质量吸收的能量高于规则蜂窝结构。研究结果可为蜂窝材料的研究和设计提供参考。Abstract: 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.
-
Key words:
- solid mechanics /
- deformation modes /
- energy absorption /
- plateau stress /
- inclusions /
- circular honeycomb
-
蜂窝及多孔泡沫材料作为具有重量轻、比强度和比刚度高、稳定性好、隔热隔声性好等一系列优点的复合材料, 已经广泛应用在航空、航天、建筑、包装等很多领域[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。
在规则蜂窝中引入随机固体填充孔, 假定结构的相对密度保持不变, 定义孔洞填充比为:
α=m/N (1) 式中:m为蜂窝结构中的随机固体填充数, N为蜂窝结构中的胞元总数。容易得到孔洞填充比为α时未填充胞孔的壁厚t满足下面的方程, 即:
ρL1L2−ms0=(N−m)[πR2−π(R−t)2] (2) 式中:s0=πR2为单个胞孔的面积, ρ为蜂窝结构相对密度, 即:ρ=ρ0/ρs=π/2(t/R)。
为简化模型, 固体填充孔的细观结构与未填充胞孔相同, 将固体填充孔质量均匀分布于胞壁上, 可以算得相应的固体填充孔胞元胞壁材料密度为:
ρe=ρss/s0=R2t(2R−t)ρs (3) 式中:s=πR2-π(R-t)2, 为胞壁的面积。不同孔洞填充比α的蜂窝结构与未填充胞孔壁厚t及胞壁材料密度ρe的对应关系如表 1所示。
表 1 含随机固体填充孔蜂窝结构参数Table 1. Parameters of honeycombs with randomly arranged solid inclusionsα t/mm ρe/(g·cm-3) 0.02 0.125 33.106 0.04 0.096 42.650 0.06 0.066 61.490 0.08 0.036 114.130 本文中所研究的二维蜂窝结构的相对密度为0.1, 其刚度和强度远低于固体填充孔。因此, 固体填充孔可通过改变胞壁材料杨氏模量和屈服强度为规则蜂窝结构胞壁材料相应值的1 000倍引入, 与C.Chen等[6]的有限元模型相同。本文中对同一孔洞填充比计算了3组不同随机样本的含随机固体填充孔圆形蜂窝结构。
2. 计算结果与分析
2.1 变形模式
限于篇幅, 图 2仅给出了α=0, 0.04时, 不同冲击速度下圆形蜂窝结构的变形模式图。当冲击速度为3 m/s时, 蜂窝结构冲击端和固定端的应力基本相同, 变形模式均为准静态模式, 在初始压缩变形阶段, 整个蜂窝试件的变形是均匀的, 各个胞元均具有横向扩张。紧接着, 变形不再是均匀的, 部分胞元开始交替平行地发生坍塌, 很多胞元被4个椭圆坍塌胞元围住, 结构总体变形呈现“分布崩塌带”模式, 结构中的填充孔对崩塌带的分布有微弱影响, 所有固体填充孔必然会被坍塌的椭圆形胞元包围。随着相对压缩量的增大, 胞元逐层压缩直至蜂窝材料密实化。
随着冲击速度的提高(v=40 m/s), 规则蜂窝结构变形由准静态模式转变为过渡模式。变形初始阶段, 惯性效应使得蜂窝结构在靠近冲击端处出现变形局部化, 结构冲击端应力处于应力应变曲线的压实应变侧, 因此高于准静态模式下的压缩应力值。随着相对压缩量的增大, 规则圆形蜂窝结构中部的胞元变形要早于蜂窝结构两边的胞元, 使得蜂窝结构整体变形出现“V”形变形带, 而固体填充孔在变形过程中具有一定的牵制效应[5], 使含填充孔蜂窝结构靠近冲击端的变形局部化现象更加明显, 且变形带向结构两侧扩展, 呈现出倒置“W”形。
冲击速度进一步增加为100 m/s时, 含填充孔的圆形蜂窝结构的变形过程与规则蜂窝材料的变形模式类似, 只产生一个呈“I”形状的局部变形带, 对应着从冲击端到固定端的逐层压溃变形模式。与含随机固体填充孔六边形蜂窝结构的动态压缩变形模式不同[7], 圆形蜂窝结构密实化时, 所有的胞元(除填充孔外)均被压实, 并没有出现渗漏效应(当蜂窝结构压缩至密实化应变时, 部分被固体填充孔包围的未完全坍塌的胞元将不再发生变形)。
研究不同填充比的蜂窝结构在不同冲击速度下的变形模式, 结果表明随机填充孔圆形蜂窝材料的变形模式包括准静态模式、过渡模式及动态模式, 即固体填充孔并没有改变圆形蜂窝结构y方向动态冲击时的变形模式。图 2则给出了这几种典型变形模式的变形轮廓图。
图 3则给出了不同填充比下蜂窝结构变形模式转变时的冲击速度值。填充孔的引入对临界速度有一定的影响, 根据图 3中的数据, 得到含填充孔蜂窝结构由准静态变形模式向过渡模式转变的临界速度v1(c)和由过渡模式向动态模式转变的临界速度v2(c)的近似公式:
v(c)1=−350α+38 (4) v(c)2=−450α+77 2.2 冲击端平台应力的速度效应
作为表征多胞材料能量吸收能力大小的重要参数, 平台应力σp定义为冲击端刚性板压缩应力从第1个应力峰值到试件密实化前最大应变所对应的压缩反力的平均名义应力[10], 即定义平台应力为:
σp=1εd−εy∫εdεyσdε (6) 式中:名义应力σ为刚性板作用在试件上的压缩反力F与初始横截面积A(L2×b)的比值, 名义应变ε定义为试件在y方向的压缩量与初始长度L1之比, εd为蜂窝结构的密实化应变值, εy为屈服应变, 本文中取εy=0.02[9]。
分析可得, 当冲击速度大于由准静态变形模式向过渡模式转变的临界速度v1(c)时, 平台应力σp与v2呈线性关系, 并且可以由以下公式拟合:
σp=a1+a2v2 (7) 式中:a1和a2为拟合系数。
图 4(a)给出了含填充孔的圆形蜂窝结构冲击端平台应力与速度参量v2之间的关系, 图 4(b)所示为局部放大图。从局部放大图中也能够看出, 当冲击速度小于临界速度v1(c), 即蜂窝结构的变形仍处于准静态模式时, 冲击端平台应力曲线表现出明显的非线性。S.R.Reid等[11]基于一维冲击波理论建立了均匀多孔材料动态冲击平台应力σp与冲击速度之间的关系, 并指出当冲击速度大于一定值时, 规则蜂窝结构的平台应力存在速度效应, 平台应力值与速度的平方成线性关系。各蜂窝结构的拟合参数值如表 2所示, 本文中的拟合结果表明, 当蜂窝结构中含有随机固体填充孔时这样的速度效应同样存在, 但此时拟合参数与填充比的大小有关。从表中可知, 随着填充比的增大, 拟合参数a2也有所增大, 这说明填充比的增加使得圆形蜂窝结构冲击端平台应力的速度效应更加明显。
表 2 平台应力拟合关系式中的系数Table 2. Coefficients in fitting equation of plateau stressα a1/MPa a2/(Pa·s2·m-2) 0.00 0.086 0.238 0.02 0.071 0.244 0.04 0.055 0.252 0.06 0.043 0.260 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时候, 相对平台应力甚至有下降的趋势。
2.4 随机填充孔对能量吸收特性的影响
作为一种高效的阻尼材料, 多孔胞元材料已被广泛地应用于各种能量吸收结构中。评价多胞材料与结构能量吸收能力大小的重要指标就是单位质量吸收的能量大小[12]:
Wm=Wvρρs (8) 式中: Wv=∫εdεyσdε为蜂窝材料单位体积吸收的能量(用名义应力应变曲线下所围成的面积来表征)。图 7给出了不同填充比的蜂窝材料单位质量吸收的能量与相对压缩量之间的关系。
由图 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) 当含随机固体填充孔的圆形蜂窝结构变形处于准静态模式时, 冲击端平台应力曲线表现出明显的非线性, 冲击速度大于由准静态变形模式向过渡模式转变的临界速度v1(c)时, 平台应力σp与v2呈线性关系, 存在明显的速度效应。
(3) 低速冲击下, 含固体填充孔的蜂窝结构平台应力随孔洞填充比的增大而显著降低, 导致含填充孔蜂窝结构的吸能能力下降, 冲击速度为3 m/s时, α=0.06的蜂窝结构平台应力约为相同相对密度的规则蜂窝结构平台应力的50%, 单位质量吸收能量大小仅为规则圆形蜂窝结构的37.5%。随着冲击速度的提高, 固体填充孔的牵制效应使得动态冲击的惯性效应影响增强, 冲击端应力应变曲线出现明显的应力尖峰, 提高了蜂窝结构的单位质量吸能能力。
本文中对含随机固体填充孔的圆形蜂窝结构的面内冲击性能的研究可为有缺陷蜂窝材料的研究提供参考, 但关于随机固体填充孔与蜂窝结构动态响应特性更深层次的探讨还有待进一步展开。
-
表 1 含随机固体填充孔蜂窝结构参数
Table 1. Parameters of honeycombs with randomly arranged solid inclusions
α t/mm ρe/(g·cm-3) 0.02 0.125 33.106 0.04 0.096 42.650 0.06 0.066 61.490 0.08 0.036 114.130 表 2 平台应力拟合关系式中的系数
Table 2. Coefficients in fitting equation of plateau stress
α a1/MPa a2/(Pa·s2·m-2) 0.00 0.086 0.238 0.02 0.071 0.244 0.04 0.055 0.252 0.06 0.043 0.260 -
[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. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=lxjz200604004Lu 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. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=lxjz200604004 [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. http://www.sciencedirect.com/science/article/pii/S1359645405002089 [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. http://d.wanfangdata.com.cn/Periodical/lxxb-e200704005 [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. doi: 10.1080/01418619608242994 [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. http://www.sciencedirect.com/science/article/pii/S0020740399001228 [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. http://www.sciencedirect.com/science/article/pii/S0734743X08000675 [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. http://www.sciencedirect.com/science/article/pii/S0734743X09000281 [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. http://www.sciencedirect.com/science/article/pii/S0263822312004837 [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. http://www.ingentaconnect.com/content/el/00225096/2005/00000053/00000010/art00002 [11] Reid S R, Peng C. Dynamic uniaxial crushing of wood[J]. International Journal of Impact Engineering, 1997, 19(5/6): 531-570. http://www.sciencedirect.com/science/article/pii/S0734743X9700016X [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. http://www.sciencedirect.com/science/article/pii/S1359645404003131 期刊类型引用(5)
1. 赵著杰,侯海量,李典,王克,姚梦雷. 部分充液多胞元结构的面内动态力学特性研究. 爆炸与冲击. 2022(03): 37-52 . 本站查看
2. 李响,曹祥斌,杨蔚华,李林. 泡沫填充类蜂窝夹层结构的耐撞性. 武汉科技大学学报. 2021(01): 20-26 . 百度学术
3. 孙玉瑾,孙德强,安兴,焦思涵. 圆形蜂窝的异面冲击性能研究. 陕西科技大学学报. 2021(05): 139-145 . 百度学术
4. 游颖,张泽涛,郭琪. 两种新型蜂窝结构的承载性能对比. 湖北工业大学学报. 2020(01): 10-13 . 百度学术
5. 张金山,乔及森,孔海勇,苗红丽. 铝圆管蜂窝材料制备及其准静态压缩性能研究. 机械工程学报. 2020(16): 78-83 . 百度学术
其他类型引用(5)
-