The calculation of the factor of safety(FOS)is an important means of slope evaluation.This paper proposed an improved double strength reductionmethod(DRM)to analyze the safety of layered slopes.The physical properties...The calculation of the factor of safety(FOS)is an important means of slope evaluation.This paper proposed an improved double strength reductionmethod(DRM)to analyze the safety of layered slopes.The physical properties of different soil layers of the slopes are different,so the single coefficient strength reduction method(SRM)is not enough to reflect the actual critical state of the slopes.Considering that the water content of the soil in the natural state is the main factor for the strength of the soil,the attenuation law of shear strength of clayey soil changing with water content is fitted.This paper also establishes the functional relationship between different reduction coefficients.Then,a USDFLD subroutine is programmed using the secondary development function of finite element software.Controlling the relationship between field variables and calculation time realizes double strength reduction applicable to the layered slope.Finally,by comparing the calculation results of different examples,it is proved that the stress and displacement distribution of the critical slope state obtained by the improved method is more realistic,and the calculated safety factor is more reliable.The newly proposedmethod considers the difference of intensity attenuation between different soil layers under natural conditions and avoids the disadvantage of the strength reduction method with uniform parameters,which provides a new idea and method for stability analysis of layered and complex slopes.展开更多
This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuit...This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuitive particle-relocation algorithm ensuring material points always remain within the computational mesh.The explicit generalized-a time scheme was implemented in MPM to enable the damping of spurious high frequency oscillations.Firstly,the MPM was verified against finite element method(FEM).Secondly,ability of the MPM in capturing the analytical transfer function was investigated.Thirdly,a symmetric embankment was adopted to investigate the effects of ground motion arias intensity(I_(a)),geometry dimensions,and constitutive models.The results show that the larger the model size,the higher the crest runout and settlement for the same ground motion.When using a Mohr-Coulomb model,the crest runout increases with increasing I_(a).However,if the strain-softening law is activated,the results are less influenced by the ground motion.Finally,the MPM results were compared with the Newmark sliding block solution.The simplified analysis herein highlights the capabilities of MPM to capture the full deformation process for earthquake engineering applications,the importance of geometry characterization,and the selection of appropriate constitutive models when simulating coseismic site response and subsequent large deformations.展开更多
To study the safety and stability of large slopes, taking the right side slope of the Yuxi’an tunnel of the Yuchu Expressway Bridge in Yunnan Province as an example, limit equilibrium and finite element analysis were...To study the safety and stability of large slopes, taking the right side slope of the Yuxi’an tunnel of the Yuchu Expressway Bridge in Yunnan Province as an example, limit equilibrium and finite element analysis were applied to engineering examples to calculate the stability coefficient of the slope before and after excavation in the natural state. After comparative analysis, it was concluded that the former had a clear mechanical model and concept, which could quickly provide stability results;the latter could accurately determine the sliding surface of the slope and simulate the stress state changes of the rock and soil mass. The stability coefficients calculated by the two methods were within the stable range, but their values were different. On this basis, combined with the calculation principles, advantages and disadvantages of the two methods, a comprehensive analysis method of slope stability based on the limit equilibrium and finite element methods was proposed, and the rationality of the stability coefficient calculated by this method was judged for a slope case.展开更多
Infiltration–runoff–slope instability mechanism of macropore slope under heavy rainfall is unclear.This paper studied its instability mechanism with an improved Green–Ampt(GA)model considering the dual-porosity(i.e...Infiltration–runoff–slope instability mechanism of macropore slope under heavy rainfall is unclear.This paper studied its instability mechanism with an improved Green–Ampt(GA)model considering the dual-porosity(i.e.,matrix and macropore)and ponding condition,and proposed the infiltration equations,infiltration–runoff coupled model,and safety factor calculation method.Results show that the infiltration processes of macropore slope can be divided into three stages,and the proposed model is rational by a comparative analysis.The wetting front depth of the traditional unsaturated slope is 17.2%larger than that of the macropore slope in the early rainfall stage and 27%smaller than that of the macropore slope in the late rainfall stage.Then,macropores benefit the slope stability in the early rainfall but not in the latter.Macropore flow does not occur initially but becomes pronounced with increasing rainfall duration.The equal depth of the wetting front in the two domains is regarded as the onset criteria of macropore flow.Parameter analysis shows that macropore flow is delayed by increasing proportion of macropore domain(ω_(f)),whereas promoted by increasing ratio of saturated permeability coefficients between the two domains(μ).The increasing trend of ponding depth is sharp at first and then grows slowly.Finally,when rainfall duration is less than 3 h,ωf andμhave no significant effect on the safety factor,whereas it decreases with increasingωf and increases with increasingμunder longer duration(≥3 h).With the increase ofω_(f),the slope maximum instability time advances by 10.5 h,and with the increase ofμ,the slope maximum instability time delays by 3.1 h.展开更多
Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model sl...Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.展开更多
The instability of slope blocks occurred frequently along traffic corridor in Southeastern Tibet(TCST),which was primarily controlled by the rock mass structures.A rapid method evaluating the control effects of rock m...The instability of slope blocks occurred frequently along traffic corridor in Southeastern Tibet(TCST),which was primarily controlled by the rock mass structures.A rapid method evaluating the control effects of rock mass structures was proposed through field statistics of the slopes and rock mass structures along TCST,which combined the stereographic projection method,modified M-JCS model,and limit equilibrium theory.The instabilities of slope blocks along TCST were then evaluated rapidly,and the different control factors of instability were analyzed.Results showed that the probabilities of toppling(5.31%),planar(16.15%),and wedge(35.37%)failure of slope blocks along TCST increased sequentially.These instability modes were respectively controlled by the anti-dip joint,the joint parallel to slope surface with a dip angle smaller than the slope angle(singlejoint),and two groups of joints inclined out of the slope(double-joints).Regarding the control effects on slope block instability,the stabilization ability of doublejoints(72.7%),anti-dip joint(67.4%),and single-joint(57.6%)decreased sequentially,resulting in different probabilities of slope block instability.Additionally,nearby regional faults significantly influenced the joints,leading to spatial heterogeneity and segmental clustering in the stabilization ability provided by joints to the slope blocks.Consequently,the stability of slope blocks gradually weakened as they approached the fault zones.This paper can provide guidance and assistance for investigating the development characteristics of rock mass structures and the stability of slope blocks.展开更多
Repetitive mining beneath bedding slopes is identified as a critical factor in geomorphic disturbances, especially landslides and surface subsidence. Prior research has largely concentrated on surface deformation in p...Repetitive mining beneath bedding slopes is identified as a critical factor in geomorphic disturbances, especially landslides and surface subsidence. Prior research has largely concentrated on surface deformation in plains due to multi-seam coal mining and the instability of natural bedding slopes, yet the cumulative impact of different mining sequences on bedding slopes has been less explored. This study combines drone surveys and geological data to construct a comprehensive three-dimensional model of bedding slopes. Utilizing FLAC3D and PFC2D models, derived from laboratory experiments, it simulates stress, deformation, and failure dynamics of slopes under various mining sequences. Incorporating fractal dimension analysis, the research evaluates the stability of slopes in relation to different mining sequences. The findings reveal that mining in an upslope direction minimizes disruption to overlying strata. Initiating extraction from lower segments increases tensile-shear stress in coal pillar overburdens, resulting in greater creep deformation towards the downslope than when starting from upper segments, potentially leading to localized landslides and widespread creep deformation in mined-out areas. The downslope upward mining sequence exhibits the least fractal dimensions, indicating minimal disturbance to both strata and surface. While all five mining scenarios maintain good slope stability under normal conditions, recalibrated stability assessments based on fractal dimensions suggest that downslope upward mining offers the highest stability under rainfall, contrasting with the lower stability and potential instability risks of upslope downward mining. These insights are pivotal for mining operations and geological hazard mitigation in multi-seam coal exploitation on bedding slopes.展开更多
The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experi...The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.展开更多
This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering a...This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.展开更多
The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have diff...The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a...Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.展开更多
The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method (DFEM). Combining with the DFEM, the vector sum analysis method (VSAM) is...The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method (DFEM). Combining with the DFEM, the vector sum analysis method (VSAM) is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.展开更多
Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlight...Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.展开更多
Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to ...Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.展开更多
The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CG...The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.展开更多
Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafros...Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment design in frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of material properties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer (i.e., mass equation) and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety (FOS), which is definitely useful for embankment design in frozen ground regions.展开更多
In this paper, the ring-type ingot of hypereutectic high Cr cast iron was obtained by slope cooling bodycentrifugal casting method (SC-CCM), and its microstructure and impact toughness were investigated, respectivel...In this paper, the ring-type ingot of hypereutectic high Cr cast iron was obtained by slope cooling bodycentrifugal casting method (SC-CCM), and its microstructure and impact toughness were investigated, respectively. The results indicated that, first, the primary carbides in the microstructure are prominently finer than those in the hypereutectic high Cr cast iron prepared by conventional casting method. Second, in the ring-type ingot, the primary carbides near radial outer field are finer than those near radial inner field; furthermore, there is dividing field in the microstructure. Finally, the impact toughness values of the specimens impacted on the radial outer face and on the radial inner face are improved respectively about 36% and 138% more than that of the hypereutectic high Cr one prepared by conventional casting method.展开更多
The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calcula...The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calculated from the phenomenological models to deviate from the actual SFT.Currently,very limited study has been conducted on how to evaluate the effect of such uncertainties on SFT prediction.In this paper,a comprehensive slope failure database was compiled.A Bayesian machine learning(BML)-based method was developed to learn the model and observational uncertainties involved in SFT prediction,through which the probabilistic distribution of the SFT can be obtained.This method was illustrated in detail with an example.Verification studies show that the BML-based method is superior to the traditional inverse velocity method(INVM)and the maximum likelihood method for predicting SFT.The proposed method in this study provides an effective tool for SFT prediction.展开更多
This study first reviews the numerical manifold method(NMM)which possesses some advantages over the traditional limit equilibrium methods(LEMs)in calculating the factors of safety(Fs)of the slopes.Then,with regard to ...This study first reviews the numerical manifold method(NMM)which possesses some advantages over the traditional limit equilibrium methods(LEMs)in calculating the factors of safety(Fs)of the slopes.Then,with regard to a trial slip surface(TSS),associated stress fields reproduced by NMM as well as the enhanced limit equilibrium method are combined to compute Fs.In order to search for the potential critical slip surface(CSS),the MAX-MIN ant colony optimization algorithm(MMACOA),one of the best performing algorithms for some optimization problems,is adopted.Procedures to obtain Fs in conjunction with the potential CSS are described.Finally,the proposed numerical model and traditional methods are compared with stability analysis of three typical slopes.The numerical results show that Fs and CSSs of the slopes can be accurately calculated with the proposed model.展开更多
基金This research was funded by the National Natural Science Foundation of China(51709194),Qinglan Project of Jiangsu University,the Priority Academic Program Development of Jiangsu Higher Education Institutions,and Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering.
文摘The calculation of the factor of safety(FOS)is an important means of slope evaluation.This paper proposed an improved double strength reductionmethod(DRM)to analyze the safety of layered slopes.The physical properties of different soil layers of the slopes are different,so the single coefficient strength reduction method(SRM)is not enough to reflect the actual critical state of the slopes.Considering that the water content of the soil in the natural state is the main factor for the strength of the soil,the attenuation law of shear strength of clayey soil changing with water content is fitted.This paper also establishes the functional relationship between different reduction coefficients.Then,a USDFLD subroutine is programmed using the secondary development function of finite element software.Controlling the relationship between field variables and calculation time realizes double strength reduction applicable to the layered slope.Finally,by comparing the calculation results of different examples,it is proved that the stress and displacement distribution of the critical slope state obtained by the improved method is more realistic,and the calculated safety factor is more reliable.The newly proposedmethod considers the difference of intensity attenuation between different soil layers under natural conditions and avoids the disadvantage of the strength reduction method with uniform parameters,which provides a new idea and method for stability analysis of layered and complex slopes.
基金funded by National Science Foundation(NSF)(Grant No.CMMI-2211002).
文摘This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuitive particle-relocation algorithm ensuring material points always remain within the computational mesh.The explicit generalized-a time scheme was implemented in MPM to enable the damping of spurious high frequency oscillations.Firstly,the MPM was verified against finite element method(FEM).Secondly,ability of the MPM in capturing the analytical transfer function was investigated.Thirdly,a symmetric embankment was adopted to investigate the effects of ground motion arias intensity(I_(a)),geometry dimensions,and constitutive models.The results show that the larger the model size,the higher the crest runout and settlement for the same ground motion.When using a Mohr-Coulomb model,the crest runout increases with increasing I_(a).However,if the strain-softening law is activated,the results are less influenced by the ground motion.Finally,the MPM results were compared with the Newmark sliding block solution.The simplified analysis herein highlights the capabilities of MPM to capture the full deformation process for earthquake engineering applications,the importance of geometry characterization,and the selection of appropriate constitutive models when simulating coseismic site response and subsequent large deformations.
文摘To study the safety and stability of large slopes, taking the right side slope of the Yuxi’an tunnel of the Yuchu Expressway Bridge in Yunnan Province as an example, limit equilibrium and finite element analysis were applied to engineering examples to calculate the stability coefficient of the slope before and after excavation in the natural state. After comparative analysis, it was concluded that the former had a clear mechanical model and concept, which could quickly provide stability results;the latter could accurately determine the sliding surface of the slope and simulate the stress state changes of the rock and soil mass. The stability coefficients calculated by the two methods were within the stable range, but their values were different. On this basis, combined with the calculation principles, advantages and disadvantages of the two methods, a comprehensive analysis method of slope stability based on the limit equilibrium and finite element methods was proposed, and the rationality of the stability coefficient calculated by this method was judged for a slope case.
基金funded by the Natural Science Foundation of Fujian Province(Grant No.2023J011133)。
文摘Infiltration–runoff–slope instability mechanism of macropore slope under heavy rainfall is unclear.This paper studied its instability mechanism with an improved Green–Ampt(GA)model considering the dual-porosity(i.e.,matrix and macropore)and ponding condition,and proposed the infiltration equations,infiltration–runoff coupled model,and safety factor calculation method.Results show that the infiltration processes of macropore slope can be divided into three stages,and the proposed model is rational by a comparative analysis.The wetting front depth of the traditional unsaturated slope is 17.2%larger than that of the macropore slope in the early rainfall stage and 27%smaller than that of the macropore slope in the late rainfall stage.Then,macropores benefit the slope stability in the early rainfall but not in the latter.Macropore flow does not occur initially but becomes pronounced with increasing rainfall duration.The equal depth of the wetting front in the two domains is regarded as the onset criteria of macropore flow.Parameter analysis shows that macropore flow is delayed by increasing proportion of macropore domain(ω_(f)),whereas promoted by increasing ratio of saturated permeability coefficients between the two domains(μ).The increasing trend of ponding depth is sharp at first and then grows slowly.Finally,when rainfall duration is less than 3 h,ωf andμhave no significant effect on the safety factor,whereas it decreases with increasingωf and increases with increasingμunder longer duration(≥3 h).With the increase ofω_(f),the slope maximum instability time advances by 10.5 h,and with the increase ofμ,the slope maximum instability time delays by 3.1 h.
基金The authors wish to thank National Key R&D Program of China(Grant No.2022YFC308100)the National Nature Science Foundation of China(Grant Nos.42107172 and 42072303)for financial support.
文摘Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.41941019,42177142)the Second Tibetan Plateau Scientific Expedition and Research(STEP)program(Grant NO.2019QZKK0904)the Fundamental Research Funds for the Central Universities,CHD(Grant No.300102212213).
文摘The instability of slope blocks occurred frequently along traffic corridor in Southeastern Tibet(TCST),which was primarily controlled by the rock mass structures.A rapid method evaluating the control effects of rock mass structures was proposed through field statistics of the slopes and rock mass structures along TCST,which combined the stereographic projection method,modified M-JCS model,and limit equilibrium theory.The instabilities of slope blocks along TCST were then evaluated rapidly,and the different control factors of instability were analyzed.Results showed that the probabilities of toppling(5.31%),planar(16.15%),and wedge(35.37%)failure of slope blocks along TCST increased sequentially.These instability modes were respectively controlled by the anti-dip joint,the joint parallel to slope surface with a dip angle smaller than the slope angle(singlejoint),and two groups of joints inclined out of the slope(double-joints).Regarding the control effects on slope block instability,the stabilization ability of doublejoints(72.7%),anti-dip joint(67.4%),and single-joint(57.6%)decreased sequentially,resulting in different probabilities of slope block instability.Additionally,nearby regional faults significantly influenced the joints,leading to spatial heterogeneity and segmental clustering in the stabilization ability provided by joints to the slope blocks.Consequently,the stability of slope blocks gradually weakened as they approached the fault zones.This paper can provide guidance and assistance for investigating the development characteristics of rock mass structures and the stability of slope blocks.
基金funded by the Sichuan Science and Technology Program (grant number 2022NSFSC1176)the open Fund for National Key Laboratory of Geological Disaster Prevention and Environmental Protection (grant number SKLGP2022K027)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project (SKLGP2022Z001)。
文摘Repetitive mining beneath bedding slopes is identified as a critical factor in geomorphic disturbances, especially landslides and surface subsidence. Prior research has largely concentrated on surface deformation in plains due to multi-seam coal mining and the instability of natural bedding slopes, yet the cumulative impact of different mining sequences on bedding slopes has been less explored. This study combines drone surveys and geological data to construct a comprehensive three-dimensional model of bedding slopes. Utilizing FLAC3D and PFC2D models, derived from laboratory experiments, it simulates stress, deformation, and failure dynamics of slopes under various mining sequences. Incorporating fractal dimension analysis, the research evaluates the stability of slopes in relation to different mining sequences. The findings reveal that mining in an upslope direction minimizes disruption to overlying strata. Initiating extraction from lower segments increases tensile-shear stress in coal pillar overburdens, resulting in greater creep deformation towards the downslope than when starting from upper segments, potentially leading to localized landslides and widespread creep deformation in mined-out areas. The downslope upward mining sequence exhibits the least fractal dimensions, indicating minimal disturbance to both strata and surface. While all five mining scenarios maintain good slope stability under normal conditions, recalibrated stability assessments based on fractal dimensions suggest that downslope upward mining offers the highest stability under rainfall, contrasting with the lower stability and potential instability risks of upslope downward mining. These insights are pivotal for mining operations and geological hazard mitigation in multi-seam coal exploitation on bedding slopes.
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.
基金funded by the National Key Research and Development Program of China(No.2021YFB2600704)the National Natural Science Foundation of China(No.52171272)the Significant Science and Technology Project of the Ministry of Water Resources of China(No.SKS-2022112).
文摘This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.
文摘The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.
基金Supported by the Program of Yunnan Provincial Institute of Communications Planning,Design and Research (2011(D)11-b)
文摘The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method (DFEM). Combining with the DFEM, the vector sum analysis method (VSAM) is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.
文摘Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.
基金NRDMS Division,Department of Science and Technology,Government of India for providing financial assistance for field investigations.
文摘Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.
基金National Natural Science Foundation of China under grant No.59839330China Postdoctoral Science Foundation
文摘The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.
基金supported by the National 973 Project of China (No. 2012CB026104)the National Natural Science Foundation of China (No. 51378057)
文摘Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment design in frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of material properties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer (i.e., mass equation) and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety (FOS), which is definitely useful for embankment design in frozen ground regions.
基金This work was supported by the National Natural Science Foundation of China under grant No.50571079.
文摘In this paper, the ring-type ingot of hypereutectic high Cr cast iron was obtained by slope cooling bodycentrifugal casting method (SC-CCM), and its microstructure and impact toughness were investigated, respectively. The results indicated that, first, the primary carbides in the microstructure are prominently finer than those in the hypereutectic high Cr cast iron prepared by conventional casting method. Second, in the ring-type ingot, the primary carbides near radial outer field are finer than those near radial inner field; furthermore, there is dividing field in the microstructure. Finally, the impact toughness values of the specimens impacted on the radial outer face and on the radial inner face are improved respectively about 36% and 138% more than that of the hypereutectic high Cr one prepared by conventional casting method.
基金substantially supported by the Shuguang Program from Shanghai Education Development FoundationShanghai Municipal Education Commission, China (Grant No. 19SG19)+1 种基金National Natural Science Foundation of China (Grant No. 42072302)Fundamental Research Funds for the Central Universities, China
文摘The data-driven phenomenological models based on deformation measurements have been widely utilized to predict the slope failure time(SFT).The observational and model uncertainties could lead the predicted SFT calculated from the phenomenological models to deviate from the actual SFT.Currently,very limited study has been conducted on how to evaluate the effect of such uncertainties on SFT prediction.In this paper,a comprehensive slope failure database was compiled.A Bayesian machine learning(BML)-based method was developed to learn the model and observational uncertainties involved in SFT prediction,through which the probabilistic distribution of the SFT can be obtained.This method was illustrated in detail with an example.Verification studies show that the BML-based method is superior to the traditional inverse velocity method(INVM)and the maximum likelihood method for predicting SFT.The proposed method in this study provides an effective tool for SFT prediction.
基金This study is supported by the Youth Innovation Promotion Association of Chinese Academy of Sciences(Grant No.2020327)the National Natural Science Foundation of China(Grant No.51609240).
文摘This study first reviews the numerical manifold method(NMM)which possesses some advantages over the traditional limit equilibrium methods(LEMs)in calculating the factors of safety(Fs)of the slopes.Then,with regard to a trial slip surface(TSS),associated stress fields reproduced by NMM as well as the enhanced limit equilibrium method are combined to compute Fs.In order to search for the potential critical slip surface(CSS),the MAX-MIN ant colony optimization algorithm(MMACOA),one of the best performing algorithms for some optimization problems,is adopted.Procedures to obtain Fs in conjunction with the potential CSS are described.Finally,the proposed numerical model and traditional methods are compared with stability analysis of three typical slopes.The numerical results show that Fs and CSSs of the slopes can be accurately calculated with the proposed model.