3、应力光学原理:平面应力与光弹性效应、应力与折射率张量的关系、应力光学常数标定方法
光学应力测量(光弹性法)的核心,在于利用透明材料在受力状态下产生的双折射现象来反推内部应力状态。要真正掌握这门技术,不能只停留在“会看条纹”的层面,而必须理解应力如何改变材料的微观光学属性,以及如何将光学响应定量地转换为应力数值。本章将围绕三个关键问题展开:平面应力状态如何与光弹性效应建立联系、应力与折射率张量之间的数学物理关系、以及如何通过实验标定获得可靠的应力光学常数。这三者构成了从“看到条纹”到“测准应力”的完整逻辑链。
上图展示了本章三个核心模块之间的逻辑关系以及从实验到定量计算的工程链路。可以看到,平面应力与光弹性效应是物理基础,应力与折射率张量的关系是理论桥梁,而应力光学常数标定则是将理论应用于实际测量的关键步骤。三者环环相扣,缺一不可。下面我们逐一深入展开。
3.1 平面应力状态与光弹性效应的物理基础
核心原理:光弹性法最典型的应用场景是平面应力问题,即构件厚度方向的应力σ₃(或σz)远小于面内应力,可近似视为零。此时,构件处于二维应力状态,任意一点的应力状态由两个主应力σ₁和σ₂(σ₁≥σ₂)及其方向角θ完全确定。当一束偏振光垂直入射到受力的透明模型时,由于应力引起的双折射效应,光波会沿两个主应力方向分解为振动方向互相垂直、传播速度不同的两束光。这两束光通过厚度为h的模型后产生光程差δ,其大小与主应力差(σ₁-σ₂)和模型厚度h成正比,这就是著名的应力-光学定律:
应力-光学定律(平面应力形式)
δ = C · h · (σ₁ - σ₂) 或 N = (δ/λ) = (C·h/λ)·(σ₁ - σ₂)
其中:δ为光程差(nm),C为应力光学常数(布儒斯特,1 Brewster = 10⁻¹² m²/N),h为模型厚度(m),λ为光源波长(nm),N为条纹级数(无量纲)。
在圆偏振光场中观察时,光程差表现为等色线(同一条线上主应力差相等),而在平面偏振光场中,还会出现等倾线(同一条线上主应力方向相同)。等色线给出的是(σ₁-σ₂)的分布,等倾线给出的是主应力方向角θ的分布。两者结合,再辅以分离技术(如剪应力差法、有限元辅助法),才能获得完整的应力分量。
判断方法或操作流程:在实际光弹实验中,判断平面应力状态是否成立,通常遵循以下流程:
- 几何检查:模型厚度h远小于面内特征尺寸(通常h < 1/10 最小面内尺寸),且厚度方向无外载荷。
- 边界检查:模型自由边界上无面外法向力,加载方式为面内加载。
- 光学验证:旋转模型观察条纹是否随旋转角度变化——若条纹不随旋转变化,说明厚度方向应力可忽略。
- 数值估算:用有限元快速估算σ₃/σ₁比值,若小于5%则可按平面应力处理。
案例拆解或实战场景:以带孔平板拉伸为例(经典光弹教学试件),平板宽100mm、厚5mm、中心孔径20mm,材料为环氧树脂(C≈2.2 Brewster),光源波长546.1nm。在拉伸载荷F=5kN作用下,孔边最大应力集中处出现第3级条纹(N=3)。根据应力-光学定律,可快速估算孔边最大主应力差:
σ₁ - σ₂ = N·λ / (C·h) = 3 × 546.1nm / (2.2×10⁻¹² m²/N × 0.005m)
= 1638.3 / (1.1×10⁻¹⁴) ≈ 1.49×10⁸ Pa = 149 MPa
该值与理论应力集中系数Kt=3的预测值(名义应力50MPa × 3 = 150MPa)高度吻合,验证了光弹测量的准确性。
经验提醒:条纹级数N的判读是光弹实验中最容易出错的环节。建议采用“零级条纹追踪法”——从模型的自由边界(应力为零处)开始数起,每经过一条完整的明暗交替即增加0.5级。务必使用白光光源辅助判断零级条纹(零级条纹呈黑色或灰色,非零级呈彩色)。
常见误区或总结复盘:初学者常犯的错误是将等色线直接当作主应力差云图,忽略了等色线只代表(σ₁-σ₂)的等值线,并不包含主应力方向信息。另外,平面应力假设不适用于厚度方向有约束的场合(如厚壁管、接触问题),此时必须采用三维光弹性或数值方法。总结:平面应力是光弹法的适用前提,应力-光学定律是定量测量的核心公式,条纹判读是实验技能的关键。
3.2 应力与折射率张量的关系:从光率体到应力-光学定律
核心原理:各向同性透明材料在无应力状态下,其折射率在各个方向相同,折射率椭球(光率体)为球体。当材料受到应力作用时,原子间的键长和键角发生变化,导致材料的介电张量改变,折射率椭球发生变形,由球体变为三轴椭球。这一现象可以用应力-光学张量关系精确描述。
在三维情况下,折射率变化与应力之间的关系为:
ΔBᵢⱼ = qᵢⱼₖₗ · σₖₗ
其中Bᵢⱼ为逆介电张量(B = 1/n²),qᵢⱼₖₗ为四阶光弹张量(压光系数),σₖₗ为应力张量。对于各向同性材料,独立的压光系数只有两个:q₁₁和q₁₂(或q₄₄ = (q₁₁-q₁₂)/2)。在平面应力状态下,沿两个主应力方向的折射率变化为:
n₁ = n₀ + q₁₁·σ₁ + q₁₂·σ₂
n₂ = n₀ + q₁₂·σ₁ + q₁₁·σ₂
两束光通过厚度h的模型后,光程差为:
δ = h·(n₁ - n₂) = h·(q₁₁ - q₁₂)·(σ₁ - σ₂) = C·h·(σ₁ - σ₂)
由此得到应力光学常数C = q₁₁ - q₁₂。这就是应力-光学定律的张量推导过程,它揭示了光弹效应的微观本质——应力通过改变材料电子云的极化率,进而改变折射率,最终表现为宏观的光程差。
判断方法或操作流程:在实际工程中,我们通常不需要直接测量q₁₁和q₁₂,而是通过标定实验直接获得C值。但理解张量关系有助于判断测量的适用性和误差来源:
- 材料各向同性验证:对未加载试件进行光学检查,确认无初始双折射(残余应力)。
- 线性范围确认:在标定实验中,改变载荷水平,绘制N-Δσ曲线,检查线性度(R²>0.999)。
- 波长依赖性检查:不同波长下C值略有差异(色散效应),高精度测量需使用单色光源。
- 温度敏感性评估:光弹材料对温度敏感,标定和使用应在相同温度条件下进行(±1°C)。
案例拆解或实战场景:某航空透明件(聚碳酸酯)需要评估其内部应力水平。已知该材料在25°C、波长546.1nm下的应力光学常数C = 3.8 Brewster。在偏振光场中观察到某区域条纹级数N = 2.5,试件厚度h = 8mm。计算主应力差:
σ₁ - σ₂ = N·λ / (C·h) = 2.5 × 546.1nm / (3.8×10⁻¹² × 0.008m)
= 1365.25 / (3.04×10⁻¹⁴) ≈ 44.9 MPa
该值接近聚碳酸酯屈服强度的60%,需要进一步评估是否在安全范围内。通过这个案例可以看出,准确的材料常数C直接决定了测量结果的可靠性。
风险说明:应力光学常数C并非常数,它随温度、波长、材料老化程度而变化。对于高精度测量(误差<5%),必须在与实际测量相同的环境条件下进行标定。切勿直接引用文献值而不加验证,特别是对于复合材料或改性材料。
常见误区或总结复盘:一个常见误区是将折射率椭球的变化简单理解为折射率各向同性变化,忽略了应力引起的双折射是差值效应——只有主应力差(σ₁-σ₂)才产生光程差,而平均应力(σ₁+σ₂)/2只改变折射率,不产生双折射。因此,光弹法天然无法直接测量平均应力,需要借助其他方法(如全息干涉、莫尔干涉)或数值方法辅助。总结:张量关系是理论基石,C值是工程应用的关键参数,理解其物理意义有助于正确设计实验和解释数据。
3.3 应力光学常数标定方法:从实验设计到数据处理
核心原理:应力光学常数C是光弹测量的“标尺”,其精度直接决定测量结果的可靠性。标定的本质是:在已知应力分布的条件下,测量光弹条纹级数,反推C值。最常用的标定方法包括对径压缩圆盘法(巴西圆盘法)和四点弯曲梁法。其中,对径压缩圆盘法因其应力分布有精确解析解而被广泛采用。
对于直径为D、厚度为h的圆盘,沿直径方向施加压缩力F时,圆盘中心处的应力状态为:
σ₁ = 2F / (π·D·h) (拉应力,水平方向)
σ₂ = -6F / (π·D·h) (压应力,垂直方向)
σ₁ - σ₂ = 8F / (π·D·h)
因此,只要在圆盘中心处读取条纹级数N₀,即可计算:
C = N₀·λ·π·D / (8·F)
判断方法或操作流程:标准标定流程如下:
| 步骤 | 操作内容 | 关键参数/注意事项 |
|---|---|---|
| 1. 试件制备 | 加工圆盘试件,直径D=40~60mm,厚度h=4~8mm | 边缘光滑无缺口,厚度均匀(±0.02mm),退火消除残余应力 |
| 2. 加载装置 | 使用万能试验机或专用加载架,加载方向沿直径 | 加载点加垫片防止局部压碎,力值精度±0.5% |
| 3. 光学系统 | 圆偏振光场,光源波长λ=546.1nm(汞灯绿光) | 确保起偏镜和检偏镜正交,1/4波片角度正确 |
| 4. 条纹读取 | 在圆盘中心区域读取条纹级数N₀ | 中心处条纹最清晰,使用补偿器可提高读数精度至±0.02级 |
| 5. 多级加载 | 从F₁到F₅逐步加载(如500N、1000N、1500N、2000N、2500N) | 每级记录N₀,绘制N₀-F曲线,检查线性度 |
| 6. 数据处理 | 线性拟合N₀ = k·F,计算C = k·λ·π·D / 8 | 拟合R²应大于0.999,取多次测量的平均值 |
案例拆解或实战场景:某实验室需要对一批环氧树脂板材进行标定。圆盘试件D=50.00mm,h=5.02mm,使用546.1nm光源。加载数据如下:
F (N) 500 1000 1500 2000 2500
N₀ 0.82 1.65 2.47 3.30 4.12
线性拟合得到斜率k = 0.00165 级/N(R²=0.9999)。计算C值:
C = k·λ·π·D / 8 = 0.00165 × 546.1×10⁻⁹ × π × 0.050 / 8
= 0.00165 × 546.1×10⁻⁹ × 0.19635 / 8
= 0.00165 × 1.072×10⁻⁸ / 8
= 2.21×10⁻¹² m²/N = 2.21 Brewster
该值与文献报道的环氧树脂典型值(2.0~2.5 Brewster)吻合。进一步分析不确定度:力值误差±0.5%、厚度误差±0.4%、条纹读数误差±0.02级(约±0.5%)、直径误差±0.1%,合成标准不确定度约为±1.2%。
经验提醒:标定实验中,圆盘中心条纹级数的读取是关键。建议使用Tardy补偿法或旋转检偏镜法将读数精度提高到±0.01级。同时,加载速度要缓慢均匀,避免粘弹性效应引起的条纹漂移。每次加载后等待30秒再读数,确保蠕变稳定。
常见误区或总结复盘:标定中常见的误区包括:①忽略温度影响——光弹材料C值随温度变化可达0.5%/°C,标定和使用必须在相同温度下进行;②加载点应力集中——圆盘加载点附近存在高应力梯度,条纹密集难以判读,应只读取中心区域;③厚度测量位置不当——应测量圆盘中心处的厚度而非边缘;④忽略初始双折射——试件制备过程中可能引入残余应力,标定前必须进行退火处理并验证无初始条纹。总结:标定不是一次性的工作,建议每次实验前对同一批材料进行抽检标定,确保C值的时效可靠性。
本章核心要点回顾:
- 平面应力状态(σ₃≈0)是光弹法应用的前提,应力-光学定律δ = C·h·(σ₁-σ₂)是定量测量的基础。
- 应力通过改变折射率椭球形状产生双折射,各向同性材料的光弹行为由两个独立压光系数决定,C = q₁₁-q₁₂。
- 对径压缩圆盘法是最常用的标定方法,利用圆盘中心精确解析解σ₁-σ₂ = 8F/(πDh)反推C值。
- 标定精度受温度、波长、加载速度、读数误差等多因素影响,合成不确定度应控制在±2%以内。