Numerical investigations of rock bridge effect on open pit slope stability

In this study,the effect of rock bridges on rock slope stability was investigated by incorporating nonpersistent joint networks in numerical models,and the critical profiles of an open pit mine were analysed.Parallel deterministic networks of infinite and finite lengths,ubiquitous joint network model and Veneziano joint network model were used in order to simulate the rock fractures.Materials were modelled based on the generalised Hoek-Brown and equivalent Mohr-Coulomb failure criteria.The parallel deterministic infinite and the ubiquitous joint network models produced lower safety factors.The introduction of rock bridges along discontinuity planes in the parallel deterministic network and Veneziano joint network models significantly contributed to the stability and strain distribution,which should be considered in stability analysis of rock mass in open pit by rock slope practitioners.The results show the significance of joints in hard rock behaviour and the joints should be included in order to attain practical and realistic simulations.

Numerical analysis of pile–slope stability and the soil arching around two adjacent piles

Seismic pile–slope stability analysis and the formation mechanism of soil arching have not been well studied. This study used a threedimensional(3D) finite difference to determine soil and pile parameter changes in the static and seismic stability of the pile–slope caused by the interaction between stabilizing piles. Pile–slope stability analysis was performed to determine the optimal design of piles along a slope and the corresponding failure mode involving the formation of soil arching around two adjacent piles. The Factor of Safety(FS) of the slope was evaluated using the shear strength reduction method for static and seismic analyses. The results of the analysis show that suitable pile spacing(S) and a suitable pile diameter(D) in the middle of a slope result in the maximum FS for the pile–slope system and the formation of soil arching around two adjacent piles. FS of the pile–slope increased negligibly in the seismic analysis of piles located at the slope crest and toe. An optimized pile diameter and installation location afforded the maximum FS for the slope that corresponded to a specified slope failure mode for different pile locations. A pile spacing S ≤ 2.5D for piles installed in the middle of the slope is suggested for increasing the static and seismic pile–slope stability.

Critical Assessment of Slope Stability: A Case Study on the Toffo-Lalo Road Project

This article systematically delves into a comprehensive analysis of the latest and most advanced techniques for the assessment of slope stability. It particularly focuses on strategies aimed at enhancing slope stability in road construction. In addition to this analysis, the article presents an illustrative case study centered on the Toffo-Lalo Road Project. The core objective of this paper is to scrutinize the stability of large embankments in road construction, with a specific emphasis on the development and asphalt overlay of the Toffo-Lalo road. This scrutiny is conducted through the utilization of stability calculation software, GEOSTUDIO2018, specifically its SLOPE/W module. Within this framework, a detailed model of the cutbank located at KP1+750-2+250 was meticulously developed. This model takes into account the physical-mechanical characteristics of the soil at the site, as well as the topographic layout. Its attributes include a cohesion value of 11.3 Kpa, a density of 16.57 KN/m3, and a friction angle of 27˚. The modeling results, employing the Morgenstern-Price method—an approach renowned for its adherence to equilibrium conditions and provision of precise results—conclude that the safety coefficient (Fs = 1.429) prior to any reinforcement signifies a critical state of slope stability. To address this, the article explores the implementation of reinforcement techniques, particularly focusing on rigid inclusions like nailing and piles. The modeling exercises reveal a noteworthy enhancement in the safety coefficient (Fs) post-reinforcement. Furthermore, the article undertakes a parametric study to optimize the reinforcement strategies. This analysis highlights that anchoring at 0˚ downward relative to the horizontal plane and employing a pile angle of 90˚ represent the most favorable approaches. These measures yield safety coefficients of 3.60 and 2.34, respectively, indicating substantially improved slope stability.

3D Identification and Stability Analysis of Key Surface Blocks of Rock Slope

