基于振动台试验的不同含石率土-石混合体边坡地震动响应差异性研究

针对土-石混合体边坡动力响应特性尚不明确的问题,基于相似原理设计开展了3组不同含石率土-石混合体边坡大型振动台试验,系统地对比分析了不同含石率土-石混合体边坡动力响应的差异。研究结果表明:土-石混合体边坡在地震作用下的加速度响应符合自由面效应,即边坡在近坡顶端处的加速度放大效应要明显强于坡体内,而不同含石率土-石混合体边坡在不同频率正弦波激振作用下的动力响应差异明显,这是因为不同含石率土-石混合体边坡结构的动力特性存在差异;地震作用下,模型边坡动土压力由边坡浅表层至坡体内部呈不断增大趋势,但由于坡体变形与损伤程度不同,不同含石率边坡整体的动土压力响应程度不同;且在地震波分级加载过程中,边坡不同部位处动土压力的突变可作为边坡产生动力破坏的依据;随着含石率的上升,地震作用下边坡整体的变形程度逐渐降低:20%含石率边坡表现为由浅层至深层的连续滑动,而40%与60%含石率的边坡变形则相对较小,40%含石率边坡仅表层岩土体发生滑移,60%含石率边坡仅浅表层土体发生剥落,含石率的提升对土-石混合体边坡起到了加固作用。

Assessing the range of blasting-induced cracks in the surrounding rock of deeply buried tunnels based on the unified strength theory

Blasting-induced cracks in the rock surrounding deeply buried tunnels can result in water gushing and rock mass collapse,posing significant safety risks.However,previous theoretical studies on the range of blasting-induced cracks often ignore the impact of the in-situ stress,especially that of the intermediate principal stress.The particle displacement−crack radius relationship was established in this paper by utilizing the blasthole cavity expansion equation,and theoretical analytical formulas of the stress−displacement relationship and the crack radius were derived with unified strength theory to accurately assess the range of cracks in deep surrounding rock under a blasting load.Parameter analysis showed that the crushing zone size was positively correlated with in-situ stress,intermediate principal stress,and detonation pressure,whereas negatively correlated with Poisson ratio and decoupling coefficient.The dilatancy angle-crushing zone size relationship exhibited nonmonotonic behavior.The relationships in the crushing zone and the fracture zone exhibited opposite trends under the influence of only in-situ stress or intermediate principal stress.As the in-situ stress increased from 0 to 70 MPa,the rate of change in the crack range and the attenuation rate of the peak vibration velocity gradually slowed.

阶梯状滑动断续节理顺层边坡稳定性分析

阶梯状滑动顺层岩质滑坡在工程中较为常见,其岩桥和节理特性对阶梯状滑动面的形成至关重要。考虑岩桥和节理面之间强度参数弱化特性,构建阶梯状滑动顺层边坡稳定性分析模型。研究表明:与平面型滑动对比验证了分析模型的准确性。裂隙连通率k_(Y)、基本摩擦角φb、粗糙度系数JRC和弱化系数K_(c)对安全系数F_(s)影响显著。岩桥倾角β_(3)越大,F_(s)减小越显著,边坡越不稳定。台阶法中开挖次数和坡角角度对边坡稳定性影响显著,通过算例分析验证了该方法与极限平衡解差异率在1%以内。某些情况下当开挖使节理面出露时,F_(s)显著降低(如β=53°时降低了22.8%),易诱发滑坡的发生。同时,随着拉裂缝倾角β_(1)和充水高度h_(w)的增大,F_(s)显著降低,出流缝被堵塞时水力效应更为显著,此时边坡最不稳定。

Limit equilibrium method(LEM) of slope stability and calculation of comprehensive factor of safety with double strength-reduction technique

