Research Article  Open Access
Hansheng Geng, Hongfa Xu, Lei Gao, Ling Chen, Bo Wang, Xiao Zhao, "Elastic Modulus and Strength of RockLike Material with LockedIn Stress", Mathematical Problems in Engineering, vol. 2018, Article ID 5320906, 17 pages, 2018. https://doi.org/10.1155/2018/5320906
Elastic Modulus and Strength of RockLike Material with LockedIn Stress
Abstract
Lockedin stress is a special stress objectively existing in the rock which has an important influence on the mechanical properties. However, there were few studies on the lockedin stress in rocks. The lockedin stress was simplified into spherical inclusions of stress. Based on the Eshelby inclusion theory, the influence equation of the lockedin stress on the elastic modulus of the rock is derived. Through experiment, the expression of elastic modulus of rocklike material with lockedin stress based on MT method is determined. Based on MohrCoulomb criterion and effective stress principle, the equations for the influence of lockedin stress on strength of rocklike material under uniaxial and conventional triaxial conditions were deduced. Through experimental verification, an empirical formula for the strength of rocklike material containing lockedin stress is obtained. This article provided a method for further studying of the lockedin stress in rocklike material and has certain reference value.
1. Introduction
Rocks with randomly distributed inclusions can be considered as natural composite material. The main component of rock can be regarded as the matrix of composite material; cracks, heterogeneous particles, and pores in the rock can be regarded as inclusions in composite material. The elastic properties of composite materials are important issues and many theoretical results have been achieved so far. For engineering, rock is generally in elastic deformation stage before yielding. In this article, the micromechanics theory was used to study the elastic properties of rock containing inclusions. So far, many achievements have been made in the study of lockedin stresses in metals, for example, Voyiadjis’s study on the determination of lockedin stresses in metals (Voyiadjis [1, 2]), but there is little research on lockedin stresses in rocks. Therefore, this article will further study the situation of inclusion with lockedin stresses.
The strength criterion of rock materials is a macroscopic concept. However, microinclusions can significantly affect the strength of the rock material, which can be derived from experiments. In fact, the influence of the mesostructure on failure criterion not only exists in the study of lockedin stress, but also exists in other respects. For example, Chen [3] studied a new equivalent statistical damage constitutive model on rock block mixed up with fluid inclusions. Chen [4] studied the influence of porosity on the strength of sandstone based on the MohrCoulomb strength criterion. Wang [5] studied the strength and constitutive relationship of rock under the joint action of pore water pressure and confining pressure based on statistical theory. Lai [6] and Li [7] studied the elastic properties and constitutive relations of frozen soils at different temperatures. Zhou [8] studied the mechanical properties of rocks under the combined action of temperature, pore water, and confining pressure. Ren [9] researched threedimensional mesodamage model of rock concrete based on the deformation and expansion of elliptic microcrack. Yang [10] studied the joint action mechanism of cyclic loading and pore water pressure on rocks. This paper intends to establish the criteria for the destruction of rock materials affected by inclusion content and lockedin stress. The MohrCoulomb criterion (Jaeger [11]) is a commonly used failure criterion for rock materials. It can describe the failure of rock under uniaxial compression and conventional triaxial compression. The cohesion force of rocklike test blocks used in this paper is provided by cement; the shear friction is provided by granular fine sand. The failure characteristics of samples are in accordance with MohrCoulomb strength theory.
In this article, the inclusion mechanics of composite materials and the MohrCoulomb strength criterion were used to analyze the mechanism of the variation of elastic modulus and peak strength of rock specimens with inclusions content and value of lockedin stress.
2. Mechanical Model of Rock with Lockedin Stress
Lockedin stress in rocks actually exists in various forms and is inhomogeneous, which is determined by the nature of geological evolution. Carrier materials of lockedin stress (inclusions) are of various shapes and sizes. The material outside the lockedin stress carrier (matrix) also has random diversity. Randomness is the difficulty in quantifying mechanical parameters of rock, as well as in the study of lockedin stress. It is only possible to obtain a simple analytical solution by simplifying the problem into ideal inclusion with limited content in the elastic matrix. With reference to the practice in the mesomechanics of composite materials, one cell element in the rock is analyzed. Let the volume of the cell be 1; the inclusion volume is the volume content of the inclusion in the test block; is the lockedin stress when no external load is applied (Figure 1). The problem of rock contain lockedin stress is reduced to the problem that the ellipsoidal heterogeneous inclusions contain initial intrinsic strain in the homogeneous matrix, and, from this, we get the relevant conclusions about the coupling of external loads and lockedin stress in ellipsoidal heterogeneous inclusions.
, , and are the three principal stresses of the external load; is the lockedin stress caused by inconsistent deformation of the inclusion and the matrix (Figure 1). There are several reasons for inconsistent deformation, such as nonrecoverable plasticity, slippage, nonuniform thermal expansion, weathering, and martensite deformation. Figure 2 can simply state the process of the generation of lockedin stress. There was a heterogeneous sphere in the isotropic material; initially, the sphere and the external medium were in close contact with each other, and the lockedin stress is zero. The heterogeneous sphere was taken out alone and expanded by for some reason (thermal expansion, etc.). Only one cavity remains in the matrix material, from which the heterogeneous spheres have been removed. Due to the different material properties, the change in the matrix under the same conditions was negligible and so did the radius of the spherical cavity in the base. At this time, a normal force was applied to the outside of the heterogeneous sphere to restore it to its original shape and reimposed it into the cavity. Then the normal applied stress was removed and the heterogeneous sphere expands, causing compression between the sphere and the matrix. Due to the constraints of the external matrix, the amount of springback of the heterogeneous inclusion cannot reach , but there is a difference from . It can be seen, the cavity radius of the base material is expanded by ; and the compression of the heterogeneous sphere relative to free expansion is . Therefore, the value of the lockedin stress can be obtained according to the size change and the mechanical parameters of each material (Figure 2).
The problem studied in this article is that the rock contains lockedin stress and is subject to external loads; i.e., the initial state is the selfbalanced state in Figure 2. An external load is applied on the basis of Figure 2, as shown in Figure 1.
3. Changes of Lockedin Stress with External Loads
3.1. Tensor Solution Based on Inclusion Theory
The lockedin stress in the inclusion is not fixed. When the rock is subjected to an external load, the deformation of the matrix and the inclusion must be linked to each other. When the external load is tensile, it will reduce the constraints of the matrix on the inclusion. If the adhesion between the inclusion body and the base body is weak, the interface will be separated and the medium will not be continuous. In this study, the applied loads are all compressive stresses. In the elastic stage, the compression of the matrix material will transfer the external load to the inclusion, so that the lockedin stress in the inclusion will increase, and react on the matrix.
Based on inclusion theory (Eshelby [12]) and MT method (Mori [13]), combined with the simplified mechanical model shown in Figure 1, the lockedin stress in the inclusions after external load was applied on the rock was analyzed.
The problem in Figure 1 is a heterogeneous inclusion with intrinsic strain. The cells are affected by the uniform stress in the far field. When there is no heterogeneous inclusion in the cell, the resulting strain is . The relationship between stress and strain is
is the elastic constant tensor of the matrix. Due to the presence of inclusions, the actual cell strain is in fact different from and will produce a perturbation strain . The average stress of the matrix part in this cell is
Obviously, the perturbation part of the stress in the matrix is
Due to the different elastic properties of the material, the average stress and strain in the inclusions are different from the corresponding values in the matrix, and their differences are and , respectively. Therefore the stress in the inclusion is
If there is a lockedin stress in the inclusion, i.e., the intrinsic strain , then there are
Subtracting the right sides of (5) and (4) yields
It can be seen that, for a given inclusion body’s intrinsic strain (or lockedin stress), the amount of change in external stress to the inclusion body’s stress is irrelevant to the original lockedin stress. Therefore, only the situation of (4) needs to be analyzed.
is the elastic constant tensor of the inclusion and is the equivalent intrinsic strain of the inclusion. According to the Eshelby inclusion theory there are
From (7), (4), and (2) we can find
The volumeaverage stress of the cell should be equal to the average stress of the farfield effect and is
According to (2) and (4) we get
Substituting (7) and (11) into (4), the solution is
Substituting (12), (11), (1), and (7) into (4) can solve the relationship between the stress in the inclusion and the far field stress of the cell.with
Equation (13) is the general solution of the lockedin stress in the inclusion with the change of the external load of the test block.
3.2. The Solution When the Inclusion Is Spherical
Isotropic fourthorder tensor in composite micromechanicsis abbreviated as
, , and in accordance with the algorithm, , , , fourthorder unit tensor can be expressed as and the elastic constant tensor can be expressed as .
The equation of the tensor can be obtained by the following calculation process:with (Shen [14]),is the bulk modulus and shear modulus of the matrix, andis the bulk modulus and shear modulus of the inclusion (Xu [15]).
LetEquation (25) is obtained by (19):
3.3. Examples
(1) The matrix material is made to have an elastic modulus of 1MPa and a Poisson’s ratio of 0.3. The inclusion material has a Poisson’s ratio of 0.3, an elastic modulus of 0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, and 5, and an external load applied to the rock is a uniaxial pressure of 6 MPa. The variation of the component of the lockedin stress parallel to the loading surface with the content of the inclusion and the ratio of the elastic modulus of the inclusion to the matrix can be obtained by (25), as shown in Figure 3(a). As the elastic modulus of the inclusion changes from small to large, the lockedin stress generated in the inclusion is also increasing.
(a) Hydrostatic pressure 6MPa
(b) Uniaxial pressure 6MPa
(2) The external load applied to the rock is a hydrostatic pressure of 6MPa, and the other conditions are the same as (1), as shown in Figure 3(b). As the elastic modulus of the inclusion changes from small to large, the lockedin stress generated in the inclusion is also increasing.
(3) The matrix material is made to have an elastic modulus of 1MPa and a Poisson’s ratio of 0.3. The inclusion material has a modulus of elasticity of 1 and a Poisson’s ratio of 0.1, 0.2, 0.3, 0.4, and 0.49, and an external load applied to the rock is a uniaxial pressure of 6MPa. As the Poisson’s ratio of the inclusion changes from small to large, the lockedin stress generated in the inclusion is also increasing.
(4) The external load applied to the rock is a hydrostatic pressure of 6MPa, and the other conditions are the same as (3), as shown in Figure 4(b). As the Poisson’s ratio of the inclusion changes from small to large, the lockedin stress generated in the inclusion is also increasing.
(a) Hydrostatic pressure 6MPa
(b) Uniaxial pressure 6MPa
When the elastic parameters of the inclusions are the same as the matrix, This model degenerates into a homogeneous medium, and the lockedin stress generated in the inclusions is independent of the inclusion content. As the inclusion content increases, the influence of elastic parameters becomes smaller and smaller.
Consider the case of uniaxial compression; there is at this time. Test data of Geng [16] is cited here. The matrix material has an elastic modulus of 2560 MPa and a Poisson’s ratio of 0.3; the elastic modulus of inclusion is , and the Poisson’s ratio is 0.49. The relation between and can be obtained from (27). Figure 5 shows the increase of the lockedin stress when the uniaxial pressure is 2MPa, 4MPa, 6MPa, and 8MPa, respectively, and it is related to the initial lockedin stress and inclusion content.
Secondly, consider the case of hydrostatic compression; i.e. . Figure 6 shows the increase of the lockedin stress when the hydrostatic pressure is 2 MPa, 4 MPa, and 6 MPa, respectively, and it is related to the initial lockedin stress and inclusion content.
From Figures 5 and 6 it can be seen that the external load will generate additional stress on the inclusions with lockedin stress, which is related to the material properties and external forces. The effect of the lockedin stress shown in the figure on the results is due to the fact that the lockedin stress of the inclusion in the object studied here is related to the temperature, and the volumetric elasticity of the inclusion is also related to the temperature, which is essentially the effect of temperature. There is no connection with the lockedin stress. When the lockedin stress is small, the increase of it is not negligible. The increase is proportional to the increase in uniaxial compression load and hydrostatic pressure load.
The derivation of (1), (2), (3), (4), (5), (7), (8), (9), (10), (11), (12), (13), (14), (15), (16), (17), (18), (19), (23), (24), and (25) is also applicable to other forms of inclusion material. In fact, this increase in the lockedin stress is different from the lockedin stress but is the additional stress in the loading. After unloading this increment disappears but the lockedin stress does not.
4. Theoretical Study on the Change of Elastic Modulus of Rocks with Lockedin Stress
4.1. Solution without Considering LockedIn Stress
(1) Micromechanical Analysis of Equivalent Elastic Modulus. The upper limit of inclusion content is 1 and the lower limit is 0. When the content of the inclusion is 0, it is obvious that the elastic modulus of the rock is the elastic modulus of the intact rock. When the content reaches 1, it becomes a sample composed of the inclusion material, and its elastic modulus is consistent with the inclusion material. Therefore, the relationship between the elastic modulus of inclusioncontaining specimens and the variation of inclusion content should satisfy these two prerequisites.
First, the relationship between elastic modulus and content is analyzed. When calculating the elastic modulus, both the matrix and inclusions are considered to be in a linear elastic state. Therefore, the problem of the elastic modulus in this paper is simplified to discuss the problem of spherical inclusions in elastic materials. For the elastic solution of inclusioncontaining materials, there are currently Voigt upper limits and Reuss lower limits:
Taylor method is
MT method iswhere is the macroscopically equivalent elastic tensor of the rock containing the inclusion, is the elastic tensor of the rock, is the elastic tensor of the rphase inclusions.
In addition, there are differential method, generalized selfconsistent method, and Hashin and Shtrikman’s upper and lower limit methods. For cases with a single inclusion, (27) to (28) can be written in the following scalar form.withThe equivalent elastic modulus can be obtained according to
The matrix material is made to have an elastic modulus of 1MPa and a Poisson’s ratio of 0.3; the inclusion material is made to have a Poisson’s ratio of 0.3 and an elastic modulus of 0.2, 0.6. The variation of the equivalent elastic modulus with the content of inclusions and the ratio of the elastic modulus of the inclusions to the matrix can be obtained by (26), (27), (28), (29), (30), and (32) (Figure 7(a)). Let the matrix material have an elastic modulus of 1MPa and a Poisson’s ratio of 0.3; let Poisson’s ratio of the inclusion be 0.3 and the elastic modulus of the inclusion be 2 and 3. The variation of the equivalent elastic modulus with the inclusion content and the ratio of the elastic modulus of inclusion to the matrix can be obtained by (26), (27), (28), (29), (30), and (32) (Figure 7(b)).
(a) Flexible inclusion
(b) Hard inclusions
When the elastic modulus of inclusion is larger than that of the matrix, the macroelastic modulus increases with the increase of the inclusion content, and the elastic modulus of inclusion increases, when the elastic modulus of inclusion is smaller than that of the matrix. The macroelastic modulus decreases with the increase of inclusion content and increases with the increase of inclusion elastic modulus. Under the same conditions, the value obtained by the Voigt method is the highest, the value obtained by the Reuss method is the lowest, and the MT method curve is slightly higher than the Taylor method curve. In the case where the elastic modulus of the inclusion is small enough, the curve of the Taylor method will be distorted (Figure 7).
There is no Poisson’s ratio in the expressions of elastic modulus of Voigt upper limit method and Reuss lower limit method, so the change of Poisson’s ratio will not affect the macroscopic elastic modulus, so it will not be discussed here.
Secondly, we discuss the influence of the Poisson’s ratio change on the elastic modulus when the Taylor method and the MT method are used. Let the elastic modulus and Poisson ratio of the matrix be 1 and 0.3, the elastic modulus of inclusion be 2 and 0.2, and Poisson’s ratio be 0.2, 0.4. Figure 8 shows that as Poisson’s ratio increases, the slope of the curve becomes larger. The macroelastic modulus at the same inclusion content increases. Under the same conditions, curve of the MT method is slightly higher than that of Taylor method, but they are very closed.
(2) Test Verification. Using rocklike materials and rubber particles to make a test block rocklike model containing inclusions, Geng [16] studied the influence of lockedin stress on the mechanical properties of rocks. By heating the test piece, the rubber particles expand, creating a closed stress in the inclusion. Uniaxial compression tests of similar model blocks with different lockedin stresses and different inclusion contents were performed with a water bath and MTS5T uniaxial compression testing machine. Results show that the peak strength and elastic modulus decrease with the increase of lockedin stress and decrease with the increase of the inclusion content. With the decrease of the inclusion content, the effect of the closed stress on the peak strength and elastic modulus gradually decreases. When the content of the inclusion is zero, the influence of the closure stress also disappears. Through the mathematical fitting analysis, the quantitative formula for calculating the relationship between the peak strength, the elastic modulus of the test piece, the lockedin stress, and the inclusion content was obtained. The expression and value curve of joint impact factor are given. Results of this study reveal the importance of the lockedin stress in the mechanical properties of the rock and provide a reference for the further study of the lockedin stress in rocks.
Substituting the mechanical parameters of the matrix material and the inclusion material in the test of Geng [16] into (26), (27), and (28), the resulting equation is plotted as a curve (Figure 9).
Figure 9 depicted the elastic moduluscontent relationship obtained by different methods when the lockedin stress is 0. It can be seen that both the MT method and the Taylor method can well match the experimental data in the 00.04 inclusion stage.
The existing research results suggest that the Taylor method is only applicable to the case where the inclusion content is less than 5% (Du [17]). This is because when the content of inclusions exceeds 5%, the interaction of inclusions will become unnegligible.
There is small but negligible deviation between calculated results and the experimental results. Equations (29) and (30) are similar in structure, where the denominator has one more parameter related to inclusion content. In general, the MT method is more suitable for this article.
In summary, when there is no lockedin stress in the inclusion, the scalar form of the macroelastic modulus of the sample can be obtained by substituting (30) into (32):
4.2. Solution with Considering LockedIn Stress
The relationship between the equivalent elastic modulus and the lockedin stress can not be found by the equation derived from Eshelby and the MT method. The change in the volumetric content of inclusions due to intrinsic strain is negligible relative to the volume content. In the MT method, the lockedin stress is an “irrelevant quantity”, because both the Eshelby theory and the MT method have a preset condition; namely, it’s continuous between media. At present, the research on the inclusion of stress in rocks mainly focuses on the inclusion of fluids, such as pore water or pore gas, but such research is rare. Liu [18] studied the effects of pore pressure changes on rock strength characteristics and found that the increase in pore pressure caused a decrease in elastic modulus. As the pore pressure increases, some cracks open and the contact between the particles decreases, resulting in a change in the matrix skeleton unit. Li [19] studied the relationship between the shear modulus of saturated sand and water pressure through a saturated sand liquefaction test under triaxial conditions.
(1) Mesoscopic Analysis of Equivalent Elastic Modulus. The lockedin stress studied in this article is essentially different from the pore pressure caused by gas or fluid which is isotropic. Due to the fluidity, the geometry of the fluid changes with the pore morphology during rock damage, while the solid inclusions remain relatively intact. Therefore, the lockedin stress field of fluid inclusions under external loads is more complicated. But as the lockedin stress increases, the trend of the decrease of the elastic modulus of the test block does not change.
The lockedin stress caused by the expansion of the inclusion body can not only squeeze the crack parallel to the spherical surface of the inclusion body to make it more compact, but also open the crack perpendicular to the spherical surface and increase the damage of the rock. This initial damage makes the mechanical properties of matrix material change, and the shear modulus and bulk modulus have undergone a certain degree of reduction. Assume that the bulk modulus and shear modulus of the matrix after reduction are and , respectively. Therefore, the macroelastic modulus of the test block iswith
Therefore, the change of the mechanical properties of the sample is the result of the joint action of the inclusion content and the lockedin stress. When studying the statistical damage constitutive model, Li [20] equated damaged rock materials with pores, and the volume content of the initial damage (i.e., additional pore) caused by the lockedin stress is . In 1961, Greertsma and Smit [17] obtained an empirical formula for the relationship between rock bulk modulus and porosity based on experimental results. This paper assumes that this empirical relationship also applies to the relationship between shear modulus and porosity, yielding
and is a function of the lockedin stress and the inclusion content . The equation is
Since a and b are undetermined coefficients, (37) and (38) can be abbreviated as
Equation (39) is the problem of expansion damage of spherical inclusions. This formula is discussed in detail below.
(2) Damage Mechanical Analysis of Spherical Inclusions. The initial damage is caused by the squeezing of the lockedin stress. The expansion of the inclusion, in addition to the strain energy created by the expansion of the inclusion in itself, also performs work on the base material. The volume of compression of the inclusion is , and the amount of expansion of the matrix body is . is the intrinsic radius increment of the inclusion (Figure 2). According to the expression in Geng [16] regarding lockedin stress,Equations (43) and (44) are obtained:
When the lockedin stress is caused by the thermal expansion of inclusions, can be expressed by (39):
is the increase in radius due to free thermal expansion when the rubber ball is unconstrained. The thermal expansion of the matrix is neglected here since the coefficient of it is extremely small relative to the rubber ball. According to (42), (43), and (39), it can be obtained that, under the constraint of the matrix, the sphere is compressed by
From (44), the radius of the inner wall of the matrix of is expanded by
Therefore, the internal energy generated in the sphere during the generation of the lockedin stress is
The work done by the expansion of individual rubber particles on the matrix is
The ratio of and can be obtained from (40) and (47):
Substituting (50) into (49) yields
The total rubber content is , and the volume of the test block is 1. The total work done by all rubbers in the test block on the matrix is
According to the uniaxial compression test results of the matrix material, the required strain energy for destruction of matrix material iswhere
is the strain energy required for the matrix material destruction, and the damage degree of the matrix material under uniaxial compression failure is 1. The damage degree of the matrix material when no external load is applied is 0. Let the relationship between strain energy and damage degree follow power function distribution:where and are distribution parameters that can be combined with the test.
indicate the degree of damage:
Combining (52) and (53) yields
Substituting (56) into (40) and (41) yields
Substituting (57) and (58) into (32) yields
(3) Test Verification. The test data of Geng [16] is cited here; the lockedin stress is caused by the uncoordinated thermal expansion of the matrix and inclusions. Substituting (57) and (58) into (59) yields an expression of the elastic modulus of the test block when the inclusion contains lockedin stress. Then, the threedimensional nonlinear fitting of the elastic modulus test data of the Geng [16] by (61) results in , , (). The fitted result is shown in Figure 10. When in (57) and (58) is 0, (59) degenerates to (33).
5. Strength Change Mechanism of Test Block
5.1. Theoretical Analysis of Peak Strength without Considering LockedIn Stress
The failure mode of the sample is shear failure. Therefore, the MohrCoulomb failure criterion is used (Figure 11), in the form ofwhere is the shear strength, is the normal stress of the shear plane, is the cohesion of the sample, and is the internal friction angle of the sample.
For the case of uniaxial compressive stress , the specimen is broken along the surface at a angle of to the crosssection and can be expressed as
In the case of softer inclusions, the shear modulus of inclusion is much smaller than that of the matrix, and basically no shear stress is assumed in the sample. For normal stress, the presence of inclusions reduces the area of the matrix on the shear plane.
Due to the small deformation of the specimen, the inclusion can be considered to be always spherical with a radius of r and a content of . Take a representative cell with a volume of 1 containing a single inclusion with a volume of as an example at the cut section, and the six faces of the cell are parallel or perpendicular to the shear plane (Figure 12). Therefore, according to the volume formula of the sphere, there is a radius of the inclusion. The part of the sheared broken section passing through the cell is a square region with an area of 1. The section of the sheared broken section passing the inclusion is a circle and its area is
(a) Force model of the specimen
(b) Representative unit body
(c) Representative unit surface
The area of the matrix material part on the broken surface is 1. Assuming that the elastic modulus and shear modulus of the inclusions are very small compared to rocklike materials, it can be considered that there is basically no ability to carry compressive and shear stresses, and the bearing capacity is close to that of the holes compared to the matrix made of rocklike materials. Therefore, the average normal stress and the average shear stress of the matrix part on the shear surface are not equal to and . According to Lemaitre’s strain equivalence assumption (Lemaitre [21]), the true stresses on the shear plane are the effective normal stress and the effective shear stress , and the corresponding and are called the nominal normal stress and the nominal shear stress here.
Figure 13 shows the MohrCoulomb model based on effective normal stress and effective shear stress. It can be seen that, as the content of inclusions increases, the nominal stress required for the shear plane to reach failure is getting smaller and smaller. According to the assumption of strain equivalence, the effective normal stress and the effective shear stress on the failure surface of the sample containing the inclusion are
According to the theory of elasticity, the relationship between the axial stress , , and of the specimen before the specimen is destroyed always conforms to (61) and (62). Substituting (61) and (62) into (64) and (65) yields
From (66), (67), and (60), is obtained:
Taking the assumptions of the distribution of the inclusions into account, the above formula is modified. Parameter is introduced in (68) and results in
Equation (69) needs to satisfy two conditions: (1) when , it degenerates to the strength of the matrix; (2) when , there is no matrix present on the failure surface and the intensity is zero. That is, in (69), needs to satisfy the following: (1) when , ; (2) when , .
Let ; it was found that when meets the following formula, it can meet the above two conditions:
From formula (69) and formula (70), is obtained:where parameter can be obtained in combination with the uniaxial compression test data.
5.2. Theoretical Analysis of Peak Strength Considering LockedIn Stress
Figure 14 is a schematic diagram of the force of a representative unit on the failure surface of the sample when the inclusion contains a lockedin stress. When the initial lockedin stress of the inclusion is not 0, the influence of the stress concentration effect of the inclusion must be taken into account as well as the influence of the initial lockedin stress. The effect of the initial lockedin stress on the failure shear plane of the specimen is to cancel out a part of the normal stress on the shear plane and reduce the maximum shear stress required for the failure. Let the normal stress and shear stress on the failure surface be and :
Substituting (72) and (73) into (60) yields
Considering the stress concentration effect around the inclusion, the lockedin stress term of (74) is modified and the correction coefficient is introduced, which is considered to be related to the inclusion content:
Equation (74) is amended to
Combined with the correction factor related to the distribution of inclusions in (71), we can get
Simplify (77) and obtainwhere is the parameter determined by the uniaxial compression test of rocklike blocks without lockedin stress and is the parameter related to inclusion content determined by uniaxial compression test of rocklike blocks with lockedin stress.
5.3. Theoretical Analysis of Peak Strength Considering Confining Pressure
When considering confining pressure, there are
Substituting (79) and (80) into (60) yields
Equation (81) can be converted into
By observing (82), it is found that when the other parameters are constant, the confining pressure is proportional to the change of the deviatoric stress intensity, which is because that the equation deduced based on the MohrCoulomb equation (60) was actually a linear equation, which does not match the actual situation. The relationship between confining pressure and peak deviatoric stress is often a nonlinear relationship in the parabolic or hyperbolic form. Therefore, it is necessary to adjust the confining pressure item in (82). Parameter is introduced combined with parameters in (78) to get the final expression of peak deviatoric stress:where , , and are parameters related to confining pressure and can be obtained by fitting the test data.
When the confining pressure is 0, (83) degenerates to (77).
5.4. Test Verification
(1) Uniaxial Compression without Considering LockedIn Stress. Using (69) to fit the test data of Geng [16] when there is no lockedin stress, can be obtained, as shown in Figure 15.
So when the inclusion does not contain a lockedin stress, the expression of the strength of the test block is
(2) Uniaxial Compression with Considering LockedIn Stress. Figure 16 shows the uniaxial and conventional triaxial compressive stressstrain curves of a standard sample of cement mortar without inclusions obtained by laboratory tests. The confining pressure and its corresponding peak strength in the figure are listed in Table 1.

