JP4344601B2 - Magnetic bearing low disturbance control device - Google Patents

Magnetic bearing low disturbance control device Download PDF

Info

Publication number
JP4344601B2
JP4344601B2 JP2003427227A JP2003427227A JP4344601B2 JP 4344601 B2 JP4344601 B2 JP 4344601B2 JP 2003427227 A JP2003427227 A JP 2003427227A JP 2003427227 A JP2003427227 A JP 2003427227A JP 4344601 B2 JP4344601 B2 JP 4344601B2
Authority
JP
Japan
Prior art keywords
rotor
bearing
displacement
radial
command
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP2003427227A
Other languages
Japanese (ja)
Other versions
JP2005188545A (en
Inventor
光伯 齊藤
正夫 井上
一彦 福島
昭彦 今城
直樹 今村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2003427227A priority Critical patent/JP4344601B2/en
Publication of JP2005188545A publication Critical patent/JP2005188545A/en
Application granted granted Critical
Publication of JP4344601B2 publication Critical patent/JP4344601B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Description

この発明は、回転体(ロータ)を非接触で支持する磁気軸受が発生する擾乱を、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減する低擾乱化制御装置に関するものである。   The present invention relates to a low disturbance control device that reduces disturbance generated by a magnetic bearing that supports a rotating body (rotor) in a non-contact manner in all frequency ranges including a servo resonance frequency range of the magnetic bearing.

従来、磁気軸受の擾乱抑制装置として、同期妨害補償装置がある。この同期妨害補償装置は、ロータ半径方向位置x,yの検出手段と、二入力加算器から出力されるロータ半径方向位置x,yに基づいて、予め設定した半径方向位置にロータを保持するように軸受電磁石コイルに対して電流を供給する制御回路とを備えた通常の磁気軸受装置に対して、ロータ半径方向位置検出手段と制御回路との間に処理回路を配置した構成となっている。この処理回路は、アキシャル方向ロータ回転速度の周波数に中心を置く帯域消去フィルタとして作用し、ロータ半径方向位置x,yから回転同期成分を除去するため、回転同期周波数域における軸受剛性を著しく減少させる。これにより、ロータは慣性主軸周りに回転するため、ロータ回転時に磁気軸受装置が発生する力の振動、つまり擾乱を抑制することが可能となる(例えば、特許文献1参照)。   Conventionally, as a disturbance suppression device for a magnetic bearing, there is a synchronous interference compensation device. The synchronous disturbance compensating apparatus holds the rotor at a preset radial position based on the rotor radial position x, y detecting means and the rotor radial position x, y output from the two-input adder. In contrast to a normal magnetic bearing device having a control circuit for supplying current to the bearing electromagnet coil, a processing circuit is arranged between the rotor radial position detecting means and the control circuit. This processing circuit acts as a band elimination filter centered on the frequency of the axial direction rotor rotational speed, and removes the rotational synchronization component from the rotor radial position x, y, thereby significantly reducing the bearing rigidity in the rotational synchronization frequency range. . Thereby, since the rotor rotates around the inertial main shaft, vibration of the force generated by the magnetic bearing device when the rotor rotates, that is, disturbance can be suppressed (for example, refer to Patent Document 1).

また、従来、磁気軸受の共振周波数域通過を容易にする装置として、臨界周波数減衰装置がある。この臨界周波数減衰装置は、通常の磁気軸受装置に対して、進相網と並列にディジタル回路を配置した構成となっている。このディジタル回路は、アキシャル方向ロータ回転速度と等しい周波数に中心を置く周波数帯域において、ロータ半径方向位置x,yの位相または進相利得の選択的増幅を行う手段として作用し、回転同期周波数域における磁気軸受制御系の位相、即ち進相利得の増大を可能にする。これにより、回転同期周波数域において軸受剛性が増加し、特に、磁気軸受の共振周波数域においてはロータ振れ回り量を抑制して共振周波数域を容易に超えることが可能となる(例えば、特許文献2参照)。   Conventionally, there is a critical frequency attenuating device as a device for facilitating the passage of the magnetic bearing in the resonance frequency range. This critical frequency attenuating device has a configuration in which a digital circuit is arranged in parallel with a phase advance network with respect to a normal magnetic bearing device. This digital circuit acts as a means for selectively amplifying the phase or phase gain of the rotor radial position x, y in a frequency band centered at a frequency equal to the axial rotor rotational speed, It is possible to increase the phase of the magnetic bearing control system, that is, the phase advance gain. As a result, the bearing rigidity is increased in the rotation synchronization frequency range, and in particular, in the resonance frequency range of the magnetic bearing, it is possible to suppress the rotor swing amount and easily exceed the resonance frequency range (for example, Patent Document 2). reference).

また、従来、ロータの擾乱抑制装置として、磁気軸受制御装置がある。この磁気軸受制御装置は、同期妨害補償装置と等価なABS制御、及び臨界周波数減衰装置と等価なNクロス制御、FF制御の適用度を調整するゲインを、アキシャル方向ロータ回転速度に応じてCPUが設定する構成となっている。従って、アキシャル方向ロータ回転速度が磁気軸受の共振周波数域以外にある場合はABS制御のみを適用して、回転同期周波数域における軸受剛性を減少させて慣性主軸回りのロータ回転による低擾乱化を実現し、ロータ回転速度が磁気軸受の共振周波数域にある場合はNクロス制御、またはFF制御を適用することにより、回転同期周波数域における軸受剛性を増加させてロータ振れ回り量の抑制による共振周波数域の容易な通過を可能とすることができる(例えば、特許文献3参照)。   Conventionally, there is a magnetic bearing control device as a rotor disturbance suppressing device. In this magnetic bearing control device, the CPU adjusts the gain for adjusting the applicability of the ABS control equivalent to the synchronous disturbance compensation device, the N-cross control equivalent to the critical frequency attenuation device, and the FF control according to the axial direction rotor rotational speed. It is configured to set. Therefore, when the rotational speed of the rotor in the axial direction is outside the resonance frequency range of the magnetic bearing, only the ABS control is applied to reduce the bearing rigidity in the rotational synchronization frequency range, thereby realizing low disturbance due to the rotor rotation around the inertia spindle. When the rotor rotation speed is in the resonance frequency range of the magnetic bearing, the N cross control or FF control is applied to increase the bearing rigidity in the rotation synchronization frequency range, thereby reducing the rotor swing amount. Can be easily passed (see, for example, Patent Document 3).

特公昭60−14929号公報(第3頁−第5頁、第4図及び8図)Japanese Examined Patent Publication No. 60-14929 (pages 3 to 5, FIGS. 4 and 8) 特開昭52−93853号公報(第3頁−第4頁、第4図及び第6図)Japanese Patent Laid-Open No. 52-93853 (pages 3 to 4, FIGS. 4 and 6) 特開平06−193633号公報(第10頁−14頁、図1)JP-A-06-193633 (pages 10-14, FIG. 1)

上記特許文献1の制御方法では、ロータ回転同期周波数域のみ軸受剛性を極度に下げることで、ロータを慣性主軸回りに回転させて低擾乱化を図っている。この従来の制御方法で低擾乱化を実現するためには、回転同期周波数域における軸受剛性の低下量を所定量以上に設定する必要があるが、磁気軸受のサーボ共振周波数域において軸受剛性の低下量を大きくするとロータの安定浮上が損なわれるために、サーボ共振周波数域では十分な擾乱抑制が実現できないという問題があった。   In the control method of Patent Document 1, the bearing rigidity is extremely lowered only in the rotor rotation synchronization frequency region, so that the rotor is rotated around the inertia main shaft to reduce disturbance. In order to achieve low disturbance with this conventional control method, it is necessary to set the amount of decrease in bearing rigidity in the rotation synchronization frequency range to a predetermined value or more. However, the decrease in bearing rigidity in the servo resonance frequency range of magnetic bearings. When the amount is increased, the stable flying of the rotor is impaired, so that there is a problem that sufficient disturbance suppression cannot be realized in the servo resonance frequency range.

また、上記特許文献2及び3の制御方法では、回転同期周波数域における軸受剛性を増加させてロータ振れ回り量の抑制による共振周波数域の容易な通過を可能とすることができ、特許文献3では、さらに、ロータ回転同期周波数域のみ軸受剛性を極度に下げることで、ロータを慣性主軸回りに回転させて低擾乱化を図ることができるが、サーボ共振周波数域では軸受剛性を増加させることによって、発生擾乱も増加してしまうという問題がある。   Further, in the control methods of Patent Documents 2 and 3, the bearing rigidity in the rotation synchronous frequency region can be increased to allow easy passage of the resonance frequency region by suppressing the rotor swing amount. Furthermore, by significantly lowering the bearing rigidity only in the rotor rotation synchronization frequency range, the rotor can be rotated around the inertia main shaft to reduce disturbance, but in the servo resonance frequency range, by increasing the bearing rigidity, There is a problem that the occurrence disturbance increases.

この発明は、上記のような問題点を解決するためになされたものであり、磁気軸受装置が発生する擾乱を、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減する低擾乱化制御装置を提供することを目的としている。   The present invention has been made to solve the above-described problems, and is a low disturbance control that reduces the disturbance generated by the magnetic bearing device in the entire frequency range including the servo resonance frequency range of the magnetic bearing. The object is to provide a device.

この発明に係る第1の磁気軸受の低擾乱化制御装置は、ロータを電磁力によって浮上させる複数の電磁石からなる磁気軸受の、ラジアル方向をX軸及びY軸、中心軸をZ軸とする固定座標系において、上記ロータの幾何中心軸が上記ラジアル方向に並進運動するラジアル方向並進変位及び上記ロータの幾何中心軸が上記固定座標系のアキシャル方向(Z軸方向)に並進運動するアキシャル方向並進変位の3自由度ロータ変位と、上記ロータの幾何中心軸が上記固定座標系のZ軸を含む平面内で回動するラジアル方向回転変位の2自由度ロータ変位からなる5自由度ロータ変位及びアキシャル方向ロータ回転角度に関係する情報(ロータ変位情報)を検出する変位検出手段と、
上記変位検出手段の出力を、上記5自由度ロータ変位に変換するロータ変位演算器と、 上記ロータ変位演算器の出力に基づいて、上記ロータを磁気軸受の所定の位置に制御するための軸受力となるラジアル方向及びアキシャル方向軸受力指令、上記ロータのラジアル方向回転変位を制御する軸受力となるラジアル方向軸受モーメント指令を出力する制御指令出力手段と、
上記ロータ変位演算器の出力に基づいて、上記電磁石各々と上記ロータとのギャップ長を出力するギャップ長演算器と
上記ギャップ長演算器の出力及び上記制御指令出力手段の出力に基づいて、上記制御指令出力手段から出力されるラジアル方向及びアキシャル方向軸受力指令、及びラジアル方向軸受モーメント指令を実現するための上記電磁石各々に対する供給電流を分配して供給する電力分配器とを備えた軸受制御系、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記軸受制御系における制御指令出力手段が出力するラジアル方向軸受力指令を、上記ロータの幾何中心軸をZ軸、上記Z軸に垂直な方向をX軸及びY軸としてロータと共にZ軸まわりに回転する回転座標系表現に変換する回転座標変換器と、
上記回転座標変換器の出力のロータ回転同期成分を抽出するLPFと、
上記LPFの出力を積分して、ラジアル方向に並進させる回転座標系のラジアル方向並進変位を出力する積分器と、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記積分器が出力するラジアル方向並進変位を上記固定座標系表現に変換して上記ロータの重心位置が上記回転座標系原点に位置するように上記ロータの幾何中心軸をラジアル方向に並進させるラジアル方向並進変位指令を算出する固定座標変換器を備えた並進変位指令制御系、
を備え、
上記固定座標系に変換したラジアル方向並進変位指令により、上記軸受制御系が上記電磁石を制御して上記ロータの浮上位置を所定位置に制御するものである。
A first disturbance control device for a magnetic bearing according to a first aspect of the present invention is a fixed magnetic bearing comprising a plurality of electromagnets that levitate a rotor by electromagnetic force, with a radial direction as an X axis and a Y axis, and a central axis as a Z axis. In the coordinate system, a radial translational displacement in which the geometric center axis of the rotor translates in the radial direction, and an axial translational displacement in which the geometric center axis of the rotor translates in the axial direction (Z-axis direction) of the fixed coordinate system. 5-degree-of-freedom rotor displacement and axial direction consisting of a 3-degree-of-freedom rotor displacement and a 2-degree-of-freedom rotor displacement in which the geometric center axis of the rotor rotates in a plane including the Z-axis of the fixed coordinate system Displacement detecting means for detecting information related to the rotor rotation angle (rotor displacement information);
A rotor displacement calculator for converting the output of the displacement detector into the 5-degree-of-freedom rotor displacement; and a bearing force for controlling the rotor to a predetermined position of the magnetic bearing based on the output of the rotor displacement calculator. A control command output means for outputting a radial direction and axial direction bearing force command, and a radial direction bearing moment command serving as a bearing force for controlling the radial rotational displacement of the rotor;
Based on the output of the rotor displacement calculator, a gap length calculator that outputs the gap length between each of the electromagnets and the rotor ;
The electromagnet for realizing the radial direction and axial direction bearing force command and the radial direction bearing moment command output from the control command output means based on the output of the gap length calculator and the output of the control command output means. A bearing control system including a power distributor that distributes and supplies a supply current to each ;
Based on the rotor rotational angle in the axial direction of the rotor, the radial bearing force command output by the control command output means in the bearing control system is expressed as follows: the geometric center axis of the rotor is the Z axis, and the direction perpendicular to the Z axis is X A rotary coordinate converter for converting into a rotary coordinate system representation that rotates about the Z axis with the rotor as an axis and a Y axis;
LPF for extracting the rotor rotation synchronization component of the output of the rotary coordinate converter;
An integrator that integrates the output of the LPF and outputs a radial translational displacement of a rotating coordinate system that translates in the radial direction;
Based on the axial rotor rotation angle of the rotor, the radial translation displacement output from the integrator is converted into the fixed coordinate system representation so that the center of gravity of the rotor is located at the origin of the rotary coordinate system. translational displacement command control system that includes a fixed coordinate converter for calculating a radial translational displacement command for translating the geometric center axis of the radial direction,
With
The bearing control system controls the electromagnet to control the flying position of the rotor to a predetermined position by a radial translational displacement command converted to the fixed coordinate system.

この発明に係る第2の磁気軸受の低擾乱化制御装置は、ロータを電磁力によって浮上させる複数の電磁石からなる磁気軸受の、ラジアル方向をX軸及びY軸、中心軸をZ軸とする固定座標系において、上記ロータの幾何中心軸が上記ラジアル方向に並進運動するラジアル方向並進変位及び上記ロータの幾何中心軸が上記固定座標系のアキシャル方向(Z軸方向)に並進運動するアキシャル方向並進変位の3自由度ロータ変位と、上記ロータの幾何中心軸が上記固定座標系のZ軸を含む平面内で回動するラジアル方向回転変位の2自由度ロータ変位からなる5自由度ロータ変位及びアキシャル方向ロータ回転角度に関係する情報(ロータ変位情報)を検出する変位検出手段と、
上記変位検出手段の出力を、上記5自由度ロータ変位に変換するロータ変位演算器と、
上記ロータ変位演算器の出力に基づいて、上記ロータを磁気軸受の所定の位置に制御するための軸受力となるラジアル方向及びアキシャル方向軸受力指令、上記ロータのラジアル方向回転変位を制御する軸受力となるラジアル方向軸受モーメント指令を出力する制御指令出力手段と、
上記ロータ変位演算器の出力に基づいて、上記電磁石各々と上記ロータとのギャップ長を出力するギャップ長演算器と
上記ギャップ長演算器の出力及び上記制御指令出力手段の出力に基づいて、上記制御指令出力手段から出力されるラジアル方向及びアキシャル方向軸受力指令、及びラジアル方向軸受モーメント指令を実現するための上記電磁石各々に対する供給電流を分配して供給する電力分配器とを備えた軸受制御系、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記軸受制御系における制御指令出力手段が出力するラジアル方向軸受モーメント指令を、上記ロータの幾何中心軸をZ軸、上記Z軸に垂直な方向をX軸及びY軸としてロータと共にZ軸周りに回転する回転座標系表現に変換する回転座標変換器と、
上記回転座標変換器の出力のロータ回転同期成分を抽出するLPFと、
上記LPFの出力を積分してラジアル方向に回転させる回転座標系のラジアル方向回転変位を出力する積分器
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記積分器から出力されるラジアル方向回転変位を上記固定座標系表現に変換して上記ロータの慣性主軸が上記磁気軸受けの中心軸と平行になるように上記ロータの幾何中心軸を上記固定座標系のZ軸を含む平面内で回動させるラジアル方向回転変位指令を算出する固定座標変換器とを備えた回転変位指令制御系、
を備え、
上記固定座標系に変換したラジアル方向回転変位指令により、上記軸受制御系が上記電磁石を制御して上記ロータの浮上位置を所定位置に制御するものである。
A second disturbance control device for a magnetic bearing according to the present invention is a magnetic bearing composed of a plurality of electromagnets that levitates a rotor by electromagnetic force, with a radial direction as an X axis and a Y axis, and a central axis as a Z axis. In the coordinate system, a radial translational displacement in which the geometric center axis of the rotor translates in the radial direction, and an axial translational displacement in which the geometric center axis of the rotor translates in the axial direction (Z-axis direction) of the fixed coordinate system. 5-degree-of-freedom rotor displacement and axial direction consisting of a 3-degree-of-freedom rotor displacement and a 2-degree-of-freedom rotor displacement in which the geometric center axis of the rotor rotates in a plane including the Z-axis of the fixed coordinate system Displacement detecting means for detecting information related to the rotor rotation angle (rotor displacement information);
A rotor displacement calculator for converting the output of the displacement detection means into the 5-degree-of-freedom rotor displacement;
Based on the output of the rotor displacement calculator, a radial and axial direction bearing force command that becomes a bearing force for controlling the rotor to a predetermined position of the magnetic bearing, and a bearing force for controlling the radial rotational displacement of the rotor Control command output means for outputting a radial bearing moment command,
Based on the output of the rotor displacement calculator, a gap length calculator that outputs the gap length between each of the electromagnets and the rotor ;
The electromagnet for realizing the radial direction and axial direction bearing force command and the radial direction bearing moment command output from the control command output means based on the output of the gap length calculator and the output of the control command output means. A bearing control system including a power distributor that distributes and supplies a supply current to each ;
Based on the rotor rotational angle in the axial direction of the rotor, the radial bearing moment command output by the control command output means in the bearing control system is expressed as follows: the geometric center axis of the rotor is the Z axis, and the direction perpendicular to the Z axis is X A rotary coordinate converter that converts the rotary coordinate system to rotate around the Z axis together with the rotor as an axis and a Y axis;
LPF for extracting the rotor rotation synchronization component of the output of the rotary coordinate converter;
An integrator for outputting a radial rotational displacement of the rotating coordinate system that rotates in the radial direction by integrating the output of the LPF,
Based on the rotor rotational angle in the axial direction of the rotor, the radial rotational displacement output from the integrator is converted into the fixed coordinate system representation so that the inertial main shaft of the rotor is parallel to the central axis of the magnetic bearing. A rotational displacement command control system comprising: a fixed coordinate converter that calculates a radial rotational displacement command for rotating the geometric center axis of the rotor in a plane including the Z axis of the fixed coordinate system;
With
The bearing control system controls the electromagnet to control the flying position of the rotor to a predetermined position by a radial rotational displacement command converted into the fixed coordinate system.

この発明に係る第1の磁気軸受の低擾乱化制御装置によれば、ロータの重心位置を上記回転座標系Z軸上に位置するようにしているので、ロータ静不釣合いに起因する擾乱力を低減することができる。   According to the first disturbance control device for a magnetic bearing according to the present invention, since the position of the center of gravity of the rotor is positioned on the rotational coordinate system Z-axis, the disturbance force caused by the rotor unbalance is reduced. Can be reduced.

この発明に係る第2の磁気軸受の低擾乱化制御装置によれば、ロータ慣性主軸を磁気軸受中心軸と平行になるようにしているので、ロータ動不釣合いに起因する擾乱モーメントを低減することができる。   According to the second disturbance control device for a magnetic bearing according to the present invention, since the rotor inertia main shaft is made parallel to the magnetic bearing central axis, the disturbance moment caused by rotor dynamic imbalance is reduced. Can do.

実施の形態1.
図1は、磁気軸受装置の構成要素であるロータ、電磁石、及び変位検出手段の配置例を示す構成図であり、(a)は正面図、(b)は平面図である。
Embodiment 1 FIG.
FIG. 1 is a configuration diagram showing an arrangement example of a rotor, an electromagnet, and a displacement detection means that are components of a magnetic bearing device, where (a) is a front view and (b) is a plan view.

図1に示したように、この実施の形態1における磁気軸受装置は、ロータ1、ロータ1を電磁力によって所定位置に浮上させる電磁石2からなる磁気軸受、ロータ1の磁気軸受中心からの変位を検出する変位検出手段3を備えており、X軸、Y軸、及びZ軸は磁気軸受(図示していないステータに対して固定されている)に固定された固定座標系の座標軸で、Z軸は磁気軸受中心軸(アキシャル方向)、X軸、Y軸は磁気軸受中心軸に垂直で、互いに直交するラジアル軸であり、固定座標系原点はロータ1の幾何中心位置となっている。   As shown in FIG. 1, the magnetic bearing device according to the first embodiment includes a rotor 1, a magnetic bearing composed of an electromagnet 2 that floats the rotor 1 to a predetermined position by electromagnetic force, and displacement of the rotor 1 from the center of the magnetic bearing. Displacement detecting means 3 for detecting is provided, and the X axis, the Y axis, and the Z axis are coordinate axes of a fixed coordinate system fixed to a magnetic bearing (fixed to a stator not shown), and the Z axis Is a magnetic bearing central axis (axial direction), X axis and Y axis are radial axes perpendicular to and perpendicular to the magnetic bearing central axis, and the origin of the fixed coordinate system is the geometric center position of the rotor 1.

また、ロータ1の3自由度並進変位X,Y,Z及びラジアル方向2自由度回転変位θx、θyを制御する磁気軸受装置において、円柱状のロータ1の周囲に10個の電磁石2a〜2j及び10個の変位検出手段3a〜3jを配置している。   Further, in the magnetic bearing device for controlling the three-degree-of-freedom translational displacements X, Y, Z and the radial two-degree-of-freedom rotational displacements θx, θy of the rotor 1, ten electromagnets 2a-2j around the cylindrical rotor 1 and Ten displacement detection means 3a to 3j are arranged.

