Fuzzy Bang-Bang Relay Control of a Rigid Rotor Supported by Active Magnetic Bearings

: Active magnetic bearings, which are open-loop and unstable, require a feedback control system to ensure stable operation of the rotating machines that they support. Proportional-integral-derivative (PID) controllers are widely used in ﬁeld applications of these bearings for this purpose. PID controllers are designed to work effectively within the linear region of operation of the rotating machines. Due to the inherent nonlinearity of the active magnetic bearings, large unbalance forces that may occur in these machines result in nonlinear vibration responses. Therefore, the PID controller’s effectiveness to control the vibration of the rotating machines is considerably reduced when the unbalance forces in these machines become large. Other control strategies, such as the fuzzy logic and the sliding mode control schemes, are more apt to deal with the nonlinear responses of the rotating machines supported by active magnetic bearings. The present work proposes an integrated fuzzy bang-bang relay controller for a rigid rotor mounted on active magnetic bearings. The effectiveness of this controller to suppress rotor vibrations is examined numerically. Performance comparison of this controller with the conventional fuzzy logic and PD controllers are made for different initial conditions, rotor imbalance magnitudes, and rotor angular speeds. At extreme operating conditions due to large rotor unbalance forces, where the magnetic bearings are highly nonlinear, the proposed integrated fuzzy bang-bang relay controller proved to be more superior over the conventional fuzzy logic and PD controllers.


