Force analysis of pile foundation in rock slope based on upper-bound theorem of limit

Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex load condition was applied and the upper-bound theorem of limit analysis was used to compute thrust of rock layers with all possible distribution shapes. The interaction of slope and pile was considered design load in terms of slope thrust,and the finite difference method was derived to calculate inner-force and displacement of bridge pile foundation in rock slope under complex load condition. The result of example shows that the distribution model of slope thrust has certain impact on displacement and inner-force of bridge pile foundation. The maximum displacement growth rate reaches 54% and the maximum moment and shear growth rates reach only 15% and 20%,respectively,but the trends of inner-force and displacement of bridge pile foundation are basically the same as those of the conventional pile foundation in the flat ground. When the piles bear the same level lateral thrust,the distribution shapes of slope thrust have different influence on inner-force of pile foundation,especially the rectangle distribution,and the triangle thrust has the smallest displacement and inner-force of pile foundation.

Upper-bound limit analysis based on the natural element method

The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.

Upper bound limit analysis of roof collapse in shallow tunnels with arbitrary cross sections under condition of seepage force

The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation principle according to Hoek-Brown failure criterion. The seepage force was included in the upper bound limit analysis, and it was computed from the gradient of excess pore pressure distribution. The seepage was regarded as a work rate of external force. The numerical results of roof collapse in square and circular tunnels with different rock parameters were derived and discussed, which proves to be valid in comparison with the previous work. The influences of different parameters on the shape of collapsing blocks were also discussed.

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.

NONSMOOTH MODEL FOR PLASTIC LIMIT ANALYSIS AND ITS SMOOTHING ALGORITHM

By means of Lagrange duality of Hill＇s maximum plastic work principle theory of the convex program, a dual problem under Mises＇ yield condition has been derived and whereby a non-differentiable convex optimization model for the limit analysis is developed. With this model, it is not necessary to linearize the yield condition and its discrete form becomes a minimization problem of the sum of Euclidean norms subject to linear constraints. Aimed at resolving the non-differentiability of Euclidean norms, a smoothing algorithm for the limit analysis of perfect-plastic continuum media is proposed. Its efficiency is demonstrated by computing the limit load factor and the collapse state for some plane stress and plain strain problems.

Upper Bound Limit Analysis of Anisotropic Soils

In this paper,a novel discretization method inσ-τspace is developed to calculate the upper bound limit loads and failure modes of anisotropic Mohr-Coulomb materials.To achieve this objective,the Mohr-Coulomb yield criterion is linearized inσ-τspace,which allows for upper bound solution of soils whose cohesion and friction coefficient varying with direction.The finite element upper limit analysis formulation using the modified anisotropic yield criterion is then developed.Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous upper bounds for anisotropic soils.

Upper bound analysis of slope stability with nonlinear failure criterion based on strength reduction technique

Based on the upper bound limit analysis theorem and the shear strength reduction technique,the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and its corresponding critical failure mechanism by means of the kinematical approach of limit analysis theory.The nonlinear shear strength parameters were treated as variable parameters and a kinematically admissible failure mechanism was considered for calculation schemes.The iterative optimization method was adopted to obtain the safety factors.Case study and comparative analysis show that solutions presented here agree with available predictions when nonlinear criterion reduces to linear criterion,and the validity of present method could be illuminated.From the numerical results,it can also be seen that nonlinear parameter m,slope foot gradient β,height of slope H,slope top gradient α and soil bulk density γ have significant effects on the safety factor of the slope.

A new failure mechanism for deep cavity and upper bound solution of supporting pressure

The investigation of supporting pressure is of great significance to the design of underground structures.Based on the kinematical approach of limit analysis,an improved failure mechanism is proposed,and the supporting pressure is investigated for deep buried cavity.Three failure mechanisms are first introduced according to the existing failure mechanisms of geotechnical structures of limit analysis.A comparison with respect to the optimal failure mechanisms and the upper bound solutions provided among these three mechanisms are then conducted in an attempt to obtain the improved failure mechanism.The results provided by the improved failure mechanism are in good agreement with those by the existing method,the numerical solution and field monitoring,which demonstrates that the proposed failure mechanism is effective for the upper bound analysis of supporting pressure.

