(57)【要約】【課題】窓の時間幅の選択に依存せず基底振動周波数の「ゆらぎ」を正確に測定して概周期振動の周期の安定性の指標とする、振動波形の周期の安定性評価方法を提供する。【解決手段】ある時間Ｔ１について、近似的周期が０．７Ｔ１〜１．４Ｔ１の範囲にある概周期性の振動波形の時系列データに対して、Ｔ１より大きいＴ２について、時刻ｔ〜ｔ＋Ｔ２の「窓」に属する時系列データに第１のＦＦＴ（高速フーリエ変換）を施し、得た周波数スペクトルの基底振動における最大振幅値をＭ（ｔ）とし、Ｔ１より小さいＴ３だけ前記窓を順次シフトして得た時系列データＭ（ｔ）に第２のＦＦＴを施し、得た周波数スペクトルの基底振動の周波数拡散幅ＤＦｃと基本周波数Ｆｃの比である周波数拡散係数ＤＦｃ／Ｆｃを指標とする。【選択図】図１
2008039499.tif 000002
【特許請求の範囲】
【請求項１】
第１の時間Ｔ１について、近似的周期がＴ１の７０％より大きく１４０％より小さい範囲にある概周期性の振動波形の時系列データに対して、
第１の時間Ｔ１より大きい第２の時間Ｔ２について、ある時刻ｔから時刻ｔ＋Ｔ２までの時間（以下、「窓」という）に属する前記時系列データに第１のＦＦＴ（高速フーリエ変換）を施し、前記第１のＦＦＴ演算により得られた周波数スペクトルの基底振動における最大振幅値をＭ（ｔ）とし、
第１の時間Ｔ１より小さい第３の時間Ｔ３について、前記窓を第３の時間Ｔ３だけシフトして同様にＭ（ｔ＋Ｔ３）を得、
第２の時間Ｔ２より大きい第４の時間Ｔ４について、時刻ｔ＝０から時刻ｔ＝Ｔ４＋Ｔ２まで、複数個の順次シフトされた窓に対して得た前記最大振幅値Ｍ（ｔ）の時系列データに第２のＦＦＴを施し、
前記第２のＦＦＴ演算により得られた周波数スペクトルの基底振動において振幅が最大となる周波数を基本周波数Ｆｃとし、
前記基底振動の周波数拡散幅ＤＦｃと前記基本周波数Ｆｃの比である周波数拡散係数ＤＦｃ／Ｆｃを指標とすることを特徴とする振動波形の周期の安定性評価方法。
【請求項２】
前記振動波形の時系列データが、機械装置の物理的振動を観測して得られるものであり、前記周波数拡散係数により前記機械装置の部品の消耗（磨耗又はキズ、または電源の消耗）の程度を評価することを特徴とする請求項１に記載の振動波形の周期の安定性評価方法。
【請求項３】
前記機械装置がドリルを含む工作機械であり、前記機械装置の物理的振動が、前記ドリル先端の被加工物への接触・貫通により発生するものであり、前記機械装置の部品の消耗が前記ドリルを駆動するモータの部品（ブラシ、ベアリングを含む）の消耗またはモータの電源の消耗であることを特徴とする請求項２に記載の振動波形の周期の安定性評価方法。
【請求項４】
前記機械装置がリレーを含む機械装置であり、前記機械装置の物理的振動が、前記リレーのチャタリングにより発生するものであり、前記機械装置の部品の消耗が前記リレーの接点の消耗であることを特徴とする請求項２に記載の振動波形の周期の安定性評価方法。
【請求項５】
前記機械装置が車輪を含む機械装置であり、前記機械装置の物理的振動が、前記車輪の軸受け摩擦により発生するものであり、前記機械装置の部品の消耗が前記車輪の軸受け部品（ベアリングを含む）の消耗であることを特徴とする請求項２に記載の振動波形の周期の安定性評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、機械装置などの振動波形の周期の安定性評価方法に係り、特に、概周期性の振動を発生する機械装置などにおいて、その構成部品の異常が振動周期の安定性に影響している場合に前記異常とその程度を検出することができる振動波形の周期の安定性評価方法に関するものである。
【背景技術】
【０００２】
工作機械、電車、ＨＤＤ、あるいはそれらに使われるモータなどの各種の機械装置を稼動させると、いずれも周期性の振動を発生し、正常動作時にはその周期性が厳密に維持されるが、ベアリング、ブラシなどの消耗（磨耗・キズ）を含む構成部品の異常が発生すると、振動の周期性が失われ、ついには故障（暴走又は停止）に至る。
故障を予防するためには、この異常を初期段階で検知することが肝要であるが、異常の初期で消耗（磨耗・キズ）などが軽微な場合は、振動の周期性が僅かに失われた、概周期性振動となる場合が多い。
【０００３】
そのため、この概周期性振動の判定、即ち周期の安定性の評価に当たっては、各種の方法が用いられている。
通常用いられる方法は、振動波形をある時間幅（以下、「窓」という）についてＦＦＴ（高速フーリエ変換）により周波数領域に変換し、その基底振動周波数を算定し、これを別の窓について繰り返し、得られた基底振動周波数の「ゆらぎ」を評価するものである。
しかしながらこの方法では、窓の時間幅が概周期に比べて、大きいと「ゆらぎ」が平均化されてしまい、逆に小さいと基底振動周波数の分解能が劣化し「ゆらぎ」が歪曲されてしまい、いずれにしても適切な窓の時間幅を選ぶのが困難であった。
さらに「ゆらぎ」の程度が経時的に変化する場合が多く、その場合は少なくとも窓の最適時間幅を経時的に変動させる必要があり、そのような窓の時間幅の選択はますます困難であった。
従って、「ゆらぎ」の程度を定量的に把握することが困難であり、個別の機械装置ごと、個別の異常ごとに、熟練者のノウハウ（例えば、振動音の微妙な音色の差異）に依存するしかなかった。
【０００４】
例えば特許文献１には、振動波形のエンベロープを検出してこれを周波数解析する技術が開示されているが、上記の困難は必ずしも回避できるものではない。
【特許文献１】特開２００６−１１３００３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上記のような振動波形の周期の安定性評価方法における諸問題を解決するためになされた本発明の目的は、窓の時間幅をどう選ぶかに実質的に依存しないようにして基底振動周波数の「ゆらぎ」を正確に測定して概周期振動の周期の安定性の指標とする、振動波形の周期の安定性評価方法を提供することにある。
【課題を解決するための手段】
【０００６】
上記の課題を解決するためになされた本発明による振動波形の周期の安定性評価方法は請求項１に示すように、第１の時間Ｔ１について、近似的周期がＴ１の７０％より大きく１４０％より小さい範囲にある概周期性の振動波形の時系列データに対して、第１の時間Ｔ１より大きい第２の時間Ｔ２について、ある時刻ｔから時刻ｔ＋Ｔ２までの時間（以下、「窓」という）に属する前記時系列データに第１のＦＦＴ（高速フーリエ変換）を施し、前記第１のＦＦＴ演算により得られた周波数スペクトルの基底振動における最大振幅値をＭ（ｔ）とし、第１の時間Ｔ１より小さい第３の時間Ｔ３について、前記窓を第３の時間Ｔ３だけシフトして同様にＭ（ｔ＋Ｔ３）を得、第２の時間Ｔ２より大きい第４の時間Ｔ４について、時刻ｔ＝０から時刻ｔ＝Ｔ４＋Ｔ２まで、複数個の順次シフトされた窓に対して得た前記最大振幅値Ｍ（ｔ）の時系列データに第２のＦＦＴを施し、前記第２のＦＦＴ演算により得られた周波数スペクトルの基底振動において振幅が最大となる周波数を基本周波数Ｆｃとし、前記基底振動の周波数拡散幅ＤＦｃと前記基本周波数Ｆｃの比である周波数拡散係数ＤＦｃ／Ｆｃを指標とすることを特徴とする。
【０００７】
また請求項２に示すように、前記振動波形の時系列データが、機械装置の物理的振動を観測して得られるものであり、前記周波数拡散係数により前記機械装置の部品の消耗（磨耗又はキズ、または電源の消耗）の程度を評価することを特徴とする。
【０００８】
また請求項３に示すように、前記機械装置がドリルを含む工作機械であり、前記機械装置の物理的振動が、前記ドリル先端の被加工物への接触・貫通により発生するものであり、前記機械装置の部品の消耗が前記ドリルを駆動するモータの部品（ブラシ、ベアリングを含む）の消耗またはモータの電源の消耗であることを特徴とする。
【０００９】
また請求項４に示すように、前記機械装置がリレーを含む機械装置であり、前記機械装置の物理的振動が、前記リレーのチャタリングにより発生するものであり、前記機械装置の部品の消耗が前記リレーの接点の消耗であることを特徴とする。
【００１０】
また請求項５に示すように、前記機械装置が車輪を含む機械装置であり、前記機械装置の物理的振動が、前記車輪の軸受け摩擦により発生するものであり、前記機械装置の部品の消耗が前記車輪の軸受け部品（ベアリングを含む）の消耗であることを特徴とする。
【発明の効果】
【００１１】
本発明の請求項１に係る振動波形の周期の安定性評価方法によれば、窓の時間幅Ｔ２を概周期Ｔ１に比べて「ゆらぎ」が平均化されない程度に小さくとって第１のＦＦＴ演算を行い、その際、分解能が劣化している基底振動周波数を使わない。
その代わりに、基底振動周波数に対応する最大振幅値を用い、窓を順次スタガーして得られる最大振幅値の時系列データを、周波数分解能が十分出せる時間Ｔ３にわたって得て、これに対して第２のＦＦＴ演算を行い、その基底振動周波数及びその周波数拡散係数を高い分解能をもって得ることができるので、精度のよい定量的な、振動波形の周期の安定性評価方法を提供することができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明に係る実施の形態と効果を、図１ないし図１０に従って具体的に説明する。
【実施例１】
【００１３】
対象とする振動波形は、概周期性、即ち、理想的な周期振動に対して、その周期に１周期ごとに少しずつ変動（「ゆらぎ」）を与えたものであり、「ゆらぎ」の範囲は、ここでは、ある概周期Ｔ１に対して、例えば７０％から１４０％の範囲にある。
図１は第１の実施例に係る、振動波形の周期の安定性評価方法による、安定性指標としての周波数拡散係数を得るためのフローチャートである。
全体のフローはフロー（１）とフロー（２）に分けられる。
フロー（１）では、最初にステップＳ１で、対象とする機械装置の振動を圧力センサにより検出し、次にステップＳ２で、圧力センサからの微弱な電気信号をアンプにより増幅する。
【００１４】
次にスレッド１のステップＳ１０で、時刻ｓの初期値を設定した後、ステップＳ１１で、アンプの出力信号をＡ／Ｄ変換し、サンプルデータＤ（ｓ）を得る。
ここで時刻ｓはｓ＝０で始まる、刻みΔＴでデジタイズされたサンプリング時刻である。
次にステップＳ１２で、サンプルデータＤ（ｓ）をリングメモリに書き込む。
次にステップＳ１３で、時刻ｓをΔＴだけシフトし、ステップＳ１３で時刻ｓが所定の時刻Ｔ４＋Ｔ２を越えない間はステップＳ１１に戻り、Ｓ１１、Ｓ１２を繰り返す。
ここで、後述するように、Ｔ４は最後の「窓」の開始時刻（最初の「窓」の開始時刻を０とする）、Ｔ２は「窓」の時間幅である。
この結果、サンプルデータはＤ（０）、Ｄ（ΔＴ）、Ｄ（２ΔＴ）・・・の順にシフトしながらリングメモリに格納され、Ｄ（Ｔ４＋Ｔ２）を格納するまで続けられる。
【００１５】
一方、フロー（２）では、スレッド２のステップＳ２０、Ｓ２１で、各々時刻ｔと窓上の時刻ｕの初期値を設定した後、ステップＳ２２で、サンプルデータＤ（ｔ）をリングメモリから読み出し、ステップ２３で、窓上の時刻ｕをΔＴだけシフトして、ステップＳ２４で時刻ｕが所定の時刻Ｔ２を越えない間はステップＳ２２に戻り、Ｓ２２、Ｓ２３を繰り返す。
この結果、時刻ｔに始まる１つの窓［ｔ、ｔ＋Ｔ２］に関するサンプルデータＤ（ｔ）〜Ｄ（ｔ＋Ｔ２）がリングメモリから順次読み出される。
窓の時間幅Ｔ２は、対象とする振動波形の概周期Ｔ１より大きく、好ましくは概周期Ｔ１の１．４倍より大きい。
【００１６】
最後のＤ（ｔ＋Ｔ２）が読み出された時点で、リングメモリ中の先頭のＤ（ｔ）はシフト先がなく、消去されてもよい。
次にステップＳ３１で、窓［ｔ、ｔ＋Ｔ２］に関する時系列データＤ（ｔ）〜Ｄ（ｔ＋Ｔ２）に対して第１のＦＦＴ演算が行われ、ステップＳ３２で、第１のＦＦＴ演算の結果である周波数特性から基底振動の最大振幅値Ｍ（ｔ）を得る。
次にステップ３３で、時刻ｔを窓のシフト量Ｔ３だけシフトして、ステップＳ３４で、時刻ｔがＴ４を越えない間はステップＳ２１に戻り、次の窓［ｔ＋Ｔ３、ｔ＋Ｔ３＋Ｔ２］に関する時系列データＤ（ｔ＋Ｔ３）〜Ｄ（ｔ＋Ｔ３＋Ｔ２）がリングメモリから読み出され、その窓に関する基底振動の最大振幅値Ｍ（ｔ＋Ｔ３）を得、これを繰り返す。
ここで、窓のシフト量Ｔ３は、求める周波数拡散係数の精度により決まるが、相次ぐ窓の一部が互いに重畳するように、少なくとも窓の時間幅Ｔ２より小さい。
【００１７】
スレッド１では、リングメモリへの格納（書き込み）だけが行われ、スレッド２では、リングメモリからの読み出しだけが行われる。
リングメモリとして、直列書き込み機能と並列読み出し機能が非同期で実行可能なＦＩＦＯ（Ｆｉｒｓｔ−Ｉｎ−Ｆｉｒｓｔ−Ｏｕｔ）メモリを使うと、さらに高速演算処理が可能になる。
【００１８】
ステップＳ４１では、上記スレッド２で得られた、データＭ（０）、Ｍ（Ｔ３）、Ｍ（２・Ｔ３）、・・・Ｍ（Ｔ４）に対して第２のＦＦＴ演算が行われる。
時間Ｔ４は、周波数の分解能を確保するために十分大きく、少なくとも「窓」の時間幅Ｔ２よりも大きくとられる。
【００１９】
最後にステップＳ４２では、上記第２のＦＦＴ演算から基底振動の基本周波数Ｆｃと、その周波数拡散幅ＤＦｃを得て、両者の比、ＤＦｃ／Ｆｃをもって求める概周期の安定性指標とする。
【００２０】
もとの概周期性の振動波形の周波数の「ゆらぎ」を平均化してしまわず反映するためには、個々の窓の時間幅Ｔ２はあまり大きくとれない。
従って、ステップＳ３１で得られる周波数特性の周波数分解能は、必ずしも高くない。
しかしながら、本実施例では周波数特性のうち、そのような分解能の低い周波数値を使わず、基底振動の最大振幅値Ｍ（ｔ）だけを使い、しかも十分長い時間Ｔ４にわたって、Ｍ（ｔ）の時系列データＭ（０）、Ｍ（Ｔ３）、Ｍ（２・Ｔ３）、・・・Ｍ（Ｔ４）を得ることができる。
従って、このような時系列データから得られる周波数特性の周波数分解能は十分に高くできるので、得られた概周期の安定性指標ＤＦｃ／Ｆｃは精度が高い。
【００２１】
本発明のさらに具体的な実施の形態を、ケースＡ、Ｂ、Ｃ、Ｄ、Ｅに分けて、各々以下の実施例２〜６により説明する。
いずれも、対象とする機械装置はドリルを備えた工作機械であり、ドリルが一定の時間間隔で被加工物である鉄板に孔を開けている場合に生じる工作機械の振動を想定する。
工作機械がドリルの駆動系にあるモータを含めて正常に稼動している場合は、振動波形の基底振動の基本周波数は上記一定の時間間隔の逆数に相当し、「ゆらぎ」がない、即ち、ＦＦＴによる周波数特性において、基底振動の基本周波数及びその倍調波の振幅だけが顕著に大きく（波形の「山がシャープ」）、基本周波数の近傍の周波数拡散幅が実質的にゼロになる。
しかしながら、工作機械になんらかの異常がある場合、特にモータに異常がある場合は、基底振動周波数に「ゆらぎ」を生じ（波形の「山がくずれる」）、周波数拡散幅と基本周波数の比である周波数拡散係数が実質的にゼロでなくなる。
（以下実施例２〜６では、周波数拡散幅を「基本周波数の両側において振幅が基底振動における最大振幅の４分の１になる２つの周波数の差」として定義するが、周波数拡散幅の定義はこれに限定されるものではない。）
本実施例２〜６では、この「ゆらぎ」を正確にかつ定量的に捉えることができることが分かる。
【実施例２】
【００２２】
本実施例はケースＡの場合であり、図２、図３を参照すると、図２は正常な場合、図３はモータのベアリングの消耗に相当する異常が生じている場合である。
図２において、図２（ａ）はもととなる振動波形を示し、約２００ｍｓの周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に「ゆらぎ」が殆どない。
即ち、概周期Ｔ１＝２００ｍｓ、であり、ここでは、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
【００２３】
図２（ｂ）は、図２（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍ（ｔ）を順次プロットした結果である。
もとの振動波形の周期性を反映して、本図のＭの波形も周期性がよいことが視認できる。
事実、図２（ｂ）をＦＦＴ演算して得られる図２（ｃ）は、楕円内に示すように、シャープな基底周波数特性を示し、周波数拡散係数は約０．２である。
（図２（ｃ）において、基底振動の基本周波数Ｆｃ＝６．８００（相対値）、基本周波数での最大振幅＝２０目盛（相対値）、４分の１値幅ＤＦｃ（５目盛での周波数拡散幅）＝１．３６（相対値）、従って周波数拡散係数ＤＦｃ／Ｆｃ＝０．２）
【００２４】
一方図３において、図３（ａ）は異常が生じている場合の振動波形を示し、約１７０ｍｓの概周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に１５０〜１８０ｍｓ程度の「ゆらぎ」がある。
即ち、概周期Ｔ１＝１７０ｍｓ、であり、ここでも、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
しかしながら、「ゆらぎ」の程度を視認によって定量的に把握するのは困難である。
【００２５】
図３（ｂ）は、図３（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍ（ｔ）を順次プロットした結果である。
もとの振動波形の概周期の「ゆらぎ」を反映して、本図のＭ（ｔ）の波形も周期性が崩れていることが視認できる。
事実、図３（ｂ）をＦＦＴ演算して得られる図３（ｃ）は、楕円内に示すように、山が崩れた基底周波数特性を示し、基底周波数での最大振幅値が１４．５目盛と低下しているだけでなく、周波数拡散係数は約０．７となり、周波数の「ゆらぎ」の程度が定量的に把握でき、周波数拡散係数が概周期の安定性指標として好適に使えることが分かる。
【実施例３】
【００２６】
本実施例はケースＢ、即ち、上記実施例２の場合よりもモータの回転数を上げた場合であり、図４、図５を参照すると、図４は正常な場合、図５はモータに過負荷をかけた場合である。
ただし本実施例では、被加工物である鉄板の配置が上記実施例１の場合と異なり、ドリルが鉄板に当たる衝撃の周期も異なる。
図４において、図４（ａ）はもととなる振動波形を示し、実質的に約１００ｍｓの周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に「ゆらぎ」が殆どない。（図４（ａ）（ｂ）では横軸が５０ｍｓ／ＤＩＶであることに留意）
ただし鉄板の配置が異なる結果として振動は実施例１の場合よりも速やかに減衰している。
本実施例でも、上記実施例１と同じく概周期Ｔ１、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、２００ｍｓ、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
【００２７】
図４（ｂ）は、図４（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の周期性を反映して、本図のＭの波形も周期性がよいことが視認できる。
事実、図４（ｂ）をＦＦＴ演算して得られる図４（ｃ）は、楕円内に示すように、シャープな基底周波数特性を示し、周波数拡散係数は約０．２である。
【００２８】
一方図５において、図５（ａ）は異常が生じている場合の振動波形を示し、約１００ｍｓの概周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に９０〜１１０ｍｓ程度の「ゆらぎ」がある。
即ち、概周期Ｔ１＝１００ｍｓ、であり、ここでも、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
しかしながら、「ゆらぎ」の程度を視認によって定量的に把握するのは困難である。
【００２９】
図５（ｂ）は、図５（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の概周期の「ゆらぎ」を反映して、本図のＭ（ｔ）の波形も周期性が崩れていることが視認できる。
事実、図３（ｂ）をＦＦＴ演算して得られる図３（ｃ）は、楕円内に示すように、山が崩れた基底周波数特性を示し、周波数拡散係数は約０．５である。
ただし本実施例では、基底周波数での最大振幅値が２０目盛であって、正常な場合と同等であり、「ゆらぎ」の程度は上記実施例１の場合ほど大きくないことが分かる。
【実施例４】
【００３０】
本実施例はケースＣの場合であり、図６、図７を参照すると、図６は正常な場合、図７は上記実施例２と同じくモータのベアリングの消耗に相当する異常が生じている場合である。
ただし本実施例では、上記実施例２と異なり、鉄板が工作機械本体に対してより強固に固定されている。
図６において、図６（ａ）はもととなる振動波形を示し、実質的に約２００ｍｓの周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に「ゆらぎ」がない。
ただし鉄板が強固に固定されている結果として振動は実施例１の場合よりも速やかに減衰している。
本実施例でも、上記実施例２と同じく概周期Ｔ１、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、２００ｍｓ、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
【００３１】
図６（ｂ）は、図６（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の周期性を反映して、本図のＭの波形も周期性がよいことが視認できる。
事実、図６（ｂ）をＦＦＴ演算して得られる図６（ｃ）は、楕円内に示すように、（最大振幅値が１８目盛と、上記実施例１の場合よりは若干低いがシャープな基底周波数特性を示し、周波数拡散係数は約０．３である。
【００３２】
一方図７において、図７（ａ）は異常が生じている場合の振動波形を示し、約２００ｍｓの概周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に１６０〜２４０ｍｓ程度の「ゆらぎ」がある。（図７（ａ）（ｂ）では横軸が２００ｍｓ／ＤＩＶであることに留意）
即ち、概周期Ｔ１＝２００ｍｓ、であり、ここでも、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
しかしながら、「ゆらぎ」の程度を視認によって定量的に把握するのは困難である。
【００３３】
図７（ｂ）は、図７（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の概周期の「ゆらぎ」を反映して、本図のＭ（ｔ）の波形の周期性が崩れていることが視認できる。
事実、図７（ｂ）をＦＦＴ演算して得られる図７（ｃ）は、楕円内に示すように、山が崩れた基底周波数特性を示し、最大振幅値は１３目盛、周波数拡散係数は約０．６である。
【実施例５】
【００３４】
本実施例はケースＤ、即ち、上記実施例３において観測時間Ｔ４を増加した場合であり、図８、図９を参照すると、図８は正常な場合、図９はモータの駆動電圧が変動している場合である。
図８において、図８（ａ）はもととなる振動波形を示し、上記実施例３の場合と同じく、実質的に約１００ｍｓの周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に「ゆらぎ」が殆どない。
本実施例でも、上記実施例１と同じく概周期Ｔ１、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、２００ｍｓ、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
【００３５】
図８（ｂ）は、図８（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の周期性を反映して、本図のＭの波形も周期性がよいことが視認できる。
事実、図８（ｂ）をＦＦＴ演算して得られる図８（ｃ）は、楕円内に示すように、最大振幅値は１３．３目盛となり、かなり低下しているが、シャープな基底周波数特性を示し、周波数拡散係数は約０．２である。
【００３６】
一方図９において、図９（ａ）は異常が生じている場合の振動波形を示し、約１００ｍｓの概周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっても周波数に８０〜１４０ｍｓ程度の「ゆらぎ」がある。
即ち、概周期Ｔ１＝１００ｍｓ、であり、ここでも、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
しかしながら、「ゆらぎ」の程度を視認によって定量的に把握するのは困難である。
【００３７】
図９（ｂ）は、図５（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍを順次プロットした結果である。
もとの振動波形の概周期の「ゆらぎ」を反映して、本図のＭ（ｔ）の波形も周期性が崩れていることが視認できる。
事実、図３（ｂ）をＦＦＴ演算して得られる図３（ｃ）は、楕円内に示すように、山が崩れた基底周波数特性を示し、最大振幅値が８目盛と低下しているだけでなく、周波数拡散係数は約０．８である。
【実施例６】
【００３８】
本実施例はケースＥの場合であり、図１０を参照すると、上記実施例２において、同一条件でさらに長時間動作させ、モータの電源が消耗した場合である。
図１０において、図１０（ａ）はもととなる振動波形を示し、約７０ｍｓの周期でドリルが鉄板に当たるたびにその衝撃による減衰振動が発生しており、視認によっては周波数に「ゆらぎ」が大きくあるようには見えない。
即ち、概周期Ｔ１＝７０ｍｓ、であり、ここでは、窓の時間幅Ｔ２、窓のシフト量Ｔ３を各々、５００ｍｓ、２０ｍｓにとっている。（Ｔ３＜Ｔ１＜Ｔ２）
【００３９】
図１０（ｂ）は、図１０（ａ）の波形を窓ごとにＦＦＴ演算し、得られた周波数特性（図示せず）から基底周波数の最大振幅値Ｍ（ｔ）を順次プロットした結果である。
図１０（ｂ）では、もとの振動波形では視認し難かった概周期の「ゆらぎ」を反映して、Ｍ（ｔ）の波形も周期性が崩れていることが視認できる。
事実、図３（ｂ）をＦＦＴ演算して得られる図３（ｃ）は、楕円内に示すように、山が崩れた基底周波数特性を示し、基底周波数での最大振幅値が１０目盛と低下しているだけでなく、周波数拡散係数は約１．０となり、基底周波数の「ゆらぎ」の程度が大変大きいことが定量的に把握でき、周波数拡散係数が概周期の安定性指標として好適に使えることが分かる。
【実施例７】
【００４０】
以上実施例１において、一般的な概周期振動に対して、本発明による振動波形の周期の安定性評価方法を適用する手順の一例を示し、実施例２〜６においては、対象とする機械装置がドリルを備えた工作機械であり、ドリルが一定の時間間隔で被加工物である鉄板に孔を開けている場合に生じる工作機械の振動から、特にモータのベアリングの消耗、モータの過負荷、又はモータの電源の消耗などの異常の存在を予知できることを示した。
しかしながら、本発明は、上記の具体例に限定されず、同様の手順に従い、リレーのチャタリングに起因する振動からリレーの接点の消耗を予知する際に適用することができる。
また、同様の手順に従い、車輪の軸受け摩擦に起因する振動から車輪の軸受け部品の磨耗を予知する際に適用することができる。
【産業上の利用可能性】
【００４１】
本発明による振動波形の周期の安定性評価方法を適用すると、色々な機械装置において観測される概周期性の振動波形から、振動波形を視認しただけでは必ずしも判別が不可能な周波数の「ゆらぎ」を定量的に検出し、その機械装置の部品の消耗を早い段階で予知できるので、機械装置の安定性評価、保守点検を効率的かつ的確に行うことができ、産業上利するところが大である。
【図面の簡単な説明】
【００４２】
【図１】実施例１に係る、振動波形の周期の安定性評価方法のフローチャートである。
【図２】実施例２に係る正常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図３】実施例２に係る異常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図４】実施例３に係る正常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図５】実施例３に係る異常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図６】実施例４に係る正常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図７】実施例４に係る異常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図８】実施例５に係る正常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図９】実施例５に係る異常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【図１０】実施例６に係る異常な場合の（ａ）原振動波形、（ｂ）最大振幅値波形、（ｃ）ＦＦＴ演算結果である。
【符号の説明】
【００４３】
Ｔ１    第１の時間（概周期）
Ｔ２    第２の時間（「窓」の時間幅）
Ｔ３    第３の時間（「窓」のシフト量）
Ｔ４    第４の時間（総観測時間）
ΔＴ    （サンプリング時間の）刻み
Ｄ（ｓ）、Ｄ（ｔ）、Ｄ（ｕ）    時刻ｓ、ｔ、ｕにおけるサンプルデータ
Ｍ（ｔ）    基底振動における最大振幅値
Ｆｃ    基底振動における基本周波数
ＤＦｃ    基底振動における基本周波数
ＤＦｃ／Ｆｃ    周波数拡散係数
【図１】
2008039499.tif 000003
【図２】
2008039499.tif 000004
【図３】
2008039499.tif 000005
【図４】
2008039499.tif 000006
【図５】
2008039499.tif 000007
【図６】
2008039499.tif 000008
【図７】
2008039499.tif 000009
【図８】
2008039499.tif 000010
【図９】
2008039499.tif 000011
【図１０】
2008039499.tif 000012