[go: up one dir, main page]

JP3824546B2 - Rotational accuracy measuring device - Google Patents

Rotational accuracy measuring device Download PDF

Info

Publication number
JP3824546B2
JP3824546B2 JP2002036777A JP2002036777A JP3824546B2 JP 3824546 B2 JP3824546 B2 JP 3824546B2 JP 2002036777 A JP2002036777 A JP 2002036777A JP 2002036777 A JP2002036777 A JP 2002036777A JP 3824546 B2 JP3824546 B2 JP 3824546B2
Authority
JP
Japan
Prior art keywords
data
displacement
sampling
rotating body
sampling point
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 - Lifetime
Application number
JP2002036777A
Other languages
Japanese (ja)
Other versions
JP2003240536A (en
Inventor
恭之 脇田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JTEKT Corp
Original Assignee
JTEKT Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by JTEKT Corp filed Critical JTEKT Corp
Priority to JP2002036777A priority Critical patent/JP3824546B2/en
Publication of JP2003240536A publication Critical patent/JP2003240536A/en
Application granted granted Critical
Publication of JP3824546B2 publication Critical patent/JP3824546B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、所定軸を中心として回転する回転体の回転精度を測定する回転精度測定装置に関する。
【0002】
【従来の技術】
転がり軸受けや流体軸受け等に支持されて回転する回転体(例えばモータによって回転するスピンドル)の回転精度を測定する回転精度測定装置では、従来、回転体の変位を示す時系列データを回転体の回転とは無関係に採取して大量の変位データを収集し、その大量の変位データに対して一括してデータ処理を行っていた。このデータ処理では、例えば、時系列的に得られた変位データ(以下「採取データ」という)に対して高速フーリエ変換(以下「FFT」という)を施すことにより、周波数スペクトルを示すデータ(以下「スペクトルデータ」という)を算出し、そのスペクトルデータにおいて同期振れ誤差(RRO:Repeatable Run Out)を除去することにより、1回転毎に繰り返されない回転体の変位に相当する非同期振れ誤差(NRRO:Non Repeatable Run Out)を求めている。
【0003】
【発明が解決しようとする課題】
しかし、従来の回転精度測定装置では、回転体を回転させるモータ等の駆動源の回転ムラが測定に大きな影響を与えていた。また、上記の採取データは、回転体の回転とは無関係に採取された時系列データであるので、上記データ処理においてFFTが施される際には、採取データが所定区間外で滑らかに0になるように、予め窓関数によって採取データに対して適当な重み付けがなされていた。一般に窓関数が使用されると、FFTによって得られるスペクトルが拡散するので、この窓関数の使用も、従来の回転精度測定装置による測定に大きな影響を与えることになる。このように従来の測定装置では、モータの回転ムラや窓関数の影響が大きいために、回転精度についての高精度な測定、すなわちRROやNRRO、真円度等の正確な算出は困難であった。
【0004】
これに対し本願出願人は、回転体の変位を示す標本化値からなる時系列データとしての変位データを回転周期毎に分離してブロック化することによりブロック化データを生成し、そのブロック化データにおける各データブロックを構成する標本化値の数をレート変換によって同一とし、レート変換後のブロック化データを用いてNRRO等を算出するという周期分離型信号処理に基づく回転精度測定装置の開発を進めている。この周期分離型信号処理によれば、回転ムラや窓関数の使用等による測定への影響が低減されて正確なNRRO等の算出が可能となる。
【0005】
上記の周期分離型信号処理では、回転体を回転させるモータ等を含む駆動部から回転周期に同期した信号(例えば1回転毎に1個パルスが現れる信号)がインデックスパルスとして出力される場合には、そのインデックスパルスを使用して変位データを回転周期毎に分離することができる。これに対し、インデックスパルスのような回転周期に同期した信号を取得できない場合には、例えば変位データから直流成分を除去した後にゼロ点検出を行うことにより、変位データを回転周期毎に分離する必要がある。しかし、ノイズ等の外乱やA/D変換器の誤差等により、正確にゼロ点を求めるのは困難であり、その結果、変位データを正確に回転周期毎に分離することができないという問題があった。
【0006】
そこで本発明では、正確な回転精度測定のために変位データを正確に回転周期毎に分離することができる回転精度測定装置を提供することを目的とする。
【0007】
【課題を解決するための手段および発明の効果】
第1の発明は、所定軸を中心として回転する回転体の回転精度を測定する回転精度測定装置であって、
前記回転体の半径方向の変位または前記所定軸方向の変位を示す信号の標本化値からなる時系列データを変位データとして取得するデータ取得手段と、
前記変位データを構成する各標本化値に対する移動平均を算出し、当該移動平均からなる時系列データである平滑化変位データを生成する移動平均手段と、
前記平滑化変位データの表す信号波形におけるゼロ点に基づき、前記データ取得手段によって取得された変位データを前記回転体の回転周期毎に分離してブロック化するブロック化手段と、
前記ブロック化手段によってブロック化された変位データにおける各データブロックを構成する標本化値の数が同一となるように、当該各データブロックに対してレート変換を行う変換手段と、
前記レート変換後の前記データブロックからなるブロック化データに基づき、前記回転精度を示す指標を算出する算出手段と
を備え、
前記ブロック化手段は、
前記平滑化変位データに基づき、前記回転体の変位を示す標本化値のサンプリング点における隣接サンプリング点対のうち、一方のサンプリング点における移動平均が正値で他方のサンプリング点における移動平均が負値である隣接サンプリング点対を検出するサンプリング点対検出手段と、
前記検出された隣接サンプリング点対の間に前記ゼロ点が存在するとみなして、前記データ取得手段によって取得された変位データを前記回転周期毎に分離する分離手段とを含むことを特徴とする。
【0008】
このような第1の発明によれば、回転体の変位を表す変位データが回転周期毎に分離されてブロック化され、ブロック化された変位データに基づき回転精度を示す指標が算出される。このような処理において、変位データの回転周期毎の分離は、変位データにおける各標本化値に対する移動平均からなる平滑化変位データに基づいて行われるので、変位データを正確に回転周期毎に分離することができる。このとき、回転体の変位を示す標本化値のサンプリング点における隣接サンプリング点対のうち、その間にゼロ点が存在するとみなせる隣接サンプリング点対が検出され、その隣接サンプリング点対に基づいて変位データが回転周期毎に分離される。したがって、正確なゼロ点の位置を求めることなく、効率よく正確に変位データを回転周期毎に分離することができる。このようにして変位データが正確に回転周期毎に分離されてブロック化され、更に、各データブロックを構成する標本化値の数が同一となるようにレート変換が行われる。そして、そのレート変換後のブロック化データに基づいて回転精度を示す指標が算出される。このため回転精度測定において、回転体を回転させるモータ等の回転ムラの影響を軽減することができる。また、変位データが回転周期毎にブロック化されると、そのブロック化された変位データに対して窓関数を使用することなくFFTを施すことにより、回転体の変位における同期成分およびその高調波成分を広がりの無い線スペクトルとして得ることができる。このため、そのFFTによって得られたスペクトルデータからNRROを求める際には、RROを他に影響を与えることなく精度よく除去することができる。
【0013】
【発明の実施の形態】
以下、添付図面を参照して本発明の実施形態について説明する。
<構成および動作>
図1は、本発明の一実施形態に係る回転精度測定装置の構成を模式的に示す構成図であり、図2は、この回転精度測定装置による回転体の測定を示す平面図である。
【0014】
この回転精度測定装置は、軸受けによって支持され所定の回転軸11を中心としてモータ(不図示)により回転する回転体10の回転精度を測定する装置であって、回転体10の半径方向の変位を検出する非接触式の変位センサ20と、変位センサ20の検出信号Sdに基づき回転体10の半径方向の変位を示す標本化値からなる時系列データを収集し、その時系列データである採取データを処理するデータ収集処理装置30とを備えている。
【0015】
変位センサ20は、回転体10の外周面近傍に配置され、その外周面と変位センサ20との距離を検出し、その検出結果を回転体10の半径方向の変位を示す検出信号Sdとして出力する。なお、回転体10の軸方向の変位を検出し、軸方向の変位についての回転精度を測定する場合には、図1において点線で示すように、回転体10の上面近傍に回転体10の軸方向の変位を検出する非接触式の変位センサ21を配置し、その変位センサ21から出力される信号すなわち回転体10の上面と変位センサ21との距離を示す信号を検出信号Sdとして使用すればよい。ただし以下では、説明の便宜上、半径方向の変位を検出する変位センサ20のみが配置されているものとする。
【0016】
データ収集処理装置30は、中央処理装置としてのCPU31、入力インターフェース部32、メモリ33、および表示制御部34をバスで接続した構成となっており、表示制御部34には表示部36が接続されている。上記変位センサ20からの検出信号Sdは、このデータ収集処理装置30における入力インターフェース部32に入力される。入力インターフェース部32はA/D変換器を有し、上記検出信号Sdは、このA/D変換器によって標本化され、デジタルデータとしてメモリ33に一時的に格納される(以下、このデジタルデータを「原データ」という)。この入力インターフェース部32は、上記の変位センサ20と共に、回転体10の変位を表す時系列データの取得手段を構成する。
【0017】
CPU31は、予めメモリ33に格納された所定プログラムを実行することにより、上記の原データに対し、後述のDCカット処理、周期分離処理、レート変換、FFT等のデータ処理を順に施す。これにより、データ収集処理装置30は、図3に示すような機能的構成の装置、すなわち、DCカット部111と、周期分離部112と、レート変換部113と、信号処理部114とを備える装置として動作する。
【0018】
DCカット部111は、変位センサ20から出力される検出信号Sdから直流成分を除去するものであり、具体的には、検出信号Sdを表すデジタルデータである原データに対する信号処理によって、検出信号Sdから直流成分を除去した信号を表すデジタルデータを採取データDaとして生成する(この信号処理を「DCカット処理」という)。なお、このDCカット部111はソフトウェア的に実現されているが、入力インターフェース部32に直流成分遮断回路またはオフセット調整回路を設け、これによって検出信号Sdから直流成分を除去した後に、A/D変換器によって上記の採取データDaを作成するようにしてもよい。この場合、DCカット部111は、ハードウェアとして実現されることになり、入力インターフェース部32の一部を構成する。
【0019】
周期分離部112は、DCカット部111によって生成された採取データDaを回転体10の回転周期毎に分離してブロック化することにより、複数のデータブロックからなるブロック化データDbを生成する(この処理を「周期分離処理」という)。具体的には、検出信号Sdを表すデジタルデータである採取データDaから、回転体10の変位量が0となる時点であるゼロ点に基づき、採取データDaを回転周期毎に分離してブロック化データDbを得る。
【0020】
図4は、上記周期分離処理の手順を示すフローチャートである。本実施形態では、CPU31がこれらの手順を実行することにより周期分離部112がソフトウェア的に実現される。
【0021】
上記周期分離処理ではCPU31は、まず、DCカット部111から出力される変位データとしての採取データDaを構成する各標本化値に対し、移動平均を算出する(ステップS10)。すなわち、採取データDaにおける各標本化値に順次注目し、注目標本化値のサンプリング点(以下「注目サンプリング点」という)を中心とする所定期間における複数サンプリング点の標本化値の平均値を求め、その平均値を注目標本化値に対する移動平均とする。このようにして、採取データDaにおける各標本化値に対する移動平均からなる時系列データを得る(以下、この時系列データを「平滑化変位データ」という)。
【0022】
回転体10の変位を示す検出信号Sdから図5に示すような信号波形を表す採取データDaが得られた場合、その信号波形の一部を拡大すると、その信号波形は、例えば図6(a)に示すように不規則に変動するノイズ成分を含み滑らかな波形とはなっていない。また、そのような波形の信号から得られるデジタルデータである採取データDaには、A/D変換器に起因する誤差(以下「A/D変換誤差」という)も含まれる。しかし、上記のような移動平均による処理を施せば、ノイズおよびA/D変換誤差が相殺されて除去され、図6(b)に示すように平滑化された波形の信号を表すデジタルデータが平滑化変位データとして得られる。なお、図5および図6において一定間隔で描かれている縦線はサンプリング点に対応しているものとする。
【0023】
このような平滑化変位データが得られると、CPU31は、その平滑化変位データの表す図6(b)に示すような信号波形においてゼロ点を検出する(ステップS12)。ここで、ゼロ点には、信号が負値から正値へと変化するときのゼロ点(以下「立上りゼロ点」という)と、信号が正値から負値へと変化するときのゼロ点(以下「立下りゼロ点」という)とがある。採取データDaを回転周期毎に分離するためには、これら2種類のゼロ点のうちいずれか一方のみを検出すればよいが、以下では、立上りゼロ点のみを検出するものとして説明する。また、採取データDaを回転周期毎に分離するためには、ゼロ点を正確に求める必要はなく、ゼロ点を挟んで隣接する2つのサンプリング点を求めればよい。そこで、本実施形態では、第1のサンプリング点とそれに続く第2のサンプリング点とからなるサンプリング点対であって、第1のサンプリング点の標本化値に対する移動平均が負値で第2のサンプリング点の標本化値に対する移動平均が正値となるサンプリング点対(以下「立上りゼロ交差サンプリング点対」という)を求める。この立上りゼロ交差サンプリング点対は、上記の平滑化変位データにおいて隣接する2つのサンプリング点での移動平均の正負を順次調べれば容易に求めることができる。
【0024】
図6(b)に示した例では、サンプリング点対(ta1,ta2)と、サンプリング点対(tc1,tc2)とが、立上りゼロ交差サンプリング点対として検出される。この場合、立上りゼロ交差サンプリング点対(ta1,ta2)の間には立上りゼロ点Zpaが存在し、立上りゼロ交差サンプリング点対(tc1,tc2)の間には次の立上りゼロ点Zpcが存在する。したがって、サンプリング点ta2からtc1までのサンプリング点(ta2およびtc1も含む)が回転体10の1回転周期に対応する。なお、サンプリング点対(tb1,tb2)は立下りゼロ交差サンプリング点対であって、その間には立下りゼロ点Zpbが存在する。このような立下りゼロ点を検出することによっても採取データDaを回転周期毎に分離できるが、既述のように本実施形態では、立上りゼロ点対のみを検出するものとしている。
【0025】
次にCPU31は、検出された立上りゼロ点により採取データDaを分離して1回転周期毎にブロック化し、複数のデータブロックからなるブロック化データDbを得る(ステップS14)。図6(b)に示した例では、立上りゼロ交差サンプリング点対(ta1,ta2)の間には立上りゼロ点Zpaが、立上りゼロ交差サンプリング点対(tc1,tc2)の間には次の立上りゼロ点Zpcが、それぞれ存在するので、サンプリング点ta1以前の標本化値とサンプリング点ta2以降の標本化値とを分離すると共に、サンプリング点tc1以前の標本化値とサンプリング点tc2以降の標本化値とを分離し、サンプリング点ta2からtc1までのサンプリング点(ta2およびtc1も含む)の標本化値を1つのデータブロックとしてブロック化することになる。
【0026】
上記のような周期分離処理を図5に示した採取データDaに施すと、図7に示すようにそれぞれが1回転周期分の変位データに相当する4個のブロックデータDb1〜Db4からなるブロック化データDbが得られる。このようにして得られるブロック化データDbにおける各データブロックを構成する標本化値の個数(これは1回転周期のサンプリング点数であり、以下「データ数」という)は、モータの回転ムラ等のため、通常、全て等しくはならない。例えば図8(a)に示すように、データブロックDb1,Db2,Db3,Db4のデータ数は、それぞれ異なるn1個,n2個,n3個,n4個となる。
【0027】
レート変換部113は、上記のようにデータ数のばらつく複数のデータブロックからなるブロック化データDbに対して補間処理を行うことにより、各データブロックのデータ数を同一にする。すなわち、レート変換によって各データブロックのサンプリング点数を同一にする。このとき、信号処理部114で実行されるFFT(高速フーリエ変換)を考慮して、各データブロックのデータ数を2のべき乗とする。例えば図8(a)に示したようなブロック化データDbに対してレート変換を施すことにより、図8(b)に示すようにデータ数が全て2mであるデータブロックDc1,Dc2,Dc3,Dc4からなるブロック化データDcを得る。なお、ここでは、ブロック化データDbを構成する4個のデータブロックDb1,Db2,Db3,Db4は、ブロック化データDcを構成する4個のデータブロックDc1,Dc2,Dc3,Dc4にそれぞれ変換されるものとする。
【0028】
信号処理部114は、回転精度を示す指標の算出手段であって、レート変換後のブロック化データDcに対して窓関数を使用することなくFFTを施すことにより、スペクトルデータを算出する。そして、このスペクトルデータに基づき、従来と同様の手法により、RROや、NRRO、真円度等を求める。
【0029】
このようにして得られたRROや、NRRO、真円度等、回転体10の回転精度を示す指標は、測定結果としてメモリ33に格納されると共に、他の所定プログラムに基づき、表示制御部34に送られて、その表示制御部34によって表示部36に表示される。
【0030】
<効果>
上記のような本実施形態によれば、採取データDaが回転周期毎に分離されることにより回転周期毎にブロック化されたブロック化データDbが得られ、更に、そのブロック化データDbに対するレート変換によって各データブロックのデータ数が同一となる。これにより、回転精度測定において、回転体10を回転させるモータの回転ムラの影響を低減することができる。このため、RROや、NRRO、真円度等、回転精度を示す指標を正確に求めることができる。また、回転周期毎に分離されてブロック化されたデータに対して窓関数を使用することなくFFTを施すことにより、回転体10の1回転周期の逆数を基本周波数とする検出信号Sdの高調波成分を広がりの無い線スペクトルとして求めることができる。このため、そのFFTを施すことにより得られるスペクトルデータから、RROを精度よく求めることができ、そのスペクトルデータからNRROを求める際には、RROを他に影響を与えることなく精度よく除去することができる。
【0031】
ところで本実施形態では、周期分離処理において、採取データDaの回転周期毎の分離はゼロ点検出に基づいて行われるが、そのときゼロ点検出には、採取データDaの各標本化値に対する移動平均からなる平滑化変位データが使用される。このため、採取データDaを正確に回転周期毎に分離することができる。すなわち、ゼロ点検出に採取データDaをそのまま用いると、ノイズやA/D変換誤差等のために1つのゼロ点(真のゼロ点)の近傍に複数のゼロ点が見かけ上存在することがあるが、平滑化変位データを用いることにより、本来のゼロ点に1対1に対応して正確にゼロ点を検出することができる。例えば、図6(a)に示した採取データDaからゼロ点を検出すると、サンプリング点対(ta1,ta2)の間とサンプリング点対(ta2,ta3)の間との2区間に立上りゼロ点が見かけ上存在し、真の立上りゼロ点を含む立上りゼロ交差サンプリング点対を検出できない場合がある。このため、採取データDaを正確に回転周期毎に分離することができない。これに対し、図6(b)に示した平滑化変位データからゼロ点を検出すると、サンプリング点対(ta1,ta2)の間に立上りゼロ点Zpaが存在することが検出され、その近傍では他にゼロ点は検出されない。このため、採取データDaを正確に回転周期毎に分離することが可能となる。
【0032】
なお、上述の周期分離処理の手順からわかるように、採取データDaを分離するためには、ゼロ点の正確な位置を求める必要はなく、ゼロ点の存在するサンプリング点対を求めれば十分である。すなわち、立上りゼロ交差サンプリング点対または立下りゼロ交差サンプリング点対のいずれか一方を検出すれば、検出されたサンプリング点対の間にゼロ点が存在するとみなして採取データDaを回転周期毎に分離することができる。また、上記周期分離処理では、採取データDaの各標本化値に対して、その標本化値のサンプリング点を中心とする所定期間のサンプリング点の標本化値の平均を算出するが、その所定期間を適切に設定する必要がある。本実施形態では、この所定期間を、回転体10の回転周期や、サンプリング周波数、ノイズのレベルおよび周波数等を考慮して、本来のゼロ点に1対1に対応するゼロ点が平滑化変位データから検出されるように適切な所定期間を選定するものとする。
【0033】
<変形例>
上記実施形態では、周期分離処理によって得られる各データブロックのデータ数を同一にするためにレート変換が行われるが、これに代えて、次のような処理を行ってもよい。すなわち、周期分離処理によって採取データDaを回転周期毎に分離して得られる各データブロックのデータ数のうち最も出現頻度の高いデータ数を求め、最も出現頻度の高いデータ数を有するデータブロックのみを採取データDaから抽出する。そして、このようにして抽出されたデータブロックからなるブロック化データを用いてNRRO等の回転精度を示す指標を算出する。このような変形例によっても、回転精度測定において、回転体10を回転させるモータの回転ムラ等の影響を低減することができる。
【図面の簡単な説明】
【図1】本発明の一実施形態に係る回転精度測定装置の構成を模式的に示す構成図である。
【図2】上記実施形態に係る回転精度測定装置による回転体の測定を示す平面図である。
【図3】上記実施形態に係る回転精度測定装置におけるデータ処理装置の機能的構成を示すブロック図である。
【図4】上記実施形態における周期分離処理の手順を示すフローチャートである。
【図5】上記実施形態における採取データを説明するための信号波形図である。
【図6】上記実施形態における周期分離処理を説明するための信号波形図である。
【図7】上記実施形態において採取データに対する周期分離処理によって得られるブロック化データを説明するための信号波形図である。
【図8】上記実施形態におけるレート変換を説明するためのデータ構成図である。
【符号の説明】
10 …回転体
11 …回転軸
20,21…変位センサ
30 …データ収集処理装置
31 …CPU
32 …入力インターフェース部
33 …メモリ
111 …DCカット部
112 …周期分離部
113 …レート変換部
114 …信号処理部
Sd …検出信号
Da …採取データ(変位データ)
Db …ブロック化データ
Dc …レート変換後のブロック化データ
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a rotational accuracy measuring device that measures rotational accuracy of a rotating body that rotates about a predetermined axis.
[0002]
[Prior art]
In a rotational accuracy measuring device that measures the rotational accuracy of a rotating body (for example, a spindle that is rotated by a motor) that is supported by a rolling bearing or a fluid bearing and rotates, time series data indicating the displacement of the rotating body is conventionally used to rotate the rotating body. A large amount of displacement data was collected regardless of the amount of data collected, and data processing was collectively performed for the large amount of displacement data. In this data processing, for example, fast Fourier transform (hereinafter referred to as “FFT”) is performed on the displacement data (hereinafter referred to as “collected data”) obtained in time series, whereby data indicating a frequency spectrum (hereinafter referred to as “FFT”). By calculating the spectral data (hereinafter referred to as “spectrum data”) and removing the synchronous shake error (RRO: Repeatable Run Out) in the spectral data, the asynchronous shake error (NRRO: Non) corresponding to the displacement of the rotating body that is not repeated every rotation. Repeatable Run Out).
[0003]
[Problems to be solved by the invention]
However, in the conventional rotational accuracy measuring device, rotation unevenness of a driving source such as a motor that rotates a rotating body has a great influence on the measurement. Further, since the collected data is time-series data collected regardless of the rotation of the rotating body, when the FFT is performed in the data processing, the collected data is smoothly set to 0 outside a predetermined interval. As shown, appropriate weighting is performed on the collected data by a window function in advance. In general, when a window function is used, the spectrum obtained by the FFT is diffused. Therefore, the use of this window function also greatly affects the measurement by the conventional rotational accuracy measuring device. As described above, in the conventional measuring apparatus, since the influence of the rotation irregularity of the motor and the window function is large, it is difficult to accurately measure the rotation accuracy, that is, to accurately calculate RRO, NRRO, roundness, and the like. .
[0004]
On the other hand, the applicant of the present application generates block data by separating the displacement data as time-series data consisting of sampled values indicating the displacement of the rotating body for each rotation period, and blocking the data. Development of a rotational accuracy measurement device based on periodic separation type signal processing, in which the number of sampling values constituting each data block in the same is made the same by rate conversion, and NRRO and the like are calculated using the blocked data after rate conversion ing. According to this period-separated signal processing, the influence on measurement due to rotation unevenness, use of a window function, etc. is reduced, and accurate calculation of NRRO and the like becomes possible.
[0005]
In the periodic separation type signal processing described above, when a signal synchronized with the rotation period (for example, a signal in which one pulse appears for each rotation) is output as an index pulse from a drive unit including a motor that rotates a rotating body. The index pulse can be used to separate the displacement data for each rotation period. On the other hand, when a signal synchronized with the rotation cycle such as an index pulse cannot be obtained, it is necessary to separate the displacement data for each rotation cycle, for example, by detecting the zero point after removing the DC component from the displacement data. There is. However, it is difficult to accurately determine the zero point due to disturbances such as noise and errors of the A / D converter, and as a result, there is a problem that the displacement data cannot be separated accurately for each rotation period. It was.
[0006]
Therefore, an object of the present invention is to provide a rotation accuracy measuring device capable of accurately separating displacement data for each rotation cycle for accurate rotation accuracy measurement.
[0007]
[Means for Solving the Problems and Effects of the Invention]
A first invention is a rotational accuracy measuring device for measuring rotational accuracy of a rotating body that rotates about a predetermined axis,
Data acquisition means for acquiring, as displacement data, time-series data consisting of sampling values of signals indicating displacement in the radial direction of the rotating body or displacement in the predetermined axis direction;
A moving average means for calculating a moving average for each sampling value constituting the displacement data, and generating smoothed displacement data that is time-series data composed of the moving average;
Blocking means for separating and blocking the displacement data acquired by the data acquisition means for each rotation period of the rotating body based on the zero point in the signal waveform represented by the smoothed displacement data;
Conversion means for performing rate conversion on each data block so that the number of sampling values constituting each data block in the displacement data blocked by the blocking means is the same;
Calculation means for calculating an index indicating the rotation accuracy based on the blocked data composed of the data blocks after the rate conversion ;
The blocking means includes
Based on the smoothed displacement data, the moving average at one sampling point is a positive value and the moving average at the other sampling point is a negative value among the adjacent sampling point pairs at the sampling point of the sampling value indicating the displacement of the rotating body. Sampling point pair detection means for detecting adjacent sampling point pairs,
Separating means for separating the displacement data acquired by the data acquiring means for each rotation period on the assumption that the zero point exists between the detected pair of adjacent sampling points .
[0008]
According to the first aspect of the invention, the displacement data representing the displacement of the rotating body is separated and blocked for each rotation period, and an index indicating the rotation accuracy is calculated based on the blocked displacement data. In such processing, since the separation of the displacement data for each rotation period is performed based on the smoothed displacement data consisting of the moving average for each sampled value in the displacement data, the displacement data is accurately separated for each rotation period. be able to. At this time, among the adjacent sampling point pairs at the sampling point of the sampling value indicating the displacement of the rotating body, an adjacent sampling point pair that can be regarded as having a zero point between them is detected, and the displacement data is based on the adjacent sampling point pair. Separated every rotation cycle. Therefore, it is possible to efficiently and accurately separate the displacement data for each rotation period without obtaining an accurate zero point position. In this way, the displacement data is accurately separated for each rotation period and is divided into blocks, and further rate conversion is performed so that the number of sampling values constituting each data block is the same. Then, an index indicating the rotation accuracy is calculated based on the blocked data after the rate conversion. Therefore, in the rotation accuracy measurement, it is possible to reduce the influence of rotation unevenness such as a motor that rotates the rotating body. Further, when the displacement data is blocked for each rotation period, the synchronous component and its harmonic component in the displacement of the rotating body are obtained by performing FFT on the blocked displacement data without using a window function. Can be obtained as a line spectrum without spread. For this reason, when obtaining the NRRO from the spectrum data obtained by the FFT, the RRO can be accurately removed without affecting the other.
[0013]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings.
<Configuration and operation>
FIG. 1 is a configuration diagram schematically showing a configuration of a rotational accuracy measuring device according to an embodiment of the present invention, and FIG. 2 is a plan view showing measurement of a rotating body by the rotational accuracy measuring device.
[0014]
This rotational accuracy measuring device is a device that measures the rotational accuracy of a rotating body 10 supported by a bearing and rotated by a motor (not shown) around a predetermined rotating shaft 11. Based on the non-contact type displacement sensor 20 to be detected and the sampled value indicating the displacement in the radial direction of the rotating body 10 based on the detection signal Sd of the displacement sensor 20, time series data is collected, and the collected data which is the time series data is collected. And a data collection processing device 30 for processing.
[0015]
The displacement sensor 20 is disposed in the vicinity of the outer peripheral surface of the rotator 10, detects the distance between the outer peripheral surface and the displacement sensor 20, and outputs the detection result as a detection signal Sd indicating the radial displacement of the rotator 10. . When detecting the axial displacement of the rotating body 10 and measuring the rotational accuracy of the axial displacement, the axis of the rotating body 10 is located near the upper surface of the rotating body 10 as shown by the dotted line in FIG. If a non-contact type displacement sensor 21 for detecting the displacement in the direction is arranged and a signal output from the displacement sensor 21, that is, a signal indicating the distance between the upper surface of the rotating body 10 and the displacement sensor 21 is used as the detection signal Sd. Good. However, in the following, for convenience of explanation, it is assumed that only the displacement sensor 20 for detecting the displacement in the radial direction is arranged.
[0016]
The data collection processing device 30 has a configuration in which a CPU 31 as a central processing unit, an input interface unit 32, a memory 33, and a display control unit 34 are connected by a bus, and a display unit 36 is connected to the display control unit 34. ing. The detection signal Sd from the displacement sensor 20 is input to the input interface unit 32 in the data collection processing device 30. The input interface unit 32 includes an A / D converter, and the detection signal Sd is sampled by the A / D converter and temporarily stored as digital data in the memory 33 (hereinafter, this digital data is referred to as digital data). "Original data"). The input interface unit 32, together with the displacement sensor 20, constitutes time series data acquisition means representing the displacement of the rotating body 10.
[0017]
The CPU 31 sequentially executes data processing such as DC cut processing, period separation processing, rate conversion, and FFT, which will be described later, on the original data by executing a predetermined program stored in the memory 33 in advance. Accordingly, the data collection processing device 30 has a functional configuration as shown in FIG. 3, that is, a device including a DC cut unit 111, a period separation unit 112, a rate conversion unit 113, and a signal processing unit 114. Works as.
[0018]
The DC cut unit 111 removes a direct current component from the detection signal Sd output from the displacement sensor 20, and specifically, the detection signal Sd is performed by signal processing on original data that is digital data representing the detection signal Sd. The digital data representing the signal from which the DC component is removed from is generated as the collected data Da (this signal processing is referred to as “DC cut processing”). Although the DC cut unit 111 is realized by software, a DC component cutoff circuit or an offset adjustment circuit is provided in the input interface unit 32, thereby removing the DC component from the detection signal Sd, and then A / D conversion. The sampling data Da may be created by a container. In this case, the DC cut unit 111 is realized as hardware, and constitutes a part of the input interface unit 32.
[0019]
The period separation unit 112 generates block data Db composed of a plurality of data blocks by separating the collection data Da generated by the DC cut unit 111 into blocks for each rotation period of the rotating body 10 (this is the block data Db). Processing is called “periodic separation processing”). Specifically, from the collection data Da that is digital data representing the detection signal Sd, the collection data Da is separated into blocks based on the zero point, which is the time when the displacement of the rotating body 10 becomes zero, and is blocked. Data Db is obtained.
[0020]
FIG. 4 is a flowchart showing the procedure of the period separation process. In this embodiment, when the CPU 31 executes these procedures, the period separation unit 112 is realized by software.
[0021]
In the periodic separation process, the CPU 31 first calculates a moving average for each sampled value constituting the collection data Da as displacement data output from the DC cut unit 111 (step S10). That is, paying attention to each sampled value in the sampling data Da in sequence, an average value of sampled values of a plurality of sampling points in a predetermined period centering on a sampling point of the sampled value of interest (hereinafter referred to as “targeted sampling point”) is obtained. The average value is a moving average for the sampled value of interest. In this way, time-series data consisting of moving averages for each sampled value in the collected data Da is obtained (hereinafter, this time-series data is referred to as “smoothed displacement data”).
[0022]
When the sampling data Da representing the signal waveform as shown in FIG. 5 is obtained from the detection signal Sd indicating the displacement of the rotating body 10, when a part of the signal waveform is enlarged, the signal waveform is, for example, FIG. As shown in (), a smooth waveform including a noise component that fluctuates irregularly is not obtained. Further, the collected data Da, which is digital data obtained from a signal having such a waveform, includes an error caused by an A / D converter (hereinafter referred to as “A / D conversion error”). However, if the above-described moving average processing is performed, noise and A / D conversion errors are canceled out and removed, and digital data representing a smoothed waveform signal as shown in FIG. 6B is smoothed. Obtained as converted displacement data. Note that vertical lines drawn at regular intervals in FIGS. 5 and 6 correspond to sampling points.
[0023]
When such smoothed displacement data is obtained, the CPU 31 detects a zero point in the signal waveform shown in FIG. 6B represented by the smoothed displacement data (step S12). Here, the zero point includes a zero point when the signal changes from a negative value to a positive value (hereinafter referred to as “rising zero point”) and a zero point when the signal changes from a positive value to a negative value ( Hereinafter referred to as “falling zero point”). In order to separate the collection data Da for each rotation period, only one of these two types of zero points may be detected. However, in the following, description will be made assuming that only the rising zero point is detected. Further, in order to separate the collection data Da for each rotation cycle, it is not necessary to accurately obtain the zero point, and two adjacent sampling points may be obtained with the zero point interposed therebetween. Therefore, in the present embodiment, the second sampling point pair is a sampling point pair including a first sampling point and a second sampling point that follows the first sampling point, and the moving average with respect to the sampling value of the first sampling point is a negative value. A sampling point pair (hereinafter referred to as “rising zero-crossing sampling point pair”) whose moving average with respect to the sampled value of the point is a positive value is obtained. The rising zero-crossing sampling point pair can be easily obtained by sequentially examining the positive and negative of the moving average at two adjacent sampling points in the smoothed displacement data.
[0024]
In the example shown in FIG. 6B, the sampling point pair (ta1, ta2) and the sampling point pair (tc1, tc2) are detected as rising zero-crossing sampling point pairs. In this case, the rising zero point Zpa exists between the rising zero crossing sampling point pair (ta1, ta2), and the next rising zero point Zpc exists between the rising zero crossing sampling point pair (tc1, tc2). . Therefore, the sampling points (including ta2 and tc1) from the sampling points ta2 to tc1 correspond to one rotation period of the rotating body 10. The sampling point pair (tb1, tb2) is a falling zero crossing sampling point pair, and a falling zero point Zpb exists between them. Although the sampling data Da can be separated for each rotation cycle by detecting such a falling zero point, as described above, only the rising zero point pair is detected in the present embodiment.
[0025]
Next, the CPU 31 separates the collected data Da by the detected rising zero point and blocks it for each rotation cycle, thereby obtaining the blocked data Db composed of a plurality of data blocks (step S14). In the example shown in FIG. 6B, the rising zero point Zpa is between the rising zero-crossing sampling point pair (ta1, ta2), and the next rising edge is between the rising zero-crossing sampling point pair (tc1, tc2). Since the zero point Zpc exists, the sampling value before the sampling point ta1 and the sampling value after the sampling point ta2 are separated, and the sampling value before the sampling point tc1 and the sampling value after the sampling point tc2 are separated. And sampling values of sampling points (including ta2 and tc1) from sampling points ta2 to tc1 are blocked as one data block.
[0026]
When the period separation process as described above is applied to the sampling data Da shown in FIG. 5, as shown in FIG. 7, each block is made up of four block data Db1 to Db4 corresponding to displacement data for one rotation period. Data Db is obtained. The number of sampled values constituting each data block in the block data Db thus obtained (this is the number of sampling points in one rotation cycle, hereinafter referred to as “data number”) is due to uneven rotation of the motor. , Usually not all equal. For example, as shown in FIG. 8A, the data numbers of the data blocks Db1, Db2, Db3, and Db4 are n1, n2, n3, and n4, respectively.
[0027]
The rate conversion unit 113 makes the number of data of each data block the same by performing interpolation processing on the blocked data Db composed of a plurality of data blocks in which the number of data varies as described above. That is, the number of sampling points in each data block is made the same by rate conversion. At this time, considering the FFT (Fast Fourier Transform) executed by the signal processing unit 114, the number of data in each data block is set to a power of two. For example, by performing rate conversion on the blocked data Db as shown in FIG. 8 (a), data blocks Dc1, Dc2, Dc3, all of which are 2 m in number as shown in FIG. 8 (b). Blocked data Dc consisting of Dc4 is obtained. Here, the four data blocks Db1, Db2, Db3, Db4 constituting the blocked data Db are converted into the four data blocks Dc1, Dc2, Dc3, Dc4 constituting the blocked data Dc, respectively. Shall.
[0028]
The signal processing unit 114 is a means for calculating an index indicating rotation accuracy, and calculates spectrum data by performing FFT on the blocked data Dc after rate conversion without using a window function. And based on this spectrum data, RRO, NRRO, roundness, etc. are calculated | required by the method similar to the past.
[0029]
Indices indicating the rotational accuracy of the rotating body 10 such as RRO, NRRO, and roundness obtained in this manner are stored in the memory 33 as measurement results, and based on other predetermined programs, the display control unit 34. And is displayed on the display unit 36 by the display control unit 34.
[0030]
<Effect>
According to the present embodiment as described above, the sampling data Da is separated for each rotation period, whereby the blocked data Db that is blocked for each rotation period is obtained, and the rate conversion for the blocked data Db is further performed. Thus, the number of data in each data block becomes the same. Thereby, in rotation accuracy measurement, the influence of the rotation nonuniformity of the motor which rotates the rotary body 10 can be reduced. For this reason, it is possible to accurately obtain an index indicating rotational accuracy, such as RRO, NRRO, and roundness. Further, by applying FFT to the data separated for each rotation period and made into a block without using a window function, the harmonics of the detection signal Sd having a fundamental frequency that is the reciprocal of one rotation period of the rotating body 10. The component can be obtained as a line spectrum having no spread. For this reason, RRO can be accurately obtained from the spectrum data obtained by performing the FFT, and when obtaining NRRO from the spectrum data, RRO can be accurately removed without affecting the other. it can.
[0031]
By the way, in the present embodiment, in the period separation process, separation of the collection data Da for each rotation period is performed based on zero point detection. At that time, the zero point detection includes moving average for each sampling value of the collection data Da. Smoothed displacement data consisting of For this reason, the collection data Da can be accurately separated for each rotation cycle. That is, if the collected data Da is used as it is for zero point detection, a plurality of zero points may apparently exist in the vicinity of one zero point (true zero point) due to noise, A / D conversion error, or the like. However, by using the smoothed displacement data, the zero point can be accurately detected in a one-to-one correspondence with the original zero point. For example, when a zero point is detected from the collected data Da shown in FIG. 6A, a rising zero point is detected in two sections between the sampling point pair (ta1, ta2) and the sampling point pair (ta2, ta3). There may be instances where a rising zero crossing sampling point pair that exists in appearance and includes a true rising zero point cannot be detected. For this reason, the collection data Da cannot be accurately separated for each rotation cycle. On the other hand, when the zero point is detected from the smoothed displacement data shown in FIG. 6B, it is detected that the rising zero point Zpa exists between the sampling point pair (ta1, ta2). No zero point is detected. For this reason, it becomes possible to isolate | separate the collection data Da correctly for every rotation period.
[0032]
As can be seen from the above-described procedure of the periodic separation process, it is not necessary to obtain the exact position of the zero point in order to separate the collection data Da, and it is sufficient to obtain the sampling point pair in which the zero point exists. . That is, if either one of the rising zero crossing sampling point pair or the falling zero crossing sampling point pair is detected, it is assumed that there is a zero point between the detected sampling point pairs, and the collected data Da is separated for each rotation period. can do. In the periodic separation process, for each sampled value of the sampled data Da, an average of the sampled values of the sampling points for a predetermined period centered on the sampling points of the sampled values is calculated. Must be set appropriately. In the present embodiment, the zero point corresponding to the original zero point is smoothed displacement data in consideration of the rotation period of the rotating body 10, the sampling frequency, the noise level, the frequency, and the like in this predetermined period. An appropriate predetermined period shall be selected so as to be detected from
[0033]
<Modification>
In the above embodiment, rate conversion is performed in order to make the number of data of each data block obtained by the period separation process the same, but instead, the following process may be performed. That is, the number of data with the highest appearance frequency is obtained from the number of data of each data block obtained by separating the collected data Da for each rotation period by the period separation process, and only the data block having the number of data with the highest appearance frequency is obtained. Extracted from the collected data Da. Then, an index indicating the rotation accuracy, such as NRRO, is calculated using the block data composed of the data blocks extracted in this way. Such a modification can also reduce the influence of rotation unevenness of the motor that rotates the rotating body 10 in the rotation accuracy measurement.
[Brief description of the drawings]
FIG. 1 is a configuration diagram schematically showing a configuration of a rotational accuracy measuring apparatus according to an embodiment of the present invention.
FIG. 2 is a plan view showing measurement of a rotating body by the rotation accuracy measuring device according to the embodiment.
FIG. 3 is a block diagram showing a functional configuration of a data processing device in the rotational accuracy measuring device according to the embodiment.
FIG. 4 is a flowchart showing a procedure of period separation processing in the embodiment.
FIG. 5 is a signal waveform diagram for explaining collected data in the embodiment.
FIG. 6 is a signal waveform diagram for explaining period separation processing in the embodiment.
FIG. 7 is a signal waveform diagram for explaining blocked data obtained by a periodic separation process on collected data in the embodiment.
FIG. 8 is a data configuration diagram for explaining rate conversion in the embodiment.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 ... Rotating body 11 ... Rotating shaft 20, 21 ... Displacement sensor 30 ... Data collection processing apparatus 31 ... CPU
32 ... Input interface unit 33 ... Memory 111 ... DC cut unit 112 ... Period separation unit 113 ... Rate conversion unit 114 ... Signal processing unit Sd ... Detection signal Da ... Collection data (displacement data)
Db ... Blocked data Dc ... Blocked data after rate conversion

