CN113012102A - Rock damage evolution process analysis method and device, storage medium and electronic equipment - Google Patents

Rock damage evolution process analysis method and device, storage medium and electronic equipment Download PDF

Info

Publication number
CN113012102A
CN113012102A CN202110170225.1A CN202110170225A CN113012102A CN 113012102 A CN113012102 A CN 113012102A CN 202110170225 A CN202110170225 A CN 202110170225A CN 113012102 A CN113012102 A CN 113012102A
Authority
CN
China
Prior art keywords
rock
damage
target
target rock
variable
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
Application number
CN202110170225.1A
Other languages
Chinese (zh)
Inventor
李曙光
张茹萍
霍润科
任少强
吴应明
谢江胜
仲维玲
王平安
申艳军
王恩龙
张良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway 20th Bureau Group Corp
Original Assignee
China Railway 20th Bureau Group Corp
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 China Railway 20th Bureau Group Corp filed Critical China Railway 20th Bureau Group Corp
Priority to CN202110170225.1A priority Critical patent/CN113012102A/en
Publication of CN113012102A publication Critical patent/CN113012102A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention discloses a rock damage evolution process analysis method, a device, a storage medium and electronic equipment, and relates to the technical field of rock analysis, wherein the method comprises the following steps: obtaining rock sample images corresponding to different stages of a target rock damage process; performing image processing and curve analysis to obtain damage variables of the target rock under the chemical action; obtaining a rock damage variable of the target rock under the combined action of chemistry and load by combining the damage variable of the target rock under the action of load; establishing a segmented statistical damage constitutive model based on a stress-strain curve compaction section of a target rock; and analyzing the damage evolution process of the target rock by utilizing the piecewise statistic damage constitutive model. The invention solves the problem that the rock corroded by the acidic environment is difficult to analyze in the damage evolution process in the prior art, realizes the technical effects of fully reflecting the damage and damage processes of the rock under different acidic corrosion effects and ensuring the accuracy of the analysis of the rock damage evolution process.

Description

