The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-e...The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.展开更多
Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with th...Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.展开更多
In this paper, a series of experiments, including atomic force microscope （AFM）, environmental scanning electron microscope （ESEM）, and core displacement tests were conducted to investigate the effect of polymer s...In this paper, a series of experiments, including atomic force microscope （AFM）, environmental scanning electron microscope （ESEM）, and core displacement tests were conducted to investigate the effect of polymer solution structure on solution properties and oil displacement efficiency. The results show that in the HPAM solution polymer coils were formed and then aggregated into a loose structure, while the HAP2010 solution formed a strong network structure, which would significantly improve the solution viscosity and flow resistance so as to upgrade the capacity of piston-like oil displacement in highly permeable porous media. Meanwhile, the retention of the HAP2010 solution at pore throats were also enhanced, which could reduce water production during subsequent water flooding and enlarge the swept volume during polymer flooding. Therefore, enhancing the interaction among polymer molecules is an effective way to improve the displacement efficiency of polymer solutions in heavy oil reservoirs with high permeability.展开更多
In this paper the method of reciprocal theorem is extended to find solutions of three-D problems of elasticity.First we give the basic solution of the cube with six surfaces fixed as the basic system and then using th...In this paper the method of reciprocal theorem is extended to find solutions of three-D problems of elasticity.First we give the basic solution of the cube with six surfaces fixed as the basic system and then using the reciprocal theorem between the basic system acted on by unit concentrated loads and the actual system with prescribed surface displacements, we find displacement solution of the actual system.展开更多
Single-phase concentrated solid solution alloys(SP-CSAs),including high-entropy alloys,have received extensive attention due to their excellent irradiation resistance.In this work,displacement cascade simulations are ...Single-phase concentrated solid solution alloys(SP-CSAs),including high-entropy alloys,have received extensive attention due to their excellent irradiation resistance.In this work,displacement cascade simulations are conducted using the molecular dynamics method to study the evolution of defects in Ni-based SP-CSAs.Compared with pure Ni,the NiCr,NiCo,and NiCu alloys exhibit a larger number of Frankel pairs(FPs)in the thermal peak stage,but a smaller number of surviving FPs.However,the NiFe alloy displays the opposite phenomenon.To explain these different observations for NiFe and other alloys,the formation energy and migration energy of interstitials/vacancies are calculated.In the NiFe alloy,both the formation energy and migration energy barrier are higher.On the other hand,in NiCr and other alloys,the formation energy of interstitials/vacancies is lower,as is the migration energy barrier of interstitials.The energy analysis agrees well with previous observations.The present work provides new insights into the mechanism behind the irradiation resistance of binary Ni-based SP-CSAs.展开更多
Recently, a new （2＋1）-dimensional shallow water wave system, the （2＋1）-dlmenslonal displacement shallow water wave system （2DDSWWS）, was constructed by applying the variational principle of the analytic mechan...Recently, a new （2＋1）-dimensional shallow water wave system, the （2＋1）-dlmenslonal displacement shallow water wave system （2DDSWWS）, was constructed by applying the variational principle of the analytic mechanics in the Lagrange coordinates. The disadvantage is that fluid viscidity is not considered in the 2DDSWWS, which is the same as the famous Kadomtsev-Petviashvili equation and Korteweg-de Vries equation. Applying dimensional analysis, we modify the 2DDSWWS and add the term related to the fluid viscidity to the 2DDSWWS. The approximate similarity solutions of the modified 2DDSWWS （M2DDSWWS） is studied and four similarity solutions are obtained. For the perfect fluids, the coefficient of kinematic viscosity is zero, then the M2DDSWWS will degenerate to the 2DDSWWS.展开更多
The sequence characteristics and focal mechanism solution of the Jiashi, Xinjiang strong earthquake swarm are analyzed and studied in this paper. The result shows that before the M S=6.6 earthquake, value h o...The sequence characteristics and focal mechanism solution of the Jiashi, Xinjiang strong earthquake swarm are analyzed and studied in this paper. The result shows that before the M S=6.6 earthquake, value h of sequence frequency attenuation coefficient was less than 1, then value h was more than 1. Before occurrence of M S6.0 earthquakes the energy is released either in a continuously strengthened way or a sharply strengthened way, and before M S5.0 earthquakes the sequence frequency shows calm. The study on the focal mechanism solution of the strong earthquake swarm shows that the source faults are mainly in a right lateral, strike slip way and the faults have characteristics of tensor shear.展开更多
The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak stren...The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock’s characteristics of constant modulus and double moduli.Then,the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli,respectively,considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship.Finally,the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics,the strength drop,the stiffness degradation,and the dilatancy.The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression.The coefficients of the expression are related to the stress field of the original rock,the elastic constant of the surrounding rock,the strength parameters,and the dilatancy angle.The strength drop,the stiffness degradation,and the dilatancy effects all have effects on the displacement.The effects can be characterized by quantitative relationships.展开更多
In this paper, the author obtains the more general displacement solutions for the isotropic plane elasticity problems. The general solution obtained in ref. [ 1 ] is merely the particular case of this paper, In compar...In this paper, the author obtains the more general displacement solutions for the isotropic plane elasticity problems. The general solution obtained in ref. [ 1 ] is merely the particular case of this paper, In comparison with ref. [1], the general solutions of this paper contain more arbitrary constants. Thus they may satisfy more boundary conditions.展开更多
It is demonstrated that when tension leg platform （TLP） moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP,...It is demonstrated that when tension leg platform （TLP） moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP, The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theo- retical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displace- ment, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force, Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.展开更多
The focal mechanism parameters of small earthquakes are determined by the maximum velocity and displacement amplitude ratio of the direct -and -waves recorded by digital stations. The displacement is obtained from...The focal mechanism parameters of small earthquakes are determined by the maximum velocity and displacement amplitude ratio of the direct -and -waves recorded by digital stations. The displacement is obtained from the velocity by emulation, and the two results are compared and analyzed. Results of theoretical analysis and practical measurement indicate that the two results of velocity and displacement are consistent, and it is feasible that the maximum displacement amplitude ratio be replaced by the maximum velocity amplitude ratio of the direct -and -waves recorded by regional seismic networks when determining focal mechanism solutions of small earthquakes.展开更多
The unified displacement function(UDF)is presented to describe the deformation behaviours of the tunnel profile along with time under the surface slope condition.Based on the discrete Fourier method,the third-order UD...The unified displacement function(UDF)is presented to describe the deformation behaviours of the tunnel profile along with time under the surface slope condition.Based on the discrete Fourier method,the third-order UDF in the physical plane is expanded to the Laurent series in the complex variable plane.The complex variable method is employed to derive the elastic analytical solution of stra-tum displacement,when the third-order UDF is taken as the displacement boundary condition of tunnel cross-section(DBCTC).The proposed elastic solution agrees well with the results of the finite element method for the consistent model,which verifies the correctness of the proposed analytical solution.Combining the corresponding principle and fractional Generalized Kelvin viscoelastic constitutive model,the fractional viscoelastic solution under the surface slope condition is determined.The time effect of stratum displacement is presented in two aspects:time-dependent DBCTC and time-dependent material parameters.The parameter analysis is performed to investigate influences of deformation modes of the third-order UDF,slope angle,tunnel radius and fractional order on the time effect of stratum vertical and horizontal displacement.展开更多
基金Projects(52378411,52208404)supported by the National Natural Science Foundation of China。
文摘The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.
基金We would like to acknowledge all the reviewers and editors and the sponsorship of National Natural Science Foundation of China(42030103)the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(2021QNLM020001-6)the Laoshan National Laboratory of Science and Technology Foundation(LSKJ202203400).
文摘Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.
基金supported by the National Science and Technology Major Project (2011ZX05024-004)National High Technology Research and Development Program of China (863 Program: 2007AA090701-3)
文摘In this paper, a series of experiments, including atomic force microscope （AFM）, environmental scanning electron microscope （ESEM）, and core displacement tests were conducted to investigate the effect of polymer solution structure on solution properties and oil displacement efficiency. The results show that in the HPAM solution polymer coils were formed and then aggregated into a loose structure, while the HAP2010 solution formed a strong network structure, which would significantly improve the solution viscosity and flow resistance so as to upgrade the capacity of piston-like oil displacement in highly permeable porous media. Meanwhile, the retention of the HAP2010 solution at pore throats were also enhanced, which could reduce water production during subsequent water flooding and enlarge the swept volume during polymer flooding. Therefore, enhancing the interaction among polymer molecules is an effective way to improve the displacement efficiency of polymer solutions in heavy oil reservoirs with high permeability.
文摘In this paper the method of reciprocal theorem is extended to find solutions of three-D problems of elasticity.First we give the basic solution of the cube with six surfaces fixed as the basic system and then using the reciprocal theorem between the basic system acted on by unit concentrated loads and the actual system with prescribed surface displacements, we find displacement solution of the actual system.
基金supported by the National Natural Science Foundation of China(12232008,12072211)Foundation of Key laboratory(2022JCJQLB05703)Sichuan Province Science and Technology Project(2023NSFSC0914,2020JDJQ0029).
文摘Single-phase concentrated solid solution alloys(SP-CSAs),including high-entropy alloys,have received extensive attention due to their excellent irradiation resistance.In this work,displacement cascade simulations are conducted using the molecular dynamics method to study the evolution of defects in Ni-based SP-CSAs.Compared with pure Ni,the NiCr,NiCo,and NiCu alloys exhibit a larger number of Frankel pairs(FPs)in the thermal peak stage,but a smaller number of surviving FPs.However,the NiFe alloy displays the opposite phenomenon.To explain these different observations for NiFe and other alloys,the formation energy and migration energy of interstitials/vacancies are calculated.In the NiFe alloy,both the formation energy and migration energy barrier are higher.On the other hand,in NiCr and other alloys,the formation energy of interstitials/vacancies is lower,as is the migration energy barrier of interstitials.The energy analysis agrees well with previous observations.The present work provides new insights into the mechanism behind the irradiation resistance of binary Ni-based SP-CSAs.
基金Project supported by the Natural Science Foundation of Guangdong Province of China (Grant No.10452840301004616)the National Natural Science Foundation of China (Grant No.61001018)the Scientific Research Foundation for the Doctors of University of Electronic Science and Technology of China Zhongshan Institute (Grant No.408YKQ09)
文摘Recently, a new （2＋1）-dimensional shallow water wave system, the （2＋1）-dlmenslonal displacement shallow water wave system （2DDSWWS）, was constructed by applying the variational principle of the analytic mechanics in the Lagrange coordinates. The disadvantage is that fluid viscidity is not considered in the 2DDSWWS, which is the same as the famous Kadomtsev-Petviashvili equation and Korteweg-de Vries equation. Applying dimensional analysis, we modify the 2DDSWWS and add the term related to the fluid viscidity to the 2DDSWWS. The approximate similarity solutions of the modified 2DDSWWS （M2DDSWWS） is studied and four similarity solutions are obtained. For the perfect fluids, the coefficient of kinematic viscosity is zero, then the M2DDSWWS will degenerate to the 2DDSWWS.
文摘The sequence characteristics and focal mechanism solution of the Jiashi, Xinjiang strong earthquake swarm are analyzed and studied in this paper. The result shows that before the M S=6.6 earthquake, value h of sequence frequency attenuation coefficient was less than 1, then value h was more than 1. Before occurrence of M S6.0 earthquakes the energy is released either in a continuously strengthened way or a sharply strengthened way, and before M S5.0 earthquakes the sequence frequency shows calm. The study on the focal mechanism solution of the strong earthquake swarm shows that the source faults are mainly in a right lateral, strike slip way and the faults have characteristics of tensor shear.
基金Project supported by the National Natural Science Foundation of China and Shandong Province Joint Program(No.U1806209)the National Natural Science Foundation of China(Nos.51774196 and 51774194)and Shandong University of Science and Technology(SDUST)Research Fund(No.2019TDJH101)。
文摘The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock’s characteristics of constant modulus and double moduli.Then,the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli,respectively,considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship.Finally,the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics,the strength drop,the stiffness degradation,and the dilatancy.The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression.The coefficients of the expression are related to the stress field of the original rock,the elastic constant of the surrounding rock,the strength parameters,and the dilatancy angle.The strength drop,the stiffness degradation,and the dilatancy effects all have effects on the displacement.The effects can be characterized by quantitative relationships.
文摘In this paper, the author obtains the more general displacement solutions for the isotropic plane elasticity problems. The general solution obtained in ref. [ 1 ] is merely the particular case of this paper, In comparison with ref. [1], the general solutions of this paper contain more arbitrary constants. Thus they may satisfy more boundary conditions.
基金Project supported by "Creativeness Project of the Tenth Five-Year Plan" of Chinese Academy of Sciences (No.KJCX2-SW-L03)the National High-Tech Research and Development Program of China (863 Program) (No.2004AA617010)
文摘It is demonstrated that when tension leg platform （TLP） moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP, The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theo- retical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displace- ment, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force, Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.
基金the"Strong Earthquake Monitoring and Prediction Research"for the 11th"Five-Year Plan"under the grant No.2006BA(01B02-01-01)subspecial subject of 2006BAC01 B02-01-01the State Key Basic Research Development Programme (2004CB418404),China
文摘The focal mechanism parameters of small earthquakes are determined by the maximum velocity and displacement amplitude ratio of the direct -and -waves recorded by digital stations. The displacement is obtained from the velocity by emulation, and the two results are compared and analyzed. Results of theoretical analysis and practical measurement indicate that the two results of velocity and displacement are consistent, and it is feasible that the maximum displacement amplitude ratio be replaced by the maximum velocity amplitude ratio of the direct -and -waves recorded by regional seismic networks when determining focal mechanism solutions of small earthquakes.
基金the financial supports from the National Natural Science Foundation of China(Grant No.52025084)the Beijing Natural Science Foundation,China(Grant No.8212007).
文摘The unified displacement function(UDF)is presented to describe the deformation behaviours of the tunnel profile along with time under the surface slope condition.Based on the discrete Fourier method,the third-order UDF in the physical plane is expanded to the Laurent series in the complex variable plane.The complex variable method is employed to derive the elastic analytical solution of stra-tum displacement,when the third-order UDF is taken as the displacement boundary condition of tunnel cross-section(DBCTC).The proposed elastic solution agrees well with the results of the finite element method for the consistent model,which verifies the correctness of the proposed analytical solution.Combining the corresponding principle and fractional Generalized Kelvin viscoelastic constitutive model,the fractional viscoelastic solution under the surface slope condition is determined.The time effect of stratum displacement is presented in two aspects:time-dependent DBCTC and time-dependent material parameters.The parameter analysis is performed to investigate influences of deformation modes of the third-order UDF,slope angle,tunnel radius and fractional order on the time effect of stratum vertical and horizontal displacement.
