JPS5822926A - 回転軸系捩り振動監視装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 生した捩り振動を少数の測定可能な位置で測定し，この
振動が線型に分解できるとすることにより回転軸糸上の
任意位置の捩り振動を推定する回転軸系捩り振動監視装
置に関するものである。

擾乱の実態を把握してお（ことは設計の技術上非常に重
要なことであり，特に実際に使用している者にとっては
その擾乱により刻々と変化する疲労寿命を知ることが重
要なことである。

しかしタービン発電機の回転軸糸は長さ数十メートルに
及ぶ場合もあり測定点は数箇所必要となるのが普通であ
る。

また回転軸糸に生じる捩り振動はその回転軸に疲労を及
ぼす可能性があり．これを測定することは重要なことで
あるが回転軸の数箇所に捩り振動の測定装置を設置する
ことは極めて不経済でもあり，物理的に設置不可能な所
もある。

振動学士の知見によれば回転軸糸に生じる捩り振動ｙｒ
ｘ．ｔ）はｔ−ド別振動Ｙｉｒｘ．ｔ）の和で表わされ
，同モード別振動は振動ｔード形ＧＬ（ｘ）と振動モー
ド成分Ｈ’ｉ　（　ｔ　）の積で表わされるとみなすこ
とができる。

Ｙ（ｘ，ｔ）＝ΣＹｉ（　ｘ，　ｔ　）　　　　・−・
・ｉｌｌ＝ΣＧｉｒｘ）・Ｈｉｒｔ）　　　−＝　＋２
１ここで振動七一ド形ＧｉＬｘ）は第１図に示すよ５な
軸糸によりあらかじめ定まる波形であってｉ−４ｏ３１
　Ｇｉ（ｘ）　ｌ　＝０ となる。

そこで、回転軸系．の定位置Ｐ（座標Ｘｐ）において、
適宜な捩り振動検出装置で振動 Ｙ（Ｘｐ．ｔ）＝ΣＹｉ　（Ｘｐ．　ｔ　）　＝ΣＧｉ
　（Ｘｐ）　ａＨｉ（ｔ）　　−−　　［３１を検出す
る。（第２Ａ図に示す。） 次に適当なフィルタを用いることにより第２Ｂ図に示す
ようなモード別振動 Ｙｌ　（Ｘｐ，　ｔ）　＝Ｇｉ　（Ｘｐ）・ｉｔ（ｔ）
　　　　　　　　　・・・・・・　（４）が求められる
。

ここで振動ｔード形Ｇｌ（Ｘｐ＞は第１図に示したよう
に，座標Ｘによって値が一意的に定まるものである。

従って軸糸の任意の点，たとえば捩り振動を監視したい
位置Ｊの座標なＸｊとすわば，点Ｊにおける七一ド別振
動は Ｙｉ　（ＸＪ，　ｔ）　＝Ｇｉ　（Ｘｊ）　・Ｈｉ　（
　ｔ）＝（Ｇｉ（Ｘｊ）／ｌｌｒｉ（Ｘｐ））・Ｇｔ（
ｘｐ）・Ｈｔ　（　ｔ）：＝　（ＧｉｒＸｊ）／Ｄｉｒ
Ｘｐ））Ｙｉ（Ｘｐ，　ｔ）二Ｇｉｊ・Ｙｉ（Ｘｐ，　
ｔ）　　　　　　　　・・・・・・（５）と表わさｎる
。Ｃ第２Ｃ図〕 すなわちあらかじめ各七ード毎における定位置Ｐと任意
位置Ｊとの換算係数Ｇｉｊ　（　Ｊは１つとは限らない
）を求めておくことによって同Ｊ点の捩り振動は Ｙ（ＸＪ．　ｔ）＝ΣｙＢｘｊ，ｔ） ＝ΣＧｔ　ｊ　ａＹｉ（　Ｘｐ．　ｔ　）　　　・−−
・・（６１として．定位置Ｐの七一ド別振動から推定す
ることができる。（第２Ｄ図） しかし、このフィルタを用いる方法は．フィルタの過渡
応答特性のため振動振幅が急激に変化する場合などでは
．応答お（れを生じて精度よ（監視点における振動を推
定できないことがある／ そ　で大きなｉに対しＧｌ（ｘ）が近似的に零に等Ｌ％
いことに着目して次のような方法が考えら／ れ／た。

