JP4312477B2 - 回転機械の診断方法、診断装置及びそのプログラム - Google Patents
回転機械の診断方法、診断装置及びそのプログラム Download PDFInfo
- Publication number
- JP4312477B2 JP4312477B2 JP2003067071A JP2003067071A JP4312477B2 JP 4312477 B2 JP4312477 B2 JP 4312477B2 JP 2003067071 A JP2003067071 A JP 2003067071A JP 2003067071 A JP2003067071 A JP 2003067071A JP 4312477 B2 JP4312477 B2 JP 4312477B2
- Authority
- JP
- Japan
- Prior art keywords
- vibration
- dimensional
- dimensionless
- principal component
- evaluation index
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、回転機械の振動波形データを解析することによって回転機に生じる異常モードを検知する技術に関し、特に異常発生初期の軽微な段階で異常を検知するための回転機械の診断方法に関する。
【０００２】
【従来の技術】
従来、回転機械の良否を判定する方法としては以下の技術が知られている。
【０００３】
第１の技術は、振動測定値を用いて回転機械の良否を判定するものである。
【０００４】
具体的には、その回転機械に発生する振動変位、振動速度、振動加速度値などを計測し、その計測値から振動の最大値や実効値などのパラメータを求め、それらのパラメータ値とＩＳＯ基準やユーザが設定した基準値とを比較することにより、正常、注意、危険の判定を行う。この技術は一般的な手法として広く知られているものである。そして、この一般的な手法をさらに発展させた、以下の技術が公知となっている。
【０００５】
第２の技術は、有次元パラメータと無次元パラメータを用いて回転機械の良否を判定するものである。
【０００６】
ここで、「有次元パラメータ」とは、測定した振動の大きさ等を所定の単位系によって表したパラメータであり、「無次元パラメータ」とは、測定した振動の波形の特徴を無単位の指標として表したパラメータである。
【０００７】
有次元パラメータは、異常の程度が重くなったときに有効である。従って、初期異常、または一部の異常モードに対しては有次元パラメータが十分に大きくならない、という欠点がある。また、異常とは無関係な運転条件（回転数や負荷など）が変化したときにも、有次元パラメータの値が影響を受けるという欠点がある。しかし、多くの機械設備では運転の継続が可能か否かは、有次元パラメータによって決定されることから有次元パラメータは必要不可欠なパラメータである。
【０００８】
一方、無次元パラメータは運転条件の影響を受けにくいという特徴を有するが、一部の異常に対しては、異常の程度が重くなった場合でも、その異常に反応しないことがある。
【０００９】
従って第２の技術では、有次元パラメータ、無次元パラメータの両方の長所短所を生かした判定方法を提案している。
【００１０】
具体的には、有次元パラメータ１個と無次元パラメータ１４個を採取する。そして、有次元パラメータはその大きさ（実効値または最大値）が正常値からどれくらい変化したかの変化量によって正常、注意、危険の判定を行い、無次元パラメータも１４個のそれぞれのパラメータに対して同様に正常値からの変化量によって正常、注意、危険の判定を行っている。そして、有次元パラメータと無次元パラメータのそれぞれの判定結果を組み合わせて総合的に正常、注意、危険の判定を行う。
【００１１】
第３の技術は、判定モデルを用いて回転機械の異常を判定するものである（例えば、特許文献１参照）。
【００１２】
この技術では、機器の正常状態において動作変量を計測してそのデータから異常状態における変量データを推定する。そして正常状態、異常状態それぞれについてＡＲモデルを適用し、ＡＲモデルをもとに、それぞれのマハラノビスの距離を求めておく。以上の準備の下に、診断対象機器について実際に変量データを計測してＡＲモデルを適用し、マハラノビスの距離を対比してＡＲモデルどうしの近さをもとに正常か異常かを判定する。
【００１３】
第４の技術は、振幅確率密度関数を用いて回転機械の異常を判定するものである（例えば、特許文献２参照）。
【００１４】
この技術は、正常状態の機械から発生する振動等の計測信号の振幅確率密度関数が正規分布であり、故障や異常状態になると正規分布から外れることを基本原理としたものである。即ち、正規化した振幅確率密度関数からグラムチャーリー級数を算出し、正規分布との差分を算出して異常か否かを判定する。本発明の従来の手法と比較した顕著な特徴は、計測信号の大きさや周波数成分の情報は利用していないことであり、上述の判定方法を用いることによって多種多様な回転機械に共通に使える判定基準値を提案しているところにある。
【００１５】
【特許文献１】
特開２０００−４６７０号公報
【００１６】
【特許文献２】
特開２０００−１７１２９１号公報
【００１７】
【発明が解決しようとする課題】
しかしながら、上述した従来の技術では、なお解決すべき以下のような問題点を有していた。
【００１８】
先ず、一般手法を用いる第１の技術においては、振動速度、加速度値の最大値や実効値といった有次元パラメータのみを用いるものであるため、回転機械に発生するアンバランスやミスアライメント、軸受の欠陥などの異常を検出することができるが、これはあくまでも損傷の程度がある程度大きな場合に限られる。即ち初期の異常などはほとんど検知できないという問題がある。
【００１９】
有次元パラメータと無次元パラメータを用いて回転機械の良否を判定する第２の技術においては、有次元パラメータを１個しか採用していないため、十分に異常の特徴を抽出しているとはいえない。このため、その有次元パラメータで検出感度が低い異常に対してはその兆候を見逃す危険性があった。特に初期の異常に対してはその恐れが大きいと考えられる。
【００２０】
さらに、どのパラメータが判定に有効かはユーザがデータを見て決定しなければならない。従って、ユーザが注目すべきパラメータを正しく選定していない場合には、異常の兆候を見逃す危険性がある。
【００２１】
判定モデルを用いて回転機械の異常を判定する第３の技術（特許文献１記載の技術）においては、機器の異常を客観的に高い確度で判断することのできる診断方法とあるが、動作変量が適正に選択されていなければ所定の異常を高い精度で検出することは困難であると考えられる。特に異常発生初期の軽微な段階で異常を検知するための技術に関しては、そのような記述もなくまた示唆もされていない。
【００２２】
振幅確率密度関数を用いて回転機械の異常を判定する第４の技術（特許文献２記載の技術）においては、機器ごとに煩わしい基準値を設定することなく、異常有無の判断ができる診断方法であるが、これはあくまでその異常が振幅確率密度関数に影響を与えることを前提とするものである。従って、本技術は計測信号がその異常状態を検出している場合の異常判定方法であり、通常の一般的な計測信号を用いて振幅確率密度関数が異常発生初期の軽微な段階で異常を検知できるとの根拠は記載されていない。
【００２３】
以上のように、従来技術では正常状態から異常状態に差しかかった異常初期の軽微な段階での異常を十分に検出できるには至っていないという問題点が指摘されていた。異常の早期検知は、回転機械の信頼性を向上させる上で重要な課題であり、特に、高い信頼性、安全性が要求される重要な回転機械に対しては必要不可欠な課題であった。
【００２４】
本発明は、上述の問題点を解決するためになされたもので、異常を初期の段階で早期に検知することができる回転機械の診断方法及び回転機械の診断プログラムを提供することを目的とする。
【００２５】
【課題を解決するための手段】
上記課題を解決するための本発明は、回転機械の作動状態における振動情報を採取し、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出し、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出し、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行い、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断方法である。
【数７】
【００２６】
また本発明は、回転機械の作動状態における振動情報を採取し、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出し、前記統合次元の主成分に基く統合次元の状態評価指数を算出し、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行い、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断方法である。
【数８】
【００２７】
また本発明は、上記記載の発明である回転機械の診断方法において、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断方法である。
【００２８】
また本発明は、上記記載の発明である回転機械の診断方法において、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する回転機械の診断方法である。
【００２９】
また本発明は、回転機械の診断プログラムにおいて、コンピュータに、回転機械の作動状態における振動情報を採取する手順、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手順、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手順、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手順、を実行させ、前記状態評価指数が以下の式により導き出されるプログラムである。
【数９】
【００３０】
また本発明は、回転機械の診断プログラムにおいて、コンピュータに、回転機械の作動状態における振動情報を採取する手順、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手順、前記統合次元の主成分に基く統合次元の状態評価指数を算出する手順、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手順、を実行させ、前記状態評価指数が以下の式により導き出されるＳであるプログラムである。
【数１０】
【００３１】
また本発明は、上記記載の発明であるプログラムにおいて、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断プログラムである。
【００３２】
また本発明は、上記記載の発明であるプログラムにおいて、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する異常原因推定手順、を備えた回転機械の診断プログラムである。
【００３３】
また本発明は、回転機械の作動状態における振動情報を採取する手段と、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手段と、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手段と、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手段とを備え、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断装置である。
【数１１】
【００３４】
また本発明は、回転機械の作動状態における振動情報を採取する手段と、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手段と、前記統合次元の主成分に基く統合次元の状態評価指数を算出する手段と、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手段とを備え、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断装置である。
【数１２】
【００３５】
また本発明は、上記記載の発明である回転機械の診断装置において、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断装置である。
【００３６】
また本発明は、上記記載の発明である回転機械の診断装置において、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する手段を更に供えた回転機械の診断装置である。
【００４０】
【発明の実施の形態】
図１は、本発明に係る診断システムの概略の構成を示すブロック図である。 本発明に係る診断システムは、異常判定を行うための諸機能を備えた判定処理装置１、振動センサ２からの電荷信号を電圧信号等に変換する信号変換器３、記録媒体に記録されている振動信号を再生する記録信号再生装置４、判定処理装置１からの情報を表示する表示装置５及び判定処理装置１からの情報を印刷する印刷装置６で構成されている。
【００４１】
そして判定処理装置１は、入出力制御部１０、波形読込部１１、ＦＦＴ処理部１２、パラメータ算出部１３、主成分分析部１４、統合パラメータ値算出部１５、良否判定処理部１６、異常原因分析処理部１７、判定結果処理部１８及びデータ部２０とで構成されている。
【００４２】
入出力制御部１０は、外部装置との信号授受を行い情報交換を行うためのインターフェースである。波形読込部１１は振動センサ２または記録信号再生装置４から振動信号を読み込むとともに必要時にはその信号から加速度信号と速度信号を生成する。ＦＦＴ処理部１２は読み込んだ振動信号の周波数分析を行う。
【００４３】
パラメータ算出部１３は、波形読込部１１とＦＦＴ処理部１２の処理結果に基づいて有次元パラメータと無次元パラメータを算出する。主成分分析部１４は、正常状態にある回転機械の振動信号に基づいて主成分分析を行い固有ベクトルを算出する。統合パラメータ値算出部１５は、有次元パラメータ、無次元パラメータ及び固有ベクトルとから回転機械の異常を判定するための指標である統合パラメータ値を算出する。
【００４４】
良否判定処理部１６は、統合パラメータ値に基づいて回転機械の良否を判定する。異常原因分析処理部１７は、回転機械が異常であると判定された場合にその異常の原因を分析して推定する。判定結果処理部１８は、良否判定結果及び異常原因推定結果を表示装置５あるいは印刷装置６に出力する。
【００４５】
データ部２０は、上述の各処理部が異常判定処理を行うために用いる各種データを記憶する。初期状態データベース２１は、回転機械の正常状態において算出した各種パラメータ値などを記憶する。測定値データベース２２は、回転機械の現在状態において算出した各種パラメータ値などを記憶する。基準値データベース２３は、回転機械の良否を判定するための基準値を記憶する。異常原因分析データベース２４は、異常原因を推定するための基準値データを記憶する。
【００４６】
尚、本発明に係る診断システムでは、振動信号を振動センサ２から直接に、あるいは記録信号再生装置４から読み込んでいるが、この形態に限定されず例えば図示しない通信回線を介して遠隔から振動信号を読み取るように構成しても良い。また波形読込部１１、ＦＦＴ処理部１２を判定処理装置の内部に設けず外部の装置を用いて構成し、その処理結果を入力するように構成しても良い。
【００４７】
次に、このように構成された診断システムの動作について説明する。
図２は、診断システムの概略の動作を示すフロー図である。
【００４８】
まず、波形読込部１１が、機械の初期状態（正常状態）における振動波形を読み込む（Ｓ１）。ここで、波形読込部１１は振動加速度波形ならびに振動速度波形を読み込むように構成しても良く、また振動加速度波形を読み込んだ後、その信号を積分して振動速度波形を算出するように構成しても良い。
【００４９】
次に、ＦＦＴ処理部１２が読み込んだ振動信号の周波数分析を実行し（Ｓ２）、各周波数毎の成分値を算出する。
【００５０】
続いて、パラメータ算出部１３が起動して、それぞれの波形データから以下に示すような複数の有次元パラメータと無次元パラメータを算出する（Ｓ３）。
【００５１】
（１）有次元パラメータ
１）振動速度
ａ.速度ピーク値 ：ＶＥＬ−Ｐ
測定した振動速度波形の振幅値ｘｉの内、｜ｘｉ｜の大きなものから数えた上位5％の｜ｘｉ｜の平均値である。
【００５２】
ｂ．速度ＲＭＳ値 ：ＶＥＬ−Ｒ
【数５】
【００５３】
ｃ．周波数成分値 ：ＶＥＬ−ｆｒ（回転周波数ｆｒ成分）
ここで回転周波数ｆｒは、Ｎ:回転数［rpm］を用いて下式で表される。
【００５４】
【数６】
【００５５】
ｄ．周波数成分値 ：ＶＥＬ−２ｆｒ（周波数２×ｆｒ成分）
ｅ．周波数成分値 ：ＶＥＬ−３ｆｒ（周波数３×ｆｒ成分）
ｆ．周波数成分値 ：ＶＥＬ−４ｆｒ（周波数４×ｆｒ成分）
ｇ．周波数成分値 ：ＶＥＬ−５ｆｒ（周波数５×ｆｒ成分）
ｈ．周波数成分値 ：ＶＥＬ−１／２ｆｒ（周波数１／２×ｆｒ成分）
ｉ．周波数成分値 ：ＶＥＬ−１／３ｆｒ（周波数１／３×ｆｒ成分）
ここで表された有次元パラメータの内、ｃの回転周波数成分は、主にアンバランス異常時に卓越する成分、ｄ〜ｇの回転周波数の高調波成分は、ミスアライメントやガタ・ゆるみ異常時に発生する成分、ｈ〜ｉの分数調波成分は、ガタ・ゆるみ発生時に表れる特徴的な成分である。
【００５６】
２）振動加速度
ａ．加速度ピーク値：ＡＣＣ−Ｐ
測定した振動加速度波形の振幅値ｘｉの内、｜ｘｉ｜の大きなものから数えた上位5％の｜ｘｉ｜の平均値である。
【００５７】
ｂ．加速度ＲＭＳ値：ＡＣＣ−Ｒ
【数７】
【００５８】
ｃ．周波数成分値 ：ＡＣＣ−ｆ０（外輪キズ周波数成分）
【数８】
【００５９】
ｄ．周波数成分値 ：ＡＣＣ−ｆｉ（内輪キズ周波数成分）
【数９】
【００６０】
ｅ．周波数成分値 ：ＡＣＣ−ｆｂ（ボールキズ周波数成分）
【数１０】
【００６１】
ｆ．周波数成分値 ：ＡＣＣ−２ｆｂ（ボールキズ周波数成分）
【数１１】
【００６２】
ｇ．周波数成分値 ：ＡＣＣ−ｆｃ（保持器キズ周波数成分）
【数１２】
【００６３】
ｈ．周波数成分値：ＡＣＣ−ｆｒ（回転周波数成分）
【数１３】
【００６４】
ここで、Ｎ:回転数［rpm］、ｆｒ：軸（内輪）の回転周波数(Hz)、Ｄ：軸受のピッチ円直径(ｍｍ)、d：転動体の直径(ｍｍ)、α：接触角(度)、z：転動体の数である。図３は一般的な軸受の主要諸元を示す図である。
【００６５】
（２）無次元パラメータ
１）振動速度
ａ．波形率Ｓｆ：ＶＥＬ−Ｓｆ
振動波形の正弦波からのずれを表すパラメータ。低周波領域のアンバランスやミスアライメント、脈動波形などの定量化に有効なパラメータ。
【００６６】
【数１４】
【００６７】
ｂ．波高率Ｃｆ：ＶＥＬ−Ｃｆ
振動波形の衝撃性を表すパラメータ。局部異常の検出に有効なパラメータ。
【００６８】
【数１５】
【００６９】
ｃ．衝撃指数Ｉｐ：ＶＥＬ−Ｉｐ
波高率Ｃｆと同様に衝撃性を表すパラメータ。
【００７０】
【数１６】
【００７１】
ｄ．スキューネス（歪み度）β１：ＶＥＬ−β１
振動波形がゼロ点を中心にしていかに非対称となっているかを示すパラメータ。摩耗系の異常が発生すると、振動波形が非対称となり、歪み度が増大する。
【００７２】
【数１７】
【００７３】
ｅ．クートシス（尖り度）β２：ＶＥＬ−β２
振動波形がゼロ点を中心にしていかに尖っているかを示すパラメータ。転がり軸受や歯車装置の異常診断に有効なパラメータ。
【００７４】
【数１８】
【００７５】
２）振動加速度
振動速度と同様に、以下のパラメータを求める。
【００７６】
ａ．波形率Ｓｆ ：ＡＣＣーＳｆ
ｂ．波高率Ｃｆ ：ＡＣＣ−Ｃｆ
ｃ．衝撃指数Ｉｐ ：ＡＣＣ−Ｉｐ
ｄ．スキューネス（歪み度）β１：ＡＣＣ−β１
ｅ．クートシス（尖り度）β２：ＡＣＣ−β２
以上の有次元パラメータと無次元パラメータを算出後、主成分分析部１４が起動して算出された有次元パラメータあるいは無次元パラメータ毎に、主成分分析を行う（Ｓ４）。
【００７７】
主成分算出までの計算手順を以下に説明する。尚、以下の説明では簡便化のために有次元パラメータのみについて説明するが、無次元パラメータについても同様の手順とする。
【００７８】
ある設備の正常状態下で上述の手順により収集されたｍ個の有次元パラメータを要素とする兆候パラメータＹｐ＝（Ｙ１，Ｙ２，Ｙ３ ・・・，Ｙｍ）のｎ組のデータからなる行列Ｙ０を以下の式で定義する。
【００７９】
【数１９】
【００８０】
次に、この行列Ｙ０の列方向、即ち縦方向の行列要素に対して以下の式を用いて演算を行い、ｙｉを算出する。
【００８１】
【数２０】
【００８２】
この演算を列毎に行うことによって、新たな標準化された兆候パラメータを要素とするデータ行列Ｙを求める。
【００８３】
【数２１】
【００８４】
この行列は各列毎に平均値＝０、分散＝１に変換された行列である。
【００８５】
そこで、このデータ行列Ｙから相関行列Ｒを算出すると以下の式となる。
【００８６】
【数２２】
【００８７】
ここに、ｒijは２つの列の相関係数である。すなわち、
【数２３】
【００８８】
次に、相関行列Ｒの固有値λを求める。即ち、以下の式を満たす固有値λを求める。
【００８９】
【数２４】
【００９０】
求められた固有値をλ１,λ２，・・・λｎとする。ただし、λ１＞λ２，・・・＞λｎである。これらのｎ個の固有値に対する固有ベクトルをａｉ（ｉ＝１、・・・ｎ）とする。
【００９１】
【数２５】
【００９２】
そうするとそれぞれの固有ベクトルに対して次の式が成立する。
【００９３】
【数２６】
【００９４】
この式に基づいて、ａｉ１、ａｉ２、・・・ａｉｍを求める。但し、この係数は以下の式を充たす値である。
【００９５】
【数２７】
【００９６】
この手順を繰り返して、夫々の固有値λ１,λ２，・・・λｎに対する固有ベクトルａ１，ａ２，・・・ａｎを求めることができる。そしてこれらのｎ個の固有ベクトルと標準化された兆候パラメータ（ｙ１、ｙ２、・・・ｙｍ）を組み合わせて、主成分Ｚ１，Ｚ２，・・・Ｚｎが、以下の式で表される。
【００９７】
【数２８】
【００９８】
ここで、Ｚ１を第１主成分、Ｚ２を第２主成分、Ｚｎを第ｎ主成分と呼ぶ。
【００９９】
次に、統合パラメータ値算出部１５が起動して、有次元統合パラメータＳｙを算出する（Ｓ５）。
【０１００】
前述の正常状態下で求めた主成分Ｚｉ（母集団）が、正規分布に従うと仮定する。この母集団から独立に取り出されたｎ個の標本で構成される統計量χ２ は、自由度ｎ−１のカイ２乗分布に従う。
【０１０１】
【数２９】
【０１０２】
ここで、母集団から取り出したｎ個の標本をＸ１，Ｘ２，・・・Ｘｎとすると、標本分散ｓ２は次の式で表される。
【０１０３】
【数３０】
【０１０４】
また、主成分Ｚｉの母分散σ２は固有値λiに等しいことから、
σ２＝ λi
これらの関係を整理すると以下の式となり、標準化されたＸ値の偏差平方和は、自由度ｎー１のカイ自乗分布に従う。
【０１０５】
【数３１】
【０１０６】
ここで、式のＸｉを主成分Ｚｉに置き換えると、
【数３２】
【０１０７】
次に、正常状態下におけるデータの主成分Ｚｉが１−αの確率で入る領域は、以下の式で表される。
【０１０８】
【数３３】
【０１０９】
よって、正常状態の状態確定領域は次の式を満たす範囲となる。
【０１１０】
【数３４】
【０１１１】
たとえば、有意水準α＝0.05、自由度φ＝3の場合は、
【数３５】
【０１１２】
となり、この正常状態確定領域に入ったときは正常、領域外のときは異常という判定が可能となる。
【０１１３】
従って、正常状態からの変化量を監視するために、下の式で表される状態量を統合パラメータとして定義する。
【０１１４】
【数３６】
【０１１５】
このＳ値が大きくなれば異常状態と判定でき、回転機械の良否を判定するための普遍的な状態量を示すものである。この式において、有次元パラメータについて集約した指標値を有次元統合パラメータ値Ｓｙとし、無次元パラメータについて集約した指標値を無次元統合パラメータ値Ｓｍとする。そして、この２つの状態量Ｓｙ、Ｓｍを監視することにより、設備の状態監視を行う。
【０１１６】
ＳｙとＳｍの２つを監視する理由は、異常モードによって、有次元パラメータに大きな変化が現れる場合と、無次元パラメータに大きな変化が現れる場合があるためである。
【０１１７】
以上の手順に従って、初期状態波形データについて有次元パラメータ、無次元パラメータ、主成分値、固有値、固有ベクトルなどが求められるが、それらの値は正常状態における基準データとして初期状態データベース２１に記憶する（Ｓ６）。
【０１１８】
以上説明した波形データ処理は正常状態における基準データを作成するための準備である。したがって、この正常状態処理が行われた時点以降に読み込まれる、回転機械の良否を判定するための振動波形の処理については、この基準データを利用する形態での処理となる。
【０１１９】
正常状態処理が行われた時点以降に読み込まれる判定用振動波形の処理においては、ステップＳ１からＳ３までは前述と同様に実行される。しかし、ステップＳ４の主成分分析においては、標準化された兆候パラメータを要素とするデータ行列Ｙを求めた後は、再度新たな固有ベクトルの算出は行わず、初期状態データベース２１に記憶されている固有ベクトルを抽出して、その固有ベクトルに基づいて主成分Ｚ１、・・・Ｚｎを算出する。そして、この主成分に基づいて前述のステップＳ５と同様の手順で有次元統合パラメータＳｙと無次元統合パラメータＳｍを算出する。このようにして算出された判定用振動波形の処理データおよび処理結果は、測定値データベース２２に記憶される。
【０１２０】
次に、良否判定処理部１６が起動して算出された有次元統合パラメータＳｙと無次元統合パラメータＳｍに基づいて回転機械の良否を判定する（Ｓ７）。良否判定処理部１６は、このＳｙまたはＳｍの値が、初期状態よりもｎ倍になれば注意、ｍ倍になれば異常、ｎ倍未満であれば良と判定する。この判定に際し、初期状態のＳｙまたはＳｍは初期状態データベース２１から取り出され、判定基準であるｎ倍またはｍ倍は基準値データベース２３に格納されている値が用いられる。
【０１２１】
判定結果で、注意または異常と判定された場合は、異常原因分析処理部１７が起動して異常原因を分析して推定する（Ｓ８）。異常原因分析処理部１７は異常原因分析データベース２４から、図４に示す異常原因分析マトリックスを取り出す。
【０１２２】
異常原因分析マトリックスは、縦方向に異常原因を列記し、横方向に有次元及び無次元パラメータを配置した構成である。そして、この異常原因分析マトリックスの行列要素には、異常度指数１（ａｉｊ）と異常度指数２（ｂｉｊ）が記載されている。この異常度指数１と異常度指数２は、該異常が発生した場合にそのパラメータが示す値を統計的に処理した指標であり、過去に発生した異常に実績データに基づいて算出し決定した値である。ここで、異常度指数１と異常度指数２の複数の指数を設けているのは、異常原因によってはいずれか一方のみの指数しか変化しない場合があるためである。尚、その異常度指数が該異常に関係しない場合は、その異常度指数には数値は設定されず、以降の計算には組み込まれない。
【０１２３】
次に、異常原因分析処理部１７は、測定値データベース２２に格納されている判定用振動波形の処理データを用いて、図５に示す異常度指数１（Ａｊ）と異常度指数（Ｂｊ）とを算出する。そして、次の式に基づいて異常原因（ｉ）についてパラメータ（ｊ）毎に異常確率指数Ｐｉｊを算出する。
【０１２４】
Ｐｉｊ ＝ αｊ＊｛（Ａｊ／ａｉｊ）＋（Ｂｊ／ｂｉｊ）｝
但し、ａｉｊ、ｂｉｊともに値が設定されているとき：αｊ＝０．５
ａｉｊ、ｂｉｊいずれか一方に値が設定されているとき：αｊ＝１
ａｉｊ、ｂｉｊともに値が設定されていないとき：αｊ＝０
そして、次の式によって異常原因（ｉ）毎に総合確率Ｔｉを算出する。
【０１２５】
Ｔｉ ＝１／Ｎ＊ΣＰｉｊ
ここで、Ｎ：ａｉｊまたはｂｉｊに値が設定されているパラメータの数
この式において、総合確率Ｔｉが大きな値をとる場合は、あらかじめ異常原因分析マトリックスに設定されてある異常原因と類似している場合であると考えられる。従って、総合確率Ｔｉが所定以上になれば、その異常が発生している可能性が高いと判断する。
【０１２６】
続いて、判定結果処理部１８が起動して、判定結果である良、注意、異常と、異常原因（ｉ）毎の総合確率Ｔｉを編集して表示装置５あるいは印刷装置６に出力する（Ｓ９）。その後、必要なデータを保存して（Ｓ１０）、異常判定処理を終了する。
【０１２７】
次に、本発明の診断システムの適用例を、図面を参照しながら説明する。図６は、振動センサ２が取り付けられた回転機械であるモータ３０の概略図である。
【０１２８】
振動センサ２によって、正常状態における振動速度波形と加速度波形を採取し、採取した２種類の波形から、有次元パラメータと無次元パラメータの値を算出する。
【０１２９】
そして、上記の有次元・無次元パラメータに対し、それぞれ主成分分析を行い、初期波形データにおける各パラメータの値、主成分値、固有値、固有ベクトル等をすべて記憶しておく。
【０１３０】
次に、モータ３０に初期異常を発生させそのデータを計測する。具体的には軸受に軽微な欠陥のある異常データ５データを計測し、初期データ解析時と同様、有次元統合パラメータと無次元統合パラメータを算出する。図７は、正常状態と異常状態における波形データの有次元統合パラメータの値を示す図である。無次元統合パラメータについても同様に算出する。
【０１３１】
そして主成分分析を行い状態量Ｓｙ、Ｓｍ値を算出して加速度ＲＭＳ値（ＡＣＣ−Ｒ）とともに記載した結果を図８に示す。図９はＳｙ値とＡＣＣ−Ｒのデータをグラフ化した図であり、図１０はＳｍ値とＡＣＣ−Ｒのデータをグラフ化した図である。
【０１３２】
Ｎｏ．１１〜１５の異常データに対し、有次元パラメータである加速度ＲＭＳ値(ＡＣＣ−Ｒ)はほとんど変化していなく、異常を検知していない。これに対して、有次元統合パラメータＳｙ値と無次元統合パラメータＳｍ値は５倍以上も大きく変化しており、Ｓｙ値とＳｍ値を用いることで初期異常を感度良く検知していることがわかる。
【０１３３】
尚、有次元統合パラメータＳｙと無次元統合パラメータＳｍ単独で判定せずに、ＳｙとＳｍを変数とする関数を考え、その関数の示す値によって判定するように構成しても良い。
【０１３４】
図１１は、Ｓｙ値とＳｍ値と積Ｓｙ＊Ｓｍのデータをグラフ化した図である。Ｎｏ．１１〜１５の異常データに対し、指標Ｓｙ＊Ｓｍが初期異常の兆候をより顕著に表していることがわかる。このように、Ｓｙ値とＳｍ値を組合わせた指標を用いることで、初期異常の発生を確実に検出することが可能となる。
【０１３５】
また、本実施の形態では、有次元統合パラメータＳｙと無次元統合パラメータＳｍの２つの指標を求めているが、この形態に限定されず有次元パラメータと無次元パラメータとを分離せずにまとめ、これから統合パラメータを算出するように構成しても良い。
【０１３６】
また、有次元パラメータと無次元パラメータの選択を異常モード毎に行い、異常モード毎の有次元統合パラメータと無次元統合パラメータを用いるように構成しても良い。
【０１３７】
更に、異常モードの初期、中期、末期の各状態毎に、有次元パラメータと無次元パラメータを選択し、その状態毎の有次元統合パラメータと無次元統合パラメータを用いて判定するように構成しても良い。
【０１３８】
以上のように、この発明によれば、振動波形から算出した２つの統合パラメータ値を監視する機能により、回転機械の劣化状態を的確に、かつ早期に検知することができる。また、初期データからの波形の大きさと形状から統合パラメータを算出するので、回転機械の型式に依存しない精度の高い異常判定が可能となる。従って従来手法では困難であった特殊回転機械の診断も可能となるなど、信頼性の高い診断が可能となり、産業上有用な効果がもたらされる。
【０１３９】
【発明の効果】
この発明における診断手法を用いて回転機械の状態監視を行うことにより、異常の早期検知が可能となり、精度の高い良否判定を行うことができる。
【図面の簡単な説明】
【図１】 本発明に係る診断システムの概略の構成を示すブロック図。
【図２】 診断システムの概略の動作を示すフロー図。
【図３】 一般的な軸受の主要諸元を示す図。
【図４】 異常原因分析マトリックスを示す図。
【図５】 異常度指数を示す図。
【図６】 振動センサが取り付けられたモータを示す図。
【図７】 正常状態と異常状態における波形データの有次元統合パラメータの値を示す図。
【図８】 Ｓｙ、Ｓｍ値、加速度ＲＭＳ値を示す図。
【図９】 Ｓｙ値とＡＣＣ−Ｒの推移を示す図。
【図１０】 Ｓｍ値とＡＣＣ−Ｒの推移を示す図。
【図１１】 Ｓｙ、Ｓｍ値、Ｓｙ＊Ｓｍ値の推移を示す図。
【符号の説明】
１…判定処理装置、２…振動センサ、１１…波形読込部、１２…ＦＦＴ処理部、１３…パラメータ算出部、１４…主成分分析部、１５…統合パラメータ値算出部、１６…良否判定処理部、１７…異常原因分析処理部、１８…判定結果処理部。
【発明の属する技術分野】
本発明は、回転機械の振動波形データを解析することによって回転機に生じる異常モードを検知する技術に関し、特に異常発生初期の軽微な段階で異常を検知するための回転機械の診断方法に関する。
【０００２】
【従来の技術】
従来、回転機械の良否を判定する方法としては以下の技術が知られている。
【０００３】
第１の技術は、振動測定値を用いて回転機械の良否を判定するものである。
【０００４】
具体的には、その回転機械に発生する振動変位、振動速度、振動加速度値などを計測し、その計測値から振動の最大値や実効値などのパラメータを求め、それらのパラメータ値とＩＳＯ基準やユーザが設定した基準値とを比較することにより、正常、注意、危険の判定を行う。この技術は一般的な手法として広く知られているものである。そして、この一般的な手法をさらに発展させた、以下の技術が公知となっている。
【０００５】
第２の技術は、有次元パラメータと無次元パラメータを用いて回転機械の良否を判定するものである。
【０００６】
ここで、「有次元パラメータ」とは、測定した振動の大きさ等を所定の単位系によって表したパラメータであり、「無次元パラメータ」とは、測定した振動の波形の特徴を無単位の指標として表したパラメータである。
【０００７】
有次元パラメータは、異常の程度が重くなったときに有効である。従って、初期異常、または一部の異常モードに対しては有次元パラメータが十分に大きくならない、という欠点がある。また、異常とは無関係な運転条件（回転数や負荷など）が変化したときにも、有次元パラメータの値が影響を受けるという欠点がある。しかし、多くの機械設備では運転の継続が可能か否かは、有次元パラメータによって決定されることから有次元パラメータは必要不可欠なパラメータである。
【０００８】
一方、無次元パラメータは運転条件の影響を受けにくいという特徴を有するが、一部の異常に対しては、異常の程度が重くなった場合でも、その異常に反応しないことがある。
【０００９】
従って第２の技術では、有次元パラメータ、無次元パラメータの両方の長所短所を生かした判定方法を提案している。
【００１０】
具体的には、有次元パラメータ１個と無次元パラメータ１４個を採取する。そして、有次元パラメータはその大きさ（実効値または最大値）が正常値からどれくらい変化したかの変化量によって正常、注意、危険の判定を行い、無次元パラメータも１４個のそれぞれのパラメータに対して同様に正常値からの変化量によって正常、注意、危険の判定を行っている。そして、有次元パラメータと無次元パラメータのそれぞれの判定結果を組み合わせて総合的に正常、注意、危険の判定を行う。
【００１１】
第３の技術は、判定モデルを用いて回転機械の異常を判定するものである（例えば、特許文献１参照）。
【００１２】
この技術では、機器の正常状態において動作変量を計測してそのデータから異常状態における変量データを推定する。そして正常状態、異常状態それぞれについてＡＲモデルを適用し、ＡＲモデルをもとに、それぞれのマハラノビスの距離を求めておく。以上の準備の下に、診断対象機器について実際に変量データを計測してＡＲモデルを適用し、マハラノビスの距離を対比してＡＲモデルどうしの近さをもとに正常か異常かを判定する。
【００１３】
第４の技術は、振幅確率密度関数を用いて回転機械の異常を判定するものである（例えば、特許文献２参照）。
【００１４】
この技術は、正常状態の機械から発生する振動等の計測信号の振幅確率密度関数が正規分布であり、故障や異常状態になると正規分布から外れることを基本原理としたものである。即ち、正規化した振幅確率密度関数からグラムチャーリー級数を算出し、正規分布との差分を算出して異常か否かを判定する。本発明の従来の手法と比較した顕著な特徴は、計測信号の大きさや周波数成分の情報は利用していないことであり、上述の判定方法を用いることによって多種多様な回転機械に共通に使える判定基準値を提案しているところにある。
【００１５】
【特許文献１】
特開２０００−４６７０号公報
【００１６】
【特許文献２】
特開２０００−１７１２９１号公報
【００１７】
【発明が解決しようとする課題】
しかしながら、上述した従来の技術では、なお解決すべき以下のような問題点を有していた。
【００１８】
先ず、一般手法を用いる第１の技術においては、振動速度、加速度値の最大値や実効値といった有次元パラメータのみを用いるものであるため、回転機械に発生するアンバランスやミスアライメント、軸受の欠陥などの異常を検出することができるが、これはあくまでも損傷の程度がある程度大きな場合に限られる。即ち初期の異常などはほとんど検知できないという問題がある。
【００１９】
有次元パラメータと無次元パラメータを用いて回転機械の良否を判定する第２の技術においては、有次元パラメータを１個しか採用していないため、十分に異常の特徴を抽出しているとはいえない。このため、その有次元パラメータで検出感度が低い異常に対してはその兆候を見逃す危険性があった。特に初期の異常に対してはその恐れが大きいと考えられる。
【００２０】
さらに、どのパラメータが判定に有効かはユーザがデータを見て決定しなければならない。従って、ユーザが注目すべきパラメータを正しく選定していない場合には、異常の兆候を見逃す危険性がある。
【００２１】
判定モデルを用いて回転機械の異常を判定する第３の技術（特許文献１記載の技術）においては、機器の異常を客観的に高い確度で判断することのできる診断方法とあるが、動作変量が適正に選択されていなければ所定の異常を高い精度で検出することは困難であると考えられる。特に異常発生初期の軽微な段階で異常を検知するための技術に関しては、そのような記述もなくまた示唆もされていない。
【００２２】
振幅確率密度関数を用いて回転機械の異常を判定する第４の技術（特許文献２記載の技術）においては、機器ごとに煩わしい基準値を設定することなく、異常有無の判断ができる診断方法であるが、これはあくまでその異常が振幅確率密度関数に影響を与えることを前提とするものである。従って、本技術は計測信号がその異常状態を検出している場合の異常判定方法であり、通常の一般的な計測信号を用いて振幅確率密度関数が異常発生初期の軽微な段階で異常を検知できるとの根拠は記載されていない。
【００２３】
以上のように、従来技術では正常状態から異常状態に差しかかった異常初期の軽微な段階での異常を十分に検出できるには至っていないという問題点が指摘されていた。異常の早期検知は、回転機械の信頼性を向上させる上で重要な課題であり、特に、高い信頼性、安全性が要求される重要な回転機械に対しては必要不可欠な課題であった。
【００２４】
本発明は、上述の問題点を解決するためになされたもので、異常を初期の段階で早期に検知することができる回転機械の診断方法及び回転機械の診断プログラムを提供することを目的とする。
【００２５】
【課題を解決するための手段】
上記課題を解決するための本発明は、回転機械の作動状態における振動情報を採取し、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出し、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出し、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行い、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断方法である。
【数７】
また本発明は、回転機械の作動状態における振動情報を採取し、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出し、前記統合次元の主成分に基く統合次元の状態評価指数を算出し、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行い、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断方法である。
【数８】
また本発明は、上記記載の発明である回転機械の診断方法において、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断方法である。
【００２８】
また本発明は、上記記載の発明である回転機械の診断方法において、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する回転機械の診断方法である。
【００２９】
また本発明は、回転機械の診断プログラムにおいて、コンピュータに、回転機械の作動状態における振動情報を採取する手順、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手順、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手順、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手順、を実行させ、前記状態評価指数が以下の式により導き出されるプログラムである。
【数９】
また本発明は、回転機械の診断プログラムにおいて、コンピュータに、回転機械の作動状態における振動情報を採取する手順、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手順、前記統合次元の主成分に基く統合次元の状態評価指数を算出する手順、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手順、を実行させ、前記状態評価指数が以下の式により導き出されるＳであるプログラムである。
【数１０】
また本発明は、上記記載の発明であるプログラムにおいて、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断プログラムである。
【００３２】
また本発明は、上記記載の発明であるプログラムにおいて、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する異常原因推定手順、を備えた回転機械の診断プログラムである。
【００３３】
また本発明は、回転機械の作動状態における振動情報を採取する手段と、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手段と、前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手段と、前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手段とを備え、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断装置である。
【数１１】
また本発明は、回転機械の作動状態における振動情報を採取する手段と、前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手段と、前記統合次元の主成分に基く統合次元の状態評価指数を算出する手段と、前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手段とを備え、前記状態評価指数が以下の式により導き出されるＳである回転機械の診断装置である。
【数１２】
また本発明は、上記記載の発明である回転機械の診断装置において、前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含む回転機械の診断装置である。
【００３６】
また本発明は、上記記載の発明である回転機械の診断装置において、正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する手段を更に供えた回転機械の診断装置である。
【００４０】
【発明の実施の形態】
図１は、本発明に係る診断システムの概略の構成を示すブロック図である。 本発明に係る診断システムは、異常判定を行うための諸機能を備えた判定処理装置１、振動センサ２からの電荷信号を電圧信号等に変換する信号変換器３、記録媒体に記録されている振動信号を再生する記録信号再生装置４、判定処理装置１からの情報を表示する表示装置５及び判定処理装置１からの情報を印刷する印刷装置６で構成されている。
【００４１】
そして判定処理装置１は、入出力制御部１０、波形読込部１１、ＦＦＴ処理部１２、パラメータ算出部１３、主成分分析部１４、統合パラメータ値算出部１５、良否判定処理部１６、異常原因分析処理部１７、判定結果処理部１８及びデータ部２０とで構成されている。
【００４２】
入出力制御部１０は、外部装置との信号授受を行い情報交換を行うためのインターフェースである。波形読込部１１は振動センサ２または記録信号再生装置４から振動信号を読み込むとともに必要時にはその信号から加速度信号と速度信号を生成する。ＦＦＴ処理部１２は読み込んだ振動信号の周波数分析を行う。
【００４３】
パラメータ算出部１３は、波形読込部１１とＦＦＴ処理部１２の処理結果に基づいて有次元パラメータと無次元パラメータを算出する。主成分分析部１４は、正常状態にある回転機械の振動信号に基づいて主成分分析を行い固有ベクトルを算出する。統合パラメータ値算出部１５は、有次元パラメータ、無次元パラメータ及び固有ベクトルとから回転機械の異常を判定するための指標である統合パラメータ値を算出する。
【００４４】
良否判定処理部１６は、統合パラメータ値に基づいて回転機械の良否を判定する。異常原因分析処理部１７は、回転機械が異常であると判定された場合にその異常の原因を分析して推定する。判定結果処理部１８は、良否判定結果及び異常原因推定結果を表示装置５あるいは印刷装置６に出力する。
【００４５】
データ部２０は、上述の各処理部が異常判定処理を行うために用いる各種データを記憶する。初期状態データベース２１は、回転機械の正常状態において算出した各種パラメータ値などを記憶する。測定値データベース２２は、回転機械の現在状態において算出した各種パラメータ値などを記憶する。基準値データベース２３は、回転機械の良否を判定するための基準値を記憶する。異常原因分析データベース２４は、異常原因を推定するための基準値データを記憶する。
【００４６】
尚、本発明に係る診断システムでは、振動信号を振動センサ２から直接に、あるいは記録信号再生装置４から読み込んでいるが、この形態に限定されず例えば図示しない通信回線を介して遠隔から振動信号を読み取るように構成しても良い。また波形読込部１１、ＦＦＴ処理部１２を判定処理装置の内部に設けず外部の装置を用いて構成し、その処理結果を入力するように構成しても良い。
【００４７】
次に、このように構成された診断システムの動作について説明する。
図２は、診断システムの概略の動作を示すフロー図である。
【００４８】
まず、波形読込部１１が、機械の初期状態（正常状態）における振動波形を読み込む（Ｓ１）。ここで、波形読込部１１は振動加速度波形ならびに振動速度波形を読み込むように構成しても良く、また振動加速度波形を読み込んだ後、その信号を積分して振動速度波形を算出するように構成しても良い。
【００４９】
次に、ＦＦＴ処理部１２が読み込んだ振動信号の周波数分析を実行し（Ｓ２）、各周波数毎の成分値を算出する。
【００５０】
続いて、パラメータ算出部１３が起動して、それぞれの波形データから以下に示すような複数の有次元パラメータと無次元パラメータを算出する（Ｓ３）。
【００５１】
（１）有次元パラメータ
１）振動速度
ａ.速度ピーク値 ：ＶＥＬ−Ｐ
測定した振動速度波形の振幅値ｘｉの内、｜ｘｉ｜の大きなものから数えた上位5％の｜ｘｉ｜の平均値である。
【００５２】
ｂ．速度ＲＭＳ値 ：ＶＥＬ−Ｒ
【数５】
ｃ．周波数成分値 ：ＶＥＬ−ｆｒ（回転周波数ｆｒ成分）
ここで回転周波数ｆｒは、Ｎ:回転数［rpm］を用いて下式で表される。
【００５４】
【数６】
ｄ．周波数成分値 ：ＶＥＬ−２ｆｒ（周波数２×ｆｒ成分）
ｅ．周波数成分値 ：ＶＥＬ−３ｆｒ（周波数３×ｆｒ成分）
ｆ．周波数成分値 ：ＶＥＬ−４ｆｒ（周波数４×ｆｒ成分）
ｇ．周波数成分値 ：ＶＥＬ−５ｆｒ（周波数５×ｆｒ成分）
ｈ．周波数成分値 ：ＶＥＬ−１／２ｆｒ（周波数１／２×ｆｒ成分）
ｉ．周波数成分値 ：ＶＥＬ−１／３ｆｒ（周波数１／３×ｆｒ成分）
ここで表された有次元パラメータの内、ｃの回転周波数成分は、主にアンバランス異常時に卓越する成分、ｄ〜ｇの回転周波数の高調波成分は、ミスアライメントやガタ・ゆるみ異常時に発生する成分、ｈ〜ｉの分数調波成分は、ガタ・ゆるみ発生時に表れる特徴的な成分である。
【００５６】
２）振動加速度
ａ．加速度ピーク値：ＡＣＣ−Ｐ
測定した振動加速度波形の振幅値ｘｉの内、｜ｘｉ｜の大きなものから数えた上位5％の｜ｘｉ｜の平均値である。
【００５７】
ｂ．加速度ＲＭＳ値：ＡＣＣ−Ｒ
【数７】
ｃ．周波数成分値 ：ＡＣＣ−ｆ０（外輪キズ周波数成分）
【数８】
ｄ．周波数成分値 ：ＡＣＣ−ｆｉ（内輪キズ周波数成分）
【数９】
ｅ．周波数成分値 ：ＡＣＣ−ｆｂ（ボールキズ周波数成分）
【数１０】
ｆ．周波数成分値 ：ＡＣＣ−２ｆｂ（ボールキズ周波数成分）
【数１１】
ｇ．周波数成分値 ：ＡＣＣ−ｆｃ（保持器キズ周波数成分）
【数１２】
ｈ．周波数成分値：ＡＣＣ−ｆｒ（回転周波数成分）
【数１３】
ここで、Ｎ:回転数［rpm］、ｆｒ：軸（内輪）の回転周波数(Hz)、Ｄ：軸受のピッチ円直径(ｍｍ)、d：転動体の直径(ｍｍ)、α：接触角(度)、z：転動体の数である。図３は一般的な軸受の主要諸元を示す図である。
【００６５】
（２）無次元パラメータ
１）振動速度
ａ．波形率Ｓｆ：ＶＥＬ−Ｓｆ
振動波形の正弦波からのずれを表すパラメータ。低周波領域のアンバランスやミスアライメント、脈動波形などの定量化に有効なパラメータ。
【００６６】
【数１４】
ｂ．波高率Ｃｆ：ＶＥＬ−Ｃｆ
振動波形の衝撃性を表すパラメータ。局部異常の検出に有効なパラメータ。
【００６８】
【数１５】
ｃ．衝撃指数Ｉｐ：ＶＥＬ−Ｉｐ
波高率Ｃｆと同様に衝撃性を表すパラメータ。
【００７０】
【数１６】
ｄ．スキューネス（歪み度）β１：ＶＥＬ−β１
振動波形がゼロ点を中心にしていかに非対称となっているかを示すパラメータ。摩耗系の異常が発生すると、振動波形が非対称となり、歪み度が増大する。
【００７２】
【数１７】
ｅ．クートシス（尖り度）β２：ＶＥＬ−β２
振動波形がゼロ点を中心にしていかに尖っているかを示すパラメータ。転がり軸受や歯車装置の異常診断に有効なパラメータ。
【００７４】
【数１８】
２）振動加速度
振動速度と同様に、以下のパラメータを求める。
【００７６】
ａ．波形率Ｓｆ ：ＡＣＣーＳｆ
ｂ．波高率Ｃｆ ：ＡＣＣ−Ｃｆ
ｃ．衝撃指数Ｉｐ ：ＡＣＣ−Ｉｐ
ｄ．スキューネス（歪み度）β１：ＡＣＣ−β１
ｅ．クートシス（尖り度）β２：ＡＣＣ−β２
以上の有次元パラメータと無次元パラメータを算出後、主成分分析部１４が起動して算出された有次元パラメータあるいは無次元パラメータ毎に、主成分分析を行う（Ｓ４）。
【００７７】
主成分算出までの計算手順を以下に説明する。尚、以下の説明では簡便化のために有次元パラメータのみについて説明するが、無次元パラメータについても同様の手順とする。
【００７８】
ある設備の正常状態下で上述の手順により収集されたｍ個の有次元パラメータを要素とする兆候パラメータＹｐ＝（Ｙ１，Ｙ２，Ｙ３ ・・・，Ｙｍ）のｎ組のデータからなる行列Ｙ０を以下の式で定義する。
【００７９】
【数１９】
次に、この行列Ｙ０の列方向、即ち縦方向の行列要素に対して以下の式を用いて演算を行い、ｙｉを算出する。
【００８１】
【数２０】
この演算を列毎に行うことによって、新たな標準化された兆候パラメータを要素とするデータ行列Ｙを求める。
【００８３】
【数２１】
この行列は各列毎に平均値＝０、分散＝１に変換された行列である。
【００８５】
そこで、このデータ行列Ｙから相関行列Ｒを算出すると以下の式となる。
【００８６】
【数２２】
ここに、ｒijは２つの列の相関係数である。すなわち、
【数２３】
次に、相関行列Ｒの固有値λを求める。即ち、以下の式を満たす固有値λを求める。
【００８９】
【数２４】
求められた固有値をλ１,λ２，・・・λｎとする。ただし、λ１＞λ２，・・・＞λｎである。これらのｎ個の固有値に対する固有ベクトルをａｉ（ｉ＝１、・・・ｎ）とする。
【００９１】
【数２５】
そうするとそれぞれの固有ベクトルに対して次の式が成立する。
【００９３】
【数２６】
この式に基づいて、ａｉ１、ａｉ２、・・・ａｉｍを求める。但し、この係数は以下の式を充たす値である。
【００９５】
【数２７】
この手順を繰り返して、夫々の固有値λ１,λ２，・・・λｎに対する固有ベクトルａ１，ａ２，・・・ａｎを求めることができる。そしてこれらのｎ個の固有ベクトルと標準化された兆候パラメータ（ｙ１、ｙ２、・・・ｙｍ）を組み合わせて、主成分Ｚ１，Ｚ２，・・・Ｚｎが、以下の式で表される。
【００９７】
【数２８】
ここで、Ｚ１を第１主成分、Ｚ２を第２主成分、Ｚｎを第ｎ主成分と呼ぶ。
【００９９】
次に、統合パラメータ値算出部１５が起動して、有次元統合パラメータＳｙを算出する（Ｓ５）。
【０１００】
前述の正常状態下で求めた主成分Ｚｉ（母集団）が、正規分布に従うと仮定する。この母集団から独立に取り出されたｎ個の標本で構成される統計量χ２ は、自由度ｎ−１のカイ２乗分布に従う。
【０１０１】
【数２９】
ここで、母集団から取り出したｎ個の標本をＸ１，Ｘ２，・・・Ｘｎとすると、標本分散ｓ２は次の式で表される。
【０１０３】
【数３０】
また、主成分Ｚｉの母分散σ２は固有値λiに等しいことから、
σ２＝ λi
これらの関係を整理すると以下の式となり、標準化されたＸ値の偏差平方和は、自由度ｎー１のカイ自乗分布に従う。
【０１０５】
【数３１】
ここで、式のＸｉを主成分Ｚｉに置き換えると、
【数３２】
次に、正常状態下におけるデータの主成分Ｚｉが１−αの確率で入る領域は、以下の式で表される。
【０１０８】
【数３３】
よって、正常状態の状態確定領域は次の式を満たす範囲となる。
【０１１０】
【数３４】
たとえば、有意水準α＝0.05、自由度φ＝3の場合は、
【数３５】
となり、この正常状態確定領域に入ったときは正常、領域外のときは異常という判定が可能となる。
【０１１３】
従って、正常状態からの変化量を監視するために、下の式で表される状態量を統合パラメータとして定義する。
【０１１４】
【数３６】
このＳ値が大きくなれば異常状態と判定でき、回転機械の良否を判定するための普遍的な状態量を示すものである。この式において、有次元パラメータについて集約した指標値を有次元統合パラメータ値Ｓｙとし、無次元パラメータについて集約した指標値を無次元統合パラメータ値Ｓｍとする。そして、この２つの状態量Ｓｙ、Ｓｍを監視することにより、設備の状態監視を行う。
【０１１６】
ＳｙとＳｍの２つを監視する理由は、異常モードによって、有次元パラメータに大きな変化が現れる場合と、無次元パラメータに大きな変化が現れる場合があるためである。
【０１１７】
以上の手順に従って、初期状態波形データについて有次元パラメータ、無次元パラメータ、主成分値、固有値、固有ベクトルなどが求められるが、それらの値は正常状態における基準データとして初期状態データベース２１に記憶する（Ｓ６）。
【０１１８】
以上説明した波形データ処理は正常状態における基準データを作成するための準備である。したがって、この正常状態処理が行われた時点以降に読み込まれる、回転機械の良否を判定するための振動波形の処理については、この基準データを利用する形態での処理となる。
【０１１９】
正常状態処理が行われた時点以降に読み込まれる判定用振動波形の処理においては、ステップＳ１からＳ３までは前述と同様に実行される。しかし、ステップＳ４の主成分分析においては、標準化された兆候パラメータを要素とするデータ行列Ｙを求めた後は、再度新たな固有ベクトルの算出は行わず、初期状態データベース２１に記憶されている固有ベクトルを抽出して、その固有ベクトルに基づいて主成分Ｚ１、・・・Ｚｎを算出する。そして、この主成分に基づいて前述のステップＳ５と同様の手順で有次元統合パラメータＳｙと無次元統合パラメータＳｍを算出する。このようにして算出された判定用振動波形の処理データおよび処理結果は、測定値データベース２２に記憶される。
【０１２０】
次に、良否判定処理部１６が起動して算出された有次元統合パラメータＳｙと無次元統合パラメータＳｍに基づいて回転機械の良否を判定する（Ｓ７）。良否判定処理部１６は、このＳｙまたはＳｍの値が、初期状態よりもｎ倍になれば注意、ｍ倍になれば異常、ｎ倍未満であれば良と判定する。この判定に際し、初期状態のＳｙまたはＳｍは初期状態データベース２１から取り出され、判定基準であるｎ倍またはｍ倍は基準値データベース２３に格納されている値が用いられる。
【０１２１】
判定結果で、注意または異常と判定された場合は、異常原因分析処理部１７が起動して異常原因を分析して推定する（Ｓ８）。異常原因分析処理部１７は異常原因分析データベース２４から、図４に示す異常原因分析マトリックスを取り出す。
【０１２２】
異常原因分析マトリックスは、縦方向に異常原因を列記し、横方向に有次元及び無次元パラメータを配置した構成である。そして、この異常原因分析マトリックスの行列要素には、異常度指数１（ａｉｊ）と異常度指数２（ｂｉｊ）が記載されている。この異常度指数１と異常度指数２は、該異常が発生した場合にそのパラメータが示す値を統計的に処理した指標であり、過去に発生した異常に実績データに基づいて算出し決定した値である。ここで、異常度指数１と異常度指数２の複数の指数を設けているのは、異常原因によってはいずれか一方のみの指数しか変化しない場合があるためである。尚、その異常度指数が該異常に関係しない場合は、その異常度指数には数値は設定されず、以降の計算には組み込まれない。
【０１２３】
次に、異常原因分析処理部１７は、測定値データベース２２に格納されている判定用振動波形の処理データを用いて、図５に示す異常度指数１（Ａｊ）と異常度指数（Ｂｊ）とを算出する。そして、次の式に基づいて異常原因（ｉ）についてパラメータ（ｊ）毎に異常確率指数Ｐｉｊを算出する。
【０１２４】
Ｐｉｊ ＝ αｊ＊｛（Ａｊ／ａｉｊ）＋（Ｂｊ／ｂｉｊ）｝
但し、ａｉｊ、ｂｉｊともに値が設定されているとき：αｊ＝０．５
ａｉｊ、ｂｉｊいずれか一方に値が設定されているとき：αｊ＝１
ａｉｊ、ｂｉｊともに値が設定されていないとき：αｊ＝０
そして、次の式によって異常原因（ｉ）毎に総合確率Ｔｉを算出する。
【０１２５】
Ｔｉ ＝１／Ｎ＊ΣＰｉｊ
ここで、Ｎ：ａｉｊまたはｂｉｊに値が設定されているパラメータの数
この式において、総合確率Ｔｉが大きな値をとる場合は、あらかじめ異常原因分析マトリックスに設定されてある異常原因と類似している場合であると考えられる。従って、総合確率Ｔｉが所定以上になれば、その異常が発生している可能性が高いと判断する。
【０１２６】
続いて、判定結果処理部１８が起動して、判定結果である良、注意、異常と、異常原因（ｉ）毎の総合確率Ｔｉを編集して表示装置５あるいは印刷装置６に出力する（Ｓ９）。その後、必要なデータを保存して（Ｓ１０）、異常判定処理を終了する。
【０１２７】
次に、本発明の診断システムの適用例を、図面を参照しながら説明する。図６は、振動センサ２が取り付けられた回転機械であるモータ３０の概略図である。
【０１２８】
振動センサ２によって、正常状態における振動速度波形と加速度波形を採取し、採取した２種類の波形から、有次元パラメータと無次元パラメータの値を算出する。
【０１２９】
そして、上記の有次元・無次元パラメータに対し、それぞれ主成分分析を行い、初期波形データにおける各パラメータの値、主成分値、固有値、固有ベクトル等をすべて記憶しておく。
【０１３０】
次に、モータ３０に初期異常を発生させそのデータを計測する。具体的には軸受に軽微な欠陥のある異常データ５データを計測し、初期データ解析時と同様、有次元統合パラメータと無次元統合パラメータを算出する。図７は、正常状態と異常状態における波形データの有次元統合パラメータの値を示す図である。無次元統合パラメータについても同様に算出する。
【０１３１】
そして主成分分析を行い状態量Ｓｙ、Ｓｍ値を算出して加速度ＲＭＳ値（ＡＣＣ−Ｒ）とともに記載した結果を図８に示す。図９はＳｙ値とＡＣＣ−Ｒのデータをグラフ化した図であり、図１０はＳｍ値とＡＣＣ−Ｒのデータをグラフ化した図である。
【０１３２】
Ｎｏ．１１〜１５の異常データに対し、有次元パラメータである加速度ＲＭＳ値(ＡＣＣ−Ｒ)はほとんど変化していなく、異常を検知していない。これに対して、有次元統合パラメータＳｙ値と無次元統合パラメータＳｍ値は５倍以上も大きく変化しており、Ｓｙ値とＳｍ値を用いることで初期異常を感度良く検知していることがわかる。
【０１３３】
尚、有次元統合パラメータＳｙと無次元統合パラメータＳｍ単独で判定せずに、ＳｙとＳｍを変数とする関数を考え、その関数の示す値によって判定するように構成しても良い。
【０１３４】
図１１は、Ｓｙ値とＳｍ値と積Ｓｙ＊Ｓｍのデータをグラフ化した図である。Ｎｏ．１１〜１５の異常データに対し、指標Ｓｙ＊Ｓｍが初期異常の兆候をより顕著に表していることがわかる。このように、Ｓｙ値とＳｍ値を組合わせた指標を用いることで、初期異常の発生を確実に検出することが可能となる。
【０１３５】
また、本実施の形態では、有次元統合パラメータＳｙと無次元統合パラメータＳｍの２つの指標を求めているが、この形態に限定されず有次元パラメータと無次元パラメータとを分離せずにまとめ、これから統合パラメータを算出するように構成しても良い。
【０１３６】
また、有次元パラメータと無次元パラメータの選択を異常モード毎に行い、異常モード毎の有次元統合パラメータと無次元統合パラメータを用いるように構成しても良い。
【０１３７】
更に、異常モードの初期、中期、末期の各状態毎に、有次元パラメータと無次元パラメータを選択し、その状態毎の有次元統合パラメータと無次元統合パラメータを用いて判定するように構成しても良い。
【０１３８】
以上のように、この発明によれば、振動波形から算出した２つの統合パラメータ値を監視する機能により、回転機械の劣化状態を的確に、かつ早期に検知することができる。また、初期データからの波形の大きさと形状から統合パラメータを算出するので、回転機械の型式に依存しない精度の高い異常判定が可能となる。従って従来手法では困難であった特殊回転機械の診断も可能となるなど、信頼性の高い診断が可能となり、産業上有用な効果がもたらされる。
【０１３９】
【発明の効果】
この発明における診断手法を用いて回転機械の状態監視を行うことにより、異常の早期検知が可能となり、精度の高い良否判定を行うことができる。
【図面の簡単な説明】
【図１】 本発明に係る診断システムの概略の構成を示すブロック図。
【図２】 診断システムの概略の動作を示すフロー図。
【図３】 一般的な軸受の主要諸元を示す図。
【図４】 異常原因分析マトリックスを示す図。
【図５】 異常度指数を示す図。
【図６】 振動センサが取り付けられたモータを示す図。
【図７】 正常状態と異常状態における波形データの有次元統合パラメータの値を示す図。
【図８】 Ｓｙ、Ｓｍ値、加速度ＲＭＳ値を示す図。
【図９】 Ｓｙ値とＡＣＣ−Ｒの推移を示す図。
【図１０】 Ｓｍ値とＡＣＣ−Ｒの推移を示す図。
【図１１】 Ｓｙ、Ｓｍ値、Ｓｙ＊Ｓｍ値の推移を示す図。
【符号の説明】
１…判定処理装置、２…振動センサ、１１…波形読込部、１２…ＦＦＴ処理部、１３…パラメータ算出部、１４…主成分分析部、１５…統合パラメータ値算出部、１６…良否判定処理部、１７…異常原因分析処理部、１８…判定結果処理部。
Claims (12)
- 回転機械の作動状態における振動情報を採取し、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、
前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出し、
前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出し、
前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行い、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とする回転機械の診断方法。
- 回転機械の作動状態における振動情報を採取し、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出し、
前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出し、
前記統合次元の主成分に基く統合次元の状態評価指数を算出し、
前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行い、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とする回転機械の診断方法。
- 前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、
前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含むことを特徴とする請求項１又は２に記載の回転機械の診断方法。 - 正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定することを特徴とする請求項１又は２に記載の回転機械の診断方法。
- 回転機械の診断プログラムにおいて、
コンピュータに、
回転機械の作動状態における振動情報を採取する手順、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、
前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手順、
前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手順、
前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手順、
を実行させ、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とするプログラム。
- 回転機械の診断プログラムにおいて、
コンピュータに、
回転機械の作動状態における振動情報を採取する手順、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手順、
前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手順、
前記統合次元の主成分に基く統合次元の状態評価指数を算出する手順、
前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手順、
を実行させ、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とするプログラム。
- 前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、
前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含むことを特徴とする請求項５又は６に記載の回転機械の診断プログラム。 - 正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する異常原因推定手順、
を備えたことを特徴とする請求項５又は６に記載の回転機械の診断プログラム。 - 回転機械の作動状態における振動情報を採取する手段と、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、
前記複数の有次元振動パラメータを標準化した後主成分分析法をもちいて有次元の主成分を抽出するとともに、前記複数の無次元振動パラメータを標準化した後主成分分析法をもちいて無次元の主成分を抽出する手段と、
前記有次元の主成分に基く有次元の状態評価指数と、前記無次元の主成分に基づく無次元の状態評価指数とを算出する手段と、
前記有次元の状態評価指数と無次元の状態評価指数とを組合わせた指数に基づき前記回転機械の良否判定を行う手段と
を備え、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とする回転機械の診断装置。
- 回転機械の作動状態における振動情報を採取する手段と、
前記振動情報に基づき振動を特徴づける複数の有次元振動パラメータと複数の無次元振動パラメータとを算出する手段と、
前記複数の有次元振動パラメータと、前記複数の無次元振動パラメータとをまとめて標準化した統合振動パラメータから主成分分析法をもちいて統合次元の主成分を抽出する手段と、
前記統合次元の主成分に基く統合次元の状態評価指数を算出する手段と、
前記統合次元の状態評価指数に基づき前記回転機械の良否判定を行う手段と
を備え、
前記状態評価指数が以下の式により導き出されるＳであることを特徴とする回転機械の診断装置。
- 前記有次元振動パラメータは、少なくとも振動のピーク値、ＲＭＳ値、回転周波数成分値、回転周波数の高調波成分値、回転周波数の分数調波成分値、及びベアリング傷周波数成分値のいずれかを含み、
前記無次元振動パラメータは、少なくとも波形率、波高率、衝撃指数、スキューネス及びクートシスのいずれかを含むことを特徴とする請求項９又は１０に記載の回転機械の診断装置。 - 正常作動状態の振動パラメータと、異常と判定した状態の振動パラメータとの変化量を算出し、予めこれらの変化量と各異常原因との相関を関連付けた判定基準と比較・判定することにより異常原因を推定する手段を更に供えたことを特徴とする請求項８又は９に記載の回転機械の診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003067071A JP4312477B2 (ja) | 2003-03-12 | 2003-03-12 | 回転機械の診断方法、診断装置及びそのプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003067071A JP4312477B2 (ja) | 2003-03-12 | 2003-03-12 | 回転機械の診断方法、診断装置及びそのプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004279056A JP2004279056A (ja) | 2004-10-07 |
| JP4312477B2 true JP4312477B2 (ja) | 2009-08-12 |
Family
ID=33284789
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003067071A Expired - Fee Related JP4312477B2 (ja) | 2003-03-12 | 2003-03-12 | 回転機械の診断方法、診断装置及びそのプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4312477B2 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4573036B2 (ja) * | 2005-03-16 | 2010-11-04 | オムロン株式会社 | 検査装置および検査方法 |
| JP3780299B1 (ja) * | 2005-06-24 | 2006-05-31 | 独立行政法人科学技術振興機構 | 対象設備の診断方法、コンピュータプログラム、及び、対象設備を診断するための装置 |
| JP2008058191A (ja) * | 2006-08-31 | 2008-03-13 | Jfe Steel Kk | 回転機械の診断方法、そのプログラム、及びその診断装置 |
| CN102798528B (zh) * | 2012-08-07 | 2015-08-12 | 洛阳轴研科技股份有限公司 | 通过能耗系数来评价轴承自身能耗的方法 |
| JP6475589B2 (ja) * | 2015-08-03 | 2019-02-27 | 一般財団法人電力中央研究所 | 人の行動の同定方法、同定装置、及び同定プログラム |
| JP6944285B2 (ja) * | 2017-06-29 | 2021-10-06 | 川崎重工業株式会社 | 回転部材を有する装置の異常原因特定システム |
| JP2023085839A (ja) * | 2021-12-09 | 2023-06-21 | 株式会社日立インダストリアルプロダクツ | 流体機械の異常原因推定装置及びその異常原因推定方法並びに流体機械の異常原因推定システム |
-
2003
- 2003-03-12 JP JP2003067071A patent/JP4312477B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004279056A (ja) | 2004-10-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3321487B2 (ja) | 機器／設備診断方法およびシステム | |
| JP4874406B2 (ja) | 軸受の診断システム | |
| Mechefske et al. | Fault detection and diagnosis in low speed rolling element bearings Part I: The use of parametric spectra | |
| US8720275B2 (en) | Detecting rotor anomalies | |
| JPH0315698B2 (ja) | ||
| JP5740208B2 (ja) | 軸受診断方法及びシステム | |
| JP7198089B2 (ja) | 電力変換装置、回転機システム、及び診断方法 | |
| CN111964909A (zh) | 滚动轴承运行状态检测方法、故障诊断方法及系统 | |
| CN113757093B (zh) | 一种闪蒸汽压缩机组故障诊断方法 | |
| JP6750644B2 (ja) | 異常監視方法および異常監視装置 | |
| CN113418730A (zh) | 一种卷烟机运行状态的在线监测方法 | |
| KR20190081933A (ko) | 제조 설비의 이상 감지 및 진단 방법 | |
| JP4312477B2 (ja) | 回転機械の診断方法、診断装置及びそのプログラム | |
| JP2010175446A (ja) | 状態診断装置 | |
| CN110553789A (zh) | 一种压阻式压力传感器的状态检测方法、装置及制动系统 | |
| JP3333494B2 (ja) | 装置の異常診断方法 | |
| JPWO2004068078A1 (ja) | 状態判定方法と状態予測方法及び装置 | |
| CN112711850A (zh) | 一种基于大数据的机组在线监测方法 | |
| JP2008058191A (ja) | 回転機械の診断方法、そのプログラム、及びその診断装置 | |
| JP2020027563A (ja) | ヘルスモニタリングシステム | |
| CN113418731A (zh) | 一种卷烟机组的在线故障诊断方法 | |
| CN112964470A (zh) | 一种基于轴承健康指标检测滚动轴承早期故障的方法 | |
| JP3103193B2 (ja) | 回転機械の診断装置 | |
| CN112857806B (zh) | 一种基于移动窗口时域特征提取的轴承故障检测方法 | |
| CN113609964A (zh) | 一种电机异常振动预警方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20040702 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060309 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070613 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090203 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090401 |
|
| 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: 20090428 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090513 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120522 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4312477 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120522 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130522 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140522 Year of fee payment: 5 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |