Abstract
Recently, the holding states of nanoindentation experiments have been widely used to analyze the time-dependent deformations of various rocks, and the dynamic mechanical analysis (DMA) method seems to be more applicable than the quasi-static mechanical analysis (QMA) method when the influence of creep deformation on mechanical properties of rocks was analyzed. However, the former method causes an abnormal behavior during the creep holding stages that was not clearly interpreted.2 Consequently, in this study, by amplifying the oscillation of the DMA method, the mechanical mechanism of this phenomenon was explained. Experimental results confirm that the rheological deformation of rocks consists of the creep deformation (depth increasing) and the elastic aftereffect deformation (depth decreasing) during the creep time with small oscillation; once the elastic aftereffect deformation exceeds the creep deformation, the abnormal behavior can be observed. Besides, some other abnormal behaviors might be found for other rock materials when the DMA method with different oscillations is used, which illustrates the complexity and limitation of applying this method. Thus, the QMA method was recommended to investigate the above questions in future studies.
1. Introduction
In the deep underground, the rheological phenomenon is the most important mechanical property of geomaterials, which negatively impacts the safety of engineering structure[1–3], coalbed methane (CBM) and shale gas production [4–6], geological disposal repository effectiveness [7], and so on. As a result, to enhance the understanding of this phenomenon is of importance for deep rock mechanics. The most popular method to investigate the creep behavior of geomaterials is based on in-scale core samples, which is time-consuming and challenging to obtain enough intact specimens [3, 8, 9].
Recently, the holding states of nanoindentation experiments have been widely used to analyze the time-dependent deformations of biomaterials [10], metals and alloys [11, 12], ceramics [13, 14], polymers and composites [15–18], types of cement [19, 20], and even crystals [21, 22]. Besides, the estimation of the macrofracture toughness of minerals and rocks is given in [23], while the analysis of the micromechanical properties of mineral crystal grains is becoming more and more extensive [24–26]. All those above mean that the experimental technique was a robust way to conduct creep tests. And it was convenient to obtain the sizeable samples from broken rock blocks since the size requirement was just about one centimeter. Consequently, the nanoindentation technique starts being adopted in the field of geomaterials’ rheology, especially for shale in the petroleum and gas industries, helping reveal some undiscovered problems at the macroscale. Taking gas shale formations as examples, the role of the organic matter (OM) in gas shale is “the dominant creeping phase,” which drives the creep rate [27]. Analogously, the hard minerals in shale collected from Bakken formation have much smaller creep deformation than that of softer ones [28]. The quasi-static mechanical analysis (QMA) method and the dynamic mechanical analysis (DMA) method that was regarded as two basic approaches to conduct nanoindentation creep tests [29] were respectively used in above the two articles. And for other geomaterials obtained in coal mines, such as coal, mudstone, and sandstone, the quasi-static mechanical analysis method was also widely used to reveal their time-dependent deformation phenomenon [25, 30].
However, this technique also causes some debate in the field of rock mechanics. Using the quasi-static nanoindentation method, prior research results presented that “as the creep time increases, the Young’s modulus, hardness, and fracture toughness decreases” [31], which means that the depth of indents increases with holding time. While, for the convenience of conducting a continuous and accurate measurement, the DMA method was used, showing that “the dynamic mechanical properties such as storage modulus and hardness were found to first decrease and then increase after the holding time exceeds a certain value” [32]. “Elastic recovery during the holding period” and “stress relaxation” were considered as the reasons for this abnormal behavior [32], which cannot explain why the quasi-static nanoindentation method yields normal results [31] and brings much debate on whether the DMA method is suitable to discuss the rheological responses of rocks. Consequently, further explanation for different results is deemed to be made to help comprehend the rheological characteristics of the DMA method and propagate the application of the nanoindentation technique in rock mechanics.
Thus, in this article, considering that the only difference between the two methods is that only the DMA method has small oscillations, the method named oscillation amplification was adopted to try to explain the abnormal behavior. This paper was structured as follows: first, the oscillation amplification method was clearly introduced (Section 2); second, in Section 3, the experimental results were given in detail to illustrate the existence of elastic aftereffect beforehand; third, the discussions in Section 4 calculated the nanoindentation rheological strain-time curves, determined the constitutive model of rock nanoindentation rheological behaviors, and analyzed several situations that might be encountered if the DMA method was used in the future study. In conclusion, the rheological deformation of rock investigated by using the DMA method is very complicated, showing the limitation of using this method.
2. Oscillation Amplification Method
In order to overcome the main limitation of the quasi-static mechanical analysis (QMA) method in which the hardness and modulus values can only be calculated at the initial part of the unloading curve, the dynamic mechanical analysis (DMA) method with quasi-static force being kept constant and small oscillation being superimposed was proposed and used [28], as shown in the left part of Figure 1(a). However, as mentioned above, prior study results found an abnormal behavior when the holding time exceeds a certain value by using this experimental method. To try to explain this in terms of rock mechanics, in the right part of Figure 1(a), the method to amplify the small oscillation was schematically illustrated.