Complicated geological structures make it difficult to analyze the stability of rock slopes, such as faults, weak intercalated layers or joint fissures. Based on 3D geological modeling and surface block identifying methods, an integrated methodology framework was proposed and realized to analyze the stability of surface blocks in rock slopes. The surface blocks cut by geological structures, fissures or free faces could be identified subjected to the four principles of closure, completeness, uniqueness and validity. The factor of safety(FOS)of single key block was calculated by the limit equilibrium method. If there were two or more connected blocks, they were defined as a block-group. The FOS of a block-group was computed by the Sarma method. The proposed approach was applied to an actual rock slope of a hydropower project, and some possible instable blocks were demonstrated and analyzed visually. The obtained results on the key blocks or block-groups provide essential information for determining potential instable region of rock slopes and designing effective support scheme in advance.

Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Rock-soil slope stability analysis by two-phase random media and finite elements

To investigate the strong random nature of the geometric interfaces between soil and rock, a rock-soil slope is considered as a two-phase random medium. A nonlinear translation of a Gaussian field is utilized to simulate the two-phase random media, such that the soil(or rock) volume fraction and the inclination of the soil layer can be examined. The finite element method with random media incorporated as the material properties is used to determine the factor of safety of the rock-soil slope. Monte-Carlo simulations are used to estimate the statistical characteristics of the factor of safety. The failure mode of the rock-soil slope is examined by observing the maximum principal plastic strain at incipient slope failure. It is found that the critical surface of a rock-soil slope is fairly irregular, and it significantly differs from that of a pure soil slope. The factor of safety is sensitive to the soil volume faction, but it is predictable. The average factor of safety could be well predicted by the weighted harmonic average between the strength of soil and rock; the prediction model is practical and simple. Parametric studies on the inclination of the soil layer demonstrate that the most instable scenario occurs when the slope angle is consistent with the inclination of the soil layer.

Stability assessment of landslide-prone road cut rock slopes in Himalayan terrain:A finite element method based approach

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.

Development of a preliminary slope stability calculation method based on internal horizontal displacements

This paper proposes a horizontal displacement-based approach to determine the potential slip surface of the slope. Firstly, a group of in-situ inclinometers with an appropriate spacing in the horizontal direction is located in the model slope.The equation of horizontal displacement with time for each in-situ inclinometer is fitted during the whole simulation process. Furthermore, the intersection of each inclinometer with potential slip surface is determined by using an optimization model. The slip surface can be obtained by using least square fitting method. Finally, the feasibility and accuracy of the method are validated by a series of numerical simulations. It is noted that the optimization model taking the maximum value of displacement increment gradient as an objective has higher accuracy when compared with other optimization models. This method employed in this study provides a preliminary approach to determine the real-time slope stability based on displacement, which can also be measured by using conventional instruments on site.

Full 3-D numerical modeling of borehole electric image logging and the evaluation model of fracture

A full 3-D finite element method numerical modeling program is written based on the principle and technical specification of borehole electric image well logging tool. The response of well logging is computed in the formation media model with a single fracture. The effect of changing fracture aperture and resistivity ratio to the logging response is discussed. The identification ability for two parallel fractures is also present. A quantitative evaluation formula of fracture aperture from borehole electric image logging data is set up. A case study of the model well is done to verify the accuracy of the for-mula. The result indicates that the formula is more accurate than the foreign one.

Exploring mechanism of hidden,steep obliquely inclined bedding landslides using a 3DEC model:A case study of the Shanyang landslide in Shaanxi Province,China

