Abstract

Aiming at the damage evolution and energy release characteristics of sandstone in failure process, uniaxial and cyclic loading experiments were conducted with the average temperature, infrared temperature field, acoustic emission (AE), and displacement field being monitored simultaneously. The results show that in the low stress stage, the initial pores and fissures are compressed and the thermoelastic effect is insignificant. At the same time, the AE events are active, but most of the events are released with low energy. Then, the surface of the specimen shows stratified displacement phenomenon in the vertical direction. In the following stage, the thermoelastic effect is significant, and the average temperature of specimens shows a linear correlation with stress, whereas the AE is relatively inactive, and vertical surface deformation shows further homogenization and saliency. Because of the top and bottom constraints, the horizontal displacement field shows axis symmetrical distribution of double half-ellipse like a drum. During the yielding phase and the rupture moment, the average temperature rises obviously because of the heat from the friction of the rupture surface, and the temperature field also appears the phenomenon of differentiation; in the meantime, a large amount of AE events occur and the proportion of high energy events increases; further, the deformation field is significantly gathered near the ruptured position. Finally, the time sequence of the multiple physical parameters (AE parameters, average temperature, temperature fields, horizontal and vertical displacement fields) is summarized, which can be a reference for the stress state and failure analysis.

1. Introduction

The failure of rock is the result of interior crack initiation, propagation, and penetration. The whole rupture process of rock is a complex physical phenomenon. With the increasing load, many physical parameters, including the force, deformation, elastic wave [1], thermal energy [2, 3], electric potential [4], and resistivity [5], will occur abnormal variations. These abnormal signals are very important for predicting the instability and quantifying the damage degree.

The energy evolution plays an important role in explaining the mechanism of rock engineering disasters. The energy input by the loading machine can transform and release in multiple forms, such as elastic deformation, plastic deformation [6], infrared heat radiation, elastic wave, and electric energy. The initiation and expansion of microcracks will cause energy release in the form of acoustic emission (AE), which has been used for the prediction of quantitative rock damage widely [7]. The friction caused by microcrack and the thermoelastic effect will increase infrared temperature [8]. Further, the relationship between the infrared radiation temperature and stress has been discussed by bending and shearing experiments [9]; the experiment makes it possible to measure the stress level by the noncontact method. On the other hand, the specimens deformation (elastic and plastic deformation) is also a result of energy transformation, and the digital image correlation (DIC) technology can be used to monitor and inverse the full-field deformation [10]. With the characteristic of high frequency and noncontact, the DIC technology is usually used in the complicated and harsh circumstances, such as in high temperature [11] or polluted environment. Compared with the traditional measure method such as using extensometer or strain gauge, the full-field monitoring of DIC has a great advantage in obtaining the full-field deformation rather than the strain at a point [12, 13].

From the above studies, it is clear that the variation characteristics of infrared radiation temperature, full-field deformation, and AE are quite meaningful for quantifying damage and explaining the failure mechanism [14, 15], and the monitoring of single physical quantity has been widely researched. But the joint monitoring of multiple physical quantities is more comprehensive for understanding and analyzing the internal damage evolution process. By now, some joint monitoring systems have been established, such as the multiparameter comprehensive monitoring system (including microseismic, charge induction, AE, microcomputed tomography, and scanning electron microscopy) [16, 1719]. The joint monitoring systems show that comprehensive use of multimethod monitoring is a potential method for analyzing the rock failure process in multiple views.

Thus, in this paper, the uniaxial and cyclic loading of sandstone specimens have been conducted. During the loading process, the average temperature, infrared radiation temperature field, AE activity, vertical displacement fields, and horizontal displacement fields have been monitored simultaneously. Based on the obtained multiple physical parameters, the chronological order of multi-physics parameters has been analyzed, and the chronological order can be a guidance for the failure warning and damage quantifying.

2. Specimens Preparation and Loading Schemes

The mesozoic sandstone drilled from Xiaojihan coal mine, located in the city of Yulin, China, was prepared for the experiments. All the specimens were cut into standard cylinder shape with 50 mm diameter and 100 mm length in accordance with the ISRM. The tests were carried out on the TYJ-500 kN loading machine with a loading speed of 0.001 mm/s in the axial direction.

