4、静力学分析基础：有限元法（FEM）入门、网格划分技巧、边界条件设定

各位，欢迎来到实战解析的第四讲。前面我们聊了反射镜的材料、受力，但光有理论公式还不够。你设计出来的支撑结构，到底能不能扛得住？变形量是 10 纳米还是 100 纳米？

这时候，就得请出我们的老朋友——有限元法（FEM）。说白了，就是把一个连续的大镜子，切成成千上万个微小的单元，然后让计算机去解每个小单元的受力，最后拼回整体结果。我当年刚入行时，觉得这东西就是个黑盒子，点一下“求解”就完事了。结果呢？算出来的变形和实测差了 3 倍。从那以后，我再也不敢轻视前处理了。

4.1 有限元法（FEM）入门：从“切西瓜”说起

有限元法的思想其实很朴素。你想想看，一个完整的西瓜，你很难精确描述它每一处的应力。但如果你把它切成很多小块（单元），每个小块内部的变形规律就简单多了——可以用一个多项式近似。把所有小块的方程联立起来，就得到了整个西瓜的变形场。

在反射镜分析中，我们通常关注的是静力学问题。也就是载荷恒定，结构达到平衡状态。FEM 求解的核心方程是：

[K]{u} = {F}

其中 [K] 是刚度矩阵，{u} 是节点位移向量，{F} 是节点力向量。计算机的工作，就是组装这个巨大的矩阵，然后求解这个线性方程组。

概念	通俗理解	在反射镜分析中的角色
单元	西瓜块	镜体被划分成的小块（如四面体、六面体）
节点	西瓜块的角点	单元之间的连接点，位移和力的传递点
自由度	每个角点能动的方向	每个节点有 3 个平移自由度（UX, UY, UZ）
刚度矩阵	西瓜块的“硬骨头”	描述单元抵抗变形的能力，由材料属性和几何决定

4.2 网格划分技巧：好网格是成功的一半

网格划分，是 FEM 分析中最耗时、也最考验经验的环节。我见过太多人，模型一导入，直接点“自动划分”，然后就去喝咖啡了。结果回来一看，应力集中区域网格太粗，变形值完全不对。

网格划分的核心原则就一句话：在应力梯度大的地方加密，在梯度小的地方放粗。

4.2.1 网格类型选择

对于反射镜支撑结构，我们主要面对两种几何：

• 镜体（实体）：推荐使用 四面体 或 六面体 单元。六面体精度更高，但划分复杂；四面体适应性好，适合复杂曲面。
• 支撑结构（如柔性铰链、筋板）：如果厚度远小于其他尺寸，可以用 壳单元 或 梁单元 来简化，大幅降低计算量。

4.2.2 网格质量检查

网格画完了，别急着求解。先做质量检查。我常用的几个指标：

指标	理想值	警告值	我的建议
长宽比	< 3	> 10	超过 10 的单元，刚度矩阵会变“病态”，结果不可信
扭曲度	< 60°	> 90°	扭曲度太大，单元会“翻过去”，导致计算不收敛
雅可比	> 0.7	< 0.3	雅可比为负，单元就“自相交”了，必须重画

4.2.3 网格收敛性验证

这是很多新手容易忽略的一步。你画了一套网格，算出一个变形值。但你怎么知道这个值准不准？

方法很简单：逐步加密网格。比如先画 5 万单元，算一次；再画 10 万单元，算一次；再画 20 万单元……直到相邻两次的结果差异小于 5%。这时候，我们就说结果“收敛”了。

# 网格收敛性验证流程（伪代码）
1. 初始网格：50,000 单元 → 计算变形 δ1
2. 加密网格：100,000 单元 → 计算变形 δ2
3. 检查：|δ2 - δ1| / δ1 < 5% ?
   - 是：结果收敛，使用 δ2
   - 否：继续加密到 200,000 单元，重复步骤 2-3

4.3 边界条件设定：别让模型“飞”起来

边界条件，就是告诉计算机：你的模型是怎么固定、怎么受力的。设定错了，结果就是“空中楼阁”。

4.3.1 约束条件

反射镜支撑结构的约束，通常来自与镜室或基座的连接面。常见的约束类型：

• 固定约束：限制所有自由度（UX=UY=UZ=0）。用于模拟螺栓连接或粘接面。
• 位移约束：只限制某个方向的位移。比如镜体底面只能沿 Z 向移动，不能沿 X/Y 向移动。
• 对称约束：如果结构对称，可以只建一半模型，在对称面上施加对称约束。能省一半计算量。

4.3.2 载荷条件

反射镜承受的载荷主要有：

• 重力：最基础也最容易被忽略的载荷。记得设置正确的重力方向（通常是 -Z 向）。
• 温度载荷：如果镜体与支撑结构材料不同（比如 SiC 镜体 + 钛合金支撑），温度变化会产生热应力。
• 装配预紧力：螺栓拧紧产生的力，会影响镜面面形。

4.3.3 接触与绑定

如果模型包含多个零件（比如镜体 + 支撑座），就需要定义零件之间的接触关系。

• 绑定接触：两个面粘在一起，没有相对滑动。适用于胶接或过盈配合。
• 摩擦接触：允许滑动，但有摩擦力。适用于螺栓连接面。
• 无摩擦接触：允许滑动，无摩擦力。适用于理想润滑面。

4.4 本章知识体系总览

为了让你更直观地理解本章的知识结构，我画了一张流程图。它展示了从几何模型到结果验证的完整路径。

这张图清晰地展示了分析流程：从几何模型出发，赋予材料属性，然后进行网格划分（同时做质量检查和收敛性验证），接着设定边界条件，最后求解并查看结果。每一步都环环相扣，缺一不可。