対数減衰率による減衰比測定シミュレーター 戻る
振動工学

対数減衰率による減衰比測定シミュレーター

物体を一度たたいて自由に振動させ、その揺れが収まっていく波形から減衰特性を測るツールです。1回目とn周期後のピーク振幅を入力すると、対数減衰率 δ・減衰比 ζ・Q値がリアルタイムで分かり、減衰の評価ができます。

パラメータ設定
1回目のピーク振幅 A₁
mm
減衰波形で最初に現れる山の高さ
n周期後のピーク振幅 Aₙ
mm
n サイクル経過後に読み取る山の高さ
経過サイクル数 n
周期
A₁ と Aₙ の間に挟まる完全な周期の数
不減衰固有振動数 fₙ
Hz
減衰がないと仮定したときの振動数
計算結果
対数減衰率 δ
減衰比 ζ
減衰固有振動数 fd (Hz)
1サイクルの振幅比
Q値(品質係数)
減衰の分類
自由振動の減衰波形 — 描画アニメーション

減衰サイン波 x(t) と、その包絡線 ±A·e^(−ζωₙt)(破線)を描きます。マーカーが波をなぞり、1回目のピーク A₁ と n 周期後のピーク Aₙ を示します。

自由振動の減衰波形 x(t)
ピーク振幅の減衰(対数縦軸)
理論・主要公式

$$\delta=\frac{1}{n}\ln\frac{A_1}{A_{n}},\qquad \zeta=\frac{\delta}{\sqrt{4\pi^2+\delta^2}}$$

対数減衰率 δ と減衰比 ζ。A₁:1回目のピーク振幅、Aₙ:n周期後のピーク振幅、n:経過サイクル数。

$$f_d=f_n\sqrt{1-\zeta^2},\qquad Q=\frac{1}{2\zeta}$$

減衰固有振動数 fd(実際の振動数)と Q値。fₙ:不減衰固有振動数。

減衰が小さい(ζ < 0.1)ときは ζ ≈ δ/(2π) の近似が成り立ち、δ から減衰比を暗算で見積もれます。

対数減衰率による減衰比測定とは

🙋
「減衰比」って、振動がどれくらい早く収まるかの指標ですよね。でも、それって実際にどうやって測るんですか?式に出てくる減衰係数 c なんて、現物を見ても分からない気がして…。
🎓
いい疑問だね。減衰係数 c を直接測るのは確かに難しい。そこで使うのが「対数減衰率」という考え方なんだ。やることはシンプルで、対象をコツンと一度たたいて自由に揺らす。あとは揺れがだんだん小さくなっていく波形を記録するだけ。波の山の高さが1周期ごとにどれだけ縮むか、それを見れば減衰が分かる、というアイデアだよ。
🙋
山の高さの「縮み具合」ですか。となりの山同士を比べればいいんですか?
🎓
そう。減衰自由振動の包絡線は指数関数 e^(−ζωₙt) で減っていくから、山の高さは1周期ごとに「同じ比率」で小さくなる。隣り合うピークの比 A₁/A₂ の自然対数をとったものが対数減衰率 δ だ。ただ、隣同士だと差が小さくて読み取り誤差が乗りやすい。だから実務では n 周期ぶん離れた山を比べて δ = (1/n)·ln(A₁/Aₙ) とする。左の「経過サイクル数 n」を増やすほど、平均化されて測定が安定するんだよ。
🙋
なるほど、n で割って1周期あたりに直すんですね。その δ から減衰比 ζ はどう出てくるんですか?
🎓
厳密には ζ = δ/√(4π²+δ²) という式で結びついている。でも現場ではもっとざっくり ζ ≈ δ/(2π) と覚えていることが多い。減衰が小さいうちは両者ほぼ同じだからね。例えばデフォルト値だと δ ≈ 0.285、ζ ≈ 0.045。鋼構造物やよく設計された機械はだいたい ζ = 0.01〜0.05 くらいに収まる。逆に ζ が 0.2 を超えるようなら、ゴムやダンパーがしっかり効いている系だと考えていい。
🙋
「Q値」っていうのも結果に出ていますね。これは何を表しているんですか?
🎓
Q値は「共振の鋭さ」を表す数字で、Q = 1/(2ζ) という単純な関係だ。減衰が小さいほど Q は大きく、共振ピークが針のように鋭くなる。ζ = 0.045 なら Q ≈ 11。スピーカーや音叉のように「よく響く」ものは Q が高く、地震対策の制振ダンパーのように「すぐ止めたい」ものは Q を意図的に低くする。同じ減衰特性を、減衰比・対数減衰率・Q値という3つの言葉で言い換えているだけ、と思っておけばいい。
🙋
最後にひとつ。減衰があると振動数も変わるって本当ですか?固有振動数って一定だと思っていました。
🎓
本当だよ。減衰があると振動はほんの少しゆっくりになる。減衰がないと仮定した不減衰固有振動数 fₙ に対して、実際に振動する減衰固有振動数は fd = fₙ·√(1−ζ²) だ。ただし ζ が小さいときは √(1−ζ²) がほぼ1で、差は1%未満。デフォルト値でも fd ≈ 9.99 Hz と、fₙ=10 Hz とほぼ変わらない。だから多くの実務計算では fd ≈ fₙ としてしまう。でも ζ が 0.3、0.5 と大きくなると差が無視できなくなるから、そこは頭の片隅に置いておこう。