Catastrophic geological disasters frequently occur on slopes with obliquely inclined bedding structures(also referred to as obliquely inclined bedding slopes),where the apparent dip sliding is not readily visible.This phenomenon has become a focal point in landslide research.Yet,there is a lack of studies on the failure modes and mechanisms of hidden,steep obliquely inclined bedding slopes.This study investigated the Shanyang landslide in Shaanxi Province,China.Using field investigations,laboratory tests of geotechnical parameters,and the 3DEC software,this study developed a numerical model of the landslide to analyze the failure process of such slopes.The findings indicate that the Shanyang landslide primarily crept along a weak interlayer under the action of gravity.The landslide,initially following a dip angle with the support of a stable inclined rock mass,shifted direction under the influence of argillization in the weak interlayer,moving towards the apparent dip angle.The slide resistance effect of the karstic dissolution zone was increasingly significant during this process,with lateral friction being the primary resistance force.A reduction in the lateral friction due to karstic dissolution made the apparent dip sliding characteristics of the Shanyang landslide more pronounced.Notably,deformations such as bending and uplift at the slope’s foot suggest that the main slide resistance shifts from lateral friction within the karstic dissolution zone to the slope foot’s resistance force,leading to the eventual buckling failure of the landslide.This study unveils a novel failure mode of apparent dip creep-buckling in the Shanyang landslide,highlighting the critical role of lateral friction from the karstic dissolution zone in its failure mechanism.These insights offer a valuable reference for mitigating risks and preventing disasters related to obliquely inclined bedding landslides.

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG)terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.

Evaluation of potential failure of rock slope at the left abutment of Jinsha River Bridge by model test and numerical method

Jinsha River Bridge is located on Tiger Leaping Gorge town,China.The left bank slope composes of moderately thick layer of slate overlain by schistose basalt,and where rocks are controlled by two sets of joint planes.To evaluate the stability of the rock slope under bridge foundation,model test and calculation model based on the geological parameters and the slope stability was simulated and analyzed using Universal Distinct Element Code(UDEC)and Finite Element Mehod(FEM).According to model test,failure mainly initiated at the toe with shear movement along the joint planes,eventually resulting in the sliding along the slope surface and formation of tension crack at the crest of the model.This result coincide with the UDEC model,which shows that slope surface will produce loosening damage and slipping expected along the joint planes under different loading conditions.Moreover,the result of FEM analysis indicates that the rock mass under the main pier has potential shear failure region.So,the bridge foundation should be strengthened to prevent the slope failure under external forces.

Effect of surface retaining elements on rock stability:laboratory investigation with sand powder 3D printing

This study aims to investigate the benefcial efects of surface retaining elements (SREs) on the mechanical behaviors of bolted rock and roadway stability. 3D printing (3DP) technology is utilized to create rock analogue prismatic specimens for conducting this investigation. Uniaxial compression tests with acoustic emission (AE) and digital image correlation techniques have been conducted on 3DP specimens bolted with diferent SREs. The results demonstrate that the strength and modulus of elasticity of the bolted specimens show a positive correlation with the area of the SRE;the AE characteristics of the bolted specimens are higher than those of the unbolted specimen, but they decrease with an increase in SRE area, thus further improving the integrity of the bolted specimens. The reinforcement efect of SREs on the surrounding rock of roadways is further analyzed using numerical modelling and feld test. The results provide a better understanding of the role of SREs in rock bolting and the optimization of rock bolting design. Furthermore, they verify the feasibility of 3DP for rock analogues in rock mechanics tests.

Numerical modeling calculation for the spatial distribution characteristics of horizontal field transfer functions

Applying 3-dimension finite difference method, the distribution characteristics of horizontal field transfer functions for rectangular conductor have been computed, and the law of distribution for Re-part and Im-part has been given. The influences of source field period, the conductivity, the buried depth and the length of the conductor on the transfer functions were studied. The extrema of transfer functions appear at the center, the four corners and around the edges of conductor, and move with the edges. This feature demonstrates that around the edges are best places for transfer functions' observation.

Simulation of subcritical flow pattern in 180° uniform and convergent open-channel bends using SSIIM 3-D model

In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel＇s centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.

3-D Modelling of the Confederation Bridge Using Data of Full Scale Tests

Long-span bridges are special structures that require advanced analysis techniques to examine their performance. This paper presents a procedure developed to model the Confederation Bridge using 3-D beam elements. The model was validated using the data collected before the opening of the bridge to the public. The bridge was instrumented to conduct fullscale static and dynamic tests. The static tests were to measure the deflection of the bridge pier while the dynamic tests to measure the free vibrations of the pier due to a sudden release of the static load. Confederation Bridge is one of the longest reinforced concrete bridges in the world. It connects the province of Prince Edward Island and the province of New Brunswick in Canada. Due to its strategic location and vital role as a transportation link between these two provinces, it was designed using higher safety factors than those for typical highway bridges. After validating the present numerical model, a procedure was developed to evaluate the performance of similar bridges subjected to traffic and seismic loads. It is of interest to note that the foundation stiffness and the modulus of elasticity of the concrete have significant effects on the structural responses of the Confederation Bridge.