When the slope is in critical limit equilibrium(LE) state, the strength parameters have different contribution to each other on maintaining slope stability. That is to say that the strength parameters are not simultaneously reduced. Hence, the LE stress method is established to analyze the slope stability by employing the double strengthreduction(DSR) technique in this work. For calculation model of slope stability under the DSR technique, the general nonlinear Mohr–Coulomb(M–C) criterion is used to describe the shear failure of slope. Meanwhile, the average and polar diameter methods via the DSR technique are both adopted to calculate the comprehensive factor of safety(FOS) of slope. To extend the application of the polar diameter method, the original method is improved in the proposed method. After comparison and analysis on some slope examples, the proposed method's feasibility is verified. Thereafter, the stability charts of slope suitable for engineering application are drawn. Moreover, the studies show that:(1) the average method yields similar results as that of the polardiameter method;(2) compared with the traditional uniform strength-reduction(USR) technique, the slope stability obtained using the DSR techniquetends to be more unsafe; and(3) for a slope in the critical LE state, the strength parameter φ, i.e., internal friction angle, has greater contribution on the slope stability than the strength parameters c, i.e., cohesion.

Combined influence of nonlinearity and dilation on slope stability evaluated by upper-bound limit analysis

The combined influence of nonlinearity and dilation on slope stability was evaluated using the upper-bound limit analysis theorem.The mechanism of slope collapse was analyzed by dividing it into arbitrary discrete soil blocks with the nonlinear Mohr–Coulomb failure criterion and nonassociated flow rule.The multipoint tangent(multi-tangent) technique was used to analyze the slope stability by linearizing the nonlinear failure criterion.A general expression for the slope safety factor was derived based on the virtual work principle and the strength reduction technique,and the global slope safety factor can be obtained by the optimization method of nonlinear sequential quadratic programming.The results show better agreement with previous research result when the nonlinear failure criterion reduces to a linear failure criterion or the non-associated flow rule reduces to an associated flow rule,which demonstrates the rationality of the presented method.Slope safety factors calculated by the multi-tangent inclined-slices technique were smaller than those obtained by the traditional single-tangent inclined-slices technique.The results show that the multi-tangent inclined-slices technique is a safe and effective method of slope stability limit analysis.The combined effect of nonlinearity and dilation on slope stability was analyzed,and the parameter analysis indicates that nonlinearity and dilation have significant influence on the result of slope stability analysis.

Limit equilibrium analysis for stability of soil nailed slope and optimum design of soil nailing parameters

Reinforcement of slopes using soil nailing can effectively improve slope stability, and it has been widely used in upgrading cut slopes. Based on the assumptions of stresses on the slip surface, a new method for analyzing the stability of a slope reinforced with soil nails was established in the limit equilibrium theory framework, by considering that slope sliding occurs owing to shear failure of the slip surface, which subjects to Mohr–Coulomb(M–C) strength criterion. Meanwhile, in order to easily analyze the stability of a soil nailed slope in actual engineering and facilitate optimum design of parameters for soil nailing, factor of safety(FOS) contour curve charts were drawn on the basis of the established linear proportional relationship between the spacing of soil nails and slope height, and the length of soil nails and slope height. Then, by analyzing and verifying the results obtained from classic examples, some conclusions can be got as follows: 1) The results obtained from the current method are close to those obtained from the traditional limit equilibrium methods, and the current method can provide a strict solution for the slope FOS as it satisfies all the static equilibrium conditions of a sliding body, thus confirming the feasibility of the current method; 2) The slope FOS contour curve charts can be used not only to reliably analyze the stability of a soil nailed slope, but also to design optimally the parameters of soil nailing for the slope with a certain safety requirement.

Upper bound analysis of ultimate pullout capacity of shallow 3-D circular plate anchors based on nonlinear Mohr-Coulomb failure criterion

Based on the nonlinear Mohr-Coulomb failure criterion and the associated flow rules,the three-dimensional(3-D)axisymmetric failure mechanism of shallow horizontal circular plate anchors that are subjected to the ultimate pullout capacity(UPC)is determined.A derivative function of the projection function for projecting the 3-D axisymmetric failure surface on plane is deduced using the variation theory.By using difference principle,the primitive function of failure surface satisfying boundary condition and numerical solution to its corresponding ultimate pullout capacity function are obtained.The influences of nonlinear Mohr-Coulomb parameters on UPC and failure mechanism are studied.The result shows that UPC decreases with dimensionless parameter m and uniaxial tensile strength increases but increases when depth and radius of plate anchor,surface overload,initial cohesion,geomaterial density and friction angle increase.The failure surface is similar to a symmetrical spatial funnel,and its shape is mainly determined by dimensionless parameter m;the surface damage range expands with the increase of radius and depth of the plate anchor as well as initial cohesion but decreases with the increase of dimensionless parameter m and uniaxial tensile strength as well as geomaterial density.As the dimensionless parameter m=2.0,the numerical solution of UPC based on the difference principle is proved to be feasible and effective through the comparison with the exact solution.In addition,the comparison between solutions of UPC computed by variation method and those computed by upper bound method indicate that variation method outperforms upper bound method.

