1 概要
有限要素法は、連続体や構造物を細かな部分に分け、各部分の振る舞いを近似して全体を解く数値計算法である。複雑な形状、材料の違い、境界条件の多様さを扱いやすく、工学分野を中心に広く普及している。
この手法では、対象領域を要素の集合として表し、節点で結んだ離散モデルを用いる。解析対象は静的問題だけでなく、時間変化を伴う現象や非線形性を含む場合にも及ぶ。
1.1 定義
有限要素法は、連続的な物理量を有限個の自由度に置き換えて解く近似解法である。未知量は通常、節点における変位、温度、電位などとして表される。各要素内では形状関数によって場の変化を補間し、全体の方程式を組み立てる。
1.2 歴史
起源は20世紀中頃の構造解析にあり、航空機や土木構造の設計需要とともに発展した。初期には平面応力やトラス解析が中心だったが、計算機の性能向上により、三次元連続体、動的問題、材料非線形まで対象が広がった。現在では、工学計算の標準的手法の一つとして位置づけられている。
1.3 他の数値解析手法との関係
有限要素法は、有限差分法や境界要素法と並ぶ代表的な離散化手法である。有限差分法が格子点上の微分近似を重視するのに対し、有限要素法は要素ごとの近似と弱形式に基づく点が特徴である。境界要素法は境界のみを離散化するため、無限領域や外部場に適する一方、有限要素法は内部の分布量を詳細に扱いやすい。
2 基本原理
有限要素法の核心は、複雑な連続問題を小さな領域の集合として捉え直すことにある。各部分の近似解を組み合わせることで、元の問題に近い解を得る。
2.1 領域分割
解析対象は、三角形、四辺形、四面体などの単純な形状の要素に分割される。分割が細かいほど局所的な変化を表現しやすくなるが、計算量は増える。要素の配置はメッシュと呼ばれ、解析精度と効率を左右する重要な要素である。
2.2 近似関数
要素内の未知量は、節点値から決まる補間関数で表される。これらの関数は形状関数と呼ばれ、要素の次数に応じて一次、二次、それ以上の表現が可能である。適切な形状関数を選ぶことで、局所的な変化を滑らかに近似できる。
2.3 変分原理
多くの有限要素定式化は、物理法則を変分原理の形に書き換えることで導かれる。これは、解を直接微分方程式で扱うのではなく、全体のエネルギーや仕事の釣り合いから求める考え方である。構造力学では、最小ポテンシャルエネルギー原理が典型例となる。
2.4 弱形式
弱形式は、支配方程式に重み関数を掛けて領域全体で積分した表現である。これにより、必要な微分の階数が下がり、近似関数を導入しやすくなる。境界条件の扱いも整理され、数値計算への適合性が高まる。
3 数学的定式化
有限要素法では、物理現象を数式として定義し、それを離散化して解く。ここでは微分方程式、境界条件、行列形式への変換が中心となる。
3.1 支配方程式
対象となる現象は、釣り合い式、熱伝導方程式、波動方程式などの支配方程式で表される。構造解析では、外力と内部応力の釣り合いが基本となる。これらを要素単位で積分し、全体の方程式へ展開する。
3.2 境界条件
境界条件は、解析領域の外縁で与える制約であり、解の一意性に不可欠である。拘束の有無や外力の指定によって、問題の性質が大きく変わる。
3.2.1 ディリクレ境界条件
ディリクレ境界条件は、変位や温度などの値そのものを指定する条件である。構造問題では固定端のような拘束に対応する。数値計算では、節点自由度を直接与えることで実装される。
3.2.2 ノイマン境界条件
ノイマン境界条件は、力、熱流束、応力のような法線方向の流束を指定する条件である。外力が作用する自由表面や、熱の出入りを表す際に用いられる。弱形式では自然に現れることが多い。
3.3 離散化
離散化は、連続問題を有限個の未知数の問題に置き換える過程である。これにより、微分方程式は代数方程式系へ変換される。
3.3.1 節点と要素
節点は未知量を定義する代表点であり、要素は節点同士を結ぶ基本単位である。各要素は局所座標系を持ち、内部の場を節点値から補間する。全体の計算は、要素ごとの情報を結合して進められる。
3.3.2 剛性行列
剛性行列は、変位に対する抵抗特性を表す中心的な行列である。各要素の剛性を組み合わせて全体剛性行列を構成し、釣り合い式を行列表現にまとめる。対称性や疎性を持つことが多く、効率的な解法が重要となる。
3.3.3 荷重ベクトル
荷重ベクトルは、外力や体積力、境界での作用を離散化したものである。集中荷重だけでなく、分布荷重や慣性力も含まれる。要素ごとの寄与を加算して全体の右辺項を作る。
4 要素と形状関数
要素の種類と形状関数の選択は、解析の精度と安定性を左右する。対象の次元や変形の特徴に応じて使い分けられる。
4.1 一次元要素
一次元要素は、棒、ばね、梁などの細長い構造に適用される。節点数が少なく、基本概念を理解するうえでわかりやすい。梁要素では、曲げ変形や回転自由度も考慮される。
4.2 二次元要素
二次元要素は、平面応力、平面ひずみ、薄板などの解析に用いられる。面内の変形を効率よく表現でき、構造解析で頻繁に利用される。
4.2.1 三角形要素
三角形要素は、複雑な領域の分割に適しており、メッシュ生成の自由度が高い。一次要素は単純だが、二次以上では曲線境界の表現力が増す。局所的に細かな形状へ追従しやすい。
4.2.2 四辺形要素
四辺形要素は、規則的なメッシュで高い精度を得やすい。辺方向の補間特性が明瞭で、応力分布の滑らかな表現に向く。変形が大きい領域では、要素の歪みに注意が必要となる。
4.3 三次元要素
三次元要素は、立体構造や内部応力の詳細解析に使われる。部材の厚みや局所的な応力集中を扱う際に有効である。
4.3.1 四面体要素
四面体要素は、複雑な三次元形状への適応性が高い。自動メッシュ生成との相性がよく、実務で広く用いられる。高精度化には、要素次数や分割密度の調整が必要になる。
4.3.2 六面体要素
六面体要素は、整った形状のモデルで高い計算効率を示す。応力やひずみの表現が安定しやすく、良質なメッシュが得られれば精度面で有利である。生成には手間がかかることが多い。
4.4 高次要素
高次要素は、要素内でより高い次数の形状関数を用いる方式である。少ない要素数でも滑らかな近似が可能になり、曲率の大きい変形や応力勾配の再現に適する。ただし、計算負荷と実装の複雑さは増す。
5 解析の流れ
有限要素解析は、モデル化から結果評価まで一連の手順で進む。各段階の品質が、最終的な信頼性に直接影響する。
5.1 前処理
前処理では、解析対象を計算用のモデルへ整える。形状の単純化、境界の整理、材料情報の入力などが含まれる。
5.1.1 モデル作成
モデル作成では、実物や設計案を解析可能な形へ抽象化する。不要な細部を省きつつ、重要な寸法や接触面は保持する。ここでの判断が、結果の解釈を左右する。
5.1.2 メッシュ生成
メッシュ生成では、対象を要素へ分割する。要素サイズ、形状の均一性、局所的な細分化が重要である。応力集中部や境界近傍は、細かい分割が必要になりやすい。
5.1.3 材料定数の設定
材料定数の設定では、ヤング率、ポアソン比、密度、降伏応力などを指定する。解析目的に応じて、線形材料モデルか非線形モデルかを選ぶ。値の妥当性は、実験データや規格情報で確認される。
5.2 求解
求解は、離散化された方程式系を計算機で解く段階である。問題の規模や性質に応じて、直接法や反復法が使い分けられる。
5.2.1 連立方程式の解法
連立方程式の解法には、ガウス消去法、LU分解、共役勾配法などがある。大規模問題では、疎行列を活用した効率的手法が重要となる。メモリ使用量と計算時間のバランスが実務上の焦点である。
5.2.2 線形解析
線形解析は、材料特性や変形が比例関係を保つと仮定する方法である。解釈が比較的容易で、初期設計や概略評価に適している。小変形範囲では十分な近似を与えることが多い。
5.2.3 非線形解析
非線形解析は、材料、幾何、接触などの非線形性を含む問題に対応する。反復計算を要するため、収束判定や初期値の設定が重要になる。現実的な挙動を表現できる一方、計算負荷は増大する。
5.3 後処理
後処理では、得られた数値解を物理量として解釈し、設計判断に結び付ける。数値結果の読み取りには、表示方法の工夫も求められる。
5.3.1 変位の評価
変位の評価では、節点の移動量や部材のたわみを確認する。最大変位の位置や変形の傾向は、使用性や干渉の判断に役立つ。全体の変形模式図は、理解の補助になる。
5.3.2 応力の評価
応力の評価では、ミーゼス応力、主応力、せん断応力などを調べる。局所的な集中が現れやすく、要素境界付近の解釈には注意が必要である。許容値との比較により安全性を判断する。
5.3.3 可視化
可視化では、等高線図、変形図、ベクトル表示などを用いて結果を示す。数値の大小だけでなく、分布の連続性や局所傾向を把握しやすくなる。実務では報告資料の基礎表現としても重視される。
6 構造工学への応用
構造工学では、有限要素法が設計・解析・評価の中心的手段となっている。部材から複雑な構造体まで、幅広い対象に適用される。
6.1 骨組構造の解析
骨組構造の解析では、梁や柱を線要素として表し、節点で接合された構造を扱う。部材力、節点変位、全体剛性を効率よく求められる。建築物や架台、橋梁の初期設計でよく利用される。
6.2 板・殻構造の解析
板・殻構造の解析は、薄い構造物の曲げや面内変形を対象とする。自動車外板、船体、屋根構造などで重要である。膜応力と曲げ応力の両方を考慮できる点が特徴である。
6.3 連続体構造の解析
連続体構造の解析では、部材を線ではなく体積をもつ領域として扱う。局所的な応力集中、穴周辺の挙動、接合部の詳細評価に向く。複雑な三次元状態を直接表現できる。
6.4 固有値解析
固有値解析は、構造が持つ固有振動数や座屈荷重を求める方法である。設計上の不安定現象を事前に把握するために重要である。
6.4.1 振動解析
振動解析では、構造の自然振動モードと固有振動数を求める。地震、回転、機械振動への応答評価に関係する。共振回避や制振設計の基礎情報となる。
6.4.2 座屈解析
座屈解析では、圧縮荷重下での安定性喪失を調べる。細長い部材や薄肉構造で特に重要である。線形座屈は初期推定に有用であり、より詳細には非線形解析が用いられる。
7 高度な解析
高度な解析では、現実の材料や構造が示す複雑な挙動を扱う。計算の難度は高いが、実際の挙動に近い結果を得やすい。
7.1 材料非線形
材料非線形は、応力とひずみの関係が単純な比例にならない現象である。塑性、粘性、損傷などが代表例である。
7.1.1 弾塑性解析
弾塑性解析では、弾性変形に加えて塑性変形を考慮する。降伏後の永久変形や荷重履歴の影響を表現できる。金属部材の評価で頻繁に用いられる。
7.1.2 クリープ解析
クリープ解析は、一定荷重の下で時間とともに進む変形を扱う。高温環境で顕著になり、長期使用の安全性に関わる。時間依存の材料モデルが必要となる。
7.2 幾何学的非線形
幾何学的非線形は、大変形や回転が無視できない場合に生じる。変形後の形状変化が剛性や応力状態に影響する。薄膜、柔軟構造、座屈後挙動の解析で重要である。
7.3 接触解析
接触解析は、物体同士が接触・分離・滑りを起こす問題を対象とする。条件設定が難しく、反復計算も複雑になりやすい。歯車、締結部、衝突前後の評価に利用される。
7.4 動解析
動解析は、慣性力を含む時間依存問題を扱う。外力の変化や運動の履歴を追跡できる。
7.4.1 陽解法
陽解法は、次時刻の状態を現在値から明示的に求める計算法である。各ステップの計算は比較的軽いが、時間刻みの制約が厳しい。衝撃や短時間現象に適している。
7.4.2 陰解法
陰解法は、次時刻の未知量を方程式として解く方式である。安定性に優れ、大きめの時間刻みを取りやすい。計算ごとの負荷は増えるが、長時間解析では有利なことが多い。
8 数値計算上の注意
有限要素法の結果は、設定や計算条件に強く左右される。数値的な性質を理解しなければ、見かけ上の精度に惑わされやすい。
8.1 収束性
収束性は、メッシュを細かくしたときに解が安定した値へ近づく性質を指す。解析目的に対して適切に収束しているかを確認することが重要である。誤差の減少傾向を追うことで、信頼度を見積もれる。
8.2 安定性
安定性は、計算誤差が増幅せずに抑えられる性質である。解法や時間積分法、要素選択によって左右される。条件が不適切だと、発散や振動的な不自然解が生じることがある。
8.3 精度評価
精度評価では、解析結果がどの程度現実に近いかを判断する。エネルギー誤差、応力誤差、比較試験などが用いられる。単一の数値だけでなく、分布の妥当性も見る必要がある。
8.4 メッシュ依存性
メッシュ依存性は、要素分割の仕方によって結果が変わる現象である。特に応力集中や局所損傷の評価では顕著である。要素サイズの変更に対する感度を確認し、必要に応じて局所的に細分化する。
9 応用分野
有限要素法は、構造力学に限らず多様な工学分野で使われる。共通する利点は、複雑な境界や不均一な物性を扱いやすいことである。
9.1 土木工学
土木工学では、橋梁、トンネル、ダム、地盤構造物の解析に用いられる。長大構造や複合的な荷重条件を扱うため、実務上の重要性が高い。耐震性能や耐久性の検討にも活用される。
9.2 機械工学
機械工学では、機械部品、筐体、締結部、回転機械の解析に使われる。強度確認だけでなく、振動低減や軽量化設計にも役立つ。試作回数を減らす手段としても有効である。
9.3 航空宇宙工学
航空宇宙工学では、機体構造、翼、ロケット部材の評価に応用される。軽量で高信頼な設計が求められるため、詳細なシミュレーションが不可欠である。熱や振動との連成解析も重要になる。
9.4 生体力学
生体力学では、骨、関節、軟部組織などの力学特性を近似する。個体差や材料の異方性が大きく、モデル化には注意が要る。医療機器設計や術前計画の補助にも利用される。
10 ソフトウェアと実務
実務での有限要素法は、理論だけでなくソフトウェア運用と品質管理によって支えられている。解析者には、数値計算と工学的判断の両方が求められる。
10.1 解析ソフトウェア
解析ソフトウェアは、前処理、求解、後処理を一体化して提供する。市販品から研究用まで種類は幅広く、用途に応じて機能が異なる。操作性だけでなく、使用できる要素、非線形機能、連成解析の範囲も選定基準となる。
10.2 モデル化の限界
モデル化には、理想化や簡略化が伴う。材料の不均一性、製造ばらつき、実際の接触状態などは完全には再現できない。したがって、解析結果は絶対値としてではなく、仮定付きの予測として扱う必要がある。
10.3 結果の検証
結果の検証では、荷重条件、拘束条件、単位系、メッシュ品質を改めて確認する。力の釣り合い、反力、変形の整合性も重要な点である。異常な値が出た場合は、設定の見直しが欠かせない。
10.4 実験との比較
実験との比較は、解析の信頼性を高める基本手段である。変位、ひずみ、固有振動数などを測定値と照合し、差の原因を検討する。解析と実験を往復させることで、材料定数や境界条件の精度が向上する。