よくある質問

対数減衰率 δ は、自由振動の減衰波形で連続するピーク振幅の比の自然対数です。1回目のピーク振幅 A1 と n 周期後のピーク振幅 An を使い、δ = (1/n)·ln(A1/An) で求めます。n 周期にわたって平均することで、1周期だけ読むより読み取り誤差が小さくなります。減衰が大きいほど振幅は速く小さくなり、δ も大きくなります。本ツールは2つのピーク振幅と経過サイクル数からこの δ を計算します。
減衰比 ζ は対数減衰率 δ から ζ = δ / √(4π² + δ²) で厳密に求められます。減衰が小さい(ζ < 0.1 程度)場合は ζ ≈ δ/(2π) という簡便な近似が使えます。例えば δ = 0.285 なら ζ ≈ 0.285/6.283 ≈ 0.0454 で、厳密式とほぼ一致します。減衰が大きくなると分母の δ² の項が無視できなくなるため、本ツールは常に厳密式で計算しています。
不減衰固有振動数 fn は減衰がまったくないと仮定したときの振動数で、減衰固有振動数 fd は実際に減衰系が振動する振動数です。両者は fd = fn·√(1−ζ²) の関係にあり、減衰があると振動は少しだけゆっくりになります。減衰比が小さい(ζ < 0.1)場合、fd と fn の差は1%未満でほぼ無視できますが、ζ が 0.3 を超えると差が目立ち始めます。減衰周期 Td = 1/fd は、減衰波形で隣り合うピークの時間間隔として実測できます。
Q値は系の鋭さを表す無次元量で、減衰比とは Q = 1/(2ζ) の関係にあります。減衰が小さいほど Q は大きく、共振が鋭くなります。例えば ζ = 0.045 なら Q ≈ 11、ζ = 0.01 なら Q = 50 です。Q値は1周期あたりに蓄えられたエネルギーと失われるエネルギーの比とも解釈でき、振動の半値幅から Q = fn/Δf として周波数応答から求めることもできます。対数減衰率からの Q と周波数応答からの Q が一致するかは、測定の妥当性チェックになります。

実世界での応用

建築・土木構造物のモーダル試験:橋梁・高層ビル・煙突などは、風や交通振動、起振機による加振で自由振動させ、収まっていく波形から対数減衰率を求めます。鋼構造物の減衰比はおおむね ζ = 0.005〜0.02、鉄筋コンクリート構造で 0.02〜0.05 が目安です。耐震・耐風設計では、この実測減衰比が応答計算の前提値として直接使われるため、設計値より小さい減衰しか得られないと想定外の大きな揺れにつながります。

機械部品・回転機械の減衰評価:工作機械の主軸、タービンブレード、配管系などでは、インパルスハンマで一度たたいて応答を記録するインパルス加振試験が広く使われます。波形のピーク減衰から対数減衰率を求め、共振時の振幅増幅率(Q値)を見積もることで、危険速度付近の運転がどれだけ危ないかを評価します。Q値が高い(減衰が小さい)部品は、共振に入った瞬間に振幅が跳ね上がります。

