一种建立岩石理论损伤演化模型的方法技术

本申请涉及一种建立岩石理论损伤演化模型的方法，涉及地下工程技术领域。方法包括：基于冈珀茨模型和待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式，确定破坏单元随应变发展的第二破坏单元数量表达式；基于第二破坏单元数量表达式和第一损伤变量表达式，确定待测岩石样本的损伤变量随应变发展的第二损伤变量表达式；对待测岩石样本进行单轴压缩试验，获取待测岩石样本的应力

【技术实现步骤摘要】
一种建立岩石理论损伤演化模型的方法


[0001]本申请涉及地下工程
，特别是涉及一种建立岩石理论损伤演化模型的方法。

技术介绍

[0002]目前，随着地下固体矿产资源开采由浅部逐步向深部空间进军，深部岩石工程面临的高地应力、高地温、高渗透压、强烈开挖扰动等复杂地质环境导致岩石内部缺陷不可逆发展，损伤程度非线性累积，极易发生岩爆、冲击地压等灾害事故。岩石内部损伤破裂过程伴随着能量耗散，而岩爆发生可能性及剧烈程度与岩石内部能量储存、耗散、释放过程密切相关。深入研究岩石损伤发展规律对深部地下工程围岩石稳定性监测、变形控制及灾害预警具有重要意义。
[0003]因此，亟需一种建立岩石理论损伤演化模型的方法。

技术实现思路

