岩土为何越压越难压？ 读懂 Cam-Clay 模型
PETROLEUM PAPER · 高级地质力学
L32 / Cam-Clay 模型（上）
岩土为何越压越难压？
读懂 Cam-Clay 模型
临界状态线、屈服面、各向同性固结线：用一组计算，读懂材料对加载历史的“记忆”。
把一团松散、含水的沉积物压实，再卸去一部分压力，它通常不会恢复原来的蓬松程度。接着对它施加剪切，同样的当前压力下，经历过强烈压密的试样与刚刚压到该压力的试样，还可能表现出不同的体积变化和强度演化。
这说明，判断岩土响应，只有“现在受到多大应力”还不够。我们还需要知道孔隙有多大、过去压到了什么程度，以及接下来的加载会让材料压密还是剪胀。Cam-Clay 把这些问题放进了同一个力学框架。
本讲范围 / 先把模型说准确
本讲采用具有椭圆屈服面的修正剑桥模型（Modified Cam-Clay，简称 MCC）。原始 Cam-Clay 的屈服面形式不同，不能混用。这里以饱和、近似各向同性的可压缩岩土骨架为对象，压应力与压缩体积应变均取正。
01 / 状态变量先用三个量，描述材料的处境
MCC 最常用的状态空间由平均有效应力 p′、偏应力强度 q 和比体积 v 构成。前两个量描述受力，第三个量描述材料有多疏松。
p′ = (σ′1 + σ′2 + σ′3) / 3
q = √(3J2)
σ′ 为有效应力；J₂ = ½sᵢⱼsᵢⱼ，s 为偏有效应力张量。p′、q 都具有应力单位，例如 MPa。p′ 不是孔隙压力，本讲将孔隙压力另记为 u。
对常规轴对称三轴压缩，σ′2 = σ′3，且 σ′1 ≥ σ′3，于是 q = σ′1 − σ′3。三个有效主应力都相等时，q = 0，但 p′ 可以很大：材料没有偏应力，也可能发生不可逆压密。
e = Vp / Vs
v = V / Vs = 1 + e
n = Vp / V = e / (1 + e)
Vₚ、Vₛ、V 分别为孔隙、固体和总体积；e 是孔隙比，v 是比体积，n 是孔隙度，三者均无量纲。本讲 v 不表示速度。
例如，孔隙比 e = 0.5 表示孔隙体积是固体体积的一半，对应孔隙度 n = 1/3，并不是 50%。把孔隙比和孔隙度混用，会连带算错压缩曲线与状态参数。
02 / 临界状态线不断剪切，却不再改变应力和体积
临界状态的核心是：材料可以继续产生剪切变形，而平均有效应力、偏应力和比体积保持稳定。应变仍在累积，稳定的是应力与体积状态。它与“刚刚开始屈服”不是同一个概念。
图 01｜蓝色虚线为临界状态线，蓝色椭圆为某一硬化状态下的屈服边界。二者有特定交点，但不能把整条虚线当作当前屈服面。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.16 局部。
q = M p′
M 为临界状态应力比，无量纲。CSL 是 Critical State Line，即临界状态线。只有满足相应的体积状态并进入持续塑性剪切，才构成完整的临界状态。
对于按三轴压缩条件匹配的摩擦型材料，可以用临界状态摩擦角 φcs 表示 M：
M = 6 sin φcs / (3 − sin φcs)
这是一条三轴压缩匹配关系，不能不加说明地用于所有应力路径。当 φcs = 30° 时，M = 1.2。
图 02｜把三轴压缩的主应力比换算为 q/p′，就得到临界状态斜率 M。板书中的 φcs = 30° 对应 M = 1.2。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.18 原图。
容易忽略的一点
某个应力点恰好满足 q/p′ = M，并不意味着它已经达到临界状态：它可能仍位于当前屈服面内部。不排水试验中总体积近似不变，也不能单独证明临界状态已经建立。
03 / 屈服椭圆同一条临界状态线，可以对应不同大小的屈服面
MCC 用预固结压力 p′c 控制当前屈服面的大小。它记录材料的压密历史，是会随塑性变形演化的内部变量；公开板书采用 p′0 记号，本讲统一写为 p′c。
f = q² + M²p′(p′ − p′c) = 0
f < 0 表示当前屈服面内部；f = 0 为边界。这里 f 的量纲为应力的平方。标准模型讨论正的平均有效压应力，不把靠近原点的形式外推为可靠的抗拉模型。
图 03｜令 q = 0 可确定椭圆在平均应力轴上的截点；再联立 q = Mp′，得到临界状态与当前屈服面的非零交点。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.19 原图。
三个位置，读懂一只椭圆
“顶部为什么在一半的位置”可以直接算出来。由屈服式得到 q² = M²p′(p′c − p′)，右侧关于 p′ 是一条开口向下的抛物线，最大值出现在 p′ = p′c/2。再代入，就得到 q = Mp′c/2。
椭圆两侧，为什么会出现压密与剪胀？
标准 MCC 采用关联流动法则：塑性应变增量的方向由屈服函数的梯度决定。对平均应力求偏导，就能判断塑性体积变化的方向：
dεvp = dΛ · M²(2p′ − p′c)
dΛ ≥ 0 为塑性乘子增量。本讲压缩为正：椭圆右半部产生塑性压缩，左半部产生塑性剪胀，顶部的塑性体积应变增量为零。此判断以正在发生塑性加载为前提。
因此，在顶部可以持续产生塑性剪切，却不再靠塑性体积压缩来扩张屈服面。临界状态、零塑性体积变化与椭圆顶部，正是通过流动法则联系起来的。单有一个椭圆图形，还不足以完整确定应变。
04 / 压密与硬化材料的“记忆”，存在哪里？
图 04｜塑性压密使预固结压力增加，屈服椭圆逐步扩张。图中的应力路径与应变曲线共同表示硬化过程，椭圆不是一成不变的边界。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.20 原图。
把材料压密后，它在再次压缩时往往需要更高的应力才继续产生塑性体积变化。在 MCC 中，这种响应通过 p′c 增大来表达。相反，模型中的塑性剪胀对应软化和屈服面的收缩。
dp′c / p′c = [v / (λ − κ)] dεvp
这是标准 MCC 的体积硬化微分关系。λ、κ 是压缩与卸载曲线的斜率，下一节给出定义；λ > κ > 0。它们不是这里的塑性乘子 Λ。
在纯各向同性压缩—卸载历史中，p′c 可理解为此前到达过的最大平均有效压力。推广到剪切加载后，它更准确的含义是当前等效预固结压力：剪切引发的塑性压密也能使它变化，不一定伴随新的历史最大 p′。
如果只是从屈服面上的各向同性状态卸载，未产生塑性变形，p′ 会减小，而 p′c 保持不变。这时材料进入超固结状态，当前压力与过去形成的压密边界便分开了。
05 / 固结线从孔隙变化，读出压缩曲线的两种斜率
图 05｜孔隙比、孔隙度和各向同性压缩之间的联系。原图以 φ 表示孔隙度；为避免与摩擦角混淆，正文统一用 n。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.21 原图。
对各向同性正常固结状态，q = 0，且 p′ = p′c。在 v–ln p′ 坐标中，标准模型将这组状态近似为一条直线，即正常固结线；在本讲的加载条件下，也称各向同性正常压缩线。
v = N − λ ln(p′ / pref)
N 为参考压力 pref 处的正常固结比体积，λ 为正常压缩斜率。对数的自变量必须无量纲，因此明确写出 p′/pref。图中出现的 ICL 指各向同性压缩线。
从线上某点 A 卸载后，材料沿较平缓的回弹线变化。用 A 点作锚点，可以避免把卸载线的截距误认为始终不变的材料常数：
v = vA − κ ln(p′ / p′A)
κ 是弹性卸载—再加载斜率。卸载时 p′ < p′A，对数为负，v 增加；再加载且尚未重新屈服时，沿该线返回。
λ 比 κ 大，意味着首次正常压缩的体积变化显著大于同一压力区间内的弹性回弹。两者之差 λ − κ 对应不可逆压密的贡献，也解释了它为什么进入硬化关系。
参数识别 / 坐标轴必须一起读
如果同一压缩曲线用 e–log₁₀p′ 表示，斜率为 Cc，换到 e–ln p′ 后有 λ = Cc/ln10；回弹斜率同理。e 与 v 只差 1，斜率相同。但是，一维侧限固结与各向同性压缩不是同一加载路径，不能仅换对数底就把两种试验参数无条件互换。
06 / 三维联系三条线，怎样描述同一个材料状态？
图 06｜把孔隙比作为第三个坐标，才能同时看见压缩历史、屈服面大小和临界状态。仅在 p′–q 平面投影，会隐藏材料的疏密信息。
配图：D. Nicolas Espinoza 公开课程讲义，Figure 4.22 原图。
到这里，可以把三者的分工说清楚：临界状态线约束持续剪切的最终状态；屈服椭圆给出当前弹性范围的边界；固结线则把压力、体积和压密历史联系起来。它们描述的是同一套模型的不同侧面。
临界状态不仅在 p′–q 平面上满足 q = Mp′，还在比体积平面上满足：
v = Γ − λ ln(p′ / pref)
Γ = N − (λ − κ) ln2
标准 MCC 中，临界状态线与正常固结线在 v–ln p′ 平面上平行。Γ 是临界状态线在参考压力处的比体积；在这些约定下，它与 N、λ、κ 存在上述约束。
这个 ln2 从哪里来？
取一个预固结压力为 p′c 的状态。在正常固结线上，它的比体积为 N − λln(p′c/pref)。同一屈服椭圆的临界点位于 p′ = p′c/2，因此可先沿弹性卸载线将平均压力减半，再在恒平均有效应力下弹性剪切到椭圆顶部。
v = N − λ ln(2p′ / pref) + κ ln2
= [N − (λ − κ) ln2] − λ ln(p′ / pref)
推导使用各向同性弹性关系中体积与偏应力的解耦，并采用标准 MCC 的体积硬化假设。ln2 来自椭圆顶部与右端点之间固定的二倍压力关系。
这一步很关键：Γ 并不是可以在同一标准模型中任意单独调节的截距。如果使用另一套归一化、修改了压缩规律或引入结构性，应先重新核对参数关系。
07 / 完整算例从 20 MPa 卸载到 10 MPa，再开始剪切
教学假设 / 数值用于理解模型
取 M = 1.2，N = 1.80，λ = 0.12，κ = 0.025，pref = 1 MPa。初始为各向同性正常固结状态：p′ = p′c = 20 MPa，q = 0。以下采用标准 MCC 假设，不是某块实测岩土的参数。
第一步：计算正常固结状态的比体积
v20 = 1.80 − 0.12 ln(20/1)
≈ 1.44051
压力比无量纲；孔隙比 e20 = v20 − 1 ≈ 0.44051。
第二步：卸载，改变当前应力而不改变压密边界
v10 = v20 + 0.025 ln(20/10)
≈ 1.45784
当前 p′ 降到 10 MPa，体积略有回弹；弹性卸载使 p′c 仍为 20 MPa。此各向同性状态的超固结比 OCR = p′c/p′ = 2。
此时 e ≈ 0.45784，孔隙度 n = e/(1 + e) ≈ 31.41%。卸载没有使材料退回到 p′ = 10 MPa 的正常固结线；材料保留了经历 20 MPa 压密的历史。
第三步：保持 p′ = 10 MPa，增大 q
qy = M√[p′(p′c − p′)]
= 1.2√[10 × (20 − 10)]
= 12 MPa
qy 为在该固定平均有效应力下首次触及当前屈服面的偏应力。这里控制的是平均有效应力恒定，不是常规三轴试验中的围压恒定。
先检查 q = 6 MPa：f = 6² + 1.2² × 10 × (10 − 20) = −108 MPa²，仍在屈服面内部。增至 q = 12 MPa 时 f = 0，且 p′ = p′c/2，恰好触及椭圆顶部。
第四步：再用体积条件核对临界状态
Γ = 1.80 − (0.12 − 0.025) ln2
≈ 1.73415
vCSL = Γ − 0.12 ln10 ≈ 1.45784
结果与卸载后算出的 v10 一致，同时 q/p′ = 12/10 = M。应力、体积与屈服面三个条件相互吻合。
因此，在这个理想模型和指定路径下，试样可以在该点进入持续塑性剪切，而不再发生塑性体积变化。这比“看到 q/p′ = M 就宣布破坏”多了一次关键核对：材料的体积状态也必须一致。
08 / 课堂自检会读图，也要会识别模型的边界
1．三个主有效应力相等，能否发生塑性变形？
能。此时 q = 0，但各向同性加载达到 p′ = p′c 后，仍可产生塑性压密与硬化。MCC 的封闭压密边界正是为描述这种行为而设置的。
2．卸载后 p′ 减小，屈服椭圆一定同步缩小吗？
不一定。纯弹性卸载不改变 p′c，因此当前应力点移动，而屈服面大小保持不变。椭圆的变化取决于硬化变量，不能只看当前平均应力。
3．到达椭圆顶部，所有体积应变都必定为零吗？
要区分塑性与总量。顶部的关联流动方向给出零塑性体积应变增量；如果同时改变平均有效应力，仍可能有弹性体积变化。持续临界状态还要求应力与比体积保持稳定。
4．临界状态意味着材料没有应变了吗？
恰好相反。剪切应变可以持续累积，只是应力与体积不再继续演化。它描述稳定流动状态，不是变形停止。
5．得到一组漂亮的椭圆，就能用于所有储层岩石吗？
不能。标准 MCC 不自动描述胶结破坏、各向异性、颗粒破碎、裂隙张开及独立抗拉强度。用于泥岩、弱胶结沉积物或压密储层前，应分别检验压缩、回弹、剪胀、峰后响应和应力路径，而不是只拟合一条曲线。
带走的认识
材料当前受到的压力由 p′ 描述，过去形成的压密边界由 p′c 描述，疏密程度由 v 描述。临界状态线、屈服椭圆与固结线必须共同满足，才能把应力判断变成一套自洽的变形解释。
课程对应主题：Advanced Geomechanics（PGE 383），L32：Cam-Clay 模型（上）——临界状态线、屈服面与各向同性固结线。本篇依据作者对应公开课程讲义组织讲解，配图取自原讲义及局部；公式用可复 制文字重新排版，数值算例为教学假设与独立演算。
图文依据：D. Nicolas Espinoza，Advanced Geomechanics · Inelasticity，Critical State Soil Mechanics，Figures 4.16、4.18—4.22。
模型关系核对：Itasca · Modified Cam-Clay Model，压缩曲线、硬化关系与临界状态体积关系。