During the experiments, the average temperature was measured by Fluke 568 Thermometer with wavelength of 8–14 µm. Infrared temperature field was collected by H2640 imager with wavelength of 8–13 µm. The AE monitoring was conducted by PAC PCI-II AE system, which includes a Nano30 piezoelectric sensor with response frequency of 125–750 kHz, and the sensor was arranged in the middle area of the cylindrical specimen. The images used for the DIC were obtained by XTDIC speckle system with 2.3 million pixels. The related testing equipment is shown in Figures 1 and 2.

In the experiment, 12 sandstone specimens were prepared for the experiment, the loading paths contained uniaxial loading and cyclic loading, in the uniaxial loading experiments, the specimens were compressed to failure directly, and in the cyclic loading process, the maximum pressure of each cyclic loading increased 10 kN until the failure. The loading paths and monitoring methods are shown in Table 1.

In the experiment, the energy of average temperature and infrared image correspond to thermal energy, AE energy corresponds to elastic wave energy (from a microfracture inside the rock), and displacement field of DIC corresponds to elastic and plastic deformation energy; all the energy come from the energy input by the press machine during the loading process.

3. The Variations of Multiple Physical Parameters during the Loading Process

3.1. The Average Temperature Variation Laws during Loading Process

Some theories and experiments have concluded that the solid material shows thermoelastic and thermoplastic effect during the loading process. Compared with thermoplastic, the study of thermoelastic theory is more systemic and deeply [20]. In a system considering heat transmitting, the relationship between temperature and elastic deformation is shown in Equation (1) [21]:where T is the absolute temperature, Cε is the specific heat capacity of the material at constant strain, Q is the heat input or output, ρ is the mass density, σij is the stress change tensor, and εij is the strain change tensor.

For isotropic elastic material, the Equation (1) is simplified to:where α is the coefficient of linear thermal expansion, Cp is the specific heat capacity at constant pressure, and ∑σij is the first stress invariant.

Aimed at the thermal phenomenon in the loading process, massive studies have been developed in the uniaxial loading test [5, 22]. However, because of the limit of the sampling frequency, few experiments are conducted to obtain the temperature data directly by a thermometer. Based on the high-frequency and high-precision temperature acquisition instrument, the average temperature of sandstone with different loading paths has been investigated.

As shown in Table 1, six specimens were conducted for the uniaxial loading experiment. In the average temperature experiment, there are four specimens showing declining trend at low stress stage, then during the elastic and plastic stage, the average temperature increases stably, and finally reaches highest temperature at the failure moment. The typical temperature variation trend is shown in Figure 3(a). However, there are two specimens showing declining trend during the whole loading process, and only an abrupt increase occurs at the failure moment.

For the cyclic loading experiments, also six specimens were conducted, as shown in Table 1. All the specimens show the same variation trend. At low stress level, the average temperature shows declining trend, and in the following cyclic loading, the average temperature shows notable linear positive correlation with the stress, namely, the thermoelastic effect is remarkable. Moreover, the temperature shows serious increasing trend at the failure moment, and the typical average temperature variation is shown in Figure 3(b).

The reason of average temperature dropping at the initial loading stage is that the fracture and cracks in the specimens are closed in the initial stage, and the specimens mainly show plastic deformation, so the thermoelastic effect is not remarkable. Moreover, the steel machine starts contacting with the specimens closely, and the steel is a better heat conduct material resulting in the heat of the specimens transferring to the machine, so the temperature of specimens decreases. After a period of contacting, the heat transferring between rock and machine reaches a balance, and also the specimens enter elastic stage, so the primary stress and the primary strain show linear relationship, and because of thermoelastic effect, the average temperature variation is consistent with the stress increasing and declining. At last, when the specimen reaches its peak strength, the deformation converges in the failure zone, and the friction effect of particles is also serious; meanwhile, the elastic wave may cause particles vibration at a high frequency; all the above factors may cause the temperature increasing, so the average temperature increases sharply at the failure moment.

3.2. The Infrared Radiation Characteristic of Temperature Field

The thermometer has the characteristic of high frequency and high sensitivity, but the temperature obtained is the average value in a small zone, and it can not reflect the temperature field variation. With the purpose of analyzing the temperature field, the infrared radiation images have been collected in the test simultaneously.