／ ′上記（２１式において．適当なｎを定めることにより ｙＣｘ，ｔ）＝ΣＧｉ（ｘ）ｅＨｌ（ｔ）　　　　　　
　　−・・・・・ｔ７１とみなすことができる。

従って、軸糸上ｎ個の定位ｆｔＰｋ（座標Ｘｐｋｋ−１
・・・・・・ｎ）で測定されるｎ組の捩り振動はとして
表わされるはずである。

すでに述べたようにＧｉ　（Ｘｐｋ　）は一定の値であ
るから上記ｎ組の式は撮動［−ド成分Ｈｔｒｔ）を未知
数とするｎ元１次の連立方程式である。

従ってｎＸｎ行列の要素Ｇｉ（Ｘｐｋ）　（ｉが列。

ｋが行を示す）からなる行列の逆行列の要素０１ｋ（ｌ
は行、ｋが列を示す）はあらかじめ求めてお（ことがで
きる。

よって１次の撮動ｔ−ド成分は Ｈｉ（ｔ）＝ΣＧ’ｉｋ＠Ｙ（Ｘｐｋ、　ｔ）　　　　
・・曲・・・　（８ンと表わすことができるので任意位
置Ｊの捩り振動は Ｙ（Ｘｊ、　ｔ）＝ΣＧｉ　（Ｘｊ　）・Ｈｉ　（ｔ　
）＝ΣＧｉ　（Ｘｊ　）　（ΣＧ’１ｋ−Ｙ（Ｘｐｋ、
ｔ））　　・・・・・・（９）で推定されることになる
。

ただしこの方法では回転軸系の捩り振動をｎ次までの七
−ド別振動の和であ・るとみなすと。

それと同数すなわちｎ個の点において軸糸に生じる捩り
振動を検出しな（てはならないことになる。

しかし一般に回転軸糸の捩り振動をピ・シクア・シづを
用いて測定する場合においては、その取付場所に物理的
な制限がある。従って、定位置とすべき箇所が少ない方
が良いものの１反面監視点の振動を推定する場合には出
来るだけ多くの七−ド別振動まで考慮した方が推定精度
が向上するわけである。

またｎ次までの振動ｔ−トド成分求める際に１次以上の
七−ド別振動を含んでいる捩り振動を用いているので振
動モード成分自体にノイズが含まれてしまうことになる
。

従ってこれらフィルタを用いる方法、連立方程式を解く
方法の夫々の特徴を生かすべくなされたものとして回転
軸糸のｎ個の定位置で同軸糸の捩り振動を検出するｎ組
の検出装置と、同装置により得られた各々の捩り振動を
上記回転軸糸固有のｎ＋１次以上の各モード別振動に分
解する分解装置と、上記検出装置で得られた各定位置の
振動から上記分解装置で得られた各定位置のｎ　＋　１
次以上のｔ−ド別振動の和を減算する減算装置と、同装
置で得られたｎ個の定位置の夫々ｎ次までの和の振動に
あらかじめ求められているｎ個ｎ組の定数を各側割に乗
じて和を取りｎ次までの振動モード成分を求める行列演
算装置と、上記定位置のうちの１つの位置のｎ　＋１次
以上のモード別振動を同位置の同次の撮動ｔ−ド形で除
するとともに上記回転軸上の任意位置の同次の振動モー
ド形を乗じ同任意位置のｎ　＋　１次以上の各々のｔ−
ド別振動を求めるＭｌのモード別演算装置と、上記行列
演算装置で得られたｎ次までの振動モード成分に上記任
意位置の同次の振動ｔ−ド形を乗じ同位置のｎ次までの
ｔ−ド別振動を求める第２のｔ−ド別演算装置と、これ
らＭｌ、第２のｔ−ド別演算装置で得られたｔ−ド別振
動を加算し上記任意位置の振り振動を求める加算装置と
からなるものが提供されている。

