# Active Magnetic Bearings Stiffness and Damping Identification from Frequency Characteristics of Control System

Authors: Chaowu Jin; Yuanping Xu; Jin Zhou; Changli Cheng
Journal: Shock and Vibration
 (2016)

---

## Abstract

At present, the stiffness and damping identification for active magnetic bearings (AMBs) are still in the stage of theoretical analysis. The theoretical analysis indicates that if the mechanical structure and system parameters are determined, AMBs stiffness and damping are only related to frequency characteristic of control system, ignoring operating condition. More importantly, few verification methods are proposed. Considering the shortcomings of the theoretical identification, this paper obtains these coefficients from the experiment by using the magnetic bearing as a sine exciter. The identification results show that AMBs stiffness and damping have a great relationship with the control system and rotating speed. Specifically, at low rotating speed, the stiffness and damping can be obtained from the rotor static suspension by adding the same excitation frequency. However, at high speed, different from the static suspension situation, the AMBs supporting coefficients are not only related to the frequency characteristics of control system, but also related to the system operating conditions.


---

## Body

## 1. Introduction

From the dynamic point, the AMBs rotor system can operate stably since the AMBs can provide appropriate supporting stiffness and damping parameters (bearing parameters) for the rotor. Now, the researches on AMB supporting parameters identification are mainly from two aspects.One idea is from the mechanical bearing parameters identification. Through the relationship between the excitation and rotor bearing system response, the bearing parameters can be identified through the identification algorithm along with experiment. Based on the unbalance response of rotor bearing system, Zhou and Ni [1] identified the AMBs stiffness and damping of rigid rotor bearing model and verified identification by finite element simulation. Zhang et al. [2] set the modal parameters of AMBs rotor system as the target response, using the complex modal model updating technology, and identified AMBs stiffness and damping under static suspension. De Santiago and San Andrés [3, 4] identified these coefficients for mechanical bearing from transient rotor dynamic response due to impacts and imbalances. Bi et al. [5] proposed a time-frequency-domain hybrid method of identifying the dynamic characteristic parameters of journal bearings and the unbalance of the rotor on field condition and without artificial excitations. Wang and Xia [6] gave a method and data processing program for extracting the stiffness parameters of the rolling element bearings from random vibration responses. Chen et al. [7] identified the dynamic stiffness and damping of bearing by the boundary element method. Tiwari [8, 9] proposed identification algorithms for the estimation of mechanical bearing parameters and residual unbalances through forced responses of the system. Tiwari and Chakravarthy [10] developed an identification algorithm for the simultaneous identification of the bearing and unbalance parameters for flexible rotor bearing systems by using run-down measurements. The disadvantage of these identification methods for mechanical bearing is that they need to apply excitation to the rotor and an accurate rotor model is indispensable. Test data must be under the rotating situation, which also puts forward higher requirements on excitation forms, system response, and stability of the system.Identifying the stiffness and damping from the theoretical analysis of AMB control system by obtaining the system transfer function is another identification idea. The AMBs system uses active control strategy; therefore, the supporting parameters are closely related to the system structure and control strategy. Wang [11] studied the effects of AMB control system frequency response characteristics and structural parameters on the AMBs stiffness and damping and obtains the frequency characteristics of the equivalent stiffness and damping. Zhao et al. [12] studied the AMB rotor system with analog PID controller based on the independent control principle; it showed that the AMBs stiffness and damping were influenced by the differential and proportional coefficient in controller. Humphris et al. [13] obtain the stiffness and damping values for a magnetic bearing at static situation and found that these coefficients of AMBs are affected by the position and velocity feedback parameters of controller. By analyzing the frequency dependent feedback controller transfer function, Williams et al. [14] obtained the equivalent stiffness and damping coefficients theoretically. Under the shaker excitation, Lim et al. [15] experimentally identified a five-degree-of-freedom hybrid magnetic bearing for blood bump applications at static suspension situation with proportional-integral-derivative (PID) controllers. Tsai et al. [16] applied the wavelet transform algorithm to identify the magnetic damping and stiffness coefficients. This identification method had to simplify the control system in order to get the transfer function. But, for some system using the complex control strategy such as adaptive control method, it was very difficult to get the transfer function accurately.Aiming at the shortages of the researches on AMBs stiffness and damping, this paper uses magnetic bearing as the exciter, tests the frequency characteristics of the control system, obtains the relationship between bearing parameters and frequency characteristics of the control system, and verifies the identified stiffness and damping values. The results show that, at low speed, the bearing parameters of AMBs are mainly related to the controller frequency characteristics, but at high speed, the bearing parameters are not only related to the controller frequency characteristics, but also related to the speed.The remainder of the paper is organized as follows. Section2 describes the theoretical basis of the relationship between bearing parameters and frequency characteristics of the control system. The test results of the frequency characteristics of the control system are presented in Section 3 using magnetic bearing as the exciter. Section 4 shows experimental verification in order to validate the proposed method. Conclusions are drawn in Section 5.


