振動減衰シミュレーター 戻る
構造解析

振動減衰シミュレーター — SDOF減衰振動・対数減衰率

質量・剛性・減衰比を変えて不足減衰・臨界減衰・過減衰の応答をリアルタイム比較。対数減衰率δと整定時間を対話的に計算しよう。

パラメータ設定
質量 m (kg)
kg
剛性 k (N/m)
N/m
減衰比 ζ
初期変位 x₀ (m)
m
初速度 v₀ (m/s)
m/s
計算結果
減衰
減衰
理論・主要公式
$$\ddot{x}+ 2\zeta\omega_n\dot{x}+ \omega_n^2 x = 0$$

対数減衰率:

$$\delta = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}}, \quad \omega_n = \sqrt{k/m}$$

振動減衰シミュレーターとは

🙋
このシミュレーターで「減衰比ζ」を変えるとグラフが全然変わるみたいだけど、これって何ですか?
🎓
大まかに言うと、振動がどれだけ早く静まるかの「ブレーキの強さ」を表すパラメータだよ。右のスライダーでζを0.05くらいにすると振動がなかなか収まらないでしょ?これが「不足減衰」。ζを1.0にすると、一気に静まる「臨界減衰」になるんだ。実務では、自動車のサスペンションや建物の耐震設計で、この値を最適化するんだよ。
🙋
え、そうなんですか!「対数減衰率δ」って表示もありますけど、これはζとどう違うんですか?
🎓
ζが理論的な設計値だとすると、δは現場で実際に測定する値なんだ。シミュレーターでζを0.1に設定してグラフを確認してみて。山の高さが順番に小さくなってるよね。その隣り合う山の高さの比の自然対数を取ったものがδなんだ。実験で振動波形を測れば、このδから逆にζを求められる。現場で多いのは、ハンマーで構造物を叩いて出た振動波形から減衰性能を評価する方法だね。
🙋
なるほど!「整定時間」ってのも出てきますね。これは何に使うんですか?
🎓
振動が「実質的に終わった」とみなす時間だよ。例えば、精密機械が衝撃を受けた後、振動がいつまで続くかは製品の信頼性に直結する。上のパラメータで「質量m」を大きくしたり「剛性k」を小さくすると、振動自体が遅くなるから整定時間も長くなるよね。制御システムの設計では、この整定時間を短くすることが重要な設計目標になるんだ。

よくある質問

質量や剛性を変えると固有振動数が変化し、振動の速さが変わります。減衰比ζが1未満(不足減衰)では振動しながら収束、1(臨界減衰)では最も速く振動せずに収束、1超(過減衰)ではゆっくりと収束します。スライダーを動かしてリアルタイムに比較できます。
対数減衰率δは、隣り合う振幅の比の自然対数で計算します。シミュレーター上で波形のピーク値を読み取り、δ = ln(x1/x2) を求めると減衰比ζを推定できます。整定時間は応答が目標値の±2%以内に収まるまでの時間で、画面上の波形から直接確認できます。
例えば、車両のサスペンションや精密機器の防振設計において、適切な減衰比や剛性を試行錯誤する際に役立ちます。臨界減衰に近づけると振動が早く収まるため、応答速度とオーバーシュートのバランスを視覚的に確認しながらパラメータを調整できます。
ζ=0は無減衰状態で、振動が永久に減衰せず持続します。現実の機械では摩擦や空気抵抗で必ず減衰が生じるため、あくまで理論上の理想状態です。シミュレーター上では振幅が一定の正弦波が続くので、設計時には実際の減衰を考慮する必要があります。

実世界での応用

自動車のサスペンション設計:乗り心地(振動の早い収束)と接地性(振動しすぎない剛性)はトレードオフの関係にあります。減衰比ζを調整することで、スポーツカーと高級セダンでは全く異なる振動特性が設計されます。シミュレーターでmを車体重量、kをバネ定数に見立てて操作してみましょう。

建築物の耐震・耐風設計:高層ビルや橋梁は風や地震で振動します。この振動を早く減衰させるため、建物内部に巨大な質量(チューニッドマスダンパー)やオイルダンパーを設置します。実構造物の典型的な減衰比は0.01〜0.05と非常に小さい値です。

精密機械・半導体製造装置:製造中の微細な振動は製品の不良率に直結します。装置の框架(フレーム)の剛性kを高め、防振材などで減衰比ζを最適化することで、外部からの振動を素早く減衰させ、整定時間を短くします。

家電製品の振動・騒音低減:洗濯機の脱水時や掃除機のモーターには大きな振動が発生します。この振動を筐体に伝えにくくし、かつ発生した振動を素早く静めるために、質量mや支持部の剛性k、ダンパーのζが詳細に設計されます。

よくある誤解と注意点

このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるんだ。まず一つ目は、「減衰比ζが大きければ大きいほど、振動は常に早く収まる」と思いがちな点。確かにζが1.0(臨界減衰)まではそうなんだけど、ζをそれ以上(過減衰)にすると、逆に平衡位置に戻るまでの「整定時間」が長くなってしまうんだ。例えばζ=2.0にすると、振動はしない代わりに、のっそりとゆっくり戻るだけ。最速で収束させるのはζ=1.0の時だから、これが「臨界」と呼ばれる理由だね。

二つ目は、パラメータを個別にいじる時の相互作用を見落とすこと。減衰比ζは、質量m、減衰係数c、剛性kの3つから計算される(ζ = c / (2√mk))よね。だから、例えば「剛性kを2倍にしたら振動は速くなるから、減衰も強くしないと」と思ってcも2倍にしたとする。実はこれ、ζの値自体は変わらないんだ(計算式に入れてみて!)。振動の速さ(ω_n)は変わるから応答は変わるけど、「ブレーキの強さ」の割合は同じまま。パラメータを変更する時は、この式を頭の片隅においておこう。

三つ目は実務的な落とし穴。シミュレーション上の「整定時間」と、実機で測定される「振動が収まった時間」は必ずしも一致しない。シミュレーションは理想的な初期条件と数学的な定義(例えば、振幅が初期値の2%以内に入る時間)で計算する。でも現場では、測定ノイズや微小な継続振動、そもそも測定開始タイミングのずれなどがある。シミュレーション結果をそのまま信じるのではなく、「理論値としてはこのくらい」という指標として使うことが大事だよ。

使い方ガイド

  1. 質量(kg)と剛性(N/m)をスライダーで設定し、固有角振動数ωnを決定する
  2. 減衰比ζ(ゼータ)を0~2の範囲で調整して、不足減衰(ζ<1)・臨界減衰(ζ=1)・過減衰(ζ>1)の3パターンを比較観察する
  3. 初期変位x0(mm)を入力後、シミュレーション実行ボタンを押して時刻歴応答波形を確認し、対数減衰率δと整定時間(2%基準)をリアルタイム出力で読み取る

具体的な計算例

鋼製フレーム構造(質量m=50kg、剛性k=10000N/m)に対して初期変位x0=10mmを与えた場合:固有角振動数ωn=14.14rad/s、固有振動数fn=2.25Hzが自動計算される。減衰比ζ=0.1(軽度減衰)を設定すると、対数減衰率δ=0.628、1周期エネルギー損失=2.3mJとなり、整定時間は約4.5秒で振動が静止する。同じ条件でζ=1.0(臨界減衰)に変更すると、整定時間は約0.28秒に短縮され、オーバーシュートなしの滑らかな応答が得られる。

実務での注意点