固有振動数解析 — トラブルシューティングガイド
固有振動数解析のトラブル
固有振動数解析でよくあるトラブルを教えてください。
固有振動数のトラブルは入力データの間違いがほとんどだ。
固有振動数がゼロ
固有振動数が0 Hzのモードが出ます(拘束しているのに)。
剛体モードが残っている。拘束が不足して構造が動ける方向がある。
確認:
- 6自由度全てが拘束されているか(並進3+回転3)
- RBE2のスレーブDOFが正しいか
- メカニズム(不安定な構造)がないか
自由境界解析($f = 0$ のモードを含む)をするつもりだったのに、ゼロモードが7個以上出るのは?
メカニズムがある。構造の一部が結合されていない(接続忘れ)。モード形状を見れば、どの部分が分離しているかわかる。
固有振動数が理論値と大きくずれる
密度とヤング率の両方が $f$ に効くんですね。
$f \propto \sqrt{E/\rho}$。密度を10倍間違えると $f$ が $1/\sqrt{10} \approx 0.32$ 倍。ヤング率を1000倍間違え(MPa→Pa)ると $f$ が $\sqrt{1000} \approx 32$ 倍。どちらも桁違いのずれになる。
モードが期待と違う
1次モードが全体曲げではなく局所モードです。
薄いパネルや軽い付属物が全体モードより低い振動数で振動することがある。FEMでは「低い順にモードを出す」ため、局所モードが先に出ると全体モードが後回しになる。
対策:
- モード数を増やして全体モードを探す
- 局所モードの原因(薄いパネル、軽い部品)を特定
- 必要に応じてモデルを簡略化(局所モードの原因を除去)
固有振動数が実験と合わない
FEMと実験で10%以上の差があります。
FEMと実験の5〜10%の差は一般的。10%以上の差がある場合:
- FEMの境界条件が実験と異なる — 実験の支持条件がピンでも固定でもない中間状態
- FEMの質量が不正確 — 非構造質量(配線、塗装、付帯物)の見落とし
- FEMの剛性が不正確 — 接合部の剛性(RBE2→RBE3で変わる)
- 減衰の影響 — 高減衰では固有振動数がずれる($f_d = f_n \sqrt{1-\zeta^2}$)
まとめ
固有振動数解析のトラブル対処、整理します。
- $f = 0$ のモード → 拘束不足(剛体モード)。メカニズムの確認
- $f$ が桁違い → $\rho, E$ の単位を確認。$f \propto \sqrt{E/\rho}$
- 局所モードが先行 → モード数を増やす。原因を特定
- 実験との不一致 → 境界条件、非構造質量、接合部剛性を確認
- 質量サマリーが全てのデバッグの出発点 — $\rho$ と集中質量の確認
「密度の設定確認」が固有振動数解析の第一歩。これだけは忘れません。
密度の設定忘れは固有振動数解析の最も多いミスだ。静解析から動的解析に移行するとき、密度の入力を確認する習慣をつけよう。
タイタニック号と安全率の教訓
「不沈」と謳われたタイタニック号は、低温でのリベット材の脆性破壊が沈没の一因とされています。現代の破壊力学CAEでは、温度依存の材料特性と応力拡大係数を計算して「その温度で本当に大丈夫か?」を事前に検証できます。技術の進歩は、過去の悲劇から学んだ結果です。
トラブル解決の考え方
デバッグのイメージ
構造解析のトラブルシューティングは「医師の問診」に似ている。「いつから症状が出たか」(どのステップでエラーが出るか)、「どこが痛いか」(どの要素で収束しないか)、「何をしたか」(直前に何を変更したか)を系統的に聞くことで原因を特定する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——固有振動数解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
構造解析の収束問題や計算コストに課題を感じていませんか? — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「固有振動数解析をもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →