土的压缩性和地基沉降计算

官方地理老师·29 页·深入(追求细节与边界)·0 次浏览·2 天前
压缩性沉降计算分层总和法固结理论

土的压缩性与地基沉降计算

从压缩机理到分层总和法,掌握地基沉降的预测方法与适用边界

按 空格/→ 演示下一步

1 / 29 页

全部页面点击任意一页,跳回舞台从这页播放

压缩性沉降计算分层总和法固结理论

土的压缩性与地基沉降计算

从压缩机理到分层总和法,掌握地基沉降的预测方法与适用边界

1第 1 页 · 土的压缩性与地基沉降计算

为什么要研究土的压缩性

从一栋楼歪了说起——比萨斜塔每年还在以肉眼看不见的速度倾斜。这不是施工失误,而是土被压实的结果。

比萨斜塔
建于1173年,因地基土压缩不均至今仍在缓慢倾斜
差异沉降危害
同一建筑不同部位下沉量不同,导致墙体开裂、管道断裂
相邻建筑影响
新建建筑加载或基坑开挖引起邻近房屋沉降与开裂
长期工后沉降
高速公路、地铁运营数年后仍持续下沉,影响行车安全
经济代价沉重
事后加固与纠偏费用往往是沉降计算投入的数十倍
湿海绵上压重物对应 →土的压缩与沉降

水被挤出体积变小;土中孔隙水排出土骨架压密下沉

s=σihiEsi——分层总和法核心公式s = \sum \frac{\sigma_i h_i}{E_{si}} ——分层总和法核心公式
2第 2 页 · 为什么要研究土的压缩性

室内压缩试验的原理

侧限压缩仪里,一块土饼如何被逐级压实、排出孔隙水——五步看清全过程。

1
土样装入环刀
原状土切成圆饼压入刚性金属环刀,建立侧限:只许上下动,不许横向胀
2
透水石就位
上下放透水石和加压盖板,让孔隙水有出路,荷载能均匀传到土面
3
分级加压排水
荷载分级施加(常按 50→100→200→400 kPa 翻倍),每级让孔隙水经透水石排出
4
骨架压缩变形
土骨架承担有效应力,颗粒滑移重排,土样变薄,孔隙比减小
5
记录并绘曲线
读各级稳定变形量,算孔隙比,绘制 e-p 或 e-log p 压缩曲线
3第 3 页 · 室内压缩试验的原理

压缩仪的构造与试验流程

沿「取样→加荷→读数」主线,看环刀、加压系统、变形量测三大部件如何接棒完成一次压缩试验。

图解渲染中…
a2环刀:定体积并提供侧限边界a4加压系统:按等级表逐级施荷a5变形量测:百分表记录沉降a6稳定判据:每小时变形≤0.01 mm
4第 4 页 · 压缩仪的构造与试验流程

试验成果:e-p曲线的形成

上一节我们把土样装进压缩仪、逐级加荷、读变形。这一页把所有数据画出来——横轴压力p、纵轴孔隙比e,每级荷载对应一个点,串成一条曲线。

坐标系建立
横轴压力p、纵轴孔隙比e,每级荷载对应一个点
数据点来源
每级加荷稳定后读取变形,代入公式换算成孔隙比
典型形态
曲线向下凹,前期陡、后期缓,不是一条直线
物理含义
压力越大孔隙比越小,且单位增量引起的孔隙比变化越来越小
分段信息
初次压缩、回弹、再压缩三条曲线对应不同加载历史
海绵挤水的过程对应 →e-p曲线形态

一压就出很多水(孔隙比骤降),越压越难再出水(曲线趋缓)

ei=e01+e0H0ΔHie_i = e_0 - \frac{1+e_0}{H_0}\sum \Delta H_i
5第 5 页 · 试验成果:e-p曲线的形成

什么是一维压缩

上一页我们看到了 e-p 曲线,它能成立靠一个关键假设:试验中土样被钢环紧紧箍住,只能竖向变形。这个'只能竖向'就是今天要拆开看的一维压缩。

侧限条件
土体水平方向被刚性边界约束,εx = εy = 0
一维竖向压缩
只在竖直方向发生变形与排水固结
K₀ 静止土压力
侧压力系数 K₀ = σ'h/σ'v,即水平与竖向有效应力之比
适用场景
宽阔基础中心点下、均质深厚土层的沉降估算
局限性
靠近边坡、基坑、挡土墙时侧向变形不可忽略
硬罐里的一摞弹簧对应 →侧限一维压缩

罐壁挡住四周,弹簧只能沿轴向被压短——对应土样被钢环箍住

$\varepsilon_x = \varepsilon_y = 0, \quad K_0 = \sigma'_h / \sigma'_v$
6第 6 页 · 什么是一维压缩