In short, the deformation response of a small oscillation can be replaced by a holding period at the maximum load followed by another holding period after unloading to a small load (minimum load) value. The whole creep time has numerous alternately occurring maximum load and minimum load-holding periods. The load-holding periods at the maximum load are seated together, first (stage b-c) and then the periods at the minimum load happen together next (stage d-e). The two holding periods are connected by a ramp unloading process (stage c-d) in practice. The experimental protocol in Figure 1(b) further amplifies the different rheological responses during stage b-c and stage d-e by setting the maximum load at 50 mN and then unloading to 10% of it (minimum load). Besides, to calculate the hardness and elastic modulus of materials, before both load-holding periods, the typical nanoindentation experiment (stage o-a) was also conducted in this study. It consists of the loading segment, holding segment of 5 s, and unloading segment with the loading and unloading rates of 100 mN/min.
Many previous study results displayed that the deformation during the load-holding period at 10% of the maximum load originates from temperature drift when the rheological properties of other materials were investigated, such as thin Cu film [33] and polycrystalline Al film [22]. However, this method implicitly assumed that the time-dependent deformation during this period was ignored [34] and nanoindentation experimental systems have low resolutions. Thus, to not only consider the time-dependent deformation during the load-holding period at 10% of the peak load but also to avoid the influence of this error of equipment-temperature drift, the thermal equilibrium time of the Aton Pharr NHT3 system equiped with Berkovich indenter was more than 1 h, and compensations were autumatically formulated [35]. Both the thermal equilibrium and compensations were done previously based on other materials and the thermal drift also reaches much less than 0.01 nm/s [35], since we found the deformation of rocks (shale, mudstone, and coal) during the holding period of 10% of the maximum load is time-dependent rather a linear increase [36], which is too large to be ignored. After compensating, three random indentations were made on each rock sample to research the time-dependent deformation characteristics.
3. Results
Taking shale, mudstone, and coal, derived from China’s coal mines as examples [37], three nanoindentation tests labeled sample-1, -2, and -3 following the protocol in Figure 1(b) were made on the random locations on each sample. The hardness and elastic modulus of three nanoindentation tests on each sample was computed by using the Oliver-Pharr method through stage o-a in Figure 1(b) [34]. Calculated results show that the mechanical properties of different positions on each sample are various. For example, the hardness of shale-1, shale-2, and shale-3 are 532.71 MPa, 523.53 MPa, and 801.05 MPa; and that for mudstone-1, mudstone-2, and mudstone-3 are 598.46 MPa, 376.08 MPa, and 217.32 MPa; and for coal-1, coal-2, and coal-3 are 388.43 MPa, 317.08 MPa, and 153.88 MPa. The elastic modulus of coal and rock increases linearly with their hardness [38]. Obviously, at the mesoscale, geomaterials are typical heterogeneous materials, consisting of various mineralogical constituents in different positions.
The deformations during stage b-c (-C with solid symbols) and that during stage d-e (-E with open symbols) of shale, mudstone, and coal are given in Figures 2(a), 2(b), and 2(c), respectively, where the deformation-time curves of all tests on different surface locations start from 0. The increasing depth is positive deformation. Besides, as shown in Figure 2, just like the hardness and elastic modulus in different positions on each sample are different, the deformations during stage b-c (-C) and that during stage d-e (-E) are different, which not only shows that the rheological properties of different positions are various but also illustrates that these displacements of stage d-e (-E) are not caused by thermal drift.

(a)

(b)

