減衰自由振動(減衰固有振動数)シミュレーター 戻る
振動工学

減衰自由振動(減衰固有振動数)シミュレーター

ばね・質量・ダンパからなる1自由度系を弾いて手を離したときの「減衰自由振動」を体験するツールです。質量・ばね定数・減衰係数を変えると、非減衰固有振動数・減衰比 ζ・減衰固有振動数・対数減衰率 δ がリアルタイムで分かり、振幅が指数関数的に減っていく振動波形を見られます。

パラメータ設定
質量 m
kg
振動する物体の等価質量
ばね定数 k
N/m
系の復元剛性。大きいほど速く振動する
減衰係数 c
N·s/m
ダッシュポット等が振動エネルギーを奪う強さ
評価サイクル数 n
何サイクル後の振幅比を評価するか
計算結果
非減衰固有振動数 (Hz)
減衰比 ζ
減衰固有振動数 (Hz)
対数減衰率 δ
振幅半減サイクル数 (回)
n サイクル後の振幅比 (%)
ばね・質量・ダンパ系と減衰振動 — アニメーション

質量を弾いて手を離すと、減衰固有振動数で振動しながら振幅が指数関数的に減衰していきます。減衰包絡線(破線)と連続するピーク、対数減衰率 δ を表示します。

減衰自由振動の波形 — 変位 vs 時間
減衰固有振動数 vs 減衰比 ζ
理論・主要公式

$$\omega_n=\sqrt{\frac{k}{m}},\qquad \zeta=\frac{c}{2\sqrt{m\,k}}$$

非減衰固有角振動数 ωₙ と減衰比 ζ。m:質量、k:ばね定数、c:減衰係数。臨界減衰係数は c_c = 2√(mk) で、ζ = c/c_c。

$$\omega_d=\omega_n\sqrt{1-\zeta^{2}},\qquad \delta=\frac{2\pi\zeta}{\sqrt{1-\zeta^{2}}}$$

減衰固有角振動数 ω_d(ζ < 1 の不足減衰で有効)と対数減衰率 δ。減衰は振動数をわずかに下げ、δ は振幅が1サイクルごとに減る速さを表す。

$$x(t)=x_0\,e^{-\zeta\omega_n t}\cos(\omega_d t),\qquad \frac{x_n}{x_0}=e^{-n\delta}$$

減衰自由振動の変位 x(t) と、n サイクル後の振幅比。包絡線 e^(−ζωₙt) が振動を閉じ込めながら指数関数的に縮む。

減衰固有振動数とは

🙋
ばねに付いたおもりを引っ張って離すと、しばらく揺れて、だんだん止まりますよね。あの「だんだん止まる」のが減衰ってやつですか?
🎓
そのとおり。摩擦も空気抵抗も何もない理想の世界なら、一度弾いたおもりは固有振動数のまま永遠に揺れ続ける。でも現実には必ず「減衰」がある。空気抵抗、接触面の摩擦、材料内部のエネルギー損失、あるいはわざと付けたダッシュポット。これらが振動のエネルギーを少しずつ奪うから、現実の自由振動は「鳴るけれど、ひと揺れごとに少しずつ小さくなって、やがて消える」減衰振動になるんだ。
🙋
なるほど。でも「減衰固有振動数」って言葉、固有振動数とは別物なんですか?
🎓
いい質問。実は減衰は2つの違うことを同時にやっていて、ここを混同しないのが大事なんだ。1つめ、減衰は振動を「ほんの少し遅くする」。減衰がない理想の固有振動数 fₙ ではなく、ちょっとだけ低い「減衰固有振動数」f_d で揺れる。式で言うと f_d = fₙ·√(1−ζ²)。ζ は減衰比だ。ただし多くの実構造の軽い減衰では、この差はごくわずか。減衰比5%でも振動数が下がるのはたった0.1%くらいだから、実務では計算がラクな非減衰の fₙ をそのまま使えることが多い。
🙋
差がそんなに小さいなら、減衰固有振動数って気にしなくてもいいんですか?
🎓
軽減衰ならね。でも減衰がどんどん大きくなると、その差は無視できなくなる。そしてある特別な値——臨界減衰、つまり減衰比がちょうど 1 のところ——で、減衰固有振動数はゼロまで落ちて、系はもう振動しなくなる。ただスーッと静止位置に戻るだけ。左の減衰係数 c を最大まで上げてみると、減衰比が 1 を超えて減衰固有振動数が 0 Hz になり、波形が「振動なし」に切り替わるのが見えるよ。
🙋
減衰のもう1つの働きって何ですか?2つあるって言ってましたよね。
🎓
2つめは、もっと目に見える方。減衰は「振幅を減らす」。そのスピードを表すのが対数減衰率 δ なんだ。連続する2つのピーク振幅の比の自然対数、δ = ln(xₙ/xₙ₊₁)。1回揺れるごとに振幅が e^(−δ) 倍になる、と読める。デフォルト値(m=10 kg、k=10000 N/m、c=40 N·s/m)だと δ≈0.398 で、10サイクル後には振幅が最初の約1.9%まで落ちる。上のグラフで指数関数的に縮む包絡線が、まさにこの δ で決まっているよ。
🙋
対数減衰率って、計算で出すだけのものなんですか?それとも実際の役に立つ場面が?
🎓
これがすごく実用的なんだ。δ は「逆算」で真価を発揮する。構造を1回コツンと叩いて、その自由振動が消えていく波形を録る。ピーク振幅がどれくらいの速さで減るかを測れば δ が分かり、そこから減衰比 ζ を逆に求められる。これが「対数減衰率法」で、振動試験の中でも最も手軽で実用的な減衰測定テクニックの1つだよ。叩いて、録って、ピークを数えるだけ。減衰という「見えない物性」を引き出せるんだ。

