Abstract

We propose a novel combined five-degrees-of-freedom (5-DOFs) hybrid magnetic bearing (HMB) with only one permanent magnet ring (PMR) used for turboexpanders. It has two radial magnetic bearing (RMB) units; each has four poles and one thrust magnetic bearing (TMB) to control 5-DOFs. Based on one PMR, the bias flux of the two radial magnetic bearing units and the one thrust magnetic bearing unit is constructed. As a result, ultra-high-speed, lower power loss, small size, and low cost can be achieved. Furthermore, the equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the combined 5-DOFs HMB. The force-current, force-position, torque-coil currents, the torque-angle position, and the stiffness models of the combined 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed in this paper. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.

1. Introduction

Magnetic bearing can provide many advantages for turbomachinery such as high-speed, long-life, higher efficiency, lower maintenance, and higher reliability. However, the greatest benefit for turboexpanders is realized by the ability of the magnetic bearings to run in the process gas, which in turn eliminates the lube oil and dry gas seals. As a result, there is no lube oil to contaminate the process. Elimination of the lubrication system also provides significant size and mass savings for off-shore applications. The turboexpanders in S2M, GE, SKF, and so forth have been equipped with the magnetic bearings [15]. Jumonville et al. provide useful photographs of actual magnetic bearing hardware, loading plots, and general design guidelines. A neon cryogenic turboexpander with magnetic bearing is developed by Hirai et al. [6]. The experimental equipment of active magnetic bearing on a 150 M3 turbo-oxygen gas expander is tested by Xiping et al. [7]. Active magnetic bearing applied on a gas expander is mentioned in the literature [8, 9]. A magnetically levitated 400 kW turbogenerator system for natural gas expansion is developed in the literature [10]. Schmied provides excellent information on magnetic bearing control system in a very understandable fashion [11]. The rotor of the turboexpander mentioned above is supported by active magnetic bearings in 5 degrees-of-freedom (DOFs) (Figure 1), which consists of two radial magnetic bearings with bias current, two thrust magnetic bearings with bias current, expander wheel, displacement sensors, shaft, touchdown bearings, and compressor wheel. The radial magnetic bearing and thrust magnetic bearing with bias current have better control ability and high stiffness, whereas they suffer from high power loss due to the bias current. To reduce the power consumption, hotness, and nonlinearity, the radial and axial hybrid magnetic bearings with permanent magnet providing bias flux for suspension and the electromagnet providing necessary control flux are developed [1215]. Most of the 5-axis hybrid magnetic bearing (HMB) systems are usually composed of two radial hybrid magnetic bearings (RMBs) and one pair of axial hybrid magnetic bearings (AMBs) [1618]. These magnetic bearing systems are easy to manufacture, but their disadvantages are high power consumption, bulky, complex process, and high cost. The reason is that each RMB has eight poles that lead to more eddy current loss. Moreover, each of RMBs and AMBs has one permanent magnet ring (PMR) that is used to provide bias flux. As many as four PMRs are required in this 5-DOFs HMB system.

In order to reduce the power consumption, cost, and size, it is important to reduce the number of the units by means of combining the RMB and AMB. Integrated axial magnetic bearing and radial magnetic bearing with the conical rotor are designed and analyzed [19, 20] and one downside of the integrated bearings is a strong coupling problem between the radial and axial degrees-of-freedom. In addition, the combined radial-axial hybrid magnetic bearings (CRAMBs) are proposed [2128]. The CRAMBs mentioned above have been applied in ultra-high-speed machinery, but it may not suit the turboexpanders.

In this paper, a novel combined 5-DOFs HMB with only one PMR is proposed to be used for turboexpanders in this paper. The combined 5-DOFs HMB consists of two RMB units each with four poles and one AMB unit. The bias flux of the two RMB units and AMB unit is supplied by only one PMR. The low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB units. The low cost is gained due to the one PMR used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk. The equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the combined 5-DOFs HMB. The force-current, force-position, torque-current, the torque-angular position, and the stiffness models of the integral 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.

This paper is organized as follows. First, in Section 2, the configuration of combined 5-DOFs HMB is given. Then, the equivalent magnetic circuit method is used to model and analyze the combined 5-DOFs HMB in Section 3. Section 4 is devoted to developing simulation and clarifying the linearized models. Finally, Section 5 concludes this paper. We focus on the configuration, modeling, and analysis of the combined 5-DOFs HMB; a future paper will compare with experimental results.

