# Research on Rock Creep Characteristics Based on the Fractional Calculus Meshless Method.

1. Introduction

The rheological properties of rock are very important to the safety and long-term stability of rock engineering. However, it is difficult to test rheological properties of rock in engineering practice under a limited engineering environment. Therefore, the laboratory rheological test is always the principal method to study the rheological properties of rock. Suitable rheological constitutive models have been established to study the creep characteristics of rock by many scholars. Rheological models based on elastic elements and viscous elements have been proposed, such as Kelvin Model, Burgers Model, and Nishihara Model, and improved models of those models have also been investigated extensively [1-5]. Some scholars have carried out compression experiments on the different stress conditions of frozen rock and established a theoretical model of low-temperature rock creep, which accurately reflects the creep properties of frozen rock [6-10]. Some scholars have used the unified creep model and its corresponding creep equation to study the creep properties of rock under the condition of graded loading , which is more accurate in response to the increasing stress in the rock mass engineering. A modified Kelvin model was used to describe attenuation creep and steady creep .

Fractional calculus that is a branch of mathematical analysis has also been introduced into the study of rheological properties of rocks [13-16], and a large amount of achievements have been obtained. Some scholars have proposed an Abel dashpot on the basis of fractional order, further extended the Nishihara model, and derived a new constitutive model of viscoelastic fractional order [17, 18]. In order to describe the rheological constitutive model of salt rock, the fractional calculus theory was introduced to construct the rheological constitutive model of salt rock .

In the large-scale rock engineering, the deformation of rock is difficult to measure directly due to the complex working conditions, e.g., water seepage [20, 21], high surrounding pressure , and dynamic disturbance [23, 24]. The use of a finite element is an effective method for the simulation. In order to improve the accuracy of numerical simulation, some scholars studied the viscoelastic finite element method based on fractional calculus . However, crack propagation, rock fracture, and excessive deformation in the process of rock rheology lead to the finite element mesh deformation and failure [26, 27]. In order to solve these problems, mesh reconstruction is needed to generate new grid information, but this has a great impact on the accuracy and speed of computation. The meshless method is a new method of the numerical analysis. Based on the node information and the geometric boundary information of the computational domain, it can avoid the mesh constraints partially or thoroughly. The process of forming approximate functions in the meshless method is based on a series of discrete nodes. The problems mentioned above can be solved easily by adding the required nodes, instead of reconstructing the mesh. This method gets away from the concept of the unit, which avoids the above weakness of the finite element .

According to the above advantages, the basic equations of the generalized Kelvin model were studied based on the elementary theory of fractional calculus in this paper. Based on the data obtained from the compression creep test of the argillaceous sandstone, the meshless method was introduced to study the rock creep, whose formula was derived from the generalized Kelvin model based on fractional calculus. The results showed that the algorithm had a good stability and accuracy.

2. Meshless Method for Generalized Kelvin Model Based on Fractional Calculus

The generalized Kelvin model is presented in Figure 1, consisting of a Kelvin body and a Hooke body in series. It reflects the gradual change in deformation with time.

[k.sub.1], [[eta].sub.1] , and [[epsilon].sub.1] are the elastic modulus, viscosity coefficient, and strain of the Kelvin body, respectively; [k.sub.2] and [[epsilon].sub.2] are the elastic modulus and dependent variable, respectively. The fractional form of the stress [sigma] and strain [epsilon] relations in the model can be expressed by component parameters:

[mathematical expression not reproducible], (1)

[mathematical expression not reproducible], (2)

According to the basic definition of the Riemann-Liouville fractional calculus, the fractional order integral of Riemann-Liouville  is given as follows.

For any plural a > 0, Equation (3) can be obtained as

[mathematical expression not reproducible], (3)

where [GAMMA] ([alpha]) = [[integral].sup.[infinity].sub.0] [t.sup.[alpha]- 1][e.sup.-t] is a gamma function.

The Riemann-Liouville fractional calculus has strong singularity, which is not easy for engineering and physical modeling. Italy geophysicist Caputo proposed a fractional differential definition with weak singularity , which can be expressed as follows:

When [mathematical expression not reproducible]. (4)

Equation (2), Equation (3), and Equation (4) are combined. Set a = [k.sub.1]/[[eta].sub.1] and b = [sigma]/[[eta].sub.1], considering an initial condition [[epsilon].sub.1] = 0 when t = 0, Equation (5) can be obtained:

[d.sup.[alpha]][[epsilon].sub.1](t)/[dt.sup.[alpha]] + [alpha][[epsilon].sub.1](t) = b. (5)

Equation (5) is transformed by Laplace as

[s.aup.[alpha]][[epsilon].sub.1](s) + a[[epsilon].sub.1](s) = b/s. (6)

After transferring, it is expressed as

[[epsilon].sub.1](s) = b/[s([s.sup.[alpha]] + a)]. (7)

Equation (7) is transformed by Laplace as

[mathematical expression not reproducible], (8)

where [E.sub.[alpha],[alpha]] (z) = [[summation].sup.[infinity].sub.n=0] ([z.sup.n]/[GAMMA] (n[alpha] + [alpha])), and then the MittagLeffler equation is obtained .

Solving Equation (8), a, b are brought in:

[mathematical expression not reproducible]. (9)

Substitute Equation (9) into Equation (2), which can be written as

[mathematical expression not reproducible]. (10)

In the computational domain Q, the integral treatment of Equation (10) is carried out:

[mathematical expression not reproducible]. (11)

Then the equation can be changed as

K(t)u(t) = P, (12)

where

K(t) = [summation] [k.sub.IJ] (t), (13)

[mathematical expression not reproducible], (14)

[mathematical expression not reproducible], (15)

[mathematical expression not reproducible], (16)

where K(t) is the total stiffness matrix of the equivalent meshless method, [K.sub.IJ] (t) is the stiffness matrix between the meshless node I and the node J, l (t) is the fractional order form of the generalized Kelvin model, u(t) is the node displacement vector, and P is the node load vector.

3. Experimental Study on Compressive Creep Properties of Rock

Specimens were derived from the unweathered and nonvibratory argillaceous sandstone collected from coal mines. A standard cylindrical sample was produced with a diameter of 55 mm and a height of 100 mm, and the processed rock samples would be maintained the state of natural moisture. During the experiment, indoor temperature and humidity conditions were kept constant as much as possible to reduce the influence of external factors on the experimental results.

The three-axis compression rheology test was carried out in the State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology. The MTS815.02 electrohydraulic servo rock mechanics test system was used, as shown in Figure 2. The test system has the advantages of high loading control precision and automatic real-time acquisition of test results. It is equipped with three sets of independent closed-loop servo control systems, such as axial compression, confining pressure, and pore water pressure. It has five basic test functions, namely, uniaxial compression test, conventional triaxial compression test, true triaxial compression test, pore water pressure test, and water penetration test [31-33].

From the results of the creep test, it is shown that, at low stress levels (when the load is less than the long-term strength of the rock mass), the deformation of the specimen consists of two parts: instantaneous deformation and stable creep deformation. According to the constitutive relation of the generalized Kelvin model of fractional calculus, there are 4 unknown parameters in the model, which are two elastic moduli [k.sub.1], [[eta].sub.1] and two viscosity coefficients [k.sub.2], [alpha]. Since the constitutive equations of the model are complex, it is difficult to obtain them directly from the experimental results. According to Equation (10), the parameters of the model can be solved by the least square method, namely,

[mathematical expression not reproducible], (17)

where [[epsilon].sub.i] and [t.sub.i] are the strain and corresponding time of the experimental test, respectively, and N is the number of data. The fitting of parameters can be obtained by general nonlinear least squares fitting. When [absolute value of (x)] is small, the Mittag-Leffler function converges faster, but when [absolute value of (x)] is large, the series converges slowly and is difficult to calculate. Therefore, the asymptotic approximation of the Mittag-Leffler function is introduced.

[mathematical expression not reproducible]. (18)

With the method in the literature , the results of various parameters can be obtained.

Under different stress conditions, the rheological parameters of rock are not constant. Theoretically, there is a corresponding function relationship between the rheological parameters of rock and the stress level. However, because of the complexity in rock engineering practice and many uncertain factors in the rock mechanics test, it is not necessary to emphasize the difference of the rheological parameters under different stresses so as to improve the calculation precision. Therefore, the rheological parameters of rock under different confining pressure are assumed to be constant, and the final rheological parameters of rock can be determined according to the statistical average of the rheological parameters under various loads. The final result of parameter fitting based on the experiment is shown in Table 1.

4. Experimental Analysis and Numerical Method Verification

