Abstract
Considering the heterogeneous and viscoelastic characteristics of the soil around the tunnel, the dynamic characteristics of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil are studied. To accurately describe the viscoelasticity of the soil around the circular tunnel, the fractional derivative viscoelastic model is used to describe the stress-strain relationship of the soil. To consider the heterogeneous characteristics of the soil around the circular tunnel caused by construction and in situ stress, the soil around the tunnel lining is divided into undisturbed area and disturbed area, and the disturbed area is divided into multiple thin concentric rings. A dynamic interaction model of fractional derivative viscoelastic radially inhomogeneous soil-circular tunnel lining is established. By solving the vibration of the soil in undisturbed area and each circle and considering the continuity conditions of the radial displacement and radial stress of the soil at the junction of each circle, the transfer matrix of the radial displacement and radial stress of the soil at the outer and inner boundaries of disturbed area is obtained. Based on the transfer matrix of the radial displacement and radial stress of soil in the disturbed area and considering the boundary conditions and continuity conditions of the problem, the solution of the dynamic response of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil is obtained in the frequency domain. The results show that the fractional derivative viscoelastic inhomogeneous soil-circular tunnel lining system also has resonance under the action of simple harmonic internal pressure, and the curves of the radial displacement and radial stress of soil varying with frequency obtained by the fractional derivative viscoelastic model fluctuate greatly. When studying the dynamic response of circular tunnel lining in inhomogeneous soil, the influence of viscoelasticity of the soil around the tunnel should not be ignored. The shear modulus ratio of soil in the affected area around the tunnel has a great influence on the dynamic characteristics of circular tunnel lining in fractional derivative viscoelastic radial inhomogeneous soil, and the inhomogeneous characteristics of soil should not be ignored. The greater the thickness and shear modulus of the tunnel lining, the greater the stiffness of tunnel lining and the smaller the radial stress and radial displacement of the soil.
1. Introduction
The mechanical properties of soil are important factors affecting the dynamic response of civil engineering structures. Due to the different engineering environment of the soil, there are great differences in the mechanical properties of the soil under different working conditions and in different regions. Usually, we need to focus on the influence of two properties of soil, one is the influence of the viscosity of soil, and another is the influence of the heterogeneity of soil, including radial heterogeneity and stratification. Under the action of long-term load, the material usually leads to creep or relaxation, that is, it shows viscoelastic properties. In the mechanical analysis of tunnel, pile foundation, and other structures, the influence of viscoelastic properties of soil cannot be ignored. At present, the classical viscoelastic constitutive model is mainly used to describe the viscoelastic properties of soil, but the classical viscoelastic model cannot be in good agreement with the experimental data at the initial stage of creep and relaxation [1]. To accurately fit the experimental data, the high-order differential term often has to be cancelled or the application scope of the constitutive model is reduced. Because the fractional derivative viscoelastic model can describe the mechanical behavior of soil in a wide frequency range, few experimental parameters are required to determine the model [2]. In recent years, the application of fractional derivative viscoelastic model in the field of geotechnical engineering has been paid more and more attention. Based on viscoelastic theory and theory of fractional calculus, Zhu et al. [3] proposed a fractional Kelvin-Voigt model to describe the time-varying characteristics of foundation under vertical load and gave the analytical solution of foundation settlement. By changing the fractional order of the differential operator and viscosity coefficient, the settlement time relationship can be accurately obtained, and the long-term settlement of foundation can be more accurately predicted. Wang et al. [4] introduced the theory of fractional calculus into the Kelvin-Voigt constitutive model to describe the consolidation characteristics of viscoelastic saturated soil, gave the semi analytical solution of the one-dimensional consolidation equation of the fractional derivative Kelvin-Voigt model for viscoelastic saturated soil under different variable loads, and obtained the analytical solution of effective stress and settlement in Laplace transform domain. Xiao et al. [5] used the three-dimensional fractional derivative viscoelastic model to describe the rheological properties of the soil around the pile and used the correspondence principle and Laplace transform technology to give the long-term settlement of a single pile in the fractional derivative viscoelastic soil under time-dependent load. Li et al. [6] analyzed the triaxial creep and shear test results of deep artificial frozen soil under different confining pressure and temperature, established the fractional derivative constitutive model of deep artificial frozen soil based on Nishihara model, developed the corresponding constitutive finite element program, and simulated the deep artificial freezing curtain of a coal mine, it is verified that the fractional derivative constitutive model can describe the nonattenuation creep deformation characteristics of deep artificial frozen soil under high confining pressure. Engineering construction and in situ stress usually cause the compression, disturbance, and remodeling of the soil near the engineering structure, which will lead to the change of the properties of the soil near the engineering structure, such as shear modulus and density, resulting in the nonuniform distribution of the soil around the structure. At present, considering the influence of soil heterogeneity, the research on the dynamic response of engineering structures is mainly carried out for the vibration of pile foundation [7–12].
With the construction of subway, railway, subsea tunnel, and other projects, the research on the dynamic stability of underground engineering structures is particularly important. The research on the dynamic response of underground engineering structures such as tunnels is of great significance for the seismic design and dynamic monitoring of tunnels, which has attracted extensive attention of scholars. Li [13] obtained the analytical solution of stress-displacement field caused by tunnel excavation in saturated soil in Laplace transform domain and obtained the solution in time domain with the help of inverse Laplace transform numerical technology. The results show that the stress field and displacement field are highly sensitive to the flow conditions imposed on the tunnel wall. Xie et al. [14] gave analytical solutions of stress, displacement, and pore pressure of partially closed circular tunnel lining in viscoelastic saturated soil under progressive loading in Laplace transform domain, the research shows that the local closure of boundary has a significant impact on stress, displacement, and pore pressure, and the radial displacement decreases with the increase of the viscoelastic damping coefficient of soil and the relative stiffness between lining and soil. Lu and Jeng [15] studied the dynamic response of circular tunnel in saturated porous media under moving axisymmetric annular load and obtained the numerical solutions of stress, displacement, and pore pressure. Yang et al. [16] assumed that the tunnel was a beam element, imposed the free field displacement of the horizontal layered viscoelastic foundation on the long tunnel through the viscoelastic foundation model, and analyzed the dynamic response of the long tunnel in the horizontal layered viscoelastic foundation under the action of inclined SH wave by Fourier transform. Yu and Wang [17] established the overall seismic response calculation model of tunnel surface buildings by using the finite element method, and studied the seismic response of dynamic characteristics of a subway tunnel structure under earthquake loading and the interaction characteristics between tunnel structure and ground buildings. Alielahi and Feizi [18] studied the interaction between underground tunnel and surface structure under the action of SV wave by finite element method. It was found that the influence of soil-tunnel-structure interaction on surface structure displacement and displacement depends on the characteristics of soil, underground tunnel, number of layers and input motion. Huang et al. [19] considered the slip effect at the interface between lining and ground lining, and obtained the simplified analytical solution of stress and deformation of circular composite lined tunnel under the action of in-plane shear wave, which can be used to approximately calculate the seismic response of composite lined circular tunnel. Lu et al. [20] simplified the blasting P-wave into a triangular pulse, and obtained the distribution functions of dynamic stress concentration coefficient and radial and circumferential dynamic velocity scale coefficient around the circular tunnel based on the Fourier Bessel expansion method. The effects of the rise duration and total duration, Poisson’s ratio on dynamic stress concentration coefficient, and distribution function are discussed. It can be seen that, although many achievements have been made in the current research on the dynamic response and dynamic characteristics of tunnel structure, with the rapid development of tunnel and underground engineering, the engineering environment, and geological environment of tunnel structures are more complex. Due to the influence of excavation disturbance or in situ stress, the heterogeneity of tunnel surrounding rock is common. To ensure the safety and stability of a tunnel construction and to guide the design and construction of tunnel in inhomogeneous soil, it is necessary to study the dynamic characteristics of a tunnel structure in inhomogeneous soil. Therefore, considering the heterogeneity and viscoelasticity of the soil around the tunnel, based on the fractional derivative viscoelastic model, the dynamic characteristics of the circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil are studied in the frequency domain by using the multicycle transfer matrix method and the initial parameter method.
2. Model of Viscoelastic Radially Inhomogeneous Soil-Circular Tunnel Lining
The dynamic characteristics of a wireless long and deep buried circular tunnel in infinite soil shown in Figure 1 will be investigated. There is a simple harmonic load q = q0eiωt with radial uniform distribution on the inner boundary of the tunnel lining, where ω is the frequency of the load, i is the imaginary unit, and q0 is the load amplitude. The soil around the tunnel is regarded as viscoelastic medium, and the lining of the tunnel is regarded as elastic solid. The outer radius and lining thickness of the tunnel lining are R and d, respectively. It is assumed that the viscoelastic soil in the range from R to RO around the circular tunnel is disturbed due to construction and other reasons, and the viscoelastic soil outside RO is not disturbed. The shear modulus of viscoelastic soil in undisturbed area is μO. The shear modulus of the soil in the disturbed area is given in reference [21, 22] (as shown in Figure 2).where μI is the shear modulus of the soil at the junction of tunnel lining and soil, RS is the length of the soil disturbance range, is the soil modulus ratio in the affected area, and m is a positive index, reflecting the variation law of soil along the radial direction.