よくある質問

非減衰固有振動数 fₙ は、減衰がまったくない理想的な系が自由振動するときの振動数で、fₙ = (1/2π)·√(k/m) で決まります。減衰固有振動数 f_d は、実際の減衰がある系が自由振動するときの振動数で、f_d = fₙ·√(1−ζ²) と少しだけ低くなります。ζ は減衰比です。減衰が大きいほど振動はゆっくりになり、ζ = 1(臨界減衰)で f_d はゼロになって振動しなくなります。軽い減衰(ζ が数パーセント)では f_d と fₙ の差はごくわずかで、実務では非減衰の式で十分なことが多いです。
対数減衰率 δ は、減衰自由振動の連続する2つのピーク振幅の比の自然対数で、δ = ln(xₙ/xₙ₊₁) = 2πζ/√(1−ζ²) で定義されます。1サイクルごとに振幅がどれだけ速く減るかを表す無次元量で、δ が大きいほど振動が速く収まります。実務で非常に役立つのは「逆算」で、構造を1回叩いて自由振動波形を録り、ピーク振幅の減り方から δ を測れば、そこから減衰比 ζ を求められます。これは対数減衰率法と呼ばれ、最も基本的な実験的減衰測定法の一つです。
減衰比 ζ が 1 に達すると臨界減衰となり、系は振動せずに最短時間でまっすぐ静止位置へ戻ります。このとき減衰固有振動数 f_d はゼロです。ζ が 1 を超えると過減衰となり、系はさらにゆっくりと、やはり振動せずに戻ります。臨界減衰係数は c_c = 2·√(m·k) で、減衰比は ζ = c/c_c と表せます。ドアクローザや車のサスペンション、計測器の指針などは、行き過ぎ(オーバーシュート)なく素早く整定させたいので、臨界減衰の少し手前(ζ ≈ 0.6〜0.8)に設計されることが多いです。
多くの実構造の減衰比は ζ = 0.01〜0.05 程度と非常に小さく、このとき f_d = fₙ·√(1−ζ²) の補正項 √(1−ζ²) はほぼ 1 になります。例えば ζ = 0.05 でも振動数の低下はわずか約0.1%で、測定誤差より小さいほどです。そのため固有振動数の計算には、計算が簡単で減衰の値を知らなくても求まる非減衰の式 fₙ = (1/2π)·√(k/m) を使うのが実務の定石です。一方、減衰が大きい防振系や制振設計、ζ が 1 に近づくケースでは、減衰固有振動数の低下が無視できなくなるため f_d を正しく扱う必要があります。

実世界での応用

実験モード解析と打撃試験:橋梁・建築・機械部品の動特性を調べるとき、インパクトハンマで構造を1回叩いて自由振動を録り、その減衰波形のピーク振幅の減り方から対数減衰率 δ を求め、減衰比 ζ を算出します。これは対数減衰率法と呼ばれる時間領域の標準的な減衰測定法で、周波数領域の半値幅法と補完関係にあります。鋭い単一モードの構造では、たった1回の打撃と数個のピークの測定だけで減衰比が得られるため、現場で最も手軽に使われる手法です。

