CN112432634B - Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation - Google Patents
Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation Download PDFInfo
- Publication number
- CN112432634B CN112432634B CN202011380266.5A CN202011380266A CN112432634B CN 112432634 B CN112432634 B CN 112432634B CN 202011380266 A CN202011380266 A CN 202011380266A CN 112432634 B CN112432634 B CN 112432634B
- Authority
- CN
- China
- Prior art keywords
- frequency
- harmonic
- rotor
- displacement
- force
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001629 suppression Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000009466 transformation Effects 0.000 title claims description 35
- 239000000725 suspension Substances 0.000 claims abstract description 41
- 238000011426 transformation method Methods 0.000 claims abstract description 14
- 238000006073 displacement reaction Methods 0.000 claims description 49
- 230000001360 synchronised effect Effects 0.000 claims description 28
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 230000000452 restraining effect Effects 0.000 abstract 2
- 238000004422 calculation algorithm Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 10
- 238000012546 transfer Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 238000003780 insertion Methods 0.000 description 4
- 230000037431 insertion Effects 0.000 description 4
- 238000002955 isolation Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000003252 repetitive effect Effects 0.000 description 3
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000005339 levitation Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
- G01C19/02—Rotary gyroscopes
- G01C19/04—Details
- G01C19/16—Suspensions; Bearings
- G01C19/24—Suspensions; Bearings using magnetic or electrostatic fields
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
The invention discloses a magnetic suspension rotor harmonic vibration force suppression method based on a multi-synchronous coordinate transformation method, which comprises the following steps of: firstly, a magnetic suspension rotor dynamic model containing mass unbalance and sensor harmonic is established, and then a magnetic suspension rotor harmonic vibration force suppression method based on MSRFT is adopted. The MSRFT can accurately restrain harmonic vibration force, and the controller is used for restraining vibration in X and Y directions at the same time, so that hardware computing resources are reduced, and the restraining dynamic process is accelerated. Meanwhile, the absolute stability of the system in a larger frequency range can be ensured by introducing a phase compensation angle. The MSRFT controller has a simple structure, is very convenient in practical application, can inhibit harmonic vibration force in the magnetic suspension rotor, and is suitable for inhibiting the harmonic vibration force of a magnetic suspension rotor system with mass unbalance and sensor harmonic.
Description
Technical Field
The invention relates to the technical field of suppression of harmonic vibration force of a magnetic suspension rotor, in particular to a suppression method of harmonic vibration force of a magnetic suspension rotor based on a multi-synchronous rotation coordinate transformation method, which is used for suppressing the harmonic vibration force of a rotor system of a magnetic suspension control moment gyroscope in a full working rotation speed range and providing technical support for application of the magnetic suspension control moment gyroscope on an ultra-static satellite platform and an ultra-stable satellite platform.
Background
The control moment gyroscope has the advantages of large output moment, high response speed and the like, and becomes a key attitude control inertia actuating mechanism of a high-performance satellite. Compared with the traditional mechanical bearing, the active magnetic bearing can realize the non-contact support of the rotor, so that the active magnetic bearing has the advantages of no friction, no need of lubrication, high rotating speed and the like, and the active vibration is controllable. Compared with the traditional mechanical gyroscope, the magnetic suspension control moment gyroscope has the characteristic of micro-vibration and is widely used on high-performance satellites.
However, some vibration is inevitable in the magnetic suspension control moment gyro due to machining errors of the rotor and measurement errors of the sensor. The two main vibration sources are rotor mass unbalance and sensor harmonic waves, the mass unbalance is represented by the misalignment of the geometric center and the mass center of the rotor, when the rotor rotates at high speed, larger centrifugal force can be generated, the amplitude of the centrifugal force is increased along with the increase of the rotating speed, and the frequency is the same as the rotating frequency; the harmonic wave of the sensor comes from the roundness error of the measured surface, so that the output signal of the displacement sensor contains harmonic interference with the same frequency and frequency multiplication of the rotating speed, and the electromagnetic coil of the magnetic bearing of the active magnetic suspension control moment gyroscope generates harmonic current with the same frequency and frequency multiplication, thereby generating harmonic vibration force. These harmonic vibration forces are transmitted to the satellite platform through the base, so that the pointing accuracy and attitude stability of the satellite are seriously affected.
In order to solve the vibration, a traditional method is to use a mechanical vibration isolation device, and a vibration source is isolated by adding a physical device, on one hand, the vibration isolation device is expensive, large in size and inconvenient to maintain, and the vibration isolation performance is inevitably influenced along with the increase of the service time; on the other hand, the vibration isolation device does not eliminate the vibration, but converts the low-frequency high-amplitude vibration into the high-frequency low-amplitude vibration, and the energy of the vibration is not reduced. Therefore, the active vibration control is carried out on the magnetic suspension control moment gyroscope, and the vibration is eliminated through a vibration control algorithm, so that the method for effectively solving the vibration is provided.
The active vibration control algorithm of the magnetic suspension control moment gyroscope mainly comprises a wave trap, a resonance controller, a repetitive controller, an LMS (least mean square) algorithm and the like, wherein the repetitive controller is an effective method for processing harmonic interference, but the repetitive control has the defect of low response speed, the LMS algorithm has large calculation amount and is difficult to realize, and although the wave trap and the resonance controller have simple structures, a plurality of controllers are required to be connected in parallel when a plurality of frequencies are processed, so that the calculation burden of a system is increased.
Disclosure of Invention
The purpose of the invention is as follows: the method can realize the complete inhibition of the vibration force by using the bearing force constructed by current and displacement according to an electromagnetic force model as the input of a control algorithm, can simultaneously inhibit the vibration force in two directions by using the orthogonal characteristics of signals in the X and Y directions, reduces the calculation resources, improves the dynamic response speed, and realizes the absolute stability of the full working frequency band by introducing a phase compensation angle to compensate the phase in different frequency bands.
The technical scheme adopted by the invention is as follows: a magnetic suspension rotor harmonic vibration force suppression method based on a multi-synchronous coordinate transformation method comprises the following steps:
step (1) establishing a full-active magnetic suspension rotor dynamic model containing mass unbalance and sensor harmonic
The invention is applied to an active magnetic bearing system in a magnetic suspension control moment gyro. Let N be the geometric center of the stator, NXY be the inertial coordinate system, C and O be the center of mass and the geometric center of the rotor, respectively, and O ε η be the rotational coordinate system. The method mainly aims at the suppression of the vibration force of the radial translation two-degree-of-freedom, so that the modeling of the translation two-degree-of-freedom is only considered. From newton's second law, the following kinetic equation is obtained:
where m is the rotor mass, X (t), Y (t) represent the translational displacements of the rotor mass center in the X and Y directions, respectively, f x (t),f y (t) resultant forces of bearing forces in X-direction and Y-direction, respectively, f ax (t),f bx (t),f ay (t),f by (t) is the bearing force of the four pairs of radial magnetic bearings. When the rotor moves in a small displacement, the nonlinear bearing force can be approximately linearized and can be expressed as follows:
wherein K i ,K h Current stiffness and displacement stiffness, i, respectively ax (t),i bx (t),i ay (t),i by (t) four pairs of radial magnetic bearing coil currents, x a (t),x b (t),y a (t),y b (t) displacement under a bearing coordinate system; subscripts a, B denote the ends of a, B of the rotor system;
because the mass unbalance makes geometric center and center of mass misalignment, displacement sensor measures as rotor geometric center displacement, has following relation:
X(t)＝x(t)+Θ x (t)
Y(t)＝y(t)+Θ y (t)
wherein X (t), Y (t) represent the displacement of the geometric center, respectively, and can be obtained by the geometric relationship:
Θ x (t),Θ y (t) represents the amount of unbalance, which has the following form:
Θ x (t)＝ecos(Ωt+χ)
Θ y (t)＝esin(Ωt+χ)
e is the magnitude of the unbalance, χ is the initial phase, Ω is the rotor rotational speed. It can be seen that the unbalanced mass can generate co-frequency interference signals, so that the active magnetic bearing generates co-frequency vibration force.
Because the displacement sensor has sensor harmonic waves, the output signals of the four pairs of displacement sensors are not the real displacement signals of the rotor, and the displacement sensors provide the following signals:
x as (t),x bs (t),y as (t),y bs (t) is the displacement sensor output signal; d is a radical of axs (t),d bxs (t),d ays (t),d bys (t) is the sensor harmonic interference signal, which can be expressed as follows:
wherein s is ai ,s bi Is the harmonic amplitude of the sensor, χ i The initial phase is I, the harmonic order is I, and the sensor harmonic can generate an interference signal with the same frequency and frequency multiplication with the rotating speed, so that the magnetic bearing system generates harmonic vibration force.
Step (2) designs a magnetic suspension rotor harmonic vibration force suppression method based on a multi-synchronous rotation coordinate transformation method
The controller takes the same-frequency vibration force and frequency multiplication current as input, is connected into an original closed-loop system in a parallel mode, and the output of the controller is fed back to the power amplifier input end of the original control system, and the design of the module mainly comprises the following two aspects:
(1) the multi-synchronous rotation coordinate transformation method comprises the following steps: according to different forms of common-frequency and frequency-doubled vibration forces generated by an actual magnetic suspension rotor system, common-frequency vibration comprises a current rigidity force and a displacement rigidity force, and frequency-doubled vibration only comprises the current rigidity force; according to a model of system electromagnetic force, current and displacement are used for constructing same-frequency vibration force to be used as input of a multiple synchronous rotation coordinate transformation same-frequency force suppression part, and frequency multiplication vibration can be suppressed by suppressing frequency multiplication current, so that the current is used as input of the frequency multiplication force suppression part;
(2) obtaining the stable condition of a system with multi-synchronous rotation coordinate transformation through theoretical analysis and proof; and designing a corresponding phase compensation angle according to the closed-loop characteristic of the actual magnetic suspension rotor system, and realizing the absolute stability of the system in a working rotating speed range through the phase compensation angle.
Further, the vibration force suppression algorithm in the step (2) is as follows:
the synchronous rotating coordinate transformation for realizing harmonic suppression mainly comprises three parts, namely, a disturbing signal is transformed from a static coordinate system to a rotating coordinate system through synchronous rotating coordinate transformation, at the moment, a same frequency/frequency multiplication disturbing signal is represented as a direct current quantity under the rotating coordinate system, then the direct current quantity is identified through low-pass filtering, and finally, the same frequency/frequency multiplication disturbing signal under the static coordinate system is obtained through synchronous rotating coordinate inverse transformation. The output is fed back to a closed loop structure formed by the output end of the power amplifier, so that the same frequency/frequency multiplication vibration force can be effectively inhibited.
Because the phases of all frequency bands of the system are inconsistent, a phase compensation angle needs to be introduced for phase compensation in order to ensure absolute stability of the system in the full working frequency band range. The transformation equation of the synchronous rotation coordinate with the phase compensation angle is as follows:
wherein u 1 (t),u 2 (t) is the input signal, u dc1 (t),u dc2 (t) is a synchronous rotation coordinate transformation output signal; omega is the rotational speed of the rotor and,
the phase compensation angle is used for ensuring the stability of a closed loop, i is a positive integer, when the harmonic frequency of an input signal is i times of omega, the input signal is output as a direct current signal after synchronous rotating coordinate transformation, harmonic components in an original signal can be extracted through low-pass filtering, and a low-pass filter is provided with the following transfer functions:
where k is the low pass filter gain factor, ω c Is the cut-off frequency;
suppose the signal after low-pass filtering is
Let u be dc1 (t),u dc2 (t),
Respectively u as the laplace transform dc1 (s),u dc2 (s),
The following equation holds:
the above formula is rewritten as follows:
inverse laplace transform on the above equation can yield the following differential equation:
let κ = k ω c The above equation is written in matrix form as follows:
the low-pass filtered signal needs to be inversely transformed to the original stationary coordinate system by the synchronous rotating coordinate, and the transformation relation is expressed as follows:
x 1 (t),x 2 (t) is the output signal after the synchronous rotation coordinate inverse transformation;
the following state space expression can be solved from the above equations:
the transfer function matrix is derived from the state space expression:
the corresponding input and output relationship is as follows:
u 1 (t),u 2 (t) are orthogonal sinusoidal signals, assuming the form:
A m and γ represents the amplitude and initial phase of the signal, respectively;
from the above formula, the following relationship exists:
iΩu 1 (s)＝su 2 (s)
an equivalent single input transfer function can be derived from the above equation:
let s = j ω, it can be obtained that when ω = i Ω, there are:
when ω is c And k, the synchronous rotation coordinate transformation can generate large gain at harmonic frequency, and the harmonic component can be effectively inhibited after the synchronous rotation coordinate transformation and the original system form a closed loop. Thereby eliminating the harmonic vibration forces of the system.
The basic principle of the invention is as follows: the magnetic suspension control moment gyroscope is supported by a magnetic suspension bearing, and for a magnetic suspension rotor, the main vibration sources are mass unbalance and sensor harmonic waves. Due to mass unbalance and the existence of sensor harmonic waves, harmonic vibration force is contained in the magnetic suspension rotor system. The harmonic vibration force is transmitted to the spacecraft through the base, and the performance of the spacecraft platform is seriously influenced. The invention provides a harmonic vibration force suppression algorithm based on a multi-synchronous rotation coordinate transformation method by establishing a magnetic suspension rotor dynamic model containing mass unbalance and sensor harmonic aiming at the magnetic suspension rotor harmonic vibration force of a magnetic suspension control moment gyroscope. The orthogonal characteristic of the output signals of the displacement sensors in the X direction and the Y direction is utilized, one controller is used for simultaneously inhibiting the vibration force in two directions, a plurality of algorithms with different frequencies are connected in parallel, the inhibition of common-frequency and frequency-doubling vibration is realized, a phase compensation angle is introduced, and the stability of the system in the full rotating speed range is ensured by changing the phase angle in different frequency sections.
Compared with the prior art, the invention has the advantages that:
1) Most of the traditional magnetic suspension rotor vibration suppression methods are zero current control, most of vibration force is suppressed by suppressing harmonic current, but the suppression of harmonic current can only suppress current stiffness force, and the participated displacement stiffness force cannot be suppressed. The invention uses current and displacement to construct vibration force as input signal, and can realize the suppression of all vibration force.
2) The invention utilizes the orthogonal characteristic of the output signal of the displacement sensor, and uses one controller to simultaneously realize the vibration suppression in two directions, thereby reducing the calculation resource of the system and improving the dynamic response speed. And different phase angles are selected in different frequency sections through the phase compensation angle, so that the stability of the system in the full rotating speed range can be ensured.
Drawings
FIG. 1 is a flow chart of the present invention;
fig. 2 is a schematic structural diagram of a magnetic suspension rotor system, wherein 1 is an active magnetic bearing, 2 is a rotor, 3 is a geometric axis of the rotor, and 4 is an inertial axis of the rotor;
FIG. 3 is a block diagram of a basic control system of a magnetic levitation rotor;
FIG. 4 is a schematic diagram of a principle of synchronous selective coordinate transformation, wherein 5 is a sensor probe and 2 is a rotor;
FIG. 5 is a schematic block diagram of a synchronous selective coordinate transformation algorithm;
FIG. 6 is a block diagram of a multi-synchronous rotational coordinate transformation method and a master controller compound control system;
FIG. 7 is a simplified block diagram of an equivalent single-input synchronous rotation coordinate transformation method and a main controller composite control system.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, rather than all embodiments, and all other embodiments obtained by a person skilled in the art based on the embodiments of the present invention belong to the protection scope of the present invention without creative efforts.
According to the embodiment of the invention, as shown in fig. 1, the implementation process of the magnetic suspension rotor harmonic vibration force suppression method based on the multi-synchronous rotation coordinate transformation method is as follows: firstly, establishing a magnetic suspension rotor dynamic model containing mass unbalance and sensor harmonic; and then designing a controller based on a multi-synchronous rotation coordinate transformation method to suppress harmonic vibration force.
Step (1) establishing a magnetic suspension rotor dynamic model containing mass unbalance and sensor harmonic
The invention is applied to an active magnetic bearing system in a magnetic suspension control moment gyroscope. The structural schematic diagram is shown in fig. 2, wherein N is the geometric center of the stator, NXY is the inertial coordinate system, C and O are the mass center and the geometric center of the rotor, respectively, and O epsilon η is the rotational coordinate system. The method mainly aims at the suppression of the vibration force of two degrees of freedom of radial translation, so that only the modeling of two degrees of freedom of translation is considered. From newton's second law, the following kinetic equation is obtained:
where m is the rotor mass, X (t), Y (t) represent the translational displacements of the rotor mass center in the X and Y directions, respectively, f x (t),f y (t) resultant forces of bearing forces in X-direction and Y-direction, respectively, f ax (t),f bx (t),f ay (t),f by (t) is the bearing force of four pairs of radial magnetic bearings (two pairs at each end of A and B, the magnetic bearings are only in section, and only one pair is visible at each end of A and B, so only one pair is shown in the figure). When the rotor moves in a small displacement, the nonlinear bearing force can be approximately linearized and can be expressed as follows:
wherein K is i ,K h Current stiffness and displacement stiffness, i ax (t),i bx (t),i ay (t),i by (t) four pairs of radial magnetic bearing coil currents, x a (t),x b (t),y a (t),y b And (t) is the displacement of the bearing coordinate system.
As can be seen from fig. 2, since the mass imbalance causes the geometric center to be misaligned with the mass center, the displacement sensor measures the displacement of the geometric center of the rotor, which has the following relationship:
X(t)＝x(t)+Θ x (t)
Y(t)＝y(t)+Θ y (t)
wherein X (t), Y (t) represent the displacement of the geometric center, respectively, and can be obtained by the geometric relationship:
Θ x (t),Θ y (t) represents the amount of unbalance, and has the following form:
Θ x (t)＝ecos(Ωt+χ)
Θ y (t)＝esin(Ωt+χ)
e is the magnitude of the unbalance, χ is the initial phase, and Ω is the rotor rotational speed. It can be seen that the unbalanced mass can generate co-frequency interference signals, so that the active magnetic bearing generates co-frequency vibration force.
Because the displacement sensor has sensor harmonic waves, the output signals of the four pairs of displacement sensors are not the real displacement signals of the rotor, and the displacement sensors provide the following signals:
x as (t),x bs (t),y as (t),y bs (t) is the displacement sensor output signal; d axs (t),d bxs (t),d ays (t),d bys (t) is the sensor harmonic interference signal, which can be expressed as follows:
wherein s is ai ,s bi Is the harmonic amplitude of the sensor, χ i The initial phase is I, the harmonic order is I, and the harmonic of the sensor can generate an interference signal with the same frequency and the frequency multiplication with the rotating speed, so that the magnetic bearing system generates harmonic vibration force. The control diagram of a magnetic levitation rotor system with mass unbalance and sensor harmonics is shown in fig. 3, where G c (s),G w (s) respective tableShowing the controller transfer function and the power amplifier transfer function, K s Representing the sensor gain factor.
Step (2): magnetic suspension rotor harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation method
The controller takes the same-frequency vibration force and frequency multiplication current as input, is connected into an original closed-loop system in a parallel mode, and the output of the controller is fed back to the power amplifier input end of the original control system, and the design of the module mainly comprises the following two aspects:
(1) Multi-Synchronous rotational coordinate Transformation method (Multi-Synchronous Rotating Frame Transformation, MSRFT): according to different forms of common-frequency and frequency-doubled vibration forces generated by an actual magnetic suspension rotor system, common-frequency vibration comprises a current rigidity force and a displacement rigidity force, and frequency-doubled vibration only comprises the current rigidity force; according to a model of system electromagnetic force, current and displacement are used for constructing same-frequency vibration force to be used as input of a multiple synchronous rotation coordinate transformation same-frequency force suppression part, and frequency multiplication vibration can be suppressed by suppressing frequency multiplication current, so that the current is used as input of the frequency multiplication force suppression part;
(2) obtaining the stable condition of a system with multi-synchronous rotation coordinate transformation through theoretical analysis and proof; and designing a corresponding phase compensation angle according to the closed-loop characteristic of the actual magnetic suspension rotor system, and realizing the absolute stability of the system in the working rotating speed range through the phase compensation angle.
Further, the vibration force suppression algorithm in the step (2) is as follows:
the synchronous rotating coordinate transformation for realizing harmonic suppression mainly comprises three parts, namely, a disturbing signal is transformed from a static coordinate system to a rotating coordinate system through synchronous rotating coordinate transformation, at the moment, a same frequency/frequency multiplication disturbing signal is represented as a direct current quantity under the rotating coordinate system, then the direct current quantity is identified through low-pass filtering, and finally, the same frequency/frequency multiplication disturbing signal under the static coordinate system is obtained through synchronous rotating coordinate inverse transformation. The schematic diagram of the principle is shown in FIG. 4, wherein M is the geometric center of the rotor, C is the center of mass, epsilon is the vibration quantity, and CX is the rotation around the center of mass due to the self-centering principle when the rotor rotates at high speed s Y s As a stationary coordinate system, CX r Y r Is a rotating coordinate system, x s ,y s Is a vibration quantity coordinate, x, of a stationary coordinate system r ,y r The vibration amount coordinate under the rotating coordinate is a constant flow amount, and the harmonic component can be processed by low-pass filtering in the rotating coordinate system and synchronous rotating coordinate inverse transformation. The output is fed back to a closed loop structure formed by the output end of the power amplifier, so that the same frequency/frequency multiplication vibration force can be effectively inhibited.
1. Synchronous rotating coordinate transformation algorithm analysis
Because the phases of all frequency bands of the system are inconsistent, a phase compensation angle needs to be introduced for phase compensation in order to ensure absolute stability of the system within the full working frequency band range. As shown in fig. 5, the transformation equation of the synchronous rotating coordinate with the phase compensation angle is as follows:
wherein u is 1 (t),u 2 (t) is the input signal, u dc1 (t),u dc2 (t) is a synchronous rotation coordinate transformation output signal; omega is the rotational speed of the rotor and,
the phase compensation angle is used for ensuring the stability of a closed loop, when the harmonic frequency of an input signal is k times of omega, the input signal is output as a direct current signal after synchronous rotation coordinate transformation, the harmonic component in the original signal can be extracted through low-pass filtering, and a low-pass filter is provided with the following transfer functions:
where k is the low pass filter gain factor, ω c Is the cut-off frequency;
suppose the signal after low-pass filtering is
Order to
Respectively u as a laplace transform dc1 (s),u dc2 (s),
The following equation holds:
the above formula is rewritten as follows:
inverse laplace transform on the above equation can yield the following differential equation:
let κ = k ω c The above equation is written in matrix form as follows:
the low-pass filtered signal needs to be reversely transformed to the original stationary coordinate system by the synchronous rotating coordinate, and the transformation relation is expressed as follows:
x 1 (t),x 2 (t) is the output signal after the synchronous rotation coordinate inverse transformation;
from the above equations, the following state space expressions can be solved:
the transfer function matrix is derived from the state space expression:
the corresponding input and output relationship is as follows:
u 1 (t),u 2 (t) are the orthogonal sinusoidal signals, assuming the form:
A m and γ represents the amplitude and initial phase of the signal, respectively;
from the above formula, the following relationship exists:
iΩu 1 (s)＝su 2 (s)
an equivalent single input transfer function can be derived from the above equation:
let s = j ω, it can be found that when ω = i Ω, there are:
when ω is c And k, the synchronous rotation coordinate transformation can generate large gain at harmonic frequency, and the harmonic component can be effectively inhibited after the synchronous rotation coordinate transformation and the original system form a closed loop.
As shown in fig. 6, the co-frequency vibration force and the frequency doubling current are used as the input of the multi-synchronous rotation coordinate transformation, the output is fed back to the input end of the power amplifier and added with the control signal of the main controller, and the harmonic vibration force of the system is eliminated.
2. Stability analysis
Because the X direction and the Y direction have symmetrical structures and the same parameters, the stability analysis is carried out by taking the X direction as an example, the stability analysis idea of the invention is that the insertion of each controller is based on a first stabilization system, namely the insertion of the same-frequency vibration force algorithm is based on the stability of the original system, and the insertion of the frequency doubling suppression algorithm is based on the system stability of the adding same-frequency suppression algorithm. The invention only carries out expansion analysis on the stability condition of the same-frequency suppression algorithm insertion system, and the frequency multiplication algorithm is similar when being inserted, and the simplified block diagram is shown in figure 7.
The closed-loop characteristic polynomial of the system obtained from fig. 7 is:
1+2K s K i G c (s)G w (s)P(s)-2K h P(s)+2K i G w (s)G SRFx (s)＝0
the above equation can be converted into:
wherein,
is a system function.
When k =0, s = - ω c + -j omega, when a closed loop pole exists at-omega when k → 0 c Moving in the ± j Ω field, to ensure a sufficient stability margin of the system, using κ as an independent variable and s as a dependent variable, the following equation is obtained by differentiating κ =0 and s = j Ω:
to ensure that the closed loop eigenfunctions follow the s left half plane after the algorithm is added to the system, the following conditions need to be satisfied:
because of omega c Omega, so
I.e. the above formula is about
The following stability conditions can be obtained:
the phase compensation angle is properly adjusted according to the phases of the system function at different frequencies, so that the system can be stable in the full working rotating speed range.
Although the illustrative embodiments of the present invention have been described in order to facilitate those skilled in the art to understand the invention, it is to be understood that the invention is not limited in scope to the specific embodiments, but rather, it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and it is intended that all matter contained in the invention and created by the inventive concept be protected.
Claims (1)
1. A magnetic suspension rotor harmonic vibration force suppression method based on a multi-synchronous coordinate transformation method is characterized by comprising the following steps: the method comprises the following steps:
step (1): establishing a magnetic suspension rotor dynamics model containing mass unbalance and sensor harmonic
For an active magnetic bearing system in a magnetic suspension control moment gyroscope, setting N as a geometric center of a stator, NXY as an inertial coordinate system, C and O as a mass center and a geometric center of a rotor respectively, and O epsilon eta as a rotating coordinate system, and carrying out modeling of two translational degrees of freedom aiming at the suppression of the vibration force of the two translational degrees of freedom in radial direction; from newton's second law, the following kinetic equation is obtained:
where m is the rotor mass, X (t), Y (t) represent the translational displacements of the rotor mass center in the X and Y directions, respectively, f x (t),f y (t) resultant forces of bearing forces in X-direction and Y-direction, respectively, f ax (t),f bx (t),f ay (t),f by (t) is the bearing force of the four pairs of radial magnetic bearings, and the nonlinear bearing force is approximately linearized when the rotor performs a displacement motion less than a predetermined threshold as follows:
wherein K i ,K h Current stiffness and displacement stiffness, i ax (t),i bx (t),i ay (t),i by (t) four pairs of radial magnetic bearing coil currents, x a (t),x b (t),y a (t),y b (t) displacement under a bearing coordinate system; subscripts a, B denote the ends of a, B of the rotor system;
because the mass unbalance makes geometric center and center of mass misalignment, displacement sensor measures as rotor geometric center displacement, has following relation:
X(t)＝x(t)+Θ x (t)
Y(t)＝y(t)+Θ y (t)
wherein X (t) and Y (t) respectively represent the displacement of the geometric center in the X direction and the Y direction, and the displacement can be obtained by the geometric relationship:
Θ x (t),Θ y (t) represents the amount of unbalance, and has the following form;
Θ x (t)＝ecos(Ωt+χ)
Θ y (t)＝esin(Ωt+χ)
e is the amplitude of unbalance, chi is the initial phase, omega is the rotation speed of the rotor, and the mass unbalance can generate a same-frequency interference signal, so that the active magnetic bearing can generate a same-frequency vibration force;
because the displacement sensor has sensor harmonic waves, the output signals of the four pairs of displacement sensors are not the real displacement signals of the rotor, and the displacement sensors provide the following signals:
x as (t),x bs (t),y as (t),y bs (t) is the displacement sensor output signal; d axs (t),d bxs (t),d ays (t),d bys (t) is the sensor harmonic interference signal, expressed as follows:
wherein s is ai ,s bi Is the harmonic amplitude of the sensor, χ i The sensor harmonic can generate an interference signal with the same frequency and frequency multiplication with the rotating speed, so that a magnetic bearing system generates harmonic vibration force;
step (2): harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation
Designing a controller, wherein the controller takes the same-frequency vibration force and frequency multiplication current as input, is connected into an original closed-loop system in a parallel mode, and outputs and feeds back the output to the power amplifier input end of the original control system, and the controller comprises the following two aspects:
(1) the multi-synchronous rotation coordinate transformation method comprises the following steps: according to different forms of the same-frequency and frequency-doubled vibration force generated by an actual magnetic suspension rotor system, the same-frequency vibration comprises a current stiffness force and a displacement stiffness force, and the frequency-doubled vibration only comprises the current stiffness force; according to a model of system electromagnetic force, current and displacement are used for constructing same-frequency vibration force to be used as input of a multiple synchronous rotation coordinate transformation same-frequency force suppression part, and frequency multiplication vibration can be suppressed by suppressing frequency multiplication current, so that the current is used as input of the frequency multiplication force suppression part;
(2) calculating to obtain the stable condition of the system of multi-synchronous rotation coordinate transformation; and designing a corresponding phase compensation angle according to the closed-loop characteristic of the actual magnetic suspension rotor system, and realizing the absolute stability of the system in a working rotating speed range through the phase compensation angle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011380266.5A CN112432634B (en) | 2020-11-30 | 2020-11-30 | Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011380266.5A CN112432634B (en) | 2020-11-30 | 2020-11-30 | Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112432634A CN112432634A (en) | 2021-03-02 |
| CN112432634B true CN112432634B (en) | 2023-03-14 |
Family
ID=74698373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011380266.5A Active CN112432634B (en) | 2020-11-30 | 2020-11-30 | Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112432634B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113125073B (en) * | 2021-04-21 | 2022-10-14 | 北京控制工程研究所 | Online dynamic balance system and method for detection control rotor fixed magnetic suspension rotary joint |
| CN113485472B (en) * | 2021-08-04 | 2022-02-11 | 北京航空航天大学 | Magnetic suspension rotor same-frequency vibration torque suppression method based on double-channel wave trap |
| CN113670288B (en) * | 2021-08-24 | 2023-05-26 | 北京航空航天大学 | Magnetic suspension rotor harmonic vibration suppression method based on multi-rate quasi-resonance controller |
| CN114114919B (en) * | 2021-11-25 | 2022-08-02 | 北京航空航天大学 | Same-frequency vibration force inhibition method and system based on active magnetic bearing control system |
| CN114322971B (en) * | 2022-01-07 | 2023-09-12 | 北京航空航天大学 | Magnetic suspension rotor same-frequency vibration force suppression method based on biquad generalized integrator |
| CN114371622B (en) * | 2022-01-07 | 2024-04-12 | 北京航空航天大学 | Magnetic suspension rotor harmonic vibration force suppression method based on multi-harmonic inverse Park transformation |
| CN114326409B (en) * | 2022-01-07 | 2024-04-12 | 北京航空航天大学 | Magnetic suspension rotor direct vibration force suppression method based on double-channel harmonic reconstruction |
| CN114967446B (en) * | 2022-05-06 | 2024-09-17 | 中国人民解放军战略支援部队航天工程大学 | Magnetic suspension rotor multi-frequency vibration suppression method based on adaptive frequency estimation |
| CN116221276B (en) * | 2023-05-04 | 2023-08-15 | 山东华东风机有限公司 | Special control system and control method for magnetic suspension motor |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002006700A1 (en) * | 2000-07-13 | 2002-01-24 | Ebara Corporation | Vibration isolating device using magnetic levitating device |
| CN104503237A (en) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | Harmonic vibration control method for magnetic suspension rotor based on Fourier transform |
| CN105783898A (en) * | 2016-04-28 | 2016-07-20 | 北京航空航天大学 | Magnetic suspension rotor harmonic vibration suppression method based on frequency-domain self-adaption LMS (least mean square) algorithm |
| CN106773673A (en) * | 2016-11-23 | 2017-05-31 | 北京航空航天大学 | A kind of magnetic suspension rotor method for inhibiting harmonic current of the fractional compensation repetitive controller based on frequency self adaptation |
| CN106873655A (en) * | 2017-04-13 | 2017-06-20 | 北京航空航天大学 | A kind of magnetic bearing system multiple-harmonic vibration suppressing method based on finite-dimensional repetitive control |
| CN108490777A (en) * | 2018-03-19 | 2018-09-04 | 北京航空航天大学 | A kind of magnetic suspension rotor harmonic vibration power suppressing method based on improvement odd times Repetitive controller |
| CN109976403A (en) * | 2019-03-31 | 2019-07-05 | 南京航空航天大学 | The magnetic suspension rotor vibration harmonics suppressing method and system of more trapper parallel connections |
-
2020
- 2020-11-30 CN CN202011380266.5A patent/CN112432634B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002006700A1 (en) * | 2000-07-13 | 2002-01-24 | Ebara Corporation | Vibration isolating device using magnetic levitating device |
| CN104503237A (en) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | Harmonic vibration control method for magnetic suspension rotor based on Fourier transform |
| CN105783898A (en) * | 2016-04-28 | 2016-07-20 | 北京航空航天大学 | Magnetic suspension rotor harmonic vibration suppression method based on frequency-domain self-adaption LMS (least mean square) algorithm |
| CN106773673A (en) * | 2016-11-23 | 2017-05-31 | 北京航空航天大学 | A kind of magnetic suspension rotor method for inhibiting harmonic current of the fractional compensation repetitive controller based on frequency self adaptation |
| CN106873655A (en) * | 2017-04-13 | 2017-06-20 | 北京航空航天大学 | A kind of magnetic bearing system multiple-harmonic vibration suppressing method based on finite-dimensional repetitive control |
| CN108490777A (en) * | 2018-03-19 | 2018-09-04 | 北京航空航天大学 | A kind of magnetic suspension rotor harmonic vibration power suppressing method based on improvement odd times Repetitive controller |
| CN109976403A (en) * | 2019-03-31 | 2019-07-05 | 南京航空航天大学 | The magnetic suspension rotor vibration harmonics suppressing method and system of more trapper parallel connections |
Non-Patent Citations (4)
| Title |
|---|
| A two-stage synchronous vibration control for magnetically suspended rotor system in the full speed range;Cong，Peng等;《IEEE Trans. Ind. Electron》;20200131;第67卷(第1期);第480-489页 * |
| Synchronous vibration suppression of magnetically suspended rotor system using improved adaptive frequency estimation;Li，Jinlei等;《IEEE SENSORS JOURNAL》;20201001;第20卷(第19期);第11212-11220页 * |
| Vibration Suppression Control for AMB-Supported Motor Driveline System Using Synchronous Rotating Frame Transformation;Zheng,Shiqiang等;《IEEE Transactions on Industrial Electronics》;20150930;第62卷(第9期);第5700-5708页 * |
| Vibration Suppression Control for MSFW with Gyroscopic Effects Using Synchronous Rotate Frame;Cai,Kaiwen等;《2018 37th Chinese Control Conference (CCC)》;20180727;第3756-3760页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112432634A (en) | 2021-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112432634B (en) | Harmonic vibration force suppression method based on multi-synchronous rotation coordinate transformation | |
| CN108490777B (en) | Magnetic suspension rotor harmonic vibration force suppression method based on improved odd-order repetitive control | |
| CN106647843B (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current based on composite score repetitive controller | |
| CN114326409B (en) | Magnetic suspension rotor direct vibration force suppression method based on double-channel harmonic reconstruction | |
| CN109058292B (en) | A kind of novel magnetically levitated direct suppressing method of bearing unbalance vibration power | |
| CN108106611B (en) | Magnetic suspension rotor harmonic current suppression method based on multiple phase shift quasi-resonance control | |
| CN106444390B (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current based on FIR filter and fractional order repetitive controller | |
| CN108897218B (en) | Magnetic suspension rotor odd harmonic vibration suppression method based on mixed odd-order repetitive controller | |
| CN107870568A (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current based on second order bimodulus Repetitive controller | |
| CN106886152B (en) | Magnetic suspension rotor odd harmonic current suppression method based on second-order odd repetitive controller | |
| CN106610586B (en) | It is a kind of based on parallel FORC and delayed phase-lead compensation magnetic suspension rotor method for inhibiting harmonic current | |
| Fang et al. | Adaptive complete suppression of imbalance vibration in AMB systems using gain phase modifier | |
| CN106773673A (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current of the fractional compensation repetitive controller based on frequency self adaptation | |
| CN108227499A (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current of harmonic oscillator parallel connection fractional order repetitive controller | |
| CN111650975B (en) | Magnetic suspension rotor harmonic current suppression method based on multi-order repetitive controller | |
| Cui et al. | Harmonic vibration control of MSCMG based on multisynchronous rotating frame transformation | |
| CN112525180B (en) | Harmonic vibration torque suppression method based on hybrid repetitive controller | |
| CN106896726A (en) | A kind of magnetic suspension rotor method for inhibiting harmonic current of the digitial controller based on parallel second order Repetitive controller DPSORC | |
| Chen et al. | Internal model control for the AMB high‐speed flywheel rotor system based on modal separation and inverse system method | |
| CN114322971B (en) | Magnetic suspension rotor same-frequency vibration force suppression method based on biquad generalized integrator | |
| CN114371622B (en) | Magnetic suspension rotor harmonic vibration force suppression method based on multi-harmonic inverse Park transformation | |
| CN113670288B (en) | Magnetic suspension rotor harmonic vibration suppression method based on multi-rate quasi-resonance controller | |
| CN112525179B (en) | Harmonic current suppression method based on virtual sampling second-order repetitive control | |
| Zhao et al. | Feedforward decoupling control for rigid rotor system of active magnetically suspended high‐speed motors | |
| Ren et al. | Rotation modes stability analysis and phase compensation for magnetically suspended flywheel systems with cross feedback controller and time delay |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |