Experimental Analysis and Full Prediction Model of a 5-DOF Motorized Spindle

The cost and power consumption of DC power amplifiers are much greater than that of AC power converters. Compared to a motorized spindle supported with DC magnetic bearings, a motorized spindle supported with AC magnetic bearings is inexpensive and more efficient. This paper studies a five-degrees-of-freedom (5-DOF) motorized spindle supported with AC hybrid magnetic bearings (HMBs). Most models of suspension forces, except a “switching model”, are quite accurate, but only in a particular operating area and not in regional coverage. If a “switching model” is applied to a 5-DOF motorized spindle, the real-time performance of the control system can be significantly decreased due to the large amount of data processing for both displacement and current. In order to solve this defect, experiments based on the “switching model” are performed, and the resulting data are analyzed. Using the data analysis results, a “full prediction model” based on the operating state is proposed to improve real-time performance and precision. Finally, comparative, verification and stiffness tests are conducted to verify the improvement of the proposed model. Results of the tests indicate that the rotor has excellent characteristics, such as good real-time performance, superior anti-interference performance with load and the accuracy of the model in full zone. The satisfactory experimental results demonstrate the effectiveness of the “full prediction model” applied to the control system under different operating stages. Therefore, the results of the experimental analysis and the proposed full prediction model can provide a control system of a 5-DOF motorized spindle with the most suitable mathematical models of the suspension force.


Introduction
Magnetically suspended high-speed motorized spindles display the advantages of being contact-free and wear-free; they also exhibit low noise, high rotational speed, compact structures, high mechanical efficiency, low vibration, and high precision, which has aided the further development of high-speed machining processes.Magnetic bearings are the core part of high-speed electrospindles and directly affect their properties.The most common magnetic bearings used to support magnetically suspended motorized spindles are DC magnetic bearings [1][2][3][4][5].In terms of cost savings and power consumption, AC magnetic bearings and hybrid magnetic bearings (HMBs) have obvious advantages compared with other magnetic bearings [6][7][8][9][10][11]. Therefore, AC hybrid magnetic bearings (AC HMBs) that combine the advantages of HMBs and AC magnetic bearings are worth attention [12][13][14][15].
This study investigates a five-degrees-of-freedom (5-DOF) motorized spindle supported with AC HMBs [16].The magnetic bearings used to support the 5-DOF motorized spindle are the AC 2-DOF HMB and the AC-DC 3-DOF HMB, both of which have perfect structures, accurate models and mature classical control technologies [17].Although these two magnetic bearing models are highly reliable, Energies 2017, 10, 75 2 of 17 not all model generations have been applied to the 5-DOF motorized spindle.Therefore, several new model generations are greatly anticipated in order to determine the most suitable and accurate motorized spindle model for a given working condition.
This paper establishes full model generations of the 5-DOF motorized spindle.The term "full" indicates that the model contains two parts: the first is the AC 2-DOF HMB model, and the other is the AC-DC 3-DOF HMB model.All the full models are quite accurate in a particular area of displacement and current.As for the full operating area of the motorized spindle, the most comprehensive and practical model is the "switching model".However, compared to the 3-DOF magnetic bearing system in [16], a "switching model" applied to a 5-DOF motorized spindle may adversely affect the real-time performance of the control system due to the large amount of data processing for displacement and current.Therefore, if a model based on a different operating state can be set in advance, the model will be more efficient and achieve better real-time performance in practice.In order to realize this goal, experiments based on the "switching model" are performed along with data analysis to be used as a basis for prediction.Using the data analysis results, a "full prediction model" based on the operating state is proposed to improve the real-time performance and precision of the motorized spindle system.
To verify the validity and accuracy of the "full prediction model", performance experiments are conducted and the results are analyzed.The results of comparative tests indicate that the control system based on the "full prediction model" results in the start-of-suspended state exhibit better real-time performance than the control system based on the "switching model".Results of the verification tests indicate that there is still good anti-interference performance, even with a heavier load for the rotor.In addition, the regulation times and real-time performances are comparable to the disturbance response experiment's results based on the "switching model".Moreover, the results of the stiffness tests show that the proposed full prediction model can provide a control system for the most suitable mathematical models of the suspension force, which is the closest to the actual condition model.Therefore, the effectiveness of the "full prediction model" applied to the control system can be verified by referring to the results of the performance test experiments under different operating stages.

