(57)【特許請求の範囲】
【請求項１】
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得する画像データ取得部と、
前記動作画像に基づいて、前記回転体の振動状態を解析する演算部と、を備え、
前記演算部は、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する、
振動解析システム。
【請求項２】
前記演算部は、
前記動作画像と前記非共振画像との類似度に基づいて、前記動作画像の前記解析位相を推定し、前記位相ロックイン画像を生成する、
請求項１に記載の振動解析システム。
【請求項３】
前記演算部は、
前記動作画像に撮影されている、前記回転体の位相を明示する情報に基づいて、前記動作画像の前記解析位相を推定し、前記位相ロックイン画像を生成する、
請求項１に記載の振動解析システム。
【請求項４】
前記回転体を撮影するカメラを備え、
前記カメラは、フレームタイミングを制御可能な可変フレームレートカメラであり、
前記画像データ取得部は、
前記カメラを制御して、前記回転体の位相が前記解析位相となるタイミングで、前記回転体を撮影し、前記位相ロックイン画像を生成する、
請求項１に記載の振動解析システム。
【請求項５】
前記画像データ取得部は、
前記回転体の回転角情報に基づいて、前記回転体の位相を推定する、
請求項４に記載の振動解析システム。
【請求項６】
前記画像データ取得部は、
動作中の前記回転体の画像から、前記回転体の位相を表す特徴量を抽出し、
前記特徴量に基づいて前記回転体の位相を推定する、
請求項４に記載の振動解析システム。
【請求項７】
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得し、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する、
振動解析方法。
【請求項８】
コンピュータを、
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得する画像データ取得部、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する演算部、
として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、振動解析システム、振動解析方法及びプログラムに関する。
【背景技術】
【０００２】
エンドミル、ドリル等の切削工具、駆動力を伝達する各種の動力軸等の回転体では、回転動作によって、強制振動、自励振動等に起因する共振現象が発生し、振動振幅が増大する場合がある。強制振動は、外力の変化の周波数と一致する振動であり、自励振動は、外力の変化の周波数とは無関係で、回転体の固有振動数と一致する振動である。共振現象が生じると、加工精度の低下、機器の損傷等を生じるおそれがあるため、回転体の振動状態を解析することが重要となる。
【０００３】
一般的に、回転体の振動状態の解析は、回転体近傍に設置された加速度センサを用いて行われる（例えば、特許文献１）。加速度センサを用いた計測では、回転体に直接加速度センサを取り付けることが難しいので、チャック、軸受け等に加速度センサを設置して振動計測が行われる。しかしながら、このような振動計測方法では、回転体を直接計測していないので、計測の感度は低いものとなる。
【０００４】
また、加速度センサを用いない振動解析方法として、画像解析による振動解析方法が開発されている（例えば、特許文献２）。画像解析による振動解析方法では、振動解析の対象物を直接解析できるので、解析の感度を高めることが可能となる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１０−２２１１６１号公報
【特許文献２】特開２０１７−２０７３３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
回転体に対して画像解析による振動解析を行う場合、回転体の外観形状が、回転軸に対してして不均一であると、回転に伴う外観の変化が振動成分として検出される。例えば、エンドミルは、らせん状の切れ刃を有するので、回転角度によって、撮影されるエンドミルの外観形状が変化する。これにより、回転体の回転周波数に係る振動成分が検出されるので、計測の目的である共振周波数付近の振動の検出が困難となる。
【０００７】
本発明は、上述の事情に鑑みてなされたものであり、回転体の回転に伴う振動成分と共振成分とを分離して、共振成分を精度よく検出できる振動解析システム、振動解析方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、この発明の第１の観点に係る振動解析システムは、
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得する画像データ取得部と、
前記動作画像に基づいて、前記回転体の振動状態を解析する演算部と、を備え、
前記演算部は、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する。
【００１０】
また、前記演算部は、
前記動作画像と前記非共振画像との類似度に基づいて、前記動作画像の前記解析位相を推定し、前記位相ロックイン画像を生成する、
こととしてもよい。
【００１１】
また、前記演算部は、
前記動作画像に撮影されている、前記回転体の位相を明示する情報に基づいて、前記動作画像の前記解析位相を推定し、前記位相ロックイン画像を生成する、
こととしてもよい。
【００１２】
また、前記回転体を撮影するカメラを備え、
前記カメラは、フレームタイミングを制御可能な可変フレームレートカメラであり、
前記画像データ取得部は、
前記カメラを制御して、前記回転体の位相が前記解析位相となるタイミングで、前記回転体を撮影し、前記位相ロックイン画像を生成する、
こととしてもよい。
【００１３】
また、前記画像データ取得部は、
前記回転体の回転角情報に基づいて、前記回転体の位相を推定する、
こととしてもよい。
【００１４】
また、前記画像データ取得部は、
動作中の前記回転体の画像から、前記回転体の位相を表す特徴量を抽出し、
前記特徴量に基づいて前記回転体の位相を推定する、
こととしてもよい。
【００１５】
また、本発明の第２の観点に係る振動解析方法では、
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得し、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する。
【００１６】
また、本発明の第３の観点に係るプログラムは、
コンピュータを、
解析対象である回転体の動作状態を表す動作画像と、前記回転体の非共振状態を表す非共振画像とを取得する画像データ取得部、
前記動作画像のうち、前記回転体の位相が前記回転体の回転周期の１／Ｎ（Ｎは自然数）となる解析位相である画像からなる時系列の位相ロックイン画像と前記非共振画像との差分画像に基づいて、前記解析位相毎の振動変位を算出し、
前記解析位相毎の振動変位を合成した合成振動変位に基づいて、前記回転体の振動状態を解析する演算部、
として機能させる。
【発明の効果】
【００１７】
本発明の振動解析システム、振動解析方法及びプログラムによれば、回転体の位相毎の時系列画像に基づいて回転体の振動を解析するので、回転体の回転に伴う振動成分と共振成分とを分離して、共振成分を精度よく検出することが可能である。
【図面の簡単な説明】
【００１８】
【図１】本発明の実施の形態１に係る振動解析システムの機能ブロック図である。
【図２】実施の形態１に係る振動解析の流れを示すフローチャートである。
【図３】回転体及び振動解析システムの構成を示す概略図である。
【図４】回転体の回転周期と位相との関係を示す概念図である。
【図５】動作画像と位相ロックイン画像との関係を示す概念図である。
【図６】回転成分と共振成分のスペクトルを示す図であり、（Ａ）は、動作画像に基づいて検出される振動のスペクトル、（Ｂ）は、位相ロックイン画像に基づいて検出される振動のスペクトルである。
【図７】回転成分と共振成分とを分離するためのマージンを示す概念図である。
【図８】解析位相毎の振動変位を示すグラフであり、（Ａ）は、補正前の振動変位を示すグラフ、（Ｂ）は、基準変位の補正後の振動変位を示すグラフである。
【図９】実施の形態２に係る振動解析の流れを示すフローチャートである。
【図１０】実施の形態３に係る振動解析の流れを示すフローチャートである。
【発明を実施するための形態】
【００１９】
以下、図を参照しつつ、本発明の実施の形態に係る振動解析システム１について説明する。
【００２０】
（実施の形態１）
図１のブロック図に示すように、本実施の形態に係る振動解析システム１は、カメラ１１、制御ユニット２０を備える。
【００２１】
カメラ１１は、振動の解析対象である回転体５０を撮影するカメラである。回転体５０は、エンドミル、ドリル等の切削工具、駆動力を伝達する各種の動力軸等、回転軸を中心として回転するものであり、特に限定されない。カメラ１１は、予め定められた所定のフレームレートＦで回転体５０を撮影し、各時刻の画像データを出力する。
【００２２】
制御ユニット２０は、例えばコンピュータ装置であり、図１に示すように、制御部２１、記憶部２２、表示部２３、入力部２４を備える。
【００２３】
制御部２１は、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、水晶発振器等から構成されており、振動解析システム１の動作を制御する。また、制御部２１は、カメラ１１で撮影された回転体５０の画像データに基づいて、回転体５０の振動状態を解析する。制御部２１は、制御部２１のＲＯＭ、記憶部２２等に記憶されている各種動作プログラム及びデータをＲＡＭに読み込んでＣＰＵ及びＧＰＵを動作させることにより、図１に示される制御部２１の各機能を実現させる。これにより、制御部２１は、画像データ取得部２１１、演算部２１２として動作する。
【００２４】
画像データ取得部２１１は、カメラ１１を制御して、回転体５０を所定のフレームレートＦで撮影する。また、画像データ取得部２１１は、撮影された回転体５０を表す時系列の画像データをカメラ１１から取得する。
【００２５】
演算部２１２は、画像データ取得部２１１で取得された画像データに基づいて、回転体５０の振動状態を解析する。より具体的には、演算部２１２は、デジタル画像相関法等の画像処理により、回転体５０の動作によって生じる共振周波数付近の振動の大きさ、変動等を検出する。
【００２６】
記憶部２２は、ハードディスク、フラッシュメモリ等の不揮発性メモリであり、時系列の画像データから回転体５０の振動状態を解析する演算アルゴリズム、演算部２１２で算出された回転体５０の振動状態を示すデータ等を記憶する。
【００２７】
表示部２３は、コンピュータ装置である制御ユニット２０に備えられた表示用デバイスであり、例えば液晶パネルである。表示部２３は、カメラ１１で撮影された回転体５０の画像、演算部２１２で算出された回転体５０の振動状態を示すグラフ等を表示する。
【００２８】
入力部２４は、解析を行う回転体５０の解析位相等のパラメータの設定、解析の開始、終了指示等を入力するための入力デバイスである。入力部２４は、制御ユニット２０に備えられたキーボード、タッチパネル、マウス等である。
【００２９】
続いて、図２のフローチャートを参照しつつ、振動解析システム１を用いた振動解析方法について説明する。本実施の形態では、図３に示すように、回転体５０であるエンドミルを撮影して、撮影された時系列の画像から、回転体５０の共振周波数付近の振動を検出する場合を例として説明する。
【００３０】
本実施の形態に係る振動解析処理では、まず、回転体５０の振動を解析するための解析位相φＮ（Ｎは１〜ｎの自然数）、すなわち回転体５０の回転角度を設定する（ステップＳ１１）。より具体的には、図４に示すように、回転体５０が１回転する間の、１又は２以上の解析位相φＮを設定し、設定された解析位相φＮ毎に振動を解析する。これにより、回転体５０の外観形状が、位相によって異なることに起因して、振動として検出される回転成分の影響を低減することができる。
【００３１】
解析位相φＮの数であるｎは、フレームレートＦ（ｆｐｓ）と回転体５０の回転周波数ｆＲ（Ｈｚ）とを用いて、例えば、Ｎ＝Ｆ／ｆＲとすることができる。
【００３２】
続いて、画像データ取得部２１１は、回転体５０の非共振画像０Ｉを取得してＬＵＴ（Lookup Table）を作成する（ステップＳ１２）。ＬＵＴに登録される非共振画像０Ｉは、回転体５０の非共振状態、例えばエンドミルである回転体５０が被削体に接触しておらず、無負荷状態で回転することにより、共振を生じない状態の画像である。画像データ取得部２１１は、カメラ１１を制御して、解析位相φＮに対応する非共振画像０Ｉを撮影する。また、画像データ取得部２１１は、撮影した非共振画像０Ｉをカメラ１１から受信してＬＵＴを作成し、記憶部２２に記憶させる。
【００３３】
続いて、画像データ取得部２１１は、回転体５０の動作画像Ｉ（ｘ，ｙ，ｔ）（以下、動作画像Ｉともいう。）を取得する（ステップＳ１３）。より具体的には、画像データ取得部２１１は、カメラ１１を制御して、予め定められた所定のフレームレートＦで、回転体５０を撮影する。動作画像Ｉは、回転体５０の動作状態の画像であり、本実施の形態では、エンドミルである回転体５０が被削材を切削している状態の画像である。すなわち、動作画像Ｉは、負荷状態で動作することによって生じる、共振成分等、回転成分以外の振動成分を含む時系列画像である。
【００３４】
演算部２１２は、画像データ取得部２１１で取得された動作画像Ｉについて、位相を推定する。本実施の形態に係る回転体５０は、一定の回転周波数ｆＲで回転しているので、位相の推定は、時刻を基準に行うことができる。すなわち、図４に示すように、回転体５０の回転周期ＴＲ＝１／ｆＲを用いて、位相φを表すことができる。例えば、時刻ｔの動作画像Ｉ（ｘ，ｙ，ｔ）の位相をφ＝０とすると、時刻ｔ＋ＴＲ／４の動作画像Ｉ（ｘ，ｙ，ｔ＋ＴＲ／４）の位相はφ＝π／２となる。
【００３５】
演算部２１２は、図５に示すように、推定された位相に対応する解析位相φＮ毎に動作画像Ｉを割り振って、時系列の位相ロックイン画像ＩφＮ（ｘ，ｙ）（以下、位相ロックイン画像ＩφＮともいう。）を生成する（ステップＳ１４）。また、演算部２１２は、生成された位相ロックイン画像ＩφＮに基づいて、振動検出を行う。振動検出方法は、特に限定されず、デジタル画像相関法等の公知の方法を用いることができる。本実施の形態では、解析位相φＮ毎に回転体５０の振動状態を解析するので、回転体５０の回転偏心成分がキャンセルされ、位相による回転体５０の外観の変化の影響を低減することができる。
【００３６】
ここで、初期データである動作画像ＩのフレームレートＦが、共振周波数ｆＮに対して十分高く設定されている場合であっても、間引きデータとなる位相ロックイン画像ＩφＮのフレームレートｆＲは、共振周波数ｆＮより低くなるおそれがある。この場合、位相ロックイン画像ＩφＮでは、共振成分の振動がアンダーサンプリングされる。
【００３７】
十分高いフレームレートＦで撮影された動作画像Ｉから検出される振動のスペクトルの例を図６（Ａ）に示す。また、図６（Ａ）に示される回転成分と共振成分の振動がアンダーサンプリングされた場合のスペクトルを図６（Ｂ）に示す。図６（Ｂ）に示すように、アンダーサンプリングされた画像の観測周波数領域は０〜ｆＲ／２であり、見かけの共振周波数ｆ’Ｎは、０≦ｆ’Ｎ≦ｆＲ／２の範囲となる。また、回転成分の振動は、ＤＣ成分であるので、図６（Ｂ）に示すように、見かけの回転周波数とその高調波成分はアンダーサンプリングにより周波数０付近の成分となる。
【００３８】
アンダーサンプリングされた画像に基づいて検出される振動から、共振成分と回転成分とを分離するためには、図７に示すように、見かけの共振周波数と見かけの回転周波数とが離れており、共振周波数ｆＮと回転周波数ｆＲの倍数とが一致しないことが求められる。したがって、分離条件は、マージンΔｆＲを用いて、以下の式で表される。
ｋｆＲ＋ΔｆＲ≦ｆＮ≦（ｋ＋１）ｆＲ−ΔｆＲ
【００３９】
共振周波数ｆＮが回転周波数ｆＲの整数倍とほぼ一致する場合、回転成分のピークがなだらかな場合等では、マージンΔｆＲを小さくできないので、回転成分と共振成分との分離は難しい。
【００４０】
したがって、アンダーサンプリングによる不具合を低減するため、解析位相φＮ毎に算出された変位の時間変化を合成して、サンプリング周波数を高める。より詳細には、各解析位相φＮにおいて算出された変位の時間変化ＶφＮ（ｔ）（図８（Ａ））について、基準変位をゼロに補正する（ステップＳ１５）。本実施の形態では、演算部２１２は、ＬＵＴを参照して、位相ロックイン画像ＩφＮと位相ロックイン画像ＩφＮに対応する解析位相φＮの非共振画像０Ｉとの差分画像を生成する。そして、演算部２１２は、生成した差分画像に基づいて、解析位相φＮ毎の振動変位Ｖ’φＮ（ｔ）を算出する。これにより、図８（Ｂ）に示すように、位相毎の補正後の振動変位Ｖ’φＮ（ｔ）は、基準変位をゼロとする波形となる。
【００４１】
補正後の位相毎の振動変位Ｖ’φＮ（ｔ）は、以下の式で表される。
【数１】
0006842218.tif 000002
ただし、０Ｖ’φｎは位相毎の共振のない状態の変位である。
【００４２】
演算部２１２は、ステップＳ１５で算出された振動変位Ｖ’φＮ（ｔ）を合成して、合成振動変位Ｖ’（ｔ）を生成する（ステップＳ１６）。合成振動変位Ｖ’（ｔ）は、以下の式で表される。
【数２】
0006842218.tif 000003
【００４３】
上述のように、補正後の振動変位Ｖ’φＮ（ｔ）を合成することにより、解析位相φＮ毎の外観形状の違いによって生じる回転成分の影響を低減することができる。また、合成により、フレームレートｆＲをＮ倍した合成振動変位Ｖ’（ｔ）に基づいて振動を解析することができる。したがって、共振周波数ｆＮの観測周波数領域が０〜ｆＲ／２から、０〜ＮｆＲ／２へと拡張され、アンダーサンプリングによる不具合を低減することができる。
【００４４】
制御部２１は、合成振動変位Ｖ’（ｔ）に基づいて、回転体５０の振動状態を解析する（ステップＳ１７）。例えば、算出された合成振動変位Ｖ’（ｔ）を高速フーリエ変換して、表示部２３に表示させる。また、所定の窓時間毎に振動状態を解析して、表示させることにより、回転体５０の振動状態の時間変化を把握することができる。
【００４５】
以上説明したように、本実施の形態に係る振動解析システム１では、回転体５０の位相毎の時系列画像に基づいて振動成分を検出するので、回転体の回転に伴う振動成分と共振成分とを分離して、共振成分の振動を精度よく検出することができる。
【００４６】
また、回転体５０の位相毎に非共振画像０Ｉと動作画像Ｉとの差分画像から変位を算出して、位相毎の変位を合成するので、位相毎の回転成分の影響を低減するとともに、アンダーサンプリングによる不具合を低減し、高周波の共振成分振動を検出することが可能となる。
【００４７】
（実施の形態２）
実施の形態１では、時刻を基準として動作画像Ｉ（ｘ，ｙ，ｔ）の解析位相φＮを推定することとしたが、撮影された回転体５０の外観形状に基づいて解析位相φＮを推定することもできる。本実施の形態では、動作画像Ｉから回転体５０の解析位相φＮを推定して、位相ロックイン画像ＩφＮを生成する場合について説明する。
【００４８】
本実施の形態に係る振動解析システム１では、演算部２１２が、動作画像Ｉについて回転体５０の解析位相φＮを推定する方法が実施の形態１と異なる。その他、カメラ１１、制御ユニット２０の構成等については、実施の形態１と同様であるので同じ符号を付して、詳細な説明は省略する。
【００４９】
以下、図９のフローチャートを参照しつつ、本実施の形態に係る振動解析方法について具体的に説明する。ステップＳ３１の解析位相φＮの設定、ステップＳ３２の非共振画像０Ｉの取得は、実施の形態１に係るステップＳ１１、ステップＳ１２と同様である。本実施の形態では、取得された非共振画像０Ｉを、予め、ステップＳ３１で設定した解析位相φＮと関連付けて、ＬＵＴを作成する。非共振画像０Ｉと解析位相φＮとの関連付けは、例えば、非共振画像０Ｉの撮影時刻に基づいて行われる。この場合、非共振画像０Ｉは、切削加工を行っていない無負荷状態で撮影されるので、切削時の摩擦等による位相のずれは生じない。したがって、動作画像Ｉと比較して、少ない誤差で非共振画像０Ｉと解析位相φＮとを関連付けることができる。
【００５０】
続いて、振動解析システム１は、実施の形態１に係るステップＳ１３と同様に、回転体５０の動作画像Ｉを撮影する（ステップＳ３３）。画像データ取得部２１１は、撮影された動作画像Ｉのデータをカメラ１１から取得する。
【００５１】
演算部２１２は、画像データ取得部２１１で取得された時系列画像データである動作画像Ｉについて、解析位相φＮを推定する（ステップＳ３４）。本実施の形態に係る解析位相φＮの推定は、動作画像Ｉに撮影されている回転体５０の外観形状に基づいて行われる。より具体的には、演算部２１２は、時系列の画像データである動作画像Ｉの各時刻の画像について、ＬＵＴに保存されている非共振画像０Ｉ（ＬＵＴ画像）を参照し、最も類似度の高いＬＵＴ画像を選択する。そして、演算部２１２は、選択されたＬＵＴ画像に関連付けられた解析位相φＮを、動作画像Ｉの解析位相φＮと推定する。
【００５２】
また、演算部２１２は、時刻ｔの動作画像Ｉを、推定した解析位相φＮの位相ロックイン画像ＩφＮを構成する画像データとして、位相ロックイン画像ＩφＮに組み込む。これにより、演算部２１２は、推定された解析位相φＮ毎に位相ロックイン画像ＩφＮを生成する（ステップＳ３５）。
【００５３】
また、演算部２１２は、生成された位相ロックイン画像ＩφＮに基づいて、振動検出を行う。具体的には、演算部２１２は、実施の形態１と同様に、各解析位相φＮにおいて算出された変位の時間変化Ｖφ１（ｔ）について、基準変位をゼロとする補正を行う（ステップＳ３６）。すなわち、演算部２１２は、ＬＵＴを参照して、位相ロックイン画像ＩφＮと位相ロックイン画像ＩφＮに対応する解析位相φＮの非共振画像０Ｉとの差分画像を生成する。そして、差分画像に基づいて振動変位を算出し、基準変位をゼロとする補正後の解析位相毎の振動変位Ｖ’φＮ（ｔ）を生成する。
【００５４】
演算部２１２は、補正後の位相毎の振動変位Ｖ’φＮ（ｔ）を合成して、合成振動変位Ｖ’（ｔ）を生成する（ステップＳ３７）。制御部２１は、実施の形態１と同様に、合成振動変位Ｖ’（ｔ）に基づいて、回転体５０の振動状態を解析する（ステップＳ３８）。例えば、算出された合成振動変位Ｖ’（ｔ）を高速フーリエ変換して、表示部２３に表示させる。また、所定の窓時間毎に振動状態を解析して、表示させることにより、回転体５０の振動状態の時間変化を把握することができる。
【００５５】
以上説明したように、本実施の形態に係る振動解析システム１では、回転体５０の外観形状に基づいて、動作画像Ｉの解析位相φＮを推定するので、回転体５０の動作によって生じる時刻と位相とのずれの影響を低減することができる。
【００５６】
本実施の形態では、回転体５０の外観形状に基づく類似度によって、動作画像Ｉの解析位相φＮを推定することとしたが、これに限られない。例えば、エンドミルの切れ刃の一部等、回転体５０の特徴量に基づいて作成された特徴量マップによって解析位相φＮを推定することとしてもよい。適切な特徴量を抽出することにより、位相推定のための演算部２１２の演算量を低減することができる。
【００５７】
また、回転体５０の一部に回転体５０の位相を明示する外観上の特徴を付し、これに基づいて動作画像Ｉの解析位相φＮを推定することとしてもよい。例えば、エンドミルのチャック部にオプティカルエンコーダで用いられる位相検出パターンを付すこととしてもよい。この場合、例えば１６ビットのパターンを用いることにより、１°以下の位相間隔で振動解析を行うことができる。
【００５８】
（実施の形態３）
上記各実施の形態では、所定のフレームレートＦで撮影された動作画像Ｉについて解析位相φＮを推定して振動解析を行うこととしたが、所定の解析位相φＮで回転体５０を撮影して動作画像Ｉを生成することもできる。本実施の形態では、撮影のタイミングであるフレームタイミングを制御可能な可変フレームレートカメラを用いて、回転体５０を所定の解析位相φＮで撮影する場合について説明する。
【００５９】
本実施の形態に係る振動解析システム１は、カメラ１１が可変フレームレートカメラであり、カメラ１１のフレームタイミングが制御部２１によって制御される点で、実施の形態１と異なる。その他、制御ユニット２０の構成等については、実施の形態１と同様であるので同じ符号を付して、詳細な説明は省略する。
【００６０】
以下、図１０のフローチャートを参照しつつ、本実施の形態に係る振動解析方法について具体的に説明する。本実施の形態に係る振動解析処理では、まず、予め定められた基準位相φ０の非共振画像０Ｉを取得する（ステップＳ５１）。本実施の形態では、回転体５０の保持部、駆動部等に、回転体５０の回転角を計測する回転角計測手段５１が配置されている。回転角計測手段５１は、例えばロータリエンコーダである。画像データ取得部２１１は、回転角計測手段５１から取得された回転体５０の回転角情報である位相データに基づいてカメラ１１を制御し、基準位相φ０の非共振画像０Ｉを取得する。また、制御部２１は、取得した非共振画像０Ｉに基づいてＬＵＴを作成し、記憶部２２に記憶させる（ステップＳ５２）。
【００６１】
続いて、画像データ取得部２１１は、回転体５０の動作画像Ｉを取得する（ステップＳ５３）。具体的には、制御部２１は、回転角計測手段５１から回転角を示すデータを取得する。そして、制御部２１は、カメラ１１のフレームタイミングを制御して、予め定められた所定の解析位相φＮで、切削動作中の回転体５０を撮影する。これにより、回転体５０のフレームタイミングと実際の回転角との誤差を修正して、精度よく解析位相φＮの動作画像Ｉを撮影することができる。回転周期をＴ＾Ｒ（ｔｋ）、時刻ｔ＝ｔｋでの回転角の推定値をφ＾（ｔｋ）、目標角をφ０とすると、次のフレームタイミングである時刻ｔ＝ｔｋ＋１は、以下の式で表される。
【数３】
0006842218.tif 000004
【００６２】
演算部２１２は、撮影された基準位相φ０の位相ロックイン画像ＩφＮを生成する（ステップＳ５４）。また、演算部２１２は、生成された位相ロックイン画像ＩφＮに基づいて、振動検出を行う。振動検出方法は、実施の形態１と同様に、画像処理を用いて変位を算出することにより行う。
【００６３】
制御部２１は、検出した振動変位に基づいて、回転体５０の振動状態を解析する（ステップＳ５５）。例えば、算出された振動情報を高速フーリエ変換して、表示部２３に表示させる。これにより、所定の窓時間毎に振動状態を解析して、動画として表示させることにより、回転体５０の振動状態の時間変化を把握することができる。
【００６４】
以上説明したように、本実施の形態に係る振動解析システム１では、回転体５０の位相、すなわち回転角を制御部２１で検知してフレームタイミングを決定するフィードバック制御によって、動作画像Ｉを取得する。したがって、動作画像Ｉの位相誤差を低減することができる。これにより、回転体５０の振動解析の精度を向上させることができる。
【００６５】
本実施の形態では、回転体５０の回転角度を計測して、フレームタイミングを決定することとしたが、これに限られない。例えば、解析位相φＮで回転体５０を撮影するためのロックイン目標となる特徴量Ｇ（φ）を用いて、回転体５０の位相を計測し、フレームタイミングを決定することとしてもよい。より具体的には、画像データ取得部２１１は、カメラ１１又は制御ユニット２０に接続された他の高フレームレート撮影可能なカメラで、取得される画像データから、特徴量Ｇ（ｔｋ）を抽出して回転体５０の位相を計測し、動作画像Ｉを撮影することとしてもよい。これにより、画像処理によりフレームタイミングを決定することが可能となる。
【００６６】
特徴量Ｇに基づくフレームタイミングは、目標特徴量をＧ０とすると、以下の式のように表すことができる。
【数４】
0006842218.tif 000005
【００６７】
上記のように、特徴量Ｇを用いてフレームタイミングを調整することで、解析位相φＮである基準位相φ０に対して誤差の小さい動作画像Ｉを撮影することができる。したがって、回転角計測手段５１がなくとも、精度よく回転体５０の振動を検出することができる。
【００６８】
また、本実施の形態では、回転体５０の回転周期ＴＲにおいて、所定の基準位相φ０のみの動作画像Ｉを取得することとしたが、これに限られない。例えば、１周期をＮ分割して、複数の解析位相φＮで回転体５０を撮影することとしてもよい。この場合、解析位相φＮ毎に、位相ロックイン画像ＩφＮを作成し、振動解析を行えばよい。また、解析位相φＮ毎に、予め非共振画像０Ｉを取得してＬＵＴを作成し、解析位相φＮ毎の位相ロックイン画像ＩφＮを補正して合成することとしてもよい。これにより、実施の形態１と同様に、アンダーサンプリングによる不具合を低減して、より精度の高い振動解析を行うことができる。
【００６９】
上記各実施の形態では、回転体５０の回転軸に対して直交する方向の画像に基づいて振動解析することを基本としているが、これに限られない。例えば、回転軸方向の画像に対して振動解析を行うこととしてもよい。このような場合であっても、回転体５０の位相毎のロックイン画像に基づいて解析を行うことにより、回転成分の影響を低減して、精度よく共振成分の振動を検出することができる。
【００７０】
また、上記各実施の形態では、カメラ１１で撮影した画像を解析することとしたが、これに限られない。例えば、予め撮影され、制御ユニット２０に接続された外部記憶装置、記憶部２２等に記憶されている動作画像Ｉ、非共振画像０Ｉ等を用いて回転体５０の振動状態を解析することとしてもよい。
【産業上の利用可能性】
【００７１】
本発明は、振動を直接計測することが難しい回転体の振動解析に好適である。特に、外観形状が回転軸に対して不均一である回転体の振動解析に好適である。
【符号の説明】
【００７２】
１  振動解析システム、１１  カメラ、２０  制御ユニット、２１  制御部、２１１  画像データ取得部、２１２  演算部、２２  記憶部、２３  表示部、２４  入力部、５０  回転体
【要約】
【課題】回転体の回転に伴う振動成分と共振成分とを分離して、共振成分を精度よく検出できる振動解析システム、振動解析方法及びプログラムを提供する。
【解決手段】振動解析システムは、解析対象である回転体５０の動作状態を表す動作画像Ｉを取得する画像データ取得部と、動作画像Ｉに基づいて、回転体５０の振動状態を解析する演算部と、を備える。また演算部は、動作画像Ｉのうち、回転体５０の位相が所定の解析位相φＮである画像からなる時系列の位相ロックイン画像ＩφＮに基づいて、回転体５０の振動状態を解析する。
【選択図】図２【図１】
0006842218.tif 000006
【図２】
0006842218.tif 000007
【図３】
0006842218.tif 000008
【図４】
0006842218.tif 000009
【図５】
0006842218.tif 000010
【図６】
0006842218.tif 000011
【図７】
0006842218.tif 000012
【図８】
0006842218.tif 000013
【図９】
0006842218.tif 000014
【図１０】
0006842218.tif 000015