第2章 有限元基础回顾：单元类型选择、网格划分、边界条件与求解器设置
各位工程师朋友，咱们直接进入正题。有限元分析这东西，说白了就是把一个连续体拆成小块块，然后算每个小块的行为。但怎么拆、用什么块、边界怎么给，这里面的门道可不少。我做了十几年光机结构，见过太多分析结果跟实测对不上的案例，十有八九都是基础没打牢。
2.1 单元类型选择：实体、壳、梁
选单元类型，就像选工具。拧螺丝用扳手，切钢板用角磨机。你非要用扳手去切钢板，那肯定不行。有限元里也一样，选错了单元，结果就是错的。
2.1.1 实体单元（Solid）
实体单元是最通用的。它有三个方向的尺寸，能模拟真实的体积行为。我个人习惯，只要结构三个方向的尺寸相差不超过5倍，优先用实体单元。
- 一阶六面体（C3D8）：计算快，但容易锁死。弯曲问题慎用。
- 二阶六面体（C3D20）：精度高，我常用的主力。代价是计算量大。
- 四面体（C3D4 / C3D10）：自动网格划分的首选。C3D10（二阶）勉强能用，C3D4（一阶）我基本不用，太硬了。
2.1.2 壳单元（Shell）
当一个方向的尺寸远小于另外两个方向（比如厚度/长度 < 1/10），壳单元就登场了。光机结构里的镜筒壁、遮光罩、薄板支撑，都是壳单元的天下。
- S4R（四节点减缩积分壳）：通用性强，我80%的壳分析都用它。
- S8R（八节点壳）：精度更高，适合有应力集中或大变形的情况。
- 注意：壳单元需要定义厚度和偏置。偏置决定了中性面位置，搞错了，连接处的弯矩传递就全乱了。
2.1.3 梁单元（Beam）
梁单元用于细长构件，长度远大于截面尺寸（长细比 > 10）。光机结构里的支撑杆、桁架、调整螺钉，用梁单元最合适。
- B31（两节点梁）：线性，适合简单连接。
- B32（三节点梁）：二次，精度更高。
- 截面定义：必须输入截面形状和尺寸。圆管、矩形、工字钢，程序会自动算惯性矩。
2.2 网格划分原则
网格是有限元的骨架。网格质量直接决定分析成败。你想想看，骨架都歪了，肌肉能长对吗？
2.2.1 网格密度
网格不是越密越好。太密算不动，太疏算不准。我有个经验法则：
| 分析类型 | 网格尺寸建议 | 说明 |
|---|---|---|
| 静力分析（无应力集中） | 特征尺寸的 1/5 ~ 1/10 | 比如镜座边长100mm，网格10~20mm |
| 静力分析（有应力集中） | 局部加密到 1/20 ~ 1/50 | 圆角、尖角处必须加密 |
| 模态分析 | 能捕捉振型即可 | 一般取最高关注频率波长的 1/6 |
| 热-结构耦合 | 兼顾温度梯度 | 温度变化剧烈处加密 |
2.2.2 网格质量指标
网格画完了，得检查质量。我常用的几个指标：
- 长宽比（Aspect Ratio）：最好 < 5，极限 < 10。太长太扁的单元，刚度矩阵会病态。
- 翘曲度（Warpage）：壳单元尤其重要。翘曲 > 15° 就要小心了。
- 雅可比（Jacobian）：反映单元畸变程度。二阶单元雅可比 < 0.7 就算差，< 0.5 必须重画。
- 最小角度：三角形单元最小角 > 20°，四边形 > 45°。
2.2.3 网格过渡
从粗网格到细网格，要平滑过渡。我建议：
- 相邻单元尺寸比 < 1.5
- 避免突然从10mm跳到1mm
- 使用过渡区域，逐步加密
2.3 边界条件与载荷施加
边界条件和载荷，是有限元分析的输入。输入错了，输出肯定错。这个道理很简单，但实际中经常出问题。
2.3.1 边界条件
边界条件模拟的是结构跟外界的连接。光机结构里常见的：
- 固定约束（Encastre）：所有自由度全锁。用于模拟安装基座、刚性支撑。
- 简支（Pinned）：锁平动，放转动。用于模拟铰链、球头连接。
- 弹簧/阻尼：模拟柔性支撑、减振垫。需要输入刚度系数。
- 对称/反对称：利用对称性减少计算量。但前提是结构和载荷都对称。
2.3.2 载荷施加
载荷类型：
- 集中力/力矩：施加在节点上。注意应力集中，实际中力是通过接触面传递的。
- 分布压力：施加在单元面上。模拟气压、水压、接触压力。
- 重力/加速度：施加在整个模型上。注意方向。
- 温度载荷：定义温度场或热膨胀系数。热-结构耦合分析必备。
载荷施加的位置和方式，要跟实际一致。比如，镜片是通过压圈固定的，那载荷应该施加在压圈接触面上，而不是直接加在镜片边缘。
2.4 求解器设置
求解器是有限元的大脑。设置不对，要么算不出来，要么算出来是错的。
2.4.1 求解器类型
| 求解器 | 适用场景 | 特点 |
|---|---|---|
| 直接法（如Abaqus/Standard） | 静力、模态、屈曲 | 精度高，内存消耗大 |
| 显式法（如Abaqus/Explicit） | 冲击、碰撞、大变形 | 适合瞬态，时间步长小 |
| 迭代法（如ANSYS PCG） | 大规模模型 | 内存效率高，收敛性依赖网格质量 |
2.4.2 关键设置参数
- 非线性开关：几何非线性（大变形）、材料非线性（塑性）、接触非线性。光机结构一般小变形，线性分析就够了。但如果有接触或大转动，必须开非线性。
- 时间步长：静力分析用自动步长。显式分析要保证稳定时间步长，跟最小单元尺寸和材料波速有关。
- 收敛准则：力收敛、位移收敛、能量收敛。默认值一般够用，但遇到接触或塑性问题，可以适当放宽。
- 输出控制：只输出关心的结果，别全输出。不然硬盘会爆。
2.5 本章知识体系
下面这张图，是我自己总结的有限元基础决策流程。每次做分析前，我都会过一遍。
这张图把本章的核心逻辑串起来了。从问题定义开始，到单元选择、网格划分、边界载荷、求解器设置，每一步都是环环相扣的。你跳过了哪一步，后面就会在哪一步出问题。