衡重式路肩挡土墙计算：土压力与稳定验算的Python实现

简介：一份面向道路桥梁工程专业学生和一线设计人员的衡重式路肩挡土墙完整计算书，系统覆盖4米墙高（不出现第二破裂面）、4米墙高（出现第二破裂面）、7米墙高、9米墙高四种典型工况，针对低墙和高墙分别处理第二破裂面、土压力分布及墙体稳定等问题。从基本资料、土压力计算、墙身截面设计，到稳定性验算、地基承载力计算、墙身强度和稳定性验算，形成完整设计链条，每一步均给出公式、参数取值与算例结果。压缩包内仅含1个doc格式文件，大小2.77MB，内容附带中英文摘要、前言、四类工况详细计算过程、结论和附录，便于直接打印或对照学习。已有335人学习下载。计算书依据《公路挡土墙设计与施工技术细则》编写，融合库仑与朗金土压力理论，并重点说明第二破裂面出现前后对设计和安全系数的影响，对课程设计、毕业设计乃至实际衡重式路肩墙选型计算具有较高参考价值。

1. 衡重式路肩挡土墙的计算书到底要算穿什么

拿到一份《公路工程衡重式路肩挡土墙设计(计算书).doc》，先不要急着翻土压力公式。这个名字已经说明两件事：一是墙型是衡重式，二是交付物是计算书而不是单纯图纸。衡重式的特点在于墙背中段伸出一个平台（衡重台），用台上填土的自重参与抵抗倾覆和滑移，因此比同等高度的普通重力式挡墙更省圬工，山区公路路肩段用得很多。但代价是计算路径变长：上墙、下墙土压力要分开算，稳定验算的竖向力多了一项台背土重，基底合力位置的判断也比直线墙背更敏感。这篇文章就顺着计算书的推进顺序，把每一段的取数原则、验算边界和容易返工的点说清楚，适合正在做路基详设或复核别人计算书的工程师参考。

2. 土压力拆成上墙和下墙：衡重式挡土墙计算的第一步

2.1 为什么库仑平面破裂面假设在这里失效

普通重力式挡土墙的墙背是直线，库仑理论的场景是：一个三角形的土楔，沿墙背和一个从墙踵发出的破裂面滑动。衡重式挡墙的墙背是一条折线，上墙段放得比较缓，中间又横着伸出一个衡重台，墙踵位置也往往比想象中靠后。此时最危险的滑动楔体不再贴着上墙背走，而是从衡重台的外端向上发展出一个新的破裂面，土力学里叫第二破裂面。墙背与第二破裂面之间的那块土，实际上已经“坐”在衡重台上，参与了稳定而不是单纯产生压力。

这就是衡重式区别于普通挡墙的核心：土压力计算时要把墙背分成两段。上墙按第二破裂面理论计算，平台以下的直线段按力多边形法处理；如果一道计算书里只有一个库仑系数算出个总的主动土压力，多半是把问题简化过头了，审图时容易被卡。

2.2 上墙第二破裂面与库仑系数的快速验证

现行《公路路基设计规范》对衡重式挡土墙的土压力计算有专门附录，思路是先假定第二破裂面的角度，用试算法找到最不利情况。手算时一般从 30° 到 50° 每隔 2° 试一次；程序化之后，这就是一个一维枚举问题。在拿到完整计算书之前，可以用库仑系数先做个量级验证：给定填土内摩擦角 φ、墙背与竖直面的夹角 α、墙背摩擦角 δ，库仑主动土压力系数可以直接这样算：

import math def coulomb_ka(phi, alpha, delta, beta=0.0): # phi 填土内摩擦角, alpha 墙背与竖直面夹角(外倾为正) # delta 墙背摩擦角, beta 填土面坡角(水平为0) pr, ar, dr, br = map(math.radians, (phi, alpha, delta, beta)) s = math.sqrt( math.sin(pr + dr) * math.sin(pr - br) / (math.cos(ar + dr) * math.cos(br - ar)) ) return math.cos(pr - ar) ** 2 / ( math.cos(ar) ** 2 * math.cos(ar + dr) * (1 + s) ** 2 ) for alpha in range(0, 26, 5): print(f"alpha={alpha:>2}° Ka={coulomb_ka(35, alpha, 20):.4f}")

