Mechanical Behavior of Shale Rock under Uniaxial Cyclic Loading and Unloading Condition

In order to investigate the mechanical behavior of shale rock under cyclic loading and unloading condition, two kinds of incremental cyclic loading tests were conducted. Based on the result of the short-term uniaxial incremental cyclic loading test, the permanent residual strain, modulus, and damage evolution were analyzed firstly. Results showed that the relationship between the residual strains and the cycle number can be expressed by an exponential function. +e deformation modulus E50 and elastic modulus ES first increased and then decreased with the peak stress under the loading condition, and both of them increased approximately linearly with the peak stress under the unloading condition. On the basis of the energy dissipation, the damage variables showed an exponential increasing with the strain at peak stress. +e creep behavior of the shale rock was also analyzed. Results showed that there are obvious instantaneous strain, decay creep, and steady creep under each stress level and the specimen appears the accelerated creep stage under the 4th stress of 51.16MPa. Based on the characteristics of the Burgers creep model, a viscoelastic-plastic creep model was proposed through viscoplastic mechanics, which agrees very well with the experimental results and can better describe the creep behavior of shale rock better than the Burgers creep model. Results can provide some mechanics reference evidence for shale gas development.


Introduction
Rock mass in underground is exposed to systematic cyclic loading during drilling, mechanical excavation, or due to mine seismology [1].It is clear that the mechanical properties of rock under cyclic loading differ dramatically from those under static loads [2,3].erefore, it is necessary to understand the response of rock mass under cyclic loading condition.
Laboratory testing is the main method to understand the mechanical behaviors of rock mass.Over the past few decades, considerable efforts have been made to assess the mechanical response of rock under cyclic loads.Liu et al. [4] investigated the deformation and fracture of coal rock under a cyclic loading and unloading condition.By proposing a new incremental cyclic loading-unloading pressure test method, Qiu et al. [5] analyzed experimentally the prepeak unloading damage evolution characteristics of marble from the Jinping II Hydropower Project.Tien et al. [6] studied deformation and fatigue behaviors of saturated sandstone subjected to cyclic loadings.Based on the load-unload response ratio (LURR) by the multilevel amplitude cyclic loading tests under the different confining pressure, Sun et al. [7] analyzed the damage evolution law and damage accumulation mode of sandstone.Fuenkajorn and Phueakphum [8] investigated the effect of cyclic loading on the mechanical properties of Maha Sarakham salt and found that the salt compressive strength decreased with the number of loading cycles according to a power function and that the elastic modulus decreased slightly during the first few cycles and kept constant until failure.
It is noteworthy that the abovementioned studies mainly focus on the mechanical behavior of rock under cyclic compressive loading and unloading condition and that while the instability of many geotechnical engineering does not appear immediately, such as underground excavation, mining tunnel, and drilling, it will last for some time after excavation.So, it is also very necessary to consider the timedependent mechanical behavior of rock mass under cyclic loading.
Shale rock may suffer cyclic loading during drilling, fracturing, and refracturing.To better understand the mechanical behavior of rock mass, the short-term and creep tests of a type of shale rock were conducted under uniaxial cyclic loading.Based on the short-term experimental test results, the deformation, modulus, and damage evolution were first analyzed.
en, in accordance with the creep experimental results, the creep behavior of the shale rock was obtained.Finally, by considering the characteristics of the Burgers creep model, a viscoelastic-plastic creep model was proposed through viscoplastic mechanics and its parameters were identified.Results can help better understand the failure mechanism of shale under cyclic loading, which has important practical significance for accurately assessing the effect of hydraulic fracturing and long-term stability of borehole.

Specimen Preparation.
e specimen used for the cyclic tests in this research was shale rock, a type of soft sedimentary rock, which was sampled from a Longmaxi formation on the east of the Sichuan Basin in China.Longmaxi Shale is the most commercially developed shale gas reservoir in China [9,10].
e natural rock sample with the average unit weight about is 2668 kg/m 3 , and with the same bedding orientation (horizontal), the natural water content of the rock sample is about 6%.According to the method suggested by the ISRM [11], the shape of all tested shale rock specimens is cylindrical with 50 mm in diameter and 100 mm in length approximately (Figure 1).