Upper bound analytical solution for surrounding rock pressure of shallow unsymmetrical loading tunnels

By combining the results of laboratory model tests with relevant flow rules, the failure mode of shallow unsymmetrical loading tunnels and the corresponding velocity field were established. According to the principle of virtual power, the upper bound solution for surrounding rock pressure of shallow unsymmetrical loading tunnel was derived and verified by an example. The results indicate that the calculated results of the derived upper bound method for surrounding rock pressure of shallow unsymmetrical loading tunnels are relatively close to those of the existing "code method" and test results, which means that the proposed method is feasible. The current code method underestimates the unsymmetrical loading feature of surrounding rock pressure of shallow unsymmetrical loading tunnels, so it is unsafe; when the burial depth is less or greater than two times of the tunnel span and the unsymmetrical loading angle is less than 45°, the upper bound method or the average value of the results calculated by the upper bound method and code method respectively, is comparatively reasonable. When the burial depth is greater than two times of the tunnel span and the unsymmetrical loading angle is greater than 45°, the code method is more suitable.

The study of slip line field and upper bound method based on associated flow and non-associated flow rules

At present,associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle,so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory,the theory of slip line field is put forward in this paper,and then the ultimate bearing capacity of strip footing is obtained based on the associated flow rule and the non-associated flow rule individually. These two results are identical since the ultimate bearing capacity is independent of flow rule. On the contrary,the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associated flow rule is incorrect.

PLASTIC LIMIT ANALYSIS OF DISCONTINUOUS LINING UNDERGROUND PRESSURE

Discontinuous lining is a special form of support in underground excavation. Based on the method of plastic limit analysis, it is found the upper and the lower bound solution of the pressure of circular discontinuous lining and discussed support parameter of discontinuous lining and its applicable conditions , which provides theoretical basis for the design and calculation of discontinuous lining.

RIGID FINITE ELEMENT AND LIMIT ANALYSIS

According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.

Stability analysis for natural slope by kinematical approach

The stability of natural slope was analyzed on the basis of limit analysis. The sliding model of a kind of natural slope was presented. A new kinematically admissible velocity field for the new sliding model was constructed. The stability factor formulation by the upper bound theorem leads to a classical nonlinear programming problem, when the external work rate and internal energy dissipation were solved, and the constraint condition of the programming problem was given. The upper bound optimization problem can be solved efficiently by applying a nonlinear SQP algorithm, and stability factor was obtained, which agrees well with previous achievements.

考虑水力效应的裂缝边坡稳定性非线性能耗分析

现有考虑水力效应的裂缝边坡稳定性分析大多基于线性破坏准则,而岩土体破坏往往呈现非线性特征,因此开展水力效应影响下的裂缝边坡稳定性非线性极限上限分析具有重要意义。基于极限分析上限定理及强度折减技术,结合“外切线法”非线性破坏准则,构建坡顶含竖直裂缝的边坡对数螺旋线破坏模式,根据虚功原理推导出裂缝边坡安全系数解析式,通过MATLAB优化计算,结合边坡工程实例探讨了典型因素对裂缝边坡稳定性、临界裂缝及滑动面位置的影响规律。研究结果表明:随着地下水位h的持续上升,边坡安全系数不断降低,临界裂缝深度逐渐增大且临界裂缝位置逐渐向坡顶缘偏移;非线性系数m明显影响边坡的稳定性,边坡安全系数随着非线性系数的增大显著减小,采用线性破坏准则会高估地下水位变化对边坡稳定性的影响;随着非线性系数的增大,临界裂缝深度随之增大,临界裂缝位置距离坡顶缘越来越远,滑坡体体积逐渐增大;裂缝与坡顶缘距离l m随着非线性系数m的增大而增大(当m增大0.2时,l m增大约1 m),且随着裂缝与坡顶缘距离的增大,边坡安全系数无明显变化,临界裂缝深度先逐渐减小后趋于稳定。