ロータ1には、ロータ1の幾何中心軸が固定座標系原点からラジアル方向に並進運動するラジアル方向並進変位X,Y、ロータ1の幾何中心軸が固定座標のアキシャル方向(Z軸方向)に並進運動するアキシャル方向並進変位Z、及びロータ1の幾何中心軸が固定座標系のZ軸を含む平面内で回動するラジアル方向回転変位θx、θyといったロータ変位が生じる。   In the rotor 1, a radial translational displacement X, Y in which the geometric center axis of the rotor 1 translates in the radial direction from the origin of the fixed coordinate system, and the geometric center axis of the rotor 1 translates in the axial direction (Z-axis direction) of the fixed coordinates. A rotor displacement such as an axial translational displacement Z that moves, and a radial rotational displacement θx and θy in which the geometric center axis of the rotor 1 rotates in a plane including the Z axis of the fixed coordinate system occurs.

この構成の磁気軸受装置では、変位検出手段3a〜3hから出力される、ロータ変位に関する情報(以下、ロータ変位情報という)によって、ロータ1のラジアル方向並進変位X,Y及びラジアル方向回転変位θx、θyを算出することができ、変位検出手段3i〜3jから出力されるロータ変位情報によってロータ1のアキシャル方向並進変位Zを算出することができる(すなわち、5自由度ロータ変位X,Y,Z、θx、θyを算出することができる)。また、電磁石2a〜2hが発生する電磁力によって、ロータ1に作用するラジアル方向軸受力及びラジアル方向軸受モーメントを制御することができ、電磁石2i〜2jが発生する電磁力によってロータ1に作用するアキシャル方向軸受力を制御することができる。   In the magnetic bearing device having this configuration, the radial translational displacements X and Y of the rotor 1 and the radial rotational displacement θx, based on the rotor displacement information (hereinafter referred to as rotor displacement information) output from the displacement detectors 3a to 3h. θy can be calculated, and the axial translational displacement Z of the rotor 1 can be calculated from the rotor displacement information output from the displacement detectors 3i to 3j (that is, the 5-degree-of-freedom rotor displacement X, Y, Z, θx and θy can be calculated). Moreover, the radial bearing force and the radial bearing moment acting on the rotor 1 can be controlled by the electromagnetic force generated by the electromagnets 2a to 2h, and the axial acting on the rotor 1 by the electromagnetic force generated by the electromagnets 2i to 2j. The directional bearing force can be controlled.

図2は、傾斜磁極リング状ロータ1の周囲に6個の電磁石2a〜2f及び6個の変位検出手段3a〜3fを配置した磁気軸受装置の正面図(a)及びA−A断面図(b)であり、変位検出手段3a〜3fから出力されるロータ変位情報によって5自由度ロータ変位X,Y,Z、θx、θyを算出することができ、また、電磁石2a〜2fが発生する電磁力によってロータに作用するラジアル方向軸受力、アキシャル方向軸受力及びラジアル方向軸受モーメントを制御することができる。   2A and 2B are a front view and a cross-sectional view taken along line AA of a magnetic bearing device in which six electromagnets 2a to 2f and six displacement detecting means 3a to 3f are arranged around the inclined magnetic pole ring-shaped rotor 1. FIG. 5 degrees of freedom rotor displacement X, Y, Z, θx, θy can be calculated from the rotor displacement information output from the displacement detection means 3a-3f, and the electromagnetic force generated by the electromagnets 2a-2f Thus, the radial bearing force, the axial bearing force and the radial bearing moment acting on the rotor can be controlled.

図3は、この発明に係る低擾乱化制御装置における実施の形態1を説明するためのブロック図、図4は、実施の形態1を説明するためのロータ半径方向断面図である。   FIG. 3 is a block diagram for explaining the first embodiment of the low disturbance control device according to the present invention, and FIG. 4 is a rotor radial cross-sectional view for explaining the first embodiment.

図3に示したように、低擾乱化制御装置は、ロータ1の幾何中心軸を磁気軸受の所定の位置に制御するための軸受力となる軸受力指令Fx,Fy,Fz及びロータ1のラジアル方向回転変位を制御する軸受力となるラジアル方向軸受モーメント指令Nx,Nyを出力する軸受制御系10及び軸受制御系10に対して軸受力指令Fx,Fyからロータ1のラジアル方向並進変位指令を算出して出力する並進変位指令制御系20を備える。   As shown in FIG. 3, the low disturbance control device includes bearing force commands Fx, Fy, Fz which are bearing forces for controlling the geometric central axis of the rotor 1 to a predetermined position of the magnetic bearing, and the radial of the rotor 1. A radial translational displacement command of the rotor 1 is calculated from the bearing force commands Fx and Fy to the bearing control system 10 and the bearing control system 10 that output radial bearing moment commands Nx and Ny that are bearing forces that control the rotational displacement of the direction. The translation displacement command control system 20 is provided.

軸受制御系10は、ロータ1、ロータ1を電磁力によって浮上させる電磁石2及びロータ1の変位を検出する変位検出手段3からなる構成要素11、n個の変位検出手段(変位センサ)3から出力されるロータ変位情報S1,S2・・・Snに基づいて、ロータ1の3自由度並進変位X,Y,Z及びラジアル方向2自由度回転変位θx、θyからなる5自由度ロータ変位を算出するロータ変位演算器12、5自由度ロータ変位X,Y,Z,θx,θyに基づいて、磁気軸受を構成するn個の電磁石2の各々におけるギャップ長G1,G2,・・・Gnを算出するギャップ長演算器13、ロータ変位の偏差(5自由度ロータ変位X,Y,Z,θx,θyとラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycとの偏差)に基づいて軸受力指令Fx,Fy,Fz及びラジアル方向軸受モーメント指令Nx,Nyを算出する、制御指令出力手段としてのPID制御器15、軸受力指令Fx,Fy,Fz、軸受モーメント指令Nx,Ny及び軸受電磁石ギャップ長G1,G2,・・・Gnに基づいて、電磁石2の各々に対して供給電流I1,I2,・・・,Inを分配して供給する電力分配器14を備える。図にはX軸方向並進運動に対するPID制御器15x及びY軸方向並進運動に対するPID制御器15yのみを図示している。 The bearing control system 10 outputs from the rotor 1, the electromagnet 2 that floats the rotor 1 by electromagnetic force, and the component 11 including the displacement detector 3 that detects the displacement of the rotor 1, and n displacement detectors (displacement sensors) 3. Based on the rotor displacement information S1, S2... Sn, the 5-degree-of-freedom rotor displacement comprising the 3-degree-of-freedom translational displacements X, Y, Z of the rotor 1 and the radial direction 2-degree-of-freedom rotational displacements θx, θy is calculated. Based on the rotor displacement calculator 12, 5-degree-of-freedom rotor displacement X, Y, Z, θx, θy, gap lengths G1, G2,... Gn in each of the n electromagnets 2 constituting the magnetic bearing are calculated. Gap length calculator 13, rotor displacement deviation (5-degree-of-freedom rotor displacement X, Y, Z, θx, θy and radial translational displacement commands Xc, Yc, axial translational displacement command Zc, radial rotation) Position command Shitaxc, bearing force command Fx on the basis of the deviation between the θyc), Fy, Fz and radial bearing moment command Nx, calculates the Ny, PID controller 15 as a control command output means, bearing force command Fx, Fy , Fz, bearing moment commands Nx, Ny and bearing electromagnet gap lengths G1, G2,... Gn, supply currents I1, I2,. A power distributor 14 is provided. Only the PID controller 15x for the X-axis direction translational movement and the PID controller 15y for the Y-axis direction translational movement are illustrated.

また、並進変位指令制御系20は、軸受制御系10において出力されたラジアル方向軸受力指令Fx,Fyを、アキシャル方向ロータ回転角度θzを用いてロータに固定された回転座標系表現Fxr,Fyrに変換する回転座標変換器21、ラジアル方向軸受力指令の回転座標系表現Fxr,Fyrに対して、通過低周波数帯域を大略1Hz程度の狭帯域とすることで、ラジアル方向軸受力指令のロータ回転同期成分のみを抽出するLPF(ローパス・フィルタ)22、LPF22の出力を積分し、ロータ1のラジアル方向並進変位指令の回転座標系表現Xcr,Ycrを出力する積分器23、積分器23の出力であるラジアル方向並進変位指令の回転座標系表現Xcr,Ycrを、アキシャル方向ロータ回転角度θzを用いて固定座標系表現Xc,Ycに変換する固定座標変換器24を備えている。   Further, the translational displacement command control system 20 converts the radial bearing force commands Fx and Fy output from the bearing control system 10 into rotational coordinate system expressions Fxr and Fyr fixed to the rotor using the axial direction rotor rotation angle θz. The rotation low-frequency band is set to a narrow band of about 1 Hz with respect to the rotary coordinate converter 21 to be converted and the rotational coordinate system representation Fxr, Fyr of the radial direction bearing force command, so that the rotor rotation synchronization of the radial direction bearing force command is performed. These are the outputs of an integrator 23 and an integrator 23 that integrate the outputs of LPF (low-pass filter) 22 and LPF 22 that extract only components, and output rotational coordinate system expressions Xcr and Ycr of the radial translational displacement command of the rotor 1. Rotational coordinate system expression Xcr, Ycr of radial direction translational displacement command is expressed in fixed coordinate system using axial direction rotor rotation angle θz c, and a fixed coordinate converter 24 which converts the Yc.

図4は、ロータ半径方向断面を示す断面図であり、図4(a)は並進変位指令制御系を適用した直後、図4(b)は並進変位指令制御系を適用後十分時間が経過した後の状態を示している。図4において、ロータ1の回転軸に垂直で互いに直交し、かつ、ロータに対して固定された回転座標系の座標軸をXr,Yrとし、この回転座標系の原点は上記固定座標系原点と一致している。また、ロータの重心位置をG、変位検出手段で検出したロータ変位情報により決定するロータの幾何中心位置をS、ロータの静不釣合いに起因する振れ回り力を30、回転座標変換器21から出力されるラジアル方向軸受力指令の回転座標系表現Fxr,Fyrを要素とするベクトルを31、ラジアル方向並進変位指令の回転座標系表現Xcr,Ycrを要素とするベクトルを32としている。   FIG. 4 is a cross-sectional view showing a cross section in the radial direction of the rotor. FIG. 4 (a) is just after applying the translational displacement command control system, and FIG. 4 (b) is enough time after applying the translational displacement command control system. The later state is shown. In FIG. 4, the coordinate axes of the rotational coordinate system perpendicular to the rotational axis of the rotor 1 and orthogonal to each other and fixed to the rotor are Xr and Yr. The origin of the rotational coordinate system is the same as the fixed coordinate system origin. I'm doing it. Further, the center of gravity position of the rotor is G, the geometrical center position of the rotor determined by the rotor displacement information detected by the displacement detecting means is S, the swinging force caused by the static imbalance of the rotor is 30, and the rotational coordinate converter 21 outputs 31 is a vector whose elements are the rotational coordinate system representations Fxr and Fyr of the radial bearing force command, and 32 is a vector whose elements are the rotational coordinate system representations Xcr and Ycr of the radial translational displacement command.

次に、動作について説明する。軸受制御系10において、ロータ変位を検出する電磁石/ロータ/変位検出手段からなる構成要素11から出力されるロータ変位情報Si(i=1、2、…、n)は、ロータ変位演算器12によってロータの3自由度並進変位X、Y、Z及びラジアル方向2自由度回転変位θx、θyからなる5自由度のロータ変位に変換される。この5自由度ロータ変位はフィードバックされ、5自由度ロータ変位とラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycとの偏差信号をPID制御器15に入力する。PID制御器15では入力された偏差信号に基づいて、磁気軸受に対する軸受力指令Fx、Fy、Fz及びラジアル方向軸受モーメント指令Nx、Nyを算出して出力する。ギャップ長演算器13では、ロータ変位演算器12から出力される5自由度ロータ変位に基づいて、磁気軸受を構成するn個の電磁石の各々におけるギャップ長G1、G2、…、Gnを算出して出力する。電力分配器14では、ギャップ長演算器13から出力される軸受電磁石ギャップ長Gi(i=1、2、…、n)において、PID制御器15から出力される軸受力指令Fx、Fy、Fz及びラジアル方向軸受モーメント指令Nx、Nyを実現するために、n個の電磁石の各々に対する供給電流I1、I2、…、Inを分配して供給する。
ただし図3では、簡単のためにアキシャル方向並進運動及びラジアル方向回転運動のフィードバック・ループは図示していない。
Next, the operation will be described. In the bearing control system 10, rotor displacement information Si (i = 1, 2,..., N) output from a component 11 comprising an electromagnet / rotor / displacement detecting means for detecting rotor displacement is transmitted by a rotor displacement calculator 12. The rotor is converted into a 5-degree-of-freedom rotor displacement consisting of a 3-degree-of-freedom translational displacement X, Y, Z of the rotor and a radial direction 2-degree-of-freedom rotational displacement θx, θy. The 5-degree-of-freedom rotor displacement is fed back, and a deviation signal between the 5-degree-of-freedom rotor displacement and the radial translational displacement commands Xc and Yc, the axial translational displacement command Zc, and the radial rotational displacement commands θxc and θyc is sent to the PID controller 15. input. The PID controller 15 calculates and outputs bearing force commands Fx, Fy, Fz and radial bearing moment commands Nx, Ny for the magnetic bearing based on the input deviation signal. The gap length calculator 13 calculates gap lengths G1, G2,..., Gn in each of the n electromagnets constituting the magnetic bearing based on the 5-degree-of-freedom rotor displacement output from the rotor displacement calculator 12. Output. In the power distributor 14, in the bearing electromagnet gap length Gi (i = 1, 2,..., N) output from the gap length calculator 13, the bearing force commands Fx, Fy, Fz output from the PID controller 15 and In order to realize the radial bearing moment commands Nx, Ny, supply currents I1, I2,..., In are distributed and supplied to each of the n electromagnets.
However, in FIG. 3, for the sake of simplicity, the feedback loop for the axial translational motion and the radial rotational motion is not shown.

このような軸受制御系10の動作により、ロータ浮上位置がラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θyc(5自由度のロータ変位指令)と一致するように制御することが可能となる。   By such operation of the bearing control system 10, the rotor floating position coincides with the radial translational displacement commands Xc and Yc, the axial translational displacement command Zc, and the radial rotational displacement commands θxc and θyc (5-degree-of-freedom rotor displacement command). It becomes possible to control to do.

しかしながら、従来の磁気軸受制御装置では、5自由度のロータ変位指令を常に零値としており、従ってアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域よりも低い場合のロータ半径方向断面図は、理想的には図4(a)に示す通りとなる。すなわち、回転座標系Xr−Yrから見た場合、ロータ幾何中心位置Sは回転座標系の原点に位置し、ロータ重心位置Gはロータ幾何中心位置Sからロータ静不釣合いに対応する距離を隔てた点に位置する。その結果、ロータ重心位置Gには静不釣合いに起因する振れ回り力(遠心力)30が作用し、振れ回り力30の作用下でロータ浮上位置、すなわちロータ幾何中心位置Sを回転座標系原点に保持するために、ラジアル方向軸受力指令の回転座標系表現Fxr、Fyrを要素とするベクトル31は振れ回り力30と逆極性ベクトルとなる。このとき、ロータに作用するラジアル方向軸受力の反作用がステータに作用し、これが磁気軸受装置のハウジングに伝達する力の振動、すなわち、擾乱の主要因となる。   However, in the conventional magnetic bearing control device, the rotor displacement command of 5 degrees of freedom is always set to zero, and therefore the rotor radial cross-sectional view when the axial direction rotor rotational speed is lower than the position control band of the bearing control system 10 is Ideally, this is as shown in FIG. That is, when viewed from the rotational coordinate system Xr-Yr, the rotor geometric center position S is located at the origin of the rotational coordinate system, and the rotor center-of-gravity position G is separated from the rotor geometric center position S by a distance corresponding to the rotor unbalance. Located at a point. As a result, a swinging force (centrifugal force) 30 caused by static imbalance acts on the rotor center of gravity position G, and the rotor floating position, that is, the rotor geometric center position S is determined as the origin of the rotational coordinate system under the action of the swinging force 30. Therefore, the vector 31 having the rotational coordinate system representations Fxr and Fyr of the radial bearing force command as the elements becomes the swinging force 30 and the reverse polarity vector. At this time, the reaction of the radial bearing force acting on the rotor acts on the stator, which becomes the main factor of vibration of the force transmitted to the housing of the magnetic bearing device, that is, disturbance.

そこで、この実施の形態1の低擾乱化制御装置では、並進変位指令制御系20を設け、並進変位指令制御系20において、ラジアル方向軸受力指令Fx、Fy及びアキシャル方向ロータ回転角度θzに基づいて、ロータ重心位置Gが回転座標系原点に位置するようにラジアル方向並進変位指令Xc、Ycを制御している。   Therefore, in the low disturbance control device according to the first embodiment, the translation displacement command control system 20 is provided, and the translation displacement command control system 20 is based on the radial direction bearing force commands Fx and Fy and the axial direction rotor rotation angle θz. The radial translational displacement commands Xc and Yc are controlled so that the rotor center-of-gravity position G is located at the origin of the rotational coordinate system.

PID制御器15から出力されるラジアル方向軸受力指令Fx、Fy及びアキシャル方向ロータ回転角度θzは並進変位指令制御系20の回転座標変換器21に入力され、回転座標変換器21ではラジアル方向軸受力指令の回転座標系表現Fxr、Fyrを次式(1)、(2)により算出する。   The radial direction bearing force commands Fx and Fy and the axial direction rotor rotation angle θz output from the PID controller 15 are input to the rotational coordinate converter 21 of the translational displacement command control system 20, and the rotational coordinate converter 21 receives the radial direction bearing force. The rotational coordinate system representations Fxr and Fyr of the command are calculated by the following expressions (1) and (2).

Figure 0004344601
Figure 0004344601

このラジアル方向軸受力指令の回転座標系表現Fxr、Fyrは、通過低周波数帯域を大略1Hz程度の狭帯域に設定したLPF22に入力され、ここでラジアル方向軸受力指令のロータ回転同期成分のみが抽出される。LPF22の出力は、図4(a)においてベクトル31で表され、これが積分器23において積分された結果、ベクトル32で表されるラジアル方向並進変位指令の回転座標系表現Xcr、Ycrが積分器23から出力される。固定座標変換器24では、積分器23の出力であるラジアル方向並進変位指令の回転座標系表現Xcr、Ycrを、アキシャル方向ロータ回転角度θzを用いて次式(3)、(4)のように固定座標系表現のラジアル方向並進変位指令Xc、Ycに変換し、軸受制御系10に対して出力する。   The rotational coordinate system representations Fxr and Fyr of the radial bearing force command are input to the LPF 22 in which the low pass frequency band is set to a narrow band of about 1 Hz. Here, only the rotor rotation synchronization component of the radial bearing force command is extracted. Is done. The output of the LPF 22 is represented by a vector 31 in FIG. 4A, and as a result of integration by the integrator 23, the rotational coordinate system representations Xcr and Ycr of the radial translational displacement command represented by the vector 32 are integrated into the integrator 23. Is output from. In the fixed coordinate converter 24, the rotational coordinate system representation Xcr, Ycr of the radial translational displacement command, which is the output of the integrator 23, is expressed by the following equations (3), (4) using the axial direction rotor rotation angle θz. It is converted to a radial translational displacement command Xc, Yc expressed in a fixed coordinate system and output to the bearing control system 10.

Figure 0004344601
Figure 0004344601

こうして求められたラジアル方向並進変位指令Xc、Ycに追従するように軸受制御系10がロータ浮上位置を制御することにより、ロータ幾何中心位置Sは回転座標系上をベクトル32に従って移動し、定常状態では図4(b)に示すようにロータ重心位置Gが回転座標系原点に位置するようになる。これにより、ロータ静不釣合いに起因する振れ回り力30及びラジアル方向軸受力指令の回転座標系表現31は共に零ベクトルに収束し、従って、ロータ1に作用するラジアル方向軸受力及びステータに作用するラジアル方向軸受力反作用が共に低減した結果、ロータ静不釣合いに起因する擾乱力のレベルが低下する。   The bearing control system 10 controls the rotor flying position so as to follow the radial translational displacement commands Xc and Yc thus obtained, so that the rotor geometric center position S moves on the rotating coordinate system according to the vector 32 and is in a steady state. Then, as shown in FIG. 4B, the rotor center-of-gravity position G is located at the origin of the rotational coordinate system. As a result, both the run-out force 30 caused by the rotor static imbalance and the rotational coordinate system representation 31 of the radial bearing force command converge to the zero vector, and thus act on the radial bearing force acting on the rotor 1 and the stator. As a result of the reduction of both radial bearing force reactions, the level of disturbance forces due to rotor static imbalance is reduced.

このように、この実施の形態1による低擾乱化制御装置により、特にアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、ロータ重心位置Gが磁気軸受回転軸上に位置するようになり、ロータ静不釣合いに起因する擾乱力を低減することができる。   As described above, the low disturbance control device according to the first embodiment allows the rotor center-of-gravity position G to be positioned on the magnetic bearing rotating shaft, particularly when the axial direction rotor rotational speed is within the position control band of the bearing control system 10. Thus, the disturbance force caused by the rotor unbalance can be reduced.

また、並進変位指令制御系20におけるLPF22によって、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, the LPF 22 in the translational displacement command control system 20 can easily remove signal components belonging to a frequency region other than the rotor rotation synchronization component.

実施の形態2.
図5は、この発明に係る低擾乱化制御装置における実施の形態2を説明するためのブロック図、図6は、実施の形態2を説明するためのロータ半径方向断面図であり、図5及び図6において、図3及び図4と同一の符号を付したものは、同一部分または相当部分を示す。
Embodiment 2. FIG.
FIG. 5 is a block diagram for explaining the second embodiment of the low disturbance control device according to the present invention, and FIG. 6 is a rotor radial direction sectional view for explaining the second embodiment. 6, the same reference numerals as those in FIGS. 3 and 4 denote the same or corresponding parts.

