?

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

2019-06-10 00:06:00 土木建筑与环境工程 2019年2期

赵国 陈建功

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

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

?#22411;?#20998;类号: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]基于墙后滑楔体整体极限平衡方程,建立了被动土压力计算的泛函极?#30340;?#22411;,计算结果表明,挡墙背后土体存在两种被动临界状态,两种状态的计算结果就是被动土压力大小的一个区间估计。该方法所建模型为填土面水平、墙背垂直的情况,不能用于一般情况。对于承受被动土压力的挡土墙(如桥台),有必要研究合力作用点位置对被动土压力大小的影响,本文考虑在一般情况下(倾斜墙背,填土表面非均匀堆载,墙后为黏性土体,坡面为曲面),利用变分法原理,提出基于合力作用点位置的被动土压力计算方法,由此可以得到挡墙被动土压力的大小和作用点位置的范围。

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

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

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

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

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

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

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

式中:ryD=σ(xD)cos φ,σ(xD)为滑裂面上D点的正应力,可由式(23)计算得到。通过变分法计算滑面上的正应力及滑面倾角,就可通过式(30)得到被动土压力强度的非线性分布。图9即为算例中?#34987;?#35010;面为平面时的被动土压力强度分布图。

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

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

力的影響

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

6 结论

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

参考文献:

[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] 侯键,夏?#25340;?#23380;祥冰,等. 基于土拱效应原理求解挡土墙被动土压力[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)

?
(function(){ var bp = document.createElement('script'); var curProtocol = window.location.protocol.split(':')[0]; if (curProtocol === 'https') { bp.src = 'https://zz.bdstatic.com/linksubmit/push.js'; } else { bp.src = 'http://push.zhanzhang.baidu.com/push.js'; } var s = document.getElementsByTagName("script")[0]; s.parentNode.insertBefore(bp, s); })();
热带动物园登陆
河南快三打法 时时彩三星直选技巧 辽宁十一选五定胆技巧 捕鱼大富翁安卓版 韩国花牌4人玩法 那些适合女生转赚钱的事业 怎样做云集才能赚钱 平特一肖高手专区 排列五乐彩网 江苏11选5开奖结果今天开奖号 北京十一选五遗漏数据 排列三杀号定胆 白天上班晚上怎样赚钱 飞艇计划软件app下载 广西快3走势图电脑版 官方吉林快三开奖直播