基于西原模型的非线性损伤蠕变模型_付腾飞.pdf

基基础础研研究究 基于西原模型的非线性损伤蠕变模型 付腾飞， 付宝杰 （安徽理工大学 能源与安全学院， 安徽 淮南 232100） ［摘 要］ 针对传统的西原模型不能较好地反映岩石蠕变非稳定阶段， 且多数学者对西原模型 的修正是通过屈服强度阈值来界定非稳定阶段， 而通过时间阈值考虑非稳定蠕变阶段界定的研究相对 较少。 在西原模型的基础上， 引入 Kachanov 损伤理论及时间损伤阈值， 建立了一种同时描述3 个阶段 的非线性损伤蠕变方程式。 通过灰岩蠕变试验得到试验数据， 对新的非线性损伤蠕变模型的参数进行 反演， 并进行曲线拟合， 分析岩石的蠕变损伤过程和对应阶段的稳定性， 证明了该模型的正确性和合 理性， 为岩石流变失稳现象提供理论依据。 ［关键词］ 西原模型； 损伤； 非线性； 加速蠕变 ［中图分类号］ TD325 ［文献标识码］ A ［文章编号］ 1006-6225 （2016） 04-0001-04 Nonlinear Damage Creep Model Based on Xiyuan Model FU Teng-fei， FU Bao-jie （Energy and Safety School， Anhui University of Science Technology， Huainan 232001， China） Abstract To the problems of instability step of rock creep deation could not be reflected by traditional Xiyuan model， and as the modified Xiyuan model for most of scholars， instability deation step was confirmed by yield strength threshold value， but instabili- ty creep deation step was confirmed by time threshold value was less. On the basis of Xiyuan model， Kachanov damage theory and time damage threshold was introduced， nonlinear damage creep ula that explained three steps at the same time was built. The ex- perimental data was obtained by creep experiment of limestone， the parameters of new nonlinear damage creep model were inverted， and curve fitting was analyzed， the creep damage process of rock and its stability of corresponding step， the results proved the model was right and reasonable， it references for rock rheology instability phenomenon theory. Key words xiyuan model； damage； nonlinear； acceleration creep ［收稿日期］ 2015-11-30［DOI］ 10. 13532/ j. cnki. cn11-3677/ td. 2016. 04. 001 ［基金项目］ 国家自然科学基金 （51174002） ［作者简介］ 付腾飞 （1992-）， 男， 山西长治人， 在读硕士研究生， 研究方向为矿山压力与岩层控制。 ［引用格式］ 付腾飞， 付宝杰 . 基于西原模型的非线性损伤蠕变模型 ［J］ . 煤矿开采， 2016， 21 （4） 1-4. 岩石流变在岩石工程 （如隧道工程、 煤矿巷 道支护） 中普遍存在， 其中许多破坏都因岩石流 变破坏而引起， 如 边坡的突然滑动、 巷道的岩爆 等， 对岩石工程的长期稳定性和安全性有显著影 响［1］。 对于岩石流变的研究， 现有的经典流变模 型都只能描述岩石的线性衰减蠕变和稳定蠕变， 不 能同时描述岩石的加速蠕变和表现岩石的破坏特 性， 因此更多的研究人员建立能反映加速蠕变阶段 的非线性流变方程， 以进一步研究岩石的流变破 坏。 可以通过增加新的非线性元件或用非线性元件 代替现有流变模型中的基本元件来建立非线性流变 方程， 如 胡波［2］等将可阶跃函数表示的开关元 件和损伤元件引入 CVISC 蠕变模型研究了岩石蠕 变特性； 陈沅江［3］等考虑裂隙塑性体和蠕变体， 并将它们和虎克体及开尔文体相结合得到了一种能 描述加速蠕变阶段的新的复合流变力学模型； 袁 林［4］等引入非线性 TN 体与传统西原模型组合为 非线性黏弹塑性流变模型， 更有效地描述了岩石的 流变特性。 也可以通过考虑经典流变模型中元件模 型参数为非常定值串并联组合成来描述岩石蠕变全 过程的新模型， 如 丁志坤［5］在考虑参数为非定 常的情况下， 采用 H-K 体的三元件模型更为准确 地反映了岩石不同应力水平下的黏弹性变形性能； 曹树刚［6］在分析岩石的全应力-应变曲线和蠕变 曲线基础上， 修正西原模型中的黏滞系数， 建立了 参数非线性蠕变型。 而在工程应用中， 岩石材料的强度与该材料的 流变特性和内部损伤机理是密切相关的［1］。 不仅 岩石材料本身非均匀性会导致加速蠕变， 在外载荷 作用下岩体内微观裂纹和空洞的发展而形成的损伤 也同样会导致加速蠕变。 事实上所谓岩石损伤就是 外荷载作用下岩石内部微观结构的不断变形和调 整， 当结构变化到一定程度就会表现为岩石宏观上 的失稳破坏。 因此， 很多学者将损伤引入蠕变研究 1 第 21 卷 第 4 期 （总第 131 期） 2016 年 8 月 煤 矿 开 采 COAL MINING TECHNOLOGY Vol. 21No. 4 （Series No. 131） August 2016 ChaoXing 中。 如 张强勇等［7］考虑岩体的流变损伤劣化效 应， 建立了变参数的蠕变损伤本构模型， 反映岩体 逐渐损伤劣化引起的蠕变参数改变对流变特性的影 响； 刘桃根［8］等应用损伤力学， 建立了 3 种蠕变 损伤本构模型， 都很好地描述了砂岩的流变特性。 岩石的非线性蠕变是其内部结构调整在蠕变过 程中的综合表现。 本文以西原模型为基础， 在蠕变 本构关系的推导中考虑蠕变损伤和时间阈值的影 响， 构造岩石非线性蠕变损伤机制解释岩石蠕变过 程中所显现的非线性变形特征。 结合试验数据， 利 用 Matlab 软件反演出新蠕变方程中的参数， 将拟 合曲线与实际试验曲线进行对比分析， 结果具有较 高的准确性， 为岩体流变的非线性研究提供一定的 理论参考。 1 流变模型的建立 1. 1 西原模型 根据应力水平的不同， 岩石蠕变曲线有 3 种情 况 较低应力水平下初始蠕变后进入速率为零的衰 减蠕变曲线； 低应力水平下初始蠕变后进入速率为 常数的定常蠕变曲线； 高应力水平下岩石进入具有 初始蠕变， 稳态蠕变和加速蠕变的蠕变全过程曲 线［9］。 传统的西原模型所含元件少、 能全面反映 岩石弹-黏弹-黏塑变形特性， 在岩石材料研究中 得到广泛应用， 其模型如图 1 所示。 但其一般只能 描述第一、 第二阶段的蠕变， 并不能反映加速蠕变 阶段。 图 1 西原模型 传统西原模型的本构、 蠕变方程为 E2ε E1E2 η1 ε σ E2 η1 E2 η2 E1 η1 ■ ■ ■ ■ σ E1E2 η1η2（σ - σs） （1） εc σ0 E1 σ0 E2 1 - e -E 2 η1 t （） σ0 - σ s η1 t（2） 式中， σ0为初始应力； σs为临界应力； εc为传统 西原模型中的应变； E 为弹性模量； η 为黏性系 数； t 为时间。 1. 2 非线性损伤蠕变模型 岩石的蠕变过程与其内在损伤机理有密切关 系， 其非线性蠕变是其内部结构变化在蠕变过程中 的综合表现， 本文引入损伤变量， 同时考虑加载时 间因素对流变特性影响， 建立非线性蠕变模型， 修 正西原体蠕变方程。 岩石材料的损伤是指在外载环境作用下， 由于 微观结构的大量微裂纹、 微孔洞的产生、 扩展和连 接， 导致该材料的力学性能的劣化过程。 其中已经 破坏材料微元与微元总体的比值即为损伤变量。 文 中设岩石材料微元总体面积为 1， 微元受损面积为 D， 则未受损部分的面积为 1-D。 假设无损岩石的 应力和应变分别为 σw和 εw， 其中 εw与式 （2） 中 εc的意义是一致的， 已破坏岩石材料的应力和应 变分别为 σf和 εf， 则此岩石的总应变［10］为 ε εw1 - D（） εfD（3） 在忽略弹性变形， 考虑损伤情况下， 蠕变率假 设为指数率耦合 ε f a σ 1 - D ■ ■ ■ ■ m （4） 取初始条件 t0 时， εf0 得 εf a σ 1 - D ■ ■ ■ ■ m t（5） 依据 Kachanov 提出的损伤变量， 其演变方程 为 D b σ 1 - D ■ ■ ■ ■ n （6） 式中， a， b， m， n 是材料常数， 均由试验确定。 取初始条件 t tb， D 0， 其中 tb为蠕变损伤 临界破坏时间 （损伤时间阈值）， 由 （6） 得 1 - 1 - D（） n1 1 n bσnt - tb（ ）（7） 对式 （7） 整理得 D 0 t ≤ tb 1 - 1 - b n 1（）σnt - tb（）［］ n1 t ＞ t b ｛ （8） 将式 （8） 代入式 （5） 得已破坏岩石材料的 应变表达式为 εf aσmt 1 - bσnt - tb（）n 1（）［］ m n1 （9） 所以， 将式 （2） 和 （9） 代入式 （3） 得到考 虑损伤的蠕变方程为 ε σ0 E1 σ0 E2 1 - e -E 2 η1 t （） σ0 - σ s η2 t ■ ■ ■ ■ 1 - D（） aσmtD 1 - bσnt - tb（）n 1（）［］ m n1 （10） 根据应力水平条件和加载时间的不同， 可以分 别得到相应的蠕变方程 （1） 当 σ＜σs时， 模型表示为广义 Kelvin 体， 可以描述衰减蠕变， 岩石处于稳定状态， 此时蠕变 方程为 2 总第 131 期煤 矿 开 采2016 年第 4 期 ChaoXing ε σ E1 σ E2 1 - e -E 2 η1 t （）（11） （2） 当 σ ≥σs， t＜tb时， 模型表达为经典的西 原模型， 可以描述稳定蠕变， 此时蠕变方程为 ε σ E1 σ E2 1 - e -E 2 η1 t （） σ - σs η2 t（12） （3） 当 σ ≥ σs， t≥tb时， 模型为非线性蠕变 模型， 可以描述蠕变全过程曲线， 反映岩石非稳定 阶段， 这一阶段将导致岩石失稳破坏， 此时蠕变方 程为 ε σ E1 σ E2 1 - e -E 2 η1 t （） σ - σs η2 t 1 - b n 1（ ） σ n t - tb（）［］ -n1（） aσm 1 - 1 - b n 1（ ） σ n t - tb（）［］ 1 n1 ［］ 1 - b n 1（ ） σ n t - tb（）［］ m n1 t（13） 1. 3 确定时间阈值 岩石材料流变破坏荷载和时间的关系如图2 所 示， 伴随时间增加， 破坏荷载呈减小趋势， 即流变 时间趋向无穷大对应的强度为长期强度 σ∞。 岩石 从开始进入加速蠕变阶段到岩石流变破坏经历的时 间可定义为加速历时。 图 2 σ 与 t 的关系 假定岩石的破坏荷载和时间的关系、 破坏荷载 和加速历时的关系用如下指数函数来表示 σ σ∞ 1 - σ0 σ∞ ■ ■ ■ ■ exp - α（）（14） σ σ∞ 1 - σ0 σ∞ ■ ■ ■ ■ exp Δtβ（）（15） 式中， σ0为岩石瞬时强度； σ∞为长期强度； t 为 破坏时间； Δt 为加速历时， α， β 为模拟参数。 则岩石加速流变的起始时间［11］可表示为 tb t - Δt β - α αβ ln σ0 - σ ∞ σ - σ∞ ■ ■ ■ ■ （16） 2 算例分析 2. 1 灰岩蠕变试验 本次试验灰岩岩样采自潘北煤矿， 其主要组成 为碳酸盐岩， 掺杂有黏土和粉砂， 密度 2210kg/ m3， 砂岩的波速为 2672 ～3570m/ s。 岩样取芯制成 高直径为 100mm50mm 的标准尺寸。 采用安徽 理工大学液压伺服加载试验系统对灰岩试件进行蠕 变试验， 将岩样分为 6 组， 分别进行加载， 得到如 图 3 所示的灰岩蠕变曲线 图 3 不同应力水平下灰岩的蠕变曲线 由图 3 可知， 分级加载岩样具有如下特征 （1） 瞬时施加应力后， 每级荷载都会产生瞬 时应变， 然后为蠕变变形， 其蠕变变形表现出非线 性， 应变率也有增加的趋势。 岩样存在一个蠕变应 力阈值 31. 2MPa （该阈值为岩样长期强度）， 阈值 右侧整体斜率小于阈值点左侧并整体偏向于应力 轴， 当应力小于阈值时蠕变很快衰减趋于零， 即只 有稳态蠕变。 （2） 岩样产生瞬时应变后， 微裂纹开始发生 随机分布， 损伤度逐渐增大， 在较高应力水平下随 时间增长开始产生非稳定蠕变变形， 低应力水平只 有瞬态蠕变和稳定蠕变阶段。 20MPa 载荷仅产生 8. 110 -4 的瞬时应变， 经过 5. 3h 后应变保持恒 定； 25MPa 载荷产生 9. 0 10 -4 的瞬时应变， 经 3. 9h 后应变发展为 1. 17 10 -3 后保持恒定不变。 随着应力水平的增大并大于起始蠕变应力阈值， 35MPa 载荷下产生 1. 110 -3 的瞬时蠕变， 然后进 入稳定蠕变阶段， 经 21h 后出现加速蠕变趋势。 应 力水平较高并大于时间阈值时将出现加速蠕变阶 段， 40MPa 载荷产生 1. 1710 -3 初始蠕变， 经过 15. 4h 后进入加速蠕变， 此载荷下加速蠕变阶段出 现相对较晚； 45MPa 载荷产生 1. 210 -3 瞬时蠕变， 经过短暂的稳定蠕变阶段， 立即进入加速蠕变阶 段， 并在较短时间内发生破坏， 此载荷下加速蠕变 阶段出现得较早。 2. 2 参数反演及验证 通过蠕变试验数据确定流变参数， 曲线拟合是 普遍适用的方法， 其中使用非线性最小二乘法进行 曲线 拟 合 应 用 得 最 为 广 泛， 即 对 给 定 数 据 点 ｛ （Xi， Yi） ｝ （i0， 1， ， m）， 在取定的函数 3 付腾飞等 基于西原模型的非线性损伤蠕变模型2016 年第 4 期 ChaoXing Φ 中， 求拟合函数 p （x） ∈Φ， 通过运用最小二 乘法使误差的平方和 ∑ ［p （Xi） - Yi］2最小。 Matlab 软件可以用于各种线性或非线性曲线拟合， 本文利用其参数反演工具箱， 读入轴向蠕变应变、 及相对应的蠕变时间， 根据公式 （11）、（12）、 （13） 对灰岩蠕变曲线进行拟合。 假定 Φ 是文中提 出的非线性蠕变曲线方程， 是自变量 t 和待定参数 B 的已知函数 Φf （t， B）， 已给出 （Φ， T） 的 n 对观测值 （εk， tk） （k1， 2， ， n）， 要求待 定参数 B， 通过使 Q ∑nk ［Φk-f （tk， B） ］2 为最小可得［12］， 得到如表 1 所示的各参数拟合值。 取外 载 荷 σ 40MPa 时， σ0 48. 3MPa， σ∞ 31. 2MPa， 由公式 （16） 可以得到 α0. 03， β 1. 04， tb15. 7h。 表 1 非线性损伤蠕变模型的拟合参数 σ/ MPa E1/ GPaE2/ GPa η1/ （GPah） η2/ （GPah） abmnR2 2015. 52175. 43182. 720. 9937 2516. 64201. 13194. 230. 9914 3016. 69223. 62199. 310. 9893 3516. 55543. 10421. 50629. 10. 9871 4017. 022214. 90973. 101987. 61. 710 -6 5. 810 -6 1. 92. 20. 9896 4516. 901841. 501236. 202510. 62. 410 -6 8. 910 -6 1. 42. 70. 9864 参数拟合结果显示出 新非线性蠕变模型的拟 合结果与实际试验结果的相关性系数 R2都较高， 为 0. 95 以上， 说明新非线性损伤蠕变模型可以较 好地描述岩石的 3 种蠕变情况， 拟合结果如图 4 所 示。 图 4 蠕变理论曲线和试验曲线 3 结 论 （1） 将岩石材料的强度及蠕变特性与内部损 伤机理相联系， 考虑微观的影响。 引入 Kachanov 损伤理论到西原模型蠕变表达式中， 较准确地解释 岩石材料的蠕变本质， 从微观结构的岩石损伤考虑 对蠕变加速的影响， 使材料蠕变特性更符合实际。 （2） 从岩石流变破坏载荷与破坏时间的关系 考虑， 引入反映加速蠕变起始点的时间阈值， 当达 到时间阈值 t＞tb时， 岩石开始出现损伤劣化， 反 映岩石蠕变的加速阶段。 （3） 通过灰岩流变试验， 运用 Matlab 软件对 建立的非线性损伤蠕变模型进行了反演， 得到岩石 非线性蠕变模型的拟合曲线。 模型能较好地描述不 同水平应力下灰岩的蠕变曲线， 尤其是描述加速蠕 变阶段， 该非线性损伤蠕变模型可以为岩石流变失 稳现象提供理论依据。 ［参考文献］ ［1］ 楼志文 . 损伤力学基础 ［M］ . 西安 西安交通大学出版社， 1990. ［2］ 胡 波， 王宗林， 梁 冰， 等 . 岩石蠕变特性试验研究 ［J］. 实验力学， 2015， 30 （4） 438-446. ［3］ 陈沅江， 潘长良， 曹 平， 等 . 软岩流变的一种新力学模型 ［J］ . 岩土力学， 2003， 24 （2） 209-214. ［4］ 丁志坤， 吕爱钟 . 岩石粘弹性非定常蠕变方程的参数识别 ［J］ . 岩土力学， 2004， 25 （ S） 37-40. ［5］ 袁 林， 高召宁， 孟祥瑞 . 岩石的非线性西原流变模型及其 应用 ［J］ . 西北地震学报， 2012， 34 （4） 331-334. ［6］ 曹树刚， 边 金， 李 鹏 . 岩石蠕变本构关系及改进的西原 正夫模型 ［J］ . 岩石力学与工程学报， 2002， 21 （ 5） 632- 634. ［7］ 张强勇， 杨文东， 张建国， 等 . 变参数蠕变损伤本构模型及 其工程应用 ［J］ . 岩石力学与工程学报， 2009， 28 （4） 732-738. ［8］ 刘桃根， 王 伟， 吴斌华， 等 . 基于损伤力学的砂岩蠕变模 型研究与参数辨识 ［ J］ . 三峡大学学报 （自然科学版）， 2010， 32 （6） 55-60. ［9］ 孙 钧 . 岩土材料流变及其工程应用 ［M］ . 北京 中国建筑 工业出版社， 1999. ［10］ 冶小平， 李 曼， 孙 强 . 巷道软岩蠕变损伤模型 ［J］ . 煤 田地质与勘探， 2009， 37 （2） 44-46. ［11］ 李良权 . 基于西原模型的非线性黏弹塑性流变模型 ［J］ . 力 学学报， 2009， 41 （5） 671-680. ［12］ 李青麒 . 软岩蠕变参数的曲线拟合计算方法 ［J］ . 岩石力学 与工程学报， 1998， 17 （5） 559-564. ［13］ 赵宝云 . 岩石拉、 压蠕变特性研究及其在地下大空间洞室施 工控制中的应用 ［D］ . 重庆 重庆大学， 2011. ［14］ 伍国军， 陈卫忠， 曹俊杰， 等 . 工程岩体非线性蠕变损伤力 学模型及其应用 ［J］ . 岩石力学与工程学报， 2010， 29 （6） 1184-1191. ［责任编辑 李宏艳］ 4 总第 131 期煤 矿 开 采2016 年第 4 期 ChaoXing