Seismic stability analysis of piled embankments: A multiphase approach

The seismic stability analysis of an embankment lying over a soft foundation soil reinforced by a group of vertical piles is performed within the framework of the upper bound kinematic approach of yield design. The analysis is based on a previously developed ‘multiphase’ model of the reinforced ground, which explicitly accounts for the shear and bending resistances of the piles. Making use of appropriate failure mechanisms involving shear zones across which the reinforcements are continuously deforming, along with ‘plastic hinge’ surfaces, upper bound estimates to the critical seismic coefficient of the structure are derived. The results, which are confirmed by the simulations obtained from a finite element elastoplastic code, give clear evidence of the key role played by the bending strength capacities of the piles in ensuring the stability of the pile reinforced embankment. Copyright © 2009 John Wiley & Sons, Ltd.     

FE analysis of grid reinforced embankment system on soft Bangkok clay

The behavior of a reinforced embankment on soft Bangkok clay has been analyzed by plane strain finite element method. The finite element analysis considers the selection of proper soil/reinforcement properties according to the relative displacement pattern of upper and lower interface elements. The large deformation phenomenon is simulated by updating the node coordinates, including those of the embankment elements above the current construction level, which ensures that the applied fill thickness simulates the actual field value. A full scale test reinforced embankment with a vertical face (wall) on Bangkok clay has been analyzed by the proposed finite element method, and the numerical results are compared with the field data. The response of a reinforced embankment on soft ground is principally controlled by the interaction between the reinforced soil mass and soft ground and the interaction between the grid reinforcement and the backfill soil. The tension in reinforcement and lateral displacement of the wall face varied during consolidation of foundation soil. The maximum tension force occurred in the reinforcement layer placed at the base of reinforced mass, due to bending of the reinforced mass resulting from differential settlements. It is considered necessary to account for the permeability variation of the soft ground foundation in the finite element analysis.     

考虑强度各向异性的边坡稳定有限元分析

天然沉积的土层总是表现出一定程度的强度各向异性,但现有的边坡稳定有限元分析极少考虑各向异性的影响。对大型有限元软件ABAQUS进行二次开发,使其能够考虑土体黏聚力随大主应力方向的变化,动态更新增量迭代求解过程中边坡不同位置处的抗剪强度,进而提出具备安全系数自动搜索功能的各向异性边坡稳定分析方法。计算结果表明,均匀边坡的有限元解与极限分析上限解相差很小。如果采用土体固结方向的黏聚力并按各向同性评价缓坡的稳定性,可能严重高估安全系数,尤其是在黏聚力较高的情况下。与极限分析不同,建立的强度更新有限元模型能够分析成层边坡的稳定性。     

A plane strain yield design approach to stability analysis of a column‐reinforced soil foundation under inclined loading

The ultimate bearing capacity problem of column‐reinforced foundations under inclined loading is investigated within the framework of static and kinematic approaches of yield design theory. The configuration of a native soft clayey soil reinforced by either a group of purely cohesive columns (lime‐column technique) or a group of purely frictional columns (stone‐column technique) is analyzed under plane strain conditions. First, lower bound estimates are derived for the ultimate bearing capacity by considering statically admissible piecewise linear stress distributions that comply with the local strength conditions of the constitutive materials. The problem is then handled by means of the yield design kinematic approach of limit analysis through the implementation of several failure mechanisms, allowing the formulation of upper bound estimates for the ultimate bearing capacity. A series of finite element limit load solutions obtained from numerical elastoplastic simulations suggests that the predictions derived from the kinematic approach appear to be more accurate than the estimates obtained from the static approach. Comparison with available results obtained in the context of yield design homogenization demonstrates the accuracy of the proposed direct analysis, which may therefore be viewed as complementary approach to homogenization‐based approaches when a small number of columns is involved. Copyright © 2016 John Wiley & Sons, Ltd.     

基于塑性极限分析上限法理论的土质边坡可靠度分析