2. Analyzed Configuration

A simplified configuration of turboexpander supported by the combined 5-DOFs HMB is shown in Figure 2. The assembly consists of one 5-DOFs HMB, expander wheel, displacement sensor, shaft, touchdown bearing, and compressor wheel. And the rigid body motion on the rotor-bearing system is considered in this paper. The configuration and magnetic circuit of the combined 5-DOFs HMB are shown in Figure 3. The force and torque produced by the RMB units and AMB unit are shown in Figure 4. It consists of two RMB units (RMB units and ), one pair of AMB unit, and one axially magnetized PMR. The bias fluxes produced by the permanent magnets are shown by the solid arrow lines and the control flux is generated by -axis (blue color) and -axis coil currents (red color), which are shown by the dotted arrow lines in Figure 3. Each stator of the RMB unit consists of two soft magnetic pieces each of which has four equally spaced magnetic poles and four coils distributed at intervals of 90 degrees around the stator iron. The control flux produced by the radial and axial coil currents will allow adding or subtracting flux to the radial air gaps and axial air gaps, and then restoring force and torque will be produced to stabilize the rotor in 5-DOFs except the spin axis in relation to the rotor position. The bias flux of the RMB units and AMB unit is produced by only one PMR. The RMB units and stabilize rotor on the two axes (-axis and -axis) orthogonal to the spin axis (-axis) and exert restoring torque on the rotor around -axis or (and) -axis when subjected to tilting disturbed torque. The AMB unit exerts restoring force on the rotor at -axis. In a word, the low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB unit. The low cost is gained due to only one PMR being used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk.

3. Modeling and Analysis of the Combined 5-DOFs HMB

3.1. The Equivalent Magnetic Circuits of the Combined 5-DOFs HMB

According to construction of the combined 5-DOFs HMB, the magnetic resistance of the iron path and flux leakage is neglected, the permeability of the iron is assumed to be infinite, and their equivalent bias and control magnetic circuits can be gained and are shown in Figure 5.

3.2. The Model of Combined 5-DOFs HMB for Translation in -Axis (or -Axis) and -Axis Direction

Figure 3 describes the operating principle of this combined 5-DOFs HMB in -axis and -axis. The rotor will depart from its center position when subject disturbance force in negative -axis direction. The control flux on the upper side of the RMB units and will add to the bias flux in upper radial air gap and subtract from the bias flux in the lower radial air gap. And the control flux of the AMB unit should also be calculated. A net restoring force on the rotor will be produced in positive -axis direction.

According to the structure of the combined 5-DOFs HMB (Figure 4) and its equivalent magnetic circuits (Figure 5), the net force of the RMB unit A in the -axis direction (, where is the number of winding turns and is the control current of -axis coil currents), , is where the constant is the permeability of free space and is equal to 4 × 10−7 H/m; is the area of the pole face of RMB unit; and are the bias flux in radial air gaps of the RMB unit ; and are the control flux in the radial air gaps of the RMB unit produced by the -axis coil currents; and are the control flux in the radial air gaps of the RMB unit produced by the -axis coil current of the AMB unit.

The bias flux in Figure 5 can be calculated, respectively, as follows: where

As shown in Figures 3 and 5, is the total reluctance of the bias circuit; and are the reluctance of the radial air gaps in -axis for the RMB unit ; and are the reluctance of the radial air gaps in -axis for the RMB unit ; and are the reluctance of the radial air gaps in -axis for the RMB unit ; and are the reluctance of the radial air gaps in -axis for the RMB unit ; and are the reluctance of the axial air gaps in -axis for the AMB unit. , where is the coercive force of the permanent magnet; is PMR thickness; is the magnetic resistance of PMR and it can be given as where is the relative permeability of permanent magnet and is the middle-section area of the PMR.

The reluctances of radial air gap are written by where is the nominal air gap of the RMB unit; is the rotor position as deviation from the operating point in -axis direction; is the deviation from the operating point in -axis direction.

The reluctances of axial air gap are written by where is the nominal air gap of the AMB unit; is the rotor position as deviation from the operating point in -axis direction; is the area of the pole face for the AMB unit.

