四边手性蜂窝动态压溃行为的数值模拟 numerical simulation on dynamic crushing behaviors of tetrachiral honeycombs.pdfVIP

第34卷第2期 爆炸与冲击 V01．34。No．2 WAVES Mar．，2014 2014年3月 EXPLOSl0NANDSHOCK 文章编号：1001—1455(2014)02—0181—07 四边手性蜂窝动态压溃行为的数值模拟。 卢子兴，李康 (北京航空航天大学固体力学研究所，北京100191) 摘要：建立了四边手性蜂窝的有限元模型，采用数值模拟方法研究了四边手性蜂窝在不同冲击速度下的 变形模式和能量吸收等动力学响应特性，并同普通六边形蜂窝的冲击行为进行了对比。计算得到了这2种蜂 窝的变形模式图、动力响应曲线和能量吸收曲线。模拟结果表明：低速冲击下，四边形手性蜂窝的变形模式为 “Z”字形；高速冲击下，四边手性蜂窝的变形模式与普通蜂窝的“I”字形模式类似；在适中速度的冲击下，四边 手性蜂窝表现出兼具高速冲击和低速冲击特征的一种过渡态变形模式；随着冲击速度的提高，局部变形带由 固定端向冲击端移动，并且能量吸收能力也随之提高；在中、低速度的冲击下，能够观察到拉胀材料压缩时特 有的“缩颈”现象。 关键词：固体力学；动态压溃；ANSYS／LS-DYNA；手性蜂窝；变形机制；能量吸收；拉胀 中图分类号：0347 国标学科代码：13015 文献标志码：A 与普通多孔材料相比，拉胀多孔材料具有更高的冲击阻抗、抗剪能力、抗凹能力和能量吸收能力，在 航空航天领域中有着广泛的应用前景，如用作结构防护材料和夹层板芯层材料等E15]。在现有材料的基 础上，寻找能量吸收效率更高的材料，并对现有的吸能材料和结构进行优化设计和改进，已经成为当今 一系列蜂窝在面内准静态压缩下的力学行为进行了实验研究，给出了压缩应变为1％或2％时的应力应 变关系曲线，得到了面内弹性模量和泊松比，并进行了有限元模拟与实验的对比。最近，A。Spadoni 等uu采用微极连续介质模型来描述拉胀手性蜂窝，试图消除由泊松比为一1引起的本构关系的不确定 性。F．Dos Reis等[12借助渐进均匀化方法对5种二维周期性拉胀蜂窝进行了分析，得到了模量和泊松 比的封闭解，认为导致负泊松比的2种主要机制是内凹机制和旋转机制。目前，仍有一些新型拉胀拓扑 结构在不断涌现[5J3。1 4|，寻找更高负泊松比的蜂窝拓扑结构已成为该领域的主题之一。拉胀蜂窝是一 种新型多孔材料。要在结构上广泛使用拉胀蜂窝，首先需要对其力学性能有充分的研究和认识。当前 的研究工作主要集中于发现新的拉胀蜂窝拓扑结构和预测其准静态力学性能及解释负泊松比变形机制 上，然而，实际应用中对其用作夹层结构或能量吸收构件材料时，遭受外来物冲击的动态行为尚缺乏了 解，因此，对于拉胀蜂窝结构在冲击载荷作用下的动力学性能展开研究就显得尤为必要，而目前对此还 未见报到。 本文中，采用数值模拟方法，对四边手性蜂窝在不同冲击速度下的变形模式和能量吸收等动力学响 应特性展开研究，以确定这种蜂窝在不同冲击速度下的变形模式和能量吸收性能，并同普通六边形蜂窝 的冲击行为进行对比，以期为进一步研究三维拉胀泡沫材料的动态冲击行为奠定基础。 1模型的建立 1．1手性蜂窝的几何结构 *收稿日期：2012—09—13；修回日期：2013—03—29 基金项目：国家自然科学基金项目11272030) 作者简介：卢子兴(1960一)，男，博士，教授，博士生导师。 万方数据 爆 炸 与 冲 击 第34卷 的几何结构参数，其中^r、占分别为支柱的长度、圆 形节点的半径及壁厚，支柱和圆形节点的壁厚相同。 所有算例中取z一20mm，艿一1．5mlTl，r一5mm。 1．2有限元模型 采用ANSYS／LS-DYNA显式动力分析有限元 软件对手性蜂窝面内冲击特性进行数值分析