A Simplified Numerical Approach for the Prediction of Rainfall-Induced Retrogressive Landslides

Retrogressive landslides are common geological phenomena in mountainous areas and on onshore and offshore slopes. The impact of retrogressive landslides is different from that of other landslide types due to the phenomenon of retrogression. The hazards caused by retrogressive landslides may be increased because retrogressive landslides usually affect housing, facilities, and infrastructure located far from the original slopes. Additionally, substantial geomorphic evidence shows that the abundant supply of loose sediment in the source area of a debris flow is usually provided by retrogressive landslides that are triggered by the undercutting of water. Moreover, according to historic case studies, some large landslides are the evolution result of retrogressive landslides. Hence the ability to understand and predict the evolution of retrogressive landslides is crucial for the purpose of hazard mitigation. This paper discusses the phenomenon of a retrogressive landslide by using a model experiment and suggests a reasonably simplified numerical approach for the prediction of rainfall-induced retrogressive landslides. The simplified numerical approach, which combines the finite element method for seepage analysis, the shear strength reduction finite element method, and the analysis criterion for the retrogression and accumulation effect, is presented and used to predict the characteristics of a retrogressive landslide. The results show that this numerical approach is capable of reasonably predicting the characteristics of retrogressive landslides under rainfall infiltration, particularly the magnitude of each landslide, the position of the slip surface, and the development processes of the retrogressive landslide. Therefore, this approach is expected to be a practical method for the mitigation of damage caused by rainfall-induced retrogressive landslides.

Approach of moving boundary and its application in 3D tidal current simulation of Changjiang (Yangtze) River Estuary based on ECOMSED model

Considering three-dimensional model ECOMSED can not simulate wetting-drying of shoal with its fixed boundary, an approach to represent moving boundary in the model is introduced here. This approach smoothly joints the internal and external mode by making use of wetting and drying technique and is verified by a numerical test which presents a good agreement with the previous test results obtained by other researchers. A three dimensional numerical model is established to simulate the hydrodynamics in spring tide in the Changjiang （Yangtze） River estuary by this modified ECOMSED model which is also validated through the observed field data, the simulation presents a good periodic tidal change. It also successfully simulates the tidal current of computational areas and reproduces the tidal flat intermittent appearance.

Numerical Si mulation of Transportation of SPMfrom the Yellow River to the Bohai Sea

A 3D diagnostic model including a suspended sediment transport module, ECOMSED, driven by the monthly discharge and sediment load of the Yellow River in 1983 ～ 1984, was applied to the Bohai Sea. The simulation of the temporal-spatial variation of SPM （Suspended Particular Matter） concentration, sediment flux, and deposition rate was carfled out. It could be seen that, the SPM concentration exhibited distinct seasonal variation, and the variation characteristic was similar for two years. Some of the fine sediment of the Yellow River extended to the western bank of the Liaodong Bay, but most of the sediment deposited around the estuary of the Yellow River, the deposition flux rapidly decreasing with the increasing distance to the estuary. The deposition rate kept at 0.05 ～ 0.1 mm/a for most area, close to the insitu data. There were two areas of high SPM concentration in the Bohai Sea, one of which could be attributed to the Yellow River discharge, and another one, located at the seas with intensive resuspension power, to the tidal current energy. The sediment from the Yellow River was transported mainly along the track of the, Lagrangian residual cireulations. In one word, the suspended load transport was basically attributed to the wind-driven and tidal-induced Lagrangian residual circulations. The resuspenslon flux was related to the bottom currents： the larger the current speed, the more the sediment resuspended.