(c)
For shale, mudstone, and coal samples, during the period of stage b-c, the indentation depth or deformation is rapidly growing at first and then slows down over time. Throughout stage d-e, indentation depth decreases with time for mudstone, coal, and shale-3. However, for shale-1 and shale-2, when the holding time exceeds a certain value, the indentation depth begins to increase. During the holding period of 50 mN, the deformations or depths of indentations are not abnormal by using the quasi-static nanoindentation method [31]; while during the holding period of 5 mN is negative and decrease with time or firstly decline then slightly increase. The creep deformation during the holding period of the maximum load was widely detected and investigated as mentioned above, but there are limited studying results to discuss the depth decreasing phenomenon after unloading to a small value, probably causing some confusion of the time-dependent deformation during the unloading period in any form, such as thermal drift [33] and small oscillation [32]. Consequently, understanding this phenomenon is of significance.
4. Discussions
4.1. Nanoindentation Rheological Strain-Time Curves
In previous nanoindentation tests, only the creep deformations were observed as an increase in depth during a hold period at the maximum load in the load-displacement data for a range of materials [29]. However, this method was out of the framework of rock rheology and limited the analysis of rheological properties of rocks to some extent. Thus, converting the nanoindentation depth-time data to the strain-time data is necessary. For the nanoindentation rheological experiments, the stress and strain rates can be expressed as [39, 40]where is the stress, MPa; is the strain rate, s−1; is the indentation depth changing with time, nm; is the time, s; and is the geological constant of indenter (for Berkovich indenter, = 0.0407 [40]).
In rock rheology, the creep deformation can be recorded during the constant load time, and after unloading to a small load value or 0, elastic aftereffect deformation happens [41, 42].
4.1.1. Nanoindentation Creep Stage (Stage b-c)
In this stage, the indentation depth () increases with time. Equation (2) was integrated to get the nanoindentation creep strain where is the nanoindentation creep strain, and is the indentation depth at the point of b, nm.
4.1.2. Nanoindentation Elastic Aftereffect Stage (Stage d-e)
In stage d-e, the indentation depth decreases with time during the initial part, meaning that the deformation is time-dependent. The elastic aftereffect is a phenomenon of delayed recovery of the elastic deformation [41, 42]. The strain expression in the elastic aftereffect stage can also be obtained by integrating equation (2). In this process, it should be noted that the elastic aftereffect strain is a negative value.where is the elastic aftereffect strain and is the indentation depth at the point of d, nm.
Therefore, according to the typical nanoindentation rheological experiment process (Figure 1(b)) and the above nanoindentation rheological theory, the depth-time data of the nanoindentation experiments can be converted into the strain-time curves consisting of creep and elastic aftereffect stages, as shown in Figure 3 taking shale-3, mudstone-3, and coal-3 as examples.

Both the nanoindentation creep and elastic aftereffect curves include two phases: the transition deformation phase and the stable deformation phase. During the transition deformation phase, at first, the creep strain increases sharply with time, and then the strain rate slows down; during the stable deformation phase, the creep strain increases linearly at an almost constant rate. After unloading to 10% of the maximum load, the elastic strain gradually recovers with the first fast strain rate and then a small constant rate during the whole period of elastic aftereffect for shale-3, mudstone-3, and coal-3. However, for shale-1, the strain increases during the stable deformation phase of elastic aftereffect, where its rheological mechanism was not explained, too. Consequently, discussing the constitutive model of nanoindentation rheological curves is of importance for making sure the deformation characteristic during both load-holding periods.
4.2. Nanoindentation Rheological Constitutive Model
Considering the viscoelastic properties of shale, mudstone, and coal samples and the rheological deformation properties of nanoindentation experiments. The simplified Burgers model was tried to be utilized to fit the nanoindentation experimental data, and this model is schematically presented in Figure 4.

For the viscous element and Kelvin model in series, the stress and strain relationship can be expressed as
And for the Kelvin model with Hooke element and viscous element in parallel, the stress and strain relationship can be expressed aswhere and are the stress and strain of the Kelvin model, and are the stress and strain of the viscous element on the left part, and are the stress and strain of the Hooke element in the Kelvin model, and and are the stress and strain of the viscous element in the Kelvin model.
The creep equation has been given previously, but the equation of elastic aftereffect without complete unloading needs to be further derived here for direct application in the future study. According to equations (5) and (6), the equations of creep curve and elastic aftereffect curve of the simplified Burgers model is presented aswhere is the constant stress during the creep stage, and is the constant stress during the stage of the elastic aftereffect.
The whole rheological process can be described by a piecewise function consisting of a creep equation (equation (7)) and an elastic aftereffect equation (equation (8)). To verify the utilization of the creep equation, first, equation (7) was used to fit the experimental data derived from stage b-c, as shown in Figure 5.

(a)

(b)

