3、应力基本概念：应力定义、应变、弹性模量、应力单位及薄膜应力与体材料应力的区别
光学薄膜的应力问题贯穿于薄膜设计、制备工艺、器件集成与可靠性评估的全流程。无论是镀膜后的翘曲变形、膜层开裂，还是长期服役中的性能漂移，其根源都指向薄膜内部存在的残余应力。要系统掌握应力控制与消除方法，首先必须建立清晰、准确的应力基本概念体系。本章从材料力学的基础定义出发，逐步引入应变、弹性模量与应力单位，并重点剖析薄膜应力与体材料应力在物理本质、约束状态和工程表现上的显著差异，为后续章节的应力测量、工艺调控与失效分析奠定理论基础。
3.1 应力与应变的核心定义
3.1.1 应力的物理本质与数学表达
应力（Stress）描述的是材料内部单位面积上所承受的内力，其定义为 σ = F / A，其中 F 为作用于截面上的内力，A 为截面积。在连续介质力学框架下，应力是一个二阶张量，完整描述需要九个分量。但在薄膜力学分析中，通常关注的是正应力（法向应力）和剪应力（切向应力）两类。正应力垂直于截面，剪应力平行于截面。对于光学薄膜，由于膜层厚度远小于面内尺寸，应力状态可简化为平面应力问题，即仅考虑面内两个正应力分量 σx 和 σy，以及面内剪应力 τxy。
应力的符号约定在工程应用中极为重要：张应力（Tensile Stress）取正值，表示材料受拉，有收缩趋势；压应力（Compressive Stress）取负值，表示材料受压，有膨胀趋势。在薄膜体系中，张应力易导致膜层开裂，压应力易导致膜层起皱或剥离，这一符号判断是后续失效分析的基础。
3.1.2 应变与弹性变形的关联
应变（Strain）是材料在应力作用下发生的相对形变，定义为 ε = ΔL / L₀，其中 ΔL 为长度变化量，L₀ 为原始长度。应变是无量纲量，工程中常用微应变（με，即 10⁻⁶）表示。应变分为弹性应变和塑性应变：弹性应变在应力卸载后完全恢复，遵循胡克定律；塑性应变则产生永久变形。对于光学薄膜，绝大多数服役工况要求膜层处于弹性变形范围内，一旦进入塑性区，膜层的光学性能和力学完整性将不可逆地劣化。
在薄膜体系中，应变与应力通过弹性模量建立联系。需要特别指出的是，薄膜的弹性模量往往与体材料不同，这源于薄膜的微结构特征——柱状晶生长、晶界密度高、孔隙率大、以及可能的非晶态结构，都会导致弹性模量偏离体材料值。例如，电子束蒸发沉积的 SiO₂ 薄膜，其弹性模量约为体材料（约 73 GPa）的 60%~80%。
3.1.3 弹性模量的物理意义与测试方法
弹性模量（Elastic Modulus）是材料抵抗弹性变形的能力指标，在薄膜力学中主要涉及以下三种模量：
| 模量类型 | 符号 | 物理意义 | 薄膜测试方法 |
|---|---|---|---|
| 杨氏模量 | E | 单轴拉伸时应力与应变之比 | 纳米压痕、悬臂梁弯曲 |
| 剪切模量 | G | 剪应力与剪应变之比 | 表面声波法、布里渊散射 |
| 体积模量 | K | 静水压力与体积应变之比 | 高压X射线衍射 |
对于各向同性材料，E、G、K 之间存在关系：E = 2G(1+ν) = 3K(1-2ν)，其中 ν 为泊松比。在薄膜应力计算中，常使用双轴模量 M = E / (1-ν)，这是因为薄膜处于双轴应力状态，面内两个方向同时受约束，其等效刚度高于单轴拉伸情形。
3.2 应力单位与换算体系
3.2.1 国际单位制与常用工程单位
应力的国际单位制（SI）单位为帕斯卡（Pa），定义为 1 N/m²。在薄膜领域，由于应力值通常在数十至数千兆帕范围，常用单位及其换算关系如下：
- 1 MPa = 10⁶ Pa = 1 N/mm²，薄膜应力最常用单位
- 1 GPa = 10⁹ Pa = 1000 MPa，用于高应力硬质薄膜
- 1 dyn/cm² = 0.1 Pa，CGS单位制，部分文献仍在使用
- 1 kgf/mm² ≈ 9.8 MPa，工程单位，需注意换算
- 1 psi ≈ 6.895 kPa，英制单位，部分进口设备规格书中出现
3.2.2 薄膜应力数值的量级感知
为了建立直观的量级概念，以下给出典型光学薄膜的应力范围：电子束蒸发 SiO₂ 薄膜通常呈现压应力，约 -100 ~ -400 MPa；离子辅助沉积（IAD）的 TiO₂ 薄膜应力可在 -200 ~ +300 MPa 之间调控；磁控溅射的金属薄膜（如 Ag、Al）应力范围较宽，从 -500 MPa 到 +500 MPa 不等。当应力超过 1 GPa 时，绝大多数光学薄膜已接近或超过其断裂强度，极易发生开裂或剥离。
-350 MPa = -0.35 GPa
-350 MPa = -350 × 10⁶ Pa = -350 × 10⁷ dyn/cm² = -3.5 × 10⁹ dyn/cm²
在撰写或与国外同行交流时，建议统一使用 MPa 或 GPa，避免单位混淆。
3.3 薄膜应力与体材料应力的本质区别
3.3.1 约束状态与应力维度的差异
体材料应力分析通常假设材料在三个方向自由变形，应力状态可以是单轴、双轴或三轴，取决于外加载荷方式。而薄膜沉积在基板上，面内方向受到基板的刚性约束，厚度方向自由。这种约束状态导致薄膜内部形成双轴应力场：面内两个方向存在应力分量，厚度方向应力近似为零（薄膜表面自由）。这一差异的工程后果是：薄膜的应力-应变关系必须使用双轴模量 M = E/(1-ν) 来描述，而非简单的杨氏模量 E。
3.3.2 应力来源与演化机制的差异
体材料的应力通常来源于外部载荷（机械力、热载荷）或加工残余应力，其分布相对均匀，可通过热处理或机械加工消除。薄膜应力则复杂得多，主要包含以下两类：
- 热应力（Thermal Stress）：由于薄膜与基板热膨胀系数（CTE）不匹配，在温度变化时产生。其表达式为 σ_thermal = E_f/(1-ν_f) × (α_s - α_f) × ΔT，其中 α_s 和 α_f 分别为基板和薄膜的 CTE。热应力可通过工艺温度控制来调节。
- 本征应力（Intrinsic Stress）：源于薄膜生长过程中的微观机制，包括晶格失配、晶界弛豫、杂质嵌入、相变、空位湮灭等。本征应力与沉积工艺参数（能量、速率、气压）密切相关，是薄膜应力调控的核心对象。
3.3.3 尺度效应与微结构影响
薄膜厚度通常在纳米至微米量级，其微结构与体材料截然不同。薄膜生长初期形成岛状结构，随后合并成连续膜，这一过程中产生大量晶界、位错和孔隙。这些缺陷在应力场中充当应力集中源，导致薄膜的力学性能呈现显著的尺度效应：强度随厚度减小而升高（Hall-Petch效应在纳米尺度可能反转），弹性模量随厚度变化，断裂韧性远低于体材料。此外，薄膜与基板的界面结合强度、界面处的应力传递效率，都直接影响薄膜应力的实际表现。
3.3.4 应力梯度与厚度依赖性
体材料应力在宏观尺度上可视为均匀分布，而薄膜应力沿厚度方向存在显著梯度。这源于薄膜生长过程中微结构随厚度的演化：靠近基板界面处，薄膜受基板晶格和表面粗糙度影响，应力状态复杂；膜层中部趋于稳态生长，应力相对均匀；表面层则受环境氧化、吸附等影响。应力梯度的存在使得简单的平均应力值不足以描述薄膜的力学行为，必须结合应力分布剖面进行分析。
3.4 薄膜应力工程分析案例：SiO₂/TiO₂ 多层膜应力失配问题
3.4.1 案例背景与初始现象
某激光薄膜项目需要制备 1064nm 高反膜，膜系结构为交替沉积的 SiO₂/TiO₂ 多层膜，共 21 层，总厚度约 3.5 μm，基板为直径 25.4 mm、厚度 1 mm 的熔石英。采用离子辅助电子束蒸发工艺，沉积温度 200°C。镀膜完成后，使用激光干涉仪测量基板面形，发现基板严重翘曲，中心凹陷深度达 12 μm（初始平面度要求小于 2 μm），且膜层表面出现微裂纹。
3.4.2 关键数据与应力分析
通过基板曲率法（Stoney 公式）计算薄膜平均应力：σ = E_s × t_s² / (6 × (1-ν_s) × t_f × R)，其中 E_s = 72 GPa，ν_s = 0.17，t_s = 1 mm，t_f = 3.5 μm，R 为基板曲率半径。实测平均应力为 -280 MPa（压应力）。进一步使用纳米压痕法分别测量单层 SiO₂ 和 TiO₂ 薄膜的弹性模量与硬度：SiO₂ 薄膜 E = 68 GPa，TiO₂ 薄膜 E = 105 GPa，均低于对应体材料值（SiO₂ 体材料 73 GPa，TiO₂ 体材料 230 GPa）。
3.4.3 分析判断与问题定位
根据膜系结构，SiO₂ 层总厚度约 2.1 μm，TiO₂ 层总厚度约 1.4 μm。若假设各层应力均匀，则平均应力应为各层应力的厚度加权平均。通过单层膜应力测试发现：SiO₂ 单层膜应力约 -150 MPa（压应力），TiO₂ 单层膜应力约 +180 MPa（张应力）。按此计算，多层膜平均应力应为：(-150 × 2.1 + 180 × 1.4) / 3.5 ≈ -18 MPa，与实测 -280 MPa 相差甚远。
进一步分析发现，问题出在界面应力贡献和应力梯度效应。多层膜中每层厚度仅约 160 nm（SiO₂）和 110 nm（TiO₂），远小于单层膜测试时的 500 nm 厚度。在如此薄的膜层中，界面效应占比显著增大，且 TiO₂ 在 SiO₂ 表面成核生长时，晶格失配导致额外应力。此外，离子辅助沉积过程中的离子轰击效应在多层膜中逐层累积，使得各层实际应力与单层膜测试值存在系统性偏差。
3.4.4 解决方案与结果对比
针对应力失配问题，采取以下措施：
- 调整沉积参数：将 TiO₂ 层的离子辅助能量从 300 eV 降至 200 eV，降低压应力贡献；同时将 SiO₂ 层的沉积速率从 0.5 nm/s 降至 0.3 nm/s，减少张应力积累。
- 引入应力补偿层：在膜系中插入一层 50 nm 的 MgF₂ 过渡层，利用其张应力特性平衡整体压应力。
- 优化退火工艺：镀膜后在 250°C 下退火 2 小时，促进应力弛豫，使应力分布更均匀。
改进后，基板翘曲从 12 μm 降至 1.8 μm，满足平面度要求；膜层表面微裂纹消失；反射率在 1064nm 处达到 99.8%，满足设计指标。该案例表明，薄膜应力控制必须基于实际工艺条件下的应力数据，而非简单套用体材料参数或单层膜测试结果。
3.4.5 案例复盘与工程启示
本案例揭示了薄膜应力分析的三个关键要点：第一，多层膜应力不是各层应力的简单叠加，界面效应和应力梯度不可忽略；第二，薄膜弹性模量必须实测，体材料参数仅可作为初值参考；第三，应力调控需要从材料选择、工艺参数、膜系结构三个维度协同优化。在实际工程中，建议建立工艺-应力-性能的数据库，通过系统实验积累本单位的薄膜应力特征数据。
3.5 常见误区与总结
| 误区类型 | 错误认知 | 正确理解 |
|---|---|---|
| 模量取值 | 直接使用体材料弹性模量 | 薄膜模量受工艺影响，需实测 |
| 应力状态 | 将薄膜应力视为三轴应力 | 薄膜为双轴应力，厚度方向自由 |
| 应力来源 | 仅考虑热应力 | 本征应力往往占主导地位 |
| 多层膜叠加 | 各层应力线性叠加 | 界面效应和应力梯度导致非线性 |
| 单位混淆 | MPa 与 GPa 混用不换算 | 统一单位，注意数量级差异 |
本章系统梳理了应力、应变、弹性模量等基本概念，明确了应力单位的换算体系，并深入剖析了薄膜应力与体材料应力的本质区别。核心要点可归纳为：薄膜应力是双轴应力状态，由热应力和本征应力共同构成，其弹性模量、应力分布和失效模式均与体材料存在显著差异。掌握这些基本概念，是后续学习应力测量方法、工艺调控策略和失效分析技术的前提。在实际工程中，建议工程师建立"实测优先、数据库支撑、系统优化"的应力管理理念，避免因概念混淆或参数误用导致的薄膜失效问题。