一维固结理论：Terzaghi 解析解与固结度

饱和软土层受到瞬时均布荷载时，沉降不会立刻完成：荷载最初几乎全部由孔隙水承担，形成超孔隙水压力；随后孔隙水从排水边界缓慢排出，超孔压消散，荷载逐步转移给土骨架，地层随之压密。这个过程就是固结。它同时是一个水力扩散问题和一个骨架变形问题——只盯着孔压或只盯着沉降，都看不到全貌。

本文给出 Terzaghi 一维固结理论的控制方程与解析解，说明时间因子和排水路径如何决定消散速度，最后用计算实验室的一维固结实验核对孔压剖面与固结度。阅读本文需要有效应力原理（$\sigma = \sigma' + u$）和达西定律的基础。

一、基本假设与符号

Terzaghi 一维固结理论建立在以下假设上：

• 土层只发生竖向变形，孔隙水只沿竖向渗流；
• 土体完全饱和，土颗粒与孔隙水不可压缩；
• 渗流服从达西定律，渗透系数 $k$在固结过程中保持常数；
• 土骨架线弹性，体积压缩系数 $m_v$保持常数；
• 荷载瞬时施加，此后保持不变。

符号约定：$z$自土层顶面（排水面）向下为正，土层厚度为 $H$；$u(z,t)$为超孔隙水压力；$\Delta\sigma$为瞬时施加的均布荷载增量。低渗透土中加载瞬间几乎无水排出，因此初始条件为

$$u(z,0) = \Delta\sigma, \qquad 0 \le z \le H.$$

排水条件决定边界：排水面处 $u=0$；不透水面处流量为零，即 $\partial u/\partial z = 0$。定义 $H_{dr}$为最远排水距离：单面排水时 $H_{dr} = H$，双面排水时 $H_{dr} = H/2$。

二、固结微分方程

取深度 $z$处厚度为 $dz$的单元体。孔隙水流速服从达西定律，超孔压梯度驱动的水力梯度为 $-\dfrac{1}{\gamma_w}\dfrac{\partial u}{\partial z}$，其中 $\gamma_w$为水的重度。水量连续要求：单位时间流出与流入的水量差，等于单元体积的压缩速率。体积压缩由有效应力增量控制，而总应力不变时 $\Delta\sigma' = -\Delta u$，即孔压消散多少，有效应力就增长多少。合并三式得到

$$\frac{\partial u}{\partial t} = c_v \frac{\partial^2 u}{\partial z^2}, \qquad c_v = \frac{k}{m_v \gamma_w} \tag{1}\label{eq:consolidation-pde}$$

其中 $c_v$为固结系数，单位 $\mathrm{m^2/s}$。式 \eqref{eq:consolidation-pde} 与热传导方程同形：超孔压以"扩散"的方式消散，快慢由 $c_v$与排水距离共同决定，而不是只取决于渗透性——土越软（$m_v$越大），同样渗透性下固结反而越慢。

三、解析解与时间因子

以单面排水（顶面排水、底面不透水）为例，式 \eqref{eq:consolidation-pde} 在上述初边值条件下的级数解为

$$\frac{u(z,t)}{\Delta\sigma} = \sum_{m=0}^{\infty} \frac{2}{M} \sin\left(\frac{M z}{H_{dr}}\right) \exp\left(-M^2 T_v\right), \qquad M = \frac{\pi}{2}(2m+1) \tag{2}\label{eq:pore-pressure-series}$$

其中无量纲时间因子

$$T_v = \frac{c_v t}{H_{dr}^2}. \tag{3}\label{eq:time-factor}$$

$T_v$把 $c_v$、真实时间 $t$与排水距离 $H_{dr}$打包成一个数：不同厚度、不同土性的土层，只要 $T_v$相同，孔压剖面形状就相同。双面排水时剖面关于土层中面镜像对称，对每一半使用同一个解、取 $H_{dr} = H/2$即可。

工程上更常用的是平均固结度 $U$，即已消散的平均超孔压占 $\Delta\sigma$的比例，也等于已完成沉降与最终沉降之比 $s_t / s_\infty$：

$$U = 1 - \sum_{m=0}^{\infty} \frac{2}{M^2} \exp\left(-M^2 T_v\right) = 1 - \sum_{m=0}^{\infty} \frac{8}{\pi^2 (2m+1)^2} \exp\left(-\frac{(2m+1)^2 \pi^2 T_v}{4}\right). \tag{4}\label{eq:average-consolidation}$$

两个便于记忆和手算核对的锚点：$T_v = 0.2$时 $U \approx 50\%$；$T_v = 0.848$时 $U = 90\%$。实用近似方面，$U < 60\%$时可用 $U \approx \sqrt{4T_v/\pi}$，$U$较大时级数取第一项（$m=0$）已足够准。

四、排水路径决定时间尺度

式 \eqref{eq:time-factor} 中 $H_{dr}$以平方出现，这是固结问题最要紧的一条 scaling：双面排水把排水距离减半，达到同一固结度所需的时间变为单面排水的 $1/4$。砂井、塑料排水板等地基处理手段的本质，就是人为缩短 $H_{dr}$。

用实验室实验的示例参数换算一遍：层厚 $H = 10\ \mathrm{m}$、$c_v = 1\times10^{-7}\ \mathrm{m^2/s}$。达到 $U \approx 50\%$（$T_v = 0.2$）：

• 双面排水：$H_{dr} = 5\ \mathrm{m}$，$t = T_v H_{dr}^2 / c_v = 0.2 \times 25 / 10^{-7} \approx 5\times10^7\ \mathrm{s}$，约 1.6 年；
• 单面排水：$H_{dr} = 10\ \mathrm{m}$，$t \approx 2\times10^8\ \mathrm{s}$，约 6.3 年，正好是 4 倍。

五、用交互实验核对

打开一维固结实验，可以这样核对本文的结论：

• 拖动时间因子 $T_v$从 0 增大：孔压剖面从初始的矩形（$u/\Delta\sigma = 1$）开始，总是从排水边界一侧先消散，不透水边界附近最后；
• 把 $T_v$停在 0.2：页面给出的平均固结度应约为 50%，荷载分担条上孔隙水与土骨架各承担约一半；
• 切换单面/双面排水：同一 $T_v$下固结度不变（$U$只由 $T_v$决定），但剖面对称性不同，"等效经过时间"变为 4 倍或 1/4，与第四节的换算一致。

六、模型边界

使用解析解时不要越过它的假设：

• $c_v$恒定只在小应变、应力增量较小时成立；软土大变形固结中 $k$与 $m_v$随有效应力显著变化，需要 Gibson 大应变固结理论或数值求解；
• 分级施工不是瞬时加载，沉降曲线需按叠加原理修正；
• 砂井地基存在井阻与涂抹效应，实际固结速率低于理想径向固结解。

即便如此，解析解仍是检验数值方法的第一道基准：任何两相耦合的 FEM 或 MPM 程序，都应先在一维固结上同时对上孔压剖面、平均固结度与沉降三条曲线——只对上其中一条，往往意味着耦合系数或应力基准处理有问题，只是尚未暴露。固结分析的起点是与自重平衡的初始应力场，其建立方法见《K₀ 初始应力场：自重应力与静止土压力系数》。

本文含有 AI 生成内容