Rock damage evolution process analysis method and device, storage medium and electronic equipment
Technical Field
The invention relates to the technical field of rock analysis, in particular to a rock damage evolution process analysis method and device, a storage medium and electronic equipment.
Background
With the rapid development of global industry, the production and living environment of human beings is increasingly worsened, and environmental problems such as air pollution and water pollution are gradually attracted by wide attention of various national scholars. Under an acidic environment, the corrosive destruction effect of materials such as soil, rock, concrete and the like also draws wide attention of numerous scholars, and in environmental geotechnics, the influence of water chemistry on the physical and mechanical properties of the rock is still in an exploration research stage and becomes one of the key contents of the current research.
The rock material has anisotropy, and the mechanical properties of the rock are complex and changeable due to external environmental conditions such as water chemical solution, freeze thawing, high temperature and the like, so that the analysis of the rock damage evolution process is difficult due to the combined action of chemistry and load under the corrosion action of an acid environment, and the damage and damage process of the rock cannot be fully reflected.
Therefore, it is an urgent technical problem to provide a rock damage evolution process analysis method capable of sufficiently reflecting the rock damage destruction process under the corrosion action of an acidic environment.
Disclosure of Invention
The main purposes of the invention are as follows: the method and the device for analyzing the rock damage evolution process, the storage medium and the electronic equipment are provided, and the technical problem that the rock corroded by the acidic environment in the prior art is difficult to analyze the damage evolution process is solved.
In order to achieve the purpose, the invention adopts the following technical scheme:
in a first aspect, the invention provides a method for analyzing a rock damage evolution process, the method comprising the following steps:
obtaining rock sample images corresponding to different stages of a target rock damage process;
performing image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the action of load;
according to the rock damage variable, a segmented statistical damage constitutive model is established based on a stress-strain curve compaction section of the target rock;
and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
Optionally, in the method for analyzing the evolution process of rock damage, the damage process includes a process of corrosion by acid;
the step of obtaining rock sample images corresponding to different stages of the target rock damage process specifically includes:
and carrying out CT scanning on the rock which is not corroded by the acid and the rock at different stages corroded by the acid to obtain rock sample images corresponding to different stages of the target rock damage process.
Optionally, in the method for analyzing a rock damage evolution process, the step of performing image processing and curve analysis on the rock sample image to obtain the damage variable of the target rock under the chemical action specifically includes:
performing image processing on the rock sample image to obtain a CT value frequency distribution curve;
carrying out CT number change analysis and peak type change analysis on the CT value frequency distribution curve to obtain the damage variable of the target rock under the chemical action, wherein the damage variable D of the target rock under the chemical actioncComprises the following steps:
Figure BDA0002939470010000021
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0The density of the target rock without damage is expressed in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
Optionally, in the method for analyzing a rock damage evolution process, the step of performing image processing on the rock sample image to obtain a CT value frequency distribution curve specifically includes:
extracting the CT number of the rock sample image, and performing statistical analysis to obtain a CT number histogram;
and converting the CT number histogram into a curve graph to obtain a CT value frequency distribution curve.
Optionally, in the method for analyzing a rock damage evolution process, the step of obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable under the load action specifically includes:
according to the damage variable of the target rock under the chemical action and the damage variable under the load action, based on the chemical-load coupling damage variable equationObtaining a rock damage variable of the target rock under the combined action of chemistry and load, wherein the damage variable of the target rock under the load is obtained based on Weibull distribution; the damage variable D of the target rock under the loadmComprises the following steps:
Figure BDA0002939470010000031
wherein F represents a randomly distributed variable of the Weibull distribution of the target rock, F0Representing a Weibull distribution parameter of the target rock, and m represents a Weibull distribution statistical parameter;
the chemical-load coupling damage variable equation is as follows:
D=Dm+Dc-DmDc
wherein D represents the rock damage variable of the target rock under the combined action of chemistry and load, and DcRepresenting damage variables of the target rock under chemical action;
the calculation formula for obtaining the rock damage variable D of the target rock under the combined action of chemistry and load is as follows:
Figure BDA0002939470010000032
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0The density of the target rock without damage is expressed in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
Optionally, in the method for analyzing a rock damage evolution process, the step of establishing a piecewise statistical damage constitutive model based on a stress-strain curve compaction section of a target rock according to the rock damage variable specifically includes:
establishing a piecewise statistical damage constitutive equation based on a stress-strain curve compaction section of the target rock;
obtaining a Weibull distribution parameter by using a fitting method according to the rock damage variable;
and substituting the Weibull distribution parameters into the segmentation statistical damage constitutive equation to obtain a segmentation statistical damage constitutive model.
Optionally, in the method for analyzing the rock damage evolution process, before the step of establishing a piecewise statistical damage constitutive equation based on a compressive section of a stress-strain curve of the target rock, the method further includes:
and carrying out uniaxial compression test on the target rock to obtain a stress-strain curve of the target rock.
In a second aspect, the present invention provides an apparatus for analyzing a rock damage evolution process, the apparatus comprising:
the image acquisition module is used for acquiring rock sample images corresponding to different stages of the target rock damage process;
the image processing module is used for carrying out image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
the damage variable acquisition module is used for acquiring a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable under the load action;
the model building module is used for building a segmental statistical damage constitutive model based on a stress-strain curve compaction section of the target rock according to the rock damage variable;
and the damage analysis module is used for analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
In a third aspect, the present invention provides a storage medium having stored thereon a computer program executable by one or more processors to implement the rock damage evolution process analysis method as described above.
In a fourth aspect, the present invention provides an electronic device, which includes a memory and a processor, wherein the memory stores a computer program, and the computer program, when executed by the processor, executes the rock damage evolution process analysis method.
One or more technical solutions provided by the present invention may have the following advantages or at least achieve the following technical effects:
according to the rock damage evolution process analysis method, the rock sample images corresponding to different stages of the target rock damage process are subjected to image processing and curve analysis, the damage variable of the target rock under the chemical action is obtained, then the rock damage variable of the target rock under the combined action of chemistry and load is obtained by combining the damage variable of the target rock under the chemical action, different environment factors of rock damage are fully considered, the damage evolution process of the target rock is analyzed by utilizing a segmented statistical damage constitutive model established based on a stress-strain curve compaction section of the target rock, the damage process of the rock under different acid corrosion actions is fully reflected, and the accuracy of rock damage evolution process analysis is guaranteed.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
Fig. 1 is a schematic flow chart of a rock damage evolution process analysis method according to an embodiment of the present invention;
fig. 2 is a rock sample image obtained in step S1 in the rock damage evolution process analysis method according to the first embodiment of the present invention;
fig. 3 is a frequency distribution graph of CT values obtained in step S2.1 in the rock damage evolution process analysis method according to the first embodiment of the present invention;
fig. 4 is a stress-strain graph obtained in step a1 in the rock damage evolution process analysis method according to the first embodiment of the present invention;
fig. 5 is a functional module schematic diagram of a rock damage evolution process analysis apparatus according to a second embodiment of the present invention.
The implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that all the directional indicators (such as up, down, left, right, front, and rear … …) in the embodiment of the present invention are only used to explain the relative position relationship between the components, the movement situation, etc. in a specific posture (as shown in the drawing), and if the specific posture is changed, the directional indicator is changed accordingly.
In the present invention, unless otherwise expressly stated or limited, the terms "connected," "secured," and the like are to be construed broadly, and for example, "secured" may be a fixed connection, a removable connection, or an integral part; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium; either internal or interactive relationship, unless expressly defined otherwise. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In addition, if there is a description of "first", "second", etc. in an embodiment of the present invention, the description of "first", "second", etc. is for descriptive purposes only and is not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In addition, the meaning of "and/or" appearing throughout includes three juxtapositions, exemplified by "A and/or B" including either A or B or both A and B. In addition, technical solutions between various embodiments may be combined with each other, but must be realized by a person skilled in the art, and when the technical solutions are contradictory or cannot be realized, such a combination of technical solutions should be considered to be absent and not be within the protection scope of the present invention.
Rock, a natural body of land, is a common building material in hydraulic and civil engineering. The deformation and destruction characteristics of rock are not only related to their own mineral components, inter-particle connections, internal microcracks and other complex structures, but also affected by external environmental factors such as confining pressure, temperature, pore water and the like. Over a long period of time, most relevant researchers have conducted extensive studies on water permeability in rock, but less on water chemistry in rock. In the chemical action of water, the research on the corrosive destruction effect of acidic environment on soil, rock, concrete and other materials has attracted much attention.
The analysis of the prior art finds that the theories in the aspects of establishment of a rock damage constitutive model, improvement of a corresponding integral algorithm, nonlinear finite element solution and the like are not mature in the analysis of the rock damage evolution process under the action of water chemistry at present, and are still in an exploration research stage.
In view of the technical problem that the analysis of the damage evolution process of the rock corroded by the acidic environment is difficult in the prior art, the invention provides an analysis method of the damage evolution process of the rock, and the general idea is as follows:
obtaining rock sample images corresponding to different stages of a target rock damage process; performing image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action; obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the action of load; according to the rock damage variable, a segmented statistical damage constitutive model is established based on a stress-strain curve compaction section of the target rock; and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
According to the technical scheme, the rock sample images corresponding to different stages of the target rock damage process are subjected to image processing and curve analysis, after the damage variable of the target rock under the chemical action is obtained, the rock damage variable of the target rock under the combined action of chemistry and load is obtained by combining the damage variable of the target rock under the chemical action, different damaged environmental factors of the rock are fully considered, the damage evolution process of the target rock is analyzed by utilizing a segmented statistical damage constitutive model established based on the stress-strain curve compaction section of the target rock, the damage process of the rock under different acidic corrosion actions is fully reflected, and the accuracy of analysis of the rock damage evolution process is ensured.
Example one
Referring to fig. 1 to 4, a first embodiment of the present invention provides a rock damage evolution process analysis method applicable to an electronic device, where the method, when applied to the electronic device, performs the following steps:
step S1: obtaining rock sample images corresponding to different stages of a target rock damage process;
step S2: performing image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
step S3: obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the action of load;
step S4: according to the rock damage variable, a segmented statistical damage constitutive model is established based on a stress-strain curve compaction section of the target rock;
step S5: and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
In the damage evolution process, an equation is generally used to describe the rule that the damage inside the structural material changes along with the change of the action of external factors, such as load, temperature change, corrosion and other external factors. Different damage variables are selected, damage evolution equations are different, and the damage evolution equations can reflect the damage change development condition of the material. The damage model established based on the damage evolution equation is a physical model used for describing the damage development and change rule. This example analyzes the evolution process of the damage of the rock corroded by acid.
The rock damage evolution process analysis method provided by this embodiment is described in detail below with reference to fig. 1, and specifically includes the following steps:
step S1: and obtaining rock sample images corresponding to different stages of the target rock damage process.
Specifically, the damage process includes a process of corrosion by acid; the step S1 may include:
step S1.1: and carrying out CT scanning on the rock which is not corroded by the acid and the rock at different stages corroded by the acid to obtain rock sample images corresponding to different stages of the target rock damage process.
In a specific embodiment, rocks corroded by different acid solutions can be used as the rock sample, and the acid solution can be hydrochloric acid solution with the pH of 1 or 3, or sulfuric acid solution with the pH of 1 or 3. The rock may be laminated rock such as sandstone, conglomerate, shale, etc., or metamorphic rock such as marble rock, slate, etc. The sandstone is a sedimentary clastic rock composed of quartz, feldspar and the like, and is formed by stacking weathered corrosion, transportation and the like in a plurality of hilly areas. The different stages of the damage process may be different damage days, e.g. 0, 30, 90 and 180 days. Carrying out CT scanning on the rock which is not corroded by acid to obtain a CT image of the original rock; and performing CT scanning on rocks corroded by acid at different stages, namely performing CT scanning on the middle layer of rocks corroded by the same acid solution but with different corrosion days to obtain rock sample CT images, wherein the original rock CT image and the rock sample CT image are used as rock sample images corresponding to different stages of the target rock damage process together.
In this embodiment, as shown in fig. 2, the rock sample image obtained in this embodiment is shown; the acid solution is sulfuric acid solution with pH 1, and sandstone corroded by the sulfuric acid solution with pH 1 is used as the rock sample. Performing CT scanning on the sandstone which is not corroded by the acid solution to obtain a CT image of the original rock shown in figure 2 (a); carrying out CT scanning on the sandstone corroded by the acid solution for 30 days to obtain a rock sample CT image shown in figure 2 (b); carrying out CT scanning on the sandstone corroded by the acid solution for 90 days to obtain a rock sample CT image shown in figure 2 (c); CT scanning was performed on sandstone corroded by acid solution for 180 days to obtain a CT image of the rock sample as shown in fig. 2 (d).
By carrying out CT scanning on the rock sample, the rock sample image is obtained, the advantages of nondestructive testing of the CT scanning are fully utilized, and the processing precision and accuracy of the subsequent steps are improved.
Step S2: and carrying out image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action.
Specifically, the step S2 may include:
step S2.1: and carrying out image processing on the rock sample image to obtain a CT value frequency distribution curve.
Specifically, the step S2.1 may include:
step S2.1.1: and extracting the CT number of the rock sample image, and performing statistical analysis to obtain a CT number histogram.
In a specific embodiment, ImageJ software is used for extracting the CT number of the rock sample image, and statistical analysis is carried out to obtain a CT number histogram.
Step S2.1.2: and converting the CT number histogram into a curve graph to obtain a CT value frequency distribution curve.
In a specific embodiment, ImageJ image processing software is used for converting the CT number histogram into a curve graph to obtain a CT value frequency distribution curve.
In this embodiment, the above processing is performed on the rock sample image of fig. 2, and then the frequency distribution graph of the CT value shown in fig. 3 is obtained. Performing CT number extraction and histogram conversion on the original rock CT image shown in FIG. 2(a) to obtain a CT value frequency distribution curve shown in FIG. 3 (a); performing CT number extraction and histogram conversion on the rock sample CT image shown in the figure 2(b) to obtain a CT value frequency distribution curve shown in the figure 3 (b); performing CT number extraction and histogram conversion on the rock sample CT image shown in the figure 2(c) to obtain a CT value frequency distribution curve shown in the figure 3 (c); CT number extraction and histogram conversion are performed on the rock sample CT image shown in fig. 2(d), and a CT value frequency distribution curve shown in fig. 3(d) is obtained.
Step S2.2: and carrying out CT number change analysis and peak type change analysis on the CT value frequency distribution curve to obtain the damage variable of the target rock under the chemical action.
Specifically, the damage evolution law of the rock material is analyzed through the change of the CT number and the change of the peak type in the CT value frequency distribution curve, and the damage variable formula of the rock material with the bimodal distribution of the CT value frequency distribution curve is deduced and determined.
Specifically, the damage variable D of the target rock under the chemical actioncComprises the following steps:
Figure BDA0002939470010000091
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0The density of the target rock without damage is expressed in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
In this embodiment, the original rock CT image and the rock sample CT image are obtained in step 1, the uncorroded sandstone is subjected to CT scanning, the CT number is extracted by ImageJ software, and then the average value is taken to obtain H1Performing CT scanning on corroded sandstone, extracting CT number by ImageJ software, and averaging to obtain H2The CT number of the matrix material of the sandstone is determined according to the type of the sandstone, and an empirical value 25003000 is generally taken.
Analyzing the CT value frequency distribution curve shown in FIG. 3, wherein the CT value frequency distribution curve of the sandstone not corroded by acid shown in FIG. 3(a) is a unimodal curve and is mainly caused by the defects of the material; the peak patterns of the CT value frequency distribution curves of the acid-corroded sandstone are gradually distributed in a double-peak manner, the sandstone mineral structure distribution is changed due to cavities caused by dissolution of sandstone mineral particles after the acid corrosion, and the cavity damage is more and more dominant in the sandstone damage along with the severity of the acid corrosion, and is reflected in the change of the peak patterns in the CT value frequency distribution curves. The change of the CT number and the change of the peak shape indicate that the damage of the sandstone corroded by acid is not only the original single damage caused by the material defect, but also the cavity damage caused by the dissolution of the mineral particles, so the bimodal distribution damage variable D of the sandstone corroded by acid in the embodiment is deduced and determinedc
Step S3: and obtaining the rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the load action.
Specifically, the step S3 may include:
and obtaining the rock damage variable of the target rock under the combined action of chemistry and load based on a chemical-load coupling damage variable equation according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the load action, wherein the damage variable of the target rock under the load action is obtained based on Weibull distribution.
Specifically, based on Weibull distribution, obtaining damage variable D of the target rock under the action of loadmComprises the following steps:
Figure BDA0002939470010000101
wherein F represents the Weibull's score of the target rock characterizing the strength of the rock infinitesimal bodyRandomly distributed variation of cloth, F0The Weibull distribution parameter represents the characteristic rock mechanical characteristics of the target rock, and m represents the Weibull distribution statistical parameter of the characteristic rock mechanical characteristics of the target rock;
in the specific implementation mode, Weibull distribution is a theoretical basis for researching material reliability and testing material service life, and the invention introduces Weibull distribution to calculate damage variable D of rock sample under loadm
Specifically, the chemical-load coupling damage variable equation is as follows:
D=Dm+Dc-DmDc
wherein D represents the rock damage variable of the target rock under the combined action of chemistry and load, and DcRepresenting damage variables of the target rock under chemical action;
the calculation formula for obtaining the rock damage variable D of the target rock under the combined action of chemistry and load is as follows:
Figure BDA0002939470010000111
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0The density of the target rock without damage is expressed in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
In a specific implementation mode, a generalized damage variable is obtained through an equivalent reaction principle, and a chemical-load coupling damage variable equation is introduced to obtain a rock damage variable D of a rock sample under the combined action of chemistry and load.
Step S4: and according to the rock damage variable, establishing a segmental statistical damage constitutive model based on the stress-strain curve compaction section of the target rock.
Specifically, the step S4 may include:
step S4.1: and establishing a piecewise statistic damage constitutive equation based on the stress-strain curve compaction section of the target rock.
In a specific embodiment, the stress-strain curve may be divided into four phases: the stress-strain curve based damage constitutive equation is established by utilizing the concave characteristic of the compression section of the stress-strain curve:
Figure BDA0002939470010000112
wherein σ1Representing the stress when the rock is undamaged, v representing the Poisson's ratio when the rock is undamaged, σ3Representing the principal stress when the rock is undamaged, i.e. the surrounding pressure of the triaxial test, E representing the modulus of elasticity in the undamaged state of the rock, ε1Indicating the strain when the rock is not damaged,
Figure BDA0002939470010000113
indicating the internal angle of friction of the rock when it is undamaged.
In this embodiment, the piecewise statistical damage constitutive equation established based on the stress-strain curve compaction section of the target rock is:
when epsilon1≤ε1cWhen the temperature of the water is higher than the set temperature,
Figure BDA0002939470010000121
when epsilon11cWhen the temperature of the water is higher than the set temperature,
Figure BDA0002939470010000122
wherein σ1Representing the stress when the target rock is undamaged, v representing the Poisson's ratio, σ, when the target rock is undamaged3Which represents the principal stress when the target rock is not damaged, i.e., the triaxial test ambient pressure, E represents the modulus of elasticity of the target rock in a non-damaged state,ε1representing the strain of the target rock when it is not damaged,
Figure BDA0002939470010000124
representing the internal friction angle, σ, of the target rock when undamaged1cRepresenting the stress at the onset of damage, epsilon, of the target rock1cRepresenting the strain at the onset of damage to the target rock.
Step S4.2: and obtaining a Weibull distribution parameter by using a fitting method according to the rock damage variable.
In a specific embodiment, first, the demonstration statistical damage constitutive equation of step S4.1 is deformed.
Specifically, logarithms are taken for two sides of the demonstration statistical damage constitutive equation in the step S4.1 to obtain:
Figure BDA0002939470010000123
and (4) taking logarithm of two sides after the term of the formula is shifted to obtain:
Figure BDA0002939470010000131
in order to ensure that the water-soluble organic acid,
A=-mln(F0),
Figure BDA0002939470010000132
Figure BDA0002939470010000133
the above-mentioned demonstration statistical damage constitutive equation can be expressed as:
Y=mX+A。
then, a fitting method is used for obtaining the Weibull distribution parameters.
Concretely, stress and strain data in a rock stress-strain curve are substituted into a demonstration statistical damage bookForming an equation, and calculating to obtain an (X, Y) sequence; performing linear regression analysis on the (X, Y) sequence to obtain a parameter A and a value of a Weibull distribution statistical parameter m; then calculating the Weibull distribution parameter F of the rock according to the following formula0
Figure BDA0002939470010000134
In this embodiment, the piecewise statistical damage constitutive equation of step S4.1 is deformed.
Specifically, when epsilon1≤ε1cThen, obtaining:
Figure BDA0002939470010000135
in order to ensure that the water-soluble organic acid,
A1=-m1 ln(F01),
Figure BDA0002939470010000136
Figure BDA0002939470010000141
the above piecewise statistical impairment constitutive equation can be expressed as:
Y=m1 X+A1
then, a fitting method is used for obtaining the Weibull distribution parameters.
Specifically, linear fitting is carried out on data of a stress-strain curve compaction section of the target rock, namely stress and strain data of the stress-strain curve compaction section of the target rock are substituted into a piecewise statistics damage constitutive equation, and an (X, Y) sequence is calculated; then carrying out linear regression analysis on the (X, Y) sequence to obtain a parameter A1And a Weibull distribution statistical parameter m1A value of (d); then calculating to obtain a Weibull distribution parameter F of the target rock according to the following formula01
Figure BDA0002939470010000142
Specifically, when epsilon11cThen, obtaining:
Figure BDA0002939470010000143
in order to ensure that the water-soluble organic acid,
A2=-m2 ln(F02),
Figure BDA0002939470010000144
Figure BDA0002939470010000145
the above piecewise statistical impairment constitutive equation can be expressed as:
Y=m2 X+A2
then, a fitting method is used for obtaining the Weibull distribution parameters.
Specifically, linear fitting is carried out on data of a stress-strain curve compaction section of the target rock, namely stress and strain data of the stress-strain curve compaction section of the target rock are substituted into a piecewise statistics damage constitutive equation, and an (X, Y) sequence is calculated; then carrying out linear regression analysis on the (X, Y) sequence to obtain a parameter A2And a Weibull distribution statistical parameter m2A value of (d); then calculating to obtain a Weibull distribution parameter F of the target rock according to the following formula02
Figure BDA0002939470010000151
Step S4.3: and substituting the Weibull distribution parameters into the segmentation statistical damage constitutive equation to obtain a segmentation statistical damage constitutive model.
In a specific embodiment, a Weibull distribution statistical parameter m and a Weibull distribution parameter F are used0And substituting the model into the demonstration statistical damage constitutive equation to obtain a demonstration statistical damage constitutive model.
In this embodiment, the weibull distribution statistical parameter m obtained in step S4.2 is used1And the Weibull distribution parameter F of the target rock01Statistical parameter m of Weibull distribution2And the Weibull distribution parameter F of the target rock02And (4) substituting the equation into the piecewise statistic damage constitutive equation in the step (S4.1) to determine a specific piecewise statistic damage constitutive equation, wherein the equation is an expression of the piecewise statistic damage constitutive model in the embodiment.
In the embodiment, a sectional statistical damage constitutive equation of the sandstone stress-strain curve compaction section in the acidic environment is considered, a sectional statistical damage constitutive model is further established, and the obtained sectional statistical damage constitutive model can fully reflect the damage and damage processes of the sandstone under the corrosion action of different acidic environments.
In one embodiment, before step S4.1, the method may further comprise a method of obtaining a stress-strain curve of the target rock, which may comprise:
step A1: and carrying out uniaxial compression test on the target rock to obtain a stress-strain curve of the target rock.
As shown in fig. 4, which is a stress-strain curve of the target rock damage process obtained in the present embodiment, the stress-strain curve can be divided into four stages: a compacting section, an elastic section, a plastic section, and a breaking and destroying section, wherein the compacting section is concave. Wherein E-C represents an experimental value, T-C represents a true value, the abscissa of the graph represents strain, the ordinate represents stress, the unit is MPa, and the solid line in the graph represents the true values of the stress-strain curves of the sandstone obtained in the step when the sandstone is damaged for 30 days, 90 days and 180 days.
In another embodiment, after step S4.3, a method for verifying the piecewise statistical impairment constitutive model may be further included, and the method may include:
step A2: and carrying out uniaxial compression test on the rock sample to obtain a stress-strain curve of the rock sample, and verifying the piecewise statistic damage constitutive model obtained by the embodiment on the basis of the stress-strain curve of the rock sample.
As shown in fig. 4, the dotted line in the graph is the experimental values of the stress-strain curves of the sandstone obtained in the step when the sandstone is damaged for 30 days, 90 days and 180 days, and the real value of the stress-strain curve in the step a1 is combined, so that the accuracy of the piecewise statistic damage constitutive model obtained in the embodiment is high, and the evolution process of rock damage can be fully reflected.
Step S5: and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
In this embodiment, the damage evolution process of the sandstone corroded by acid is analyzed by using the piecewise statistical damage constitutive model obtained in the above steps.
According to the rock damage evolution process analysis method provided by the embodiment, image processing and curve analysis are performed on rock sample images corresponding to different stages of a target rock damage process, rock damage variables of the target rock under the combined action of chemistry and load are obtained based on a chemical-load coupling damage variable equation, different environmental factors of rock damage are fully considered, the damage evolution process of the target rock is analyzed through a segmented statistical damage constitutive model established based on a stress-strain curve compaction section of the target rock, the damage process of the rock under different acidic corrosion actions is fully reflected, and the accuracy of rock damage evolution process analysis is guaranteed.
Example two
Based on the same inventive concept, referring to fig. 5, a second embodiment of the present invention provides an apparatus for analyzing a rock damage evolution process, which is described in detail with reference to the schematic diagram of functional modules of fig. 5, and includes:
the image acquisition module is used for acquiring rock sample images corresponding to different stages of the target rock damage process;
the image processing module is used for carrying out image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
the damage variable acquisition module is used for acquiring a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable under the load action;
the model building module is used for building a segmental statistical damage constitutive model based on a stress-strain curve compaction section of the target rock according to the rock damage variable;
and the damage analysis module is used for analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
For the sake of brevity of the description, repeated descriptions are not repeated in this embodiment.
In the device for analyzing the rock damage evolution process provided by this embodiment, after the image obtaining module obtains rock sample images corresponding to different stages of the target rock damage process, the image processing module performs image processing and curve analysis on the rock sample images to obtain damage variables of the target rock under chemical action, the damage variable obtaining module combines the damage variables of the target rock under chemical action to obtain rock damage variables of the target rock under combined action of chemistry and load, different environmental factors of rock damage are fully considered, the model building module builds a segmented statistical damage constitutive model based on a stress-strain curve compaction section of the target rock, and the damage analysis module analyzes the damage evolution process of the target rock by using the segmented statistical damage constitutive model to fully reflect the damage destruction process of the rock under different acidic corrosion actions, the accuracy of rock damage evolution process analysis is guaranteed.
EXAMPLE III
Based on the same inventive concept, the present embodiment provides a computer-readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card-type memory (e.g., SD or DX memory, etc.), a Random Access Memory (RAM), a Static Random Access Memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, a server, an App, etc., on which a computer program is stored, which when executed by a processor, may implement the following method steps:
obtaining rock sample images corresponding to different stages of a target rock damage process;
performing image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the action of load;
according to the rock damage variable, a segmented statistical damage constitutive model is established based on a stress-strain curve compaction section of the target rock;
and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
The specific embodiment process of the above method steps can be referred to as embodiment one, and the detailed description of this embodiment is not repeated herein.
Example four
Based on the same inventive concept, the present embodiment provides an electronic device, which may be a mobile phone, a computer, or a tablet computer, and the electronic device includes a memory and a processor, where the memory stores a computer program, and the computer program is executed by the processor to implement the rock damage evolution process analysis method as described in the above embodiment.
It is understood that the electronic device may also include multimedia components, input/output (I/O) interfaces, and communication components.
Wherein the processor is used for executing all or part of the steps in the rock damage evolution process analysis method as described in the first embodiment. The memory is used to store various types of data, which may include, for example, instructions for any application or method in the electronic device, as well as application-related data.
The Processor may be an Application Specific Integrated Circuit (ASIC), a Digital Signal Processor (DSP), a Digital Signal Processing Device (DSPD), a Programmable Logic Device (PLD), a Field Programmable Gate Array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic components, and is configured to perform all or part of the steps of the rock damage evolution process analysis method according to the first embodiment.
The Memory may be implemented by any type of volatile or non-volatile Memory device or combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk or optical disk.
The multimedia components may include a screen, which may be a touch screen, and an audio component for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may further be stored in a memory or transmitted through a communication component. The audio assembly also includes at least one speaker for outputting audio signals.
The I/O interface provides an interface between the processor and other interface modules, such as a keyboard, a mouse, buttons, etc. These buttons may be virtual buttons or physical buttons.
The communication component is used for carrying out wired or wireless communication between the electronic equipment and other equipment. Wireless Communication, such as Wi-Fi, bluetooth, Near Field Communication (NFC), 2G, 3G or 4G, or a combination of one or more of them, so that the corresponding Communication component may include: Wi-Fi module, bluetooth module, wireless communication modules such as NFC module.
It should be noted that, since the drawings in the specification should not be colored or modified, it is difficult to display the parts of the drawings in the present invention where the parts are clearly distinguished from each other, and if necessary, a color picture can be provided.
The above description is only an alternative embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes performed by the present specification and drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.

