4. 二维问题网格划分：平面应力/应变问题的三角形与四边形网格划分、网格密度控制

各位同学，咱们今天聊聊二维问题的网格划分。说实话，这是有限元分析里最基础也最容易被忽视的一环。我见过太多人一上来就闷头算，结果网格画得乱七八糟，算出来的结果自己都不敢信。

平面应力和平面应变问题，说白了就是二维简化。你想想看，实际结构都是三维的，但有些情况我们可以偷个懒——比如薄板受面内力，或者厚坝体受横向力。这时候用二维模型，效率能翻好几倍。

4.1 三角形网格 vs 四边形网格

先说说两种最常用的单元。三角形单元和四边形单元，各有各的脾气。

三角形单元：

• 适应性极强。再复杂的几何边界，它都能贴上去。我早期做汽车零件分析时，有些铸件形状特别怪，三角形网格一上去就服服帖帖。
• 但精度嘛...说实话，线性三角形单元（3节点）的精度很一般。它内部是常应变场，说白了就是一个单元里应变不变。应力梯度大的地方，你得画很密才行。
• 二次三角形单元（6节点）会好很多，但计算量也上去了。

四边形单元：

• 我个人习惯，能用四边形尽量用四边形。为什么？精度高啊！同样节点数，四边形的计算精度比三角形高一个档次。
• 但四边形有个毛病——怕畸变。你把它拉得太歪，或者内角太离谱，结果就会飘。我记得有次做悬臂梁分析，网格稍微歪了点，应力结果差了30%多。
• 对于复杂边界，四边形网格生成难度大。有时候你得手动切分几何区域。

4.2 网格密度控制——避坑指南

网格密度怎么定？这是门学问。太粗算不准，太细算不动。我曾经吃过这个亏——有次做压力容器分析，网格画得太粗，应力集中区的结果完全不对，差点导致设计变更。

基本原则：

1. 应力梯度大的地方加密。比如孔边、缺口、拐角、接触区域。这些地方应力变化剧烈，粗网格根本捕捉不到。
2. 应力平缓的地方可以放粗。远离边界、载荷的区域，网格可以大一些。没必要均匀加密，那是浪费计算资源。
3. 过渡要平滑。从密网格到疏网格，尺寸变化率不要超过1.5倍。突然跳变会导致数值误差。

具体怎么控制？ 我常用的方法：

• 先画粗网格算一遍，看应力分布。找到高应力区。
• 在高应力区局部加密，再算。对比两次结果。
• 如果应力变化超过5%，继续加密。直到变化小于5%为止。

这就是所谓的网格收敛性分析。说白了，就是不断加密网格，看结果还变不变。不变了，就说明网格够密了。

4.3 平面应力与平面应变的区别

这里有个容易混淆的点。平面应力和平面应变，虽然都是二维，但物理本质不同。

类型	适用场景	厚度方向	典型例子
平面应力	薄板（厚度远小于其他尺寸）	应力为零，有应变	薄板、壳体、膜结构
平面应变	长结构（长度远大于截面尺寸）	应变为零，有应力	坝体、隧道、长轴

网格划分时，两者没有本质区别。但要注意：平面应变问题通常刚度更大，应力集中更明显。所以网格密度要求往往更高。

4.4 知识体系框架

下面这张图，是我自己总结的二维网格划分决策流程。你照着走，基本不会出错。

4.5 实战要点总结

最后，我把这些年踩过的坑总结成几条：

• 网格质量比数量重要。一个畸变的四边形单元，不如两个规整的三角形单元。
• 不要迷信自动网格划分。软件自动生成的网格，往往在关键区域不够密。你得手动干预。
• 收敛性分析必须做。不做收敛性分析的结果，我不敢用。这是底线。
• 平面应变问题网格要更密。因为厚度方向有约束，应力集中更严重。

嗯，二维网格划分就讲这么多。你回去找几个简单例子练练手，比如带圆孔的平板拉伸。先画粗网格算，再加密，看看应力集中系数怎么变。实践出真知。

---

公众号:蓝海资料掘金营，微信deep3321