一种岩石渗流-应力耦合损伤本构模型的构建方法

本发明涉及岩石力学与工程领域，尤其涉及一种岩石渗流-应力耦合损伤本构模型的构建方法。

背景技术：

1、岩石的渗流-应力耦合特性是岩石力学与工程领域的研究热点。近年来，针对岩石的渗流-应力耦合特性的研究层出不穷，提出了一些新的耦合模型、经验公式等，用来描述岩石的渗流-应力耦合特性。渗流-应力耦合作用下的本构模型是研究岩石在水环境耦合特性的有效手段之一。已有的模型中渗流-应力耦合作用的研究方法主要是经验公式拟合、渗透系数法等，这些方法简单、直观，与试验结果吻合较好，但普适性不强。因此，通过引入渗流-应力耦合变量，模拟岩石内渗流与应力的相互作用，可以较为有效的实现渗透水对于岩石内部应力场的影响和应力作用下的变形行为对岩石内部的渗透演化规律的影响，从而实现岩石的渗流－应力耦合效应。

技术实现思路

1、发明目的：本发明旨在提供一种岩石渗流-应力耦合损伤本构模型的构建方法，该构建方法通过引入渗流-应力耦合变量，建立耦合变量与损伤和塑性变形之间的关系，得到岩石渗流-应力耦合损伤本构模型。

2、技术方案：本发明所述的岩石渗流-应力耦合损伤本构模型的构建方法，包括以下步骤：

3、(1)引入渗流-应力耦合变量，用于表征水压作用下的岩石渗流-应力耦合特性，得到耦合变量ξ；

4、(2)引入加载阶段的损伤演化准则，用于表征岩石加载过程中的损伤演化，得到损伤变量d；

5、(3)将耦合变量和损伤变量引入塑性屈服面方程和塑性势能方程中，得到渗流-应力耦合本构的屈服面方程和塑性势能方程；

6、(4)依据一致性原理，根据塑性屈服面、塑性势能函数，列出塑性和损伤的一致性条件，基于塑性流动法则，计算损伤乘子dd、渗流-应力耦合乘子dξ和耦合作用下的塑性乘子λs；

7、(5)依据岩石渗流-应力耦合作用下模量退化原理，引入渗流-应力耦合的耦合变量、损伤变量以及损伤乘子dd、渗流-应力耦合乘子dξ和耦合作用下的塑性乘子λs，计算渗流-应力耦合作用下的刚度张量，得到渗流-应力耦合损伤本构模型；

8、(6)根据岩石三轴压缩试验结果和岩石物理特性测试结果，确定岩石渗流-应力耦合本构模型参数。

9、进一步的，步骤(1)中，渗流-应力耦合变量ξ为

10、

11、式中，as为岩样的横截面面积，单位为m2；ls为岩样的高度，单位为m；φ为孔隙率；δvl为瞬时时刻渗透水浸入岩样的体积，单位为m3；

12、瞬时时刻渗透水浸入岩样的体积δvl为

13、

14、式中，pl为岩样上下两端施加的渗透压差，单位为pa；δt为时间间隔，单位为s；μl为水的动力粘滞系数，单位为pa·s，优选20℃时水的动力粘滞系数为1×10-3pa·s；k为岩石的渗透率，单位为m2；

15、岩石的渗透率k为

16、

17、式中，k0为初始渗透率，φ0为初始孔隙率，φ为孔隙率；

18、孔隙率φ为

19、

20、

21、式中，d为损伤变量，md为表征岩石扩容速率的参数，p、q分别为岩石的平均应力与偏应力，为扩容孔隙率的最大阈值，φdil为扩容孔隙率，dεp为塑性应变增量。

22、进一步的，步骤(2)中，瞬时损伤fd为

23、

24、式中，ydp为岩石损伤演化的驱动力，dc为损伤变量的最大临界值，bd为控制损伤演化速率的参数。

