回転機械の転がり軸受診断方法
【課題】振動センサの共振特性の影響を含んだ信号でも、従来の有線方式のセンサ出力で診断するＱ値と同等の値を確保することのできる転がり軸受診断方法を提供する。
【解決手段】振動センサ１から送信されてきたデジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した振動センサの周波数応答倍率を補正係数αとしてＱ値を補正し、この補正されたＱａ値に基づいて診断を行なう。
【解決手段】振動センサ１から送信されてきたデジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した振動センサの周波数応答倍率を補正係数αとしてＱ値を補正し、この補正されたＱａ値に基づいて診断を行なう。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発電機や電動機など機械装置の軸受け部分に発生する異常振動を機械装置を運転したままの状態で測定し、電波として出力する振動センサを用いた回転機械の転がり軸受診断方法に関する。
【背景技術】
【０００２】
従来、多くの回転機械で構成されるプラントでは、設備の異常による操業停止を防止するため、主要な回転機械の振動監視が行われており、その方法も振動の変化をオンラインで常時監視する方法やオフラインで定期的に監視する方法により、コンピュータ等によるデータ収集装置でFFT(Fast Fourier Transform)等の解析処理が行われている（特許文献1参照）。
回転機械の軸受部で検出される振動には、例えばロータやファンなどの回転部の不釣合いに起因する回転数成分の１倍に相当する周波数成分や、電動機や発電機等ではロータの極数Ｎに起因する回転周波数のＮ倍に相当する周波数成分が計測された加速度信号の中に含まれている。さらに、転がり軸受の場合には、回転周波数ｆと軸受の玉数ｍの整数倍（ｎ*ｆ*ｍ ここでｎ=1,2・・・・）や転動面に発生した傷に起因して励振される高調波成分が含まれるため、振動センサで測定された振動信号の中には、約10Hz〜10kHzの広帯域の振動周波数成分を含んでいる。
この中で、転がり軸受の振動は1kHz〜10kHzの高周波領域の振動加速度となるため、1KHz〜10kHzを通過するバンドパスフィルタ処理後に解析処理が行われる。
圧電型の加速度センサを用いた従来の振動センサでは、加速度センサの信号をアンプ回路で増幅し、そのアナログ信号をＡ／Ｄ変換器によりデジタル信号に変換し、伝送回路によりセンサ周囲に設置した受信機で収集し、収集したデジタル信号を用いて解析処理により診断が行なわれる。
【特許文献１】特開平11-337400号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前記で説明したような振動センサは、検出した加速度信号の増幅、Ａ／Ｄ変換およびメモリ機能をもつアンプ回路基板や送信回路基板、それらを駆動するための電池を内蔵し、さらに送受信用アンテナが取付けられる。また、これらの電気部品は振動センサの取付け環境、例えば粉塵・油・湿気などに対応するため樹脂モールドされる場合が一般的である。
このように、無線方式の振動センサでは、有線方式のセンサと比較して必然的にサイズが大きくなり重くなるため、振動センサの固有振動数を転がり軸受の対象周波数1kHz〜10kHzに影響しない数10KHzにすることが困難になる。
また、固有振動数が10kHzよりも十分に高い市販の振動センサは、重さが数ｇ程度であり、100ｇを超えるような重さ（大きさ）のものはない。このことから分かるように、伝送方式の振動センサ、例えば振動センサや無線基板が組み込まれた一体型のRFIDタグは、重さ（大きさ）の関係で10kHz以内にタグの共振特性による影響を避けることは困難であった。
このため、無線方式の振動センサを用いる場合には、振動センサの固有振動数を転がり軸受の対象周波数の上限以上にしても、10kHzから十分離れていないために共振特性の影響が現れて周波数応答特性が平坦にならず、正確な加速度振幅値が検出できないといった問題が生じる。
【０００４】
本発明は、上記に鑑み、振動センサの共振特性の影響を含んだ信号でも従来の有線方式のセンサ出力で診断するＱ値と同等の値を確保することのできる転がり軸受診断方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記の課題を解決するために、本発明は、回転機械の軸受部に設けた振動センサで振動加速度を検出して電気信号に変換し、該電気信号から所定帯域の高周波数領域のアナログ信号を抽出した後、前記アナログ信号をデジタル化して外部に電波として出力し、前記デジタル化された出力信号を包絡線化処理して波高率に変換し、前記包絡線化処理波形と波高率よりＱ値を算出し、該Ｑ値の大きさにより回転機械の転がり軸受を診断する方法において、前記デジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した補正係数αで前記Ｑ値を補正し、この補正されたＱａ値に基づいて診断を行なうことを特徴とする。
また、上記の発明において、前記補正係数αを振動センサの周波数応答倍率とすることを特徴とする。
また、上記の発明において、前記Ｑ値を次式（１）にて算出し、
【０００６】
[数１]
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
前記Ｑａ値を次式（２）にて算出することを特徴とする。
【０００７】
[数２]
Ｑa＝Ｑ/α2・・・（２）
さらに、上記の各発明において、前記所定帯域の高周波領域は、1ｋHz〜10ｋHzの領域であることを特徴とする。
【発明の効果】
【０００８】
出力信号の主要周波数に対応した応答倍率で補正するため、振動センサの共振特性による応答レベルの増大を削減することができ、固有振動数が10KHｚに近い振動センサでも精度よいQ値を求めることができる。
また、無線方式の振動センサの固有振動数を数10KHzにする必要がないために、比較的形状を大きくすることができ、振動センサの回路基板やバッテリーの大きさに余裕ができる。このため、基板設計の簡素化や電池の長寿命化が可能になる。 さらに、振動センサの形状・構造が簡素化でき低コスト化が可能になる。
【発明を実施するための最良の形態】
【０００９】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明の転がり軸受診断方法を実施するための振動センサおよび診断器の構成を示すブロック図であり、振動センサ１と診断器２の内部構成を示している。振動センサ１は、圧電素子を使用した加速度検出素子１１と、その加速度検出素子１１の電荷出力を電圧に変換するチャージアンプ１２と、チャージアンプ１２の出力信号から例えば1kHz〜10KHzの高周波帯域を通過するバンドパスフィルタ１３と、その出力信号を増幅するアンプ１４と、アンプ１４からのアナログ信号をデジタル信号に変換するためのA/D変換器１５と、変換されたデジタル信号を一旦記憶するデータ記憶部１６と、記憶されたデジタル信号を外部からの指令に基づいて無線送信するための無線送信部１７と、各部に電源を供給するバッテリーからなる電源供給部１８とで構成されている。
また、診断器２は、データ記憶部１６に記憶されているデジタル信号を無線送信部１７を介して伝送させるための無線受信部２１と、受信したデジタル信号の周波数分析を行なって周波数スペクトルを算出する周波数分析処理部２２と、その周波数スペクトルの中の最大値における周波数を検出する最大ピーク周波数検出部２３と、最大ピーク周波数検出部２３で検出された周波数に対する補正値を予め記憶された補正値から抽出、あるいは演算式から算出する補正値設定部２４と、受信したデジタル信号を包絡線の波形の波高率に変換する包絡線化処理部２５と、包絡線波形の自乗平均平方根と尖頭値と、補正値よりＱ値を算出するＱ値算出処理部２６とで構成されている。
【００１０】
このような構成において、軸受部の振動を例えば圧電素子などの加速度検出素子１１により検出し、チャージアンプ１２で電気信号に変換し、バンドパスフィルタ１３で１ｋHz〜10KHzの帯域成分を抽出しアンプ１４で増幅する。増幅された電気信号はＡ／Ｄ変換器１５によりデジタル変換し、データ記憶部１６、無線送信部１７を介して無線受信部２１に向けて電波として出力する。
診断器２は、無線受信部２１でデジタルデータを受信し、その受信信号を周波数分析器２２で周波数スペクトルを算出し、最大ピーク周波数検出部２３で周波数スペクトルの中から最大値を示す周波数を検出し、その検出周波数に対応したＱ値の補正値を補正値設定部２４で設定する。
Ｑ値は、一例として、特開昭56-135140号公報に記載されている次の（１）式により算出される。
【００１１】
[数３]
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
この（１）式の意味は、特開昭56-135140号公報に記載されているが、次の通りである。
すなわち、Ｖrms×Ｖpeakを変形すると、（２）式のように表すことができる。
【００１２】
[数４]
Ｖrms×Ｖpeak＝（Ｖrms）2×（Ｖpeak／Ｖrms）・・・（２）
（２）式の右辺において、（Ｖrms）2はランダム波形のパワーを示しており、通常ランダム波のエネルギーを表すのに用いられる値である。一方、（Ｖpeak／Ｖrms）はクレストファクタリーまたは波形率と称して、波形のおおよその形の定義として用いられるものである。ここで、（Ｖrms）2の値について考えてみると、軸受にキズ、うねり、電食等の損傷がある場合にはランダム波形は大きくなるため、当然のことながら（Ｖrms）2の値は大きくなる。
これに対して、（Ｖpeak／Ｖrms）の値はランダム波よりもピーク値が大きく出るような波形の場合、すなわち、軸受にキズがある場合には５よりも大きくなるため、図２において（Ｖpeak／Ｖrms）＝５の線より上方に位置する。
なお、図２は、（Ｖpeak／Ｖrms）、Ｖrms×Ｖpeakを演算し、縦軸に（Ｖpeak／Ｖrms）、横軸にＶrms×Ｖpeakをとって両対数グラフ上にプロットしたものである。図２からも明らかなとおり、（Ｖpeak／Ｖrms）＝５に集中しており、軸受を分解した結果、まだ十分使用に耐えるものは○印、かなり損傷があるが、まだ使用にたえなくもないものは△印、損傷がひどく使用できないものは×印でそれぞれプロットされている。
【００１３】
また、ランダム波よりも低周波成分の多い波形の場合、すなわち、軸受にうねり、電食等がある場合には５よりも小さくなるため、図２において、（Ｖpeak／Ｖrms）の値を掛け算することにより、（Ｖpeak／Ｖrms）＝５を掛け算する場合に比べて（Ｖrms）2×（Ｖpeak／Ｖrms）の値は小さくなる。しかしながら、損傷がある場合には、（Ｖrms）2の値のほうが大きな変化を示すこと、損傷としてはうねりや電食よりもキズのほうが重大であることを考えれば、Ｖrms×Ｖpeakの値により良否の判別を行なうことができるのである。
次に（１）式の右辺において、[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]について説明する。前述において、Ｖrms×Ｖpeakの値により良否の判別を行なうことができることを説明したが、（Ｖrms）2の値に対して（Ｖpeak／Ｖrms）の値の変化のほうが小さいため、Ｖrms×Ｖpeakの値だけの良否の判別では（Ｖpeak／Ｖrms）＝５の線から離れた損傷に対しては不十分なことがある。
そのため、[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]なる補正係数を掛け算し、Ｖrms×Ｖpeakの値を補正してＱ値を算出しているのである。
この補正の効果は、図２を用いて説明すると、Ｖrms×Ｖpeakの値にて良否の判別を行なう場合には、破線ｂが境界となりＶrms×Ｖpeakの値が小さいが、（Ｖpeak／Ｖrms）の値が５より離れているものも含むことができるようになる。なお、軸受の種類によって定数ａ，ｂの値を設定することにより、図２において境界ｂの傾きを変えることができる。
【００１４】
本発明においては、このＱ値を補正値設定部２４における補正値αで補正する。この補正値αは、最大ピーク周波数検出部２３で検出した周波数に基づいて、図３に示すような振動センサの周波数応答特性による応答倍率の数値データを記憶したテーブルから選択したり、応答倍率の近似式から求めたりすることができる。
近似式を用いる場合には、一例として次の（３）式を用いることができる。
【００１５】
[数５]
α＝ａ×Ｆmax4−ｂ×Ｆmax3＋ｃ×Ｆmax2−ｄ×Ｆmax ＋ｅ・・・（３）
但し、ａ〜ｅは近似のための係数
この（３）式による求められる応答倍率（補正値α）は、各係数ａ〜ｅを、例えばａ＝0.0122，ｂ＝0.1872，ｃ＝1.0779，ｄ＝2.3607，ｅ＝2.5442とした場合には、図４に破線で示すとおりであり、実線で記載されている図３の応答倍率とほぼ一致する。
この補正値αを用いての補正は、ＶrmsとＶpeakとをそれぞれαで割るため、結果的にＱ値をαの２乗で割ることと同じになる。従って、補正後の軸受判定値Ｑaは、（４）式で示される。
【００１６】
[数６]
Ｑa＝Ｑ/α2・・・（４）
この軸受判定値Ｑaに応じて、図５に示すように良好、注意、危険の判定を行なう。図５は、転動体（軸受）の公転速度の指標であるＤｍ×Ｎに対してのＱａ値の判定図表を示したものである。但し、Ｄｍは転動体の公転直径（ｍｍ）、Ｎは機械の回転速度（ｒ．ｐ．ｍ）である。
全周波数帯域をまずフーリエ変換し、各周波数スペクトル毎に振幅補正し、その後に逆フーリエ変換して波形に戻し、波形のピーク値や実行値を求めてＱ値を演算することもできるが、ソフトウェアが複雑となり、コストが高くなる。
しかし、軸受の傷による振動は、軸受部の固有振動数で励振されるため、主要成分が特定の周波数成分になることから、本発明のようにピーク周波数成分で補正することが有効になる。
【図面の簡単な説明】
【００１７】
【図１】本発明による転がり軸受の診断方法の処理を説明するためのブロック図。
【図２】軸受を分解点検した結果をプロットした図。
【図３】振動センサの周波数応答特性による応答倍率を示す図。
【図４】近似式により求めた振動センサの周波数応答特性による応答倍率を示す図。
【図５】Ｑａ値の判定図表。
【符号の説明】
【００１８】
１・・・・振動センサ
２・・・・診断器
１１・・・加速度検出素子
１５・・・Ｑ値算出処理
１２・・・チャージアンプ
１３・・・バンドローパスフィルタ
１５・・・Ａ／Ｄ変換器
１６・・・データ記憶部
１７・・・無線送信部
２１・・・無線受信部
２２・・・周波数分析処理部
２３・・・最大ピーク周波数検出部
２４・・・補正値設定部
２５・・・包絡線化処理部
２６・・・Ｑ値算出処理部
【技術分野】
【０００１】
本発明は、発電機や電動機など機械装置の軸受け部分に発生する異常振動を機械装置を運転したままの状態で測定し、電波として出力する振動センサを用いた回転機械の転がり軸受診断方法に関する。
【背景技術】
【０００２】
従来、多くの回転機械で構成されるプラントでは、設備の異常による操業停止を防止するため、主要な回転機械の振動監視が行われており、その方法も振動の変化をオンラインで常時監視する方法やオフラインで定期的に監視する方法により、コンピュータ等によるデータ収集装置でFFT(Fast Fourier Transform)等の解析処理が行われている（特許文献1参照）。
回転機械の軸受部で検出される振動には、例えばロータやファンなどの回転部の不釣合いに起因する回転数成分の１倍に相当する周波数成分や、電動機や発電機等ではロータの極数Ｎに起因する回転周波数のＮ倍に相当する周波数成分が計測された加速度信号の中に含まれている。さらに、転がり軸受の場合には、回転周波数ｆと軸受の玉数ｍの整数倍（ｎ*ｆ*ｍ ここでｎ=1,2・・・・）や転動面に発生した傷に起因して励振される高調波成分が含まれるため、振動センサで測定された振動信号の中には、約10Hz〜10kHzの広帯域の振動周波数成分を含んでいる。
この中で、転がり軸受の振動は1kHz〜10kHzの高周波領域の振動加速度となるため、1KHz〜10kHzを通過するバンドパスフィルタ処理後に解析処理が行われる。
圧電型の加速度センサを用いた従来の振動センサでは、加速度センサの信号をアンプ回路で増幅し、そのアナログ信号をＡ／Ｄ変換器によりデジタル信号に変換し、伝送回路によりセンサ周囲に設置した受信機で収集し、収集したデジタル信号を用いて解析処理により診断が行なわれる。
【特許文献１】特開平11-337400号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前記で説明したような振動センサは、検出した加速度信号の増幅、Ａ／Ｄ変換およびメモリ機能をもつアンプ回路基板や送信回路基板、それらを駆動するための電池を内蔵し、さらに送受信用アンテナが取付けられる。また、これらの電気部品は振動センサの取付け環境、例えば粉塵・油・湿気などに対応するため樹脂モールドされる場合が一般的である。
このように、無線方式の振動センサでは、有線方式のセンサと比較して必然的にサイズが大きくなり重くなるため、振動センサの固有振動数を転がり軸受の対象周波数1kHz〜10kHzに影響しない数10KHzにすることが困難になる。
また、固有振動数が10kHzよりも十分に高い市販の振動センサは、重さが数ｇ程度であり、100ｇを超えるような重さ（大きさ）のものはない。このことから分かるように、伝送方式の振動センサ、例えば振動センサや無線基板が組み込まれた一体型のRFIDタグは、重さ（大きさ）の関係で10kHz以内にタグの共振特性による影響を避けることは困難であった。
このため、無線方式の振動センサを用いる場合には、振動センサの固有振動数を転がり軸受の対象周波数の上限以上にしても、10kHzから十分離れていないために共振特性の影響が現れて周波数応答特性が平坦にならず、正確な加速度振幅値が検出できないといった問題が生じる。
【０００４】
本発明は、上記に鑑み、振動センサの共振特性の影響を含んだ信号でも従来の有線方式のセンサ出力で診断するＱ値と同等の値を確保することのできる転がり軸受診断方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記の課題を解決するために、本発明は、回転機械の軸受部に設けた振動センサで振動加速度を検出して電気信号に変換し、該電気信号から所定帯域の高周波数領域のアナログ信号を抽出した後、前記アナログ信号をデジタル化して外部に電波として出力し、前記デジタル化された出力信号を包絡線化処理して波高率に変換し、前記包絡線化処理波形と波高率よりＱ値を算出し、該Ｑ値の大きさにより回転機械の転がり軸受を診断する方法において、前記デジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した補正係数αで前記Ｑ値を補正し、この補正されたＱａ値に基づいて診断を行なうことを特徴とする。
また、上記の発明において、前記補正係数αを振動センサの周波数応答倍率とすることを特徴とする。
また、上記の発明において、前記Ｑ値を次式（１）にて算出し、
【０００６】
[数１]
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
前記Ｑａ値を次式（２）にて算出することを特徴とする。
【０００７】
[数２]
Ｑa＝Ｑ/α2・・・（２）
さらに、上記の各発明において、前記所定帯域の高周波領域は、1ｋHz〜10ｋHzの領域であることを特徴とする。
【発明の効果】
【０００８】
出力信号の主要周波数に対応した応答倍率で補正するため、振動センサの共振特性による応答レベルの増大を削減することができ、固有振動数が10KHｚに近い振動センサでも精度よいQ値を求めることができる。
また、無線方式の振動センサの固有振動数を数10KHzにする必要がないために、比較的形状を大きくすることができ、振動センサの回路基板やバッテリーの大きさに余裕ができる。このため、基板設計の簡素化や電池の長寿命化が可能になる。 さらに、振動センサの形状・構造が簡素化でき低コスト化が可能になる。
【発明を実施するための最良の形態】
【０００９】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明の転がり軸受診断方法を実施するための振動センサおよび診断器の構成を示すブロック図であり、振動センサ１と診断器２の内部構成を示している。振動センサ１は、圧電素子を使用した加速度検出素子１１と、その加速度検出素子１１の電荷出力を電圧に変換するチャージアンプ１２と、チャージアンプ１２の出力信号から例えば1kHz〜10KHzの高周波帯域を通過するバンドパスフィルタ１３と、その出力信号を増幅するアンプ１４と、アンプ１４からのアナログ信号をデジタル信号に変換するためのA/D変換器１５と、変換されたデジタル信号を一旦記憶するデータ記憶部１６と、記憶されたデジタル信号を外部からの指令に基づいて無線送信するための無線送信部１７と、各部に電源を供給するバッテリーからなる電源供給部１８とで構成されている。
また、診断器２は、データ記憶部１６に記憶されているデジタル信号を無線送信部１７を介して伝送させるための無線受信部２１と、受信したデジタル信号の周波数分析を行なって周波数スペクトルを算出する周波数分析処理部２２と、その周波数スペクトルの中の最大値における周波数を検出する最大ピーク周波数検出部２３と、最大ピーク周波数検出部２３で検出された周波数に対する補正値を予め記憶された補正値から抽出、あるいは演算式から算出する補正値設定部２４と、受信したデジタル信号を包絡線の波形の波高率に変換する包絡線化処理部２５と、包絡線波形の自乗平均平方根と尖頭値と、補正値よりＱ値を算出するＱ値算出処理部２６とで構成されている。
【００１０】
このような構成において、軸受部の振動を例えば圧電素子などの加速度検出素子１１により検出し、チャージアンプ１２で電気信号に変換し、バンドパスフィルタ１３で１ｋHz〜10KHzの帯域成分を抽出しアンプ１４で増幅する。増幅された電気信号はＡ／Ｄ変換器１５によりデジタル変換し、データ記憶部１６、無線送信部１７を介して無線受信部２１に向けて電波として出力する。
診断器２は、無線受信部２１でデジタルデータを受信し、その受信信号を周波数分析器２２で周波数スペクトルを算出し、最大ピーク周波数検出部２３で周波数スペクトルの中から最大値を示す周波数を検出し、その検出周波数に対応したＱ値の補正値を補正値設定部２４で設定する。
Ｑ値は、一例として、特開昭56-135140号公報に記載されている次の（１）式により算出される。
【００１１】
[数３]
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
この（１）式の意味は、特開昭56-135140号公報に記載されているが、次の通りである。
すなわち、Ｖrms×Ｖpeakを変形すると、（２）式のように表すことができる。
【００１２】
[数４]
Ｖrms×Ｖpeak＝（Ｖrms）2×（Ｖpeak／Ｖrms）・・・（２）
（２）式の右辺において、（Ｖrms）2はランダム波形のパワーを示しており、通常ランダム波のエネルギーを表すのに用いられる値である。一方、（Ｖpeak／Ｖrms）はクレストファクタリーまたは波形率と称して、波形のおおよその形の定義として用いられるものである。ここで、（Ｖrms）2の値について考えてみると、軸受にキズ、うねり、電食等の損傷がある場合にはランダム波形は大きくなるため、当然のことながら（Ｖrms）2の値は大きくなる。
これに対して、（Ｖpeak／Ｖrms）の値はランダム波よりもピーク値が大きく出るような波形の場合、すなわち、軸受にキズがある場合には５よりも大きくなるため、図２において（Ｖpeak／Ｖrms）＝５の線より上方に位置する。
なお、図２は、（Ｖpeak／Ｖrms）、Ｖrms×Ｖpeakを演算し、縦軸に（Ｖpeak／Ｖrms）、横軸にＶrms×Ｖpeakをとって両対数グラフ上にプロットしたものである。図２からも明らかなとおり、（Ｖpeak／Ｖrms）＝５に集中しており、軸受を分解した結果、まだ十分使用に耐えるものは○印、かなり損傷があるが、まだ使用にたえなくもないものは△印、損傷がひどく使用できないものは×印でそれぞれプロットされている。
【００１３】
また、ランダム波よりも低周波成分の多い波形の場合、すなわち、軸受にうねり、電食等がある場合には５よりも小さくなるため、図２において、（Ｖpeak／Ｖrms）の値を掛け算することにより、（Ｖpeak／Ｖrms）＝５を掛け算する場合に比べて（Ｖrms）2×（Ｖpeak／Ｖrms）の値は小さくなる。しかしながら、損傷がある場合には、（Ｖrms）2の値のほうが大きな変化を示すこと、損傷としてはうねりや電食よりもキズのほうが重大であることを考えれば、Ｖrms×Ｖpeakの値により良否の判別を行なうことができるのである。
次に（１）式の右辺において、[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]について説明する。前述において、Ｖrms×Ｖpeakの値により良否の判別を行なうことができることを説明したが、（Ｖrms）2の値に対して（Ｖpeak／Ｖrms）の値の変化のほうが小さいため、Ｖrms×Ｖpeakの値だけの良否の判別では（Ｖpeak／Ｖrms）＝５の線から離れた損傷に対しては不十分なことがある。
そのため、[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]なる補正係数を掛け算し、Ｖrms×Ｖpeakの値を補正してＱ値を算出しているのである。
この補正の効果は、図２を用いて説明すると、Ｖrms×Ｖpeakの値にて良否の判別を行なう場合には、破線ｂが境界となりＶrms×Ｖpeakの値が小さいが、（Ｖpeak／Ｖrms）の値が５より離れているものも含むことができるようになる。なお、軸受の種類によって定数ａ，ｂの値を設定することにより、図２において境界ｂの傾きを変えることができる。
【００１４】
本発明においては、このＱ値を補正値設定部２４における補正値αで補正する。この補正値αは、最大ピーク周波数検出部２３で検出した周波数に基づいて、図３に示すような振動センサの周波数応答特性による応答倍率の数値データを記憶したテーブルから選択したり、応答倍率の近似式から求めたりすることができる。
近似式を用いる場合には、一例として次の（３）式を用いることができる。
【００１５】
[数５]
α＝ａ×Ｆmax4−ｂ×Ｆmax3＋ｃ×Ｆmax2−ｄ×Ｆmax ＋ｅ・・・（３）
但し、ａ〜ｅは近似のための係数
この（３）式による求められる応答倍率（補正値α）は、各係数ａ〜ｅを、例えばａ＝0.0122，ｂ＝0.1872，ｃ＝1.0779，ｄ＝2.3607，ｅ＝2.5442とした場合には、図４に破線で示すとおりであり、実線で記載されている図３の応答倍率とほぼ一致する。
この補正値αを用いての補正は、ＶrmsとＶpeakとをそれぞれαで割るため、結果的にＱ値をαの２乗で割ることと同じになる。従って、補正後の軸受判定値Ｑaは、（４）式で示される。
【００１６】
[数６]
Ｑa＝Ｑ/α2・・・（４）
この軸受判定値Ｑaに応じて、図５に示すように良好、注意、危険の判定を行なう。図５は、転動体（軸受）の公転速度の指標であるＤｍ×Ｎに対してのＱａ値の判定図表を示したものである。但し、Ｄｍは転動体の公転直径（ｍｍ）、Ｎは機械の回転速度（ｒ．ｐ．ｍ）である。
全周波数帯域をまずフーリエ変換し、各周波数スペクトル毎に振幅補正し、その後に逆フーリエ変換して波形に戻し、波形のピーク値や実行値を求めてＱ値を演算することもできるが、ソフトウェアが複雑となり、コストが高くなる。
しかし、軸受の傷による振動は、軸受部の固有振動数で励振されるため、主要成分が特定の周波数成分になることから、本発明のようにピーク周波数成分で補正することが有効になる。
【図面の簡単な説明】
【００１７】
【図１】本発明による転がり軸受の診断方法の処理を説明するためのブロック図。
【図２】軸受を分解点検した結果をプロットした図。
【図３】振動センサの周波数応答特性による応答倍率を示す図。
【図４】近似式により求めた振動センサの周波数応答特性による応答倍率を示す図。
【図５】Ｑａ値の判定図表。
【符号の説明】
【００１８】
１・・・・振動センサ
２・・・・診断器
１１・・・加速度検出素子
１５・・・Ｑ値算出処理
１２・・・チャージアンプ
１３・・・バンドローパスフィルタ
１５・・・Ａ／Ｄ変換器
１６・・・データ記憶部
１７・・・無線送信部
２１・・・無線受信部
２２・・・周波数分析処理部
２３・・・最大ピーク周波数検出部
２４・・・補正値設定部
２５・・・包絡線化処理部
２６・・・Ｑ値算出処理部
【特許請求の範囲】
【請求項１】
回転機械の軸受部に設けた振動センサで振動加速度を検出して電気信号に変換し、該電気信号から所定帯域の高周波数領域のアナログ信号を抽出した後、前記アナログ信号をデジタル化して外部に電波として出力し、前記デジタル化された出力信号を包絡線化処理して波高率に変換し、前記包絡線化処理波形と波高率とよりＱ値を算出し、該Ｑ値の大きさにより回転機械の転がり軸受を診断する方法において、
前記デジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した補正係数αで前記Ｑ値を補正し、この補正されたＱａ値に基づいて診断を行なうことを特徴とする回転機械の転がり軸受診断方法。
【請求項２】
請求項１に記載の回転機械の転がり軸受診断方法において、前記補正係数αを振動センサの周波数応答倍率とすることを特徴とする回転機械の転がり軸受診断方法。
【請求項３】
請求項２に記載の回転機械の転がり軸受診断方法において、前記Ｑ値を次式（１）にて算出し、
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
前記Ｑａ値を次式（２）にて算出することを特徴とする回転機械の転がり軸受診断方法。
Ｑa＝Ｑ/α2・・・（２）
【請求項４】
請求項１ないし３の何れか１項に記載の回転機械の転がり軸受診断方法において、前記所定帯域の高周波領域は、1ｋHz〜10ｋHzの領域であることを特徴とする回転機械の転がり軸受診断方法。
【請求項１】
回転機械の軸受部に設けた振動センサで振動加速度を検出して電気信号に変換し、該電気信号から所定帯域の高周波数領域のアナログ信号を抽出した後、前記アナログ信号をデジタル化して外部に電波として出力し、前記デジタル化された出力信号を包絡線化処理して波高率に変換し、前記包絡線化処理波形と波高率とよりＱ値を算出し、該Ｑ値の大きさにより回転機械の転がり軸受を診断する方法において、
前記デジタル化された出力信号を分析して周波数スペクトルを求め、その最大ピーク周波数に対応した補正係数αで前記Ｑ値を補正し、この補正されたＱａ値に基づいて診断を行なうことを特徴とする回転機械の転がり軸受診断方法。
【請求項２】
請求項１に記載の回転機械の転がり軸受診断方法において、前記補正係数αを振動センサの周波数応答倍率とすることを特徴とする回転機械の転がり軸受診断方法。
【請求項３】
請求項２に記載の回転機械の転がり軸受診断方法において、前記Ｑ値を次式（１）にて算出し、
Ｑ＝Ｖrms×Ｖpeak×[a｜（Ｖpeak／Ｖrms）-５｜＋ｂ]・・・（１）
但し、Ｖrms ：包絡線波形の自乗平均平方根値
Ｖpeak：包絡線波形のピーク値
ａ、ｂ：定数
前記Ｑａ値を次式（２）にて算出することを特徴とする回転機械の転がり軸受診断方法。
Ｑa＝Ｑ/α2・・・（２）
【請求項４】
請求項１ないし３の何れか１項に記載の回転機械の転がり軸受診断方法において、前記所定帯域の高周波領域は、1ｋHz〜10ｋHzの領域であることを特徴とする回転機械の転がり軸受診断方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２００９−２２９０８９（Ｐ２００９−２２９０８９Ａ）
【公開日】平成２１年１０月８日（２００９．１０．８）
【国際特許分類】
【出願番号】特願２００８−７１３４３（Ｐ２００８−７１３４３）
【出願日】平成２０年３月１９日（２００８．３．１９）
【出願人】（５９１０８３２４４）富士電機システムズ株式会社 (1,717)
【Ｆターム（参考）】
【公開日】平成２１年１０月８日（２００９．１０．８）
【国際特許分類】
【出願日】平成２０年３月１９日（２００８．３．１９）
【出願人】（５９１０８３２４４）富士電機システムズ株式会社 (1,717)
【Ｆターム（参考）】
[ Back to top ]