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

赵国 陈建功

摘 要：基于滑楔体整体极限平衡方程，根据变分法原理推导了被动土压力泛函极值的变分模型，并引入拉格朗日乘子，将等周变分模型转化为含有两个函数自变量的泛函极值模型。依据欧拉方程、边界条件和横截条件，得到了滑裂面函数和滑裂面上的应力函数，函数泛函极值模型转化为两个未知量的函数优化模型。算例表明，对于一般土体，在作用点位置系数下界限处，滑裂面呈现对数螺旋曲面，此时被动土压力最小;当作用点位置上移时，被动土压力呈非线性增长，在作用点位置系数上界限处，滑裂面为平面，被动土压力达到最大，与库仑土压力理论解完全一致，但作用点在墙体的相对位置并非在墙高的1/3处。结果表明，被动土压力大小和作用点位置受坡面的起伏和坡面超载的不均匀性影响比较明显。

关键词：刚性挡墙;被动土压力;变分极限平衡法;滑裂面

中图分类号：TU476

文献标志码：A  文章编号：2096-6717（2019）02-0029-07

Abstract：Considering an inclined rough retaining wall under the general conditions such as curvilinear fill， cohesive soil and uneven surface load， the functional extreme-value isoperimetric model about passive earth pressure is deduced based on the force equilibrium equations of the sliding mass. Then， the model can be transferred into a functional extreme-value problem with two undetermined functions by introducing Lagrange undetermined multiplier. According to Euler equations， Logarithmic spiral slip surface and normal stress distribution along the slip surface are obtained. Combined with the boundary conditions and transversality conditions， the conditional functional extremum problem of passive earth pressure  involves searching the minimum of unconstrained optimizations of function with two unknown Lagrange multiplier. Results show that the passive earth pressure resultant force is minimal when the point of resultant force is on the lower bound and it increase nonlinearly as the point of resultant force moved up to the upper bound for general soil. Accordingly， the slip face evolves from logarithmic spiral face to plan. Although the magnitude of passive earth pressure reaches maximal that is the same with result calculated from Coulombs theory， the application point of earth pressure is not at 1/3 height of the retaining wall. In addition， curvilinear fill and uneven surface load have significant effect on both the magnitude and  location of application point of passive earth pressure.

Keywords：rigid retaining wall; passive earth pressure; variational limit equilibrium method; slip surface

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

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

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

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

设挡土墙墙高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为滑裂面随作用点位置系数ξ的变化情况。随着作用点位置系数的增大，被动土压力呈非线性增长，在作用点位置下限处，被动土压力最小，滑裂面为对数螺旋曲面，但随着作用点位置的上移，曲率半径逐渐增大，到上限处变为平面，此时，被动土压力值达到最大。

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

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

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

式中：ryD=σ（xD）cos φ，σ（xD）为滑裂面上D点的正应力，可由式（23）计算得到。通过变分法计算滑面上的正应力及滑面倾角，就可通过式（30）得到被动土压力强度的非线性分布。图9即为算例中当滑裂面为平面时的被动土压力强度分布图。

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

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

力的影響

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

6 结论

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

参考文献：

[1] 李永刚. 挡土墙被动土压力研究[J]. 岩土力学，2003，24（2）： 273-276.

LI Y G. Study on passive earth pressure of retaining wall [J]. Rock and Soil Mechanics， 2003， 24（2）： 273-276.（in Chinese）

[2] 王立国， 王志凯. 考虑土拱效应的挡土墙被动土压力分析[J]. 水利与建筑工程学报， 2010， 8（4）： 146-149.

WANG L G， WANG Z K. Analysis on passive earth pressure for retaining walls considering soil arching [J]. Journal of Water Resources and Architectural Engineering， 2010， 8（4）： 146-149. （in Chinese）

[3] 侯键，夏唐代，孔祥冰，等. 基于土拱效应原理求解挡土墙被动土压力[J]. 岩土力学， 2012， 33（10）： 2296-3000.

HOU J ，XIA T D，KONG X B， et al. Passive earth pressure on retaining walls calculated by principle of soil arching effect [J]. Rock and Soil Mechanics， 2012， 33（10）： 2296-3000. （in Chinese）