Generations of Mathematical Models
Figure 1 shows the prototype and its exploded view of machinery parts of the 5-DOF AC HMB-supported motorized spindle.The AC 2-DOF HMB and AC-DC 3-DOF HMB models are responsible for supporting the shaft, as shown in Figures 2 and 3.
Energies 2017, 10, 75 2 of 17 and mature classical control technologies [17].Although these two magnetic bearing models are highly reliable, not all model generations have been applied to the 5-DOF motorized spindle.Therefore, several new model generations are greatly anticipated in order to determine the most suitable and accurate motorized spindle model for a given working condition.This paper establishes full model generations of the 5-DOF motorized spindle.The term "full" indicates that the model contains two parts: the first is the AC 2-DOF HMB model, and the other is the AC-DC 3-DOF HMB model.All the full models are quite accurate in a particular area of displacement and current.As for the full operating area of the motorized spindle, the most comprehensive and practical model is the "switching model".However, compared to the 3-DOF magnetic bearing system in [16], a "switching model" applied to a 5-DOF motorized spindle may adversely affect the real-time performance of the control system due to the large amount of data processing for displacement and current.Therefore, if a model based on a different operating state can be set in advance, the model will be more efficient and achieve better real-time performance in practice.In order to realize this goal, experiments based on the "switching model" are performed along with data analysis to be used as a basis for prediction.Using the data analysis results, a "full prediction model" based on the operating state is proposed to improve the real-time performance and precision of the motorized spindle system.
To verify the validity and accuracy of the "full prediction model", performance experiments are conducted and the results are analyzed.The results of comparative tests indicate that the control system based on the "full prediction model" results in the start-of-suspended state exhibit better real-time performance than the control system based on the "switching model".Results of the verification tests indicate that there is still good anti-interference performance, even with a heavier load for the rotor.In addition, the regulation times and real-time performances are comparable to the disturbance response experiment's results based on the "switching model".Moreover, the results of the stiffness tests show that the proposed full prediction model can provide a control system for the most suitable mathematical models of the suspension force, which is the closest to the actual condition model.Therefore, the effectiveness of the "full prediction model" applied to the control system can be verified by referring to the results of the performance test experiments under different operating stages.

Generations of Mathematical Models
Figure 1 shows the prototype and its exploded view of machinery parts of the 5-DOF AC HMB-supported motorized spindle.The AC 2-DOF HMB and AC-DC 3-DOF HMB models are responsible for supporting the shaft, as shown in Figures 2 and 3.

ClassicalFull Model
In summarizing the previous model's shortcomings for the AC 2-DOF HMB and AC-DC 3-DOF HMB, [12,14] respectively established comparatively accurate models.Combining the two models results in a full model covering the mathematical models of the 5-DOF HMBs, which are primarily covered by the classical full linear and nonlinear models.

Classical Full Nonlinear Model
The classical full nonlinear model of radial and axial suspension forces can be written as follows where

ClassicalFull Model
In summarizing the previous model's shortcomings for the AC 2-DOF HMB and AC-DC 3-DOF HMB, [12,14] respectively established comparatively accurate models.Combining the two models results in a full model covering the mathematical models of the 5-DOF HMBs, which are primarily covered by the classical full linear and nonlinear models.

Classical Full Nonlinear Model
The classical full nonlinear model of radial and axial suspension forces can be written as follows where

ClassicalFull Model
In summarizing the previous model's shortcomings for the AC 2-DOF HMB and AC-DC 3-DOF HMB, [12,14] respectively established comparatively accurate models.Combining the two models results in a full model covering the mathematical models of the 5-DOF HMBs, which are primarily covered by the classical full linear and nonlinear models.

