板の座屈 — トラブルシューティングガイド
板座屈解析のトラブル
板座屈のFEM解析で、特に注意すべきトラブルを教えてください。
板座屈は柱座屈に比べて、FEM特有のトラブルが多い。代表的なものを見ていこう。
座屈荷重が理論値と合わない
四辺単純支持の板で $k = 4.0$ のはずが、FEMでは $k = 4.3$ になりました。
考えられる原因を順にチェック:
1. メッシュが粗い — 板座屈は最低6×6要素(2次シェル)が必要。粗いメッシュは剛性を過大評価し、座屈荷重が高くなる方向にずれる。
2. 境界条件が厳しすぎる — 「単純支持」のつもりが回転を拘束してしまっている。シェル要素の場合、面外変位($w = 0$)だけ拘束すべきところ、面外回転($\theta = 0$)まで拘束すると固定支持になる。
3. 面内拘束の影響 — 面内変位の拘束条件が理論の仮定と合っていない。面内自由な辺で面内を拘束すると剛性が上がる。
境界条件のミスが一番多いですか?
圧倒的に多い。板座屈の「単純支持」は:
| DOF | 支持辺 | 自由辺 |
|---|---|---|
| 面外変位 $w$ | 0(拘束) | 自由 |
| 面外回転 $\theta$ | 自由 | 自由 |
| 面内(荷重方向) | 荷重条件に依存 | 自由 |
| 面内(幅方向) | ポアソン膨張に注意 | 自由 |
「ポアソン膨張に注意」が盲点だ。荷重辺で幅方向の変位を拘束すると、ポアソン効果による横膨張が抑制されて座屈荷重が上がる。理論の「単純支持」はポアソン膨張を許容する条件だから、幅方向は自由にすべきだ。
座屈モードが期待と違う
$a/b = 3$ の板で3半波の座屈を期待していたのに、1半波のモードが出ました。
荷重辺の拘束条件を確認。荷重辺が面内に完全に拘束(一様変位)されていると、短辺方向の変形も拘束されて高次モードが抑制される。
もう一つの可能性はアスペクト比に対してモード数が求められている固有値の中に入っていないこと。$a/b = 3$ では $m = 3$ のモードの $k$ 値は $m = 1$ よりわずかに高いが、$m = 2$ のモードはさらに高い。固有値を多めに(20モード程度)求めて、各モードの半波数を確認しよう。
メッシュ細分化で座屈荷重が低下し続ける
メッシュを細かくするたびに座屈荷重が下がっていきます。収束する気配がありません。
応力集中部やコーナー部で局所座屈モードが出現している可能性がある。穴の周辺、荷重集中点、形状の不連続部では、メッシュを細かくするほど局所的なモードが低い固有値で出る。
対策:
- モード形状を確認して、目的の座屈モード(全体座屈)が何番目にあるかを追跡
- 応力集中部はサブモデリングで別途評価し、全体モデルでは粗いメッシュのままにする
- NAFEMSベンチマーク問題と比較して、要素タイプとメッシュ密度の妥当性を確認
複合荷重の相互作用
圧縮とせん断を同時に受ける板の座屈はどう評価しますか?
設計基準では相互作用式を使う:
$\alpha = 1, \beta = 1$ が線形相互作用(保守的)。$\alpha = 2, \beta = 2$ が二次相互作用(一般的)。
FEMでは複合荷重を直接与えて固有値解析すればいいですか?
そう。FEMの固有値座屈解析で圧縮とせん断を同時に参照荷重として与えれば、相互作用を自動的に考慮した座屈荷重が得られる。ただし荷重比が固定される点に注意。荷重比が変わる場合は複数ケースの解析が必要だ。
よくある設計ミス
板座屈の設計でよくあるミスをまとめてもらえますか?
| ミス | 結果 | 対策 |
|---|---|---|
| 座屈係数 $k$ の境界条件を間違える | 座屈応力が数倍ずれる | 支持・自由の区別を明確に |
| フランジの突出板を内部板として扱う | $k$ が過大(0.425 → 4.0) | 自由辺の有無を確認 |
| 有効幅を適用すべきClass 4断面を見逃す | 耐力過大評価(危険側) | $\bar{\lambda}_p$ を計算して判定 |
| スティフナーの偏心効果を無視 | 歪み座屈の見落とし | スティフナーの剛性と偏心を正しくモデル化 |
| せん断座屈を検討しない | ウェブの座屈を見落とす | $h_w / t_w$ が大きいI桁は必ず検討 |
「フランジの突出板を内部板として扱う」ミスが特に危険ですね。$k$ が10倍近く違う。
これは実務で本当に起きるミスだ。H形鋼のフランジは片側がウェブに接続(支持辺)、もう片側が自由端。$k = 0.425$ で評価すべきところを $k = 4.0$ にしたら、座屈応力が10倍になってしまう。断面のどの部分が「内部板」でどの部分が「突出板」かを正しく認識することが板座屈設計の第一歩だ。
タイタニック号と安全率の教訓
「不沈」と謳われたタイタニック号は、低温でのリベット材の脆性破壊が沈没の一因とされています。現代の破壊力学CAEでは、温度依存の材料特性と応力拡大係数を計算して「その温度で本当に大丈夫か?」を事前に検証できます。技術の進歩は、過去の悲劇から学んだ結果です。
トラブル解決の考え方
デバッグのイメージ
構造解析のトラブルシューティングは「医師の問診」に似ている。「いつから症状が出たか」(どのステップでエラーが出るか)、「どこが痛いか」(どの要素で収束しないか)、「何をしたか」(直前に何を変更したか)を系統的に聞くことで原因を特定する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——板の座屈の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
構造解析の収束問題や計算コストに課題を感じていませんか? — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「板の座屈をもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →