共查询到20条相似文献，搜索用时 62 毫秒
1.
S. W. Sloan 《国际地质力学数值与分析法杂志》1989,13(3):263-282
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. 相似文献
2.
Wang Jianfeng 《中国地质大学学报(英文版)》1998,9(2):2
ＩＮＴＲＯＤＵＣＴＩＯＮＤｅｓｐｉｔｅｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｅｌａｓｔｉｃ－ｐｌａｓｔｉｃ－ｖｉｓｃｏｕｓｎｕｍｅｒｉ－ｃａｌｍｅｔｈｏｄｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｓｌｏｐｅｓ,ｔｒａｄｉｔｉｏｎａｌｔｅｃｈｎｉｑｕｅｓｂａｓｅｄｏｎ... 相似文献
3.
In this paper, a nonlinear numerical technique is developed to calculate the limit load and failure mode of structures obeying an ellipsoid yield criterion by means of the kinematic limit theorem, nonlinear programming theory and displacement-based finite element method. Using an associated flow rule, a general yield criterion expressed by an ellipsoid equation can be directly introduced into the kinematic theorem of limit analysis. The yield surface is not linearized and instead a nonlinear purely kinematic formulation is obtained. The nonlinear formulation has a smaller number of constraints and requires less computational effort than a linear formulation. By applying the finite element method, the kinematic limit analysis with an ellipsoid yield criterion is formulated as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the dissipation power which is to be minimized and an upper bound to the plastic limit load of a structure can then be calculated by solving the minimum optimization problem. An effective, direct iterative algorithm has been developed to solve the resulting nonlinear programming formulation. The calculation is based purely on kinematically admissible velocities. The stress field does not need to be calculated and the failure mode of structures can be obtained. The proposed method can be used to calculate the bearing capacity of clay soils in a direct way. Some examples are given to illustrate the validity and effectiveness of the proposed method. 相似文献
4.
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. 相似文献
5.
将极限分析的上限定理、有限元离散思想、随机规划理论和蒙特卡洛方法这四者结合起来,提出了一种土质边坡可靠度分析的上限数值方法。首先采用三节点有限单元离散土质边坡,然后将土体的抗剪参数设为随机变量,根据上限定理构建同时满足三角形单元的塑性流动约束条件、单元公共边的塑性流动约束条件和单元速度边界条件的机动许可速度场,并根据内功功率等于外功功率条件建立目标函数,构建土质边坡可靠度分析的上限法随机规划模型。采用蒙特卡洛方法求解上限法随机规划模型,同时提出了一种基于上限法速度场的边坡失效风险系数估算方法,该方法特别适用于具有多种失效模式的边坡风险分析。对2个经典算例进行了深入分析,验证了方法的正确性。 相似文献
6.
极限分析是岩土工程稳定性评价的重要方法之一。传统的有限元极限分析方法,采用低阶三角形单元时需要引入速度间断面并采用特殊网格布局,或者采用高阶三角形单元等措施来克服体积锁定问题和提高数值精度。在光滑有限元法（smoothed finite element method,简称SFEM）的基础上,提出了一种基于新型混合常应力−光滑应变单元的极限分析方法（mixed constant stress-smoothed strain element limit analysis,简称MCSE-LA方法）。在服从关联流动法则和Mohr-Coulomb屈服准则的基础上,MCSE-LA方法最终将数值极限分析转化为以应力和极限荷载乘子为基本未知量的二阶锥规划（second order cone programming,简称SOCP）问题。MCSE-LA方法具有形式简单、优化变量相对较少和无需显式的写出塑性内能耗散函数的优点,并且根据凸锥优化的对偶理论,可以从对偶问题中获得速度场和塑性乘子等信息。此外,还采用基于最大塑性剪应变率的网格自适应加密算法,该算法在塑性区细化网格,显著提高了新数值极限分析方法的计算效率和精度。最后通过边坡稳定分析的结果对比,验证了MCSE-LA方法的计算精度和效率均高于传统的有限元极限分析方法。 相似文献
7.
A plane strain yield design approach to stability analysis of a column‐reinforced soil foundation under inclined loading 
下载免费PDF全文
下载免费PDF全文 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. 相似文献
8.
对屈服条件与屈服函数的极值条件进行了讨论,在不需要流动法则的情况下建立了速度方程。平衡方程、屈服条件与屈服函数的极值条件、速度方程就是土体极限分析的基本方程。这样,对荷载、速度边界条件（包括荷载与速度边界同时存在）,均构成了完备的极限平衡问题。对边坡稳定问题的上、下限定理给出了严密的证明,并建立了求解极限平衡问题的广义极限平衡法,均质土的计算结果表明,圆弧滑动面的广义极限平衡法与对数螺旋面的上限解法基本一致。 相似文献
9.
地基承载力问题是岩土工程领域最重要的研究课题之一。针对圆形浅基础地基,从空间问题着手,根据极限分析上限分析理论,选择合适的地基破坏模式及机动位移速度场,首次考虑了单元土体所受的侧向土压力对地基极限承载力的影响,同时考虑到土体自重及地面超载等因素的作用,推导出理论上更为合理的三维圆形浅基础地基极限承载力上限解。结合工程实例,用相应Matlab计算程序计算出上限解,并与实测值进行了对比,证明了其方法的合理性与实用性。 相似文献
10.
In this study, upper bound finite element (FE) limit analysis is applied to stability problems of slopes using a nonlinear criterion. After formulating the upper bound analysis as the dual form of a second-order cone programming (SOCP) problem, the stress field and corresponding shear strength parameters can be determined iteratively. Thus, the nonlinear failure criterion is represented by the shear strength parameters associated with stress so that the analysis of slope stability using a nonlinear failure criterion can be transformed into the traditional upper bound method with a linear Mohr–Coulomb failure criterion. Comparison with published solutions illustrates the accuracy and feasibility of the proposed method for a simple homogeneous slope stability problem. The proposed approach is also applied to a seismic stability problem for a rockfill dam to study the influence of different failure criterions on the upper bound solutions. The results show that the seismic stability coefficients obtained using two different nonlinear failure criteria are similar but that the convergence differs significantly. 相似文献
11.
S. W. Sloan 《国际地质力学数值与分析法杂志》1988,12(1):61-77
This paper describes a technique for computing lower bound limit loads in soil mechanics under conditions of plane strain. In order to invoke the lower bound theorem of classical plasticity theory, a perfectly plastic soil model is assumed, which may be either purely cohesive or cohesive-frictional, together with an associated flow rule. Using a suitable linear approximation of the yield surface, the procedure computes a statically admissible stress field via finite elements and linear programming. The stress field is modelled using linear 3-noded traingles and statically admissible stress discontinuities may occur at the edges of each triangle. Imposition of the stress-boundary, equilibrium and yield conditions leads to an expression for the collapse load which is maximized subject to a set of linear constraints on the nodal stresses. Since all of the requirements for a statically admissible solution are satisfied exactly (except for small round-off errors in the optimization computations), the solution obtained is a strict lower bound on the true collapse load and is therefore ‘safe’. A major drawback of the technique, as first described by Lysmer,1 is the large amount of computer time required to solve the linear programming problem. This paper shows that this limitation may be avoided by using an active set algorithm, rather than the traditional simplex or revised simplex strategies, to solve the resulting optimization problem. This is due to the nature of the constraint matrix, which is always very sparse and typically has many more rows that columns. It also proved that the procedure can, without modification, be used to derive strict lower bounds for a purely cohesive soil which has increasing strength with depth. This important class of problem is difficult to tackle using conventional methods. A number of examples are given to illustrate the effectiveness of the procedure. 相似文献
12.
This paper introduces sequential limit analysis (SLA) as a method for modelling large plastic deformations of purely cohesive materials such as undrained clay. The method involves solving a series of consecutive small‐deformation plastic collapse problems using finite element limit analysis, thus ensuring high levels of accuracy, efficiency, and robustness. The techniques needed to develop an SLA implementation for two‐dimensional (plane strain) problems are described in detail, including model geometry updating routines, treatment of rigid bodies, interfaces and boundaries, and periodic remeshing and interpolation of field variables. A simple total stress‐based constitutive model is used to account for strain softening and strain rate effects. Extensive verifications and validations are performed using analytical solutions and physical model test results, comparing both collapse loads and failure mechanisms, to demonstrate the effectiveness of the SLA approach. Additional solution quality checks on the bracketing discrepancy between lower‐bound and upper‐bound limit analysis solutions, and on the incompressibility of the rigid‐plastic material, are also presented. 相似文献
13.
为明确盾构施工掌子面滑移破坏机制并确定掌子面支护力的合理范围,基于滑移线理论和极限分析上限定理,利用空间离散技术提出了一种盾构施工掌子面三维滑移破裂模型。依据大主应力拱理论计算滑移区顶部竖向土压力值,并以此作为滑移破坏区上部的竖向荷载计算掌子面极限支护力。研究表明,土拱效应显著影响掌子面前方土体竖向应力的大小及分布规律;将本模型与已有研究方法进行比较,验证了本模型获取的掌子面极限支护力极限分析上限解在黏性土地层以及摩擦土地层中的适用性。同时本模型构建的掌子面破坏区域形态更加贴近离心试验结果与数值计算结果。 相似文献
14.
极限分析方法是土边坡稳定性分析的重要方法之一。刚体有限元上限法是其中的一类,此类方法仍旧存在一些关键问题需要完善。由于单元的刚性假设,系统的塑性变形内能耗散仅发生在单元间的界面上,故此类方法的性能主要取决于界面的布局,即采用非结构化三角形单元计算往往精度较差。为此,提出了基于滑动面摄动的刚体有限元上限法及临界滑动面的搜索方法。首先,在考虑刚体转动的基础上构造刚体有限元上限法的二阶锥规划模型,用于确定在给定试滑动面条件下的运动许可速度场。其次,将试滑动面的控制参数视为决策变量,建立搜索临界滑动面的非线性非凸优化问题模型,并采用非线性单纯形方法和粒子群方法求解此优化问题找出临界滑动面。通过经典边坡稳定问题的分析求解,验证了所提出的新方法,进一步证实了网格类型（即界面的布局）是影响刚体有限元上限法计算精度的主要因素。经过计算结果的对比发现,在刚体有限元上限法中考虑刚体转动是非常必要的,不仅可以提高刚体有限元上限法的计算精度,还可以克服此方法对界面布局的依赖性。 相似文献
15.
16.
This paper presents an inelastic element for the analysis of beams resting on two-parameter foundations. The element is derived from a two field mixed formulation with independent approximation of forces and displacements. The values for the two parameters of the foundation are derived through an iterative technique that is based on an assumption of plane strain for the soil medium. This iterative behavior is repeated at each time step of the nonlinear solution algorithm. The nonlinear response of structures resting on this improved two-parameter foundation model is analyzed following both a Vlasov and a Pasternak approach. Numerical examples that clarify the advantage of the newly developed model are conducted. These studies confirmed the importance of accounting for the foundation second parameter, and the efficiency and accuracy of the proposed model. 相似文献
17.
刚性块体极限分析上限法常应用于岩土工程稳定性研究,然而应用时需假定刚性块体破坏模式并递推繁琐的几何关系。为此,提出一种适应性更广的基于非线性规划模型的刚体平动运动单元上限有限元法,并解决了其优化模型初始值的确定问题。通过引入有限单元思想,将计算区域离散成刚体单元,同时以单元速度和节点坐标作为决策变量,由上限定理建立非线性规划模型获得上限解。利用编制的上限有限元程序进行边坡和浅埋隧道稳定性算例验证,表明运动单元上限有限元法能调整速度间断线至较优方位,所得破坏模式特征鲜明,上限解精度高,可广泛应用于边坡、隧道等稳定性分析研究。 相似文献
18.
针对地表超载作用下隧道稳定性和破坏模式问题,基于刚体平动运动单元上限有限元理论编程并计算分析,获得了浅埋隧道失稳临界超载系数上限解和刚性运动块体体系破坏模式。通过与现有的刚性块体极限分析上限法以及极限分析上、下限有限元法计算结果的对比分析,验证了上限解的可靠性。研究结果表明,（1）临界超载系数 黏聚力c之比 随土体内摩擦角 和隧道埋深C与直径D之比( )的增大而相应增大,随土体重度与黏聚力参数 的增大而减小;（2） 和 对隧道破坏模式的影响较明显; 增大,则隧道破坏范围增加;内摩擦角 增大,刚性运动块体破坏模式相互错动更加显著,相比而言, 对破坏模式的影响并不显著;（3）刚体平动运动单元上限有限元上限解精度高,所得刚性运动块体破坏模式具有滑移线形态,能精细地反映隧道失稳破坏特征。 相似文献
19.
考虑地基土体的非均质特性,采用非线性Mohr-Coulomb强度准则及其关联流动法则构造了浅埋水平条形锚板的曲线型破裂机制与机动许可速度场,根据极限分析上限定理推导了条形锚板抗拔承载力的表达式。利用变分极值原理求得了锚板抗拔承载力及其上方土体破裂面的上限解,分析了锚板埋深、土体非均质和非线性强度特性对锚板抗拔承载特性的影响,并将该上限解与已有计算方法进行了对比。结果表明:锚板埋深、土体非均质和非线性强度特性对其抗拔承载力与破裂面特征具有明显的影响。锚板埋深和土体非均质系数越大以及土体非线性强度系数越小,锚板抗拔承载力和土体破裂面深度、宽度均是越大。该上限解与极限平衡和极限分析有限元方法的计算结果基本一致,验证了所采用的曲线型破裂机制和地基非均质变化规律有效性,为条形锚板设计提供了一定的参考。 相似文献