In the previous study, it is concluded that the infrared image shows alienation phenomenon and the average infrared radiation temperature (AIRT) shows abnormal trend. Some studies show that the failure pattern may also affect the temperature, the shearing fracture may cause the temperature increase, and the splitting fracture may cause the temperature decline [3, 23].

In this part, the H2640 infrared radiation imager was used to monitor the temperature field of the specimens surface. Due to the effect of the environment variation (including the light and human body radiation), during the test, the windows and curtains were closed, and the unnecessary activities were forbidden.

The typical results are shown in Figure 4; it can be found that the temperature is not evenly distributed on the specimens surface. When the rock reaches its peak strength, the infrared images show differentiation characteristic, the temperature of nonfailure zone keeps stable, and that of failure zone releases massive heat because of shearing crack and the friction of particles; the released energy leads to a seriously increase in the temperature.

Further, with the purpose of studying the temperature variation in different zones, a measuring line (white line) is located on each picture, as shown in Figure 4, and the temperature value on the measuring line has been extracted. The temperature of the measuring line with the length of the measuring line (pixel points) is shown in Figure 5. It is clear that the temperature increases seriously in the failure zone, and the highest temperature (as shown in Figures 5(a) and 5(b)) of typical specimens of A16Q and H21 increases to 1 and 5°C, respectively, which is similar with the previous research [24]. In the infrared image, it can be found that the energy produced by the shearing crack and friction is much more than the thermoelastic effect.

The variations of average temperature and infrared images show that in the initial loading stage, the specimens show more plastic deformation, and the thermoelastic effect is not obvious; due to the heat transforming between the rock and machine, the average temperature shows declining trend. In the elastic stage, thermoelastic effect is notable, and the temperature shows linear relationship with the stress. At the failure moment, massive heat is produced in the failure zone because of shearing cracks and friction, and the infrared image shows differentiation phenomenon, the average temperature increases sharply.

3.3. Characteristic of Acoustic Emission Parameters in the Loading Process

In the rock failure process, the radiation energy contains the thermo energy, which is expressed in the temperature increasing and elastic wave energy, which release in the form of AE. AE is the energy released by cracks initiation and propagation, and it is widely discussed in the field of crack location [25, 26], damage evaluation [27], and energy release [28] The elastic wave parameters include amplitude, energy, ring counts, rising time, and frequency; these parameters are meaningful for the damage analysis [29], and also these AE parameters can be calculated from the AE waveform. The part-related parameters are shown in Figure 6.

In AE studies, the ringing counts, energy, and b value based on the AE wave are chosen for studying the energy release and damage evolution. The b value is a universal parameter in the earthquake domain; it means that the earthquake frequency declines in the exponential form with the earthquake magnitude [30], and the b value decreases to a minimum at the time of failure [31]. There are many similarities between the earthquake and AE event except the magnitude and frequency, and so many methods and theory of earthquake are cited in the AE studies. In the AE monitoring system, the Nano30 sensor receives the voltage signal and transmits the signal into the acquisition card. So, compared with the earthquake signal, the AE number and the amplitude obey the following equation [32] as follows:where AdB is the signal maximum amplitude, N is number of events larger than magnitude AdB, and a and b are fitting parameters.

The fitting parameters present the proportion of the high amplitude events in the all the events. A higher b value will be obtained if most events are small events with less energy, indicating less cracks or cracks growing slowly, whereas a lower b value represents most AE events that are large events with high energy, indicating many cracks or cracks growing fast. A undulate b value meaning unstable cracks expand [6, 33], so the b value is related with the damage degree closely.

In the loading process, the typical temporal variation of b value and AE energy is shown in Figure 7. It shows that in the compaction stage, the initial pore and cracks are closed, and massive AE events occur; however, the b value is also at a high level; it means that most of the AE events produced in the initial stage are low energy. During the elastic stage, the AE events are relatively inactive, and the b value shows stable value with little variation. However, in the yield stage, the AE energy increases seriously; it means that the inner damage evolves violently; meanwhile, the b value vibrates acutely meaning that the AE events with high and low energy occur alternately. At the failure moment, the AE energy reaches its maximum, and the maximum energy value is often hundredfolds of previous value (a break is located at the counts axial, as shown in Figure 7). Further, the b value declines clearly; it means that cracks with higher energy appear in the rock, and then the failure comes. It can be seen that the AE parameters (AE energy and b value) show different variations in different stress stage, which can be used to provide early warning for fracture instability.

