Thermal Behavior of a Magnetically Levitated Spindle for Fatigue Testing of Fiber Reinforced Plastic

This article discusses the critical thermal behavior of a magnetically levitated spindle for fatigue testing of cylinders made of fiber reinforced plastic. These cylinders represent the outer-rotor of a kinetic energy storage. The system operates under vacuum conditions. Hence, even small power losses in the rotor can lead to a high rotor temperature. To find the most effective way to keep the rotor temperature under a critical limit in the existing system, first, transient electromagnetic finite element simulations are evaluated for the active magnetic bearings and the electric machine. Using these simulations, the power losses of the active components in the rotor can be derived. Second, a finite element simulation characterizes the thermal behavior of the rotor. Using the power losses calculated in the electromagnetic simulation, the thermal simulation provides the temperature of the rotor. These results are compared with measurements from an experimental spindle. One effective way to reduce rotational losses without major changes in the hardware is to reduce the bias current of the magnetic bearings. Since this also changes the characteristics of the magnetic bearings, the dynamic behavior of the rotor is also considered.


Introduction
Flywheels store energy as kinetic energy of the rotor and can provide a cost efficient solution for short-term energy storage and load smoothening services in electricity grids (e.g., [1,2]).Their advantages lie in the high possible number of cycles and the low initial cost per unit of power.The main drawbacks are their high stand-by-losses and the relatively low energy density compared to other storage technologies.In order to utilize the advantages, the energy density of one system should be increased while decreasing the stand-by losses at the same time.

Outer-Rotor Fywheel Design
One possible flywheel design is an outer-rotor setup (e.g., [3][4][5]).It promises high energy densities due to the large radii of the rotor and high rotational speeds.A realized full-scale system is described in [6].A model of the system is shown in Figure 1a.The rotor is a magnetically levitated hollow cylinder.The outside of the rotor is made out of fiber reinforced plastic (FRP) with a circumferential fiber orientation, using carbon fibers.The inside of the rotor consists of different rotor parts of the active and passive components.Active magnetic bearings (AMBs) are used for radial levitation and a passive magnetic bearing is used for axial levitation.A permanent magnet synchronous machine (PMSM) accelerates and decelerates the rotor.All rotating components of the magnetic bearings and the PMSM are integrated on the inner circumference of the FRP rotor.Under rotation, these components press against the FRP, resulting in radial compressive stress transversal to the fiber orientation, which is superimposed by circumferential stress in the fiber direction (see Figure 1b).The energy density of the system increases with the radii of the rotor and its rotational speed, both factors lead to an increased stress in the FRP [7].Consequently, high stress in the material is necessary to reach high energy densities.Charging and discharging the flywheel leads to cyclically varying mechanical stresses.To investigate the cyclic stability and lifetime of the FRP rotor, cyclic material tests, such as cyclic transverse tensile tests (derived from [8]), cyclic transverse compressive tests (derived from [9]), and cyclic four-point bending tests are performed on material samples.These samples are thin walled and relatively easy to fabricate with a high quality, whereas the rotor of the flywheel is much thicker and harder to produce in an industrial fiber winding process.Furthermore, the thermal expansion of carbon fiber and the epoxy plastic matrix used differ, leading to inner stress when a thick walled FRP structure cools down from the curing temperature.Hence, there is a much higher probability of imperfections and defects in the thick walled FRP structure of an outer-rotor flywheel than in the thin material samples.Therefore, the applicability of test results obtained with thin material samples to the rotor of the flywheel has to be investigated by tests with thick walled rotors.A specialized test rig was designed and set up to perform cyclic tests on thick-walled FRP specimens, which represent the outer-rotor of kinetic energy storage.These specimens are smaller, cheaper and have a lower energy content, reducing the danger during destructive testing.
hollow cylinder.The outside of the rotor is made out of fiber reinforced plastic (FRP) with a circumferential fiber orientation, using carbon fibers.The inside of the rotor consists of different rotor parts of the active and passive components.Active magnetic bearings (AMBs) are used for radial levitation and a passive magnetic bearing is used for axial levitation.A permanent magnet synchronous machine (PMSM) accelerates and decelerates the rotor.All rotating components of the magnetic bearings and the PMSM are integrated on the inner circumference of the FRP rotor.Under rotation, these components press against the FRP, resulting in radial compressive stress transversal to the fiber orientation, which is superimposed by circumferential stress in the fiber direction (see Figure 1b).The energy density of the system increases with the radii of the rotor and its rotational speed, both factors lead to an increased stress in the FRP [7].Consequently, high stress in the material is necessary to reach high energy densities.Charging and discharging the flywheel leads to cyclically varying mechanical stresses.To investigate the cyclic stability and lifetime of the FRP rotor, cyclic material tests, such as cyclic transverse tensile tests (derived from [8]), cyclic transverse compressive tests (derived from [9]), and cyclic four-point bending tests are performed on material samples.These samples are thin walled and relatively easy to fabricate with a high quality, whereas the rotor of the flywheel is much thicker and harder to produce in an industrial fiber winding process.Furthermore, the thermal expansion of carbon fiber and the epoxy plastic matrix used differ, leading to inner stress when a thick walled FRP structure cools down from the curing temperature.Hence, there is a much higher probability of imperfections and defects in the thick walled FRP structure of an outer-rotor flywheel than in the thin material samples.Therefore, the applicability of test results obtained with thin material samples to the rotor of the flywheel has to be investigated by tests with thick walled rotors.A specialized test rig was designed and set up to perform cyclic tests on thick-walled FRP specimens, which represent the outer-rotor of kinetic energy storage.These specimens are smaller, cheaper and have a lower energy content, reducing the danger during destructive testing.

Testing Procedure
The ratio between the circumferential and the radial stress in a test specimen should be as close as possible to the one present in a full-scale flywheel.In the system described in [6], the radial transverse compressive stress in the FRP is 60 MPa and the circumferential longitudinal stress is 389 MPa at its maximum speed of 15,000 rpm.In the design of the flywheel, the maximum stress was limited to half of the expected strength of the FRP, to account for the high data uncertainty.To test the limits of the FRP, the stress in the specimen is doubled compared to the flywheel.To create this

Testing Procedure
The ratio between the circumferential and the radial stress in a test specimen should be as close as possible to the one present in a full-scale flywheel.In the system described in [6], the radial transverse compressive stress in the FRP is 60 MPa and the circumferential longitudinal stress is 389 MPa at its maximum speed of 15,000 rpm.In the design of the flywheel, the maximum stress was limited to half of the expected strength of the FRP, to account for the high data uncertainty.To test the limits of the FRP, the stress in the specimen is doubled compared to the flywheel.To create this state of stress in the test specimen, a circumferentially segmented steel ring is placed inside the FRP cylinder of the specimen, and rotated at 30,000 rpm.This results in a transverse compressive stress of 100 MPa and a longitudinal tension of 775 MPa in the FRP.With an outer diameter of 190 mm, the surface speed of the specimen reaches approximately 298 m/s.The fatigue test will be conducted by alternating between a maximum speed of 30,000 rpm and a minimum speed of 15,000 rpm.At minimum speed the transversal compressive stress is 26 MPa and the longitudinal tension is 200 MPa.Because of the well-known high tensile strength of carbon fiber, the chance of fiber fracture is low at this state of stress.The probability of matrix fracture is expected to be much higher [10].It is planned to perform up to 200,000 cycles per specimen.With a cycle time of 30 s this will take about 70 days.Additionally, overload tests are planned, where the specimen is accelerated to 40,000 rpm, resulting in a surface speed of nearly 398 m/s, a circumferential longitudinal tension of 1,400 MPa and a radial compressive stress of 180 MPa.The latter should lead to failure of the matrix.In the longitudinal direction, the FRP can withstand static tensions over 2,000 MPa.To reduce air drag and the consequential heating of the specimen, all tests are performed under vacuum conditions, with a pressure below 0.05 Pa.

Test Rig Description
Figure 2a shows a section view of the designed test rig, which was introduced in [11].A hub and a shaft coupling connected the specimen to the driving spindle.This configuration, rather than an outer-rotor setup, was chosen to protect the active parts of the spindle in case of a failure of the specimen at high speed.Each of the eight segments of the steel ring inside the FRP weighed 1.027 kg, and the distance between their center of gravity (CoG) and the rotation axis is about 55 mm.Hence, at 40,000 rpm one segment contained around 30 kJ of kinetic energy.The surrounding containment was designed to absorb this energy in case of specimen failure.The design criteria for the containment were derived from [12].The containment also serves as a vacuum chamber.A detailed view of the spindle is shown in Figure 2b.The motor in the middle is a water cooled PMSM with four poles, a maximum torque of 9.6 Nm, and a maximum power of 30 kW.In order to avoid excessive wear of the high-speed drive, the rotor is supported by AMBs.The radial AMBs were designed in a heteropolar configuration with laminated cores on rotor and stator, both made out of 0.2 mm thick sheets of non-oriented silicon steel (NO20).The radial position of the rotor is measured with four eddy current sensors at each bearing, two for each radial direction.Because of the well-known high tensile strength of carbon fiber, the chance of fiber fracture is low at this state of stress.The probability of matrix fracture is expected to be much higher [10].It is planned to perform up to 200,000 cycles per specimen.With a cycle time of 30 s this will take about 70 days.Additionally, overload tests are planned, where the specimen is accelerated to 40,000 rpm, resulting in a surface speed of nearly 398 m/s, a circumferential longitudinal tension of 1,400 MPa and a radial compressive stress of 180 MPa.The latter should lead to failure of the matrix.In the longitudinal direction, the FRP can withstand static tensions over 2,000 MPa.To reduce air drag and the consequential heating of the specimen, all tests are performed under vacuum conditions, with a pressure below 0.05 Pa.

Test Rig Description
Figure 2a shows a section view of the designed test rig, which was introduced in [11].A hub and a shaft coupling connected the specimen to the driving spindle.This configuration, rather than an outer-rotor setup, was chosen to protect the active parts of the spindle in case of a failure of the specimen at high speed.Each of the eight segments of the steel ring inside the FRP weighed 1.027 kg, and the distance between their center of gravity (CoG) and the rotation axis is about 55 mm.Hence, at 40,000 rpm one segment contained around 30 kJ of kinetic energy.The surrounding containment was designed to absorb this energy in case of specimen failure.The design criteria for the containment were derived from [12].The containment also serves as a vacuum chamber.A detailed view of the spindle is shown in Figure 2b.The motor in the middle is a water cooled PMSM with four poles, a maximum torque of 9.6 Nm, and a maximum power of 30 kW.In order to avoid excessive wear of the high-speed drive, the rotor is supported by AMBs.The radial AMBs were designed in a heteropolar configuration with laminated cores on rotor and stator, both made out of 0.2 mm thick sheets of non-oriented silicon steel (NO20).The radial position of the rotor is measured with four eddy current sensors at each bearing, two for each radial direction.To lift the combined mass of the spindle rotor and the specimen of roughly 20 kg, the axial AMB needs a large pole area, which leads to a big outer diameter of the thrust disk.Because of the high stress at high speed, a monolithic design was chosen (see [13]).The stator of the axial AMB was fabricated of the soft magnetic composite Siron S400b of the PMG Füssen GmbH [14].An inductive sensor at the upper end of the rotor detects the axial position.Two-rowed hybrid spindle bearings with a static load rating of 20 kN are used as backup bearings.The rotor temperature is measured with two infrared sensors-one between the upper radial AMB and the thrust disk and one next to the lower radial position sensors, below the lower AMB.The moment of inertia of the specimen is about ten times higher than that of the rotor, and the surface speed of the specimen is more than two times higher.To protect the rotor and the spindle in case of a specimen failure, a predetermined breaking point was included in the hub (see Figure 3).This tapering also uncouples the dynamics of the rotor from the specimen to some extent.To lift the combined mass of the spindle rotor and the specimen of roughly 20 kg, the axial AMB needs a large pole area, which leads to a big outer diameter of the thrust disk.Because of the high stress at high speed, a monolithic design was chosen (see [13]).The stator of the axial AMB was fabricated of the soft magnetic composite Siron S400b of the PMG Füssen GmbH [14].An inductive sensor at the upper end of the rotor detects the axial position.Two-rowed hybrid spindle bearings with a static load rating of 20 kN are used as backup bearings.The rotor temperature is measured with two infrared sensors-one between the upper radial AMB and the thrust disk and one next to the lower radial position sensors, below the lower AMB.The moment of inertia of the specimen is about ten times higher than that of the rotor, and the surface speed of the specimen is more than two times higher.To protect the rotor and the spindle in case of a specimen failure, a predetermined breaking point was included in the hub (see Figure 3).This tapering also uncouples the dynamics of the rotor from the specimen to some extent.One goal of the test rig operation is to minimize the total test duration for one specimen, maximizing the number of specimens that can be tested in a given amount of time.The test rig should thus be operated with as much power as possible.High power normally leads to high losses, as well as subsequent heating.Losses on the stator can be effectively cooled with water, but without a medium for convection or thermal conduction through physical contact, the rotor can only transfer heat to the stator via radiation.Hence, even small power losses in the rotor can lead to a high rotor temperature.The temperature of the rotor is crucial for the feasibility of the cyclic fatigue testing, since the temperature of the specimen affects the strength of the FRP and can change the test result.Furthermore, the magnets of the PMSM should not be operated over a certain temperature.Thus, the overall goal of this paper is to minimize the test duration under the boundary condition of a maximum rotor and specimen temperature.For this purpose, the losses in the rotor during operation and their influence on the rotor temperature have to be identified.Furthermore, feasible methods of minimizing the most crucial losses in an existing test rig are derived.

Method
With the goal of finding the shortest total test duration, this paper proposes a structured approach to evaluate the heating of a magnetically levitated rotor in a vacuum, and deduces measures to avoid critical temperatures.An overview on the approach is shown in Figure 4. First, losses in the rotor of both radial AMBs, the axial AMB, and the PMSM were calculated for different rotational speeds, using transient electromagnetic 2D finite element (FE) models implemented in ANSYS Maxwell 2015.The highest and lowest operating speeds of the spindle followed from the demanded stress in the FRP, as described in Section 1.2.For each component, a polynomial was fitted to the related losses, to obtain a mean loss value, which was used for the calculation of the steady state One goal of the test rig operation is to minimize the total test duration for one specimen, maximizing the number of specimens that can be tested in a given amount of time.The test rig should thus be operated with as much power as possible.High power normally leads to high losses, as well as subsequent heating.Losses on the stator can be effectively cooled with water, but without a medium for convection or thermal conduction through physical contact, the rotor can only transfer heat to the stator via radiation.Hence, even small power losses in the rotor can lead to a high rotor temperature.The temperature of the rotor is crucial for the feasibility of the cyclic fatigue testing, since the temperature of the specimen affects the strength of the FRP and can change the test result.Furthermore, the magnets of the PMSM should not be operated over a certain temperature.Thus, the overall goal of this paper is to minimize the test duration under the boundary condition of a maximum rotor and specimen temperature.For this purpose, the losses in the rotor during operation and their influence on the rotor temperature have to be identified.Furthermore, feasible methods of minimizing the most crucial losses in an existing test rig are derived.