The control flux in the radial air gaps of -axis, and , generated by -axis coil currents of RMB unit (, where is the number of winding turns and is the control current of -axis coil currents) can be calculated as where is the control current of -axis coil currents for the same displacement of RMB unit .

The control flux in the radial air gaps, and , which are produced by the -axis coil currents of the AMB unit can be calculated as where the control flux in the axial air gap is calculated by where is the number of winding turns of the AMB and is the control current of -axis coil current. And .

Substituting (2), (7), and (8) into the right side of (1) results in the net force on -axis direction. Using the Taylor series expansion for small values of and , we can get the following attractive force with linear terms in -axis direction:

Since the net force is zero at the center position in radial and axial directions, (10) can be written by where is the current stiffness of the RMB unit and is the negative position stiffness of the RMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to magnetic bearing force (MBF) is in the form of hyperbolic functions indicating that grows with time. They can be calculated as follows:

According to the structure of the combined 5-DOFs HMB (Figure 4) and its equivalent magnetic circuits (Figure 5), the net force of the AMB unit in the -axis direction is where and are the bias flux in axial air gaps of the AMB unit and is the control flux in the axial air gaps of the AMB unit by the -axis coil current.

Substituting (2) and (9) for the right side of (13) results in the net force on -axis direction. Using the Taylor series expansion for small values of and , we can get the following attractive force with linear terms in -axis direction:

Since the net force is zero at the center, (14) can be written by where is the current stiffness of the AMB unit and is the negative position stiffness of the AMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to MBF is in the form of hyperbolic functions indicating that grows with time. They can be calculated as follows:

3.3. The Model of Combined 5-DOFs HMB for Tilting

The RMB units and of the combined 5-DOFs HMB also control 2-DOFs around -axis and -axis which are orthogonal to the spin axis (-axis) when the rotor is subject to disturbance torque around -axis and -axis. The RMB units and create tilting stiffness to restrain the rotor when subjected to the disturbed torques as shown in Figure 4. According to the structure of the combined 5-DOFs HMB (Figure 3) and its equivalent magnetic circuits (Figure 5), the torque, , acting on the rotor can be given by (17) ( and ). Consider where is the distance from center of the magnetic pole of the RMB units and to the center of rotor as shown in Figure 4. The for small value of the tilting angle, . Equation (17) can be simplified as

The change in radial air gap is for small value of the tilting angle, . The radial air gap resistance can be written as

Using the Taylor series expansion for small values of   and we can get the following attractive torque with linear terms: where is the tilting current stiffness and is the negative tilting displacement stiffness. They can be calculated as follows:

3.4. The Coupling Model of the Combined 5-DOFs HMB
3.4.1. The Current Stiffness and Position of the RMB Unit Affected by the Axial Air Gap

The current stiffness and position stiffness of the RMB unit in (12) do not consider the effect of the rotor in -axis direction, namely, the axial air gap. However, the bias flux generated from the PMR flows across the radial air gaps and the axial air gaps (Figure 5(a)); the alternation of the axial air gap will affect the current stiffness and the position stiffness of the RMB and will be calculated by

3.4.2. The Current Stiffness and Displacement Stiffness of the RMB Unit Affected by the Coil Current in AMB Unit

From Figure 5(d), the control flux generated by the coil current in AMB unit will flow through radial air gaps of the RMB unit, and the current stiffness and displacement stiffness of the RMB unit could be affected and given by

3.4.3. The Current Stiffness and Displacement Stiffness of the AMB Unit Affected by the Radial Air Gap

The alternation of the radial air gap will affect the current stiffness and the displacement stiffness of the AMB unit and will be calculated by

4. The Performance Analysis of the Combined 5-DOFs HMB Using Field-Circuit Method

4.1. The Combined 5-DOFs HMB Parameters and the 3D Finite Element Model

The combined 5-DOFs HMB parameters are given in Table 1. The example aims to validate the design and analysis method based on linearized model and 3D finite element method (FEM).

The linearized model is given by where and are the current stiffness of the system in -axis and -axis, respectively; is the displacement stiffness of the system in -axis and -axis; is the displacement stiffness in -axis; and are the tilting displacement stiffness around -axis and -axis; and are the tilting current stiffness around -axis and -axis. Moreover, the rotor can be controlled in five degrees-of-freedom.

