有限要素法(FEM)
有限要素法（FEM）は、連続体を多数の要素に分割し、場の量（変位、温度、電位など）を要素内で近似して偏微分方程式を離散化する数値解析手法である。複雑形状・異材接合・非均質材料にも適用しやすく、構造力学、熱伝導、流体、電磁場、音響など多分野の境界値問題を統一的な枠組みで解けることが強みである。弱形式（エネルギー原理・仮想仕事式）に基づき、連立一次方程式を構成して解を得るため、厳密解が難しい工学問題でも実用的な精度を達成する。
基本概念
連続体領域を三角形・四辺形・四面体・六面体などの要素に分割し、節点で未知量を定義する。要素内の場は形状関数で近似し、重み付き残差法（ガラーキン法）により弱形式を導出して要素剛性行列と要素荷重ベクトルを得る。要素をアセンブルして全体剛性方程式を構築し、境界条件を課して解くのが基本手順である。
数学的基礎
支配方程式は多くの場合、楕円型や放物型の偏微分方程式である。強形式を積分恒等式に変換し、試験関数と近似関数を同一空間にとるガラーキン法を用いることで、近似解の安定性と収束性を確保する。境界条件は本質境界（Dirichlet）と自然境界（Neumann）に分類され、前者は未知量の既知化、後者は弱形式中の境界積分として扱う。
要素と形状関数
要素は1次（線形）から高次（2次以上）まで用意され、幾何と場の近似を一致させる等写像（isoparametric）定式化が一般的である。数値積分にはガウスの数値積分を用い、参照要素への座標変換で効率的に評価する。
- 1D：ばね・棒・梁要素
- 2D：三角形（一次/二次）、四辺形（双一次/二次）
- 3D：四面体、六面体、プリズム要素
メッシュ生成と品質
解の勾配が大きい領域では要素を細分化し、滑らかな領域では粗くする適応分割が有効である。要素品質はアスペクト比、ねじれ、ヤコビアンの符号などで評価し、劣悪な要素は誤差と発散の原因となる。h-法（細分化）、p-法（次数上げ）、hp-法（併用）により精度と計算コストのバランスを最適化する。
剛性方程式と境界条件
全体方程式は一般に対称正定値の形 K u = f をとる。拘束条件はペナルティ法、ラグランジュ乗数法、多点拘束（MPC）などで実装する。接触や摩擦のような不等式制約は非線形性を生み、増分反復で解く必要がある。
解法と前処理
線形問題では直接法（Cholesky、LU）と反復法（CG、BiCGSTAB、GMRES）が用いられる。大規模疎行列には反復法と前処理（ILU、IC、AMG）が有効で、スケーリングや再番号付けで条件数を改善する。並列計算（MPI、スレッド）により大規模メッシュでも実用時間内に解ける。
非線形・動的解析
材料非線形（塑性、粘弾性、超弾性）、幾何非線形（大変形）、接触非線形に対してはNewton–Raphson法を基礎とする増分反復を行う。動的問題は時間積分で扱い、Newmarkやα法などの陰解法、中央差分などの陽解法を選択する。減衰の扱い（Rayleighなど）と時間刻みの安定条件が精度を左右する。
誤差評価と検証
誤差はエネルギーノルムや後験誤差推定で評価し、メッシュ収束性（要素を細かくしたときの解の安定）を確認する。パッチテストの合格は要素定式化の健全性を示す。モデル検証では、解析と理論解または実験の整合性を確認し、材料データ（JISやISO規格に準拠）と境界条件の妥当性を点検する。
FDM・FVMとの比較
有限差分法（FDM）は直交格子上での離散化が簡潔で離散化誤差の解析が容易である一方、複雑形状の表現が苦手である。有限体積法（FVM）は保存則をセル境界で厳密に満たしやすく、流体で強みがある。有限要素法（FEM）は任意形状への適用性と高次近似の柔軟性に優れ、多物理連成に適するが、要素選択とメッシュ設計に専門性を要する。
適用例とモデリング上の注意
構造ではシャーシ、フレーム、板殻の座屈・固有値解析、接合部の局所応力評価に用いられる。たとえばボルト・ナット接合のプリテンションやすべり限界の評価、熱では非定常伝熱と放射・対流の併用、電磁場では渦電流や静電容量の見積りが典型例である。モデル化では荷重・拘束の実在性、対称性の活用、単位系の統一、材料の温度依存や塑性則の選定、接触剛性や摩擦係数のキャリブレーションが重要である。
ポストプロセスと結果の解釈
応力特異点（切欠き先端や境界拘束点）では公称応力がメッシュ依存となるため、経路積分、平均化、応力拡大係数や評価点の規約に基づいて判断する。反力・エネルギー・質量保存の整合、メッシュ独立性、感度解析による設計変数の影響評価を併用し、設計意思決定に資する指標（安全率、許容たわみ、熱抵抗、インダクタンスなど）へ適切に落とし込む。
コメント（β版）