## 2. Fundamental

The structure of radial AMB with single degree of freedom used in this paper is shown in Figure1. Table 1 shows its specific parameters.Table 1
The related parameters of radial AMB.
Parameter
Value
Coil number of signal-poleN
120
Magnetic pole areaA/m2
2 × 10−4
Bias currentI0/A
1
Air gapx0/mm
0.3
Permeability of vacuumμ0
4
π
×
10
-
7
Pole angleβ/°
22.5Figure 1
The structure principle diagram of radial AMB with single degree of freedom.The electromagnetic force generated by a pair of differential magnetic poles in the figure can be expressed as(1)Ft=μ0AN24cos⁡βI0-it2x0-xt2-I0+it2x0+xt2.As shown in (1), the electromagnetic force of differential arrangement AMB is nonlinear related to the displacement x(t) and the control current i(t). The nonlinearity of the electromagnetic force is detrimental for the application of linear control theory, so linear analysis of the formula is necessary, which can be expressed as(2)Ft=-kx0·xt+ki0·it,where kx0, ki0 are displacement stiffness and current stiffness of radial AMBs at the equilibrium position, respectively, which can be written as follows:(3)kx0=μ0AN2I02cos⁡βx03,ki0=-μ0AN2I0cos⁡βx02.Because of the complete symmetry of the AMB radial structure, one degree of freedom of the radial AMB system can be the research object; control block diagram of a one degree-of-freedom system is shown in Figure2.Figure 2
Control block diagram of one degree-of-freedom radial AMB system.I
0 and i(t) are the bias current and the control current into the electromagnet coil, respectively; vs is the voltage signal converted from the rotor displacement by the displacement sensor, which has a linear relationship with the rotor vibration displacement x(t); vc is the input of controller; cr is the gap between the catcher bearing and rotor, named protective gap, which is usually half of air gap x0; v0 is reference voltage of controller, which is to multiply displacement gain by protective gap; Ga, Gs, and Gc are the transfer function of the power amplifier, the displacement sensor, and the PID controller, respectively.According to the block diagram of control system, the relationship between control current and rotor vibration displacement is obtained:(4)it=vc·Gc·Ga=v0-vs·Gc·Ga=v0-cr+xt·Gs·Gc·Ga=-xt·Gs·Gc·Gav0=cr·Gs.Take (4) into (2): (5)Ft=-kx0·xt-ki0·Gs·Gc·Ga·xt.Laplace transform (5):(6)Fs=-kx0-ki0·Gss·Gcs·Gas·xs.AMBs generalized stiffness is the ratio of electromagnetic force Laplace transformation and rotor displacement Laplace transformation; that is,(7)Ks=Fsxs=-kx0-ki0·Gss·Gcs·Gas.One degree-of-freedom bearing force with the spring damping characteristics can be expressed as(8)ft=k·xt+c·x˙t.Transform (8) by Laplace transform and get the generalized stiffness:(9)K′s=fsxs=k+c·s.Considering the frequency characteristics of the generalized stiffness, takes=jω, so (7) and (9) are (10)Kjω=-kx0-ki0·Gsjω·Gcjω·Gajω,K′jω=k+c·jω.From (10), the magnetic bearing equivalent stiffness and damping can be written as(11)k=-kx0-ki0·Re⁡Gsjω·Gcjω·Gajω,c=-ki0·Im⁡Gsjω·Gcjω·Gajωω,where Re[Gs(jω)·Gc(jω)·Ga(jω)] and Im[Gs(jω)·Gc(jω)·Ga(jω)] in the equation represent the real and imaginary parts of Gs(jω)·Gc(jω)·Ga(jω), respectively, where ω is the external excitation frequency, rather than the speed of rotor. Seen by (11) in certain control parameters, and ignoring the influence of eddying effect, the equivalent stiffness and damping of the radial AMB are only related to the external excitation frequency ω, which has nothing to do with the rotor speed Ω unless the excitation is due to unbalance rotor, in which the excitation frequency and rotor speed are clearly related.By formula transformation, (4) can be expressed as(12)Gsjω·Gcjω·Gajω=-ijωxjω=-iAxAeφi-φx,where xA and iA are the amplitude of the rotor displacement and control current, respectively, and φx and φi are the phase of displacement and control current, respectively. If the amplitude and phase of control current i(t) and rotor displacement x(t) can be obtained by testing method, the frequency characteristics of Gs(jω)·Gc(jω)·Ga(jω) and AMBs equivalent stiffness and damping can also be obtained. According to the characteristics of AMB system, it can test the frequency characteristics of the control system loop using AMB as an incentive. Schematic of the excitation applied way is shown in Figure 3.Figure 3
Schematic of the excitation applied way.Using signal generator to generate a signal with frequencyω and amplitude 1 V, put this signal to the input of the power amplifier 1 corresponding to the 1+ magnetic pole. At this time, the coil current of 1+ magnetic pole includes bias current I0, feedback control current i(t), and external excitation current ie(t). The rotor displacement signal x(t) and the control current i(t) of 1− pole coil will be collected, respectively. The data will be substituted into (6) to get the amplitude ratio and phase difference. Change the excitation frequency and repeat the testing process, the frequency characteristic curve of Gs(jω)·Gc(jω)·Ga(jω) about ω can be drawn, and the equivalent stiffness and damping of AMBs can be obtained with the combination of (11).