Claims (1)

所定軸を中心として回転する回転体の回転精度を測定する回転精度測定装置であって、
前記回転体の半径方向の変位または前記所定軸方向の変位を示す信号の標本化値からなる時系列データを変位データとして取得するデータ取得手段と、
前記変位データを構成する各標本化値に対する移動平均を算出し、当該移動平均からなる時系列データである平滑化変位データを生成する移動平均手段と、
前記平滑化変位データの表す信号波形におけるゼロ点に基づき、前記データ取得手段によって取得された変位データを前記回転体の回転周期毎に分離してブロック化するブロック化手段と、
前記ブロック化手段によってブロック化された変位データにおける各データブロックを構成する標本化値の数が同一となるように、当該各データブロックに対してレート変換を行う変換手段と、
前記レート変換後の前記データブロックからなるブロック化データに基づき、前記回転精度を示す指標を算出する算出手段と
を備え、
前記ブロック化手段は、
前記平滑化変位データに基づき、前記回転体の変位を示す標本化値のサンプリング点における隣接サンプリング点対のうち、一方のサンプリング点における移動平均が正値で他方のサンプリング点における移動平均が負値である隣接サンプリング点対を検出するサンプリング点対検出手段と、
前記検出された隣接サンプリング点対の間に前記ゼロ点が存在するとみなして、前記データ取得手段によって取得された変位データを前記回転周期毎に分離する分離手段とを含むことを特徴とする回転精度測定装置。
A rotational accuracy measuring device that measures rotational accuracy of a rotating body that rotates about a predetermined axis,
Data acquisition means for acquiring, as displacement data, time-series data consisting of sampling values of signals indicating the radial displacement of the rotating body or the displacement in the predetermined axial direction;
A moving average means for calculating a moving average for each sampling value constituting the displacement data, and generating smoothed displacement data that is time-series data composed of the moving average;
Blocking means for separating and blocking the displacement data acquired by the data acquisition means for each rotation period of the rotating body based on the zero point in the signal waveform represented by the smoothed displacement data;
Conversion means for performing rate conversion on each data block so that the number of sampling values constituting each data block in the displacement data blocked by the blocking means is the same;
Calculation means for calculating an index indicating the rotation accuracy based on the blocked data composed of the data blocks after the rate conversion ;
The blocking means includes
Based on the smoothed displacement data, the moving average at one sampling point is a positive value and the moving average at the other sampling point is a negative value among the adjacent sampling point pairs at the sampling point of the sampling value indicating the displacement of the rotating body. Sampling point pair detection means for detecting adjacent sampling point pairs,
Rotational accuracy , comprising separation means for separating the displacement data acquired by the data acquisition means for each rotation period on the assumption that the zero point exists between the detected pair of adjacent sampling points measuring device.
JP2002036777A 2002-02-14 2002-02-14 Rotational accuracy measuring device Expired - Lifetime JP3824546B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002036777A JP3824546B2 (en) 2002-02-14 2002-02-14 Rotational accuracy measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002036777A JP3824546B2 (en) 2002-02-14 2002-02-14 Rotational accuracy measuring device

