闭孔泡沫金属变形模式的有限元分析 Finite element analysis on deformation modes of closed-cell metallic foam.pdfVIP

第34卷第4期 爆炸与冲击 V01．34，No．4 2014年7月 EXPLCISIONANDSHOCKWAVES Jul．，2014 1455(2014)04—0464—07 文章编号：1001 闭孔泡沫金属变形模式的有限元分析。 李妍妍，郑志军，虞吉林，王长峰 (中国科学技术大学中国科学院材料力学行为和设计重点实验室，安徽合肥230026) 形行为。随着冲击速度的提高，得到了3种变形模式：准静态均匀模式、过渡模式和冲击模式，并以相对密度 和冲击速度为坐标建立了变形模式图。引入应力均匀性指标和变形局部化指标，确定了模式转化的临界速 度，并与已有的冲击速度预测公式进行了比较。根据临界速度的数值和理论结果，提出了一种确定锁定应变 的方案，结果介于压实应变和完全密实应变之间。 关键词：固体力学；变形模式；有限元方法；Voronoi模型；临界速度；锁定应变 中图分类号：0347．1 国标学科代码：1301515文献标志码：A 泡沫金属已广泛用作能量吸收材料，很有必要了解其动态力学行为。通过细观有限元模拟能有效 地研究随机性或缺陷对泡沫金属力学性能的影响。 已有的有限元模拟多建立在二维情况下。D．Ruan等Ⅲ研究了不同冲击速度下规则蜂窝的变形模 研究随机蜂窝的动态冲击中从变形机理的角度提出了变形模式的分类，即准静态均匀模式、过渡模式和 冲击模式。G．W．Ma等H1则采用随机模式、过渡模式和渐进模式来分类。在蜂窝面内冲击性能的进一 步研究中，刘颖等口1和胡玲玲等卟3分别考虑了缺陷分布不均匀性、胞壁夹角的影响。 态压缩下变形出现的方式。宋延泽等研究了规则正十四面体泡沫的能量吸收性能[8]，并进一步研究了 随机性对闭孔Voronoi模型的平台应力的影响口]。王鹏飞等[103在改进的Hopkison压杆上，开展闭孔 泡沫金属的动态冲击中，发现变形模式对实验影响很大，高速冲击下得到的结果不能真实反映材料的动 态力学性能。因此，有必要对泡沫金属的变形模式和模式转变的临界速度进行深入研究。 本文中，采用闭孑L泡沫金属的细观有限元模型，研究相对密度和随机性对变形模式的影响，给出模 式转变的临界冲击速度的计算方法，并与文献中存在的理论预测公式进行比较和讨论。 1数值计算模型 1．1三维随机Voronoi模型 基于三维随机Voronoi技术构建闭孔泡沫金属的计算模型。在指定的区域内随机生成N个点作 为成核点，各成核点之间的最小距离d…限制为： d…一(1一k)d。一(1一是)压(Vcell／4)∽ (1) 式中：d。和V叫。分别为规则正十四面体模型中成核点间的最小距离和单胞的体积，k为不规则因子。将 mm×12mm×12 留指定区域内的构型作为计算模型。本文中，指定区域大小为16 mm，成核点305 个，不规则因子取0．3，一个计算样本如图1所示。不考虑胞内气压的影响，并采用均一的胞壁厚度，由 给定的相对密度确定胞壁厚度。计算中，选用了4种相对密度p一0．05，0．10，0．15，0．20。 *收稿日期：2012—12—14；修回日期：2013-03—21 基金项目：国家自然科学基金项目90916026) 作者简介：李妍妍(1989)，女，硕士研究生。 万方数据 第4期 李妍妍等：闭孔泡沫金属变形模式的有限元分析 1．2有限元模拟 个端面分别与2个刚性板接触。一个刚性板固定，另一个刚性板以恒定的速度冲击泡沫金属模型。采 mm，并调 用S3R和S4R壳单元对泡沫金属模型进行网格划分，经过收敛性分析，网格尺寸为约0．03 整较小的单元以确保计算效率，共划分成近10000个S3R壳单元和38000个S4R壳单元。基体材料