Method
With the goal of finding the shortest total test duration, this paper proposes a structured approach to evaluate the heating of a magnetically levitated rotor in a vacuum, and deduces measures to avoid critical temperatures.An overview on the approach is shown in Figure 4. First, losses in the rotor of both radial AMBs, the axial AMB, and the PMSM were calculated for different rotational speeds, using transient electromagnetic 2D finite element (FE) models implemented in ANSYS Maxwell 2015.The highest and lowest operating speeds of the spindle followed from the demanded stress in the FRP, as described in Section 1.2.For each component, a polynomial was fitted to the related losses, to obtain a mean loss value, which was used for the calculation of the steady state temperature.Next, air friction losses were calculated analytically, based on the kinetic theory of gases.Finally, to calculate the rotor temperature, a 3D-FE-model of a quarter of the test rig was set up and evaluated using ANSYS Workbench 19.1.Because of the water cooling, losses on the stator were not considered; instead the stator temperature in the simulation was predefined with measured data.A transient evaluation of the model was compared to measured temperatures on the test rig, to adjust thermal and loss calculation.To compute the rotor temperature after a long time of cycling, a steady state evaluation of the adjusted model was used.From this model, the most influential components on the rotor temperature can be derived.The focus was on components that allow for an improvement of the thermal behavior without major changes in the hardware.One option, which will be discussed as an example, is the reduction of the bias current in the radial AMBs [15].Since this also changed the behavior of the rotor in the AMBs, a brief analysis of the rotor dynamics was performed using measurements and a mechanical 3D-FE-model of the rotor implemented in ANSYS Workbench 19.1.Another loss reduction can be achieved by reducing the control activity of the AMB [15].Most of these strategies, like unbalance compensation (e.g., [16]) or a linear-quadratic-Gaussian-control (e.g., [17]), were not feasible at this point due to limitations of the AMB controller hardware.If no further loss reduction can be achieved in the components, the next step is to adjust the cycle time.Since rotational losses increase with speed and losses in the PMSM also increase with acceleration, the cycle time has a major influence on the mean losses and therefore on the rotor temperature.If the rotor temperature is below the critical values, the cycle time can be reduced.If the critical values are exceeded, the cycle time has to be increased.These options will be discussed further on.
Actuators 2019, 8, x FOR PEER REVIEW 5 of 18 temperature.Next, air friction losses were calculated analytically, based on the kinetic theory of gases.Finally, to calculate the rotor temperature, a 3D-FE-model of a quarter of the test rig was set up and evaluated using ANSYS Workbench 19.1.Because of the water cooling, losses on the stator were not considered; instead the stator temperature in the simulation was predefined with measured data.A transient evaluation of the model was compared to measured temperatures on the test rig, to adjust thermal and loss calculation.To compute the rotor temperature after a long time of cycling, a steady state evaluation of the adjusted model was used.From this model, the most influential components on the rotor temperature can be derived.The focus was on components that allow for an improvement of the thermal behavior without major changes in the hardware.One option, which will be discussed as an example, is the reduction of the bias current in the radial AMBs [15].Since this also changed the behavior of the rotor in the AMBs, a brief analysis of the rotor dynamics was performed using measurements and a mechanical 3D-FE-model of the rotor implemented in ANSYS Workbench 19.1.Another loss reduction can be achieved by reducing the control activity of the AMB [15].Most of these strategies, like unbalance compensation (e.g., [16]) or a linear-quadratic-Gaussiancontrol (e.g., [17]), were not feasible at this point due to limitations of the AMB controller hardware.
If no further loss reduction can be achieved in the components, the next step is to adjust the cycle time.Since rotational losses increase with speed and losses in the PMSM also increase with acceleration, the cycle time has a major influence on the mean losses and therefore on the rotor temperature.If the rotor temperature is below the critical values, the cycle time can be reduced.If the critical values are exceeded, the cycle time has to be increased.These options will be discussed further on.

Modeling of the Thermal Stability of the Rotor
The following sections discuss the loss calculation of the active magnetic bearings for the radial and axial direction and the electric machine, as well as losses caused by air friction.Subsequently, the thermal model focusing on the rotor temperature is described.For the model validation and potential adjustment, the calculation results are compared with measurements.

Loss Calculation of the Radial Active Magnetic Bearings
Both radial AMBs have eight poles, which are evenly distributed on the inner circumference of the stator and arranged in the sequence N-S-S-N-N-S-S-N.Flux leakage between two poles with equal flux direction-that is, N-N or S-S-is small, and thus two neighboring poles with differing flux direction, called a pole pair, function as one independent electromagnet.One pole pair of the upper radial AMB with an idealized average path of the magnetic flux φ is shown in Figure 5a.The magnetic flux of the AMB passes through the NO20-metal sheets that surrounded the solid rotor, which is shown in the bottom left corner.The nominal air gap between the rotor and the surrounding stator in the centered position is 0.4 mm.To linearize the actuator characteristic, a differential winding

Modeling of the Thermal Stability of the Rotor
The following sections discuss the loss calculation of the active magnetic bearings for the radial and axial direction and the electric machine, as well as losses caused by air friction.Subsequently, the thermal model focusing on the rotor temperature is described.For the model validation and potential adjustment, the calculation results are compared with measurements.