3. Dynamic Equations of Fractional Derivative Viscoelastic Soil and the Solution by Multicycle Matrix Transfer Method
According to the theory of elastic dynamics, ignoring the circumferential displacement, the motion equation of soil under polar scale iswhere σr and σθ are the radial and circumferential stresses of soil, respectively, ur is the radial displacement of soil, and ρ is the density of soil mass.
Here, it is assumed that the soil around the tunnel is a viscoelastic medium, and the fractional derivative viscoelastic constitutive relationship is used to describe the relationship between stress and strain of soil [3], that is,where τ1 and τ2 are the model parameter of soil fractional derivative viscoelastic model, are the lame constant of soil, , ν is Poisson’s ratio of soil mass, is the Riemann Liouville fractional derivative, α is the order of fractional derivative, and . Riemann Liouville fractional derivative is defined as [23].where is the gamma function.
From equations (1)–(3), the motion control equation of fractional derivative viscoelastic soil expressed in displacement can be obtained as
Here, we establish the mathematical and physical model of fractional derivative viscoelastic radially inhomogeneous soil-circular tunnel lining, as shown in Figure 1. The viscoelastic soil around the tunnel is divided into disturbed area viscoelastic soil and undisturbed area viscoelastic soil. The viscoelastic soil in disturbed area is divided into n thin-layers concentric ring areas (as shown in Figure 1), and the radial thickness of each thin-layer concentric ring area is small, which can be regarded as homogeneous viscoelastic soil. The fractional derivative viscoelastic model is used to describe the stress-strain relationship of both the viscoelastic soil in the disturbed area and the viscoelastic soil in the undisturbed area. According to equation (5), the dynamic equations of the viscoelastic soil in the undisturbed area and the soil in the ith circle in the disturbed area can be established, respectively, which arewhere uOr, ρO, νO, and μO are the radial displacement, density, Poisson’s ratio, and shear modulus of viscoelastic soil in the undisturbed area, respectively, τ1O, τ2O, and siαO are the model parameters and the order of the fractional derivative of the fractional derivative viscoelastic model of viscoelastic soil in the undisturbed area, respectively. uOir, ρOi, νOi, and μOi are the radial displacement, density, Poisson’s ratio, and shear modulus of the viscoelastic soil in the ith circle of the disturbed area, τ1i, τ2i, and αi are the model parameters and the order of fractional derivative of the fractional derivative viscoelastic model of the viscoelastic soil in the ith circle in the disturbed area.
Because the circular tunnel will vibrate stably under the action of simple harmonic internal pressure, the parameters such as radial displacement and radial stress of soil in the undisturbed area and the ith circle of disturbed area meet , , , , where , , , are the amplitudes of radial displacement and radial stress of soil in the undisturbed area and the ith circle of disturbed area, respectively. Then, the parameters are substituted into equations (6) and (7), and let , , , , , , , , , , , , , . We havein which , , , and . By solving Bessel equations (8) and (9), the dimensionless radial displacements of undisturbed fractional derivative viscoelastic soil and ith circle of fractional derivative viscoelastic soil are, respectively, obtained aswhere K1(·) and I1(·) are first-order second-class and first-class deformed Bessel functions, respectively, and A1, A2, A1i, and A2i are undetermined coefficients.
Since the displacement of the undisturbed fractional derivative viscoelastic soil at infinity is zero and considering the properties of first-order first-class deformation Bessel function, we know A2 = 0, then, we have
Considering the fractional derivative viscoelastic stress-strain relationship (6), the dimensionless radial stress of the undisturbed fractional derivative viscoelastic soil is obtained:in which . Similarly, considering the fractional derivative viscoelastic stress-strain relationship (7), it can be obtained that the dimensionless radial stress of ith layer of fractional derivative viscoelastic soil in the disturbed area iswhere .
At the junction of ith layer soil and i − 1th layer soil (i.e., ), the radial displacement and radial stress of fractional derivative viscoelastic soil are, respectivelywhere K0(·) and I0(·) are zero-order second-class and first-class deformed Bessel functions, respectively. The undetermined coefficients A1i and A2i can be obtained by using the initial parameter method [24] from equations (14) and (15), and then substituted into equations (10) and (13) to obtain the following recurrence relationship:
If (16) and (17) are written in matrix form, there arein which,
Consider the continuity boundary conditions of the radial stress and the radial displacement of the soil at the outer boundary of i − 1th circle and the inner boundary of ith, we have
It can be obtained that the transfer relationship of the radial displacement and the radial stress of the soil at the inner boundary of the first circle and the outer boundary of the nth circle iswhere is the transfer matrix of radial displacement and radial stress of soil at the inner boundary of first circle and the outer boundary of the nth circle.
4. Solution of Dynamic Response of Circular Tunnel Lining in Fractional Derivative Viscoelastic Radially Inhomogeneous Soil
The lining is regarded as an elastic body, and its dynamic governing equation in polar coordinates iswhere uLr, νL, GL, and ρL are the radial displacement, Poisson’s ratio, shear modulus, and density of lining, respectively.
Similarly, under the action of simple harmonic load, the lining also makes simple harmonic vibration. The radial displacement and stress of the lining meet , , where and are the radial displacement and radial stress amplitudes of the lining, respectively, which are substituted into equation (22), at the same time, let , , , then equation (22) is dimensionless aswhere . We can obtain the radial displacement of the lining by solving Bessel equation (23)where AL1 and AL2 are undetermined coefficients. According to equation (24), the dimensionless radial stress of the lining is
If the simple harmonic internal pressure acts at the inner boundary of the lining, there iswhere , .
Considering the continuity conditions at the interface between the outer boundary of the lining and the disturbed area viscoelastic soil, and at the interface between the disturbed area viscoelastic soil and the undisturbed area viscoelastic soil, the following continuity boundary conditions can be obtained
Considering equations (12), (13), (21), and (24)–(27), the following equations can be obtained:
The undetermined coefficients A1, AL1, and AL2 can be obtained by solving equations (28)–(30), and then the expressions of the radial displacement and the radial stress of the soil at the contact surface between lining and soil can be obtained as follows:where
5. Numerical Examples and Discussion
The dynamic response characteristics of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil are discussed in the form of numerical examples. In the example, the region of fractional derivative viscoelastic inhomogeneous soil is divided into 10 circles, i.e., n = 10. If no description is given, the parameter values in the calculation example are as follows: m = 2, α1 = α2 = α3 = α4 = α5 = α6 = α7 = α8 = α9 = α10 = 0.5, αO = 0.5, β11 = β12 = β13 = β14 = β15 = β16 = β17 = β18 = β19 = β10 = 1.0, β21 = β22 = β23 = β24 = β25 = β26 = β27 = β28 = β29 = β20 = 1.0, T1O = 10, T2O = 5, = 0.33, = 0.33, = = = = = = = = = = 0.33, μ = 0.1, μL = 2.0, .
5.1. Comparative Analysis of Models
Figures 3 and 4 show the curves of the radial displacement and the radial stress of soil at the contact surface between lining and soil varying with dimensionless frequency obtained by using elastic model, classical viscoelastic model, and fractional derivative viscoelastic model, respectively. No matter what model is used, there are peaks and troughs in the curves of the radial displacement and radial stress varying with frequency. Under the action of harmonic internal pressure, there is resonance in the inhomogeneous soil-circular tunnel lining system. At the same time, it can be seen that the fractional derivative viscoelastic model can degenerate to the classical viscoelastic model and elastic model, which shows the correctness of the analysis method in this paper. When the elastic model is used, the stiffness of the system is relatively large, the radial displacement and radial stress of the soil are the smallest and the fluctuation is relatively stable. When the classical viscoelastic model and fractional derivative viscoelastic model are used, the curves of the radial displacement and radial stress of soil varying with frequency fluctuate more, and the radial displacement and radial stress of the classical viscoelastic model are the largest. Therefore, the influence of viscoelasticity of soil should not be ignored when analyzing the dynamic response of circular tunnel lining in fractional derivative viscoelastic radial inhomogeneous soil.