Introduction
The utilization of active magnetic bearings in rotating machines has increased over the past years.They are used in applications where they have advantages over the rollingelement bearing and fluid-film bearing types.Some of the main advantages of the active magnetic bearings are their oil-free operation and their higher mechanical efficiency.Lubrication is not needed for these bearings.Contact between the rotor and stator also do not occur during the operation of these bearings.Another major advantage of the magnetic bearings over the conventional bearing types, from the rotor dynamics viewpoint, is the possibility of the dynamic characteristics to be varied without changing the dimensions of the rotor or bearing.The dynamic characteristics of the magnetic bearings are governed by the parameters of their controllers.By varying these parameters, the dynamic characteristics of these bearings can be changed over a considerably wide operating range.A major disadvantage of the active magnetic bearings is their highly nonlinear characteristics, which are attributed to several inherent sources.
The main causes of nonlinear behavior in magnetic bearings are the relationship between the magnetic-bearing force with the coil current and the air gap between the rotor and the stator.The force is inversely proportional to the square of the air gap and proportional to the square of the current in the coil.Another cause of nonlinearity in these bearings is the cross-coupling that arises from the magnetic forces that act in the two orthogonal directions.Cross-coupling forces may occur in these bearings due to the geometry of the magnetic-bearing actuators, the effect of gyroscopic motion, and the effect of the eddy current in the magnetic actuator.Other sources of nonlinearities in magnetic bearings include hysteresis and saturation of the magnetic core material.The nonlinear nature of the active magnetic bearings usually results in non-synchronous vibration in the response of the rotor supported by them for certain ranges of the design and operating parameter values.Non-synchronous vibrations, which induce cyclic stresses, should be prevented in the rotating machines' operations.Cyclic or fluctuating stresses may shorten the fatigue life of the machine components and may eventually result in their failure.
Due to being open-loop unstable, magnetic bearings require feedback control for stable operation.The most common control strategy used for field applications of magnetic bearings is the linear proportional-integral-derivative (PID) controller [1].The effective working operation of the PID controller is however limited to the proximity of the equilibrium point where the linearization of its dynamics is performed.Deviation of the magnetic-bearing operation from the equilibrium point may render the PID controller to be ineffective.Such deviation may occur in the field operation of rotating machines mounted in magnetic bearings due to the occurrence of large forces acting on these machines.Some of the common causes of these forces include rotor imbalance and flow-induced vibrations, amongst others.Works in the past have proposed various nonlinear control techniques to address the occurrence of nonlinear vibrations in magnetic bearings.A brief review of the applications of these techniques, in particular the fuzzy logic and sliding mode control (SMC) schemes, to the control of magnetic bearings follows.The recent paper by Srinivas et al. [2] may be referred to for a more exhaustive review of the application of magnetic bearings in rotating machines.Huang et al. [3], on the other hand, provides a focused review on the control technology for magnetic bearings.
An adaptive fuzzy scheme was employed by Liebert [4] in magnetic-bearing application.In this scheme, the gains of the linear PID controller were varied using fuzzy logic.Measurements undertaken on a magnetic-bearing-test rig demonstrated that the proposed scheme had a better step response when compared to the use of the conventional linear PID controller.Hung [5] proposed the use of a variable structure scheme to control the motion of a single-axis magnetic-bearing system.This scheme made use of two different controllers.For operation in the vicinity of the equilibrium point, a linear PID controller was used.A nonlinear feedback-linearizing controller was utilized for the operation of the bearing away from the equlibrium point.In order to have a seamless transition between these two controllers, a fuzzy control scheme was used.The fuzzy control scheme proposed in this work was also utilized to prevent chattering from occuring in the response of the system.The Takagi-Sugeno-Kang fuzzy model was utilized by Hong et al. [6] in the development of a scheme to control nonlinear vibrations in magnetic bearings.The scheme was found to perform better compared to the single operating point linear controller based on the outcomes of numerical studies.The magnetic bearing's regime of stable operation was also found to be at a maximum with the use of the proposed scheme.Experimental work undertaken by Hong and Langari [7] also confirmed the superior performance of the proposed control scheme [7].Hong and Langari [8] further extended their proposed control scheme and demonstrated that this scheme was robust to parameter uncertainties and disturbances that were harmonic in nature.The effective performance of this control scheme was further demonstrated through numerical simulation as well as measurement on an experimental setup.Kosaki et al. [9] proposed a multivariable fuzzy controller for rotors mounted on magnetic bearings subjected to cross-coupled gyroscopic motion.The proposed control scheme was found to be able to maintain the system's performance when the speed of the rotor was varied and when external disturbance was imposed on it.The numerical work of Lei et al. [10] on the utilization of a magnetic bearing as an exciter for an active rotating-stall control test rig showed that the proposed fuzzy logic controller used for this purpose was able to provide a whirl orbit of prescribed magnitudes.The work further showed that the whirl orbit radius that the fuzzy logic controller was able to generate was much higher than that obtained using a conventional PD controller.Du et al. [11] has proposed a fuzzy-control strategy based on the Takegi-Sugeno model [12] to stabilize a single-axis magnetic-bearing system with fast-response speed subjected to parameter uncertainties and saturation of the control voltage.Numerical simulations, used to ascertain the effectiveness of this control scheme, showed the robustness of this scheme to uncertainties in the active magnetic-bearing-model parameters.Shi and Lee [13] compared the performance of a fuzzy logic controller with a model-based controller for the stabilization of an experimental magnetic-bearing setup.They found that the fuzzy logic controller was more robust and had a better performance in rejecting disturbances compared to the model-based controller.Noshadi et al. [14] proposed time-domain objective functions, which were used to tune the design parameters of a PID-type fuzzy.They further performed experiments on a laboratory-based active magnetic-bearing system and showed that the PID-type fuzzy controllers, which maintained lower profiles of control signals, had better performance as compared to the linear on-board controllers.Nguyen et al. [15] reported work on a fuzzy logic controller for magnetic-bearing application.The magnetic-bearing system's mathematical model was the basis for the fuzzy controller's membership functions and rule design.The proposed control algorithm was numerically evaluated, and the outcome revealed its effectiveness in controlling the response of the rotor.Specifically, the controller demonstrated its effectiveness in suppressing steady-state and transient responses of the rotor.A novel design procedure for a neuro-fuzzy controller of a magnetically supported rotor based on optimal robust control was proposed by Carvalho et al. [16].The tuning of this controller considered the optimal balance between vibration attenuation and robustness, which was carried out using the robust optimization procedure.Evaluation of the unbalance response of the rotor revealed the effectiveness of the proposed procedure.
For nonlinear systems, robust control techniques, such as the SMC, have also been used for magnetic-bearing applications.This control scheme accommodates the variability and un-modeled dynamics of the system.In this control scheme, the control problem of a multi-state control system is reduced to a control problem of a first-order control system [17].For a simple control structure, the performance of this scheme is satisfactory.This control scheme has been utilized by Yeh et al. [18] and Rundell et al. [19] to address issues of nonlinear vibrations in magnetic-bearing systems.Maslen and Montie [20] have contended that there is no added advantage of using the sliding mode scheme over the the PID controller for vibration control in the magnetic-bearing system.On the other hand, Allaire and Sinha [21] have shown the robustness of SMC schemes to plant uncertainities.Tian et al. [22] have further demonstrated the robustness of this control scheme to address disturbances in the flexible rotor model.For first-order systems, Foust et al. [23] showed that the SMC scheme is equivalent to the bang-bang control scheme.The control activities are regulated by the SMC through a bang-bang controller to obtain time-optimal control.By switching between the two opposite extremes, the bang-bang controller provides minimumtime control [24].Yao and Chen [25] proposed a deep convolutional neural network-an SMC scheme for the trajectory control of a rotor mounted on magnetic bearings.The sliding mode scheme was used to obtain the tracking control with high robustness.The deep convolutional neural network, on the other hand, was used for the compensation of the system's uncertainties.Numerical simulation undertaken on a five-degree-of-freedom rotor-magnetic-bearing system model showed the effectiveness of the proposed control scheme.Magnetic-bearing controllers based on the sliding mode scheme and neural network algorithm was also presented by Cao et al. [26].The SMC scheme was utilized to increase the robustness and response speed of the system.The neural network algorithm, on the other hand, provided online tuning for compensation of unmodeled uncertainties and external disturbances.Numerical simulation was carried out to examine the effectiveness of the proposed control scheme against the PID controller.Numerical results showed superior performance of the proposed control scheme over the PID controller.In order to overcome the problem of high-frequency chattering of the conventional SMC scheme, which affects the operational stability of rotors mounted in active magnetic bearings, Wang et al. [27] proposed a second-order SMC algorithm.Lyapunov criterion [28] was used to prove the stability of the proposed SMC scheme.Numerical simulation results revealed that the high-frequency chattering was absent in the response of the rotor with the use of the proposed control scheme.Experimental results from a test rig that employed this control scheme further showed that the stability time of the magnetic bearing was decreased by 50%, and the axis trajectory range was reduced by 41.5% as compared to the use of the SMC with exponential-reaching law.
Past numerical and experimental works on the control of magnetic bearings using fuzzy-based and SMC schemes have shown the effectiveness of these schemes in controlling non-linear vibrations in the response of rotors supported by magnetic bearings.These works further showed the superior performance of these control schemes over the conventionally used PID controllers.In an earlier paper [29], the authors have shown the effectiveness of a fuzzy bang-bang relay controller (FBBRC) for the control of a single-axis active magnetic-bearing system.A two-level bang-bang crisp output is directly produced by the FBBRC.Consequently, it does not need further saturation or any hard-limiting devices.This output is based on the largest of maxima (LOM) defuzzification method.This controller is also structurally simple since its fuzzy output set and rule matrix need only two membership functions.The present paper numerically examines the effectiveness of the same controller for the control of vibration of a rigid rotor supported by two-axes active magnetic bearings.The mathematical model of the nonlinear magnetic-bearing forces considers the coil current and air gap relationship, as well as the cross-coupling due to the geometry of the magnetic bearing.The performance of the FBBRC is compared with the conventional fuzzy logic controller (FLC) and the linear PD controller for controlling the vibration of the rotor-bearing system, particularly, for severe operations due to large rotor unbalance force and high rotor angular speed.
The subsequent sections of this paper are organized as follows.Section 2 presents the formulation of the governing equations of the rigid rotor-magnetic-bearing system.The design of the fuzzy controller is described in Section 3. Results of the numerical simulation are presented and discussed in Section 4. Finally, Section 5 presents the conclusions.