図6(a)、(c)及び(e)は、アキシャル方向ロータ回転速度を磁気軸受のサーボ共振周波数域以上まで上昇させた場合において、固定座標変換器/位相変換器25における進相量を各々零値、最適値及び過大値として並進変位指令制御系を適用した直後を示し、図6(b)、(d)及び(f)は、各々図6(a)、(c)及び(e)において並進変位指令制御系を適用後十分時間が経過した後を示している。   6 (a), 6 (c) and 6 (e) show the amount of phase advance in the fixed coordinate converter / phase converter 25 when the rotational speed of the rotor in the axial direction is increased to the servo resonance frequency range of the magnetic bearing or higher. FIGS. 6 (b), (d), and (f) are respectively shown in FIGS. 6 (a), (c), and (e). ) After a sufficient time has elapsed after the translation displacement command control system is applied.

図5に示したように、この実施の形態2の低擾乱化制御装置は、位置制御ループにおけるPI制御器16、速度制御ループにおいてゲインKを施すP制御器17、ロータ変位に基づいてロータ速度を算出する微分器18、積分器23の出力であるXcr’、Ycr’を固定座標系表現に変換し、かつ、軸受制御系の位相特性を考慮してアキシャル方向ロータ回転速度に応じて所定の位相量を進めた信号を軸受制御系10のラジアル方向並進変位指令Xc、Ycとして出力する固定座標変換器/位相変換器25を備えている。 As shown in FIG. 5, lower perturbation of the control device of the second embodiment, PI controller 16 in the position control loop, P controller 17 to perform the gain K D in the speed control loop, based on the rotor displacement rotor Xcr ′ and Ycr ′, which are the outputs of the differentiator 18 and integrator 23 for calculating the speed, are converted into a fixed coordinate system expression, and predetermined according to the axial direction rotor rotational speed in consideration of the phase characteristics of the bearing control system. The fixed coordinate converter / phase converter 25 is provided which outputs a signal obtained by advancing the phase amount as radial translational displacement commands Xc and Yc of the bearing control system 10.

図6において、固定座標変換器/位相変換器25でアキシャル方向ロータ回転速度に応じて進める所定の位相量をα、ロータ回転方向を33、ロータ幾何中心位置Sの移動方向ベクトルを34とする。   In FIG. 6, a predetermined phase amount advanced according to the axial direction rotor rotational speed by the fixed coordinate converter / phase converter 25 is α, the rotor rotational direction is 33, and the moving direction vector of the rotor geometric center position S is 34.

次に、動作について説明する。軸受制御系10において、電磁石/ロータ/変位検出手段11から出力されるロータ変位情報Si(i=1、2、…、n)は、ロータ変位演算器12によってロータの3自由度並進変位X、Y、Z及びラジアル方向2自由度回転変位θx、θyからなる5自由度ロータ変位に変換される。   Next, the operation will be described. In the bearing control system 10, the rotor displacement information Si (i = 1, 2,..., N) output from the electromagnet / rotor / displacement detecting means 11 is converted by the rotor displacement calculator 12 into the three-degree-of-freedom translational displacement X, It is converted into a 5-degree-of-freedom rotor displacement comprising Y, Z and radial direction 2-degree-of-freedom rotational displacements θx, θy.

この5自由度ロータ変位はフィードバックされ、5自由度ロータ変位とラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycとの偏差信号が、位置制御ループにおけるPI制御器16に入力される。PI制御器16では入力された偏差信号に基づいて、磁気軸受に対する5自由度の速度指令を算出して出力する。   The 5-degree-of-freedom rotor displacement is fed back, and deviation signals between the 5-degree-of-freedom rotor displacement and the radial direction translational displacement commands Xc and Yc, the axial direction translational displacement command Zc, and the radial direction rotational displacement commands θxc and θyc are output in the position control loop. Input to the PI controller 16. The PI controller 16 calculates and outputs a 5-degree-of-freedom speed command for the magnetic bearing based on the input deviation signal.

ロータ変位演算器12から出力される5自由度ロータ変位は、更に、微分器18において微分され、5自由度ロータ速度に変換され、PI制御器16の出力である5自由度速度指令と微分器18の出力である5自由度のロータ速度から、速度制御ループにおける速度偏差信号が得られる。この速度偏差信号は速度制御ループにおけるP制御器17に入力され、速度制御ループゲインKを施した信号を磁気軸受に対する軸受力指令Fx、Fy、Fz及びラジアル方向軸受モーメント指令Nx、Nyとして出力する。 The 5-degree-of-freedom rotor displacement output from the rotor displacement calculator 12 is further differentiated by the differentiator 18 and converted into a 5-degree-of-freedom rotor speed, and the 5-degree-of-freedom speed command and the differentiator that are the outputs of the PI controller 16. A speed deviation signal in the speed control loop is obtained from the rotor speed of 5 degrees of freedom which is 18 outputs. This speed deviation signal is input to the P controller 17 in the speed control loop, outputs the signal subjected to the speed control loop gain K D bearing force command Fx with respect to the magnetic bearing, Fy, Fz and radial bearing moment command Nx, as Ny To do.

ギャップ長演算器13では、ロータ変位演算器12から出力される5自由度ロータ変位に基づいて、磁気軸受を構成するn個の電磁石の各々におけるギャップ長G1、G2、…、Gnを算出して出力し、電力分配器14では、ギャップ長演算器13から出力される軸受電磁石ギャップ長Gi(i=1、2、…、n)において、速度制御ループにおけるP制御器17から出力される軸受力指令Fx、Fy、Fz及びラジアル方向軸受モーメント指令Nx、Nyにより、n個の電磁石の各々に対する供給電流I1、I2、…、Inが分配して供給される。   The gap length calculator 13 calculates gap lengths G1, G2,..., Gn in each of the n electromagnets constituting the magnetic bearing based on the 5-degree-of-freedom rotor displacement output from the rotor displacement calculator 12. In the power distributor 14, the bearing force output from the P controller 17 in the speed control loop at the bearing electromagnet gap length Gi (i = 1, 2,..., N) output from the gap length calculator 13. Supply currents I1, I2,..., In are distributed and supplied to each of the n electromagnets according to the commands Fx, Fy, Fz and the radial bearing moment commands Nx, Ny.

ただし、図5では、簡単のためにアキシャル方向並進運動及びラジアル方向回転運動のフィードバック・ループは図示していない。   However, in FIG. 5, for the sake of simplicity, the feedback loop of the axial translational motion and the radial rotational motion is not shown.

このような軸受制御系10の動作により、ロータ浮上位置がラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycと一致するように制御することが可能となり、かつ軸受制御系10を位置制御ループと速度制御ループの2重ループ制御としたことで、図3に示した単一フィードバック・ループのPID制御器に比べて、位置制御閉ループゲイン特性が平坦となるように軸受制御系10を設計することができ、これによって磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱を抑制することができる。   By such an operation of the bearing control system 10, it is possible to control the rotor levitation position so as to coincide with the radial translational displacement commands Xc and Yc, the axial translational displacement command Zc, and the radial rotational displacement commands θxc and θyc. In addition, since the bearing control system 10 is a double loop control of a position control loop and a speed control loop, the position control closed loop gain characteristic is flatter than the single feedback loop PID controller shown in FIG. Thus, the bearing control system 10 can be designed so that the rotor swing amount in the servo resonance frequency region of the magnetic bearing and the disturbance of the magnetic bearing device can be suppressed.

また、並進変位指令制御系20の作用により、固定座標変換器/位相変換器25における進相量を零値としても、実施の形態1と同様にアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、ロータ重心位置Gが磁気軸受の回転軸上に位置するようになり、ロータ静不釣合いに起因する擾乱力を低減することができる。   Even if the phase advance amount in the fixed coordinate converter / phase converter 25 is set to zero by the action of the translational displacement command control system 20, the axial direction rotor rotational speed is the position of the bearing control system 10 as in the first embodiment. In the case of being within the control band, the rotor center-of-gravity position G is positioned on the rotating shaft of the magnetic bearing, so that the disturbance force due to the rotor unbalance can be reduced.

しかしながら、従来は、固定座標変換器/位相変換器25における進相量を零値にしており、固定座標変換器/位相変換器25における進相量を零値に保持したままアキシャル方向ロータ回転速度をサーボ共振周波数域まで上昇させると、ロータ幾何中心位置Sは、図6(a)に示したように、ベクトル32で表されるラジアル方向並進変位指令の回転座標系表現Xcr、Ycrに対して、軸受制御系10の位置制御閉ループ位相特性の遅れに対応した位相角だけロータ回転方向33と逆方向に回転したベクトル34に沿って、回転座標系上を移動するようになる。その結果、所定の時間が経過した後のロータ重心位置G及びロータ幾何中心位置Sは、図6(b)に示す配置となり、ロータ静不釣合いに起因する振れ回り力及びラジアル方向軸受力指令の回転座標系表現が各々ベクトル30’及び31’となる結果、ラジアル方向軸受力指令の回転座標系表現ベクトル31はロータ回転方向33と同一方向に回転し、ロータ重心位置Gと磁気軸受回転軸との距離が零値に収束しないために、ロータ静不釣合いに起因する擾乱力を低減することができなくなる。   However, conventionally, the phase advance amount in the fixed coordinate converter / phase converter 25 is set to zero, and the axial direction rotor rotation speed is maintained while the phase advance amount in the fixed coordinate converter / phase converter 25 is held at zero value. Is increased to the servo resonance frequency range, the rotor geometric center position S is obtained with respect to the rotational coordinate system expressions Xcr and Ycr of the radial translational displacement command represented by the vector 32 as shown in FIG. The bearing control system 10 moves on the rotational coordinate system along the vector 34 rotated in the direction opposite to the rotor rotation direction 33 by a phase angle corresponding to the delay of the position control closed loop phase characteristic. As a result, the rotor center-of-gravity position G and the rotor geometric center position S after a lapse of a predetermined time are arranged as shown in FIG. 6B, and the whirling force and radial direction bearing force command caused by the rotor unbalance are obtained. As a result of the rotational coordinate system representations being vectors 30 ′ and 31 ′, the rotational coordinate system representation vector 31 of the radial bearing force command rotates in the same direction as the rotor rotational direction 33, and the rotor center-of-gravity position G, the magnetic bearing rotational axis, Therefore, the disturbance force due to the rotor unbalance cannot be reduced.

そこで、この実施の形態2の低擾乱化制御装置では、並進変位指令制御系20において、ラジアル方向軸受力指令Fx、Fy、アキシャル方向ロータ回転角度θz及び軸受制御系10の位置制御閉ループ位相特性に基づいて、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ重心位置Gが回転座標系原点に位置するようにラジアル方向並進変位指令Xc、Ycを制御している。   Thus, in the low disturbance control device of the second embodiment, in the translational displacement command control system 20, the radial direction bearing force commands Fx and Fy, the axial direction rotor rotation angle θz, and the position control closed loop phase characteristic of the bearing control system 10 are used. Based on this, the rotor center-of-gravity position G is located at the origin of the rotational coordinate system regardless of the axial direction rotor rotational speed in any frequency region within the position control band of the bearing control system 10, the servo resonance frequency region, and the position control band. The radial translation displacement commands Xc and Yc are controlled.

P制御器17から出力されるラジアル方向軸受力指令Fx、Fyは、アキシャル方向ロータ回転角度θzに基づいて、回転座標変換器21により、上記式(1)及び(2)に示した回転座標系表現Fxr、Fyrに変換される。このラジアル方向軸受力指令の回転座標系表現Fxr、Fyrは、通過低周波数帯域を大略1Hz程度の狭帯域に設定したLPF22によってロータ回転同期成分のみが抽出される。LPF22の出力は、図6(c)においてベクトル31で表され、これが積分器23において積分され、積分されて得られた信号Xcr’、Ycr’が固定座標変換器/位相変換器25に入力される。   The radial direction bearing force commands Fx and Fy output from the P controller 17 are rotated by the rotation coordinate converter 21 based on the axial direction rotor rotation angle θz and the rotation coordinate system shown in the above formulas (1) and (2). The expressions Fxr and Fyr are converted. In the rotational coordinate system representations Fxr and Fyr of the radial bearing force command, only the rotor rotation synchronization component is extracted by the LPF 22 in which the low pass frequency band is set to a narrow band of about 1 Hz. The output of the LPF 22 is represented by a vector 31 in FIG. 6C. This is integrated by the integrator 23, and the signals Xcr ′ and Ycr ′ obtained by the integration are input to the fixed coordinate converter / phase converter 25. The

固定座標変換器/位相変換器25では、積分器23の出力Xcr’、Ycr’を、アキシャル方向ロータ回転角度θzを用いて固定座標系表現に変換するが、このとき、軸受制御系10の位置制御閉ループ位相特性を考慮して、次式(5)及び(6)に示すようにアキシャル方向ロータ回転速度に応じて所定の位相量αを進めた上で固定座標系表現に変換した信号を、ラジアル方向並進変位指令Xc、Ycとして軸受制御系10に対して出力する。   In the fixed coordinate converter / phase converter 25, the outputs Xcr ′ and Ycr ′ of the integrator 23 are converted into a fixed coordinate system expression using the axial direction rotor rotation angle θz. At this time, the position of the bearing control system 10 is changed. Considering the control closed-loop phase characteristics, a signal converted into a fixed coordinate system expression after a predetermined phase amount α is advanced according to the axial direction rotor rotational speed as shown in the following equations (5) and (6), It outputs to the bearing control system 10 as radial direction translational displacement command Xc, Yc.

Figure 0004344601
Figure 0004344601

このラジアル方向並進変位指令Xc、Ycの回転座標系表現Xcr、Ycrは、図6(c)においてベクトル32で表され、ロータ幾何中心位置Sはベクトル32に対して、軸受制御系10の位置制御閉ループ位相特性の遅れに対応した位相角だけ、ロータ回転方向33と逆方向に回転したベクトル34に沿って回転座標系上を移動するため、定常状態では図6(d)に示すようにロータ重心位置Gが回転座標系原点に位置するようになる。   The rotational coordinate system representations Xcr and Ycr of the radial translational displacement commands Xc and Yc are represented by a vector 32 in FIG. 6C, and the rotor geometric center position S is controlled by the position control of the bearing control system 10 with respect to the vector 32. Since it moves on the rotating coordinate system along the vector 34 rotated in the direction opposite to the rotor rotation direction 33 by the phase angle corresponding to the delay of the closed loop phase characteristic, in the steady state, as shown in FIG. The position G is positioned at the origin of the rotating coordinate system.

これにより、ロータ静不釣合いに起因する振れ回り力30及びラジアル方向軸受力指令の回転座標系表現31は共に零ベクトルに収束し、従って、ロータに作用するラジアル方向軸受力及びステータに作用するラジアル方向軸受力反作用が共に低減した結果、ロータ静不釣合いに起因する擾乱力のレベルが低下する。   As a result, both the run-out force 30 caused by the rotor static imbalance and the rotational coordinate system representation 31 of the radial bearing force command converge to the zero vector, and therefore the radial bearing force acting on the rotor and the radial acting on the stator. As a result of the reduction in both directional bearing force reactions, the level of disturbance forces due to rotor static imbalance is reduced.

ここで、固定座標変換器/位相変換器25における進相量αを、図6(e)に示したように過大値に設定した場合、ロータ幾何中心位置Sの移動方向ベクトル34は、ラジアル方向軸受力指令ベクトル31に対してロータ回転方向33に所定位相角だけ回転した配置となる。   Here, when the phase advance amount α in the fixed coordinate converter / phase converter 25 is set to an excessive value as shown in FIG. 6E, the moving direction vector 34 of the rotor geometric center position S is expressed in the radial direction. The arrangement is such that the bearing force command vector 31 is rotated by a predetermined phase angle in the rotor rotation direction 33.

その結果、所定の時間が経過した後のロータ重心位置G及びロータ幾何中心位置Sは、図6(f)に示す配置となり、ロータ静不釣合いに起因する振れ回り力及びラジアル方向軸受力指令の回転座標系表現が各々30’及び31’で表されるベクトルとなる結果、ラジアル方向軸受力指令の回転座標系表現ベクトル31はロータ回転方向33と逆方向に回転する。   As a result, the rotor center-of-gravity position G and the rotor geometric center position S after a lapse of a predetermined time are arranged as shown in FIG. 6 (f), and the whirling force and radial direction bearing force command caused by the rotor unbalance are obtained. As a result of the rotational coordinate system representations being vectors represented by 30 ′ and 31 ′, the rotational coordinate system representation vector 31 of the radial bearing force command rotates in the direction opposite to the rotor rotational direction 33.

すなわち、ラジアル方向軸受力指令の回転座標系表現ベクトル31は、固定座標変換器/位相変換器25における進相量αが最適値に対して小さい場合(図6(a)、(b))はロータ回転方向33と同一方向に回転し、最適値に対して大きい場合(図6(e)、(f))はロータ回転方向33と逆方向に回転し、最適値の場合(図6(c)、(d))はベクトルの方向は変化せずに零ベクトルに収束する。   That is, the rotational coordinate system expression vector 31 of the radial bearing force command is obtained when the phase advance amount α in the fixed coordinate converter / phase converter 25 is smaller than the optimum value (FIGS. 6A and 6B). When rotating in the same direction as the rotor rotation direction 33 and larger than the optimum value (FIGS. 6 (e) and (f)), it rotates in the opposite direction to the rotor rotation direction 33, and in the case of the optimum value (FIG. 6 (c) ), (D)) converges to a zero vector without changing the direction of the vector.

従って、どのようなアキシャル方向ロータ回転速度においても、固定座標変換器/位相変換器25における進相量αをラジアル方向軸受力指令の回転座標系表現ベクトル31が回転せずに零ベクトルに収束するような値に設定することにより、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ静不釣合いに起因する振れ回り力30及びラジアル方向軸受力指令の回転座標系表現31を共に零ベクトルに収束させることができ、ロータに作用するラジアル方向軸受力及びステータに作用するラジアル方向軸受力反作用が共に低減するため、ロータ静不釣合いに起因する擾乱力のレベルがサーボ共振周波数域を含めた全周波数域において低下するようになる。   Accordingly, the phase advance amount α in the fixed coordinate converter / phase converter 25 converges to the zero vector without rotating the rotational coordinate system expression vector 31 of the radial bearing force command at any axial rotor rotational speed. By setting to such a value, the axial direction rotor rotational speed is caused by the rotor static imbalance regardless of the frequency range within the position control band, the servo resonance frequency range, and the position control band of the bearing control system 10. The swinging force 30 and the rotational coordinate system representation 31 of the radial bearing force command can be converged to the zero vector, and both the radial bearing force acting on the rotor and the radial bearing force reaction acting on the stator are reduced. , The level of disturbance force due to rotor unbalance decreases in the entire frequency range including the servo resonance frequency range. So as to.

このように、この実施の形態2の低擾乱化制御装置により、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ重心位置Gが磁気軸受の回転軸上に位置するようになり、ロータ静不釣合いに起因する擾乱力を、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減することができる。   As described above, the low disturbance control device according to the second embodiment allows the axial direction rotor rotational speed to be in any frequency region within the position control band, the servo resonance frequency range, and the position control band of the bearing control system 10. The rotor center-of-gravity position G is positioned on the rotating shaft of the magnetic bearing, and the disturbance force caused by the rotor unbalance can be reduced in the entire frequency region including the servo resonance frequency region of the magnetic bearing.

また、並進変位指令制御系20におけるLPF22がロータ回転同期成分のみを抽出することで、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, since the LPF 22 in the translational displacement command control system 20 extracts only the rotor rotation synchronization component, signal components belonging to a frequency region other than the rotor rotation synchronization component can be easily removed.

更に、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御とし、位置制御ループゲイン特性が平坦となるように軸受制御系10を設計することにより、磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱の更なる抑制を図ることができる。   Furthermore, the bearing control system 10 is configured as a double loop control of a position control loop and a speed control loop, and the bearing control system 10 is designed so that the position control loop gain characteristic becomes flat, so that in the servo resonance frequency region of the magnetic bearing. It is possible to further suppress the amount of rotor swing and disturbance of the magnetic bearing device.

実施の形態3.
図7は、この発明に係る実施の形態3の低擾乱化制御装置を説明するためのブロック図、図8は、閉ループゲイン特性を示したボード線図であり、図9、図10及び図11は、この実施の形態3の低擾乱化制御装置を説明するためのロータ断面図であり、図7において、図3及び図5と同一の符号を付したものは、同一部分または相当部分を示し、図9、図10及び図11において、図2、図4及び図6と同一の符号を付したものは、同一部分または相当部分を示す。
Embodiment 3 FIG.
FIG. 7 is a block diagram for explaining the low disturbance control device according to the third embodiment of the present invention, and FIG. 8 is a Bode diagram showing the closed-loop gain characteristics, and FIGS. 9, 10, and 11 are used. These are rotor sectional drawings for demonstrating the low disturbance control apparatus of this Embodiment 3, In FIG. 7, what attached | subjected the code | symbol same as FIG.3 and FIG.5 shows the same part or an equivalent part. 9, 10, and 11, the same reference numerals as those in FIGS. 2, 4, and 6 denote the same or corresponding parts.

図7に示したように、この実施の形態3の低擾乱化制御装置は、軸受制御系10におけるラジアル方向軸受モーメント指令Nx、Nyから、軸受制御系10に対するラジアル方向回転変位指令を算出して出力する回転変位指令制御系40を備える。   As shown in FIG. 7, the low disturbance control device according to the third embodiment calculates a radial rotational displacement command for the bearing control system 10 from the radial bearing moment commands Nx and Ny in the bearing control system 10. A rotational displacement command control system 40 for outputting is provided.

回転変位指令制御系40は、軸受制御系10におけるラジアル方向軸受モーメント指令Nx、Nyを、アキシャル方向ロータ回転角度θzを用いてロータに固定された回転座標系表現Nxr、Nyrに変換する回転座標変換器41、軸受モーメント指令の回転座標系表現Nxr、Nyrに対して、通過低周波数帯域を大略1Hz程度の狭帯域とすることで、軸受モーメント指令のロータ回転同期成分のみを抽出するLPF(ローパス・フィルタ)42、LPF42の出力を積分してラジアル方向回転変位の信号θxcr’、θycr’を出力する積分器43、信号θxcr’、θycr’を固定座標系表現に変換し、かつ軸受制御系10の位相特性を考慮してアキシャル方向ロータ回転速度に応じて所定の位相量を進めた信号を軸受制御系10のラジアル方向回転変位指令θxc、θycとして出力する固定座標変換器/位相変換器45を備える。   The rotational displacement command control system 40 converts the radial bearing moment commands Nx, Ny in the bearing control system 10 into rotational coordinate system representations Nxr, Nyr fixed to the rotor using the axial rotor rotational angle θz. LPF (low-pass filter), which extracts only the rotor rotation synchronization component of the bearing moment command, by setting the low pass frequency band to a narrow band of about 1 Hz with respect to the rotating coordinate system representation Nxr, Nyr of the bearing 41 Filter) 42, integrator 43 that integrates the outputs of LPF 42 and outputs radial rotational displacement signals θxcr ′ and θycr ′, converts signals θxcr ′ and θycr ′ into a fixed coordinate system representation, and In consideration of the phase characteristics, a signal obtained by advancing a predetermined phase amount in accordance with the axial rotor rotational speed is transmitted to the bearing control system 10. A fixed coordinate converter / phase converter 45 that outputs radial rotational displacement commands θxc and θyc is provided.