5.2. The Influence of the Order of the Fractional Derivative
The influence of the order of the fractional derivative on the radially inhomogeneous soil-circle tunnel lining system is shown in Figures 5 and 6. The larger the order of the fractional derivative, the greater the fluctuation of curves the radial displacement and radial stress of soil varying with the frequency. The fluctuation of curves the radial displacement and radial stress of soil varying with the frequency will be greater when the order of the fractional derivative is larger. The larger the order of the fractional derivative, the larger the peak and valley values of the curves, and the smaller the frequency corresponding to the peak and valley values. We can fit different properties of soil by changing the order of fractional derivative. It can be seen that the application range of fractional derivative viscoelastic model is wider than that of classical viscoelastic model.


5.3. The Influence of the Shear Modulus Ratio of Affected Area
Figures 7 and 8 are the influence of the shear modulus ratio of affected area on radial displacement and radial stress of soil. Obviously, when μ = 1.0, the fractional derivative viscoelastic radial heterogeneous soil can degenerate to the fractional derivative radial homogeneous soil. At low frequency, the influence of heterogeneity of soil around the tunnel on the radial displacement and radial stress of soil is greater. The smaller the shear modulus ratio of affected area, the smaller the rigidity of the system, the larger the radial displacement of soil, and the smaller the radial stress of soil. At high frequency, the influence of shear modulus ratio of affected area on the radial displacement and radial stress of soil is relatively small, but it has a great influence on the fluctuation of the curves. It can be seen that when studying the dynamic characteristics of tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil, the influence of soil heterogeneity around the tunnel should be considered.