Governing Equations of the Rotor-Bearing System
Several assumptions were made in the derivation of the equations of motion for the rigid rotor mounted on active magnetic bearings: (i) the rotor is symmetric, (ii) the rotor unbalance is static and located at the middle of the bearing span, (iii) the axial motion of the rotor and gyroscopic effect are neglected, (iv) the magnetic flux leakage and fringing effect are neglected, and (v) the operation of the magnetic iron is below the level of saturation.
Figure 1 shows the schematics of an active magnetic bearing and a rigid rotor supported by two active magnetic bearings.The coordinates, X and Y, represent the directions perpendicular to the rotor's axis of rotation.The coordinate, Z, on the other hand, is along the direction of the rotor's axis of rotation.The displacement of the geometric center of the rotor in the X and Y directions are represented by x and y, respectively.The equations of motion of the rotor's geometric center, accounting for the magnetic-bearing forces, rotor unbalance force, and gravity, are given by [30].
The rotor's half-mass is denoted by m , the eccentricity of the rotor's center of mass is denoted by ρ , the rotor's angular speed is denoted by ω , and the gravity constant is represented by g .The magnetic bearing forces in the X and Y directions are respec- tively represented by x F and y F .
The electromagnetic force, F, generated by a single-acting magnetic actuator is presented in Equation (2).The permeability of free space is denoted by 0 μ , the number of magnetic actuator coil turns is denoted by N, and the area of one magnetic pole is represented by g A .The coil current and the length of one air gap are respectively denoted by i and 0 In the operation of these bearings, the differential driving mode configuration is utilized, where a pair of actuators are arranged to counteract each other thus enabling the generation of positive and negative forces.).
The details on the derivation of the magnetic-bearing forces in the X and Y directions for the differential-driving-mode configuration of the magnetic actuators can be found in [32].For a nominal air gap, 0 g , the current in the decreasing air gap side, for the magnetic force in the X direction, and the current in the increasing air gap side, ( The rotor's half-mass is denoted by m, the eccentricity of the rotor's center of mass is denoted by ρ, the rotor's angular speed is denoted by ω, and the gravity constant is represented by g.The magnetic bearing forces in the X and Y directions are respectively represented by F x and F y . The electromagnetic force, F, generated by a single-acting magnetic actuator is presented in Equation ( 2).The permeability of free space is denoted by µ 0 , the number of magnetic actuator coil turns is denoted by N, and the area of one magnetic pole is represented by A g .The coil current and the length of one air gap are respectively denoted by i and g 0 [31].
In the operation of these bearings, the differential driving mode configuration is utilized, where a pair of actuators are arranged to counteract each other thus enabling the generation of positive and negative forces.This arrangement requires the bias flux, B b , and the perturbation flux, B p , which are respectivley induced by the bias steady-state current, i b , and the control current, i p , in the coil.The superposition of B b and B p yields the total flux.The maximum dynamic range of the magnetic actuator is obtained by setting i b to about half the saturation current, with one actuator driven by (i b + i p ) and the other by The details on the derivation of the magnetic-bearing forces in the X and Y directions for the differential-driving-mode configuration of the magnetic actuators can be found in [32].For a nominal air gap, g 0 , the current in the decreasing air gap side, (g 0 − x), is (i b + i x ) for the magnetic force in the X direction, and the current in the increasing air gap side, (g 0 + x), is (i b − i x ).Substituting these values into Equation (2) yields the nonlinear magnetic-bearing forces for the differential-driving-mode configuration, where the positive and negative forces in the X direction are respectively denoted by F x+ and F x− Energies 2022, 15, 3975 6 of 20 The positive and negative forces in the Y direction, F y+ and F y− , can be similarly derived to yield Equation (4).
Accounting for the effect of the geometrical coordinate coupling, α, which couples the governing equations in the X and Y directions, the magnetic-bearing forces in both these directions can be expressed by Equation ( 5) [29].The geometrical curvature of the magnetic pole causes the air gap to vary over the entire pole area, and this, in turn, results in the occurrence of a normal force, which is 90 degrees away from the principal force during motion in the principal axis.The ratio of the principal force to the normal force is known as the geometrical coordinate coupling, α [33].
Inserting Equation ( 5) into Equation ( 1) gives the governing equation of a rigid rotor supported by active magnetic bearings.
The proportional, K P , and derivative, K D , feedback gains for a linear PD controller are related to i x , i y , x, y, .
x and .y by Equation ( 7) [32].The present work makes use of the system parameters presented by Delprete et al. [34].These values are reproduced in Table 1.The geometrical coupling coordinate, which was not considered in [34], has a value of 0.2.This value was selected based on the work of Knight et al. [35], who had presented a range of geometrical coupling coordinate values between 0.05 to 0.25.

Fuzzy Controller Design
In this work, the FBBRC and the conventional FLC are used to control the motion of the rotor in magnetic bearings.The FBBRC is made up of fuzzy logic and a hard limiter relay combined in one structure.The same fuzzy set is given as input to both controllers.However, the controllers give different output fuzzy sets.The LOM defuzzification technique used in the FBBRC results in a bang-bang output [36].The centroid defuzzification technique is used by the FLC.

Controller Fuzzy Set
For the control of the active magnetic bearing system in both the X − Y axes, two identical fuzzy controllers are utilized, i.e., one in each direction.Hence, construction details are discussed for one controller only.In the design of these fuzzy controllers, ranges of their input and output variables need to be determined since they reasonably represent all the operating conditions that the rotor-bearing system controllers may encounter to yield stable operation of the system.These limiting values are derived in the following sections.

Fuzzy Membership Sets
The nominal air gap, g 0 , of the active magnetic-bearing system, given in Table 1, is 0.5 mm.For the control in the X axis, the displacement x(t), which is measured in the positive direction from the equilibrium position, x(t) = 0, is indexed as input i = 1.The error position is e(t) = x d (t) − x(t), and the range of the universe of discourse is χ 1 ∈ [−0.01, 0.01] m, which is much larger than g 0 to cater for large offsets from the operating conditions.x d (t) = ω is the desired input frequency.Numerical simulations of the model are performed in closed loop for different initial displacements (x 0 , y 0 ), eccentricity of the rotor's center of mass, ρ, and angular speed, ω.Two extremes of output level are needed for the control of the rotor-bearing system's motion with bang-bang action.Independent on/off switching of i x± or i y± produces the bang-bang action required to keep the rotor in its equilibrium position.The choice of applied forces are determined by the control law.The actuation of F x± or F y± is achieved from Equations ( 3) and ( 4), respectively, by switching on/off the current, i x (t).The required bang-bang action, which sets the output universe of discourse to I bb = [−1, +1] A for the bang-bang controller, Figure 3, is obtained by switching the i x (t) between ±1 A for the case of i b (t) = 0.5 A. A small dead band is introduced between the bang-bang action to reduce chattering, which can potentially excite high frequencies in the rotor-bearing system model.The output of the FLC, on the other hand, is shown in Figure 4.The linguistic set A k i defines the k th membership function of i th input variables, x i=1 = 'error position' and x i=2 = 'error velocity'.The error position and error velocity are subsequently defined over the universe of discourse χ i .For input x 1 , the following linguistic values: , see Figure 5.
are the on/offfiring command for actuator forces 1 F and 2 F , respectively.

Fuzzy Rules
The rotor is reset to its equilibrium position using the fuzzy linguistic rules proposed herein.Two input variables, each having three linguistic values, form the basis of these rules.Therefore, there are at most nine possible rules.Table 2 describes these rules in matrix form.Wong et al. [38] presented a direct Lyapunov's stability method for individual firing rules, which form the basis of the rules given in Table 2.The diagonal entries in the rules of Table 2 represent the dead band and act as dampeners to slow down the switching action.The bang-bang action is produced by the off-diagonal entries in Table 2, which represent the inputs' unequal strength.The smooth control of the rotor over the universes of discourse are made possible by heuristically choosing the rules partitions.The symmetry of the system's dynamics give rise to the symmetry of the rules matrix.Equation ( 8) gives the decomposition of linguistic rules from the inputs of the FBBRC to its output.The th j fuzzy rule that contributes to the output is presented by , in this equation.

Fuzzy Rules
The rotor is reset to its equilibrium position using the fuzzy linguistic rules proposed herein.Two input variables, each having three linguistic values, form the basis of these rules.Therefore, there are at most nine possible rules.Table 2 describes these rules in matrix form.Wong et al. [38] presented a direct Lyapunov's stability method for individual firing rules, which form the basis of the rules given in Table 2.The diagonal entries in the rules of Table 2 represent the dead band and act as dampeners to slow down the switching action.The bang-bang action is produced by the off-diagonal entries in Table 2, which represent the inputs' unequal strength.The smooth control of the rotor over the universes of discourse are made possible by heuristically choosing the rules partitions.The symmetry of the system's dynamics give rise to the symmetry of the rules matrix.Equation ( 8) gives the decomposition of linguistic rules from the inputs of the FBBRC to its output.The j th fuzzy rule that contributes to the output is presented by the index, j = 1, 2, . . ., J, in this equation.
The standard decomposition and centroid defuzzification technique are used by the conventional FLC [33].The shape of the input and output membership functions for the conventional FLC are the same as in Figures 2 and 3, but the defuzzification method is 'centroid', and the rules are shown in Table 3.
Table 3. Fuzzy rules used in the conventional FLC.P = positive, Z = 0, N = negative.The formulation of the fuzzy bang-bang relay controller is based on the output membership functions, shown in Figure 3, and the LOM aggregation, depicted in Figure 5.The rule matrix in Table 2 defines the output membership functions, which are acted upon by any disturbance of the rotor from its equilibrium position.The maximum value of the degree of membership function, µ B j k , shown in Figure 3, governs the overall output of the FBBRC.By taking the maximum of aggregation given in Equation ( 9), µ ol , which represents the membership function of the overall implied fuzzy rules, j = 1, 2, . . ., J, is determined.
Equation (10) gives the expression for the defuzzified crisp output, u, which is related to µ ol in Equation (9).
The LOM value, which occurs at the output universe of the discourse's extremes, I bb = [−I x , I x ], is represented by the supremum in Equation (10).The argument, arg sup [µ ol (i x )], in Equation ( 10) returns u = [−1, +1]. Figure 3b presents the bang-bang firing action of the membership functions, I2 and I1.

PD Controller
The performance of FBBRC, proposed in this work, was compared with the performance of the PD controller.The PD controller feedback gains values, K P and K D , are based on the work reported by Delprete et al. [34], which are reproduced in Table 1.Numerical simulations of the rotor-bearing system with the linear PD controller were undertaken using the same values of initial conditions, rotor unbalance parameter, and rotor angular speed, as those used for the FLC and FBBRC controllers.

Numerical Results and Discussion
Numerical simulations were performed using MATLAB and Simulink Release 2021b (The MathWorks, Inc., Natick, MA, USA) computer package to determine the strengths and drawbacks of the three controllers, PD, FLC, and FBBRC, when used to control the vibration of the rotor-bearing system.The values of the parameters used in these simulations are given in Table 1.For the control of the rotor-bearing system using the linear PD scheme, the fuzzy controller blocks were replaced by PD controller blocks.Since the magnetic bearing has a nominal air gap of 0.50 mm, numerical simulations were terminated when the displacements, x or y, reached this value, and the simulations were reported to have failed.The simulations were run for a duration equal to 1500 cycles of the rotor's motion based on the angular speed, ω, and the time increment was chosen to be 0.001 s.
In the first set of simulations, the initial conditions were chosen to be the same along both axes, i.e., x 0 = y 0 , and .
x 0 = .y 0 = 0.The system was then simulated for different values of angular speed ω and eccentricity of the rotor center of mass ρ chosen over the range 50 to 500 rad/s in steps of 50, and 0.03 to 0.08 mm in steps of 0.01, respectively.It was found that none of the controllers was able to bring the rotor to its equilibrium position (x = y = 0) when (x 0 = y 0 ≥ 0.30 mm).In other words, the simulation failed with x or y having exceeded 0.50 mm sometime during the simulation.For x 0 = y 0 = 0.20 mm, the PD controller failed for all combinations of (ω,ρ) while the FLC failed for all combinations of (ω,ρ) except ω = 50 rad/s, ρ = 0.03, 0.04 mm.However, the FBBRC gave better results with the simulations being successful for ω = 50, 100 rad/s and all ρ, for ω = 150 rad/s, ρ =0.03, 0.04, 0.05, 0.06 mm, and for ω = 200 rad/s, ρ = 0.03 mm.
The response of the rotor-bearing system for 0.01 ≤ x 0 = y 0 ≤ 0.10 mm and with the above values of ω and ρ is shown in Figure 6.A value of 10 along the vertical axis implies that the system simulation for all 10 values of initial conditions is successful.We refer to this value as the number of convergent (initial) conditions.A value less than 10 implies that the simulation fails for one of the initial conditions.It can be seen from Figure 6a that the presence of the PD controller improves the system performance for all initial conditions when ω is small (50 and 100 rad/s).The presence of the FLC controller makes the system behavior erratic with simulations being successful only for small ω (50 rad/s) and small ρ (0.03 and 0.04 mm), Figure 6b.The FBBRC controller gives superior performance for a larger range of (ω, ρ) pairs, as can be seen from the flat region at the top of Figure 6c.Figures 7 and 8 show the last 10 cycles of the system responses for two different (ω, ρ) pairs and x 0 = y 0 = 0.05 mm.These responses are typical for initial conditions, x 0 = y 0 , that were simulated.The left column shows the time response for x (time response for y is similar), and the right column shows the rotor response on the X − Y plane.The X − Y plane axes have been intentionally selected to represent the air gap length of 0.50 mm.
Figure 7 shows successful simulations when a PD controller or an FBBRC controller is used.The FLC controller fails for this case.Note that the response to a PD controller contains a single frequency, and hence, the rotor position makes a perfect circle.However, the response due to the FBBRC contains multiple frequencies, and hence, the rotor response forms a non-uniform shape.Figure 8 shows the same responses for higher ω and ρ values.Here, both the PD controller and FLC responses failed while the FBBRC produced a bounded response, albeit with a non-uniform rotor response.
Figures 9 and 10 simulate the system for a fixed (ω, ρ) pair but different initial displacements.The marks on the graph show the initial positions of the rotor (x 0 , y 0 ) for which the system was simulated.The '×' represents a simulation that failed, while an 'o' represents a successful simulation.The initial positions were chosen on radial circles with the equilibrium point (x 0 = y 0 = 0) as the center, and radii 0.01 to 0.10 mm in steps of 0.01 mm and 0.20, 0.30, 0.40 mm.Angular increment with 15 • step totaling 24 points was chosen on each circle.As expected, the simulations failed for large initial displacements in most cases for all three controllers.For ρ = 0.08 mm and ω = 100 rad/s, Figure 9, the PD controller response was successful for all initial displacements up to 0.10 mm.The behavior of the FLC controller was, as seen earlier, erratic with the '×' and 'o' not showing any discernable pattern.The FBBRC showed better performance with successful simulations up to 0.20 mm.For ρ = 0.06 mm and ω = 250 rad/s, Figure 10, the PD controller response showed an interesting pattern.The simulation was successful for initial displacements up to 0.10 mm only in the first quadrant and part of the second quadrant.The erratic behavior of the FLC responses did not come as a surprise.However, the FBBRC showed a truly superior behavior with the responses being successful for all initial displacements up to 0.10 mm, most initial displacements of 0.20 mm, and a few initial displacements of 0.30 mm.It is noted that the robustness of the proposed FBBRC was not addressed in the present work.The robustness of the proposed bang-bang control scheme to external disturbances and variations of the rigid rotor-magnetic-bearing system parameters can be achieved with the introduction of SMC to the control scheme.SMC offers the most well-known robust strategy, which works under discontinuous on/off (bang-bang) operation [39].This controller can handle external disturbances and variations in the system parameters to give robust performance.To obtain time-optimal control, the SMC directs the control activity through bang-bang controllers, which are entrusted for minimum time response.This robust and minimum time property of the control scheme for the rigid rotor-magneticbearing system can be achieved with the combination of SMC and bang-bang controller.It is noted that the robustness of the proposed FBBRC was not addressed in the present work.The robustness of the proposed bang-bang control scheme to external disturbances and variations of the rigid rotor-magnetic-bearing system parameters can be achieved with the introduction of SMC to the control scheme.SMC offers the most well-

Conclusions
The fuzzy bang-bang relay controller, which integrates into one entity the fuzzy logic intelligence and the hard-limiting relay, has been employed for the control of a magnetically supported rigid rotor.This controller is inherently stable because of its bang-bang action properties.Its configuration is also less complex as opposed to that of the conventional fuzzy controller as its output membership function only uses two choices, which results in a relatively simpler rule matrix.Numerical simulations of the rotor-bearing system were undertaken for a range of rotor imbalance magnitude and angular speed to evaluate the performance of the FBBRC.A comparison of the FBBRC with the conventional FLC and PD controller for similar operating conditions of the rotor-bearing system demonstrated its superiority, in particular at extreme operating conditions.The improved performance of the FBBRC, however, comes at the expense of the non-synchronous rotor response.The considerably small amplitude of the non-synchronous response is nevertheless within the acceptable limit of vibrations for rotating machinery that falls under the category of rigid rotor.
The present work can be extended to the control of flexible rotors supported by magnetic bearings.Unlike the rigid rotors considered in the present work, the dynamics of flexible rotors are more complex since the flexibility of the shaft introduces not only two additional degrees of freedom but also a bending critical speed and associated mode shape.For the case of a symmetrical flexible rotor, the magnetic-bearing controller should be able to control the vibration at the mid-span of the rotor and at the location of the magnetic bearing.

Figure 1 .
Figure 1.(a) Schematic of an active magnetic bearing.(b) Schematic of a rigid rotor supported by two active magnetic bearings.
This arrangement requires the bias flux, b B , and the perturbation flux, p B , which are respectivley induced by the bias steady-state current, b i , and the control current, p i , in the coil.The superposition of b B and p B yields the total flux.The maximum dynamic range of the magnetic actuator is obtained by setting b i to about half the saturation current, with one actuator driven by (

Figure 1 .
Figure 1.(a) Schematic of an active magnetic bearing.(b) Schematic of a rigid rotor supported by two active magnetic bearings.

Figure 2
shows the input membership function of the fuzzy set for the FLC and FBBRC.The numbers 1, 2, and 3 in Figure 2 represent the three different input membership functions of the fuzzy set.The smaller values of the universe of discourse for χ 1 , shown in Figure 2, ensures lower sensitivity of the controller near the dead band zone [37].The error velocity, .e(t) = .x d (t) − .x(t), is indexed as i = 1, and the range of the universe of discourse is χ 2 ∈ [−10, 10] mm/s.The error velocity is given more derivative action for fast response to the tracking frequencies.Energies 2022, 15, x FOR PEER REVIEW 8 mm/s.The error velocity is given more derivative action for fast response to the tracking frequencies.

Figure 2 .Figure 2 .
Figure 2. Input membership function of the fuzzy set for FBBRC and FLC.Two extremes of output level are needed for the control of the rotor-bearing system's motion with bang-bang action.Independent on/off switching of ± x i or ± y i produces the bang-bang action required to keep the rotor in its equilibrium position.The choice of applied forces are determined by the control law.The actuation of ± x F or

Figure 4 .=
Figure 4. Output of FLC.The linguistic set k i A defines the th k membership function of th i input variables,

1 =
P] can adequately define the control level of the system operation.Similar linguistic values are selected for input x 2 , A k 2 ≡ A k 1 .The linguistic set B k defines the output membership function, y f is the individual rules output, when defuzzified with LOM gives the aggregated out membership functions

Figure 5 .
I1 = +i x (t) and I2 = −i x (t) are the on/off-firing command for actuator forces F 1 and F 2 , respectively.Energies 2022, 15, x FOR PEER REVIEW 11 of 22 The linguistic set k B defines the output membership function, f y is the individual rules output, when defuzzified with LOM gives the aggregated out membership functions ]

Figure 7 .
Figure 7.Last 10 cycles of rotor response for

Figure 8 .
Figure 8.Last 10 cycles of rotor response for

Figures 9
Figures 9 and 10 simulate the system for a fixed (  ,  ) pair but different initial displacements.The marks on the graph show the initial positions of the rotor ( 0 x , 0 y ) for

Table 1 .
Parameters for the numerical simulation model.

Table 2 .
Fuzzy rules used in FBBRC.