挡土墙上被动土压力的变分求解方法

赵国，陈建功

(重庆大学 土木工程学院；山地城镇建设与新技术教育部重点实验室，重庆 400045)

被动土压力的计算是挡土墙上土压力理论研究内容之一，除了对经典的朗肯土压力理论和库仑土压力理论进行扩展、补充外，学者们在考虑其他诸如土拱效应、强度理论、位移模式等方面也对被动土压力的计算进行了深入探讨。李永刚[1]、王立国等[2]考虑了土拱效应原理，按库仑理论得到土压力强度的计算方法。侯键等[3]对平移模式下刚性挡土墙和滑裂面间处于被动极限平衡状态的楔形土体进行应力分析，推导了基于土拱原理的被动土压力系数和滑裂面水平倾角,并提出被动土压力分布、土压力合力及其作用位置公式。Chang[4]、蒋波等[5]和杨泰华等[6-9]对不同变位模式下的非极限主动土压力进行计算，提出考虑位移效应的被动土压力计算理论。程康等[10]推导了平动模式下非饱和土填土挡墙的被动土压力系数及被动土压力解析解。赵均海等[11]考虑基质吸力等因素的影响，结合双剪强度理论，建立了非饱和土库仑被动土压力的统一解。

变分极限平衡法最早是由Kopácsy提出的[12-14]。李兴高、刘维宁等[15-17]基于墙后滑楔体整体极限平衡方程，建立了被动土压力计算的泛函极值模型，计算结果表明,挡墙背后土体存在两种被动临界状态，两种状态的计算结果就是被动土压力大小的一个区间估计。该方法所建模型为填土面水平、墙背垂直的情况，不能用于一般情况。对于承受被动土压力的挡土墙(如桥台)，有必要研究合力作用点位置对被动土压力大小的影响，本文考虑在一般情况下(倾斜墙背，填土表面非均匀堆载,墙后为黏性土体，坡面为曲面)，利用变分法原理，提出基于合力作用点位置的被动土压力计算方法，由此可以得到挡墙被动土压力的大小和作用点位置的范围。

1 被动土压力变分计算模型

设挡土墙高为H，坡面形状为y=g(x)，倾斜角为α，墙后土体容重为γ，墙土间的摩擦角为δ，黏聚力为c、内摩擦角为φ，q(x)为作用在坡面上的竖向压力集度。考虑极限平衡条件下墙后土体形成滑动楔体，假定其滑裂面通过挡墙墙踵，用曲线y=s(x)表示，Pp为挡土墙对滑动土体的作用合力，其值等于被动土压力合力。σ(x)、τ(x)为滑裂面上的法向和切向应力，点B为滑裂面的上端的点，其X坐标为x1，A为墙背的上端点，其X坐标为x2点，x2=-tanα，被动土压力作用点位置至墙踵的垂直距离与挡土墙高度的比值称为作用点位置系数ξ，见图1。滑楔体OAB的静力平衡方程为

图1 计算模型示意图Fig.1 Sketch map of computation

由∑X=0得

(1)

由∑Y=0得

(2)

由∑MO=0得

qx+(g-s)γx]dx=0

(3)

式中：s′=ds/dx；k=tanα，滑裂面上的应力遵循Mohr-Coulomb 破坏准则，即

τ=n1σ+c

(4)

式中：n1=tanφ。

由式(1)～式(3)，可建立相应的泛函极值模型。

由式(3)得泛函

(5)

式中：

F0=(x+ss′+n1s-n1s′x)σ-(g-s)γx-s′xc+sc-qx；

由式(1)得约束条件

(6)

式中：

F1=(-n1-s′)σ-c；

z2=z0z1=const。

由式(2)得约束条件

(7)

式中：

F2=(-n1s′+1)σ-cs′-q-γ(g-s)；

上述模型是变分学中的等周模型，可转化为无约束泛函极值模型，构造泛函J*。

(8)

F=F0+λ1F1+λ2F2

式中：λ1、λ2为拉格朗日乘子。变分模型的欧拉微分方程、边界条件及可动边界处的横截条件为

1)Euler微分方程

(9)

(10)

2)约束方程

同式(6)、式(7)。

3)边界条件

s(0)=0

(11)

s(x1)=g(x1)

(12)

4)横截条件

(13)

2 被动土压力变分计算模型的求解

由式(9)，可得到

(14)

进行坐标平移，引入新坐标

u=x+λ2，v=s-λ1

(15)

坐标平移后的模型见图2。

图2 坐标变换后的计算模型Fig.2 Computation model after coordinate

令w=v/u，则微分方程式(14)变为

(16)

分离变量得通解为

ln[u2(1+w2)]=2n1arctanw+z5

(17)

平移后的坐标换成极坐标

u=rcosθ，v=rsinθ

(18)

式(17)变为

r=z6en1θ

(19)

式中：z5、z6为任意积分常数。

由s(0)=0得原坐标原点O在新坐标中为u0=λ2，v0=-λ1，新坐标中O点的极坐标为(r0,θ0)，则滑裂面方程为

r=r0e-n1(θ0-θ)

(20)

滑裂面为对数螺旋面。

由式(10)得

-2n1σ+(-n1x+s-λ1-n1λ2)σ′-

γx-λ2γ-2c=0

(21)

式(21)在新坐标下的极坐标形式为

(22)

式(22)的通解为

式中：z7为积分常数；θ1为任意角度，可取θ1=0。