压缩曲线的分类与判读

从同一种土样出发,看单调加荷与中途卸荷再加载两条路径,各得到什么曲线、工程意义有何不同。

图解渲染中…
D分支点:是否做卸荷再加载试验G门槛:历史最大压力Pc,超固结判据Hp<Pc段:斜率为Cs,沉降小Ip>Pc段:斜率变Cc,沉降剧增
7第 7 页 · 压缩曲线的分类与判读

压缩系数 vs 压缩指数

都是衡量压缩性的指标,但来源曲线、是否为常数、适用场景都不同。

压缩系数 a
  • 来自 e-p 普通坐标曲线
  • 单位压力引起的孔隙比变化
  • 随压力区间变化,不是常数
  • 规范分层总和法推荐使用
压缩指数 Cc
  • 来自 e-lg p 半对数曲线
  • 孔隙比对数压力的变化速率
  • 高压段斜率近似为常数
  • 重大工程或高压沉降估算
常规沉降用 a 走规范法;高压力或重要工程用 Cc 更精确。
8第 8 页 · 压缩系数 vs 压缩指数

压缩系数的定义与应用

上一节 e-p 曲线刻画了孔隙比随压力变化的全过程,但整条曲线不便直接用于工程判别。这一节问:能否抽出一个'代表数字',像弹簧柔软度那样刻画土的压缩灵敏?这个数字就是压缩系数 a。

压缩系数定义
a = Δe/Δp,即 e-p 曲线某压力段的割线斜率,单位 MPa⁻¹
物理含义
单位压力增量引起的孔隙比减量——土对压力的灵敏程度
代表值 a₁₋₂
取 p=100~200 kPa 段斜率,对应一般浅基应力增量范围
压缩性分级
a₁₋₂<0.1 低、0.1~0.5 中、≥0.5 高压缩性土
为何取代表段
e-p 非线性,全段斜率不均;取代表段才能跨场地横向比较
弹簧的软硬对应 →土的压缩系数 a

硬弹簧压力增一点形变小(低 a=低压缩性),软弹簧同压力形变大(高 a=高压缩性)

a=e1e2p2p1=ΔeΔpa = \frac{e_1 - e_2}{p_2 - p_1} = -\frac{\Delta e}{\Delta p}
9第 9 页 · 压缩系数的定义与应用

压缩指数的计算方法

上一节我们看到,e-p 曲线随压力增大越来越平缓,所以压缩系数 a 不是定值。如果把横坐标换成 log p,曲线在常规荷载段就近似为一条直线——这条直线的斜率,就是压缩指数 Cc。

计算公式
Cc = -(e1-e2)/(log p2 - log p1),直线段上两点求斜率
取点位置
必须选在直线段内,常取 p=100、200 kPa 两点
几何意义
e-log p 曲线直线段的斜率,表征压缩性强弱
与 a 的换算
Cc = a·p/ln10 ≈ 0.434·a·p,Cc 不随压力变化
一段直坡路的坡度对应 →Cc 作为直线斜率

直路坡度是个固定值,对应 e-log p 直线段的稳定斜率 Cc

Cc=e1e2log10p2log10p1C_c = -\dfrac{e_1 - e_2}{\log_{10} p_2 - \log_{10} p_1}
10第 10 页 · 压缩指数的计算方法

先期固结压力的确定

上节我们从 e-log p 曲线算出了压缩指数,但这条曲线还藏着更关键的信息——土体历史上的最大压力。Casagrande(1936)用一个几何作图法把它「读」了出来,这就是确定先期固结压力的经典方法。

找最弯点
在 e-log p 曲线上找出曲率半径最小(弯得最厉害)的那个点
作两线
过该点画一条水平线,再画一条切线,二者构成夹角
作角平分线
对这个夹角作角平分线
延长直线段
将曲线尾部的直线段(正常压缩段)向后延长
交点即 Pc
角平分线与延长线的交点对应横坐标,就是先期固结压力
照片中下蹲的膝盖对应 →e-log p 曲线上的拐点

站立到下蹲,膝盖处曲率最大;曲线在 Pc 处也最弯,这正是作图法的几何根基

