有限要素法・有限差分法・有限体積法 — それぞれの原理・適用分野・精度・計算コストをエンジニア視点で解説
CAE・数値解析の分野ではFEM(有限要素法)・FDM(有限差分法)・FVM(有限体積法)の3手法が広く使われています。それぞれ適した問題領域が異なり、使い分けが重要です。本記事では原理・特徴・代表ソフトウェアを比較します。
解析領域を有限個の要素(element)に分割し、各要素内で変位・温度などの物理量を形状関数で近似する手法。弱形式(重み付き残差法)に基づき、Galerkin法・変分法で定式化されます。
強み:複雑形状に対応可能、境界条件の適用が容易、固体力学との相性が良い。
代表ソフト:ANSYS Mechanical、Abaqus、Nastran、LS-DYNA
偏微分方程式の微分演算子を、格子点上の差分商で置き換える手法。テイラー展開による局所的近似が基礎です。実装がシンプルで教育・研究用途に多く使われます。
強み:実装が容易、規則格子との相性が良い、計算コスト低。
弱点:複雑形状への対応が難しい(境界フィッティングが必要)。
代表ソフト:FDTD(電磁界)、独自コード、Matlab PDE Toolbox
解析領域を制御体積(コントロールボリューム)に分割し、各体積内で保存則の積分形を離散化する手法。質量・運動量・エネルギーの保存が離散レベルで厳密に成立します。
強み:保存則の厳密な満足、非構造格子対応、流体解析に最適。
代表ソフト:OpenFOAM、FLUENT(Ansys)、STAR-CCM+、Code_Saturne
| 項目 | FEM | FDM | FVM |
|---|---|---|---|
| 離散化の基礎 | 弱形式・変分原理 | 差分近似(強形式) | 積分保存則 |
| 格子タイプ | 非構造・構造両対応 | 主に構造(規則)格子 | 非構造・構造両対応 |
| 複雑形状 | ◎ 得意 | △ 難しい | ◎ 得意 |
| 保存則の保証 | △ 要素レベルで近似 | △ 大局的には成立 | ◎ 離散的に厳密 |
| 主な適用 | 構造・固体力学・振動 | 熱拡散・波動・電磁 | 流体・熱流体・燃焼 |
| 代表ソフト | ANSYS, Abaqus, Nastran | FDTD系、研究コード | OpenFOAM, FLUENT |
| 実装難易度 | 中〜高 | 低〜中 | 中 |
梁・シェル・ソリッド要素を用いた応力・変位解析、モーダル解析、座屈解析にはFEMが最適です。複雑形状の自動メッシュ生成ツールも充実しています。
OpenFOAMを筆頭に、主要なCFDコードはほぼすべてFVMベースです。圧力・速度の連成(SIMPLE法など)も保存則と相性が良く、乱流モデルも豊富です。
規則格子で良い場合(2D熱伝導、FDTD電磁波解析など)はFDMが最もシンプルに実装でき、教育目的にも適しています。
熱応力連成にはFEM(ANSYS Mechanical + Fluent連成など)、熱流体単独ならFVMが主流です。
精度はいずれもメッシュ依存性が高く、十分に細かいメッシュであれば3手法とも同等の精度を達成できます。計算コストはFDMが一般に最も低く、FEMは行列アセンブリ・ソルバコストが高い傾向があります。FVMは陽解法(非定常)で効率的、陰解法(定常)ではFEMと同程度です。
目的・形状・精度要求に合わせて手法を選択することが、効率的なCAE解析の第一歩です。
厚さ10mmのアルミニウム板(熱伝導率237 W/m·K)に片側100℃、他側0℃の定常熱伝導解析を実施。FEMは2次要素4,096個で中央部の温度50℃(誤差0.2%)、計算時間2.3秒。FDMは20×20格子点で49℃(誤差2.0%)、計算時間0.4秒。FVMは同格子で50.5℃(誤差1.0%)、計算時間0.8秒。複雑な曲面や不規則領域ではFEMが精度向上、単純矩形領域ではFDMが高速。