(a) The compression curve of the matrix under different confining pressure
(b) The relationship between matrix material strength and confining pressure
Using the test data in Table 1, the cohesion c and the friction angle of the cement mortar matrix can be obtained by linear fitting according to MohrCoulomb equation (85). The fitting curve is shown in Figure 16(b) in the solid line.
The cohesive force of the matrix material is 1MPa; the friction angle is by fitting.
Equation (77) is used to fit the strength test data when there is lockedin stress; it is found that when meets , (77) can be consistent with the experimental data, and the fitting result is shown in Figure 17, with , ().
Therefore, the strength expression of the specimen containing lockedin stress inclusions is
When in (86) is 0, it degenerates to (69).
It can be seen that the inclusion content and the value of lockedin stress are significant for the change of the peak stress and elastic modulus of rocklike materials, and with the increase of lockedin stress, the effect on elastic modulus and peak stress will also grow bigger.
(3) Conventional Triaxial Compression Considering LockedIn Stress. From (83), firstly determine the parameter based on the experimental values of the pure matrix material under different confining pressures. Then on the basis of the test values without considering the lockedin stress under each confining pressure, the parameter is determined by fitting. On this basis, the parameters , are determined based on all the data under the corresponding confining pressure. See Table 2 for various parameters under different confining pressures. The result of the fitting is shown in Figure 18.