制振・防振設計の効果確認:ゴムマウント、オイルダンパ、制振合金などを追加する前後で対数減衰率を測れば、その対策で減衰比がどれだけ増えたかを定量的に示せます。「なんとなく揺れが減った」ではなく「δ が 0.05 から 0.25 に増えた」と数字で言えるのが対数減衰率法の強みです。自動車のサスペンションやエンジンマウントの調整でも同じ手法が使われます。

CAE減衰モデルの検証:有限要素法による振動解析では、レイリー減衰やモード減衰として減衰比を入力する必要があります。この入力値を当てずっぽうで決めると応答計算が大きく外れます。実機をインパルス加振して対数減衰率から ζ を実測し、その値を解析モデルに反映する「モデルキャリブレーション」が、信頼できる振動シミュレーションには不可欠です。

よくある誤解と注意点

まず多いのが、「ピーク振幅をどこから読むかが曖昧」という問題です。対数減衰率は同じ符号・同じ向きのピーク同士で比べる必要があります。プラス側の山とマイナス側の谷を混ぜて数えたり、減衰が大きい系で1周期を半周期と取り違えたりすると、δ が2倍・半分にずれます。さらに、たたいた直後の最初の数周期はインパルスの過渡的な乱れや高次モードの混入が残っていることが多く、ここを A₁ にすると誤差が乗ります。波形が安定したきれいな1サイクルを A₁ に選び、十分に減衰したが計測ノイズに埋もれていない山を Aₙ にするのがコツです。

次に、「対数減衰率は1自由度系を前提にしている」という点を忘れがちです。この方法が成り立つのは、減衰波形がきれいな単一周波数の指数減衰サイン波である場合だけです。複数の固有振動数が近接して同時に励起される多自由度系では、波形に「うなり」が生じてピーク高さが単調に減らず、対数減衰率を素直に適用すると誤った値になります。実機では帯域通過フィルタで対象モードだけを抽出してから δ を求める、あるいは曲線あてはめ(カーブフィット)で指数包絡線を推定するのが正攻法です。本ツールは1自由度・粘性減衰の理想系を扱っていることに注意してください。

最後に、「減衰は振幅によらず一定」だと思い込むことです。本ツールが前提とする粘性減衰では、減衰比は振幅に関係なく一定で、対数減衰率も1周期ごとに同じ値になります。しかし現実の減衰には、ボルト継手のすべりや材料内部摩擦に由来する「クーロン摩擦減衰(一定減衰)」や「履歴減衰」が混ざっていることが多く、これらは振幅が大きいほど減衰の効き方が変わります。摩擦減衰が支配的だと、ピークの減り方が指数関数ではなく直線的になり、振幅が小さくなるほど見かけの対数減衰率が大きくなります。減衰波形が片対数プロットで直線に乗らない場合は、粘性減衰の仮定そのものを疑うべきサインです。

使い方ガイド

  1. 初期振幅a1と最初の計測範囲a1Rangeを入力。鋼構造物や機械装置の自由振動における初期ピーク値(例:10mm)を設定します
  2. n番目の振幅anと計測範囲anRangeを指定。周期経過後の振幅値(例:5周期後に2.5mm)を入力して対数減衰率δを自動計算
  3. 減衰固有振動数fnと周期数ncを設定。減衰比ζ・Q値・1サイクル振幅比が算出され、减衰特性が自動分類されます

具体的な計算例

コンクリート橋脚の耐震診断で、初期ピーク振幅a1=15mm、5周期後の振幅a5=3.2mm、固有振動数fn=2.1Hzの場合、対数減衰率δ≒0.498、減衰比ζ≒7.9%、Q値≒6.3が得られます。1サイクル当たり振幅比≒0.78で、低減衰領域(ζ<10%)と判定。アルミ合金機械部品でfn=8.5Hz、初期10mmから3周期で2.8mmの場合、δ≒0.443、ζ≒6.9%となります

実務での注意点