多孔材料内冲击导致特征结构及其演化.pdfVIP

中国核科学技术进展报告（第一卷） 计算物理分卷 Progress Report on China Nuclear Science Technology （Vol.1 ） 2009 年11 月 多孔材料内冲击导致的特征结构及其演化 许爱国，张广财，蔚喜军，朱建士 （北京应用物理与计算数学研究所国家级计算物理重点实验室，北京 100088 ） 摘要：使用颗粒接触物质点方法对冲击波与多孔材料的相互作用过程进行模拟， 对其特征结构进行形态学分析，引 入“结构相似”和“过程相似” 的概念，对不同温度阈值条件下的特征结构极其演化过程进行比较分析。发现：在同一材 料中，如果冲击强度和温度阈值选取合适，那么高温区域及其演化过程会表现出极其相似的特征。当孔隙度较高时， 过程相似要求的温度阈值和冲击速度成二次抛物线关系；随着孔隙度降低，温度阈值随冲击速度升高的速率有所加 快。在不同材料中，在冲击强度一定的情况下，过程相似要求的温度阈值和孔隙度在一定范围内表现出幂律行为。 其它类型的特征结构可以按同样方式进行研究分析。 关键词：多孔材料；冲击波；特征结构；物质点法；形态学分析 非均质材料在国防与民用等工业领域有着极其广泛的应用[1] [2-3] 。除了颗粒材料 ，多孔材料很早 就被用来保护精密设备和仪器免受猛烈震动和冲击。材料的非均匀性能够极大地改变其冲击动力学 特征。当物体受到弱冲击作用时，其中空隙的存在可能使两碰撞物体出现二次碰撞；当冲击作用较 强时，材料内部空隙处可能会出现射流[4，5] 。从整体上看，在冲击波作用下，多孔材料内会出现各类 特征结构。例如“高压区”、“高温区”、“粒子速度较大区” ，“偏应力较大区”等。这些特征结构是实际 应用过程中容易出“ 问题” （相变、化学反应、粒子微喷、塑性变形等）的区域，其各类特征和规律 值得关注。虽然多孔材料冲击动力学的研究已有较长的历史，但以前的研究主要集中在冲击 Hugoniots[6] [7] 和状态方程 。为了研究多孔材料中复杂的时空耗散结构，许爱国等将积分几何和数值 图像处理中的Minkowski 泛函引入冲击波物理，从形态学的角度来分析温度、密度、压强、粒子速 度等特征量在材料中的分布及其演化[8] 。 多孔材料的冲击响应过程涉及到大变形和很强的非线性，其理论描述非常复杂。在相关实验中， 一些快速过程和细节的测量在技术上也具有极强的挑战性。所以，数值模拟对于提供完整动力学图 像、对于推动本学科的发展都起着不可替代的作用。从模拟方面来说，分子动力学可以揭示出空隙 塌缩的一些微观机制[9,10] ，但其所涉及的时间和空间尺度尚局限在纳秒和几十纳米量级，远不能满足 实验要求。纯拉氏和纯欧氏算法在处理这类涉及到材料大变形的问题时均遇到了一些不易克服的困 难：拉氏方法能够很好地跟踪物质的界面，但会遇到网格畸变的问题；在欧氏方法中，网格不随物 体的变形而变化，但是它不易跟踪物质的界面。为充分发挥欧氏和拉氏算法的优点、避免其缺点， 本文将使用最近发展起来的一种混合方法（物质点法）[11,12]来模拟多孔材料的冲击响应过程。 1 材料模型和数值方法 1.1 材料模型 在我们的模拟研究中，多孔材料的初态是由一块完整的固体内嵌一定数量的空隙来构成。固体材料 特性用带线性硬化的von Mises 塑性模型来描述[4] 。材料在屈服条件满足之前表现出线弹性响应； 作者简介：许爱国（1970—），男，汉族，山东人，研究员，博士，理论物理 基金项目：国家自然科学基金资助项目10775018,、中物院发展基金和实验室基金资助项目 21 2 屈服判据为3J 2 −σ 0 ，其中σ 是塑性屈服应力， Y Y