(a) The confining pressure is 2MPa
(b) The confining pressure is 4MPa
(c) The confining pressure is 6MPa
From Table 2, it can be seen that the parameters , , , of (83) change with the confining pressure. Among them, decreases with the increase of confining pressure; increases with the increase of confining pressure; decreases with the increase of confining pressure; decreases with the increase of confining pressure and then increases. The fitting results of these parameters with the confining pressure are shown in (87). The fitting curve is shown in Figure 19.
(a)
(b)
(c)
(d)
Substituting the value in (87) into (83) yields the strength expression of the specimen containing the lockedin stress when considering confining pressure.
6. Conclusion
The change mechanism of the elastic modulus and peak strength of the sample along with the variation of the inclusion content and the lockedin stress were analyzed, and the following results were obtained.
(1) According to Eshelby’s equivalent inclusion theory, the variation law of lockedin stress with external load in rock was deduced. Results showed that, with the increase of the external load, the lockedin stress would also increase; the increase range was related to the ratio of the elastic modulus of the inclusion to the matrix, the ratio of Poisson’s ratio of the inclusion and the matrix, and the content of inclusions. When the inclusions and the matrix had the same parameters except the elastic modulus, if the ratio of the elastic modulus of the inclusions to the matrix is greater than 1, the lockedin stress decreases with the increase of the content; while the ratio of the elastic modulus of the inclusion to the matrix is less than 1, the lockedin stress increases with the increase of the content.
(2) Composite micromechanics theory was introduced. The relationship between the elastic modulus of the rock and the content of inclusions without considering the lockedin stress was obtained, based on Voigt upper limit and Reuss lower limit method, Taylor method, and MT method, respectively. Through verification and comparison with the previous test results, the MT method was used to describe the equivalent elastic modulus of rocklike material containing inclusions.
(3) The concept of strain energy and damage was introduced, and the relationship between strain energy and damage degree were established through a statistical probability model. Based on the MT method, the theoretical solution of the equivalent elastic modulus of rock with lockedin stress is deduced. Through the threedimensional nonlinear fitting, the empirical formula of the change law of the elastic modulus of the sample with the content of the inclusion and the lockedin stress value is obtained, which is in good agreement with the test result.
(4) MohrCoulomb strength theory was introduced, and the relationship between the peak stress and the inclusion content of the test block without considering the lockedin stress was obtained by combining the effective stress principle. The empirical formula of the rule that the peak stress of the test block containing lockedin stress inclusions varies with the inclusion content and the lockedin stress value was obtained by the threedimensional nonlinear fitting method.
(5) The concept of effective stress was extended to the conventional threeaxis case. Based on the MohrCoulomb model, the empirical formula of the variation of the peak deviatoric stress of the test block containing lockedin stress under different confining pressures with the inclusion body content and the lockedin stress value is obtained. The parameters corresponding to the empirical formulas for different confining pressures were obtained by the test data. Through polynomial fitting, the relationship between these parameters and confining pressure was obtained.
The spherical rubber particles used to simulate stress inclusions can only represent a certain type of stress inclusions. There are many types of stress inclusions in rocks. To conduct systematic research, more advanced instruments are needed to conduct a large number of tests.
Data Availability
Previously reported experiment data were used to support this study. The prior study is cited at relevant places within the text as [16].
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors would like to thank Key Laboratory of Ministry of Education for Geomechanics and Geomechanics and Embankment Engineering, Hohai University (no. 201707) and Natural Science Foundation of Jiangsu Province of China (BK20141067) for their financial support.
References
 G. Z. Voyiadjis, P. D. Kiousis, and C. S. Hartley, “Analysis of residual stresses in cylindrically anisotropic materials,” Experimental Mechanics, vol. 25, no. 2, pp. 145–147, 1985. View at: Publisher Site  Google Scholar
 G. Z. Voyiadjis and C. S. Hartley, “Residualstress determination of concentric layers of cylindrically orthotropic materials,” Experimental Mechanics, vol. 27, no. 3, pp. 290–297, 1987. View at: Publisher Site  Google Scholar
 X. Chen, H. Xu, H. Geng, L. Dong, and J. Zhang, “A New Equivalent Statistical Damage Constitutive Model on Rock Block Mixed Up with Fluid Inclusions,” Mathematical Problems in Engineering, vol. 2018, Article ID 3080173, 11 pages, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 Y. M. Chen, “Study on damage mechanics of effect of porosity on sandstone strength,” Sichuan Hydropower, no. 2, pp. 39–43, 1997. View at: Google Scholar
 W. Wang, Z. Y. Tain, and Q. Z. Zhu, “Study of statistical damage constitutive model for rock considering pore water pressure,” Chinese Journal of Rock Mechanics and Engneering, no. s2, Article ID 10.13722/j.cnki.jrme.2014.1293, pp. 3676–3682, 2015. View at: Google Scholar
 Y. Lai, S. Li, and J. Qi, “Strength distributions of warm frozen clay and its stochastic damage constitutive model,” Cold Regions Science & Technology, vol. 53, no. 2, pp. 200–215, 2008. View at: Publisher Site  Google Scholar
 S. Li, Y. Lai, and S. Zhang, “An improved statistical damage constitutive model for warm frozen clay based on Mohr–Coulomb criterion,” Cold Regions Science & Technology, vol. 57, no. No, pp. 154–159, 2009. View at: Publisher Site  Google Scholar
 Q. C. Zhou, “Study on the mechanical property of a sandstone under geothermalmechanical and hydraulicmechanical coupling,” Chinese academy of science (Institute of rock & soil mechanics), vol. 2006, 2006. View at: Google Scholar
 Z. J. Ren, A Micromechanical damage model for rocks and concretes based on the deformation and growth of elliptic microcracks, Chongqing University, 2008.
 H. W. Yang, Study on coupling mechanism of rock and pore water under cyclic loading, Chongqing University, 2011.
 J. C. Jaeger, N. G. W. Cook, and R. W. Zimmerman, Fundamentals of rock mechanics[M]., , Maklen, Blackwell Publishing, Maklen, 4th edition, 2007.
 J. D. Eshelby, “The determination of the elastic field of an ellipsoidal inclusion, and related problems,” Proceedings of the Royal Society London A: Mathematical, Physical and Engineering Sciences, vol. 241, pp. 376–396, 1957. View at: Publisher Site  Google Scholar  MathSciNet
 T. Mori and K. Tanaka, “Average stress in matrix and average elastic energy of materials with misfitting inclusions,” Acta Metallurgica et Materialia, vol. 21, no. 5, pp. 571–574, 1973. View at: Publisher Site  Google Scholar
 G. L. Shen, G. K. Hu, and B. Liu, Mechanics of Composite Materials, Tsinghua University Press, 2013.
 Z. L. Xu, A brief tutorial on elastic mechanics, Higher Education Press, 2002.
 H. S. Geng, Influence of lockedin stress on compressive properties of rock, Army Engineering University of PLA, 2018.
 S. Y. Du, Micromechanics of composite materials, Science Press, 1998.
 X. J. Liu, J. K. Shen, L. X. Liang, and al. et, “Effects of pore pressure on rock strength properties,” Chinese Journal of Rock Mechanics and Engneering, vol. 30, no. s2, pp. 3457–3463, 2011. View at: Google Scholar
 W. Y. Li and H. S. Liu, “Influence of Pore Water Pressure on Shear Modulus and Damping Ratio of Saturated Sands,” Chinese Journal of Geotechnical Engineering, vol. 5, no. 4, pp. 56–67, 1983. View at: Google Scholar
 X. Li, W.G. Cao, and Y.H. Su, “A statistical damage constitutive model for softening behavior of rocks,” Engineering Geology, vol. 143144, pp. 1–17, 2012. View at: Publisher Site  Google Scholar
 J. A. Lemaitre, “How to use damage mechanics,” Nuclear Engineering and Design, vol. 80, no. 2, pp. 233–245, 1984. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Hansheng Geng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.