(57)【要約】【課題】この発明の実施の形態は、イベント間の時間間隔の確率的な決定のための方法を開示するものであり、イベントは周期的に信号を妨害する。【解決手段】本方法は、イベントの発生確率をジッタされる信号の値に基づいて、時間の関数として、決定するものであり、イベント間の時間間隔を決定するのに適する１組の可能な時間間隔の確率を生成する、イベント間の時間間隔に対する１組の可能な時間間隔の一致の確率を、イベントの発生の確率に基づいて決定する。【選択図】図２
2011013214.tif 000002
【特許請求の範囲】
【請求項１】
信号を周期的に妨害するイベント間の時間間隔を確率的に決定するための方法であって、その方法の工程を実行するためのプロセッサを備え、
前記イベントの発生確率をジッタされる前記信号の値に基づいて、時間の関数として、決定する工程と、
前記イベント間の時間間隔を決定するのに適する１組の可能な時間間隔の確率を生成する、前記イベント間の時間間隔に対する１組の可能な時間間隔の一致の確率を、前記イベントの発生確率に基づいて決定する工程と、
を含む方法。
【請求項２】
前記信号は、前記イベントによりベアリングによって生成された振動信号であり、また前記時間間隔は前記ベアリングの診断に適している、請求項１の方法。
【請求項３】
前記信号は前記ジッタを引き起こす白色ノイズによって達成される、請求項１の方法。
【請求項４】
前記信号を照合フィルタによって濾波することをさらに含む、請求項１の方法。
【請求項５】
前記イベントの発生確率の関数を濾波することをさらに含む、請求項１の方法。
【請求項６】
前記一致の確率を決定することは、前記イベント間の可能な時間間隔の確率と前記信号を妨害する前記イベントの発生確率との間の関数関係に基づく、請求項１の方法。
【請求項７】
前記イベント間の時間間隔を、前記可能な時間間隔の確率に基づいて前記可能な時間間隔の組から選択することをさらに含む、請求項１の方法。
【請求項８】
前記可能な時間間隔の組を前記ベアリングの幾何学的配列に基づいて決定することをさらに含む、請求項２の方法。
【請求項９】
前記イベントの発生確率を決定することは、
対応する時間で前記信号の１組の極大値の値を決定すること、
前記極大値の値に基づいて前記信号の妨害確率特性を決定すること、
前記妨害確率特性に基づいて前記イベントの発生の事後確率関数を決定すること、
をさらに含む、請求項１の方法。
【請求項１０】
前記極大値の値の組を決定することは半波高全幅値（ＦＷＨＭ）方法に基づく、請求項９の方法。
【請求項１１】
前記妨害確率特性は、前記イベントの発生の事前確率、イベントの発生の条件付き確率、およびイベントの不発生の条件付き確率を含む、請求項９の方法。
【請求項１２】
前記妨害確率特性を決定することは、前記イベントの発生の条件付き確率を、前記極大値の値の組への対数正規分布の最大尤度適合度として決定することをさらに含む、請求項１１の方法。
【請求項１３】
前記イベントの発生の条件付き確率の期待値最大化（ＥＭ）を使用して、前記イベントの不発生の条件付き確率を決定することをさらに含む、請求項１２の方法。
【請求項１４】
前記イベントの発生の条件付き確率の期待値最大化（ＥＭ）を使用して、前記イベントの発生の事前確率を決定することをさらに含む、請求項１２の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、一般的に、信号を妨害する周期的なイベント間の時間間隔を決定するための方法に関し、特に信号を妨害するイベント間の時間間隔を実時間で確率的に推定するための方法に関する。
【背景技術】
【０００２】
ベアリングは、滑り（平）軸受け、ボール（玉）軸受け、ローラ（ころ）軸受け、ニードル軸受け、テーパ（円錐）軸受け、スヘリカル（球面）軸受け、あるいはスラスト軸受けであろうが、遍在している。ベアリングは、世界を動かしており、モータ（電動機）、ジェネレータ（発電機）、ホイール（車輪）、タービン、ディスクドライブおよびジェットエンジンなどのすべてのタイプ（型）の設備で見つかる。ベアリングのデザイン（設計）は非常に簡単であるが、ベアリングの故障により大惨事に至る場合がある。したがって、ベアリングの故障をリアルタイムに検知することが望まれる。
【０００３】
図１は、ボールベアリング１１０を示し、それは内レース１１５、外レース１１７、およびそれらのレース間の軸受エレメント、たとえばボール１１９、を含む。ベアリングは、回転するにつれて、故障により、ベアリングの角加速度に比例する周波数で共振する。共振は、信号１２２として経時的にサンプリングすることができる。
【０００４】
典型的に、故障は、単一のポイントにあり、ベアリングによって生成された信号１２２を妨害して信号中にピーク１２０を生じさせる周期的なイベントを引き起こす。理想的には、ピークは周期的であり、時間間隔１２５によって分離される。図示された信号は時間（ｔ）の関数としての加速度（ｇ）である。
【０００５】
ベアリングおよび相互イベント間隔の幾何学的配列に基づいて、異なる故障を識別することが可能である。しかしながら、ノイズは故障関連の周期的なイベントを隠すことがある。たとえば、ベアリング・スリップによって引き起こされた機械騒音すなわちボールとレースとの間の欠陥による接点の変動は、振幅１３０および時間間隔１３１が変化するように、信号をジッタ（小刻みに変動）させる。
【０００６】
信号を妨害するイベント間の時間間隔を測定するための１つの従来手法は、覆い隠されている信号を自動的に関連させて、固有振動数に対応する遅れでピークを検知することである。ここで、信号１４０は、遅れの関数として自己相関ｒ（ｔ）として示される。しかしながら、自己相関は少数の擬似的な最大の振幅ピークによって支配されるので、その方法はジッタ（小刻みに変動）させられた信号の場合には失敗する。
【０００７】
別の方法は、同時に信号を復調し、かつ固有振動数でエネルギを検知するために、バイスペクトル分析を使用する。しかしながら、そのバイスペクトルに基づく方法は、信号の周期性に関する仮定を行う。
【０００８】
別の方法は、信号に由来した特徴ベクトルに機械学習技法を適用することにより、信号の周期性に変動がある状態で故障を検出する。しかしながら、その方法は訓練データを必要とし、また訓練データに表われない状況へ上手く一般化できない場合があり、その方法をリアルタイムの故障検出のために役立たなくしている。
【発明の概要】
【発明が解決しようとする課題】
【０００９】
したがって、トレーニングを必要とせずに、ジッタ（小刻みに変動）させられた信号を妨害する周期的イベント間の時間間隔を決定することが望まれる。
【００１０】
この発明の目的は、トレーニングを必要とせずに、ジッタ（小刻みに変動）させられた信号を妨害する周期的イベント間の時間間隔を決定することである。
【課題を解決するための手段】
【００１１】
この発明の実施の形態は、イベント間の可能な時間間隔の事後確率が、時間の関数として、信号を妨害するイベントの発生確率の関数であるという認識に基づく。さらに、その確率は、たとえば信号の加速度の振幅などの値に基づいて決定することができる。
【００１２】
この発明の実施の形態は、イベント間の時間間隔の確率的な決定のための方法を開示するものであり、イベントは周期的に信号を妨害する。本方法は、イベントの発生確率をジッタされる信号の値に基づいて、時間の関数として、決定するものであり、イベント間の時間間隔を決定するのに適する１組の可能な時間間隔の確率を生成する、イベント間の時間間隔に対する１組の可能な時間間隔の一致の確率を、イベントの発生の確率に基づいて決定する。
【発明の効果】
【００１３】
この発明の実施の形態は、イベント間の時間間隔を決定するための方法を記述するものであり、イベントは周期的に信号を妨害する。イベント間の時間間隔は、単一点の欠陥を有するベアリングの故障関連の妨害（外乱）の特性を検出するために使用される。これらの実施の形態は、故障関連の振動妨害（外乱）の確率モデルを使用し、リアルタイムに実行することができる。
【図面の簡単な説明】
【００１４】
【図１】この発明の実施の形態１によって解決されるベアリング故障診断問題の概略図である。
【図２】この発明の実施の形態１による、信号を妨害するイベント間の時間間隔を決定する方法のブロックダイアグラムである。
【図３Ａ】イベントによって部分的に妨害される信号のグラフである。
【図３Ｂ】イベントによって部分的に妨害される信号のグラフである。
【図４】イベントの発生確率のヒストグラムである。
【図５】可能な時間間隔の確率のヒストグラムである。
【図６】この発明の実施の形態１の実際的なセッティングの例である。
【発明を実施するための形態】
【００１５】
実施の形態１．
図２は、イベント２０５の時間間隔２９０を決定する方法２００を示し、イベントは周期的に信号ｘ（ｔ）２１０を妨害する。ノイズおよび他のランダム（無作為）な影響により、信号２１０はジッタ（小刻みに変動）させられる。ジッタは、信号の連続するパルスの周波数、信号の振幅、および信号の位相を含む。本方法のステップ（工程）はプロセッサ２０１によって実行される。プロセッサ２０１は、当該技術において知られているように、メモリ、Ｉ／Ｏインタフェース、および信号プロセッサを含む。
【００１６】
図３Ａは信号２１０の一例３１０を示す。信号はジッタさせられた期間（以下、ジッタ期間と称す）３２０を有する。ジッタ期間は、信号の特性の最大あるいは最小イフェクト（効力）の２つのイベント間の時間間隔である。さらに、信号の振幅３２０もジッタさせられる、すなわち、経時的に変化する。
【００１７】
１つの実施の形態では、信号２１０は、イベント２０５により故障するベアリング１１０によって生成される。故障が局地化される場合、イベントはイベント間の時間間隔２９０で周期的である。時間間隔２９０の長さは故障の診断２８０のための重要な要因である。しかしながら、他の実施の形態では、信号２１０は、ジッタを蒙りやすい任意の周期的信号、たとえば任意の電磁気信号、である。信号のジッタは、たとえば、指数関数的に分布した白色加法的バックグランドノイズによって引き起こされる場合がある。
【００１８】
この発明の実施の形態は、イベント間の可能な時間間隔の事後確率ｐｉ（τ）が、時間の関数として、信号を妨害するイベントの発生確率ｐｄ（ｔ）の関数であるという認識に基づく。さらに、その確率ｐｄ（ｔ）は、たとえば信号の加速度の振幅などの値に基づいて決定することができる。
【００１９】
したがって、本方法２００は、信号の値に基づいてイベントの発生確率ｐｄ（ｔ）２６５を決定２６０する。１つの実施の形態では、信号からノイズを除去するために信号を濾波２５０して濾波信号ｍ（ｔ）２５５を生成する。この実施の形態では、濾波工程は、たとえば、信号のＮ個の最大ピークを含む窓で信号を平均することにより作成された照合フィルタを使用する。この記述のために、信号ｘ（ｔ）およびｍ（ｔ）は区別なく同義で用いられている。図３Ｂは、フィルタリング（濾波）後の信号の例を示す。この濾波は、低振幅イベントを明らかにするためにバックグランドノイズを排除する。
【００２０】
工程２６０は点別変換ｐｄ（ｔ）＝ｇ（ｍ（ｔ））を適用する。ここで、ｇ（ｍ（ｔ））は、次式（１）によって信号ｍ（ｔ）の値を観察した後に、時間ｔでのイベントの発生の事後確率である。
【００２１】
【数１】
2011013214.tif 000003
【００２２】
ここで、Ｐｄはイベントの発生の事前確率であり、ｐ（ｍ｜ｄ＝１）はイベントの発生の条件付き確率であり、すなわち、イベントが時間ｔで発生するとして、信号のｍ（ｔ）の値を観察することであり、また、ｐ（ｍ｜ｄ＝０）はイベントの不発生の条件付き確率であり、すなわち、イベントが時間ｔで生じないとして、信号のｍ（ｔ）の値を観察することである。
【００２３】
ここに定義されたように、イベントの発生の事前確率、イベントの発生の条件付き確率、およびイベントの不発生の条件付き確率は、イベントの確率特性である。
【００２４】
いくつかの実施の形態は、信号を妨害するイベントが、イベントの（発生）時間に近接する時間ｔの範囲に対して信号ｍ（ｔ）の値を増大させるという観察結果告を使用する。これらの実施の形態では、確率ｐｄ（ｔ）の比較的大きな極大値のみを保存することによりイベントの発生確率をフィルタし、また、個々の大きな極大値に近い確率ｐｄ（ｔ）を零に設定する。したがって、信号の中間から大きなピークの時に、妨害（外乱）確率ｐｄ（ｔ）は１に近く、また明白な妨害（外乱）の無い時には、確率ｐｄ（ｔ）は零（０）近い。信号の値の中の小さなピークは、不確実な照合フィルタ出力を生成し、図４で示されるように、ｐｄ（ｔ）に対して中間値に帰着する。
【００２５】
確率ｐｄ（ｔ）２６５に基づいて、本方法２００は、１組の可能な時間間隔２３０の確率を生成する時間間隔２９０に対する可能な時間間隔２２０の一致の確率を決定２７０する。その組２３０はイベント間の時間間隔を決定するために使用される。たとえば、１つの実施の形態では、その組２３０から、最も高い確率の可能な時間間隔を選択する。１つの実施の形態では、可能な時間間隔２２０はベアリング１１０の幾何学的配列に基づいて選択される。
【００２６】
その可能な時間間隔の確率の組は、必ずしも、時間間隔２９０に対する実際の一致の確率を含まない。１つの実施の形態では、確率は、正規化され、すなわち、互いに関連して、最も有望な可能な時間間隔を決定することを可能にする。
【００２７】
これらの実施の形態では、次の式（２）によりその可能な時間間隔の確率の組を決定する。
【００２８】
【数２】
2011013214.tif 000004
【００２９】
ここで、Ｐｉ（τ）は可能な時間間隔に関する事前確率であり、それは一定であり、かつその可能な時間間隔の値である。
【００３０】
図５は、可能な時間間隔の確率ｐｉ（τ）の正規化されたダイアグラムを示す。たとえば、イベント間の時間間隔が１０ミリセカンド（ｍｓ）である確率は０であり、また、時間間隔が８ミリセカンド（ｍｓ）であるの確率は５である。したがって、１つの実施の形態１では、最も高い確率値を有する可能な時間間隔を選択することにより、１組の確率２３０に基づいて期間２９０を決定する。
【００３１】
イベントの発生確率
典型的には、信号を妨害する殆どのイベントは信号ｍ（ｔ）において顕著なピークとなる。したがって、恐らくイベントの発生に対応する時間｛ｔ０、…、（ｔＭ｝）で信号の１組の極大値｛ｍ（ｔ０）、…、ｍ（ｔＭ）｝を決定する。
【００３２】
１つの実施の形態では、半波高全幅値（ＦＷＨＭ）方法を使用して、その組の中の極大値の数を制限する。ＦＷＨＭとは、従属変数がその極大値の半分と等しくなる独立変数の２つの極値間の差により与えられる、関数の範囲の表現である。
【００３３】
最初に、信号のｍ（ｔ）の大域的な（全体）最大を決定して、ｍ（ｔ０）選択する。そして、追加のイベントが、時間ｔ０の近くの領域、ｍ（ｔ）＞０．５ｍ（ｔ０）、すなわち、ｔ０の付近の半波高全幅値の領域  に存在しないと仮定する。
【００３４】
次に、時間ｔ１で、その領域の外部におけるｍ（ｔ）の最大値を求め、また極大値の値の組の中へｍ（ｔ１）を選択する。続いて、その信号を「半波高全幅値」領域へ分割し、極大値の値の組および対応する時間を決定する。
【００３５】
次いで、極大値の値の組に基づいてイベントの確率特性を決定する。１つの実施の形態では、極大値の値の組に対する対数正規分布の最大尤度適合度として確率ｐ（ｍ｜ｄ＝１）を決定する。他の実施の形態では、正数のデータをモデル化するための異なるパラメータ形式、たとえば正規分布、を使用する。
【００３６】
条件付き確率ｐ（ｍ｜ｄ＝１）を与えられることで、この分配を、複数の対数正規分布の２成分混合物の１つの成分として固定し、また、混合分布モデルが単に極大値の値に対してだけではなくすべての時間に対して完全な組（フルセット）の値ｍ（ｔ）をモデル化するように、別の（第２の）成分のパラメータおよび第２の成分の事前確率を決定するために、期待値最大化（ＥＭ）手続きを使用する。
【００３７】
期待値最大化（ＥＭ）を行なった後では、１つの混合成分をｐ（ｍ｜ｄ＝１）の最大尤度適合度に固定したので、第２の成分は非妨害（妨害のない）時の信号を表わすであろう。したがって、第２の成分はｐ（ｍ｜ｄ＝０）に対応する。ｐ（ｍ｜ｄ＝１）に対応する成分に対する事前確率は、確率Ｐｄの推定値を生み出す。
【００３８】
この発明の実施の形態は、イベント間の時間間隔を決定するための方法を記述するものであり、イベントは周期的に信号を妨害する。イベント間の時間間隔は、単一点の欠陥を有するベアリングの故障関連の妨害（外乱）の特性を検出するために使用される。これらの実施の形態は、故障関連の振動妨害（外乱）の確率モデルを使用し、リアルタイムに実行することができる。
【００３９】
図６は、この発明の実施の形態を使用した機械故障診断方法のためのセッティングの非制限（制限されない）例を示す。モータ６１０の振動信号２１０は加速度計６２０によって検知され、本方法２００を実行するプロセッサ２０１に向けられる６３０。本方法２００の実行の結果は解析器６４０に導かれる。さらに、信号２１０は、更なる調査のために振動監視モジュール６５０に転送される。
【００４０】
この発明は好ましい実施の形態を例として記述されたが、この発明の趣旨および範囲内で様々な他の改変および変更を行ってもよいことが理解されるべきである。したがって、この発明の真実の趣旨および範囲内に入るような、すべての変更例および変形例をカバーすることが、添付の特許請求の範囲の目的である。
【図１】
2011013214.tif 000005
【図２】
2011013214.tif 000006
【図３Ａ】
2011013214.tif 000007
【図３Ｂ】
2011013214.tif 000008
【図４】
2011013214.tif 000009
【図５】
2011013214.tif 000010
【図６】
2011013214.tif 000011
【外国語明細書】
2011013214.tif 000012
2011013214.tif 000013
2011013214.tif 000014
2011013214.tif 000015
2011013214.tif 000016
2011013214.tif 000017
2011013214.tif 000018
2011013214.tif 000019
2011013214.tif 000020
2011013214.tif 000021
2011013214.tif 000022
2011013214.tif 000023
2011013214.tif 000024
2011013214.tif 000025
2011013214.tif 000026
2011013214.tif 000027
2011013214.tif 000028
2011013214.tif 000029
2011013214.tif 000030
2011013214.tif 000031