将极限分析的上限定理、有限元离散思想、随机规划理论和蒙特卡洛方法这四者结合起来,提出了一种土质边坡可靠度分析的上限数值方法。首先采用三节点有限单元离散土质边坡,然后将土体的抗剪参数设为随机变量,根据上限定理构建同时满足三角形单元的塑性流动约束条件、单元公共边的塑性流动约束条件和单元速度边界条件的机动许可速度场,并根据内功功率等于外功功率条件建立目标函数,构建土质边坡可靠度分析的上限法随机规划模型。采用蒙特卡洛方法求解上限法随机规划模型,同时提出了一种基于上限法速度场的边坡失效风险系数估算方法,该方法特别适用于具有多种失效模式的边坡风险分析。对2个经典算例进行了深入分析,验证了方法的正确性。     

非均质土坡的有限元塑性极限分析

在评价现有边坡稳定分析方法的基础上,提出了边坡稳定的有限元塑性极限分析方法。借助有限单元思想和线性规划,建立了边坡稳定的数学规划模型,由此可以求出安全系数的上下限解,从而,界定了边坡的安全范围,同时,可以给出下限状态下的应力场和上限状态下的速度场。以一个经典非均质土坡的边坡稳定作为算例,比较了多种方法的分析结果,论证了本方法的正确性。     

不排水双层粘性土地基极限承载力的数值分析

采用弹塑性有限元分析了条形基础作用下不排水条件的双层粘性土地基极限承载力性状。采用修正的地基承载力系数表征,并将不同的几何与土层参数条件下的数值解与上下限解和经典的经验解进行比较。表明弹塑性位移有限元法可以很好地求解地基的极限承载力问题,其求解得到的修正地基承载力系数与基于下限原理的有限元解很接近,而上限解高估了地基的极限荷载值,传统的经验解在某些条件下却偏小。     

Slope stability analysis by means of finite element limit analysis and finite element strength reduction techniques. Part II: Back analyses of a case history

This paper deals with a back analysis of a slope failure. The case history investigated is located in an alpine environment in central Europe and is characterized by a very steep original terrain, indicating in situ soil with high strength. To study the factor of safety, two different approaches applying the so-called φ′/c′ reduction are used, namely finite element limit analysis and strength reduction finite element analysis. Comparison of a strength reduction technique with rigorous finite element limit analysis confirms that the factors of safety (FoS) obtained are very similar for associated plasticity, an intrinsic assumption of limit analysis. For non-associated plasticity, a modified version of the so-called Davis approach has been applied because it has been shown that the original formulation proposed by Davis works well when the FoS is defined in terms of loads but is not appropriate when the FoS is defined in terms of soil strength. The results show that, with the modified Davis parameter, both strength reduction finite element analyses and finite element limit analyses provide very similar factors of safety. The key advantage of limit analysis, however, is that the value of the FoS can be bracketed from above and below with upper and lower bound calculations.     

Upper bound limit analysis using finite elements and linear programming

This paper describes a technique for computing rigorous upper bounds on limit loads under conditions of plane strain. The method assumes a perfectly plastic soil model, which is either purely cohesive or cohesive-frictional, and employs finite elements in conjunction with the upper bound theorem of classical plasticity theory. The computational procedure uses three-noded triangular elements with the unknown velocities as the nodal variables. An additional set of unknowns, the plastic multiplier rates, is associated with each element. Kinematically admissible velocity discontinuities are permitted along specified planes within the grid. The finite element formulation of the upper bound theorem leads to a classical linear programming problem where the objective function, which is to be minimized, corresponds to the dissipated power and is expressed in terms of the velocities and plastic multiplier rates. The unknowns are subject to a set of linear constraints arising from the imposition of the flow rulé and velocity boundary conditions. It is shown that the upper bound optimization problem may be solved efficiently by applying an active set algorithm to the dual linear programming problem. Since the computed velocity field satisfies all the conditions of the upper bound theorem, the corresponding limit load is a strict upper bound on the true limit load. Other advantages include the ability to deal with complicated loading, complex geometry and a variety of boundary conditions. Several examples are given to illustrate the effectiveness of the procedure.     

地表超载作用下隧道失稳破坏的上限有限元分析