Classical Full Nonlinear Model
The classical full nonlinear model of radial and axial suspension forces can be written as follows where Energies 2017, 10, 75 4 of 17 where µ 0 is the permeability in vacuum.In this model, one part of the parameters (i.e., the structure and electromagnetic parameters with subscript (1) belongs to the AC 2-DOF HMB.F cn1j is the classical nonlinear model of the radial suspension forces of the AC 2-DOF HMB; x 1 and y 1 are positive displacements of the rotor in the xand y-directions; ϕ 1j is the radial resultant magnetic fluxes in the air-gap corresponding to each pole; S r1 is the face area of the radial magnetic pole; δ r1 is the length of the uniform air-gap without rotor eccentricity; F m1 is the magnetomotive force provided to the outer circuit of the AC 2-DOF HMB; i j1 is the radial three-phase control current; j = A, B, C; N r1 is the number of turns that the radial control coils.The other part of the parameters (i.e., the structure and electromagnetic parameters with subscript (2) belongs to the AC-DC 3-DOF HMB.F cn2j and F cnz are the classical nonlinear model of the radial and axial suspension forces of the AC-DC 3-DOF HMB; x 2 and y 2 are positive displacements in the xand y-directions and z is right displacement of the rotor; ϕ 2j is the radial resultant magnetic fluxes in the air-gap corresponding to each pole; F m2 is the magnetomotive force provided to the outer circuit of the AC-DC 3-DOF HMB; i z is the axial control current and i j2 is the radial three-phase control current; j = A, B, C. F z1 is the right magnetic suspension force acting on the rotor; F z2 is the left magnetic suspension force acting on the rotor; N r2 is the number of turns that the radial control coils; N z is the number of turns that the axial control coils; ϕ z1 and ϕ z2 are the axial resultant magnetic-biased fluxes generated by the permanent magnet and axial control coils in the right and left axial air-gaps, respectively; G z1 and G z2 describe the right and the left air-gap permeances, respectively; G A , G B , and G C describe the radial air-gap permeances.

Classical Full Linear Model
The classical full linear models of radial and axial suspension forces (F cl1x , F cl1y , F cl2x , F cl2y and F clz ) can be written as follows: Similarly, in this model, one part of parameters with subscript1 belongs to the AC 2-DOF HMB: i x1 , i y1 are the current components in the xand y-axis transformed from the three-phased currents by Clark coordinate transformation; k cl1xy is the radial force-displacement coefficient; k cl1ir is the radial force-current coefficient.The other part of the parameters with subscript2 belongs to the AC-DC 3-DOF HMB: i x2 , i y2 are the current components in the xand y-axis transformed from the three-phased currents by Clark coordinate transformation; k cl2xy is the radial force-displacement coefficient; k cl2ir is the radial force-current coefficient; k clz is the axial force-displacement coefficient, and k cliz is the axial force-current coefficient.
where S r2 is the face area of the radial magnetic pole; δ r2 is the length of the uniform air-gap without rotor eccentricity and δ z is the axial air-gap length of the AC-DC 3-DOF HMB; S z is the face area of the axial magnetic pole of the AC-DC 3-DOF HMB; and S mag is the average area of the inner and outer faces of the permanent magnet.

Improved Full Model
In seeking further improvement, the magnetic circuit analysis process focuses on error capture and analysis.In [14], the negligence of permanent magnetic reluctance in analyzing the equivalent magnetic fluxes path limited accurate assessments in the modeling of AC-DC 3-DOF HMB.Therefore, [15] proposed an improved model for AC-DC 3-DOF HMB to address these defects.
Other than improving the magnetic circuit, a breakthrough in principle analysis that perfected the model was also achieved.In considering the similarities of the working principle of radial AC magnetic bearings and the radial suspension subsystem of a bearing-free motor, a new modeling method based on the Maxwell tensor method for the radial suspension force of an AC magnetic bearing was presented in [13], and the AC 2-DOF HMB was modeled by this method.The experimental results showed that the invention designed especially for AC magnetic bearings overcame inaccuracies and possessed the advantages of directness and universality compared to the previous method.
In combining the two models, a full model covering the mathematical models of the 5-DOF HMBs can be obtained that is primarily covered by improved full linear and nonlinear models.

Improved Full Nonlinear Model
The improved full nonlinear model of the radial and axial suspension forces can be written as follows: Energies 2017, 10, 75 Similarly, in this model, one part of the parameters with subscript1 belongs to the AC 2-DOF HMB: F in1j is the improved full nonlinear model of the radial suspension forces generated by the resultant magnetic fluxes of each radial air-gap; l is the equivalent length of rotor; r is radius of rotor; dS r1 is per unit area; and θ is the dimensional mechanical angle.With the rotor eccentricity at time t, the expression of the magnetic flux density in air-gap B (θ, t) is determined by the permanent magnet and the control current.The other part of the parameters with subscript2 belongs to the AC-DC 3-DOF HMB: F in2j is the improved full nonlinear model of the radial suspension forces generated by the resultant magnetic fluxes of each radial air-gap; F inz is the improved full nonlinear model of the axial suspension force; G z1 and G z2 describe the right and the left air-gap permeances, respectively; G A , G B , and G C describe the radial air-gap permeances; G mag describes the permanent magnetic permeance; ϕ cz is the axial control magnetic fluxes generated by the control coils; ϕ pz1 and ϕ pz2 are the axial magnetic biased fluxes generated by the permanent magnet in right and left axial air-gaps, respectively.

Improved Full Linear Model
The improved linear model of the radial and axial suspension forces (F il1x , F il1y , F il2x , F il2y and F ilz ) can be written as follows: where S mag is the average area of the inner and outer faces of the permanent magnet; and δ mag is the vertical distance from the outer surface to inner surface of the permanent magnet.
In the model, one part of the parameters with subscript1 belongs to the AC 2-DOF HMB: k il1xy is the radial force-displacement coefficient; k il1ir is the radial force-current coefficient.The other part of the parameters with subscript2 belongs to the AC-DC 3-DOF HMB: k il2xy is the radial force-displacement coefficient; k il2ir is the radial force-current coefficient; k ilz is the axial force-displacement coefficient; and k iliz is the axial force-current coefficient.

SwitchingFullModel
The suspension force models mentioned above are quite accurate, but only in a particular operating area and not in regional coverage; thus, the magnetic bearing model must gradually improve and expand.In [16], a "switching model" was proposed for the AC-DC 3-DOF HMB that involved a wide range of displacements and currents for the rotor.The concept of a "switching model" can then be applied to a 5-DOF HMB to form a full operation-domain model; this will make the mathematical model of suspension forces for the 5-DOF motorized spindle more complete and accurate.Referring to the existing analysis result of "switching model" for the AC-DC 3-DOF HMB in [16], the "full switching model" of the 5-DOF AC HMB-supported motorized spindle selection table Energies 2017, 10, 75 7 of 17 can be derived, as shown in Table 1.There are still four submodels in the composite "full switching model", which are called Model (1), Model (2), Model (3) and Model (4).The four submodels in the different current and displacement regions consist of combinations of the model type (improved full nonlinear model and classical full nonlinear model) and the overlapping combination mode.It is noteworthy that the model selection in Zone (3) can be further divided according to the last model state for the different displacements and current regions.The selection algorithm is as follows: (a) If the last model is Model (1), Model (3) will be selected for efficient and simplified control.Generally, this type of model selection will occur in the cases of a suspension experiment or disturbance response experiment; (b) If the last model is Model (2) or Model (4), Model (3) will be selected for efficient and simplified control.Generally, this type of model selection will occur in the cases in which the rotor is in the start-of-suspension state or in the case of the middle transition process for a rotor going back and forth between the equilibrium position and a distance.
Using different rotor displacements and control currents, the mathematical model of the suspension forces used in the fore-current transformation link of the control system can be switched at any time.However, theoretically, more complex than the 3-DOF HMB system, the real-time performance of the control system for the 5-DOF motorized spindle will be decreased due to the extra 2-DOF control inputs and outputs.Therefore, a large amount of displacement and current measurement time reduction is demanded because of the constant judging and selection for switching models to achieve the stable and rapid control effect.Thus, the 5-DOF motorized spindle model still needs improvement by analyzing the switching full model based on experiments.

Switching Full Model Analysis Based on Experiments
Unlike the "switching model" for the AC-DC 3-DOF HMB, a full model integrating the most accurate AC 2-DOF HMB and AC-DC 3-DOF HMB system model is established.Although the "switching model" has been proven accurate in an AC-DC 3-DOF HMB system in [16], whether the "full switching model" is useful in 5-DOF HMBs is unproven and the real-time performance should be emphatically investigated.Therefore, the related experimental data must be selected to determine a more effective modeling method.Therefore, the following experiments are performed to enable data analysis.The experimental platform of the 5-DOF motorized spindle supported with AC HMBs is shown in Figure 4.The block diagram of the control algorithm for the control system adopting the "full switching model" for the spindle is presented in Figure 5. is shown in Figure 4.The block diagram of the control algorithm for the control system adopting the "full switching model" for the spindle is presented in Figure 5. (2) Switching model (2)

Radial force-current transformation
Switching model decision table (2)

Radial force-current transformation
Switching model is shown in Figure 4.The block diagram of the control algorithm for the control system adopting the "full switching model" for the spindle is presented in Figure 5. (2) Switching model (2)

Radial force-current transformation
Switching model decision table (2)

Radial force-current transformation
Switching model The basic block of control algorithm for control system.
Figure 5 illustrates the position and current feedback signals that are processed using the decision table when the rotor is displaced from an equilibrium position by an external disturbance.Meanwhile, the position deviation signals (i.e., the difference between the equilibrium position reference signals and the position feedback signals) are sent to the PID controllers.
For the radial and axial control parts of the magnetic bearings, the PID controllers transform radial and axial deviation signals into force signals (F xx1 , F yy1 , F xx2 , F yy2 and F z ) using the force-current transformation modules.Notably, the radial and axial force-current transformation modules reflect the efforts of the mathematical suspension force models.Therefore, the mathematical models of the suspension forces are important to obtain precise control of these magnetic bearings.The key mathematical expressions (i.e., mathematical models of the suspension forces) in the algorithm of the control system of the 5-DOF motorized spindle are as mentioned above.
These modules transform radial and axial force signals (F xx1 , F yy1 , F xx2 , F yy2 , and F zz ) into radial and axial control current reference signals (i A1 *, i B1 *, i C1 *, i A2 *, i B2 *, i C2 *, and i z *) for the 5-DOF HMBs.The control current reference signals function as the reference signals of the internal closed loop current control link.Moreover, the exciting currents can be obtained using the current-controlled voltage-source inverter and axial switching power amplifier.Through this process, control currents are adjusted using the control system with a negative position and current feedbacks to readjust flux distribution.Thus, the accuracy of a suspension forces model directly affects the effectiveness of the control system.

Start-of-Suspension Response Experiment
The dynamic stability of the magnetic-rotor system can be observed in the start-of-suspension response curve of the rotor; its purpose is to realize the preparation for the high-speed rotor test and test the effectiveness of the full models across the multi-zone.Figure 6a shows the mass center orbit of the start-of-suspension characteristics of the rotor under the support of the AC 2-DOF HMB. Figure 6b shows the mass orbit of the start-of-suspension characteristics of the rotor under the support of the AC-DC 3-DOF HMB. Figure 6c shows the displacement waveforms of the rotor of the start-of-suspension characteristics of the rotor under the support of the AC-DC 3-DOF HMB.
The first observation is the original position for the rotor.Prior to introducing electricity, the rotor stops on the backup bearing at a random point.In order to gather as much different data and make the rotor cross as many displacement and current regions as possible, the lengths of the radial and axial air-gaps are both set to 0.49 mm-slightly lower than the maximum air-gap length-for the safety and limit test.
At one end of the rotor supported by the AC 2-DOF HMB, the rotor starts at a position with initial coordinates of x = 0.20 mm and y = 0.46 mm, as shown in Figure 6a.At the other end of the rotor supported by the AC-DC 3-DOF HMB, the rotor starts at a position with initial coordinates of x = −0.24mm, y = 0.43 mm, and z = 0.067 mm, as shown in Figure 6b.Both radial and axial eccentricity deviations can be seen from the equilibrium position of the rotor.Thus, the torque coupling must be considered when designing the controller and regulating parameters.
In the experiment, as can be seen in Figure 6a, for the AC 2-DOF HMB, the initial displacement in x-direction is 0.20 mm (belong to Zone (1)) and the displacement in y-direction is 0.46 mm (belong to Zone (4)).Model (1) and Model (4) are selected by referring to Table 1.As can be seen in Figure 6b, the initial displacement in the x-direction is −0.24 mm, the displacement in y-direction is 0.43 mm and the displacement in z-direction is 0.067 mm for the AC-DC 3-DOF HMB.Model (1) and Model (3) are then selected by referring to Table 1.Based on the different initial positions in the experiment, Model (1), Model (3) and Model (4) were selected as controls.
The next observation is the start-of-suspension process for the rotor.It can be seen in Figure 6 that each end of the rotor starts simultaneously from a stationary position controlled by the two magnetic bearings and then rapidly and smoothly moves to an equilibrium position after a slight oscillation under current control.Model (1) was selected as the control in the experiment.Therefore, as the rotor approaches equilibrium and enters Zone (1), the displacement and current will be more specifically suited to Model (1), which is designated as the control experience.
The final observation is the return process for the rotor.Based on Figure 6a,b, it can be seen that the rotor oscillation amplitude near the equilibrium position in Figure 6b is greater than that in Figure 6a.The reason for this is that Figure 6b,c shows the radial and axial start-of-suspension situation separately for the same AC-DC 3-DOF HMB-only the view in angle is different.Because the axial sensor is Energies 2017, 10, 75 10 of 17 installed near one side of the AC-DC 3-DOF HMB and the influence produced by the installation position of the sensor on the offset of the rotor near the sensor is greater, the influence of distance to the sensor is slight.In addition, Figure 6c shows that the rotor realizes its axial start-of-suspension under the support of the AC-DC 3-DOF HMB with a slight oscillation when it has almost reached the equilibrium position; therefore, the rotor is also affected in the radial direction.However, it can be seen that the axial disturbance of the rotor shows less of an effect on the other end of the rotor supported by the AC 2-DOF HMB in the radial direction due to the installation location factors of the axial sensor.
Energies 2017, 10, 75 10 of 17 the axial sensor is installed near one side of the AC-DC 3-DOF HMB and the influence produced by the installation position of the sensor on the offset of the rotor near the sensor is greater, the influence of distance to the sensor is slight.In addition, Figure 6c shows that the rotor realizes its axial start-of-suspension under the support of the AC-DC 3-DOF HMB with a slight oscillation when it has almost reached the equilibrium position; therefore, the rotor is also affected in the radial direction.However, it can be seen that the axial disturbance of the rotor shows less of an effect on the other end of the rotor supported by the AC 2-DOF HMB in the radial direction due to the installation location factors of the axial sensor.In conclusion, the experimental results show that the control system based on the "full switching model" crossing multi-zone has good start-of-suspension response characteristics.The model-conversion mode during the start-of-suspension process can be a useful control experience in forming a new model.

Suspension Experiment
The static stability of the magnetic-rotor system can be observed by the characteristic curve of the static suspension of the rotor.In order to operate the spindle at the maximum speed of 60,000 rpm, the indispensable modal frequencies must be steadily crossed.The backup bearing is installed to demonstrate the effectiveness during the test.In addition, the rotor is operated in Zone (1), so the displacement and current will be more specifically suited to Model (1) as referenced in Table 1.
Figure 7a shows the mass center orbit of the suspended rotor under the support of the AC 2-DOF HMB. Figure 7b shows the mass center orbit of the suspended rotor under the support of the AC-DC 3-DOF HMB. Figure 7c shows the displacement waveforms of the suspended rotor under the support of the AC-DC 3-DOF HMB.The purpose of the test is to observe the suspension characteristics of the rotor.Based on the rotor trajectories, it can be seen that the rotor in the In conclusion, the experimental results show that the control system based on the "full switching model" crossing multi-zone has good start-of-suspension response characteristics.The model-conversion mode during the start-of-suspension process can be a useful control experience in forming a new model.

Suspension Experiment
The static stability of the magnetic-rotor system can be observed by the characteristic curve of the static suspension of the rotor.In order to operate the spindle at the maximum speed of 60,000 rpm, the indispensable modal frequencies must be steadily crossed.The backup bearing is installed to demonstrate the effectiveness during the test.In addition, the rotor is operated in Zone (1), so the displacement and current will be more specifically suited to Model (1) as referenced in Table 1.
Figure 7a shows the mass center orbit of the suspended rotor under the support of the AC 2-DOF HMB. Figure 7b shows the mass center orbit of the suspended rotor under the support of the AC-DC 3-DOF HMB. Figure 7c shows the displacement waveforms of the suspended rotor under the support of the AC-DC 3-DOF HMB.The purpose of the test is to observe the suspension characteristics of the rotor.Based on the rotor trajectories, it can be seen that the rotor in the equilibrium position realizes the stability of suspension in 5-DOF directions with small fluctuations around the equilibrium position; this shows that the system achieves a stable five-degrees-of-freedom suspension.
Energies 2017, 10, 75 11 of 17 equilibrium position realizes the stability of suspension in 5-DOF directions with small fluctuations around the equilibrium position; this shows that the system achieves a stable five-degrees-of-freedom suspension.
x displacement (mm) y displacement (mm) In conclusion, the experimental results show that the control system based on the "full switching model" crossing the multi-zone (there is only one model here) demonstrates good suspension response characteristics.The model-conversion mode during the suspension process can be used as a control experience when forming a new model.

Disturbance Response Experiment
The dynamic stability of the magnetic-rotor system can also be observed through the anti-turbulence response characteristic curves of the rotor.When the rotor realizes stable suspension, an external disturbance suddenly imposed on the rotor in a certain direction can test the dynamic performance of the system by observing the changes to the displacement and control current.The disturbance test is designated at the end of the rotor supported by the AC-DC 3-DOF HMB.In the test, the backup bearing (the radius of the backup bearing is 0.25 mm) is installed to demonstrate the effectiveness.In addition, the rotor is operated and disturbed in Zone (1); therefore, the displacement and current will be more specifically suited to Model (1) as referenced in Table 1.In conclusion, the experimental results show that the control system based on the "full switching model" crossing the multi-zone (there is only one model here) demonstrates good suspension response characteristics.The model-conversion mode during the suspension process can be used as a control experience when forming a new model.

Disturbance Response Experiment
The dynamic stability of the magnetic-rotor system can also be observed through the anti-turbulence response characteristic curves of the rotor.When the rotor realizes stable suspension, an external disturbance suddenly imposed on the rotor in a certain direction can test the dynamic performance of the system by observing the changes to the displacement and control current.The disturbance test is designated at the end of the rotor supported by the AC-DC 3-DOF HMB.In the test, the backup bearing (the radius of the backup bearing is 0.25 mm) is installed to demonstrate the effectiveness.In addition, the rotor is operated and disturbed in Zone (1); therefore, the displacement and current will be more specifically suited to Model (1) as referenced in Table 1.
Figure 8 shows the wave forms of the shock test.In Figure 8, a sudden knock by a wooden hammer is imposed on one end of the shaft in the axial direction; it can be seen that the rotor requires 40 ms to reacquire stable suspension.In Figure 9, interference by a constant force of 100 N in the negative y-direction causes the rotor to deviate from the equilibrium position and then revert to the original seat after a period of adjustment by the control system.Simultaneously, the disturbance imparted in the y-direction has little effect on the performance of the rotor in the x-direction.In addition, once the disturbance disappears, the displacements in both directions quickly revert to the equilibrium position.Based on the experimental results, it can be seen that the coupling of the two degrees of freedom (in xand y-direction) is not obvious; the scope of rotor activity remains close to the equilibrium position due to the limit of the backup bearing.The experimental results show that the designed control system provides the rotor with good anti-interference performance.
Figure 8 shows the wave forms of the shock test.In Figure 8, a sudden knock by a wooden hammer is imposed on one end of the shaft in the axial direction; it can be seen that the rotor requires 40 ms to reacquire stable suspension.In Figure 9, interference by a constant force of 100 N in the negative y-direction causes the rotor to deviate from the equilibrium position and then revert to the original seat after a period of adjustment by the control system.Simultaneously, the disturbance imparted in the y-direction has little effect on the performance of the rotor in the x-direction.In addition, once the disturbance disappears, the displacements in both directions quickly revert to the equilibrium position.Based on the experimental results, it can be seen that the coupling of the two degrees of freedom (in x-and y-direction) is not obvious; the scope of rotor activity remains close to the equilibrium position due to the limit of the backup bearing.The experimental results show that the designed control system provides the rotor with good anti-interference performance.In conclusion, the instantaneous and delayed disturbance response experiments show that a rotor control based on the "switching model" crossing a multi-zone (there is only one model here) can reacquire stable suspension in a short time.The model-conversion mode during the disturbance response process can be used as a control experience when forming a new model.

Full Prediction Model of Suspension Force according to Operating State
Based on the results of the above experiments and the analysis, certain useful control experiences can be summarized.Other similar performance experiments have been performed to collect more information.Based on the experimental results, it can be seen that the model Figure 8 shows the wave forms of the shock test.In Figure 8, a sudden knock by a wooden hammer is imposed on one end of the shaft in the axial direction; it can be seen that the rotor requires 40 ms to reacquire stable suspension.In Figure 9, interference by a constant force of 100 N in the negative y-direction causes the rotor to deviate from the equilibrium position and then revert to the original seat after a period of adjustment by the control system.Simultaneously, the disturbance imparted in the y-direction has little effect on the performance of the rotor in the x-direction.In addition, once the disturbance disappears, the displacements in both directions quickly revert to the equilibrium position.Based on the experimental results, it can be seen that the coupling of the two degrees of freedom (in x-and y-direction) is not obvious; the scope of rotor activity remains close to the equilibrium position due to the limit of the backup bearing.The experimental results show that the designed control system provides the rotor with good anti-interference performance.In conclusion, the instantaneous and delayed disturbance response experiments show that a rotor control based on the "switching model" crossing a multi-zone (there is only one model here) can reacquire stable suspension in a short time.The model-conversion mode during the disturbance response process can be used as a control experience when forming a new model.

Full Prediction Model of Suspension Force according to Operating State
Based on the results of the above experiments and the analysis, certain useful control experiences can be summarized.Other similar performance experiments have been performed to collect more information.Based on the experimental results, it can be seen that the model In conclusion, the instantaneous and delayed disturbance response experiments show that a rotor control based on the "switching model" crossing a multi-zone (there is only one model here) can reacquire stable suspension in a short time.The model-conversion mode during the disturbance response process can be used as a control experience when forming a new model.

Full Prediction Model of Suspension Force according to Operating State
Based on the results of the above experiments and the analysis, certain useful control experiences can be summarized.Other similar performance experiments have been performed to collect more information.Based on the experimental results, it can be seen that the model transformation time and response time-even in milliseconds-can play a significant role in the control process.Therefore, a composite model such as a "full prediction model" of a suspension force based on an operating state is an ideal way to simplify the model.The mathematical model of suspension forces used in oscillations and less adjustment time when almost reaching the equilibrium position; thus, the real-time performance is improved.
Thus, it is revealed by the comparison between Figure 6 and 10 that the control system based on the "full prediction model" results in the start-of-suspended rotor having a better performance of real-time processing than the control system based on the "switching model."

Verification Experiment
In the test, the "full prediction model" of suspension force according to operating state is used in the control system, where the models according to the different operating state were previously set.By checking Table 2 and detecting the displacement and current detection, the initial model is set as Model (1), and the maximum deviation position is set as Model (4), and then the Model variation is set as Model (1)→Model (2-4)→Model (1).
In the disturbance response test, the backup bearing (the radius of the backup bearing is 0.49 mm) is installed to demonstrate the effectiveness of the model analysis results that are applied to the control strategies.Unlike the previous section, the activity range of the rotor extends to the nonlinear zone because of the change to the backup bearing.The contrasting experimental conditions remain the same except for the different mathematical models.In the test, a constant external force of 600 N is instantaneously applied in the negative y-direction.In Figure 11, the rotor x-direction almost steps away from its equilibrium position instantaneously with the advent of interference.Meanwhile, the performance of the rotor in the x-direction is relatively unscathed by the disturbance imparted in the negative y-direction.Therefore, for the rotor in Figure 11, superior anti-interference performance is demonstrated, even with a heavier load.In the test, the regulation time and real-time performance of the test results obtained with a load disturbance of 600 N in Figure 11 are comparable to the disturbance response experiment results obtained with a load disturbance of 100 N in Figure 9.

Verification Experiment
In the test, the "full prediction model" of suspension force according to operating state is used in the control system, where the models according to the different operating state were previously set.By checking Table 2 and detecting the displacement and current detection, the initial model is set as Model (1), and the maximum deviation position is set as Model (4), and then the Model variation is set as Model (1)→Model (2-4)→Model (1).
In the disturbance response test, the backup bearing (the radius of the backup bearing is 0.49 mm) is installed to demonstrate the effectiveness of the model analysis results that are applied to the control strategies.Unlike the previous section, the activity range of the rotor extends to the nonlinear zone because of the change to the backup bearing.The contrasting experimental conditions remain the same except for the different mathematical models.In the test, a constant external force of 600 N is instantaneously applied in the negative y-direction.In Figure 11, the rotor x-direction almost steps away from its equilibrium position instantaneously with the advent of interference.Meanwhile, the performance of the rotor in the x-direction is relatively unscathed by the disturbance imparted in the negative y-direction.Therefore, for the rotor in Figure 11, superior anti-interference performance is demonstrated, even with a heavier load.In the test, the regulation time and real-time performance of the test results obtained with a load disturbance of 600 N in Figure 11 are comparable to the disturbance response experiment results obtained with a load disturbance of 100 N in Figure 9.
Figure 12a shows that the force errors between the experimental and full prediction models are near zero in the control current range of 0 A to 1 A. By contrast, the difference in force errors between the experimental and traditional models (classical and improved full models) become larger over the same linear zone (0 A to 1 A).With a further increase in the control current range to between 1 A and 2 A, the calculation results of the full prediction model are still concluded to be closer to the experiment results than the other two traditional models within the nonlinear zone of the current (1 A to 2 A).
Figure 12b shows that the force error between the experimental and full prediction models is nearly zero within the displacement range from 0.0 mm to 0.25 mm (in the range of linear model), especially from 0.0 mm to 0.1 mm.By contrast, the difference in force errors between the experimental and traditional models (classical and improved full models) becomes larger in the same linear zone.With a further increase in the displacement range of 0.25 mm to 0.49 mm, the calculation results of the full prediction model are concluded to be much closer to the experiment results compared to the other two traditional models within the nonlinear zone of displacement (0.25 mm to 0.49 mm).
Therefore, the results of the model analysis show the proposed full prediction model can provide a control system for the most suitable mathematical models of the suspension force in the full zone.

Conclusions
In the paper, the full models of the five-degrees-of-freedom (5-DOF) motorized spindle are established, which include the most classical model, an improved model and the "full switching model".Most suspension force models, except for the "switching model", are quite accurate, but only in a particular operating area and not in regional coverage.If a "switching model" is applied to a 5-DOF motorized spindle, the real-time performance of the control system can be significantly decreased due to the large amount of data processing for both displacement and current.In order to solve this defect above, the performance experiments based on "switching model" have been undertaken for data analysis.Subsequently, based on the data analysis results, a "full prediction model" according to operating state is proposed to improve real-time performance and precision of the motorized spindle system, where the model can be previously set according to the different operating state.In addition, some comparative, verification and stiffness tests are conducted to verify the improvement of the proposed model.The satisfactory comparative experimental results demonstrate the "full prediction model" applied to the control system under different operating stages has good real-time performance.The fine verification experimental results demonstrate the "full prediction model" has strong robustness and real-time performance can still be guaranteed.The stiffness test results verify the accuracy of the proposed model.Therefore, the results of the experimental analysis and the proposed full prediction model can provide a control system of a 5-DOF motorized spindle with the most suitable mathematical models of the suspension force.

Figure 5 .
Figure 5.The basic block of control algorithm for control system.

Figure 5
Figure 5 illustrates the position and current feedback signals that are processed using the decision table when the rotor is displaced from an equilibrium position by an external disturbance.Meanwhile, the position deviation signals (i.e., the difference between the equilibrium position reference signals and the position feedback signals) are sent to the PID controllers.

Figure 5 .
Figure 5.The basic block of control algorithm for control system.

Figure 5
Figure 5 illustrates the position and current feedback signals that are processed using the decision table when the rotor is displaced from an equilibrium position by an external disturbance.Meanwhile, the position deviation signals (i.e., the difference between the equilibrium position reference signals and the position feedback signals) are sent to the PID controllers.

Figure 6 .
Figure 6.Trajectories and displacement waves of the rotor adopting "switching model" at the start-of-suspension state.(a) Mass center orbit of the rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Figure 6 .
Figure 6.Trajectories and displacement waves of the rotor adopting "switching model" at the start-of-suspension state.(a) Mass center orbit of the rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Figure 7 .
Figure 7. Radial trajectories of the rotor at the stable suspension state.(a) Mass center orbit of the suspended rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the suspended rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Figure 7 .
Figure 7. Radial trajectories of the rotor at the stable suspension state.(a) Mass center orbit of the suspended rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the suspended rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Figure 8 .
Figure 8.Waveforms of impact test in axial direction.

Figure 9 .
Figure 9. Waveforms of disturbance test in radial direction.

Figure 8 .
Figure 8. Waveforms of impact test in axial direction.

Figure 8 .
Figure 8.Waveforms of impact test in axial direction.

Figure 9 .
Figure 9. Waveforms of disturbance test in radial direction.

Figure 9 .
Figure 9. Waveforms of disturbance test in radial direction.

Figure 10 .
Figure 10.Trajectories and displacement waves of the rotor adopting "full prediction model" at the start-of=suspension state.(a) Mass center orbit of the rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Figure 10 .
Figure 10.Trajectories and displacement waves of the rotor adopting "full prediction model" at the start-of=suspension state.(a) Mass center orbit of the rotor under the support of the AC 2-DOF HMB; (b) Mass center orbit of the rotor under the support of the AC-DC 3-DOF HMB; (c) Displacement waveforms of the rotor under the support of the AC-DC 3-DOF HMB.

Table 1 .
The "switching model" across multi-zone selection table.