Claims (10)

1. A rock damage evolution process analysis method, characterized in that the method comprises the following steps:
obtaining rock sample images corresponding to different stages of a target rock damage process;
performing image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
obtaining a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the action of load;
according to the rock damage variable, a segmented statistical damage constitutive model is established based on a stress-strain curve compaction section of the target rock;
and analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
2. The method for analyzing the evolution process of rock damage according to claim 1, wherein the damage process comprises a process of acid corrosion;
the step of obtaining rock sample images corresponding to different stages of the target rock damage process specifically includes:
and carrying out CT scanning on the rock which is not corroded by the acid and the rock at different stages corroded by the acid to obtain rock sample images corresponding to different stages of the target rock damage process.
3. The rock damage evolution process analysis method of claim 1, wherein the step of performing image processing and curve analysis on the rock sample image to obtain the damage variable of the target rock under the chemical action specifically comprises:
performing image processing on the rock sample image to obtain a CT value frequency distribution curve;
carrying out CT number change analysis and peak type change analysis on the CT value frequency distribution curve to obtain the damage variable of the target rock under the chemical action, wherein the damage variable D of the target rock under the chemical actioncComprises the following steps:
Figure FDA0002939467000000011
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0The density of the target rock without damage is expressed in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
4. The method for analyzing a rock damage evolution process according to claim 3, wherein the step of performing image processing on the rock sample image to obtain the CT value frequency distribution curve specifically comprises:
extracting the CT number of the rock sample image, and performing statistical analysis to obtain a CT number histogram;
and converting the CT number histogram into a curve graph to obtain a CT value frequency distribution curve.
5. The method for analyzing the evolution process of rock damage according to claim 1, wherein the step of obtaining the rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the action of chemistry and load comprises:
obtaining a rock damage variable of the target rock under the combined action of chemistry and load based on a chemical-load coupling damage variable equation according to the damage variable of the target rock under the chemical action and the damage variable of the target rock under the load action, wherein the damage variable of the target rock under the load action is obtained based on Weibull distribution; the damage variable D of the target rock under the loadmComprises the following steps:
Figure FDA0002939467000000021
wherein F represents a randomly distributed variable of the Weibull distribution of the target rock, F0Representing a Weibull distribution parameter of the target rock, and m represents a Weibull distribution statistical parameter;
the chemical-load coupling damage variable equation is as follows:
D=Dm+Dc-DmDc
wherein D represents the rock damage variable of the target rock under the combined action of chemistry and load, and DcRepresenting damage variables of the target rock under chemical action;
the calculation formula for obtaining the rock damage variable D of the target rock under the combined action of chemistry and load is as follows:
Figure FDA0002939467000000022
where ρ isrRepresents the density of the target rock matrix material in g/cm3,ρ0Indicating that the target rock is freeDensity at damage in g/cm3,H1Represents the CT number, H, of the target rock when damaged2Represents the CT number, H, of the target rock when undamagedrRepresents the CT number of the target rock matrix material.
6. The method for analyzing the evolution process of rock damage according to claim 1, wherein the step of establishing a piecewise statistical damage constitutive model based on the stress-strain curve compaction section of the target rock according to the rock damage variable specifically comprises:
establishing a piecewise statistical damage constitutive equation based on a stress-strain curve compaction section of the target rock;
obtaining a Weibull distribution parameter by using a fitting method according to the rock damage variable;
and substituting the Weibull distribution parameters into the segmentation statistical damage constitutive equation to obtain a segmentation statistical damage constitutive model.
7. The rock damage evolution process analysis method of claim 6, wherein prior to the step of establishing a piecewise statistical damage constitutive equation based on the stress-strain curve compaction section of the target rock, the method further comprises:
and carrying out uniaxial compression test on the target rock to obtain a stress-strain curve of the target rock.
8. An apparatus for analyzing a rock damage evolution process, the apparatus comprising:
the image acquisition module is used for acquiring rock sample images corresponding to different stages of the target rock damage process;
the image processing module is used for carrying out image processing and curve analysis on the rock sample image to obtain a damage variable of the target rock under the chemical action;
the damage variable acquisition module is used for acquiring a rock damage variable of the target rock under the combined action of chemistry and load according to the damage variable of the target rock under the chemical action and the damage variable under the load action;
the model building module is used for building a segmental statistical damage constitutive model based on a stress-strain curve compaction section of the target rock according to the rock damage variable;
and the damage analysis module is used for analyzing the damage evolution process of the target rock by utilizing the segmented statistical damage constitutive model.
9. A storage medium having stored thereon a computer program executable by one or more processors to implement the rock damage evolution process analysis method of any one of claims 1 to 7.
10. An electronic device, characterized in that the electronic device comprises a memory and a processor, the memory having stored thereon a computer program which, when executed by the processor, implements the rock damage evolution process analysis method according to any one of claims 1 to 7.
CN202110170225.1A 2021-02-08 2021-02-08 Rock damage evolution process analysis method and device, storage medium and electronic equipment Pending CN113012102A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110170225.1A CN113012102A (en) 2021-02-08 2021-02-08 Rock damage evolution process analysis method and device, storage medium and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110170225.1A CN113012102A (en) 2021-02-08 2021-02-08 Rock damage evolution process analysis method and device, storage medium and electronic equipment

Publications (1)

Publication Number Publication Date
CN113012102A true CN113012102A (en) 2021-06-22

Family

ID=76385255

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110170225.1A Pending CN113012102A (en) 2021-02-08 2021-02-08 Rock damage evolution process analysis method and device, storage medium and electronic equipment

Country Status (1)

Country Link
CN (1) CN113012102A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113484499A (en) * 2021-07-26 2021-10-08 东北大学 Rock deformation stage identification method based on temperature change characteristics
CN113610914A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone pore structure analysis method, device, equipment and storage medium
CN113610758A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Rock damage analysis method, device and equipment based on medical CT scanning
CN113607764A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone sample scanning test method and system
CN113742958A (en) * 2021-08-05 2021-12-03 华北理工大学 Calculation method of rock digital representation model based on matter element theory

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SHUGUANG LI: "Experimental Study on Physicomechanical Properties of Sandstone under Acidic Environment", 《DAMAGE AND FRACTURE BEHAVIOR OF ROCK》, pages 5 *
SHUGUANG LI: "Mechanical Properties of Acid-corroded Sandstone Under Uniaxial Compression", 《 ROCK MECHANICS AND ROCK ENGINEERING》, pages 2 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113610914A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone pore structure analysis method, device, equipment and storage medium
CN113610758A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Rock damage analysis method, device and equipment based on medical CT scanning
CN113607764A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone sample scanning test method and system
CN113484499A (en) * 2021-07-26 2021-10-08 东北大学 Rock deformation stage identification method based on temperature change characteristics
CN113484499B (en) * 2021-07-26 2022-03-25 东北大学 Rock deformation stage identification method based on temperature change characteristics
CN113742958A (en) * 2021-08-05 2021-12-03 华北理工大学 Calculation method of rock digital representation model based on matter element theory