Loss Calculation of the Radial Active Magnetic Bearings
Both radial AMBs have eight poles, which are evenly distributed on the inner circumference of the stator and arranged in the sequence N-S-S-N-N-S-S-N.Flux leakage between two poles with equal flux direction-that is, N-N or S-S-is small, and thus two neighboring poles with differing flux direction, called a pole pair, function as one independent electromagnet.One pole pair of the upper radial AMB with an idealized average path of the magnetic flux ϕ is shown in Figure 5a.The magnetic flux of the AMB passes through the NO20-metal sheets that surrounded the solid rotor, which is shown in the bottom left corner.The nominal air gap between the rotor and the surrounding stator in the centered position is 0.4 mm.To linearize the actuator characteristic, a differential winding design is used, where a common bias coil and a control coil, with associated currents, are used for two counteracting pole pairs [18].In the test rig, two coils are mounted on each pole of the radial AMBs, the inner for bias current I B and the outer for the control current I Y .Both coils have the same number of turns.All bias coils are connected in series for one AMB and the control coils for each direction in one AMB.The magnetic flux in a specific point in the rotor changes direction each time it passed a pole with a different flux direction.The flux also drops between two poles with equal flux direction.This remagnetization leads to hysteresis losses P h and excess losses P e , as well as eddy currents and subsequent losses P c .To calculate the losses, a transient electromagnetic 2D-FE-model of the rotor was set up and evaluated.Figure 5b The corresponding loss constants  ,  , and  were derived by fitting data from the manufacturer of the electric steel obtained by standardized measurements [20] to Equation (1), which yields: T 2 m 3 and  , 193.6 Ws T 2 m 3 .For the NO20 sheets used in the radial AMBs,  , could not be significantly determined and therefore was set to zero.The calculated total iron losses for the upper and lower radial AMB are shown in Table 1.Additionally, the switching of the amplifiers leads to further high frequency changes in the magnetic flux.The resulting rotor losses can be calculated with the same model and yielded 0.11 W for the upper AMB and 0.09 W for the lower AMB.For the loss calculation, the Bertotti formula [19] in the form of Equation (1) was applied, in which the three loss mechanisms P c , P h , and P e depend on the amplitude of the flux density B m and the frequency f with which the magnetic flux changes: m dV, and P e = k e ( f B m ) 1.5 dV. ( The corresponding loss constants k c , k h , and k e were derived by fitting data from the manufacturer of the electric steel obtained by standardized measurements [20] to Equation (1), which yields: For the NO20 sheets used in the radial AMBs, k e,AMB could not be significantly determined and therefore was set to zero.The calculated total iron losses for the upper and lower radial AMB are shown in Table 1.Additionally, the switching of the amplifiers leads to further high frequency changes in the magnetic flux.The resulting rotor losses can be calculated with the same model and yielded 0.11 W for the upper AMB and 0.09 W for the lower AMB.

Loss Calculation the of Axial Active Magnetic Bearing
The magnetic flux in the axial AMB is mostly symmetrical to the rotational axis, so no significant remagnetization losses should occur due to rotation.Nevertheless, the switching of the amplifier leads to changes in the magnetic flux and associated rotor losses.Equation ( 1) is only applicable for thin sheets, but in the solid thrust disk of the axial AMB, without the usage of a laminated core or soft-magnetic composites, eddy current losses dominate.Therefore, the other loss mechanisms were neglected.To calculate eddy current losses a transient electromagnetic 2D-FE-model was analyzed.Figure 6a shows the magnetic flux density B in the axial AMB, with an air gap of 0.4 mm and a coil current of 3.15 A. The resulting force is 200 N, which is required to lift the rotor with the specimen, with a total weight of 20.3 kg.Changes in the magnetic flux in the rotor result in eddy currents and therefore in a non-zero current density J in the rotor.Eddy current losses P c,J were calculated as the integral of the square of J over the rotor volume V, divided by the conductivity of the material σ el (see Equation ( 2)).The conductivity of X14CrMoS17 is σ el,rot = 1.43 × 10 6 1  Ωm .

Loss Calculation the of Axial Active Magnetic Bearing
The magnetic flux in the axial AMB is mostly symmetrical to the rotational axis, so no significant remagnetization losses should occur due to rotation.Nevertheless, the switching of the amplifier leads to changes in the magnetic flux and associated rotor losses.Equation ( 1) is only applicable for thin sheets, but in the solid thrust disk of the axial AMB, without the usage of a laminated core or soft-magnetic composites, eddy current losses dominate.Therefore, the other loss mechanisms were neglected.To calculate eddy current losses a transient electromagnetic 2D-FE-model was analyzed.Figure 6a shows the magnetic flux density B in the axial AMB, with an air gap of 0.4 mm and a coil current of 3.15 A. The resulting force is 200 N, which is required to lift the rotor with the specimen, with a total weight of 20.3 kg.Changes in the magnetic flux in the rotor result in eddy currents and therefore in a non-zero current density  in the rotor.Eddy current losses  , were calculated as the integral of the square of  over the rotor volume , divided by the conductivity of the material  (see Equation ( 2)).The conductivity of X14CrMoS17 is  , 1.43 • 10 .
Since eddy currents oppose their provoking field, high frequency changes of the magnetic field only affect the outer layers of the conducting material.Hence, eddy currents are also limited to a thin region on the surface of the rotor (see Figure 6b).The losses of the axial AMB in the rotor due to the amplifier switching, derived by Equation ( 2

Loss Calculation of the Permanent Magnet Synchronous Machine
The last active component that induces losses in the rotor is the PMSM.Here, the magnetic field rotates synchronously with the rotor.Remagnetization of the rotor occurs because of flux drops at slots in the stator and non-harmonic changes of the motor current due to switching of the inverter.Losses were calculated via a transient electromagnetic 2D-FE model.In the solid permanent magnets, losses were calculated using Equation ( 2), with an electric conductivity of the magnets of  , 1.1 • 10 . The rotor is laminated underneath the permanent magnets to further reduce losses.
Remagnetization losses in these sheets were calculated according to Equation (1).The corresponding coefficients are: T 1.5 m 3 .Results for different phase currents and rotational speeds are shown in Figure 7.The losses increased nearly linearly with the speed above 10,000 rpm.For root mean square values of the phase current (rms) below 10 A, losses were nearly independent of the phase current.Since eddy currents oppose their provoking field, high frequency changes of the magnetic field only affect the outer layers of the conducting material.Hence, eddy currents are also limited to a thin region on the surface of the rotor (see Figure 6b).The losses of the axial AMB in the rotor due to the amplifier switching, derived by Equation (2), did not exceed 0.1 W.

Loss Calculation of the Permanent Magnet Synchronous Machine
The last active component that induces losses in the rotor is the PMSM.Here, the magnetic field rotates synchronously with the rotor.Remagnetization of the rotor occurs because of flux drops at slots in the stator and non-harmonic changes of the motor current due to switching of the inverter.Losses were calculated via a transient electromagnetic 2D-FE model.In the solid permanent magnets, losses were calculated using Equation (2), with an electric conductivity of the magnets of σ el,mag = 1.1 × 10 6 1  Ωm .The rotor is laminated underneath the permanent magnets to further reduce losses.Remagnetization losses in these sheets were calculated according to Equation (1).The corresponding coefficients are:

Air Friction Losses
For the calculation of air friction losses, it has to be evaluated if the gas in the system can be modeled with continuum dynamics.This was done by calculating the mean free path  ̅ of the gas. ̅ is the avarage distance a molecule moves before it collides with another molecule of the gas.For a gas with pressure  , temperature  , and a mean atom diameter  , the mean free path  ̅ is calculated by Equation (3) (see [21]), where  1.381 10 The mean molecule diameter of air is  3.559 10 m [21].During rotation, the maximum pressure in the test rig is  0.05 Pa and the minimum temperature  290 K, resulting in a minimum mean free path of  ̅ 0.144 .Since  ̅ is much bigger than most of the air gaps in the test rig, an air molecule is more likely to collide with the walls of the test rig than with other air molecules.Hence, continuum dynamics are not applicable for this system.Instead, the kinetic theory of gases was utilized, where molecules are modeled as randomly moving elastic spheres.In [22], air friction losses  were derived by means of the momentum excange beween a spinning rotor and air molecules.For an axial disc with an outer radius  and an inner radius  this yields [22]: for a radial cylindrical surface with a radius  and a height ℎ [22]: and for a truncated cone with a height ℎ , a lower radius  and an upper radius  : where  is the rotational speed of the rotor and  is the mean molecule mass of the gas.Hence, air friction losses increase linearly with , quadratically with , and decrease with .For dry air, the mean atom mass is approximately  4.81 10 kg [21].The losses were calculated for each surface of the rotor and the specimen and were subsequently superimposed.For the worst case, with a maximum pressure of  0.05 Pa and a minimum temperature of  290 K, the sums of the calculated losses for the rotor and the specimen at different speeds are listed in Table 2.

Air Friction Losses
For the calculation of air friction losses, it has to be evaluated if the gas in the system can be modeled with continuum dynamics.This was done by calculating the mean free path l of the gas.l is the avarage distance a molecule moves before it collides with another molecule of the gas.For a gas with pressure p, temperature T, and a mean atom diameter d m , the mean free path l is calculated by Equation (3) (see [21]), where k = 1.381 × 10 −23 J K is the Boltzmann constant, The mean molecule diameter of air is d m = 3.559 × 10 −10 m [21].During rotation, the maximum pressure in the test rig is p = 0.05 Pa and the minimum temperature T = 290 K, resulting in a minimum mean free path of l = 0.144 m.Since l is much bigger than most of the air gaps in the test rig, an air molecule is more likely to collide with the walls of the test rig than with other air molecules.Hence, continuum dynamics are not applicable for this system.Instead, the kinetic theory of gases was utilized, where molecules are modeled as randomly moving elastic spheres.In [22], air friction losses P a were derived by means of the momentum excange beween a spinning rotor and air molecules.For an axial disc with an outer radius r o and an inner radius r i this yields [22]: for a radial cylindrical surface with a radius r r and a height h r [22]: and for a truncated cone with a height h c , a lower radius r l and an upper radius r u : where ω is the rotational speed of the rotor and m W is the mean molecule mass of the gas.Hence, air friction losses increase linearly with p, quadratically with ω, and decrease with T. For dry air, the mean atom mass is approximately m W ≈ 4.81 × 10 −26 kg [21].The losses were calculated for each surface of the rotor and the specimen and were subsequently superimposed.For the worst case, with a maximum pressure of p = 0.05 Pa and a minimum temperature of T = 290 K, the sums of the calculated losses for the rotor and the specimen at different speeds are listed in Table 2.

Simulation of Thermal Rotor Behavior
The thermal simulation was performed using a 3D-FE model of a quarter of the test rig, both with and without the specimen.The rotor transmits thermal energy to the stator via radiation, and vice versa.It was assumed that the emissivity and absorptivity of the rotor and stator surfaces do not depend on the wavelength or direction of the radiation.Under this assumption, the emissivity and absorptivity of the surface are equal [23].For real materials, only a part of the radiation that hits a surface is absorbed while the rest is reflected, hence, the coefficient of emission is smaller than one.When considering the heat exchange between two surfaces through radiation, both surfaces have to be taken into account.Consequently, for the simulation of the rotor heating, the inner surface of the stator has to be included in the model.The stator and rotor were painted with lacquer on the surface of the AMBs and the PMSM.The coefficient of emission was assumed to be 0.9 for all painted parts and the FRP.For the unpainted, blank parts of the steel rotor, a coefficient of emission of 0.3 was assumed, and for the aluminum parts of the stator and the specimen a value of 0.05 was assumed [23].The temporal stator temperature profile of the active components was roughly approximated with three measured temperatures, which were linearly interpolated.The temperature values can be found in Appendix A. This transient temperature change has only a small impact on the simulation results; using a constant mean temperature value for all stator components increases the calculated rotor temperature only about 1 • C. For the rest of the stator, a constant temperature of 25 • C was predefined.To obtain the temperature distribution in the rotor, heat conduction in the rotor was also taken into account in the model, with a conductivity of 25 W/Km.Calculating the steady state rotor temperature during the test cycles requires the mean loss value per cycle for each actuator.As a baseline, a cycling time of 30 s was defined-i.e., during one cycle, the rotor constantly accelerates for 15 s from 15,000 rpm to 30,000 rpm, and then constantly decelerates with the same slope back to 15,000 rpm.The rotational speed and the speed dependent losses of the active components are illustrated in Figure 8a.The mean values of all calculated losses are shown in Figure 8b.

Simulation of Thermal Rotor Behavior
The thermal simulation was performed using a 3D-FE model of a quarter of the test rig, both with and without the specimen.The rotor transmits thermal energy to the stator via radiation, and vice versa.It was assumed that the emissivity and absorptivity of the rotor and stator surfaces do not depend on the wavelength or direction of the radiation.Under this assumption, the emissivity and absorptivity of the surface are equal [23].For real materials, only a part of the radiation that hits a surface is absorbed while the rest is reflected, hence, the coefficient of emission is smaller than one.When considering the heat exchange between two surfaces through radiation, both surfaces have to be taken into account.Consequently, for the simulation of the rotor heating, the inner surface of the stator has to be included in the model.The stator and rotor were painted with lacquer on the surface of the AMBs and the PMSM.The coefficient of emission was assumed to be 0.9 for all painted parts and the FRP.For the unpainted, blank parts of the steel rotor, a coefficient of emission of 0.3 was assumed, and for the aluminum parts of the stator and the specimen a value of 0.05 was assumed [23].The temporal stator temperature profile of the active components was roughly approximated with three measured temperatures, which were linearly interpolated.The temperature values can be found in Appendix A. This transient temperature change has only a small impact on the simulation results; using a constant mean temperature value for all stator components increases the calculated rotor temperature only about 1 °C.For the rest of the stator, a constant temperature of 25 °C was predefined.To obtain the temperature distribution in the rotor, heat conduction in the rotor was also taken into account in the model, with a conductivity of 25 W/Km.Calculating the steady state rotor temperature during the test cycles requires the mean loss value per cycle for each actuator.As a baseline, a cycling time of 30 s was defined-i.e., during one cycle, the rotor constantly accelerates for 15 s from 15,000 rpm to 30,000 rpm, and then constantly decelerates with the same slope back to 15,000 rpm.The rotational speed and the speed dependent losses of the active components are illustrated in Figure 8a.The mean values of all calculated losses are shown in Figure 8b.

Temperature Measurements
To validate the described models, the rotor temperature was measured during cycling without a specimen.Two infrared sensors measured the rotor temperature: one between the upper radial AMB and the thrust disk and one next to the lower radial position sensors (see Figure 2b).The utilized

Temperature Measurements
To validate the described models, the rotor temperature was measured during cycling without a specimen.Two infrared sensors measured the rotor temperature: one between the upper radial AMB and the thrust disk and one next to the lower radial position sensors (see Figure 2b).The utilized sensors can measure the temperature of metallic surfaces above 50 • C, with an uncertainty of ±2 • C. First, the cooling of the levitated but not spinning rotor from a known temperature was analyzed.This excluded the rotational losses and the switching losses in the PMSM.Hence, only the small switching losses of the AMBs were present.The results of a transient simulation and the measurement are shown in Figure 9a.With a maximum deviation of 1.85 • C, the simulation results lie in the range of the measurement uncertainty.Consequently, the switching losses of the AMBs and the modeling of the heat transfer are sufficiently accurate.Next, the thermal behavior of the rotor during rotation was analyzed.Before the measurement, the system was turned off for 12 h, to ensure the rotor was cooled down to room temperature.During the measurement, the PMSM continuously accelerated and decelerated the rotor without a specimen from 15,000 rpm to 30,000 rpm.In the transient simulation, the mean losses from Figure 8b were used.The initial speed up to 15,000 rpm, which takes about 15 s, was neglected.Figure 9b compares the measured temperatures at both sensor positions with the simulation.While the rotor temperature at both positions was nearly the same in the simulation, it differed considerably in the measurement.The model underestimated the rotor heating at both measurement positions.At the lower position, the deviation was 5.2 • C, while it was 22.4 • C at the upper position.Due to the good agreement of the rotor cooling, the modeling error is expected to lie in the loss calculation rather than the thermal model.the modeling of the heat transfer are sufficiently accurate.Next, the thermal behavior of the rotor during rotation was analyzed.Before the measurement, the system was turned off for 12 h, to ensure the rotor was cooled down to room temperature.During the measurement, the PMSM continuously accelerated and decelerated the rotor without a specimen from 15,000 rpm to 30,000 rpm.In the transient simulation, the mean losses from Figure 8b were used.The initial speed up to 15,000 rpm, which takes about 15 s, was neglected.Figure 9b compares the measured temperatures at both sensor positions with the simulation.While the rotor temperature at both positions was nearly the same in the simulation, it differed considerably in the measurement.The model underestimated the rotor heating at both measurement positions.At the lower position, the deviation was 5.2 °C, while it was 22.4 °C at the upper position.Due to the good agreement of the rotor cooling, the modeling error is expected to lie in the loss calculation rather than the thermal model.One reason why the losses could be underestimated can be found in the loss coefficients of the NO20 sheets in the AMBs and PMSM.The loss measurements done by the manufacturer that were used to identify the loss coefficients are normally performed on unmachined sheets after a precise annealing process.However, the material behavior can change during milling and heating.For example, a damaged sheet isolation can increase eddy current losses by up to 30% [24] and hysteresis losses can vary more than 50% depending on their heat treatment [25].Because these changes can hardly be quantified, the calculated losses are adjusted to fit the temperature measurement.At the lower AMB, the calculated losses were increased by 10% to fit the temperature at the lower temperature sensor.Since the sheets at the upper and lower AMB were milled and heated similarly, the losses in the upper AMB increased by 10% as well.To further match the simulation with the measured temperature at the upper sensor, 8 W of additional losses in the axial AMB were required in the simulation.For the loss calculation, it was assumed that the axial AMB is perfectly symmetrical One reason why the losses could be underestimated can be found in the loss coefficients of the NO20 sheets in the AMBs and PMSM.The loss measurements done by the manufacturer that were used to identify the loss coefficients are normally performed on unmachined sheets after a precise annealing process.However, the material behavior can change during milling and heating.For example, a damaged sheet isolation can increase eddy current losses by up to 30% [24] and hysteresis losses can vary more than 50% depending on their heat treatment [25].Because these changes can hardly be quantified, the calculated losses are adjusted to fit the temperature measurement.At the lower AMB, the calculated losses were increased by 10% to fit the temperature at the lower temperature sensor.Since the sheets at the upper and lower AMB were milled and heated similarly, the losses in the upper AMB increased by 10% as well.To further match the simulation with the measured temperature at the upper sensor, 8 W of additional losses in the axial AMB were required in the simulation.For the loss calculation, it was assumed that the axial AMB is perfectly symmetrical with respect to the rotational axis.In reality, small circumferential variations in the material of stator and rotor lead to flux changes during rotation and hence to higher rotor losses.Furthermore, control activities in the AMBs were neglected in the simulation, but during the measurement, considerable movement in the axial direction was observed.The adjusted losses are summarized in Table 3.A transient thermal simulation with the adjusted loss values showed good accordance with the measured data (see Figure 10a).With the adjusted loss values, a steady state thermal simulation of the test rig with and without the specimen was performed.Without the specimen, the motor needs an rms phase current of 3.8 A to accelerate the rotor during the described cycle, but with the specimen an rms phase current of about 37 A is required.The losses of the PMSM were considered accordingly in the simulation.The other losses in the rotor were assumed to be equal for both cases.The simulated temperature distribution in the center of the rotor over the axial length is shown in Figure 10b.In the specimen, the temperature between the FRP and the steel segments is displayed.The maximum temperature without a specimen was 176 °C, located at the upper AMB.Even though higher rotor losses occur with the specimen, the rotor temperature was 10 °C lower than without it, as the large surface area of the FRP benefits heat transfer to the surrounding containment.The thin tapering at the hub led to a large temperature drop between rotor and specimen, preventing the temperature of the FRP from exceeding 40 °C.The operation temperature of the FRP in the flywheel must not exceed 80 °C.The magnets of the PMSM should not be operated over 120 °C.In order to have a safety margin, the magnet temperature should not exceed 110 °C.The simulation results, however, show a temperature of the magnets of over 140 °C.Thus, to avoid reducing the cycle time, rotor losses have to be reduced.

Loss Reduction
As seen in the previous section, the highest power losses in the rotor were due to rotational losses in the radial AMBs and PMSM, as well as switching losses in the PMSM, which all a have similar order of magnitude.Losses due to control activities in the axial AMB were slightly lower.Switching in the AMBs and air friction can be neglected.Hence, the main focus should be on the reduction of rotational losses in the radial AMBs and the PMSM as well as switching losses in the PMSM.In this paper only the radial AMBs will be further discussed.With the adjusted loss values, a steady state thermal simulation of the test rig with and without the specimen was performed.Without the specimen, the motor needs an rms phase current of 3.8 A to accelerate the rotor during the described cycle, but with the specimen an rms phase current of about 37 A is required.The losses of the PMSM were considered accordingly in the simulation.The other losses in the rotor were assumed to be equal for both cases.The simulated temperature distribution in the center of the rotor over the axial length is shown in Figure 10b.In the specimen, the temperature between the FRP and the steel segments is displayed.The maximum temperature without a specimen was 176 • C, located at the upper AMB.Even though higher rotor losses occur with the specimen, the rotor temperature was 10 • C lower than without it, as the large surface area of the FRP benefits heat transfer to the surrounding containment.The thin tapering at the hub led to a large temperature drop between rotor and specimen, preventing the temperature of the FRP from exceeding 40 • C. The operation temperature of the FRP in the flywheel must not exceed 80 • C. The magnets of the PMSM should not be operated over 120 • C. In order to have a safety margin, the magnet temperature should not exceed 110 • C. The simulation results, however, show a temperature of the magnets of over 140 • C. Thus, to avoid reducing the cycle time, rotor losses have to be reduced.

Loss Reduction
As seen in the previous section, the highest power losses in the rotor were due to rotational losses in the radial AMBs and PMSM, as well as switching losses in the PMSM, which all a have similar order of magnitude.Losses due to control activities in the axial AMB were slightly lower.Switching in the AMBs and air friction can be neglected.Hence, the main focus should be on the reduction of rotational losses in the radial AMBs and the PMSM as well as switching losses in the PMSM.In this paper only the radial AMBs will be further discussed.

Reducing Losses in the Radial Bearings
One effective way to reduce the rotational losses of an AMB is to reduce its bias current I B .For example, reducing I B from 5.67 A to 4 A reduces the rotational losses in the rotor by more than 50%, to 8.2 W in the upper AMB and 6.8 W in the lower AMB.A thermal simulation and measurement with I B = 4 A was performed, where the rotor was again cycled according to Figure 8a.The results are shown in Figure 11.The transient simulation and the measurement of the test rig without a specimen again show sufficient correspondence (see Figure 11a).The steady state simulation predicted a maximum temperature of 129 • C at the axial AMB and 123 • C at the PMSM without a specimen, as well as 127 • C and 117 • C with a specimen (see Figure 11b).Changing I B also changes the rotor dynamics in the AMBs.To evaluate the possibility of a bias reduction, the dynamics of the levitated rotor will be investigated in the next section.One effective way to reduce the rotational losses of an AMB is to reduce its bias current  .For example, reducing  from 5.67 A to 4 A reduces the rotational losses in the rotor by more than 50%, to 8.2 W in the upper AMB and 6.8 W in the lower AMB.A thermal simulation and measurement with  4 A was performed, where the rotor was again cycled according to Figure 8a.The results are shown in Figure 11.The transient simulation and the measurement of the test rig without a specimen again show sufficient correspondence (see Figure 11a).The steady state simulation predicted a maximum temperature of 129 °C at the axial AMB and 123 °C at the PMSM without a specimen, as well as 127 °C and 117 °C with a specimen (see Figure 11b).Changing  also changes the rotor dynamics in the AMBs.To evaluate the possibility of a bias reduction, the dynamics of the levitated rotor will be investigated in the next section.

Rotor Dynamics
A mechanical 3D-FE model of the rotor was created and evaluated.For a modal analysis of the model, the AMBs were modeled as springs with a constant stiffness of 2.5 × 10 5 N/m.The first bending mode of the rotor without a specimen lay at 1,030 Hz, thus 353 Hz above the highest rotational frequency.The calculated second bending mode was at 2,170 Hz.Measured Campbell diagrams of the rotor without the specimen are shown in Figure 12.They were derived from the signals of the radial position sensors of the AMBs during an acceleration of the rotor from 0 rpm to 40,000 rpm in 110 s.The signals were filtered with a 5 kHz low-pass filter and sampled with 10 kHz.To create the Campbell diagram, the measurement was divided into 4096 sample long windows, and for each window a Fourier transformation was performed.Figure 12a

Rotor Dynamics
A mechanical 3D-FE model of the rotor was created and evaluated.For a modal analysis of the model, the AMBs were modeled as springs with a constant stiffness of 2.5 × 10 5 N/m.The first bending mode of the rotor without a specimen lay at 1,030 Hz, thus 353 Hz above the highest rotational frequency.The calculated second bending mode was at 2,170 Hz.Measured Campbell diagrams of the rotor without the specimen are shown in Figure 12.They were derived from the signals of the radial position sensors of the AMBs during an acceleration of the rotor from 0 rpm to 40,000 rpm in 110 s.The signals were filtered with a 5 kHz low-pass filter and sampled with 10 kHz.To create the Campbell diagram, the measurement was divided into 4096 sample long windows, and for each window a Fourier transformation was performed.Figure 12a shows the Campbell diagram of the position data of the upper AMB with I B = 5.67 A. The lines starting at the lower left corner are harmonics of the rotational speed, the thickest one being the first harmonic.Odd harmonics of the rotational speed are clearly visible, whereas even harmonics are barely detectable.The approximately horizontal lines are eigenfrequencies (EFs) of the system.At about 90 Hz lies the rigid body tilting mode of the rotor in the AMBs.The translational rigid body mode is not visible in the diagram.The first bending mode of the rotor can be seen at 1057 Hz, the second at 2063 Hz.The second splits into a forward and backward mode.While the calculated first bending mode only differed by 27 Hz, the second bending mode was about 100 Hz lower than the simulation predicted.Underneath each Campbell diagram, the mean deflection of the rotor for each window is shown.The deflections were less than 20% of the backup bearing clearance of 200 µm.Peaks can be seen when the rotational speed or its harmonics hit an EF.The first two peaks were due to the rigid body modes.The next peak at 25,000 rpm was caused by the fifth harmonic reaching the forward mode of the second bending EF, and the last peak at 39,000 rpm was caused by the third harmonic reaching the backward mode of the second bending EF.Smaller peaks can also be seen when higher harmonics reach the second bending EF, where every second odd harmonic excited its forward mode and every other odd harmonic excited the backward mode.A notch filter was placed at 1 kHz in the AMB control, so there was no additional excitation of this EF due to harmonics of the rotational speed.Figure 12b shows the same evaluations for I B = 4 A. The harmonics were not influenced by I B and the EFs only changed slightly.However, the amplitude of the deflection differed noticeable.The peak where the first harmonic crossed the second rigid body mode was about 30% higher with a reduced I B , whereas the peaks at higher speeds were now below 20 µm.Because the normal speed range starts at 15,000 rpm, the rigid body mode has to be crossed only once per test, while the peak at 25,000 rpm is reached twice every cycle.Consequently, without a specimen, reducing I B does not only reduce the rotational losses, but also reduces the control activity of the AMB during cycling.
Actuators 2019, 8, x FOR PEER REVIEW 13 of 18 25,000 rpm was caused by the fifth harmonic reaching the forward mode of the second bending EF, and the last peak at 39,000 rpm was caused by the third harmonic reaching the backward mode of the second bending EF.Smaller peaks can also be seen when higher harmonics reach the second bending EF, where every second odd harmonic excited its forward mode and every other odd harmonic excited the backward mode.A notch filter was placed at 1 kHz in the AMB control, so there was no additional excitation of this EF due to harmonics of the rotational speed.Figure 12b shows the same evaluations for  4 A. The harmonics were not influenced by  and the EFs only changed slightly.However, the amplitude of the deflection differed noticeable.The peak where the first harmonic crossed the second rigid body mode was about 30% higher with a reduced  , whereas the peaks at higher speeds were now below 20 μm.Because the normal speed range starts at 15,000 rpm, the rigid body mode has to be crossed only once per test, while the peak at 25,000 rpm is reached twice every cycle.Consequently, without a specimen, reducing  does not only reduce the rotational losses, but also reduces the control activity of the AMB during cycling.Measurements with a specimen were not performed yet, but Figure 13a shows the calculated EFs and Figure 13Error!Reference source not found.bthe simulated Campbell diagram of the rotor with the specimen.The first elastic mode was at 9 Hz, where only the specimen tilts at the tapering.At speeds higher than 9 Hz, the CoG of the specimen did not move translationally.In the second, third, and fourth modes, the specimen was tilting around its CoG and in the fifth mode it does not Measurements with a specimen were not performed yet, but Figure 13a shows the calculated EFs and Figure 13b the simulated Campbell diagram of the rotor with the specimen.The first elastic mode was at 9 Hz, where only the specimen tilts at the tapering.At speeds higher than 9 Hz, the CoG of the specimen did not move translationally.In the second, third, and fourth modes, the specimen was tilting around its CoG and in the fifth mode it does not move at all.The rotor showed elastic deformations beginning with the fourth mode.Due to the high inertia ratio of the specimen, the higher EFs showed a strong dependency on the rotational speed, which can be seen in the Campbell diagram in Figure 13b.The first three EFs were below 250 Hz, where the fatigue testing begins, and have to be passed while starting a test.The fourth EF increased in such a way that it is not reached by speed-synchronous excitations within the operation range.However, due to the close proximity of the fourth mode, the bandwidth of the AMBs has to include this frequency.The fifth EF had a minimum distance of 290 Hz to the rotational speed and therefore it will have to be filtered out from the position signal.For the differential winding design, which is used for the radial AMBs, the control current has to be smaller than  .To allow control currents bigger than  the AMBs would have to be rebuilt for a differential control design and a nonlinear control might be needed (see [15,[26][27][28]).Both would require major changes in the hardware.Hence, with the existing hardware, reducing  also reduces the maximum force the actuators can generate.In the centered position, one AMB can generate a force of 278 N with  5.67 A, but only 144 N with  4 A. For the rotor without a specimen this is unproblematic, because there are nearly no external loads in the test rig and the rotor can be precisely balanced to minimize unbalance forces.However, the initial unbalance of an assembled specimen is big and balancing it precisely is difficult.The force generated by the AMB has to be big enough to counteract the unbalance force while passing the first three EFs.This has to be further investigated.

Increasing the Cycle Time
When no further loss reduction can be achieved, the cycle time has to be changed.In the upper example, even with a bias current of 4 A, the rotor temperature were above 110 °C but below 120 °C, as shown in Section 4.1.If further bias reduction or other loss reduction strategies are not applicable, a less favorable possibility to reduce the rotor temperature is to increase the cycle time.Here, two variations were considered, which are illustrated in Figure 14a.In the first variation, the acceleration of the rotor is reduced to lower load dependent losses in the PMSM, in the second variation, after a cycle, the rotational speed is held constant at 15,000 rpm for a while to allow the rotor to cool down.How long the cycle has to be increased depends on the associated loss reduction and the maximum allowable rotor temperature.For example, the cycle time could be increased from 30 s to 37 s to reduce the rms motor phase current from 37 A to 30 A, and thereby the rotor losses in the PMSM by 1.7 W to 17.4 W. Losses in the AMBs stay the same.Alternatively, after a 30 s cycle, the speed is held constant at 15,000 rpm for 7 s.Here, the mean losses can be calculated by averaging the mean losses for the first 30 s, as shown in Sections 3.5 and 3.6, and the losses at a constant speed.To keep the rotor without a specimen at a constant speed of 15,000 rpm, an rms phase current of 1.8 A is needed in the PMSM.Hence, it was assumed that an rms phase current of less than 10 A is needed to keep the rotor with a specimen at constant speed.Below 10 A, losses were nearly independent of the phase current (see Section 3.3).The resulting mean losses are summarized in Table 4.For the axial AMB constant, rotor For the differential winding design, which is used for the radial AMBs, the control current has to be smaller than I B .To allow control currents bigger than I B the AMBs would have to be rebuilt for a differential control design and a nonlinear control might be needed (see [15,[26][27][28]).Both would require major changes in the hardware.Hence, with the existing hardware, reducing I B also reduces the maximum force the actuators can generate.In the centered position, one AMB can generate a force of 278 N with I B = 5.67 A, but only 144 N with I B = 4 A. For the rotor without a specimen this is unproblematic, because there are nearly no external loads in the test rig and the rotor can be precisely balanced to minimize unbalance forces.However, the initial unbalance of an assembled specimen is big and balancing it precisely is difficult.The force generated by the AMB has to be big enough to counteract the unbalance force while passing the first three EFs.This has to be further investigated.

Increasing the Cycle Time
When no further loss reduction can be achieved, the cycle time has to be changed.In the upper example, even with a bias current of 4 A, the rotor temperature were above 110 • C but below 120 • C, as shown in Section 4.1.If further bias reduction or other loss reduction strategies are not applicable, a less favorable possibility to reduce the rotor temperature is to increase the cycle time.Here, two variations were considered, which are illustrated in Figure 14a.In the first variation, the acceleration of the rotor is reduced to lower load dependent losses in the PMSM, in the second variation, after a cycle, the rotational speed is held constant at 15,000 rpm for a while to allow the rotor to cool down.How long the cycle has to be increased depends on the associated loss reduction and the maximum allowable rotor temperature.For example, the cycle time could be increased from 30 s to 37 s to reduce the rms motor phase current from 37 A to 30 A, and thereby the rotor losses in the PMSM by 1.7 W to 17.4 W. Losses in the AMBs stay the same.Alternatively, after a 30 s cycle, the speed is held constant at 15,000 rpm for 7 s.Here, the mean losses can be calculated by averaging the mean losses for the first 30 s, as shown in Sections 3.5 and 3.6, and the losses at a constant speed.To keep the rotor without a specimen at a constant speed of 15,000 rpm, an rms phase current of 1.8 A is needed in the PMSM.Hence, it was assumed that an rms phase current of less than 10 A is needed to keep the rotor with a specimen at constant speed.Below 10 A, losses were nearly independent of the phase current (see Section 3.3).The resulting mean losses are summarized in Table 4.For the axial AMB constant, rotor losses of 8.1 W were assumed for both variations.The total cycle time in both variations was 37 s and the whole test would take the same amount of time-86 days to perform 200,000 cycles.losses of 8.1 W were assumed for both variations.The total cycle time in both variations was 37 s and the whole test would take the same amount of time-86 days to perform 200,000 cycles.Since losses in the PMSM showed only a small dependency on the load, the total mean losses in variation 2 were smaller than in variation 1.The calculated rotor temperature is shown in Figure 14b.For variation 1, the maximum calculated temperature in the permanent magnets of the PMSM was 115.6 °C, and for variation 2 it was 111.0 °C.Thus, with a slightly increased cooling time at 15.000 rpm, variation 2 should not reach the critical magnet temperature, whereas variation 1 has no big benefit for the rotor temperature.However, in the model the thermal conduction between rotor and magnets of the PMSM is assumed as ideal.In reality, small gaps could decrease the conductivity and therefore increase the dependency of the magnet temperature on the internal losses of the PMSM.This would make the loss reduction in the PMSM critical compared to the other components, giving variation 1 higher relevance.

Discussion
This paper proposed a structured approach to evaluate the heating of a magnetically levitated rotor in a vacuum, and deduced measures to avoid critical temperatures.First losses and the subsequent rotor heating were calculated.The rotational losses in the radial AMBs and the PMSM were identified as the main rotor loss mechanisms in the test rig.Losses in the axial AMB during rotation are smaller but also in the same order of magnitude.Switching losses in the AMBs and air friction losses only have a very small contribution to the total losses.One way to achieve loss reduction in the radial AMBs, without major changes in the spindle hardware, is to reduce the bias current.While doing so, the controllability of the rotor has to be ensured.The applicability of this approach was shown for the rotor without a specimen, but was not yet shown with a specimen.Stability problems could arise at low rotational speeds where many EFs can be excited.An adjustment of the bias with the rotational speed might be necessary.Further loss reduction can be accomplished by reducing the control activity.A centralized control of the radial AMBs should reduce the activities in the upper radial AMB.To realize a centralized control, however, the amplifier needs to be replaced.An overview of loss reduction strategies in AMBs can be found, for example, in [15].Losses in the PMSM can be reduced by increasing the switching frequency of its inverter or by employing an output filter.Both options also require changes in the hardware.If no further loss reduction can be achieved, the cycle time has to be adjusted, so that the critical rotor temperature is exactly reached.In the present case, the cycle time had to be increased.With a cycle time of 37 s and a bias current of  Since losses in the PMSM showed only a small dependency on the load, the total mean losses in variation 2 were smaller than in variation 1.The calculated rotor temperature is shown in Figure 14b.For variation 1, the maximum calculated temperature in the permanent magnets of the PMSM was 115.6 • C, and for variation 2 it was 111.0 • C. Thus, with a slightly increased cooling time at 15.000 rpm, variation 2 should not reach the critical magnet temperature, whereas variation 1 has no big benefit for the rotor temperature.However, in the model the thermal conduction between rotor and magnets of the PMSM is assumed as ideal.In reality, small gaps could decrease the conductivity and therefore increase the dependency of the magnet temperature on the internal losses of the PMSM.This would make the loss reduction in the PMSM critical compared to the other components, giving variation 1 higher relevance.

Discussion
This paper proposed a structured approach to evaluate the heating of a magnetically levitated rotor in a vacuum, and deduced measures to avoid critical temperatures.First losses and the subsequent rotor heating were calculated.The rotational losses in the radial AMBs and the PMSM were identified as the main rotor loss mechanisms in the test rig.Losses in the axial AMB during rotation are smaller but also in the same order of magnitude.Switching losses in the AMBs and air friction losses only have a very small contribution to the total losses.One way to achieve loss reduction in the radial AMBs, without major changes in the spindle hardware, is to reduce the bias current.While doing so, the controllability of the rotor has to be ensured.The applicability of this approach was shown for the rotor without a specimen, but was not yet shown with a specimen.Stability problems could arise at low rotational speeds where many EFs can be excited.An adjustment of the bias with the rotational speed might be necessary.Further loss reduction can be accomplished by reducing the control activity.A centralized control of the radial AMBs should reduce the activities in the upper radial AMB.To realize a centralized control, however, the amplifier needs to be replaced.An overview of loss reduction strategies in AMBs can be found, for example, in [15].Losses in the PMSM can be reduced by increasing the switching frequency of its inverter or by employing an output filter.Both options also require changes in the hardware.If no further loss reduction can be achieved, the cycle time has to be adjusted, so that the critical rotor temperature is exactly reached.In the present case, the cycle time had to be increased.With a cycle time of 37 s and a bias current of 4 A in both radial AMBs, a calculated steady state temperature of 111 • C can be reached, however, this has the drawback of an increased test duration.To account for modeling errors, the cycle time can further be adjusted during operation, depending on the rotor temperature, in a closed control loop.For the calculations, a constant thermal conductivity in the rotor was assumed.This assumption is especially critical between the magnets and the rotor.In reality, the thermal coupling will be poorer, resulting in a higher magnet temperature.

Figure 1 .
Figure 1.(a) Model the outer-rotor Flywheel described in [6] with a halved rotor; (b) Top view on the rotor showing the segments pressing against the fiber reinforced plastic (FRP) under rotation.

Figure 1 .
Figure 1.(a) Model the outer-rotor Flywheel described in [6] with a halved rotor; (b) Top view on the rotor showing the segments pressing against the fiber reinforced plastic (FRP) under rotation.

Figure 2 .
Figure 2. (a) Cross section of the test rig; (b) Cross section of the spindle.Figure 2. (a) Cross section of the test rig; (b) Cross section of the spindle.

Figure 2 .
Figure 2. (a) Cross section of the test rig; (b) Cross section of the spindle.Figure 2. (a) Cross section of the test rig; (b) Cross section of the spindle.

Figure 4 .
Figure 4. Applied approach to minimize the test duration by calculating and reducing rotor losses for the existing test rig.

Figure 4 .
Figure 4. Applied approach to minimize the test duration by calculating and reducing rotor losses for the existing test rig.

Figure 5 .
Figure 5. (a) Quarter of the cross section of the upper radial active magnetic bearings (AMB); (b) Magnetic flux density in the upper radial AMB with only the bias current as excitation.The rotor turns clockwise.For the loss calculation, the Bertotti formula[19] in the form of Equation (1) was applied, in which the three loss mechanisms  ,  , and  depend on the amplitude of the flux density  and the frequency f with which the magnetic flux changes:     , with     ,     , and    ..(1)

Figure 5 .
Figure 5. (a) Quarter of the cross section of the upper radial active magnetic bearings (AMB); (b) Magnetic flux density in the upper radial AMB with only the bias current as excitation.The rotor turns clockwise.

Figure 6 .
Figure 6.(a) Magnetic flux density in the axial active magnetic bearing (AMB); (b) Current density in the surface area of the thrust disk during switching of the axial AMB.

Figure 6 .
Figure 6.(a) Magnetic flux density in the axial active magnetic bearing (AMB); (b) Current density in the surface area of the thrust disk during switching of the axial AMB.

Figure 7 .
Figure 7. Rotor losses of the permanent magnet synchronous machine.

Figure 7 .
Figure 7. Rotor losses of the permanent magnet synchronous machine.

Figure 8 .
Figure 8.(a) Rotational speed and calculated rotational losses of the active components during one cycle; (b) Resulting mean losses of all calculated loss mechanisms in the rotor during cycling.

Figure 8 .
Figure 8.(a) Rotational speed and calculated rotational losses of the active components during one cycle; (b) Resulting mean losses of all calculated loss mechanisms in the rotor during cycling.

Figure 9 .
Figure 9. (a) Simulated and measured rotor cooling without the specimen; (b) Simulated and measured rotor temperatures without the specimen during cycling.

Figure 9 .
Figure 9. (a) Simulated and measured rotor cooling without the specimen; (b) Simulated and measured rotor temperatures without the specimen during cycling.

Figure 10 .
Figure 10.(a) Simulated and measured rotor temperatures with adjusted AMB losses without the specimen during cycling; (b) Simulated steady state temperature of the rotor with and without the specimen.

Figure 10 .
Figure 10.(a) Simulated and measured rotor temperatures with adjusted AMB losses without the specimen during cycling; (b) Simulated steady state temperature of the rotor with and without the specimen.

Figure 11 .
Figure 11.(a) Simulated and measured rotor temperatures with a reduced bias current of 4 A in both AMBs without the specimen during cycling; (b) Steady state rotor temperature with a reduced bias current of 4 A in both AMBs with and without the specimen.
shows the Campbell diagram of the position data of the upper AMB with  5.67 A. The lines starting at the lower left corner are harmonics of the rotational speed, the thickest one being the first harmonic.Odd harmonics of the rotational speed are clearly visible, whereas even harmonics are barely detectable.The approximately horizontal lines are eigenfrequencies (EFs) of the system.At about 90 Hz lies the rigid body tilting mode of the rotor in the AMBs.The translational rigid body mode is not visible in the diagram.The first bending mode of the rotor can be seen at 1057 Hz, the second at 2063 Hz.The second splits into

Figure 11 .
Figure 11.(a) Simulated and measured rotor temperatures with a reduced bias current of 4 A in both AMBs without the specimen during cycling; (b) Steady state rotor temperature with a reduced bias current of 4 A in both AMBs with and without the specimen.

Figure 12 .
Figure 12.(a) Measured Campbell diagram and deflection of the rotor with a bias current of 5.67 A without specimen; (b) Measured Campbell diagram and deflection of the rotor with a reduced bias current of 4 A without specimen.

Figure 12 .
Figure 12.(a) Measured Campbell diagram and deflection of the rotor with a bias current of 5.67 A without specimen; (b) Measured Campbell diagram and deflection of the rotor with a reduced bias current of 4 A without specimen.

Figure 13 .
Figure 13.(a) Calculated bending eigenmodes below 1 kHz of the rotor with sample; (b) Calculated Campbell diagram of the rotor with the specimen.Each eigenmode shown in (a) splits in a forward and a backward mode.

Figure 13 .
Figure 13.(a) Calculated bending eigenmodes below 1 kHz of the rotor with sample; (b) Calculated Campbell diagram of the rotor with the specimen.Each eigenmode shown in (a) splits in a forward and a backward mode.

Figure 14 .
Figure 14.(a) Variations of increasing the cycle time to reduce the rotor temperature; (b) Steady state rotor temperature of both cycle variations with a reduced bias current of 4 A in both AMBs with a specimen.

Figure 14 .
Figure 14.(a) Variations of increasing the cycle time to reduce the rotor temperature; (b) Steady state rotor temperature of both cycle variations with a reduced bias current of 4 A in both AMBs with a specimen.
Actuators 2019, 8, x FOR PEER REVIEW 3 of 18 state of stress in the test specimen, a circumferentially segmented steel ring is placed inside the FRP cylinder of the specimen, and rotated at 30,000 rpm.This results in a transverse compressive stress of 100 MPa and a longitudinal tension of 775 MPa in the FRP.With an outer diameter of 190 mm, the surface speed of the specimen reaches approximately 298 m/s.The fatigue test will be conducted by alternating between a maximum speed of 30,000 rpm and a minimum speed of 15,000 rpm.At minimum speed the transversal compressive stress is 26 MPa and the longitudinal tension is 200 MPa.

Table 1 .
Rotational losses in the radial AMBs.

Table 1 .
Rotational losses in the radial AMBs.

Table 2 .
Air friction losses in the test rig with  0.05 Pa and  290 K.

Table 2 .
Air friction losses in the test rig with p = 0.05 Pa and T = 290 K.

Table 3 .
Adjusted mean losses in the rotor during cycling.

Table 3 .
Adjusted mean losses in the rotor during cycling.

Table 4 .
Mean rotor losses during cycling variation 1 and 2 with a reduced bias current of 4 A.

Table 4 .
Mean rotor losses during cycling variation 1 and 2 with a reduced bias current of 4 A.

Table A3 .
Measured stator temperature during cycling with a bias current of 4 A.