3.4. Deformation Field Characteristic with the Method of Digital Image Correlation

Part of the work done by the press machine is used for radiation in the form of thermo and elastic wave, and the rest of the work is used for changing the shape of the specimens. The deformation contains elastic and plastic deformation. In the traditional deformation monitoring experiment, the deformation is usually measured by extensometer or strain gauge, which is a quite small zone deformation with the probability of losing whole deformation characteristic. In this part, a high-function digital camera and image correlation technology are used to study the displacement field characteristic in the loading process.

DIC is a new optical method, which captures images of the rock surface and obtains full-field deformation by comparing the current images with the reference image. Here, the Ncorr code is used for the DIC, Ncorr is an open-source 2D DIC program compiled by Matlab and C++ language. With the purpose of subset displacement tracking, an improved DIC based on the conventional Newton–Raphson algorithm is compiled in the code [34]. By processing the displacement field with the method of differentiation, the Green–Lagrangian strain field can be obtained [11, 35]. The image processing mainly contains setting reference image, defining deformation image and target zone, arranging reference seeds, and calculating the displacement and deformation field. Further, the code flow can map the physical dimension with the pixels, so the displacement size in the object zone is the real deformation scale. Some studies show that the deformation result of Ncorr is consistent with the commercial 2D DIC packages [36]; it means the deformation processed by the Ncorr is reliable and reasonable.

Figure 8 shows the vertical displacement variation at different stress level. According to the loading machine principle, in the axial loading process, the upper platen is fixed, and the lower platen moves toward to the upper one as the setting rate (0.001 mm/s). It can be found that at low stress level (Figure 8(a)), the vertical displacement field shows disordered pattern. With the increasing stress, the displacement field shows evenly stratified characteristic. The lower and upper faces have the biggest and smallest displacement, respectively (Figure 8(b)). When the stress reaches the peak strength, the displacement image shows the concentration and heterogeneity, and the zone with concentrated displacement contour is also the failure zone in the next step.

In the uniaxial loading process, there is no limit in the radial direction around the specimens, so the expansion effect exists. Figure 9 shows the horizontal displacement characteristic. In the low stress level (Figure 9(a)), the displacement shows no regular pattern, similar with vertical displacement. As the stress increases (Figure 9(b)), the horizontal displacement shows particular characteristic. In the axial line, the horizontal displacement is quite small, and the zone far from the axial line shows obvious expanding trend. Moreover, because of the friction between the loading platen and specimens, the expanding of the top and bottom of specimens is limited, so the displacement in the horizontal direction is much smaller than that of other zone; thus, the middle of the specimen shows biggest displacement. The whole horizontal displacement field shows axis symmetrical distribution as double half-ellipse. In the failure moment, the horizontal displacement also shows heterogeneity characteristic clearly.

In above DIC analysis, it can be summarized that the length of specimens in the axial direction is evenly compressed, as shown in Figure 10(c), and the radial size shows the phenomenon of expanding, but there is a friction force between indenter of the machine and the top (bottom) end faces of the specimen, which restricts the radial expansion of the two end faces, so the middle part bulges freely and end faces (top and bottom) restrict bulging, as shown in Figure 10(d). The displacement in vertical direction shows uniformly-layered characteristic and the displacement in horizontal direction shows axis symmetrical distribution as double half-ellipse. So, the cylindrical specimens exhibit deformation of drum shaped, as shown in Figure 10(b).

4. The Time Sequence Characteristics of Multiple Physical Parameters Variation in the Process of Rock Failure

It can be seen from the above analysis that the characteristic, including surface temperature field, AE activity, and surface displacement field of the specimen, shows significant regular variation during the loading process. Further, the change characteristics and sequence of different physical quantities have certain guiding significance for explaining the damage mechanism.

