含柔性垫层的刚性挡土墙土压力计算方法

梁利生 郭俊源 王慧芳

摘要：在剛性挡土墙与填土之间设置柔性垫层能减小作用于挡土墙的土压力，但目前仍缺乏针对设置聚苯乙烯土工泡沫（EPS）柔性垫层的刚性挡土墙土压力计算方法。将EPS柔性垫层的压缩量视为墙后填土的位移量，考虑挡土墙后土拱效应，基于挡土墙土压力-位移的关系曲线，引入迭代法进行收敛计算，得到设置EPS柔性垫层的刚性挡土墙土压力计算方法。该计算方法的优势是可在EPS柔性垫层压缩量未知的情况下求解土压力，即可应用工程设计阶段。建立FLAC3D有限差分数值模型，对推导的理论解进行验证，并对EPS柔性垫层减载效果进行分析。结果表明：基于土压力-位移关系曲线并采用迭代法得出的墙后设置EPS柔性垫层的刚性挡土墙土压力理论解具有较好的合理性。在EPS柔性垫层弹性模量不变的情况下，EPS柔性垫层减小土压力的效果随着EPS柔性垫层厚度的增加而增强；在EPS柔性垫层厚度一定的情况下，随着EPS柔性垫层弹性模量的增加，其减小土压力的效果逐渐减弱。

关键词：挡土墙；EPS柔性垫层；土压力；减载；土拱效应

中图分类号：TU476.4     文献标志码：A     文章编号：2096-6717（2023）06-0158-07

Calculation for lateral earth pressure on rigid retaining wall with flexible cushion

LIANG Lisheng1， GUO Junyuan2， WANG Huifang1

（1. Department of Civil Engineering; The Cultivation Base of Shanxi Key Laboratory of Mining Area Ecological Restoration and Solid Wastes Utilization， Shanxi Institute of Technology， Yangquan 045000， Shanxi， P. R. China；2.College of Civil Engineering， Taiyuan University of Technology， Taiyuan 030024， P. R. China）

Abstract： The installation of polystyrene geofoam （EPS） inclusions between the rigid retaining wall and backfill can reduce the earth pressure acting on the retaining wall. To date， there is a lack of the calculation method for the earth pressure on rigid retaining walls with EPS inclusions. In this study， the compression from geofoam inclusions is equated to the displacement of the backfill， and the soil arching behind the retaining wall is considered. Based on the displacement-dependent earth-pressure model （DDEP）， the iterative method is introduced for the convergence calculation， and then the theoretical solution is obtained to calculate the earth pressure on the rigid retaining wall with EPS inclusion. The advantage of the proposed method is that it can be used to solve the earth pressure when the EPS compression  modulus is not provided. A series of FLAC3D numerical models is established to verify the proposed method and analyze the load reduction effect of EPS inclusions. The proposed theoretical solutions based on the DDEP model and iterative method are in preferably consistency with the numerical results. Effect of EPS on reducing the earth pressure increases with increase of the thickness of EPS inclusions when EPS elastic modulus is constant， and decreases with increase of the elastic modulus of EPS inclusions while EPS thickness is constant.

Keywords： retaining wall； EPS flexible inclusion； earth pressure； load reduction； soil arching

挡土墙广泛应用于土木工程领域中，作用于挡土墙的侧向土压力计算是其设计的关键指标。聚苯乙烯土工泡沫（EPS）是一种具有密度小、压缩变形大、强度高和减振性能好等优点的柔性材料。将EPS垫层设置于刚性挡土墙和填土之间后，EPS垫层受到的土压力将会使其产生一定的压缩量，该部分压缩量即转化为填土的位移。在刚性挡土墙自身位移量很小的情况下，设置EPS垫层后，挡土墙受到的土压力将会明显减小[1]。

