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

State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China Hunan Provincial Key Laboratory of Safe Mining Techniques of Coal Mines, Hunan University of Science and Technology, Xiangtan, Hunan 411201, China 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 HuaiYin Institute of Technology, Huai’an, Jiangsu 223003, China


Introduction
e 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.
erefore, 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][2][3][4][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][7][8][9][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 [11], 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 [12].
Fractional calculus that is a branch of mathematical analysis has also been introduced into the study of rheological properties of rocks [13][14][15][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 [19].
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 [22], and dynamic disturbance [23,24].e 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 [25].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.e 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.e process of forming approximate functions in the meshless method is based on a series of discrete nodes.e problems mentioned above can be solved easily by adding the required nodes, instead of reconstructing the mesh. is method gets away from the concept of the unit, which avoids the above weakness of the finite element [28].
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.e results showed that the algorithm had a good stability and accuracy.

Meshless Method for Generalized Kelvin Model Based on Fractional Calculus
e 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 1 , η 1 , and ε 1 are the elastic modulus, viscosity coefficient, and strain of the Kelvin body, respectively; k 2 and ε 2 are the elastic modulus and dependent variable, respectively.e fractional form of the stress σ and strain ε relations in the model can be expressed by component parameters: Equation ( 2) can be obtained by Equation (1): According to the basic definition of the Riemann-Liouville fractional calculus, the fractional order integral of Riemann-Liouville [29] is given as follows.
Solving Equation ( 8), a, b are brought in: Substitute Equation ( 9) into Equation ( 2), which can be written as In the computational domain Ω, the integral treatment of Equation ( 10) is carried out: en the equation can be changed as where where K(t) is the total stiffness matrix of the equivalent meshless method, K 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.

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.e 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.
e MTS815.02electrohydraulic servo rock mechanics test system was used, as shown in Figure 2. e 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][32][33].
For rock rheological tests, there are usually two loading methods, namely, separate loading method and multistage loading method.
e multistage loading method was adopted in this test.e rheological tests were carried out with a confining pressure of 0 MPa and 10 MPa, respectively, the axial pressure of the specimen was changed, and the multistage loading was carried out.e test was terminated when the creep deformation at all levels of load tended to be stable or had a constant steady-state creeping state, and then it was transferred to the next level of loading.e failure samples are shown in Figure 3.No RS01 was tested under confining pressure 0 MPa, and No RS02 was tested under confining pressure 10 MPa. e results of the experiment are presented in Figures 4 and 5.
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 1 , η 1 and two viscosity coefficients k 2 , α.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, where ε i and t i are the strain and corresponding time of the experimental test, respectively, and N is the number of data.e fitting of parameters can be obtained by general nonlinear least squares fitting.When |x| is small, the Mittag-Leffler function converges faster, but when |x| is large, the series converges slowly and is difficult to calculate.erefore, the asymptotic approximation of the Mittag-Leffler function is introduced.
With the method in the literature [19], the results of various parameters can be obtained.
Under different stress conditions, the rheological parameters of rock are not constant.
eoretically, 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.erefore, 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.e final result of parameter fitting based on the experiment is shown in Table 1.

Experimental Analysis and Numerical Method Verification
According to the parameters determined in Table 1, σ � 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).e 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).
ere 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 δ 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.
e 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.e 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.Advances in Civil Engineering (4) Using the newly obtained numerical method, the model was calculated and compared with the fractional model results.e 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.Advances in Civil Engineering

4
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.

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