Testing Equipment.
Both conventional uniaxial compression experiments and cyclical experiments of specimens were carried out on a TFD-2000 rheology testing system, as shown in Figure 2. e test machine comprises axial loading, confining pressure loading, hydraulic pressure loading, numerical control, and measurement systems.
e maximum axial load is 2000 kN, and the measurement resolution is 200 N; the maximum axial deformation is 10 mm, and measurement accuracy is ±0.5%. is machine can perform the controlled tests and data acquisition/analysis by computer and robotized operations, which ensures that the tests are analyzed safely, timely, and precisely.

Testing
Procedure.Due to the shale gas pressure and fracturing fluid pressure, the surrounding rock is always subjected to cyclic loading in the development of shale gas field [12].Two different kinds of uniaxial cyclic loading and unloading tests were conducted in this research, which were short-term uniaxial cyclic compressive loading and unloading test and long-term compressive cyclic creep test, as shown in Figures 3(a) and 3(b), respectively.Before setting up the testing procedure, the uniaxial compression tests were conducted firstly to obtain the uniaxial compressive strength and the typical stress-strain curve of the shale rock specimen as shown in Figure 4.It can be obtained that the uniaxial compressive strength of the shale rock is 54.32 MPa.
e short-term cyclic loading and unloading test (Figure 3(a)) consisted of the following steps.e axial and lateral displacement sensor LVDTs were first adjusted to set the initial value, and then the first loading stress was applied to the designed value; meanwhile, the data of stress, displacement, and strain were recorded by computer software.Next, the stress was unloaded to zero.e preceding steps were repeated until the specimen failed in the last stress level loading process.A constant axial stress control rate of 0.5 MPa/s was selected during the each loading and unloading stage.In the whole procedure, the first loading stress is 27.16 MPa (50% of the uniaxial compressive strength), and then the stress level was increased by intervals of 8 MPa.
at is to say, the second to fourth stress levels are 35.16MPa, 43.16 MPa, and 51.16 MPa, respectively.
e procedure for the uniaxial creep test under cyclic loading condition (Figure 3(b)) was as follows.
e first loading stress of 27.16 MPa (same as the short-term test) was applied to the specimen at a rate of 0.5 MPa/s after the displacement sensor LVDTs were first adjusted to set the initial value.Next, the stress was kept constant for a time interval of about 60 minutes, during which the axial deformation was continuously monitored.And then, the axial stress was unloaded to zero at same rate of 0.5 MPa/s.e zero stress level was lasted for about 60 minutes, and then, the second stress will be applied to the specimen.e preceding steps were   Advances in Civil Engineering repeated until the fourth stress level, and the stress level was also increased by intervals of 8 MPa.Finally, at the fth stress level, the stress was maintained until creep failure took place.
To prevent the e ects of variation in diurnal temperature, a constant indoor temperature was maintained at 20 ± 0.1 °C during the entire tests.

Mechanical Behavior under Short-Term Cyclic Loading
3.1.1.Deformation Characteristic.Figure 5 shows the stressstrain curve of the shale rock specimen for a typical cyclic uniaxial test.It can be seen from Figure 5 that the curve of stress-strain formed hysteretic loops after each stress path of compressive loading-compressive unloading, and the area of the hysteretic loops increased with the increasing of the loading stress; this phenomenon is the characteristic of a material where viscosity might play an important role in its stress-strain behavior.It was also found that the stress-strain curve in the compressive loading stages was almost overlapped, which indicated that the stress-strain curve of shale rock has obvious memory e ect in the compressive loading stages.
As shown in Figure 5, considerable irreversible deformation of shale rock is observed.Figure 6 shows the relationship between the residual strains and the cycle number.With the increase in the stress levels and cyclic Advances in Civil Engineering loading numbers, the irreversible deformation becomes more signi cant, and the relationship between the residual strains and the cycle number can be expressed by an exponential function.Meanwhile, the in uence of stress levels on the irreversible deformation appears more evident than the cyclic loading numbers.
It can be also obtained that the rupture stress is 53.18MPa, which is 97.91% of the uniaxial compressive strength.Figure 7 shows a schematic view of the shale rock specimen after failure under uniaxial compressive loading and unloading condition, the main feature associated with the failure is the tension destruction.

