4. 光机结构有限元建模要点
光机结构有限元建模的目标，不是把三维CAD模型原封不动搬进求解器，而是用尽可能少的计算代价，得到对随机振动与冲击响应足够准确的刚度、质量与阻尼分布。几何简化、单元选择、网格划分、连接与接触、材料参数、边界条件这六项，本质上都在回答同一个问题：哪些细节会显著改变结构的模态频率、振型与应力传递路径，哪些细节可以安全地忽略。下面按建模的实际推进顺序展开，并给出可复用的判断依据与示例参数。
一、几何简化：先判断“传力路径”再动手
光机结构常见对象包括镜筒、镜座、柔性支撑、压圈、隔圈、基板与减振器。简化时优先保留主传力路径上的连续实体，删除或等效处理不影响整体刚度的局部特征。典型处理原则如下：
- 可删除：倒角、小圆角（R<1 mm）、螺纹牙型、标识刻字、工艺孔、线缆走线槽。它们对整体模态影响通常小于1%，但会显著增加网格量。
- 可等效：螺纹连接用等效圆柱或刚性单元代替；薄壁翻边用壳单元中面代替；蜂窝夹层板用等效正交各向异性板代替。
- 必须保留：柔性铰链的细颈、减振器的安装接口、镜座与镜筒的配合止口、大开口附近的加强筋。这些位置是刚度突变点，直接决定局部应力与传力路径。
一个实用的判断方法是：先做一次粗网格整体模态分析，记录前10阶频率；再对可疑特征做保留/删除的对比计算。若某特征删除后前10阶频率变化均小于2%，且关注点应力变化小于5%，则可安全简化。示例：某Φ200 mm镜筒，删除全部M4螺纹孔后一阶频率由412 Hz变为414 Hz，变化0.5%，可删除；而删除镜座止口后一阶频率由412 Hz降至356 Hz，变化13.6%，必须保留。
二、单元选择：按“厚度方向应力状态”决定
光机结构以实体和板壳为主，单元选择的核心依据是结构厚度方向是否存在明显应力梯度。常用选择如下：
| 结构对象 | 推荐单元 | 理由 |
|---|---|---|
| 镜筒、基板、镜座实体 | 二阶四面体或六面体 | 厚度方向有弯曲应力梯度，需三维实体 |
| 薄壁遮光罩、翻边 | 壳单元（S4/S3） | 厚度远小于面内尺寸，弯曲为主 |
| 柔性铰链细颈 | 六面体或高阶四面体 | 应力集中，需高精度 |
| 减振器、隔振垫 | 弹簧-阻尼单元或超单元 | 只关心接口阻抗，不关心内部细节 |
| 螺钉连接 | 梁单元或刚性单元 | 简化传力，避免局部网格畸变 |
四面体单元对复杂几何适应性强，但一阶四面体刚度过硬，必须用二阶（C3D10）或六面体。若结构以弯曲为主且厚度方向只有2~3层单元，六面体精度明显优于四面体。示例：某镜座用一阶四面体时一阶频率为520 Hz，改用二阶四面体后降至438 Hz，与试验值445 Hz更接近，说明一阶四面体高估了刚度。
三、网格划分：控制“每波长单元数”与“厚度层数”
随机振动分析关注到较高频率（通常到500~2000 Hz），网格密度必须与最高关注频率对应的弯曲波长相匹配。经验规则：
- 弯曲波每波长至少6~8个二阶单元，或4~6个线性单元。
- 厚度方向至少2层单元（实体），若关注弯曲应力则建议3层。
- 应力集中区（柔性铰链、止口根部）局部加密，过渡区网格尺寸变化不超过2倍。
以铝合金镜筒为例，弹性模量70 GPa，密度2700 kg/m³，板厚3 mm，最高关注频率1500 Hz。弯曲波速约 c ≈ 1500 m/s量级，对应波长约1 m量级，实际按经验取单元尺寸5~8 mm即可。若镜筒直径200 mm，圆周方向约80~120个单元，轴向按5 mm划分。网格质量检查重点：雅可比>0.6，翘曲角<15°，长宽比<5（应力区<3）。
四、连接与接触：区分“线性连接”与“非线性接触”
光机结构中的连接方式直接决定阻尼与刚度。随机振动分析通常先做线性模态，因此连接应尽量线性化：
- 螺栓连接：用刚性单元（RBE2）或梁单元模拟，耦合被连接件对应节点。若需考虑连接刚度，可用CBUSH单元，示例轴向刚度1×10⁶ N/mm，剪切刚度5×10⁵ N/mm。
- 胶接：用薄层实体或弹簧单元，胶层厚度0.05~0.2 mm，弹性模量按胶种取1~10 MPa。
- 接触：仅当存在明显分离、滑移或间隙时才用非线性接触。随机振动若必须用接触，需先做预紧模态，再在预紧状态下做线性化。
五、材料参数：随机振动最敏感的是阻尼
光机结构常用材料参数如下（示例值，实际以材料证书为准）：
| 材料 | 弹性模量 E (GPa) | 密度 ρ (kg/m³) | 泊松比 ν | 结构阻尼比 ζ |
|---|---|---|---|---|
| 铝合金 6061 | 69 | 2700 | 0.33 | 0.01~0.03 |
| 钛合金 TC4 | 110 | 4430 | 0.34 | 0.01~0.02 |
| 铟钢 | 140 | 8100 | 0.30 | 0.005~0.01 |
| 碳纤维复合材料 | 各向异性 | 1600 | 0.30 | 0.02~0.05 |
阻尼比是随机振动响应幅值最敏感的参数。若缺少实测值，金属结构可取0.01~0.02，胶接结构取0.03~0.05，减振器按厂家提供的损耗因子。示例：某镜筒一阶频率420 Hz，阻尼比从0.01提高到0.03，共振点加速度响应均方根值下降约40%。因此阻尼参数必须谨慎取值，并在报告中注明来源。
六、边界条件：约束“真实安装界面”而非“方便位置”
边界条件应模拟结构在整机中的实际安装状态。常见做法：
- 固定约束：仅用于安装面确实与刚性基础相连的情况，约束安装孔周围节点的全部平动自由度。
- 弹性约束：用弹簧单元模拟减振器或柔性安装，示例三向刚度分别为X:200 N/mm、Y:200 N/mm、Z:500 N/mm。
- 自由-自由：用于模态试验对比，前6阶为刚体模态，频率应接近0。
随机振动输入通常以基础加速度PSD给出，方向按发射或运输环境谱。示例：某星载光机随机振动条件为10~2000 Hz，总均方根加速度6.0 g，施加于安装面三个方向。边界条件一旦确定，全文与图示必须保持一致，不得中途切换参照系。
六项要点并非孤立，而是相互制约：几何简化过度会导致刚度失真，网格过密会放大计算量，连接线性化不当会丢失阻尼，边界条件错误则全盘皆输。实际建模时建议按“几何简化→单元选择→网格划分→连接与接触→材料参数→边界条件”的顺序推进，每完成一步做一次模态自检，记录一阶频率与振型，与上一版对比。若一阶频率变化超过5%，应回溯上一步检查。最终模型应能通过自由-自由模态试验或整机模态试验验证，前若干阶频率误差控制在5%~10%以内，方可进入随机振动与冲击响应分析。