振動解析・固有値解析シミュレーター 戻る
振動・構造動解析

振動解析・固有値解析シミュレーター(多自由度)

1〜3自由度の質量バネ系で固有振動数と振動モードをリアルタイム計算。各モードのアニメーションと周波数応答グラフで振動挙動を直感的に理解できる。

自由度選択
質量 [kg]
m₁
kg
m₂
kg
m₃
kg
ばね定数 [N/m]
k₁(固定端〜m₁)
N/m
k₂(m₁〜m₂)
N/m
k₃(m₂〜m₃)
N/m
k₄(m₃〜自由端)
N/m
固有振動数
f₁
f₂
f₃
モード表示
モード形アニメーション
周波数応答関数(FRF)— m₁加振・m₁応答
理論・主要公式
$$(\mathbf{K}- \omega^2 \mathbf{M})\mathbf{X}= \mathbf{0}$$ $$\det(\mathbf{K}- \omega^2 \mathbf{M}) = 0$$ $$f_n = \frac{\omega_n}{2\pi}\text{ [Hz]}$$

振動解析・固有値解析とは

🙋
このシミュレーターで「固有振動数」って計算してるけど、それが分かると何がいいんですか?
🎓
大まかに言うと、その構造物が「揺れやすい周波数」が分かるんだ。例えば自動車のエンジンが特定の回転数でガタガタ振動するのは、その回転数が車体の固有振動数に近づいて「共振」するからだよ。上のスライダーでm1の質量を大きくしてみて。どうなる?
🙋
え、全部の周波数が下がりました!重くすると揺れにくくなるということですか?
🎓
その通り!質量が大きいと動きにくくなるから、固有振動数は低くなる。逆に、右のk1のばね定数を強く(大きく)してみて。今度は周波数が上がるだろ?実務では、重い機械を軽量化すると振動問題が起きやすくなったりするんだ。
🙋
「振動モード」ってグラフで出てくる変な形は何ですか?第1モードと第2モードで全然違いますね。
🎓
あれは、その周波数で自由に揺らした時に、各質点がどういう「形」で動くかを表しているんだ。第1モードはみんな同じ方向に揺れる。第2モードは真ん中の質点がほとんど動かず、端の2つが逆方向に動く。実際の橋やビルでも、低い周波数では全体が同じ方向に、高い周波数では複雑にねじれるモードが出てくるよ。シミュレーターのアニメーションを再生して、動きを確認してみて。

よくある質問

解析したい実際の構造物の動きを単純化して考えてください。例えば、1自由度は1つの質量の上下振動、2自由度は2つの質量が連成する振動、3自由度はより複雑な連成を表現します。自由度が多いほど実際の挙動に近づきますが、計算の解釈は複雑になります。
通常は低い次数(1次、2次、3次)の固有振動数ほど実現象で励起されやすく、設計上も重要です。このツールでは各モードのアニメーションと周波数応答グラフで、どの振動数でどのような変形が起きるかを直感的に確認できます。
単位系は統一されていれば任意ですが、一般的には質量を[kg]、ばね定数を[N/m]とすると、固有振動数は[Hz]で表示されます。数値の大小関係が正しければ、相対的な振動特性を理解する上で問題ありません。
ブラウザのハードウェアアクセラレーションが有効かご確認ください。また、タブやウィンドウを多数開いていると処理が重くなることがあります。それでも改善しない場合は、ページをリロードするか、別のブラウザ(Chrome推奨)をお試しください。

実世界での応用

自動車・航空機の設計:エンジンやプロペラの回転数が車体・機体の固有振動数と一致すると、共振による大振幅振動や騒音、疲労破壊が発生します。設計段階でCAE(モーダル解析)を用いて固有値を把握し、危険な周波数域を避けるように質量や剛性を調整します。

建築・橋梁の耐震設計:地震動は広い周波数帯域の振動を含みます。建物の固有振動数が地震動の主要な周波数と一致すると共振し、倒壊の危険が高まります。固有値解析により建物の揺れやすいモードを把握し、制振装置の設置や構造の見直しを行います。

工作機械・精密機器:マシニングセンターなどで工具が回転する際、その回転数が機械構造や工具自体の固有振動数に近づくと「びびり振動」が発生し、加工精度が著しく低下します。工具の長さや保持方法を変えて固有振動数をシフトさせる対策が取られます。

家電製品の静粛化:洗濯機の脱水時や冷蔵庫のコンプレッサー動作時に発生する振動や騒音は、内部構造の固有振動が励起されることが一因です。CAEを用いて振動モードを可視化し、防振ゴムの配置や質量バランスを最適化することで静粛性を向上させます。

よくある誤解と注意点

このツールで遊んでいると、いくつか勘違いしやすいポイントがあるんだ。まず、「固有振動数が低い(重い・ばねが柔らかい)ほど危険」と思いがちだけど、それは状況による。確かに低い周波数は日常的に発生しやすいけど、問題は「加振力の周波数」と「固有振動数」がどれだけ近いかだ。例えば、高速回転するファン(加振周波数100Hz)に、固有振動数2Hzの筐体はほとんど影響を受けない。逆に、エンジンのアイドリング振動(20Hz)と固有振動数が20Hzのマウント部品は、たとえ周波数が高くても共振してしまう。

次に、ツール上では「減衰」を考慮していない点に注意だ。現実の構造物には必ず減衰(振動をエネルギーを熱などに変えて消す効果)がある。減衰が大きいと、共振ピークの応答振幅は理論値より大幅に低くなる。このシミュレーターの周波数応答グラフは「非減衰系」の理想的な結果だから、実務ではここに減衰項を加えたモデルで再評価するのが次のステップだ。

最後に、パラメータ設定の落とし穴。質量やばね定数を「直列」や「並列」で複数設定する場合、その合成値の計算を間違えることが多い。例えば、2本のばねが質点間に「直列」なら、合成ばね定数kは $1/k = 1/k_1 + 1/k_2$ で求められる。ツールでk1とk2を別々にいじる時、この関係を頭に入れておかないと、意図しない剛性分布になってしまうよ。

使い方ガイド

  1. 質量パラメータ(vm1, vm2, vm3)をスライダーで設定します。例えば、1段目に2kg、2段目に1.5kg、3段目に1kgを入力し、多層構造の振動特性をモデル化します。
  2. ばね定数(sk1など)を調整して、支持条件を決定します。鋼製フレームの場合、典型値として10kN/m~50kN/mの範囲で設定してください。
  3. 「解析実行」ボタンで固有振動数と各モードの振動パターンを計算し、周波数応答グラフとアニメーション表示を確認します。

具体的な計算例

3自由度系(機械基礎の振動解析):m1=3kg、m2=2kg、m3=1kg、k1=15kN/m、k2=12kN/m、k3=8kN/mの場合、第1固有振動数は約8.2Hz、第2固有振動数は約15.7Hz、第3固有振動数は約22.1Hzとなります。第1モードでは全質量が同位相で振動し、第2モードでは中間層が逆位相となる特性を示します。

実務での注意点