Characteristics of Modulus.
Figure 8 shows the relationship between peak stress at each cycle loop and elastic modulus E S .Figure 9 shows the relationship between peak stress and deformation modulus E 50 .It should be claimed that, in Figures 8 and 9, the elastic modulus E S is the slope of the linear segment in the loading or unloading curve [12].e deformation modulus E 50 is the slope of the line between the point at which the stress value is 50% of each loading at the peak strength and the minimum stress point of the coordinates.
e deformation modulus E 50 can be expressed as follows: where E 50 is the deformation modulus; σ (1/2)p and ε (1/2)p are the stress and strain at the point of half peak strength under each cycle loading or unloading section, respectively; and σ min and ε min are the stress and strain at the minimum stress point under each cycle loading or unloading section, respectively.Advances in Civil Engineering e elastic modulus E S can be expressed as follows: where Δσ is the stress increment on the linear section under each cycle loading or unloading and Δε is the strain increment which is corresponded to Δσ.
From Figures 8 and 9, the deformation modulus E 50 and elastic modulus E S rst increase and then decrease with increase of peak stress under the loading condition, while the both of them increase approximately linearly with increase of the peak stress under the unloading condition.

Damage Evolution of Shale Rock under Cyclic Loading.
e existed initial cracks or bedding planes within shale make the quantity of damaged parts increase, and the area expands when the shale is under cyclic loading situation, and the increased damage causes the e ective stress on a microscale to increase.All sorts of the mesoscopic aws in shale rock are randomly distributed because rock is a product of long geological history [13].erefore, damage mechanics is a powerful tool to study the occurrence, propagation, and coalescence processes of these mesoscopic aws and their e ect on rock mechanical behaviors.Damage mechanics consists of the study of initiation and propagation of internal microcracks and defects under di erent loading condition, which may lead to deterioration of the mechanical properties of the material and eventually to failure [14].To describe the damage evolution process, the calculation of the damage variable is essential.
Based on the stress-strain curve of the rock mass under cyclic loading, we can calculate the constitutive energy and the dissipated energy at a given stress state, as shown in Figure 10.According to the previous studies [14 -16], the constitutive energy of the rocks is represented by the area between the complete stress-strain curve and the strain axis, and the area of OA 1 A 1 ′ is surrounded by the rst loading curve, and the strain axis can be de ned as absorbed energy; the area of O ′ A 1 A 1 ′ which is surrounded by the rst unloading curve and the strain axis can be de ned as released energy.So, the area of OAO ′ which is surrounded by the rst loading curve, the rst unloading curve, and the strain axis can be de ned as the dissipated energy.en the damage variable [14] at point A can be obtained as follows: where U OAO′ is the dissipated energy under rst cycle loads; U OA 1 A 1 ′ is the absorbed energy under rst cycle loads; is the released energy under rst cycle loads; and U is the constitutive energy.
According to the results from cyclic loading test on the shale rock, the damage variable can be calculated by using (3). Figure 11 shows the relationship between the calculated damage variables and the strain at peak stress.An exponential function was used to t the damage variables of the shale rock, with a coe cient of determination greater than 0.98.

Cyclic Creep Test Results.
In accordance with the testing procedure shown in Figure 3(b), the cyclic loading creep test result of shale rock was obtained, as shown in Figure 12. e stress levels were the same as these of shortterm cyclic loading tests, and it was lasted for about 2 hours at each cycle.
It can be deduced from Figure 12 that under stress levels of 27.16 MPa, 35.16 MPa, and 43.16 MPa, the axial and lateral strain rate increases rst and then gradually decreases to a constant value after a period time, and the similar variance trend was found in the unloading stage.According to the evolution of the creep strain, the creep curve under each stress level can be divided into decaying and steady stages except 4th stress of 51.16 MPa.Both axial instantaneous strain and steady strain of the specimen were increased with compression stress; the residual strains were also appeared Advances in Civil Engineering increase trend with increasing of former stress levels.e rock specimen fails eventually at 4th stress of 51.16 MPa which is 94.18% of the uniaxial compressive strength, also lower than the short-term rupture stress (53.18 MPa) indicated that time-dependent cycle loads have a signi cant impact on the damage evolution of the shale rock, the creep curve at this stress level stage includes three stages of the decaying, steady, and accelerating creep, and the whole creep only lasted for several minutes.

Creep Model and Parameter Identi cation.
Form Figure 12, it can be easily observed that there are obvious instantaneous strain, decay creep, and steady creep under each stress level, and steady creep strain approaches a certain constant value with the time when stress less than 51.16 MPa.While under the 4th stress of 51.16 MPa, the specimen shows the accelerated creep stage.erefore, it is necessary to simulate the creep behavior of rock by using a creep model.
Previous studies have shown that the viscosity coe cient of rock will decrease when stress exceeds the long-term ultimate strength [17][18][19]; since an exponential function can better re ect the relationship between the damage variables and the strain (Figure 10), we assume that the damage variables in the cycle creep stage can be expressed as an exponential function (as shown in ( 4)), and based on the characteristic of the Burgers creep model [20], a new viscoelastic-plastic creep model was proposed, as shown in Figure 13.
where α and β are the parameters controlling the damage evolution with time, which are related to the stress levels and cycle times; t is the creep time; and t 0 is the unit time, which is 1 minute in this work.
e creep equation of the viscoelastic-plastic creep model can be written as follows: where σ 0 and ε are creep stress and creep strain, respectively; E 1 is the elastic modulus of the Maxwell body; and E 2 and η 2 are the elastic modulus and viscosity coe cient of the Kelvin body, η 0 is the initial value of the viscosity coe cient and t is the creep time.By using the viscoelastic-plastic creep model (5) to t the creep curves, the relevant parameters of the viscoelasticplastic creep model were identi ed from data processing, as shown in Table 1; meanwhile, the Burgers creep model was also used to make a comparison with this of the proposed viscoelastic-plastic creep model.Figure 14 shows the calculated results of creep curves and the tested ones which were below stress level of 51.16 MPa, and the tting curve and tested curve under stress level of 51.16 MPa are shown in Figure 15.From the two gures, it was found that the calculated results of the proposed model are in good agreement with the tested results than the Burgers creep model under each stress level, and the proposed model also can better re ect accelerated creep stage of the shale rock specimen.
e correlation coe cient square (R 2 ) of four curves is more than 0.959, which indicates that the prediction precision of

Conclusions
Two di erent kinds of cyclic loading and unloading tests were conducted in this study for characterizing the mechanical behavior of shale rock, which were short-term uniaxial incremental compressive cyclic loading and unloading test and long-term uniaxial incremental compressive cyclic loading and unloading creep test.From the test results, the following conclusions can be drawn: (1) e relationship between the residual strains and the cycle number can be expressed by an exponential function.e deformation modulus and elastic modulus rst increase and then decrease with increase of peak stress under the loading condition, and both of them increase approximately linearly with increase of the peak stress under the unloading condition.
(2) e cyclic loads formed hysteretic loops, and the area of the hysteretic loops increased with the increasing of the loading stress, which can be used to interpret energy dissipation of the specimen.Based on the energy dissipation, the damage variables under short-term cyclic loading condition were analyzed, and it shows an exponential increasing with the strain at peak stress.
(3) ere exist obvious instantaneous strain, decay creep, and steady creep under each loading or unloading stage, and under the stress of 51.16 MPa, the specimen shows the accelerated creep stage.(4) A viscoelastic-plastic creep model was proposed.
Comparison of the test results with the model predictions shows that the proposed viscoelastic-plastic creep model is capable of describing the creep behavior of shale rock subjected to cyclic loading and unloading condition.
Finally, it should be noted that further studies both in triaxial tests and in theory are necessary to gain more relative data for better understanding the cycle load mechanical behavior of shale rock.Numerical simulation of engineering projects associated with shale rock should also be taken into account in the future.
is study was partially supported by the National Natural Science Foundation of China (Grant no.41302223), Science and Technology Plan Projects of Chongqing Administration

Figure 3 :Figure 4 :
Figure 3: e sketch of (a) short-term cyclic loading test and (b) cyclic loading creep test.

Figure 7 :Figure 8 :
Figure 7: Schematic view of the shale rock specimen after failure.

Figure 10 :Figure 9 :
Figure 10: Calculation model of dissipated energy and constitutive energy.

Figure 12 :Figure 13 :Figure 11 :
Figure 12: Curves of axial strain versus time of shale rock specimens under cyclic loading.

Figure 14 :
Figure 14: Comparison between calculated results and experimental results below the stress level of 51.16 MPa.

Figure 15 :
Figure 15: Test results and model calculation curves of specimen under stress levels of 51.16 MPa.

Table 1 :
Parameters of the viscoelastic-plastic creep model.