また、図8において、2重ループ制御で構成した軸受制御系10の位置制御閉ループゲイン特性を50、軸受制御系10の速度制御閉ループゲイン特性を51、位置制御周波数帯域をf、速度制御周波数帯域をf、アキシャル方向ロータ回転速度をf、速度制御ループゲイン切換閾値をfとしている。 Further, in FIG. 8, the position control loop gain characteristic 50 of the bearing control system 10 constructed in double loop control, the speed control loop gain characteristic of the bearing control system 10 51, position control frequency band f p, speed control frequency The band is f D , the axial rotor rotational speed is f, and the speed control loop gain switching threshold is f 1 .

図9(a)、(b)は円柱状ロータ(図1)の磁気軸受装置でアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、回転変位指令制御系40を適用した直後を示し、図9(c)、(d)は回転変位指令制御系40を適用後十分時間が経過した後を示し、図10(a)、(b)は円板状ロータ(図2)の磁気軸受装置でアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、回転変位指令制御系40を適用した直後を示し、図10(c)、(d)は回転変位指令制御系40を適用後十分時間が経過した後を示している。   9A and 9B show a magnetic bearing device of a cylindrical rotor (FIG. 1), and the rotational displacement command control system 40 is applied when the axial direction rotor rotational speed is within the position control band of the bearing control system 10. 9 (c) and 9 (d) show a state after sufficient time has elapsed after applying the rotational displacement command control system 40, and FIGS. 10 (a) and 10 (b) show a disk-shaped rotor (FIG. 2). ) In the case where the axial direction rotor rotational speed is within the position control band of the bearing control system 10, immediately after applying the rotational displacement command control system 40, FIGS. 10C and 10D show the rotation. It shows after a sufficient time has elapsed after applying the displacement command control system 40.

また、図9〜図11において、アキシャル方向のロータ慣性主軸をP、アキシャル方向ロータ慣性主軸P上の固定点をa、ロータを2分割して各部分を質点で置き換えた場合に、各質点に作用する遠心力をF、ロータの動不釣合いに起因する振れ回りモーメントを60、回転座標変換器41から出力されるラジアル方向軸受モーメント指令の回転座標系表現Nxr、Nyrを要素とするモーメントを61、ラジアル方向回転変位指令の回転座標系表現を62としている。   Further, in FIGS. 9 to 11, when the axial rotor inertia main shaft is P, the fixed point on the axial rotor inertia main shaft P is a, the rotor is divided into two parts and each part is replaced with a mass point, The acting centrifugal force is F, the swinging moment due to the rotor dynamic imbalance is 60, and the moment having the rotational coordinate system expression Nxr, Nyr of the radial bearing moment command output from the rotational coordinate converter 41 as an element is 61. The rotational coordinate system representation of the radial direction rotational displacement command is 62.

図11(a)、(c)及び(e)は、アキシャル方向ロータ回転速度を磁気軸受のサーボ共振周波数域以上まで上昇させた場合において、固定座標変換器/位相変換器45における進相量を各々零値、最適値及び過大値として回転変位指令制御系40を適用した直後を示し、図11(b)、(d)及び(f)は、各々図11(a)、(c)及び(e)において回転変位指令制御系40を適用後十分時間が経過した後のロータ半径方向断面図である。   11 (a), 11 (c) and 11 (e) show the amount of phase advance in the fixed coordinate converter / phase converter 45 when the axial rotor rotational speed is increased to the servo resonance frequency range of the magnetic bearing or higher. FIGS. 11 (b), (d) and (f) show FIGS. 11 (a), 11 (c) and (f), respectively, immediately after applying the rotational displacement command control system 40 as a zero value, an optimum value and an excessive value. 7 is a radial cross-sectional view of the rotor after a sufficient time has elapsed after applying the rotational displacement command control system 40 in e). FIG.

更に、図11において、固定座標変換器/位相変換器45における、アキシャル方向ロータ回転速度に応じて進める所定の位相量をβ、ロータ回転方向を63、アキシャル方向ロータ慣性主軸Pの回転方向ベクトルを64としている。   Further, in FIG. 11, in the fixed coordinate converter / phase converter 45, a predetermined phase amount advanced according to the axial direction rotor rotational speed is β, the rotor rotational direction is 63, and the rotational direction vector of the axial direction rotor inertial main axis P is 64.

次に、動作について説明する。図7では、簡単のためにラジアル方向並進運動及びアキシャル方向並進運動のフィードバック・ループは図示していないが、軸受制御系10の動作は実施の形態2と同一である。このとき、軸受制御系10の位置制御閉ループゲイン特性50及び速度制御閉ループゲイン特性51は、図8(a)に示した通りであり、速度制御周波数帯域fが位置制御周波数帯域fよりも十分高くなるように軸受制御系10を設計する。このような設計とすることで位置制御ループゲイン特性が平坦となり、磁気軸受のサーボ共振周波数域(位置制御周波数帯域f)におけるロータ振れ回り量及び磁気軸受装置の発生擾乱を抑制することができるが、その反面、速度制御系ゲイン交点周波数域(速度制御周波数帯域f)における磁気軸受装置の発生擾乱が多少増加する傾向を示す。 Next, the operation will be described. In FIG. 7, for the sake of simplicity, the feedback loop of the radial translation and the axial translation is not shown, but the operation of the bearing control system 10 is the same as that of the second embodiment. At this time, the position control loop gain characteristic 50 and a speed control loop gain characteristic 51 of the bearing control system 10 is as shown in FIG. 8 (a), the speed control frequency band f D than the position control frequency band f P The bearing control system 10 is designed to be sufficiently high. By adopting such a design, the position control loop gain characteristic becomes flat, and the rotor swing amount in the servo resonance frequency region (position control frequency band f P ) of the magnetic bearing and the disturbance generated in the magnetic bearing device can be suppressed. On the other hand, however, the disturbance generated by the magnetic bearing device in the speed control system gain intersection frequency range (speed control frequency band f D ) tends to increase somewhat.

そこで、この実施の形態3の低擾乱化制御装置では、速度制御ループゲイン切換閾値fを設定し、アキシャル方向ロータ回転速度fがfよりも低い場合は、図8(a)に示したように、ロータ回転速度f<f< 速度制御周波数帯域f、となるように軸受制御系10の速度制御ループゲインKを増加させ、アキシャル方向ロータ回転速度fがfよりも高い場合は、図8(b)に示したように、速度制御周波数帯域f < f < ロータ回転速度fとなるように軸受制御系10の速度制御ループゲインKを減少させることで、速度制御周波数帯域fがアキシャル方向ロータ回転速度fの近傍に属することが無い構成とし、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域(速度制御周波数帯域f)における磁気軸受装置の発生擾乱を抑制することができる。 Therefore, in the low disturbance control device according to the third embodiment, when the speed control loop gain switching threshold f 1 is set and the axial direction rotor rotational speed f is lower than f 1 , FIG. As described above, when the speed control loop gain K D of the bearing control system 10 is increased so that the rotor rotational speed f <f 1 <speed control frequency band f D , and the axial direction rotor rotational speed f is higher than f 1. As shown in FIG. 8B, the speed control loop gain K D of the bearing control system 10 is decreased so that the speed control frequency band f D <f 1 <rotor rotational speed f is satisfied, thereby speed control. frequency band f D is configured not to belong to the vicinity of the axial rotor rotational speed f, without impairing the stability of the bearing control system 10, the speed control system gain crossover frequency range (speed control frequency band it is possible to suppress the occurrence disturbance of the magnetic bearing apparatus in f D).

このような軸受制御系10の動作により、ロータ浮上位置がラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycと一致するように制御することが可能となり、かつ、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御としたことで、図3に示した単一フィードバック・ループのPID制御器に比べて、位置制御閉ループゲイン特性が平坦となるように軸受制御系10を設計することができ、これによって磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱を抑制することができる。   By such an operation of the bearing control system 10, it is possible to control the rotor levitation position so as to coincide with the radial translational displacement commands Xc and Yc, the axial translational displacement command Zc, and the radial rotational displacement commands θxc and θyc. In addition, since the bearing control system 10 is a double loop control of a position control loop and a speed control loop, the position control closed loop gain characteristic is flatter than the single feedback loop PID controller shown in FIG. Thus, the bearing control system 10 can be designed so that the rotor swing amount in the servo resonance frequency region of the magnetic bearing and the disturbance of the magnetic bearing device can be suppressed.

加えて、軸受制御系10において速度制御周波数帯域fがアキシャル方向ロータ回転速度fの近傍に属することが無いように、ロータ回転速度fに応じて速度制御ループゲインKを変化させることにより、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域における磁気軸受装置の発生擾乱の更なる抑制を図ることができる。 In addition, by changing the speed control loop gain K D according to the rotor rotational speed f so that the speed control frequency band f D does not belong to the vicinity of the axial direction rotor rotational speed f in the bearing control system 10, Without impairing the stability of the bearing control system 10, it is possible to further suppress the disturbance of the magnetic bearing device in the speed control system gain intersection frequency range.

従来の磁気軸受制御装置では、軸受制御系10に対するラジアル方向回転変位指令θxc、θycを常に零値としており、従ってアキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域よりも低い場合のロータ断面図は、理想的には円柱状ロータ(図1)で図9(a)、(b)に示した通りとなり、円板状ロータ(図2)で図10(a)、(b)に示した通りとなる。つまり、ロータ回転軸(Z軸)に垂直に直交し、かつ、ロータに対して固定された回転座標系Xr−Yrから見た場合、ロータ変位演算器12で算出されるラジアル方向2自由度回転変位θx、θyは零値に保持され、アキシャル方向ロータ慣性主軸Pはロータ回転軸(Z軸)に対してロータ動不釣合いに対応する角度だけ傾いた配置となる。ただし図9及び10では簡単のために、アキシャル方向ロータ慣性主軸PがZr−Xr平面内に位置する場合を記載している。   In the conventional magnetic bearing control device, the radial direction rotational displacement commands θxc and θyc for the bearing control system 10 are always zero, and therefore the rotor cross section when the axial direction rotor rotational speed is lower than the position control band of the bearing control system 10. The figure is ideally shown in FIGS. 9 (a) and 9 (b) for a cylindrical rotor (FIG. 1), and shown in FIGS. 10 (a) and 10 (b) for a disk-like rotor (FIG. 2). It will be true. That is, when viewed from the rotational coordinate system Xr-Yr perpendicular to the rotor rotation axis (Z axis) and fixed to the rotor, the rotation in the radial direction calculated by the rotor displacement calculator 12 is two degrees of freedom. The displacements θx and θy are held at zero values, and the axial direction rotor inertia main shaft P is inclined with respect to the rotor rotation axis (Z axis) by an angle corresponding to the rotor motion imbalance. However, in FIGS. 9 and 10, for the sake of simplicity, the case where the axial direction rotor inertia main axis P is located in the Zr-Xr plane is described.

このとき、ロータを2分割して各部分を質点で置き換えた場合の、各質点に作用する遠心力Fは、図9(a)、図10(a)に示した通りとなり、その結果、ロータには動不釣合いに起因する振れ回りモーメント60が作用する。振れ回りモーメント60の作用下でロータ浮上位置、すなわちラジアル方向2自由度回転変位θx、θyを零値に保持するために、ラジアル方向軸受モーメント指令の回転座標系表現Nxr、Nyrを要素とするモーメント61は、振れ回りモーメント60と逆極性となる。このとき、ロータに作用するラジアル方向軸受モーメントの反作用がステータに作用し、これが磁気軸受装置のハウジングに伝達するモーメントの振動、すなわち擾乱の主要因となる。   At this time, when the rotor is divided into two parts and each part is replaced with a mass point, the centrifugal force F acting on each mass point is as shown in FIG. 9 (a) and FIG. 10 (a). The swinging moment 60 caused by the dynamic imbalance acts on. In order to maintain the rotor floating position, that is, the rotational displacement θx, θy in the radial direction with two degrees of freedom under the action of the swinging moment 60 at a zero value, the moment having the rotational coordinate system expression Nxr, Nyr of the radial direction bearing moment command as an element. 61 has a reverse polarity to the swinging moment 60. At this time, the reaction of the radial bearing moment acting on the rotor acts on the stator, and this becomes the main factor of vibration of the moment transmitted to the housing of the magnetic bearing device, that is, disturbance.

そこで、この実施の形態3の低擾乱化制御装置では、回転変位指令制御系40において、第1に、ラジアル方向軸受モーメント指令Nx、Ny及びアキシャル方向ロータ回転角度θzに基づいて、ロータ慣性主軸Pが磁気軸受回転軸と平行になるようにラジアル方向回転変位指令θxc、θycを制御している。   Thus, in the low disturbance control device of the third embodiment, in the rotational displacement command control system 40, first, the rotor inertia main shaft P is based on the radial direction bearing moment commands Nx, Ny and the axial direction rotor rotation angle θz. Are controlled so as to be parallel to the rotation axis of the magnetic bearing so that the radial rotational displacement commands θxc and θyc are controlled.

P制御器17から出力されるラジアル方向軸受モーメント指令Nx、Nyは、アキシャル方向ロータ回転角度θzに基づいて、回転座標変換器41により次式(7)、(8)のように回転座標系表現Nxr、Nyrに変換される。   The radial bearing moment commands Nx and Ny output from the P controller 17 are expressed in the rotational coordinate system as shown in the following equations (7) and (8) by the rotational coordinate converter 41 based on the axial rotor rotational angle θz. Converted to Nxr and Nyr.

Figure 0004344601
Figure 0004344601

このラジアル方向軸受モーメント指令の回転座標系表現Nxr、Nyrは、通過低周波数帯域を大略1Hz程度の狭帯域に設定したLPF42によってロータ回転同期成分のみが抽出される。LPF42の出力は、図9(a)、(b)及び図10(a)、(b)においてモーメント61で表され、このモーメント61が円柱状ロータでは同極性(図9(b))で、円板状ロータでは逆極性(図10(b))で積分器43において積分された結果、回転ベクトル62で表されるラジアル方向回転変位指令の回転座標系表現θxcr’、θycr’が積分器43から出力される。固定座標変換器/位相変換器45では、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、積分器43の出力θxcr’、θycr’を、アキシャル方向ロータ回転角度θzを用いて次式(9)、(10)のように固定座標系表現θxc、θycに変換し、軸受制御系10に対して出力する。   In the rotational coordinate system representations Nxr and Nyr of the radial bearing moment command, only the rotor rotation synchronization component is extracted by the LPF 42 in which the low pass frequency band is set to a narrow band of about 1 Hz. The output of the LPF 42 is represented by a moment 61 in FIGS. 9A and 9B and FIGS. 10A and 10B, and this moment 61 has the same polarity (FIG. 9B) in the cylindrical rotor. In the disk-shaped rotor, as a result of integration in the integrator 43 with the reverse polarity (FIG. 10B), the rotational coordinate system representations θxcr ′ and θycr ′ of the radial direction rotational displacement command represented by the rotation vector 62 are integrated in the integrator 43. Is output from. In the fixed coordinate converter / phase converter 45, when the axial rotor rotational speed is within the position control band of the bearing control system 10, the outputs θxcr ′ and θycr ′ of the integrator 43 are set to the axial rotor rotational angle θz. By using them, they are converted into fixed coordinate system expressions θxc and θyc as in the following equations (9) and (10), and output to the bearing control system 10.

Figure 0004344601
Figure 0004344601

このようにして求められたラジアル方向回転変位指令θxc、θycに追従するように軸受制御系10がロータ浮上位置を制御することにより、アキシャル方向ロータ慣性主軸Pは回転座標系上を回転ベクトル62に従って回転し、定常状態では図9(c)、(d)及び図10(c)、(d)に示すようにロータ慣性主軸Pが磁気軸受中心軸に対して平行となる。これにより、ロータ動不釣合いに起因する振れ回りモーメント60及びラジアル方向軸受モーメント指令の回転座標系表現61は共に零値に収束し、従ってロータに作用するラジアル方向軸受モーメント及びステータに作用するラジアル方向軸受モーメント反作用が共に低減した結果、ロータ動不釣合いに起因する擾乱モーメントのレベルが低下する。   The bearing control system 10 controls the rotor floating position so as to follow the radial direction rotational displacement commands θxc and θyc thus determined, whereby the axial direction rotor inertia main shaft P follows the rotational vector 62 on the rotational coordinate system. In a steady state, the rotor inertia main shaft P is parallel to the magnetic bearing central axis as shown in FIGS. 9C, 9D, 10C, and 10D. Thus, the rotational coordinate system representation 61 of the run-out moment 60 and the radial bearing moment command due to the rotor dynamic imbalance converges to zero values, and thus the radial bearing moment acting on the rotor and the radial direction acting on the stator. As a result of both bearing moment reactions being reduced, the level of disturbance moments due to rotor dynamic imbalance is reduced.

ラジアル方向回転運動では上述のように、ロータ動不釣合いに起因する振れ回りモーメント60、従ってラジアル方向軸受モーメント指令61の極性が、ロータのアキシャル方向慣性モーメントとラジアル方向慣性モーメントとの大小関係により反転する。従って、ラジアル方向回転変位指令の回転座標系表現62を算出するために積分器43で実行する積分動作の極性が、ラジアル方向慣性モーメント>アキシャル方向慣性モーメントの場合(円柱状ロータ)では同極性、アキシャル方向慣性モーメント>ラジアル方向慣性モーメントの場合(円板状ロータ)では逆極性となることに注意を要する。   As described above, in the rotational motion in the radial direction, the swinging moment 60 caused by the rotor dynamic imbalance, and hence the polarity of the radial bearing moment command 61, is reversed by the magnitude relationship between the axial inertia moment and the radial inertia moment of the rotor. To do. Accordingly, when the polarity of the integral operation executed by the integrator 43 to calculate the rotational coordinate system representation 62 of the radial direction rotational displacement command is radial moment of inertia> axial moment of inertia (cylindrical rotor), the polarity is Note that if the moment of inertia in the axial direction is greater than the moment of inertia in the radial direction (disk-shaped rotor), the polarity is reversed.

このように、回転変位指令制御系40の作用により、固定座標変換器/位相変換器45における進相量を零値としても、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内にある場合において、アキシャル方向ロータ慣性主軸Pが磁気軸受回転軸に対して平行となり、ロータ動不釣合いに起因する擾乱モーメントを低減することができる。   As described above, due to the action of the rotational displacement command control system 40, the axial direction rotor rotational speed is within the position control band of the bearing control system 10 even if the phase advance amount in the fixed coordinate converter / phase converter 45 is zero. In this case, the axial direction rotor inertia main shaft P is parallel to the magnetic bearing rotating shaft, and the disturbance moment caused by rotor dynamic imbalance can be reduced.

しかしながら、固定座標変換器/位相変換器45における進相量を零値に保持したままアキシャル方向ロータ回転速度をサーボ共振周波数域まで上昇させると、ロータ慣性主軸Pは、図11(a)に示すように、回転ベクトル62で表されるラジアル方向回転変位指令の回転座標系表現θxcr’、θycr’に対して、軸受制御系10の位置制御閉ループ位相特性の遅れに対応した位相角だけロータ回転方向63と逆方向に回転した回転ベクトル64に従って、回転座標系上を回転するようになる。その結果、所定の時間が経過した後のロータ慣性主軸P上の固定点aは図11(b)に示す配置となり、ロータ動不釣合いに起因する振れ回りモーメント及びラジアル方向軸受モーメント指令の回転座標系表現が各々60’及び61’で表されるモーメントとなる結果、ラジアル方向軸受モーメント指令の回転座標系表現61はロータ回転方向63と同一方向に回転し、ロータ慣性主軸Pと磁気軸受回転軸との角度が零値に収束しないために、ロータ動不釣合いに起因する擾乱モーメントを低減することができなくなる。   However, when the axial direction rotor rotation speed is increased to the servo resonance frequency range while the phase advance amount in the fixed coordinate converter / phase converter 45 is kept at zero, the rotor inertia main axis P is shown in FIG. Thus, with respect to the rotational coordinate system representation θxcr ′, θycr ′ of the radial direction rotational displacement command represented by the rotation vector 62, the rotor rotational direction is equal to the phase angle corresponding to the delay of the position control closed loop phase characteristic of the bearing control system 10. According to the rotation vector 64 rotated in the opposite direction to 63, the rotation coordinate system is rotated. As a result, the fixed points a on the rotor inertia main shaft P after a predetermined time have passed are arranged as shown in FIG. 11B, and the rotational coordinates of the swinging moment and the radial bearing moment command caused by the rotor dynamic imbalance. As a result of the system representations being 60 ′ and 61 ′, the rotational coordinate system representation 61 of the radial bearing moment command rotates in the same direction as the rotor rotational direction 63, and the rotor inertia main shaft P and the magnetic bearing rotational shaft Since the angle with the angle does not converge to zero, the disturbance moment due to rotor dynamic imbalance cannot be reduced.

そこで、この実施の形態3の低擾乱化制御装置では、回転変位指令制御系40において、軸受制御系10の位置制御閉ループ位相特性に基づいて、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ慣性主軸Pが磁気軸受回転軸と平行になるように固定座標変換器/位相変換器45における進相量を調整し、ラジアル方向回転変位指令θxc、θycを制御している。   Therefore, in the low disturbance control device according to the third embodiment, in the rotational displacement command control system 40, the axial direction rotor rotational speed is controlled based on the position control closed loop phase characteristic of the bearing control system 10 according to the position control of the bearing control system 10. The phase advance amount in the fixed coordinate converter / phase converter 45 is set so that the rotor inertial main shaft P is parallel to the magnetic bearing rotation axis in any frequency region outside the band, the servo resonance frequency region, and the position control band. It adjusts and controls radial direction rotational displacement commands θxc and θyc.

