3. 应力光学定律：应力-光学系数、平面应力与双折射的关系、三维应力场的双折射表达

好，咱们进入第三章。这一章是整个应力双折射标注的“心脏”。

前面两章我们聊了光弹性的现象和基本概念。现在要回答一个核心问题：应力到底是怎么变成双折射的？ 这个转化关系，就是应力光学定律。

我个人习惯把这条定律看作一座桥。桥的一头是力学量（应力），另一头是光学量（折射率变化）。没有这座桥，我们测出来的条纹就只是好看的图案，跟应力值对不上号。

3.1 应力光学定律的核心：应力-光学系数

先看最简单的形式。对于各向同性的透明材料（比如玻璃、环氧树脂），在没有应力时，它的折射率在各个方向都一样，记作 n₀。

一旦施加应力，材料内部就变得“各向异性”了。折射率会沿着主应力方向发生变化。这个变化量，跟主应力的大小成正比。

数学上，应力光学定律可以写成：

n₁ - n₀ = C₁ · σ₁ + C₂ · (σ₂ + σ₃)
n₂ - n₀ = C₁ · σ₂ + C₂ · (σ₁ + σ₃)
n₃ - n₀ = C₁ · σ₃ + C₂ · (σ₁ + σ₂)

这里：

• n₁, n₂, n₃ —— 三个主应力方向上的折射率
• σ₁, σ₂, σ₃ —— 三个主应力
• C₁, C₂ —— 应力-光学系数（材料常数）

你想想看，这个公式其实挺直观的。某个方向上的折射率变化，不仅由该方向的应力决定，还受到另外两个方向应力的影响。这就是泊松效应的光学表现。

3.2 平面应力状态：从三维到二维的简化

实际工程中，很多构件处于平面应力状态。比如薄板、壳体、压力容器的壁面。这种情况下，垂直于板面的应力 σ₃ = 0。

好，简化来了。当 σ₃ = 0 时，光沿着板厚方向（即 σ₃ 方向）穿过，我们测到的双折射只跟 σ₁ 和 σ₂ 有关。

光程差 δ 的表达式变得很简单：

δ = C · d · (σ₁ - σ₂)

其中：

• δ —— 光程差（单位：nm）
• C —— 相对应力-光学系数
• d —— 试件厚度
• σ₁ - σ₂ —— 主应力差

说白了，在平面应力下，双折射的大小直接等于主应力差乘以一个常数。这就是为什么光弹条纹图上的等色线，其实就是主应力差的等值线。

3.3 三维应力场的双折射表达

到了三维情况，事情就复杂了。光弹实验通常只能给出二维信息（沿光路积分的平均值）。但理论上，三维应力场的双折射关系需要完整描述。

我们回到最一般的表达式。对于任意应力状态，折射率椭球（又称光率体）的变化由下式描述：

ΔBᵢⱼ = Cᵢⱼₖₗ · σₖₗ

这里：

• ΔBᵢⱼ —— 介电不透张量的变化（与折射率相关）
• Cᵢⱼₖₗ —— 四阶应力-光学张量
• σₖₗ —— 应力张量

对于各向同性材料，这个四阶张量只有两个独立分量。写成矩阵形式（Voigt 记号）就是：

其中 C₄₄ = (C₁₁ - C₁₂)/2。

嗯，这里要注意。三维情况下的双折射不再只是主应力差的函数。切应力也会引起折射率变化，而且方向是旋转的。这就是为什么三维光弹需要“冻结切片法”——把三维应力状态“冻结”在模型里，然后切成薄片，逐片分析。

Δn = C · √[(σₓ - σᵧ)² + 4τₓᵧ²]

3.4 实际应用中的注意事项

讲完理论，说点实际的。我在做应力双折射标注时，有几个坑是反复踩过的：

1. 温度影响：应力-光学系数 C 是温度的函数。室温下测的 C 值，不能直接用到高温环境。我曾经用加热到 80°C 的试件做实验，结果条纹级数比室温下少了 30%。后来查资料才发现，C 值随温度升高而下降。
2. 材料非线性：高应力下，C 值可能不再是常数。对于某些聚合物材料，当应力超过 20 MPa 时，C 值会开始变化。这时候再用线性公式，误差就大了。
3. 光路积分效应：三维应力场中，光沿着路径穿过，测到的是积分效果。如果应力沿厚度方向变化剧烈，测得的双折射是“平均”结果，不代表某一点的局部应力。