这段代码的用处不是替代规范附录，而是快速感受参数敏感度：同一组 φ=35°、δ=20° 下，墙背从竖直渐变到外倾 25°，Ka 大约从 0.25 升到 0.44，土压力直接放大七成。换句话说，第二破裂面法的上墙土压力通常低于用坦墙背直接套库仑公式的结果，而很多计算书的“偏安全”恰恰是这里偏过头，导致断面浪费。墙背摩擦角 δ 的取值也很关键：浆砌片石与填土之间一般取 (1/3～2/3)φ，混凝土墙面略低；这个值越大，土压力的竖直分量越大，对抗滑是有利的，但抗倾覆力矩里该分量的力臂不能想当然。

2.3 下墙算法与衡重台上土重的归并

下墙一般按力多边形法计算，它和直接套库仑公式的差别在于：上墙传下来的一部分土压力会改变下墙背后土楔的平衡。工程上也有简化做法，即把下墙背视为一个新的墙背单独算，但结果通常比力多边形法小 10% 左右，偏不安全。我做复核时，如果计算书用的是简化库仑，会重点看它的稳定系数是不是贴着限值，贴着限值的必须按力多边形法重算。

衡重台上的土重是衡重式最划算的一项“免费”竖向力。这部分土直接压在平台上，墙越高、台越宽，它贡献的抗滑力与抗倾覆力矩越可观。计算截面时注意：台面上是一个宽度略大于台宽的梯形土体，而不是矩形；墙背有倾角时，上墙背后的三角区域也有一部分土要归到这一项里。常见的错误是只按台宽乘以平台以上墙高取一个矩形面积，这会少算 10%～20% 的竖向力。上墙、下墙主要计算参数的参考范围整理如下：

3. 衡重式路肩挡土墙稳定性验算的参数边界：抗滑、抗倾覆与基底应力

3.1 抗滑稳定：Kc 的构成与基底摩擦系数 f

衡重式路肩挡土墙的抗滑验算公式是所有稳定验算里最直观的：Kc = (ΣN·f + Ep)/ΣT。ΣN 是作用于墙底的竖向力之和，包括墙身自重、衡重台上土重、上墙与下墙土压力的竖向分力；ΣT 是水平土压力之和；Ep 是墙前被动土压力，计算书里一般不计或只计一半。需要明确，规范规定的主力组合下 Kc 不应小于 1.3，这是硬指标，不能拿 f 取大值来凑。

f 的取值与基底和地基土的摩擦特性直接相关，也是最容易被审图追问的参数。常见做法是按基底所处土层的性质取：硬塑黏土 0.30～0.35，密实砂土与碎石土 0.40～0.50，软质岩 0.40～0.50，硬质岩 0.60～0.70。如果地基是黏性土且浸水，f 还要乘折减系数。抗滑不足时的处理顺序有讲究：先看能不能把基底做成向内倾斜（倾斜角一般不大于 10°），再看能不能加宽踵板让衡重台土重更多压到基底上，最后才是加深基础或用齿坎增加被动抗力。直接提高 f 去通过验算，计算书拿到现场就是风险。

3.2 抗倾覆验算的力臂容易错在三个地方

K0 = ΣMy / ΣM0，力矩参考点取墙趾（基底最前端）。抗倾覆力矩来源是墙身自重、衡重台土重以及土压力竖向分力，倾覆力矩主要来自上墙和下墙的水平土压力。力臂计算有三个高频错误：

第一个，衡重台前上方那块土漏算。上墙墙背是斜的，墙背与第二破裂面之间的土体形状不是矩形，很多计算书只按衡重台宽度乘台面以上填土高度取矩形面积。正确做法是把这个区域按梯形分解，重心位置也要随之调整，不然竖向力 N 和它的力臂同时出错。

第二个，上墙土压力的作用点被误放到总墙高的 1/3 处。上墙的土压力分布按三角形，作用点应在上墙段内，即距上墙顶 H1/3，折算到墙趾的竖直力臂是 H−H1/3。直接套 H/3 会把倾覆力矩高估不少，稳定系数看起来偏保守，但断面的经济性就很差了。

第三个，墙身自重合并成一个整体重心。衡重式墙身截面由上下两段不同坡度的梯形组成，有的还带基础台阶，合成一个梯形求重心会与实际相差很大。程序验算时我会把截面拆成多边形分别求面积和形心，手算时至少分上墙梯形、下墙梯形和基础底板三块。