固定座標変換器/位相変換器45では、積分器43の出力θxcr’、θycr’を、アキシャル方向ロータ回転角度θzを用いて固定座標系表現に変換するが、このとき、軸受制御系10の位置制御閉ループ位相特性を考慮して、次式(11)、(12)に示すようにアキシャル方向ロータ回転速度に応じて所定の位相量βを進めた上で固定座標系表現に変換した信号を、ラジアル方向回転変位指令θxc、θycとして軸受制御系10に対して出力する。   The fixed coordinate converter / phase converter 45 converts the outputs θxcr ′ and θycr ′ of the integrator 43 into a fixed coordinate system expression using the axial direction rotor rotation angle θz. At this time, the position of the bearing control system 10 is changed. Considering the control closed-loop phase characteristics, a signal converted into a fixed coordinate system expression after a predetermined phase amount β is advanced according to the axial direction rotor rotational speed as shown in the following equations (11) and (12), Output to the bearing control system 10 as radial direction rotational displacement commands θxc, θyc.

Figure 0004344601
Figure 0004344601

このラジアル方向回転変位指令θxc、θycの回転座標系表現θxcr、θycrは、図11(c)において回転ベクトル62で表され、ロータ慣性主軸Pは回転ベクトル62に対して、軸受制御系10の位置制御閉ループ位相特性の遅れに対応した位相角だけ、ロータ回転方向63と逆方向に回転した回転ベクトル64に従って回転座標系上を回転するため、定常状態では、図11(d)に示したように、ロータ慣性主軸Pが磁気軸受回転軸に対して平行となる。これにより、ロータ動不釣合いに起因する振れ回りモーメント60及びラジアル方向軸受モーメント指令の回転座標系表現61は共に零値に収束し、従ってロータに作用するラジアル方向軸受モーメント及びステータに作用するラジアル方向軸受モーメント反作用が共に低減した結果、ロータ動不釣合いに起因する擾乱モーメントのレベルが低下する。   The rotational coordinate system representations θxcr and θycr of the radial direction rotational displacement commands θxc and θyc are represented by the rotation vector 62 in FIG. 11C, and the rotor inertia main axis P is the position of the bearing control system 10 with respect to the rotation vector 62. As shown in FIG. 11 (d), in the steady state, the rotation is performed on the rotating coordinate system according to the rotation vector 64 rotated in the direction opposite to the rotor rotation direction 63 by the phase angle corresponding to the delay of the control closed loop phase characteristic. The rotor inertia main shaft P is parallel to the magnetic bearing rotating shaft. Thus, the rotational coordinate system representation 61 of the run-out moment 60 and the radial bearing moment command due to the rotor dynamic imbalance converges to zero values, and thus the radial bearing moment acting on the rotor and the radial direction acting on the stator. As a result of both bearing moment reactions being reduced, the level of disturbance moments due to rotor dynamic imbalance is reduced.

ここで、固定座標変換器/位相変換器45における進相量βを図11(e)に示したように過大値に設定した場合、ロータ慣性主軸Pの回転方向ベクトル64は、ラジアル方向軸受モーメント指令ベクトル61に対してロータ回転方向63に所定位相角だけ回転した配置となる。その結果、所定の時間が経過した後のロータ慣性主軸P上の固定点aは、図11(f)に示す配置となり、ロータ動不釣合いに起因する振れ回りモーメント及びラジアル方向軸受モーメント指令の回転座標系表現が各々60’、および61’で表されるモーメントとなる結果、ラジアル方向軸受モーメント指令の回転座標系表現61はロータ回転方向63と逆方向に回転する。   Here, when the phase advance amount β in the fixed coordinate converter / phase converter 45 is set to an excessive value as shown in FIG. 11 (e), the rotational direction vector 64 of the rotor inertia spindle P is represented by a radial bearing moment. The arrangement is rotated by a predetermined phase angle in the rotor rotation direction 63 with respect to the command vector 61. As a result, the fixed points a on the rotor inertia main shaft P after a predetermined time have elapsed are arranged as shown in FIG. 11 (f), and the rotation moment and radial direction bearing moment command due to rotor dynamic imbalance are rotated. As a result of the coordinate system representations being moments represented by 60 ′ and 61 ′, the rotational coordinate system representation 61 of the radial bearing moment command rotates in the direction opposite to the rotor rotational direction 63.

すなわち、ラジアル方向軸受モーメント指令の回転座標系表現61は、固定座標変換器/位相変換器45における進相量βが最適値に対して小さい場合(図11(a)、(b))はロータ回転方向63と同一方向に回転し、最適値に対して大きい場合(図11(e)、(f))はロータ回転方向63と逆方向に回転し、最適値の場合(図11(c)、(d))は方向は変化せずに零値に収束する。従って、どのようなアキシャル方向ロータ回転速度においても、固定座標変換器/位相変換器45における進相量βをラジアル方向軸受モーメント指令の回転座標系表現61が回転せずに零値に収束するような値に設定することにより、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ慣性主軸Pが磁気軸受中心軸と平行になることで、ロータ動不釣合いに起因する振れ回りモーメント60及びラジアル方向軸受モーメント指令の回転座標系表現61を共に零値に収束させることがでる。これにより、ロータに作用するラジアル方向軸受モーメント及びステータに作用するラジアル方向軸受モーメント反作用が共に低減するため、ロータ動不釣合いに起因する擾乱モーメントのレベルがサーボ共振周波数域を含めた全周波数域において低下するようになる。   In other words, the rotational coordinate system representation 61 of the radial bearing moment command indicates that the rotor is in the case where the phase advance amount β in the fixed coordinate converter / phase converter 45 is smaller than the optimum value (FIGS. 11A and 11B). When rotating in the same direction as the rotation direction 63 and larger than the optimum value (FIGS. 11E and 11F), the rotation is in the direction opposite to the rotor rotation direction 63, and the optimum value is obtained (FIG. 11C). , (D)) converges to a zero value without changing the direction. Accordingly, the phase advance amount β in the fixed coordinate converter / phase converter 45 is converged to a zero value without rotating the rotational coordinate system representation 61 of the radial bearing moment command at any axial rotor rotational speed. By setting the value to a proper value, the rotor inertia main shaft P is located at the center of the magnetic bearing regardless of the axial direction rotor rotational speed in any position within the position control band of the bearing control system 10, the servo resonance frequency area, and the position control band. By being parallel to the shaft, it is possible to converge both the run-out moment 60 caused by rotor dynamic imbalance and the rotational coordinate system representation 61 of the radial bearing moment command to zero values. As a result, the radial bearing moment acting on the rotor and the radial bearing moment reaction acting on the stator are both reduced, so that the level of disturbance moment caused by rotor dynamic imbalance is in the entire frequency range including the servo resonance frequency range. It begins to decline.

このように、この実施の形態3の低擾乱化制御装置により、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ慣性主軸が磁気軸受の回転軸と平行になり、ロータ動不釣合いに起因する擾乱モーメントを、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減することができる。   As described above, by the low disturbance control device of the third embodiment, the axial direction rotor rotational speed is in any frequency region within the position control band, the servo resonance frequency range, and the position control band of the bearing control system 10. The rotor inertia main shaft is parallel to the rotating shaft of the magnetic bearing, and the disturbance moment due to rotor dynamic imbalance can be reduced in the entire frequency region including the servo resonance frequency region of the magnetic bearing.

また、回転変位指令制御系40におけるLPF42がロータ回転同期成分のみを抽出することで、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, since the LPF 42 in the rotational displacement command control system 40 extracts only the rotor rotation synchronization component, signal components belonging to a frequency region other than the rotor rotation synchronization component can be easily removed.

更に、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御とし、位置制御ループゲイン特性が平坦となるように軸受制御系10を設計することにより、磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱の更なる抑制を図ることができる。   Furthermore, the bearing control system 10 is configured as a double loop control of a position control loop and a speed control loop, and the bearing control system 10 is designed so that the position control loop gain characteristic becomes flat, so that in the servo resonance frequency region of the magnetic bearing. It is possible to further suppress the amount of rotor swing and disturbance of the magnetic bearing device.

加えて、軸受制御系10において速度制御帯域51がアキシャル方向ロータ回転速度fの近傍に属することが無いように、ロータ回転速度fに応じて速度制御ループゲインKを変化させることにより、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域における磁気軸受装置の発生擾乱の更なる抑制を図ることができる。 In addition, as the speed control bandwidth 51 in the bearing control system 10 that is not belonging to the vicinity of the axial direction the rotor rotational speed f, by varying the speed control loop gain K D according to the rotor rotational speed f, bearing control Without impairing the stability of the system 10, it is possible to further suppress the disturbance of the magnetic bearing device in the speed control system gain intersection frequency range.

また、この実施の形態3では、速度制御ループゲインKの切換閾値が単一周波数fである場合について説明したが、この切換閾値はfを中心として所定のヒステリシス幅を持たせてもよい。この場合、上記実施の形態3と同様の効果があることに加えて、更にロータ回転速度信号fにノイズが重畳している場合、またはロータ回転速度fが速度制御ループゲイン切換閾値fの近傍で増減する場合において、より安定して速度制御ループゲインKの切換を実現することができる。 Further, in the third embodiment, although the switching threshold of the speed control loop gain K D of been described is a single frequency f 1, be the switching threshold is to have a predetermined hysteresis width around the f 1 Good. In this case, in addition to the same effects as in the third embodiment, noise is further superimposed on the rotor rotational speed signal f, or the rotor rotational speed f is close to the speed control loop gain switching threshold f 1 . in the case of increasing or decreasing, it is possible to realize the switching of a more stable speed control loop gain K D.

実施の形態4.
図12は、この発明に係る低擾乱化制御装置の実施の形態4を説明するためのブロック図であり、図3及び図5と同一の符号を付したものは、同一部分または相当部分を示す。
Embodiment 4 FIG.
FIG. 12 is a block diagram for explaining the fourth embodiment of the low disturbance control device according to the present invention, where the same reference numerals as those in FIGS. 3 and 5 designate the same or corresponding parts. .

次に動作について説明する。図12では、簡単のためにアキシャル方向並進運動及びラジアル方向回転運動のフィードバック・ループは図示していないが、軸受制御系10の動作はロータ変位演算器12を除いて実施の形態3と同一である。また、並進変位指令制御系20の動作は実施の形態2と同一である。   Next, the operation will be described. In FIG. 12, for the sake of simplicity, the feedback loop of the axial translational motion and the radial rotational motion is not shown, but the operation of the bearing control system 10 is the same as that of the third embodiment except for the rotor displacement calculator 12. is there. The operation of the translational displacement command control system 20 is the same as that in the second embodiment.

この実施の形態4において、ロータ変位演算器12では、電磁石/ロータ/変位検出手段11から出力されるロータ変位情報Si(i=1、2、…、n)に対して、例えばロータ1が幾何中心に位置する場合のロータ変位情報Si0(i = 1、2、…、n)及び変換係数γを用いて、次式(13)に示すようなロータが幾何中心に位置する場合を基準とした変位検出手段3とロータ1とのギャップ長変動量ΔXi(i=1、2、…、n)に変換する。   In the fourth embodiment, the rotor displacement calculator 12 uses, for example, the geometry of the rotor 1 for the rotor displacement information Si (i = 1, 2,..., N) output from the electromagnet / rotor / displacement detecting means 11. Using the rotor displacement information Si0 (i = 1, 2,..., N) and the conversion coefficient γ when located at the center, the case where the rotor as shown in the following equation (13) is located at the geometric center is used as a reference. It is converted into a gap length variation ΔXi (i = 1, 2,..., N) between the displacement detection means 3 and the rotor 1.

Figure 0004344601
Figure 0004344601

次に、上記ギャップ長変動量ΔXi(i=1、2、…、n)に対して、ロータ幾何形状及び変位検出手段3の配置を考慮して、次式(14)に示すようにロータの3自由度並進変位X、Y、Z及びラジアル方向2自由度回転変位θx、θyからなる5自由度ロータ変位を算出する。   Next, with respect to the gap length variation ΔXi (i = 1, 2,..., N), the rotor geometry and the arrangement of the displacement detection means 3 are taken into consideration, as shown in the following equation (14). A 5-degree-of-freedom rotor displacement consisting of 3-degree-of-freedom translational displacements X, Y, Z and radial direction 2-degree-of-freedom rotational displacements θx, θy is calculated.

Figure 0004344601
Figure 0004344601

上記式(14)において、Aは5行n列の変換行列である。
変位検出手段3がターゲットとするロータ1の変位検出面が真円であり、かつ変位検出手段3がロータ回転軸を中心とする同心円上に配置されているという前提のもとでは、ロータ1が幾何中心に位置する場合のロータ変位情報Si0(i=1、2、…、n)は、アキシャル方向ロータ回転角度θzに依存せずに全て同一の値(DC)となる。通常、ロータが幾何中心に位置する場合のロータ変位情報Si0(i=1、2、…、n)は、上記前提が成り立つと仮定してアキシャル方向ロータ回転角度θzと独立に全て同一のDC値を用いる場合が一般的であるが、この場合、ロータ変位検出面が真円に対して歪をもっている場合、または変位検出手段3の配置がロータ回転軸を中心とする同心円上にない場合において、ロータ変位演算器12の出力と、真のロータ変位との誤差が拡大する傾向を示す。このとき、軸受制御系10では誤差を含んだロータ5自由度変位に基づいて軸受力指令、ラジアル方向軸受モーメント指令及び軸受電磁石ギャップ長を算出し、この軸受力指令、ラジアル方向軸受モーメント指令及び軸受電磁石ギャップ長に基づいて電力分配器14において軸受電磁石に対して供給電流Ii(i=1、2、…、n)を分配して供給する。
In the above equation (14), A is a 5 × n conversion matrix.
Under the premise that the displacement detection surface of the rotor 1 targeted by the displacement detection means 3 is a perfect circle and the displacement detection means 3 is arranged on a concentric circle centered on the rotor rotation axis, the rotor 1 The rotor displacement information Si0 (i = 1, 2,..., N) when located at the geometric center all have the same value (DC) regardless of the axial direction rotor rotation angle θz. Normally, the rotor displacement information Si0 (i = 1, 2,..., N) when the rotor is located at the geometric center is assumed to be satisfied by the above assumption, and is all the same DC value independently of the axial direction rotor rotation angle θz. In this case, in the case where the rotor displacement detection surface is distorted with respect to a perfect circle, or the arrangement of the displacement detection means 3 is not on a concentric circle centered on the rotor rotation axis, An error between the output of the rotor displacement calculator 12 and the true rotor displacement tends to increase. At this time, the bearing control system 10 calculates a bearing force command, a radial direction bearing moment command and a bearing electromagnet gap length based on the rotor 5 degree-of-freedom displacement including an error, and this bearing force command, radial direction bearing moment command and bearing Based on the electromagnet gap length, the power distributor 14 distributes and supplies the supply current Ii (i = 1, 2,..., N) to the bearing electromagnet.

従って、ロータ変位演算器12の出力が誤差を含んだ結果、軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が拡大し、並進変位指令制御系20及び回転変位指令制御系40においてラジアル方向軸受力指令及びラジアル方向軸受モーメント指令が零ベクトルとなるように制御しても、静不釣合いに起因する擾乱力及び動不釣合いに起因する擾乱モーメントは所定の大きさで磁気軸受装置から発生し続けることとなる。   Accordingly, as a result of the output of the rotor displacement calculator 12 including an error, the error between the bearing force command and the radial bearing moment command and the true bearing force and the true radial bearing moment acting on the rotor is increased, and the translation is performed. In the displacement command control system 20 and the rotational displacement command control system 40, even if the radial direction bearing force command and the radial direction bearing moment command are controlled to be zero vectors, it is caused by disturbance force and dynamic imbalance caused by static unbalance. The disturbance moment to be generated continues to be generated from the magnetic bearing device with a predetermined magnitude.

そこで、この実施の形態4の低擾乱化制御装置では、ロータ変位演算器12の出力と、真のロータ変位との誤差が低減するように、ロータ変位演算器12に対してアキシャル方向ロータ回転角度θzに基づいた補正関数を導入する。上述のように、ロータ変位演算器12の出力が誤差を含む要因としては、ロータ変位検出面の真円に対する歪及び変位検出手段3配置のロータ回転軸を中心とする同心円上からのずれが挙げられ、ロータ1が幾何中心に位置する場合のロータ変位情報Si0(i=1、2、…、n)に対して、前者はアキシャル方向ロータ回転角度θzに応じて変動するAC成分誤差、後者は変位検出手段3の各々で個別にばらつくDC成分誤差として作用する。   Therefore, in the low disturbance control device of the fourth embodiment, the axial rotor rotation angle with respect to the rotor displacement calculator 12 is reduced so that an error between the output of the rotor displacement calculator 12 and the true rotor displacement is reduced. A correction function based on θz is introduced. As described above, the factors that include errors in the output of the rotor displacement calculator 12 include the distortion of the rotor displacement detection surface with respect to the true circle and the deviation of the displacement detection means 3 from the concentric circle centered on the rotor rotation axis. For the rotor displacement information Si0 (i = 1, 2,..., N) when the rotor 1 is located at the geometric center, the former is an AC component error that varies according to the axial direction rotor rotation angle θz, and the latter is Each of the displacement detection means 3 acts as a DC component error that varies individually.

従って、幾何中心におけるロータ変位情報Si0(i=1、2、…、n)に対して、上述のAC成分誤差及びDC成分誤差を除去するようにアキシャル方向ロータ回転角度θzに基づいた補正関数を導入し、次式(15)に示すようにギャップ長変動量ΔXi(i=1、2、…、n)を算出する。   Therefore, a correction function based on the axial direction rotor rotation angle θz is applied to the rotor displacement information Si0 (i = 1, 2,..., N) at the geometric center so as to remove the AC component error and the DC component error. Then, the gap length variation ΔXi (i = 1, 2,..., N) is calculated as shown in the following equation (15).

Figure 0004344601
Figure 0004344601

上記式(15)におけるSi0(θz)(i=1、2、…、n)は、ロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出することができる。図12に示したブロック図において、並進変位指令制御系20を停止させ、軸受制御系10に対してラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycとして全て零値を入力した場合、ロータが回転せずに静止浮上しているならば、PI制御器16における積分器の作用により指令に対する定常偏差が零となり、ロータ変位演算器12の出力である5自由度ロータ変位X、Y、Z、θx、θyは全て零値に保持される。このとき、ロータ重心の真の3自由度並進変位をXt、Yt、Zt、ロータ慣性主軸Pの真のラジアル方向2自由度回転変位をθxt、θyt、ロータ質量をM、重力加速度をg、磁気軸受によりロータに作用する真の3自由度軸受力をFxt、Fyt、Fzt、真のラジアル方向2自由度軸受モーメントをNxt、Nytとすると、図2に示した固定座標系において次式(16)、(17)の静力学が成り立つ。

Figure 0004344601
Si0 (θz) (i = 1, 2,..., N) in the above formula (15) is the statics when the rotor is statically levitated without rotating around the rotation axis, and the bearing electromagnet supply current measurement at this time It can be derived based on the value. In the block diagram shown in FIG. 12, the translational displacement command control system 20 is stopped, and the radial translational displacement commands Xc and Yc, the axial translational displacement command Zc, and the radial rotational displacement commands θxc and θyc with respect to the bearing control system 10. When all the zero values are input, the steady deviation with respect to the command becomes zero due to the action of the integrator in the PI controller 16 if the rotor is still floating without rotating, which is the output of the rotor displacement calculator 12. The 5-degree-of-freedom rotor displacements X, Y, Z, θx, and θy are all held at zero values. At this time, true three-degree-of-freedom translational displacement of the rotor center of gravity is Xt, Yt, Zt, true radial two-degree-of-freedom rotational displacement of the rotor inertial spindle P is θxt, θyt, rotor mass is M, gravity acceleration is g, magnetic When the true three-degree-of-freedom bearing force acting on the rotor by the bearing is Fxt, Fyt, Fzt, and the true radial two-degree-of-freedom bearing moment is Nxt, Nyt, in the fixed coordinate system shown in FIG. , (17) holds.
Figure 0004344601

このとき、軸受電磁石で発生している真の磁気吸引力をFit(i=1、2、…、n)とすると、真の軸受力Fxt、Fyt、Fzt、真のラジアル方向軸受モーメントNxt、Nytと真の磁気吸引力Fit(i=1、2、…、n)との関係は、Bを5行n列の変換行列として次式(18)のように表される。   At this time, if the true magnetic attractive force generated by the bearing electromagnet is Fit (i = 1, 2,..., N), the true bearing forces Fxt, Fyt, Fzt, the true radial bearing moments Nxt, Nyt And the true magnetic attraction force Fit (i = 1, 2,..., N) is expressed by the following equation (18), where B is a 5 × n conversion matrix.

Figure 0004344601
Figure 0004344601

また、軸受電磁石における真の磁気吸引力Fit(i=1、2、…、n)は、各軸受電磁石に対する供給電流計測値をIi(i=1、2、…、n)、各軸受電磁石とロータ吸引面との真のギャップ長をGit(i=1、2、…、n)、変換係数をkとして次式(19)のように表される。   In addition, the true magnetic attractive force Fit (i = 1, 2,..., N) in the bearing electromagnet is the measured supply current value for each bearing electromagnet Ii (i = 1, 2,..., N), The true gap length with the rotor suction surface is expressed by the following equation (19), where Git (i = 1, 2,..., N) and the conversion coefficient are k.

Figure 0004344601
Figure 0004344601

更に、真の軸受電磁石ギャップ長Git(i=1、2、…、n)は、ロータが幾何中心に位置する場合の軸受電磁石ギャップ長をGi0(i=1、2、…、n)、Cをn行5列の変換行列として次式(20)で与えられる。   Further, the true bearing electromagnet gap length Git (i = 1, 2,..., N) is equal to Gi0 (i = 1, 2,..., N), C when the rotor is located at the geometric center. Is given by the following equation (20) as an n × 5 conversion matrix.

Figure 0004344601
Figure 0004344601

上記式(18)〜(20)を上記式(16)、(17)に代入し、所定のアキシャル方向ロータ回転角度θzで静止浮上させた場合の軸受電磁石供給電流計測値Ii(i=1、2、…、n)を適用することで、当該回転角度θzにおける真の5自由度ロータ変位Xt、Yt、Zt、θxt、θytに関する非線形連立方程式を得る。この非線形連立方程式は、Newton法等の計算手法により数値解を導出することができ、更にこのとき、ロータ変位演算器12の出力である5自由度ロータ変位X、Y、Z、θx、θyは全て零値に保持されているので、この数値解(真の5自由度ロータ変位Xt、Yt、Zt、θxt、θyt)はロータ変位演算器12における出力誤差と一致する。   The above formulas (18) to (20) are substituted into the above formulas (16) and (17), and the bearing electromagnet supply current measurement value Ii (i = 1, i.e., when statically levitated at a predetermined axial direction rotor rotation angle θz. 2,..., N) is applied to obtain a simultaneous simultaneous equation relating to true 5-degree-of-freedom rotor displacement Xt, Yt, Zt, θxt, θyt at the rotation angle θz. The nonlinear simultaneous equations can be derived numerically by a calculation method such as Newton's method, and at this time, the 5-degree-of-freedom rotor displacement X, Y, Z, θx, θy, which is the output of the rotor displacement calculator 12, is Since all of them are held at zero values, this numerical solution (true five-degree-of-freedom rotor displacement Xt, Yt, Zt, θxt, θyt) matches the output error in the rotor displacement calculator 12.