Publications (2)

Publication Number Publication Date
JP2003240536A JP2003240536A (en) 2003-08-27
JP3824546B2 true JP3824546B2 (en) 2006-09-20

Family

ID=27778567

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002036777A Expired - Lifetime JP3824546B2 (en) 2002-02-14 2002-02-14 Rotational accuracy measuring device

Country Status (1)

Country Link
JP (1) JP3824546B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004027800B4 (en) * 2004-06-08 2006-04-06 Fag Kugelfischer Ag & Co. Ohg Method and computer program for determining operating parameters in a rolling bearing and evaluable rolling bearing hereby
JP4699714B2 (en) * 2004-06-15 2011-06-15 富士フイルム株式会社 Eccentricity measuring apparatus and eccentricity measuring method
JP4667186B2 (en) * 2005-09-26 2011-04-06 学校法人慶應義塾 Rotational accuracy measurement method

Also Published As

Publication number Publication date
JP2003240536A (en) 2003-08-27

Similar Documents

Publication Publication Date Title
Shi et al. Bearing fault diagnosis under variable rotational speed via the joint application of windowed fractal dimension transform and generalized demodulation: A method free from prefiltering and resampling
US5501105A (en) Digital signal processing of encoder signals to detect resonances in rotating machines
CN110163190B (en) A kind of rolling bearing fault diagnosis method and device
US20080177485A1 (en) System and method for providing vibration detection in turbomachinery
US9791856B2 (en) Fault frequency set detection system and method
CN102353500B (en) Extraction method of unbalanced signal for dynamic balance measurement
JPH1026580A (en) Diagnosis method and apparatus for variable speed rotary machine equipment
CN102175412A (en) Rotor torsion oscillation testing device and testing method thereof
WO2007052550A1 (en) Method for measuring s/n ratio in eddy current scratch on inner surface of tube
CN104227501A (en) Testing and analyzing system for rotating errors of main shaft
BR112014007556B1 (en) method of determining stationary signals for the diagnosis of an electromechanical system
CN111458122A (en) Rotary machine fault diagnosis method based on matching enhancement time-frequency representation
Zhang et al. Research on the identification of asynchronous vibration parameters of rotating blades based on blade tip timing vibration measurement theory
JP3824546B2 (en) Rotational accuracy measuring device
EP2672240B1 (en) Campbell diagram displays and methods and systems for implementing same
CN110988676B (en) Mechanical equipment fault diagnosis method and device and storage medium
JP3765480B2 (en) Rotational accuracy measuring device
JP4096768B2 (en) Rotational accuracy measuring method and rotational accuracy measuring device
JP3997936B2 (en) Rotational accuracy measuring method and rotational accuracy measuring device
JP3997935B2 (en) Rotational accuracy measuring method and rotational accuracy measuring device
JPH02309212A (en) Method and apparatus for monitoring torsion of shaft
JP4096766B2 (en) Rotational accuracy measuring method and rotational accuracy measuring device
JP2000146771A (en) Simple abnormality diagnosis method for acceleration / deceleration type rotating machinery
CN111880095A (en) Induction motor fault analysis method and system based on the sum of squares of three-phase current signals
US10473535B2 (en) Methods and systems for non-contact magnetostrictive sensor runout compensation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040609

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20051017

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051101

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051213

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: 20060606

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060627

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090707

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100707

Year of fee payment: 4