第3章 二维网格划分策略：壳单元选择与网格控制

二维网格划分，说白了就是把薄壁结构用壳单元来模拟。光机结构里，反射镜支撑板、镜筒壁、遮光罩，这些厚度远小于长宽尺寸的零件，都适合用壳单元。我做了这么多年光机仿真，发现很多新手在二维网格上栽跟头——不是单元选错，就是网格质量太差，算出来的变形结果根本没法用。

这一章，咱们就聊聊二维网格划分的核心策略。我会把壳单元选择、四边形主导网格、三角形过渡技巧、网格偏置与局部加密这几个关键点，结合我实际项目中的经验，给你讲透。

3.1 壳单元选择：SHELL181 vs SHELL63

ANSYS里最常用的壳单元就两个：SHELL63和SHELL181。很多人问我，到底该用哪个？我的回答很直接：能用SHELL181就别用SHELL63。

为什么？咱们看个对比：

对比项	SHELL63	SHELL181
单元类型	弹性壳（4节点）	有限应变壳（4节点）
厚度变化	不支持	支持渐变厚度
复合材料	不支持	支持多层铺层
大变形	小变形假设	支持大变形、大应变
应力刚化	需额外设置	默认包含
输出精度	线性位移场	二次位移场（更准）

我记得有一次，一个同事用SHELL63算镜筒的变形，结果跟实验差了30%。我帮他换成SHELL181，重新划分网格，误差直接降到5%以内。你想想看，光机系统里微米级的变形都会影响成像质量，30%的误差谁能接受？

3.2 四边形主导网格算法

二维网格划分，首选四边形网格。为什么？因为四边形网格在相同节点数下，计算精度比三角形高，而且收敛性更好。我个人的习惯是，能用四边形的地方绝不用三角形。

四边形主导网格的生成策略，我总结为三步：

1. 几何清理：去掉小倒角、小孔、短边，避免生成畸形网格
2. 分区映射：把复杂几何分割成多个可映射的四边形区域
3. 网格生成：用Mapped或Sweep方法生成规则四边形

举个例子，一个圆形反射镜的支撑板，如果直接自由划分，会生成大量三角形。正确的做法是：

! ANSYS APDL 示例：圆形板四边形网格
/prep7
et,1,181        ! 定义SHELL181单元
rectng,0,100,0,100  ! 先画一个正方形
cyl4,50,50,50   ! 在正方形内画圆
asba,1,2        ! 布尔减，得到圆环区域
lesize,all,5    ! 设置单元尺寸5mm
mshape,0,2d     ! 指定四边形单元
mshkey,1        ! 映射网格
amesh,all       ! 划分网格

3.3 三角形网格过渡技巧

说实话，完全避免三角形网格是不现实的。复杂几何里，总有一些区域四边形网格生成不了。这时候，三角形网格的过渡技巧就派上用场了。

三角形网格的核心问题是什么？刚度偏硬。同样的载荷下，三角形单元比四边形单元变形小，这会导致局部应力失真。所以，三角形网格只能用在非关键区域，而且必须控制数量。

我常用的过渡策略：

• 金字塔过渡：在四边形区域和三角形区域之间，用一层金字塔单元（五面体）过渡
• 密度渐变：从四边形区域到三角形区域，网格尺寸逐渐变化，避免突变
• 局部替换：如果三角形区域出现在关键位置，我会手动调整几何，重新分区

3.4 网格偏置与局部加密

网格偏置，说白了就是让网格在某个方向上疏密变化。这在光机结构里非常实用。比如，镜筒的根部应力集中，需要加密；远离根部的区域应力平缓，可以稀疏。

我常用的偏置方法：

偏置类型	适用场景	偏置比建议
线性偏置	悬臂梁、支撑臂	1:5 ~ 1:10
指数偏置	应力集中区域	1:10 ~ 1:20
双向往复偏置	对称结构	1:3 ~ 1:5

局部加密就更直接了——在关键区域把网格尺寸缩小。比如，螺栓孔周围、支撑点、焊缝附近，这些地方应力梯度大，必须加密。

3.5 知识体系总览

为了让你更直观地理解本章的知识结构，我画了一张流程图：

3.6 实战要点总结

最后，我把本章的核心要点整理成一张速查表，方便你以后查阅：

场景	推荐做法	注意事项
薄壁支撑板	SHELL181 + 四边形映射网格	厚度方向至少2层单元
镜筒壁	SHELL181 + 四边形主导 + 局部加密	法兰连接处加密到1mm
复杂曲面	四边形为主 + 三角形过渡	三角形比例<5%
应力集中区	指数偏置 + 局部加密	偏置比不超过1:20
初步估算	SHELL63 + 粗四边形网格	仅用于趋势判断