Dynamic Stiffness Analysis and Measurement of Radial Active Magnetic Bearing in Magnetically Suspended Molecular Pump

Current stiffness and displacement stiffness are two important parameters of radial active magnetic bearing (RAMB) that are generally considered as constants in a control system. However, such presumption may lead to the probable degradation of the control performance of the RAMB system when the perspective that the current and displacement stiffness should be variable due to variations in the speed. On this regard, a structure analysis and stiffness measurement method based on non-variable stiffness would not be feasible for RAMB under high-speed conditions. This paper presented an analysis of the dynamic stiffness characteristics of RAMB by means of a dynamic equivalent magnetic circuit method and laid out a comparison between the results of the theoretical analysis and the simulation results of finite element method. In particular, a novel dynamic stiffness measurement method of RAMB under high-rotation frequency was introduced. Results of experiments on the dynamic stiffness of RAMB demonstrated an excellent agreement with the theoretical research and the finite element analysis results, thereby verifying the rationality of the discussed dynamic stiffness characteristics. Practically, the proposed measurement method of RAMB dynamic stiffness provides an accurate analysis for the dynamic stiffness and contributes inspiring research significance for the dynamic properties of RAMB.

gimbal magnetically suspended control moment gyroscope. Moreover, literature 21 proposed a new stiffness measurement method for the magnetically suspended flywheel for determining the current stiffness and displacement stiffness of permanent magnet biased radial magnetic bearings. Further, a detailed stiffness measurement method for a radial hybrid magnetic bearing was described in 22 . In all these studies, the stiffness of magnetic bearings has been regarded as a constant (not affected by rotation frequency) and the measurement experiments were performed when the rotor was suspended steadily at zero-frequency. For the low speed applications in [20][21][22] , the stiffness at zero speed can be used as a constant in the whole speed range. However, some researchers have found that the stiffness of radial active magnetic bearing (RAMB) would decline in a discernible manner with an increase in rotor frequency and would influence the stability of the radial magnetic bearing system in the process 23,24 . The literature presented herein also depict that the frequency response of the current stiffness could be obtained by addition of a balancing mass; nonetheless, the method is complicated and no explicit description of measuring current stiffness is provided in this research. This aside, the frequency response of the displacement stiffness is not analysed nor measured. Therefore, it is necessary implement the meterage of dynamic stiffness of RAMB for the magnetically suspended molecular pump (MSMP).
In this paper, an accurate dynamic stiffness measurement method of RAMB is put forward for MSMP. A build-up of the characteristics of the RAMB dynamic stiffness is first established by means of an equivalent magnetic circuit with varying frequency. Diagrams of magnetic line distribution and change trend of the stiffness is then obtained using finite element method (FEM). Finally, the proposed dynamic stiffness measurement method of RAMB is verified by prototyped experiments via comparison of results of the theoretical and FEM analyses.

