3. 网格类型选择与单元阶次决策

在二维有限元建模中，网格类型与单元阶次是两个相互独立又彼此耦合的决策。网格类型决定单元几何形状（三角形或四边形），单元阶次决定插值多项式的完备程度（一阶线性或二阶二次）。这两个选择共同决定了离散误差、计算规模与结果可信度。本节从力学特性出发，给出可操作的决策依据。

三角形与四边形单元的力学差异

三节点常应变三角形（CST）是最基础的单元，其位移场为线性插值，单元内应变为常数。这一特性带来两个后果：一是弯曲变形需要大量单元才能逼近，因为常应变无法表达单元内部的应变梯度；二是存在明显的剪切自锁倾向，在弯曲主导问题中刚度偏大。四节点四边形（Q4）采用双线性插值，应变在单元内可随坐标变化，弯曲适应能力优于CST，但同样存在剪切自锁，且在畸变严重时精度下降明显。

二阶单元显著改善上述问题。六节点二次三角形（T6）和八节点二次四边形（Q8）的位移场为二次完备，能够自然表达线性应变场，对弯曲和应力集中区域的逼近效率远高于一阶单元。代价是每个单元的自由度数增加，带宽和求解时间上升。工程经验是：在应力梯度平缓的区域，一阶四边形已足够；在弯曲主导或存在孔洞、圆角、裂纹尖端的区域，二阶单元能以更少的单元数达到同等精度。

阶次选择的量化依据

阶次并非越高越好。以悬臂梁弯曲为例，若用一阶四边形沿厚度方向仅布置一层单元，剪切自锁会使挠度显著偏小；布置四层以上则逐步收敛。若改用二阶四边形，两层即可获得接近理论解的挠度。这说明阶次选择与网格密度是同一问题的两个侧面：提高阶次等价于在单元内部增加自由度，其效果类似于加密网格，但自由度增长更经济。

判断是否需要用二阶单元，可参考以下准则：当关注量为位移且结构以弯曲为主时，二阶单元的优势明显；当关注量为应力且峰值位于几何不连续处时，二阶单元能更准确地捕捉应力梯度，但需注意节点应力与积分点应力的外推误差。对于接触、大变形等非线性问题，一阶单元因积分点少、计算稳定，常被优先采用，此时通过局部加密补偿精度。

网格类型的几何适应性

四边形单元要求区域可映射为规则拓扑，对于矩形、圆环等规则区域，结构化四边形网格质量高、精度好。三角形单元对任意几何的适应性更强，自动划分成功率高，适合复杂边界。工程实践中常见的策略是：主体规则区域用四边形，过渡区和复杂边界附近用三角形，形成混合网格。需注意三角形与四边形交界处应保证节点协调，避免出现悬挂节点。

单元畸变是比类型更关键的精度影响因素。四边形单元的长宽比宜控制在3以内，内角避免接近0°或180°；三角形单元应避免过小的内角。畸变严重时，即使采用二阶单元，精度也会急剧下降。因此，在完成网格划分后，应检查单元质量指标（如雅可比行列式、长宽比、翘曲度），对超标区域进行局部重划。

决策流程与示例

以一个带圆孔矩形板的平面应力分析为例，说明决策过程。分析目标为获取孔边应力集中系数。孔边存在明显的应力梯度，属于应力集中主导问题，应选用二阶单元；板的主体区域几何规则，可采用四边形；孔边曲边用二阶三角形过渡更为自然。最终方案为：主体区域Q8单元，孔边一圈T6单元，两者节点协调。若计算资源紧张，可退化为Q4主体加孔边局部加密，但需通过网格收敛性验证确认应力峰值误差在可接受范围内。

网格收敛性验证是决策的最终裁判。建议对同一模型采用两套不同密度的网格（如单元尺寸减半），比较关注量的变化。若两次结果的相对偏差小于5%，可认为网格已收敛；若偏差较大，应优先考虑提高阶次或局部加密，而非全局细化。这一验证过程本身也是判断网格类型与阶次是否合理的直接证据。