The finite element model using Ansoft Maxwell V13 is built and shown in Figure 6(a) except the air. Its meshes are shown in Figures 6(b) and 6(c), and the total number of mesh is 286036.

4.2. The Performance Analysis for Radial and Axial Translation of the Combined 5-DOFs HMB
4.2.1. The Analysis Results of the Force-Displacement and Displacement Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm in -axis (or -axis), and hold in center position in -axial direction, the relationship between net force and the rotor position in -axis (or -axis) direction is shown in Figure 7(a) using linearized method, nonlinear method, and FEM. It is obvious that the force is nonlinear. The calculated results are approximately equal for small values of the rotor in -axis (or -axis) and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffnesses using the three calculated methods are shown in Figure 7(b). The displacement stiffness in rotor center is using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is ; the calculated error between them is 2.6%.

4.2.2. The Analysis Results of the Force-Coil Current and Current Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current in -axis or-axis of the RMB unit changes from −1 A to 1 A (the coil current of the AMB unit is 0 A), the relationship between net force of the RMB unit and coil current is shown in Figure 8 using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the RMB unit using the three calculated methods is 284.4 N/A (linearized method and nonlinear method) and 284.13 N/A (FEM); the calculated error is between linearized method and FEM 1.5%.

4.2.3. The Analysis Results of the Force-Coil Current and Rotor Position in -Axis (or -Axis) of the RMB Unit (Using Nonlinear Method)

When the rotor in -axis locates its center position and the coil current of the AMB unit is 0 A, the curves of the net force-coil current-rotor position of the RMB unit in -axis (or -axis) using nonlinear method are shown in Figure 9. It is obvious that it is approximately linearity for small displacement and coil current and nonlinear for large displacement and coil current. The displacement stiffness and current stiffness affected by the rotor position and the coil current in -axis (or -axis) of the RMB unit are shown in Figures 10 and 11, respectively. They are approximately linear for small displacement and coil current.

4.2.4. The Analysis Results of the Force-Displacement and Displacement Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current all are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm in -axis, and hold in center position in -axial and -axis direction, the relationship between net force and the rotor position in -axis direction is shown in Figure 12(a) using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the force is also nonlinear. The calculated results are approximately equal for small values of the rotor in -axis and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffness using the three calculated methods is shown in Figure 12(b). The calculated results are approximately equal for small value of the displacement in -axis between the linearized method and FEM. The displacement stiffness in rotor center is using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is  N/m; the calculated error between them is 4.4%.

