冲击损伤应力波传播过程的数值流形方法模拟研究.pdfVIP

第九届全国冲击动力学学术会议论文集 冲击损伤应力波传播过程的数值流形方法模拟 刘红岩 王贵和 赵乐华 周辉峰马孝春 (中国地质人学(北京)工程技术学院，北京，100083) 摘要：根据岩石在冲击载荷作用下的损伤与破坏特征，对现有的仅适用于模拟基于线性本构关系和恒定载荷作用下 物体变形规律的数值流形程序进行了扩展。程序的扩展内容主要集中在两个方面：把载荷形式由原来的恒定载荷扩 展为任意的三角载荷，改变原来的单一线弹性本构关系为由冲击损伤本构关系和线弹性本构关系共同组成的复合本 构关系。最后利用扩展后的程序对一维应力波的传播过程进行了模拟，结果很好地再现了应力波的传播过程。 关键词：一维应力波，冲击损伤，数值流形方法，损伤演化 中图分类号： TGl56 文献标志码： 0 引言 岩石材料在国防与民用建筑中都得到了广泛的应用，它们在使用过程都不可避免地要承受来自 外界冲击载荷的作用，由于冲击载荷的作用将在岩囱．内部产生应力波，由于岩石都是天然的损伤材 料，在应力波的作用下，其内部的损伤也将会加剧，从而加速材料的破坏。因此，对冲击载荷作用 下岩石内部应力波的传播进行研究具有十分重要的理论意义和实用价值。 由于数值模拟具有可重复性和低成本，现已在科学研究领域中发挥着越来越重要的作用。数值 提出的一种新的数值分析方法【l儿21。由于其具有明显的优越性，目前已在岩土工程分析中得到了较为 广泛的应用【3J。然而由于其提出的时间比较短，其应用范围还不象有限元和离散元那么广泛。但是 由丁其在计算的过程中是采用的两套网格——物理网格和数学网格，因此在计算材料变形过程中裂 纹的产生、扩展等方面具有明显的优越性。目前断裂和损伤理论也都被成功地应用于岩石冲击动力 学的研究中，并建立了相应的断裂和损伤模犁。因此如何更好地利用数值流形方法来对此类问题进 行研究也是一个重要的课题，本文就在石根华博士开发的线弹性数值流形方法程序的基础上加入以 应力波衰减为基础的冲击损伤模型，来模拟一维岩石长杆在冲击载荷作用下的损伤演化及其内部的 应力波传播过程。 1 岩石冲击损伤应力波衰减模型 1．1损伤变量的定义 量表达式14J： D=l一(v／％)2 (1) 式中，’，为受载损伤后岩石的纵波波速，1，o为受载前岩石的纵波波速。由于(∥％)2表示受载损伤前 后岩石弹性模量之比，所以D表述的是岩石的相关模量的变化。 2．2损伤演化方程及本构方程 岩石冲击实验和超声波实验研究揭示了应力波在损伤岩石中传播时的衰减特性，即随着损伤的 不断增加，应力波在其内传播的衰减愈人，表现为衰减系数的不断增加，损伤耗散能增加，因此可 以认为衰减系数口和损伤耗散能K之间的关系式为【5】： 口=口o+K。％ (2) 反映了损伤演化的基本规律。其中K。是拟合常数，％可以认为是岩石介质的初始衰减系数， ·基金项目：中国博十后基金资助(20060400498) 328 第九届全国冲击动力学学术会议论文集 随着应力波的传播和损伤的演化，介质的衰减系数增加，其损伤程度越来越大。损伤参数与衰减系 数间的关系为线性关系，即： 伐=A七BD (3) 将式(2)和式(3)写成导数形式为： 西=K。圪 (4) 和 D=C磊 (5) 式中C是材料肯量(C=1／B)。式(4)和式(5)为体积拉伸下的损伤演化方程。由于岩石在压 缩状态下，损伤演化一般很小，所以一般不考虑岩石在压缩状态下的损伤累计。 损伤材料的本构关系为： 盯，，=3K(1一D)哦+2G(1一D)P盯