CN113591392A - Rock damage assessment method and device and computer readable storage medium - Google Patents
Rock damage assessment method and device and computer readable storage medium Download PDFInfo
- Publication number
- CN113591392A CN113591392A CN202110910255.1A CN202110910255A CN113591392A CN 113591392 A CN113591392 A CN 113591392A CN 202110910255 A CN202110910255 A CN 202110910255A CN 113591392 A CN113591392 A CN 113591392A
- Authority
- CN
- China
- Prior art keywords
- rock
- damage
- evaluated
- ratio
- under different
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Mathematical Optimization (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Economics (AREA)
- Computational Mathematics (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- General Business, Economics & Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Tourism & Hospitality (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a rock damage assessment method, a rock damage assessment device and a computer readable storage medium.
Description
Technical Field
The invention belongs to the technical field of rock damage assessment, and particularly relates to a rock damage assessment method and device and a computer-readable storage medium.
Background
The damage research of rock materials is always a problem which is paid great attention to in the geotechnical engineering world, and the analysis and the test to rocks at present are generally as follows: the rock is assumed to be homogeneous and homogeneous materials, and the damage degree of the rock is evaluated through macroscopic mechanical parameters such as shear strength, elastic modulus, deformation modulus and the like, and actually, a plurality of fine and complex defects such as cracks and fissures exist in the rock, and the fine defects greatly influence the parameters of the rock, so that the damage analysis precision strength of the rock is not high, and the reliability is poor.
The influence of the defects such as fine cracks on the mechanical properties of the rock can be reflected more accurately by utilizing the damage mechanics theory to research the relation between the strength and the deformation of the rock; therefore, many scholars set up corresponding rock statistical damage constitutive models by defining damage variables with real physical backgrounds from the microscopic structure of the material; however, since the damage variable is used as a parameter in the statistical damage constitutive model of the rock, the damage degree of the rock cannot be evaluated, and therefore, how to accurately evaluate the rock damage becomes an urgent problem to be solved.
Disclosure of Invention
The invention aims to provide a rock damage evaluation method, a rock damage evaluation device and a computer-readable storage medium, which are used for solving the problem that the rock damage degree cannot be evaluated based on damage variables in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a rock damage assessment method, which comprises the following steps:
acquiring experiment data of a rock to be evaluated after a mechanical experiment is carried out, wherein the mechanical experiment comprises a single-axis compression experiment and a three-axis compression experiment;
obtaining distribution parameters in a rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data;
obtaining a critical damage ratio of the rock to be evaluated under different confining pressure conditions by utilizing the rock statistical damage constitutive model, the distribution parameters and the infinitesimal body strength under different confining pressure conditions, wherein the critical damage ratio is a corresponding damage ratio when stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated;
and obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions, so as to obtain the damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
Based on the disclosure, the invention determines the distribution parameters in the rock statistical constitutive model and the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions through the mechanical experiment on the rock; secondly, the damage ratio of the rock to be evaluated is defined based on the damage area and the lossless area of the rock, so that the critical damage ratio of the rock to be evaluated is extended through the damage ratio, and the rock statistical constitutive model, the distribution parameters and the infinitesimal strength are calculated and calculated; finally, the damage index of the rock to be evaluated is obtained through the critical damage ratios under different confining pressure conditions, and is compared with the damage indication threshold value, so that the damage evaluation result of the rock to be evaluated can be obtained according to the comparison result.
Through the design, the damage ratio of the rock is defined, so that the rock critical damage ratio is calculated by using the rock statistical constitutive model, the damage index of the rock is obtained through the critical damage ratio, and the damage index is used as the evaluation index of the rock damage degree.
In one possible design, the experimental data includes an elastic modulus of the rock to be evaluated, and the obtaining of the infinitesimal strength of the rock to be evaluated under different confining pressure conditions in the mechanical experiment by using the experimental data includes:
calculating the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions according to the following formula;
in the formula, F*Alpha and K are D-P failure criterion parameters, I1Is the first invariant of the stress tensor, J2Is the second variable of the stress deflection, E is the elastic modulus, mu is the Poisson's ratio of the rock to be evaluated, sigma1And σ2Respectively the maximum principal stress and the middle principal stress, epsilon, on the regular hexahedral microelements in the rock to be evaluated under different confining pressure conditions1The strain corresponding to the maximum principal stress.
Based on the disclosure, the invention discloses a specific calculation formula of the infinitesimal body strength, namely the infinitesimal body strength under different confining pressure conditions can be obtained by calculating according to a D-P failure criterion and substituting the maximum main stress and the middle main stress of the rock to be evaluated under different confining pressure into the formula.
In one possible design, the experimental data further includes secant modulus, uniaxial compressive strength and strain corresponding to the uniaxial compressive strength of the rock to be evaluated, and the distribution parameters include a first parameter and a second parameter;
the method for obtaining the distribution parameters in the statistical damage constitutive model of the rock corresponding to the rock to be evaluated by using the experimental data comprises the following steps:
obtaining the first parameter according to the elastic modulus, the secant modulus, the uniaxial compressive strength and the strain corresponding to the uniaxial compressive strength;
and obtaining second parameters under different confining pressure conditions according to the first parameters and the infinitesimal body strength under different confining pressure conditions.
In one possible design, obtaining the critical damage ratio of the rock to be evaluated under different confining pressures by using the statistical damage constitutive model of the rock, the distribution parameters and the infinitesimal body strengths under different confining pressure conditions includes:
obtaining a critical damage ratio expression by utilizing the rock statistical damage constitutive model;
and substituting the distribution parameters and the infinitesimal body strength under different confining pressure conditions into the critical damage ratio expression to obtain the critical damage ratio of the rock to be evaluated under different confining pressure conditions.
Based on the above disclosure, the invention obtains a critical damage ratio expression through a rock statistical damage constitutive model, which is substantially as follows: the critical damage ratio is indirectly expressed by using a damage variable, so that the influence of defects such as cracks and fissures in the rock on damage evaluation is considered; and finally, substituting the distribution parameters and the infinitesimal body strength under different confining pressure conditions into an expression to obtain the critical damage ratio under different confining pressure conditions.
In one possible design, obtaining a damage index of the rock to be evaluated by using the critical damage ratio of the rock to be evaluated under different confining pressure conditions includes:
performing linear fitting on the critical damage ratio of the rock to be evaluated under different confining pressure conditions to obtain a linear fitting curve of the critical damage ratio and the confining pressure;
obtaining a linear expression of the critical damage ratio and the confining pressure according to the linear fitting curve;
and calculating the slope of the linear fitting curve according to the linear expression to serve as the damage index.
Based on the above disclosure, the essence of the above steps is: performing linear fitting on the confining pressure and the corresponding critical damage ratio to obtain a linear fitting curve of the critical damage ratio based on the confining pressure; and the slope of the linear fitting curve is used as the damage index.
In one possible design, the damage indication threshold includes a first threshold, a second threshold, and a third threshold, where the first threshold is greater than the second threshold, and the second threshold is greater than the third threshold;
correspondingly, obtaining the damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value, wherein the method comprises the following steps:
if the damage index is larger than the first threshold value, judging that the rock to be evaluated is a seriously damaged rock;
if the damage index is between the third threshold and the second threshold, judging that the rock to be evaluated is medium damage rock;
and if the damage index is smaller than the third threshold value, judging that the rock to be evaluated is slightly damaged rock.
In a second aspect, the present invention provides an apparatus for assessing rock damage, comprising: the system comprises an acquisition unit, a parameter calculation unit, a critical damage ratio generation unit, a damage index generation unit and a damage evaluation unit;
the acquisition unit is used for acquiring experiment data of a rock to be evaluated after a mechanical experiment is carried out on the rock, wherein the mechanical experiment comprises a single-axis compression experiment and a three-axis compression experiment;
the parameter calculation unit is used for obtaining distribution parameters in a rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data;
the critical damage ratio generating unit is configured to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions by using the rock statistical damage constitutive model, the distribution parameters, and the infinitesimal body strengths under different confining pressure conditions, where the critical damage ratio is a corresponding damage ratio when stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated;
the damage index generation unit is used for obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions;
and the damage evaluation unit is used for obtaining a damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
In one possible design, the parameter calculation unit includes: a infinitesimal strength meter operator unit;
the infinitesimal strength meter operator unit is used for calculating and obtaining the infinitesimal strength of the rock to be evaluated under different confining pressure conditions according to the following formula;
in the formula, F*Alpha and K are D-P failure criterion parameters, I1Is the first invariant of the stress tensor, J2Is the second variable of the stress deflection, E is the elastic modulus, mu is the Poisson's ratio of the rock to be evaluated, sigma1And σ2Respectively the maximum principal stress and the middle principal stress, epsilon, on the regular hexahedral microelements in the rock to be evaluated under different confining pressure conditions1The strain corresponding to the maximum principal stress.
In one possible design, the parameter calculation unit further includes: a distribution parameter calculation subunit;
the distribution parameter calculation subunit is configured to obtain the first parameter according to the elastic modulus, the secant modulus, the uniaxial compressive strength, and a strain corresponding to the uniaxial compressive strength;
and the distribution parameter calculating subunit is further configured to obtain a second parameter under different confining pressure conditions according to the first parameter and the infinitesimal body strengths under different confining pressure conditions.
In one possible design:
the critical damage ratio generating unit is specifically configured to obtain a critical damage ratio expression by using the rock statistical damage constitutive model;
the critical damage ratio generating unit is further specifically configured to substitute the distribution parameter and the infinitesimal strength under different confining pressure conditions into the critical damage ratio expression to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions.
In one possible design:
the damage index generation unit is specifically configured to perform linear fitting on the critical damage ratio of the rock to be evaluated under different confining pressure conditions to obtain a linear fitting curve of the critical damage ratio and the confining pressure;
the damage index generation unit is specifically configured to obtain a linear expression of the critical damage ratio and the confining pressure according to the linear fitting curve;
the damage index generating unit is further specifically configured to calculate a slope of the linear fitting curve according to the linear expression as the damage index.
In a third aspect, the present invention provides another rock damage assessment apparatus, which is a computer-based device, and includes a memory, a processor and a transceiver, which are communicatively connected in sequence, where the memory is used to store a computer program, the transceiver is used to transmit and receive messages, and the processor is used to read the computer program and execute the rock damage assessment method as described in the first aspect or any one of the possible designs in the first aspect.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon instructions which, when run on a computer, perform the method of assessing rock damage as described in the first aspect or any one of the possible designs of the first aspect.
In a fifth aspect, the present invention provides a computer program product comprising instructions which, when run on a computer, cause the computer to perform the method of assessing rock damage as set forth in the first aspect or any one of the possible designs of the first aspect.
Drawings
FIG. 1 is a schematic flow chart of the steps of the rock damage assessment method provided by the present invention;
FIG. 2 is an equivalent schematic diagram of a nondestructive area and a damage area in a rock to be evaluated provided by the present invention;
FIG. 3 is a schematic structural diagram of an apparatus for evaluating rock damage according to the present invention;
fig. 4 is a schematic structural diagram of a computer main device provided in the present invention.
Detailed Description
The invention is further described with reference to the following figures and specific embodiments. It should be noted that the description of the embodiments is provided to help understanding of the present invention, but the present invention is not limited thereto. Specific structural and functional details disclosed herein are merely illustrative of example embodiments of the invention. This invention may, however, be embodied in many alternate forms and should not be construed as limited to the embodiments set forth herein.
It will be understood that, although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first element could be termed a second element, and, similarly, a second element could be termed a first element, without departing from the scope of example embodiments of the present invention.
It should be understood that, for the term "and/or" as may appear herein, it is merely an associative relationship that describes an associated object, meaning that three relationships may exist, e.g., a and/or B may mean: a exists alone, B exists alone, and A and B exist at the same time; for the term "/and" as may appear herein, which describes another associative object relationship, it means that two relationships may exist, e.g., a/and B, may mean: a exists independently, and A and B exist independently; in addition, for the character "/" that may appear herein, it generally means that the former and latter associated objects are in an "or" relationship.
Examples
At present, in the damage evaluation of rocks, if the rocks are assumed to be homogeneous isotropic materials and are evaluated through macroscopic mechanical parameters such as shear strength, elastic modulus, deformation modulus and the like, the parameters are greatly influenced by defects such as cracks and cracks in the rocks, so that the damage analysis precision strength is not high, and the reliability is poor; if the damage constitutive model of the rock statistical damage based on the damage variable is adopted for analysis, the problem that the damage degree of the rock cannot be evaluated because the damage variable is only one parameter in the model is faced; therefore, according to the rock damage assessment method provided by this embodiment, the damage ratio of the rock is defined, so that the critical damage ratio of the rock is extended, and the damage ratio is calculated by combining the damage ratio with the rock statistical damage constitutive model (the damage ratio is actually indirectly expressed by the damage variables), so that the rock damage assessment based on the damage ratio is realized by defining the damage ratio of the rock in this embodiment, and not only is the accuracy high, but also the reliability is improved.
The method for evaluating rock damage provided by the first aspect of this embodiment defines a ratio of a damaged area to an undamaged area of rock (i.e., a damage ratio of rock), so as to realize rock damage degree evaluation based on the damage ratio, which may include, but is not limited to, the following steps S101 to S104.
S101, acquiring experimental data of the rock to be evaluated after a mechanical experiment is carried out.
Step S101 is to perform a mechanical experiment on the rock to obtain various experimental data, so as to provide basic data for the subsequent calculation of the model parameters.
In this embodiment, the mechanical experiment may include, but is not limited to: uniaxial compression experiments and triaxial compression experiments; the uniaxial compression experiment of the rock measures the strength, deformation and failure characteristics of the rock subjected to uniaxial compression; obtaining index data such as uniaxial compressive strength, elastic modulus, Poisson's ratio and the like of the rock to be evaluated through uniaxial compression; the triaxial compression test refers to a confined compression and shear test, and can determine data such as an internal friction angle and cohesive force in the shear strength index of the rock to be evaluated.
After acquiring the experimental data of the uniaxial compression experiment and the triaxial compression experiment of the rock to be evaluated, calculating the distribution parameters in the corresponding rock statistical damage constitutive model and the infinitesimal body strength of the rock under different confining pressure conditions by using the experimental data; as shown in step S102 below.
Of course, before calculating the distribution parameters and the strength of the infinitesimal body, a rock statistical damage constitutive model of the rock to be evaluated needs to be constructed, which may include, but is not limited to, the following steps:
a. and acquiring a damage variable of the rock to be evaluated, wherein the damage variable is the ratio of the damage area to the whole area of the rock to be evaluated.
b. And carrying out statistical processing on the rock strength of the rock to be evaluated by utilizing a Weibull distribution function and the damage variable to obtain the statistical damage constitutive model of the rock.
The construction process of the aforementioned model is illustrated below as an example:
suppose that the damage area of the rock to be evaluated is defined as S*The total area of the glass fiber is S when the glass fiber is not damagedmThen the damage variable D*It can be expressed as:
and according to the infinitesimal strength in the rock, the logarithm positive Tai distribution rule is obeyed, so that the rock strength of the rock to be evaluated is subjected to statistical processing by using a Weibull distribution function, and the damage variable D can be obtained*And performing re-expression based on the rock strength, and using the obtained expression as a rock statistical damage constitutive model of the rock to be evaluated, wherein the expression can be, but is not limited to, as follows:
in the above formula, m and F0Respectively, distribution parameters in the model, which reflect the mechanical properties of the brittle material (i.e. the rock to be evaluated), F representing a random distribution variable of the infinitesimal strength, and F*It represents the infinitesimal strength of the rock to be evaluated.
After the statistical damage constitutive model of the rock to be evaluated is constructed, the distribution parameters and the infinitesimal body strength in the model can be calculated by using the experimental data, as shown in the following step S102.
S102, obtaining distribution parameters in the rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data.
In this embodiment, example experimental data may include, but is not limited to: the elastic modulus, cohesive force, internal friction angle, secant modulus, uniaxial compressive strength and strain corresponding to the uniaxial compressive strength of the rock to be evaluated.
In this embodiment, for example, but not limited to, the calculation of the infinitesimal strength may be performed according to Drucker-Prager destruction criteria, and the calculation formula may be, but is not limited to, the following:
in the above formula, F*Alpha and K are D-P failure criterion parameters, I1Is the first invariant of the stress tensor, J2Is the second variable of the stress deflection, E is the elastic modulus, mu is the Poisson's ratio of the rock to be evaluated, sigma1And σ2Respectively the maximum principal stress and the middle principal stress, epsilon, on the regular hexahedral microelements in the rock to be evaluated under different confining pressure conditions1The strain corresponding to the maximum principal stress.
Namely, because the mechanical experiment is carried out under different confining pressure conditions (such as 3MPa, 6MPa, 12MPa and 24MPa), the maximum principal stress, the middle principal stress and the strain corresponding to the maximum principal stress on the regular hexahedral microelements in the rock to be evaluated are different under different confining pressure conditions; therefore, the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions can be obtained by substituting the parameter data such as the maximum principal stress under different confining pressure conditions into the formula.
In this example, the poisson's ratio of the rock to be evaluated is: the ratio of the absolute value of the transverse positive strain to the absolute value of the axial positive strain of the rock to be evaluated under unidirectional tension or compression can be obtained through measurement.
For the same reason, as already explained in the foregoing, the distribution parameter has two, i.e., m and F0Therefore, the first parameter m needs to be calculated first, and then the second parameter F can be calculated through the first parameter m0。
In the present embodiment, the calculation of the distribution parameter may include, but is not limited to, the following step S102a and step S102b.
S102a, obtaining the first parameter according to the elastic modulus, the secant modulus, the uniaxial compressive strength and the strain corresponding to the uniaxial compressive strength.
In this embodiment, the first parameter m is calculated by using the following formula:
in the above formula, E represents an elastic modulus, EsRepresenting secant model, σcDenotes uniaxial compressive strength, andcit represents strain corresponding to uniaxial compressive strength.
After the first parameter is obtained by calculation, the second parameter can be obtained by using the first parameter and the infinitesimal intensity, as shown in the following step S102 b:
s102b, obtaining second parameters under different confining pressure conditions according to the first parameters and the infinitesimal body strength under different confining pressure conditions.
In the present embodiment, the following formula is used to calculate the second parameter F by way of example0。
According to the formula, the infinitesimal body strength under different confining pressure conditions is substituted into the formula, and then the second parameters under different confining pressure conditions can be obtained.
After the distribution parameters of the statistical damage constitutive model of the rock and the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions are obtained, the critical damage ratio of the rock to be evaluated under different confining pressure conditions can be calculated by using the model, so that the damage index of the rock to be evaluated can be calculated according to the critical damage ratio.
In the present embodiment, the critical damage ratio is obtained by extending the damage ratio; in the embodiment, the damage ratio of the rock is defined, so that the damage ratio is combined with a rock statistical damage constitutive model to realize rock damage assessment based on the damage ratio.
As shown in fig. 2, the definition of the damage ratio, i.e. the ratio of the damaged area to the undamaged area of the rock to be evaluated, is first explained; σ in FIG. 20Representing the stress, σ, acting on the rock to be evaluated0' denotes the stress acting on the undamaged face in the rock to be evaluated.
Referring to FIG. 2, it can be seen from the foregoing examples that the damage area of the rock to be evaluated under a certain confining pressure condition (e.g. 3 MPa) is S*(ii) a The total area of the rock to be evaluated is S when the rock is not damagedmThen the damage-free area S of the rock to be evaluatede＝Sm-S*。
Then damage ratio DeCan be expressed as follows:
therefore, the damage ratio can be used as an index of the structural damage characteristic of the rock to reflect the evolution characteristic of the microcracks in the rock under the action of the space stress.
Meanwhile, in the rock mechanics experiment process, when the confining pressure is fixed (for example, under the condition of 3 MPa), the internal damage is gradually accumulated along with the increase of the load, when the load cannot be continuously increased, namely the stress applied to the rock to be evaluated reaches the peak value, the internal damage is also accumulated to a certain value, and at the moment, the corresponding damage ratio is taken as the critical damage ratio.
Therefore, the present embodiment may calculate the critical damage ratio of the rock to be evaluated under different confining pressure conditions by combining the statistical damage constitutive model of the rock, as shown in step S103.
S103, obtaining a critical damage ratio of the rock to be evaluated under different confining pressure conditions by utilizing the rock statistical damage constitutive model, the distribution parameters and the infinitesimal body strength under different confining pressure conditions, wherein the critical damage ratio is a corresponding damage ratio when stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated.
In this embodiment, the example of calculating the critical damage ratio under different ambient pressure conditions in step S103 may include, but is not limited to, the following step S103a and step S103b.
And S103a, obtaining a critical damage ratio expression by using the rock statistical damage constitutive model.
And S103b, substituting the distribution parameters and the infinitesimal body strengths under different confining pressure conditions into the critical damage ratio expression to obtain the critical damage ratio of the rock to be evaluated under different confining pressure conditions.
The essence of step S103a is: and characterizing the critical damage ratio by using the damage variable so as to obtain a critical damage ratio expression based on the damage variable.
Since the foregoing has disclosed the expression of damage ratio, the use of damage variables to express damage ratio may be, but is not limited to, the following:
meanwhile, since it has been stated in the foregoing that when the stress applied to the rock to be evaluated reaches a peak value, the corresponding damage ratio is a critical damage ratio, so that only the infinitesimal strength in the foregoing formula needs to be replaced by: the corresponding infinitesimal body strength when the stress reaches the peak value can obtain the expression of the critical damage ratio, as shown in the following:
in the above-mentioned formula, the compound of formula,
representing the corresponding infinitesimal body strength when the stress on the rock to be evaluated in the application process reaches the peak value; referring to the above formula for calculating the strength of the infinitesimal body, it is only necessary to determine the corresponding sigma when the stress reaches the peak value under a certain confining pressure1、σ2And epsilon1Substituting the obtained force into a infinitesimal body strength calculation formula to obtain the corresponding infinitesimal body strength when the stress reaches a peak value, thereby calculating the critical damage ratio of the rock to be evaluated under a certain confining pressure condition.
Similarly, the calculation of the critical damage ratio of the rock to be evaluated under the other confining pressure conditions is consistent with the foregoing example, and details are not repeated herein.
After the critical damage ratio of the rock to be evaluated under different confining pressure conditions is obtained, the damage index of the rock to be evaluated can be calculated according to the critical damage ratio, as shown in the following step S104.
S104, obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions, so as to obtain the damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
In step S104, the calculated damage index may be, but is not limited to, derived by linear fitting of critical damage ratios under different confining pressure conditions, as shown in steps S104a, S104b, and S104c below.
S104a, performing linear fitting on the critical damage ratio of the rock to be evaluated under different confining pressure conditions to obtain a linear fitting curve of the critical damage ratio and the confining pressure.
And S104b, obtaining a linear expression of the critical damage ratio and the confining pressure according to the linear fitting curve.
S104c, calculating the slope of the linear fitting curve according to the linear expression to serve as the damage index.
Namely, the principle of steps S104a to S104c is: the critical damage ratio is obtained by performing linear fitting on different confining pressure data and the critical damage ratio corresponding to the confining pressure data, and a linear fitting curve of the critical damage ratio based on the confining pressure is obtained, and the slope of the linear fitting curve can be used as the damage index of the rock to be evaluated.
In the present embodiment, the critical damage ratio linear fit curve based on the confining pressure can be characterized by, but is not limited to, the following linear expression:
Der＝Alnσ3+B
wherein A represents the slope and B represents a constant, σ3Representing a variable, namely the confining pressure.
From the above linear expression, the slope a is expressed as:
after the damage index of the rock to be evaluated is calculated, the damage index can be used for judgment, namely, the damage evaluation result of the rock to be evaluated is obtained by comparing the size relation between the damage index and the damage indication threshold value.
In this embodiment, the example damage indication threshold includes a first threshold, a second threshold and a third threshold, where the first threshold is greater than the second threshold, and the second threshold is greater than the third threshold; and the three threshold values are obtained by mechanical experiments.
In this embodiment, the evaluation criterion is:
if the damage index is larger than the first threshold value, judging that the rock to be evaluated is a seriously damaged rock; if the damage index is between the third threshold and the second threshold, judging that the rock to be evaluated is medium damage rock; and if the damage index is smaller than the third threshold value, judging that the rock to be evaluated is slightly damaged rock.
Therefore, by the rock damage assessment method described in detail in the foregoing steps S101 to S104, the damage ratio of the rock is defined, so that the rock critical damage ratio is calculated by using the rock statistical constitutive model, the damage index of the rock is obtained by using the critical damage ratio, and the damage index is used as the assessment index of the rock damage degree.
As shown in fig. 3, a second aspect of the present embodiment provides a hardware device for implementing the method for evaluating rock damage described in the first aspect of the embodiment, including: the device comprises an acquisition unit, a parameter calculation unit, a critical damage ratio generation unit, a damage index generation unit and a damage evaluation unit.
The acquisition unit is used for acquiring experimental data of the rock to be evaluated after mechanical experiment, wherein the mechanical experiment comprises a single-axis compression experiment and a three-axis compression experiment.
And the parameter calculation unit is used for obtaining distribution parameters in the rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data.
The critical damage ratio generating unit is configured to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions by using the rock statistical damage constitutive model, the distribution parameters, and the infinitesimal body strengths under different confining pressure conditions, where the critical damage ratio is a corresponding damage ratio when a stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated.
The damage index generation unit is used for obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions.
And the damage evaluation unit is used for obtaining a damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
In one possible design, the parameter calculation unit includes: and a infinitesimal strength measuring operator unit.
And the infinitesimal strength meter operator unit is used for calculating and obtaining the infinitesimal strength of the rock to be evaluated under different confining pressure conditions according to the following formula.
In the formula, F*Alpha and K are D-P failure criterion parameters, I1Is the first invariant of the stress tensor, J2Is the second variable of the stress deflection, E is the elastic modulus, mu is the Poisson's ratio of the rock to be evaluated, sigma1And σ2Respectively the maximum principal stress and the middle principal stress, epsilon, on the regular hexahedral microelements in the rock to be evaluated under different confining pressure conditions1The strain corresponding to the maximum principal stress.
In one possible design, the parameter calculation unit further includes: and a distribution parameter calculation subunit.
And the distribution parameter calculation subunit is used for obtaining the first parameter according to the elastic modulus, the secant modulus, the uniaxial compressive strength and the strain corresponding to the uniaxial compressive strength.
And the distribution parameter calculating subunit is further configured to obtain a second parameter under different confining pressure conditions according to the first parameter and the infinitesimal body strengths under different confining pressure conditions.
In one possible design:
the critical damage ratio generating unit is specifically configured to obtain a critical damage ratio expression by using the rock statistical damage constitutive model.
The critical damage ratio generating unit is further specifically configured to substitute the distribution parameter and the infinitesimal strength under different confining pressure conditions into the critical damage ratio expression to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions.
In one possible design:
the damage index generation unit is specifically configured to perform linear fitting on the critical damage ratio of the rock to be evaluated under different confining pressure conditions to obtain a linear fitting curve of the critical damage ratio and the confining pressure.
And the damage index generation unit is specifically configured to obtain a linear expression of the critical damage ratio and the confining pressure according to the linear fitting curve.
The damage index generating unit is further specifically configured to calculate a slope of the linear fitting curve according to the linear expression as the damage index.
For the working process, the working details, and the technical effects of the hardware apparatus provided in this embodiment, reference may be made to the first aspect of the embodiment, which is not described herein again.
As shown in fig. 4, a third aspect of the present embodiment provides a computer main apparatus, including: a memory, a processor and a transceiver, which are in communication with each other in that order, wherein the memory is used for storing a computer program, the transceiver is used for transmitting and receiving messages, and the processor is used for reading the computer program and executing the method for assessing rock damage according to the first aspect of the embodiments.
For example, the Memory may include, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Flash Memory (Flash Memory), a First In First Out (FIFO), and/or a First In Last Out (FILO), and the like; the processor may not be limited to a microprocessor of a model number STM32F105 series, a reduced instruction set computer (RSIC) microprocessor, an architecture processor such as X86, or a processor integrated with a neural-Network Processing Unit (NPU); the transceiver may be, but is not limited to, a wireless fidelity (WIFI) wireless transceiver, a bluetooth wireless transceiver, a General Packet Radio Service (GPRS) wireless transceiver, a ZigBee wireless transceiver (ieee802.15.4 standard-based low power local area network protocol), a 3G transceiver, a 4G transceiver, and/or a 5G transceiver, etc. In addition, the device may also include, but is not limited to, a power module, a display screen, and other necessary components.
For the working process, the working details, and the technical effects of the computer main device provided in this embodiment, reference may be made to the first aspect of the embodiment, which is not described herein again.
A fourth aspect of the present embodiment provides a computer-readable storage medium storing instructions including the method for evaluating rock damage according to the first aspect of the present embodiment, that is, the computer-readable storage medium storing instructions that, when executed on a computer, perform the method for evaluating rock damage according to the first aspect.
The computer-readable storage medium refers to a carrier for storing data, and may include, but is not limited to, floppy disks, optical disks, hard disks, flash memories, flash disks and/or Memory sticks (Memory sticks), etc., and the computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices.
For the working process, the working details, and the technical effects of the computer-readable storage medium provided in this embodiment, reference may be made to the first aspect of the embodiment, which is not described herein again.
A fifth aspect of the present embodiments provides a computer program product comprising instructions which, when run on a computer, which may be a general purpose computer, a special purpose computer, a computer network, or other programmable device, cause the computer to perform the method of assessing rock damage according to the first aspect of the embodiments.
Finally, it should be noted that: the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method of assessing rock damage, comprising:
acquiring experiment data of a rock to be evaluated after a mechanical experiment is carried out, wherein the mechanical experiment comprises a single-axis compression experiment and a three-axis compression experiment;
obtaining distribution parameters in a rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data;
obtaining a critical damage ratio of the rock to be evaluated under different confining pressure conditions by utilizing the rock statistical damage constitutive model, the distribution parameters and the infinitesimal body strength under different confining pressure conditions, wherein the critical damage ratio is a corresponding damage ratio when stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated;
and obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions, so as to obtain the damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
2. The evaluation method according to claim 1, wherein the experimental data comprises an elastic modulus of the rock to be evaluated, and wherein the obtaining of the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in the mechanical experiment by using the experimental data comprises:
calculating the infinitesimal body strength of the rock to be evaluated under different confining pressure conditions according to the following formula;
in the formula, F*Alpha and K are D-P failure criterion parameters, I1Is the first invariant of the stress tensor, J2Is the second variable of the stress deflection, E is the elastic modulus, mu is the Poisson's ratio of the rock to be evaluated, sigma1And σ2Respectively the maximum principal stress and the middle principal stress, epsilon, on the regular hexahedral microelements in the rock to be evaluated under different confining pressure conditions1The strain corresponding to the maximum principal stress.
3. The evaluation method according to claim 2, wherein the experimental data further includes secant modulus, uniaxial compressive strength, and strain corresponding to the uniaxial compressive strength of the rock to be evaluated, and the distribution parameters include a first parameter and a second parameter;
the method for obtaining the distribution parameters in the statistical damage constitutive model of the rock corresponding to the rock to be evaluated by using the experimental data comprises the following steps:
obtaining the first parameter according to the elastic modulus, the secant modulus, the uniaxial compressive strength and the strain corresponding to the uniaxial compressive strength;
and obtaining second parameters under different confining pressure conditions according to the first parameters and the infinitesimal body strength under different confining pressure conditions.
4. The evaluation method according to claim 1, wherein the obtaining the critical damage ratio of the rock to be evaluated under different confining pressures by using the statistical damage constitutive model of the rock, the distribution parameters and the infinitesimal body strengths under different confining pressures comprises:
obtaining a critical damage ratio expression by utilizing the rock statistical damage constitutive model;
and substituting the distribution parameters and the infinitesimal body strength under different confining pressure conditions into the critical damage ratio expression to obtain the critical damage ratio of the rock to be evaluated under different confining pressure conditions.
5. The evaluation method according to claim 1, wherein obtaining the damage index of the rock to be evaluated by using the critical damage ratio of the rock to be evaluated under different confining pressure conditions comprises:
performing linear fitting on the critical damage ratio of the rock to be evaluated under different confining pressure conditions to obtain a linear fitting curve of the critical damage ratio and the confining pressure;
obtaining a linear expression of the critical damage ratio and the confining pressure according to the linear fitting curve;
and calculating the slope of the linear fitting curve according to the linear expression to serve as the damage index.
6. The evaluation method of claim 1, wherein the damage indication threshold comprises a first threshold, a second threshold, and a third threshold, wherein the first threshold is greater than the second threshold, and the second threshold is greater than the third threshold;
correspondingly, obtaining the damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value, wherein the method comprises the following steps:
if the damage index is larger than the first threshold value, judging that the rock to be evaluated is a seriously damaged rock;
if the damage index is between the third threshold and the second threshold, judging that the rock to be evaluated is medium damage rock;
and if the damage index is smaller than the third threshold value, judging that the rock to be evaluated is slightly damaged rock.
7. An apparatus for evaluating rock damage, comprising: the system comprises an acquisition unit, a parameter calculation unit, a critical damage ratio generation unit, a damage index generation unit and a damage evaluation unit;
the acquisition unit is used for acquiring experiment data of a rock to be evaluated after a mechanical experiment is carried out on the rock, wherein the mechanical experiment comprises a single-axis compression experiment and a three-axis compression experiment;
the parameter calculation unit is used for obtaining distribution parameters in a rock statistical damage constitutive model corresponding to the rock to be evaluated and infinitesimal body strength of the rock to be evaluated under different confining pressure conditions in a mechanical experiment by using the experimental data;
the critical damage ratio generating unit is configured to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions by using the rock statistical damage constitutive model, the distribution parameters, and the infinitesimal body strengths under different confining pressure conditions, where the critical damage ratio is a corresponding damage ratio when stress applied to the rock to be evaluated reaches a peak value, and the damage ratio is a ratio of a damaged area to a non-damaged area in the rock to be evaluated;
the damage index generation unit is used for obtaining the damage index of the rock to be evaluated by utilizing the critical damage ratio of the rock to be evaluated under different confining pressure conditions;
and the damage evaluation unit is used for obtaining a damage evaluation result of the rock to be evaluated according to the size relation between the damage index and the damage indication threshold value.
8. The evaluation apparatus according to claim 7, wherein the critical damage ratio generation unit is specifically configured to derive a critical damage ratio expression using the statistical damage constitutive model of the rock; the critical damage ratio generating unit is further specifically configured to substitute the distribution parameter and the infinitesimal strength under different confining pressure conditions into the critical damage ratio expression to obtain a critical damage ratio of the rock to be evaluated under different confining pressure conditions.
9. An apparatus for evaluating rock damage, comprising: a memory, a processor and a transceiver connected in sequence, wherein the memory is used for storing a computer program, the transceiver is used for transmitting and receiving messages, and the processor is used for reading the computer program and executing the rock damage assessment method according to any one of claims 1-6.
10. A computer-readable storage medium having stored thereon instructions for performing the method of assessing rock damage according to any one of claims 1 to 6 when run on a computer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110910255.1A CN113591392A (en) | 2021-08-09 | 2021-08-09 | Rock damage assessment method and device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110910255.1A CN113591392A (en) | 2021-08-09 | 2021-08-09 | Rock damage assessment method and device and computer readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113591392A true CN113591392A (en) | 2021-11-02 |
Family
ID=78256543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110910255.1A Pending CN113591392A (en) | 2021-08-09 | 2021-08-09 | Rock damage assessment method and device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113591392A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114397432A (en) * | 2022-03-24 | 2022-04-26 | 四川省水利水电勘测设计研究院有限公司 | Concrete structure surface crack detection device |
| CN116611264A (en) * | 2023-07-14 | 2023-08-18 | 中国矿业大学(北京) | Rock damage evolution model considering initial damage recovery and construction method thereof |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN107515291A (en) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | A kind of construction method of the lower rock freezing-thawing damage constitutive model of confined pressure effect such as |
| CN108287112A (en) * | 2018-01-31 | 2018-07-17 | 成都理工大学 | A method of damage of rock parameter is measured based on triaxial compression test |
-
2021
- 2021-08-09 CN CN202110910255.1A patent/CN113591392A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN107515291A (en) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | A kind of construction method of the lower rock freezing-thawing damage constitutive model of confined pressure effect such as |
| CN108287112A (en) * | 2018-01-31 | 2018-07-17 | 成都理工大学 | A method of damage of rock parameter is measured based on triaxial compression test |
Non-Patent Citations (3)
| Title |
|---|
| 康亚明，等: "岩石的统计损伤本构模型及临界损伤度研究", 《四川大学学报（工程科学版）》, vol. 41, no. 4, pages 42 - 47 * |
| 张毅, 廖华林, 李根生: "岩石连续损伤统计本构模型", 石油大学学报(自然科学版), vol. 28, no. 03, pages 37 - 39 * |
| 曹文贵, 赵明华, 刘成学: "基于统计损伤理论的德鲁克-普拉格岩石强度准则的修正", 水利学报, no. 09, pages 18 - 23 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114397432A (en) * | 2022-03-24 | 2022-04-26 | 四川省水利水电勘测设计研究院有限公司 | Concrete structure surface crack detection device |
| CN114397432B (en) * | 2022-03-24 | 2022-07-08 | 四川省水利水电勘测设计研究院有限公司 | Concrete structure surface crack detection device |
| CN116611264A (en) * | 2023-07-14 | 2023-08-18 | 中国矿业大学(北京) | Rock damage evolution model considering initial damage recovery and construction method thereof |
| CN116611264B (en) * | 2023-07-14 | 2023-10-13 | 中国矿业大学(北京) | Rock damage evolution model considering initial damage recovery and construction method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113591392A (en) | Rock damage assessment method and device and computer readable storage medium | |
| CN110296892B (en) | Method for determining characteristic stress in rock damage evolution process based on energy analysis | |
| CN103197042B (en) | Computing method for representative elementary volume of jointed rock | |
| Munoz et al. | Fracture energy-based brittleness index development and brittleness quantification by pre-peak strength parameters in rock uniaxial compression | |
| CN103115820B (en) | Method for confirming anisotropism of jointed rock mass | |
| CN111060396B (en) | Material mechanical property calibration method based on Ls-Dyna No. 187 material card | |
| CN111238931B (en) | Shale brittleness index evaluation method based on energy evolution | |
| CN110414140A (en) | The model training and test method of wood damage, device, equipment and storage medium | |
| Ren et al. | Numerical investigation of the dynamic increase factor of ultra-high performance concrete based on SHPB technology | |
| CN106769456A (en) | Fully graded concrete long-term behaviour test device and method under a kind of load of long duration | |
| Ai et al. | Crack detection and evolution law for rock mass under SHPB impact tests | |
| CN116662920B (en) | Abnormal data identification method, system, equipment and medium for drilling and blasting method construction equipment | |
| CN112730057A (en) | Method for evaluating long-term service life of brittle solid material under action of compressive load | |
| Gao et al. | Investigation on the expression ability of a developed constitutive model for rocks based on statistical damage theory | |
| CN202903630U (en) | Dynamic indentation testing device of material | |
| Shen et al. | Impact of cyclic wetting and drying on slate properties in the Miaowei Reservoir area, southwest China | |
| CN114492914A (en) | Rock mass characterization unit volume value prediction method and system based on machine learning | |
| CN105424554A (en) | Method for determining fatigue damage degree of metal material based on specific gravity measurement | |
| CN115270490A (en) | Construction method of layered filling body damage dynamic constitutive model | |
| Yang et al. | Prediction method of impact deformation mode based on multimodal fusion with point cloud sequences: Applied to thin-walled structures | |
| Hu et al. | Fatigue life analysis of specimens subjected to infrequent severe loading using a nonlinear kinematic hardening cyclic plasticity model | |
| CN201497684U (en) | Dynamic calibration system of fatigue testing machine | |
| Xie et al. | An automatic pixel-level crack identification method for coals experiencing SHPB impact tests | |
| CN112800583A (en) | Medium and high strain rate dynamic damage constitutive model suitable for hard rock and application method thereof | |
| CN116167147B (en) | Coal rock impact tendency direct index evaluation method based on multi-layer perceptron algorithm |
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 |