Bearing capacity of two interfering strip footings from lower bound finite elements limit analysis

A rigorous lower bound solution, with the usage of the finite elements limit analysis, has been obtained for finding the ultimate bearing capacity of two interfering strip footings placed on a sandy medium. Smooth as well as rough footing–soil interfaces are considered in the analysis. The failure load for an interfering footing becomes always greater than that for a single isolated footing. The effect of the interference on the failure load (i) for rough footings becomes greater than that for smooth footings, (ii) increases with an increase in ?, and (iii) becomes almost negligible beyond S/B > 3. Compared with various theoretical and experimental results reported in literature, the present analysis generally provides the lowest magnitude of the collapse load. Copyright © 2011 John Wiley & Sons, Ltd.     

无厚度界面单元边坡浅基础极限承载力上限值的数值分析

采用研发的有限单元分析程序,并同时考虑以库仑摩擦界面为组成模式的多段式无厚度界面单元,藉以探讨边坡坡角与坡缘距离对边坡坡顶浅基础极限承载力的影响,并应用于不同几何破坏机制条件,计算边坡基础在凝聚性土体的最小荷载上限值,模拟土体间的滑动破坏趋势与渐进式破坏行为。研究分析内容包含：（1）浅基础位于半无限空间;（2）在不同边坡坡角条件下,浅基础位于坡顶边缘处;（3）在不同边坡坡角条件下,浅基础位于不同坡缘距离等。经由数值计算分析结果显示,坡缘距离值愈小或边坡坡角愈陡峭时,边坡浅基础极限承载力值则愈小。换言之,当浅基础位于坡顶边缘处且坡角为垂直开挖情况时,则导致浅基础极限承载力降低约50 %。     

基于线性规划模型的极限分析上限有限元的实现

极限分析上限有限元法常利用三角形常应变率单元和摩尔-库仑屈服函数线性化的方法,以形成较易求解的线性规划模型。然而为满足计算精度,需引入大量的优化变量,增加了计算和存储的难度。为此,基于线性规划模型,利用MATLAB编制极限分析上限有限元程序,针对线性规划模型等式约束矩阵的高度稀疏性的特点,以稀疏矩阵的方式存储,从而在一定程度上解决了上述问题,使得上限有限元法能处理较大规模的岩土工程稳定性问题。以条形基础地基承载力课题为例进行算例分析,验证该方法的有效性,同时讨论了模型网格单元和塑性乘子数目对计算结果精度的影响。     

Bearing capacity of non-associative coaxial granular materials by upper bound limit analysis and finite elements

In this paper, an upper bound estimate of the limit load on non-associative coaxial granular materials is presented. The kinematic approach of the upper bound limit analysis has been utilised. The failure mechanism is assumed to coincide with the direction of the shear bands at every point throughout the body. The shear band orientation in non-associative coaxial materials, i.e. those with the same major principal stress and major principal strain increment directions, can be found based on the angle of dilation and the major principal stress direction. Therefore, having known the stress field at limiting equilibrium, the orientation of the shear bands and hence, the failure mechanism can be obtained. In this study, the stress field is first determined by the method of stress characteristics. Then, the finite element interpolation technique is used to interpolate the stress field and to find the orientation of the shear bands at every point within the field. Once the failure mechanism and the stress state at every point along velocity discontinuities have been found, the upper bound limit analysis has been performed to estimate the limit load.     

Collapse loads for rectangular foundations by three-dimensional upper bound limit analysis using radial point interpolation method

A three-dimensional kinematic limit analysis approach based on the radial point interpolation method (RPIM) has been used to compute collapse loads for rectangular foundations. The analysis is based on the Mohr-Coulomb yield criterion and the associated flow rule. It is understood that the internal plastic power dissipation function and flow rule constraints can be expressed entirely in terms of plastic strain rates without involving stresses. The optimization problem has been solved on basis of the semidefinite programming (SDP) by using highly efficient primal-dual interior point solver MOSEK in MATLAB. The results have been presented in terms of the variation of the shape factors with changes in the aspect ratio (L/B) of the footing for different values of soil internal friction angle (ϕ). Computations have revealed that the shape factors, s_c and s_q, due to effects of cohesion and surcharge increase continuously with (1) decrease in L/B and (2) increase in ϕ. On the other hand, the shape factor s_γ, due to the effect of soil unit weight, increases very marginally with an increase in L/B up to (1) ϕ = 25° for a rough footing and (2) ϕ = 35° for a smooth footing. Thereafter, for greater values of ϕ, the variation of s_γ with L/B has been found to be quite similar to that of the factors s_c and s_q. The variations of (1) nodal velocity patterns, (2) plastic power dissipation, and (3) maximum plastic shear strain rates have also been examined to interpret the associated failure mechanism.     

Bearing capacity factors of circular foundations for a general c–ϕ soil using lower bound finite elements limit analysis

By using the lower bound limit analysis in conjunction with finite elements and linear programming, the bearing capacity factors due to cohesion, surcharge and unit weight, respectively, have been computed for a circular footing with different values of ?. The recent axisymmetric formulation proposed by the authors under ?=0 condition, which is based on the concept that the magnitude of the hoop stress (σ_θ) remains closer to the least compressive normal stress (σ₃), is extended for a general c–? soil. The computational results are found to compare quite well with the available numerical results from literature. It is expected that the study will be useful for solving various axisymmetric geotechnical stability problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Stability determination for layered soil slopes using the upper bound limit analysis

By using the upper bound limit analysis, a method has been proposed in this article to compute stability numbers for layered soil slopes. The rupture surface was assumed to be a combination of different logarithmic spiral arcs with a common focus. The analysis ensures the kinematic admissibility of the collapse mechanism with respect to the rigid rotation of the bounded soil mass about the focus of the logarithmic spiral. For the sake of illustration, the computations, however, were exclusively performed only for a two layered soil slope. The effect of the pore water pressure and horizontal earthquake body forces was also incorporated in the computations. The computational procedure was validated by making a comparison of the obtained stability numbers with those reported in literature.     

边坡稳定极限分析的单元集成法

提出的单元集成法为采用塑性力学上限定理分析边坡的稳定性提供了一种通用性较强的手段。它采用类似于有限元网格划分的方式离散边坡计算区域,并设定一个机动许可的滑动机构（包括滑裂面位置和速度场,它与网格划分方式无关）,在此滑动机构下可以很方便地计算每个单元贡献的外力功率和内能耗散率。把所有单元的能量相加就是滑坡体的总能量。然后,根据上限定理求得安全系数,并通过非线性数学规划方法找到最小值。对于直线滑裂面或对数螺旋滑裂面的单滑块机制,几个典型算例说明了该方法的有效性。     

Upper bound limit analysis using linear finite elements and non‐linear programming

A new method for computing rigorous upper bounds on the limit loads for one‐, two‐ and three‐dimensional continua is described. The formulation is based on linear finite elements, permits kinematically admissible velocity discontinuities at all interelement boundaries, and furnishes a kinematically admissible velocity field by solving a non‐linear programming problem. In the latter, the objective function corresponds to the dissipated power (which is minimized) and the unknowns are subject to linear equality constraints as well as linear and non‐linear inequality constraints. Provided the yield surface is convex, the optimization problem generated by the upper bound method is also convex and can be solved efficiently by applying a two‐stage, quasi‐Newton scheme to the corresponding Kuhn–Tucker optimality conditions. A key advantage of this strategy is that its iteration count is largely independent of the mesh size. Since the formulation permits non‐linear constraints on the unknowns, no linearization of the yield surface is necessary and the modelling of three‐dimensional geometries presents no special difficulties. The utility of the proposed upper bound method is illustrated by applying it to a number of two‐ and three‐dimensional boundary value problems. For a variety of two‐dimensional cases, the new scheme is up to two orders of magnitude faster than an equivalent linear programming scheme which uses yield surface linearization. Copyright © 2001 John Wiley & Sons, Ltd.     