上記手順により、所定の回転角度θzにおけるロータ変位演算器12の出力誤差、従って上記式(15)におけるSi0(θz)(i=1、2、…、n)を導出することができる。そして、上記式(15)によって得られたギャップ長変動量ΔXi(i=1、2、…、n)を上記式(14)に適用することにより、ロータ変位演算器12の出力と、真のロータ変位との誤差を低減することができ、従って軸受制御系10における軸受力指令、およびラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力、および真のラジアル方向軸受モーメントとの誤差が減少し、磁気軸受装置の発生擾乱をより一層低減することができる。   By the above procedure, it is possible to derive the output error of the rotor displacement calculator 12 at a predetermined rotation angle θz, and thus Si0 (θz) (i = 1, 2,..., N) in the above equation (15). Then, by applying the gap length variation ΔXi (i = 1, 2,..., N) obtained by the above equation (15) to the above equation (14), the output of the rotor displacement calculator 12 and the true value The error with the rotor displacement can be reduced, and therefore the error between the bearing force command and the radial bearing moment command in the bearing control system 10 and the true bearing force and the true radial bearing moment acting on the rotor can be reduced. This can reduce the disturbance of the magnetic bearing device.

このように、この発明の実施の形態4の低擾乱化制御装置により、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ重心位置Gが磁気軸受の回転軸上に位置するようになり、ロータ静不釣合いに起因する擾乱力を、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減することができる。   Thus, with the low disturbance control device according to the fourth embodiment of the present invention, the axial direction rotor rotational speed is in any frequency region within the position control band of the bearing control system 10, the servo resonance frequency region, and outside the position control band. Even in this case, the rotor center-of-gravity position G is positioned on the rotating shaft of the magnetic bearing, and the disturbance force caused by the rotor static imbalance can be reduced in the entire frequency region including the servo resonance frequency region of the magnetic bearing. it can.

また、並進変位指令制御系20におけるLPF22がロータ回転同期成分のみを抽出することで、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, since the LPF 22 in the translational displacement command control system 20 extracts only the rotor rotation synchronization component, signal components belonging to a frequency region other than the rotor rotation synchronization component can be easily removed.

更に、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御とし、位置制御ループゲイン特性が平坦となるように軸受制御系10を設計することにより、磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱の更なる抑制を図ることができるばかりか、軸受制御系10において速度制御帯域51がアキシャル方向ロータ回転速度fの近傍に属することが無いように、ロータ回転速度fに応じて速度制御ループゲインKを変化させることにより、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域における磁気軸受装置の発生擾乱の更なる抑制を図ることができる。 Furthermore, the bearing control system 10 is configured as a double loop control of a position control loop and a speed control loop, and the bearing control system 10 is designed so that the position control loop gain characteristic becomes flat, so that in the servo resonance frequency region of the magnetic bearing. In addition to being able to further suppress the amount of rotor swing and disturbance generated by the magnetic bearing device, the rotor is controlled so that the speed control band 51 does not belong to the vicinity of the axial rotor rotational speed f in the bearing control system 10. by varying the speed control loop gain K D according to the rotational speed f, without impairing the stability of the bearing control system 10, achieve a further reduction of generation disturbance of the magnetic bearing apparatus in the speed control system gain crossover frequency range be able to.

加えて、ロータ変位演算器12の出力と、真のロータ変位との誤差を低減することにより、軸受制御系10における軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が低減し、磁気軸受装置の発生擾乱をより一層低減することができる。   In addition, by reducing an error between the output of the rotor displacement calculator 12 and the true rotor displacement, the bearing force command and the radial bearing moment command in the bearing control system 10, the true bearing force acting on the rotor, and The error from the true radial bearing moment is reduced, and the disturbance generated in the magnetic bearing device can be further reduced.

実施の形態5.
図13は、この発明に係る低擾乱化制御装置の実施の形態5を説明するためのブロック図、図14は、補正関数を導出するためのデータ計測における変位検出手段及びロータ変位検出面の配置を示した図であり、図13において、図3及び図5と同一の符号を付したものは、同一部分または相当部分を示す。
Embodiment 5 FIG.
FIG. 13 is a block diagram for explaining the fifth embodiment of the low disturbance control device according to the present invention, and FIG. 14 is an arrangement of displacement detecting means and rotor displacement detecting surfaces in data measurement for deriving a correction function. In FIG. 13, the same reference numerals as those in FIGS. 3 and 5 designate the same or corresponding parts.

図14に示したように、この実施の形態5の低擾乱化制御装置は、ロータ変位検出方向71を検出する任意の変位検出手段70、変位検出手段70により、ロータ変位検出方向71のロータ変位が計測されるロータ変位検出面72、ロータ回転軸(Z軸)73を備える。   As shown in FIG. 14, the low disturbance control device according to the fifth embodiment has an arbitrary displacement detection means 70 for detecting the rotor displacement detection direction 71, and a rotor displacement in the rotor displacement detection direction 71 by the displacement detection means 70. Is provided with a rotor displacement detection surface 72 and a rotor rotation axis (Z-axis) 73.

次に動作について説明する。図13では、簡単のためにアキシャル方向並進運動及びラジアル方向回転運動のフィードバック・ループは図示していないが、軸受制御系10の動作はロータ変位演算器12及びギャップ長演算器13を除いて実施の形態3と同一である。また、並進変位指令制御系20の動作は実施の形態2と同一である。   Next, the operation will be described. In FIG. 13, for the sake of simplicity, the feedback loop of the axial translational motion and the radial rotational motion is not shown, but the operation of the bearing control system 10 is performed except for the rotor displacement calculator 12 and the gap length calculator 13. This is the same as Form 3. The operation of the translational displacement command control system 20 is the same as that in the second embodiment.

この実施の形態5においては、ロータ変位演算器12では、電磁石/ロータ/変位検出手段11から出力されるロータ変位情報Si(i=1、2、…、n)に対して、アキシャル方向ロータ回転角度θzに基づいた補正関数を導入したロータ幾何中心時のロータ変位情報Si0(i=1、2、…、n)及び変換係数γを用いて、上記式(15)により変位検出手段70とロータとのギャップ長変動量ΔXi(i=1、2、…、n)を算出し、次いでギャップ長変動量ΔXi(i=1、2、…、n)を用いて上記式(14)により5自由度ロータ変位X、Y、Z、θx、θy算出することにより、ロータ変位演算器12の出力と、真のロータ変位との誤差を低減することができる。   In the fifth embodiment, the rotor displacement calculator 12 rotates the rotor in the axial direction with respect to the rotor displacement information Si (i = 1, 2,..., N) output from the electromagnet / rotor / displacement detecting means 11. Using the rotor displacement information Si0 (i = 1, 2,..., N) and the conversion coefficient γ at the rotor geometric center with the correction function based on the angle θz and the conversion coefficient γ, the displacement detection means 70 and the rotor are expressed by the above equation (15). The gap length variation amount ΔXi (i = 1, 2,..., N) is calculated, and then the gap length variation amount ΔXi (i = 1, 2,..., N) is used to calculate 5 degrees according to the above equation (14). By calculating the degree of rotor displacement X, Y, Z, θx, θy, the error between the output of the rotor displacement calculator 12 and the true rotor displacement can be reduced.

また、ギャップ長演算器13では、ロータが幾何中心に位置する場合の軸受電磁石ギャップ長Gi0(i=1、2、…、n)、n行5列変換行列Cを用いて、ロータ変位演算器12の出力である5自由度ロータ変位X、Y、Z、θx、θyから次式(21)に示すように軸受電磁石とロータ吸引面とのギャップ長Gi(i=1、2、…、n)を算出する。   Further, the gap length calculator 13 uses a bearing electromagnet gap length Gi0 (i = 1, 2,..., N) and an n-row / 5-column conversion matrix C when the rotor is located at the geometric center, to calculate the rotor displacement calculator. From the 5-degree-of-freedom rotor displacements X, Y, Z, θx, and θy that are 12 outputs, the gap length Gi (i = 1, 2,..., N) between the bearing electromagnet and the rotor attraction surface as shown in the following equation (21): ) Is calculated.

Figure 0004344601
Figure 0004344601

ここで、軸受電磁石がターゲットとするロータ吸引面が真円であり、かつ軸受電磁石がロータ回転軸を中心とする同心円上に配置されているという前提のもとでは、ロータが幾何中心に位置する場合の軸受電磁石ギャップ長Gi0(i=1、2、…、n)は、アキシャル方向ロータ回転角度θzに依存せずに全て同一のDC値となる。通常、Gi0(i=1、2、…、n)は、上記前提が成り立つと仮定してアキシャル方向ロータ回転角度θzと独立に全て同一のDC値を用いる場合が一般的であるが、ロータ吸引面が真円に対して歪をもっている場合、または、軸受電磁石の配置がロータ回転軸を中心とする同心円上にない場合において、ギャップ長演算器13の出力と、真の軸受電磁石ギャップ長との誤差が拡大する傾向を示す。このとき、軸受制御系10における電力分配器14では誤差を含んだ軸受電磁石ギャップ長に基づいて、軸受力指令、およびラジアル方向軸受モーメント指令を実現するための供給電流Ii(i=1、2、…、n)を軸受電磁石に対して分配して供給する。   Here, the rotor is positioned at the geometric center under the assumption that the rotor attracting surface targeted by the bearing electromagnet is a perfect circle and the bearing electromagnet is arranged on a concentric circle centered on the rotor rotation axis. The bearing electromagnet gap lengths Gi0 (i = 1, 2,..., N) in this case all have the same DC value without depending on the axial rotor rotation angle θz. In general, Gi0 (i = 1, 2,..., N) is generally assumed to use the same DC value independently of the axial-direction rotor rotation angle θz on the assumption that the above assumption is satisfied. When the surface is distorted with respect to the perfect circle, or when the arrangement of the bearing electromagnet is not on a concentric circle centered on the rotor rotation axis, the output of the gap length calculator 13 and the true bearing electromagnet gap length Shows a tendency for errors to expand. At this time, the power distributor 14 in the bearing control system 10 supplies a supply current Ii (i = 1, 2,) for realizing a bearing force command and a radial bearing moment command based on the bearing electromagnet gap length including an error. ..., n) are distributed and supplied to the bearing electromagnets.

従って、ギャップ長演算器13の出力が誤差を含んだ結果、軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が拡大し、並進変位指令制御系20及び回転変位指令制御系40においてラジアル方向軸受力指令及びラジアル方向軸受モーメント指令が零ベクトルとなるように制御しても、静不釣合いに起因する擾乱力及び動不釣合いに起因する擾乱モーメントは所定の大きさで磁気軸受装置から発生し続けることとなる。   Therefore, as a result of the output of the gap length calculator 13 including an error, errors between the bearing force command and the radial bearing moment command and the true bearing force and the true radial bearing moment acting on the rotor are expanded, and the translation is performed. In the displacement command control system 20 and the rotational displacement command control system 40, even if the radial direction bearing force command and the radial direction bearing moment command are controlled to be zero vectors, it is caused by disturbance force and dynamic imbalance caused by static unbalance. The disturbance moment to be generated continues to be generated from the magnetic bearing device with a predetermined magnitude.

そこで、この実施の形態5の低擾乱化制御装置では、ギャップ長演算器13の出力と、真の軸受電磁石ギャップ長との誤差が低減するように、ギャップ長演算器13に対してアキシャル方向ロータ回転角度θzに基づいた補正関数を導入する。上述のように、ギャップ長演算器13の出力が誤差を含む要因としては、ロータ吸引面の真円に対する歪及び軸受電磁石配置のロータ回転軸を中心とする同心円上からのずれが挙げられ、ロータ幾何中心時の軸受電磁石ギャップ長Gi0(i=1、2、…、n)に対して、前者はアキシャル方向ロータ回転角度θzに応じて変動するAC成分誤差、後者は軸受電磁石の各々で個別にばらつくDC成分誤差として作用する。   Therefore, in the low disturbance control device of the fifth embodiment, the axial direction rotor with respect to the gap length calculator 13 is reduced so that the error between the output of the gap length calculator 13 and the true bearing electromagnet gap length is reduced. A correction function based on the rotation angle θz is introduced. As described above, the factors that include an error in the output of the gap length calculator 13 include distortion of the rotor attracting surface with respect to a perfect circle and deviation of the bearing electromagnet arrangement from a concentric circle around the rotor rotation axis. With respect to the bearing electromagnet gap length Gi0 (i = 1, 2,..., N) at the geometric center, the former is an AC component error that varies according to the axial rotor rotation angle θz, and the latter is individually determined for each bearing electromagnet. It acts as a DC component error that varies.

従って、幾何中心における軸受電磁石ギャップ長Gi0(i=1、2、…、n)に対して、上述のAC成分誤差及びDC成分誤差を除去するようにアキシャル方向ロータ回転角度θzに基づいた補正関数を導入し、次式(22)に示すように軸受電磁石ギャップ長Gi(i=1、2、…、n)を算出する。   Therefore, the correction function based on the axial direction rotor rotation angle θz so as to remove the above-described AC component error and DC component error with respect to the bearing electromagnet gap length Gi0 (i = 1, 2,..., N) at the geometric center. And the bearing electromagnet gap length Gi (i = 1, 2,..., N) is calculated as shown in the following equation (22).

Figure 0004344601
Figure 0004344601

このように、この実施の形態5の低擾乱化制御装置では、アキシャル方向ロータ回転角度θzに基づいた補正関数を、ロータ変位演算器12における上記式(15)及びギャップ長演算器13における上記式(22)に導入しており、上記式(15)におけるSi0(θz)(i=1、2、…、n)は磁気軸受の保護を目的としたタッチダウン軸受によってロータを機械的に拘束した場合のロータ変位情報に基づいて、上記式(22)におけるGi0(θz)(i=1、2、…、n)はロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出することができる。   As described above, in the low disturbance control device of the fifth embodiment, the correction function based on the axial direction rotor rotation angle θz is expressed by the above equation (15) in the rotor displacement calculator 12 and the above equation in the gap length calculator 13. (22), and Si0 (θz) (i = 1, 2,..., N) in the above formula (15) mechanically constrains the rotor by a touchdown bearing for the purpose of protecting the magnetic bearing. Gi0 (θz) (i = 1, 2,..., N) in the above equation (22) is based on the rotor displacement information in the case, and the statics when the rotor is statically levitated without rotating around the rotation axis. It can derive | lead-out based on the bearing electromagnet supply current measured value at this time.

一般に、磁気軸受装置には、軸受電磁石及び変位検出手段とロータとの接触を回避して磁気軸受装置を保護するために、ロータがその幾何学中心から所定量だけ変位した場合に、ロータを機械的に拘束できるようなタッチダウン軸受が装備されている。そこで、軸受制御系10を動作させずにタッチダウン軸受でロータを機械的に拘束し、特に変位検出手段のロータ変位検出方向にロータ1をタッチダウンさせた場合のロータ変位情報を用いることで、ロータ幾何中心時のロータ変位情報Si0(i=1、2、…、n)を得ることができる。   In general, in a magnetic bearing device, in order to protect the magnetic bearing device by avoiding contact between the bearing electromagnet and the displacement detecting means and the rotor, the rotor is mechanically moved when the rotor is displaced from the geometric center by a predetermined amount. It is equipped with touch-down bearings that can be restrained. Therefore, by using the rotor displacement information when the rotor 1 is mechanically restrained by the touchdown bearing without operating the bearing control system 10, and the rotor 1 is touched down in the rotor displacement detection direction of the displacement detection means, in particular, Rotor displacement information Si0 (i = 1, 2,..., N) at the rotor geometric center can be obtained.

例えば、円柱形状ロータ(図1)のラジアル方向並進変位を計測する変位検出手段と、ロータ変位検出面の関係は、図14(a)のように示すことができる。ロータ回転軸73回りに回転する円柱形状ロータでは、ロータのラジアル方向並進変位を計測する変位検出手段70に対して、ロータ側面がロータ変位検出面72となる。ロータが幾何中心に位置する場合のロータ変位検出面72に対して、ロータ変位検出方向71のうち、変位検出手段70とのギャップ長が小さくなる向きに、ロータをタッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’、ロータ変位検出方向71のうち、変位検出手段70とのギャップ長が大きくなる向きに、ロータをタッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’’である。通常、ロータ変位検出方向71に対するロータ変位検出面72とロータ変位検出面72’との距離及びロータ変位検出面72とロータ変位検出面72’’との距離は同一であるので、ロータ変位検出面72’におけるロータ変位情報及びロータ変位検出面72’’におけるロータ変位情報の平均により、ロータ幾何中心時(ロータ変位検出面72)のロータ変位情報Si0(i=1、2、…、n)を推定することができる。   For example, the relationship between the displacement detection means for measuring the radial translational displacement of the cylindrical rotor (FIG. 1) and the rotor displacement detection surface can be shown as shown in FIG. In the cylindrical rotor rotating around the rotor rotation shaft 73, the rotor side surface becomes the rotor displacement detection surface 72 with respect to the displacement detection means 70 for measuring the radial translational displacement of the rotor. The rotor is mechanically constrained by a touch-down bearing in a direction in which the gap length with the displacement detection means 70 becomes smaller in the rotor displacement detection direction 71 with respect to the rotor displacement detection surface 72 when the rotor is located at the geometric center. The rotor displacement detection surface when the rotor is mechanically constrained by a touch-down bearing in a direction in which the gap length with the displacement detection means 70 becomes larger in the rotor displacement detection direction 71 is 72 ′. Is 72 ″. Usually, the distance between the rotor displacement detection surface 72 and the rotor displacement detection surface 72 ′ and the distance between the rotor displacement detection surface 72 and the rotor displacement detection surface 72 ″ with respect to the rotor displacement detection direction 71 are the same. Based on the average of the rotor displacement information at 72 ′ and the rotor displacement information at the rotor displacement detection surface 72 ″, the rotor displacement information Si0 (i = 1, 2,..., N) at the rotor geometric center (rotor displacement detection surface 72) is obtained. Can be estimated.

また、ロータ変位検出面72、従ってロータ変位検出方向71がロータ回転軸73に対して傾斜するような形状のロータ(図2)において、変位検出手段とロータ変位検出面の関係は図14(b)の如く示すことができる。ロータが幾何中心に位置する場合のロータ変位検出面72に対して、ロータ配置をアキシャル方向(Z軸方向)には正方向、ラジアル方向には変位検出手段70とのギャップ長が小さくなる向きに、タッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’、ロータ配置をアキシャル方向(Z軸方向)には正方向、ラジアル方向には変位検出手段70とのギャップ長が大きくなる向きに、タッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’’、ロータ配置をアキシャル方向(Z軸方向)には負方向、ラジアル方向には変位検出手段70とのギャップ長が大きくなる向きに、タッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’’’、ロータ配置をアキシャル方向(Z軸方向)には負方向、ラジアル方向には変位検出手段70とのギャップ長が小さくなる向きに、タッチダウン軸受により機械的に拘束した場合のロータ変位検出面が72’’’’である。   Further, in the rotor having a shape in which the rotor displacement detection surface 72, and hence the rotor displacement detection direction 71 is inclined with respect to the rotor rotation shaft 73 (FIG. 2), the relationship between the displacement detection means and the rotor displacement detection surface is as shown in FIG. ). With respect to the rotor displacement detection surface 72 when the rotor is located at the geometric center, the rotor is arranged in the positive direction in the axial direction (Z-axis direction) and in the direction in which the gap length with the displacement detection means 70 becomes smaller in the radial direction. When the rotor is mechanically restrained by the touchdown bearing, the rotor displacement detection surface is 72 ', the rotor arrangement is positive in the axial direction (Z-axis direction), and the gap length with the displacement detection means 70 is large in the radial direction. When the rotor is mechanically constrained by a touchdown bearing, the rotor displacement detection surface is 72 ″, the rotor arrangement is negative in the axial direction (Z-axis direction), and the gap length with the displacement detection means 70 in the radial direction. When the rotor is mechanically constrained by a touchdown bearing, the rotor displacement detection surface is 72 '' 'and the rotor arrangement is in the axial direction (Z-axis direction). In the negative direction, facing the gap length between the displacement detection means 70 in the radial direction becomes smaller, the rotor displacement detection surface in the case of mechanically constrained by touchdown bearing is a 72 '' ''.

このとき、ロータ変位検出面72’〜72’’’’におけるロータ変位情報の平均により、ロータ幾何中心時(ロータ変位検出面72)のロータ変位情報Si0(i=1、2、…、n)を推定することができる。   At this time, the rotor displacement information Si0 (i = 1, 2,..., N) at the time of the rotor geometric center (rotor displacement detection surface 72) is obtained by averaging the rotor displacement information on the rotor displacement detection surfaces 72 ′ to 72 ″ ″. Can be estimated.

上記手順により、所定の回転角度θzにおけるロータ幾何中心時のロータ変位情報Si0(i=1、2、…、n)、従って上記式(15)におけるSi0(θz)(i=1、2、…、n)を導出することができる。   According to the above procedure, the rotor displacement information Si0 (i = 1, 2,..., N) at the rotor geometric center at the predetermined rotation angle θz, and thus Si0 (θz) (i = 1, 2,..., In the above equation (15). , N) can be derived.

次に、上記式(22)におけるGi0(θz)(i=1、2、…、n)は、ロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出することができる。図13に示すブロック図において、ロータ変位演算器12に対して上述のアキシャル方向ロータ回転角度θzに基づいた補正関数を導入し、かつ、並進変位指令制御系20を停止させた状態で軸受制御系10に対してラジアル方向並進変位指令Xc、Yc、アキシャル方向並進変位指令Zc、ラジアル方向回転変位指令θxc、θycとして全て零値を入力した場合、ロータが回転せずに静止浮上しているならば、PI制御器16における積分器の作用により指令に対する定常偏差が零となり、ロータ変位演算器12の出力である5自由度ロータ変位X、Y、Z、θx、θyは全て零値に保持される。   Next, Gi0 (θz) (i = 1, 2,..., N) in the above formula (22) is the statics when the rotor is floated without rotating around the rotation axis and the bearing electromagnet at this time. It can be derived based on the measured supply current. In the block diagram shown in FIG. 13, the bearing control system is introduced in a state where the correction function based on the axial rotor rotation angle θz is introduced to the rotor displacement calculator 12 and the translational displacement command control system 20 is stopped. If all zero values are input as radial translation commands Xc and Yc, axial translational displacement commands Zc, radial rotational displacement commands θxc and θyc with respect to 10, the rotor does not rotate and floats statically. The steady deviation with respect to the command becomes zero by the action of the integrator in the PI controller 16, and the 5-degree-of-freedom rotor displacements X, Y, Z, θx, and θy that are the outputs of the rotor displacement calculator 12 are all held at zero values. .

