B3 · 出版卷 8
应力、应变与变形机制
Stress tensor intuition、strain、elastic/brittle/ductile
学习目标
完成本课后，学习者应能区分力、牵引力、应力、位移和应变；分解一个面上的正应力和剪应力；从概念上区分平均应力和偏应力；比较弹性、脆性、摩擦和晶体塑性响应；解释温度、压力、流体压力、应变速率、粒度和矿物组成为何会改变岩石保存的变形机制。
从野外问题开始
一条陡倾破碎带切过能干砂岩，向下延伸到片理化片岩中。砂岩段含角砾和张开脉体，片岩段则变宽为动态重结晶颗粒和连续弯曲面理组成的带状体。这是同一构造在不同条件下形成的连续记录、两个后来连接的构造，还是露头造成的假象？
“岩石受到挤压”并不足够。挤压可以是静水式，也可以具有差应力；对一个面是正向载荷，对另一个面可能是斜向剪切。同一套区域边界条件可以使某一方向张裂、另一方向剪切，并让高温或弱层发生韧性流动。现存产物还整合了完整路径：破裂、滑动、愈合、埋藏、升温和再次变形都可能发生在同一继承带上。
应先记录接触、宽度、产状、错距、碎屑形态、脉体、粒度梯度以及断层岩与面理的关系，再区分哪些证据约束几何、位移、机制和时间。应力通常通过模型反演；应变和显微构造更接近直接观测。
核心过程模型
应力是在足够小面积极限上的单位面积力，并以二阶张量表示。对于单位法向量 \mathbf{n}，该面上的牵引力为
\mathbf{t}(\mathbf{n})=\boldsymbol{\sigma}\mathbf{n}.
正应力分量为 \sigma_n=\mathbf{n}\cdot\mathbf{t}，剪应力矢量为 \boldsymbol{\tau}=\mathbf{t}-\sigma_n\mathbf{n}。因此，同一个应力张量会在不同方向的面上产生不同正应力和剪应力。主方向上的剪应力为零。平均应力主要影响体积和孔隙，偏应力推动形状改变，但天然材料的本构响应会耦合二者。
应变描述相对位移。小变形条件下，正应变表示相对长度变化，剪应变表示角度变化。有限变形需要变形梯度；不同增量的先后顺序可能改变结果。刚体平移和刚体旋转改变位置或方向，但本身不是应变。
弹性应变在适用范围内卸载后可以恢复。脆性破裂形成或扩展不连续面；摩擦滑动把位移集中于已有面；碎裂流动结合颗粒破碎、旋转和滑移；压溶通过化学势梯度迁移物质；位错蠕变、扩散蠕变和晶界滑动形成晶体塑性或黏性变形。多种机制可同时工作，最终保存的显微构造未必对应位移最大的阶段。
围压抑制张开并促进分布式变形；升温通常加速热活化过程；连通孔隙中的流体压力降低有效正应力；应变速率改变蠕变和反应可用时间；粒度影响扩散和晶界作用；矿物组成控制强度、各向异性和反应。因此，“脆—韧转换”不存在一个普适固定深度。
证据与测量
露头尺度可测量裂隙间距、断层带宽度、标志层错距、脉体张开量、褶皱波长、面理强度和横穿构造带的梯度。手标本和薄片尺度应记录颗粒破裂、粉碎基质、波状消光、亚晶、新生颗粒、双晶、压力影、溶蚀缝、矿物纤维和切割关系。所有运动学图像都必须带方向和观察面说明。
颗粒形状本身具有多解性。长形颗粒可能是沉积、岩浆、变质或变形产物；细粒基质可由碎裂、反应、蚀变或原始粒度造成。动态重结晶判据取决于矿物、温度、水活度和后期退火。应把显微证据与地图尺度梯度及矿物化学结合，而不能只凭一张照片命名机制。
实验室强度、摩擦或流动律只适用于其边界条件。与天然岩石比较前，应记录岩性、孔隙度、饱和度、围压、温度、加载路径、应变速率、样品方向和尺度。跨越多个时间或长度数量级的外推是带不确定性的模型，而不是直接测量。
还应检查同一构造带内部的材料差异。石英、云母、长石和黏土矿物可以在相同边界条件下采用不同机制，使应变在矿物和层间分配。显微尺度上的弱相连通可能控制宏观强度，而少量强颗粒仍保存早期应变。把机制图与实测矿物组成、含水状态和粒度范围逐项对应，比用单一“脆性”或“韧性”标签更可靠。
计算示例
设二维压应力张量单位为兆帕：
\boldsymbol{\sigma}=
\begin{bmatrix}
100 & 30\\
30 & 60
\end{bmatrix}.
若某面的单位法向量为 \mathbf{n}=(0.707,0.707)，牵引力约为 \mathbf{t}=(91.9,63.6) MPa。正应力为 \sigma_n=\mathbf{n}\cdot\mathbf{t}=110 MPa，扣除 \sigma_n\mathbf{n} 后的剪应力大小约为 20 MPa。这个计算说明：张量的某个坐标分量不能直接当作地图裂隙面上的应力。
若连通流体压力为 70 MPa，则简单有效应力近似给出 \sigma'_n=40 MPa。设黏聚力 c=5 MPa、摩擦系数 \mu=0.6，摩擦阻力为
\tau_f=c+\mu\sigma'_n=5+0.6(40)=29\ \text{MPa}.
剪应力 20 MPa 低于这个模型的阈值。提高流体压力、降低黏聚力、旋转裂隙面或改变应力张量都可能改变结论。粗糙度、尺度、瞬态压力、化学弱化和动态作用被简化，因此计算不能证明该面一定不活动。
误判与不确定性
不能直接根据褶皱或断层画出应力箭头。构造形成后可能旋转，可能受早期各向异性控制，也可能在非共轴流动中形成或被后期应力场再活化。现代原位应力不等于古应力。
不能简单等同“韧性—深部”“脆性—浅部”，也不能从碎裂岩唯一推断滑动速率或深度。有效压力、温度、矿物、流体、粒度和速率共同控制响应。大的韧性带中可以包含脆性增量，裂隙附近也可发生韧性蠕变。“能干”和“软弱”是比较性、尺度相关描述，不是固定常数。
应力大小通常比构造几何约束更弱。用断层擦痕或地震数据反演应力依赖均一应力、滑动方向和断层相互作用等假设。报告时应分别标明测量量、计算量和假设量。
实践调查
画一个方形应力单元，选择法线为 0°、30°和 60°的三个面。对一个对称二维应力张量分别计算牵引力、正应力和剪应力，并在正应力—剪应力图上表示，再用圆图法复核。改变两档流体压力后重复计算。
随后对一个露头、岩芯图像或薄片组进行三个尺度的观察，建立“观测—机制假设—检验”三张独立表。每次赋予机制名称前至少需要两项证据，并为每个诊断纹理写出一个可能的模拟或擦除过程。
掌握检查
- 为什么一个区域箭头不能描述任意面上的应力？
- 计算应变时排除了哪些位移分量？
- 流体压力如何影响简单摩擦破坏计算？
- 为什么脆—韧转换不是固定深度？
- 哪些观测可区分碎裂造成的细粒化与原生细粒？