25、进一步的，步骤(3)中，渗流-应力耦合损伤本构的屈服面方程f和塑性势能方程g为

26、f＝q2+a0(1-ξ)αp(p-c0)p0＝0            (7)

27、

28、式中，p为平均应力，q为偏应力，a0和c0是根据平均应力p与偏应力q的平面上的轨迹而获得的模型参数，p0为归一化参量，η为岩石压缩和扩容区域边界线的斜率，i0为塑性势面与平均应力的交点值。

29、进一步的，步骤(4)中，损伤乘子dd为

30、

31、式中，εp为塑性应变，γp为塑性剪切变形值，σ为应力张量；

32、渗流-应力耦合乘子dξ为

33、

34、耦合作用下塑性乘子λs为

35、

36、渗流-应力耦合的塑性硬化模量hξd为

37、

38、进一步的，步骤(5)中，岩石渗流-应力耦合损伤本构模型表示为

39、

40、式中，为岩石初始的四阶弹性张量，εe为应变张量，为渗流-应力耦合弹性刚度矩阵；依据岩石渗流-应力耦合作用下模量退化原理，岩石渗流-应力耦合损伤本构增量方程为

41、

42、式中，为四阶渗流-应力耦合柔度矩阵，为四阶损伤柔度矩阵。

43、有益效果：本发明与现有技术相比，其显著优点是：本发明将渗流-应力耦合变量ξ和损伤变量d引入到本构模型中，不仅可以更好的描述岩石的的塑性变形、损伤演化、对压力依赖、压缩到扩容的转变以及峰前塑性硬化与峰后的应变软化等现象，还可以准确描述渗流-应力耦合引起岩样变形和强度特性的改变，能够更好的诠释岩石在渗透水压力下的渗流-应力耦合特性，且不同渗透水压下的数值模拟结果与试验结果基本一致，克服了传统模型无法准确模拟渗流-应力耦合特性等问题，对准确评价岩石工程在渗透水压作用下的安全与稳定性问题具有重要的科学意义和参考价值。

技术特征：

1.一种岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，包括以下步骤：

2.根据权利要求1所述岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，步骤(1)中，渗流-应力耦合变量ξ为

3.根据权利要求2所述岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，所述瞬时时刻渗透水浸入岩样的体积δvl为

4.根据权利要求3所述岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，所述岩石的渗透率k为

5.根据权利要求4所述岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，所述孔隙率φ为

6.根据权利要求5所述岩石三轴流变损伤本构模型的构建方法，其特征在于，步骤(2)中，瞬时损伤fd为

7.根据权利要求6所述岩石三轴流变损伤本构模型的构建方法，其特征在于，步骤(3)中，渗流-应力耦合损伤本构的屈服面方程f和塑性势能方程g为

8.根据权利要求7所述岩石三轴流变损伤本构模型的构建方法，其特征在于，步骤(4)中，损伤乘子dd为

9.根据权利要求8所述岩石三轴流变损伤本构模型的构建方法，其特征在于，所述渗流-应力耦合的塑性硬化模量hξd为

10.根据权利要求9所述的一种岩石渗流-应力耦合损伤本构模型的构建方法，其特征在于，步骤(5)中，岩石渗流-应力耦合损伤本构模型表示为

技术总结
本发明公开了一种岩石渗流‑应力耦合损伤本构模型的构建方法，该构建方法基于岩石渗流‑应力耦合的基本理论，通过对岩石渗流‑应力耦合中内部结构的孔隙率变化特征的分析，提出渗流‑应力耦合变量的概念；通过建立耦合变量与损伤和塑性变形之间的关系，将耦合变量引入到弹塑性损伤本构框架中来表征渗流‑应力耦合的特性，对准确评价重大工程中岩石的安全与稳定性具有参考价值。

技术研发人员：刘琳
受保护的技术使用者：河海大学
技术研发日：
技术公布日：2024/1/15