According to the parameters determined in Table 1, [sigma] = 6.35 (MPa) is taken. According to the actual situation of the experiment, the strain curve of 160 h is determined by using Equation (10). The result can be seen as an analytical solution of the experiment.

In order to further verify the meshless method of the generalized Kelvin model of fractional calculus, numerical analysis is carried out by using Equation (13). There are many ways to build meshless normal form functions. To facilitate the application of essential boundary conditions, the meshless radial basis interpolation (RPIM) method is used to construct the shape function. In the RPIM formula, due to the use of the Galerkin weak form, the surface force boundary condition has been naturally integrated into the discrete system equations formed during the formulation process. Because the RPIM shape function has the property of the Kronecker 8 function, its essential boundary condition can be applied as easily as the finite element method. Since the result determined by using Equation (10) is the strain and the result calculated by using Equation (13) is displacement, the numerical method determining the unique value is converted to strain for ease of comparison.

The fractional model solution determined by using Equation (10) and the numerical method solution determined by using Equation (13) are shown in Figure 6, and some results are shown in Table 2.

As shown in Figure 6 and Table 2, the results of the fractional model solution and numerical method solution are very close. The results show that the algorithm based on the generalized Kelvin model of fractional calculus is reliable and stable.

A study on creep behavior in the field of rock mechanics and mining science has been performed extensively in the literature, involving theoretical, experimental, and numerical methods. In general, the rheological behavior of rocks is typically characteristic of the viscosity, elasticity, and plasticity. Many creep models have been proposed by researchers to describe the visco-elastoplastic rheological behavior. Although the classical Burgers models are widely adopted in practical engineering due to their simplicity, the fractal derivative model shows the flexibility in the description of the rheological property of soft rock and how to interpret the physical meaning of the fractal derivative order requires further research.

5. Conclusion

(1) In this paper, the common forms of fractional calculus were introduced. The constitutive relation of the generalized Kelvin model based on fractional calculus was studied. This constitutive relation was also introduced into the meshless method, deriving a new formula for numerical calculation.

(2) Based on the elementary theory of rock rheology, using the MTS815.02 hydraulic servo rock mechanics test system, the creep test of the argillaceous sandstone was carried out, and the related data of the creep process in soft rock were obtained. Based on the generalized Kelvin model of fractional calculus, the related parameters were fitted.

(3) Using the constitutive relation of the generalized Kelvin model of fractional calculus and the obtained parameters, the analytical curves of argillaceous sandstone creep were solved.

(4) Using the newly obtained numerical method, the model was calculated and compared with the fractional model results. The results showed that the curves of the two methods were very close with a minor difference. It is shown that the meshless method based on fractional calculus had a good stability and accuracy for the generalized Kelvin model.

https://doi.org/ 10.1155/2018/1472840

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Nos. 51434006, 51774131, and 51274097) and the Open Research Fund Program of Hunan Province Key Laboratory of Safe Mining Techniques of Coal Mines (Hunan University of Science and Technology) (No. 201206).