针对地表超载作用下隧道稳定性和破坏模式问题,基于刚体平动运动单元上限有限元理论编程并计算分析,获得了浅埋隧道失稳临界超载系数上限解和刚性运动块体体系破坏模式。通过与现有的刚性块体极限分析上限法以及极限分析上、下限有限元法计算结果的对比分析,验证了上限解的可靠性。研究结果表明,（1）临界超载系数 黏聚力c之比 随土体内摩擦角 和隧道埋深C与直径D之比( )的增大而相应增大,随土体重度与黏聚力参数 的增大而减小;（2） 和 对隧道破坏模式的影响较明显; 增大,则隧道破坏范围增加;内摩擦角 增大,刚性运动块体破坏模式相互错动更加显著,相比而言, 对破坏模式的影响并不显著;（3）刚体平动运动单元上限有限元上限解精度高,所得刚性运动块体破坏模式具有滑移线形态,能精细地反映隧道失稳破坏特征。     

Undrained stability of wide rectangular tunnels

Upper bound rigid block methods and finite element limit analysis are applied to investigate the undrained stability of rectangular tunnels. Solutions for a range of soil properties and tunnel geometries are obtained for situations in which the tunnel is wider than it is high. The upper and lower solutions for tunnel stability, which accurately bracket the true collapse load, are summarised in a series of stability charts for use in tunnel design.     

A homogenization method for estimating the bearing capacity of soils reinforced by columns

The ultimate bearing capacity problem of a strip foundation resting on a soil reinforced by a group of regularly spaced columns is investigated in the situation when both the native soil and reinforcing material are purely cohesive. Making use of the yield design homogenization approach, it is shown that such a problem may be dealt with as a plane strain yield design problem, provided that the reinforced soil macroscopic strength condition has been previously determined. Lower and upper bound estimates for such a macroscopic criterion are obtained, thus giving evidence of the reinforced soil strong anisotropy. Performing the upper bound kinematic approach on the homogenized bearing capacity problem, by using the classical Prandtl's failure mechanism, makes it then possible to derive analytical upper bound estimates for the reinforced foundation bearing capacity, as a function of the reinforced soil parameters (volume fraction and cohesion ratio), as well as of the relative extension of the reinforced area. It is shown in particular that such an estimate is closer to the exact value of the ultimate bearing capacity, than that derived from a direct analysis which implicitly assumes that the reinforced soil is an isotropic material. Copyright © 2005 John Wiley & Sons, Ltd.     

基于刚体平动运动单元的上限有限元研究

刚性块体极限分析上限法常应用于岩土工程稳定性研究,然而应用时需假定刚性块体破坏模式并递推繁琐的几何关系。为此,提出一种适应性更广的基于非线性规划模型的刚体平动运动单元上限有限元法,并解决了其优化模型初始值的确定问题。通过引入有限单元思想,将计算区域离散成刚体单元,同时以单元速度和节点坐标作为决策变量,由上限定理建立非线性规划模型获得上限解。利用编制的上限有限元程序进行边坡和浅埋隧道稳定性算例验证,表明运动单元上限有限元法能调整速度间断线至较优方位,所得破坏模式特征鲜明,上限解精度高,可广泛应用于边坡、隧道等稳定性分析研究。     

边坡稳定分析刚体有限元上限法的锥规划模型

极限分析方法是土边坡稳定性分析的重要方法之一。刚体有限元上限法是其中的一类,此类方法仍旧存在一些关键问题需要完善。由于单元的刚性假设,系统的塑性变形内能耗散仅发生在单元间的界面上,故此类方法的性能主要取决于界面的布局,即采用非结构化三角形单元计算往往精度较差。为此,提出了基于滑动面摄动的刚体有限元上限法及临界滑动面的搜索方法。首先,在考虑刚体转动的基础上构造刚体有限元上限法的二阶锥规划模型,用于确定在给定试滑动面条件下的运动许可速度场。其次,将试滑动面的控制参数视为决策变量,建立搜索临界滑动面的非线性非凸优化问题模型,并采用非线性单纯形方法和粒子群方法求解此优化问题找出临界滑动面。通过经典边坡稳定问题的分析求解,验证了所提出的新方法,进一步证实了网格类型（即界面的布局）是影响刚体有限元上限法计算精度的主要因素。经过计算结果的对比发现,在刚体有限元上限法中考虑刚体转动是非常必要的,不仅可以提高刚体有限元上限法的计算精度,还可以克服此方法对界面布局的依赖性。     

