基礎からの周波数分析（25）－「振動計測の基礎－4」

これまで、3 回にわたって「振動計測の基礎」に関してお話してきました。今回はその 4 として、減衰比についてのお話となります。

減衰比ζは、振動の挙動に重要な役目を果たしており、自由振動（1 回目）では振動するか否かまたは振動した場合での振動が止まる速さに影響を与え、強制振動(2 回目)では共振周波数におけるピークの大きさに影響を与え、振動伝達率（3 回目）では共振周波数におけるピークの大きさや共振周波数を超えた周波数帯域での伝達率に影響を与えます。
今回（4 回目）では、FFTアナライザで実際に振動計測を行いその結果から減衰比ζを求める方法について説明します。

減衰を表現するパラメータは応用分野によっていろいろありますが、まず、対数減衰率から求める方法を紹介します。
対数減衰率δは、自由振動の減衰波形の隣り合う振幅比の自然対数として定義されます。
なぜ、自然対数なのか？ それはもちろん、図1 にあるように、自由振動の振幅は指数関数的に減衰するからです。

周期T

図1 減衰自由振動波形

図1 の減衰自由振動波形 x (t)は、

となります。
図1 において、対数減衰率δは

ここで、ζ<<1 ならば

次に、対数減衰率δの求め方ですが、まず、ヒルベルト変換の技術を利用して減衰波形の包絡線を求め、横軸を時間、縦軸は振幅の対数表示のグラフを描きます。(図2)

対数減衰率 δ は、減衰周期 T_d に対する傾きだから、dB 表示が常用対数だということに注意して、図2 の直線の傾きより、

ここで、 だから

FFT アナライザ(DS-3000 シリーズ)のデルタサーチ機能を使って、(5)式から、対数減衰率δを算出します。
さらに、(3)式から

実測したデータと(5)式から対数減衰率δを算出して、(6)式から減衰比ζを求めることが出来ます。
図 3 は、ハンマリングして自由振動の減衰時間波形を計測して、それから対数減衰率を求めさらに減衰比を求めた例です。この図は、ヒルベルト変換の帯域制限機能を利用して、1つの共振周波数の減衰波形にしていることに注意して下さい。

File(ハンマリング　時間_0008.txt):ﾋﾙﾍﾞﾙﾄ変換データ Mag

File(ハンマリング　時間_0008.txt):時間軸波形 Real

Log.d	：	対数減衰率
Damp	：	減衰比
Freq	：	減衰固有振動数

図3 減衰時間波形と包絡線データから対数減衰率と減衰比を求めた例

次に、周波数応答関数の共振ピークから半値幅法により、減衰比を求めます。

周波数(Hz)

図4 半値幅法を説明する図

1 自由度系の周波数応答関数(コンプライアンス)G(f )とすると、

ここで、 f_n :固有周波数　　ζ :減衰比

図4において、共振周波数 fn の値 G(f)max から－3 dB(パワーで 1/2、振幅で 1/√2)となる点を f₁ と f₂、またΔf=f₂-f₁とします。
f₁ と f₂は近代的に f_n ± ∆f / 2 とおけるので、 f_nとの比は

この値を (7) 式に代入して逆数をとると、

また、最大値は近代的に (7) 式よりとなるから (9) 式に代入すると

(10) 式を変形して

パワーが 1/2 の帯域を使っているので、上記の(11)式で減衰比を求める方法を半値幅法(half-power bandwidth method)と呼びます。
半値幅法以外に、mdB 法も考えることが出来ます。この時は、

ここで、

ちなみに、1 dB の時は、k=1.97　2 dB の時は、k=1.31 です。

半値幅法で、制振材料の評価パラメータで用いられる損失係数(loss factor)も求められます。履歴減衰を伴う減衰振動では、振動変位に比例する復元力を生じさせる複素弾性率 K’を定義して、実数部を K₁、虚数部を K₂とすると

ここで、

となります。(図5) このηを損失係数と呼びます。複素弾性率 K’は、1 自由度系モデルでばね定数 k と粘性減衰係数 c の両特性を含んだものです。図 5 からわかるように、損失係数は、タンデルタ(タンデル)とも呼ばれます。

半値幅法から、損失係数ηは以下の式で求められます。

CH1-2 周波数応答関数 Mag

Damp	：	減衰比
Lose F	：	損失係数
Freq	：	減衰固有振動数

図6 半値幅法による減衰比と損失係数の実測例

図 6 は、FFT アナライザを使ってハンマリング法で計測した FRF(周波数応答関数)から、半値幅法により減衰比と損失係数を求めた例です。

いろんなパラメータのまとめとして、減衰特性を表す各種パラメータ間の関係の一覧表が、下記の表 1 です。

表1 各種パラメータ間の相互関係

注意1：本表は、ζ ≪ 1  とみなせる場合です。
注意2：Q値は、振動系が共振する場合、この共振の鋭さを表す量で、損失係数の逆数です。

【キーワード】

減衰比、自由振動、強制振動、共振周波数、振動伝達率、対数減衰率、固有角振動数、減衰固有角振動数、減衰周期、ヒルベルト変換、包絡線、デルタサーチ機能、半値幅法、half‐power band width method、損失係数、loss factor、複素弾性率、ばね定数、粘性減衰係数、タンデルタ、タンデル、ハンマリング法、Q 値

【参考】

「モード解析入門」長松昭男著 コロナ社(1994 年)
小野測器技術レポート「減衰特性をあらわす係数」

(2016年1月21日発行メールマガジンより抜粋)