This application claims priority to Chinese Patent Application No. 202110549984.9, filed on May 20, 2021, the contents of which are hereby incorporated by reference.
The application relates to a method for determining rock deformation and failure, and in particular to a method for determining a whole macro-micro process of rock deformation and failure based on a four-parameter test.
At present, the deformation and failure of rocks is the main factor that induces the disasters in practical engineering such as mines, underground spaces, tunnels and dams. As a service life of rock engineering increases, accumulation of rock deformation and damage, which may result in the collapse of engineering structure corresponding to the rock mass, threatens people's lives and engineering safety.
This deformation process may associate with the macroscopic damage and microscopic structural changes of rock. Rock deformation and failure have attracted much attention in mines, underground spaces, tunnels, dams and other projects. Because of the complexity of rock deformation and failure, to reveal its mechanism, it is necessary to analyze it from both macro and micro perspectives, and build a bridge between macro failure process and micro structural change. However, at present, macro and micro parameters are measured and analyzed independently, but no connection has been established. Therefore, how to provide a method to establish a quantitative relationship between macro and micro in the whole process of rock deformation and failure is a research direction of this industry, which is important to provide theoretical support for follow-up researches.
Aiming at the problems existing in the prior art, the present disclosure provides method for determining a whole macro-micro process of rock deformation and failure based on a four-parameter test, which establishes the quantitative relationship between macro and micro of the whole process of rock deformation and failure, and provides theoretical support for subsequent research.
In order to achieve the above objectives, the technical scheme adopted by the disclosure is as follows: a method for determining a whole macro-micro process of rock deformation and failure based on a four-parameter test which specifically comprises following steps:
S1: sampling an engineering rock mass to be tested, and processing a sampled rock mass into a cylindrical specimen;
S2: placing the cylindrical specimen on a testing machine in a compression test system, sticking a deformation sensor and an acoustic emission probe on a surface of the cylindrical specimen, then starting a compression test, and collecting acoustic emission data through the acoustic emission probe and deformation data on the surface of the cylindrical specimen through the deformation sensor while the compression test is being carried out;
S3: calculating the deformation data collected in S2 according to a finite deformation theory, and obtaining a parameter-mean rotation angle θ, which characterizes a macroscopic deformation characteristic of materials at each stress level and is specifically:
F
j
i
=S
j
i
+R
j
i (1),
where Fji is a deformation gradient; orthogonal transformation Rji is a rotation tensor; symmetric transformation Sji is a strain tensor, and an expression of Sji is as follows:
during a test measurement, calculating out a strain component based on a small deformation theory, namely:
where ui|j is a covariant derivative of displacement and εji is a small deformation strain;
combined a small deformation strain component with a finite deformation strain component to get:
S
j
i=εji−(1−cos θ)LkiLji (4),
where Ljk is azimuth tensor of a rotation axis;
according to Hooke's law, one-dimensional elastic lossless constitutive formula is
σ=ES (5),
from formula (4) and formula (5), getting
σ=Eεji−E(1−cos θ)LkiLjk (6),
where σ is a stress;
extending formula (6) to a three-dimensional state and writing the formula (6) as
in a triaxial test, σ22=σ33=σcon, combined with formula (10), getting:
in a triaxial compression test, there being an assumption as follows:
(L21)2=(L32)2=(L13)2 (9),
writing formula (8) as
obtaining a formula for calculating the mean rotation angle θ from formula (12):
so as to calculate the mean rotation angle θ;
S4: using Grassberger-Procaccia (G-P) algorithm on the acoustic emission data collected in S2 to calculate a fractal dimension of a temporal distribution DT of an acoustic emission signal and calculating a fractal dimension of a spatial distribution DS according to a spatial projection method; specifically:
taking time series of the acoustic emission signal as a research object, then corresponding each time series to a series set with a capacity of n:
X={x
1
,x
2
, . . . ,x
n} (12),
constructing formula (12) as a m-dimensional phase space (m<n), firstly, taking m numbers as a vector of m-dimensional space
X
1
={x
1
,x
2
,x
3
, . . . ,x
m} (13),
then shifting one data to the right and taking m numbers again to form another vector, and so on to form N=n−m+1 vectors; the corresponding correlation function is:
where H is Heaviside function; r is a given scale; when assigning a value to scale r, making r=kr0 in order to avoid dispersion, where k is taken as a scale coefficient and
obtaining n points in a double logarithmic coordinate system, and performing data fitting on T1 points. If the result is a straight line, it shows that the acoustic emission series has fractal characteristics in a given scale range, and a slope of the straight line is the fractal dimension of the temporal distribution DT of the acoustic emission parameter, namely
D
T=1 gW(r)/1 g(r) (15).
For DS, the space box dimension is used. The box dimension is defined as:
N(r)=Cr−D
where N(r) is the number of discrete bodies whose characteristic size is greater than r, C is a material constant, and the other form of the above formula is the number-radius relation as follows:
M(r)=Cr−D
where r is different radii covering natural discrete bodies, and M(r) is the number of discrete bodies covered in a circle with radius of r. Taking the logarithm on both sides to get:
1 gM(r)=1 gC+DS1 g(r) (18),
where DS is the fractal dimension of the spatial distribution;
S5: carrying out a scanning electron microscope (SEM) test on a fracture surface of the sample after the compression test is completed, to obtain a microscopic morphology of the fracture surface, observing the morphology of the fracture surface and calculating the fractal dimension DA of the fracture surface.
The number of units needed to cover an image in units of δ is N(δ), DA=−log (N(δ))/log δ.
S6: because of a correspondence of a change of the mean rotation angle θ to each process of rock deformation, including a compaction stage, a linear stage and a plastic yield stage in a compression process, finally obtaining a mathematical trend relationship between θ and DT, DS and DA through comprehensively analyzing the obtained fractal dimension of the temporal distribution DT of the acoustic emission, the fractal dimension of the spatial distribution DS of the acoustic emission and the fractal dimension DA of the fracture surface at each stress level (prior to a peak strength) and the mean rotation angle θ at a same stress level, as shown in a following formula,
θ=a*DT+b*DS+c*DA (19),
eventually, obtaining values of a, b and c, so as to establish the quantitative relationship between macro and micro in the whole process of rock deformation and failure.
Further, a height of the cylindrical specimen is 100 mm and a diameter of the cylindrical specimen is 50 mm.
Further, the deformation data includes axial deformation and circumferential deformation.
Compared with the prior art, the disclosure firstly obtains acoustic emission data and deformation data of the sample through the deformation sensor and the acoustic emission probe in the process of the sample compression test, and then calculates the deformation data according to the finite deformation theory to obtain a parameter-mean rotation angle θ which represents the macroscopic deformation characteristics of the material at each stress level; G-P algorithm is used to calculate the acoustic emission data, and the fractal dimension of the temporal distribution DT of acoustic emission signal is obtained, and the fractal dimension of the spatial distribution DS is calculated according to the spatial projection method. After the compression test, the microscopic morphology of the fracture surface is obtained by a scanning electron microscope (SEM) test, and the fractal dimension DA of the fracture surface is calculated out. Finally, the mathematical trend relationship between θ and DT, DS and DA is obtained through comprehensively analysing the obtained fractal dimension of the temporal distribution DT of the acoustic emission, the fractal dimension of the spatial distribution DS of the acoustic emission and fractal dimension DA of the fracture surface at each stress level (prior to a peak strength) and the mean rotation angle θ at a same stress level, thus establishing a quantitative relationship between macro and micro in the whole process of rock deformation and failure and providing theoretical support for follow-up researches.
The technical solutions in the embodiments of the present disclosure will be clearly and completely described below. Obviously, the described embodiments are only part of the embodiments in this disclosure, but not all of them. Based on the embodiment in this disclosure, all other embodiments obtained by ordinary technicians in this field without creative effort are within the scope of protection in this disclosure.
As shown in
S1: firstly, sampling an engineering rock mass to be tested, and processing a sampled rock mass into a cylindrical specimen whose height is 100 mm high and diameter is 50 mm.
S2: placing the cylindrical specimen on a testing machine in a compression test system, sticking a deformation sensor and an acoustic emission probe on a surface of the cylindrical specimen, then starting a compression test, and collecting acoustic emission data through the acoustic emission probe and deformation data on the surface of the cylindrical specimen through the deformation sensor while the compression test is being carried out; with the deformation data comprising axial deformation and circumferential deformation.
S3: calculating the deformation data collected in S2 according to a finite deformation theory, and obtaining a parameter—mean rotation angle θ, which characterizes a macroscopic deformation characteristic of materials at each stress level and is specifically:
F
j
i
=S
j
i
+R
j
i (1),
where Fji is a deformation gradient, orthogonal transformation Rji is a rotation tensor while symmetric transformation Sji is a strain tensor, and an expression of Sji is as follows:
during a test measurement, calculating out a strain component based on a small deformation theory below:
where ui|j is a covariant derivative of displacement and εji is a small deformation strain;
combined a small deformation strain component with a finite deformation strain component, getting:
S
j
i=εji−(1−cos θ)LkiLji (4),
where Ljk, is an azimuth tensor of a rotation axis;
according to Hooke's law, one-dimensional elastic lossless constitutive formula is
σ=ES (5),
from formula (4) and formula (5), getting
σ=Eεji−E(1−cos θ)LkiLjk (6),
where σ is a stress;
extending formula (6) to a three-dimensional state and writing formula (6) as
in a triaxial test, σ22=σ33=σcon, combined with formula (10), getting:
in a triaxial compression test, there being an assumption as follows:
(L21)2=(L32)2=(L13)2 (9),
writing formula (8) as
obtaining a formula for calculating the mean rotation angle θ from formula (12):
so as to calculate the mean rotation angle θ;
S4: using Grassberger-Procaccia (G-P) algorithm on the acoustic emission data collected in S2 to calculate a fractal dimension of a temporal distribution DT of an acoustic emission signal and calculating a fractal dimension of a spatial distribution DS according to a spatial projection method; specifically:
taking time series of the acoustic emission signal as a research object, then corresponding each time series to a series set with a capacity of n:
X={x
1
,x
2
, . . . ,x
n} (12),
constructing formula (12) as a m-dimensional phase space (m<n), firstly, taking m numbers as a vector of m-dimensional space
X
1
={x
1
,x
2
,x
3
, . . . ,x
m} (13),
then shifting one data to the right and taking m numbers again to form another vector, and so on to form N=n−m+1 vectors; the corresponding correlation function is:
where H is Heaviside function; r is a given scale; when assigning a value to scale r, making r=kr0 in order to avoid dispersion, where k is taken as a scale coefficient and
obtaining n points in a double logarithmic coordinate system, and performing data fitting on n points. If the result is a straight line, it shows that the acoustic emission series has fractal characteristics in a given scale range, and a slope of the straight line is the fractal dimension of the temporal distribution DT of the acoustic emission parameter, namely
D
T=1 gW(r)/1 g(r) (15).
For DS, the space box dimension is used. The box dimension is defined as:
N(r)=Cr−D
where N(r) is the number of discrete bodies whose characteristic size is greater than r, C is a material constant, and the other form of the above formula is the number-radius relation as follows:
M(r)=Cr−D
where r is different radii covering natural discrete bodies, and M(r) is the number of discrete bodies covered in a circle with radius of r. Taking the logarithm on both sides to get:
1 gM(r)=1 gC+DS1 g(r) (18),
where DS is the fractal dimension of the spatial distribution;
S5: carrying out a scanning electron microscope (SEM) test on a fracture surface of the sample after the compression test is completed, to obtain a microscopic morphology of the fracture surface, observing the morphology of the fracture surface and calculating the fractal dimension DA of the fracture surface.
The number of units needed to cover an image in units of δ is N(δ), DA=−log(N(δ))/log δ.
S6: because of a correspondence of a change of the mean rotation angle θ to each process of rock deformation, including a compaction stage, a linear stage and a plastic yield stage in a compression process, finally obtaining a mathematical trend relationship between θ and DT, DS and DA through comprehensively analyzing the obtained the fractal dimension of the temporal distribution DT of the acoustic emission, the fractal dimension of the spatial distribution DS of the acoustic emission and the fractal dimension DA of the fracture surface at each stress level (prior to a peak strength) and the mean rotation angle θ at a same stress level, as shown in a following formula,
θ=a*DT+b*DS+c*DA (19),
eventually obtaining values of a, b and c, so as to establish the quantitative relationship between macro and micro in the whole process of rock deformation and failure.
In this application, specific examples are used to explain the principle and implementation of this application. The explanations of the above examples are only used to help understand the methods and core ideas of this application. At the same time, according to the ideas in this application, there will be some changes in the specific implementation and application scope for ordinary technicians in this field. To sum up, the contents of this specification should not be construed as a limitation to this application.
Number | Date | Country | Kind |
---|---|---|---|
202110549984.9 | May 2021 | CN | national |