第二章 有限元法入门：基本原理、网格划分、边界条件与材料本构
各位同行，大家好。这一章我们聊聊有限元法。说实话，很多搞铸造的朋友一听到“有限元”三个字就头疼，觉得那是搞力学的人的事。其实不然。你想想看，我们做应力场模拟，本质上就是在解一个偏微分方程。有限元法，说白了就是把一个连续体切成很多小块，每块单独算，再拼回去。
我在刚入行那会儿，也踩过不少坑。有一次算一个大型铸钢件的热应力，网格画得太粗，结果裂纹位置完全算偏了。嗯，从那以后，我对网格划分就再也不敢马虎了。
2.1 有限元法的基本原理
有限元法的核心思想，我总结成八个字：离散近似，分片插值。什么意思呢？
- 离散：把连续的铸件几何体，划分成有限个单元（比如四面体、六面体）。
- 近似：在每个单元内部，用简单的形函数（通常是多项式）来近似描述位移场或温度场。
- 组装：把所有单元的方程组装成一个全局方程组，[K]{u} = {F}。
- 求解：解这个方程组，得到每个节点的位移，进而算出应变和应力。
我个人习惯把有限元法比作“搭积木”。每个单元就是一块积木，形状简单，但拼在一起就能模拟出复杂的结构。你想想看，是不是这个道理？
核心公式（静力分析）：
[K]{u} = {F}其中 [K] 是整体刚度矩阵，{u} 是节点位移向量，{F} 是节点力向量。对于热应力分析，我们还要加上温度载荷项。
我建议初学者不要一上来就啃理论推导。先理解这个“离散-近似-组装-求解”的流程，后面再慢慢深入。
2.2 网格划分策略
网格划分是有限元分析中最耗时、也最考验经验的一步。我见过太多人，模型建得漂亮，但网格一塌糊涂，结果自然不可信。
2.2.1 单元类型选择
| 单元类型 | 适用场景 | 我的建议 |
|---|---|---|
| 一阶四面体（C3D4） | 复杂几何、快速试算 | 慎用！太刚硬，精度差 |
| 二阶四面体（C3D10） | 复杂几何、中等精度 | 我常用，适合铸造件 |
| 一阶六面体（C3D8） | 简单几何、高精度 | 精度好，但网格难画 |
| 二阶六面体（C3D20） | 关键区域、弯曲问题 | 计算量大，慎用 |
我在项目中遇到过这样的情况：一个大型壳体，用四面体算了一周，结果应力集中区算出来偏小。后来换成六面体，虽然网格花了三天时间，但结果跟实测吻合得很好。所以，网格质量直接决定分析成败。
2.2.2 网格密度控制
网格不是越密越好。太密，计算时间成倍增加；太疏，结果失真。我的经验是：
- 应力集中区（如圆角、台阶、壁厚突变处）：网格加密，至少3层单元。
- 温度梯度大的区域（如铸件与铸型界面）：网格加密。
- 远离热节和应力峰的区域：网格可以适当稀疏。
我的小技巧：先用粗网格算一遍，找到应力热点。然后在热点区域局部加密，再算一遍。这样既省时间，又保证精度。
2.2.3 网格质量检查
网格画完后，一定要检查质量。我常用的指标有：
- 长宽比（Aspect Ratio）：最好小于5，最大不超过10。
- 扭曲度（Skewness）：最好小于0.8，最大不超过0.95。
- 雅可比（Jacobian）：必须大于0，最好大于0.7。
警告：我曾经遇到一个案例，网格雅可比出现负值，结果算出来的应力场完全乱套。检查后发现是网格畸变太严重。所以，负雅可比的单元必须重画，不能凑合。
2.3 边界条件设定
边界条件设定，说白了就是告诉软件“这个零件是怎么被约束的、受什么力”。设定错了，结果就是错的。
2.3.1 位移边界条件
在铸造应力分析中，我们通常需要约束刚体位移。我的做法是：
- 在铸件上选三个不共线的点，施加适当的位移约束。
- 注意不要过约束，否则会产生虚假应力。
举个例子：一个长条形的铸件，我通常会在左端约束X、Y、Z三个方向，右端只约束Y、Z方向，让它可以自由伸缩。这样模拟出来的热应力才真实。
2.3.2 热边界条件
热应力分析离不开温度场。热边界条件包括：
- 初始温度：浇注温度、铸型预热温度。
- 换热系数：铸件与铸型之间、铸件与空气之间。
- 环境温度：车间温度。
我记得有一次，一个同事把换热系数设错了，结果算出来的冷却速度比实际快了一倍。后来我们实测了铸型温度，修正了换热系数，结果才对了。所以，边界条件一定要有实测数据支撑。
2.4 材料本构模型选择
材料本构模型，就是描述材料“受力后怎么变形”的数学关系。选错了模型，应力场模拟就是空中楼阁。
2.4.1 常用本构模型
| 模型名称 | 适用材料 | 特点 |
|---|---|---|
| 线弹性模型 | 脆性材料、低温阶段 | 简单，但忽略塑性 |
| 理想弹塑性模型 | 低碳钢、铝合金 | 屈服后应力不变 |
| 线性硬化弹塑性模型 | 合金钢、铸铁 | 屈服后应力继续增加 |
| 蠕变模型 | 高温下的铸件 | 考虑时间效应 |
2.4.2 我的选择建议
对于铸造应力场模拟，我个人习惯这样选：
- 凝固阶段：用蠕变模型，因为高温下材料会蠕变松弛。
- 冷却阶段（高温段）：用弹塑性模型，考虑屈服和硬化。
- 冷却阶段（低温段）：用线弹性模型，因为材料已经变脆。
重要提醒：材料参数（弹性模量、屈服强度、热膨胀系数等）必须是随温度变化的。我见过有人用室温参数算高温应力，结果完全失真。记住，没有温度相关的材料参数，就不要做热应力分析。
2.5 本章知识体系
下面这张图，是我自己整理的有限元法入门知识框架。你可以把它当作一张“地图”，随时回来对照。
好了，这一章的内容就到这里。有限元法入门，说白了就是理解“离散近似”这个核心思想，然后掌握网格、边界、材料这三个实操要点。下一章我们会深入讨论热-力耦合分析的具体实现。各位回去可以拿一个简单的铸件练练手，先画网格，再设边界，最后跑一遍。遇到问题，随时回来翻翻这一章。
公众号:蓝海资料掘金营，微信deep3321