泡沫杆撞击刚性壁的动态压溃模型 dynamic crushing models for a foam rod striking a rigid wall.pdfVIP

第33卷第6期 爆 炸 与 冲 击 V01．33·No．6 WAVES Nov．·2013 2013年11月 EXPI。OSIONANDSHOCK 文章编号：i001—1455(2013)06—0587—07 泡沫杆撞击刚性壁的动态压溃模型” 王长峰，郑志军，虞吉林 (中国科学技术大学中国科学院材料力学行为与设计重点实验室，安徽合肥230026) 摘要：针对泡沫杆撞击刚性壁的情形建立了2类动态压溃模型：一维冲击波模型和三维细观有限元模 型。以连续介质框架下的应力波理论为基础，并假定了刚性一非线性塑性硬化的加载和刚性卸载的本构关系， 建立了一维冲击波模型．给出了冲击波波后应变与冲击时间的隐式表达式。利用随机Voronoi技术构建了闭 孔泡沫金属结构的三维细观有限元模型，使用ABAQus／Explicit有限元软件模拟了泡沫材料的动态压溃过 程，并基于最小二乘法计算局部变形梯度和局部应变得到了三维泡沫结构的应变场。通过理论解和数值解的 比较，发现该理论模型能够较好地预测泡沫金属杆撞击刚性壁的力学行为，得到了较为精确的结果。 关键词：固体力学；冲击波模型；Voronoi结构；泡沫金属；局部应变场 中图分类号：0347．3 国标学科代码：13015 文献标志码：A 泡沫金属具有稳定、持续的能量吸收平台，因而广泛用作能量吸收材料。在动态冲击下，泡沫金属 表现出变形局部化和应力增强等特点，但由于实验技术的限制，关于其变形机理和本构的应变率敏感性 问题仍然存在很大分歧_。已经发展了一些动态压溃模型，如冲击波模型、弹簧一质量模型和有限元模 型，但多数研究局限于对泡沫金属宏观响应的预测，缺乏对细观变形机制的分析和比较。 S．Reid等12。认为泡沫金属在动态冲击下的局部化变形行为是结构冲击波传播导致的，并率先提出 了率无关、刚性理想塑性一锁定(R—P—P—I，)一维冲击波模型，该模型可以很好地估计泡沫的冲击增强∞]。 Z．J．Zheng等“进一步发展了刚性一线性硬化塑性一锁定(R—LHP—I。)模型，用于预测泡沫金属冲击模式 和过渡模式下的动态行为。S．I。．Lopatnikov等n1考虑了弹性行为的影响。实际上，泡沫材料的动态压 实应变随着冲击速度的提高而增加…。不同学者在锁定应变的选取上并未达成一致。然而模型的预测结 果对压实应变的选取非常敏感一?。为此，发展了一些具有非线性塑性硬化段的理论模型，如弹性一理想 塑性一应变硬化(E—P—P—H)模型8。S．Patto{atto等口。将R—P—PI。模型的锁定硬化段改为幂次硬化段， 但只预测了恒速加载下的冲击应变和冲击波波速。 V．P．W．Shim等9首先提出了一维弹簧一质量模型，Q．M．Li等。1 01进一步研究了泡沫金属中应力 衰减和强化的问题。王志华等“。12发展了二维弹簧一质量一连杆模型，研究了泡沫材料中应力增强现象 以及非均匀性对冲击波传播特性的影响。值得注意的是，J．J．Harrigan等心】指出弹簧一质量模型并不适 合用来预测泡沫材料内压缩产生的冲击波的不连续性。细观有限元模型能够克服这一缺点，但已有的 研究较多地局限在二维的情况，如Z．J．Zheng等I川研究了随机蜂窝的动态冲击行为。三维结构能够较 为真实地模拟泡沫金属的力学性能，宋延泽等。1“嘲研究了规则正十四面体泡沫的能量吸收特性以及三 维闭孔Voronoi模型的随机性对平台应力的影响。 本文中着重研究泡沫杆撞击刚性壁时冲击波波后的应变和应力随时间的变化，发展2类动态压溃 模型：一维冲击波模型和三维细观有限元模型。一维冲击波模型给出对于任意塑性硬化情况下应变与 时间的关系式，三维细观有限元模型给出局部应变场。通过冲击波模型与细观有限元模型的比较，验证 冲击波模型的有效性。 *收稿日期：2012—06—29；修回日期：2012一11—19 基金项目：国家自然科学基金项目90916