Variation analysis of ultimate pullout capacity of shallow horizontal strip anchor plate with 2-layer overlying soil based on nonlinear M-C failure criterion

Based on the nonlinear Mohr-Coulomb failure criterion and an associated flow rule,a kinematic admissible velocity field of failure mechanism of the 2-layer soil above a shallow horizontal strip anchor plate is constructed.The ultimate pull-out force and its corresponding failure mechanism through the upper bound limit analysis according to a variation principle are deduced.When the 2-layer overlying soil is degraded into single-layer soil,the model of ultimate pullout force could also be degraded into the model of single-layer soil.And the comparison between results of single-layer soil variation method and those calculated by rigid limit analysis method proves the correctness of our method.Based on that,the influence of changes of geotechnical parameters on ultimate pullout forces and failure mechanism of a shallow horizontal strip anchor with the 2-layer soil above are analyzed.The results show that the ultimate pull-out force and failure mechanism of a shallow horizontal strip anchor with the 2-layer soil above are affected by the nonlinear geotechnical parameters greatly.Thus,it is very important to obtain the accurate geotechnical parameters of 2-layer soil for the evaluation of the ultimate pullout capacity of the anchor plate.

Influence of non-dimensional strength parameters on the seismic stability of cracked slopes

Cracks in rock or soil slopes influence the stability and durability of the slopes. Seismic forces can trigger slope disasters, particularly in the cracked slopes. Considering the nonlinear characteristics of materials, the more generalized nonlinear failure criterion proposed by Baker is adopted. The influence of non-dimensional strength parameters on the stability of cracked slopes under earthquakes is performed using the upper bound limit analysis. The seismic displacement is calculated by adopting the logarithmic spiral failure surface according to the sliding rigid block model. Based on the existing studies, two methods for the stability analysis of cracked slopes under earthquakes are introduced: the pseudo-static method(with the factor of safety(Fs) as an evaluation index), and the displacement-based method(with the seismic displacement as an evaluation index). The pseudo-static method can only determine the instantaneous stability state of the cracked slope, yet the displacement-based methodreflects the stability variation of cracked slopes during earthquakes. The results indicate that the nondimensional strength parameters affect the factor of safety and seismic displacement of slopes significantly. The non-dimensional strength parameter(n) controlling the curvature of strength function shapes on the slope stability is affected by other parameters. Owing to cracks, the effect of non-dimensional strength parameters on seismic displacement becomes more significant.

Investigation on seismic response of a three-stage soil slope supported by anchor frame structure

Based on a typical prototype of a soil slope in engineering practice, a numerical model of a three-stage soil slope supported by the anchor frame structure was established by means of FLAC3D code. The dynamic responses of three-stage soil slope and frame structure were studied by performing a series of bidirectional Wenchuan motions in terms of the failure mode of three-stage structure, the acceleration of soil slope, the displacement of frame structure, and the anchor stress of frame structure. The response accelerations in both horizontal and vertical directions are the most largely amplified at the slope top of each stage subjected to different shaking cases. The platforms among the stages reduce the amplification effect of response acceleration. The residual displacement of frame structure increases significantly as the intensity of shaking case increases. The frame structure at each stage presents a combined displacement mode consisting of a translation and a rotation around the vertex. The anchor stress of frame structure is mainly increased by the first intense pulse of Wenchuan seismic wave, and it is sensitive to the intensity of shaking case. The anchor stress of frame structure at the first stage is the most considerably enlarged by earthquake loading.

Limit equilibrium stability analysis of slopes under external loads

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.

System reliability analysis of seismic pseudo-static stability of rock wedge based on nonlinear Barton–Bandis criterion