この5自由度ロータ変位は、ロータ変位演算器12においてアキシャル方向ロータ回転角度θzに基づいた補正関数を導入していることから、真の5自由度ロータ変位と見なすことができ、従って、ロータ質量をM、重力加速度をg、磁気軸受によりロータに作用する真の3自由度軸受力をFxt、Fyt、Fzt、真のラジアル方向2自由度軸受モーメントをNxt、Nytとすると、図1及び図2に示した固定座標系において次式(23)、(24)のような静力学が成り立つ。   This 5-degree-of-freedom rotor displacement can be regarded as a true 5-degree-of-freedom rotor displacement because the rotor displacement calculator 12 introduces a correction function based on the axial-direction rotor rotation angle θz. 1 and 2, where M is gravity acceleration, g is gravitational acceleration, Fxt, Fyt, Fzt are true three-degree-of-freedom bearing forces acting on the rotor by magnetic bearings, and Nxt, Nyt are true radial two-degree-of-freedom bearing moments. In the fixed coordinate system shown in (1), the following static equations (23) and (24) hold.

Figure 0004344601
Figure 0004344601

このとき、軸受電磁石で発生している真の磁気吸引力をFit(i=1、2、…、n)とすると、真の軸受力Fxt、Fyt、Fzt、真のラジアル方向軸受モーメントNxt、Nytと真の磁気吸引力Fit(i=1、2、…、n)との関係は、Bを5行n列の変換行列として上記式(18)のように表される。   At this time, if the true magnetic attractive force generated by the bearing electromagnet is Fit (i = 1, 2,..., N), the true bearing forces Fxt, Fyt, Fzt, the true radial bearing moments Nxt, Nyt And the true magnetic attraction force Fit (i = 1, 2,..., N) is expressed by the above equation (18), where B is a 5 × n conversion matrix.

また、軸受電磁石における真の磁気吸引力Fit(i=1、2、…、n)は、各軸受電磁石に対する供給電流計測値をIi(i=1、2、…、n)、各軸受電磁石とロータ吸引面との真のギャップ長をGit(i=1、2、…、n)、変換係数をkとして上記式(19)のように表される。   In addition, the true magnetic attractive force Fit (i = 1, 2,..., N) in the bearing electromagnet is the measured supply current value for each bearing electromagnet Ii (i = 1, 2,..., N), The true gap length with the rotor suction surface is represented by the above equation (19) with Git (i = 1, 2,..., N) and the conversion coefficient k.

更に、真の軸受電磁石ギャップ長Git(i=1、2、…、n)は、ロータが幾何中心に位置する場合の軸受電磁石ギャップ長をGi0(i=1、2、…、n)として、上記式(20)において真の5自由度ロータ変位が全て零値であるとして、次式(25)で与えられる。   Further, the true bearing electromagnet gap length Git (i = 1, 2,..., N) is defined as Gi0 (i = 1, 2,..., N) when the rotor is located at the geometric center. In the above equation (20), assuming that the true 5-degree-of-freedom rotor displacement is all zero, the following equation (25) is given.

Figure 0004344601
Figure 0004344601

上記式(18)、(19)、(25)を上記式(23)、(24)に代入し、所定のアキシャル方向ロータ回転角度θzで静止浮上させた場合の軸受電磁石供給電流計測値Ii(i=1、2、…、n)を適用することで、当該回転角度θzにおける真のロータ幾何中心時の軸受電磁石ギャップ長Gi0(i=1、2、…、n)に関する非線形連立方程式を得る。この非線形連立方程式は、Newton法等の計算手法により数値解を導出することができ、上記手順を所定の回転角度θzにおいて実施することにより、上記式(22)におけるGi0(θz)(i=1、2、…、n)を導出することができる。上記式(22)によってギャップ長演算器13を構成することにより、ギャップ長演算器13の出力と、真の軸受電磁石ギャップ長との誤差を低減することができ、従って、軸受制御系10における軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が低減し、磁気軸受装置の発生擾乱をより一層低減することができる。   The above formulas (18), (19), and (25) are substituted into the above formulas (23) and (24), and the measured value of the bearing electromagnet supply current Ii (when statically levitated at a predetermined axial direction rotor rotation angle θz) By applying i = 1, 2,..., n), a nonlinear simultaneous equation relating to the bearing electromagnet gap length Gi0 (i = 1, 2,..., n) at the true rotor geometric center at the rotation angle θz is obtained. . This nonlinear simultaneous equation can derive a numerical solution by a calculation method such as Newton's method, and by performing the above procedure at a predetermined rotation angle θz, Gi0 (θz) (i = 1) in the above equation (22). 2, ..., n) can be derived. By configuring the gap length calculator 13 by the above equation (22), an error between the output of the gap length calculator 13 and the true bearing electromagnet gap length can be reduced. Therefore, the bearing in the bearing control system 10 can be reduced. The error between the force command and the radial direction bearing moment command and the true bearing force and the true radial direction moment acting on the rotor is reduced, and the disturbance generated in the magnetic bearing device can be further reduced.

このように、この発明における低擾乱化制御装置の実施の形態5により、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ重心位置Gが磁気軸受の回転軸上に位置するようになり、磁気軸受のサーボ共振周波数域を含めた全周波数領域において、ロータ静不釣合いに起因する擾乱力を低減することができる。   Thus, according to the fifth embodiment of the low disturbance control device of the present invention, the axial direction rotor rotational speed is in any frequency region within the position control band of the bearing control system 10, the servo resonance frequency region, and outside the position control band. Even in this case, the rotor center-of-gravity position G is positioned on the rotating shaft of the magnetic bearing, and the disturbance force due to the rotor unbalance can be reduced in all frequency regions including the servo resonance frequency region of the magnetic bearing. it can.

また、並進変位指令制御系20におけるLPF22がロータ回転同期成分のみを抽出することで、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, since the LPF 22 in the translational displacement command control system 20 extracts only the rotor rotation synchronization component, signal components belonging to a frequency region other than the rotor rotation synchronization component can be easily removed.

更に、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御とし、位置制御ループゲイン特性が平坦となるように軸受制御系10を設計することにより、磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱の更なる抑制を図ることができるばかりか、軸受制御系10において速度制御帯域51がアキシャル方向ロータ回転速度fの近傍に属することが無いように、ロータ回転速度fに応じて速度制御ループゲインKを変化させることにより、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域における磁気軸受装置の発生擾乱の更なる抑制を図ることができる。 Furthermore, the bearing control system 10 is configured as a double loop control of a position control loop and a speed control loop, and the bearing control system 10 is designed so that the position control loop gain characteristic becomes flat, so that in the servo resonance frequency region of the magnetic bearing. In addition to being able to further suppress the amount of rotor swing and disturbance generated by the magnetic bearing device, the rotor is controlled so that the speed control band 51 does not belong to the vicinity of the axial rotor rotational speed f in the bearing control system 10. by varying the speed control loop gain K D according to the rotational speed f, without impairing the stability of the bearing control system 10, achieve a further reduction of generation disturbance of the magnetic bearing apparatus in the speed control system gain crossover frequency range be able to.

加えて、ロータ変位演算器12の出力と真のロータ変位との誤差及びギャップ長演算器13の出力と真の軸受電磁石ギャップ長との誤差を低減することにより、軸受制御系10における軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が低減し、磁気軸受装置の発生擾乱をより一層低減することができる。   In addition, the bearing force command in the bearing control system 10 is reduced by reducing the error between the output of the rotor displacement calculator 12 and the true rotor displacement and the error between the output of the gap length calculator 13 and the true bearing electromagnet gap length. In addition, the error between the radial bearing moment command and the true bearing force and true radial bearing moment acting on the rotor is reduced, and the occurrence of disturbance in the magnetic bearing device can be further reduced.

実施の形態6.
図15は、この発明に係る低擾乱化制御装置の実施の形態6を説明するためのブロック図であり、図7と同一の符号を付したものは、同一部分または相当部分を示す。
Embodiment 6 FIG.
FIG. 15 is a block diagram for explaining the sixth embodiment of the low disturbance control device according to the present invention, where the same reference numerals as those in FIG. 7 denote the same or corresponding parts.

次に、動作について説明する。図15では、簡単のためにラジアル方向並進運動及びアキシャル方向並進運動のフィードバック・ループは図示していないが、軸受制御系10の動作は実施の形態5と同一である。また、回転変位指令制御系40の動作は実施の形態3と同一である。   Next, the operation will be described. In FIG. 15, for the sake of simplicity, the feedback loop of radial translation and axial translation is not shown, but the operation of the bearing control system 10 is the same as that of the fifth embodiment. The operation of the rotational displacement command control system 40 is the same as that of the third embodiment.

上記実施の形態4における低擾乱化制御装置では、アキシャル方向ロータ回転角度θzに基づいた補正関数を、ロータ変位演算器12における上記式(15)に導入しており、上記式(15)におけるSi0(θz)(i=1、2、…、n)はロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出している。   In the low disturbance control device according to the fourth embodiment, a correction function based on the axial direction rotor rotation angle θz is introduced into the equation (15) in the rotor displacement calculator 12, and Si0 in the equation (15) is calculated. (Θz) (i = 1, 2,..., N) is derived based on the statics when the rotor is statically levitated without rotating around the rotation axis and the measured value of the supply current of the bearing electromagnet at this time. .

一方、上記実施の形態5による低擾乱化制御装置では、アキシャル方向ロータ回転角度θzに基づいた補正関数を、ロータ変位演算器12における上記式(15)及びギャップ長演算器13における上記式(22)に導入しており、上記式(15)におけるSi0(θz)(i=1、2、…、n)は磁気軸受の保護を目的としたタッチダウン軸受によってロータを機械的に拘束した場合のロータ変位情報に基づいて、上記式(22)におけるGi0(θz)(i=1、2、…、n)はロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出している。   On the other hand, in the low disturbance control device according to the fifth embodiment, the correction function based on the axial rotor rotation angle θz is expressed by the above equation (15) in the rotor displacement calculator 12 and the above equation (22) in the gap length calculator 13. ), And Si0 (θz) (i = 1, 2,..., N) in the above formula (15) is obtained when the rotor is mechanically constrained by a touchdown bearing for the purpose of protecting the magnetic bearing. Based on the rotor displacement information, Gi0 (θz) (i = 1, 2,..., N) in the above equation (22) represents the statics when the rotor is floated without rotating about the rotation axis, and at this time This is derived based on the measured value of the current supplied to the bearing magnet.

このように、ロータ変位演算器12における補正関数の導出方法には、タッチダウン軸受によってロータを機械的に拘束した場合のロータ変位情報に基づく方法と、ロータを静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づく方法の2通りがある。このうち、ロータ変位演算器12の出力と真のロータ変位との誤差をより低減できる方法は後者であると考えられる。   As described above, the correction function derivation method in the rotor displacement calculator 12 includes a method based on rotor displacement information when the rotor is mechanically constrained by a touchdown bearing, statics when the rotor is statically levitated, and There are two methods based on the measured value of the supply current of the bearing electromagnet at this time. Of these, the latter is considered to be the method that can further reduce the error between the output of the rotor displacement calculator 12 and the true rotor displacement.

そこで、この実施の形態6の低擾乱化制御装置では、アキシャル方向ロータ回転角度θzに基づいた補正関数を、ロータ変位演算器12における上記式(15)及びギャップ長演算器13における上記式(22)に導入し、上記式(15)におけるSi0(θz)(i=1、2、…、n)及び上記式(22)におけるGi0(θz)(i=1、2、…、n)の双方を、ロータを回転軸回りに回転させることなく静止浮上させた場合の静力学及びこのときの軸受電磁石供給電流計測値に基づいて導出する。   Therefore, in the low disturbance control device of the sixth embodiment, the correction function based on the axial direction rotor rotation angle θz is expressed by the above equation (15) in the rotor displacement calculator 12 and the above equation (22) in the gap length calculator 13. ) And both Si0 (θz) (i = 1, 2,..., N) in the above formula (15) and Gi0 (θz) (i = 1, 2,..., N) in the above formula (22). Is derived based on the statics when the rotor is floated without rotating around the rotation axis and the measured value of the supply current of the bearing electromagnet at this time.

このときの導出手順としては、第1に、実施の形態5と同様、タッチダウン軸受によってロータを機械的に拘束した場合のロータ変位情報に基づいて、ロータ変位演算器12における補正関数を導出する。   As a derivation procedure at this time, first, as in the fifth embodiment, a correction function in the rotor displacement calculator 12 is derived based on the rotor displacement information when the rotor is mechanically restrained by the touchdown bearing. .

第2に、これも実施の形態5と同様、ロータを静止浮上させた場合の静力学及びこのときの軸受電磁石への供給電流計測値に基づいて、ギャップ長演算器13における補正関数を導出する。実施の形態6による低擾乱化制御装置では更に、実施の形態4と同様にロータを静止浮上させた場合の静力学である上記式(16)、(17)に対して、真の軸受力及びラジアル方向軸受モーメントと、真の電磁石磁気吸引力の関係式である上記式(18)、真の電磁石磁気吸引力と、各電磁石に対する供給電流計測値及び各電磁石とロータ吸引面との真のギャップ長との関係式である上記式(19)、アキシャル方向ロータ回転角度θzに基づいた補正関数を導入し、真の軸受電磁石ギャップ長を算出する上記式(22)を代入し、所定のアキシャル方向ロータ回転角度θzで静止浮上させた場合の軸受電磁石供給電流計測値Ii(i=1、2、…、n)を適用することで、当該回転角度θzにおける真の5自由度ロータ変位Xt、Yt、Zt、θx、θyに関する非線形連立方程式を得る。   Second, as in the fifth embodiment, the correction function in the gap length calculator 13 is derived based on the statics when the rotor is statically levitated and the measured current supplied to the bearing electromagnet at this time. . Further, in the low disturbance control device according to the sixth embodiment, the true bearing force and the above formulas (16) and (17), which are the statics when the rotor is statically levitated as in the fourth embodiment, The above equation (18), which is a relational expression between the radial bearing moment and the true electromagnet magnetic attraction force, the true electromagnet magnetic attraction force, the measured current supplied to each electromagnet, and the true gap between each electromagnet and the rotor attraction surface The above formula (19), which is a relational expression with respect to the length, and a correction function based on the axial direction rotor rotation angle θz are introduced, and the above formula (22) for calculating the true bearing electromagnet gap length is substituted, and a predetermined axial direction By applying the bearing electromagnet supply current measurement value Ii (i = 1, 2,..., N) in the case of static levitation at the rotor rotation angle θz, the true 5-DOF rotor displacement Xt, Yt at the rotation angle θz is applied. , Z , You get θx, the non-linear simultaneous equations related to θy.

この非線形連立方程式の数値解を算出することで、所定の回転角度θzにおけるロータ変位演算器12の出力誤差、従って上記式(15)におけるSi0(θz)(i=1、2、…、n)をより高精度に導出することができる。   By calculating a numerical solution of this nonlinear simultaneous equation, the output error of the rotor displacement calculator 12 at a predetermined rotation angle θz, and thus Si0 (θz) (i = 1, 2,..., N) in the above equation (15). Can be derived with higher accuracy.

これにより、ロータ変位演算器12の出力と、真のロータ変位との誤差を更に低減することができ、従って、軸受制御系10における軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が減少し、磁気軸受装置の発生擾乱を更に一層低減することができる。   As a result, the error between the output of the rotor displacement calculator 12 and the true rotor displacement can be further reduced. Therefore, the bearing force command and the radial bearing moment command in the bearing control system 10 and the true force acting on the rotor can be reduced. The error between the bearing force and the true radial bearing moment is reduced, and the magnetic disturbance generated in the magnetic bearing device can be further reduced.

このように、この発明の実施の形態6による低擾乱化制御装置により、アキシャル方向ロータ回転速度が軸受制御系10の位置制御帯域内、サーボ共振周波数域及び位置制御帯域外のどの周波数領域にある場合でも、ロータ慣性主軸Pが磁気軸受回転軸と平行になり、ロータ動不釣合いに起因する擾乱モーメントを、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減することができる。   Thus, with the low disturbance control device according to Embodiment 6 of the present invention, the axial direction rotor rotational speed is in any frequency region within the position control band of the bearing control system 10, the servo resonance frequency region, and outside the position control band. Even in this case, the rotor inertia main shaft P is parallel to the magnetic bearing rotating shaft, and the disturbance moment caused by the rotor dynamic imbalance can be reduced in the entire frequency region including the servo resonance frequency region of the magnetic bearing.

また、回転変位指令制御系40におけるLPF42がロータ回転同期成分のみを抽出することで、ロータ回転同期成分以外の周波数域に属する信号成分を容易に除去することができる。   Further, since the LPF 42 in the rotational displacement command control system 40 extracts only the rotor rotation synchronization component, signal components belonging to a frequency region other than the rotor rotation synchronization component can be easily removed.

更に、軸受制御系10を位置制御ループと速度制御ループの2重ループ制御とし、位置制御ループゲイン特性が平坦となるように軸受制御系10を設計することにより、磁気軸受のサーボ共振周波数域におけるロータ振れ回り量及び磁気軸受装置の発生擾乱の更なる抑制を図ることができるのみならず、軸受制御系10において速度制御帯域51がアキシャル方向ロータ回転速度fの近傍に属することが無いように、ロータ回転速度fに応じて速度制御ループゲインKを変化させることにより、軸受制御系10の安定性を損なうことなく、速度制御系ゲイン交点周波数域における磁気軸受装置の発生擾乱の更なる抑制を図ることができる。 Furthermore, the bearing control system 10 is configured as a double loop control of a position control loop and a speed control loop, and the bearing control system 10 is designed so that the position control loop gain characteristic becomes flat, so that in the servo resonance frequency region of the magnetic bearing. In addition to being able to further suppress the rotor swing amount and the disturbance generated in the magnetic bearing device, in the bearing control system 10, the speed control band 51 does not belong to the vicinity of the axial rotor rotational speed f. by varying the speed control loop gain K D according to the rotor rotational speed f, without impairing the stability of the bearing control system 10, a further reduction of generation disturbance of the magnetic bearing apparatus in the speed control system gain crossover frequency range Can be planned.

加えて、ロータ変位演算器12の出力と真のロータ変位との誤差及びギャップ長演算器13の出力と真の軸受電磁石ギャップ長との誤差を低減することにより、軸受制御系10における軸受力指令及びラジアル方向軸受モーメント指令と、ロータに作用する真の軸受力及び真のラジアル方向軸受モーメントとの誤差が低減し、磁気軸受装置の発生擾乱をより一層低減することができる。   In addition, the bearing force command in the bearing control system 10 is reduced by reducing the error between the output of the rotor displacement calculator 12 and the true rotor displacement and the error between the output of the gap length calculator 13 and the true bearing electromagnet gap length. In addition, the error between the radial bearing moment command and the true bearing force and true radial bearing moment acting on the rotor is reduced, and the occurrence of disturbance in the magnetic bearing device can be further reduced.

この発明は、回転電機の回転体(ロータ)を非接触で支持する磁気軸受が発生する擾乱を、磁気軸受のサーボ共振周波数域を含めた全周波数領域において低減するのに利用される。   The present invention is used to reduce disturbance generated by a magnetic bearing that supports a rotating body (rotor) of a rotating electrical machine in a non-contact manner in all frequency ranges including a servo resonance frequency range of the magnetic bearing.

磁気軸受装置の構成要素であるロータ、電磁石、および変位検出手段の配置例を説明するための構成図である。It is a block diagram for demonstrating the example of arrangement | positioning of the rotor which is a component of a magnetic bearing apparatus, an electromagnet, and a displacement detection means. 傾斜磁極リング状ロータ1の周囲に6個の電磁石2a〜2f及び6個の変位検出手段3a〜3fを配置した磁気軸受装置の正面図(a)及びA−A断面図(b)である。They are a front view (a) and AA sectional view (b) of a magnetic bearing device in which six electromagnets 2a to 2f and six displacement detection means 3a to 3f are arranged around the inclined magnetic pole ring rotor 1. この発明に係る低擾乱化制御装置の実施の形態1を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 1 of the low disturbance control apparatus which concerns on this invention. 実施の形態1の低擾乱化制御装置を説明するためのロータ半径方向断面図である。FIG. 4 is a rotor radial direction cross-sectional view for explaining the low disturbance control device of the first embodiment. この発明に係る低擾乱化制御装置の実施の形態2を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 2 of the low disturbance control apparatus which concerns on this invention. 実施の形態2の低擾乱化制御装置を説明するためのロータ半径方向断面図である。FIG. 5 is a rotor radial direction cross-sectional view for explaining a low disturbance control device of a second embodiment. この発明に係る低擾乱化制御装置の実施の形態3を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 3 of the low disturbance control apparatus which concerns on this invention. 実施の形態3の低擾乱化制御装置を説明するための閉ループゲイン特性を示したボード線図である。FIG. 10 is a Bode diagram showing closed loop gain characteristics for explaining the low disturbance control device according to the third embodiment. 実施の形態3を説明するためのロータ断面図である。FIG. 6 is a rotor cross-sectional view for explaining a third embodiment. 実施の形態3を説明するためのロータ断面図である。FIG. 6 is a rotor cross-sectional view for explaining a third embodiment. 実施の形態3を説明するためのロータ断面図である。FIG. 6 is a rotor cross-sectional view for explaining a third embodiment. この発明に係る低擾乱化制御装置の実施の形態4を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 4 of the low disturbance control apparatus which concerns on this invention. この発明に係る低擾乱化制御装置の実施の形態5を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 5 of the low disturbance control apparatus which concerns on this invention. 実施の形態5の補正関数を導出するためのデータ計測における変位検出手段及びロータ変位検出面の配置を示した図である。FIG. 10 is a diagram showing an arrangement of displacement detection means and rotor displacement detection surfaces in data measurement for deriving a correction function according to the fifth embodiment. この発明に係る低擾乱化制御装置の実施の形態6を説明するためのブロック図である。It is a block diagram for demonstrating Embodiment 6 of the low disturbance control apparatus which concerns on this invention.