大截面矩形顶管施工对既有管线影响研究

顶管施工会扰动邻近建筑物,危及工程安全。该文以苏州市城北路综合管廊元和塘段顶管工程为背景,通过ABAQUS有限元分析软件建立三维模型,研究顶管下穿地下管线施工对管周土及既有管线影响,并基于极限分析方法探讨了管线失稳模式。结果表明:大直径矩形顶管顶进方向的影响范围约为2倍顶管宽度,两侧的影响范围约为1倍顶管宽度,地下管线到顶管轴线的水平距离为0.5倍~1倍顶管宽度时,应注意管线受到的水平方向附加荷载是否会导致管线破坏;管线与顶管相对位置对管线破坏影响较为显著,越靠近顶管的管线破坏越严重,滑移线由顶管两侧经由管线向地面发展,破坏部位由管线顶部两侧延伸至底部;除受与顶管相对位置影响外,管线破坏模式还受矩形顶管形状影响,管线顶部及靠近顶管顶角区域破坏风险最大。

考虑基质吸力影响的非饱和路堤三维稳定性上限分析

采用极限上限分析定理构建了路堤失稳的三维旋转破坏机构,继而根据Bishop非饱和土抗剪强度理论,考虑非饱路堤内部基质吸力的空间分布及其随地下水位的变化,建立了路堤破坏土体外力功率与内部能量耗散功率的能量守恒方程,并利用遗传算法编写了非饱和路堤最小上限解的高效搜索算法。通过将非饱和土路堤基底破坏模式退化为边坡破坏模式,并与现有非饱和土边坡稳定性计算结果对比,验证了所提上限解的正确性和遗传搜寻算法的准确性。进一步地,对路堤三维稳定性的关键影响因素展开了系统分析,研究了路堤填土孔径分布、基质吸力、进气值倒数、路堤倾角以及有效内摩擦角等因素对路堤三维稳定性的影响规律。研究表明,非饱和土质路堤的稳定性不仅取决于路堤填土性质,而且依赖于影响土体吸力大小与分布的填土孔径参数和进气值等因素。地下水位升降引起的基质吸力变化对路堤稳定性存在显著影响。研究结果为路堤稳定性精细化分析提供了重要的理论依据。

海啸作用下滨水挡土墙抗震滑动稳定性上限分析

针对沿海地区滨水挡土墙在海啸和地震联合作用下易发生被动破坏的问题,结合极限分析上限法和拟静力法确定地震荷载,研究了海啸作用下滨水挡土墙的抗震滑动稳定性。通过建立墙体滑动与破坏区土楔体整体滑动的墙-土系统及速度容许场,推导出水平地震加速度系数(kh)的表达式。分析了海啸、海平面与地下水高度及墙-地和墙-土摩擦角对地震屈服加速度系数的影响。结果表明,控制地下水位在适当高度并增大墙-土摩擦角可以有效提高滨水挡土墙的抗震滑动稳定性。与经典极限平衡理论相比,该方法考虑了土体的塑性流动及挡土墙的位移模式,且无需计算地震被动土压力的值。

地震作用下双层土质边坡稳定性上限分析

为探究地震作用对双层土质边坡稳定性的影响,基于极限分析上限理论,建立了地震作用下双层土质边坡旋转破坏机构。采用拟静力法简化地震荷载作用,推导了双层土质边坡安全系数表达式。分析土体黏聚力比、边坡坡角以及地震作用对双层土质边坡稳定性及其滑动面参数的影响。研究结果表明:在水平地震作用下,随着黏聚力比的增大,边坡安全系数呈非线性增大且潜在滑动面变浅;当水平地震作用增强时,边坡失稳破坏范围显著加大;随着竖向地震系数及黏聚力比的增大,边坡安全系数略微降低;水平地震作用增强时,竖向地震作用对边坡稳定性的影响逐渐增大,此时需要考虑竖向地震的影响。