Based on the nonlinear Barton–Bandis(B–B)failure criterion,this study considers the system reliability of rock wedge stability under the pseudo-static seismic load.The failure probability(Pf)of the system is calculated based on the Monte−Carlo method when considering parameter correlation and variability.Parameter analysis and sensitivity analysis are carried out to explore the influence of parameters on reliability.The relationships among the failure probability,safety factor(Fs),and variation coefficient are explored,and then stability probability curves of the rock wedge under the pseudo-static seismic load are drawn.The results show that the parameter correlation of the B–B failure criterion has a significant influence on the failure probability,but correlation increases system reliability or decreases system reliability affected by other parameters.Under the pseudo-static seismic action,sliding on both planes is the main failure mode of wedge system.In addition,the parameters with relatively high sensitivity are two angles related to the joint dip.When the coefficient of variation is consistent,the probability of system failure is a function of the safety factor.

Barton-Bandis criterion-based system reliability analysis of rock slopes

Several potential failure modes generally exist in rock slopes because of the existence of massive structural planes in rock masses. A system reliability analyses method for rock slopes with multiple failure modes based on nonlinear Barton-Bandis failure criterion is proposed. The factors of safety associated with the sliding and overturning failure modes are derived, respectively. The validity of this method is verified through a planar rock slope with an inclined slope top and tension crack. Several sensitivity analyses are adopted to study the influences of structural-plane parameters, geometric parameters, anchoring parameters and fracture morphology on the rock slopes system reliability.

Water inflow forecasting for tunnel considering nonlinear variation of permeability coefficient

To assess the water inflow which is more suitable to the actual conditions of tunnel,an empirical correlation about the permeability coefficient changing with depth is introduced.Supposing that the surrounding rock is heterogeneous isotropy,the formula for calculating water inflow of tunnel with the nonlinear variation of permeability coefficient is deduced.By the contrast analysis with the existing formulas,the presented method has the similar value to them;moreover,the presented method has more simple form and easy to use.Due to parameter analysis,the water inflow decreases after considering the nonlinear variation of permeability coefficient.When the attenuation coefficient a>0,the water inflow increases first till reaches the maximum at a certain depth,then decreases and is close to 0 finally if deep enough.Thus,it is better to keep away from the certain depth where it is with the maximum water inflow for safe operation and economical construction,and reduce the water damage.Based on the analysis,the radius of tunnel has less impact on the amount of water inflow,and the water inflow just increases by 6.7% when the radius of tunnel increases by 1 m.

Limit variation analysis of shallow rectangular tunnels collapsing with double-layer rock mass based on a three-dimensional failure mechanism

In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.

一种可预测非饱和土壤中入渗和汇水的水力耦合模型

非饱和土壤在自然和工程实践中广泛存在,进一步了解非饱和多孔介质渗透过程中的水力耦合效应,对岩土工程和水文学界有着重要意义。本文提出一个非线性的完全水力耦合模型,用于控制可变形多孔介质中非饱和水的流动。该模型依赖于有限元(FE)求解器来计算受瞬态流作用下土壤中的水力耦合响应。使用已有的水力耦合解析解和积水入渗试验数据作为参考,对模型模拟结果的准确性进行验证。模拟结果表明,该模型在非稳定降雨条件和地表积水条件下能够准确模拟出降雨过程中的土壤入渗量和地表汇水量。此外,分别考虑吸力弹性模量、渗透模型的控制参数、饱和渗透系数以及边坡角度四个因素的变化而进行入渗和汇水过程的模拟结果进一步表明,土壤的膨胀性延迟了孔隙水压力的增长,而土壤的湿陷性则促进了孔隙水压力的增长;土壤孔径分布越粗糙,水力耦合效应越明显;土壤饱和渗透系数增大会导致更大的降雨入渗量和土体变形,并延缓地表开始积水的时间。边坡坡度对入渗过程的影响则存在争议,在本研究中,边坡坡度增大,土壤中的降雨入渗量减少,地表流失水增加,有待进一步探索。长远来看,在水力耦合系统中进一步考虑弹塑性耦合法则是更全面了解多孔介质中降雨入渗和地表汇水过程有用的方法。

Sensitivity analysis of regional rainfall-induced landslide based on UAV photogrammetry and LSTM neural network