水平受荷桩的弹性有限元虚拟加载上限分析

基于将塑性上限分析等效为弹性迭代计算的总量虚拟加载上限分析理论,在商业化有限元软件ABAQUS中实现了弹性有限元虚拟加载上限方法（弹性有限元T-EMSD）。应用弹性有限元T-EMSD法分析了不排水黏土中的二维水平受荷桩,其获得的荷载-位移曲线与弹塑性有限元分析结果一致,其极限承载力与塑性解相近。在极限位移加载量下弹性有限元T-EMSD法对应的上限机构从弹性始速度场开始随迭代逐渐演化,迭代收敛后的速度场和解析塑性破坏机构相似。与其他基于可变强度概念（MSD）的方法相比,弹性有限元T-EMSD法对水平受荷桩桩身的分析具有更高的精度。弹性有限元T-EMSD法最大的优势在于可在计算中自然地获得塑性机构,因而可被用于研究一些塑性机构难以构造的复杂问题,并对弹塑性数值方法进行验证。     

Parametric studies of disturbed rock slope stability based on finite element limit analysis methods

This study employs the finite element upper bound and lower bound limit analysis methods to investigate the stability of inhomogeneous rock slopes. The differences in the stability numbers of the upper and lower bound solutions are bracketed within ±10.5% or better, and the stability numbers obtained are presented in rock slope stability charts. These stability charts can provide a convenient tool for preliminary stability designs of inhomogeneous rock slopes. Various recommended blasting damage zones are considered, and disturbance factors are used to represent damage levels. Results showed that rock mass disturbance could significantly influence the evaluation of rock slope stability.     

土体极限分析的基本方程与广义极限平衡法

对屈服条件与屈服函数的极值条件进行了讨论,在不需要流动法则的情况下建立了速度方程。平衡方程、屈服条件与屈服函数的极值条件、速度方程就是土体极限分析的基本方程。这样,对荷载、速度边界条件（包括荷载与速度边界同时存在）,均构成了完备的极限平衡问题。对边坡稳定问题的上、下限定理给出了严密的证明,并建立了求解极限平衡问题的广义极限平衡法,均质土的计算结果表明,圆弧滑动面的广义极限平衡法与对数螺旋面的上限解法基本一致。     

三维塑性极限分析下限法原理及应用

介绍了三维塑性极限分析下限法的基本原理。将计算区域离散为有限个三维空间单元,利用极限分析中的下限定理,借助有限元思想建立数学规划模型,并引入非线性规划的数值方法寻求问题的下限解。以均匀土质边坡的极限荷载作为算例,比较了所得数值解与经典塑性力学理论解,论证了该方法的正确性。     

横观各向同性软土地基上大圆筒防波堤稳定性研究

对于横观各向同性软土地基上沉入式大圆筒防波堤结构提出了一种准三维上限极限分析方法,所假设的破坏机制为大圆筒结构绕筒体内中轴线上某点发生转动失稳,泥面处形成楔体破坏,而筒底部形成圆弧滑裂面。本方法可以考虑土的三轴压缩、拉伸强度与直剪强度的差异。利用ABAQUS分别进行平面应变以及三维有限元分析,软土采用Hill本构模型,所得到的破坏模式以及大圆筒结构水平承载力与上限极限分析方法吻合较好,同时可以得出考虑地基土各向异性的大圆筒结构承载力比不考虑时有较大降低。     

Analysis of geosynthetic reinforced soil structures with orthogonal anisotropy

Reinforced soil with geosynthetics as a composite material represents significant orthogonally anisotropic properties. However, current analytical methods usually treat the soil and reinforcement separately, which is not true of practical situations. Therefore, it is difficult to use these methods to study the real effects of the reinforcement. This paper presents an analytical model based on the theory of elasticity for orthogonally anisotropic materials that can be used in analyzing reinforced soil structures with geosynthetics. The stresses and deformations at any point within the reinforced soil structure can be determined by the proposed model. The capabilities of the model have been illustrated by application to an example problem involving a physical model test of a geosynthetic reinforced soil structure. The results of the model prediction are compared with those obtained from the model tests as well as finite element analysis. It is shown that the results of the analytical solution are in good agreement with those of the physical model tests and the finite element analysis.