岩石损伤时效模型资料.pptVIP

岩石动态破坏的时效损伤本构模型 单 仁 亮 （中国矿大 北京，土木工程系） 1 概述 1.1 General：弹模和强度（峰前） With d? /dt ? ，E?，?f?，但在103左右突变 强度（Before 1970’s） 对数方程类[?f ? log(d? /dt) ]和指数方程类[?f ? (d? /dt )n )] 1.2 本构方程（After 1970’s） (1) 过应力模型 本构方程为： 优点 能够反映强度随应变率而提高 缺陷 不能反映弹性模量随加载率而变化的特征。 (2) 粘弹性连续损伤本构模型（郑永来等） 首先将岩石看成某种粘弹性体，再假设整个粘弹性体在变形过程中收到了连续损伤，也就是将粘弹性模型推得的应力σ (t)用σ (t)/(1- D)取代，σ (t)由粘弹性模型确定，D为损伤参量，D = 1 - [(ε /ε 0) m + 1]exp[-(ε /ε 0) m]，ε为应变，ε 0和m是与材料性质和形状相关的Weibull分布参数。他们分别对Ⅱ型标准线性粘弹性模型和广义流变模型考虑了损伤效应，得到了不同的粘弹性连续损伤模型。 粘弹性连续损伤模型比较全面反映了岩石变形的应变率效应。但模型的参数很多，参数确定颇为费事，且需要大量的不同应变率的实验。此外，该模型在物理概念上还有一些模糊，因为粘弹性体是由弹性元件和粘性元件组成的，弹性固体元件可以受到损伤，而粘性液体元件不会受到损伤，让整个粘弹性体的应力σ (t)用σ (t)/(1- D)取代，意味着在变形中粘性元件也受到了损伤，即粘缸在变形中断面也缩小了。所以该模型尚需商榷。 近年来我们在证明了利用SHPB装置测试岩石动态全应力—应变曲线的可行性[12]的基础上，测试了大量的花岗岩和大理岩的动态全应力—应变曲线，对它们在受到冲击后从变形至破坏全过程的机制进行了深入的讨论[12-15]。本文进一步提出了一个简明的岩石冲击损伤时效模型，并得到了SHPB实验验证，模拟效果良好。 2 实测的本构曲线 3 时效损伤模型 3.1 基本假设 1) 岩石单元同时具有统计损伤特性和粘性液体的特性，因而可以把岩石试件看成损伤体Da和粘缸? b的并联体，如下图所示。 2) 损伤体Da在损伤之前是线弹性的，平均弹性模量为E，强度服从参数为(m, ?)的Weibull分布。其概率密度? (? )、损伤参数D以及本构关系? - ? 可以分别由公式(7)，(8)，(9)求得。 3) 粘缸没有损伤特性，遵循的本构关系为 并联组合体中损伤体的应变等于粘缸的应变，因而组合体的应变等于两个分体的应变，而应力等于两个分体的应力之和，即： 将损伤体和粘缸的本构关系代入(10)式，容易得到并联体的本构关系： 对上述模型的本构关系进行数值计算需要确定四个参数E、m、?、? ，它们的确定往往需要分析实测的应变波形，并且不可避免的需要一定的试算。通常E与岩石的冲击应力应变曲线的初始上升斜率相近，说明E的确能够表示未损岩石的初始弹性模量，m表示weibull分布中分布曲线的形状系数，一般它在1附近变化，? 一般位于峰值应力对应的应变与平均应变之间，?的变化范围一般在0.1~0.5之间；而应变? 和应变速率d? / dt应当用实测数组，下面两个算例中的应变和应变率数组都是在SHPB实验得到的。 4 结论 根据对实测岩石动态应力应变曲线的分析和现有固体本构模型，本文提出了一个简明的岩石冲击破坏时效损伤模型，模型的物理概念清晰。模型的数值计算结果与SHPB实验结果对比表明：该模型在总体上很好的模拟了大理岩的冲击破坏特性，较好的模拟了花岗岩的冲击破坏特性。 * * σ S E P P ε Fig.1 Excess stress model Fig.2 Analytical chart for equation of excess stress model Fig.5 Dynamic stress-strain curves of granite Fig.6 Dynamic stress-strain curves of marble Da ? b Fig.7 Time dependent damage model ? ≥0 (7) ? ≥0 (8) ? ≥0 (9) 10 11 3.2 实验验证 (1) 21th granite (2) 89th marble 花岗岩的应力应变曲线在最大值之前大多具有跃进特性，但在第一极大值之前具有较好的线弹性，且弹性模量E与应变率没有明显的相关性，在最大值之后形态具有较大的差异： *