JP2014114806A5 - - Google Patents
Download PDFInfo
- Publication number
- JP2014114806A5 JP2014114806A5 JP2013245535A JP2013245535A JP2014114806A5 JP 2014114806 A5 JP2014114806 A5 JP 2014114806A5 JP 2013245535 A JP2013245535 A JP 2013245535A JP 2013245535 A JP2013245535 A JP 2013245535A JP 2014114806 A5 JP2014114806 A5 JP 2014114806A5
- Authority
- JP
- Japan
- Prior art keywords
- blade
- blades
- determining
- arrival time
- coefficients
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 claims description 86
- 230000000875 corresponding Effects 0.000 claims description 80
- 230000003068 static Effects 0.000 claims description 64
- 238000000034 method Methods 0.000 claims description 46
- 230000000694 effects Effects 0.000 claims description 14
- 238000000513 principal component analysis Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000009434 installation Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 20
- 239000012530 fluid Substances 0.000 description 4
- 210000001503 Joints Anatomy 0.000 description 3
- 238000005452 bending Methods 0.000 description 3
- 241000699666 Mus <mouse, genus> Species 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000006356 dehydrogenation reaction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000004231 fluid catalytic cracking Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006011 modification reaction Methods 0.000 description 1
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 description 1
- 239000001294 propane Substances 0.000 description 1
- 230000000717 retained Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Description
開示内容の実施形態は一般的に、ロータブレード又は翼形部の健全性をモニタするためのシステム及び方法に関する。
ロータブレード又は翼形部は、多くの装置で重要な役割を演じている。いくつかの例を挙げれば、軸流式圧縮機、タービン、エンジン、ターボ機械などである。例えば、軸流式圧縮機には一連の段があり、各段は、ロータブレード又は翼形部の列と、それに続く静翼又は静止翼形部の列とを備えている。すなわち、各段は、一対のロータブレード又は翼形部と静止翼形部とを備えている。典型的には、ロータブレード又は翼形部によって、入口を通って軸流式圧縮機に入る流体の運動エネルギーが増加する。さらに、静翼又は静止翼形部によって一般的に、増加した流体の運動エネルギーが拡散を通して静圧に変換される。したがって、ロータブレード又は翼形部と静止翼形部とによって、流体の圧力が増加する。
さらに、ロータブレード又は翼形部と静止翼形部とを備える軸流式圧縮機の応用例は、多種多様である。軸流式圧縮機を、例えば、複数の装置で用いる場合がある。例えば、陸用ガスタービン、ジェットエンジン、高速船エンジン、小規模発電所などである。加えて、軸流式圧縮機には他の応用例がある場合がある。例えば、大規模空気分離装置、高炉空気、流動接触分解空気、プロパン脱水素などである。
翼形部は、翼形部の健全性に影響する極端かつ変動する動作条件、例えば、高速、圧力及び温度の下で長期間動作する。極端かつ変動する動作条件に加えて、他の特定の因子によって、翼形部の疲労及び応力を招くこともある。そのような因子としては、例えば、遠心力を含む慣性力、圧力、翼形部の共振周波数、翼形部内の振動、振動応力、温度応力、翼形部の再着座、ガスその他の流体の負荷などが挙げられる場合がある。応力及び疲労の長期的増加がある時間に渡って起こると、欠陥及び割れが翼形部内に生じる。割れの1以上が時間とともに広がって、翼形部又は翼形部の一部が遊離する場合がある。翼形部の遊離は、翼形部を備える装置にとって危険な場合があり、そのため、莫大な金銭的損失に至る場合がある。加えて、装置の付近にいる人達にとって安全ではない場合がある。
したがって、翼形部の健全性をリアルタイムで予測し得るシステム及び方法を開発することは非常に望ましい。より詳細には、割れ又は破砕をリアルタイムで検出及び予測し得るシステム及び方法を開発することは望ましい。
簡潔に述べると、本技術の一態様では、方法を提供する。本方法は、複数のブレードに対する実際の到着時間に基づいて複数のブレードに対する正規化Δ到着時間を決定するステップと、複数のブレードに対する正規化Δ到着時間から1以上の共通因子の影響を取り除くことによって複数のブレードの静的たわみを決定するステップとを含む。
一態様では、処理サブシステムを含むシステムを提供する。処理サブシステムは、複数のブレードに対する実際の到着時間に基づいて複数のブレードに対する正規化Δ到着時間を決定し、また複数のブレードに対する正規化Δ到着時間から1以上の共通因子の影響を取り除くことによって複数のブレードの静的たわみを生成する。
本システムの別の態様では、処理サブシステムを提供する。処理サブシステムは、複数ブレードに対する正規化Δ到着時間に基づいて複数のブレードに対する複数のモードを決定し、正規化Δ到着時間に基づいて複数のモードと複数のブレードとに対応する複数のブレード係数を決定し、複数のブレード係数のうち複数のモードの共通モードに対応する1以上のブレード係数を特定し、複数のブレード係数における1以上のブレード係数をゼロに等しくして再構成行列を生成し、正規化Δ到着時間と複数のモードとに基づいてブレードに対する静的たわみを決定する。
本発明の上記その他の特徴、態様及び利点は、以下の詳細な説明を、添付図面を参照して読むと理解を深めることができよう。図面の全体に渡って同様の文字は同様の部分を表す。
以下に詳細に説明するように、本システム及び技術の実施形態によって、1以上のブレード又は翼形部の健全性が評価される。より詳細には、本システム及び技術によって、ブレード又は翼形部の静的たわみが決定される。例えば、ブレードの静的たわみを用いてブレードの健全性をモニタしてもよい。以下、「翼形部」及び「ブレード」という用語は互換的に用いられる。静的たわみは、例えば、ブレードの予測又は原位置からのブレードの原位置又は予測位置の定常変化を示すのに用いることができる。
動作時には、ブレードの回転中に、ブレード内の1以上の割れ又は欠陥に起因して基準位置でのブレードの到着時間(TOA)(以下、実TOAという)が予想TOAから変化する場合がある。したがって、ブレードのTOAの変化を用いて、ブレードの静的たわみを決定することができる。本明細書で用いる「予想TOA」という用語は、ブレード内に欠陥も割れも存在せずにブレードが理想的な状況で動作し、負荷状態が最適で、ブレード内の振動が最小であるときの基準位置でのブレードのTOAをいう。
ブレードの割れ又は欠陥に加えて、実TOAは、1以上の共通因子の影響に起因して変化する場合もある。本明細書で用いる「共通因子」という用語は、装置内のブレードに共通する理由であって、ブレードに対する実TOAに影響(例えば、早着又は遅延)する理由をいう。共通因子としては、例えば、動作パラメータ、ブレードの再着座などが挙げられる。例えば、動作パラメータとしては、入口案内翼(IGV)角度、負荷変化、ブレードの再着座、速度変化、温度、速度などが挙げられる。
本明細書で用いる「ブレードの再着座」という用語は、継手(例えば、ダブテール継手)内のブレードの原位置又は予想位置とは異なる位置にブレードがロックされることをいう。典型的に、ブレードはロータに、1以上の継手(例えばダブテール継手)を介して固定されている。ブレードを含む装置の起動中に、ブレードは、継手内のその原位置からシフトする場合があり、また継手内で、ブレードの原位置とは異なる位置にロックされる場合がある。例として、装置には、ガスタービン、圧縮機などが含まれていてもよい。ブレードの原位置とは異なる位置の継手内にブレードがロックされることを、ブレードの再着座という。ブレードの位置が変わることによって、ブレードの実TOAが変わる場合がある。
その結果、ブレードに対する実TOAに対応する共通因子及びブレード内の割れ又は欠陥の影響に起因して、ブレードの静的たわみは、正確又は精密な静的たわみから変化する。したがって、ブレードの健全性のモニタリング又はブレード内の割れ又は欠陥の決定に際しては、ブレードに対する実TOAに対する共通因子の影響を無効にすることが望ましい。本システム及び技術の特定の実施形態では、共通因子の影響を取り除いて、ブレードの静的たわみを決定する。
図1は、本システムの実施形態によるブレード健全性モニタリングシステム10の概略図である。図1に示すように、システム10は、1以上のブレード又は翼形部12を備え、これらは、ブレードの静的たわみ12を決定するためにシステム10によってモニタされている。さらに、システム10は、ブレードの静的たわみ12に基づいてブレード12の健全性を決定する。現在想定される構成で示すように、システム10は1以上のセンサ14、16を備えている。各センサ14、16は、ブレード通過信号(BPS)18、20をそれぞれ発生させる。これらは、基準点へのブレード12の実際の到着時間(TOA)を表す。一実施形態では、センサ14、16は、基準点への1以上のブレード12の到着を検知して、BPS18、20を発生させる。例えば、基準点は、センサ14、16の真下であってもよいし、或いはセンサ14、16に隣接していてもよい。一実施形態では、各BPS18及び20は所定の時間サンプリング及び/又は測定され、ブレードの実TOAの決定に用いられる。例えば、実TOAは、時間単位又は度の単位で測定してもよい。
一実施形態では、センサ14、16は、1以上のブレード12の前縁の到着を検知して、BPS18、20を発生させてもよい。別の実施形態では、センサ14、16は、1以上のブレード12の後縁の到着を検知して、BPS18、20を発生させてもよい。さらに他の実施形態では、センサ14は、1以上のブレード12の前縁の到着を検知して、BPS18を発生させてもよく、センサ16は、1以上のブレード12の後縁の到着を検知してBPS20を発生させてもよいし、その逆であってもよい。例えば、センサ14、16を、1以上のブレード12に隣接して、静止物体上に、1以上のブレード12の到着が効率的に検知され得るような位置に取り付けてもよい。一実施形態では、センサ14、16の少なくとも一方が、1以上のブレード12のケーシング(図示せず)上に取り付けられている。非限定的な例として、センサ14、16は、磁気センサ、容量センサ、渦電流センサなどであってもよい。
現在想定される構成で例示するように、BPS18、20は、処理サブシステム22によって受信される。処理サブシステム22は、BPS18、20に基づいて1以上のブレード12の実TOAを決定する。さらに、処理サブシステム22は、1以上のブレード12の実TOAに基づいて1以上のブレード12の静的たわみを決定する。より詳細には、処理サブシステム22は、1以上のブレード12の実TOAを処理することによってブレード12の1以上の静的たわみを決定するように構成されている。ブレード12の実TOAは、1以上の共通因子に起因して影響を受ける場合がある。本明細書で用いる「共通因子」という用語は、装置内のすべてのブレードに共通する理由であって、ブレードに対する実TOAに影響(例えば、早着又は遅延)する理由をいう。共通因子としては、例えば、動作パラメータ、ブレードの再着座などが挙げられる。
ブレードの静的たわみ12を、共通因子の影響を取り除かずにこのような実TOAに基づいて決定した場合には、たとえブレード12内に割れ又は欠陥が存在していなくても、このような静的たわみがブレード12の1以上における割れを誤って示唆することがある。したがって、現在想定される技術では、処理サブシステム22は、例えば、ブレード12に対する実TOAから共通因子の影響を取り除くことによってブレードの静的たわみ12を決定する。一実施形態では、例えば、処理サブシステム22は、ブレード12に対する実TOAに基づいて決定される正規化ΔTOAから共通因子の影響を取り除くことによってブレード12に対する静的たわみを決定することができる。例えば、主成分分析技術、特異値分解技術、独立成分分析技術又はそれらの組合せを含む技術を適用することによって共通因子の影響を取り除いてもよい。本明細書で用いる「正規化ΔTOA」という用語は、複数のブレードにおけるあるブレードの実TOAに対する数値であって、複数のブレードに対する実TOAと、ブレード間隔パラメータとに基づいて決定される数値をいう。正規化ΔTOAとブレード間隔パラメータとの決定について、図2及び図3を参照して詳細に説明する。
一実施形態では、処理サブシステム22は、ブレード12に対する静的たわみを、主成分分析技術(PCA)を実TOA又は正規化ΔTOAに適用することによって決定する。したがって、処理サブシステム22は、主成分分析技術を実TOAに適用することによって、実TOAから共通因子の影響を取り除いてもよい。一実施形態では、処理サブシステムは、主成分分析技術、特異値分解技術、独立成分分析技術又はそれらの組合せを含む技術を用いて実TOAから共通因子の影響を取り除く。静的たわみの決定について、図2〜図6を参照してさらに詳しく説明する。一実施形態では、処理サブシステム22は、データリポジトリ24を有していてもよい。データリポジトリ24は、データ、例えば静的たわみ、動的たわみ、TOA、ΔTOA、任意の中間データなどを記憶する。
なお、現在想定される構成では、共通因子(例えば動作パラメータ)の影響は、共通因子(例えば動作パラメータ、ブレードの再着座など)に関係する任意のデータに依拠せずに、実TOAから除去される。したがって、現在想定される技術では、外部装置(例えばオンサイトのモニタリング装置その他の装置)からのデータは、実TOA又は正規化ΔTOAから共通因子の影響を除去するのに必要とされない。動作パラメータとしては、例えば、入口案内翼(IGV)角度、負荷変化、ブレードの再着座、速度変化、温度、速度などが挙げられる。
次に、図2を参照すると、本技術の実施形態により複数のブレードの静的たわみを決定するための典型的な方法200を表すフローチャートが示してある。複数のブレードは、例えば、ブレード12(図1参照)であってもよい。以下、理解を促すため、ブレード12を参照して静的たわみの決定について説明する。ステップ202において、ブレード12に対するブレード通過信号(BPS)を、処理サブシステム、例えば処理サブシステム22(図1参照)で受信してもよい。図1を参照して既に説明したように、BPSは、センサ、例えばセンサ14、16(図1参照)によって発生させてもよい。BPSは、例えば、BPS18、20(図1参照)であってもよい。
さらに、ステップ204において、ブレード12の実際の到着時間(TOA)206を処理サブシステムによって決定する。処理サブシステムは、BPSのサンプリングによって実TOAを決定する。特に、処理サブシステムは、ブレードに対するBPSを用いて、ブレードに対する1以上の実TOAを決定する。ステップ208において、ブレード12に対する正規化ΔTOA210を決定してもよい。本明細書で用いる「正規化ΔTOA」という用語は、複数のブレード内のブレードの実TOAに対応する数値であって、複数のブレードに対する実TOAとブレード間隔パラメータとに基づいて決定される数値をいう。正規化ΔTOA210は、例えば、処理サブシステムによって決定してもよい。一実施形態では、正規化ΔTOA210は、ロバスト最小二乗法技術又は加重最小二乗法技術を実TOA206に適用することによって決定される。ロバスト最小二乗法技術を用いて正規化ΔTOA210を決定すると、明確な正規化ステップがなくなる。ロバスト最小二乗法技術を用いて正規化ΔTOAを決定すると、1以上の異常値実TOA(1以上の異常値実TOAが存在するとき)例えば割れたブレードの実TOAに対する極端な感度が小さくなる。一実施形態では、ロバスト最小二乗法技術を用いて正規化ΔTOA210を決定することによって、負荷データに依拠せずに、負荷の影響に対して正規化された正規化ΔTOAが得られる。正規化ΔTOAの決定について、図3及び図4を参照してさらに詳しく説明する。ステップ212において、モード行列M又は複数のモード214を決定してもよい。さらに、ステップ212において、係数行列U又は複数のモード214に対応する複数のブレード係数215を決定してもよい。モード行列M214及び係数行列U215は、例えば、主成分分析技術を正規化ΔTOA210に適用することによって決定される。なお、現在想定される技術は、主成分分析技術を用いて静的たわみ及び/又はモード行列を決定するための一実施形態について説明しているが、他の技術、例えば、特異値分解技術、独立成分分析技術又はそれらの組合せを用いてもよい。特に、モード行列M214を、例えば、正規化ΔTOA210に基づいて決定してもよく、また固有ベクトル行列Vを、正規化ΔTOA210に基づいて決定してもよい。固有ベクトル行列Vの決定について、図5を参照してさらに詳しく説明する。モード行列M214は、ブレード12に対する複数のモード214を表す。一実施形態では、モード行列M214内の各行が、複数のモード214におけるモードを表す。別の実施形態では、モード行列M214内の各列が、複数のモード214におけるモードを表す。一実施形態では、複数のモード214又はモード行列Mを、以下の等式を用いて決定してもよい。
M=X×V (1)
式中、Mはモード行列であり、Xは正規化ΔTOAの行列であり、Vは固有ベクトル行列である。
式中、Mはモード行列であり、Xは正規化ΔTOAの行列であり、Vは固有ベクトル行列である。
さらに、ステップ212において、複数のブレード係数又は係数行列U215を決定してもよい。係数行列U215は、例えば、固有ベクトル行列Vに基づいて決定することができる。係数行列U215は、複数のモード214に対応するブレード12のブレード係数を表す。一実施形態では、モード行列M214内の行が複数のモード214を表すとき、係数行列U215内の列が、複数のモード214に対応するブレード12のブレード係数を表すが、その逆であってもよい。例えば、係数行列U215内の行が、モード行列215内の列によって示されるモードに対応するブレード係数を表すこともあり、その逆であってもよい。例えば、係数行列Uの第1行が、モードに対応するブレード12のブレード係数を表す場合がある。モードは、モード行列M214内の第1列によって示される。図7(a)及び図7(b)はそれぞれ係数行列Uの2つの行の典型的なグラフ表示を示す。
さらに、ステップ216において、複数のモード214の共通モードに対応する係数行列U215内のブレード係数を決定する。本明細書で用いる「共通モード」という用語は、複数のブレードのブレード係数に対応するモードをいう。ブレード係数はある範囲に含まれる。一実施形態では、モードに対応するブレードのブレード係数に基づいて範囲をリアルタイムで決定する。別の実施形態では、共通モード選択閾値に基づいて範囲を決定する。別の実施形態では、例えば、対応する共通モード選択閾値に基づいて、共通モードに対応する係数行列U215内のブレード係数を特定する。本技術の一実施形態による共通モードに対応するブレード係数の特定及び共通モード選択閾値の決定について図6、図7(a)及び図7(b)を参照してさらに詳しく説明する。
さらに、ステップ218において、係数行列U215内の共通モードに対応するブレード係数をゼロに等しくすることによって再構成行列U 1 を生成させる。その後、ステップ220において、再構成行列U1 及びモード行列Mを用いてブレード12に対する静的たわみ222を決定する。静的たわみは、例えば以下の等式(2)を用いて決定することができる。
Y=M×U 1 (2)
式中、Yは静的たわみを表し、Mはモード行列を表し、U1は再構成行列を表す。その後、ステップ224において、静的たわみ222に基づいてブレード12の健全性を分析してもよい。一実施形態では、ブレード12の1以上に対する静的たわみ222の1以上が所定の閾値を超えるとき、ブレード12の1以上における障害、欠陥又は割れを認定してもよい。なお、一実施形態では、モード行列M214を再構成行列U1に乗じることによって、静的たわみ222の共通モードの影響が除去される。特に、モード行列Mと再構成行列とを乗じることによって、正規化ΔTOA210から共通因子の影響を除去して、ブレード12の静的たわみ22が決定される。
式中、Yは静的たわみを表し、Mはモード行列を表し、U1は再構成行列を表す。その後、ステップ224において、静的たわみ222に基づいてブレード12の健全性を分析してもよい。一実施形態では、ブレード12の1以上に対する静的たわみ222の1以上が所定の閾値を超えるとき、ブレード12の1以上における障害、欠陥又は割れを認定してもよい。なお、一実施形態では、モード行列M214を再構成行列U1に乗じることによって、静的たわみ222の共通モードの影響が除去される。特に、モード行列Mと再構成行列とを乗じることによって、正規化ΔTOA210から共通因子の影響を除去して、ブレード12の静的たわみ22が決定される。
次に、図3を参照すると、本技術の実施形態により、ブレード12に対する正規化ΔTOA210を決定するための典型的な方法300を表すフローチャートが示してある。一実施形態では、図2におけるステップ208を図3でさらに詳しく説明する。前述したように、符号206は、ブレード12に対する実際の到着時間(TOA)を表す。ステップ302において、ラインを、ロバスト最小二乗法技術を用いて実TOA206上にフィッティングしてもよい。図4に、ロバスト最小二乗法技術を用いた実TOA上へのラインの典型的なフィッティングを示す。
ステップ304において、ブレード間間隔パラメータ及び負荷パラメータを決定してもよい。実TOA206上にフィッティングされたラインを用いて、ブレード間間隔パラメータと負荷パラメータとを求めることができる。ブレード間間隔パラメータと負荷パラメータとの典型的な決定について、図4を参照して説明する。ステップ306において、ブレード間間隔パラメータ、実TOA206及び負荷パラメータに基づいて、ブレード12の1以上に対する正規化ΔTOA210を決定することができる。特に、対応する実際の到着時間(TOA)、対応するブレード間間隔パラメータ及び負荷パラメータに基づいて、ブレードに対する正規化ΔTOAを決定することができる。一実施形態では、以下の等式(3)を用いて、正規化ΔTOAを決定することができる。
Norm.Δθ j (k)=θ A (k)−[λ(k)(j−1)+ζ(k)] (3)
式中、Δθjは、ブレードjに対する正規化ΔTOAであり、θ Aは、ブレードjに対する実際の到着時間であり、jはブレードの識別番号であり、λはブレード間間隔パラメータであり、ζは負荷パラメータであり、kはタイムスタンプである。
式中、Δθjは、ブレードjに対する正規化ΔTOAであり、θ Aは、ブレードjに対する実際の到着時間であり、jはブレードの識別番号であり、λはブレード間間隔パラメータであり、ζは負荷パラメータであり、kはタイムスタンプである。
次に、図4を参照して、本技術の一実施形態によりロバスト最小二乗法技術を用いてライン406上にフィッティングされた実TOA401の典型的なグラフ表示400を示す。図4に示すように、X軸402はブレードの識別番号を表し、Y軸404は、ブレードに対する実際の到着時間(TOA)を表す。実TOAは、例えば、実TOA206であってもよい。現在想定される構成で示すように、ライン406を実TOA401上に、ロバスト最小二乗法技術を用いてフィッティングする。さらに、図4に示すように、符号408はライン406の切片である。現在想定される構成では、切片408は、ブレードに対する負荷パラメータである。さらに、符号410は、ブレードに対するブレード間間隔パラメータである。
次に、図5を参照して、本技術の実施形態により、図2におけるモード又はモード行列214とブレード係数又は係数行列U215とを決定するための典型的な方法500を表すフローチャートを示す。一実施形態では、図2におけるステップ212を図5でさらに詳しく説明する。図2を参照して既に説明したように、符号210は、ブレード12に対する正規化ΔTOAを表す。一実施形態では、ステップ502において、正規化ΔTOA210を用いて、正規化ΔTOA行列Xを決定してもよい。ステップ504において、正規化ΔTOA行列Xの共分散を決定して、共分散行列を生成してもよい。さらに、ステップ506において、固有ベクトル行列Vを共分散行列に対して決定する。
ステップ508において、正規化ΔTOA行列Xと固有ベクトル行列Vとに基づいて、モード行列M214を決定する。一実施形態では、行列Xと固有ベクトル行列Vとを乗じることによってモード行列Mを決定する。モード行列M214が決定される結果、複数のモード214が決定される。特に、モード行列214内の各列は、複数のモード214におけるモードを表す。さらに、ステップ510において、固有ベクトル行列Vに基づいて係数行列U215を決定する。係数行列U215は、例えば、固有ベクトル行列Vの逆行列である。なお、現在想定される構成では、係数行列U215は、複数のモード214に対応する複数のブレード係数215を表す。特に、係数行列U215内の各行又は各列のいずれかが、複数のモード214におけるモードに対応するブレード係数を表す。なお、現在想定される構成では、複数のモード214と複数の係数215とを、正規化ΔTOA210を行列として表現及び使用することによって決定するが、ある実施形態では、正規化ΔTOA210の他の表現を用いて、複数のモード214と複数の係数215とを決定してもよい。
次に、図6を参照すると、本技術の実施形態により、モード214の共通モードに対応するブレード係数を決定するための典型的な方法600を表すフローチャートが示してあ。一実施形態では、図2におけるステップ216を図6でさらに詳しく説明する。ステップ602において、1以上の共通モード選択閾値を決定してもよい。本明細書で用いる「共通モード選択閾値」という用語は、モードが共通モードであるか否かを判定するために、モードに対応するブレード係数に基づいて決定される数値をいう。一実施形態では、1以上の共通モード選択閾値は、モードに対応するブレード係数の標準偏差又は中央値を計算することによって決定することができる。
共通モード選択閾値を決定した後に、ステップ604において、複数のブレード係数215の共通モードに対応するブレード係数を特定してもよい。共通モードに対応するブレード係数は、例えば、共通モード選択閾値を用いて特定することができる。図7(a)及び図7(b)に、共通モードに対応するブレード係数を特定することを説明するための複数のブレード係数の典型的なグラフ表示700を示す。特に、本技術の一実施形態により、図6におけるステップ602及び604を図7(a)及び7(b)を用いてさらに詳しく説明する。現在想定される構成では、バープロット702は、モードAに対応するブレード係数を表し、バープロット704は、モードBに対応するブレード係数を表す。一実施形態では、モードAに対応するブレード係数702及びモードBに対応するブレード係数704は、複数のブレード係数215(図2及び図5参照)であってもよい。特に、モードAに対応するブレード係数は、複数のブレード係数又は係数行列215内の行であってもよく、モードBに対応するブレード係数704は、複数のブレード係数又は係数行列215内の別の行であってもよい(図2参照)。
図7(a)及び図7(b)のそれぞれにおいて、X軸706は34個のブレードのブレード番号を表し、Y軸708は、対応するモード702、704における34個のブレードのブレード係数を表す。したがって、現在想定される構成では、バープロット702において、X軸706は34個のブレードのブレード番号を表し、Y軸708は、モードAに対応する34個のブレードのブレード係数を表す。同様に、バープロット704において、X軸706は34個のブレードのブレード番号を表し、Y軸708は、モードBに対応する34個のブレードのブレード係数を表す。例えば、バープロット702におけるバー710は、モードAに対応する第1のブレードのブレード係数を表す。
一実施形態では、共通モードに対応するブレード係数を、1以上の共通モード選択閾値に基づいて特定する。現在想定される構成では、モードA702に対応する34個のブレードのブレード係数には、2つの共通モード選択閾値712、714がある。さらに、モードBに対応するブレード係数704には、2つの共通モード選択閾値716、718がある。一実施形態では、共通モード選択閾値712、714は、例えば、モードAに対応する34個のブレードのブレード係数の標準偏差又は中央値を計算することによって決定される。同様に、共通モード選択閾値716、718は、モードBに対応する34個のブレードのブレード係数の標準偏差又は中央値を計算することによって決定される。
さらに、一実施形態では、モードに対応するブレード係数を、対応する1以上の共通モード選択閾値と対比することによって、共通モードに対応するブレード係数を特定する。したがって、共通モードに対応するブレード係数を決定するために、モードAに対応するブレード係数702を共通モード選択閾値712、714と対比する。同様に、モードBに対応するブレード係数704を、対応する共通モード選択閾値716、718と対比する。
図7(a)に示すように、モードAにおける34個のブレードのブレード係数は、共通モード選択閾値712、714を超えていない。したがって、モードAに対応する各ブレード係数702は、共通モードに対応するブレード係数として特定される。すなわち、モードAは共通モードであると認定できる。モードAは共通モードであるので、モードA702に対応する34個のブレードのブレード係数をゼロに等しくして、再構成行列(例えば、図2におけるステップ218において参照される再構成行列U1)内の行又は列を生成する。しかし、図7(b)では、モードBに対応する17番目のブレードのブレード係数720が、対応する共通モード選択閾値718を超えている。したがって、一実施形態では、モードBは共通モードではないと認定できる。一実施形態では、モードBは共通モードではないと認定した場合、モードB704に対応する34個のブレードのブレード係数を保持して、ゼロには等しくしない。したがって、このような実施形態では、モードBに対応するブレード係数704は、図2におけるステップ218において参照される再構成行列U1内の行又は列を形成する。別の実施形態では、モードB720に対応する17番目のブレード以外のすべてのブレードのブレード係数704を、ゼロに等しくしてもよい。したがって、この実施形態では、ブレード係数720を維持して、再構成行列U1内の行又は列を生成する。
図8は、本技術の一実施形態により、複数のブレードの静的たわみを表す信号、正規化ΔTOAを表す信号、共通モードを表す信号のグラフ表示である。図8に示すように、グラフ802は、30個のブレードに対する正規化ΔTOAを表す信号801のグラフ表示である。グラフ802におけるX軸804はタイムスタンプを表し、グラフ802におけるY軸806は、正規化ΔTOAを表す。グラフ802に示す各ライン801は、正規化ΔTOAを表す信号を表す。グラフ802は、30個のブレードに対する正規化ΔTOAを表す信号を、複数のタイムスタンプにおいて示している。グラフ802における正規化ΔTOA801は、例えば、正規化ΔTOA210(図2及び図3参照)であってもよい。
一実施形態では、グラフ802における正規化ΔTOA801を表す信号を処理して、グラフ808に示す複数のモードの共通モード810、812を表す信号を特定する。共通モード810、812を表す信号は、例えば、図2〜図7を参照して説明した技術を用いて特定してもよい。ステップ218及び220において既に説明したように、正規化ΔTOAから共通モードの影響を取り除くことによって複数のブレードに対する静的たわみを決定する。現在想定される構成では、共通モード810、812を表す信号の影響を、正規化ΔTOA801を表す各信号から除去して、グラフ816に示す静的たわみ814を表す信号を発生させる。共通モード810、812を表す信号の影響を除去することによって、正規化ΔTOA801を表す信号から共通因子の影響が除去されて、信号814が生成される。共通モード810、812を表す信号の影響の除去は、例えば、図2のステップ218を参照して説明した技術を用いることによって行ってもよい。現在想定される構成では、静的たわみ814を表す信号は、30個のブレード内のどんな欠陥及び障害も示してはいない。現在想定される構成では、静的たわみ814を表す信号は、どんな欠陥も示してはいない。なぜならば、30個のブレードすべてに対する信号814が同様の傾向に従っているからである。典型的な目的に対して、図9に、30個のブレードの組に対する静的たわみ902を表す信号のグラフ表示を示して、30個のブレードのうちの1つにおける欠陥を示す。図9において、各ラインは、30個のブレードにおけるブレードに対する静的たわみを表す信号を表す。現在想定される構成で示すように、30個のブレードは、ブレード(以下、ブレードAという)に対する信号904以外は、それぞれ同様の傾向に従っている。言い換えれば、ブレードAに対する静的たわみ904を表す信号は、残りの29個のブレードに対する静的たわみを表す信号と比べて、固有の傾向を示している。ブレードAに対する静的たわみ904を表す信号をさらに抽出及び処理して、欠陥又は割れがブレードA内に存在するか否かを判定してもよい。図2を参照して既に説明したように、ブレードに対する静的たわみを、対応する静的たわみ閾値と対比することによってブレード内の割れ又は欠陥を決定してもよい。
本技術の実施形態によって、ブレードの静的たわみがリアルタイムで決定される。ブレードの静的たわみを、例えば、ブレード内の障害又は欠陥を検出する際に用いてもよい。加えて、本技術では、実TOA又は正規化ΔTOAから共通因子の影響を差し引いて、静的たわみを決定する。なお、本システム及び技術では、共通因子に関する任意のデータ(例えば、動作パラメータ)に依拠せずに、ブレードに対する共通因子の影響を取り除くことによって静的たわみを決定する。本技術によって、静的たわみが、人が介入することなく自動化方式で及びリアルタイムで決定される。ブレードの静的たわみを用いてブレードの健全性を決定する。例えば、静的たわみを用いて、ブレード内の割れ、曲げその他の障害をリアルタイムで特定してもよい。
当然のことながら、必ずしも前述したこのような目的又は利点のすべてが、任意の特定の実施形態により実現されるわけではない場合がある。したがって、例えば、当業者であれば分かるように、本明細書で説明したシステム及び技術の具体化又は実施を、本明細書で教示した1つの優位点又は優位点の群を達成又は最適化するように行うことを、本明細書で教示又は示唆される場合がある他の目的又は優位点を必ずしも実現することなく、実行してもよい。
本発明を限られた数の実施形態に関してのみ詳細に説明してきたが、本発明はこのような開示された実施形態に限定されないことが容易に理解されるはずである。むしろ、これまで説明してはいないが本発明の趣旨及び範囲に見合う任意の数の変形、変更、置換又は均等な配置を取り入れるように、本発明を変更することができる。さらに加えて、本発明の種々の実施形態について説明してきたが、本発明の態様には、説明した実施形態の一部のみが含まれる場合があることを理解されたい。したがって本発明は、前述の説明によって限定されると考えるべきではなく、添付の請求項の範囲のみによって限定される。
Claims (18)
- 複数のブレードに対応する規格化デルタ到着時間を、前記複数のブレードに対応する実際の到着時間に基づいて決定するステップと、
前記複数のブレードの静的たわみを、1または複数の共通因子の影響を前記複数のブレードに対応する前記規格化デルタ到着時間から取り除くことによって決定するステップと、
を含み、
前記複数のブレードの前記静的たわみを決定するステップが、
複数のモードと複数のブレード係数とを、前記規格化デルタ到着時間に基づいて決定するステップと、
前記複数のブレード係数のうち前記複数のモードにおける共通モードに対応するブレード係数を、前記複数のブレード係数に基づいて特定するステップと、
再構成行列を、前記共通モードに対応する前記ブレード係数をゼロに等しくすることによって生成するステップと、
前記複数のブレードの前記静的たわみを、前記複数のモードと前記再構成行列とに基づいて決定するステップと、
を含む、
方法。 - 前記規格化デルタ到着時間を決定するステップが、
ブレード間間隔パラメータと負荷パラメータとを、前記実際の到着時間に基づいて決定するステップと、
前記規格化デルタ到着時間を、前記実際の到着時間、前記ブレード間間隔パラメータ、および前記負荷パラメータに基づいて決定するステップと、
を含む、請求項1に記載の方法。 - 前記複数のブレードの前記静的たわみを決定するステップが、前記1または複数の共通因子の前記影響を、前記1または複数の共通因子に関するデータを用いることなく取り除くステップを含む、請求項1に記載の方法。
- 前記複数のブレードの前記静的たわみを分析して前記複数のブレードの調子を決定するステップをさらに含む、請求項1に記載の方法。
- 前記規格化デルタ到着時間を決定するステップが、ロバスト最小二乗法技術、加重最小二乗法技術、またはそれらの組み合わせを前記実際の到着時間に適用するステップを含む、請求項1に記載の方法。
- 前記複数のモードを決定するステップが、
共分散行列を、前記規格化デルタ到着時間に基づいて決定するステップと、
前記共分散行列に対する固有ベクトル行列を決定するステップと、
前記複数のモードを、前記規格化デルタ到着時間と前記固有ベクトル行列とに基づいて決定するステップと、
を含む、請求項1に記載の方法。 - 前記複数のブレード係数を決定するステップが、前記固有ベクトル行列の逆行列を決定するステップを含む、請求項6に記載の方法。
- 前記複数のブレード係数のうち前記共通モードに対応するブレード係数を特定するステップが、
前記複数のモードにおけるモードに対応する個別のブレード係数に基づいて1または複数の共通モード選択閾値を決定するステップと、
前記ブレード係数を、前記1または複数の共通モード選択閾値に基づいて特定するステップと、
を含む、請求項1に記載の方法。 - 複数のブレードに対応する規格化デルタ到着時間を、前記複数のブレードに対応する実際の到着時間に基づいて決定し、
前記複数のブレードの静的たわみを、1または複数の共通因子の影響を前記複数のブレードに対応する前記規格化デルタ到着時間から取り除くことによって生成する、
処理サブシステムを含み、
前記生成は、
複数のモードと複数のブレード係数とを、前記規格化デルタ到着時間に基づいて決定し、
前記複数のブレード係数のうち前記複数のモードにおける共通モードに対応するブレード係数を、前記複数のブレード係数に基づいて特定し、
再構成行列を、前記共通モードに対応する前記ブレード係数をゼロに等しくすることによって生成し、
前記複数のブレードの前記静的たわみを、前記複数のモードと前記再構成行列とに基づいて決定する
ことにより、生成される、
システム。 - 前記処理サブシステムは、前記規格化デルタ到着時間を、最小二乗法技術、ロバスト最小二乗法技術、加重最小二乗法技術、またはそれらの組み合わせを前記実際の到着時間に適用することによって決定する、請求項9に記載のシステム。
- 前記1または複数の共通因子には、動作パラメータ、前記複数のブレードの再設置、またはそれらの組み合わせが含まれる、請求項9に記載のシステム。
- 前記処理サブシステムは、前記1または複数の共通因子の前記影響を、前記1または複数の共通因子に関するデータを用いずに取り除く、請求項9に記載のシステム。
- 前記処理サブシステムと動作通信状態にある複数の検知装置をさらに備え、前記複数の検知装置は、前記複数のブレードの前記実際の到着時間を表わすブレード通過信号を発生させる、請求項9に記載のシステム。
- 前記処理サブシステムは、前記1または複数の共通因子の前記影響を前記規格化デルタ到着時間から取り除き、主成分分析技術、特異値分解技術、独立成分分析技術、またはそれらの組み合わせを含む技術を適用することを含む、請求項9に記載のシステム。
- 複数のブレードに対応する複数のモードを、前記複数ブレードに対応する規格化デルタ到着時間に基づいて決定し、
前記複数のモードと前記複数のブレードとに対応する複数のブレード係数を、前記規格化デルタ到着時間に基づいて決定し、
前記複数のブレード係数のうち前記複数のモードにおける共通モードに対応する1または複数のブレード係数を特定し、
前記複数のブレード係数における前記1または複数のブレード係数をゼロに等しくして再構成行列を生成し、
前記複数のブレードに対応する静的たわみを、前記規格化デルタ到着時間と前記複数のモードとに基づいて決定する
ように構成された処理サブシステムを含む、システム。 - 前記処理サブシステムはさらに、前記複数のブレードに対応する前記規格化デルタ到着時間を、ロバスト最小二乗法技術を実際の到着時間に適用することによって決定するように構成された、請求項15に記載のシステム。
- 前記処理サブシステムはさらに、前記複数のブレードに対応する前記規格化デルタ到着時間の決定を、
ブレード間間隔パラメータと負荷パラメータとを、前記実際の到着時間に基づいて決定し、
前記規格化デルタ到着時間を、前記実際の到着時間、前記ブレード間間隔パラメータ、および前記負荷パラメータに基づいて決定する
ことによって行なうように構成された、請求項15に記載のシステム。 - 前記処理サブシステムは、前記複数のブレードに対応する前記複数のモードの決定を、
規格化デルタ到着時間行列を、前記複数のブレードに対応する前記規格化デルタ到着時間に基づいて決定し、
共分散行列を、前記規格化デルタ到着時間行列の共分散を決定することによって決定し、
前記共分散行列に対する固有ベクトル行列を決定し、
前記複数のモードを、前記規格化デルタ到着時間行列と前記固有ベクトル行列とに基づいて決定する、
ことによって行なうように構成された、請求項15に記載のシステム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IN5014CH2012 | 2012-11-30 | ||
IN5014/CHE/2012 | 2012-11-30 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2014114806A JP2014114806A (ja) | 2014-06-26 |
JP2014114806A5 true JP2014114806A5 (ja) | 2018-01-25 |
JP6280354B2 JP6280354B2 (ja) | 2018-02-14 |
Family
ID=49626878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013245535A Active JP6280354B2 (ja) | 2012-11-30 | 2013-11-28 | 翼形部の健全性をモニタするためのシステム及び方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9404386B2 (ja) |
EP (1) | EP2749740B1 (ja) |
JP (1) | JP6280354B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111353129B (zh) * | 2020-02-10 | 2021-12-28 | 西安交通大学 | 叶端定时数据存储矩阵化处理方法 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9618096D0 (en) | 1996-08-29 | 1996-10-09 | Rolls Royce Plc | Identification of resonant frequencies of vibration of rotating blades |
US6594619B1 (en) | 1999-08-02 | 2003-07-15 | Hood Technology Corporation | Apparatus and method for predicting failures of spinning disks in turbo-machinery |
US7043348B2 (en) | 2003-03-03 | 2006-05-09 | Honeywell International, Inc. | Transient fault detection system and method |
US7693643B2 (en) | 2005-02-14 | 2010-04-06 | Honeywell International Inc. | Fault detection system and method for turbine engine fuel systems |
WO2007047375A1 (en) | 2005-10-14 | 2007-04-26 | Fisher-Rosemount Systems, Inc. | Statistical signatures used with multivariate statistical analysis for fault detection and isolation and abnormal condition prevention in a process |
US7243048B2 (en) | 2005-11-28 | 2007-07-10 | Honeywell International, Inc. | Fault detection system and method using multiway principal component analysis |
US7421351B2 (en) | 2006-12-21 | 2008-09-02 | Honeywell International Inc. | Monitoring and fault detection in dynamic systems |
JP4560110B2 (ja) * | 2008-09-17 | 2010-10-13 | ジヤトコ株式会社 | 異常診断装置及び異常診断方法 |
US8543341B2 (en) * | 2010-06-29 | 2013-09-24 | General Electric Company | System and method for monitoring health of airfoils |
US7941281B2 (en) * | 2008-12-22 | 2011-05-10 | General Electric Company | System and method for rotor blade health monitoring |
US8532939B2 (en) * | 2008-10-31 | 2013-09-10 | General Electric Company | System and method for monitoring health of airfoils |
US8676514B2 (en) * | 2010-06-29 | 2014-03-18 | General Electric Company | System and method for monitoring health of airfoils |
US20110313726A1 (en) | 2009-03-05 | 2011-12-22 | Honeywell International Inc. | Condition-based maintenance system for wind turbines |
US9045999B2 (en) * | 2010-05-28 | 2015-06-02 | General Electric Company | Blade monitoring system |
US9103741B2 (en) * | 2010-08-27 | 2015-08-11 | General Electric Company | Methods and systems for assessing residual life of turbomachine airfoils |
US9046000B2 (en) * | 2011-06-18 | 2015-06-02 | Prime Photonics, Lc | Method for detecting foreign object damage in turbomachinery |
US20130082833A1 (en) * | 2011-09-30 | 2013-04-04 | General Electric Company | System and method for monitoring health of airfoils |
-
2013
- 2013-11-26 EP EP13194494.4A patent/EP2749740B1/en active Active
- 2013-11-26 US US14/089,995 patent/US9404386B2/en active Active
- 2013-11-28 JP JP2013245535A patent/JP6280354B2/ja active Active
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Madhavan et al. | Vibration based damage detection of rotor blades in a gas turbine engine | |
KR101718251B1 (ko) | 회전 블레이드 강건성 모니터링 방법 및 시스템 | |
Chen et al. | New step to improve the accuracy of blade tip timing method without once per revolution | |
JP6302152B2 (ja) | エーロフォイルの健全性を監視するためのシステムおよび方法 | |
Megerle et al. | Numerical and experimental investigation of the aerodynamic excitation of a model low-pressure steam turbine stage operating under low volume flow | |
JP5898865B2 (ja) | エアフォイルの健全性を監視するシステム及び方法 | |
JP5879055B2 (ja) | エアフォイルの健全性を監視するシステム及び方法 | |
JP2010144727A (ja) | ロータブレード健全性を監視するためのシステム及び方法 | |
US9250153B2 (en) | Methods and systems for monitoring health of blades | |
JP2013231721A (ja) | ステータベーンの健全性を監視するためのシステムおよび方法 | |
Krause et al. | Asynchronous response analysis of non-contact vibration measurements on compressor rotor blades | |
Wang et al. | Time-dependent vibration frequency reliability analysis of blade vibration of compressor wheel of turbocharger for vehicle application | |
JP6280354B2 (ja) | 翼形部の健全性をモニタするためのシステム及び方法 | |
Kielb et al. | Experimental and numerical study of forced response in a full-scale rotating turbine | |
JP2014114806A5 (ja) | ||
Grądzki et al. | Rotor blades diagnosis method based on differences in phase shifts | |
Bhattacharya et al. | System to monitor blade health in axial flow compressors | |
Rao et al. | In situ detection of turbine blade vibration and prevention | |
Hanachi et al. | Bladed disc crack diagnostics using blade passage signals | |
Rajagopalan et al. | Estimation of static deflection under operational conditions for blade health monitoring | |
Rozman et al. | Characterizing Flow Instabilities During Transient Events in the Turbine Rim Seal Cavity | |
Przysowa et al. | Health monitoring of turbomachinery based on blade tip-timing and tip-clearance | |
Civardi | Development and validation of algorithms for vibration monitoring of aeronautical engines | |
Kotambkar et al. | Analytical and numerical investigation of lacing wire damage induced mistuning in turbine blade packet | |
Ferrari | Experimental characterization of the stable operating range of a highly-loaded low-pressure compressor stage |