References

 Y. J. Chen, C. L. Pan, P. Cao, and W. X. Wang, "A new mechanical model for soft rock rheology," Rock and Soil Mechanics, vol. 24, no. 2, pp. 209-214, 2003.

 Y. Zhao, Y. Wang, W. Wang, W. Wan, and J. Tang, "Modeling of non-linear rheological behavior of hard rock using triaxial rheological experiment," International Journal of Rock Mechanics and Mining Sciences, vol. 93, pp. 66-75, 2017.

 S. G. Cao, J. Bian, and P. Li, "Comparison research of creep experiments and the analysis of the theoretical model of soft rocks," Journal of Chongqing University, vol. 25, no. 7, pp. 96-98, 2002.

 Y. L. Zhao, L. Y. Zhang, W. J. Wang, W. Wan, and W. H. Ma, "Separation of elastoviscoplastic strains of rock and a nonlinear creep model," International Journal of Geomechanics, vol. 18, no. 1, article 04017129, 2018.

 W. Y. Xu, J. W. Zhou, S. Q. Yang, and C. Shi, "Study on creep damage constitutive relation of greenschist specimen," Chinese Journal of Rock Mechanics and Engineering, vol. 25, pp. 3093-3097, 2006.

 R. L. Shan, L. W. Song, D. Y. Li, B. L. Huang, N. Liu, and W. F. Zhao, "Study of nonlinear creep model of frozen red sandstone," Rock and Soil Mechanics, vol. 35, no. 6, pp. 1541-1546, 2016.

 J. F. Zhu, Z. M. Yao, and J. H. Chen, "The uniaxial mechanical performance test and creep model of artificial frozen deep soft rock," Mechanics in Engineering, vol. 38, no. 6, pp. 651-657, 2016.

 J. Zhu, Y. Xu, and Y. Q. Qi, "Research on the creep characteristics of frozen red--Bed Soft Rock," Journal of Anhui University of Science and Technology (Natural Science), vol. 33, no. 4, pp. 26-31, 2013.

 X. J. Mao, Z. M. Yao, and F. Mao, "Kelvin creep model of artificial frozen deep soft rock based on particle swarm plasticity strengthening," Chinese Journal of Underground Space and Engineering, vol. 10, no. 1, pp. 1602-1605, 2014.

 D. W. Li, R. H. Wang, and J. H. Fan, "Nonlinear rheological model for frozen soft rock during cretaceous period," Chinese Journal of Geotechnical Engineering, vol. 33, no. 3, pp. 398-403, 2011.

 Z. T. Zhang and J. J. Luo, "Study on creep properties of rock under step load," Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 2, pp. 218-222, 2004.

 Y. G. Kang and X. E. Zhang, "An improved constitutive model for rock creep," Rock and Soil Mechanics, vol. 35, no. 4, pp. 1049-1055, 2014.

 R. L. Bagley and P. J. Torvik, "Fractional calculus-a different approach to the analysis of viscoelastically damped structures," AIAA Journal, vol. 21, no. 5, pp. 741-748, 1983.

 R. L. Bagley and P. J. Torvik, "Fractional calculus in the transient analysis of viscoelastically damped structures," AIAA Journal, vol. 23, no. 6, pp. 918-925, 1985.

 R. C. Koeller, "Application of fractional calculus to the theory of viscoelasticity," Journal of Applied Mechanics, vol. 51, no. 2, pp. 299-307, 1984.

 A. D. Drozdov, "Fractional differential models in finite viscoelasticity," Acta Mechanica, vol. 124, no. 1-4, pp. 155-180, 1997.

 V. Kiryakova and A. Sauabi, "Explicit solutions to hyper-Bessel integral equations of second kind," Computers and Mathematics with Applications, vol. 37, no. 1, pp. 75-86,1999.

 M. D. Paola, A. Pirrotta, and A. Valenza, "Visco-elastic behavior through fractional calculus: an easier method for best fitting experimental results," Mechanics of Materials, vol. 43, no. 12, pp. 799-806, 2011.

 H. W. Zhou, C. P. Wang, Z. Q. Duan, M. Zhang, and J. F. Liu, "Time-based fractional derivative approach to creep constitutive model of salt rock," Scientia Sinica (Physica, Mechanica & Astronomica, vol. 42, no. 3, pp. 310-318, 2012.

 Y. L. Zhao, S. L. Luo, Y. X. Wang, W. Wang, L. Zhang, and W. Wan, "Numerical analysis of karst cave water inrush and a criterion for the setting width of water-resistant rock pillar," Mine Water and Environment, vol. 36, no. 4, pp. 508-519, 2017.

 D. Ma, Z. L. Zhou, J. Y. Wu, Q. Li, and H. B. Bai, "Grain size distribution effect on the hydraulic properties of disintegrated coal mixtures," Energies, vol. 10, no. 5, p. 612, 2017.

 X. Fan, K. Li, H. Lai, Y. Xie, R. Cao, and J. Zheng, "Internal stress distribution and cracking around flaws and openings of rock block under uniaxial compression: a particle mechanics approach," Computers and Geotechnics, vol. 102, no. 10, pp. 28-38, 2018.

 D. Ma, Q. Li, M. R. Hall, and Y. Wu, "Experimental investigation of stress rate and grain size on gas seepage characteristics of granular coal," Energies, vol. 10, no. 4, p. 527, 2017.

 Z. L. Zhou, X. Cai, D. Ma, W. Z. Cao, C. Lu, and Z. Jing, "Effects of water content on fracture and mechanical behavior of sandstone with a low clay mineral content," Engineering Fracture Mechanics, vol. 193, pp. 47-65, 2018.

 H. Yin and N. Chen, "Finite element method for viscoelastic fractional derivative model," Chinese Journal of Computational Mechanics, vol. 29, no. 6, pp. 966-971, 2012.

 X. Zhang, Y. Liu, and S. Mang, "Meshfree methods and their applications," Advances in mechanics, vol. 39, no. 1, pp. 1-26, 2009.

 Z. L. Zhou, X. Cai, D. Ma, C. Lu, S. F. Wang, and L. H. Tan, "Dynamic tensile properties of sandstone subjected to wetting and drying cycles," Construction and Building Materials, vol. 82, pp. 215-232, 2018.

 Y. Zhang, Z. K. Shong, and W. M. Lu, "Research progress and application of meshless method," Chinese Journal of Computational Mechanics, vol. 20, no. 6, pp. 730-742, 2003.

 G. S. F. Frederico and D. F. M. Torres, "Fractional isoperimetric Noether's theorem in the Riemann-Liouville sense," Reports On Mathematical Physics, vol. 71, no. 3, pp. 291-304, 2013.

 I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, CA, USA, 1999.

 Y. L. Zhao, L. Y. Zhang, W. J. Wang, J. Z. Tang, H. Lin, and W. Wan, "Transient pulse test and morphological analysis of single rock fractures," International Journal of Rock Mechanics and Mining Sciences, vol. 91, pp. 139-154, 2017.

 D. Ma, X. Cai, Z. L. Zhou, and X. B. Li, "Experimental investigation on hydraulic properties of granular sandstone and mudstone mixtures," Geofluids, vol. 2018, Article ID 9216578, 13 pages, 2018.

 Y. X. Wang, P. P. Guo, W. X. Ren et al., "Laboratory investigation on strength characteristics of expansive soil treated with jute fiber reinforcement," International Journal of Geomechanics, vol. 17, no. 11, article 04017101, 2017.

Gang Peng [ID], (1,2,3) Zhanqing Chen [ID], (1) and Jiarui Chen [ID] (4)

(1) State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China

(2) Hunan Provincial Key Laboratory of Safe Mining Techniques of Coal Mines, Hunan University of Science and Technology, Xiangtan, Hunan 411201, China

(3) Work Safety Key Lab on Prevention and Control of Gas and Roof Disasters for Southern Coal Mines, Hunan University of Science and Technology, Xiangtan, Hunan 411201, China

(4) HuaiYin Institute of Technology, Huai'an, Jiangsu 223003, China

Correspondence should be addressed to Gang Peng; gpengcn@126.com

Received 7 June 2018; Accepted 28 August 2018; Published 26 September 2018

Caption: FIGURE 1: Generalized Kelvin model.

Caption: FIGURE 2: MTS815.02 electrohydraulic servo rock mechanics test system.

Caption: FIGURE 3: Failure specimens after the test.

Caption: FIGURE 4: Creep curve: the confining pressure is 0 MPa.

Caption: FIGURE 5: Creep curve: the confining pressure is 10 MPa.

Caption: FIGURE 6: Fractional model solution and numerical method solution.
```TABLE 1: Parameters determined by fitting analysis based on creep
tests.

Confining       [k.sub.1]      [[eta].sub.1]      [k.sub.2]   [alpha]
pressure          (GPa)     ([10.sup.14])Pa x s     (GPa)
(MPa)

0                 4.43             2.42             6.24       0.90
10                4.64             2.06             5.60       0.91
Average value     4.54             2.24             5.92       0.91

TABLE 2: Partial data of the fractional model solution and numerical
method solution.

Time (h)                     0        40       80      120      160

Fractional model           10.727   22.255   24.047   24.511   24.680
([epsilon]([10.sup.-3]))
Numerical method           10.729   22.259   24.052   24.516   24.685
([epsilon]([10.sup.-3]))
```
Title Annotation: Printer friendly Cite/link Email Feedback Research Article Peng, Gang; Chen, Zhanqing; Chen, Jiarui Advances in Civil Engineering Report 9CHIN Jan 1, 2018 3734 A New Mining Scheme for Hanging-Wall Ore-Body during the Transition from Open Pit to Underground Mining: A Numerical Study. Coupled Geomechanical-Flow Assessment of C[O.sub.2] Leakage through Heterogeneous Caprock during CCS. Calculus Calculus (Mathematics) Creep Materials Rock mechanics Testing equipment