## 3. Experimental Investigations

According to the structural parameters in Table1, displacement stiffness kx and current stiffness ki at different rotor displacement x and control current i can be obtained by using the finite element analysis method, as shown in Figure 4.Figure 4
Current stiffness (a) and displacement stiffness (b) of finite element analysis.
(a)
(b)Current stiffnesski0 and displacement stiffness kx0 of the rotor at the equilibrium position can be obtained by taking the structure parameters of Table 1 into (3). The value at x=0 and i=0 in Figure 4 is the current stiffness ki0 and displacement stiffness kx0 of the rotor at the equilibrium position calculated by finite element method. List kx0 and ki0 at the equilibrium position in Table 2, which are obtained by finite element method and theoretical calculation.Table 2
k
x
0 and ki0 at the equilibrium position of radial AMB.
Name
Values of theoretical calculation
Values of finite element analysis
Displacement stiffnesskx0/N/m
4.639
×
10
5
4.507
×
10
5
Current stiffnesski0/N/A
-
1.5063
×
10
2
-
1.4906
×
10
2From Table2, we can see that values of finite element analysis and theoretical calculation are very close, but the displacement stiffness and current stiffness calculated by finite element method are less than the theoretical calculation values. The reason is that, in the theoretical calculation, the magnetic reluctance, magnetic leakage, magnetic saturation, and other factors of silicon steel sheets are ignored. In order to get accuracy stiffness and damping, we adopt kx0 and ki0 calculated from finite element analysis in the following.When the external excitation frequency is less than the cutoff frequency of displacement sensor and power amplifier, it can be seen as a proportional component. In this paper, the transfer function of displacement sensor can be written asGs=2.5/cr, the transfer function of power amplifier can be written as Ga=2I0/10, and the transfer function for PID controller can be written as(13)Gc=Kp+Kis+Kds1+τds,where Kp, Ki, and Kd are the proportional coefficient, the integral coefficient, and the differential coefficient, respectively, and τd is the filtering time constant. The parameters of the controller are shown in Table 3; the frequency characteristic of the controller can be calculated by taking the controller parameters into (13).Table 3
Parameters of the controller by testing.
Direction
K
p
K
i
K
d
τ
d/s
1#
7.15
32.93
2.7
e
-
3
5.69
e
-
6
2#
6.85
32.06
2.6
e
-
3
7.79
e
-
6
3#
8.19
39.69
3.2
e
-
3
6.42
e
-
6
4#
8.21
39.12
3.1
e
-
3
6.33
e
-
6Figure5 shows the studied AMBs system structure, the left radial AMBs are 1# and 2#, and the right radial AMBs are 3# and 4#. The experimental test rig is pictured in Figure 6.Figure 5
Structure of AMB system.Figure 6
AMB unbalance excitation test system.Theoretical values ofGs(jω)·Gc(jω)·Ga(jω) can be obtained by multiplying the transfer function of displacement sensor, power amplifier, and controller and analyzing frequency characteristics. In order to acquire the experimental data, signal generator is employed to drive the external current and NI data acquisition card is employed to record the current and displacement values according to the principle and method of Figure 3. Then use one-order Fourier series based on the least square method to fit the filtered signals to get the amplitude and phase data. The initial frequency is 20 Hz and the stop frequency is 400 Hz, tested once every 5 Hz.Experimental and theoretical values of the control system frequency characteristics are shown in Figure7, from which it can be observed that theoretical and experimental values trends of the frequency characteristics in all radial bearings of the control system are in good agreement but not perfectly matched with the change of frequency. This may be due to the error of system model and experiment measurement. However, the error at low frequency is relatively small, and the error increases with the increase of frequency. The reason maybe is that the displacement sensor and power amplifier cannot always be considered as a proportional part in high frequency.Figure 7
Experimental values and theoretical values of each frequency characteristic.The equivalent stiffness and damping of each radial AMB can be obtained by substituting testing frequency characteristics of each AMB control system into (11). Figures 8 and 9 are, respectively, the equivalent AMBs stiffness and damping with the change of vibration frequency in 1#, 2#, 3#, and 4#.Figure 8
The stiffness and damping in 1#, 2#.Figure 9
The stiffness and damping in 3#, 4#.According to Figures7 and 8, in the range of the exciting frequency, the equivalent stiffness and damping in 1# are less than those in 2#; the equivalent stiffness and damping in 3# are less than those in 4#. And with the change of frequency, the equivalent stiffness of each bearing increases, while the equivalent damping of each bearing increases firstly and then decreases slowly, which are broadly consistent with the stiffness and damping of theory analysis in the trend of change.


## 4. Experimental Verification

Experimental verification is aimed to test the accuracy of the bearing parameter identified in this paper. Equation (11) shows that AMB parameters are only related to the excitation frequency ω, which have nothing to do with the rotor speed Ω when the system structure and control parameters are determined. When the rotor is in unbalanced excitation, the excitation frequency is equal to the rotor speed, ω=Ω. It means that the stiffness and damping with the rotor speed of Ω are the same as the stiffness and damping coefficient with the rotor in static suspension but the excitation frequency ω=Ω.In order to verify the accuracy of the bearing parameters identified in this paper. Firstly, we add unbalance mass on the rotor as the excitation and the rotor unbalance response measurements are conducted for rotor speeds from 20 Hz to 300 Hz. Secondly, the AMB rotor finite element model is built in Patran (a commercial software for Finite Element Analysis). Previous work [17] by us has already detailed the rotor modeling, model updating, and model verification of the updated model. By providing the stiffness and damping values to the bush cell and adding the same unbalance excitation for finite element model, the rotor unbalance response is simulated to obtain the unbalance displacement amplitude. Finally, we will compare the simulation amplitudes and testing unbalance response amplitudes. Flow chart of verification is shown in Figure 10.Figure 10
Testing flow chart.Before the unbalance response experiment, the rotor has been adjusted to nearly balanced condition under the dynamic balancing machine to avoid the residual unbalance influence from the rotor itself. When the rotational frequency is greater than 300 Hz, the rotor displacement response amplitude is bigger under the unbalanced excitation. In view of the stability and safety of system, the rotational frequency in this paper is below 300 Hz. Figure11 is the comparison charts of each finite element simulation amplitude and testing unbalance response amplitudes.Figure 11
Simulation and measured rotor unbalance responses.As can be seen from Figure11, when the rotational frequency is below 200 Hz, the finite element model of rotor bearing can be established based on the testing stiffness and damping in the rotor static suspension. The unbalance response amplitudes obtained by the simulation of rotor operation are basically consistent with the testing unbalance response amplitudes. This shows that the impact of speed on the radial AMBs parameters is very small, which means that the bearing parameters are mainly related to the frequency characteristics of the control system. That is, when the rotor speed is Ω, bearing parameters are basically consistent with the bearing parameters corresponding to the static suspension and the excitation frequency ω=Ω.However, when the rotational frequency is above 200 Hz, there is large difference between the finite element simulation amplitudes and the testing amplitudes in 1# and 2#. And the finite element simulation amplitudes and the testing amplitudes in 3# and 4# also have certain difference. Parameters are not only related to the frequency characteristics of the control system but also related to the speed factor. This is mainly because the eddying effect in high speed influences the bearing characteristics of system, which leads to the difference of the amplitudes of two kinds of model. The difference between the finite element model and the actual rotor model and the difference between the unbalanced incentive of finite element model and the unbalanced mass of actual test are also the cause of the difference of the amplitudes of two kinds of model.


## 5. Conclusions

(1)
This paper studies the relationship between the AMBs bearing parameters and the frequency characteristics of control system using the experimental method. By employing the magnetic bearings as a magnetic exciter, the AMBs stiffness and damping coefficients are obtained under the frequency domain, which is a simple and convenient method.
(2)
This paper verifies the rationality of the testing AMBs parameters through the comparison of the finite element simulation amplitudes and the testing amplitudes using the unbalance excitation.
(3)
The experimental results show that at low speed the bearing parameters are mainly determined by the controller characteristics. While at high speed, the bearing parameters are not only related to the control rule but also related to the speed. This may be due to the influence of eddying effect.



---

*Source: 1067506-2016-01-10.xml*