OCR=Pcσ0OCR = \dfrac{P_c}{\sigma'_0}
11第 11 页 · 先期固结压力的确定

压缩模量的概念

压缩系数 a 描述的是孔隙比随压力的变化率,但工程师算沉降更习惯用'应力÷应变'这种模量形式。土里对应的就是压缩模量 Es。

压缩模量 Es
侧限条件下,应力增量 Δp 与竖向应变 Δεz 的比值
应变与孔隙比
一维压缩、土粒不变时,εz = -Δe/(1+e0)
公式推导
把 a = -Δe/Δp 与 Es = Δp/Δε 联立得到
局限性
Es 随压力段变化,并非材料常数;仅适用侧限
应用要点
计算沉降须取实际应力范围对应的 Es
弹簧的劲度系数 k对应 →压缩模量 Es

都反映'压得难不难';k 是常数,Es 随压力段走

εz=Δe1+e0,a=ΔeΔp,Es=ΔpΔεz=1+e0a\varepsilon_z = -\frac{\Delta e}{1+e_0},\quad a = -\frac{\Delta e}{\Delta p},\quad E_s = \frac{\Delta p}{\Delta\varepsilon_z}=\frac{1+e_0}{a}
12第 12 页 · 压缩模量的概念

沉降的三个阶段

压缩模量衡量的是土被压密的能力,但工程里更关心地基加载后沉降怎么随时间展开——它分三个阶段依次发生。

瞬时沉降 Si
加载瞬间发生,形状畸变与侧向挤出,饱和土体积不变,按弹性理论估算
主固结沉降 Sc
超静孔压消散,骨架逐步承荷,体积减小,是沉降最主要的部分
次固结沉降 Ss
孔压消散后骨架蠕变与结构调整,持续时间长,软土与有机土显著
总沉降叠加
三阶段机理不同,需分别计算后叠加得到总沉降
挤压饱水海绵对应 →三阶段沉降

瞬时=形状先变形水未排出;主固结=水慢慢渗出;次固结=水排完纤维仍在缓慢调整

S=Si+Sc+SsS = S_i + S_c + S_s
13第 13 页 · 沉降的三个阶段

分层总和法的原理

从地基与荷载出发,逐层取中心点应力,把每层压缩量累加即得总沉降。

图解渲染中…
D1σc:上覆土自重产生的有效应力E1σz:基底压力在地基中引起的附加应力F1Es:压缩模量,由e-p曲线对应段求取H1累加只在压缩层下限zn以内进行
14第 14 页 · 分层总和法的原理

应力计算:角点法

用角点法求矩形荷载下任意点的附加应力,需四步完成叠加。

1
划分矩形
把荷载面积拆成以计算点为公共角点的若干子矩形
2
求影响系数
按布氏解算出每个子矩形的角点应力影响系数αc
3
叠加求和
点内子矩形取正号、外侧补出的小矩形取负号再求代数和
4
得到σz
把叠加结果乘以荷载q,得到该深度处的附加应力
15第 15 页 · 应力计算:角点法

应力分布:应力影响线

地基中σz等值线呈「压力泡」形态,从荷载中心向深处和两侧扩散衰减。

图解渲染中…
a4σz相等的点连成的曲线,类似地形等高线a5等值线围成的灯泡形区域,工程上称「压力泡」
16第 16 页 · 应力分布:应力影响线

分层沉降量的计算步骤

分层总和法把地基切成薄层,逐层算压缩再求和,共六步。

1
地基分层
按 h_i≤0.4b 切分,保证层内应力近似均匀
2
计算自重应力
求各分层中点的 σ_cz,作为初始应力基准
3
计算附加应力
用角点法或影响线求中点处附加应力 σ_z
4
确定受压层下限
以 σ_z=0.2σ_cz 处为底,更深即忽略
5
计算分层压缩量
由 e-p 曲线得 e1、e2,套 s_i 公式
6
累加得总沉降
s=Σs_i,得到地基最终沉降量
17第 17 页 · 分层沉降量的计算步骤

应力面积法 vs 规范法

两种都用分层求和,但曲线形式与参数选取不同,选错会显著低估沉降。

应力面积法
  • 基于 e-log p 曲线,半对数坐标
  • 关键参数:压缩指数 Cc
  • 公式:S = Σ CcH/(1+e₀)·lg(p₁/p₀)
  • 直接适用超固结土,可分段处理 OCR
规范法
  • 基于 e-p 曲线,普通坐标
  • 关键参数:压缩模量 Es(或 a)
  • 公式:S = ψs·Σ(Δp/Es)·H
  • 主要适用于正常固结土
超固结土或重点工程优先用应力面积法;常规工程按规范取 ψs 修正即可。
18第 18 页 · 应力面积法 vs 规范法

沉降计算的工程实例

python

用分层总和法计算方形基础沉降:角点法求附加应力,按 si=aiσh/(1+e0) 累加各层沉降得到总沉降量。

代码高亮加载中…

L20 用四点角点叠加把基础中心应力放大 4 倍;L24 是分层总和法核心公式;L28 把各层沉降累加即工程报告里的总沉降量。

19第 19 页 · 沉降计算的工程实例

渗流固结的物理过程

饱和黏土受荷瞬间,水来不及排出,外力全由孔隙水承担——超静孔压由此而生。排水固结,就是它随时间消散的过程。

超静孔压 u₀
荷载瞬时施加时,孔隙水承担了全部新增应力
有效应力增长
水从边界排出,σ' 逐步替代 u 承担外荷
总应力守恒
σ = σ' + u 不变,内部两者此消彼长
时间依赖
固结快慢取决于 k、mv 及排水路径长度
湿海绵垫重物对应 →饱和土渗流固结

挤出的水≈消散的孔压;海绵纤维逐渐接住的力≈增长的σ'

σ=σu\sigma' = \sigma - u
20第 20 页 · 渗流固结的物理过程

太沙基固结模型

上下两栏对照,左列机械模型,右列土体原型;同位节点用虚线相连,标注力学对应关系。

图解渲染中…
a3小孔↔渗透性k,决定排水快慢a4弹簧↔土骨架刚度,决定沉降量b2瞬时u=P,骨架尚未受力b4有效应力原理σ'=σ-u 在此体现
21第 21 页 · 太沙基固结模型

一维固结微分方程

上一页我们搭好了太沙基物理模型:弹簧、水、土骨架。现在要把它翻译成数学语言——推导思路其实就三块拼图。

连续性方程
流出水量 = 土骨架压缩量,水量守恒
达西定律
v = k·i,把流速和水力坡降挂上钩
有效应力原理
σ = σ' + u,把孔压 u 引进应力方程
联立消元
三个方程嵌套,把 v 和 σ' 消掉,只剩 u
热在金属棒中扩散对应 →孔压在土层中消散

数学形式完全相同——都是扩散方程,分离变量、傅里叶级数的解法也通用

ut=Cv2uz2\frac{\partial u}{\partial t} = C_v \frac{\partial^2 u}{\partial z^2}
22第 22 页 · 一维固结微分方程

固结系数的确定

时间平方根法与时间对数法的应用

固结系数的确定
时间平方根法与时间对数法的应用
23第 23 页 · 固结系数的确定

固结方程的求解

解一维固结方程分五步:从已知条件出发,最后拿到固结度公式。

1
初始条件
t=0时超静孔压u₀沿深度的分布
2
边界条件
顶面/底面排水或不透水的设定
3
分离变量
把时间与深度拆为两个函数相乘
4
傅氏级数
把任意u₀分布展为正弦级数之和
5
解析解
导出固结度U关于时间因数Tv的公式
24第 24 页 · 固结方程的求解

时间因子Tv与固结度U的关系

顺查:已知t求U;反查:已知U求t。两条路最终都汇到沉降计算。

图解渲染中…
a1解一维固结微分方程得到U(Tv)解析式a2Tv→∞时U→100%,曲线呈S形渐近a4U<60%用Tv≈π/4×(U/100)²反算a5U>60%用对数近似Tv=-0.933log(1-U/100)-0.085
25第 25 页 · 时间因子Tv与固结度U的关系

固结度的工程应用

任意时刻沉降量与超静孔隙水压力的推算

固结度的工程应用
任意时刻沉降量与超静孔隙水压力的推算
26第 26 页 · 固结度的工程应用

核心知识要点回顾

  • a、Cc、Es 都源自同一条 e-p 曲线,数学上可互算。
  • 最终沉降=各分层应力增量×厚度÷压缩模量之和。
  • 沉降随时间=超静孔压消散驱动有效应力增长。
  • 主线贯穿应力、变形、时间三轴,分立又耦合。
  • 小试样指标到工程应用,须经规范经验修正。
延伸主题:次固结(蠕变)沉降三维应力路径压缩性非饱和土固结理论
27第 27 页 · 核心知识要点回顾

自测检验

点击作答

太沙基一维固结理论中,若排水距离 Hd 加倍而其他条件不变,达到同一固结度所需的时间 t 如何变化?

28第 28 页 · 自测检验

课后思考

先自己想一想,再翻参考答案——三个问题分别对应核心概念、工程应用、理论边界。

1为何在正常固结范围内,e-p曲线越压越平,而e-log p曲线却近似一条直线?

参考答案提示——想想孔隙比与土骨架体积的关系,以及对数坐标如何"展开"高压段的微小变化。直线斜率即压缩指数Cc。

2一栋建筑荷载分两次施加,中间停工半年,沉降曲线与一次性加载相比会有什么不同?

参考答案停工期间已完成部分固结,二次加载时沉降起点不再归零,而是从已有固结度接续——这正是应力历史与时间因子Tv的工程意义。

3太沙基一维固结方程假设土体饱和、渗流符合达西定律,非饱和土还适用吗?为什么?

参考答案不再适用。方程依赖水不可压缩与达西定律;非饱和土含气、渗流偏离达西,需要更一般的非饱和固结理论。

29第 29 页 · 课后思考