Seismic uplift capacity of inclined strip anchors in sand using upper bound limit analysis

Pseudo-static approach is adopted in this paper to determine the seismic uplift capacity of an inclined strip anchor using upper bound limit analysis. Two different failure mechanisms are considered to obtain the magnitudes of unit weight component of uplift factor f_γE for different values of soil friction angle, interface friction of anchor plate, anchor inclination, embedment ratio and horizontal seismic acceleration coefficient. The failure mechanism 1 consists of a triangular and quadrilateral rigid blocks; whereas the failure mechanism 2 comprises a logarithmic spiral failure zone with varied focus, sandwiched between a triangular and quadrilateral rigid blocks. It is observed that the magnitude of uplift factor f_γE decreases significantly with the increase in seismic acceleration but increases with the increase in embedment ratio and roughness of the anchor surface. However, a mixed trend in the values of f_γE can be observed for different inclination of the anchor, which is clearly discussed in this paper. The results are compared with the existing values in the literature and the significance of the present methodology for designing the inclined strip anchor is discussed.     

基于MATLAB的边坡稳定性极限上限分析程序开发

基于系统能耗最小原理的极限分析上限定理,利用MATLAB强大的数据处理功能,对边坡稳定系数Ncr采用序列二次规划法（SQP）进行优化计算,并与枚举法相结合,从而解决SQP算法求解优化解容易陷入局部最优解的问题。同时开发了一款简单实用的边坡稳定性分析软件,实现了计算结果与临界滑裂破坏面图形的可视化输出。与经典算例对比分析,计算结果具有较好的一致性,表明了程序的可靠性。分析了内摩擦角?、边坡倾斜角?、? 3个参数对边坡稳定系数Ncr的影响。计算结果表明：边坡稳定性系数Ncr对?、? 较为敏感,而对? 的敏感次之。其次,当?、? 比较接近或在小范围内波动时,对稳定性系数Ncr的影响特别敏感;当?、?值较小时,滑动面通过坡趾下方,而当?值比较大时,滑动面通过坡趾。     

Upper bound limit analysis using simplex strain elements and second‐order cone programming

In geomechanics, limit analysis provides a useful method for assessing the capacity of structures such as footings and retaining walls, and the stability of slopes and excavations. This paper presents a finite element implementation of the kinematic (or upper bound) theorem that is novel in two main respects. First, it is shown that conventional linear strain elements (6‐node triangle, 10‐node tetrahedron) are suitable for obtaining strict upper bounds even in the case of cohesive‐frictional materials, provided that the element sides are straight (or the faces planar) such that the strain field varies as a simplex. This is important because until now, the only way to obtain rigorous upper bounds has been to use constant strain elements combined with a discontinuous displacement field. It is well known (and confirmed here) that the accuracy of the latter approach is highly dependent on the alignment of the discontinuities, such that it can perform poorly if an unstructured mesh is employed. Second, the optimization of the displacement field is formulated as a standard second‐order cone programming (SOCP) problem. Using a state‐of‐the‐art SOCP code developed by researchers in mathematical programming, very large example problems are solved with outstanding speed. The examples concern plane strain and the Mohr–Coulomb criterion, but the same approach can be used in 3D with the Drucker–Prager criterion, and can readily be extended to other yield criteria having a similar conic quadratic form. Copyright © 2006 John Wiley & Sons, Ltd.     

基于极限分析的圆形浅基础地基承载力上限解

地基承载力问题是岩土工程领域最重要的研究课题之一。针对圆形浅基础地基,从空间问题着手,根据极限分析上限分析理论,选择合适的地基破坏模式及机动位移速度场,首次考虑了单元土体所受的侧向土压力对地基极限承载力的影响,同时考虑到土体自重及地面超载等因素的作用,推导出理论上更为合理的三维圆形浅基础地基极限承载力上限解。结合工程实例,用相应Matlab计算程序计算出上限解,并与实测值进行了对比,证明了其方法的合理性与实用性。     

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

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