[0004]基于此，有必要针对上述技术问题，提供一种建立岩石理论损伤演化模型的方法。所述方法包括：基于冈珀茨模型和待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式，确定所述待测岩石样本中破坏单元随应变发展的第二破坏单元数量表达式；所述第一破坏单元数量表达式是基于所述待测岩石样本中破坏单元随应变发展的增长速度、所述待测岩石样本中破坏单元的增长率和所述待测岩石样本中上一个应变时刻的破坏单元数量构建的；基于所述第二破坏单元数量表达式和第一损伤变量表达式，确定待测岩石样本的损伤变量随应变发展的第二损伤变量表达式；所述第一损伤变量表达式是基于所述待测岩石样本中的破坏单元数量和所述待测岩石样本中的岩石单元总数量构建的；对所述待测岩石样本进行单轴压缩试验，获取所述待测岩石样本的应力r/>‑
应变曲线，并根据所述待测岩石样本的应力
‑
应变曲线，确定第三损伤变量表达式；根据所述第三损伤变量表达式，确定所述待测岩石样本的每个应变对应的损伤变量，并按照所述第二损伤变量表达式进行线性拟合，确定岩石理论损伤演化模型。
[0005]作为一种可选地实施方式，所述待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式为：dv/dε=r(v)v；其中，dv/dε为待测岩石样本中的破坏单元数量随应变发展的增长速度，v为待测岩石样本中上一个应变时刻的破坏单元数量，r(v)为待测岩石样本中上一个应变时刻的破坏单元数量v的函数，表示待测岩石样本中的破坏单元的增长率。
[0006]作为一种可选地实施方式，所述破坏单元随应变发展的第二破坏单元数量表达式为：
dv/dε=
‑
aln(v/k)v；其中，dv/dε为待测岩石样本中的破坏单元随着待测岩石样本的应变发展的增长速度，a为第一模型参数，v为待测岩石样本中的破坏单元数量，k为待测岩石样本中的岩石单元总数量。
[0007]作为一种可选地实施方式，所述第一损伤变量表达式为：D=v/k；其中，D为待测岩石样本的损伤变量，v为待测岩石样本中的破坏单元数量，k为待测岩石样本中的岩石单元总数量。
[0008]作为一种可选地实施方式，所述第二损伤变量表达式为：D=exp[
‑
exp(b
‑
aε)]；其中，D为待测岩石样本的损伤变量，a为第一模型参数，b为第二模型参数，ε为待测岩石样本的应变。
[0009]作为一种可选地实施方式，所述根据所述待测岩石样本的应力
‑
应变曲线，确定第三损伤变量表达式，包括：将所述应力
‑
应变曲线中的线弹性阶段的反向延长线与所述应力
‑
应变曲线中应变轴的交点处的应变确定为岩石裂纹闭合应变；将所述应力
‑
应变曲线中线弹性阶段的斜率确定为基准弹性模量；根据所述岩石裂纹闭合应变和所述基准弹性模量，确定所述第三损伤变量表达式。
[0010]作为一种可选地实施方式，所述第三损伤变量表达式为：D=1
‑
σ/[E0(ε
‑
ε
cc
)]；其中，D为待测岩石样本的损伤变量，σ为待测岩石样本的应力，E0为基准弹性模量，ε
cc
为岩石裂纹闭合应变，ε为待测岩石样本的应变。
[0011]第二方面，提供了一种计算机设备，包括存储器及处理器，所述存储器上存储有可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如第一方面任一项所述的方法步骤。
[0012]第三方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如第一方面任一项所述的方法步骤。
[0013]本申请提供了一种建立岩石理论损伤演化模型的方法，本申请的实施例提供的技术方案至少带来以下有益效果：基于冈珀茨模型和待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式，确定所述待测岩石样本中破坏单元随应变发展的第二破坏单元数量表达式；所述第一破坏单元数量表达式是基于所述待测岩石样本中破坏单元随应变发展的增长速度、所述待测岩石样本中破坏单元的增长率和所述待测岩石样本中上一个应变时刻的破坏单元数量构建的；基于所述第二破坏单元数量表达式和第一损伤变量表达式，确定待测岩石样本的损伤变量随应变发展的第二损伤变量表达式；所述第一损伤变量表达式是基于所述待测岩石样本中的破坏单元数量和所述待测岩石样本中的岩石单元总数量构建的；对所述待测岩石样本进行单轴压缩试验，获取所述待测岩石样本的应力
‑
应变曲线，并根据所述待测岩石样本的应力
‑
应变曲线，确定第三损伤变量表达式；根据所述第三损伤变量表达式，确定所述待测岩石样本的每个应变对应的损伤变量，并按照所述第二损
伤变量表达式进行线性拟合，确定岩石理论损伤演化模型。本申请通过单轴压缩实验，获取应力
‑
应变曲线，通过应力
‑
应变曲线可以确定待测岩石样本的岩石损伤过程具有非对称且发展速度先快后慢的特点。因此，基于描述单一物种演化过程的冈珀茨模型，确定岩石理论损伤演化模型，可以描述和预测岩石损伤演化全过程。
[0014]应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。
附图说明
[0015]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0016]图1为本申请实施例提供的一种建立岩石理论损伤演化模型的方法的流程图；图2为本申请实施例提供的一种待测岩石样本的结构示意图；图3为本申请实施例提供的一种应力
‑
应变曲线的示意图；图4为本申请实施例提供的一种损伤演化曲线的示意图；图5为本申请实施例提供的另一种建立岩石理论损伤演化模型的方法的流程图；图6为本申请实施例提供的另一种应力
‑
应变曲线的示意图；图7为本申请实施例提供的三种不同待测岩石样本的应力
‑
应变曲线的示意图；图8a为本申请实本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种建立岩石理论损伤演化模型的方法，其特征在于，所述方法包括：基于冈珀茨模型和待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式，确定所述待测岩石样本中破坏单元随应变发展的第二破坏单元数量表达式；所述第一破坏单元数量表达式是基于所述待测岩石样本中破坏单元随应变发展的增长速度、所述待测岩石样本中破坏单元的增长率和所述待测岩石样本中上一个应变时刻的破坏单元数量构建的；基于所述第二破坏单元数量表达式和第一损伤变量表达式，确定待测岩石样本的损伤变量随应变发展的第二损伤变量表达式；所述第一损伤变量表达式是基于所述待测岩石样本中的破坏单元数量和所述待测岩石样本中的岩石单元总数量构建的；对所述待测岩石样本进行单轴压缩试验，获取所述待测岩石样本的应力
‑
应变曲线，并根据所述待测岩石样本的应力
‑
应变曲线，确定第三损伤变量表达式；根据所述第三损伤变量表达式，确定所述待测岩石样本的每个应变对应的损伤变量，并按照所述第二损伤变量表达式进行线性拟合，确定岩石理论损伤演化模型。2.根据权利要求1所述的方法，其特征在于，所述待测岩石样本中破坏单元随应变发展的第一破坏单元数量表达式为：dv/dε=r(v)v；其中，dv/dε为待测岩石样本中的破坏单元数量随应变发展的增长速度，v为待测岩石样本中上一个应变时刻的破坏单元数量，r(v)为待测岩石样本中上一个应变时刻的破坏单元数量v的函数，表示待测岩石样本中的破坏单元的增长率。3.根据权利要求1所述的方法，其特征在于，所述待测岩石样本中破坏单元随应变发展的第二破坏单元数量表达式为：dv/dε=
‑
aln(v/k)v；其...

【专利技术属性】
技术研发人员：郭允朋，刘冬桥，李杰宇，凌凯，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：