すなわちｎ個の定位置Ｐｋ（座標Ｘｐｋ　）で測定した
捩り振動は Ｙ　（Ｘｐｋ、　ｔ　）　＝　Ｘ：　Ｙｉ　（Ｘｐｋ、
　ｔ　）　　　　・・−−ｔｌ（１である。

そこでｎ　千１次からｎ十ｍ次までの夫々のモード別振
動を摘出するフィルタにより該モード別振動Ｙｉ　（Ｘ
ｐｋ、　ｔ）　（ただしｎ　＜　ｉ　＜、ｎ　＋　ｍ　
）を摘出するとともに（１（１式で示される捩り振動よ
り除去する。

それをＦ　（Ｘｐｋ、　ｔ　）で定義すると＝ΣＹｉ　
（Ｘｐｋ、　ｔ　）＋ｇ　（ｘ、　ｔ　）となりε（ｘ
、ｔ）はｎ　＋　ｍ次以上のｔ−ド別振動の和であり０
とみなし、あらためてｐ　（Ｘｐｋ、　ｔＪ　”Σｙｔ
　（Ｘｐｋ、　ｔ　）二ΣＧｉ　（Ｘｐｋ　）　ａＨ（
ｔ　）　　　　　　−−・・　　αＯとおく。

従って６υ式は１元１次の連立方程式となり。

あらかじめｎＸｎの要素Ｇｉ（Ｘｐｋ）からなるマトリ
・シクスの逆行列の！！素Ｇ′庚を求めてお（ことによ
り（８２式と同様 Ｈｌ（ｔ）＝ΣＧ’ｉｋ＠Ｆ（Ｘｐｋ、　ｔ）　　　　
　　・・−・−（１１）’が得られ任意位＆Ｊのｎ次ま
での七−ド別振動Ｙｉ（Ｘｊ、ｔ）＝Ｇｉ（Ｘｊ）”（
ＪＧ’１ｋ−Ｆ（ｐｋ、ｔ））　　・＝−””−Ｑ３と
して求めることができる。

一方フィルタで摘出されたｎ個の定位置のｎ＋１次以上
のｔ−ド別振動を各法において１つづつ取り出す。

たとえばそれをＹｉ　（Ｘｐｋ、　ｔ　）とし、上記（
５）式でｋｌ 示した位置ＰｋとＪとの換算係数Ｇｉｊ　をあらかじめ
求めておくことによりｎ　＋　１次以上の任意５点にお
けるモード別−動は ｋｌ Ｙｌ　ｒｘｊ、　ｔ）＝Ｇｉｊ　１ＩＹｉ　（Ｘｐｋ、
　ｔ　）　　　　　　・・・・・・・・・　０と表わさ
れる。〔ただしとのｋを選択するかは全（自由である。

） 上記αり、０式で夫々表わされるモード別撮動を加算す
ることにより任意の点Ｊでの捩り振動は として推定することができる。

このようにすることによって フィルタ方式の欠点であるフィルタ過渡応答がおくれる
点、および 連立方程式方式の欠点である多（の計測点を必要とする
点や、高次の七−ド別振動などによる誤差を生じる点 などの欠点を解消できるわけである。

しかしこのような監視装置では、検出装置で得られた各
々の捩り振動を用いて常時演算を行ない回転軸系を監視
していることになる。それが為に送電系統に係わる各種
擾乱によるいわゆる必要な振動以外の信号、たとえば入
力系にノイズが混入した場合にも見掛上振動が生じたも
っと見なされてしまい、それを用いて演算を行なうと誤
まった結果が得られることになる。

本発明はこれらの欠点を排除するものであって、上述の
装置に次の各装置を付加した。

すなわち、検出装置に得られた各々の捩り振動の信号を
待機させるシフトレジスタと、上記捩り振動を各モード
別撮動に分解する帯域フィルタと：同フィルタで分解さ
れた各モード別撮動の一定時間内の二乗の平均値を求め
る平均演算装置と、同装置で得られた平均値と所定の値
とを比較する比較器と、上記平均値が所定の値より大き
い時に作動し上記シフトレジスタで待機させられた捩り
撮動の信号を分解装置に入力させる切換器とを具備する
回転軸系捩り振動監視製蓋である。

従って、適宜位置に配設された少数の検出装置で得られ
た捩り振動の信号を、送電系統に加わる高速三相再閉路
などの各種擾乱によるもので、装置の誤作動を防止でき
る。

以下不発明を第６因ないし第５図に示す−実施例につい
て説明する。

この実施例ではｅ　　ｎ　””　２　＋　　ｒｎ　”　
８とし、また推定された任意位置での捩り振動を用いて
同位置の疲労消費寿命を求める為の装置および本発明の
装置では推定することが好ましくないような異常な振動
が発生した場合の監視装置をも付加しである。

そこでまず疲労消費寿命の求め方について説明する。

回転軸糸に生じる捩り振動による応力は、その振動の振
幅に比例するから点ＸＪでの比例定数をｏｊとすれば、
上述の製電により得られた点Ｘｊでの捩り振動ｖｒｘｊ
、ｔ）により同点Ｘｊの応力σｊは ｏｊ　＝ａｊ１１Ｙ　（ｘｊ、　ｔ　）＝Σαｊｃｔｔ
ｘｊ）・Ｈｌ　（ｔ　）　　　　　・・・・曲・　　（
１４９（第２Ｅ図に示す） で表わされる。

さて１回転軸糸の任意の点ｊでの応力σｊが判れば、既
知の疲労寿命推定法に基づいて材料のＳ−Ｎ線図（応力
振巾−疲労（り返し数線図）を用いて計算することが可
能であり９本実施例ではそれをレンジ法により行なう。

ここでし：、Ｉじ法について説明する。第２Ｅｔｇに示
すよ５な０！９式によるｊ点の応力が得られたとする。

上記０！９式における極値をｔの小さい方から順番にｏ
ｊ　ｌｔｌ　、　　ｏｊ（２）・・・・・・とすると、
に番目の極値か１ｋ）　　　（ｋ） Δａｊ＝’１２１ａｊ−σｊ（ｋｌ１）Ｉ　　　　　＋
＋＋＋＋＋＋　　［１６１−（ｋ）　　　（ｋ） σｊ−怪１σｊ　＋σｊ（ｋｌｔ　）　ｌ　　　　　・
・曲・・・　ａηで表わすことができる。

また平均応力が零の場合と等価な変動応力ｔｋ）　　　
　　　　　　　　　　　　ＢΔσｅｑｊはｊ点の回転軸
の引張強さをｏｊ　　とすると、修正Ｇｏｏｄｍａｎ　
１１１図の考えにより。

６４Ｆｅｑｊｋ）”−６ｊＢ”　、、、、６ｊ（ｋす（
ａＪＢ−ａｌ”）　　　＝−００８で表わせる。

このようにして回転軸系の任意の位置ｊの応力波形の極
値に番目とに＋１番目との間の応力缶） の等価変動応力Δσ。、ｉが判ると、ｊ点の材料の変動
応力と繰り返し数との関係を示す５−ＮＩｉＩ図（第３
図）から、このｊ点に等価変動応力（ｋ） での繰り返し数Ｎｊ　　が判る。

（ｋｌ 寿命の消貢量へＤｊは、応力の変化が応力波のサイクル
の半分であるとみなして で表わされる。

従って回転軸糸の任意点ｊでの変動応力を０１式により
演算しその極値をカウントし、ａｅ式な（ｋｌ いし０式により、寿命の消費量ムＤｊを算出し累積する
。

すなわち は１回転軸糸の任意の位置ｊ点で応力の極値がｎ　＋　
１回生じた時点までの寿命の消費量を表わしていること
になる。

さて第４図において１は発電機、２および３はそれぞれ
発電機１を駆動する低圧タービシおよび高圧タービンで
ある。４はこれらを連結する回転軸であり一直線上に配
置されて回転軸系２を構成する。５は上記回転軸４の２
箇所定ケ置ＰＫ　のうちの１つＰＩ（座標゛を均、とす
る）に取り付（すられたターニング＋セ、６は同ター二
シク＋セを利用して上記回転軸４の捩り振動を検出する
ピックアップであり、検出装置７を構成する。

もちろん定位１１Ｐ１点はここに限らず９回転−軸４の
どこにあってもよい。

人魚ないしＨ点は上記回転軸４の捩り振動を監視したい
任意の位置で夫々座標ＸＡないし別を有′するものであ
り、一般にＸｊとして表わされる。

従ってＸｊの値は任意に取れるので監視位置は自由であ
る。

８は低域フィルタであって、軸糸２がどこでも同一に振
動する剛体ｔ−ドε′ｒｘｐｋ、ｔ）を分離する約１〜
２１（ｚ以下の成分しか通さない特性を持ｌ）。

９、ないし籟　Ｃ代表して９ｉとする）は１次の七−ド
別振動Ｙｌ　（Ｘｐｋ、　ｔ　）を分離する帯域フィル
タからなる分解装置である。

１０は減算装置で、上記検出装置７で検出された捩り振
動ｙｒｘｐｋ、ｔ）から低域フィルタ８゜分解装置９１
の出力ε’（Ｘｐｋ、　ｔ）、　ｙｔ　（Ｘｐｋ、　ｔ
　）を減算し００式で示されるＦ（Ｘｐｋ、　ｔ）’を
求めるものである。

１１は行列演算装置であって、上記減算装置１０で得ら
れたＦ（Ｘｐｌ、　ｔ）、　’Ｆ（Ｘｐ！、　ｉ）　　
に２×２の逆行列要素の２個２組の値’　”Ｉ＋　＋　
”＋２　）＋　（Ｇ’２１＋ＧＺ□）を各側割に乗じて
和を取り上記（１１つの式位置Ｘｊの四次の振動モード
形Ｇｉ　（Ｘｊ）を乗じて同位置の２次までのｔ−ド別
振動Ｙｉ　（Ｘｊ、　ｔ　）　（１２式）を求める第２
のｔ−ド別演算装置ひある。

１３は第１のモード別演算装置であって、第３次、第４
次に対しては定位置ｘｐ、　、第５次ないし第１０次に
対して定位置Ｘ％　の捩り振動を分解装置９１で分解し
得られた各ｔ−ド別振動Ｙｉ（ｋ） （Ｘｐｋ、ｔ）に換算係数Ｇｉｊを乗じて、３点のｔ−
ド別振動Ｙｉ（Ｘｊ、　ｔ）　（１３式）を求めるもの
である。

１４は上記第１１第２のｔ−ド別演算装置１３゜１２で
得られた任意位置Ｘｊ（ｊ：Ａ・・・・・・Ｈ）のｔ−
ド別振動Ｙ１（Ｘｊ、ｔ）を加算し上記任意位置の捩９
撮動Ｙ（Ｘｊ、　ｔ）（１４式）を求めるものである。

なお、この実施例では比例定数αｊを乗じて崗式で示す
応力σｊを求める。

１５は応力σｊの極値σ７”（ｋ＝１．２．・・・・・
・・・・）を検出する極値検出器であって、先に示した
（１６１ないし翰式で示された疲労寿命消費を算出する
たンごユータ１７に送る切換器、１８および１９は上記
］ンじユータ１７で算出された各点における消費寿命を
指示、記録する指示器および記録器である。従ってｊ点
の消費寿命ＤＪが１になると寿命が尽きたことが判断で
きる。

２０は電力などの回転軸糸２に加わる外力を検出する外
力検出装置、２１は記録装置２２がスタートしてから回
転が整定し正常の記録が出来るようになる迄の時間を補
正するための遅延回路、２６は設定値σＳに応じた電圧
を発生する設定電圧発生器、２４は加算゛装置１４で得
られた任意の位置の応力σｊと上記設定値σ８とを比較
し、σ８〈σｊの時上記記録装置２２をスタートさせる
比較器である。

記録装置２２は上記外力検出装置２０および検出装置７
の出力を記録するものである。

従って、任意の位置ｊ点における応力が太き（なりすぎ
た場合には、定位置ＰＫ点の捩り振動Ｙ（Ｘｐｋ、　ｔ
）　、およびその時に回転軸２に生じた外力を記録装置
２２に記録させる。

２５は帯域フィルタであって、上記検出装置７で得られ
た捩り振動ｙｒｘｐｋ、ｔ）の信号を各七−ド別振動ｙ
＋ｒｘｐｋ、　ｔ）に分解するものであり。

帯域フィルタ２５の出力は平均演算装置２６に入力され
る。同平均演算装置２６では各七−ドｉ毎に一定時間内
の２乗平均値χｌが求められ。

予め所定の値ｙを設定した設定器２７と対になり比較器
２８で比較される。

一方、２９は検出装置７で得られた捩り振動を待機させ
るシフトレジスタであって、切換器６０の切換によすて
待機させていた捩り振動な分解装置９等に入力すること
が出きる。上記切換器３０は比較器２８の出力によって
作動されるものであって、比較器２８は２乗平均値χｔ
が設定値ｙより大きい時に信号を出し切換器６０を作動
させる。

とは明確に区別できる。

第６図には軸糸に擾乱が生じた場合、第７図にはスパイ
クノイズが混入した場合の波形を示しである。

図（ＡＪは検出装置７で得られる捩り振動の波形であっ
て、帯域フィルタ２５により図（Ｂ）に示すようなｉ次
のｔ−ド別振動に分解され、平均演算器２６で図ｔｃ）
に示す２乗平均値が求められる。

設定値をｙとすると第６図（Ｃ１の時間ｔ、で設定値を
超えるので、第６図ｆＤｌに示すように切換信号が切換
器３０を作動するが、＠７図（Ｃ１に示す場合では設定
値ｙより小さいので切換信号は発生しない。

このように本発明の一実施例では回転軸４に発生した捩
り振動を２箇所の定位置ＰＫにおいて検出装置７で検出
し任意の各点ｊの捩り振動を計算した。更に応力を求め
、任意の点の疲労寿命消費を計算させた。

すなわち回転軸糸２の監視点毎にごツクアラづを取り付
けて応力を求めることにより監視するのではなく、定位
置の監視から任意の位置を推定し監視するようにしたの
で監視点に物理的な制限は何もない。

また外力検出装置２０．比較器２４．記録装置２２によ
り応力が設定値以上すなわち線型において推定すること
が好ましくない様な振動が発生した場合、その生のＰ点
における振動を記録するとともに、その時に回転軸系に
生じた外力を記録するようにしたので後工程において大
型の計算機を使用して十分な解析を行なうことが可能で
ある。

更に帯域フィルタ２５．平均演算装置２６゜設定器２７
．比較器２８および切替器３０により、検出装置７の出
力が回転軸系に加わった各種擾乱によるものか、ノイズ
によるもの力）を判別し、擾乱による場合だけについて
上述の演算を行なうようにした。

