第4章 三维四面体网格划分：从Tetra10到质量优化
四面体网格，说白了就是光机结构仿真里最常用的网格类型之一。尤其是Tetra10节点单元，我做了这么多年光机仿真，几乎每个项目都离不开它。今天咱们就把它掰开揉碎了讲清楚。
4.1 Tetra10节点单元详解
先说说什么是Tetra10。普通的四面体单元有4个节点，每个节点3个自由度，算起来快但精度一般。Tetra10呢？它在每条边的中点多加了一个节点，一共10个节点。
核心区别：Tetra4是线性单元，Tetra10是二次单元。二次单元能更好地模拟弯曲变形，这对光机结构来说太重要了——镜面变形哪怕差几个纳米，成像质量就完全不一样。
我个人的习惯是：做光机结构分析，至少用Tetra10起步。Tetra4我只有在做概念设计、快速迭代时才用。为什么？你看这个对比就明白了：
| 特性 | Tetra4（4节点） | Tetra10（10节点） |
|---|---|---|
| 形函数阶次 | 线性 | 二次 |
| 边中节点 | 无 | 有 |
| 弯曲精度 | 较差，需加密 | 较好，效率高 |
| 计算量 | 低 | 中等 |
| 光机适用场景 | 初步估算 | 详细分析、优化 |
我记得有一次做空间相机的主镜支撑分析，甲方要求变形精度到纳米级。我一开始偷懒用了Tetra4，结果算出来的镜面变形云图全是锯齿状的，根本没法用。换成Tetra10之后，结果平滑多了，和实验数据对得上。
小技巧：在Abaqus里，Tetra10对应C3D10单元；在ANSYS里是SOLID187。设置时注意检查单元阶次，别选成线性单元了。
4.2 基于曲率的自适应网格
做网格划分最头疼的是什么？就是不知道哪里该加密、哪里可以粗一点。基于曲率的自适应网格，就是让软件帮你判断。
原理其实很简单：曲率大的地方（比如圆角、尖角、薄壁），网格自动加密；曲率小的地方（比如平面、大曲面），网格自动变粗。说白了，就是让网格密度跟着几何特征走。
我一般在HyperMesh或ANSA里设置这几个参数：
- 最小曲率半径：设定一个阈值，小于这个值的区域加密
- 曲率法向角：每段弧上允许的最大角度变化，通常设15°~30°
- 全局种子尺寸：作为基准尺寸，曲率加密在此基础上细化
注意：自适应网格不是万能的。我曾经遇到过一个案例，镜筒内壁有微小的螺纹特征，自适应算法把整个内壁都加密了，网格量暴增3倍。后来我手动在螺纹区域做了局部控制，才把网格量降下来。
为什么会这样？因为自适应算法只认几何曲率，不认工程判断。你想想看，螺纹虽然曲率大，但如果你不关心螺纹本身的应力，完全可以粗化处理。所以我的建议是：自适应网格做初稿，人工调整做精修。
4.3 局部尺寸控制与增长比
局部尺寸控制，就是你在关键区域手动指定网格大小。增长比，则是控制网格从密到疏的过渡速度。
这两个参数是绑在一起的。增长比设得太小，过渡太慢，网格量爆炸；增长比设得太大，过渡太陡，单元质量下降。
我个人的经验值：
| 场景 | 推荐增长比 | 说明 |
|---|---|---|
| 一般光机结构 | 1.3 ~ 1.5 | 平衡质量和效率 |
| 高精度镜面区域 | 1.2 ~ 1.3 | 保证变形精度 |
| 非关键支撑区域 | 1.5 ~ 2.0 | 快速过渡，减少单元 |
| 薄壁/尖角过渡 | 1.1 ~ 1.2 | 防止单元畸变 |
举个例子。做镜面网格时，我会在镜面表面设一个0.5mm的局部尺寸，然后以1.3的增长比向外过渡到镜体内部。这样镜面附近有足够密的网格来捕捉变形，而镜体内部可以用较粗的网格来节省计算资源。
避坑指南：我曾经在镜面边缘设了0.2mm的局部尺寸，增长比用了2.0，结果从镜面到镜体只过渡了3层网格，单元长宽比直接飙到20以上，算出来的应力全是假的。后来我把增长比降到1.3，过渡了8层，结果就对了。
4.4 四面体网格质量优化
网格质量不好，算出来的结果你敢信吗？反正我不敢。四面体网格的质量指标，我重点关注这几个：
- 雅可比（Jacobian）：反映单元畸变程度，一般要求>0.7，至少>0.5
- 长宽比（Aspect Ratio）：单元最长边与最短边之比，一般<5
- 扭曲度（Skewness）：偏离正四面体的程度，一般<0.8
- 最小角/最大角：三角形面的角度，20°~120°为佳
优化方法我总结了三板斧：
- 重新划分：调整种子尺寸或增长比，重新生成网格
- 局部重划：只对质量差的区域重新划分，保留好区域
- 网格平滑：用Laplacian平滑或优化平滑，调整节点位置
我常用的流程是这样的：先全局划分，然后检查质量，把质量差的单元标记出来。如果坏单元集中在某个区域，就局部调整那里的种子尺寸；如果坏单元分散，就做一两轮平滑。
警告：平滑不是万能的。过度平滑会把单元拉变形，反而降低质量。我一般最多做3轮平滑，每轮检查一次质量变化。
下面这张图总结了本章的核心逻辑：
嗯，这张图把咱们今天讲的内容串起来了。从Tetra10单元出发，经过曲率自适应、局部尺寸和增长比的控制，最后落到质量优化上。每一步都环环相扣。
最后说一句：网格划分没有银弹。每个项目都有自己的特点，多试几次，多检查质量，慢慢就有感觉了。