Taking the H20 as an example, the stress curve, the surface average temperature curve, and the AE characteristic curves are integrated, as shown in Figure 11(a), and three marked points (I, II, III) are arranged on the stress curve. It can be seen that the AE ring counts increase rapidly with a high b value before marked point I (stress level 15.8%), so the AE energy value is relatively low. It means that the initial fissures inside the rock are compacted. Between marked points I and III (stress level 65.84%), the AE ringing counts increase slowly, and the b value is still at a high level, indicating that a small amount of low energy AE events is generated in the rock, and the rock is at elastic deformation stage. After marked point III, the AE ring counts increase continuously, and the b value continues to decrease; it means that the high energy ruptures in the rock occur. At the moment of the rupture, the AE events suddenly increase, and the b value suddenly decreases. Before marked point II (stress level 28.5%), the thermoelastic effect in the compaction stage is not significant, and the specimen releases heat to the outside, causing the average temperature to decrease continuously. After marked point II, the average temperature begins to rise steadily, it shows a positive correlation with the stress level, and the average temperature increases about 0.2°C before the peak strength. At the moment of failure, the average temperature increases sharply due to the sliding friction of the mineral particles and the release of fracture energy.

Moreover, the representative pictures of the displacement fields in the vertical and horizontal directions of H20 are shown in Figures 11(b) and 11(c), and the serial numbers are correspond with the marked points (A. B. C… H) on the stress curve, as shown in Figure 11(a); it can be seen that at point A (stress level 2.13%), both the vertical and horizontal displacement fields appear irregular and random characteristic. At point B (stress level 9.04%), the vertical displacement field begins to be stratified, and the horizontal displacement field is still irregular. Then, as the stress increases, the stratification of vertical displacement field becomes more obvious. Until point D (stress level 70.21%), the horizontal displacement field begins to appear axisymmetric semielliptical characteristics, which means that the specimen begins to appear regular expansion in the radial direction. With the further increase of the stress, at point G (stress level 99.14%) and H (stress level 99.78%), the vertical and horizontal displacement fields begin to show the characteristics of deformation accumulation, respectively, which mean that the rock is about to appear instability.

The AE, infrared temperature field, average temperature, and vertical and horizontal displacement fields of all specimens in the uniaxial loading experiment are processed as the above method. To avoid errors and contingencies, the stress levels of each physical quantity changing characteristic points are averaged. The comprehensive analysis of the time sequence characteristics of multiple physical parameters during the loading process is shown in Figure 12. As shown in Figure 12, the time sequence and variation laws of AE, infrared temperature, vertical displacement, and horizontal displacement with the stress level are clearly presented; the sequence and the sensitivity of physical quantities can be an effective reference for analyzing the rock failure mechanism and determining the degree of damage.

5. Conclusion

In the uniaxial and cyclic experiments, the multiple physical parameters, including AE, infrared temperature field, average temperature, and vertical and horizontal displacement field, have been monitored simultaneously, and the variation and the time sequence of physical parameters have been obtained. The following conclusions can be drawn:(1)In the initial loading stage, the AE events are active, but most of the AE events are released with low energy. In the elastic stage, the AE events are not active, and the energy of most AE events is lower. In the high stress level, the AE events increase rapidly, and the b value declines constantly, which mean the proportion of higher AE energy is increasing and the failure moment is coming.(2)In the initial loading stage, the temperature shows declining trend due to the specimens releasing energy to the outside. In the elastic stage, the specimens show obvious thermoelastic effect, and the temperature and the stress show linear relationship. At the failure moment, the temperature increases suddenly, and, further, the infrared images show differentiation and concentration phenomenon; the temperature in failure zone increases clearly because of particles friction and cracks propagation.(3)In the initial loading stage, the vertical and horizontal displacement fields both appear irregular and random phenomenon. In the elastic stage, the vertical displacement shows stratified evenly characteristic, and at high stress level, the horizontal displacement shows axis symmetrical distribution as double half-ellipse. Near the failure moment, the displacement fields both show clearly inhomogeneous characteristic.(4)The multiple physical parameters of all specimens are comprehensively analyzed, and the time sequence characteristics of AE activity, average temperature, temperature field, and digital displacement fields (horizontal and vertical direction) are obtained. This multiple physical sequence feature can be used for analyzing rock failure mechanism and early warning of rock instability.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work is funded by Science Research Project of Hebei Education Department (QN2020260), High-level Talents Research Start-up Project of Hebei University (521100221025), National Natural Science Foundation of China (2013CB227902), and Natural Science Foundation of Jiangsu Province (BK20200993).