Results
Analysis of dynamic stiffness characteristics of RAMB for MSMP. Figure 1 shows the magnetically suspended molecular pump model, mainly composed of a single rotor shaft, two RAMBs, an axial active magnetic bearing, a high-speed brushless direct current motor, two integrated displacement sensors and stator blades. For the high pumping-speed MSMP CXF400-4000, designed by professors and doctors in our academy and applied in civilian equipment, the configuration of the RAMB is pictured in Fig. 2. In the figure, the materials of the stator and rotor core are silicon steels with corresponding dimensions as described in Table 1.
While the rotor operates at or near the rated speed, eddy currents caused by time-varying control currents will affect the resistance of the rotor and the stator, which is an important factor influencing the stiffness 23 . Such phenomenon results in substantial reduction in stiffness compared to the 0 Hz stiffness value that is considered  constant in the RAMB system. As such, the traditional magnetic circuit model, which ignores the rotor frequency variation and defaults the control current and magnetic resistances to a constant value, is not suitable for the RAMB stiffness analysis. Through an analysis of the dynamic characteristics when the rotor speed is variable while other parameters remain unchanged, studies in 23,24 have proved that rotor frequency variation and eddy current effect leave a marked impact on magnetic bearing stiffness. Accordingly, these studies established an analytical method for a hybrid radial magnetic bearing system based on an equivalent circuit model whose parameters were frequency-dependent. Although it has been generally recognised that eddy current effects could be mitigated by replacing the solid structure with laminated cores and by designing modifications, the magnetic bearings stiffness was still considerably influenced by eddy currents when the rotation frequency was sufficiently high [25][26][27] . Therefore, it is necessary and is of great significance to carry out a quantitative analysis on the dynamic stiffness of RAMB at different rotation speeds.
Specifically, the rotor speed change will lead to a time-varying control current, which will induce the eddy current effect in the stator core and rotor core, leading to the modification of their magnetic permeability. Based on the dynamic permeability of the stator and rotor, an analysis of accurate quantitative dynamic stiffness characteristics for the varying frequencies can be performed by establishment of an equivalent magnetic circuit.
The relative permeability of stator and rotor with different rotation frequencies has been analysed in 25 . In particular, the relative permeability of stator (μ rs ) could be deduced using Eq. (1) while the relative permeability for rotor could be determined by Fourier transform and the n-th relative permeability (μ nr ) can be expressed through Eq. (2). Here, k n is a coefficient determined by the stator structure; the other parameters are assigned in Table 1 Figure 3 shows the dynamic equivalent magnetic circuit model of RAMB, in which R s1 and R s2 are the equivalent magnetic reluctances of the stator; R x1+ , R x2+ , R x1− , R x2− , R y1+ , R y2+ , R y1− and R y2− are utilised to depict the magnetic reluctances of the air gap in x 1+ , x 2+ , x 1− , x 2− , y 1+ , y 2+ , y 1− and y 2− direction, respectively and; R rn is the dynamic reluctances of the n-th harmonic of the rotor. According to the distribution of magnetic flux density, whereas the rotor is in the central position, the reluctances can be calculated based on Eqs. (1) and (2) as follows: x x x x y y y y a p Model of dynamic current and displacement stiffness. Dynamic current stiffness. Neglecting the coupling between two directions x and y, as the rotor rotates at the balance location, the magnetic fluxes in the x+ and x− directions could be denoted as φ x+ and φ kx− , as depicted based on Ohm's law of magnetic resistance as x tx where ρ is the magnetic flux leakage coefficient; I is the bias current; i is the control current and; R tx+ and R tx− are respectively, the total reluctance of the magnetic circuit in the x+ and x− directions. According to the formula for calculating electromagnetic force, the magnitude of the electromagnetic force in the x+ and x− directions could be respectively expressed as: x x l l a p tx 2 0 cos( ) In the differential mode, the force in the x-direction of the rotor represents the difference between the electromagnetic forces in the x+ and x− directions. The magnitude is denoted as f x and the calculation equation is , i 2 can be ignored due to the fact that the control current i is a tiny value compared to the bias current I. Here, the dynamic current stiffness k i can be defined as   Table 1. Note that the value of current stiffness declined with an increase in rotor frequency. Compared at 0 Hz, the current stiffness at 500 Hz could be seen to have decreased by 4.9%. Such distinct changes could be expected to inevitably affect the dynamic performance of RAMB.
Dynamic displacement stiffness. Let d x denote the displacement of the rotor in the x-direction that is much smaller than the length of air gap g. The expressions of air gap reluctance in x+ and x− directions can be expressed as gx gx a p gx gx a p In a similar manner as Eqs. (7) through (12), the RAMB electromagnetic force can be calculated with Eqs. (15) through (17); consequently, the dynamic displacement stiffness when the higher order of d x is ignored can be obtained from Eq. (18): x x a p a p tx x x x a p tx tx www.nature.com/scientificreports www.nature.com/scientificreports/  Figure 5 illustrates the relationship between RAMB dynamic displacement stiffness and frequency, which demonstrated, in particular, that the variation tendencies of displacement stiffness and current stiffness are consistent with an increase in frequency.

Influences of control parameters and structural parameters on stiffness.
For the RAMB structure shown in Fig. 2, the equation of electromagnetic force in RAMB can be linearised at the magnetic centre as 22 s i with f as a representation of the linear function of current and displacement and k s and k i being the respective displacement stiffness and current stiffness of magnetic bearings that are closely related to magnetic circuit characteristics and structural parameters. On such basis, the stiffness of the RAMB system can be obtained using Eq. (20), which defines it as the ratio of the electromagnetic force f and the rotor displacement that is dependent on current stiffness, displacement stiffness and the current. Furthermore, the current is a function of rotor displacement, which is dependent on the controller design of the RAMB system. Accordingly, the control and structure parameters have great influences on the system stiffness and need to be analysed.
i s Stiffness analysis with different control parameters. When the rotor is acted upon by an external force, the result of force analysis in the x-direction can be obtained with Eq. (21). Essentially, Eq. (21) employs the Laplace transform, whereas Eq. (22) defines the transfer function with the external force as input and displacement as output: with G(s) as the total transfer function of the displacement sensor, PID controller and power amplifier. According to the theory of vibration, the equation of motion for second-order system can be expressed by Eq. (23), where D represents damping and K the stiffness of the RAMB system. Applying Eqs. (22) and (23), s = 2πfj could be set and an equal stiffness of the RAMB system can be obtained through Eq. (24), with φ(2πf) for the phase angle of the transfer function G(s), The substantial influence of the control parameters on the equal stiffness for RAMB system has been clarified; moreover, the incomplete differential PID controller 28 , whose transfer function is presented as Eq. (25), is generally utilised to overcome high-frequency noises: where K p , K i , T d and K d are the typical parameters of the incomplete differential PID controller, and K i is finally designed according to the general actual effect for its decisive role in controlling the steady-state error. K p , K d and T d , thus, become the main research objects in this part, whereas the other parameters such as the sample AD converter and current controller amplification coefficients are kept stable, as listed in Table 2.
Equal system stiffness with different K d . Figure 6 shows the stiffness characteristics of the close loop system for RAMB with different K d values. Apparently, there was almost no difference in the equal stiffness of the system at the low frequency bands. However, in the medium-and high-frequency areas, the equivalent stiffness not only rose with higher frequency, but the magnitudes of the growth increased with K d as well. Take for instance the case where K d was set to 0.004, 0.006 and 0.008. Here, the increased equal stiffness from 400 to 1000 Hz were 3.1 × 10 6 , 4.8 × 10 6 and 6.3 × 10 6 N/m, respectively.    www.nature.com/scientificreports www.nature.com/scientificreports/ Equal system stiffness with different K p . Figure 7 displays equal stiffness for the close loop system of RAMB with different K p values. In this system, the equal stiffness would remain stable and the equivalent stiffness would be proportional to K p at the low-frequency band. However, the equivalent stiffness increased with frequency in the medium-and high-frequency areas, and their differences as a result of different K p values were invariable. For instance, the respective differences for K p from 1 and between 5 and 10 were 1.2 × 10 6 and 1.6 × 10 6 N/m.
Equal system stiffness with different T d . Figure 8 shows the stiffness characteristics of the RAMB system for different values of T d . The almost coincided and horizontal curves in the low-frequency band indicated that the system stiffness was not affected by the frequency and T d value. However, in the medium-and high-frequency areas, the equivalent stiffness became higher with an increase in frequency, although the magnitude of such growth increased with a decline in T d . For example, setting T d to 0.0002, 0.0003 and 0.0005 would increase equal stiffness from 400 to 1000 Hz at 5.1 × 10 6 , 3.3 × 10 6 and 1.2 × 10 6 N/m, respectively.
Stiffness analysis with different structural parameters. Equations (12) and (18) clearly depict that the values of current and displacement stiffness are influenced by many structural coefficients such as air gap length, pole width of the stator and coil turns. On this regard, this study carried out an analysis of the influence of structural parameters on dynamic stiffness.
Stiffness with air gap length. The variation of air gap length leaves a considerable influence on the magnetic reluctance of the air gap, leading to a corresponding variation in current and displacement stiffness of magnetic bearings. Herein, several air gap values were selected from 0 to 1 mm, for an investigation of the change of current and displacement stiffness at 1000 Hz. Figure 9 presents the results. The figure validated that current and displacement stiffness increase with air gap. Moreover, the rate of stiffness growth slowed down until zero as the air gap increased. Considering such decrease of magnetic force being directly caused by a large air gap, value of the air gap should not be too large; the 0.2-0.4-mm range would be perfect for the RAMB discussed in this paper.
Stiffness with pole width of stator. From the dynamic stiffness characteristics obtained in (12) and (18), it can be seen that the pole width of stator is also the key factor affecting the dynamic stiffness of magnetic bearing. Adopting the same analytical method as air gap length and the results are shown in Fig. 10. The analytical results  Stiffness with pole coil turns. As an important parameter of the proposed dynamic stiffness model shown in Eqs. (12) and (18), the value of coil turns has substantial impact on current and displacement stiffness. To study the change of dynamic stiffness with varying coil turns, several coil turn values were selected from 0 to 1000. Results of the analysis are presented in Fig. 11. In the figure, the exponentially rising curves signify that the number of turns of the coil should be set to a large value. However, it was necessary to consider that the increase of coil turns   www.nature.com/scientificreports www.nature.com/scientificreports/ would lead to large inductance, which reduces the stability of the system in the actual RAMB design. Thus, a value from 200 to 400 is ideal for the RAMB discussed in this paper.
Finite element analysis of dynamic stiffness. By the electromagnetic field analysis model in ANSYS software, the 2D-FEM mesh model of RAMB introduced in Fig. 2 was obtained as Fig. 12. For this model, the plane53 magnetic structure unit was used for its perfect performance of transient magnetic field analysis. Additionally, the rotor part was meshed by the mapped grid mode while the remaining part was meshed by the free mode. In particular, areas of air gap between the rotor surface and the fixed rotor where the electromagnetic field varies on a great scale, were finely meshed. In the FEM simulations, the bias current was set at 0.38 A, which is expected to generate a corresponding bias magnetic flux density, while the different rotor frequency values from 0 to 800 Hz were set for the rotor in the ANSYS software. Subsequently, a FEM simulation was carried out for each speed while the distributional patterns of the magnetic line of the magnetic bearing presented in Fig. 13 was obtained via the finite element calculation.
Based on Fig. 13(a), the magnetic lines of force were uniformly distributed when rotor was at a standstill. However, when the rotor rotated, the generated alternating magnetic field induced an alternating electric field that generated eddy currents and additional magnetic fields, thereby weakening the main magnetic field. As the rotor frequency increased, the theoretical analysis in 26 indicated that the rotor magnetic induction line would no longer be evenly distributed but would tend to the rotor surface due to the skin effect caused by the eddy currents. For the RAMB structure and FEM simulation conditions described in this paper, the FEM simulation results were consistent with the theoretical analysis as depicted in Fig. 13(b), in which the flux lines distorted severely and were squeezed to the surface of the rotor. This indicates that the magnetic forces would decrease rapidly and that such phenomenon could be hardly observed in AMB rotor experiments at a rotation speed of 800 Hz. To demonstrate the decrease in the variation of stiffness and magnetic force for high rotor frequency, the RAMB experiments were performed in such a way that would ensure that the magnetic force remains unchanged and the control current adjusts adaptively with a change in rotor speed. The experimental results showed that the control current would significantly rise with higher rotation frequency, which can indirectly illustrate the sharp drop of the magnetic force for RAMB.
Accordingly, the corresponding values of dynamic stiffness with varying rotor frequencies were calculated herein using 2D FEM and the calculated values were compared with those of the analytical results, as shown in Fig. 14. Here, note that the variation trends of dynamic stiffness were consistent with the theoretical analysis results. The maximum errors for current and displacement stiffness were 1.6% and 0.9%, respectively. Dynamic stiffness measurement method. It is clear from the above discussions that current and displacement stiffness are important parameters for characterising the performance of magnetic bearing and thus, stiffness measurement is a very essential research topic. This, however, most accounts on such a measurement method were limited at 0 Hz and have regarded the measured values as constants, which is obviously not feasible for high-speed conditions. To adequately verify the correctness of the proposed dynamic stiffness on this research as well as obtain the actual stiffness values for RAMB with varying rotation frequencies, a dynamic stiffness measurement method is proposed under this section.
In general, stiffness is measured on account of Eq. (19), which in turn, can be transformed to f = k x ·(h − h 0 ) + k i ·i. For this expression, f represents the electromagnetic force of the measured channel, h donates the rotor displacement, h 0 donates the magnetic centre position of RAMB, i donates the control current in the RAMB coils and k x and k i donate the measurement dynamic current and displacement stiffness, respectively. As the electromagnetic force of RAMBs in MSMP balances the gravity of the rotor with certain frequency, the electromagnetic force of each RAMB in the measured channel can be calculated with the solution f = f(mg) = G/2, with m indicating the rotor mass and G donating the rotor gravity. Consequently, the control current can be represented as: www.nature.com/scientificreports www.nature.com/scientificreports/ Under the assumption that the magnetic centre position h 0 is known, measurement of different displacements near h 0 and the corresponding control currents at any rotation speed will result to the finding the current and displacement stiffness at the rotation speed through Eq. (26). Particularly, the magnetic centre position needs only to be determined under static levitation according to its independence from the rotor speed.
Accordingly, there are two main stages in the proposed measurement method. As in the other methods described in [20][21][22] , the primary stage is determining the magnetic centre. If say, for the x-channel in RAMB the rotor is suspended stably along the x+ direction at 0 Hz, one coil current value of the RAMB can be obtained and then by inverting the rotor to stably along the x-direction, another coil current value can be obtained. If the suspended position of the rotor is adjusted until the two current values are equal, the rotor position that results could be recorded as the magnetic centre.
The secondary stage involves determining the dynamic current and displacement stiffness, which is much different from the literature [20][21][22] , by adjustment of the rotor frequency. When the rotor is suspended at the magnetic centre position, the corresponding control current for certain frequency ω is recorded as i 0 . Under this condition, the rotor is consequently adjusted to a position h j near h 0 and the corresponding control current is recorded as i j , where j = 1…n, that is, n equilibrium positions around h 0 are transformed through the adjustment of the reference value of the displacement sensor in the control programme for RAMB. According to Eq. (26), the corresponding current and displacement stiffness for ω can be expressed through Eqs. (27) and (28). Apparently, on the basis of the measured h 0 , h j , i 0 (ω) and i j (ω), the current and displacement stiffness of RAMB under different frequency can be obtained via changing the rotor speed and repeating the steps.
x j n i j 1 0 0 Measurement and experiments. The feasibility of the proposed dynamic stiffness characteristics and measurement method, the dynamic stiffness test experiment of RAMB was completed in a MSMP platform with a rated speed of 500 Hz (30000 r/min), as illustrated in Fig. 15. The MSMP model is as shown in Fig. 1 and the measured RAMB prototype and the centre of rotor mass are as depicted in Figs. 16 and 17, respectively. Figure 18 shows the principle block diagram of the RAMB control system, in which the control algorithm adopted in this work is an incomplete differential PID control strategy whose transfer function is reflected in Eq. (25). Upon application of the proposed dynamic stiffness model in the control model, the power amplifier can transform the output voltage signal of the controller into the current signal of the coil, in order to produce the required electromagnetic force. Table 2 provides a list of the main parameters. For the RAMB control system, the rotor displacement is directly measured with a displacement sensor and then compared to the reference signal of displacement set in the programme, after which the PID operation is performed so as to control the rotor to suspend stably at the reference position. Moreover, the reference signal can be regulated online via the control programme, which means that at a certain frequency, the change of rotor position can be achieved by adjusting the displacement reference signal online. Practically speaking, the measurement of dynamic stiffness for RAMB would be convenient and efficient with utilisation of the proposed measurement method.   www.nature.com/scientificreports www.nature.com/scientificreports/ During the experiment process, some factors as the interference caused by rotor unbalance will influence the performance of the RAMB system. To better simulate the system, two noise signals including a unit impulse signal and a white noise signal that result in a signal-to-noise ratio of 14 dB, were introduced for the simulation conducted in Matlab/Simulink. Table 3 provides the control parameters utilised while the representative simulation results obtained by online tuning are shown in Fig. 19. Both indicate that the RAMB system with the chosen control parameters can remain stable and is less affected by rotor unbalance and noise.
Dynamic current measurement and results. In the experiment, the magnetic centre was determined according to the primary stage described in the proposed measurement method, with the bias current set at 0.38 A. As a consequence, the control current can be measured when the rotor rotates at the magnetic centre of RAMB with different frequencies and the corresponding dynamic current stiffness under each speed can be obtained through Figure 18. Model of the RAMB control system.