3.3 基底合力偏心距与压应力控制

墙底受力不只是看平均压应力，还要看合力作用点的位置。先算出对墙趾的合力矩 M=ΣMy−ΣM0，得到合力的作用点距墙趾 c=M/ΣN，偏心距 e=B/2−c。控制条件有两个：一是 e 不能大于 B/6，否则基底出现拉应力区，黏性土地基上通常不允许；二是最大压应力 σmax=ΣN/B·(1+6e/B) 不得超过地基承载力基本容许值 [σ0]。

这里有个容易被忽略的点：作用点 c 不是越大越好。c 太大意味着合力靠近墙踵，虽然 e 满足要求但墙踵侧压应力会超；c 太小表明合力靠近墙趾，倾覆安全系数虽然高但基底前缘可能压溃。调整断面时两者要同时看。下面这张表是三个验算项目的速查：

4. 把衡重式路肩挡土墙验算写成 Python 脚本

4.1 结构设计：几何参数与岩土参数分开建模

写验算脚本的第一步不是写公式，而是把输入分好类。断面几何参数属于结构自身，岩土参数属于外部条件，分开建模的好处是出一个参数化断面后可以批量试算不同地质组合。我用 dataclass 组织这两组参数，默认值给一组典型的 8 米高路肩墙，实际使用时直接改字段即可：

from dataclasses import dataclass import math @dataclass class WallGeometry: H: float = 8.0 # 总墙高 m H1: float = 3.2 # 上墙高 m B: float = 5.0 # 基底宽 m b_top: float = 0.6 # 墙顶宽 m gamma_c: float = 24.0 # 墙身重度 kN/m3 @dataclass class SoilParams: gamma: float = 19.0 # 填土重度 kN/m3 phi: float = 35.0 # 填土内摩擦角 deg delta: float = 23.0 # 墙背摩擦角 deg f: float = 0.42 # 基底摩擦系数 sigma0: float = 280.0 # 地基承载力 kPa

这里把 H 和 H1 分开是因为衡重式的上墙段和下墙段要分别算土压力，后续所有计算函数都从 wall 和 soil 两个实例取数。alpha1 这类墙背倾角没有放进 dataclass，是为了让土压力函数显式接收它，试算时可以直接变化角度而不改断面对象。

4.2 核心计算函数：从库仑系数到三组稳定系数

接下来实现三个函数：coulomb_ka 计算主动土压力系数，earth_pressure 把系数变成水平、竖向分力，check_stability 汇总验算。墙体自重按梯形截面近似，重心距前趾取 0.42B，衡重台上土重作为竖直力单独计入：

def coulomb_ka(phi, alpha, delta, beta=0.0): pr, ar, dr, br = map(math.radians, (phi, alpha, delta, beta)) s = math.sqrt( math.sin(pr + dr) * math.sin(pr - br) / (math.cos(ar + dr) * math.cos(br - ar)) ) return (math.cos(pr - ar) ** 2 / (math.cos(ar) ** 2 * math.cos(ar + dr) * (1 + s) ** 2)) def earth_pressure(gamma, H, phi, alpha, delta): ka = coulomb_ka(phi, alpha, delta) E = 0.5 * gamma * H * H * ka return E * math.cos(math.radians(delta)), E * math.sin(math.radians(delta)) def check_stability(wall, soil, alpha1=20.0, alpha2=10.0): h2 = wall.H - wall.H1 e1x, e1y = earth_pressure(soil.gamma, wall.H1, soil.phi, alpha1, soil.delta) e2x, e2y = earth_pressure(soil.gamma, h2, soil.phi, alpha2, soil.delta) aw = (wall.b_top + wall.B) * wall.H / 2.0 # 墙体梯形面积 ww = aw * wall.gamma_c # 墙身自重 wt = soil.gamma * wall.H1 * (0.25 * wall.H) # 衡重台土重(近似) mw = ww * 0.42 * wall.B mt = wt * 0.78 * wall.B me1 = e1x * (wall.H - wall.H1 / 3.0) me2 = e2x * (h2 / 3.0) me1y = e1y * 0.55 * wall.B me2y = e2y * 0.70 * wall.B n = ww + wt + e1y + e2y t = e1x + e2x mr = mw + mt + me1y + me2y mo = me1 + me2 kc = soil.f * n / t k0 = mr / mo c = (mr - mo) / n eb = wall.B / 2.0 - c pmax = n / wall.B * (1 + 6 * eb / wall.B) return kc, k0, eb, pmax, n