4.2.5. The Analysis Results of the Force-Coil Current and Current Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current in -axis of the AMB unit changes from −1 A to 1 A (the coil currents of the RMB unit are 0 A), the relationship between net force of AMB unit and its coil current is shown in Figure 13 using linearized method (see (15)), nonlinear method (see (13)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the AMB unit using the three calculated methods is 290.5 N/A (linearized method and nonlinear method) and 280.3 N/A (FEM); the calculated error is 3.5% between linearized method and FEM.

4.2.6. The Analysis Results of the Force-Coil Current and Rotor Position in -Axis of the AMB Unit (Using Nonlinear Method)

When the rotor locates its center position in -axis and -axis and the coil currents of the RMB units are 0 A, the curves of the net force-coil current-rotor position of the AMB unit in -axis using nonlinear method are shown in Figure 14. It is obvious that the relationship of the net force-displacement in -axis is approximately linear and the relationship of the net force-versus coil current is approximately linear for small displacement and nonlinear for large coil current. The displacement stiffness affected by the rotor position and the coil current in -axis of the AMB unit is shown in Figure 15(a). The current stiffness affected by the rotor position and the coil current in -axis of the AMB unit is shown in Figure 15(b). In addition, the current stiffness is only affected by the coil current of the AMB unit.

4.3. The Performance Analysis for Tilting of the Combined 5-DOFs HMB
4.3.1. The Analysis Results of the Torque-Angular Position and Angular Displacement Stiffness of the Combined 5-DOFs HMB

The torque can be produced by the RMB units and . When the control currents in radial coil and axial coil all are 0 A and the rotor tilts within −5 mrad~5 mrad around -axis or -axis, the relationship between torque and angular position of rotor is shown in Figure 16(a) using linearized method, nonlinear method, and FEM. The calculated results are approximately equal for small values of the angular position of rotor and there is little difference between the nonlinear model and FEM. In addition, the tilting position stiffness is shown in Figure 16(b). The calculated results are approximately equal for small value of the angular position between the linearized method and FEM. The angular displacement stiffness in rotor center is −5255.6 Nm/rad.

4.3.2. The Analysis Results of the Torque-Coil Current and Tilting Current Stiffness of the Combined 5-DOFs HMB (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil currents of the RMB units and in -axis or -axis change from −1 A to 1 A (the coil current of AMB unit is 0 A), the relationship between torque and coil current of the RMB units is shown in Figure 17 using linearized method (see (15)), nonlinear method (see (13)) and FEM. It is obvious that the torque is linear, and the calculated results are approximately equal. In addition, the tilting current stiffness using the three calculated methods is 39.25 Nm/A (linearized method and nonlinear method) and 38.07 Nm/A (FEM); the calculated error between the linearized method and FEM is 3%.

4.3.3. The Analysis Results of the Torque-Coil Current and Angular Position in -Axis or -Axis of the Combined 5-DOFs HMB (Using Linearized Method and Nonlinear Method)

When the rotor locates its center position in -axis and the coil currents of the AMB unit are 0 A, the curves of the torque produced by the RMB units and versus their coil current and rotor angular position around -axis or -axis using nonlinear method are shown in Figure 18 using linearized method and nonlinear method. It is obvious that it is approximately linear for small values of the angular position and coil current. Moreover, the curves of the angular position stiffness versus angular position and coil currents in RMB units and are shown in Figure 19(a). And the curves of the tilting current stiffness versus angular position and coil currents in RMB units and are shown in Figure 19(b). They all are affected by the angular position and coil currents in RMB units and .

4.4. The Coupling Performance Analysis of the Combined 5-DOFs HMB
4.4.1. The Coupling Performance Analysis for Translation of the Combined 5-DOFs HMB

Equation (22) shows that the displacement stiffness and current stiffness of the RMB unit will be affected by the rotor position in -axis and the coil current of the AMB unit. The effect of the rotor position in -axis on the displacement stiffness and current stiffness of the RMB units and is shown in Figures 20(a) and 20(b), respectively. They all are affected by the rotor position in -axis and their change trend is opposite.

The effect of the coil current of the AMB unit on the displacement stiffness and current stiffness of the RMB units and is shown in Figures 21(a) and 21(b), respectively. They all are affected by the coil current of the AMB unit and their change trend is the same.

The 3D curves of the displacement stiffness and current stiffness of the RMB units and affected by the rotor position in -axis and the coil current of the AMB unit are shown in Figures 22(a) and 22(b).

Figure 23 shows that the displacement stiffness and current stiffness of the AMB unit will be affected by the rotor position in -axis or -axis of the RMB unit. The effect of the rotor position in -axis or -axis on the displacement stiffness and current stiffness of the AMB unit is shown in Figures 23(a) and 23(b), respectively. They all are affected by the rotor position in -axis or -axis.

4.4.2. The Coupling Performance Analysis for Tilting of the Combined 5-DOFs HMB

The tilting displacement stiffness and tilting current stiffness of the RMB units and will be affected by the rotor position in -axis and coil current of the AMB unit, and the influence curves of RMB units and are shown in Figures 24(a), 24(b), 25(a), and 25(b), respectively. It is obvious that the tilting displacement stiffness and tilting current stiffness of the RMB units and all are affected by the rotor position in -axis and coil current of the AMB unit.

5. Conclusion

A novel combined 5-DOFs HMB with only one PMR is proposed, used for turboexpander, in this paper. The combined 5-DOFs HMB consists of two RMB units and one AMB unit. The bias flux of the two RMB units and AMB unit is supplied by only one PMR. The low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB. The low cost is gained due to the one PMR used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk. The equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the integral 5-DOFs HMB. The force-current, force-position, torque-coil currents, torque-angular position, and the stiffness models of the integral 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.

Conflict of Interests

The authors declare no conflict of interests.

Acknowledgments

This work was supported by the Aviation Science Fund of China under Grant 2012ZB51019 and by the Cultivation and Development Project of Science and Technology Innovation Base of Beijing under Grant Z131104002813105.