Rainfall stands out as a critical trigger for landslides,particularly given the intense summer rainfall experienced in Zheduotang,a transitional zone from the southwest edge of Sichuan Basin to Qinghai Tibet Plateau.This area is characterized by adverse geological conditions such as rock piles,debris slopes and unstable slopes.Furthermore,due to the absence of historical rainfall records and landslide inventories,empirical methods are not applicable for the analysis of rainfall-induced landslides.Thus we employ a physically based landslide susceptibility analysis model by using highprecision unmanned aerial vehicle(UAV)photogrammetry,field boreholes and long short term memory(LSTM)neural network to obtain regional topography,soil properties,and rainfall parameters.We applied the Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability(TRIGRS)model to simulate the distribution of shallow landslides and variations in porewater pressure across the region under different rainfall intensities and three rainfall patterns(advanced,uniform,and delayed).The landslides caused by advanced rainfall pattern mostly occurred in the first 12 hours,but the landslides caused by delayed rainfall pattern mostly occurred in the last 12 hours.However,all the three rainfall patterns yielded landslide susceptibility zones categorized as high(1.16%),medium(8.06%),and low(90.78%).Furthermore,total precipitation with a rainfall intensity of 35 mm/h for 1 hour was less than that with a rainfall intensity of 1.775 mm/h for 24hours,but the areas with high and medium susceptibility increased by 3.1%.This study combines UAV photogrammetry and LSTM neural networks to obtain more accurate input data for the TRIGRS model,offering an effective approach for predicting rainfall-induced shallow landslides in regions lacking historical rainfall records and landslide inventories.

Fourier-shape-based reconstruction of rock joint profile with realistic unevenness and waviness features

Rock joint shape characteristics,waviness and unevenness play essential but distinct roles in shear mechanism of rock joints.This study presents a novel method to generate virtual rock joint profiles with realistic waviness and unevenness features.Firstly,joint profiles are obtained by 3D laser scanning device.Secondly,quantification of waviness and unevenness is conducted by traditional method,including digital filtering technique and roughness parameter RL.Thirdly,the discrete Fourier transform(DFT)method is employed to analyze the joint outlines.Two representative Fourier shape descriptors(D3,D8)for characterization of waviness and unevenness are suggested.Then,the inverse discrete Fourier transform(IDFT)is adopted to reconstruct the joint profiles with random values of phase angles but prescribed amplitudes controlled by D3 and D8.The traditional method is then applied to the reconstructed joint profiles to examine statistically the relationships between D3 and D8 and parameters RL of waviness and unevenness,respectively.The results show that larger D8 tends to result in larger waviness while higher D3 tends to increase unevenness.Reference charts for estimation of waviness and unevenness with different pairs of D3 and D8 are also provided to facilitate implementation of random joint reconstruction.

Seismic dynamic stability of double-slider rock slopes containing tension cracks

Earthquakes have significant impact on rock slopes,thus studying the seismic stability of double-slider rock slopes containing tension cracks is crucial.We proposed an analysis method on the seismic dynamic slope stability.This method utilizes discrete Fourier transform to decompose real earthquake waves into a combination of harmonic waves.These waves are then used in conjunction with the pseudo-dynamic method and safety factor calculation formula to compute the safety factor.This approach accurately captures the influence of seismic time history characteristics on the dynamic stability of double-slider rock slopes containing tension cracks.The minimum safety factor in the obtained time history curves of the safety factor reflects the most unfavorable state of the slopes under seismic effects.Quantitative analysis is conducted using six sets of actual earthquake ground motion data obtained from the Pacific Earthquake Engineering Research Center’s NGAWest2 ground-shaking record database.The conclusions are as follows:(1)There is an inverse correlation between the average seismic acceleration amplitude and the minimum safety factor.Conversely,the seismic acceleration amplitude standard deviation shows a positive correlation with the minimum safety factor.The global sensitivity of geometric parameters in the slope model is higher than other influencing factors.(2)The proposed dynamic stability analysis method can capture the dynamic characteristics of earthquakes,emphasizing the minimum safety factor of the slope in the seismic time history as a stability indicator.In contrast,the pseudo-static method may yield unsafe results.(3)A safety factor expression considering hydrostatic pressure is proposed.A negative correlation was observed between the height of the water level line and the minimum safety factor.