[go: up one dir, main page]

CN114894616A - Rock creep model viscous and elastic parameter obtaining method based on deformation modulus - Google Patents

Rock creep model viscous and elastic parameter obtaining method based on deformation modulus Download PDF

Info

Publication number
CN114894616A
CN114894616A CN202210399524.7A CN202210399524A CN114894616A CN 114894616 A CN114894616 A CN 114894616A CN 202210399524 A CN202210399524 A CN 202210399524A CN 114894616 A CN114894616 A CN 114894616A
Authority
CN
China
Prior art keywords
creep
rock
deformation modulus
confining pressure
confining
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
CN202210399524.7A
Other languages
Chinese (zh)
Other versions
CN114894616B (en
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.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science 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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN202210399524.7A priority Critical patent/CN114894616B/en
Publication of CN114894616A publication Critical patent/CN114894616A/en
Priority to PCT/CN2023/082097 priority patent/WO2023197824A1/en
Application granted granted Critical
Publication of CN114894616B publication Critical patent/CN114894616B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02Details
    • G01N3/06Special adaptations of indicating or recording means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001Type of application of the stress
    • G01N2203/0003Steady
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0071Creep
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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

The invention discloses a deformation modulus-based method for acquiring viscous and elastic parameters of a rock creep model, which comprises the following steps of: carrying out triaxial compression tests under different confining pressures to obtain the strength and deformation modulus of the rock under different confining pressures; carrying out creep tests on rocks under different confining pressures by adopting a graded loading mode; and obtaining the deformation modulus and the long-term strength under the confining pressure value according to the given confining pressure value and the established confining pressure-deformation modulus relation and confining pressure-long-term strength relation, and obtaining the viscosity and elasticity parameters and the creep-time curve of the model by combining the established relationship between the viscosity and elasticity parameters and the deformation modulus of the burgers creep model. According to the creep model parameter-deformation modulus relation method, only a few rock creep test data under confining pressure are needed to establish the creep model parameter-deformation modulus relation, and the creep model parameters under confining pressure and long-term strength and the test curve reflecting the whole process of accelerated creep of the material can be obtained only through the deformation modulus of the rock under certain confining pressure in the later period, so that the method is provided for the creep rupture behavior research of the material.

Description

一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法A method for obtaining viscoelastic and elastic parameters of rock creep model based on deformation modulus

技术领域technical field

本发明涉及一种能够通过岩石常规变形模量即可直接计算得到岩石蠕变模型参数的方法,属于岩石工程技术领域。The invention relates to a method for obtaining rock creep model parameters by directly calculating the conventional rock deformation modulus, and belongs to the technical field of rock engineering.

背景技术Background technique

作为岩石工程设计及长期稳定性评价的基础,岩石蠕变模型一直以来都是岩石工程领域科研人员关注的焦点。截止目前,人们针对不同种类、不同赋存状态及不同应力路径下的岩石蠕变模型做了诸多的开发,现存的各种蠕变本构模型基本能够满足工程需要,但如何确定蠕变模型粘、弹性参数仍是关键。以往岩石的蠕变模型参数获取,通常首先通过室内试验获得蠕变与时间的关系,然后将试验数据通过本构关系拟合得到蠕变模型参数。这种方法虽然精确度高,但费时费力。As the basis of rock engineering design and long-term stability evaluation, rock creep model has always been the focus of researchers in the field of rock engineering. Up to now, people have done a lot of development for rock creep models of different types, different occurrence states and different stress paths. The existing creep constitutive models can basically meet the needs of engineering, but how to determine the viscosity of the creep model? , the elasticity parameter is still the key. In the past, to obtain the parameters of the creep model of the rock, the relationship between creep and time is usually obtained through laboratory tests, and then the parameters of the creep model are obtained by fitting the test data through the constitutive relationship. Although this method is highly accurate, it is time-consuming and labor-intensive.

发明内容SUMMARY OF THE INVENTION

针对现有蠕变模型参数获取方法的不足,本发明提供一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法。Aiming at the deficiencies of the existing creep model parameter acquisition methods, the present invention provides a deformation modulus-based method for acquiring viscoelastic and elastic parameters of rock creep model.

本发明所采用的技术解决方案是:The technical solution adopted by the present invention is:

一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,包括以下步骤:A method for obtaining viscoelastic and elastic parameters of rock creep model based on deformation modulus, comprising the following steps:

(1)开展不同围压下的三轴压缩试验,获得不同围压下岩石的强度和变形模量;(1) Carry out triaxial compression tests under different confining pressures to obtain the strength and deformation modulus of rocks under different confining pressures;

(2)以围压为横坐标、变形模量为纵坐标,将不同围压下岩石的变形模量值绘制在坐标系中,并通过函数拟合两者之间的关系;(2) Take the confining pressure as the abscissa and the deformation modulus as the ordinate, draw the deformation modulus values of the rock under different confining pressures in the coordinate system, and fit the relationship between the two through a function;

(3)采用分级加载方式,对不同围压下的岩石开展蠕变试验;(3) The creep test is carried out on rocks under different confining pressures by adopting the hierarchical loading method;

(4)根据不同围压下岩石分级加载蠕变变形数据,获得岩石蠕变与时间关系;(4) According to the creep deformation data of rock under different confining pressures under different confining pressures, the relationship between rock creep and time is obtained;

(5)根据不同轴向加载水平下的蠕变曲线,绘制岩石等时应力-应变曲线,根据岩石等时应力-应变曲线确定岩石长期强度;(5) According to the creep curve under different axial loading levels, draw the isochronous stress-strain curve of the rock, and determine the long-term strength of the rock according to the isochronous stress-strain curve of the rock;

(6)以岩石围压为横坐标、长期强度为纵坐标,绘制围压-长期强度曲线,通过函数拟合得到围压-长期强度之间的关系;(6) Taking the confining pressure of the rock as the abscissa and the long-term strength as the ordinate, draw the confining pressure-long-term strength curve, and obtain the relationship between the confining pressure and the long-term strength through function fitting;

(7)根据步骤(4)获得的岩石蠕变与时间关系,进一步得到加速蠕变曲线;采用损伤burgers蠕变模型对加速蠕变曲线进行拟合,得到损伤burgers蠕变模型的粘弹性参数

Figure BDA0003599128540000011
Figure BDA0003599128540000012
(7) According to the relationship between rock creep and time obtained in step (4), the accelerated creep curve is further obtained; the accelerated creep curve is fitted by the damage burgers creep model, and the viscoelastic parameters of the damaged burgers creep model are obtained.
Figure BDA0003599128540000011
Figure BDA0003599128540000012

(8)通过步骤(2)得到给定围压下的变形模量,再结合步骤(7)中得到的损伤burgers蠕变模型的粘弹性参数,绘制蠕变模型参数-变形模量曲线,变形模量为横坐标、蠕变模型参数为纵坐标,通过函数拟合两者之间的关系;(8) Obtain the deformation modulus under a given confining pressure through step (2), and then combine the viscoelastic parameters of the damaged burgers creep model obtained in step (7) to draw the creep model parameter-deformation modulus curve, and the deformation The modulus is the abscissa and the creep model parameter is the ordinate, and the relationship between the two is fitted by a function;

(9)根据给定的某一围压值和步骤(2)建立的围压-变形模量关系、步骤(6)建立的围压-长期强度关系,得到该围压值下的变形模量和长期强度,再根据步骤(7)建立的burgers蠕变模型的粘、弹性参数与变形模量的关系,最终得到该围压下burgers蠕变模型的粘、弹性参数以及蠕变-时间曲线。(9) According to a given confining pressure value, the confining pressure-deformation modulus relationship established in step (2), and the confining pressure-long-term strength relationship established in step (6), the deformation modulus under the confining pressure value is obtained and long-term strength, and then according to the relationship between the viscous and elastic parameters and the deformation modulus of the Burgers creep model established in step (7), the viscosity and elastic parameters and the creep-time curve of the Burgers creep model under the confining pressure are finally obtained.

优选的,步骤(3)中:分级加载从对应围压下抗压强度的60%开始,其后逐级加压,每次加压幅度5%~10%,直至出现加速蠕变破坏为止。Preferably, in step (3): the step-by-step loading starts from 60% of the compressive strength under the corresponding confining pressure, and then pressurizes step by step, each time by 5% to 10%, until accelerated creep failure occurs.

优选的,步骤(7)中损伤burgers蠕变模型的表达式如下:Preferably, the expression of the damaged burgers creep model in step (7) is as follows:

Figure BDA0003599128540000021
Figure BDA0003599128540000021

式中:ε为应变,σ为应力,m为反映岩石蠕变损伤的程度的参数,t为时间,εc为粘性应变。where ε is the strain, σ is the stress, m is the parameter reflecting the degree of rock creep damage, t is the time, and εc is the viscous strain.

本发明的有益技术效果是:The beneficial technical effects of the present invention are:

本发明提出的蠕变模型粘、弹性参数获取方法只需要少数围压下岩石蠕变试验数据便可建立围压-变形模量、围压-长期强度以及蠕变模型参数-变形模量关系,后期只通过某围压下岩石的变形模量即可得到该围压与长期强度下的蠕变模型参数和反映材料加速蠕变全部过程试验曲线,为材料蠕变破坏行为研究提供方法。本发明方法可为岩石工程设计与评价提供准确、便捷的本构方程。The method for obtaining the viscous and elastic parameters of the creep model proposed by the invention only needs a few rock creep test data under confining pressure to establish the relationship between confining pressure-deformation modulus, confining pressure-long-term strength and creep model parameter-deformation modulus. In the later stage, only the deformation modulus of the rock under a certain confining pressure can be used to obtain the creep model parameters under the confining pressure and long-term strength and the test curve reflecting the whole process of accelerated creep of the material, which provides a method for the study of material creep failure behavior. The method of the invention can provide accurate and convenient constitutive equation for rock engineering design and evaluation.

附图说明Description of drawings

图1为本发明的实施流程图,图中给出了蠕变模型参数及基于变形模量的蠕变曲线获取流程;Fig. 1 is the implementation flow chart of the present invention, and the creep model parameter and the creep curve acquisition process based on deformation modulus are provided in the figure;

图2为中粒砂岩在围压1.5MPa、3.0MPa、4.5MPa和6.0MPa下的应力应变曲线;Figure 2 shows the stress-strain curves of medium-grained sandstone under confining pressures of 1.5MPa, 3.0MPa, 4.5MPa and 6.0MPa;

图3为中粒砂岩变形模量-围压坐标点及拟合曲线;Fig. 3 is the deformation modulus-confining pressure coordinate point and fitting curve of medium-grained sandstone;

图4为中粒砂岩在围压1.5MPa下6种应力水平的蠕变曲线;Figure 4 shows the creep curves of 6 stress levels of medium-grained sandstone under a confining pressure of 1.5MPa;

图5为中粒砂岩在围压3.0MPa下8种应力水平的蠕变曲线;Figure 5 shows the creep curves of 8 stress levels of medium-grained sandstone under confining pressure of 3.0MPa;

图6为中粒砂岩在围压4.5MPa下7种应力水平的蠕变曲线;Figure 6 shows the creep curves of 7 stress levels of medium-grained sandstone under confining pressure of 4.5MPa;

图7为中粒砂岩在围压6.0MPa下9种应力水平的蠕变曲线;Figure 7 shows the creep curves of 9 stress levels of medium-grained sandstone under confining pressure of 6.0 MPa;

图8为中粒砂岩长期强度-围压坐标点及拟合曲线;Fig. 8 is the long-term strength-confining pressure coordinate point and fitting curve of medium-grained sandstone;

图9为中粒砂岩蠕变模型参数-变形模量坐标点及拟合曲线,图中

Figure BDA0003599128540000022
Figure BDA0003599128540000023
的两条线近乎重合;Figure 9 shows the parameters of the medium-grained sandstone creep model-deformation modulus coordinates and fitting curves.
Figure BDA0003599128540000022
and
Figure BDA0003599128540000023
The two lines are almost coincident;

图10为按照本发明方法通过变形模量获得的蠕变模型参数与实验数据的对比图。FIG. 10 is a graph comparing the parameters of the creep model obtained by the deformation modulus according to the method of the present invention and the experimental data.

具体实施方式Detailed ways

如图1所示,一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,包括以下步骤:As shown in Figure 1, a method for obtaining viscoelastic and elastic parameters of rock creep model based on deformation modulus includes the following steps:

(1)开展不同围压下的常规三轴压缩试验,获得不同围压下岩石的抗压强度和变形模量。(1) Carry out conventional triaxial compression tests under different confining pressures to obtain the compressive strength and deformation modulus of rocks under different confining pressures.

(2)以围压为横坐标、变形模量为纵坐标,将不同围压下岩石的变形模量值绘制在坐标系中,并通过合适的函数拟合两者之间的关系,即得到围压-变形模量关系。(2) Take the confining pressure as the abscissa and the deformation modulus as the ordinate, draw the deformation modulus values of the rock under different confining pressures in the coordinate system, and fit the relationship between the two through a suitable function, that is, get Confining pressure-deformation modulus relationship.

(3)采用分级加载方式,对不同围压下的岩石开展蠕变试验,即开展不同围压下三轴分级加载蠕变试验。分级加载从对应围压下抗压强度的60%左右开始,其后逐级加压,每次加压幅度5%~10%,直至出现加速蠕变破坏为止。(3) The creep test is carried out on rocks under different confining pressures by adopting the hierarchical loading method, that is, the triaxial hierarchical loading creep test under different confining pressures is carried out. The graded loading starts from about 60% of the compressive strength under the corresponding confining pressure, and then is pressurized step by step, with a pressurization range of 5% to 10% each time, until accelerated creep failure occurs.

(4)根据不同围压下岩石分级加载蠕变变形数据,绘制岩石蠕变-时间曲线,获得岩石蠕变与时间关系。(4) According to the graded loading creep deformation data of rock under different confining pressures, draw the rock creep-time curve to obtain the relationship between rock creep and time.

(5)根据不同轴向加载水平下的蠕变曲线,绘制岩石等时应力-应变曲线,根据岩石等时应力-应变曲线确定岩石长期强度。(5) According to the creep curve under different axial loading levels, draw the isochronous stress-strain curve of the rock, and determine the long-term strength of the rock according to the isochronous stress-strain curve of the rock.

(6)在步骤(5)的基础上,以岩石围压为横坐标、长期强度为纵坐标,绘制围压-长期强度曲线,通过函数拟合得到围压-长期强度之间的关系。(6) On the basis of step (5), take the confining pressure of the rock as the abscissa and the long-term strength as the ordinate, draw the confining pressure-long-term strength curve, and obtain the relationship between the confining pressure and the long-term strength through function fitting.

(7)根据步骤(4)获得的岩石蠕变与时间关系,进一步得到加速蠕变曲线。采用损伤burgers蠕变模型对加速蠕变曲线进行拟合,得到损伤burgers蠕变模型的粘弹性参数

Figure BDA0003599128540000031
Figure BDA0003599128540000032
(7) According to the relationship between rock creep and time obtained in step (4), an accelerated creep curve is further obtained. The accelerated creep curve is fitted by the damage burgers creep model, and the viscoelastic parameters of the damaged burgers creep model are obtained.
Figure BDA0003599128540000031
Figure BDA0003599128540000032

损伤burgers蠕变模型的表达式如下:The expression of the damage burgers creep model is as follows:

Figure BDA0003599128540000033
Figure BDA0003599128540000033

式中:ε为应变,σ为应力,m为反映岩石蠕变损伤的程度的参数,t为时间,εc为粘性应变。where ε is the strain, σ is the stress, m is the parameter reflecting the degree of rock creep damage, t is the time, and εc is the viscous strain.

(8)通过步骤(2)得到给定围压下的变形模量,再结合步骤(7)中得到的损伤burgers蠕变模型的粘弹性参数,绘制蠕变模型参数-变形模量曲线,岩石变形模量为横坐标、蠕变模型参数为纵坐标,通过合适的函数拟合两者之间的关系,从而可进一步得出给定围压下的蠕变模型参数。(8) Obtain the deformation modulus under a given confining pressure through step (2), and then combine the viscoelastic parameters of the damaged burgers creep model obtained in step (7) to draw the creep model parameter-deformation modulus curve. The deformation modulus is the abscissa and the creep model parameter is the ordinate. The relationship between the two can be fitted by a suitable function, so that the creep model parameters under a given confining pressure can be further obtained.

(9)根据给定的某一围压值和步骤(2)建立的围压-变形模量关系、步骤(6)建立的围压-长期强度关系,得到该围压值下的变形模量和长期强度,再根据步骤(7)建立的burgers蠕变模型的粘、弹性参数与变形模量的关系,最终得到该围压下burgers蠕变模型的粘、弹性参数以及蠕变-时间曲线。(9) According to a given confining pressure value, the confining pressure-deformation modulus relationship established in step (2), and the confining pressure-long-term strength relationship established in step (6), the deformation modulus under the confining pressure value is obtained and long-term strength, and then according to the relationship between the viscous and elastic parameters and the deformation modulus of the Burgers creep model established in step (7), the viscosity and elastic parameters and the creep-time curve of the Burgers creep model under the confining pressure are finally obtained.

本发明提出的蠕变模型粘、弹性参数获取方法只需要少数围压下岩石蠕变试验数据便可建立围压-变形模量、围压-长期强度以及蠕变模型参数-变形模量关系,后期只通过某围压下岩石的变形模量即可得到该围压与长期强度下的蠕变模型参数和反映材料加速蠕变全部过程试验曲线,为材料蠕变破坏行为研究提供方法。The method for obtaining the viscous and elastic parameters of the creep model proposed by the invention only needs a few rock creep test data under confining pressure to establish the relationship between confining pressure-deformation modulus, confining pressure-long-term strength and creep model parameter-deformation modulus. In the later stage, only the deformation modulus of the rock under a certain confining pressure can be used to obtain the creep model parameters under the confining pressure and long-term strength and the test curve reflecting the whole process of accelerated creep of the material, which provides a method for the study of material creep failure behavior.

下面通过具体应用实例对本发明作进一步说明:The present invention is further described below by specific application example:

本发明一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,其实施流程如图1所示,其在中粒砂岩的蠕变模型粘、弹性参数及某给定围压下岩石加速蠕变曲线获得中的运用包括以下步骤:The present invention is a method for obtaining viscoelastic and elastic parameters of rock creep model based on deformation modulus, and its implementation process is shown in Figure 1. The use of accelerated creep curve acquisition includes the following steps:

(1)开展不同围压下(1.5MPa、3.0MPa、4.5MPa和6.0MPa)中粒砂岩的常规三轴压缩试验,获得中粒砂岩应力应变曲线(如图2所示)以及四种围压下中粒砂岩的强度和变形模量。(1) Carry out conventional triaxial compression test of medium-grained sandstone under different confining pressures (1.5MPa, 3.0MPa, 4.5MPa and 6.0MPa), and obtain the stress-strain curve of medium-grained sandstone (as shown in Figure 2) and four kinds of confining pressures Strength and deformation modulus of lower-medium-grained sandstone.

(2)以中粒砂岩围压为横坐标、变形模量为纵坐标,绘制中粒砂岩变形模量-围压坐标点,拟合等得到两者之间的关系式,如图3所示。(2) Taking the confining pressure of the medium-grained sandstone as the abscissa and the deformation modulus as the ordinate, draw the coordinate points of the deformation modulus and the confining pressure of the medium-grained sandstone, and obtain the relationship between the two by fitting, as shown in Figure 3 .

(3)采用分级加载方式,对四种围压下的中粒砂岩开展蠕变试验,加载从对应围压下抗压强度的60%左右开始,行第一级加载,此后按2.0MPa的应力梯度进行分级加载,每级加载6个小时,直至出现加速蠕变破坏为止。(3) The graded loading method was used to carry out the creep test on the medium-grained sandstone under four confining pressures. The loading started from about 60% of the compressive strength under the corresponding confining pressure, and the first-level loading was performed, and then the stress was 2.0MPa. Gradient loading is carried out in stages, each stage is loaded for 6 hours, until accelerated creep failure occurs.

(4)根据四种围压下中粒砂岩分级加载蠕变试验数据,绘制得到四种围压下中粒砂岩分级加载蠕变-时间曲线,如图4~7所示。(4) According to the graded loading creep test data of medium-grained sandstone under four confining pressures, the graded loading creep-time curves of medium-grained sandstone under four confining pressures are drawn, as shown in Figures 4-7.

(5)根据不同轴向加载水平下的蠕变曲线,绘制得到中粒砂岩等时应力-应变曲线,根据求取长期强度的等时应力-应变法得到四种围压下中粒砂岩长期强度分别为19.5MPa、26.5MPa、29.0MPa、34MPa。(5) According to the creep curves under different axial loading levels, the isochronous stress-strain curves of medium-grained sandstone are drawn, and the long-term strengths of medium-grained sandstones under four confining pressures are obtained according to the isochronous stress-strain method for obtaining long-term strength. They are 19.5MPa, 26.5MPa, 29.0MPa and 34MPa respectively.

(6)以围压为横坐标、长期强度为纵坐标,绘制中粒砂岩围压-长期强度坐标点,拟合得到两者之间的关系,如图8所示。(6) Taking the confining pressure as the abscissa and the long-term strength as the ordinate, draw the coordinate points of the confining pressure and the long-term strength of the medium-grained sandstone, and get the relationship between the two by fitting, as shown in Figure 8.

(7)采用损伤burgers蠕变模型(如下式)对加速蠕变曲线进行拟合,得到损伤burgers蠕变模型的粘弹性参数

Figure BDA0003599128540000041
(7) Use the damage burgers creep model (the following formula) to fit the accelerated creep curve to obtain the viscoelastic parameters of the damaged burgers creep model
Figure BDA0003599128540000041

Figure BDA0003599128540000042
Figure BDA0003599128540000042

式中:ε为应变,σ为应力,m为反映岩石蠕变损伤的程度的参数,t为时间,εc为粘性应变。where ε is the strain, σ is the stress, m is the parameter reflecting the degree of rock creep damage, t is the time, and εc is the viscous strain.

(8)以中粒砂岩变形模量为横坐标、蠕变模型参数

Figure BDA0003599128540000043
为纵坐标,绘制蠕变模型参数-变形模量数据坐标点,通过合适的函数拟合两者之间的关系,如图9所示。(8) Taking the deformation modulus of medium-grained sandstone as the abscissa, the creep model parameters
Figure BDA0003599128540000043
As the ordinate, draw the creep model parameter-deformation modulus data coordinate points, and fit the relationship between the two through a suitable function, as shown in Figure 9.

(9)给定某一围压值4MPa,根据围压-变形模量、围压-长期强度拟合关系,得到围压为4MPa时中粒砂岩变形模量为5.89GPa、长期强度为28.02MPa,再根据得到的burgers模型粘弹性参数与变形模量的关系,得到该围压下burgers模型的粘性和弹性参数

Figure BDA0003599128540000051
Figure BDA0003599128540000052
分别为14.26GPa、13.38GPa*h、129.83GPa、172.47GPa*h。(9) Given a certain confining pressure value of 4MPa, according to the fitting relationship between confining pressure-deformation modulus and confining pressure-long-term strength, it is obtained that the deformation modulus of medium-grained sandstone is 5.89GPa and the long-term strength is 28.02MPa when the confining pressure is 4MPa. , and then according to the obtained relationship between the viscoelastic parameters of the Burgers model and the deformation modulus, the viscous and elastic parameters of the Burgers model under the confining pressure are obtained
Figure BDA0003599128540000051
Figure BDA0003599128540000052
They are 14.26GPa, 13.38GPa*h, 129.83GPa, and 172.47GPa*h, respectively.

(10)通过上述获得的围压4MPa下产生加速蠕变时的模型参数

Figure BDA0003599128540000053
可得加载应力为28.02MPa时的蠕变-时间曲线,如图10所示,通过与实验数据对比,可知通过变形模量获得的蠕变模型参数可较好的描述岩石的加速蠕变行为。(10) Model parameters when accelerated creep occurs under the confining pressure 4MPa obtained above
Figure BDA0003599128540000053
The creep-time curve when the loading stress is 28.02MPa can be obtained, as shown in Figure 10. By comparing with the experimental data, it can be seen that the creep model parameters obtained by the deformation modulus can better describe the accelerated creep behavior of rock.

上述方式中未述及的部分采取或借鉴已有技术即可实现。The parts not mentioned in the above manner can be realized by adopting or learning from the existing technology.

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above are only the preferred embodiments of the present invention. It should be pointed out that for those skilled in the art, without departing from the principles of the present invention, several improvements and modifications can be made. It should be regarded as the protection scope of the present invention.

Claims (3)

1.一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,其特征在于包括以下步骤:1. a method for obtaining rock creep model viscosity and elasticity parameters based on deformation modulus, is characterized in that comprising the following steps: (1)开展不同围压下的三轴压缩试验,获得不同围压下岩石的强度和变形模量;(1) Carry out triaxial compression tests under different confining pressures to obtain the strength and deformation modulus of rocks under different confining pressures; (2)以围压为横坐标、变形模量为纵坐标,将不同围压下岩石的变形模量值绘制在坐标系中,并通过函数拟合两者之间的关系;(2) Take the confining pressure as the abscissa and the deformation modulus as the ordinate, draw the deformation modulus values of the rock under different confining pressures in the coordinate system, and fit the relationship between the two through a function; (3)采用分级加载方式,对不同围压下的岩石开展蠕变试验;(3) The creep test is carried out on rocks under different confining pressures by adopting the hierarchical loading method; (4)根据不同围压下岩石分级加载蠕变变形数据,获得岩石蠕变与时间关系;(4) According to the creep deformation data of rock under different confining pressures under different confining pressures, the relationship between rock creep and time is obtained; (5)根据不同轴向加载水平下的蠕变曲线,绘制岩石等时应力-应变曲线,根据岩石等时应力-应变曲线确定岩石长期强度;(5) According to the creep curve under different axial loading levels, draw the isochronous stress-strain curve of the rock, and determine the long-term strength of the rock according to the isochronous stress-strain curve of the rock; (6)以岩石围压为横坐标、长期强度为纵坐标,绘制围压-长期强度曲线,通过函数拟合得到围压-长期强度之间的关系;(6) Taking the confining pressure of the rock as the abscissa and the long-term strength as the ordinate, draw the confining pressure-long-term strength curve, and obtain the relationship between the confining pressure and the long-term strength through function fitting; (7)根据步骤(4)获得的岩石蠕变与时间关系,进一步得到加速蠕变曲线;采用损伤burgers蠕变模型对加速蠕变曲线进行拟合,得到损伤burgers蠕变模型的粘弹性参数
Figure FDA0003599128530000012
Figure FDA0003599128530000013
(7) According to the relationship between rock creep and time obtained in step (4), the accelerated creep curve is further obtained; the accelerated creep curve is fitted by the damage burgers creep model, and the viscoelastic parameters of the damaged burgers creep model are obtained.
Figure FDA0003599128530000012
Figure FDA0003599128530000013
(8)通过步骤(2)得到给定围压下的变形模量,再结合步骤(7)中得到的损伤burgers蠕变模型的粘弹性参数,绘制蠕变模型参数-变形模量曲线,变形模量为横坐标、蠕变模型参数为纵坐标,通过函数拟合两者之间的关系;(8) Obtain the deformation modulus under a given confining pressure through step (2), and then combine the viscoelastic parameters of the damaged burgers creep model obtained in step (7) to draw the creep model parameter-deformation modulus curve, and the deformation The modulus is the abscissa and the creep model parameter is the ordinate, and the relationship between the two is fitted by a function; (9)根据给定的某一围压值和步骤(2)建立的围压-变形模量关系、步骤(6)建立的围压-长期强度关系,得到该围压值下的变形模量和长期强度,再根据步骤(7)建立的burgers蠕变模型的粘、弹性参数与变形模量的关系,最终得到该围压下burgers蠕变模型的粘、弹性参数以及蠕变-时间曲线。(9) According to a given confining pressure value, the confining pressure-deformation modulus relationship established in step (2), and the confining pressure-long-term strength relationship established in step (6), the deformation modulus under the confining pressure value is obtained and long-term strength, and then according to the relationship between the viscous and elastic parameters and the deformation modulus of the Burgers creep model established in step (7), the viscosity and elastic parameters and the creep-time curve of the Burgers creep model under the confining pressure are finally obtained.
2.根据权利要求1所述的一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,其特征在于,步骤(3)中:分级加载从对应围压下抗压强度的60%开始,其后逐级加压,每次加压幅度5%~10%,直至出现加速蠕变破坏为止。2. A method for obtaining viscoelastic and elastic parameters of a rock creep model based on deformation modulus according to claim 1, wherein in step (3): grading loading from 60% of the compressive strength under corresponding confining pressure Begin, and then pressurize step by step, with a pressurization range of 5% to 10% each time, until accelerated creep failure occurs. 3.根据权利要求1所述的一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法,其特征在于,步骤(7)中损伤burgers蠕变模型的表达式如下:3. a kind of rock creep model viscous and elastic parameter acquisition method based on deformation modulus according to claim 1, is characterized in that, in step (7), the expression of damage burgers creep model is as follows:
Figure FDA0003599128530000011
Figure FDA0003599128530000011
式中:ε为应变,σ为应力,m为反映岩石蠕变损伤程度的参数,t为时间,εc为粘性应变。where ε is the strain, σ is the stress, m is the parameter reflecting the degree of rock creep damage, t is the time, and εc is the viscous strain.
CN202210399524.7A 2022-04-15 2022-04-15 A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus Active CN114894616B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210399524.7A CN114894616B (en) 2022-04-15 2022-04-15 A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus
PCT/CN2023/082097 WO2023197824A1 (en) 2022-04-15 2023-03-17 Rock creep model viscoelasticity parameter acquisition method based on deformation modulus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210399524.7A CN114894616B (en) 2022-04-15 2022-04-15 A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus

Publications (2)

Publication Number Publication Date
CN114894616A true CN114894616A (en) 2022-08-12
CN114894616B CN114894616B (en) 2023-06-06

Family

ID=82716683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210399524.7A Active CN114894616B (en) 2022-04-15 2022-04-15 A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus

Country Status (2)

Country Link
CN (1) CN114894616B (en)
WO (1) WO2023197824A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023197824A1 (en) * 2022-04-15 2023-10-19 安徽理工大学 Rock creep model viscoelasticity parameter acquisition method based on deformation modulus
CN116929927A (en) * 2023-07-19 2023-10-24 青岛理工大学 A prediction method, system, equipment and medium for creep crack initiation and expansion
CN117909624A (en) * 2024-01-05 2024-04-19 西安建筑科技大学 A method for establishing rock damage creep model

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117538154B (en) * 2023-11-15 2024-04-30 东北石油大学 3D printing proportioning cement stone ultra-early short-term creep test method based on pressing method
CN118861618A (en) * 2024-06-14 2024-10-29 中煤第三建设(集团)有限责任公司 A deep surrounding rock stability evaluation support method and device

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008209262A (en) * 2007-02-27 2008-09-11 Akita Univ Rapid evaluation method for elastic, plastic and creep properties
JP2012202908A (en) * 2011-03-28 2012-10-22 Hiroyuki Sato Prediction method of creep curve and creep lifetime
CN104598669A (en) * 2014-12-22 2015-05-06 重庆交通大学 Method for forecasting permanent deformation of bituminous mixture pavement
CN104915574A (en) * 2015-07-03 2015-09-16 长江水利委员会长江科学院 Establishing method for soft rock creep constitutive model suitable for whole loading and unloading process
CN106885733A (en) * 2017-03-15 2017-06-23 中国电建集团华东勘测设计研究院有限公司 A kind of method that rock mass long-term strength is directly determined by fast creep experiment
CN110174307A (en) * 2019-04-25 2019-08-27 华北水利水电大学 Parameter identification method based on transverse isotropic rockmass three-dimensional creep model
CN110610062A (en) * 2019-10-21 2019-12-24 安徽理工大学 A Calculation Method of Rock Elastic Modulus Under Long-term Static Load
CN110631908A (en) * 2019-09-17 2019-12-31 西安建筑科技大学 A Construction Method of Creep Damage Model in the Whole Process of Rock Uniaxial Compression
CN113053473A (en) * 2021-04-09 2021-06-29 中国工程物理研究院研究生院 Method for constructing full-stage creep model of polymer bonded composite material
CN113919148A (en) * 2021-09-29 2022-01-11 中南大学 A method for building nonlinear creep model of rock

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1479846A1 (en) * 1987-04-06 1989-05-15 Всесоюзный научно-исследовательский институт горной геомеханики и маркшейдерского дела Method for determining indexes of long-term strength of rocks
CN108152137B (en) * 2017-12-28 2020-02-14 辽宁工程技术大学 Method for determining long-term strength in rock creep test
CN109060549B (en) * 2018-10-15 2019-07-16 西安建筑科技大学 A method for establishing the creep model of rock in the whole process of uniaxial compression based on sigmoid function
CN114894616B (en) * 2022-04-15 2023-06-06 安徽理工大学 A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008209262A (en) * 2007-02-27 2008-09-11 Akita Univ Rapid evaluation method for elastic, plastic and creep properties
JP2012202908A (en) * 2011-03-28 2012-10-22 Hiroyuki Sato Prediction method of creep curve and creep lifetime
CN104598669A (en) * 2014-12-22 2015-05-06 重庆交通大学 Method for forecasting permanent deformation of bituminous mixture pavement
CN104915574A (en) * 2015-07-03 2015-09-16 长江水利委员会长江科学院 Establishing method for soft rock creep constitutive model suitable for whole loading and unloading process
CN106885733A (en) * 2017-03-15 2017-06-23 中国电建集团华东勘测设计研究院有限公司 A kind of method that rock mass long-term strength is directly determined by fast creep experiment
CN110174307A (en) * 2019-04-25 2019-08-27 华北水利水电大学 Parameter identification method based on transverse isotropic rockmass three-dimensional creep model
CN110631908A (en) * 2019-09-17 2019-12-31 西安建筑科技大学 A Construction Method of Creep Damage Model in the Whole Process of Rock Uniaxial Compression
CN110610062A (en) * 2019-10-21 2019-12-24 安徽理工大学 A Calculation Method of Rock Elastic Modulus Under Long-term Static Load
CN113053473A (en) * 2021-04-09 2021-06-29 中国工程物理研究院研究生院 Method for constructing full-stage creep model of polymer bonded composite material
CN113919148A (en) * 2021-09-29 2022-01-11 中南大学 A method for building nonlinear creep model of rock

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
丁炎志;张强勇;张龙云;任明洋;尹相杰;王斌;余光远;: "深部地质处置地下实验室花岗岩的加卸荷三轴蠕变试验与分析", 中南大学学报(自然科学版), no. 04 *
刘文博;张树光;: "基于应力和时间双重影响下岩石蠕变模型研究", 中南大学学报(自然科学版), no. 08 *
张树光;刘文博;王有涛;: "花岗岩非定常蠕变及其参数确定", 防灾减灾工程学报, no. 03 *
杨尚阳;张龙云;: "硬脆性岩石卸荷流变特性及本构模型研究", 水利水电技术, no. 10 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023197824A1 (en) * 2022-04-15 2023-10-19 安徽理工大学 Rock creep model viscoelasticity parameter acquisition method based on deformation modulus
CN116929927A (en) * 2023-07-19 2023-10-24 青岛理工大学 A prediction method, system, equipment and medium for creep crack initiation and expansion
CN117909624A (en) * 2024-01-05 2024-04-19 西安建筑科技大学 A method for establishing rock damage creep model

Also Published As

Publication number Publication date
CN114894616B (en) 2023-06-06
WO2023197824A1 (en) 2023-10-19

Similar Documents

Publication Publication Date Title
CN114894616B (en) A Method for Obtaining Viscosity and Elasticity Parameters of Rock Creep Model Based on Deformation Modulus
CN110296892B (en) Method for determining characteristic stress in rock damage evolution process based on energy analysis
CN105181435B (en) A Construction Method of Elastoplastic Constitutive Model of Rock Material
CN113032955B (en) A Construction Method of Rock Dynamic Constitutive Model Applicable to Earthquake Loads
CN110031307B (en) A method for determining the crack initiation stress index of damaged rock
CN108166963B (en) Method for evaluating fracturing effect of offshore oil and gas well
CN103216725B (en) Design method of composite material pressure container
CN108469386A (en) Method is determined based on the material of rock and soil feature of Mohr-Coulomb criterion of strength
CN114662356B (en) Fatigue micro-crack extension prediction method based on EBSD characterization and crystal plasticity
CN103760008A (en) Method for determining fracture closure stress of rock under uniaxial compression condition
CN107907409A (en) A kind of method, equipment and the storage device of definite rock crack initiation stress
CN208206682U (en) A kind of system that rock crack initiation stress is determined based on opposite compression strain
CN105606513B (en) Determine the method for crack rock permeability and porosity to effective stress sensitivity
CN112051141A (en) Construction method of composite material strain rate related compression constitutive model
CN108918396A (en) A method of compact reservoir rock permeability is tested with triaxial stress
CN109268004A (en) A kind of coupling of shale gas reservoir dielectric structure and seam net form state sentence knowledge method
Yu et al. Measuring the damage evolution of granite under different quasi-static load rates through acoustic emission time–frequency characteristics and moment tensor analysis
CN115200990B (en) Rock brittleness evaluation method and system based on constitutive parameters
CN113435017B (en) Optimization method for JH-2 constitutive integrity strength equation parameters
CN106338372A (en) Offshore platform damage positioning method based on residual strain energy and system thereof
RU2830048C2 (en) Method of obtaining viscoelastic characteristics of rock creep model based on deformation modulus
CN116305811A (en) Rapid determination method for macroscopic strength parameters of layered rock mass
CN111678787B (en) Rock long-term strength determination method based on numerical fitting
CN102252949B (en) Method for determining loose core porosity under uniaxial stress
CN114626697A (en) Method and system for evaluating cementing quality of two interfaces of high-temperature well cementation cement sheath

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
GR01 Patent grant