有限要素法とは
有限要素法(Finite Element Method:FEM)は、複雑な形状や荷重条件を持つ構造物の変形・応力・温度分布などを数値的に解析する手法です。1960年代にNASAの宇宙開発プログラムで実用化され、現在は機械・土木・電磁気・流体など工学全般で不可欠なツールになっています。
FEMの基本的な考え方は「連続体を有限個の要素(エレメント)に分割し、各要素内では変位を多項式で近似する」というものです。これにより、解析的に解けない複雑な問題を連立方程式に変換して数値的に解くことができます。
FEMの基本手順
1. モデル化とメッシュ分割
解析対象をCADデータや手作業で形状モデル化し、三角形・四辺形(2D)または四面体・六面体(3D)の要素に分割します。応力集中部や境界付近はメッシュを細かく、単純な部分は粗くする「適応メッシュ」が効果的です。
2. 要素剛性行列の作成
各要素について、変位とひずみの関係(ひずみ-変位行列 [B])と材料の応力-ひずみ関係(弾性係数行列 [D])から、要素剛性行列 [k] = ∫[B]ᵀ[D][B]dV を計算します。
3. 全体剛性行列の組み立て
各要素の剛性行列を「アセンブル」して全体剛性行列 [K] を作ります。節点番号に従って各要素の寄与を加算していくプロセスです。
4. 境界条件の設定
固定端(変位ゼロ)や荷重(力・圧力)などの境界条件を設定します。これにより不定となっている剛性方程式 [K]{u} = {F} が解けるようになります。
5. 連立方程式を解く
ガウス消去法・LU分解・共役勾配法などで {u} = [K]⁻¹{F} を解き、節点変位を求めます。変位からひずみ・応力を計算して結果を可視化します。
FEM要素の種類
| 要素タイプ | 形状 | 主な用途 |
|---|---|---|
| 梁要素(Beam) | 1次元 | 棒・梁・フレーム構造 |
| 平面応力要素(CST/Q4) | 三角形・四辺形 | 薄板・平面問題 |
| シェル要素 | 曲面 | 薄肉構造・自動車ボディ |
| ソリッド要素(Tet/Hex) | 三角錐・六面体 | ブロック形状・3D全般 |
| アクシシメトリック要素 | 軸対称 | タイヤ・配管・圧力容器 |
無料でFEMを体験できるツール
以下のツールで、ブラウザ上で実際にFEMの概念を体験できます:
FEM解析の精度を上げるポイント
- メッシュ収束確認:メッシュを細かくしても結果がほとんど変わらない状態が「収束」。2〜3段階メッシュを変えて確認する。
- 要素次数の選択:線形要素よりも二次要素(節点数が多い)の方が精度が高いが計算コストも高い。
- 材料非線形:降伏後の塑性変形を扱う場合は弾塑性解析が必要。小変形仮定が崩れる場合は幾何学的非線形も考慮。
- 境界条件の妥当性:実機の拘束状態を正確にモデル化する。完全固定と単純支持で結果が大きく異なる。
FEM vs 他の数値解析手法
| 手法 | 得意分野 | 特徴 |
|---|---|---|
| 有限要素法(FEM) | 構造・熱・電磁気 | 複雑形状に対応、標準的 |
| 有限差分法(FDM) | 流体・波動 | 実装が簡単、規則格子に最適 |
| 有限体積法(FVM) | CFD全般 | 保存則を厳密に満たす |
| 境界要素法(BEM) | 電磁気・弾性 | 表面のみ離散化、無限領域に強い |
使い方ガイド
- 解析対象の形状をCADモデルで定義し、材料特性(ヤング率E、ポアソン比ν)を入力
- メッシュ生成で四面体要素(テトラヘドロン)または六面体要素を選択し、要素サイズを0.5~2mm程度に設定
- 境界条件として固定端と荷重条件を指定し、ソルバーで剛性マトリックスを組立・求解
- 応力分布(von Mises応力)と変位量を可視化し、許容応力と比較して安全率を確認
具体的な計算例
SS400鋼製の片持ちはり(長さ1000mm、矩形断面100×50mm)に先端荷重10kNを加える場合。ヤング率E=200GPa、降伏応力σy=235MPa。理論値δ=FL³/3EI≒3.2mm、最大応力σmax≒120MPaに対し、FEM解析で四面体要素2000個のメッシュを使用すると理論値と2%以内で一致。支持点近辺の応力集中は局所メッシュ細分化により正確に評価可能。
実務での注意点
- 接触問題では非線形解析が必須。ボルト締結部の接触圧力計算では摩擦係数μ=0.15~0.25を設定し反復計算回数を50回以上確保
- 熱応力解析では温度勾配が急な領域でメッシュを5倍以上細分化。アルミ合金(熱膨張係数α=23×10⁻⁶/K)の場合、温度差ΔT=100Kで熱ひずみε_th=2.3×10⁻³を導出
- 動的解析で固有振動数を求める際は、質量マトリックスの形式(consistent mass vs lumped mass)による誤差が5~10%発生するため検証必須