3. 有限元方法入门
各位同学，欢迎来到第三章。
前两章我们聊了光机结构分析到底在干什么，也认识了有限元这个工具。但说句实在话，如果你不懂材料力学，那有限元算出来的结果就是一堆数字——你根本不知道它是对是错。
我见过不少新手，模型建得漂漂亮亮，网格剖得整整齐齐，结果一算，应力云图红得发紫，他却问我：“老师，这个值大不大？”
嗯，这就是我们今天要解决的问题。
3.1 应力与应变：最基础的两个概念
先问个问题：一根杆子，你拉它，它内部发生了什么？
简单说，应力是单位面积上承受的内力，应变是单位长度的变形量。
应力 σ = F / A（F 是力，A 是截面积）
应变 ε = ΔL / L（ΔL 是伸长量，L 是原长）
单位：应力常用 MPa（兆帕），应变无量纲。
我个人习惯把应力想象成“内部的压强”，应变想象成“被拉长的比例”。你想想看，一根橡皮筋，拉长一倍，应变就是 1.0（100%）。但光机结构里，应变通常很小，千分之一甚至万分之一。
这里有个坑，我年轻时踩过：应力分正应力和切应力。正应力垂直于截面，切应力平行于截面。有限元分析里，我们经常看的是 von Mises 应力——它把正应力和切应力综合成一个等效应力，用来判断材料是否屈服。
小技巧：在 ANSYS Workbench 或 COMSOL 里，默认输出的应力云图通常就是 von Mises 应力。但如果你关心的是某个方向的拉裂风险，记得切到最大主应力视图。
3.2 胡克定律：线弹性材料的“本构关系”
胡克定律，说白了就是：在弹性范围内，应力和应变成正比。
公式很简单：σ = E · ε
这个 E 就是弹性模量（杨氏模量），单位也是 MPa 或 GPa。它衡量材料抵抗变形的能力。E 越大，材料越“硬”。
举个例子：铝合金的 E 大约是 70 GPa，钢大约是 210 GPa。同样一根杆子，同样拉它，铝合金的变形是钢的 3 倍。
我在做一款激光雷达的转镜支架时，就吃过这个亏。一开始用了铝合金，结果模态频率不够，一振动就偏。后来换成殷钢（E 约 140 GPa），刚度上去了，问题解决。
注意：胡克定律只适用于线弹性阶段。一旦应力超过屈服强度，应力-应变曲线就不再是直线了。有限元分析里，如果你用的是“线性静力分析”，软件默认材料是线弹性的——它不会告诉你材料已经屈服了。你得自己判断。
3.3 泊松比：一个容易被忽略的参数
拉一根杆子，它会变长，同时会变细。这个“变细”的程度，就是泊松比 ν。
ν = - (横向应变) / (纵向应变)
大多数金属的泊松比在 0.3 左右。铝合金约 0.33，钢约 0.3，殷钢约 0.25。
泊松比在有限元分析里有什么用？
- 影响横向变形：如果你做的是薄壁件或梁结构，泊松比会影响弯曲时的截面畸变。
- 影响接触分析：两个零件压在一起，泊松比不同，接触压力分布会不一样。
- 影响热应力：热膨胀时，泊松比会约束横向变形，产生附加应力。
我曾经遇到一个案例：一个光学镜片的压圈，用铝合金做的，拧紧后镜片变形了。我一开始以为是预紧力太大，后来一查，是铝合金的泊松比导致压圈在径向膨胀时产生了额外的轴向力。换成殷钢后，问题消失。
3.4 屈服强度：材料的“底线”
屈服强度 σs 是材料从弹性变形转为塑性变形的临界点。超过这个值，材料就永久变形了。
光机结构里，我们通常要求：最大应力 ≤ 屈服强度 / 安全系数
安全系数怎么取？
- 一般静载：1.5 ~ 2.0
- 振动/冲击：2.0 ~ 3.0
- 光学元件固定：3.0 ~ 5.0（镜片变形直接影响成像质量）
记住：有限元分析算出来的应力，不是真实应力，而是“模型应力”。网格越粗，应力越不准。所以安全系数不是拍脑袋，而是给模型误差留的余量。
3.5 常用光机材料的力学性能
做光机结构，你打交道最多的就是这三种材料：铝合金、殷钢、微晶玻璃。我整理了一张表，建议你收藏。
| 材料 | 弹性模量 E (GPa) | 泊松比 ν | 屈服强度 σs (MPa) | 密度 ρ (g/cm³) | 热膨胀系数 α (10⁻⁶/K) | 主要用途 |
|---|---|---|---|---|---|---|
| 铝合金 6061-T6 | 68.9 | 0.33 | 276 | 2.70 | 23.6 | 结构件、镜筒、支架 |
| 殷钢 4J36 (Invar) | 141 | 0.25 | 280 | 8.10 | 1.2 | 精密镜座、基准件 |
| 微晶玻璃 (Zerodur) | 90.3 | 0.24 | —（脆性材料） | 2.53 | 0.02 | 反射镜基体、标准尺 |
说说我的个人经验：
- 铝合金：便宜、好加工、轻。但热膨胀大，温度一变，尺寸就变。适合做非精密的结构件。
- 殷钢：热膨胀极小，是精密光机结构的“定海神针”。但重、贵、难加工。我只在关键位置用它——比如镜片支撑环、基准面。
- 微晶玻璃：热膨胀几乎为零，而且硬度高、稳定性好。但它是脆性材料，不能承受拉应力。有限元分析里，微晶玻璃要用最大主应力准则来判断是否失效，而不是 von Mises。
避坑指南：我曾经用微晶玻璃做了一块反射镜，有限元算出来 von Mises 应力只有 20 MPa，我觉得没问题。结果装配时一拧螺钉，镜片裂了。后来才明白——微晶玻璃的抗拉强度只有 10~20 MPa，而抗压强度有 200 MPa。脆性材料，永远不要用 von Mises 应力去判断！
3.6 知识结构图：材料力学基础
下面这张图，帮你把今天的内容串起来。
3.7 小结：你该带走什么？
这一章内容不少，但核心就几条：
- 应力、应变、胡克定律是有限元分析的“语言”，不懂它们就看不懂结果。
- 泊松比虽然小，但影响横向变形和接触分析，别忽略。
- 屈服强度是材料的底线，安全系数是给模型误差留的余量。
- 铝合金、殷钢、微晶玻璃各有各的脾气，选材时一定要看热膨胀和脆性。