Researches on Damage Evolution and Acoustic Emission Characteristics of Rocks

Mechanical parameters of the rock are important for the design of geotechnical, mining engineering, and petroleum reservoir projects. Many researches have suggested that the mechanical variables of rock specimens, such as compressive strength and elastic modulus, do not have a single fixed value. Uncertainty in the basic mechanical variables of the rock material can significantly affect the structural performance and safety. In this study, a series of compression experiments with acoustic emission have been performed on rock specimens. (e damage evolution characteristics of the rock in the process of loading were studied, and the macromechanical behaviors were obtained at the same time. Distribution characteristics of the strength and elastic modulus as random variables are illustrated, and the statistical damage model is presented by the authors to formulate analytical constitutive relations for deformation behavior. (e comparisons between predicted results and experimental data show that the statistical damage constitutive model could well reproduce the deformation process of rock materials.


Introduction
Mechanical property of the rock is important for the design of geotechnical, mining engineering, and petroleum reservoir projects.A large number of experiments suggested that the mechanical variables of rock specimens, such as compressive strength and elastic modulus, do not have a single fixed value [1][2][3][4].
e measured mechanical parameters for a given rock usually show a large scatter characteristic.ere is no way of predicting exactly what the values of these parameters are.Uncertainty in the basic mechanical variables of the rock material can significantly affect the structural performance and safety.Understanding this variation, or uncertainty, in the design and analysis of geotechnical engineering requires the use of structural reliability-based design and assessment methodologies [5].Specimens with the same nominal material parameters and tested under the same environmental conditions may exhibit different behaviors with diversified strength and should be analyzed by a distribution function [6].Hence, regarding these mechanical parameters as random variables, the statistical mechanical approach, which has been applied successfully to deal with the mechanical behavior of the rock, becomes a quite attractive tool for investigating deformation processes and failure mechanism of the geomaterial [7][8][9].In fact, the diversity of mechanical behavior for specific samples is usually attributed to the heterogeneity of the material at the mesoscopic level.e most important aspect of the engineering material is its microstructure, which contains many randomly distributed microcracks.e rock material is initially statistically homogeneous but becomes heterogeneous due to the propagation and clustering of microcracks [10].However, the effect of mesoscopic heterogeneity on macroscopic failure is still not clear.Furthermore, considering that mechanical variables of the rock material exhibit great uncertainty, it is unreasonable to take some deterministic values as the mechanic parameters, and the reliability analysis of the mechanical variable for the rock material is necessary.
In this paper, a series of acoustic emission (AE) experiments with uniaxial compression have been performed, and the macromechanic properties and microdamage evolution have been investigated.Distributions of the strength and elastic modulus as random variables are illustrated by several samples of rock specimens, respectively.e normal, logarithmic-normal, and Weibull distributions are used to describe the statistical characteristics.
e damage mechanics is introduced to deal with the deformation characteristic of sandstone.

Experiments
e tests were performed with MTS815, and the acoustic emission testing system was illustrated in Figure 1. e sandstone specimens were prepared with the dimension of the cylindrical specimen 50 × 100 mm.
e AE detection model PCI-2 of the DISP series manufactured by PAC was employed to detect AE signals from the specimens.e AE sensors were attached to the outer wall of the rock specimen for good signal detection.Data logging was controlled by a computer during testing.After the specimen was placed on the xed lower head of the testing machine, the compressive load was applied to the specimens until failure with the strain rate of 1.5 × 10 −5 /s. e upper moving head was brought down until it just touched the specimen.e stress-strain curves of rocks are shown in Figure 2. From Figure 2, it can be seen that the shapes of stress-strain curves are similar in the process of loading.e stress-strain curve can be approximately divided into four typical stages: (1) the compaction stage: the stress-strain curve of the sample shows the downward concave and the initial nonlinear deformation at low stress levels.At the initial process of loading, the pores, voids, and microcracks of the rock sample are compacted with the increase of compression strain; (2) the linear elastic stage: after the initial compaction, stress increases approximately linearly with the further increase of axial strain, and there is little residual strain in this stage; (3) the crack propagation stage: at this stage, the slope of the stress-strain curve gradually decreases with the further increase of axial strain, the evolution of damage during compression loading causes a small amount of residual deformation, and the rigidity of the rock sample reduces; (4) the failure stage: when the microcrack develops to a certain extent, the macroscopic crack in the rock sample comes out rapidly and the stress-strain curves shows a rapid drop in the postpeak region, and the failure behavior of the rock presents an obvious brittle failure characteristic under uniaxial compression conditions.From Figure 2, it can also be seen that the stress-strain curves are signi cantly di erent from each other under the same compression conditions.Specimens with the same nominal material parameters and environmental conditions exhibit di erent behaviors with diversi ed strength and should be analyzed by a distribution function.It is necessary to adopt the probabilistic approach to study the mechanical behaviors, which takes the discreteness of mechanical parameters into account.
Acoustic emission data were analyzed to describe the damage evolution.ese correlations are illustrated in Figure 3. Test results show that the rock emits a small number of AE events before entering the failure stage.AE activity gradually increases with the increase of strain.It is shown that this AE parameter is a reasonable indicator for damage occurring within the rock specimen, except for a short time delay.e stress and strain lag behind AE activity.In the compaction stage, a few events appear suggesting that little damage of the rock happens.
e elastic stage, corresponding to a linear variation of strain with stress, exhibits a slow increase of AE activity.As the strain increases to a critical value, the AE activity increases rapidly [11].e regularities of AE activities are di erent in the loading process for di erent failure types.
e strenuous activity interval of the AE activities is narrow in the failure stage under the condition of tensile failure; for the shear failure, the strenuous activity interval of the AE activities becomes wide.

Probability Distribution Characteristics of the Rock Material
e variations of mechanical parameters of rocks are very large; it is unsafe and unreasonable for engineering design to use the average strength.In this study, the reliability of the mechanical variables, such as strength and elastic modulus, was investigated and discussed.e cumulative probabilistic distribution of data with N number of samples was calculated  2 Advances in Civil Engineering by ranking the data from the lowest value to the highest value using (1).e probability of the specimen is given by [12] where i is the ranking of the data.As it is customary for brittle materials, the strength data are interpreted statistically by assuming that specimen strength obeys a probability density function.In this section, we brie y describe di erent probabilistic models considered here and mention the estimation procedures of the unknown parameters from a given sample dataset.e performance of the models was then compared by matching the calibrated models with the individual test results.

Normal Distribution.
e probability density function of a normal distribution is given as follows [13]: where x is the value of the continuous random variable.e cumulative distribution function for the normal distribution is given as follows: Maximum-likelihood estimation is a method of estimating the parameters of a statistical model.When applied to a dataset and given a statistical model, maximumlikelihood estimation provides estimates for the model's parameters.For normal distribution, the likelihood for the whole sample is given as follows: Taking log, the log-likelihood is given as follows: To nd the critical points of the log-likelihood function, set the rst derivative with respect to each equal to zero.In di erentiating (5), note that

Advances in Civil Engineering
e maximum-likelihood estimates for μ and σ can be found by setting ( 6) and ( 7) equal to zero, and solving for each, μ and σ can be obtained as follows: 3.2.Logarithmic-Normal Distribution.In probability theory, a log-normal distribution is a continuous probability distribution of a random variable whose logarithm is normally distributed.e probability density function of a log-normal distribution is given as follows [14]: e cumulative distribution function for the normal distribution is given as follows: For determining the maximum-likelihood estimators of the log-normal distribution parameters μ and σ, we can use the same procedure as used in the normal distribution.
erefore, using the same indices to denote distributions, we can write the log-likelihood function as follows: Now taking log, the log likelihood is obtained as follows: By differentiating (12), we have Using the formulas of the normal distribution maximum-likelihood parameter method, we deduce that, for the log-normal distribution, it holds that (14)

Weibull Distribution.
e Weibull distribution is one of the most widely used lifetime distributions in reliability engineering.It is a versatile distribution that can take on the characteristics of other types of distributions.e probability density function of the Weibull distribution is given as follows [15]: where k is the shape parameter, λ is the scale parameter of the distribution, and θ is the location parameter of the distribution.e cumulative distribution function for the Weibull distribution is Considering the Weibull distribution given in (15), the likelihood function will be On taking the logarithms of ( 17), we obtain the estimating equation as follows: In differentiating (18) with respect to λ, k, and θ in turn, we have e parameters λ, k, and θ can be accomplished by the use of iterative procedures.Figures 4 and 5 show the theoretical distribution relationships together with the test results for comparison.It can be seen that the normal and 4 Advances in Civil Engineering logarithmic-normal distribution models are very similar, but neither of these models matches with the test results better than the Weibull distribution model.To analyze exactly the precision of the three kinds of probability distributions, the maximum di erence between the experimental and theoretical data for the three distributions is presented in Table 1.From Table 1, it can be seen that the maximum di erences of the strength between the theoretical results and experimental data are di erent, and their values are 0.192, 0.146, and 0.116, respectively.e maximum di erences of the elastic modulus are 0.138, 0.146, and 0.053, respectively.e Weibull distribution has less di erence between the experimental data and theoretical results either in strength or elastic modules.It can be concluded that the probability distribution law is satis ed slightly better by the Weibull distribution rather than the other distributions.
In order to ensure safety and save cost in engineering construction, the choice of strength design is important.In statistics, a con dence interval is a type of an interval estimate of the population parameter and is used to indicate the reliability of an estimate.e con dence levels for the sandstone are listed in Table 2.It can be seen that when the uniaxial compressive strength is taken as 61.50 MPa, the safe reliability is 90%, and when the strength is taken as 57.38 MPa, the undamaged probability is 95%.A con dence level of 95% means that there is a probability of at least 95% that the result is reliable.

Damage Constitutive Model of Rocks
So far, a variety of constitutive models for brittle or ductile materials have also been developed and implemented within the damage mechanic framework incorporating statistical variations [16][17][18].e statistical damage model presented by the authors in the previous article is used to formulate analytical constitutive relations for deformation behavior.Based on the linear damage mechanics [19], the stress-strain relation for the rock under uniaxial conditions can be expressed as follows: where σ is the stress applied to a damaged material, σ is the stress acting on the undamaged material, E denotes Young's modulus of the undamaged material, and D is a damage variable which takes a value between 0 and 1 corresponding to intact and fully damaged states, respectively.Equation ( 20) can be rewritten as   ( In probability theory and statistics, let N denote the number of total mesoscopic elements and N t denote the number of failed mesoscopic elements.e damage variable D can be directly de ned as e Weibull distribution is one of the most important continuous probability distributions.It was rst introduced to study the issue of structural strength and life data analysis and has been successfully applied to the mechanical failures of geomaterials [20][21][22].In this study, it is con rmed by experiment that the Weibull distribution is also valid.After the strain is regarded as a random variable, the following probability density function is used to describe the mesoscopic element failure behavior: where ε 0 is the threshold strain for the damage element; ε 0 0 for sandstone in this study; η is the mean value; and β is the shape parameter.When the axial strain reaches the value of ε, the number of damaged elements N t may be written as follows: Substituting ( 23) and ( 24) into ( 22), we have en, substituting (25) into ( 21), the following formula can be obtained: In order to demonstrate the rationality and accuracy of the statistical damage constitutive model proposed in this paper, some typical experimental curves are selected and compared with the theoretical model.e statistical distribution parameters in the constitutive model can be obtained by using the back analysis method.From Figure 6, we nd that the theoretical curves almost coincide with the experimental curves, which implies that the statistical damage constitutive model deduced in this paper could well describe the whole stress-strain process of the sandstone rock.

Conclusions
In this paper, the statistical mechanic behavior of the rock is investigated through a series of uniaxial compression experiments with acoustic emission.e following conclusions can be obtained.
e stress-strain curve can be approximately divided into four typical stages: the compaction stage, the linear elastic stage, the crack propagation stage, and the failure stage.A few events appear suggesting that a little damage of the rock happens in the compaction and elastic stages.As the strain increases to a critical value, the AE activity increases rapidly, and then, AE activity decreases with the further increase of strain.
e mechanical parameters were regressed with the Weibull distribution, normal distribution, and log-normal distribution models.
e results show that the Weibull distribution can more closely re ect the strength and elastic modulus distribution laws of the rock under uniaxial compression conditions.
A stochastic damage constitutive model for the rock was developed based on the continuous damage and probability 6 Advances in Civil Engineering theories.Comparing theoretical results of the model with experimental data, it is found that the stress-strain relationship could be well described by the damage constitutive model.

Figure 1 :
Figure 1: Collected and analyzed system of the acoustic emission.

Figure 2 :
Figure 2: e stress-strain curves of sandstone under uniaxial compression.

Figure 3 :
Figure 3: Acoustic emission characteristics for di erent failure types: (a) acoustic emission for tensile failure, (b) acoustic emission for combination failure, and (c) acoustic emission for shear failure.

Figure 4 :
Figure 4: Probability distribution characteristic of the strength under compression.

Figure 5 :
Figure 5: Probability distribution characteristic of the elastic modulus under compression.

Table 1 :
Maximum di erences between experimental and theoretical distributions.

Table 2 :
Con dence levels for the strength and elastic modulus of sandstone.