JP3404129B2 - 光学式振動センサーによるうなり周波数測定方法 - Google Patents
光学式振動センサーによるうなり周波数測定方法Info
- Publication number
- JP3404129B2 JP3404129B2 JP14397994A JP14397994A JP3404129B2 JP 3404129 B2 JP3404129 B2 JP 3404129B2 JP 14397994 A JP14397994 A JP 14397994A JP 14397994 A JP14397994 A JP 14397994A JP 3404129 B2 JP3404129 B2 JP 3404129B2
- Authority
- JP
- Japan
- Prior art keywords
- vibration
- block
- time interval
- limit point
- beat frequency
- 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
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02K—DYNAMO-ELECTRIC MACHINES
- H02K11/00—Structural association of dynamo-electric machines with electric components or with devices for shielding, monitoring or protection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
- G01H1/14—Frequency
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Description
【０００１】
【産業上の利用分野】本発明は発電機、特に発電機の振
動をモニターするシステムに係わる。
動をモニターするシステムに係わる。
【０００２】
【従来の技術】１９９２年９月１５日付け米国特許第
５，１４６，７７６号“Method for Continuously Cali
brating an Optical Vibration Sensor"は発電機に片持
ばり式に取付けた格子を利用して光ファイバーまたは光
学式振動モニター（ＦＯＶＭ）を自動的に校正するシス
テムを開示している。格子は、ただ１つの駆動周波数
（即ち、１２０Ｈｚ）で振動するセンサーの振動振幅に
正比例する周波数で光線を遮断する。この特許は発電機
端コイルの振動振幅を測定することによって発電機の異
常状態を早い段階で検知する方法を教示している。これ
によって発電機の損傷を回避し、運転停止時間を極力短
くするように保守計画を立案することが可能になる。
５，１４６，７７６号“Method for Continuously Cali
brating an Optical Vibration Sensor"は発電機に片持
ばり式に取付けた格子を利用して光ファイバーまたは光
学式振動モニター（ＦＯＶＭ）を自動的に校正するシス
テムを開示している。格子は、ただ１つの駆動周波数
（即ち、１２０Ｈｚ）で振動するセンサーの振動振幅に
正比例する周波数で光線を遮断する。この特許は発電機
端コイルの振動振幅を測定することによって発電機の異
常状態を早い段階で検知する方法を教示している。これ
によって発電機の損傷を回避し、運転停止時間を極力短
くするように保守計画を立案することが可能になる。
【０００３】要約すると、このシステムは“極限点”に
おける信号の振幅を用いてうなり周波数を求める。ここ
にいう極限点とは平衡位置を中心とする振動の過程にお
ける格子の最遠移動点を意味する。極限点は格子が反対
方向へ移動する前に瞬間的に休止する点であるから、周
波数変調出力信号における最大波長（即ち、波形中でゼ
ロ交差点が最も離れる点）はこの極限点で生じる。
おける信号の振幅を用いてうなり周波数を求める。ここ
にいう極限点とは平衡位置を中心とする振動の過程にお
ける格子の最遠移動点を意味する。極限点は格子が反対
方向へ移動する前に瞬間的に休止する点であるから、周
波数変調出力信号における最大波長（即ち、波形中でゼ
ロ交差点が最も離れる点）はこの極限点で生じる。
【０００４】本発明は、うなり振幅がフォールドオーバ
ーを発生させるほど大きくなってうなり周波数を歪ませ
ると発生する問題に着目する。この問題は、極限点にお
ける信号の振幅がピーク信号値に近い場合は小さいうな
り振幅でも生じる。フォールドオーバーが発生すると、
極限点に基づいて求めるうなり信号が不正確になり、う
なり周波数の測定が極めて困難になって、事実上、セン
サーの共振周波数の正確な測定が不可能になる。しかも
共振周波数は温度及びセンサーの使用期間と共に変化す
る。従って、センサーの共振周波数が設計通りであると
想定することはできず、発電機の稼働中、共振周波数を
現場で測定しなければならない。
ーを発生させるほど大きくなってうなり周波数を歪ませ
ると発生する問題に着目する。この問題は、極限点にお
ける信号の振幅がピーク信号値に近い場合は小さいうな
り振幅でも生じる。フォールドオーバーが発生すると、
極限点に基づいて求めるうなり信号が不正確になり、う
なり周波数の測定が極めて困難になって、事実上、セン
サーの共振周波数の正確な測定が不可能になる。しかも
共振周波数は温度及びセンサーの使用期間と共に変化す
る。従って、センサーの共振周波数が設計通りであると
想定することはできず、発電機の稼働中、共振周波数を
現場で測定しなければならない。
【０００５】
【発明が解決しようとする課題】従って、本発明の主要
目的は、うなり振幅が大きい環境においてＦＯＶＭセン
サーによりうなり周波数を測定する方法を提供すること
にある。本発明の他の目的は、うなり振幅が小さくても
極限点における信号の振幅がピーク信号値に近い環境に
おいてＦＯＶＭセンサーによりうなり周波数を測定する
方法を提供することにある。
目的は、うなり振幅が大きい環境においてＦＯＶＭセン
サーによりうなり周波数を測定する方法を提供すること
にある。本発明の他の目的は、うなり振幅が小さくても
極限点における信号の振幅がピーク信号値に近い環境に
おいてＦＯＶＭセンサーによりうなり周波数を測定する
方法を提供することにある。
【０００６】この目的にかんがみ、本発明は、装置に取
付けられた振動感知システムにおけるうなり周波数を測
定する方法であって、装置の振動周波数及び振幅を示す
振動信号を発生させ；振動信号のうなり周波数に対応す
る点を含む波形部分の２つのゼロ交差間の時間を表わす
時間データを記憶し、処理することによって振動信号か
らうなり周波数を得るステップより成る、うなり周波数
測定方法に係わる。
付けられた振動感知システムにおけるうなり周波数を測
定する方法であって、装置の振動周波数及び振幅を示す
振動信号を発生させ；振動信号のうなり周波数に対応す
る点を含む波形部分の２つのゼロ交差間の時間を表わす
時間データを記憶し、処理することによって振動信号か
らうなり周波数を得るステップより成る、うなり周波数
測定方法に係わる。
【０００７】本発明の詳細は、添付図面に沿って以下に
述べる好ましい実施例の説明からさらに明らかになるで
あろう。
述べる好ましい実施例の説明からさらに明らかになるで
あろう。
【０００８】
【実施例】本発明は、極限点での格子の運動を測定する
ＦＯＶＭシステム用の方法を提供する。以下に述べる方
法を利用するシステムを図１に示した。発電機１０、光
学式振動センサー１２、及びコンピュータ１４が振動感
知システム１６を構成する。要約すれば、このシステム
を次のように説明することができる。即ち、光学式振動
センサー１２が発電機１０の端巻きコイル１７に直接取
付けてある。ローター領域からシステムへ伝達される１
２０Ｈｚの振動の有害な作用を回避するため、発電機の
大きい励磁機端部及びタービン端部の端巻きが一体化さ
れて半剛性のバスケットになっている。センサーは端巻
きコイルの振動をモニターして、破壊的レベルの振動が
存在するかまたは振動レベルが増大傾向にあると警報信
号を発する。その場合には、発電機の運転停止を予め決
定できるようになるまで、負荷を調整するか或いは冷却
ガス温度を変化させることによって振動を制御すればよ
い。
ＦＯＶＭシステム用の方法を提供する。以下に述べる方
法を利用するシステムを図１に示した。発電機１０、光
学式振動センサー１２、及びコンピュータ１４が振動感
知システム１６を構成する。要約すれば、このシステム
を次のように説明することができる。即ち、光学式振動
センサー１２が発電機１０の端巻きコイル１７に直接取
付けてある。ローター領域からシステムへ伝達される１
２０Ｈｚの振動の有害な作用を回避するため、発電機の
大きい励磁機端部及びタービン端部の端巻きが一体化さ
れて半剛性のバスケットになっている。センサーは端巻
きコイルの振動をモニターして、破壊的レベルの振動が
存在するかまたは振動レベルが増大傾向にあると警報信
号を発する。その場合には、発電機の運転停止を予め決
定できるようになるまで、負荷を調整するか或いは冷却
ガス温度を変化させることによって振動を制御すればよ
い。
【０００９】図２は光学式振動センサー１２をさらに詳
しく示している。光学式振動センサー１２は光ファイバ
ーケーブル１８を介して受光する。センサーはハウジン
グ２０及び光／デジタルまたは光電変換手段２２を含
む。ハウジング２０は内部リード２４及び格子組立体２
６を含む。内部リード２４及び格子組立体２６は固有共
振周波数が１２０Ｈｚ以上となるように設計されてい
る。６０Ｈｚ発電機に用いる場合、共振周波数は約１３
２Ｈｚであることが好ましい。
しく示している。光学式振動センサー１２は光ファイバ
ーケーブル１８を介して受光する。センサーはハウジン
グ２０及び光／デジタルまたは光電変換手段２２を含
む。ハウジング２０は内部リード２４及び格子組立体２
６を含む。内部リード２４及び格子組立体２６は固有共
振周波数が１２０Ｈｚ以上となるように設計されてい
る。６０Ｈｚ発電機に用いる場合、共振周波数は約１３
２Ｈｚであることが好ましい。
【００１０】以下の説明は発電機が周波数６０Ｈｚの電
力を供給する場合を想定しているが、５０Ｈｚの発電機
の場合でも原理は同じである。
力を供給する場合を想定しているが、５０Ｈｚの発電機
の場合でも原理は同じである。
【００１１】内部リードが振動すると、格子組立体２６
が上下動して光パルスを発生させる。所与の時間に亘っ
て発生する光パルスの数は測定すべき１２０Ｈｚの（ヨ
ーロッパでは１００Ｈｚ）振動の振幅に比例する。格子
組立体２６は１０ミルの等間隔に配列された格子開口を
有する。即ち、所与の時間に亘って発生する光パルスの
数はセンサーの共振周波数及び格子の平衡位置からの振
動距離に応じて異なる。光パルスはハウジング２０から
光ファイバーケーブル１８を介して光電変換手段２２へ
出力される。光電変換手段２２は公知の方法で光パルス
をデジタル信号に変換する。例えば、フォトダイオード
を利用して光パルスを電気信号に変換し、次いでこれを
デジタル周波数出力信号に変換することができる。出力
信号波形は周波数変調正弦波の形態を取る。さらにこの
信号は１２０Ｈｚの励振信号をセンサーの共振周波数と
混合することで僅かに周波数変調される。
が上下動して光パルスを発生させる。所与の時間に亘っ
て発生する光パルスの数は測定すべき１２０Ｈｚの（ヨ
ーロッパでは１００Ｈｚ）振動の振幅に比例する。格子
組立体２６は１０ミルの等間隔に配列された格子開口を
有する。即ち、所与の時間に亘って発生する光パルスの
数はセンサーの共振周波数及び格子の平衡位置からの振
動距離に応じて異なる。光パルスはハウジング２０から
光ファイバーケーブル１８を介して光電変換手段２２へ
出力される。光電変換手段２２は公知の方法で光パルス
をデジタル信号に変換する。例えば、フォトダイオード
を利用して光パルスを電気信号に変換し、次いでこれを
デジタル周波数出力信号に変換することができる。出力
信号波形は周波数変調正弦波の形態を取る。さらにこの
信号は１２０Ｈｚの励振信号をセンサーの共振周波数と
混合することで僅かに周波数変調される。
【００１２】本発明のシステムは、うなり信号のピーク
と三角関数ｓｉｎ（２πｆBｔ）との間で曲線のあては
めを行うことによってうなり周波数ｆBを測定する。次
いでこのうなり周波数を利用してシステムを校正する。
具体的には、本発明のシステムは下式の増幅率を計算す
る。
と三角関数ｓｉｎ（２πｆBｔ）との間で曲線のあては
めを行うことによってうなり周波数ｆBを測定する。次
いでこのうなり周波数を利用してシステムを校正する。
具体的には、本発明のシステムは下式の増幅率を計算す
る。
【００１３】
Ｍ0 ＝（１２０／ｆ0 ）2 ／（１−（１２０／ｆ0 ）2 ）
ただし、センサーの共振周波数は
ｆ0 ＝１２０Ｈｚ ＋ｆB
によって与えられる。即ち、光学式振動センサーの共振
周波数ｆ0 が増幅率Ｍ0を決定する。１２０Ｈｚの振動
によって発電機が変位する実位置を求めるためには測定
振幅値（即ち、光パルス信号によって測定される振幅
値）を増幅率によって除算しなければならない。なお、
Ｍ0 を求める方程式は光パルス周波数と格子振幅との相
関関係から成立する式であり、格子−リードの幾何形状
が出願人が採用するものである場合、下記の式で表わす
ことができる： 振動の振幅＝ｆLP × １ミル／１８０Ｈｚ ただし、ｆLPは光パルスの周波数（Ｈｚ）であり、この
方程式は格子間隔が１０ミルの格子組立体にもあてはま
る。
周波数ｆ0 が増幅率Ｍ0を決定する。１２０Ｈｚの振動
によって発電機が変位する実位置を求めるためには測定
振幅値（即ち、光パルス信号によって測定される振幅
値）を増幅率によって除算しなければならない。なお、
Ｍ0 を求める方程式は光パルス周波数と格子振幅との相
関関係から成立する式であり、格子−リードの幾何形状
が出願人が採用するものである場合、下記の式で表わす
ことができる： 振動の振幅＝ｆLP × １ミル／１８０Ｈｚ ただし、ｆLPは光パルスの周波数（Ｈｚ）であり、この
方程式は格子間隔が１０ミルの格子組立体にもあてはま
る。
【００１４】要約すると、本発明のシステムは“極限
点”における信号の振幅を利用することによってうなり
周波数を求める。極限点とは平衡位置を起点として格子
が最も遠くまで移動した点である。周波数変調出力信号
の最大波長（即ち、波形中でゼロ交差点が最も離れる
点）はこの極限点で生じる。即ち、極限点は格子が移動
方向を反転させる前に瞬間的に休止する点を意味するか
らである。
点”における信号の振幅を利用することによってうなり
周波数を求める。極限点とは平衡位置を起点として格子
が最も遠くまで移動した点である。周波数変調出力信号
の最大波長（即ち、波形中でゼロ交差点が最も離れる
点）はこの極限点で生じる。即ち、極限点は格子が移動
方向を反転させる前に瞬間的に休止する点を意味するか
らである。
【００１５】小さいフォールドオーバー現象を図４に図
示した。極めて大きいフォールドオーバーが起こること
も珍しくないが、このような大きいフォールドオーバー
の波形の極限点を視認することは困難であるから図示し
なかった。上述したとおり、１２０Ｈｚの振動に起因す
る発電機の変位の実位置を測定するには増幅率Ｍ0 を求
めねばならない。増幅率を求めるにはセンサーの共振周
波数（ｆ0 ＋ｆB ）を正確に測定しなければならない。
示した。極めて大きいフォールドオーバーが起こること
も珍しくないが、このような大きいフォールドオーバー
の波形の極限点を視認することは困難であるから図示し
なかった。上述したとおり、１２０Ｈｚの振動に起因す
る発電機の変位の実位置を測定するには増幅率Ｍ0 を求
めねばならない。増幅率を求めるにはセンサーの共振周
波数（ｆ0 ＋ｆB ）を正確に測定しなければならない。
【００１６】図５はセンサー波形（２サイクルに亘る電
圧振幅）の一例を示す。上方極限点及び下方極限点に対
応する極限点での時間インターバルも示した。上方極限
点を“ｘ”（極限点の位置）または線“Ａ”（極限点に
おける時間インターバル）で表わし、下方極限点を
“０”または線“Ｂ”で表わしている。なお、図５から
明らかなように、上方極限点を距離Ｘ（即ち、０からｘ
までの距離）または時間インターバルＡで表わすことも
できる。同様に、下方極限点を距離０または時間インタ
ーバルＢで表わすことができる。このように時間インタ
ーバルで表現することにしたのは下記の理由による。も
し上方及び下方極限点を位置で表わすとすれば、ＸがＰ
（グラフの頂部）に達し、反射する時、ＸはＡよりも大
きくなれないから、フォールドオーバーは起こらない。
もし上方及び下方極限点を時間で表わせば、Ｘ＝０で、
Ａが大きい値から小さい値へ、あるいはその逆の方向に
急変するとフォールドオーバーが起こる。この２つの表
現方式は上方及び下方極限点を定義する別々の方法では
あるが、ほとんど等価である。ただし、位置表現は線形
に挙動しないのに対して、時間表現は線形に挙動する。
従って、もしうなり振幅がゆっくりと増大すれば、Ｘが
増大してＰに達し、次いで縮小する。Ｘのこの変化はう
なり振幅に比例しない。即ち、うなり振幅が増大してい
る間もＸがＰに“停滞”したままになり易い。これは光
の通過位置を中心とする格子の位置を感知する光学式セ
ンサーが全光量を受光することによって起こる状態であ
る。下方極限点では２つの隣接スリット間の遮光部分と
の関連でこの現象が起こる。即ち、遮光部分は移動でき
るが格子が小さく振動する時には光を遮断する。
圧振幅）の一例を示す。上方極限点及び下方極限点に対
応する極限点での時間インターバルも示した。上方極限
点を“ｘ”（極限点の位置）または線“Ａ”（極限点に
おける時間インターバル）で表わし、下方極限点を
“０”または線“Ｂ”で表わしている。なお、図５から
明らかなように、上方極限点を距離Ｘ（即ち、０からｘ
までの距離）または時間インターバルＡで表わすことも
できる。同様に、下方極限点を距離０または時間インタ
ーバルＢで表わすことができる。このように時間インタ
ーバルで表現することにしたのは下記の理由による。も
し上方及び下方極限点を位置で表わすとすれば、ＸがＰ
（グラフの頂部）に達し、反射する時、ＸはＡよりも大
きくなれないから、フォールドオーバーは起こらない。
もし上方及び下方極限点を時間で表わせば、Ｘ＝０で、
Ａが大きい値から小さい値へ、あるいはその逆の方向に
急変するとフォールドオーバーが起こる。この２つの表
現方式は上方及び下方極限点を定義する別々の方法では
あるが、ほとんど等価である。ただし、位置表現は線形
に挙動しないのに対して、時間表現は線形に挙動する。
従って、もしうなり振幅がゆっくりと増大すれば、Ｘが
増大してＰに達し、次いで縮小する。Ｘのこの変化はう
なり振幅に比例しない。即ち、うなり振幅が増大してい
る間もＸがＰに“停滞”したままになり易い。これは光
の通過位置を中心とする格子の位置を感知する光学式セ
ンサーが全光量を受光することによって起こる状態であ
る。下方極限点では２つの隣接スリット間の遮光部分と
の関連でこの現象が起こる。即ち、遮光部分は移動でき
るが格子が小さく振動する時には光を遮断する。
【００１７】従って、うなりが格子を上方極限点または
下方極限点における平均位置を中心に振動させると、こ
の時間インターバルは極限点における信号振幅の変化と
同様に変化する。うなり振幅または振動振幅によりセン
サー信号（即ち、波形）が０振幅軸と交差して非連続的
な態様で新しい極限点の時間インターバルが画定される
と、フォールドオーバーが起こる。これと同様の不連続
性は、米国特許第５，１４６，７７６号に開示されてい
る方法では、信号の上下ピークにおいてフォールドオー
バーが起こった場合に現われる。
下方極限点における平均位置を中心に振動させると、こ
の時間インターバルは極限点における信号振幅の変化と
同様に変化する。うなり振幅または振動振幅によりセン
サー信号（即ち、波形）が０振幅軸と交差して非連続的
な態様で新しい極限点の時間インターバルが画定される
と、フォールドオーバーが起こる。これと同様の不連続
性は、米国特許第５，１４６，７７６号に開示されてい
る方法では、信号の上下ピークにおいてフォールドオー
バーが起こった場合に現われる。
【００１８】図６は、上方極限点の時間インターバルが
ＡからＢへ移行する過程で極限点の時間インターバルに
現われる非連続的な変化を示す。このような不連続性は
センサーのうなり周波数を測定する上で問題を生じさせ
る。
ＡからＢへ移行する過程で極限点の時間インターバルに
現われる非連続的な変化を示す。このような不連続性は
センサーのうなり周波数を測定する上で問題を生じさせ
る。
【００１９】本発明の好ましい実施態様は、信号のゼロ
交差時だけをサンプリングして処理するが、完全なセン
サー信号をディジタル化するものではない。従って、デ
ータ収集条件が著しく軽減され、しかも極めて正確な測
定が可能になる。この測定には既存の羽根振動モニター
（ＢＶＭ）３２ＭＨｚカード（ゼロ交差カード）を使用
すればよい。ＢＶＭカードはゼロ交差時を３０／１０9
ｓｅｃ以内の精度で測定することができ、従って、うな
り周波数を極めて正確に測定できる。データファイルの
典型的な長さは１２，０００個の記述項に相当する。同
じ精度を達成するのに、上記特許に開示されている直接
ディジタル化法では３２０，０００個の記述項に相当す
る長さのデータファイルが必要である。このデータの大
部分は捨てられることになるが、そのためには多大の時
間と高価なハードウエアが必要である。従って、本発明
は必要データを最小限にとどめてハードウエアのコスト
及びコンピュータによる処理時間を節減する。
交差時だけをサンプリングして処理するが、完全なセン
サー信号をディジタル化するものではない。従って、デ
ータ収集条件が著しく軽減され、しかも極めて正確な測
定が可能になる。この測定には既存の羽根振動モニター
（ＢＶＭ）３２ＭＨｚカード（ゼロ交差カード）を使用
すればよい。ＢＶＭカードはゼロ交差時を３０／１０9
ｓｅｃ以内の精度で測定することができ、従って、うな
り周波数を極めて正確に測定できる。データファイルの
典型的な長さは１２，０００個の記述項に相当する。同
じ精度を達成するのに、上記特許に開示されている直接
ディジタル化法では３２０，０００個の記述項に相当す
る長さのデータファイルが必要である。このデータの大
部分は捨てられることになるが、そのためには多大の時
間と高価なハードウエアが必要である。従って、本発明
は必要データを最小限にとどめてハードウエアのコスト
及びコンピュータによる処理時間を節減する。
【００２０】以下に述べるアルゴリズムは、ＦＯＶＭセ
ンサーにおいて、高いうなり周波数振幅によりそれまで
よりも１本または数本多い格子スリットが光を通過させ
ることにより、うなり周波数を求めるものである。この
現象がフォールドオーバーを発生させる。センサーサイ
クルごとにセンサー信号中に２つのタイプの極限点、即
ち、上方極限点及び下方極限点が発生する。極限点の時
間インターバルは各極限点の両側における信号のゼロ交
差によって与えられる。上方極限点の時間インターバル
は格子が瞬間的にその振動経路の上端で休止する時に現
われ、同様に、下方極限点の時間インターバルは格子が
瞬間的に振動経路の下端で休止する時に現われる。
ンサーにおいて、高いうなり周波数振幅によりそれまで
よりも１本または数本多い格子スリットが光を通過させ
ることにより、うなり周波数を求めるものである。この
現象がフォールドオーバーを発生させる。センサーサイ
クルごとにセンサー信号中に２つのタイプの極限点、即
ち、上方極限点及び下方極限点が発生する。極限点の時
間インターバルは各極限点の両側における信号のゼロ交
差によって与えられる。上方極限点の時間インターバル
は格子が瞬間的にその振動経路の上端で休止する時に現
われ、同様に、下方極限点の時間インターバルは格子が
瞬間的に振動経路の下端で休止する時に現われる。
【００２１】うなり振幅が小さい時、うなりにより極限
点の時間インターバルに僅かな変調が起こる。下方極限
点の時間インターバルとセンサーサイクルとの関係を作
図すると、うなりと同じ周波数の正弦波が得られる。う
なり振幅が大きくなって（または小さくなって）２つの
ゼロ交差が追加されると（または失われると）、このう
なり周波数抽出手順は予期しない小さい（または大き
い）極限点の時間インターバルが突如現われてうなり信
号に顕著な不連続性が生ずることで混乱させられる。
点の時間インターバルに僅かな変調が起こる。下方極限
点の時間インターバルとセンサーサイクルとの関係を作
図すると、うなりと同じ周波数の正弦波が得られる。う
なり振幅が大きくなって（または小さくなって）２つの
ゼロ交差が追加されると（または失われると）、このう
なり周波数抽出手順は予期しない小さい（または大き
い）極限点の時間インターバルが突如現われてうなり信
号に顕著な不連続性が生ずることで混乱させられる。
【００２２】図７は４５０個のセンサーサイクルの各セ
ンサーサイクルごとに一度現われる下方極限点の時間イ
ンターバル”ＤＥＬＴＡ（）”を示す。この極めて混沌
とした非連続的な波形は調和のとれたうなり信号とはい
えない。図８はこのデータのフーリエ変換を示す。１
２．６Ｈｚで予想されるうなり周波数は、極限点におけ
るフォールドオーバー（追加されたゼロ交差）により発
生されるノイズに混入して検出できなくなる。
ンサーサイクルごとに一度現われる下方極限点の時間イ
ンターバル”ＤＥＬＴＡ（）”を示す。この極めて混沌
とした非連続的な波形は調和のとれたうなり信号とはい
えない。図８はこのデータのフーリエ変換を示す。１
２．６Ｈｚで予想されるうなり周波数は、極限点におけ
るフォールドオーバー（追加されたゼロ交差）により発
生されるノイズに混入して検出できなくなる。
【００２３】本発明では、フォールドオーバーが発生す
るとうなり振幅の増大または縮小を補償するため極限点
の時間インターバルＤＥＬＴＡ（）を隣接する適当個数
の時間インターバルに亘って加算することにより拡大す
る。下記アルゴリズムはフォールドオーバーが起こると
これを検出し、時間インターバルＤＥＬＴＡ（）の両側
に適当個数の時間インターバルを加算することによっ
て、測定されるこの時間インターバルを拡大する。アル
ゴリズムは例えばＱＵＩＣＫ ＢＡＳＩＳ（by Microso
ft）プログラミング語で記載されているコンピュータソ
フトウエアによって行うことができる。最小のＤＥＬＴ
Ａ（）は中心にある１個の時間インターバルであるが、
この中心時間インターバルの両側に１個、２個、３個ま
たは４個以上の時間インターバルを加算したものであっ
てもよい。極限点の時間インターバルの両側にいくつの
時間インターバルを加算すべきかの判断は、下記の事情
によって複雑になる： − 隣接する極限点の間で２つのゼロ交差が追加された
（または失われた）ことが分かっていても、波形に関す
る具体的な知識がなければ、上下いずれの極限点に起因
するものかは知るすべがない。本発明はＦＯＶＭが信号
のゼロ交差だけをサンプリングし、処理するものであ
り、完全なセンサー信号のディジタル化は行わない。こ
れによってデータ収集条件が著しく軽減され、従来より
もはるかに正確な測定が可能になる。このためには３２
ＭＨｚゼロ交差ボードを使用すればよい。
るとうなり振幅の増大または縮小を補償するため極限点
の時間インターバルＤＥＬＴＡ（）を隣接する適当個数
の時間インターバルに亘って加算することにより拡大す
る。下記アルゴリズムはフォールドオーバーが起こると
これを検出し、時間インターバルＤＥＬＴＡ（）の両側
に適当個数の時間インターバルを加算することによっ
て、測定されるこの時間インターバルを拡大する。アル
ゴリズムは例えばＱＵＩＣＫ ＢＡＳＩＳ（by Microso
ft）プログラミング語で記載されているコンピュータソ
フトウエアによって行うことができる。最小のＤＥＬＴ
Ａ（）は中心にある１個の時間インターバルであるが、
この中心時間インターバルの両側に１個、２個、３個ま
たは４個以上の時間インターバルを加算したものであっ
てもよい。極限点の時間インターバルの両側にいくつの
時間インターバルを加算すべきかの判断は、下記の事情
によって複雑になる： − 隣接する極限点の間で２つのゼロ交差が追加された
（または失われた）ことが分かっていても、波形に関す
る具体的な知識がなければ、上下いずれの極限点に起因
するものかは知るすべがない。本発明はＦＯＶＭが信号
のゼロ交差だけをサンプリングし、処理するものであ
り、完全なセンサー信号のディジタル化は行わない。こ
れによってデータ収集条件が著しく軽減され、従来より
もはるかに正確な測定が可能になる。このためには３２
ＭＨｚゼロ交差ボードを使用すればよい。
【００２４】− 時間インターバルの加算には適正な開
始パラメーターを知る必要がある。各極限点の状態は唯
１つのパラメーター”Ｊ％（）”によって判定すること
ができる。Ｊ％（）は問題の極限点と先行の極限点との
間の時間インターバルの個数（ゼロ交差回数＋１）であ
る。加算パラメーター”Ｓ１（）”は各極限点において
加算すべき時間インターバルの個数を規定する。加算す
べき時間インターバルの個数は２＊Ｓ１（）＋１に等し
い。加算パラメーターの開始値（即ち、Ｓ１（１）及び
Ｓ１（２）の値）が不適正であれば憂慮すべき結果を招
くことになる。
始パラメーターを知る必要がある。各極限点の状態は唯
１つのパラメーター”Ｊ％（）”によって判定すること
ができる。Ｊ％（）は問題の極限点と先行の極限点との
間の時間インターバルの個数（ゼロ交差回数＋１）であ
る。加算パラメーター”Ｓ１（）”は各極限点において
加算すべき時間インターバルの個数を規定する。加算す
べき時間インターバルの個数は２＊Ｓ１（）＋１に等し
い。加算パラメーターの開始値（即ち、Ｓ１（１）及び
Ｓ１（２）の値）が不適正であれば憂慮すべき結果を招
くことになる。
【００２５】奇数番目のＤＥＬＴＡ（）は各上方極限点
に対応する適正な拡大時間インターバルであり、偶数番
目のＤＥＬＴＡ（）は各下方極限点に対応する適正拡大
時間インターバルである。（なお、“ＤＥＬＴＡ（）”
はＤＥＬＴＡ（１）、ＤＥＬＴＡ（２）、、、ＤＥＬＴ
Ａ（Ｎ）を表わし、図７に示す例ではＮ＝４５０であ
る。）図９は図７で使用されているのと同じデータに関
連して下方極限点に対応するＤＥＬＴＡ（）を作図した
ものである。ただし、図７と異なり、極限点の時間イン
ターバルＤＥＬＴＡ（）は適正な個数の隣接時間インタ
ーバルを加算することによって拡大されている。この図
ではうなりが明瞭かつ鮮明に見える。図１０は下方極限
点うなり信号のフーリエ変換を示す。１２．６Ｈｚのう
なり周波数がはっきり識別される。従って、センサーの
固有周波数は１２０Ｈｚ＋１２．６Ｈｚ＝１３２．６Ｈ
ｚであり、センサーの補正率は容易に見出すことができ
る。上限に関して行われる同じ手順で全く同じ結果が得
られる（これは例えばチェックとして利用できる）。
に対応する適正な拡大時間インターバルであり、偶数番
目のＤＥＬＴＡ（）は各下方極限点に対応する適正拡大
時間インターバルである。（なお、“ＤＥＬＴＡ（）”
はＤＥＬＴＡ（１）、ＤＥＬＴＡ（２）、、、ＤＥＬＴ
Ａ（Ｎ）を表わし、図７に示す例ではＮ＝４５０であ
る。）図９は図７で使用されているのと同じデータに関
連して下方極限点に対応するＤＥＬＴＡ（）を作図した
ものである。ただし、図７と異なり、極限点の時間イン
ターバルＤＥＬＴＡ（）は適正な個数の隣接時間インタ
ーバルを加算することによって拡大されている。この図
ではうなりが明瞭かつ鮮明に見える。図１０は下方極限
点うなり信号のフーリエ変換を示す。１２．６Ｈｚのう
なり周波数がはっきり識別される。従って、センサーの
固有周波数は１２０Ｈｚ＋１２．６Ｈｚ＝１３２．６Ｈ
ｚであり、センサーの補正率は容易に見出すことができ
る。上限に関して行われる同じ手順で全く同じ結果が得
られる（これは例えばチェックとして利用できる）。
【００２６】コンピュータによって行われる本発明のア
ルゴリズムの好ましい実施例を図１１〜図１４のフロー
チャートに沿って説明する。ゼロ交差Ｉ／ＯカードはＦ
ＯＶＭセンサー電圧が０ボルトレベルと交差した絶対時
間（クロックカウント）を表わす１次元データアレイＺ
（）（例えば約１０，０００個の記述項を有するデータ
アレイ）を発生させる。従って、Ｚ（）はゼロ交差回数
をクロックカウントとして表わす一本調子に増大する数
列である。本発明の一実施例ではＺ（）アレイに１０，
０００個の記憶場所がある。
ルゴリズムの好ましい実施例を図１１〜図１４のフロー
チャートに沿って説明する。ゼロ交差Ｉ／ＯカードはＦ
ＯＶＭセンサー電圧が０ボルトレベルと交差した絶対時
間（クロックカウント）を表わす１次元データアレイＺ
（）（例えば約１０，０００個の記述項を有するデータ
アレイ）を発生させる。従って、Ｚ（）はゼロ交差回数
をクロックカウントとして表わす一本調子に増大する数
列である。本発明の一実施例ではＺ（）アレイに１０，
０００個の記憶場所がある。
【００２７】フローチャートのブロックは下記のように
ステップ１−５に対応する： ステップ−１＝ブロック１００−１０６ ステップ−２＝ブロック１０８−１３８ ステップ−３＝ブロック１４０−１４４ ステップ−４＝ブロック１４６−１７４ ステップ−５＝ブロック１７６−２００ 要約すると、ステップ１−５は次のような機能を行う：
ステップ−１は隣接ゼロ交差間の時間インターバル（ク
ロックカウント）を計算する。それぞれのクロックカウ
ントはＸ（）アレイに記憶される。ステップ−２は
Ｘ（）アレイ中の極限点の時間インターバルの位置と、
隣接する極限点間に現われる時間インターバルの個数
（Ｊ％（）で表わされる）と、２つの隣接する極限点間
に現われる時間インターバルの最小個数を見出す。この
最小個数を変数ＭＩＮで表わされる。うなり振幅が増大
すれば極限点の時間インターバルの両側において加算す
べき時間インターバルの個数が増える。ステップ−３は
一連の極限点のうちから、うなり振幅が最も小さく、
（ＤＥＬＴＡ（）で表わされる）補正時間インターバル
データを形成するのにあと１個の時間インターバルだけ
でよい、即ち、フォールドオーバーが存在しない第１の
点（ＳＴＡＲＴ）を見出だす。ステップ−４は点ＳＴＡ
ＲＴからＫまでに亘って加算すべき各極限点の時間イン
ターバルを中心とする時間インターバルの個数を求める
順方向チェーンから成る。なお、（フォールドオーバー
が存在しないから）Ｓ１（ＳＴＡＲＴ）＝Ｓ１（ＳＴＡ
ＲＴ＋１）＝０である。ステップ−５は加算すべき時間
インターバルの個数がこのタイプの最後の極限点、即
ち、上方極限点または下方極限点のあとに現われる極限
点間のクロックカウント数Ｊ％（）に基づいて変化すべ
きかどうかを判断するＫ−１から０までの逆方向チェー
ンから成る。
ステップ１−５に対応する： ステップ−１＝ブロック１００−１０６ ステップ−２＝ブロック１０８−１３８ ステップ−３＝ブロック１４０−１４４ ステップ−４＝ブロック１４６−１７４ ステップ−５＝ブロック１７６−２００ 要約すると、ステップ１−５は次のような機能を行う：
ステップ−１は隣接ゼロ交差間の時間インターバル（ク
ロックカウント）を計算する。それぞれのクロックカウ
ントはＸ（）アレイに記憶される。ステップ−２は
Ｘ（）アレイ中の極限点の時間インターバルの位置と、
隣接する極限点間に現われる時間インターバルの個数
（Ｊ％（）で表わされる）と、２つの隣接する極限点間
に現われる時間インターバルの最小個数を見出す。この
最小個数を変数ＭＩＮで表わされる。うなり振幅が増大
すれば極限点の時間インターバルの両側において加算す
べき時間インターバルの個数が増える。ステップ−３は
一連の極限点のうちから、うなり振幅が最も小さく、
（ＤＥＬＴＡ（）で表わされる）補正時間インターバル
データを形成するのにあと１個の時間インターバルだけ
でよい、即ち、フォールドオーバーが存在しない第１の
点（ＳＴＡＲＴ）を見出だす。ステップ−４は点ＳＴＡ
ＲＴからＫまでに亘って加算すべき各極限点の時間イン
ターバルを中心とする時間インターバルの個数を求める
順方向チェーンから成る。なお、（フォールドオーバー
が存在しないから）Ｓ１（ＳＴＡＲＴ）＝Ｓ１（ＳＴＡ
ＲＴ＋１）＝０である。ステップ−５は加算すべき時間
インターバルの個数がこのタイプの最後の極限点、即
ち、上方極限点または下方極限点のあとに現われる極限
点間のクロックカウント数Ｊ％（）に基づいて変化すべ
きかどうかを判断するＫ−１から０までの逆方向チェー
ンから成る。
【００２８】ステップ−１
図１１のブロック１００において、変数“Ｉ“が１にセ
ットされる。ブロック１０２において、Ｘ（Ｉ）がＺ
（Ｉ＋１）−Ｚ（Ｉ）にセットされる。ブロック１０４
において、Ｉが１０，０００と比較される。もしＩが１
０，０００に等しくなければ、プログラムはブロック１
０６へ分岐し；さもなければブロック１０８へ進む。ブ
ロック１０６において、ＩがＩ＋１にセットされる。
ットされる。ブロック１０２において、Ｘ（Ｉ）がＺ
（Ｉ＋１）−Ｚ（Ｉ）にセットされる。ブロック１０４
において、Ｉが１０，０００と比較される。もしＩが１
０，０００に等しくなければ、プログラムはブロック１
０６へ分岐し；さもなければブロック１０８へ進む。ブ
ロック１０６において、ＩがＩ＋１にセットされる。
【００２９】ステップ−２
ブロック１０８において、Ｉが５にセットされ、ＴＯＴ
が０にセットされ、Ｋが１にセットされ、ＭＩＮが１０
０にセットされる。ブロック１１０において、“ＳＫＩ
Ｐ１”がＸ（Ｉ）＋Ｘ（Ｉ＋１）にセットされ、ＳＫＩ
Ｐ２がＸ（Ｉ＋３）＋Ｘ（Ｉ＋４）にセットされる。ブ
ロック１１２において、ＳＫＩＰ１がＳＫＩＰ２と比較
される。もしＳＫＩＰ１がＳＫＩＰ２よりも小さけれ
ば、プログラムはブロック１１４に進み；さもなければ
ブロック１１６に進む。ブロック１１４においてＩがＩ
＋１にセットされ、プログラムはブロック１１０へ戻
る。即ち、ブロック１１０−１１４は時間インターバル
の長さが増大しつつあることを見極める。従って、極限
点間に現われる時間インターバルアレイＸ（）中に対称
な点は検出されない。
が０にセットされ、Ｋが１にセットされ、ＭＩＮが１０
０にセットされる。ブロック１１０において、“ＳＫＩ
Ｐ１”がＸ（Ｉ）＋Ｘ（Ｉ＋１）にセットされ、ＳＫＩ
Ｐ２がＸ（Ｉ＋３）＋Ｘ（Ｉ＋４）にセットされる。ブ
ロック１１２において、ＳＫＩＰ１がＳＫＩＰ２と比較
される。もしＳＫＩＰ１がＳＫＩＰ２よりも小さけれ
ば、プログラムはブロック１１４に進み；さもなければ
ブロック１１６に進む。ブロック１１４においてＩがＩ
＋１にセットされ、プログラムはブロック１１０へ戻
る。即ち、ブロック１１０−１１４は時間インターバル
の長さが増大しつつあることを見極める。従って、極限
点間に現われる時間インターバルアレイＸ（）中に対称
な点は検出されない。
【００３０】ブロック１１６において、Ｊが１にセット
される。ブロック１１８において、Ｔ１が（Ｘ（Ｉ＋
Ｊ）−Ｘ（Ｉ−Ｊ）÷（Ｘ（Ｉ＋Ｊ）＋Ｘ（Ｉ−Ｊ））
の絶対値にセットされる。ブロック１２０においてＴ１
が好ましい実施例では０，１に等しいトリガー変数“Ｔ
Ｒ”と比較される。もしＴ１がＴＲよりも小さければプ
ログラムはブロック１２２へ分岐し；さもなければブロ
ック１２４へ進む。ブロック１２２において、ＴＯＴが
ＴＯＴ＋１にセットされる。ブロック１２４において、
ＪがＪ＋１にセットされる。ブロック１２６においてＪ
が数５と比較される。もしＪが５以下なら、プログラム
はブロック１１８へ戻り；さもなければブロック１２８
へ進む。ブロック１２８において、ＴＯＴが数３と比較
される。もしＴＯＴが３以下ならプログラムはブロック
１１４（図１１）に戻り；さもなければブロック１３０
へ進む。ブロック１３０において、変数ＥＸ（Ｋ）がＩ
にセットされる。ブロック１３２において、Ｊ％（Ｋ）
がＥＸＴ（Ｋ）−ＥＸＴ（Ｋ−１）−１にセットされ
る。ブロック１３４においてＪ％（Ｋ）がＭＩＮと比較
される。もしＪ％（Ｋ）がＭＩＮよりも小さければ、プ
ログラムはブロック１３６へ分岐し；さもなければブロ
ック１３８へ進む。ブロック１３６において、ＭＩＮが
Ｊ％（Ｋ）にセットされる。ブロック１３８において、
インデックスＩが（Ｉ＋Ｊ％（Ｋ）／２＋２）の整数に
セットされる。もしＩが１０，０００以下ならプログラ
ムはブロック１１６（図１１）に戻る（Ｉが１０，００
０以下であるかどうかを判定するテストは図示されてい
ない）。
される。ブロック１１８において、Ｔ１が（Ｘ（Ｉ＋
Ｊ）−Ｘ（Ｉ−Ｊ）÷（Ｘ（Ｉ＋Ｊ）＋Ｘ（Ｉ−Ｊ））
の絶対値にセットされる。ブロック１２０においてＴ１
が好ましい実施例では０，１に等しいトリガー変数“Ｔ
Ｒ”と比較される。もしＴ１がＴＲよりも小さければプ
ログラムはブロック１２２へ分岐し；さもなければブロ
ック１２４へ進む。ブロック１２２において、ＴＯＴが
ＴＯＴ＋１にセットされる。ブロック１２４において、
ＪがＪ＋１にセットされる。ブロック１２６においてＪ
が数５と比較される。もしＪが５以下なら、プログラム
はブロック１１８へ戻り；さもなければブロック１２８
へ進む。ブロック１２８において、ＴＯＴが数３と比較
される。もしＴＯＴが３以下ならプログラムはブロック
１１４（図１１）に戻り；さもなければブロック１３０
へ進む。ブロック１３０において、変数ＥＸ（Ｋ）がＩ
にセットされる。ブロック１３２において、Ｊ％（Ｋ）
がＥＸＴ（Ｋ）−ＥＸＴ（Ｋ−１）−１にセットされ
る。ブロック１３４においてＪ％（Ｋ）がＭＩＮと比較
される。もしＪ％（Ｋ）がＭＩＮよりも小さければ、プ
ログラムはブロック１３６へ分岐し；さもなければブロ
ック１３８へ進む。ブロック１３６において、ＭＩＮが
Ｊ％（Ｋ）にセットされる。ブロック１３８において、
インデックスＩが（Ｉ＋Ｊ％（Ｋ）／２＋２）の整数に
セットされる。もしＩが１０，０００以下ならプログラ
ムはブロック１１６（図１１）に戻る（Ｉが１０，００
０以下であるかどうかを判定するテストは図示されてい
ない）。
【００３１】ステップ−２（比較ブロック１０８−１３
８）に関する上記説明との関連で次の点に留意された
い：即ち、ブロック１１８はシンメトリーパラメーター
Ｔ１を計算する。判断ブロック１２０でＩ番目の時間イ
ンターバルの両側で互いに対応する時間インターバルが
対称であると判断される。ブロック１２２で、プログラ
ムはＩ番目の時間の両側で互いに対称に対応する時間イ
ンターバルの個数をカウントする。ブロック１２６にお
いて、プログラムはＩ番目の時間インターバルの両側に
あって互いに対応する４個の時間インターバルが対称な
関係にあるかどうかをチェックし、ブロック１２８にお
いて、プログラムは４個の対応時間インターバルのうち
３個が極限点において判断すべきＩ番目の時間インター
バルを挟んで対称であることを確認する。ブロック１３
０において、プログラムは発見したばかりの極限点（即
ち、Ｋ番目の極限点）の時間インターバルアレイＸ（）
中に占める場所を識別する。これが上方極限点と呼ばれ
る場所である。ブロック１３２において、プログラムは
極限点Ｋと先行の極限点Ｋ−１との間に現われる時間イ
ンターバルの最小個数を求める。ブロック１３４−１３
６において、プログラムは極限点間で測定された時間イ
ンターバルの個数を記録する。この個数が可変ＭＩＮに
割当てられる。ブロック１３８において、極限点間に現
われる次の対称な時間インターバルを飛ばしてインデッ
クスＩが増分される。ステップ１及び２の完了で時間イ
ンターバルアレイＸ（Ｉ）中のすべてのＫ極限点の位置
ＥＸＴ（Ｋ）が判明する。奇数番目のＫを上方極限点、
偶数番目のＫを下方極限点と呼んでもよい。極限点間に
現われる時間インターバルの個数も計算され、可変アレ
イＪ％（Ｋ）に記憶される。Ｊ％（）アレイに記憶され
る最小値は可変ＭＩＮに記憶される。
８）に関する上記説明との関連で次の点に留意された
い：即ち、ブロック１１８はシンメトリーパラメーター
Ｔ１を計算する。判断ブロック１２０でＩ番目の時間イ
ンターバルの両側で互いに対応する時間インターバルが
対称であると判断される。ブロック１２２で、プログラ
ムはＩ番目の時間の両側で互いに対称に対応する時間イ
ンターバルの個数をカウントする。ブロック１２６にお
いて、プログラムはＩ番目の時間インターバルの両側に
あって互いに対応する４個の時間インターバルが対称な
関係にあるかどうかをチェックし、ブロック１２８にお
いて、プログラムは４個の対応時間インターバルのうち
３個が極限点において判断すべきＩ番目の時間インター
バルを挟んで対称であることを確認する。ブロック１３
０において、プログラムは発見したばかりの極限点（即
ち、Ｋ番目の極限点）の時間インターバルアレイＸ（）
中に占める場所を識別する。これが上方極限点と呼ばれ
る場所である。ブロック１３２において、プログラムは
極限点Ｋと先行の極限点Ｋ−１との間に現われる時間イ
ンターバルの最小個数を求める。ブロック１３４−１３
６において、プログラムは極限点間で測定された時間イ
ンターバルの個数を記録する。この個数が可変ＭＩＮに
割当てられる。ブロック１３８において、極限点間に現
われる次の対称な時間インターバルを飛ばしてインデッ
クスＩが増分される。ステップ１及び２の完了で時間イ
ンターバルアレイＸ（Ｉ）中のすべてのＫ極限点の位置
ＥＸＴ（Ｋ）が判明する。奇数番目のＫを上方極限点、
偶数番目のＫを下方極限点と呼んでもよい。極限点間に
現われる時間インターバルの個数も計算され、可変アレ
イＪ％（Ｋ）に記憶される。Ｊ％（）アレイに記憶され
る最小値は可変ＭＩＮに記憶される。
【００３２】ステップ−３
再び図１２において、ステップ−３はＳＴＡＲＴを０に
セットするブロック１４０から始まる。ブロック１４２
において、ＳＴＡＲＴがＳＴＡＲＴ＋１にセットされ
る。ブロック１４４において、Ｊ％（ＳＴＡＲＴ）、Ｊ
％（ＳＴＡＲＴ＋１）及びＭＩＮが比較される。もしこ
の３つが互いに等しくなければプログラムはブロック１
４２に戻り；さもなければブロック１４６へ進む。この
ようにしてプログラムは先ず２つの隣接する極限点がＭ
ＩＮに等しいことを発見する。
セットするブロック１４０から始まる。ブロック１４２
において、ＳＴＡＲＴがＳＴＡＲＴ＋１にセットされ
る。ブロック１４４において、Ｊ％（ＳＴＡＲＴ）、Ｊ
％（ＳＴＡＲＴ＋１）及びＭＩＮが比較される。もしこ
の３つが互いに等しくなければプログラムはブロック１
４２に戻り；さもなければブロック１４６へ進む。この
ようにしてプログラムは先ず２つの隣接する極限点がＭ
ＩＮに等しいことを発見する。
【００３３】ステップ−４
ブロック１４６において、インデックスＩがＳＴＡＲＴ
にセットされ、ＤＥＬＴＡ（Ｉ）がＸ（Ｉ）にセットさ
れる。ブロック１４８において、ＩがＳＴＡＲＴ＋１に
セットされ、ＤＥＬＴＡ（Ｉ）がＸ（Ｉ）にセットされ
る。このようにして最初の２個の修正時間インターバル
（ＤＥＬＴＡ（））は０に等しい合計インデックスＳ１
（）を持つことになる。プログラムはこの２つの極限点
については中心の時間インターバルを合計するだけでよ
い。ブロック１５０において、ＩがＩ＋１にセットさ
れ、Ｑが０にセットされ、ＤＥＬＴＡ（Ｉ）が０にセッ
トされる。ブロック１５２（図１３）において、Ｊ％
（Ｉ−１）がＪ％（Ｉ）と比較される。もし前者が後者
よりも大きければ、プログラムはブロック１５４へ分岐
し；さもなければブロック１５６へ進む。ブロック１５
４において、合計インデックスＳ１（Ｉ）がＳ１（Ｉ−
２）−１にセットされる。ブロック１５６において、Ｊ
％（Ｉ−１）がＪ％（Ｉ）と比較され、もし前者が後者
よりも小さければ、プログラムはブロック１５８へ分岐
し；さもなければブロック１６０へ進む。ブロック１５
８において、Ｓ１（Ｉ）がＳ１（Ｉ−２）＋１にセット
される。ブロック１６０において、合計インデックスＳ
１（Ｉ）がＳ１（Ｉ−２）にセットされる。即ち、ブロ
ック１５０−１６０は合計すべき時間インターバルの個
数（即ち、合計インデックスＳ１（））をこのタイプの
最後の極限点、即ち、上方極限点または下方極限点の前
に現われる極限点間のクロックカウント数Ｊ％（）に基
づいて変更すべきかどうかを判断する。もしクロックカ
ウントＪ％（）が減少すれば、ブロック１５４が合計イ
ンデックスＳ（１）を１だけ減分する。もしクロックカ
ウントＪ％（）が増大すれば、ブロック１５８が合計イ
ンデックスＳ１（）を１だけ増分する。もしクロックカ
ウントが増減しなければ、ブロック１６０が合計インデ
ックスＳ１（）をそのまま変化させない。
にセットされ、ＤＥＬＴＡ（Ｉ）がＸ（Ｉ）にセットさ
れる。ブロック１４８において、ＩがＳＴＡＲＴ＋１に
セットされ、ＤＥＬＴＡ（Ｉ）がＸ（Ｉ）にセットされ
る。このようにして最初の２個の修正時間インターバル
（ＤＥＬＴＡ（））は０に等しい合計インデックスＳ１
（）を持つことになる。プログラムはこの２つの極限点
については中心の時間インターバルを合計するだけでよ
い。ブロック１５０において、ＩがＩ＋１にセットさ
れ、Ｑが０にセットされ、ＤＥＬＴＡ（Ｉ）が０にセッ
トされる。ブロック１５２（図１３）において、Ｊ％
（Ｉ−１）がＪ％（Ｉ）と比較される。もし前者が後者
よりも大きければ、プログラムはブロック１５４へ分岐
し；さもなければブロック１５６へ進む。ブロック１５
４において、合計インデックスＳ１（Ｉ）がＳ１（Ｉ−
２）−１にセットされる。ブロック１５６において、Ｊ
％（Ｉ−１）がＪ％（Ｉ）と比較され、もし前者が後者
よりも小さければ、プログラムはブロック１５８へ分岐
し；さもなければブロック１６０へ進む。ブロック１５
８において、Ｓ１（Ｉ）がＳ１（Ｉ−２）＋１にセット
される。ブロック１６０において、合計インデックスＳ
１（Ｉ）がＳ１（Ｉ−２）にセットされる。即ち、ブロ
ック１５０−１６０は合計すべき時間インターバルの個
数（即ち、合計インデックスＳ１（））をこのタイプの
最後の極限点、即ち、上方極限点または下方極限点の前
に現われる極限点間のクロックカウント数Ｊ％（）に基
づいて変更すべきかどうかを判断する。もしクロックカ
ウントＪ％（）が減少すれば、ブロック１５４が合計イ
ンデックスＳ（１）を１だけ減分する。もしクロックカ
ウントＪ％（）が増大すれば、ブロック１５８が合計イ
ンデックスＳ１（）を１だけ増分する。もしクロックカ
ウントが増減しなければ、ブロック１６０が合計インデ
ックスＳ１（）をそのまま変化させない。
【００３４】ブロック１６４において、Ｓ１（Ｉ）が０
と比較される。もし０に等しければプログラムはブロッ
ク１７０へ分岐し；さもなければブロック１６６に進
む。ブロック１６６において、変数ＱがＱ＋１にセット
される。ブロック１６８において、ＤＥＬＴＡ（Ｉ）が
ＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ＋Ｑ）＋Ｘ（Ｉ−Ｑ）にセッ
トされる。ブロック１７０において、ＱがＳ１（Ｉ）と
比較される。もし双方が等しければ、プログラムはブロ
ック１７２に進み；さもなければブロック１６６に戻
る。従って、ブロック１６４−１７０は中心の時間イン
ターバルの両側の２つのＳ１（）時間インターバルを合
計する（多重フォールドオーバーの場合にはこれを２回
以上行えばよい）。ブロック１７２において、ＤＥＬＴ
Ａ（Ｉ）がＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ）にセットされ
る。ブロック１７２は中心の時間インターバルを加え
る。ブロック１７４において、ＩがＫと比較される。も
しＩがＫよりも小さければ、プログラムはブロック１５
０（図１２）に戻り；さもなければブロック１７６（図
１４）に進む。即ち、もしＩがＫよりも小さければ、プ
ログラムは次の上限を分析するためスターとに戻り；さ
もなければ、ルーチンを出て逆方向チェーンを開始す
る。
と比較される。もし０に等しければプログラムはブロッ
ク１７０へ分岐し；さもなければブロック１６６に進
む。ブロック１６６において、変数ＱがＱ＋１にセット
される。ブロック１６８において、ＤＥＬＴＡ（Ｉ）が
ＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ＋Ｑ）＋Ｘ（Ｉ−Ｑ）にセッ
トされる。ブロック１７０において、ＱがＳ１（Ｉ）と
比較される。もし双方が等しければ、プログラムはブロ
ック１７２に進み；さもなければブロック１６６に戻
る。従って、ブロック１６４−１７０は中心の時間イン
ターバルの両側の２つのＳ１（）時間インターバルを合
計する（多重フォールドオーバーの場合にはこれを２回
以上行えばよい）。ブロック１７２において、ＤＥＬＴ
Ａ（Ｉ）がＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ）にセットされ
る。ブロック１７２は中心の時間インターバルを加え
る。ブロック１７４において、ＩがＫと比較される。も
しＩがＫよりも小さければ、プログラムはブロック１５
０（図１２）に戻り；さもなければブロック１７６（図
１４）に進む。即ち、もしＩがＫよりも小さければ、プ
ログラムは次の上限を分析するためスターとに戻り；さ
もなければ、ルーチンを出て逆方向チェーンを開始す
る。
【００３５】図１４のブロック１７６において、プログ
ラムはインデックスＩをＳＴＡＲＴにセットする。ブロ
ック１７８において、ＩがＩ−１にセットされ、Ｑが０
にセットされ、ＤＥＬＴＡ（Ｉ）が０にセットされる。
ブロック１８０において、Ｊ％（Ｉ＋２）がＪ％（Ｉ＋
１）と比較される。もし前者が後者よりも大きければ、
プログラムはブロック１８２へ分岐し；さもなければブ
ロック１８４へ進む。ブロック１８２において、Ｓ１
（Ｉ）がＳ１（Ｉ＋２）−１にセットされる。ブロック
１８４において、Ｊ％（Ｉ＋２）がＪ％（Ｉ＋１）と比
較される。もし前者が後者よりも小さければ、プログラ
ムはブロック１８６に分岐し；さもなければブロック１
８８に進む。ブロック１８６において、Ｓ１（Ｉ）がＳ
１（Ｉ＋２）＋１にセットされ、ブロック１８８におい
て、Ｓ１（Ｉ）がＳ１（Ｉ＋２）にセットされる。即
ち、ブロック１８０−１８８はこのタイプの最後の時間
インターバル、即ち、上方極限点または下方極限点のあ
とに現われる極限点間のクロックカウント数に基づいて
合計すべき時間インターバルの個数を変更すべきかどう
かを判定する。もしクロックカウントＪ％（）が増大す
ればブロック１８２において、プログラムは合計インデ
ックスＳ１（）を減分し、もしクロックカウントＪ
％（）が減少すれば、ブロック１８６において、プログ
ラムは合計インデックスＳ１（）を増分する。クロック
カウントに増減がなければブロック１８８において合計
インデックスＳ１（）はそのまま変化しない。
ラムはインデックスＩをＳＴＡＲＴにセットする。ブロ
ック１７８において、ＩがＩ−１にセットされ、Ｑが０
にセットされ、ＤＥＬＴＡ（Ｉ）が０にセットされる。
ブロック１８０において、Ｊ％（Ｉ＋２）がＪ％（Ｉ＋
１）と比較される。もし前者が後者よりも大きければ、
プログラムはブロック１８２へ分岐し；さもなければブ
ロック１８４へ進む。ブロック１８２において、Ｓ１
（Ｉ）がＳ１（Ｉ＋２）−１にセットされる。ブロック
１８４において、Ｊ％（Ｉ＋２）がＪ％（Ｉ＋１）と比
較される。もし前者が後者よりも小さければ、プログラ
ムはブロック１８６に分岐し；さもなければブロック１
８８に進む。ブロック１８６において、Ｓ１（Ｉ）がＳ
１（Ｉ＋２）＋１にセットされ、ブロック１８８におい
て、Ｓ１（Ｉ）がＳ１（Ｉ＋２）にセットされる。即
ち、ブロック１８０−１８８はこのタイプの最後の時間
インターバル、即ち、上方極限点または下方極限点のあ
とに現われる極限点間のクロックカウント数に基づいて
合計すべき時間インターバルの個数を変更すべきかどう
かを判定する。もしクロックカウントＪ％（）が増大す
ればブロック１８２において、プログラムは合計インデ
ックスＳ１（）を減分し、もしクロックカウントＪ
％（）が減少すれば、ブロック１８６において、プログ
ラムは合計インデックスＳ１（）を増分する。クロック
カウントに増減がなければブロック１８８において合計
インデックスＳ１（）はそのまま変化しない。
【００３６】ブロック１９０において、Ｓ１（Ｉ）が０
と比較される。もし０に等しければ、プログラムはブロ
ック１９２へ進み；さもなければ、ブロック１９４に分
岐する。ブロック１９２において、ＤＥＬＴＡ（Ｉ）が
ＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ）にセットされる。ブロック
１９４において、ＱがＱ＋１にセットされる。ブロック
１９６において、ＤＥＬＴＡ（Ｉ）がＤＥＬＴＡ（Ｉ）
＋Ｘ（Ｉ＋Ｑ）＋Ｘ（Ｉ−Ｑ）にセットされる。ブロッ
ク１９８において、ＱがＳ１（Ｉ）と比較される。もし
両者がで等しければ、プログラムはブロック１９２へ分
岐し；さもなければブロック１９４へ戻る。即ち、ブロ
ック１９０−１９８は中心の時間インターバルＩの両側
の２つのＳ１（Ｉ）時間インターバルを合計し、多重フ
ォールドオーバーが存在するならこのプロセスを繰り返
す。ブロック１９２は中心の時間インターバルＩを加え
る。ブロック２００において、Ｉが１と比較される。も
しＩが１よりも大きければ、プログラムはブロック１７
８に戻り；さもなければプログラムは完了する。即ち、
プログラムはＩ＝１となるフロントエンドに達するまで
スターとに戻って次の下方極限点を分析する。
と比較される。もし０に等しければ、プログラムはブロ
ック１９２へ進み；さもなければ、ブロック１９４に分
岐する。ブロック１９２において、ＤＥＬＴＡ（Ｉ）が
ＤＥＬＴＡ（Ｉ）＋Ｘ（Ｉ）にセットされる。ブロック
１９４において、ＱがＱ＋１にセットされる。ブロック
１９６において、ＤＥＬＴＡ（Ｉ）がＤＥＬＴＡ（Ｉ）
＋Ｘ（Ｉ＋Ｑ）＋Ｘ（Ｉ−Ｑ）にセットされる。ブロッ
ク１９８において、ＱがＳ１（Ｉ）と比較される。もし
両者がで等しければ、プログラムはブロック１９２へ分
岐し；さもなければブロック１９４へ戻る。即ち、ブロ
ック１９０−１９８は中心の時間インターバルＩの両側
の２つのＳ１（Ｉ）時間インターバルを合計し、多重フ
ォールドオーバーが存在するならこのプロセスを繰り返
す。ブロック１９２は中心の時間インターバルＩを加え
る。ブロック２００において、Ｉが１と比較される。も
しＩが１よりも大きければ、プログラムはブロック１７
８に戻り；さもなければプログラムは完了する。即ち、
プログラムはＩ＝１となるフロントエンドに達するまで
スターとに戻って次の下方極限点を分析する。
【００３７】フーリエ変換または高速フーリエ変換を行
なうため、プログラムは捕捉されるセンサーサイクル数
を８．５３ｓｅｃのデータに相当するＫ＝１０２４とす
る。従属変数ＤＥＬＴＡ（２Ｉ）、Ｉ＝１−５１２は単
位時間を持つが、この信号の時間（周波数）変動だけが
利用されるのであるから、上記単位時間を持つことはさ
して重要ではない。従って、変数２Ｉは時間をも表わ
す。２Ｉ＝２から２Ｉ＝４までの時間インターバルは１
／１２０ｓｅｃである。これは下方極限点の場合であ
る。
なうため、プログラムは捕捉されるセンサーサイクル数
を８．５３ｓｅｃのデータに相当するＫ＝１０２４とす
る。従属変数ＤＥＬＴＡ（２Ｉ）、Ｉ＝１−５１２は単
位時間を持つが、この信号の時間（周波数）変動だけが
利用されるのであるから、上記単位時間を持つことはさ
して重要ではない。従って、変数２Ｉは時間をも表わ
す。２Ｉ＝２から２Ｉ＝４までの時間インターバルは１
／１２０ｓｅｃである。これは下方極限点の場合であ
る。
【００３８】上方極限点については、従属変数ＤＥＬＴ
Ａ（２Ｉ−１）、Ｉ＝１−５１２はデータ点２Ｉ−１＝
１から２Ｉ−１＝３までの時間インターバルであり、こ
こでも１／１２０ｓｅｃに相当する。ＤＥＬＴＡ（２
Ｉ）及びＤＥＬＴＡ（２Ｉ−１）に対するフーリエ変換
またはＦＦＴにより、分解能が０．１１７Ｈｚの全く同
じセンサーうなり周波数ＦB が得られる。センサー固有
周波数ＦN は１２０Ｈｚ＋ＦB に等しい。センサーの補
正率“ＣＦ”は下式によって与えられる。
Ａ（２Ｉ−１）、Ｉ＝１−５１２はデータ点２Ｉ−１＝
１から２Ｉ−１＝３までの時間インターバルであり、こ
こでも１／１２０ｓｅｃに相当する。ＤＥＬＴＡ（２
Ｉ）及びＤＥＬＴＡ（２Ｉ−１）に対するフーリエ変換
またはＦＦＴにより、分解能が０．１１７Ｈｚの全く同
じセンサーうなり周波数ＦB が得られる。センサー固有
周波数ＦN は１２０Ｈｚ＋ＦB に等しい。センサーの補
正率“ＣＦ”は下式によって与えられる。
【００３９】
ＣF ＝（ＦN 2 −ＦO 2) ／（ＦI 2 −ＦO 2)
ただし、ＦO はセンサーに加えられる励振または駆動周
波数（典型的には１２０Ｈｚ）、ＦI はセンサーの設計
上周波数（典型的には１３２．５Ｈｚ）である。
波数（典型的には１２０Ｈｚ）、ＦI はセンサーの設計
上周波数（典型的には１３２．５Ｈｚ）である。
【００４０】以上に本発明の多様な特徴及び利点を説明
したが、本発明の思想及び範囲に含まれるこのような特
徴及び利点のすべては頭書した特許請求の範囲によって
包含されるものとする。
したが、本発明の思想及び範囲に含まれるこのような特
徴及び利点のすべては頭書した特許請求の範囲によって
包含されるものとする。
【００４１】
【図１】発電機に取付けられるＦＯＶＭの簡略図。
【図２】光学式振動センサー１２の詳細図。
【図３】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図４】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図５】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図６】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図７】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図８】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図９】ＦＯＶＭセンサーのうなり周波数測定における
本発明システムのすぐれた性能を示す波形図。
本発明システムのすぐれた性能を示す波形図。
【図１０】ＦＯＶＭセンサーのうなり周波数測定におけ
る本発明システムのすぐれた性能を示す波形図。
る本発明システムのすぐれた性能を示す波形図。
【図１１】本発明の好ましい実施例のフローチャート。
【図１２】本発明の好ましい実施例のフローチャート。
【図１３】本発明の好ましい実施例のフローチャート。
【図１４】本発明の好ましい実施例のフローチャート。
１０ 発電機
１２ 光学式振動センサー
１４ コンピュター
１６ 振動感知システム
２０ ハウジング
２２ 光電変換手段
２４ 内部リード
２６ 格子組立体
フロントページの続き
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01H 9/00
G01H 13/00
G01M 19/00
Claims (6)
- 【請求項１】 装置に取付けられた振動感知システムに
おけるうなり周波数を測定する方法において、 （ａ）装置の振動周波数及び振幅を示す振動信号を発生
させ； （ｂ）振動信号のうなり周波数に対応する極限点を含む
波形部分の２つのゼロ交差間の時間インターバルを表わ
す時間インターバルデータを記憶し、処理することによ
って振動信号からうなり周波数を得るステップより成る
ことを特徴とするうなり周波数を測定する方法。 - 【請求項２】 ステップ（ｂ）が、うなり周波数に対応
する極限点でピーク値を有する変換済み時間インターバ
ルデータを得るために時間インターバルデータをフーリ
エ変換することを含む請求項１に記載の方法。 - 【請求項３】 ステップ（ｂ）が、振動信号の連続する
ゼロ交差間のクロックカウントを合計することによって
時間インターバルデータを計算することを含む請求項１
または２に記載の方法。 - 【請求項４】 ステップ（ｂ）が、時間インターバルデ
ータのフォールドオーバーに対応する極限点の時間イン
ターバルの位置を見出し、前記フォールドオーバーによ
り影響を受ける極限点間で合計されるクロックカウント
数を変更することによって前記フォールドオーバーを補
償することを含む請求項３に記載の方法。 - 【請求項５】 前記装置が端巻きコイルの振動が１２０
Ｈｚの発電機であり、発電機に取付けられている振動セ
ンサーによって振動信号が得られ、振動信号が発電機の
端巻きコイルの１２０Ｈｚの振動及び振動センサーの振
動周波数に起因する機械的振動信号を含むことを特徴と
する上記請求項のうち任意の請求項に記載の方法。 - 【請求項６】 前記装置が端巻きコイルの振動が１２０
Ｈｚの発電機であり、発電機に取付けられている振動セ
ンサーによって振動信号が得られ、振動信号が発電機の
端巻きコイルの１２０Ｈｚの振動及び振動センサーの共
振周波数に起因する機械的振動信号を含み；ステップ
（ｂ）が振動信号の連続するゼロ交差間のクロックカウ
ントを合計することによって時間インターバルデータを
計算し、時間インターバルデータのフォールドオーバー
に対応する極限点の時間インターバルの位置を見出し、
前記フォールドオーバーにより影響を受ける極限点間で
合計されるクロックカウント数を変更することによって
前記フォールドオーバーを補償し、うなり周波数に対応
する点でピーク値を有する変換済み時間インターバルデ
ータを得るために時間インターバルデータをフーリエ変
換することから成ることを特徴とする請求項１に記載の
方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/069,495 US5469745A (en) | 1993-06-01 | 1993-06-01 | System for determining FOVM sensor beat frequency |
| US08/069495 | 1993-06-01 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06347318A JPH06347318A (ja) | 1994-12-22 |
| JP3404129B2 true JP3404129B2 (ja) | 2003-05-06 |
Family
ID=22089375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP14397994A Expired - Lifetime JP3404129B2 (ja) | 1993-06-01 | 1994-06-01 | 光学式振動センサーによるうなり周波数測定方法 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US5469745A (ja) |
| EP (1) | EP0627616A3 (ja) |
| JP (1) | JP3404129B2 (ja) |
| KR (1) | KR100317014B1 (ja) |
| CN (1) | CN1043588C (ja) |
| CA (1) | CA2124761A1 (ja) |
| MX (1) | MX9404100A (ja) |
| SI (1) | SI9400245A (ja) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL9401949A (nl) * | 1994-11-22 | 1996-07-01 | Skf Ind Trading & Dev | Werkwijze voor het analyseren van regelmatig geëxciteerde mechanische trillingen. |
| DE19705922A1 (de) * | 1997-02-17 | 1998-08-20 | Asea Brown Boveri | Verfahren zur Bestimmung des Festigkeitszustandes von Wickelköpfen elektrischer Maschinen und Anordnung zur Durchführung des Verfahrens |
| KR100477504B1 (ko) * | 2002-12-09 | 2005-03-17 | 한국전기연구원 | 에프에프티를 이용한 실시간 전력계통 진동 감시장치 |
| US20040212353A1 (en) * | 2003-04-25 | 2004-10-28 | Siemens Westinghouse Power Corporation | Use of a closing impedance to minimize the adverse impact of out-of-phase generator synchronization |
| DE10323063A1 (de) * | 2003-05-20 | 2004-12-09 | Endress + Hauser Gmbh + Co. Kg | Verfahren zur Füllstandsmessung |
| US7533572B2 (en) * | 2006-08-15 | 2009-05-19 | Siemens Energy, Inc. | High bandwidth fiber optic vibration sensor |
| KR100883394B1 (ko) * | 2007-04-03 | 2009-02-11 | 주식회사 대웅컨설턴트 | 계측 관리형 진동현의 고유 주파수 측정 방법 |
| CN101299635B (zh) * | 2008-05-09 | 2010-09-15 | 东南大学 | 宽域全光纤扰动信号前-后判别装置及方法 |
| US20110018483A1 (en) * | 2009-07-21 | 2011-01-27 | General Electric Company | Stator end-winding component monitoring system |
| US8379225B2 (en) * | 2010-03-31 | 2013-02-19 | General Electric Company | System for monitoring a relative displacement of components |
| EP2634898A1 (de) * | 2012-03-01 | 2013-09-04 | Siemens Aktiengesellschaft | Verfahren zum Betreiben einer elektrischen Maschine |
| US10254773B2 (en) * | 2012-09-25 | 2019-04-09 | Lord Corporation | Devices, systems, and methods for balancing closely coupled rotating machinery |
| CN103558008B (zh) * | 2013-11-20 | 2017-02-08 | 无锡波汇光电科技有限公司 | 一种风洞试验中光纤传感振动信号的采集装置 |
| US9551598B2 (en) | 2014-05-12 | 2017-01-24 | Siemens Energy, Inc. | Fiber optic sensing apparatus with an improved fiber-affixing device |
| US9939458B2 (en) | 2015-08-27 | 2018-04-10 | General Electric Company | Insulated accelerometer assembly for high voltage environment |
| US11709244B2 (en) * | 2019-10-21 | 2023-07-25 | Banner Engineering Corp. | Near range radar |
| CN111397877B (zh) * | 2020-04-02 | 2021-07-27 | 西安建筑科技大学 | 一种旋转机械拍振故障检测与诊断方法 |
| CN113091879B (zh) * | 2021-04-01 | 2022-09-06 | 陕西电力建设集团有限公司 | 一种高压输电线路导线风振检测装置 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4283634A (en) * | 1971-06-23 | 1981-08-11 | Westinghouse Electric Corp. | System and method for monitoring and controlling operation of industrial gas turbine apparatus and gas turbine electric power plants preferably with a digital computer control system |
| US4218614A (en) * | 1978-06-08 | 1980-08-19 | Westinghouse Electric Corp. | Device for measuring vibration amplitude |
| US4321464A (en) * | 1978-06-08 | 1982-03-23 | Westinghouse Electric Corp. | Device for measuring vibration phase and amplitude |
| DE3112188A1 (de) * | 1980-03-31 | 1982-01-14 | Hitachi, Ltd., Tokyo | Verfahren und vorrichtung zur ueberwachung der schwingung einer sich drehenden welle |
| US4429609A (en) * | 1981-12-14 | 1984-02-07 | Warrender David J | Pitch analyzer |
| JPS61207195A (ja) * | 1985-03-08 | 1986-09-13 | Mitsubishi Electric Corp | モ−タ制御装置 |
| JPH0364294A (ja) * | 1989-08-02 | 1991-03-19 | Nec Corp | 周波数検出回路 |
| US5109700A (en) * | 1990-07-13 | 1992-05-05 | Life Systems, Inc. | Method and apparatus for analyzing rotating machines |
| US5069071A (en) * | 1990-08-27 | 1991-12-03 | United Technologies Corporation | Vibration monitoring in the frequency domain with a capacitive accelerometer |
| US5146776A (en) * | 1990-11-26 | 1992-09-15 | Westinghouse Electric Corp. | Method for continuously calibrating an optical vibration sensor |
-
1993
- 1993-06-01 US US08/069,495 patent/US5469745A/en not_active Expired - Lifetime
-
1994
- 1994-05-30 CN CN94106304A patent/CN1043588C/zh not_active Expired - Fee Related
- 1994-05-31 KR KR1019940012086A patent/KR100317014B1/ko not_active IP Right Cessation
- 1994-05-31 CA CA002124761A patent/CA2124761A1/en not_active Abandoned
- 1994-06-01 EP EP94303943A patent/EP0627616A3/en not_active Withdrawn
- 1994-06-01 SI SI9400245A patent/SI9400245A/sl unknown
- 1994-06-01 MX MX9404100A patent/MX9404100A/es not_active IP Right Cessation
- 1994-06-01 JP JP14397994A patent/JP3404129B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| SI9400245A (en) | 1994-12-31 |
| CN1043588C (zh) | 1999-06-09 |
| KR950002171A (ko) | 1995-01-04 |
| US5469745A (en) | 1995-11-28 |
| CA2124761A1 (en) | 1994-12-02 |
| KR100317014B1 (ko) | 2002-02-28 |
| JPH06347318A (ja) | 1994-12-22 |
| EP0627616A3 (en) | 1997-07-09 |
| MX9404100A (es) | 1995-01-31 |
| EP0627616A2 (en) | 1994-12-07 |
| CN1097872A (zh) | 1995-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3404129B2 (ja) | 光学式振動センサーによるうなり周波数測定方法 | |
| JPH0612296B2 (ja) | 振動センサの較正方法及び振動検知装置 | |
| EP2111535B1 (en) | Vibrating wire sensor using spectral analysis | |
| US4671659A (en) | Fiber optic displacement sensor | |
| US6851313B2 (en) | Vibratory level sensor | |
| US8671758B2 (en) | System and method for measuring the frequency of a vibrating object | |
| US20100324852A1 (en) | Determining the resonance parameters for mechanical oscillators | |
| JP3680678B2 (ja) | 振動測定方法及び振動数測定装置 | |
| JPH10253339A (ja) | 音波利用計測方法及び計測装置 | |
| JP4976924B2 (ja) | 光式ガス検知方法及び光式ガス検知装置 | |
| US7134344B2 (en) | Control method and apparatus | |
| JP2005257570A (ja) | 光ファイバによる構造物の診断方法及びシステム | |
| JP4058609B2 (ja) | 光ファイバリング干渉計型振動検知センサにおける振動位置推定装置および振動位置推定方法 | |
| JP2005037210A (ja) | 振動エネルギー損失に関する係数の測定方法及び装置 | |
| SU798185A1 (ru) | Способ контрол процесса вибро-ОбРАбОТКи КОНСТРуКций и уСТРОйСТВОдл ЕгО ОСущЕСТВлЕНи | |
| JP3272232B2 (ja) | 振動変位検出装置 | |
| RU2047835C1 (ru) | Способ определения угла наклона | |
| JPH09159521A (ja) | アクチュエータの共振周波数検出回路 | |
| SU721678A1 (ru) | Способ определени двух компонент механических колебаний конструкций и устройство дл его осуществлени | |
| SU1732179A1 (ru) | Способ определени геометрических размеров микропроволоки | |
| RU2101686C1 (ru) | Способ измерения периодического закона движения вибрирующей поверхности | |
| KR20010105982A (ko) | 가속도 센서를 이용한 휴대용 동적 변위 측정 장치 | |
| RU2017297C1 (ru) | Способ определения пляски проводов | |
| SU619863A1 (ru) | Установка дл определени чувствительности пьезоаксеерометров к деформации объекта измерени | |
| SU1294828A1 (ru) | Устройство дл измерени плотности сгущаемых растворов сахарного производства |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20030204 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| 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: 20080229 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090228 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090228 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100228 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110228 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110228 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120229 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130228 Year of fee payment: 10 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140228 Year of fee payment: 11 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |