JPH03229125A - Accuracy improving method for rotary machine state monitoring system - Google Patents
Accuracy improving method for rotary machine state monitoring systemInfo
- Publication number
- JPH03229125A JPH03229125A JP2392290A JP2392290A JPH03229125A JP H03229125 A JPH03229125 A JP H03229125A JP 2392290 A JP2392290 A JP 2392290A JP 2392290 A JP2392290 A JP 2392290A JP H03229125 A JPH03229125 A JP H03229125A
- Authority
- JP
- Japan
- Prior art keywords
- positions
- bearing
- axis
- force
- supporting force
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Description
【発明の詳細な説明】
ア)産業上の利用分野
タービン等の油膜軸受で支持した回転機械の回転軸心の
位置を監視して、回転機械の運転に支障をきたす異常状
態の有無を知るシステムを設計・利用する分野に属する
。Detailed Description of the Invention A) Industrial Application Field A system that monitors the position of the rotational axis of a rotating machine such as a turbine supported by an oil film bearing to determine whether there is an abnormal condition that may impede the operation of the rotating machine. It belongs to the field of designing and using.
イ)従来技術
渦流式変位計利用などのような単なる間隙測定法が用い
られうるだけであったから、例えば日本機械学会第68
1回講習会「やさしい回転体」教材p34表2の記載に
見受けうるように、軸心位置の監視においては、その欠
点として「測定器(特に構出器)の信鎖性が劣る」状態
にあった。従って、熱的・電圧的経時変化等によるドリ
フトによって測定値の誤差が増す結果、回転機械の状態
変化のトレンドの監視に支障を来していた。b) Conventional technology Since only a simple gap measurement method such as the use of an eddy current displacement meter could be used, for example, the Japan Society of Mechanical Engineers 68
As can be seen in the description in Table 2 on page 34 of the one-time training course "Easy Rotating Objects", the disadvantage of monitoring the axis position is that "the reliability of the measuring device (especially the measuring device) is poor". there were. Therefore, as a result of increased errors in measured values due to drift due to thermal and voltage changes over time, etc., it has been difficult to monitor trends in state changes of rotating machines.
つ)発明が解決しようとする課題
回転軸心の位置の計測値に、ドリフト等の時間的変化率
の比較的小さい計測誤差が混入していても簡便に除去し
えて、回転機械の状態変化の長期的トレンドの監視を行
いうる方法の提案を行う。(1) Problem to be Solved by the Invention Even if measurement errors with a relatively small rate of change over time, such as drift, are mixed into the measured value of the position of the rotational axis, it can be easily removed and Suggest ways to monitor long-term trends.
工)課題を解決するための手段
既設の渦流式変位計などを利用した間隙測定値を基にし
て得た軸受中心に対する回転軸心の偏心率(=偏心量を
有効半径クリアランスで除したもの)から、その偏心率
に対する容量係数(=ゾンマーフェルト数の逆数)を用
いて当該軸受の支持力を算出するに際して、鉛直方向の
支持力が最小最大になる回転軸の周方向位置を割出し、
両該位置を基準にして回転軸の周方向位置に対する位相
が同じ所の支持力のベクトル和を位相を変えて少なくと
も2回求め、それらベクトル和間の差を最小にするよう
に間隙測定値を演算補正する。Eccentricity of the rotating shaft relative to the bearing center (= eccentricity divided by effective radius clearance) obtained based on gap measurements using an existing eddy current displacement meter, etc. When calculating the supporting force of the bearing using the capacity coefficient (= reciprocal of the Sommerfeld number) for the eccentricity, determine the circumferential position of the rotating shaft where the vertical supporting force is minimum and maximum,
Using both positions as a reference, calculate the vector sum of the supporting force at a point with the same phase with respect to the circumferential position of the rotating shaft at least twice with different phases, and calculate the gap measurement value so as to minimize the difference between these vector sums. Correct the calculation.
オ)作用
軸受が支持している力は通常、回転体の自重の該軸受に
対する分力、ミスアライメントによる力、不釣合いによ
って生じる回転に伴う変動力などの合成力である。本発
明の構成に従えば、該ベクトル和を求める過程で不釣合
い力が互いに打ち消しあい、その後にとる該ベクトル和
の差により回転体の自重の分力とミスアライメント力が
打ち消しあうことになるから、求めた結果は原理的に通
常は零とならなければならない、ところが、回転軸心の
計測値に偏差的誤差があるとそれを基にして算出した力
が実際の値と異なるものになっているから、該計測値の
生じている誤差の大きさに応じて打ち消し得ない力の成
分が残ってしまうことになる。従って、求める該結果が
零に近ずくように回転軸心の計測値に対し、例えば制御
分野で適性値を求める際に用いられる公知の「山登り法
」等により、適正な補正値を求めて取り除けば該誤差は
補正されてしまうことになる。E) Action The force supported by a bearing is usually a resultant force such as a component of the weight of the rotating body against the bearing, a force due to misalignment, and a fluctuating force due to rotation caused by unbalance. According to the configuration of the present invention, the unbalanced forces cancel each other in the process of calculating the vector sum, and the component force of the rotating body's own weight and the misalignment force cancel each other due to the difference in the vector sum obtained afterwards. In principle, the obtained result should normally be zero, but if there is a deviation error in the measured value of the rotation axis, the force calculated based on it will differ from the actual value. Therefore, a force component that cannot be canceled remains depending on the magnitude of the error in the measured value. Therefore, in order to bring the desired result close to zero, it is necessary to find and remove an appropriate correction value for the measured value of the axis of rotation, for example by using the well-known "hill climbing method" used to find appropriate values in the control field. If so, the error will be corrected.
力)実施例 第1図に示した1実施例に基づいて、詳細に説明する。Power) Example This will be explained in detail based on one embodiment shown in FIG.
ジャーナル軸1の軸心の位置を、例えば渦流式変位計な
どの検出器2−1.2−2によって、軸受中心を基準に
する形で測るように配置する。本実施例では、鉛直方向
に検出器2−1を水平方向に検出器2−2を配置した。A detector 2-1.2-2, such as an eddy current displacement meter, is arranged so that the position of the axis of the journal shaft 1 is measured with the center of the bearing as a reference. In this embodiment, the detector 2-1 is arranged vertically and the detector 2-2 is arranged horizontally.
(検出器の配置方向を変えても、又検出器の数を増やし
ても本発明の基本思想から逸脱はしない、)該検出器の
出力を測定器3によって、ジャーナル軸lの軸心の位置
の値に変換した後、演算器4において、以下に記す演算
をして、誤差補正を行った修正軸心位置を監視器5へ入
力して状態監視をするものである。(Even if the arrangement direction of the detectors is changed or the number of detectors is increased, the basic idea of the present invention will not be deviated from.) After converting into the value, the calculation unit 4 performs the calculations described below, and the corrected axis center position with error correction is input to the monitor 5 to monitor the state.
本発明に従って演算する際に対象となるジャーナル軸l
の軸心に加わる力が、該軸心の各位置状態において、ど
のような状態にあるかを第2図に示す。ジャーナル軸l
の軸心に加わる力の鉛直方向成分が最小な状態のジャー
ナル軸lの軸心の軸受内における位置をCJtl、同様
に最大の状態におけるそれをCJDとする。一般に、油
膜軸受の支持力は、軸受の中心C3に対するジャーナル
軸1の偏心率によって決る。偏心率と油膜軸受の支持力
の間の関係については、公知のシンマーフェルト数の逆
数である容量係数から与えられる。偏心率は、軸受中心
C6とジャーナル軸lの細心間の距離と軸受により決る
有効半径クリアランスとの比で与え得る。該距離は測定
器3の出力として、有効半径クリアランスは軸受の固有
値として得ることができる。したがって、回転軸の回転
に同期した例えば64パルス/回転のサンプリング・パ
ルスにより測定器3の出力としての該距離を求めれば、
偏心率を仲介として回転軸の1回転内における軸受支持
力の変化が解るから、その鉛直分力が最小・最大になる
軸心の位置の形で、容易にCJU+ CJlは割出せる
ことになる。また、CJul CJD基準で所定の位
相を有する位置における軸受支持力′も割出せることに
なる。(なお、各側出しが可能なようなサンプリング・
パルス数/回転にしておくことは当然であるから、本明
細書ではこれ以上詳しくは触れない事にする。) CJ
a、 CJIを基準にしてジャーナル軸lの回転位相
差α1.α2の状態のときのジャーナル軸1の軸心の位
置を、第2図に示すように各々CJII CJl、G
Jit CJ4とする。Journal axis l that is the target when calculating according to the present invention
FIG. 2 shows the state of the force applied to the axial center at each position of the axial center. journal shaft l
Let CJtl be the position in the bearing of the axis of journal shaft l when the vertical component of the force applied to the axis is minimum, and CJD be the position when it is maximum. Generally, the supporting force of an oil film bearing is determined by the eccentricity of the journal shaft 1 with respect to the center C3 of the bearing. The relationship between the eccentricity and the supporting force of an oil film bearing is given by the capacity coefficient, which is the reciprocal of the known Simmerfeld number. The eccentricity can be given by the ratio of the distance between the bearing center C6 and the center of the journal axis l to the effective radial clearance determined by the bearing. The distance can be obtained as the output of the measuring device 3, and the effective radial clearance can be obtained as the characteristic value of the bearing. Therefore, if the distance is determined as the output of the measuring device 3 using sampling pulses of, for example, 64 pulses/rotation synchronized with the rotation of the rotating shaft, then
Since the change in bearing support force within one revolution of the rotating shaft is known through eccentricity, CJU+CJl can be easily determined in terms of the position of the shaft center where the vertical component force is minimum and maximum. Furthermore, the bearing support force' at a position having a predetermined phase based on the CJul CJD standard can also be determined. (Please note that sampling is possible so that each side can be
Since it is natural to set the number of pulses/rotation as the number of pulses, it will not be discussed in further detail in this specification. ) C.J.
a, rotational phase difference α1 of journal axis l with reference to CJI; The positions of the axis of the journal shaft 1 in the state of α2 are CJII, CJl and G, respectively, as shown in Fig. 2.
Let's call it Jit CJ4.
又、ジャーナル軸lを支えている図には示していない軸
受が支持している各々の力について、回転体の自重の分
力をf Ilnミスアライメント力をf□不釣合いによ
る力をf8とすると、CJll CJllCJ3+
cJ4といった特徴点で、軸心に働く力は次に示すよう
な関係にある。Also, regarding the respective forces supported by the bearings (not shown in the figure) that support the journal shaft l, let the component force of the rotating body's own weight be f Iln The misalignment force be f □ The force due to unbalance be f8 , CJll CJllCJ3+
At a feature point such as cJ4, the forces acting on the axis have the following relationship.
但、ジャーナル軸lの回転角は第2図中のX輪を基準に
してY軸の方向に採っている。However, the rotation angle of the journal shaft l is taken in the direction of the Y axis with the X wheel in FIG. 2 as a reference.
CJl:水平方向の力Flx=
f、・CO3(β)+f、・COS <π/2−α、)
CJI;垂直方向の力Fly=
f、+f、・5IN(β)+f1SIN(π/2−α1
)C1寞;水平方向の力F2x=
f 1I−CO5(β) + f 、−CO3(π/2
+αχ)CJz:垂直方向の力F2y=
f 、+ f 、−5IN(β)+fII−SIN(π
/2+α2)C1,;水平方向の力F3x=
f 、−CO5(β)+fu−CO5(−π/2−α1
)C1,;垂直方向の力F3y=
f、十f、・5IN(β) + f tI−SIN(−
π/2−α1)CJa;水平方向の力F4x=
fl・CO5(β)+f1・CO5(−π/2+α:)
CJa;垂直方向の力F4y=
f、+f、・5IN(β)十 f @・5IN(−π/
2+ α2)CJI+ C10の所の力のベクトル和を
採り、その結果をF13とすると、
F 13 x=2 ・f、−COS(β)IF 13
y=2・ (f 、+ f ll−5IN(β))同様
に、C,!、C,4の所の力のベクトル和を採り、その
結果をF24とすると、
F24 x=2 ・ f、−CO3(β)。CJl: Horizontal force Flx = f, ・CO3(β)+f, ・COS <π/2−α,)
CJI; Vertical force Fly = f, +f, ・5IN(β)+f1SIN(π/2-α1
)C1寞;Horizontal force F2x= f 1I−CO5(β) + f ,−CO3(π/2
+αχ)CJz: Vertical force F2y=f, +f, -5IN(β)+fII-SIN(π
/2+α2)C1,;Horizontal force F3x=f, -CO5(β)+fu-CO5(-π/2-α1
) C1,; Vertical force F3y = f, 10f, ・5IN(β) + f tI-SIN(-
π/2-α1) CJa; Horizontal force F4x= fl・CO5(β)+f1・CO5(-π/2+α:)
CJa; Vertical force F4y = f, +f, ・5IN (β) 10 f @・5IN (−π/
2+ α2) CJI+ Take the vector sum of the forces at C10 and let the result be F13, then F 13 x=2 ・f, -COS(β)IF 13
y=2・(f,+f ll-5IN(β)) Similarly, C,! , C, 4, and let the result be F24, then F24 x=2 · f, -CO3(β).
F24y=2 ・ (f、十f、・5IN(β))但、
F13.F24の添字x、yは、x、 y方向成分であ
ることを示す。F24y=2 ・(f, 10f, ・5IN(β)) However,
F13. The subscripts x and y of F24 indicate components in the x and y directions.
従って、これら2つの力のベクトル和F13゜F24の
差FSをとると、明かにFSは零にならなければならな
い。一方、油膜軸受の支持力は、軸受の中心C1に対す
るジャーナル軸lの偏心率によって決るから該偏心率に
誤差があると、偏心率から求めた上記FSは零にならな
くなる。Therefore, if we take the difference FS between the vector sum F13°F24 of these two forces, it is clear that FS must become zero. On the other hand, the supporting force of the oil film bearing is determined by the eccentricity of the journal shaft l with respect to the center C1 of the bearing, so if there is an error in the eccentricity, the FS determined from the eccentricity will not be zero.
模擬的に測定器3の出力に誤差を生ぜしめて、前記FS
の値を見てみたのが第3図である。By simulating an error in the output of the measuring device 3, the FS
Figure 3 shows the values of .
第3図は、ジャーナル径200(閣〕、軸受長150〔
閣〕、f 、−5000X 9.8 (N )、f@−
,02[kgl x、5[M] x (2π60)”[
N]、有効半径クリアランス、15[謹1.粘性係敞、
07[N −s /M町、7ミスアライメントの力と方
向が200x9.8(N]、30 (deg) 、αt
=45 (deg) 、αg=45 (deg)の条件
において、第4図に示すような大きさr、、・方向T、
の誤差を与えて、ジャーナル軸lの軸心に対して軸受の
中心C,があたかもCI’にあるかのごとき状態を呈す
る時のFSの値を図示したものである。第3図から解る
ように、公知の「山登り法」などによって、FSが最小
となるようにしてゆけば通常内在している大きさrlj
+方向γ、の誤差を求めることができる。従って、測定
器3の出力から、水平方向でr、・cosTa。Figure 3 shows a journal with a diameter of 200 mm and a bearing length of 150 mm.
], f, -5000X 9.8 (N), f@-
,02[kgl x, 5[M] x (2π60)”[
N], effective radius clearance, 15 [1. viscous modulus,
07[N-s/M town, 7 Misalignment force and direction are 200x9.8 (N], 30 (deg), αt
=45 (deg) and αg=45 (deg), the size r, ・direction T, as shown in FIG.
This figure shows the value of FS when the center C of the bearing appears to be at CI' with respect to the axis of the journal shaft l, with an error of . As can be seen from Figure 3, if FS is minimized using the well-known "hill climbing method", the inherent size rlj
The error in the + direction γ can be found. Therefore, from the output of the measuring device 3, r, ·cosTa in the horizontal direction.
鉛直方向でr、・SIN r *の補正量を差し引けば
よい。It is sufficient to subtract the correction amount of r, ·SIN r * in the vertical direction.
キ)発明の効果
本発明に従えば、僅かな演算機能の追加により回転機械
の監視システムの信頼性を向上できて回転機械の状態の
長期トレンドの把握が可能になって、状態基準保全のた
めの監視・診断が行えるようになる。G) Effects of the Invention According to the present invention, the reliability of the monitoring system for rotating machinery can be improved by adding a small amount of calculation function, and long-term trends in the condition of rotating machinery can be grasped, making it possible to maintain condition-based maintenance. monitoring and diagnosis.
第1図は本発明の実施例の機器構成に関する図、第2図
は本発明のアルゴリズムを説明するための参考図、第3
図は計測器の出力誤差の大きさによる影響の例示図、第
4図は誤差の説明図である。
:ジャーナル軸、
2−1゜
2−2:検出器、
:測定器、
:演算器、
:監視器。
出
願
人
新
日本製鐵株式会社FIG. 1 is a diagram related to the equipment configuration of an embodiment of the present invention, FIG. 2 is a reference diagram for explaining the algorithm of the present invention, and FIG.
The figure is an example of the influence of the magnitude of the output error of the measuring instrument, and FIG. 4 is an explanatory diagram of the error. : Journal axis, 2-1゜2-2: Detector, : Measuring device, : Computing unit, : Monitoring device. Applicant Nippon Steel Corporation
Claims (1)
率に対する容量係数を用いて当該軸受の支持力を算出す
るに際して、鉛直方向の支持力が最小・最大になる回転
軸の周方向位置を割出し、両該位置を基準にして回転軸
の周方向位置に対する位相が同じ所の支持力のベクトル
和を位相を変えて少なくとも2回求め、それらベクトル
和間の差を最小にするように間隙測定値を演算補正する
ことを特徴とする回転機械状態監視システムの確度向上
方法。1. When calculating the supporting force of the bearing from the eccentricity of the rotating shaft with respect to the bearing center using the capacity coefficient for that eccentricity, find the circumferential position of the rotating shaft where the vertical supporting force is minimum and maximum. Index, calculate the vector sum of the supporting force at the same position with respect to the circumferential position of the rotating shaft at least twice with the same phase relative to the circumferential position of the rotating shaft based on both positions, and set the gap so as to minimize the difference between these vector sums. A method for improving the accuracy of a rotating machine condition monitoring system, characterized by calculating and correcting measured values.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2392290A JPH03229125A (en) | 1990-02-02 | 1990-02-02 | Accuracy improving method for rotary machine state monitoring system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2392290A JPH03229125A (en) | 1990-02-02 | 1990-02-02 | Accuracy improving method for rotary machine state monitoring system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH03229125A true JPH03229125A (en) | 1991-10-11 |
| JPH0571898B2 JPH0571898B2 (en) | 1993-10-08 |
Family
ID=12124009
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2392290A Granted JPH03229125A (en) | 1990-02-02 | 1990-02-02 | Accuracy improving method for rotary machine state monitoring system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH03229125A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2685114A1 (en) | 2012-07-12 | 2014-01-15 | Institut Für Verbundwerkstoffe GmbH | Sliding bearings and method for determining the wear of a sliding bearing |
| JP2014059234A (en) * | 2012-09-18 | 2014-04-03 | Mitsubishi Heavy Ind Ltd | Shaft bearing monitoring system, rotation machine and shaft bearing monitoring method |
-
1990
- 1990-02-02 JP JP2392290A patent/JPH03229125A/en active Granted
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2685114A1 (en) | 2012-07-12 | 2014-01-15 | Institut Für Verbundwerkstoffe GmbH | Sliding bearings and method for determining the wear of a sliding bearing |
| DE102012106295A1 (en) | 2012-07-12 | 2014-10-30 | Institut Für Verbundwerkstoffe Gmbh | Slide bearing and method for determining the wear of a plain bearing |
| JP2014059234A (en) * | 2012-09-18 | 2014-04-03 | Mitsubishi Heavy Ind Ltd | Shaft bearing monitoring system, rotation machine and shaft bearing monitoring method |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0571898B2 (en) | 1993-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Glienicke | Paper 13: Experimental investigation of the stiffness and damping coefficients of turbine bearings and their application to instability prediction | |
| US4494400A (en) | Wheel balancer two plane calibration apparatus and method | |
| US4464935A (en) | Shaft vibration evaluation | |
| US4644270A (en) | Apparatus for monitoring housed turbine blading to obtain blading-to-housing distance | |
| US4033042A (en) | Shaft alignment apparatus and method | |
| US20180327111A1 (en) | Propeller health monitoring | |
| CN106153256A (en) | A kind of magnetic suspension rotor high accuracy spot dynamic balance method | |
| US20090293613A1 (en) | System and Method for Active Detection of Asymmetry In Rotating Structures | |
| Brancati et al. | Non-linear stability analysis of a rigid rotor on tilting pad journal bearings | |
| JPH03229125A (en) | Accuracy improving method for rotary machine state monitoring system | |
| Manshin et al. | About the dynamic error of strain gauge torque measuring devices | |
| JP3379170B2 (en) | Non-repetitive rotation accuracy measuring device | |
| Wegener et al. | Traceability of torque transducers under rotating and dynamic operating conditions | |
| WO2015170296A1 (en) | Method and auxiliary apparatus for balancing a rotor of a gas turbine | |
| Taylor et al. | Measured performance of a highly preloaded three-lobe journal bearing–-part II: dynamic characteristics | |
| Foiles et al. | Balancing with phase only (single-plane and multiplane) | |
| Shawki et al. | Journal bearing performance under sinusoidally alternating and fluctuating loads | |
| Ozoegwu et al. | zPure analytical approach to rotational balancing | |
| Liska et al. | Measurement and evaluation of shaft torsional vibrations using shaft instantaneous angular velocity | |
| JPH02162241A (en) | How to measure bearing friction coefficient | |
| Middleton et al. | Paper 7: An Investigation into the Performance of Dynamically Loaded Journal Bearings: Experiment | |
| Ahn et al. | Optimal multi-segment cylindrical capacitive sensor | |
| JPS6057529B2 (en) | Unbalanced position evaluation system for rotating machines | |
| Titus et al. | Establishment and development of torque metrology in CSIR-NPL for providing the traceability in torque measurements to user industries | |
| JPH11142275A (en) | Balance testing machine |