In this paper,we develop a multi-scalar auxiliary variables(MSAV)scheme for the Cahn-Hilliard Magnetohydrodynamics system by introducing two scalar auxiliary variables(SAV).This scheme is linear,fully decoupled and un...In this paper,we develop a multi-scalar auxiliary variables(MSAV)scheme for the Cahn-Hilliard Magnetohydrodynamics system by introducing two scalar auxiliary variables(SAV).This scheme is linear,fully decoupled and unconditionally stable in energy.Subsequently,we provide a detailed implementation procedure for full decoupling.Thus,at each time step,only a series of linear differential equations with constant coefficients need to be solved.To validate the effectiveness of our approach,we conduct an error analysis for this first-order scheme.Finally,some numerical experiments are provided to verify the energy dissipation of the system and the convergence of the proposed approach.展开更多
In the practical slope engineering,the stability of lower sliding mass(region A)with back tensile cracks of the jointed rock slope attracts more attentions,but the upper rock mass(region B)may also be unstable.Therefo...In the practical slope engineering,the stability of lower sliding mass(region A)with back tensile cracks of the jointed rock slope attracts more attentions,but the upper rock mass(region B)may also be unstable.Therefore,in this study,based on the stepped failure mode of bedding jointed rock slopes,considering the influence of the upper rock mass on the lower stepped sliding mass,the improved failure model for analyzing the interaction force(F_(AB))between two regions is constructed,and the safety factors(F_(S))of two regions and whole region are derived.In addition,this paper proposes a method to determine the existence of F_(AB) using their respective acceleration values(a_(A) and a_(B))when regions A and B are unstable.The influences of key parameters on two regions and the whole region are analyzed.The results show that the variation of the F_(AB) and F_(S) of two regions can be obtained accurately based on the improved failure model.The accuracy of the improved failure model is verified by comparative analysis.The research results can explain the interaction mechanism of two regions and the natural phenomenon of slope failure caused by the development of cracks.展开更多
Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensi...Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.展开更多
The traditional stability analysis method of geogrid reinforced slopes does not consider the effect of lateral swelling,so it is not applicable to reinforced expansive soil slopes.This paper reports a new stability an...The traditional stability analysis method of geogrid reinforced slopes does not consider the effect of lateral swelling,so it is not applicable to reinforced expansive soil slopes.This paper reports a new stability analysis method for geogrid reinforced expansive soil slopes.The additional pullout force of the free zone due to the lateral swelling and the anti-pullout safety factor of each geogrid layer were obtained by ensuring the overall stability of the reinforced slope.The optimum design was carried out to treat an expansive soil cut slope in Hubei Province,China,by changing the spacing and length of geogrid reinforcement.Calculation results show that the additional pullout force caused by lateral swelling has a great influence on the anti-pullout stability of geogrids,and the local stability of the reinforced slope will be overestimated if the swelling effect of soil in the free zone is not considered.展开更多
Rainfall infiltration depth and mode can severely influence slope stability.With the sustained rainfall,the influenced region of slope gradually expands.By using the Green-Ampt model to the soil slope,infiltration reg...Rainfall infiltration depth and mode can severely influence slope stability.With the sustained rainfall,the influenced region of slope gradually expands.By using the Green-Ampt model to the soil slope,infiltration regulation was discussed under sustained and small intensity rainfall.And the infiltration rate of unsaturated soil was proposed according to the saturated infiltration theory.Because of the changing of initial moisture content in depth of slope,the saturated or unsaturated infiltration rate and depth could also be changeable with the sustained rainfall infiltration.Based on the principle of strength reduction,the calculation model of slope safety factor was established under different initial moisture contents and infiltration modes.Then,the slope stability was quantitatively analyzed through software FLAC3D.The calculation results of soil slope engineering show that there is a shorter period for slope stability under different initial moisture contents and unsaturated infiltration ways at the slope wetting front.The stability period of slope is 33.3%according to different initial moisture contents of wetting front less than that of the same initial moisture content of wetting front.And the slope is easier to fail under the unsaturated infiltration.The results agree well with the actual situation under sustained and small intensity rainfall.展开更多
A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modific...A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.展开更多
To investigate the stability of rock mass in high geostress underground powerhouse caverns subjected to excavation,a microseismic(MS)monitoring system was established and the discrete element method(DEM)-based numeric...To investigate the stability of rock mass in high geostress underground powerhouse caverns subjected to excavation,a microseismic(MS)monitoring system was established and the discrete element method(DEM)-based numerical simulation was carried out.The tempo-spatial damage characteristics of rock mass were analyzed.The evolution laws of MS source parameters during the formation of a rock collapse controlled by high geostress and geological structure were investigated.Additionally,a three-dimensional DEM model of the underground powerhouse caverns was built to reveal the deformation characteristics of rock mass.The results indicated that the MS events induced by excavation of high geostress underground powerhouse caverns occurred frequently.The large-stake crown of the main powerhouse was the main damage area.Prior to the rock collapse,the MS event count and accumulated energy release increased rapidly,while the apparent stress sharply increased and then decreased.The amount and proportion of shear and mixed MS events remarkably increased.The maximum displacement was generally located near the spandrel areas.The MS monitoring data and numerical simulation were in good agreement,which can provide significant references for damage evaluation and disaster forecasting in high geostress underground powerhouse caverns.展开更多
The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a ...The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a detailed calculation procedure and a definition of FOS for slope stability was developed based on the understanding of SRM. When constructing the new definition of FOS, efforts were made to make sure that it has concise physical meanings and fully reflects the shear strength of the slope. Two examples, slopes A and B with the slope angles of 63° and 34° respectively, were given to verify the method presented. It is found that, for these two slopes, the FOSs from original strength reduction method are respectively 1.5% and 38% higher than those from double reduction method. It is also found that the double reduction method predicts a deeper potential slide line and a larger slide mass. These results show that on one hand, the double reduction method is comparative to the traditional methods and is reasonable, and on the other hand, the original strength reduction method may overestimate the safety of a slope. The method presented is advised to be considered as an additional option in the practical slope stability evaluations although more useful experience is required.展开更多
Based on strength reduction theory,the stability numbers of shallow tunnels were investigated within the framework of upper and lower bound theorems of limit analysis. Stability solutions taking into account of water ...Based on strength reduction theory,the stability numbers of shallow tunnels were investigated within the framework of upper and lower bound theorems of limit analysis. Stability solutions taking into account of water seepage were presented and compared with those without considering seepage. The comparisons indicate that the maximum difference does not exceed 3.7%,which proves the present method credible. The results show that stability numbers of shallow tunnels considering seepage are much less than those without considering seepage,and that the difference of stability numbers between considering seepage and without considering seepage increase with increasing the depth ratio. The stability numbers decrease with increasing permeability coefficient and groundwater depth. Seepage has significant effects on the stability numbers of shallow tunnels.展开更多
In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of th...In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of the particle micro parameters on the macroscopic characteristics of cohesive soil and calibrate the micro parameters of DEM model on this basis. Embankment slope stability analysis was carried out by strength reduction and gravity increase method, it is shown that the safety factor obtained by strength reduction method is more conservative, and the arc-shaped feature of the sliding surface under the gravity increase method is more obvious. Throughout the progressive failure process, the failure trends, maximum displacements, and velocity changes obtained by the two methods were consistent. When slope was destroyed, the upper part was cracked, the middle part was sheared, and the lower part was destroyed by extrusion. The conclusions of this paper can be applied to the safety factor calculation of cohesive soil slopes and the analysis of the instability process.展开更多
Two calculation modes for the effect of external load on slope stability, i.e., mode I in which the external load is thought to act on slope surface, and mode II in which the external load is thought to act on slip su...Two calculation modes for the effect of external load on slope stability, i.e., mode I in which the external load is thought to act on slope surface, and mode II in which the external load is thought to act on slip surface along the force action line, were considered. Meanwhile, four basic distribution patterns of external load were used, of which complex external loads could be composed. In analysis process, several limit equilibrium methods, such as Swedish method, simplified Bishop method, simplified Janbu method, Spencer method, Morgenstern-Price(M-P) method, Sarma method, and unbalanced thrust method, were also adopted to contrast their differences in slope stability under the external load. According to parametric analysis, some conclusions can be obtained as follows:(1) The external load, with the large magnitude, small inclination angle, and acting position close to the slope toe,has more positive effect on slope stability;(2) The results calculated using modes I and II of external load are similar, indicating that the calculation mode of external load has little influence on slope stability;(3) If different patterns of external loads are equivalent to each other, their slope stability under these external loads are the same, and if not, the external load leads to the better slope stability,as action position of the resultant force for external load is closer to the lower sliding point of slip surface.展开更多
基金Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘In this paper,we develop a multi-scalar auxiliary variables(MSAV)scheme for the Cahn-Hilliard Magnetohydrodynamics system by introducing two scalar auxiliary variables(SAV).This scheme is linear,fully decoupled and unconditionally stable in energy.Subsequently,we provide a detailed implementation procedure for full decoupling.Thus,at each time step,only a series of linear differential equations with constant coefficients need to be solved.To validate the effectiveness of our approach,we conduct an error analysis for this first-order scheme.Finally,some numerical experiments are provided to verify the energy dissipation of the system and the convergence of the proposed approach.
基金Projects(52208369,52309138,52108320)supported by the National Natural Science Foundation of ChinaProjects(2023NSFSC0284,2025ZNSFSC0409)supported by the Sichuan Science and Technology Program,ChinaProject(U22468214)supported by the Joint Fund Project for Railway Basic Research by the National Natural Science Foundation of China and China State Railway Group Co.,Ltd.
文摘In the practical slope engineering,the stability of lower sliding mass(region A)with back tensile cracks of the jointed rock slope attracts more attentions,but the upper rock mass(region B)may also be unstable.Therefore,in this study,based on the stepped failure mode of bedding jointed rock slopes,considering the influence of the upper rock mass on the lower stepped sliding mass,the improved failure model for analyzing the interaction force(F_(AB))between two regions is constructed,and the safety factors(F_(S))of two regions and whole region are derived.In addition,this paper proposes a method to determine the existence of F_(AB) using their respective acceleration values(a_(A) and a_(B))when regions A and B are unstable.The influences of key parameters on two regions and the whole region are analyzed.The results show that the variation of the F_(AB) and F_(S) of two regions can be obtained accurately based on the improved failure model.The accuracy of the improved failure model is verified by comparative analysis.The research results can explain the interaction mechanism of two regions and the natural phenomenon of slope failure caused by the development of cracks.
文摘Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.
基金Project(51978085)supported by the National Natural Science Foundation of ChinaProject(201808430102)supported by the China Scholarship Council+1 种基金Project(JTG-201507)supported by the Highway Industry Standard Compilation Project of Ministry of Transportation,ChinaProject(kfj180102)supported by the Open Fund of Changsha University of Science&Technology,China。
文摘The traditional stability analysis method of geogrid reinforced slopes does not consider the effect of lateral swelling,so it is not applicable to reinforced expansive soil slopes.This paper reports a new stability analysis method for geogrid reinforced expansive soil slopes.The additional pullout force of the free zone due to the lateral swelling and the anti-pullout safety factor of each geogrid layer were obtained by ensuring the overall stability of the reinforced slope.The optimum design was carried out to treat an expansive soil cut slope in Hubei Province,China,by changing the spacing and length of geogrid reinforcement.Calculation results show that the additional pullout force caused by lateral swelling has a great influence on the anti-pullout stability of geogrids,and the local stability of the reinforced slope will be overestimated if the swelling effect of soil in the free zone is not considered.
基金Project(2010(A)06-b)supported by the Scientific Research Fund of Yunnan Provincial Transport Department of ChinaProject(51108293)supported by the National Natural Science Foundation of China+1 种基金Project(2013PY37)supported by the Cultivated Foundation of Taizhou University of ChinaProject(LY13E080008)supported by the Zhejiang Provincial Natural Science Foundation of China
文摘Rainfall infiltration depth and mode can severely influence slope stability.With the sustained rainfall,the influenced region of slope gradually expands.By using the Green-Ampt model to the soil slope,infiltration regulation was discussed under sustained and small intensity rainfall.And the infiltration rate of unsaturated soil was proposed according to the saturated infiltration theory.Because of the changing of initial moisture content in depth of slope,the saturated or unsaturated infiltration rate and depth could also be changeable with the sustained rainfall infiltration.Based on the principle of strength reduction,the calculation model of slope safety factor was established under different initial moisture contents and infiltration modes.Then,the slope stability was quantitatively analyzed through software FLAC3D.The calculation results of soil slope engineering show that there is a shorter period for slope stability under different initial moisture contents and unsaturated infiltration ways at the slope wetting front.The stability period of slope is 33.3%according to different initial moisture contents of wetting front less than that of the same initial moisture content of wetting front.And the slope is easier to fail under the unsaturated infiltration.The results agree well with the actual situation under sustained and small intensity rainfall.
基金Projects(41002095,41172251,41272317) supported by the National Natural Science Foundation of China
文摘A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.
基金Project(2017YFC1501100)supported by the National Key R&D Program of ChinaProjects(51809221,51679158)supported by the National Natural Science Foundation of China。
文摘To investigate the stability of rock mass in high geostress underground powerhouse caverns subjected to excavation,a microseismic(MS)monitoring system was established and the discrete element method(DEM)-based numerical simulation was carried out.The tempo-spatial damage characteristics of rock mass were analyzed.The evolution laws of MS source parameters during the formation of a rock collapse controlled by high geostress and geological structure were investigated.Additionally,a three-dimensional DEM model of the underground powerhouse caverns was built to reveal the deformation characteristics of rock mass.The results indicated that the MS events induced by excavation of high geostress underground powerhouse caverns occurred frequently.The large-stake crown of the main powerhouse was the main damage area.Prior to the rock collapse,the MS event count and accumulated energy release increased rapidly,while the apparent stress sharply increased and then decreased.The amount and proportion of shear and mixed MS events remarkably increased.The maximum displacement was generally located near the spandrel areas.The MS monitoring data and numerical simulation were in good agreement,which can provide significant references for damage evaluation and disaster forecasting in high geostress underground powerhouse caverns.
基金Project(11102218) supported by the National Natural Science Foundation of China
文摘The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a detailed calculation procedure and a definition of FOS for slope stability was developed based on the understanding of SRM. When constructing the new definition of FOS, efforts were made to make sure that it has concise physical meanings and fully reflects the shear strength of the slope. Two examples, slopes A and B with the slope angles of 63° and 34° respectively, were given to verify the method presented. It is found that, for these two slopes, the FOSs from original strength reduction method are respectively 1.5% and 38% higher than those from double reduction method. It is also found that the double reduction method predicts a deeper potential slide line and a larger slide mass. These results show that on one hand, the double reduction method is comparative to the traditional methods and is reasonable, and on the other hand, the original strength reduction method may overestimate the safety of a slope. The method presented is advised to be considered as an additional option in the practical slope stability evaluations although more useful experience is required.
基金Project(200550) supported by the Foundation for the Author of National Excellent Doctoral Dissertation of ChinaProject(09JJ1008) supported by Hunan Provincial Natural Science Foundation of ChinaProject(200631878557) supported by West Traffic of Science and Technology of China
文摘Based on strength reduction theory,the stability numbers of shallow tunnels were investigated within the framework of upper and lower bound theorems of limit analysis. Stability solutions taking into account of water seepage were presented and compared with those without considering seepage. The comparisons indicate that the maximum difference does not exceed 3.7%,which proves the present method credible. The results show that stability numbers of shallow tunnels considering seepage are much less than those without considering seepage,and that the difference of stability numbers between considering seepage and without considering seepage increase with increasing the depth ratio. The stability numbers decrease with increasing permeability coefficient and groundwater depth. Seepage has significant effects on the stability numbers of shallow tunnels.
基金Project(51808116) supported by the National Natural Science Foundation of ChinaProject(BK20180404) supported by the Natural Science Foundation of Jiangsu Province, China+1 种基金Project(KFJ170106) supported by the Changsha University of Science & Technology via Open Fund of National Engineering Laboratory of Highway Maintenance Technology,ChinaProject(242020R40133) supported by Fundamental Research Funds for the Central Universities, China。
文摘In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of the particle micro parameters on the macroscopic characteristics of cohesive soil and calibrate the micro parameters of DEM model on this basis. Embankment slope stability analysis was carried out by strength reduction and gravity increase method, it is shown that the safety factor obtained by strength reduction method is more conservative, and the arc-shaped feature of the sliding surface under the gravity increase method is more obvious. Throughout the progressive failure process, the failure trends, maximum displacements, and velocity changes obtained by the two methods were consistent. When slope was destroyed, the upper part was cracked, the middle part was sheared, and the lower part was destroyed by extrusion. The conclusions of this paper can be applied to the safety factor calculation of cohesive soil slopes and the analysis of the instability process.
基金Project(2015M580702)supported by the China Postdoctoral Science FoundationProject(51608541)supported by the National Natural Science Foundation of ChinaProject(2014122006)supported by the Guizhou Provincial Department of Transportation Foundation,China
文摘Two calculation modes for the effect of external load on slope stability, i.e., mode I in which the external load is thought to act on slope surface, and mode II in which the external load is thought to act on slip surface along the force action line, were considered. Meanwhile, four basic distribution patterns of external load were used, of which complex external loads could be composed. In analysis process, several limit equilibrium methods, such as Swedish method, simplified Bishop method, simplified Janbu method, Spencer method, Morgenstern-Price(M-P) method, Sarma method, and unbalanced thrust method, were also adopted to contrast their differences in slope stability under the external load. According to parametric analysis, some conclusions can be obtained as follows:(1) The external load, with the large magnitude, small inclination angle, and acting position close to the slope toe,has more positive effect on slope stability;(2) The results calculated using modes I and II of external load are similar, indicating that the calculation mode of external load has little influence on slope stability;(3) If different patterns of external loads are equivalent to each other, their slope stability under these external loads are the same, and if not, the external load leads to the better slope stability,as action position of the resultant force for external load is closer to the lower sliding point of slip surface.