制振・防振設計:制振材・ダイナミックダンパ・防振ゴムを追加したとき、狙いどおりに減衰が増えたかを、自由振動の収まり方で確認します。対策前後で打撃応答を録り、振幅の減衰が速くなり対数減衰率 δ が大きくなっていれば制振が効いている証拠です。逆に振動がいつまでも尾を引くなら、ダンパの取り付け位置や減衰量の見直しが必要だと分かります。

臨界減衰を狙った機構設計:ドアクローザ、自動車のショックアブソーバ、アナログ計器の指針、ハードディスクのヘッド機構などは、オーバーシュートなく素早く静止させたいので、臨界減衰の少し手前(減衰比 ζ ≈ 0.6〜0.8)に設計されます。減衰が小さすぎると行き過ぎて何度も揺れ、大きすぎると戻りが遅い。このツールで減衰係数を動かすと、減衰比が 1 に近づくにつれ減衰固有振動数がゼロへ落ちていく様子が確認できます。

機械の状態監視・異常診断:回転機械や構造物では、ボルトのゆるみ・亀裂・支持剛性の変化が減衰特性を変えます。定期的に打撃応答を取り、対数減衰率から減衰比を追跡すると、減衰の急変が構造劣化の早期サインとして検出できます。減衰の変化は固有振動数のずれより敏感に現れることが多く、予知保全のパラメータとして重宝されます。

よくある誤解と注意点

まず多いのが、「減衰固有振動数と非減衰固有振動数を取り違える」誤りです。両者は f_d = fₙ·√(1−ζ²) で結ばれますが、軽減衰(ζ が数パーセント)では補正項 √(1−ζ²) がほぼ 1 で差がごくわずかなため、実務では区別せず非減衰の fₙ を使うことが多くあります。一方、減衰が大きい防振系や ζ が 1 に近づくケースでは差が無視できず、減衰固有振動数を正しく扱う必要があります。さらに紛らわしいのは、自由振動を観測して得られるのは「減衰固有振動数 f_d」だという点です。打撃試験のピーク間隔から求めた振動数は f_d であって fₙ そのものではないため、減衰が大きいときは ζ を使って fₙ に補正する必要があります。

次に、「対数減衰率を隣り合う1組のピークだけで測る」こと。δ は本来 δ = ln(xₙ/xₙ₊₁) と1サイクル分の振幅比で定義されますが、実測ではノイズの影響で1組のピークだけから求めると誤差が大きくなります。より正確には、n サイクル離れた2つのピークを使い δ = (1/n)·ln(x₀/xₙ) とします。減衰が小さいほど1サイクルの振幅減少が小さく、ノイズに埋もれやすいため、十分離れたピークを使うか、複数のピークを対数軸で直線フィットするのが実務の定石です。1組のピークだけで得た減衰比は、しばしば真値から大きく外れます。

最後に、「減衰は必ず固有振動数を下げると思い込む」こと。確かに減衰固有振動数 f_d は非減衰固有振動数 fₙ より低くなりますが、その低下は √(1−ζ²) の項によるもので、軽減衰では実質的にゼロに近いほど小さいものです。一方で、強制振動の「共振周波数」は減衰によってさらに別の量だけずれます。自由振動の減衰固有振動数、強制振動の共振周波数、非減衰固有振動数の3つは厳密には別物で、減衰が大きいときに混同すると設計を誤ります。このツールが扱うのは「自由振動の減衰固有振動数」であり、強制加振の共振とは区別して理解してください。

使い方ガイド

  1. 質量(kg)、ばね定数(N/m)、減衰係数(N·s/m)をスライダーで設定します。鋼製フレーム構造の場合、質量500kg、ばね定数50000N/m、減衰係数2000N·s/mが典型値です。
  2. 非減衰固有振動数ωnと減衰比ζが自動計算されます。ζ < 1で減衰自由振動、ζ ≥ 1で過減衰となります。
  3. 減衰固有振動数ωd、対数減衰率δ、振幅半減サイクル数を確認し、任意のサイクル数nを指定してn回後の振幅残存率(%)を観察します。

具体的な計算例

建築物制震装置で、質量1000kg、ばね定数100000N/m、減衰係数5000N·s/mの1自由度系を考えます。非減衰固有振動数ωn≈1.59Hz、減衰比ζ≈0.79となり、減衰固有振動数ωd≈0.94Hz、対数減衰率δ≈0.436です。振幅は約1.6サイクルで半減し、10サイクル後の残存振幅は初期値の約1.8%に減少します。

実務での注意点