有限要素法とは、連続体力学・熱工学・電磁気学などの物理現象を支配する偏微分方程式を数値的に解くための計算手法であり、現代の工学設計・解析において最も広く使用されている数値解析技術の一つです。
航空機・自動車・橋梁・建築構造物・医療機器など、あらゆる分野の設計・強度解析・熱解析・振動解析に有限要素法が活用されており、CAE（Computer Aided Engineering）の中核技術として欠かせない存在となっています。
この記事では、有限要素法の基本概念・支配方程式の離散化・メッシュ分割・要素の種類・解析の流れについて、数学的背景とともにわかりやすく解説していきます。
有限要素法を初めて学ぶ方にも、基礎を改めて整理したい方にも、理解が深まる内容を目指していますのでぜひ最後までお読みください。
有限要素法とは何か？基本概念と数値解析における位置づけ
それではまず、有限要素法の基本的な定義と数値解析手法の中での位置づけについて解説していきます。
有限要素法（FEM：Finite Element Method）は、解析対象の連続体を有限個の小さな要素（エレメント）に分割し、各要素内で物理量を近似関数（形状関数）で表現することで、全体の偏微分方程式を連立代数方程式に変換して数値的に解く手法です。
解析対象が単純な形状・均質な材料であれば解析解（厳密解）が存在する場合がありますが、複雑な形状・非均質材料・複雑な境界条件を持つ実際の工学問題では解析解を求めることが極めて困難です。
有限要素法はこのような複雑な問題に対して近似的な数値解を効率よく求めることができる強力なツールです。
数値解析手法には有限要素法以外に、有限差分法（FDM）・有限体積法（FVM）・境界要素法（BEM）などがありますが、複雑形状への対応力・様々な物理現象への汎用性という点で有限要素法が最も広く採用されています。
有限要素法の歴史と発展
有限要素法の起源は1950年代のターナー・クラフ・マーティン・トップらによる航空機構造解析への応用に求められます。
1960年にクラフ（Clough）が「有限要素法」という名称を初めて使用し、その後急速に理論整備と適用分野の拡大が進みました。
1970〜80年代にコンピュータの性能向上とともにFEM専用ソフトウェア（NASTRAN・ANSYS・ABAQUS等）が開発・普及し、工業界での実用化が加速しました。
現在ではオープンソースのFEMソルバー（OpenFOAM・Code_Aster・FEniCSなど）も利用可能で、研究・教育・産業の幅広い場面で活用されています。
有限要素法は現代工学設計において「バーチャル試験台」として機能しており、実物試作前の設計最適化・破損予測・性能評価に欠かせない技術です。
有限要素法の適用分野
有限要素法が適用される主な物理分野は非常に広範です。
| 分野 | 解析の種類 | 代表的な応用 |
|---|---|---|
| 構造力学 | 静的・動的・座屈・疲労解析 | 橋梁・建築・航空機・自動車の強度設計 |
| 熱工学 | 定常・非定常熱伝導解析 | エンジン冷却・電子機器熱設計 |
| 流体力学 | 層流・乱流・多相流解析 | 翼型設計・管内流動・CFD |
| 電磁気学 | 静電・静磁・電波解析 | モーター・トランス設計・アンテナ設計 |
| 地盤工学 | 土圧・浸透流・液状化解析 | トンネル・基礎構造解析 |
| 生体力学 | 骨・軟組織・血流解析 | 人工関節・医療機器設計 |
偏微分方程式と弱形式（ガラーキン法）
有限要素法の数学的基盤は偏微分方程式（PDE）の弱形式（変分形式）への変換にあります。
支配偏微分方程式（強形式）に重み関数を乗じて積分することで、より弱い条件（弱形式）に変換します。
ガラーキン法では重み関数と試験関数（形状関数）を同一に選ぶことで、系統的な離散化が可能になります。
この弱形式への変換が、不連続な材料特性・複雑な境界条件を持つ問題の取り扱いを可能にする有限要素法の数学的強みです。
有限要素法の基本手順：メッシュ・形状関数・剛性マトリックス
続いては、有限要素法の解析を進めるための基本的な手順について確認していきます。
有限要素解析は前処理・求解・後処理の三段階から構成されます。
メッシュ分割（離散化）
有限要素解析の最初の重要ステップは解析対象をメッシュ（要素の集合体）に分割することです。
メッシュの品質（要素の形状比・サイズ・分布）が解析精度と計算コストに直接影響します。
要素が細かいほど（細分化するほど）解析精度は向上しますが、計算時間・メモリも増大するため、精度とコストのバランスを考慮したメッシュ設計が重要です。
応力集中部位・境界条件の不連続点では局所的にメッシュを細密化する「適応メッシュ細化（AMR）」が有効です。
形状関数と補間
各要素内の物理量（変位・温度・圧力など）は形状関数（シェープファンクション）によって節点値から補間されます。
一次（線形）形状関数は最も単純で計算コストが低いですが、高次の形状関数は少ない要素数でより高い精度が得られます。
形状関数は要素の種類（1次元バー要素・2次元三角形・四角形要素・3次元四面体・六面体要素）によって定義されます。
アイソパラメトリック要素（形状写像と補間に同じ形状関数を使用）は曲面形状の正確な表現と数値積分の効率化を両立した現代FEMの標準的要素定式化です。
剛性マトリックスの組み立てと境界条件の適用
各要素の剛性マトリックス（要素剛性行列）を組み立てて全体剛性マトリックスKを構築します。
全体剛性マトリックスKと節点変位ベクトルu・外力ベクトルfの関係は連立方程式 Ku = f として表されます。
境界条件（変位拘束・荷重条件）を適用してこの連立方程式を解くことで、すべての節点での変位が得られます。
節点変位から各要素の応力・ひずみを計算するのが後処理（ポスト処理）のステップです。
有限要素法の精度と収束性
続いては、有限要素解析の精度に影響する因子と解の収束性の考え方を確認していきます。
有限要素解析の結果を正しく評価するためには、解の精度と収束性の概念を理解することが重要です。
h収束とp収束
有限要素解の精度を向上させる方法には大きく二通りあります。
h収束はメッシュを細分化（要素サイズhを小さく）することで精度を向上させるアプローチです。
p収束は要素サイズを維持したまま形状関数の次数p（多項式の次数）を上げることで精度を向上させるアプローチです。
hp法はhとpの両方を適応的に最適化する高精度手法で、特異点（応力集中・亀裂先端）を含む問題で効率的です。
メッシュ依存性と収束確認
有限要素解析では必ずメッシュ依存性の確認（メッシュ収束解析）を実施することが品質保証の基本です。
メッシュを段階的に細分化しながら関心量（最大応力・最大変位など）の変化を確認し、値が一定値に収束したことを確認します。
メッシュを細かくしても解が収束しない場合は、材料モデル・境界条件・要素タイプの見直しが必要です。
メッシュ収束確認を怠った解析結果は信頼性が低く、設計判断に誤りをもたらす危険があるため必ず実施すべきステップです。
数値積分とガウス求積法
有限要素法では要素剛性行列の計算に数値積分が必要で、最も広く使用されているのがガウス求積法（Gauss-Legendre quadrature）です。
ガウス積分点の数（次数）が多いほど精度が高くなりますが計算コストも増大します。
縮減積分（Reduced Integration）は積分点数を減らして計算を高速化する手法ですが、砂時計モード（Hourglass mode）と呼ばれる数値不安定モードが発生するリスクがあります。
適切な積分次数の選択は要素タイプ・問題の種類・計算コストのバランスによって決定されます。
まとめ
この記事では、有限要素法の基本概念・数値解析における位置づけ・適用分野・基本手順・精度と収束性について解説しました。
有限要素法は複雑な形状・材料・境界条件を持つ工学問題の偏微分方程式を数値的に解くための最も強力な計算ツールの一つです。
メッシュ分割・形状関数・剛性マトリックス構築・連立方程式求解という基本的な流れを理解することで、FEMソフトウェアの結果を正しく解釈・活用する力が身につきます。
構造・熱・流体・電磁気など幅広い物理現象に適用可能な有限要素法は、これからも工学設計・研究開発の中核技術として発展し続けるでしょう。
