有限元法·有限差分法·有限体积法 — 从工程师视角解说各方法的原理·适用领域·精度·计算成本
在CAE·数值分析领域，FEM（有限元法）·FDM（有限差分法）·FVM（有限体积法）这三种方法被广泛应用。各方法适合的问题领域不同，如何合理选用至关重要。本文对三种方法的原理·特点·代表软件进行对比介绍。
将分析域划分为有限个单元（element），在每个单元内用形函数对位移·温度等物理量进行近似的方法。基于弱形式（加权残差法），通过Galerkin法·变分法进行离散化。
优势：可处理复杂形状，边界条件施加方便，与固体力学契合度高。
代表软件：ANSYS Mechanical、Abaqus、Nastran、LS-DYNA
将偏微分方程中的微分算子用格点上的差商代替的方法。以泰勒展开为基础进行局部近似。实现简单，常用于教学·研究用途。
优势：实现容易，与规则格点配合良好，计算成本低。
弱点：难以处理复杂形状（需要边界拟合）。
代表软件：FDTD（电磁场）、自研代码、Matlab PDE Toolbox
将分析域划分为控制体积（Control Volume），在每个体积内对守恒律积分形式进行离散化的方法。质量·动量·能量守恒在离散层面严格成立。
优势：严格满足守恒律，支持非结构网格，最适合流体分析。
代表软件：OpenFOAM、FLUENT（Ansys）、STAR-CCM+、Code_Saturne
| 项目 | FEM | FDM | FVM |
|---|---|---|---|
| 离散化基础 | 弱形式·变分原理 | 差分近似（强形式） | 积分守恒律 |
| 网格类型 | 非结构·结构均支持 | 主要为结构（规则）网格 | 非结构·结构均支持 |
| 复杂形状 | ◎ 擅长 | △ 困难 | ◎ 擅长 |
| 守恒律保证 | △ 单元级近似 | △ 整体上成立 | ◎ 离散层面严格成立 |
| 主要应用 | 结构·固体力学·振动 | 热扩散·波动·电磁 | 流体·热流体·燃烧 |
| 代表软件 | ANSYS, Abaqus, Nastran | FDTD系、研究代码 | OpenFOAM, FLUENT |
| 实现难度 | 中～高 | 低～中 | 中 |
使用梁·壳·实体单元进行应力·位移分析、模态分析、屈曲分析时，FEM是最佳选择。复杂形状的自动网格生成工具也十分完善。
以OpenFOAM为代表，主要CFD代码几乎全部基于FVM。压力·速度耦合（SIMPLE法等）与守恒律契合度高，湍流模型也十分丰富。
适用规则网格的情况（二维热传导、FDTD电磁波分析等），FDM实现最为简单，也适合教学目的。
热应力耦合推荐FEM（如ANSYS Mechanical + Fluent联合仿真），纯热流体分析则以FVM为主流。
三种方法的精度均与网格依赖性密切相关，只要网格足够细密，三种方法均可达到相近的精度。计算成本方面，FDM通常最低，FEM因矩阵装配·求解器成本较高。FVM在显式（非稳态）方法下效率较高，隐式（稳态）方法下与FEM相当。
根据目的·形状·精度要求选择合适的方法，是高效CAE分析的第一步。
钢铁厂热轧钢板冷却过程：板厚100mm，初温900°C，冷却时间60s。FVM在CFX中划分20万四面体网格，迭代步长0.1s，计算温度场分布耗时约3小时；同工况用FDM在100×100×10网格上仅需12分钟但难以准确捕捉边界急冷的温度梯度（误差>8%）。FEM用Abaqus耦合热-应力分析，追踪残余应力分布，需40万单元，计算时间8小时，精度±2°C。