線形重ね合わせ — CAE用語解説
線形重ね合わせ
橋の設計で「死荷重＋活荷重＋風荷重」みたいに荷重を組み合わせるとき、それぞれ別々に解析して結果を足し合わせていいんですか？
理論と物理
線形重ね合わせの基本概念
「線形重ね合わせ」という言葉をよく聞きますが、CAEの文脈では具体的に何を意味しているんですか？ただ足し合わせるだけなら、なぜ特別な名前がついているんでしょう。
良い質問だ。CAEでは、支配方程式が線形であり、かつ境界条件も線形である場合に成立する非常に強力な性質を指す。具体的には、複数の荷重ケースが個別に生じる変位や応力を計算し、それらを単純に足し合わせても、全荷重が同時に作用した場合の結果と一致する、という原理だ。例えば、トラスの節点に
「支配方程式が線形」というのは、あの
その通り。
なるほど。では、もし材料が非線形、例えば塑性変形を起こしたら、この原理は完全に成り立たなくなるんですか？
完全に崩壊する。塑性域では応力-ひずみ関係が非線形となり、剛性$\mathbf{K}$が現在の応力状態や変形履歴に依存して変化する。例えば、降伏応力$355$ MPaの鋼材で、ケースAで$300$ MPa、ケースBで$100$ MPaの応力を別々に生じさせても、両方を同時に適用すると合計応力が$400$ MPaで降伏を超え、塑性変形が生じる。この時、ケースAとBの結果を単純加算したものは、実際の非線形解析結果とは一致しない。これが線形解析と非線形解析の決定的な違いの一つだ。
数値解法と実装
FEMにおける重ね合わせの利用
FEMソルバーは、実際にどのように線形重ね合わせを利用しているんですか？ 複数の荷重ケースを効率的に解くための裏技みたいなものですか？
まさに「効率化の要」と言える技術だ。ソルバーは全体剛性マトリックス$\mathbf{K}$の分解（LU分解やコレスキー分解）に最も計算コストがかかる。重ね合わせが使える線形問題では、$\mathbf{K}$は全荷重ケースで共通だから、この高コストな分解を1回だけ行えばいい。後は、異なる右辺ベクトル（荷重ベクトル$\mathbf{F}_1, \mathbf{F}_2, ...$）に対して、前進代入と後退代入という比較的軽い計算で解$\mathbf{u}_1, \mathbf{u}_2, ...$を求められる。100ケースの荷重を解く場合、分解を100回行うのと1回で済むのとでは、計算時間が桁違いだ。
ユーザーが後で「この荷重とあの荷重を組み合わせた結果が見たい」と思った時、ソルバーは再計算せずに済む、ということですか？
その通り。これが「コンビネーション」や「エンベロープ」解析の基礎だ。例えばAnsys Mechanicalでは、静構造解析で「Solution Combination」機能を使い、荷重ケース1（自重）の結果とケース2（風圧）の結果に係数（例えば1.0と1.5）を掛けて足し合わせ、設計荷重ケースを事後的に作成できる。ソルバーは再計算せず、既にある結果ファイル（.rst等）から瞬時に新しい結果を合成する。この際、応力の合成には注意が必要で、ミーゼス応力のようなスカラー量は単純加算できないが、応力テンソルの成分（$\sigma_{xx}, \tau_{xy}$など）ならば線形加算が可能だ。
実践ガイド
解析ワークフローでの適用判断
実際の解析で、今解こうとしている問題に線形重ね合わせが適用できるか、どう見極めればいいですか？ チェックすべきポイントは？
以下の4点を厳密にチェックするリストを使うといい。1. **材料の線形性**： 材料が線形弾性（フックの法則）か。塑性、超弾性、クリープはNG。2. **幾何非線形性の有無**： 変形が微小で、釣り合い方程式が変形前の形状で立てて良いか。大変形・座屈はNG。3. **接触条件**： 接触状態が荷重によって変化しないか。接触は本質的に非線形だ。4. **境界条件の線形性**： 支持条件が変位に依存しないか。例えば、非線形のばねはNG。
接触が少しでもあると完全にアウトなんですか？ 最初から接触している「面と面の固着」のような状態でもダメですか？
そこが難しいところだ。初期から完全に固着（Bonded）され、解析中に剥離や滑りが絶対に起こらないと保証できるなら、それは単なる節点の結合と同じだから線形とみなせる。しかし、多くの実務的な接触（Even Bonded）では、接触圧力の分布が荷重によって非線形に変化する可能性がある。保守的に判断するなら、Abaqus/Standardで「Linear perturbation step」を使うか、Ansysで「Linear Contact」オプションを明示的に有効にしない限り、接触を含む解析は非線形として扱うべきだ。規格ASME BPVC Section VIII Div.2の設計評価でも、接触領域の応力評価には特別な注意が喚起されている。
ソフトウェア比較
各ソルバーの実装と機能
Ansys、Abaqus、COMSOLといった主要ソフトでは、線形重ね合わせに関して機能や呼び名に違いはありますか？
コンセプトは同じだが、実装とユーザーインターフェースに違いがある。
2. **Abaqus/Standard**: 「線形摂動ステップ（Linear perturbation step）」がこれに該当する。ベースステップ（非線形でも可）の上に、複数の線形摂動ステップを重ねられ、各ステップの結果は独立で重ね合わせ可能。これが非常に明確なモデリング体系だ。
3. **COMSOL Multiphysics**: 「ステディステディ」研究で、複数の「パラメータ化されたソルバー」を定義するか、「補助掃引」機能を使う。複数の荷重やパラメータ値に対する線形解を一度の計算で効率的に得られる。
Abaqusの「線形摂動ステップ」と、通常の「静解析ステップ」の根本的な違いは何ですか？ なぜ分けているんですか？
計算アルゴリズムが全く異なるからだ。通常の「静解析ステップ」（General）は非線形方程式をニュートン・ラフソン法で反復求解する。一方、「線形摂動ステップ」は、直前のステップ（ベース状態）の剛性マトリックス$\mathbf{K}$を**固定して**、新しい荷重に対する線形方程式
トラブルシューティング
よくある誤りと対策
線形重ね合わせを使って結果を合成したら、明らかにおかしい値（例えば、圧縮応力なのに引張りになっている）が出てきました。考えられる原因は？
まず疑うべきは「合成対象の物理量の種類」だ。先ほども触れたが、**ミーゼス等価応力**のような正のスカラー量は絶対に加算してはいけない。ケースAで$100$ MPa、ケースBで$80$ MPaのミーゼス応力が生じる場所でも、合成結果が$180$ MPaになるとは限らない（多くの場合、それより小さくなる）。合成できるのは、**応力テンソル成分**（SX, SY, SXY...）、**変位**、**反力**といったベクトルまたはテンソル成分だけだ。おかしな結果が出たら、まず合成している物理量が成分データか確認せよ。後処理で「Component」表示に切り替える必要がある。
ソフトが「線形解析」と表示しているのに、二つの荷重ケースの結果を足し合わせたら、別途「両方の荷重を同時に載せた単一の解析」の結果と微妙に違いました。これはバグですか？
バグではなく、境界条件や要素特性に潜む「隠れ非線形性」が原因であることがほとんどだ。具体的には：1. **大変形効果（幾何非線形）が有効になっていないか**： ソフトのデフォルト設定を確認。Ansysでは`NLGEOM,OFF`、Abaqusでは`NLGEOM=NO`。2. **圧力荷重の方向**： 圧力荷重が「要素の法線方向」に作用する設定の場合、変形によって圧力の向きが変わり、非線形性を生む。これを「追随荷重」という。3. **張力のみを支持するケーブル要素やGap要素**： これらは荷重の向きで剛性が0か有限か変わる、明らかな非線形要素だ。これらのいずれかが存在すると、原理上、重ね合わせは成立しない。
追随荷重の影響を定量的に見積もる方法はありますか？ 無視できる基準は？
一つの目安は、荷重面の回転角度だ。変形による面の回転角$\theta$が非常に小さければ（例えば、$\theta < 1^\circ$）、$\cos\theta \approx 1$、$\sin\theta \approx 0$とみなせる。具体的には、圧力$p$が作用する面で、変位$u$、代表長さ$L$に対して、
関連トピック
なった
詳しく
報告