已有诸多学者对设置EPS的刚性挡土墙开展了研究。Bathurst等[2-3]开展了室内振动台试验研究不同密度的EPS缓冲层减小刚性挡土墙在动荷载作用下的墙背土压力的效果。Ikizler等[4]开展针对墙后填土为膨胀土的EPS垫层减载试验，得出EPS垫层可减小水平压力和竖向压力。Ertugrul等[5-6]开展了缩尺模型试验，得出EPS垫层的厚度和刚度对于土压力的减小有主要影响。Ni等[7]开展缩尺模型试验，分析EPS垫层的厚度与刚度对于减压效果的影响。Kim等[8]基于模型试验建立FLAC数值模型，分析EPS垫层的材料类型和厚度对土压力的影响，并提出刚性挡土墙后EPS垫层的分段设计方案。郑俊杰等[9-11]開展一系列模型试验，对EPS减载挡土墙土压力的效果及挡土墙后EPS垫层的设计方案进行了研究。谢明星等[12]建立设置EPS垫层的刚性挡土墙数值模型，分别对挡土墙后填土的静止、主动及被动3种位移状态进行了研究。许晓亮等[13]通过开展不同级别荷载下的模型试验，分析墙后铺设EPS柔性垫层、回填泡沫轻量土及其同时施加时的墙后土压力大小及分布特征。

目前，仍缺乏设置EPS垫层后刚性挡土墙土压力的计算方法。笔者将EPS垫层的压缩量视为墙后填土的位移量，基于Mei等[14]提出的挡土墙土压力-位移的关系曲线，推导设置EPS垫层的刚性挡土墙土压力的理论解。首先，考虑挡土墙后土拱效应，推导出考虑土拱效应的主动和被动土压力系数；其次，将所推导的主动和被动土压力系数代入土压力-位移关系曲线，并将静止土压力作为初始状态，引入迭代法进行收敛计算，求解EPS的压缩量和对应的土压力；最后，建立FLAC3D有限差分数值模型，验证所得出理论解的合理性，并对EPS垫层减载效果进行分析。

1 挡土墙土压力-位移的关系曲线

Mei等[14]将挡墙背离土体方向的位移记为负值、朝向土体方向的位移记为正值，假定存在一个函数p（s）可表示土压力p和位移s之间的关系（图1），则该函数p（s）满足以下条件：

1）p（s）为单调递增的函数。

2）p（s）函数存在上下限，下限为主动土压力pa及其对应的位移sa，上限为被动土压力pp及其对应的位移sp。

3）当s=0时，p为静止土压力p0。

进而基于上述条件提出土压力-位移的关系曲线

其中

式中：k0、ka、kp分别为静止、主动、被动土压力系数。

2 考虑土拱效应的土压力系数

挡土墙无位移的情况下，目前最常用的土压力系数计算方法为Jaky[15]针对粗粒土得出的半经验半理论的静止土压力系数公式。

k0 = 1 - sin φ （4）

式中：φ为墙后填土的内摩擦角。

在墙后填土达到主动或被动状态过程中，由于挡土墙墙背并非绝对光滑，挡土墙与填土之间存在的摩擦作用必将引起土体应力偏转，形成土拱效应，从而对土压力的分布产生影响。为便于计算，研究基于以下假定：

1）挡土墙为刚性挡土墙，且挡土墙墙背垂直；

2）挡土墙与基岩之间的填土为无黏性土，符合Mohr-Coulomb破坏准则，且不考虑孔隙水作用；

3）土压力求解视为平面应变条件下的二维问题；

4）对于刚性挡土墙，其自身的位移量通常非常小。因此，忽略挡土墙的位移量，将EPS垫层的压缩量视为墙后填土的位移；

5）考虑土拱效应时，主动状态下填土的小主应力轨迹线和被动状态下的大主应力轨迹线均假定为圆弧形。

朗肯土压力理论中假定主动和被动状态下填土内滑动面的倾角θa、θp分别为π/4+φ/2、π/4-φ/2。图2为主动和被动状态下填土的小主应力轨迹线。为便于计算，Paik等[16]将小主应力轨迹假定为圆弧形。笔者将小主应力轨迹和大主应力轨迹均假定为圆弧形。图3为主动和被动状态下墙背上和填土内滑动面上任意一点土体的应力莫尔圆。

可求出主动状态下两处土体对应的主应力偏转角αaw和αas为

被动状态下的主应力偏转角αpw和αps为