Parameter
Value Parameter Value Table 3. Parameters of the PID controller. www.nature.com/scientificreports www.nature.com/scientificreports/ force analysis of the MSMP system and by using Eq. (27). The relationship between the measured dynamic current stiffness and rotor frequency relative to the analytical and FEM results is depicted in Fig. 12.
By adoption of the measurement method presented in literature [20][21][22] , the RAMB current stiffness can be obtained in the measurement test system described herein as 133.6 N/A, with no variation in the rotor frequency. Nonetheless, Fig. 20 clearly emphasises that the current stiffness measured by the proposed method would change with rotor speed, with the RAMB current stiffness in the rated speed of 500 Hz being 120.5 N/A. Moreover, comparing the measured values with the analytical values and FEM results, the maximum error between the measured values and analytical values was 5.0% (at 550 Hz) whereas the maximum error between the measured results and FEM results was 3.1% (at 550 Hz). This indicated the good agreement between the results of the experiment and those in both the analytical analysis and FEM results.
Dynamic displacement stiffness measurement and results. In order to obtain a convenient measurement of the dynamic displacement stiffness, the attitude of the prototype of the measured MSMP was adjusted so that the x-axis faces upward in the vertical direction and the magnetic centre of RAMB was set as the coordinate origin. Under such condition, the magnetic force in the RAMB would be zero when the rotor is suspended stably, such that the control current is i 0 = 0 with h 0 = 0. Consequently, Eq. (28) can be inverted into (29) yielding x j n i j 1 1 Finally, the dynamic displacement stiffness can be obtained by the recorded control current values via a series of set equilibrium positions under different rotor frequencies in accordance to Eq. (29). Figure 21 shows a comparison between the measured dynamic displacement stiffness against rotor frequency and the analytical and FEM results.
By adopting of the measurement method presented in literature [20][21][22] , the RAMB displacement stiffness obtained through the proposed measurement test system was 1.44 × 10 5 N/m, which is regarded as a constant and is unrelated with the rotor frequency. Nonetheless, it is discerned that the displacement stiffness of RAMB obtained using the proposed measurement method varied with rotor frequency, where the displacement stiffness of RAMB in the rated speed  www.nature.com/scientificreports www.nature.com/scientificreports/ of 500 Hz was 1.34 × 10 5 N/m. Furthermore, the experimental results were consistent with the analytical data and FEM results and produced low maximum errors of 2.3% at 550 Hz and 1.5% at 550 Hz. Figures 20 and 21 illustrate that the variation trends as rotor frequency increases for theoretical analysis, FEM and the measured results were consistent. Compared to the analytical and FEM results, the maximum error of measurement was 5.0% and 2.3%, respectively, for current stiffness and displacement stiffness, which is considerably within the acceptable range. This result supports that the proposed dynamic stiffness model and measurement method can perform well and is extremely significant for improving the dynamic performance of RAMB. Additionally, the errors for the three curves in Figs. 20 and 21 could be attributed to two aspects, namely: (a) The experiments were performed on a prototyped MSMP, which suggests that interference from the electric motor was unavoidable. To describe the scenario, the electric motor's position in Fig. 16 was quite close to the measured RAMB, by which the electromagnetic field generated by RAMB is expected to be superimposed with the electromagnetic field of the motor. Additionally, with high rotation speed, the motor would have a complex influence on RAMB for a varying magnetic field under eddy current effects. (b) Some ideal hypotheses were aimed to promote calculation convenience in the theoretical analysis part, say, the laminated material of the rotor was linear and the magnetic flux was a spatial function varying with time in the Z-direction.
The errors between the measured results and theoretical analysis can be minimise in future works by focusing on the separation design of the electric motor and measured RAMB, for elimination of the influence generated by the motor.