σ=z7e-2n1θ-

(23)

由式(13)可得A点处正应力为

σ(x1)=σ(θ1)=

(24)

将式(24)代入式(23)，得

z7= e2n1θ1·

上述变分问题只包含两个未知的拉格朗日常数λ1、λ2，可由式(6)和式(7)联立求出，相当于求式(25)函数Φ的零值问题。

(25)

式(25)的解可通过求解函数Φ的极小值且极小值为0得到，可采用Matlab提供的Fminsearch函数进行求解。

3 算例分析

设挡土墙墙高6 m，墙后为砂土，重度γ=18 kN/m3，内聚力c=0，内摩擦角φ=30°，挡土墙倾角α=70°，挡墙与土体的摩擦角δ=10°，坡面倾角β=20°，上部堆载q=10 kPa。

计算结果见图3，图中Φ曲线存在一Φ=0的水平直线段，此直线段的ξ系数变化范围就是被动土压力合力作用点的位置系数范围，存在上下界限值(此例中下限值ξd=0.336 9，上限值ξu=0.535 5)，在此范围内，Φ=0，也即是说，当土压力的合力作用点位于距墙锺2.02～3.21 m范围内时，墙后土体可以达到极限平衡状态。图4为被动土压力合力随作用点位置系数ξ的变化曲线，图5为滑裂面随作用点位置系数ξ的变化情况。随着作用点位置系数的增大，被动土压力呈非线性增长，在作用点位置下限处，被动土压力最小，滑裂面为对数螺旋曲面，但随着作用点位置的上移，曲率半径逐渐增大，到上限处变为平面，此时，被动土压力值达到最大。

图3 Φ函数计算结果Fig.3 Calculated result of function Φ

图4 被动土压力随系数ξ的变化曲线Fig.4 The curves of varied passive pressure according the coefficient ξ

图5 不同ξ对应的不同滑裂面曲线Fig.5 The different curves of slip face according to different values of

4 与库仑土压力理论结果对比

按变分法的计算结果与按库仑公式的计算结果见表1。可以看出，当作用点位置在上界限时，即滑裂面为平面时，被动土压力值与按库仑公式得到的土压力值是一致的，但合力作用点位置却不一致，这主要是因为库仑理论假定土压力线性分布，且不考虑土楔体力矩平衡条件。按变分法的计算结果，当滑裂面为平面时，滑裂面上的应力分布呈非线性，相应的土压力分布也应是非线性的。

表1 变分法与库仑土压力理论的计算结果比较Table 1 Result comparing of passive earth pressure using variational method and using Coulomb theory

图6所示为库仑理论的挡土墙模型，取墙后部分三角形滑动土楔体ABD作为研究对象，如图7所示。Pp(y)为墙背对滑楔体的支撑反力合力，也即被动土压力合力；R(y)为不动土体对滑楔体的反力合力；G(y)为土体ABD的重量；Q(y)为作用在AD面上竖向压力。

图6 库仑刚性挡土墙模型图Fig.6 COULOMB’s Model of retaining

图7 土楔体ABD受力分析Fig.7 mechanical analysis on sliding wedge

土楔体ABD的力平衡条件见图8，由三角形正弦定理可得

(26)

令pp(y)为墙背对滑楔体的支撑反力强度，也即被动土压力强度，r(y)为不动土体对滑楔体的反力强度，即

(27)

将式(27)代式(26)，得

(28)

考虑几何关系

(29)

式(28)两边对y求导，得

(30)

图8 土楔体ABD的力平衡条件Fig.8 Force equilibrium of sliding wedge

图9 土压力强度分布图Fig.9 Distribution of earth

表1中的安全系数一栏是指按库仑公式进行挡墙的抗倾覆设计，取1.5的安全系数时所对应的按本文方法计算的抗倾覆安全系数范围。若按库仑理论进行设计，实际的被动土压力合力作用点位置比库仑理论设定的要高，其安全系数达不到1.5的要求，所以，应优化挡墙的结构，控制挡墙的变位模式，使合力作用点降低，提高抗倾覆安全系数。

5 坡面起伏和非均匀超载对被动土压力的影响

设挡墙高6 m，墙后土重度γ=18 kN/m3，内聚力c=10 kPa，内摩擦角φ=20°，挡土墙倾角α=70°， 坡面形状函数g(x)和非均匀超载q(x)均用三角函数模拟坡面的起伏和非均匀超载。表2为计算结果，可以看出，被动土压力合力和作用点受坡面的起伏以及坡面超载的不均匀性影响不能忽视。

表2 考虑坡面起伏及非均匀超载的计算结果Table 2 Computed results considering surface undulation of earth filling and uneven surface load

6 结论

基于墙后滑楔体整体极限平衡方程，利用变分法原理推导了被动土压力泛函极值的变分模型，提出了相应的计算方法，当作用点位置已知时，可有效确定被动土压力大小和滑裂面形状位置。算例分析表明，墙后土体的极限平衡状态对应于作用点位置系数范围，随着作用点位置系数的增大，被动土压力呈非线性增长，在作用点位置下限处，被动土压力最小，滑裂面为对数螺旋曲面，但随着作用点位置的上移，滑裂面曲率半径逐渐增大，到上限处变为平面。此时，被动土压力值达到最大，且与按库仑理论计算的被动土压力一致，但作用点位置位于墙背中点以上，与库仑理论所假定的合力作用点位置不一致。