兰州交通大学博文学院毕业设计（论文）
collapse loads for earth slopes, where soils are not heavily constrained, are quite insensitive to whether the flow rule is associated or non-associated. Principle of Virtual Work
Both the lower-and upper –bound theorems are based on the principle of virtual work. The virtual work equation is applicable, given the assumption of small deformations before collapse, and can be expressed as either
?TvdS??XvdV?????(???p)?dVsVVV'ijAijBijABii?ABiiABijij?dV （3）
Or
?TSiABi'BvdS??XiAviBdV???ij?ijdV （4）
VVWhere Ti= boundary loadings; XijA= body forces not including seepage and buoyancy forces; XijA= body forces including seepage and buoyancy forces;?ijA= total stress tensor in
'Aequilibrium with Ti and XijA; ?ij= effective stress tensor in equilibrium with TiA
A?A?Band XijA; ?ij= Kronecker delta; p = pore-water pressure; and ?ij = strain rate tensor
compatible with the velocity field viB.
B There is no need for ?, Ti, and XijA to be related to ?ij and viB in any particular
AijA?way for (3) or (4) represent the rate of the external work, while the right-hand sides represent the rate of the internal power dissipation done by the assumed stress field and external loads on the assumed strain and velocity fields. The difference between (3) and (4) is the way to incorporate the effect of pore-water pressure: the pore-water pressures are considered as internal force, reducing the internal power dissipation, in (3), while they are considered external force in (4). By taking advantage of the normality condition , it can be easily shown that elastic stress and strain have no influence on the collapse load; that is, only plastic
Pdeformation occurs during plastic flow, and ?ij= ?ij.This makes limit analysis a simple
method to solve stability problems, without loss of rigor, assuming rigid perfect plasticity. Lower-bound Theorem
If the stress field within the soil mass is stable and statically admissible, then collapse does not occur; that is, the true collapse load is definitely greater than the applied load. This can be written in the form of the virtual work equation, using (3), as
?TSLiL'vidS??XiLvidV???ij?ijdV??(?ij??ijp)L?ijdV
VVV? ??D(?ij)dV???ij?ijdV (5)
VV31
兰州交通大学博文学院毕业设计（论文）
LWhere ?ij= statically admissible stress field in equilibrium with the traction TiL and body
force XijA not including the seepage and buoyancy force; ?ij = actual stress; ?ij = actual stain rate; and vi= velocity fields.
In (5), the inequality is due to the principle of maximum plastic dissipation, according to which the actual strain rate field is always larger than the rate of work done on the actual strain rate field by a stress field not causing collapse. In (5), only the equilibrium condition and the stress boundary conditions not taken into account. The best lower bound to the true collapse load can be found by analyzing various trial statically admissible stress fields.
中文对照翻译：
孔隙水压力作用下土坡的极限分析
摘要：极限平衡法一般用于土坡的稳定性分析。然而，众所周知的是，从极限分析法中获得的解是不严密的，因为它既不满足静态的允许条件，又不满足动态的允许条件。极限分析法充分利用了塑性体的上下边界原理，在求真实解中提供了一个相对简单但又严密的边界。在这篇文章中，三点确定的三角形三边有限元法被利用与构造在下边界分析中的静态允许应力场和上边界分析中的速度场。通过假设三角形顶点的线变量和元素变量，真实解应该是一个线形的约束问题。在静态和动态的条件都满足的基础上，真实解应该处在上下边界所得的解之间。在有限元公式中，要考虑包括了孔隙压力的影响，以便使饱和土坡的有效应力分析可以得出。作者对从不同地下水形式下简单土坡的极限分析所得的结果与极限平衡法中所得的结果作了比较。
概述：稳定性和变形问题在全球技术工程领域是一个边界值问题。微分方程必须用给定的边界条件来解决。通过解决由平衡协调条件以及沙土的本构关系推出的微分方程，从而得到边界条件下的解。按照传统的说法，在土力学中，弹性理论是用来建立变形微分方程的，就象塑性理论是用来建立稳定性问题的微分方程一样 。为了获得这个解，荷载由小到大变化，直到足够大引起部分土体的滑坡。作为土体破坏的力学行为，完整的弹塑性分析以为是一个可能的方法。然而，这样一个弹塑性分析方法很少应用于实际问题当中，因为他的计算机太过复杂。站在实践的立场上，稳定性的最基本关注点应该是土体破坏条件。因此，真实解应该是通过关注即将发生的破坏条件的一个简单的方法中得来。
自然土坡、填方土坡和挖方土坡的稳定性问题是土木工程领域碰到的最常见的问题。求解通常建立在滑移线方法上，极限平衡方法或极限分析法的基础上。由于它的简单，在实践中，极限平衡法是最被广泛使用的。极限平衡法大部分建立在分块理论的基础上，在这种理论中，假设有一个破坏的滑动面，而且在此之上的土体被划分为若干垂
32
?兰州交通大学博文学院毕业设计（论文）
直土条，整个静态平衡条件下假设的失稳表面是被确定的，一个临界的滑动破裂面必须要找到，因为它的安全因数最小。在极限平衡法的发展过程中，要努力去做的是怎么降低通过内力假设的不确定性。但是，没有一种解的得来是建立在这样的极限分析法的基础上，甚至在严格的力学意义上讲，它都不算一个严密的解。在极限平衡法中，平衡方程并不是对土体的每一点都适用的。另外，在典型的假定滑动面方法中，流动法则是不满足的，同样，协调性条件和破坏前的本构关系也是不满足的。
极限分析法充分利用了边界理论，得出了相对应用于上下边界的两个严密的解。极限分析法在以下两种意义上是严密的，一是土体在外加荷载作用下的平衡，下边界解所对应的应力场；二是与外加位移相协调，上边界所对应的速度场。就简单而言，下边界荷载作用下，滑动不会发生，但是如果下边界受到外加荷载的作用，则滑动可能立即发生。同样，在上边界作用外加荷载，滑坡也会立即发生。通过寻找下边界的最大可能解和上边界的最小可能解，真实解存在于他们之间的范围内。对于土坡稳定性问题，给定土体的性质或几何尺寸的基础上，知道土坡发生滑动的临界高度和发生部分滑坡的临界荷载才能得出解来。在过去，对于土坡稳定性的应用，大多数研究工作都集中在上边界法上，这是因为通过求解下边界适合静态允许应力场方程的解是一项很困难的任务。大多数先前的工作都基于总应力之上。对于有效应力的分析，考虑孔隙水压力的作用是很有必要的。在极限平衡法中，孔隙水压力是通过限定一个地下水表面和一个可能的流动网或者通过一个孔隙水压力比率模拟地下水条件推测出来的。相似的方法可以用于明确说明孔隙水压力作用下的极限分析。
在大量的实践中，孔隙水压力的影响被看作集中考虑在土坡稳定性问题的上边界解上。Miller和Hamilton两人实验出了两种力学破坏类型：（1）刚体旋转；（2）刚体旋转和持续变形相结合。孔隙水压力被假定成流体静力学下的一个抛物线型的自由水表面，尽管他们的研究得出了正确的答案，但是，从物理学上解释他们的研究，在能量消散上是有争议的。在他们的理论中，孔隙水压力被看作是内力，在给定的滑坡机理下，它对降低内部能量消散是有影响的。孔隙水压力也可以看作是一种外力。在Michalowsk的研究中，假设刚体是沿着螺旋线破坏的。孔隙水压力被考虑用孔隙水压力比来表示：
ru?u/?z这里，u是孔隙水压力；?是沙土的比重；z是土体表面以下的深度。它表示
孔隙水压力在内部摩擦力等于零时的分析没有影响，这就证实了用总应力分析时??0。在另一项研究中，除了用不同形状的破裂面结合分块分析法时，Michalowski秉承了相同的方法。当这种努力结合孔隙水压力作用在土坡上边界的影响中，作者没有意识到就有效应力而言有任何的下边界极限分析要做。这可能因为在考虑孔隙水压力的情况下，构造静态允许应力场的难度增大。
这篇文章的目的有两个：（1）就有效应力而言，为土坡在孔隙水压力作用下的极限分析提出了一个有限元的公式；（2）通过比较Bishop的上下边界解来核实土坡稳定性分析方法在被Bishop简化的极限平衡法所得的解与简单坡中不考虑孔隙水压力作用，上下边界所得的解相比较。在这篇文章中，在平面应变条件下，上下边界的极限分析是要考虑孔隙水压力影响的。Slon和Kleenman在考虑了上边界和下边界的情况下，利用了三点构成的三边线性代数的方法将孔隙水压力计算出来的。为了模仿应力场和速度场，由三点组成的三边线性元素就要被利用。公认的力学理论包括刚体的转动和持续变形。就分点的应力和孔隙水压力或者速度而言，用平衡方程、协调条件、流动法则、屈服准则、边界条件的线性代数等式来表达，那么，求解最佳上下边界是在几个简单土坡
33
兰州交通大学博文学院毕业设计（论文）
的构造和地下水形式下被考虑的，解是以诺莫图的形式给出。 孔隙水压力作用下的极限分析：
极限分析用了屈服准则和应力应变关系：土体在流动法则下假设成一个理想的塑性体。在这个理想的塑性体假设表明了可能的应力状态形式：
')?0 (1) F(?ij'这里，F是应变函数，?ij是有效应力张量。
通过假设应变函数F配合塑性体潜在的应变函数G，用伴随的流动法则来定义塑性体应变率，塑性应变率可以从中得出：
P?ij???G?F?? (2) ''??ij??ij这里，?是非负的塑性比值，也就是当发生塑性变形时的正直。
等式（2）通常被认为是常态条件，就表面当量而言，塑性应变率的方向关系是垂直的。理想的塑性体伴随着流动法则是一个合理的假设，如果考虑荷载下伴随非常大的位移。另外，理论研究表明，不管有没有流动法则的存在，当土体没有受到严重的受压，土坡的坍塌荷载是很不敏感的。 内部作用法则：
上边界和下边界原理都建立在内部作用原则的基础上，在发生坍塌前，给定了假设的小变形，内部作用等式是可以运用的，并且可以用以下的关系来表达：
ABABTvdS?XvdV????ii?ij?ijdVsVV'B??(?ij??ijp)A?ijdVVABii? (3)
或者
?S'BTiAviBdS??XiAviBdV???ij?ijdV (4)
VV这里，Ti是边界荷载；XijA是不包括渗透量和浮力的自重XijA是包括渗透量和浮力的自
A'A重；?是在Ti和XijA平衡下的总应力张量；?ij是在TiA和XijA平衡下的有效应力张量；
?AijA??ij是Kronecker增量；P是孔隙水压力；?ijB是与速度场viB一致的应变率张量。
不需要用任何特殊的方式通过（3）式和（4）式把?,Ti和XijA联系在一起。在（3）式和（4）式的左边代表外部作用比率，而（3）式（4）式右边代表外部能量消散比率，通过假设外加荷载的应力场和假设应变的速度场，（3）式和（4）式之间的区别在于组成有效孔隙水压力的形式不同。（3）式中，孔隙水压力被认为是内力，它减少了内部的能量消散，而在（4）式中，孔隙水压力被认为是外力。
通过充分利用常态条件，显而易见的是弹性应力和应变对坍塌荷载没有影响，即：
AijA??ij= ?ijP.也就是说，只有在塑性体流动时，塑性变形才发生。假设刚体是理想的塑性体，
34
兰州交通大学博文学院毕业设计（论文）
用一个简单的极限分析去解决稳定性问题是没有严格损失的。 下边界理论：
如果土体内的应力场是稳定的，静态允许的，那么坍塌不会发生，也就是说，真实的坍塌荷载一定比计算荷载要大。这里可以用内部作用等式的形式来表示，由（3）式得：
?TSLiL'vidS??XiLvidV???ij?ijdV??(?ij??ijp)L?ijdV
VVV? ??D(?ij)dV???ij?ijdV (5)
VVL这里?ij是不包括渗透量和浮力的自重与摩擦力平衡时的静态允许应力。?ij是真实的应
力，?ij是真实的应变率，vi是速度场。
在（5）式中，不等式的成立是由于最大塑性消散原则，根据在真实应变率场下的真实坍塌应力场得出的内部作用比率总是比在真实应变率场下而引起坍塌的应力场得出的比率要大而得出。在（5）式中， 只有平衡条件和应力边界条件是满足的，流动法则和速度边界条件是不计算在内的，对于真实的坍塌荷载，最佳的下边界是通过分析不同的 静态允许应力场试验而得出的。
35