边坡稳定分析的极限分析下限解有限元法

研制开发了极限分析下限解有限元程序.在程序的开发过程中,着重解决了下限解有限元法转化为标准线性规划问题和线性优化方法的选择问题,同时,通过典型算例的分析对比,对程序的合理性和可行性进行了验证.     

岩石边坡稳定的块体单元极限分析上限法

将块体单元法与极限分析上限法相结合,提出了岩石边坡稳定的块体单元极限分析上限法。借助于块体单元法离散计算区域,通过块体系统速度场在结构面上满足Mohr-Coulomb屈服条件、关联流动法则、速度边界条件以及块体系统虚功原理,构成求解边坡强度储备安全系数上限解的非线性规划模型,模型的求解采用复合形法。最后研究了一个岩质楔形体边坡算例和小湾水电站进水口边坡实例,通过将块体单元法上限分析计算结果与刚体极限平衡法相比较,验证了该方法的正确性。     

水平加筋地基极限承载力的极限上限分析法

水平加筋是广泛应用的软土地基处理方法之一,而水平加筋地基极限承载力的确定是其地基处理设计的重要依据。因此,首先结合条形基础下水平加筋地基的工程特点,在探讨其承载机制和破坏特点的基础上,考虑破坏间断面上筋材与地基土体的变形协调,构建出反映加筋参数变化影响的可变破坏模式及机动允许速度场;然后,在此基础上,通过重点研究筋材的能量耗散分析方法,并引入上限极限分析理论,建立出条形基础下水平加筋地基极限承载力分析模型,再引入序列二次规划优化分析理论,建立出条形基础下水平加筋地基极限承载力确定方法,它能充分反映加筋设计参数对地基破坏模式及承载力的影响;最后,通过试验结果与该方法及现有同类分析模型结果的比较分析,表明了该模型与方法的合理性与可行性。     

浅埋隧道掌子面稳定性二维自适应上限有限元分析

为研究浅埋隧道掌子面稳定性及获取精细化的破坏模式,提出了一种上限有限元非结构化网格自适应加密策略。以单元耗散能权重指标作为网格自适应加密评判准则,该策略同时兼顾了单元尺度与塑性应变。应用高阶的6节点三角形单元并建立上限有限元线性规划模型,以多次反复计算和网格加密的方式实现了二维自适应上限有限元分析并编制了计算程序。利用条形基础地基极限承载力课题,从上限解精度和网格加密形态方面验证了该程序的有效性。针对浅埋隧道掌子面稳定性问题,展开多参数条件下的自适应上限有限元计算,分析了网格加密过程中单元总数与上限解精度的关系,列出不同隧道埋深和内摩擦角对应的隧道掌子面稳定性临界值的上限解,揭示出掌子面稳定性变化规律及精细化的破坏模式。     

非线性破坏准则与岩土材料地基承载力研究

根据线性Mohr-Coulomb破坏准则,前人运用极限平衡法、滑移线理论或极限分析等方法求解地基承载力问题,但实际上岩土在剪切破坏过程中破坏准则具有非线性。因此,为了研究非线性破坏准则对地基承载力的影响,基于上限定理,通过“切线法”引进变量,根据能量耗散情况,将承载力问题转变为非线性规划问题,运用“序列二次规划算法”求出地基承载力的最优解。数值计算结果表明,当非线性破坏准则转变为线性破坏准则时,非线性参数对地基承载力有重要影响。     

上限法在加筋土结构物稳定性分析中的应用

将极限分析上限法用于加筋土结构物稳定性分析的关键是滑动面上内力功如何计算的问题。在总结了基于传统塑性理论上限法内力功的计算方法的基础上,将基于广义塑性理论的上限法应用到杨俊杰的内力功计算方法中。对比基于传统和广义塑性理论的上限解及离心试验结果,认为基于传统塑性理论的上限法应用简单,计算结果基本可靠;基于广义塑性理论的上限法虽然理论上比较完善,但当滑动面为曲线时,应用繁琐。