5.4. The Influence of the Thickness of Tunnel Lining
The influence of the thickness of tunnel lining on the radial displacement and radial stress of soil is shown in Figures 9 and 10. At low frequency, the thickness of the tunnel lining has a great influence on the radial displacement of the soil, but it has little influence at high frequency. Moreover, at low frequency, the greater the thickness of tunnel lining, the smaller the radial displacement of soil. The thickness of lining has a great influence on the radial stress of the soil. The greater the thickness of lining, the smaller the radial stress of soil. This may be because the lining stiffness will be greater when the thickness of the lining is greater, so the lining will share more internal load and resulting in the smaller radial displacement and radial stress of the soil.


5.5. The Influence of the Shear Modulus Ratio of Tunnel Lining
Figures 11 and 12 show the influence of the shear modulus ratio of tunnellining on the radial displacement and radial stress of soil. The shear modulus of tunnel lining has a great influence on the dynamic characteristics of circular tunnel lining in fractional derivative radially inhomogeneous soil, but its influence will decrease with the increase of the shear modulus of lining. It can be seen that when the shear modulus of lining is large, changing the shear modulus of lining has little effect on improving the dynamic characteristics of soil-tunnel lining system.


6. Conclusion
In this paper, considering the viscoelasticity and radial heterogeneity of the soil around the tunnel, the transfer matrix at the outer boundary and inner boundary of the disturbed area is obtained by using the multicycle matrix transfer method. On this basis, the dynamic response of the circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil is studied, and the influence of relevant parameters on the dynamic characteristics of the system is numerically analyzed.
The main conclusions are as follows:(1)Whether it is an elastic model, classical viscoelastic model, or fractional derivative viscoelastic model, there is resonance in a radially inhomogeneous soil-circular tunnel lining system under simple harmonic internal pressure and the fractional derivative viscoelastic model can degenerate to classical viscoelastic model and elastic model.(2)The curves of the radial displacement and radial stress of the soil varying with frequency obtained by classical viscoelastic model and fractional derivative viscoelastic model fluctuate more than that obtained by an elastic model. When analyzing the dynamic response of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil, the influence of the viscoelasticity of soil should not be ignored.(3)By changing the order of fractional derivative, the soil with different properties under different working conditions can be fitted. The application range of fractional derivative viscoelastic model is wider than that of classical viscoelastic model.(4)The influence law of shear modulus ratio of affected area on the dynamic characteristics of circular tunnel lining in radially inhomogeneous soil has a great relationship with the frequency domain. The heterogeneous characteristics of soil around the tunnel have a great influence on the dynamic response of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil. The heterogeneous characteristics of soil around the tunnel should not be ignored.(5)The thickness of the tunnel lining and shear modulus of the tunnel lining have a great influence on the dynamic characteristics of circular tunnel lining in fractional derivative viscoelastic radially inhomogeneous soil. Increasing the thickness and shear modulus of the tunnel lining can increase the stiffness of the tunnel lining and reduce the radial stress and displacement of soil, but changing the thickness of lining at high frequency has little effect on the radial displacement of soil. When the shear modulus of the lining is large, the influence of the shear modulus of the lining on the radial displacement and radial stress of the soil will be reduced.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest to report regarding the present study.
Acknowledgments
This research was supported by National Natural Science Foundation of China (No. U1504505).