線形座屈（固有値座屈）解析 — トラブルシューティングガイド
より充実した内容を linear-buckling.html でご覧いただけます。
固有値座屈解析の典型的トラブル
固有値座屈解析を回していて、結果がおかしいときの対処法を教えてください。
固有値座屈特有のトラブルを整理しよう。オイラー座屈のトラブルシューティングとは違う、固有値解析ならではの問題がある。
固有値がゼロまたはゼロ付近に大量に出る
固有値が 0.001 とか 0.0002 とか、ほぼゼロの値が大量に出ました…。
これは前段の静解析で有意な応力が発生していない場合に起きる。$[K_\sigma]$ がほぼゼロ行列なので、固有値が極端に大きいか小さい値になる。
原因チェック：
1. 参照荷重がゼロまたは極小 — 荷重の値や方向を確認
2. 荷重が全て拘束点に作用している — 反力だけで平衡し、構造内部に応力が生じない
3. 剛体移動が拘束されていない — 静解析が特異になっている
4. Ansys APDLで PSTRES, ON を忘れている — 応力剛性が計算されない
4番は盲点ですね。Ansys固有の問題ですか？
そう。Nastranの SOL 105 やAbaqusの *BUCKLE は自動的に応力剛性を計算するが、AnsysのAPDLでは PSTRES, ON を明示的に指定しないと幾何剛性マトリクスが計算されない。Workbenchでは自動だが、APDLスクリプトでは忘れがちだ。
固有値が全て負になる
全モードの固有値がマイナスです。これは何を意味しますか？
参照荷重の方向が意図と逆になっている可能性が高い。例えば圧縮したいのに引張荷重を与えている場合だ。
荷重の符号を反転すれば正の固有値が出ますか？
そう。あるいは、参照荷重はそのままで $\lambda$ を負で読むこともできる。$\lambda = -3.0$ は「反対方向の3倍の荷重で座屈」という意味だ。ただし解釈が紛らわしいので、最初から荷重方向を合わせる方が安全だ。
1次モードと2次モードの固有値が非常に近い
$\lambda_1 = 3.42$、$\lambda_2 = 3.44$ のように、ほとんど同じ値のモードが出ました。
これが固有値のクラスタリングだ。物理的には「ほぼ同時に2つの座屈モードが発生する」状況で、モード相互作用による耐力低下が懸念される。
対処法：
- モード形状を比較 — 同じ形ならメッシュの対称性による重複（問題なし）。異なる形なら要注意
- 対称モードと反対称モード — 対称構造では対称/反対称がペアで出やすい
- 非線形解析で検証 — 両モードを初期不整として重ね合わせ、相互作用の影響を確認
非線形では、2つのモードをどう重ねるんですか？
$\{u_{impf}\} = \alpha_1 \{\phi_1\} + \alpha_2 \{\phi_2\}$ として、$\alpha_1, \alpha_2$ の比率を変えた複数ケースを実施する。例えば $\alpha_1 : \alpha_2 = 1:0, 1:1, 0:1$ の3ケースが最低限。最も低い崩壊荷重を設計値とする。
局所モードが大量に出て全体モードが見つからない
固有値が低い順に並べると、最初の20モードが全部フランジの局所座屈で、全体座屈が見つかりません。
これは薄肉構造でよくある問題だ。フランジやウェブの板要素が局所的に座屈するモードが、全体座屈より低い固有値を持つことがある。
対策：
1. 求めるモード数を大幅に増やす（50〜100モード） — 全体モードが上位に埋もれている
2. シフト-反転法で特定の固有値域を探索 — 全体座屈の概算値付近にシフトを設定
3. サブモデリング — 全体モデル（梁要素）で全体座屈を評価、局所モデル（シェル要素）で局所座屈を別途評価
サブモデリングが最もクリーンな方法に見えます。
その通り。大規模な構造（ビル全体、橋梁全体など）では、全体を詳細なシェルモデルにするのは現実的でない。梁モデルで全体座屈を確認し、クリティカルな部材だけシェルで局所座屈を確認するマルチスケールアプローチが実務的だ。
プリロードの有無で結果が大幅に変わる
プリロード（自重）を入れると座屈荷重が2割も下がりました。なぜですか？
プリロードによる応力が $[K_\sigma]$ に反映されるからだ。自重で柱に圧縮が入っていれば、追加荷重に対する座屈余裕は当然小さくなる。
逆に、プリロードが引張（例：プレストレス）だと座屈荷重が上がる。プレストレスコンクリート柱が座屈に強い理由の一つはこれだ。
プリロードを無視して座屈解析すると、危険側の過大評価になるんですね。
必ずしも過大評価とは限らない。プリロードが引張方向に効く部材もあるから。ただし自重のような永続荷重を無視する理由はないので、常にプリロードを含めるのが安全だ。
ベンチマーク検証
自分の設定が正しいか検証するためのベンチマーク問題はありますか？
NAFEMSが公開しているベンチマーク問題が最も信頼性が高い：
| ベンチマーク | 構造 | 理論解 | 出典 |
|---|---|---|---|
| NAFEMS BM-1 | 正方形板の面内圧縮座屈 | $\sigma_{cr} = 4.0 \times \frac{\pi^2 D}{b^2 t}$ | NAFEMS R0015 |
| NAFEMS BM-3 | 円筒シェルの外圧座屈 | $p_{cr}$ の理論式 | NAFEMS R0015 |
| Euler柱 | 両端ピン柱 | $P_{cr} = \pi^2 EI / L^2$ | 古典理論 |
まずEuler柱で合わせて、次にNAFEMSの板座屈で確認するのが手堅いですね。
そう。新しいソルバーや新しい要素タイプを使い始めるときは、必ずベンチマーク検証をしてから実問題に入ること。座屈解析は設定のミスが結果に直結する分野だから、検証のステップを省いてはいけない。
座屈解析のトラブルシューティング、ようやく体系的に理解できました。ポイントは「固有値の意味を物理的に解釈する」ことですね。
その通り。数字だけ見て判断するのではなく、なぜその固有値が出たのかを常に考える習慣をつけてほしい。そうすればトラブルの原因は自然と見えてくる。
線形座屈解析の結果が非保守的な場合
薄肉シェル（円筒・球殻）の線形座屈荷重は実験値の2〜8倍と大きく非保守的になることがある（ノック-ダウンファクター問題）。初期不整・残留応力の影響が極めて大きいためで、NASA SP-8007の経験式ノックダウンファクター（円筒では0.6〜0.9）を乗じる補正が必要だ。計算値だけを使った設計は危険で、不整の影響を考慮した非線形解析が原則だ。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——線形座屈（固有値座屈）解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告