主动状态下微单元体在墙背处的侧向应力paw和微单元中任意点的垂直向应力paz表示为

式中：αa为主动状态下该点的主应力偏转角。

被动状态下微单元体在墙背处侧向应力ppw和微单元任意点的垂直向应力ppz可表示为

式中：αp为该点在被动状态下的主应力偏转角。

根据图2的几何关系可求得小主应力轨迹线半径ra和大主应力轨迹线半径rp。

式中：Ba和Bp分别为主动和被动状态下挡土墙与填土内滑动面之间土体微分单元的宽度。将作用于微单元体的垂直方向的合力除以微单元体的宽度，

微单元体在挡土墙上的侧向应力除以微单元体的垂直应力平均值即为侧向应力比[17]，则侧向主动应力比Ka和侧向被动应力比Kp为

根据作用在水平微单元片上的应力平衡条件可得出

其中

C = tan θa （tan θa + tan δw ） （21）

D = tan θp （tan θp - tan δw ） （22）

不考虑填土表面荷载的情况下，边界条件为paz│z=0=ppz│z=0=0。则可求出

挡土墙上任意一点主动土压力系数ka和被动土压力系数kp为

3 迭代法的引入

在挡土墙与填土之间设置EPS垫层后，EPS的压缩是一个逐渐变化的过程。在这一过程中挡土墙的土压力p随着EPS压缩量s的增大而减小，直至p和s的值稳定。假定EPS为理想线弹性，将这一过程分为n步、每一步分为两个独立求解的部分：1）EPS受到一个已知的土压力发生弹性变形产生一个压缩量，在这一部分的求解过程中该已知的土压力保持不变；2）通过式（1）可求得在该压缩量下的土压力，在这一部分的求解过程中该压缩量保持不变，得出的土压力将作为下一步中已知的土压力。然后，引入迭代法，计算每一步中EPS的压缩量和作用于挡土墙的土压力，具体计算步骤如下。

初始状态：EPS未发生变形，即s0=0，挡土墙受到的土压力为静止土压力p0。

p0 = k0 γz （27）

第1步：EPS在p0的作用下产生压缩量s1，进而求得在s1的位移下挡土墙所受到的土压力p（s1）。

第n步：EPS在p（sn-1）的作用下产生压缩量sn，进而求得在sn的位移下挡土墙所受到的土压力p（sn）。

当p（sn）约等于p（sn-1）时，停止迭代，sn和p（sn）分别为EPS的最终压缩量和设置EPS后挡土墙的最终土压力。需要注意的是，上述计算过程中，需用土与EPS之间的界面摩擦角δE代替土与挡土墙之间的界面摩擦角δw。

4 计算结果验证及分析

目前，EPS减载挡土墙土压力在实际工程中的应用较少，缺乏现场试验数据。既有研究大多为室内缩尺试验研究，受限于挡土墙的高度，所得到的土压力远小于实际工程中挡土墙所受到的土压力。而EPS的减载效果随着作用于挡土墙的土压力的减小而减小，这将导致室内缩尺试验不能如实反映EPS的减载效果。因此，笔者采用有限差分数值软件FLAC3D建立与实际工程中相符合的挡土墙模型（图3），对所提出计算方法的合理性进行验证。以EPS厚度和弹性模量为变量建立多个数值模型。各模型中EPS的弹性模量参考文献[18]取值，且EPS垫层横断面尺寸均与挡土墙一致。此外，为验证数值模拟的准确性，建立墙后未设置EPS板的工况，以便于模拟结果与Jaky静止土压力进行对比。数值模拟方案见表1。

填土为砂土，采用Mohr-Coulomb模型，黏聚力c =0 kPa，内摩擦角φ=32°。挡土墙为刚性挡土墙，墙高H为5 m，完全限制其侧向位移。地基采用线弹性模型，并取较大的弹性模量值，以减弱地基土沉降对填土应力状态的影响。填土与挡土墙之间（墙后未设置EPS板的工况）以及填土与EPS板之间（墙后设置EPS板的工况）均设置接触面，接触面摩擦角分别为φ/2和φ/3[12]。数值模型中各材料的参数取值见表2。