符号の説明Explanation of symbols

1 ロータ、2 電磁石、3 変位検出手段、10 軸受制御系、
11 電磁石、ロータ及び変位検出手段からなる構成要素、12 ロータ変位演算器、
13 ギャップ長演算器、14 電力分配器、15 PID制御器、16 PI制御器、17 P制御器、18 微分器、20 並進変位指令制御系、
21,41 回転座標変換器、22,42 LPF、23,43 積分器、
24 固定座標変換器、25,45 固定座標変換器/位相変換器、
30 ロータ静不釣合いに起因する振れ回り力、
31 ラジアル方向軸受力指令の回転座標系表現、
32 ラジアル方向並進変位指令の回転座標系表現、33,63 ロータ回転方向、
34 ロータ幾何中心位置の移動方向ベクトル、40 回転変位指令制御系、
50 位置制御閉ループゲイン特性、51 速度制御閉ループゲイン特性、
60 ロータの動不釣合いに起因する振れ回りモーメント、
61 ラジアル方向軸受モーメント指令の回転座標系表現、
62 ラジアル方向回転変位指令の回転座標系表現、
64 ロータ慣性主軸の回転方向ベクトル、70 変位検出手段、
71 ロータ変位検出方向、72 ロータ変位検出面、73 ロータ回転軸(Z軸)。
1 rotor, 2 electromagnets, 3 displacement detection means, 10 bearing control system,
11 components comprising an electromagnet, a rotor and displacement detection means, 12 a rotor displacement calculator,
13 Gap length calculator, 14 Power distributor, 15 PID controller, 16 PI controller, 17 P controller, 18 Differentiator, 20 Translational displacement command control system,
21, 41 Rotary coordinate converter, 22, 42 LPF, 23, 43 integrator,
24 fixed coordinate converter, 25, 45 fixed coordinate converter / phase converter,
30 Swinging force due to rotor unbalance,
31 Rotary coordinate system representation of radial direction bearing force command,
32 Rotary coordinate system representation of radial direction translational displacement command, 33, 63 Rotor rotation direction,
34 Rotor geometric center position moving direction vector, 40 rotational displacement command control system,
50 Position control closed loop gain characteristics, 51 Speed control closed loop gain characteristics,
60 Swinging moment due to rotor dynamic imbalance,
61 Rotating coordinate system representation of radial direction bearing moment command,
62 Rotational coordinate system representation of radial direction rotational displacement command,
64 rotational direction vector of the rotor inertial spindle, 70 displacement detecting means,
71 Rotor displacement detection direction, 72 Rotor displacement detection surface, 73 Rotor rotation axis (Z axis).

Claims (12)

ロータを電磁力によって浮上させる複数の電磁石からなる磁気軸受の、ラジアル方向をX軸及びY軸、中心軸をZ軸とする固定座標系において、上記ロータの幾何中心軸が上記ラジアル方向に並進運動するラジアル方向並進変位及び上記ロータの幾何中心軸が上記固定座標系のアキシャル方向(Z軸方向)に並進運動するアキシャル方向並進変位の3自由度ロータ変位と、上記ロータの幾何中心軸が上記固定座標系のZ軸を含む平面内で回動するラジアル方向回転変位の2自由度ロータ変位からなる5自由度ロータ変位及びアキシャル方向ロータ回転角度に関係する情報(ロータ変位情報)を検出する変位検出手段と、
上記変位検出手段の出力を、上記5自由度ロータ変位に変換するロータ変位演算器と、 上記ロータ変位演算器の出力に基づいて、上記ロータを磁気軸受の所定の位置に制御するための軸受力となるラジアル方向及びアキシャル方向軸受力指令、上記ロータのラジアル方向回転変位を制御する軸受力となるラジアル方向軸受モーメント指令を出力する制御指令出力手段と、
上記ロータ変位演算器の出力に基づいて、上記電磁石各々と上記ロータとのギャップ長を出力するギャップ長演算器と
上記ギャップ長演算器の出力及び上記制御指令出力手段の出力に基づいて、上記制御指令出力手段から出力されるラジアル方向及びアキシャル方向軸受力指令、及びラジアル方向軸受モーメント指令を実現するための上記電磁石各々に対する供給電流を分配して供給する電力分配器とを備えた軸受制御系、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記軸受制御系における制御指令出力手段が出力するラジアル方向軸受力指令を、上記ロータの幾何中心軸をZ軸、上記Z軸に垂直な方向をX軸及びY軸としてロータと共にZ軸まわりに回転する回転座標系表現に変換する回転座標変換器と、
上記回転座標変換器の出力のロータ回転同期成分を抽出するLPFと、
上記LPFの出力を積分して、ラジアル方向に並進させる回転座標系のラジアル方向並進変位を出力する積分器と、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記積分器が出力するラジアル方向並進変位を上記固定座標系表現に変換して上記ロータの重心位置が上記回転座標系原点に位置するように上記ロータの幾何中心軸をラジアル方向に並進させるラジアル方向並進変位指令を算出する固定座標変換器を備えた並進変位指令制御系、
を備え、
上記固定座標系に変換したラジアル方向並進変位指令により、上記軸受制御系が上記電磁石を制御して上記ロータの浮上位置を所定位置に制御することを特徴とする磁気軸受の低擾乱化制御装置。
In a fixed coordinate system in which a radial direction is an X axis and a Y axis and a central axis is a Z axis of a magnetic bearing composed of a plurality of electromagnets that levitate a rotor by electromagnetic force, the geometric central axis of the rotor translates in the radial direction. 3 degree-of-freedom rotor displacement in which the translational displacement in the radial direction and the geometrical center axis of the rotor translate in the axial direction (Z-axis direction) of the fixed coordinate system, and the geometrical center axis of the rotor are fixed. Displacement detection for detecting information (rotor displacement information) related to 5-degree-of-freedom rotor displacement and axial-direction rotor-rotation angle consisting of 2-degree-of-freedom rotor displacement of radial-direction rotational displacement rotating in a plane including the Z-axis of the coordinate system Means,
A rotor displacement calculator for converting the output of the displacement detection means into the 5-degree-of-freedom rotor displacement; and a bearing force for controlling the rotor to a predetermined position of the magnetic bearing based on the output of the rotor displacement calculator. A control command output means for outputting a radial direction and axial direction bearing force command, and a radial direction bearing moment command serving as a bearing force for controlling the radial rotational displacement of the rotor;
Based on the output of the rotor displacement calculator, a gap length calculator that outputs the gap length between each of the electromagnets and the rotor ;
The electromagnet for realizing the radial direction and axial direction bearing force command and the radial direction bearing moment command output from the control command output means based on the output of the gap length calculator and the output of the control command output means. A bearing control system including a power distributor that distributes and supplies a supply current to each ;
Based on the rotor rotation angle in the axial direction of the rotor, the radial bearing force command output by the control command output means in the bearing control system is expressed as follows: the geometric center axis of the rotor is the Z axis, and the direction perpendicular to the Z axis is X A rotary coordinate converter for converting into a rotary coordinate system representation that rotates about the Z axis with the rotor as an axis and a Y axis;
LPF for extracting the rotor rotation synchronization component of the output of the rotary coordinate converter;
An integrator that integrates the output of the LPF and outputs a radial translational displacement of a rotating coordinate system that translates in the radial direction;
Based on the axial rotor rotation angle of the rotor, the radial translation displacement output from the integrator is converted into the fixed coordinate system representation so that the center of gravity of the rotor is located at the origin of the rotary coordinate system. translational displacement command control system that includes a fixed coordinate converter for calculating a radial translational displacement command for translating the geometric center axis of the radial direction,
With
A magnetic bearing low disturbance control device, wherein the bearing control system controls the electromagnet to control the flying position of the rotor to a predetermined position by a radial translational displacement command converted into the fixed coordinate system.
ロータを電磁力によって浮上させる複数の電磁石からなる磁気軸受の、ラジアル方向をX軸及びY軸、中心軸をZ軸とする固定座標系において、上記ロータの幾何中心軸が上記ラジアル方向に並進運動するラジアル方向並進変位及び上記ロータの幾何中心軸が上記固定座標系のアキシャル方向(Z軸方向)に並進運動するアキシャル方向並進変位の3自由度ロータ変位と、上記ロータの幾何中心軸が上記固定座標系のZ軸を含む平面内で回動するラジアル方向回転変位の2自由度ロータ変位からなる5自由度ロータ変位及びアキシャル方向ロータ回転角度に関係する情報(ロータ変位情報)を検出する変位検出手段と、
上記変位検出手段の出力を、上記5自由度ロータ変位に変換するロータ変位演算器と、
上記ロータ変位演算器の出力に基づいて、上記ロータを磁気軸受の所定の位置に制御するための軸受力となるラジアル方向及びアキシャル方向軸受力指令、上記ロータのラジアル方向回転変位を制御する軸受力となるラジアル方向軸受モーメント指令を出力する制御指令出力手段と、
上記ロータ変位演算器の出力に基づいて、上記電磁石各々と上記ロータとのギャップ長を出力するギャップ長演算器と
上記ギャップ長演算器の出力及び上記制御指令出力手段の出力に基づいて、上記制御指令出力手段から出力されるラジアル方向及びアキシャル方向軸受力指令、及びラジアル方向軸受モーメント指令を実現するための上記電磁石各々に対する供給電流を分配して供給する電力分配器とを備えた軸受制御系、
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記軸受制御系における制御指令出力手段が出力するラジアル方向軸受モーメント指令を、上記ロータの幾何中心軸をZ軸、上記Z軸に垂直な方向をX軸及びY軸としてロータと共にZ軸周りに回転する回転座標系表現に変換する回転座標変換器と、
上記回転座標変換器の出力のロータ回転同期成分を抽出するLPFと、
上記LPFの出力を積分してラジアル方向に回転させる回転座標系のラジアル方向回転変位を出力する積分器
上記ロータのアキシャル方向ロータ回転角度に基づいて、上記積分器から出力されるラジアル方向回転変位を上記固定座標系表現に変換して上記ロータの慣性主軸が上記磁気軸受けの中心軸と平行になるように上記ロータの幾何中心軸を上記固定座標系のZ軸を含む平面内で回動させるラジアル方向回転変位指令を算出する固定座標変換器とを備えた回転変位指令制御系、
を備え、
上記固定座標系に変換したラジアル方向回転変位指令により、上記軸受制御系が上記電磁石を制御して上記ロータの浮上位置を所定位置に制御することを特徴とする磁気軸受の低擾乱化制御装置。
In a fixed coordinate system in which a radial direction is an X axis and a Y axis and a central axis is a Z axis of a magnetic bearing composed of a plurality of electromagnets that levitate a rotor by electromagnetic force, the geometric central axis of the rotor translates in the radial direction. 3 degree-of-freedom rotor displacement in which the translational displacement in the radial direction and the geometrical center axis of the rotor translate in the axial direction (Z-axis direction) of the fixed coordinate system, and the geometrical center axis of the rotor are fixed. Displacement detection that detects information (rotor displacement information) related to the 5-degree-of-freedom rotor displacement and the axial-direction rotor-rotation angle consisting of a 2-degree-of-freedom rotor displacement of a radial-direction rotational displacement that rotates in a plane including the Z-axis of the coordinate system Means,
A rotor displacement calculator for converting the output of the displacement detection means into the 5-degree-of-freedom rotor displacement;
Based on the output of the rotor displacement calculator, a radial and axial direction bearing force command that becomes a bearing force for controlling the rotor to a predetermined position of the magnetic bearing, and a bearing force for controlling the radial rotational displacement of the rotor Control command output means for outputting a radial bearing moment command,
Based on the output of the rotor displacement calculator, a gap length calculator that outputs the gap length between each of the electromagnets and the rotor ;
The electromagnet for realizing the radial direction and axial direction bearing force command and the radial direction bearing moment command output from the control command output means based on the output of the gap length calculator and the output of the control command output means. A bearing control system including a power distributor that distributes and supplies a supply current to each ;
Based on the rotor rotational angle in the axial direction of the rotor, the radial bearing moment command output by the control command output means in the bearing control system is expressed as follows: the geometric center axis of the rotor is the Z axis, and the direction perpendicular to the Z axis is X A rotary coordinate converter that converts the rotary coordinate system to rotate around the Z axis together with the rotor as an axis and a Y axis;
LPF for extracting the rotor rotation synchronization component of the output of the rotary coordinate converter;
An integrator for outputting a radial rotational displacement of the rotating coordinate system that rotates in the radial direction by integrating the output of the LPF,
Based on the rotor rotational angle in the axial direction of the rotor, the radial rotational displacement output from the integrator is converted into the fixed coordinate system representation so that the inertial main shaft of the rotor is parallel to the central axis of the magnetic bearing. A rotational displacement command control system comprising: a fixed coordinate converter that calculates a radial rotational displacement command for rotating the geometric center axis of the rotor in a plane including the Z axis of the fixed coordinate system;
With
A magnetic bearing low disturbance control device, wherein the bearing control system controls the electromagnet to control the flying position of the rotor to a predetermined position by a radial rotational displacement command converted into the fixed coordinate system.
上記LPFの通過周波数帯域を大略1Hz程度の狭帯域としたことを特徴とする請求項1または2記載の磁気軸受の低擾乱化制御装置。 3. A low disturbance control device for a magnetic bearing according to claim 1, wherein the LPF has a narrow frequency band of about 1 Hz. 上記並進変位指令制御系の固定座標変換器が、上記積分器の出力を固定座標系表現に変換する際に、上記軸受制御系の位相特性を考慮して、アキシャル方向ロータ回転速度に応じて所定の位相量を進めて固定座標系表現に変換した信号を出力することを特徴とする請求項1記載の磁気軸受の低擾乱化制御装置。 When the fixed coordinate converter of the translational displacement command control system converts the output of the integrator into a fixed coordinate system expression, the phase characteristic of the bearing control system is taken into account and the predetermined value is determined according to the axial direction rotor rotational speed. The magnetic bearing low disturbance control device according to claim 1, wherein a signal obtained by advancing the phase amount is converted into a fixed coordinate system representation and output. 上記回転変位指令制御系の固定座標変換器が、上記積分器の出力を固定座標系表現に変換する際に、上記軸受制御系の位相特性を考慮して、アキシャル方向ロータ回転速度に応じて所定の位相量を進めて固定座標系表現に変換した信号を出力することを特徴とする請求項2記載の磁気軸受の低擾乱化制御装置。 When the fixed coordinate converter of the rotational displacement command control system converts the output of the integrator into a fixed coordinate system expression, the phase characteristic of the bearing control system is taken into consideration and the predetermined value is determined according to the axial direction rotor rotational speed. 3. A low disturbance control device for a magnetic bearing according to claim 2, wherein a signal obtained by advancing the phase amount is converted into a fixed coordinate system representation. 上記軸受制御系の制御指令出力手段は、5自由度変位指令と5自由度ロータ変位の偏差に基づいて、5自由度の速度指令を出力する位置制御ループと上記5自由度ロータ変位の微分値と5自由度の速度指令との偏差に基づいて速度制御を行う速度制御ループの2重ループ制御からなり、位置制御閉ループゲイン特性が平坦となるように制御することを特徴とする請求項1または2記載の磁気軸受の低擾乱化制御装置。 The control command output means of the bearing control system includes a position control loop for outputting a speed command of 5 degrees of freedom based on a deviation between the 5 degrees of freedom displacement command and the 5 degrees of freedom rotor displacement and a differential value of the 5 degrees of freedom rotor displacement 2. A double loop control of a speed control loop that performs speed control based on a deviation between the speed command of 5 degrees of freedom and a position control closed loop gain characteristic is controlled to be flat. 2. A low disturbance control device for a magnetic bearing according to 2. 上記速度制御の帯域がアキシャル方向ロータ回転速度の近傍に属することが無いように、上記アキシャル方向ロータ回転速度に応じて速度制御ループゲインを変化させることを特徴とする請求項6記載の磁気軸受の低擾乱化制御装置。 7. The magnetic bearing according to claim 6, wherein a speed control loop gain is changed in accordance with the axial rotor rotational speed so that the speed control band does not belong to the vicinity of the axial rotor rotational speed. Low disturbance control device. 上記ロータ変位演算器は、真のロータ変位に対する出力5自由度ロータ変位の誤差を低減する補正関数を有することを特徴とする請求項1または2記載の磁気軸受の低擾乱化制御装置。 3. The magnetic bearing low disturbance control device according to claim 1, wherein the rotor displacement calculator has a correction function for reducing an error of the rotor displacement with an output of 5 degrees of freedom with respect to the true rotor displacement. 上記磁気軸受の保護を目的としたタッチダウン軸受を設け、上記タッチダウン軸受によって上記ロータを機械的に拘束した場合のロータ変位情報に基づいて、上記補正関数を導出することを特徴とする請求項8記載の磁気軸受の低擾乱化制御装置。 A touchdown bearing for protecting the magnetic bearing is provided, and the correction function is derived based on rotor displacement information when the rotor is mechanically constrained by the touchdown bearing. 8. A low disturbance control device for a magnetic bearing according to claim 8. 上記ロータを静止浮上させた場合の静力学及び軸受電磁石供給電流計測値に基づいて、上記補正関数を導出したことを特徴とする請求項8記載の磁気軸受の低擾乱化制御装置。 The magnetic bearing low disturbance control device according to claim 8, wherein the correction function is derived based on statics when the rotor is statically levitated and a measured value of the supply current of the bearing electromagnet. 上記ギャップ長演算器は、真のギャップ長に対する出力ギャップ長の誤差を低減する補正関数を有することを特徴とする請求項1または2記載の磁気軸受の低擾乱化制御装置。 3. The magnetic bearing low disturbance control device according to claim 1, wherein the gap length calculator has a correction function for reducing an error of an output gap length with respect to a true gap length. 上記ロータを静止浮上させた場合の静力学及び軸受電磁石供給電流計測値に基づいて、上記補正関数を導出することを特徴とする請求項11記載の磁気軸受の低擾乱化制御装置。 The magnetic bearing low disturbance control device according to claim 11, wherein the correction function is derived based on statics when the rotor is statically levitated and a measured value of the supply current of the bearing electromagnet.
JP2003427227A 2003-12-24 2003-12-24 Magnetic bearing low disturbance control device Expired - Lifetime JP4344601B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003427227A JP4344601B2 (en) 2003-12-24 2003-12-24 Magnetic bearing low disturbance control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003427227A JP4344601B2 (en) 2003-12-24 2003-12-24 Magnetic bearing low disturbance control device

Publications (2)

Publication Number Publication Date
JP2005188545A JP2005188545A (en) 2005-07-14
JP4344601B2 true JP4344601B2 (en) 2009-10-14

Family

ID=34786563

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003427227A Expired - Lifetime JP4344601B2 (en) 2003-12-24 2003-12-24 Magnetic bearing low disturbance control device

Country Status (1)

Country Link
JP (1) JP4344601B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102425555A (en) * 2011-11-11 2012-04-25 北京中科科仪技术发展有限责任公司 Method for obtaining radial suspension centers of rotor of magnetic molecular pump
CN105588507A (en) * 2016-03-03 2016-05-18 武汉理工大学 Magnetic bearing position measurement error calculation method under coordinate transformation

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007263143A (en) * 2006-03-27 2007-10-11 Jtekt Corp Magnetic bearing device for machine tool
JP5673086B2 (en) * 2010-12-27 2015-02-18 株式会社明電舎 Motor control device
JP5673087B2 (en) * 2010-12-27 2015-02-18 株式会社明電舎 Motor control device
CN112347717A (en) * 2020-11-05 2021-02-09 南京工业大学 Novel centralized same-frequency disturbance suppression method for magnetic bearing

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102425555A (en) * 2011-11-11 2012-04-25 北京中科科仪技术发展有限责任公司 Method for obtaining radial suspension centers of rotor of magnetic molecular pump
CN105588507A (en) * 2016-03-03 2016-05-18 武汉理工大学 Magnetic bearing position measurement error calculation method under coordinate transformation
CN105588507B (en) * 2016-03-03 2018-06-22 武汉理工大学 Magnetic suspension bearing errors in position measurement computational methods under a kind of coordinate transform

Also Published As

Publication number Publication date
JP2005188545A (en) 2005-07-14

Similar Documents

Publication Publication Date Title
Fang et al. Decoupling control of magnetically suspended rotor system in control moment gyros based on an inverse system method
US6353273B1 (en) Hybrid foil-magnetic bearing
US5313399A (en) Adaptive synchronous vibration suppression apparatus
US6965181B1 (en) Hybrid foil-magnetic bearing with improved load sharing
Hsu et al. Nonlinear control of a 3-pole active magnetic bearing system
JP6351400B2 (en) Improved active magnetic bearing control system
JPH0242125B2 (en)
CN107102554B (en) Method for suppressing unbalanced vibration of magnetic suspension spherical flywheel
Chen et al. Experimental validation of a current-controlled three-pole magnetic rotor-bearing system
Yoon et al. Unbalance compensation for AMB systems with input delay: An output regulation approach
US8853905B2 (en) Radial direction controller and magnetic bearing apparatus utilizing same
JP4344601B2 (en) Magnetic bearing low disturbance control device
Xu et al. Dynamic modeling of the active magnetic bearing system operating in base motion condition
Zheng et al. Improving dynamic response of AMB systems in control moment gyros based on a modified integral feedforward method
Jang et al. Improvement of high-speed stability of an aerostatic bearing-rotor system using an active magnetic bearing
JP2004132441A (en) Magnetic bearing device, through-flow fan device for excimer laser using this magnetic bearing device, program for making computer execute feedback control of magnetic bearing and computer-readable recording medium with program for making computer execute feedback control of magnetic bearing recorded thereon
US5376871A (en) Method of controlling position of rotary shaft in magnetic bearing
JP4500953B2 (en) Magnetic bearing control device
KR100575556B1 (en) Method for controlling vibration of a magnetically suspended spindle
Bian et al. Synchronous disturbance suppression of active magnetic bearing rotor systems using variable period adaptive control algorithm
Yoon et al. An output regulation approach to rotor autobalancing in active magnetic bearing systems with input delay
JP2004324895A (en) Static pressure magnetic combined bearing
JP7055591B2 (en) Magnetic bearing controller
JPH10116120A (en) Fine movement positioning controller
CN114322971B (en) Magnetic suspension rotor same-frequency vibration force suppression method based on biquad generalized integrator

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051110

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20051110

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20081125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090407

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090604

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090707

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090713

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

Free format text: PAYMENT UNTIL: 20120717

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4344601

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20120717

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20130717

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term