Discussion
The content of this work highlighted the proposal of an accurate model of dynamic stiffness, along with a measurement method of dynamic stiffness for RAMB. It conducted an investigation of dynamic stiffness characteristics influenced by rotor frequency for RAMB by application of the dynamic magnetic circuit method and FEM. Based on the results, dynamic stiffness exhibited remarkable variation with an increase in frequency. The proposed measurement method is a two-stage process. The magnetic centre is first determined so that the dynamic stiffness under a certain frequency can be obtained through measurement of different current values when the rotor is adjusted to different equilibrium positions near the magnetic centre. The method was validated by an experiment on a prototyped MSMP. Upon comparison of the experimental values to the theoretical results and the values obtained by FEM, as the rotor rotates at the same frequency, the maximum error of current stiffness was 5.0% and 3.1%, respectively, while the maximum error of the displacement stiffness was 2.3% and 1.5%.
The findings herein have shown great progress in the research of the dynamic properties of RAMB and have as well provided an evaluation standard for dynamic performance improvement and structure optimisation of RAMB.

Method
Theoretical analysis for dynamic stiffness. The theoretical analysis for the dynamic stiffness model was based on the RAMB structure in MSMP, as shown in Fig. 1. The corresponding dimensions are described in Table 1 in the main text.
FEM analysis for dynamic stiffness. The parameters in the FEM magnetic bearing model were consistent with the parameters in Table 1. The 2D FEM mesh model of RAMB in Fig. 2 was obtained as Fig. 12 in the main text, in which the plane53 magnetic structure unit was used for a perfect performance of transient magnetic field analysis. In addition, in the model the rotor part was meshed by mapped grid mode while the remaining part was meshed by free mode. In particular, the areas of air gap between the rotor surface and the fixed rotor, where the electromagnetic field exhibits such a great variation, were finely meshed. In the FEM simulations, the bias current was set at 0.38 A, which is expected to generate a corresponding bias magnetic flux density. The rotor was set with different rotor frequency values from 0 to 800 Hz, in the ANSYS software, followed by FEM simulation for each respective speed.

Dynamic stiffness measurement and experiments.
The experimental conditions and control system are shown in Fig. 15 through 18. The main coefficients are defined in Tables 2 and 3 in the main text.