従って、ノイズ等による装置の誤動作を防止することが
可能となった。

【図面の簡単な説明】

第１図は回転軸系の振動ｔ−ド形の説明図。 第２図は回転軸系に生じている一各種の波形を示し、第
２Ａ図は定位置Ｐ点の捩り振動波形、第２Ｂ図は同Ｐ点
の波形をｔ−ド毎に分解した波形、第２Ｃ図は第２Ｂ図
の波形により演算される任意位置ｊ点のｔ−ド毎の捩り
振動波形、第２Ｄ図は第２Ｃ図を合成したｊ点の捩り振
動波形、第２Ｅ図は同ｊ点の応力波形、第３図はＳ−Ｎ
線図、第４図は本発明の測定方法をタービン発電機に用
いたー実施例、第５図はそのづ０・シフ線図、第６図お
よび第７図は、帯域フィルタ平均演算装置の作用の説明
図である。 に発電機、２：低圧タービシ、３：高圧タービン、４：
回転軸、５：ターニシグギセ、６：ピ・シフアラづ、７
：検出装置、９１：分解装置。 １０：減算装置、１１：行列演算装置、１２：第２の℃
−ド別演′算装置、１３：第１のｔ−ド別演算装置、１
４二加算装置、１５：極値検出器、１６：切換器、１７
：フンピユータ、１８：指示器、１９：記録器、２０：
外力検出器。 ２１：遅延回路、２２：記録装置、２３：設定電圧発生
器、２４：比較器、２５：帯域フィルタ、２６：平均演
算装置、２７：設定器、２８；比較器、２９：Ｌ／フト
レジスタ、５０：切換器Ｏ 第２Ａ図 第２Ｄ図 第す困 第７図

Claims (1)

1. 【特許請求の範囲】 回転軸糸のｎ個の定位置で同軸糸の捩り振動を検出する
   ｎ組の検出装置と、同−置により得られた各々の捩り振
   動を上記回転軸糸固有のれ＋１次以上の各モード別振動
   に分解する分解装置と、上記検出装置で得られた各定位
   置の振動から上記分解装置で得られた各定位置のｎ　＋
   　１次以上のモード別振動の和を減算する減算装置と、
   同装置で得られたｎ個の定位置の夫々ｎ次までの和の振
   動にあらかじめ求められているｎ２１　ｎ組の定数を各
   側割に乗じて和を取りｎ次までの振動Ｌ−ド成分を求め
   る行列演算装置と。 上記定位置のうちの１つの位置のｎ　＋　１次以上のモ
   ード別振動を同位置の同次の振動モード形で除するとと
   もに上記回転軸上の任意位置の同次の振動モード形を乗
   じ同任意位置のｎ＋１次以上の各々の七−ド別蚕動を求
   めるＭｌのｔ−ド別演算装置と、上記行列演算装置で得
   られたｎ次までの振動上−ド成分に上記任意位置の同次
   の振動モード形を乗じ同位置のｎ次までのモード別振動
   を求める第２のモード別演算装置と。 これら第１．第２のｔ−ド別演算装置で得られた七−ド
   別振動を加算し上記任意位置の捩り振動を求める加算装
   置とからなる回転軸糸捩り振動監視装置において、上記
   検出装置により得られた各々の捩り振動の信号を待機さ
   せるシフトレジスタと、上記捩り振動の信号を各七−ド
   別の撮動に分解する帯域フィルタと、同フィルタで分解
   さ牡た各モード別振動の一定時間内の二乗の平均値を求
   める平均演算装置と、同装置で得られた平均値と所定の
   値とを比較する比較器と、上記平均値が所定の値より大
   きい時に作動し上記シフトレジスタで待機させらｎた挟
   り振動の信号を上記分解装置に入力させる切換器とを具
   備したことを特徴とする回転軸糸捩り振動監視装置。