这段代码里的重心位置系数 0.42B、0.78B 等是简化近似，正式计算书应该用截面多边形积分，但脚本足够用来做方案比选和敏感性分析。注意 e1x、e2x 计算时用了各自墙段的计算高度：上墙是 H1，下墙是 H2=H−H1，不能混用。衡重台土重的 0.25H 是取台宽的常见比例，对应的重心力臂取 0.78B，是把台面以上梯形土体的中心和墙踵方向偏移考虑进去的结果。

4.3 用判断格式输出，结果直接能贴进计算书

验算函数返回四组数之后，最好直接打印成判断语句而不是只给数值。这样在方案比选时能一眼看出哪个断面卡在哪个指标上：

wall = WallGeometry() soil = SoilParams() kc, k0, eb, pmax, n = check_stability(wall, soil) print(f"抗滑稳定 Kc = {kc:.2f} {'通过' if kc >= 1.3 else '不通过'}") print(f"抗倾覆 K0 = {k0:.2f} {'通过' if k0 >= 1.5 else '不通过'}") print(f"偏心距 e = {eb:.2f} m (B/6={wall.B/6:.2f} m) " f"{'通过' if abs(eb) <= wall.B/6 else '不通过'}") print(f"基底最大压应力 = {pmax:.0f} kPa " f"({'通过' if pmax <= soil.sigma0 else '不通过'})")

输出形如：Kc=1.84，K0=2.12，e=0.58，Pmax=236kPa，这组数据对应的是 8 米高、底宽 5 米的典型断面。实际调整时，如果 Kc 不够，优先把 f 换回现场实测值复核而不是直接改大；如果 e 超限，说明截面形心太靠前，可以考虑把下墙背放缓或加宽踵板。脚本里所有近似参数（重心系数、台宽比例）都标注了来源，避免三个月后回来看不懂当时怎么取的。

5. 衡重式路肩挡土墙计算书收尾时的四个检查点

5.1 构造要求不是“满足规范”四个字

计算书最后通常会写一排构造说明，但审图会逐条看数据。泄水孔间距 2～3 米，最下一排出水口要高出地面不小于 0.3 米；墙背回填必须设反滤层，否则泄水孔堵住后墙背水压力会直接毁掉稳定验算的前提。伸缩缝和沉降缝按 10～15 米设置，缝宽 2～3 厘米，缝内沿墙的内、外、顶三边填塞沥青麻筋。基础埋深要超过冻结线且不小于 1 米，墙趾前襟边宽度也要和地基承载力匹配。

5.2 复核时优先看三个量

拿别人做的衡重式路肩挡土墙计算书时，我一般先看三个输出：基底合力偏心距 e 是否压着 B/6 的线，抗滑稳定系数是否刚好卡 1.3，以及上墙土压力是按第二破裂面还是按简单库仑算的。前两个体现设计余量，第三个直接反映计算路径是否正确。e 如果恰好等于 B/6，说明断面尺寸是精心凑出来的，需要怀疑是否调整过重心位置来通过验算。

5.3 让参数来源变成“可引用”的批注

真正好用的计算书，每一个关键参数旁边都有来源：填土重度来自地勘报告哪个孔，基底摩擦系数来自哪本规范或现场试验，地基承载力是深宽修正前还是修正后的值。建议在计算书正文里放一张参数来源表，把 γ、φ、δ、f 四项的取值依据全部列出来。这样审图和施工复核都不用回头翻地勘报告，也大幅减少“这个 f 为什么取这么大”这类重复质询。

5.4 用脚本批量复核一个工点的全部挡墙

一套工点的衡重式路肩挡土墙可能几十段，每段墙高和地质条件都不同。把上面的 WallGeometry 和 SoilParams 改成从 CSV 读取，循环调用 check_stability，一次就能输出所有墙段的验算汇总表，把 e、Kc、Pmax 超限的段落标红。给计算书配这样一个验证脚本，比在 Word 里反复改公式要可靠得多。

本文还有配套的精品资源，点击获取