图4为未设置EPS时（基本工况）的静止和主动模拟值与Jaky静止土压力和库仑主动土压力理论解的对比。静止土压力的理论解与数值解能够吻合；主动土压力的理论解与数值解在0～4.5 m范围内也能吻合，但在挡土墙底部位置的数值解明显小于理论解。这是由于，为减小地基土沉降对土压力的影响，数值模型中地基土弹性模量的取值远大于墙后填土。当墙后填土产生位移的过程中，填土与地基土之间存在一定的边界摩擦效应，从而会导致土压力模拟值的減小[8]。综上，不考虑边界效应的情况下，数值模型具有较好的准确性。

图5对比了设置EPS后各工况的土压力数值解与理论解（式（29）迭代计算的结果）。可以看出，在距墙顶高度0～4 m范围内模拟值略大于理论解。土压力迭代计算中任意点的EPS压缩量均按照离散点受压的弹性变形计算，尽管EPS材料的泊松比较小，但EPS垫层受压方向面积较大，各点的弹性变形并非都发生在一维应力状态下，可能导致EPS垫层压缩量的模拟值略小于理论值，从而使土压力的模拟值略大于理论解。同样受数值模型边界条件的影响，在挡土墙底部位置各工况的模拟值均明显小于理论解。因此，忽略边界条件的影响，提出的理论解总体上与数值解能较好地吻合，验证了提出方法的合理性。

图6、图7对比了设置EPS后各工况的土压力减小效果。随着距墙顶高度的增加，各工况的土压力减小效果均越来越明显。对比工况1～工况3，即在EPS弹性模量不变的情况下（E=2.34 MPa），EPS减小土压力的效果随着EPS的厚度的增加而增强；对比工况3～工况6，在EPS厚度一定的情况下（t=0.5 m），EPS减小土压力的效果随着EPS的弹性模量的增加而减弱。EPS的弹性模量随着其密度的减小而减小，但当采用低密度的EPS时，其屈服强度也较低，若挡土墙土压力值较大，会导致EPS处于屈服或非线性硬化阶段。同时，EPS厚度若过大，将会导致紧邻挡土墙区域路面出现不均匀沉降。因此，在实际工程中，需综合考虑EPS的设计参数和减载效果，避免为追求减载效果而导致工程隐患。

5 结论

1）将EPS垫层的压缩量视为墙后填土的位移量，推导了考虑土拱效应的挡土墙主动和被动土压力系数解析解，将其代入挡土墙土压力-位移的关系方程，并引入迭代法进行收敛计算，解决了设计阶段EPS实际压缩量未知这一关键问题，最终得出设置EPS垫层的挡土墙土压力解析解。

2）建立一系列有限差分数值模型，得到的模拟值与理论解总体上能较好地吻合，验证了提出方法的合理性。

3）在EPS弹性模量不变的情况下，EPS减小土压力的效果随着EPS的厚度的增加而增强；而在EPS厚度一定的情况下，随着EPS的弹性模量的增加，EPS减小土压力的效果逐渐减弱。

参考文献

[1]  HORVATH J S. The compressible inclusion function of EPS geofoam [J]. Geotextiles and Geomembranes， 1997， 15（1/2/3）： 77-120.

[2]  BATHURST R J， ZARNANI S， GASKIN A. Shaking table testing of geofoam seismic buffers [J]. Soil Dynamics and Earthquake Engineering， 2007， 27（4）： 324-332.

[3]  BATHURST R J， ZARNANI S. Earthquake load attenuation using EPS geofoam buffers in rigid wall applications [J]. Indian Geotechnical Journal， 2013， 43（4）： 283-291.

[4]  IKIZLER S B， AYTEKIN M， NAS E. Laboratory study of expanded polystyrene （EPS） geofoam used with expansive soils [J]. Geotextiles and Geomembranes， 2008， 26（2）： 189-195.

[5]  ERTUGRUL O L， TRANDAFIR A C. Reduction of lateral earth forces acting on rigid nonyielding retaining walls by EPS geofoam inclusions [J]. Journal of Materials in Civil Engineering， 2011， 23（12）： 1711-1718.

[6]  ERTUGRUL O L， TRANDAFIR A C. Lateral earth pressures on flexible cantilever retaining walls with deformable geofoam inclusions [J]. Engineering Geology， 2013， 158： 23-33.

[7]  NI P P， MEI G X， ZHAO Y L. Displacement-dependent earth pressures on rigid retaining walls with compressible geofoam inclusions： Physical modeling and analytical solutions [J]. International Journal of Geomechanics， 2017， 17（6）： 04016132.

[8]  KIM H， WITTHOEFT A F， KIM D. Numerical study of earth pressure reduction on rigid walls using EPS geofoam inclusions [J]. Geosynthetics International， 2018， 25（2）： 180-199.

[9]  郑俊杰， 景丹， 谢明星， 等. 刚性挡土墙后设置EPS垫层的设计方法研究[J]. 华中科技大学学报（自然科学版）， 2020， 48（12）： 96-101.

ZHENG J J， JING D， XIE M X， et al. Study on design scheme of installing EPS cushion behind the rigid retaining wall[J]. Journal of Huazhong University of Science and Technology （Natural Science Edition）， 2020， 48（12）： 96-101. （in Chinese）

[10]  鄭俊杰， 邵安迪， 谢明星， 等. 挡土墙后EPS板减压性能模型试验研究[J]. 土木与环境工程学报（中英文）， 2020， 42（4）： 1-9.

ZHENG J J， SHAO A D， XIE M X， et al. Model test study on decompression performance of EPS boards behind retaining walls [J]. Journal of Civil and Environmental Engineering， 2020， 42（4）： 1-9. （in Chinese）

[11]  郑俊杰， 邵安迪， 谢明星， 等. 不同填土宽度下设置EPS垫层挡土墙试验研究[J]. 岩土力学， 2021， 42（2）： 324-332.

ZHENG J J， SHAO A D， XIE M X， et al. Experimental study on retaining wall with EPS cushion under different backfill widths [J]. Rock and Soil Mechanics， 2021， 42（2）： 324-332. （in Chinese）

[12]  谢明星， 郑俊杰， 邵安迪， 等. 设置EPS柔性垫层的刚性挡土墙土压力研究[J]. 岩土工程学报， 2019， 41（Sup1）： 65-68.

XIE M X， ZHENG J J， SHAO A D， et al. Lateral earth pressure on rigid retaining walls with EPS deformable geofoam inclusions [J]. Chinese Journal of Geotechnical Engineering， 2019， 41（Sup1）： 65-68. （in Chinese）

[13]  許晓亮， 王鑫佩， 曾林风， 等. 含EPS垫层及轻量填土挡墙的墙后土压力及其减载效应[J/OL]. 土木与环境工程学报（中英文）. http：//kns.cnki.net/kcms/detail/50.1218.TU.20221202.1110.001.html.

XU X L， WANG X P， ZENG L F， et al. Lateral earth pressure and load-reducing effect behind retaining wall with EPS cushion and light-weight soil [J/OL]. Journal of Civil and Environmental Engineering. http：//kns.cnki.net/kcms/detail/50.1218.TU.20221202.1110.001.html.（in Chinese）

[14]  MEI G X， CHEN Q M， SONG L H. Model for predicting displacement-dependent lateral earth pressure[J]. Canadian Geotechnical Journal， 2009， 46（8）： 969-975.

[15]  JAKY J. The coefficient of earth pressure at rest [J]. Journal of the Society of Hungarian Architects and Engineers， 1944： 355-358.

[16]  PAIK K H， SALGADO R. Estimation of active earth pressure against rigid retaining walls considering arching effects [J]. Géotechnique， 2003， 53（7）： 643-653.

[17]  HANDY R L. The arch in soil arching [J]. Journal of Geotechnical Engineering， 1985， 111（3）： 302-318.

[18]  XIE M X， ZHENG J J， SHAO A D， et al. Study of lateral earth pressures on nonyielding retaining walls with deformable geofoam inclusions [J]. Geotextiles and Geomembranes， 2020， 48（5）： 684-690.