(c)
In Figure 5, the fitting black, red, and blue solid curves for (a) shale, (b) mudstone, and (c) coal have coefficients of determination (R) larger than 0.9863, proving that the creep equation of the simplified Burger’s model can perfectly characterize the creep behavior of shale, mudstone, and coal samples. Then, the elastic aftereffect equation (equation (8)) was used to fit the experimental data derived from stage d-e. Since nanoindentation elastic aftereffect curves have different initial strain values, the fitting results are separately analyzed. And Figures 6–8 give the fitting results of shale, mudstone, and coal, respectively.

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)
The experimental data of shale-1 and shale-2 can also be perfectly fitted by the elastic aftereffect equation of the simplified Burgers model even though their deformations or depths decrease at first and then increase with time, differing with other experimental results. During the period of elastic aftereffect time, deformations consist of time-dependent elastic recovery of the Kelvin model (right part in Figure 4) in the simplified Burgers model and the creep deformation of the viscous element (left part in Figure 4) in it. Thus, the different behaviors of shale-1 and shale-2 are caused by the creep strain rate exceeding the elastic recovery rate during this period.
4.3. Possible Situations of Using the Dynamic Mechanical Analysis (DMA) Method
The time-dependent characteristics of rocks during the period of the constant peak load time and the period of unloading to 10% of the peak loads were analyzed in section 4.1 and 4.2 by using the theories of rock rheology. During the period of creep time in the DMA method (Figure 1(b)), the deformation of small oscillation consists of creep deformation to increase the indentation depth and elastic aftereffect deformation causing time-dependent elastic recovery with small creep deformation. Comparison of the magnitude of creep deformations and elastic recovery deformations during this whole period can reveal the mechanism of abnormal behavior during the nanoindentation holding stage if the creep deformation and elastic aftereffect deformation happen at the same time. Taking shale-2, mudstone-2, and coal-1 as examples, the black lines with squares show the positive creep deformation, and the red lines with circles show the elastic aftereffect deformation with small creep deformation under the assumption that both of them are starting at 0 and develop simultaneously. Thus, the deformation throughout this period can be obtained by adding them to represent the creep time with small oscillations (Figure 1(a)), shown by the blue line with a triangle. Figure 9 gives additional results.

(a)

(b)

(c)
For an experiment named shale-2, the creep strain is much larger than that of elastic aftereffect strain during the whole process, which illustrates that the total deformation increases with holding time. For mudstone-2, during the first 10 seconds, the elastic aftereffect deformation is slightly larger than that of the creep deformation; the middle 130 seconds record a contrary result, and during the last 40 seconds, the elastic aftereffect deformation is slightly greater than that of the creep deformation again. And for Coal-1, the total deformation of the first 87 seconds is positive while the last 93 seconds have a negative total deformation. Consequently, the abnormal behavior that the creep deformation increases with time followed by a reduction of deformation with time in the article entitled “Abnormal behavior during nanoindentation holding stage: Characterization and explanation” published in 2019 (173) in the Journal of Petroleum Science and Engineering [32] is consistent with the characteristic of coal-1, which is just one situation of DMA method to discuss the influence of holding time on the mechanical properties of materials. For other materials, as shale and mudstone studied in this article, other situations may be encountered during creep time with small oscillation.
5. Conclusions
One of the highlights in the article is the dedicated discussion and explanation of the influence of nanoindentation elastic aftereffect on creep deformation when the dynamic mechanical analysis (DMA) method is used by adopting the method of oscillation amplification. Drawing the conclusions as follows:(1)Through analyzing the nanoindentation strain-time data, the shale has apparent rheological behavior while the abnormal situation appears that the strain increases during the stable deformation phase of the elastic aftereffect.(2)The simplified Burgers model was able to fit the nanoindentation experimental data well. After that, the reason of the different behaviors of shale could be drawn as the creep strain rate exceeding the elastic recovery rate during elastic aftereffect.(3)Under the influence of holding time, the mechanical properties of shale, mudstone, and coal are different. The elastic aftereffect deformation and creep aftereffect deformation of different materials have different performances, resulting in the small oscillation that appears in the creep stage.
Research results explain the abnormal behavior observed in an article entitled “Abnormal behavior during nanoindentation holding stage: Characterization and explanation” and published in 2019 (173) in the Journal of Petroleum Science and Engineering, which is that the elastic aftereffect deformation exceeds the creep deformation during the creep time in Figure 1(a). In turn, researchers will consider the practicability of the DMA method to discuss the influence of creep time on mechanical properties in future studies, and the quasi-static mechanical analysis (QMA) method is recommended for accuracy.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
The work presented in this paper was supported by the projects of “the Postgraduate Research and Practice Innovation Program of Jiangsu Province (KYCX21_2384),” “National Natural Science Foundation of China (52174089),” and “the Fundamental Research Funds for the Central Universities (2020ZDPY0221)”.