Similar Documents

Publication Publication Date Title
CN113012102A (en) Rock damage evolution process analysis method and device, storage medium and electronic equipment
Fichant et al. Isotropic and anisotropic descriptions of damage in concrete structures
Lv et al. Effects of internal pores on the mechanical properties of marine calcareous sand particles
Zhou et al. An improved grain-based numerical manifold method to simulate deformation, damage and fracturing of rocks at the grain size level
CN103076347B (en) Measurement method for mechanical injury of brittle material based on in-situ X-ray tomography
CN102081030A (en) Geotechnical model test system based on macro-micro mechanics as well as elaboration test method
Askaripour et al. A review of relationship between texture characteristic and mechanical properties of rock
Liu et al. Analysis of the compression of structured soils using the disturbed state concept
Zhu et al. Maximum entropy-based stochastic micromechanical model for a two-phase composite considering the inter-particle interaction effect
Szeląg Evaluation of cracking patterns in cement composites—From basics to advances: A review
Huang et al. A plastic-damage model for rock-like materials focused on damage mechanisms under high pressure
CN115719017A (en) Seawater sea sand concrete multi-physical field coupling analysis and construction quality control method
CN114358091B (en) Pile damage identification method, equipment and medium based on convolutional neural network
Ali et al. Analytical damage model for predicting coal failure stresses by utilizing acoustic emission
Zhou et al. A damage-based numerical manifold approach to crack propagation in rocks
Shu et al. Stress–Strain Strength Characteristics of Undisturbed Granite Residual Soil Considering Different Patterns of Variation of Mean Effective Stress
Albadri et al. The importance of incorporating hysteresis effect in determining shear strength of unsaturated soil
CN114459657B (en) Impact load automatic identification method, electronic device and storage medium
CN116556718A (en) Repairing method, repairing device, repairing equipment and repairing storage medium for building structural component
CN106934729B (en) Building detection and identification method and device
Fu et al. Slope stability analysis based on big data and convolutional neural network
Huie et al. Damage evolution process of cement-stabilized soil based on deformation and microstructure analysis
Xia et al. Strength Deterioration Model of Soft Rock Considering Mesoscopic Bonding–Expansion Coupling Mechanism under Freeze–Thaw Cycles
Rezvani The effect of particle breakage on the dilatancy and shear behavior of marine calcareous and siliceous deposits
CN113607618A (en) Rock damage analysis method, device and equipment based on industrial CT scanning

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