JP7146683B2 - 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム - Google Patents
転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム Download PDFInfo
- Publication number
- JP7146683B2 JP7146683B2 JP2019061225A JP2019061225A JP7146683B2 JP 7146683 B2 JP7146683 B2 JP 7146683B2 JP 2019061225 A JP2019061225 A JP 2019061225A JP 2019061225 A JP2019061225 A JP 2019061225A JP 7146683 B2 JP7146683 B2 JP 7146683B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- specific
- vibration
- rotation
- rolling bearing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、工作機械等に用いられて主軸等の回転体を支持する転がり軸受の異常を診断する方法及び装置、プログラムに関するものである。
回転体を支持する転がり軸受に内輪の損傷などの異常が生じると振動が発生する。軸受の異常によって発生する力は単純な正弦波状ではないため、高調波の周波数成分の振動が同時に観測される。この際に発生する振動の周波数（特徴周波数）は、回転周波数に比例しており、回転体の回転周波数と軸受諸元から算出することが可能である。
例えば特許文献１では、振動加速度を測定してエンベロープ処理および周波数分析等を行った波形を、所定の関係式に基づいて算出した内輪、外輪、転動体などの特徴周波数とともに表示し、特徴周波数に対応する振動ピークがあるか否かを判断できるようにすることで、転がり軸受の損傷部位（内輪、外輪、転動体など）を特定する手法が示されている。
特許文献２では、振動を測定してエンベロープ処理および周波数分析し、所定の関係式に基づいて算出した内輪、外輪、転動体のそれぞれの特徴周波数の値を抽出して、損傷部位（内輪、外輪、転動体）でそれぞれ異なるしきい値を用いて診断する手法が示されている。
例えば特許文献１では、振動加速度を測定してエンベロープ処理および周波数分析等を行った波形を、所定の関係式に基づいて算出した内輪、外輪、転動体などの特徴周波数とともに表示し、特徴周波数に対応する振動ピークがあるか否かを判断できるようにすることで、転がり軸受の損傷部位（内輪、外輪、転動体など）を特定する手法が示されている。
特許文献２では、振動を測定してエンベロープ処理および周波数分析し、所定の関係式に基づいて算出した内輪、外輪、転動体のそれぞれの特徴周波数の値を抽出して、損傷部位（内輪、外輪、転動体）でそれぞれ異なるしきい値を用いて診断する手法が示されている。
しかしながら、損傷部位の特定や損傷部位に応じたしきい値を用いるためには、基本周波数や特徴周波数を算出するための軸受諸元が必要となる。このため、軸受メーカの販売する軸受診断装置においては、その軸受メーカが軸受諸元を把握している自社の軸受については型式を入力することで診断ができるが、他社の軸受については軸受諸元の手入力を要求する実装となっていることが一般的である。この場合、軸受の諸元が入手できない場合には軸受を診断することができないといった課題がある。また、機械の製造者を除くと、その機械で使用されている軸受の型式を把握することさえ困難であるという課題がある。
そこで、本発明は、上記問題に鑑みなされたものであって、診断対象の機械で使用されている軸受の軸受諸元が把握できない場合であっても、軸受の損傷の種類を推定することができる転がり軸受の異常診断方法及び異常診断装置、異常診断プログラムを提供することを目的とする。
上記目的を達成するために、請求項１に記載の発明は、回転体を支持する転がり軸受の異常を診断する方法であって、
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別ステップと、を実行することを特徴とする。
請求項２に記載の発明は、請求項１の構成において、前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする。
請求項３に記載の発明は、請求項１の構成において、前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする。
請求項４に記載の発明は、請求項１乃至３の何れかの構成において、前記軸受損傷種類判別ステップは、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定する処理を含むものであり、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする。
上記目的を達成するために、請求項５に記載の発明は、回転体を支持する転がり軸受の異常を診断する方法であって、
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、前記回転同期成分強調波形を入力して前記転がり軸受の損傷の種類を判別する軸受損傷種類判別ステップと、を実行することを特徴とする。
請求項６に記載の発明は、請求項１乃至５の何れかの構成において、前記振動測定ステップでは、異なる複数の回転周波数で振動を測定することを特徴とする。
請求項７に記載の発明は、請求項１乃至６の何れかの構成において、前記回転同期成分強調波形算出ステップでは、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出する際に、ある比周波数の振動の大きさとして、その比周波数の前後に所定の幅を持たせて当該幅内の振動の大きさの最大値を採用することを特徴とする。
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別ステップと、を実行することを特徴とする。
請求項２に記載の発明は、請求項１の構成において、前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする。
請求項３に記載の発明は、請求項１の構成において、前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする。
請求項４に記載の発明は、請求項１乃至３の何れかの構成において、前記軸受損傷種類判別ステップは、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定する処理を含むものであり、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする。
上記目的を達成するために、請求項５に記載の発明は、回転体を支持する転がり軸受の異常を診断する方法であって、
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、前記回転同期成分強調波形を入力して前記転がり軸受の損傷の種類を判別する軸受損傷種類判別ステップと、を実行することを特徴とする。
請求項６に記載の発明は、請求項１乃至５の何れかの構成において、前記振動測定ステップでは、異なる複数の回転周波数で振動を測定することを特徴とする。
請求項７に記載の発明は、請求項１乃至６の何れかの構成において、前記回転同期成分強調波形算出ステップでは、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出する際に、ある比周波数の振動の大きさとして、その比周波数の前後に所定の幅を持たせて当該幅内の振動の大きさの最大値を採用することを特徴とする。
上記目的を達成するために、請求項８に記載の発明は、回転体を支持する転がり軸受の異常を診断する装置であって、
前記回転体の振動を複数回測定する振動測定手段と、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出手段と、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別手段と、を備えることを特徴とする。
請求項９に記載の発明は、請求項８の構成において、前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする。
請求項１０に記載の発明は、請求項８の構成において、前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする。
請求項１１に記載の発明は、請求項８乃至１０の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定すると共に、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする。
請求項１２に記載の発明は、請求項８乃至１１の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、前記比基本周波数と前記比振動源回転周波数とを合わせて表示する表示手段を備えることを特徴とする。
請求項１３に記載の発明は、請求項８乃至１１の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、前記特定の比周波数の位置と前記回転同期成分強調波形とを合わせて表示する表示手段を備えることを特徴とする。
上記目的を達成するために、請求項１４に記載の発明は、転がり軸受の異常診断プログラムであって、
所定の回転周波数で測定された回転体の振動が前記回転周波数と共に入力されたコンピュータに、請求項１乃至７の何れかに記載の転がり軸受の異常診断方法における回転同期成分強調波形算出ステップと軸受損傷種類判別ステップとを実行させることを特徴とする。
前記回転体の振動を複数回測定する振動測定手段と、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出手段と、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別手段と、を備えることを特徴とする。
請求項９に記載の発明は、請求項８の構成において、前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする。
請求項１０に記載の発明は、請求項８の構成において、前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする。
請求項１１に記載の発明は、請求項８乃至１０の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定すると共に、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする。
請求項１２に記載の発明は、請求項８乃至１１の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、前記比基本周波数と前記比振動源回転周波数とを合わせて表示する表示手段を備えることを特徴とする。
請求項１３に記載の発明は、請求項８乃至１１の何れかの構成において、前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、前記特定の比周波数の位置と前記回転同期成分強調波形とを合わせて表示する表示手段を備えることを特徴とする。
上記目的を達成するために、請求項１４に記載の発明は、転がり軸受の異常診断プログラムであって、
所定の回転周波数で測定された回転体の振動が前記回転周波数と共に入力されたコンピュータに、請求項１乃至７の何れかに記載の転がり軸受の異常診断方法における回転同期成分強調波形算出ステップと軸受損傷種類判別ステップとを実行させることを特徴とする。
本発明によれば、測定した振動に対し、回転体の回転周波数に比例した周波数の振動成分を強調する回転同期成分強調波形を決定し、回転同期成分強調波形の振動ピーク位置の規則性に基づいて転がり軸受における損傷の種類を判別するので、軸受諸元を用いることなく振動ピークの規則性から軸受異常の種類を判断することが可能となる。
特に、請求項４及び１１の発明によれば、上記効果に加えて、軸受諸元の値を必要とすることなく特徴周波数の値が算出できるため、従来特徴周波数の値を算出するために軸受諸元の値が必要であった診断手法が実行可能となる。よって、軸受異常の種類だけでなく異常の程度も精度よく求めることができる。
特に、請求項５の発明によれば、上記効果に加えて、軸受損傷種類判別ステップでは、回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、回転同期成分強調波形を入力して転がり軸受の損傷の種類を判別するので、人間が考える処理よりも汎用的に振動ピークの規則性から軸受異常の種類を判断することが可能となる。
特に、請求項６の発明によれば、上記効果に加えて、回転同期成分強調波形算出ステップでは、異なる複数の回転周波数で測定された振動から回転同期成分強調波形を算出するので、さらに外乱の影響が低減されて、振動ピークの規則性を捉え易くなるため、軸受異常の推定精度が向上する。
特に、請求項７の発明によれば、上記効果に加えて、振動ピークの比周波数がばらついても確実に軸受損傷に起因する振動ピークを強調した回転同期成分強調波形を算出することが可能となる。
特に、請求項１２の発明によれば、上記効果に加えて、比基本周波数と比振動源回転周波数とを合わせて表示する表示手段を備えることで、数字として表示される周波数の値を用いて他の軸受診断処理を行うことが可能となる。
特に、請求項１３の発明によれば、上記効果に加えて、特定の比周波数の位置と回転同期成分強調波形とを合わせて表示する表示手段を備えることで、振動ピークの規則性を正しく抽出できているか否かを視覚的に把握できるため、診断結果の妥当性を容易に検証できる。また、各損傷の種類において着目すべき比周波数がわかりやすくなる。
特に、請求項４及び１１の発明によれば、上記効果に加えて、軸受諸元の値を必要とすることなく特徴周波数の値が算出できるため、従来特徴周波数の値を算出するために軸受諸元の値が必要であった診断手法が実行可能となる。よって、軸受異常の種類だけでなく異常の程度も精度よく求めることができる。
特に、請求項５の発明によれば、上記効果に加えて、軸受損傷種類判別ステップでは、回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、回転同期成分強調波形を入力して転がり軸受の損傷の種類を判別するので、人間が考える処理よりも汎用的に振動ピークの規則性から軸受異常の種類を判断することが可能となる。
特に、請求項６の発明によれば、上記効果に加えて、回転同期成分強調波形算出ステップでは、異なる複数の回転周波数で測定された振動から回転同期成分強調波形を算出するので、さらに外乱の影響が低減されて、振動ピークの規則性を捉え易くなるため、軸受異常の推定精度が向上する。
特に、請求項７の発明によれば、上記効果に加えて、振動ピークの比周波数がばらついても確実に軸受損傷に起因する振動ピークを強調した回転同期成分強調波形を算出することが可能となる。
特に、請求項１２の発明によれば、上記効果に加えて、比基本周波数と比振動源回転周波数とを合わせて表示する表示手段を備えることで、数字として表示される周波数の値を用いて他の軸受診断処理を行うことが可能となる。
特に、請求項１３の発明によれば、上記効果に加えて、特定の比周波数の位置と回転同期成分強調波形とを合わせて表示する表示手段を備えることで、振動ピークの規則性を正しく抽出できているか否かを視覚的に把握できるため、診断結果の妥当性を容易に検証できる。また、各損傷の種類において着目すべき比周波数がわかりやすくなる。
以下、本発明の実施の形態を図面に基づいて説明する。
図１は転がり軸受の異常診断装置を工作機械の主軸に対して取り付けて使用する場合の構成を示した機能ブロック図で、この図に基づいて具体的に説明する。
工作機械１００において、主軸１が転がり軸受である軸受７を介して主軸ハウジング２に対して回転可能に取り付けられており、主軸１には加工を行うための工具３が固定されている。モータ４は主軸１を駆動する。モータ４には速度検出器５が設けられて、測定されたモータ４の回転周波数が制御装置６に入力されるようになっている。制御装置６は、速度検出器５で測定されたモータ４の回転周波数を測定者によって制御装置６に入力された指令回転周波数に保つようにモータ４へ供給する電流の制御を行っている。表示部８には、速度検出器５で測定されたモータ４の現在の回転周波数が表示される。
図１は転がり軸受の異常診断装置を工作機械の主軸に対して取り付けて使用する場合の構成を示した機能ブロック図で、この図に基づいて具体的に説明する。
工作機械１００において、主軸１が転がり軸受である軸受７を介して主軸ハウジング２に対して回転可能に取り付けられており、主軸１には加工を行うための工具３が固定されている。モータ４は主軸１を駆動する。モータ４には速度検出器５が設けられて、測定されたモータ４の回転周波数が制御装置６に入力されるようになっている。制御装置６は、速度検出器５で測定されたモータ４の回転周波数を測定者によって制御装置６に入力された指令回転周波数に保つようにモータ４へ供給する電流の制御を行っている。表示部８には、速度検出器５で測定されたモータ４の現在の回転周波数が表示される。
コンピュータを含む異常診断装置２００において、振動測定手段としての振動センサ９は、工作機械１００の軸受７の損傷によって生じる振動を測定可能な位置に磁力により取り付けられている。
表示手段としての表示・操作部１３は、測定者が制御装置６へ入力すべき指令回転周波数や診断結果を表示するとともに、表示部８に表示されたモータ４の現在の回転周波数が表示・操作部１３に表示された値と一致したか否かを測定者が判断した結果を入力する操作が可能となっている。振動センサ９で測定される振動加速度は、Ａ／Ｄ変換部１０でデジタル値に変換される。記憶部１１は、表示部８に表示されたモータ４の現在の回転周波数が表示・操作部１３に表示された値と一致したと測定者が判断したときから既定のデータ長となるときまでの振動加速度のデジタル値を振動測定時の回転周波数とともに記憶する。
表示手段としての表示・操作部１３は、測定者が制御装置６へ入力すべき指令回転周波数や診断結果を表示するとともに、表示部８に表示されたモータ４の現在の回転周波数が表示・操作部１３に表示された値と一致したか否かを測定者が判断した結果を入力する操作が可能となっている。振動センサ９で測定される振動加速度は、Ａ／Ｄ変換部１０でデジタル値に変換される。記憶部１１は、表示部８に表示されたモータ４の現在の回転周波数が表示・操作部１３に表示された値と一致したと測定者が判断したときから既定のデータ長となるときまでの振動加速度のデジタル値を振動測定時の回転周波数とともに記憶する。
演算部１２は、記憶部１１に記憶された振動測定時の回転周波数と振動加速度より、予め記憶された異常診断プログラムに基づき、乗算、フーリエ変換、絶対値の算出、内挿処理を行って、比周波数（周波数を回転周波数で割って算出される無次元量）ごとの振動の大きさを算出し、複数の振動測定時の回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する。さらに、回転同期成分強調波形のピーク位置の規則性より、後述する比基本周波数と比振動源回転周波数とを推定し、比振動源回転周波数の値に応じて、損傷が内輪・転動体・外輪の何れに生じているかを判断する。すなわち、演算部１２は、回転同期成分強調波形算出手段と軸受損傷種類判別手段として機能する。
軸受７において、内輪傷が存在する場合、転動体が内輪傷を通過する周波数（基本周波数）とその自然数倍の周波数の力が発生する。内輪傷の方向は主軸１の回転とともに回転するため、力の発生する方向は回転周波数で変化する。
軸受７において、転動体傷が存在する場合、転動体傷が内輪または外輪に接する周波数（基本周波数）とその自然数倍の周波数の力が発生する。傷のある転動体の位置は転動体の保持器が自転する（転動体が公転する）周波数で変化するため、力の発生する方向は転動体の保持器が自転する（転動体が公転する）周波数で変化する。保持器の自転（転動体の公転）は、主軸の回転よりも必ず遅くなるので、力の発生する方向が変化する周波数は０Ｈｚより大きく回転周波数より小さい値となる。
軸受７において、外輪傷が存在する場合、転動体が外輪傷を通過する周波数（基本周波数）とその自然数倍の周波数の力が発生する。外輪傷の方向は主軸１の回転に依らず一定であるため、力の発生する方向は変化しない。
ここでは力の発生する方向が変化する周波数を振動源回転周波数と呼称する。力の発生する方向の変化しない場合は、振動源回転周波数が便宜上０Ｈｚであるとみなしてよい。
軸受７において、転動体傷が存在する場合、転動体傷が内輪または外輪に接する周波数（基本周波数）とその自然数倍の周波数の力が発生する。傷のある転動体の位置は転動体の保持器が自転する（転動体が公転する）周波数で変化するため、力の発生する方向は転動体の保持器が自転する（転動体が公転する）周波数で変化する。保持器の自転（転動体の公転）は、主軸の回転よりも必ず遅くなるので、力の発生する方向が変化する周波数は０Ｈｚより大きく回転周波数より小さい値となる。
軸受７において、外輪傷が存在する場合、転動体が外輪傷を通過する周波数（基本周波数）とその自然数倍の周波数の力が発生する。外輪傷の方向は主軸１の回転に依らず一定であるため、力の発生する方向は変化しない。
ここでは力の発生する方向が変化する周波数を振動源回転周波数と呼称する。力の発生する方向の変化しない場合は、振動源回転周波数が便宜上０Ｈｚであるとみなしてよい。
ところで、力の発生する方向が振動源回転周波数で変化する場合、固定された座標系でその力を見ると、力の周波数（基本周波数の自然数倍）に対して、振動源回転周波数だけ前後した２つの周波数の力が生じていることになる。線形とみなせるシステムでは、力と同じ周波数の振動加速度が生じるため、固定された座標系で振動加速度を測定し、周波数分析すると、軸受が損傷している場合には、理想的には振動源回転周波数の２倍だけ離れた周波数に振動ピークの対が観測され、その２つの振動ピークの平均の周波数が、基本周波数の自然数倍となるような位置に複数観測される。なお、外輪傷の場合は振動源回転周波数を０Ｈｚとみなすことで、同一の周波数の振動ピークを２回カウントして２つの振動ピークとみなすこととする。この規則性を捉え、振動源回転周波数を推定することで異常の種類を判断することが可能となる。ただし、現実的には、振動を測定し周波数分析するだけでは、振動しにくい周波数の振動ピークはノイズに埋もれてしまい、この規則性を確実に捉えることは難しい。
回転に同期する軸受損傷に起因する振動は、繰り返し測定したり異なる回転周波数で測定したりしたとしても、振動の生じる周波数は常に回転周波数の定数倍である。振動の周波数の回転周波数に対する比の無次元量を比周波数と呼称することにすると、回転に同期する軸受損傷に起因する振動は、回転周波数に依らず同じ比周波数に振動ピークが観測されることになる。このため、複数回測定して比周波数ごとの振動の大きさを算出し、複数回の測定で得られた振動の大きさについて係数が全て正の線型和を算出してホワイトノイズを低減したり、複数の回転周波数で測定して比周波数ごとの振動の大きさを算出し、複数回の測定で得られた振動の大きさについて係数が全て正の線型和を算出してホワイトノイズや回転に同期しない外乱成分を低減したりすることで、回転に同期した振動を強調した波形（回転同期成分強調波形）を算出することができ、軸受損傷に起因する振動ピークの規則性を捉えることが容易となる。
なお、係数に異符号の値がある場合には、軸受損傷に起因する振動ピークの振動の大きさを打ち消しあうことになるため、回転に同期した軸受損傷に起因する振動を強調する効果が得られない。全ての係数が１の線型和をとった場合、回転同期成分強調波形は比周波数毎の振動の大きさの合計値となる。全ての係数を測定回数の逆数として線型和をとった場合、回転同期成分強調波形は比周波数毎の振動の大きさの平均値となる。合計値であっても平均値であっても回転同期成分強調波形のピーク位置は同一となるためピーク位置の規則性を議論する上ではいずれであってもよい。また、係数を同一に限らなくとも概ねピーク位置の規則性について同様の議論が成立する。例えば、回転周波数が大きいほど軸受損傷により生じる力の大きさが大きくなる傾向があるため、測定する回転周波数範囲が広い場合には低い回転周波数で測定したデータの影響度合いが小さくなってしまう。このため、回転周波数のべき乗の逆数を係数として線型和をとることで、異なる回転周波数で得られた振動の大きさを平等に扱うことも可能である。
複数の回転周波数で測定する場合、振動の基本周波数と振動源回転周波数とをそれぞれ回転周波数で除算して得られる無次元量の比基本周波数と比振動源回転周波数とを用いた方が議論しやすい。但し、比基本周波数は、軸受の構造によって様々な値をとるため、例えばある軸受の内輪傷の比基本周波数が別の軸受の外輪傷の比基本周波数と同じ大きさとなることもある。このため、比基本周波数の値に着目しても軸受の異常の種類を判別することは不可能である。しかし、前述したように、振動源回転周波数は、軸受損傷の種類によって回転周波数の何倍となるかが異なり、その大小関係に逆転が起こることはない。よって、比振動源回転周波数に換算すると、内輪損傷の場合は１、転動体損傷の場合は、後述する回転周波数に対する転動体公転の比の取りうる範囲の値（工作機械主軸に用いられる軸受の場合には０．４５前後の値が多い）、外輪損傷の場合は０となる。
図２に外輪および内輪が損傷している場合の回転同期成分強調波形、図３に転動体が損傷している場合の回転同期成分強調波形を示す。これらの図に基づいて、軸受損傷の種類毎の振動ピーク位置の規則性を説明する。なお、図２において、内輪損傷時の振動ピークに関連する箇所は破線で、外輪損傷時の振動ピークに関連する箇所は実線で、注釈を記した。図３において、転動体損傷時の振動ピークに関連する箇所は実線で注釈を記した。
内輪損傷があると、図２に示すように、同一の比振動源回転周波数（＝１）の２倍だけ離れた２本の振動ピーク（ピーク対）が複数箇所に存在する。さらに、それぞれのピーク対の平均の比周波数は必ず比基本周波数の自然数倍となるため、比基本周波数が未知であっても、ピーク対の平均の比周波数が自然数比となっていることを確認することで内輪損傷による振動ピークであると確信することができる。
外輪損傷があると、図２に示すように、比周波数の比が自然数比となるような振動ピークが複数箇所に存在する。軸受の比基本周波数は自然数とならないように設計されることが一般的であるため、さらにそれらの振動ピークの比周波数のいくつかが自然数でなければ、比基本周波数が自然数でないと推測できるため、外輪損傷による振動ピークであると確信することができる。なお、外輪損傷の場合、対でない１本の振動ピークが等間隔に存在するだけであるが、便宜上、比振動源回転周波数０だけ離れた２本の振動ピーク（ピーク対）が存在するとみなすことで他の種類の異常と同じ計算処理を用いることができる。
転動体損傷があると、図３に示すように、同一の比振動源回転周波数の２倍だけ離れた２本の振動ピーク（ピーク対）が複数箇所に存在する。さらに、それぞれのピーク対の平均の比周波数は必ず比基本周波数の自然数倍となるため、比基本周波数が未知であっても、ピーク対の平均の比周波数が自然数比となっていることを確認することで転動体損傷による振動ピークであると確信することができる。なお、比振動源回転周波数の値は必ず１より小さい値となるため、内輪損傷との識別が可能である。
また、回転同期成分強調波形の振動ピーク位置の規則性（回転同期成分強調波形の比基本周波数×自然数±比振動源回転周波数の位置に振動ピークが存在する）を捉える方法として、ここではピーク対の平均（外輪の場合はピーク）の比周波数が自然数比となっていることを確認する手法を例に挙げたが他の手法であっても構わない。
外輪損傷があると、図２に示すように、比周波数の比が自然数比となるような振動ピークが複数箇所に存在する。軸受の比基本周波数は自然数とならないように設計されることが一般的であるため、さらにそれらの振動ピークの比周波数のいくつかが自然数でなければ、比基本周波数が自然数でないと推測できるため、外輪損傷による振動ピークであると確信することができる。なお、外輪損傷の場合、対でない１本の振動ピークが等間隔に存在するだけであるが、便宜上、比振動源回転周波数０だけ離れた２本の振動ピーク（ピーク対）が存在するとみなすことで他の種類の異常と同じ計算処理を用いることができる。
転動体損傷があると、図３に示すように、同一の比振動源回転周波数の２倍だけ離れた２本の振動ピーク（ピーク対）が複数箇所に存在する。さらに、それぞれのピーク対の平均の比周波数は必ず比基本周波数の自然数倍となるため、比基本周波数が未知であっても、ピーク対の平均の比周波数が自然数比となっていることを確認することで転動体損傷による振動ピークであると確信することができる。なお、比振動源回転周波数の値は必ず１より小さい値となるため、内輪損傷との識別が可能である。
また、回転同期成分強調波形の振動ピーク位置の規則性（回転同期成分強調波形の比基本周波数×自然数±比振動源回転周波数の位置に振動ピークが存在する）を捉える方法として、ここではピーク対の平均（外輪の場合はピーク）の比周波数が自然数比となっていることを確認する手法を例に挙げたが他の手法であっても構わない。
図４は、転がり軸受の異常診断を行う方法のフローチャートを示したものであり、このフローチャートに基づいて具体的に説明する。
まず、予め登録されている診断回転周波数の全条件のうちの１つが表示・操作部１３に図７のように表示される（Ｓ１）。
次に、診断装置の使用者は、表示された診断回転周波数を制御装置６へ指令する（Ｓ２）。
次に、診断装置の使用者は、表示部８を確認して、診断対象の軸受に支持された回転体が表示・操作部１３に表示された回転周波数で回転しているかを判断し、一致している場合には表示・操作部１３より測定開始を指令する（Ｓ３）。
すると、異常診断装置２００では、振動加速度を測定して記録し（Ｓ４）、必要なデータ長となったら測定完了と判断してＳ６へ移行する（Ｓ５）。Ｓ６では、振動測定時の回転周波数と振動加速度とを対応付けて記録する。
Ｓ７で、予め登録されている診断回転周波数の全条件の測定が完了しているかを判断し、全条件の測定が完了していればＳ８へ移行する。完了してない場合にはＳ１に戻り、測定が完了していない別の診断回転周波数を表示する。ここでのＳ４～Ｓ７までが振動測定ステップとなる。
まず、予め登録されている診断回転周波数の全条件のうちの１つが表示・操作部１３に図７のように表示される（Ｓ１）。
次に、診断装置の使用者は、表示された診断回転周波数を制御装置６へ指令する（Ｓ２）。
次に、診断装置の使用者は、表示部８を確認して、診断対象の軸受に支持された回転体が表示・操作部１３に表示された回転周波数で回転しているかを判断し、一致している場合には表示・操作部１３より測定開始を指令する（Ｓ３）。
すると、異常診断装置２００では、振動加速度を測定して記録し（Ｓ４）、必要なデータ長となったら測定完了と判断してＳ６へ移行する（Ｓ５）。Ｓ６では、振動測定時の回転周波数と振動加速度とを対応付けて記録する。
Ｓ７で、予め登録されている診断回転周波数の全条件の測定が完了しているかを判断し、全条件の測定が完了していればＳ８へ移行する。完了してない場合にはＳ１に戻り、測定が完了していない別の診断回転周波数を表示する。ここでのＳ４～Ｓ７までが振動測定ステップとなる。
Ｓ８では、まず振動加速度をフーリエ変換し振幅スペクトル密度を算出する。予め決めておいた全条件の測定に対して共通の比周波数分解能となるように振幅スペクトル密度の絶対値について補間処理を行い比周波数毎の振動の大きさを算出する。それぞれの診断回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する（回転同期成分強調波形算出ステップ）。全Ｎ条件のうちｉ番目の診断回転周波数をｆＲＯＴ(i)、線型和を取る際のｉ番目の条件の振動の大きさＶ(i)に対する係数をＡ(i)、回転同期成分強調波形をＷとする。このとき回転同期成分強調波形は以下の数１で得られる。振動の大きさの合計を使う場合は、係数Ａ(i)は以下の数２を用いればよい。振動の大きさの平均を使う場合は、係数Ａ(i)は以下の数３を用いればよい。回転周波数により軸受損傷により生じる力がｆＲＯＴ(i)のＭ乗に比例すると近似できる場合に全ての回転周波数のデータを平等に扱うためには、係数Ａ(i)を以下の数４とすればよい。回転周波数がｆＲＯＴｂｅｓｔのとき最も測定精度が良いことがわかっている場合には、回転周波数がｆＲＯＴｂｅｓｔに近いときの振動の大きさの影響度合いが大きくなるように以下の数５のような係数Ａ(i)を用いてもよい。すなわち、数１～５の係数Ａ(i)は任意の一部又は全部を選択できる。ここで、振動の大きさＶ(i)と回転同期成分強調波形Ｗは、比周波数毎に値を持つベクトル、Ａ(i)はスカラーである。
次に、Ｓ９では、比基本周波数と比振動源回転周波数とを仮定し、回転同期成分強調波形の比基本周波数×自然数±比振動源回転周波数の位置に振動ピークが存在するかそれぞれ判断する。ある比周波数に振動ピークが存在するか否かの判断は、回転同期成分強調波形と回転同期成分強調波形の移動平均を比較し、回転同期成分強調波形の値の方が大きければ、その比周波数には振動ピークが存在すると判断すればよい。そして、比基本周波数×自然数±比振動源回転周波数の位置に振動ピークが存在した割合（一致率）を算出する。
次に、Ｓ１０では、既定された探索範囲の比基本周波数と比振動源回転周波数との組み合わせについて全て探索完了しているか判断し、完了していればＳ１１へ移行する。
Ｓ１１では、一致率が既定のしきい値以上である比基本周波数と比振動源回転周波数について、比振動源回転周波数の値に応じて軸受異常の種類（内輪損傷、転動体損傷、外輪損傷）を判断する。
さらに、Ｓ１２では、Ｓ１１で軸受異常の種類が特定された比基本周波数と比振動源回転周波数の組と診断回転周波数とから、基本周波数と特徴周波数の値を算出し、損傷の程度（修理が必要なレベルか否か）を診断する。例えば先の特許文献２に開示されているように、算出した基本周波数及び特徴周波数を、各周波数ごとに予め設定したしきい値と比較すれば、損傷の程度を診断することができる。ここでのＳ９～Ｓ１２までが軸受損傷種類判別ステップとなる。
そして、Ｓ１３では、図８に示すような診断結果を表示・操作部１３に表示する。ここでは比振動源回転周波数が外輪に対応する０、内輪に対応する１の場合に一致率が高く算出され損傷していると判断されたが、Ｓ１２の診断では異常（修理が必要）と判定するしきい値を超過しなかったため、損傷の程度では「正常範囲」と表示され、「外輪と内輪が損傷していますが、修理が必要な水準ではありません」との文字も表示される。
Ｓ１１では、一致率が既定のしきい値以上である比基本周波数と比振動源回転周波数について、比振動源回転周波数の値に応じて軸受異常の種類（内輪損傷、転動体損傷、外輪損傷）を判断する。
さらに、Ｓ１２では、Ｓ１１で軸受異常の種類が特定された比基本周波数と比振動源回転周波数の組と診断回転周波数とから、基本周波数と特徴周波数の値を算出し、損傷の程度（修理が必要なレベルか否か）を診断する。例えば先の特許文献２に開示されているように、算出した基本周波数及び特徴周波数を、各周波数ごとに予め設定したしきい値と比較すれば、損傷の程度を診断することができる。ここでのＳ９～Ｓ１２までが軸受損傷種類判別ステップとなる。
そして、Ｓ１３では、図８に示すような診断結果を表示・操作部１３に表示する。ここでは比振動源回転周波数が外輪に対応する０、内輪に対応する１の場合に一致率が高く算出され損傷していると判断されたが、Ｓ１２の診断では異常（修理が必要）と判定するしきい値を超過しなかったため、損傷の程度では「正常範囲」と表示され、「外輪と内輪が損傷していますが、修理が必要な水準ではありません」との文字も表示される。
図５は、転がり軸受の異常診断を行う他の方法のフローチャートを示したものであり、このフローチャートに基づいて具体的に説明する。
Ｓ１～Ｓ８までは図４の方法と同じである。まず、予め登録されている診断回転周波数の全条件のうちの１つが表示・操作部１３に図７のように表示される（Ｓ１）。診断装置の使用者は、表示された診断回転周波数を制御装置６へ指令する（Ｓ２）。診断装置の使用者は、表示部８を確認して、診断対象の軸受に支持された回転体が表示・操作部１３に表示された回転周波数で回転しているかを判断し、一致している場合には表示・操作部１３より測定開始を指令する（Ｓ３）。振動加速度を測定して記録し（Ｓ４）、必要なデータ長となったら測定完了と判断してＳ６へ移行する（Ｓ５）。振動測定時の回転周波数と振動加速度を対応付けて記録する（Ｓ６）。予め登録されている診断回転周波数の全条件の測定が完了しているかを判断し、全条件の測定が完了していればＳ８へ移行する。完了してない場合にはＳ１に戻り、測定が完了していない別の診断回転周波数を表示する。
Ｓ８では、まず振動加速度をフーリエ変換し振幅スペクトル密度を算出する。予め決めておいた全条件の測定に対して共通の比周波数分解能となるように振幅スペクトル密度の絶対値について補間処理を行い比周波数毎の振動の大きさを算出する。それぞれの診断回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する（回転同期成分強調波形算出ステップ）。
Ｓ１～Ｓ８までは図４の方法と同じである。まず、予め登録されている診断回転周波数の全条件のうちの１つが表示・操作部１３に図７のように表示される（Ｓ１）。診断装置の使用者は、表示された診断回転周波数を制御装置６へ指令する（Ｓ２）。診断装置の使用者は、表示部８を確認して、診断対象の軸受に支持された回転体が表示・操作部１３に表示された回転周波数で回転しているかを判断し、一致している場合には表示・操作部１３より測定開始を指令する（Ｓ３）。振動加速度を測定して記録し（Ｓ４）、必要なデータ長となったら測定完了と判断してＳ６へ移行する（Ｓ５）。振動測定時の回転周波数と振動加速度を対応付けて記録する（Ｓ６）。予め登録されている診断回転周波数の全条件の測定が完了しているかを判断し、全条件の測定が完了していればＳ８へ移行する。完了してない場合にはＳ１に戻り、測定が完了していない別の診断回転周波数を表示する。
Ｓ８では、まず振動加速度をフーリエ変換し振幅スペクトル密度を算出する。予め決めておいた全条件の測定に対して共通の比周波数分解能となるように振幅スペクトル密度の絶対値について補間処理を行い比周波数毎の振動の大きさを算出する。それぞれの診断回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する（回転同期成分強調波形算出ステップ）。
そして、Ｓ１０９では、まず振動ピークが存在すると判断された比周波数を４つ選択する。ある比周波数に振動ピークが存在するか否かの判断は、回転同期成分強調波形と回転同期成分強調波形の移動平均とを比較し、回転同期成分強調波形の値の方が大きければ、その比周波数には振動ピークが存在すると判断すればよい。
次に、Ｓ１１０では、選択された４つの振動ピークの比周波数Ｆ１、Ｆ２、Ｆ３、Ｆ４と、公約数を持たない異なる自然数Ｎ１、Ｎ２との組み合わせが以下の数６および数７の関係式を満たすか否かを判断する。なお、Ｆ１、Ｆ２、Ｆ３、Ｆ４には同一の比周波数が選ばれても良い（Ｆ１とＦ２、Ｆ３とＦ４がそれぞれ同一の比周波数のときに外輪損傷の振動ピークを検出できる）。
次に、Ｓ１１０では、選択された４つの振動ピークの比周波数Ｆ１、Ｆ２、Ｆ３、Ｆ４と、公約数を持たない異なる自然数Ｎ１、Ｎ２との組み合わせが以下の数６および数７の関係式を満たすか否かを判断する。なお、Ｆ１、Ｆ２、Ｆ３、Ｆ４には同一の比周波数が選ばれても良い（Ｆ１とＦ２、Ｆ３とＦ４がそれぞれ同一の比周波数のときに外輪損傷の振動ピークを検出できる）。
Ｓ１１０の判別で、４つの振動ピークの比周波数Ｆ１、Ｆ２、Ｆ３、Ｆ４と、公約数を持たない異なる自然数Ｎ１、Ｎ２の組み合わせが数６および数７の関係式を満たす場合、Ｓ１１１へ移行し、関係式を満たさない場合、Ｓ１０９に戻って異なる比周波数を４つ選択してＳ１１０の判別を再度行う。
Ｓ１１１では、まず数６および数７の関係式を満たす振動ピークの比周波数の組み合わせについて、以下の数８、数９より比基本周波数Ｒ０、比振動源回転周波数Ｒ１を算出する。
Ｓ１１１では、まず数６および数７の関係式を満たす振動ピークの比周波数の組み合わせについて、以下の数８、数９より比基本周波数Ｒ０、比振動源回転周波数Ｒ１を算出する。
そして、算出した比振動源回転周波数Ｒ１の値に応じて軸受異常の種類（内輪損傷、転動体損傷、外輪損傷）を判断し、Ｓ１１２で、図９に示すような診断結果を表示・操作部１３に表示する。ここでは比振動源回転周波数が転動体に対応する０．４５と算出されたため、転動体が損傷している可能性ありと表示される。ここでのＳ１０９～Ｓ１１１が軸受損傷種類判別ステップとなる。
図６は、機械学習を用いて転がり軸受の異常診断を行う他の方法のフローチャートを示したものであり、このフローチャートに基づいて具体的に説明する。
ここでは軸受異常が既知の軸受について後述する方法で算出される回転同期成分強調波形を入力とし、内輪損傷がある（１）か否か（０）、転動体損傷がある（１）か否か（０）、外輪損傷がある（１）か否か（０）、を出力とする教師データを用いて外部の学習器（図１に図示しない）において、予め機械学習させておき、学習済み数理モデルを図１の演算部１２に備えている。数理モデルの形態としては、例えば、多層ニューラルネットワークなどを用いることができる。軸受異常が既知の軸受について評価用データ（回転同期成分強調波形）の値を入力した場合の出力の分布より、内輪、転動体、外輪が異常か否かを判断するしきい値を予め決定しておき、演算部１２に備えている。
ここでは軸受異常が既知の軸受について後述する方法で算出される回転同期成分強調波形を入力とし、内輪損傷がある（１）か否か（０）、転動体損傷がある（１）か否か（０）、外輪損傷がある（１）か否か（０）、を出力とする教師データを用いて外部の学習器（図１に図示しない）において、予め機械学習させておき、学習済み数理モデルを図１の演算部１２に備えている。数理モデルの形態としては、例えば、多層ニューラルネットワークなどを用いることができる。軸受異常が既知の軸受について評価用データ（回転同期成分強調波形）の値を入力した場合の出力の分布より、内輪、転動体、外輪が異常か否かを判断するしきい値を予め決定しておき、演算部１２に備えている。
Ｓ１～Ｓ８までは図４の方法と同じである。予め登録されている診断回転周波数の全条件のうちの１つが表示・操作部１３に図７のように表示される（Ｓ１）。診断装置の使用者は、表示された診断回転周波数を制御装置６へ指令する（Ｓ２）。診断装置の使用者は、表示部８を確認して、診断対象の軸受に支持された回転体が表示・操作部１３に表示された回転周波数で回転しているかを判断し、一致している場合には表示・操作部１３より測定開始を指令する（Ｓ３）。振動加速度を測定して記録し（Ｓ４）、必要なデータ長となったら測定完了と判断してＳ６へ移行する（Ｓ５）。振動測定時の回転周波数と振動加速度を対応付けて記録する（Ｓ６）。予め登録されている診断回転周波数の全条件の測定が完了しているかを判断し、全条件の測定が完了していればＳ８へ移行する。完了してない場合にはＳ１に戻り、測定が完了していない別の診断回転周波数を表示する。
Ｓ８では、まず振動加速度をフーリエ変換し振幅スペクトル密度を算出する。予め決めておいた全条件の測定に対して共通の比周波数分解能となるように振幅スペクトル密度の絶対値について補間処理を行い比周波数毎の振動の大きさを算出する。それぞれの診断回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する（回転同期成分強調波形算出ステップ）。
Ｓ８では、まず振動加速度をフーリエ変換し振幅スペクトル密度を算出する。予め決めておいた全条件の測定に対して共通の比周波数分解能となるように振幅スペクトル密度の絶対値について補間処理を行い比周波数毎の振動の大きさを算出する。それぞれの診断回転周波数の振動の大きさについて係数が全て正の線型和を算出することで回転同期成分強調波形を算出する（回転同期成分強調波形算出ステップ）。
そして、Ｓ２０９では、回転同期成分強調波形を学習済み数理モデルに入力し、内輪損傷があるか否か、転動体損傷があるか否か、外輪損傷があるか否か、の出力値をそれぞれ算出する（軸受損傷種類判別ステップ）。
Ｓ２０９で、それぞれの出力値の値が予め設定されたしきい値を超過していれば異常と判断し、Ｓ２１０で、図１０に示すような診断結果を表示・操作部１３に表示する。図１０には、学習済み数理モデルに入力する回転同期成分強調波形がグラフとして、学習済み数理モデルからの出力値、異常か否かを判断するしきい値、異常か否かを判断した判定結果が表示されている。
Ｓ２０９で、それぞれの出力値の値が予め設定されたしきい値を超過していれば異常と判断し、Ｓ２１０で、図１０に示すような診断結果を表示・操作部１３に表示する。図１０には、学習済み数理モデルに入力する回転同期成分強調波形がグラフとして、学習済み数理モデルからの出力値、異常か否かを判断するしきい値、異常か否かを判断した判定結果が表示されている。
このように、上記各形態の異常診断方法及び異常診断装置２００、異常診断プログラムによれば、測定した振動に対し、主軸１の回転周波数に比例した周波数の振動成分を強調する回転同期成分強調波形を算出し、回転同期成分強調波形の振動ピーク位置の規則性に基づいて軸受７における損傷の種類を判別するので、軸受諸元を用いることなく振動ピークの規則性から軸受異常の種類を判断することが可能となる。
さらに、比基本周波数と比振動源回転周波数を推定する処理の場合、基本周波数や特徴周波数の算出が行える。このため、基本周波数や特徴周波数を算出するために軸受諸元を利用する必要があった従来の診断技術（例えば特許文献２）を採用して損傷レベルの診断をさらに行うことも可能となる。ある回転周波数における基本周波数が必要な場合は、比基本周波数の値に回転周波数を乗算することで算出することが可能である。ある回転周波数における、転動体の保持器が自転する（転動体が公転する）周波数が必要な場合には、転動体損傷と判断された比基本周波数と比振動源回転周波数の組み合わせの比振動源回転周波数の値に回転周波数を乗算することで算出することが可能である。この損傷レベルの診断処理は、図５の異常診断方法においても採用できる。
さらに、比基本周波数と比振動源回転周波数を推定する処理の場合、基本周波数や特徴周波数の算出が行える。このため、基本周波数や特徴周波数を算出するために軸受諸元を利用する必要があった従来の診断技術（例えば特許文献２）を採用して損傷レベルの診断をさらに行うことも可能となる。ある回転周波数における基本周波数が必要な場合は、比基本周波数の値に回転周波数を乗算することで算出することが可能である。ある回転周波数における、転動体の保持器が自転する（転動体が公転する）周波数が必要な場合には、転動体損傷と判断された比基本周波数と比振動源回転周波数の組み合わせの比振動源回転周波数の値に回転周波数を乗算することで算出することが可能である。この損傷レベルの診断処理は、図５の異常診断方法においても採用できる。
特にここでは、振動測定ステップでは、異なる複数の回転周波数で振動を測定し、回転同期成分強調波形算出ステップでは、振動を周波数分析して共通の比周波数に対応する振動の大きさをそれぞれ算出し、係数が全て正の線型和をとるようにしている（Ｓ８）ので、さらに外乱の影響が低減されて、振動ピークの規則性を捉え易くなるため、軸受異常の推定精度が向上する。
また、図６に示す異常診断方法では、回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、回転同期成分強調波形を入力して軸受損傷の種類を判別するので、人間が考える処理よりも汎用的に振動ピークの規則性から軸受異常の種類を判断することが可能となる。
そして、異常診断装置２００では、軸受損傷に対応する比基本周波数と比振動源回転周波数とを合わせて表示・操作部１３に表示するので、数字として表示される周波数の値を用いて異常診断装置２００とは別の異常診断装置による軸受診断処理を行うことが可能となる。
さらにここでは、特定の比周波数の位置と回転同期成分強調波形とを合わせて表示・操作部１３に表示するので、振動ピークの規則性を正しく抽出できているか否かを視覚的に把握できるため、診断結果の妥当性を容易に検証できる。また、各損傷の種類において着目すべき比周波数がわかりやすくなる。
そして、異常診断装置２００では、軸受損傷に対応する比基本周波数と比振動源回転周波数とを合わせて表示・操作部１３に表示するので、数字として表示される周波数の値を用いて異常診断装置２００とは別の異常診断装置による軸受診断処理を行うことが可能となる。
さらにここでは、特定の比周波数の位置と回転同期成分強調波形とを合わせて表示・操作部１３に表示するので、振動ピークの規則性を正しく抽出できているか否かを視覚的に把握できるため、診断結果の妥当性を容易に検証できる。また、各損傷の種類において着目すべき比周波数がわかりやすくなる。
なお、上記各例では、診断対象の軸受に支持された回転体の回転周波数を直接指令可能な例を示したが、ベルトや歯車を介して工作機械主軸を回転させるモータの支持軸受を診断する場合は、減速比を考慮して、モータの支持軸受の回転周波数が診断周波数となるように指定してやればよい。
また、表示・操作部１３に表示された診断回転周波数に診断対象の軸受の回転周波数を合わせることで、振動加速度と回転周波数を対応付けて記憶する方法を示したが、診断対象の軸受の回転周波数を表示・操作部１３より診断装置の使用者が入力する方法や、異常診断装置２００が速度検出器５の情報を直接取り込めるようにする方法でもよい。
さらに、回転周波数（単位はＨｚを用いることが多い）と回転速度（工作機械主軸では、単位ｍｉｎ－１が用いられることが多い）は、１Ｈｚ＝６０ｍｉｎ－１の関係がある同一の量であり、どちらを用いてもよい。
また、表示・操作部１３に表示された診断回転周波数に診断対象の軸受の回転周波数を合わせることで、振動加速度と回転周波数を対応付けて記憶する方法を示したが、診断対象の軸受の回転周波数を表示・操作部１３より診断装置の使用者が入力する方法や、異常診断装置２００が速度検出器５の情報を直接取り込めるようにする方法でもよい。
さらに、回転周波数（単位はＨｚを用いることが多い）と回転速度（工作機械主軸では、単位ｍｉｎ－１が用いられることが多い）は、１Ｈｚ＝６０ｍｉｎ－１の関係がある同一の量であり、どちらを用いてもよい。
一方、すべりのない場合の、転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）は、転動体直径ｄ、ピッチ円直径Ｄ、接触角αを用いて、以下の数１０で表現される。
数１０より、転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）はピッチ円直径に対する転動体直径の比が大きいほど小さくなることがわかる。転がり軸受を構成する上で最低３個の転動体が必要となるので、３個の転動体を密に並べた構造の軸受においてピッチ円直径に対する転動体直径の比が最大となるため、幾何学的に決定される転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）の下限値は、（１－√３÷２）÷２≒０．０６６９である。数１０より、転動体直径を無限小にした場合が、転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）が最大となることがわかる。このため、幾何学的には、転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）の上限値は０．５である。よって、転がり軸受の転動体の公転周波数を回転周波数で除算した値（転動体異常の場合の比振動源回転周波数Ｒ１）がとり得る範囲は、最も広く見積もっても（１－√３÷２）÷２以上から０．５以下である。この範囲は、内輪異常の場合の比振動源回転周波数である１、外輪異常の場合の比振動源回転周波数である０を含んでいないため、内輪異常・転動体異常・外輪異常のいずれに分類すればよいか判別できないという状況は生じない。
工作機械主軸に用いられる軸受の場合、転動体異常の場合の比振動源回転周波数は０．４５前後であるため、仮定する比振動源回転周波数は０、０．４以上０．５以下、１とするなどさらに限定して探索を行ってもよい。
工作機械主軸に用いられる軸受の場合、転動体異常の場合の比振動源回転周波数は０．４５前後であるため、仮定する比振動源回転周波数は０、０．４以上０．５以下、１とするなどさらに限定して探索を行ってもよい。
そして、各形態では、測定毎に比周波数分解能が一致しない場合に比周波数分解能を統一する方法として補間処理を行う例を示したが、軸受の接触角の変化や回転周波数の誤差により軸受損傷に起因する振動ピークの周波数がＸ％程度前後する場合は、比周波数分解能統一後のある比周波数の振動の大きさとして採用する値はその比周波数の±Ｘ％の範囲の最大値を採用するとよい。このように比周波数の前後に所定の幅を持たせて当該幅内の最大値を採用することにより、軸受の接触角の変化や回転周波数の誤差によって振動ピークの位置がＸ％ずれたとしても、必ずある比周波数において軸受損傷に起因する振動ピークの振動の大きさを合算した値が得られることになる。なお、比周波数分解能が荒いため、ある比周波数の±Ｘ％の範囲にデータが存在しない場合には、±比周波数分解能分の範囲の最大値を採用する。なお、Ｘの値の決定は、例えば実際に複数の回転周波数において、軸受諸元が既知の異常軸受に支持された主軸の振動を測定し、軸受諸元から算出される振動ピークの比周波数に最も近い振動ピークの比周波数を求め、複数の回転周波数において求めた比周波数の標準偏差の数倍（例えば３～６倍）を、複数の回転周波数において求めた比周波数の平均で割って算出するといった方法がある。このように、一旦軸受諸元が既知の異常軸受でばらつきの度合いＸ％を把握しておけば、このＸ％を軸受諸元が未知の軸受に対しても用いることができる。
一方、各形態では、係数が全て正の線型和をとることで回転同期成分強調波形を算出し、回転同期成分強調波形と回転同期成分強調波形の移動平均を比較し、回転同期成分強調波形の値の方が大きければ、その比周波数には振動ピークが存在すると判断するような局所的に最大となる振動ピークを抽出する例を示したが、係数が全て負の線型和をとることで回転同期成分強調波形を算出し、回転同期成分強調波形と回転同期成分強調波形の移動平均を比較し、回転同期成分強調波形の値の方が小さければ、その比周波数には振動ピークが存在すると判断する局所的に最小となる振動ピークを抽出する手法でも全く同様の議論が成立する。
また、軸受損傷種類判別ステップで用いる回転同期成分強調波形の比周波数の範囲は、算出可能な全ての比周波数の範囲とする必要はなく、軸受損傷に起因する振動ピークの規則性を捉え易い範囲にすることが望ましい。
さらに、機械学習を用いて転がり軸受の異常診断を行う方法の教師データの出力として、０と１の２値を与える手法を示したが、損傷度合いを定量的に表現する連続的な値を教師データとすることも可能である。
加えて、比周波数、比基本周波数、比振動源回転周波数、これらの値と比較される値、の全てに同一の値を乗算しても、全く同じ議論が成立する。つまり、例えば比振動源回転周波数が１と一致するか否かを判断することは、比振動源回転周波数×３０Ｈｚが１×３０Ｈｚと一致するか否かを判断することと同義である。
また、軸受損傷種類判別ステップで用いる回転同期成分強調波形の比周波数の範囲は、算出可能な全ての比周波数の範囲とする必要はなく、軸受損傷に起因する振動ピークの規則性を捉え易い範囲にすることが望ましい。
さらに、機械学習を用いて転がり軸受の異常診断を行う方法の教師データの出力として、０と１の２値を与える手法を示したが、損傷度合いを定量的に表現する連続的な値を教師データとすることも可能である。
加えて、比周波数、比基本周波数、比振動源回転周波数、これらの値と比較される値、の全てに同一の値を乗算しても、全く同じ議論が成立する。つまり、例えば比振動源回転周波数が１と一致するか否かを判断することは、比振動源回転周波数×３０Ｈｚが１×３０Ｈｚと一致するか否かを判断することと同義である。
そして、振動測定ステップは、一定の回転周波数のときに実行する場合に限らない。例えば工作機械の主軸のような回転負荷が小さい回転体の場合、惰性回転をさせた場合の単位時間当たりの回転周波数の変化は非常に小さい。このため、回転体の振動を複数回測定する振動測定ステップは、予め設定した回転周波数に主軸を制御した後惰性回転させ、停止若しくはある回転数に至るまでに繰り返し振動測定を行うようにしてもよい。このとき、振動測定を行ったそれぞれの時間帯では、各時間帯の平均の回転周波数で回転していると見なして比周波数を算出し、同様の処理を行えばよい。
また、比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出した値に対して、定数を加算・減算・乗算・除算したり、比周波数の変化に対して単調増加又は減少する値を加算・減算したり、比周波数の変化に対して単調増加又は減少する値でさらに全て同符号である値を乗算・除算したりしても差し支えない。すなわち、振動ピーク位置を変化させない処理であれば本発明に包含できる。
また、比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出した値に対して、定数を加算・減算・乗算・除算したり、比周波数の変化に対して単調増加又は減少する値を加算・減算したり、比周波数の変化に対して単調増加又は減少する値でさらに全て同符号である値を乗算・除算したりしても差し支えない。すなわち、振動ピーク位置を変化させない処理であれば本発明に包含できる。
その他、本実施例では、工作機械と異常診断装置とを別体で示しているが、異常診断装置を制御装置に内蔵しても良い。
また、異常診断装置を複数の工作機械と有線或いは無線で通信可能とし、工作機械側で振動を測定してデータを取得しつつ、異常診断プログラムに基づいて異常診断方法を工作機械毎に実行するようにしてもよい。さらに、振動の測定から異常診断までを連続して行う場合に限らず、振動の測定データを記憶部に保存しておき、所定のタイミングで異常診断プログラムによる異常診断を行うようにしても差し支えない。
また、異常診断装置を複数の工作機械と有線或いは無線で通信可能とし、工作機械側で振動を測定してデータを取得しつつ、異常診断プログラムに基づいて異常診断方法を工作機械毎に実行するようにしてもよい。さらに、振動の測定から異常診断までを連続して行う場合に限らず、振動の測定データを記憶部に保存しておき、所定のタイミングで異常診断プログラムによる異常診断を行うようにしても差し支えない。
１・・主軸、２・・主軸ハウジング、３・・工具、４・・モータ、５・・速度検出器、６・・制御装置、７・・軸受、８・・表示部、９・・振動センサ、１０・・Ａ／Ｄ変換器、１１・・記憶部、１２・・演算部、１３・・表示・操作部、１００・・工作機械、２００・・異常診断装置。
Claims (14)
- 回転体を支持する転がり軸受の異常を診断する方法であって、
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別ステップと、
を実行することを特徴とする転がり軸受の異常診断方法。 - 前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする請求項１に記載の転がり軸受の異常診断方法。
- 前記比振動源回転周波数の推定は、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする請求項１に記載の転がり軸受の異常診断方法。
- 前記軸受損傷種類判別ステップは、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定する処理を含むものであり、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする請求項１乃至３の何れかに記載の転がり軸受の異常診断方法。
- 回転体を支持する転がり軸受の異常を診断する方法であって、
前記回転体の振動を複数回測定する振動測定ステップと、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出ステップと、
前記回転同期成分強調波形を入力とし、内輪傷の有無、転動体傷の有無、外輪傷の有無の少なくともひとつを出力とする教師データを用いて学習した機械学習モデルを用い、前記回転同期成分強調波形を入力して前記転がり軸受の損傷の種類を判別する軸受損傷種類判別ステップと、
を実行することを特徴とする転がり軸受の異常診断方法。 - 前記振動測定ステップでは、異なる複数の回転周波数で振動を測定することを特徴とする請求項１乃至５の何れかに記載の転がり軸受の異常診断方法。
- 前記回転同期成分強調波形算出ステップでは、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出する際に、ある比周波数の振動の大きさとして、その比周波数の前後に所定の幅を持たせて当該幅内の振動の大きさの最大値を採用することを特徴とする請求項１乃至６の何れかに記載の転がり軸受の異常診断方法。
- 回転体を支持する転がり軸受の異常を診断する装置であって、
前記回転体の振動を複数回測定する振動測定手段と、
前記振動を周波数分析して周波数ごとの振動の大きさをそれぞれ求め、前記振動の周波数を回転周波数で除算して得られる無次元量を比周波数として、各回の測定で得られた前記振動の大きさを共通の前記比周波数ごとに算出し、前記比周波数ごとの振動の大きさについて係数が全て同符号の線型和を算出し、算出した値に基づいて回転同期成分強調波形を決定する回転同期成分強調波形算出手段と、
前記回転同期成分強調波形の振動ピーク位置の規則性に基づいて、前記転がり軸受において力の発生する方向が変化する周波数である振動源回転周波数と前記回転周波数との比である比振動源回転周波数を推定し、前記比振動源回転周波数の値が１の場合には、前記転がり軸受の内輪損傷と判断し、前記比振動源回転周波数の値が前記転がり軸受の転動体の公転周波数を前記回転周波数で除算した値のとり得る範囲である場合には、前記転動体の損傷と判断し、前記比振動源回転周波数の値が０の場合には、前記転がり軸受の外輪損傷と判断する軸受損傷種類判別手段と、
を備えることを特徴とする転がり軸受の異常診断装置。 - 前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを仮定して算出される特定の比周波数に振動ピークが存在する割合である一致率を算出し、前記一致率が所定のしきい値を超過する場合には、仮定した前記比基本周波数と仮定した前記比振動源回転周波数とにより表現される規則性が前記回転同期成分強調波形にあると判断し、仮定した前記比振動源回転周波数を採用することを特徴とする請求項８に記載の転がり軸受の異常診断装置。
- 前記軸受損傷種類判別手段は、前記比振動源回転周波数の推定に当たり、前記回転同期成分強調波形において、ピーク対を構成する２つの振動ピークの間隔をピーク間距離とし、前記ピーク間距離が等しい２組のピーク対について算出したそれぞれのピーク対を構成する２つの前記振動ピークの前記比周波数の平均の比が自然数比の場合に、前記ピーク間距離の２分の１を前記比振動源回転周波数とする、又はある２つの振動ピークの前記比周波数の比が自然数比である場合に、前記比振動源回転周波数を０とすることを特徴とする請求項８に記載の転がり軸受の異常診断装置。
- 前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と前記比振動源回転周波数とを推定すると共に、前記転がり軸受の損傷の種類が判別された前記比基本周波数と前記比振動源回転周波数との組を用いて、判別された損傷の程度をさらに診断することを特徴とする請求項８乃至１０の何れかに記載の転がり軸受の異常診断装置。
- 前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、
前記比基本周波数と前記比振動源回転周波数とを合わせて表示する表示手段を備えることを特徴とする請求項８乃至１１の何れかに記載の転がり軸受の異常診断装置。 - 前記軸受損傷種類判別手段は、前記転がり軸受の損傷に対応する比基本周波数と、前記比振動源回転周波数と、前記比基本周波数と前記比振動源回転周波数とから算出される特定の比周波数とを用いて前記転がり軸受の損傷の種類ごとの前記振動ピーク位置の存在を確認するものであり、
前記特定の比周波数の位置と前記回転同期成分強調波形とを合わせて表示する表示手段を備えることを特徴とする請求項８乃至１１の何れかに記載の転がり軸受の異常診断装置。 - 所定の回転周波数で測定された回転体の振動が前記回転周波数と共に入力されたコンピュータに、請求項１乃至７の何れかに記載の転がり軸受の異常診断方法における回転同期成分強調波形算出ステップと軸受損傷種類判別ステップとを実行させることを特徴とする転がり軸受の異常診断プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019061225A JP7146683B2 (ja) | 2019-03-27 | 2019-03-27 | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019061225A JP7146683B2 (ja) | 2019-03-27 | 2019-03-27 | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2020159945A JP2020159945A (ja) | 2020-10-01 |
| JP7146683B2 true JP7146683B2 (ja) | 2022-10-04 |
Family
ID=72642984
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019061225A Active JP7146683B2 (ja) | 2019-03-27 | 2019-03-27 | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7146683B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112798280B (zh) * | 2021-02-05 | 2022-01-04 | 山东大学 | 一种滚动轴承故障诊断方法及系统 |
| CN113049252B (zh) * | 2021-03-25 | 2023-04-14 | 成都天佑路航轨道交通科技有限公司 | 一种列车轴承箱的故障检测方法 |
| WO2023037882A1 (ja) * | 2021-09-09 | 2023-03-16 | ローム株式会社 | 機械学習装置 |
| JP2023096572A (ja) * | 2021-12-27 | 2023-07-07 | 株式会社神戸製鋼所 | 転がり軸受異常検出装置および転がり軸受異常検出方法 |
| CN116026598B (zh) * | 2023-03-30 | 2023-07-14 | 山东梁轴科创有限公司 | 一种轴承振动检测系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001021453A (ja) | 1999-07-09 | 2001-01-26 | Nsk Ltd | 軸受の異常診断方法および異常診断装置 |
| JP2011259624A (ja) | 2010-06-09 | 2011-12-22 | Fuji Electric Co Ltd | 転がり軸受部振動データの高周波電磁振動成分除去方法および高周波電磁振動成分除去装置、回転機械の転がりの軸受診断方法および軸受診断装置 |
| JP2012026481A (ja) | 2010-07-21 | 2012-02-09 | Hitachi Constr Mach Co Ltd | 軸受の損傷検出方法 |
| JP2017101954A (ja) | 2015-11-30 | 2017-06-08 | 日本精工株式会社 | 機械設備の評価方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3373718B2 (ja) * | 1996-03-21 | 2003-02-04 | 松下電工株式会社 | 軸受の傷検出方法 |
-
2019
- 2019-03-27 JP JP2019061225A patent/JP7146683B2/ja active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001021453A (ja) | 1999-07-09 | 2001-01-26 | Nsk Ltd | 軸受の異常診断方法および異常診断装置 |
| JP2011259624A (ja) | 2010-06-09 | 2011-12-22 | Fuji Electric Co Ltd | 転がり軸受部振動データの高周波電磁振動成分除去方法および高周波電磁振動成分除去装置、回転機械の転がりの軸受診断方法および軸受診断装置 |
| JP2012026481A (ja) | 2010-07-21 | 2012-02-09 | Hitachi Constr Mach Co Ltd | 軸受の損傷検出方法 |
| JP2017101954A (ja) | 2015-11-30 | 2017-06-08 | 日本精工株式会社 | 機械設備の評価方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2020159945A (ja) | 2020-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7146683B2 (ja) | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム | |
| JP6507297B1 (ja) | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム | |
| Antoni et al. | Feedback on the Surveillance 8 challenge: Vibration-based diagnosis of a Safran aircraft engine | |
| Su et al. | Induction machine condition monitoring using neural network modeling | |
| US11099101B2 (en) | Method for estimating bearing fault severity for induction motors | |
| JP5930789B2 (ja) | 異常音診断装置 | |
| Wang et al. | An online tacholess order tracking technique based on generalized demodulation for rolling bearing fault detection | |
| Choqueuse et al. | Current frequency spectral subtraction and its contribution to induction machines’ bearings condition monitoring | |
| CN107076155B (zh) | 用于通过手持通信装置检测泵组件中的故障的方法和系统 | |
| Schmidt et al. | A discrepancy analysis methodology for rolling element bearing diagnostics under variable speed conditions | |
| JP2008292288A (ja) | 減速機の軸受診断装置 | |
| EP3147681B1 (en) | Rotary machine diagnostic system | |
| JP6906985B2 (ja) | 振動診断システム、振動診断方法及びパラメータ設定方法 | |
| JPWO2020208743A1 (ja) | 電動機設備の異常診断装置、電動機設備の異常診断方法、および電動機設備の異常診断システム | |
| KR20180073500A (ko) | 잔여 공간 내의 진동 역학의 모델링 및 시각화 | |
| JP6511573B1 (ja) | 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム | |
| JP2020183939A (ja) | ギアボックスの異常検出方法及び情報処理装置 | |
| Bouaouiche et al. | Detection of defects in a bearing by analysis of vibration signals | |
| JP6601376B2 (ja) | 回転体の検査装置、および、回転体の検査方法 | |
| WO2024075345A1 (en) | System and method for eccentricity severity estimation of induction machines using a sparsity-driven regression model | |
| WO2019220542A1 (ja) | 圧延設備の診断装置及び診断方法 | |
| WO2021054165A1 (ja) | 設備診断装置、設備診断方法、及び設備診断プログラム | |
| RU2451299C1 (ru) | Устройство диагностирования в реальном времени системы электродвижения судна | |
| CN112557036B (zh) | 诊断装置与方法以及计算机可读取存储媒介 | |
| Abdisa et al. | A neural network based motor bearing fault diagnosis algorithm and its implementation on programmable logic controller |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211012 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220810 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220823 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220921 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7146683 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |