CN114062132A - 一种岩石单轴压缩加速蠕变起始时间的预测方法 - Google Patents

一种岩石单轴压缩加速蠕变起始时间的预测方法 Download PDF

Info

Publication number
CN114062132A
CN114062132A CN202111327972.8A CN202111327972A CN114062132A CN 114062132 A CN114062132 A CN 114062132A CN 202111327972 A CN202111327972 A CN 202111327972A CN 114062132 A CN114062132 A CN 114062132A
Authority
CN
China
Prior art keywords
creep
rock
time
point
viscosity coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111327972.8A
Other languages
English (en)
Other versions
CN114062132B (zh
Inventor
王军保
王晓鹏
张玉伟
宋战平
冯世进
王彤
姜恒坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Architecture and Technology
Original Assignee
Xian University of Architecture and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Architecture and Technology filed Critical Xian University of Architecture and Technology
Priority to CN202111327972.8A priority Critical patent/CN114062132B/zh
Publication of CN114062132A publication Critical patent/CN114062132A/zh
Application granted granted Critical
Publication of CN114062132B publication Critical patent/CN114062132B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001—Type of application of the stress
    • G01N2203/0003—Steady
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014—Type of force applied
    • G01N2203/0016—Tensile or compressive
    • G01N2203/0019—Compressive
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0071—Creep
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/025—Geometry of the test
    • G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种岩石单轴压缩加速蠕变起始时间的预测方法，本发明根据岩石单轴压缩蠕变试验结果，确定改进Maxwell模型的参数；将粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线绘制在双纵轴坐标系中，确定衰减蠕变结束点A，过A点作垂直线，垂直线与粘滞系数随时间的变化规律曲线相交于B点；过B点作水平线，水平线与粘滞系数随时间的变化规律曲线相交于C点，C点对应的时间即为岩石加速蠕变起始时间。本方法通过改进Maxwell模型能够准确预测岩石单轴压缩加速蠕变起始时间，实现了对岩石加速蠕变起始时间的准确预测，能够在岩石加速蠕变时提前预警，大大降低了了围岩坍塌对人民生命财产危害的发生。

Description

一种岩石单轴压缩加速蠕变起始时间的预测方法
技术领域
本发明属于岩石工程技术领域，具体涉及一种岩石单轴压缩加速蠕变起始时间的预测方法。
背景技术
岩石完整的蠕变过程通常包含衰减蠕变、稳态蠕变和加速蠕变等三个阶段。尽管根据Kachanov蠕变损伤定律可以预测岩石蠕变破坏时间，但对于地下工程而言，为了确保人民生命和财产安全，不允许围岩发生坍塌破坏。一旦围岩变形进入加速蠕变阶段，须立即对其进行强支护，以确保围岩稳定。因此，准确的确定和预测岩石加速蠕变起始时间具有重要意义。而目前，关于岩石加速蠕变起始时间的确定和预测方法的研究较少。
发明内容
本发明的目的在于克服上述不足，提供一种岩石单轴压缩加速蠕变起始时间的预测方法，以实现对岩石加速蠕变起始时间的准确确定和预测。
为了达到上述目的，本发明包括以下步骤：
步骤一，构建改进Maxwell模型，并确定改进Maxwell模型的参数；
步骤二，将岩石单轴压缩蠕变试验的蠕变参数代入改进Maxwell模型中粘滞系数随时间变化的函数表达式，并根据该函数表达式将粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线绘制在双纵轴坐标系中；
步骤三，在岩石蠕变曲线上确定衰减蠕变结束点A，过A点作一条垂直线，该垂直线与粘滞系数随时间的变化规律曲线的左半部分相交于B点；
步骤四，过B点作一条水平线，该水平线与粘滞系数随时间的变化规律曲线的右半部分相交于C点，C点对应的时间即为岩石加速蠕变起始时间。
步骤二中，绘制粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线的双纵轴坐标系中有一条横轴和两条纵轴，横轴为时间轴，纵轴为应变轴和粘滞系数轴。
步骤三中，衰减蠕变结束点A为岩石衰减蠕变阶段和稳态蠕变阶段的分界点，A点之前岩石处于衰减蠕变阶段，A点之后岩石处于稳态蠕变阶段。
步骤四中，B点和C点的粘滞系数相同，B点和C点对应的时间不同。
步骤一中，构建改进Maxwell模型的具体方法如下：
第一步，采用单级加载方式，对岩石试件开展不同轴向应力σ1、σ2…σi-1、σi作用下的单轴压缩蠕变试验，获得岩石在对应轴向应力下的蠕变曲线，其中σ1＜σ2＜…＜σi-1＜σi；
第二步，根据出现加速蠕变的岩石单轴压缩全过程蠕变曲线，计算岩石在不同时间时的蠕变速率，并绘制岩石蠕变速率随时间的变化规律曲线；
第三步，假设岩石蠕变速率符合元件组合蠕变模型中粘性元件的本构方程，根据岩石蠕变速率随时间的变化规律曲线和粘性元件的本构方程，推测岩石蠕变全过程中的粘滞系数随时间的变化规律，并确定粘滞系数随时间变化的函数表达式；
第四步，以粘滞系数随时间变化的函数表达式建立改进Maxwell模型的本构方程；
第五步，对改进Maxwell模型的本构方程进行积分，并结合Kachanov蠕变损伤定律，获得改进Maxwell模型的蠕变方程。
步骤一中，根据岩石单轴压缩蠕变试验结果，确定改进Maxwell模型的参数。
第三步中，粘滞系数随时间变化的函数表达式为：
Figure BDA0003347628350000021
式中：t为时间，ηFN(t)为t时刻时的岩石粘滞系数，tF为岩石蠕变破坏时间，m、n、ηFN均为材料参数。
步骤四中，改进Maxwell模型的本构方程为：
Figure BDA0003347628350000031
其中，
Figure BDA0003347628350000032
为轴向应力随时间的变化速率，E为岩石的弹性模量，t为时间，ηFN(t)为t时刻时的岩石粘滞系数，tF为岩石蠕变破坏时间，m、n和ηFN均为材料参数。
岩石的弹性模量E、材料参数m、材料参数n和材料参数ηFN的变化规律分别满足下式：
E(σ)＝a1exp(-b1σ)
Figure BDA0003347628350000033
m(σ)＝a3exp(b3σ)
n(σ)＝a4exp(b4σ)
其中，a1、b1、a2、b2、a3、b3、a4和b4均为拟合参数。
步骤五中，Kachanov蠕变损伤定律为：
Figure BDA0003347628350000034
改进Maxwell模型的蠕变方程为：
Figure BDA0003347628350000035
其中，A和v均为模型参数。
与现有技术相比，本发明根据岩石单轴压缩蠕变试验结果，确定改进Maxwell模型的参数；将粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线绘制在双纵轴坐标系中，确定衰减蠕变结束点A，过A点作一条垂直线，该垂直线与粘滞系数随时间的变化规律曲线的左半部分相交于B点；过B点作一条水平线，该水平线与粘滞系数随时间的变化规律曲线的右半部分相交于C点，C点对应的时间即为岩石加速蠕变起始时间。本方法通过改进Maxwell模型能够准确预测岩石单轴压缩加速蠕变起始时间，实现了对岩石加速蠕变起始时间的准确预测，能够在岩石加速蠕变时提前预警，大大降低了了围岩坍塌对人民生命财产危害的发生。
附图说明
图1为本发明中不同轴向应力下的岩石单轴压缩蠕变曲线图；
图2为本发明中岩石蠕变全过程的蠕变速率随时间的变化规律曲线图；
图3为本发明中岩石蠕变全过程的粘滞系数随时间的变化规律示意图；
图4为本发明中岩石弹性模量E随轴向应力的变化规律曲线图；
图5为本发明中参数ηFN随轴向应力的变化规律曲线图；
图6为本发明中参数m随轴向应力的变化规律曲线图；
图7为本发明中参数n随轴向应力的变化规律曲线图；
图8为本发明中不同轴向应力下岩石加速蠕变起始时间的确定过程示意图；其中：(a)为26MPa；(b)为24MPa；(c)为17.5MPa；(d)为9.5MPa。
具体实施方式
下面结合附图对本发明做进一步说明。
本发明包括开展岩石单轴压缩蠕变试验，绘制岩石蠕变速率随时间的变化规律曲线，确定粘滞系数随时间变化的函数表达式，建立改进Maxwell模型的本构方程，获得改进Maxwell模型的蠕变方程，确定改进Maxwell模型的参数，在双纵轴坐标系中绘制粘滞系数随时间的变化规律曲线和对应轴向应力下的岩石蠕变曲线，确定岩石加速蠕变起始时间八大步骤。
步骤一，开展岩石单轴压缩蠕变试验；
采用单级加载方式，对直径为50mm、高度为100mm的圆柱形盐岩试件开展不同轴向应力σ1＝9.5、σ2＝17.5、σ3＝24和σ4＝26MPa下的单轴压缩蠕变试验，获得盐岩在对应轴向应力下的蠕变曲线，参见图1。
参见图1，当轴向应力σ1＝9.5MPa和σ2＝17.5MPa时，盐岩总应变由加载过程产生的与时间无关的瞬时弹性应变和随时间逐渐增大的蠕变应变两部分组成，且盐岩在试验时间内均未发生加速蠕变，蠕变曲线仅包含衰减蠕变和稳态蠕变两个阶段；当轴向应力σ3＝24MPa和σ4＝26MPa时，盐岩总应变也由加载过程产生的与时间无关的瞬时弹性应变和随时间逐渐增大的蠕变应变两部分组成，但盐岩在试验时间内发生了加速蠕变，蠕变曲线包含衰减蠕变、稳态蠕变和加速蠕变三个阶段，包含三个阶段的蠕变曲线即为全过程蠕变曲线。
步骤二，绘制岩石蠕变速率随时间的变化规律曲线；
根据轴向应力σ3＝24MPa和σ4＝26MPa时的盐岩单轴压缩全过程蠕变曲线，计算盐岩在不同时间时的蠕变速率，并绘制盐岩蠕变速率随时间的变化规律曲线。参见图2，为轴向应力σ4＝26MPa时盐岩蠕变速率随时间的变化规律曲线。
参见图2，盐岩蠕变速率随时间的变化规律可划分为快速减小到基本不变到—快速增大三个阶段。
步骤三，确定粘滞系数随时间变化的函数表达式；
岩石元件组合蠕变模型中粘性元件的本构方程为
Figure BDA0003347628350000051
式中：η为粘滞系数，
Figure BDA0003347628350000052
为蠕变速率。
假设岩石蠕变速率符合元件组合蠕变模型中粘性元件的本构方程，且由于岩石蠕变过程中轴向应力σ保持恒定，则根据公式(1)可知，岩石蠕变速率和粘滞系数呈倒数关系。
同时，结合盐岩蠕变全过程中其蠕变速率随时间变化呈现出快速减小到基本不变到快速增大的三阶段性，参见图2，可推测盐岩蠕变全过程中其粘滞系数随时间的变化规律可划分为快速增大到基本不变到快速减小三个阶段。
在此基础上，经过分析，确定了盐岩蠕变过程中粘滞系数随时间变化的函数表达式为
Figure BDA0003347628350000053
式中：t为时间，ηFN(t)为t时刻时的岩石粘滞系数，tF为岩石蠕变破坏时间，m、n、ηFN均为材料参数。
进一步的，取ηFN＝500MPa·h，tF＝80h，m＝0.9和n＝0.02，参见图3，为根据公式(2)做出的盐岩蠕变全过程的粘滞系数随时间的变化规律示意图。
参见图3，根据公式(2)做出的盐岩蠕变全过程的粘滞系数随时间的变化规律可划分为快速增大到基本不变到快速减小三个阶段，与推测结果一致。
步骤四，建立改进Maxwell模型的本构方程；
岩石元件组合蠕变模型中的常规Maxwell模型的本构方程为
Figure BDA0003347628350000061
式中：
Figure BDA0003347628350000062
为轴向应力随时间的变化速率，E为岩石的弹性模量。
将盐岩蠕变过程中粘滞系数随时间变化的函数表达式(公式(2))代替常规Maxwell模型本构方程(公式(3))中的粘滞系数η，则可建立改进Maxwell模型的本构方程为
Figure BDA0003347628350000063
步骤五，获得改进Maxwell模型的蠕变方程；
对改进Maxwell模型的本构方程(公式(4))进行积分，可得改进Maxwell模型的蠕变方程为
Figure BDA0003347628350000064
公式(5)中，蠕变破坏时间tF可根据Kachanov蠕变损伤定律确定，即
Figure BDA0003347628350000065
将公式(6)代入公式(5)，可得改进Maxwell模型蠕变方程的最终表达式为
Figure BDA0003347628350000066
步骤六，确定改进Maxwell模型的参数；
改进Maxwell模型中共有岩石的弹性模量E、模型参数A、模型参数ν、材料参数ηFN、材料参数m和材料参数n六个参数需要确定。
岩石的弹性模量E根据岩石加载完成后产生的瞬时弹性应变，利用公式(7)等号右边第一项确定；
模型参数A和模型参数ν根据最后两级轴向应力σi-1和σi作用下岩石发生蠕变破坏的时间tF(i-1)和tF(i)，利用公式(6)确定；
参见图1，在最后两级轴向应力σ3＝24MPa和σ4＝26MPa作用下，盐岩在试验时间内发生了蠕变破坏，对应的破坏时间分别为tF(3)＝98.60h和tF(4)＝28.62h。将(24，98.60)和(26，28.62)两组数据代入公式(6)，可建立如下二元一次方程组：
Figure BDA0003347628350000071
解该方程组，可得A＝2.89×10-25，ν＝15.45。
岩石的弹性模量E、模型参数A和模型参数ν确定后，材料参数ηFN、材料参数m和材料参数n根据不同轴向应力下的蠕变试验数据利用曲线拟合法反演确定。
进一步的，根据不同轴向应力下的蠕变试验数据利用曲线拟合法反演确定材料参数ηFN、材料参数m和材料参数n的方法如下：
以待反演的材料参数ηFN、材料参数m和材料参数n作为设计变量X，即X＝{ηFN，m，n}；
建立目标函数Y，取
Figure BDA0003347628350000072
式中：N为试验组数，
Figure BDA0003347628350000073
为t时刻计算蠕变应变值，
Figure BDA0003347628350000074
为t时刻试验实测蠕变应变值。
设定目标函数的控制精度并进行参数迭代求解，若目标函数满足精度要求，则停止迭代，输出计算结果；若不满足，则继续迭代，直到满足精度要求为止。
按照上述方法，基于盐岩单轴压缩蠕变试验结果，利用数学优化分析软件，采用曲线拟合法反演确定盐岩在不同轴向应力下的材料参数ηFN、材料参数m和材料参数n。
表1为不同轴向应力下蠕变参数的反演确定结果。
表1不同轴向应力下蠕变参数反演确定结果
Figure BDA0003347628350000081
可以看出，岩石的弹性模量E及材料参数ηFN、材料参数m和材料参数n随轴向应力变化而变化。参见图4，为岩石的弹性模量E随轴向应力的变化规律；参见图5，为材料参数ηFN随轴向应力的变化规律；参见图6，为材料参数m随轴向应力的变化规律；参见图7，为材料参数n随轴向应力的变化规律。总体而言，随轴向应力增大，岩石的弹性模量E和参数ηFN逐渐减小，而材料参数m和材料参数n逐渐增大。为了预测盐岩在不同轴向应力下的蠕变行为，需要建立岩石的弹性模量E及材料参数ηFN、材料参数m和材料参数n随轴向应力变化的定量关系表达式。
经过拟合分析，岩石的弹性模量E及材料参数ηFN、材料参数m和材料参数n随轴向应力的变化规律分别满足公式(8)－公式(11)：
E(σ)＝15711.46exp(-0.0972σ) (8)
ηFN(σ)＝114675σ-1.36 (9)
m(σ)＝0.1239exp(0.0302σ) (10)
n(σ)＝0.0019exp(0.1377σ) (11)
步骤七，在双纵轴坐标系中绘制粘滞系数随时间的变化规律曲线和对应轴向应力下的岩石蠕变曲线；
将A＝2.89×10-25、ν＝15.45代入公式(6)，并将公式(6)和公式(9)、公式(10)、公式(11)代入公式(2)，即可预测任意轴向应力下岩石蠕变过程中粘滞系数随时间的变化规律。
将不同轴向应力下岩石蠕变过程中粘滞系数随时间的变化规律曲线和对应轴向应力下的岩石蠕变曲线绘制在双纵轴坐标系中。参见图8，为轴向应力分别为26、24、17.5和9.5MPa时盐岩加速蠕变起始时间的确定过程示意图。
参见图8，双纵轴坐标系只有一条横轴，即时间轴；有两条纵轴，分别为应变轴和粘滞系数轴。这表明两条曲线共用一条横轴，即时间轴。
步骤八，确定岩石加速蠕变起始时间；
参见图8，在不同轴向应力下的盐岩蠕变曲线上确定衰减蠕变结束点A。衰减蠕变结束点A为岩石衰减蠕变阶段和稳态蠕变阶段的分界点，即A点之前岩石处于衰减蠕变阶段，A点之后岩石处于稳态蠕变阶段。
过A点作一条垂直线，该垂直线与粘滞系数随时间的变化规律曲线的左半部分相交于B点。
过B点作一条水平线，该水平线与粘滞系数随时间的变化规律曲线的右半部分相交于C点，C点对应的时间即为岩石加速蠕变起始时间。
可以看出，当轴向应力分别为26、24、17.5和9.5MPa时，预测的盐岩加速蠕变起始时间分别大约为25h、87.5h、13100h和164856661h。同时，当轴向应力为26和24MPa时，由于试验持续时间内盐岩发生了加速蠕变，故可以利用这两组试验结果验证本发明方法预测结果的准确性。很显然，轴向应力为26MPa时，盐岩约在25h时进入加速蠕变阶段；而轴向应力为24MPa时，盐岩约在87.5h时进入加速蠕变阶段。因此，本发明提供的方法可以准确的预测岩石加速蠕变起始时间。
此外，衰减蠕变结束时和加速蠕变开始时，岩石具有相同的粘滞系数，这是本发明方法的关键。即B点和C点具有相同的粘滞系数，但该粘滞系数对应于两个不同的时间，分别为衰减蠕变结束的时间和加速蠕变开始的时间。因此，只要能够通过试验获得岩石衰减蠕变结束的时间，即可利用本发明提供的方法预测岩石的加速蠕变起始时间。

Claims (10)

1.一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，包括以下步骤：
步骤一，构建改进Maxwell模型，并确定改进Maxwell模型的参数；
步骤二，将岩石单轴压缩蠕变试验的蠕变参数代入改进Maxwell模型中粘滞系数随时间变化的函数表达式，并根据该函数表达式将粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线绘制在双纵轴坐标系中；
步骤三，在岩石蠕变曲线上确定衰减蠕变结束点A，过A点作一条垂直线，该垂直线与粘滞系数随时间的变化规律曲线的左半部分相交于B点；
步骤四，过B点作一条水平线，该水平线与粘滞系数随时间的变化规律曲线的右半部分相交于C点，C点对应的时间即为岩石加速蠕变起始时间。
2.根据权利要求1所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤二中，绘制粘滞系数随时间的变化规律曲线与对应轴向应力下的岩石蠕变曲线的双纵轴坐标系中有一条横轴和两条纵轴，横轴为时间轴，纵轴为应变轴和粘滞系数轴。
3.根据权利要求1所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤三中，衰减蠕变结束点A为岩石衰减蠕变阶段和稳态蠕变阶段的分界点，A点之前岩石处于衰减蠕变阶段，A点之后岩石处于稳态蠕变阶段。
4.根据权利要求1所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤四中，B点和C点的粘滞系数相同，B点和C点对应的时间不同。
5.根据权利要求1所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤一中，构建改进Maxwell模型的具体方法如下：
第一步，采用单级加载方式，对岩石试件开展不同轴向应力σ1、σ2…σi-1、σi作用下的单轴压缩蠕变试验，获得岩石在对应轴向应力下的蠕变曲线，其中σ1＜σ2＜…＜σi-1＜σi；
第二步，根据出现加速蠕变的岩石单轴压缩全过程蠕变曲线，计算岩石在不同时间时的蠕变速率，并绘制岩石蠕变速率随时间的变化规律曲线；
第三步，假设岩石蠕变速率符合元件组合蠕变模型中粘性元件的本构方程，根据岩石蠕变速率随时间的变化规律曲线和粘性元件的本构方程，推测岩石蠕变全过程中的粘滞系数随时间的变化规律，并确定粘滞系数随时间变化的函数表达式；
第四步，以粘滞系数随时间变化的函数表达式建立改进Maxwell模型的本构方程；
第五步，对改进Maxwell模型的本构方程进行积分，并结合Kachanov蠕变损伤定律，获得改进Maxwell模型的蠕变方程。
6.根据权利要求1所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤一中，根据岩石单轴压缩蠕变试验结果，确定改进Maxwell模型的参数。
7.根据权利要求5所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，第三步中，粘滞系数随时间变化的函数表达式为：
Figure FDA0003347628340000021
式中：t为时间，ηFN(t)为t时刻时的岩石粘滞系数，tF为岩石蠕变破坏时间，m、n、ηFN均为材料参数。
8.根据权利要求5所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤四中，改进Maxwell模型的本构方程为：
Figure FDA0003347628340000022
其中，
Figure FDA0003347628340000023
为轴向应力随时间的变化速率，E为岩石的弹性模量，t为时间，ηFN(t)为t时刻时的岩石粘滞系数，tF为岩石蠕变破坏时间，m、n和ηFN均为材料参数。
9.根据权利要求8所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，岩石的弹性模量E、材料参数m、材料参数n和材料参数ηFN的变化规律分别满足下式：
E(σ)＝a1exp(-b1σ)
Figure FDA0003347628340000031
m(σ)＝a3exp(b3σ)
n(σ)＝a4exp(b4σ)
其中，a1、b1、a2、b2、a3、b3、a4和b4均为拟合参数。
10.根据权利要求5所述的一种岩石单轴压缩加速蠕变起始时间的预测方法，其特征在于，步骤五中，Kachanov蠕变损伤定律为：
Figure FDA0003347628340000032
改进Maxwell模型的蠕变方程为：
Figure FDA0003347628340000033
其中，A和v均为模型参数。
CN202111327972.8A 2021-11-10 2021-11-10 一种岩石单轴压缩加速蠕变起始时间的预测方法 Active CN114062132B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111327972.8A CN114062132B (zh) 2021-11-10 2021-11-10 一种岩石单轴压缩加速蠕变起始时间的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111327972.8A CN114062132B (zh) 2021-11-10 2021-11-10 一种岩石单轴压缩加速蠕变起始时间的预测方法

Publications (2)

Publication Number Publication Date
CN114062132A true CN114062132A (zh) 2022-02-18
CN114062132B CN114062132B (zh) 2024-10-18

Family

ID=80274836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111327972.8A Active CN114062132B (zh) 2021-11-10 2021-11-10 一种岩石单轴压缩加速蠕变起始时间的预测方法

Country Status (1)

Country Link
CN (1) CN114062132B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116011191A (zh) * 2022-12-13 2023-04-25 广西大学 一种表征真三轴下岩石蠕变启动和加速的模型构建方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012202908A (ja) * 2011-03-28 2012-10-22 Hiroyuki Sato クリープ曲線およびクリープ寿命の予測方法
JP2017223464A (ja) * 2016-06-13 2017-12-21 株式会社日立パワーソリューションズ クリープ損傷評価方法
CN110274835A (zh) * 2019-07-12 2019-09-24 中南大学 一种改进Burgers岩石剪切蠕变模型的方法
CN110631908A (zh) * 2019-09-17 2019-12-31 西安建筑科技大学 一种岩石单轴压缩全过程蠕变损伤模型的构建方法
CN112485114A (zh) * 2020-11-23 2021-03-12 河海大学 一种预测混凝土特慢蠕变的方法
CN112946086A (zh) * 2021-02-02 2021-06-11 江西理工大学 一种岩石蠕变破坏微裂纹成核起始时间辨识方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012202908A (ja) * 2011-03-28 2012-10-22 Hiroyuki Sato クリープ曲線およびクリープ寿命の予測方法
JP2017223464A (ja) * 2016-06-13 2017-12-21 株式会社日立パワーソリューションズ クリープ損傷評価方法
CN110274835A (zh) * 2019-07-12 2019-09-24 中南大学 一种改进Burgers岩石剪切蠕变模型的方法
CN110631908A (zh) * 2019-09-17 2019-12-31 西安建筑科技大学 一种岩石单轴压缩全过程蠕变损伤模型的构建方法
CN112485114A (zh) * 2020-11-23 2021-03-12 河海大学 一种预测混凝土特慢蠕变的方法
CN112946086A (zh) * 2021-02-02 2021-06-11 江西理工大学 一种岩石蠕变破坏微裂纹成核起始时间辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MINGZHENG WANG等: "A grain-based time-to-failure creep model for brittle rocks", COMPUTERS AND GEOTECHNICS, 29 November 2019 (2019-11-29), pages 1 - 13 *
王宇等: "砂岩单轴压缩蠕变试验及对西原本构模型的改进", 兰州工业学院学报, 31 October 2020 (2020-10-31), pages 10 - 14 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116011191A (zh) * 2022-12-13 2023-04-25 广西大学 一种表征真三轴下岩石蠕变启动和加速的模型构建方法
CN116011191B (zh) * 2022-12-13 2024-05-10 广西大学 一种表征真三轴下岩石蠕变启动和加速的模型构建方法

Also Published As

Publication number Publication date
CN114062132B (zh) 2024-10-18

Similar Documents

Publication Publication Date Title
CN110631908B (zh) 一种岩石单轴压缩全过程蠕变损伤模型的构建方法
Li et al. Numerical simulation of incremental sheet forming based on GTN damage model
Malag et al. Problems determining of the mechanical properties of metallic materials from the tensile test in the aspect of numerical calculations of the technological processes
Ahmadi Brooghani et al. Loading path optimization of a hydroformed part using multilevel response surface method
CN106777457B (zh) 一种固体发动机药柱结构可靠性评估软件系统
KR101817343B1 (ko) 수정된 h-m 모델에 기반한 파괴곡선 산출 방법
CN111339703A (zh) 大应变条件下材料硬化行为的虚拟预测方法
CN114062132A (zh) 一种岩石单轴压缩加速蠕变起始时间的预测方法
CN107908917A (zh) 一种高强板冲压成形回弹预测方法
CN108562495B (zh) 一种完好钢制管道极限内压的计算方法
Gao et al. Instrumented Taylor impact test for measuring stress-strain curve through single trial
CN111931309A (zh) 一种锚杆/锚索结构拉伸硬化软化断裂模拟的方法
JP7103514B2 (ja) 鋼管圧潰強度予測モデルの生成方法、鋼管の圧潰強度予測方法、鋼管の製造特性決定方法、及び鋼管の製造方法
CN114065499B (zh) 岩石单轴蠕变全过程改进Maxwell模型及构建方法
CN112651153A (zh) 一种确定晶体塑性有限元模型材料参数的方法
CN116629048A (zh) 铝合金高温各向异性本构模型和断裂模型的建立方法
Joun et al. Tensile test based material identification program afdex/mat and its application to two new pre-heat treated steels and a conventional cr-mo steel
Shahbazi Karami et al. Experimental and numerical assessment of mechanical properties of thin-walled aluminum parts produced by liquid impact forming
CN115630544A (zh) 用于仿真分析的材料曲线处理方法
Jafari et al. Comparison of ductile fracture models on load bearing capacity of a dented aluminum pipe subjected to internal pressure
Mohammed et al. Evaluating damage scale model of concrete materials using test data
JP7292755B2 (ja) 弾塑性解析方法及び弾塑性解析プログラム
Jin et al. An incremental die forging process for producing helical tubes
Naeimi et al. Optimum designing of forging preform die for the H-shaped parts using backward deformation method and neural networks algorithm
JP2003279455A (ja) 衝撃解析用物性の取得方法および装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant