Research on Ceramic Spindle Loss Based on Ceramic Reverse Magnetic Effect

Ceramic motorized spindle is a multifield, nonlinear, and strong coupling system.3e present model ignores the reverse magnetic effect and has a poor accuracy on the loss of the spindle system. In the paper, Park transformation was used to establish the electromagnetic physical model of ceramic motorized spindle. By combining with the Jiles–Atherton theory, an eccentric magnetization model of a ceramic motorized spindle considering reverse magnetic properties was established. Energy consumption parameters were calculated under various working conditions.3e influence of ceramic reverse magnetic characteristics on the motorized spindle loss was analyzed and verified by experiments. 3e results show that the simulation of the ceramic motorized spindle loss model is in good agreement with the experimental results. Compared with the experimental results, the average loss error was 2.1%. Due to the reverse magnetic characteristics, the application of ceramic motorized spindle can help reduce the system loss. 3e ceramic motorized spindle model provides a theoretical basis for the development of ceramic spindle energy loss and efficiency.


Introduction
As the core component of computerized numerical control (CNC) machine tools, the high-speed motorized spindle has a complex coupling system. Due to the reverse magnetic properties, the ceramic motorized spindle can work better under some special conditions, such as the nonmagnetic field of petrochemical and strong chemical corrosion environment [1][2][3]. Energy development and management are crucial to further improve the machining efficiency of motorized spindles in high-speed operation [4,5]. Scholars have undertaken a significant amount of research. Many smart energy management algorithms are presented through Newton-Raphson and event-triggered-based distributed [6,7]. In general, some loss minimization algorithms are proposed, which can predict the load torque and loss of induction motor. However, the effect of torque variations on efficiency was not deeply analyzed in higher order harmonic environments [8,9]. Ismagilov et al. [10]calculated the loss of the high-speed spindle stator yoke. Wang et al. [11] established a novel small-signal model based on the characteristic equation, which can assess the system lowfrequency stability, but the field of rotor magnetic density was not studied. e common methods in motor are generally researched through air-gap flux density distribution, but in the spindle working process, the eccentric effect of the air gap was not considered [12,13]. Some researchers [14][15][16] proposed a dynamic adjustable induction motor (AIM) considering the rotor skin effect on current, which improved the energy efficiency of the manufacturing process. However, in some cases, energy efficiency did not increase significantly in a large-load situation. Other researchers reported that harmonic had great impact on the magnitude of current. However, the reverse magnetic properties of the rotating shaft materials were not taken into account [17][18][19]. Many researchers [20][21][22][23] proposed some new energy-saving management, which can enhance efficiency through optimizing the motor current and torque, but the research studies did not involve a better fit for the ceramic spindle in some special cases. Nowadays, little investigation has been conducted to improve the loss, harmonics, efficiency, energy management, and other aspects of the spindle. erefore, it is essential to pay enough attention to the magnetic field and energy loss of the ceramic motorized spindle.
In this paper, a general electromagnetic coupling model is established, which is suitable for a ceramic motorized spindle. e model takes into account the characteristics of the ceramic reverse magnetic effect. According to Park transformation, considering the deep analysis of magnetic field, the motor model of the ceramic motorized spindle is established. Based on the Jiles-Atherton theory, the mutual influence between the air gap, magnetic density, and loss of the ceramic motorized spindle is analyzed under different working conditions. e loss model of the ceramic motorized spindle is studied, revealing the relationship between reverse magnetic, harmonics, and loss. e electromagnetic model of the ceramic motorized spindle has high accuracy, which has been verified by experiments. A ceramic motorized spindle is of great significance to improve the energy output performance and reduce loss in extreme environments.

Electromagnetic Model of a Ceramic
Motorized Spindle

Mechanism Model of a Ceramic Motorized Spindle.
A ceramic motorized spindle directly connects induction motor with spindle. Driven by electromagnetic torque through an inverter, the ceramic spindle can rotate at high speed. It is a complex mixed system. ere are complex electromagnetic coupling relationships in the ceramic motorized spindle. e shaft of the motorized spindle is made of a ceramic material. When the motorized spindle rotor rotates, the stator and rotor are magnetically connected. According to Faraday's law of electromagnetic induction, the mathematical model of magnetic flux is established. e internal space vector of the motorized spindle is shown in Figure 1(a). e axis of rotation and space complex plane are represented by complex coordinate. e axis of winding is served as the real space reference axis. G is the movement trajectory of space trajectory r. e motorized spindle vector parameters of space trajectory, real axis, and imaginary axis are shown as follows: where θ r is the electrical angle of the space vector and R is the amplitude of the space vector. e relationship between current and the space vector in stationary and rotating axes is shown in Figure 1 e rotating magnetic field is generated by three-phase windings with high-speed rotation. e stator and rotor are projected into d and q axes through Park transformation. Park transformation can reduce the variables between electromagnetic parameters. e three-phase motor model of the motorized spindle is converted into a two-phase α and β static coordinate system and a twophase d and q rotating coordinate system. e two-phase α and β stationary coordinate system rotates at ω angular velocity. e angle θ of rotor and stator flux direction is relative to the slip ratio s. A backward angle allows the flux to produce a constant current. According to Figure 1(b), the transformation relationship between the current's changes in coordinate axes can be obtained by applying Park transformation, as shown in the following equation [24]: e relationship between flux and circuit in dq0 coordinates is shown in the following equation: where L is the self-induction and M is the mutual induction. e magnetic field density of the air gap is not only related to the material permeability but also has a direct relationship with the length of the air gap [25]. Due to the eccentricity of the stator and rotor, the length of the air gap is subject to a nonlinear variation rule. e eccentricity of the air gap is mainly caused by the rotor vibration in high-speed rotation. e air gap eccentricity and magnetic potential are shown in Figure 2.

Eccentric Magnetization Model of a Ceramic
In the eccentric state, the relationship between the air gap and eccentric position is shown in the following equation [26]: In the eccentric state, the element shape in the air gap is distorted. e Jacobi derivative method is used to calculate the unbalanced magnetic density in the eccentric state. e harmonic loss of the ceramic motorized spindle cannot be ignored. e eccentric air gap magnetic density of the stator and rotor is shown in the following equation [27]: where ε is the relative eccentricity, ε � d sr /R s −r r , v is the stator harmonic frequency, u is the rotor harmonic frequency, p is the pole logarithm of the motor, and φ 0 is the power factor angle.
Air gap magnetic common energy is a function of stator, rotor angle, and magnetic energy. It is directly related to the magnetic density of the motorized spindle, as shown in the following equation: where D is the inner diameter of the stator, l is the effective length of the stator and rotor, θ is the angle between the stator and rotor magnetic potential, μ 0 is the permeability of the air gap, and μ c is the permeability of ceramic.

Magnetization Model of a Motorized Ceramic Spindle.
Generally, the internal magnetic circuit of the motorized spindle is divided into five sections: stator yoke, stator tooth, air gap, rotor tooth, and rotor yoke. e rotating shaft is made of the ZrO 2 ceramic material. Rotor magnetization is mainly affected by rotational magnetization. e rotor yoke is affected by the ceramic reverse magnetic material, and the magnetic inductance line is distorted. ere are magnetization alternating and magnetization rotating modes in the motorized spindle. e rotating magnetic field breaks down into two alternating magnetic fields, differing by 90 degrees, as shown in the following equation: e magnetization intensity of the ceramic motorized spindle rotor is orderly arranged by the magnetic moment of silicon molecules. e magnetic moment vector sum of all silicon molecules is shown in the following equation:

Mathematical Problems in Engineering
where △m is the additional magnetic moment of the molecule, m is the inherent magnetic moment of the molecule, △l is the stator and rotor contact length, S is the stator and rotor contact area, △V is the volume of magnetized medium, j m is the current density of magnetized surface, and n is the number of molecules magnetized in a medium volume. e magnetic field intensity of the ceramic rotating shaft is shown in the following equation: where α is the curve parameters of ceramic magnetization (α < 1) and H r is the magnetic field intensity of the inner rotor domain. A magnetization formula is obtained with the Langevin function, as shown in the following equation: According to the ampere loop theorem in ceramic and metal magnetic media, the magnetization relationship between B, M, and H is obtained as shown in the following equation: where B m is the magnetic medium density of the spindle.

Research on Electromagnetic Characteristics of a Ceramic Motorized
Spindle. e influence of the complex motorized spindle electromagnetic characteristics is studied under different working conditions. e electromagnetic parameters of the motorized spindle eccentric magnetization model are calculated based on the ceramic reverse magnetic characteristics. e simulation condition is that the shaft material of 170SD30 ceramic motorized spindle is ZrO 2 . e metal motorized spindle shaft material of the 170HT30 is 40Cr. e rotational speed was 9000 rpm, and the air gap was 0.3 mm. Under loading conditions, the speed difference was increased from 250 rpm to 1000 rpm, increasing by 250 rpm each time. Simulation results show the current and magnetic field density of the metal and ceramic motorized spindle changes under different working conditions. e changes in current under different working conditions are shown in Figure 3. e magnetic field density distribution of the spindle at different conditions is shown in Figure 4. According to Figure 3(a), the metal and ceramic motorized spindles can be seen under the no-load condition, in which the current of the ceramic motorized spindle reaches 7.71 A and the metal motorized spindle reaches 9.36 A. e current of the metal motorized spindle is slightly higher than the ceramic motorized spindle. e main reason is the ZrO 2 shaft which has reverse magnetic properties. e magnetization moment generated by the rotor cannot make the large number of ceramic material molecules arranged in an orderly manner. e reverse magnetic properties of ceramics can effectively reduce the generation of induced current. At the same time, the alternating frequency of the magnetic field is higher. e distribution of magnetic flux density is not uniform due to the influence of eddy current. e magnetic flux is relatively large on the surface of the rotor. e skin effect on the metal motorized spindle current is obvious, which is the main reason that the current fluctuation of the metal motorized spindle is larger than the ceramic motorized spindle. Figure 3(b) shows that the motorized spindle current increases with the increase in air gap eccentricity. Because of the eccentricity of the air gap, the rotating flux of the stator and rotor is not uniform, which leads to the increase in current. Meanwhile, the current fluctuation with an air gap eccentricity of 0.2 mm is greater than that without eccentricity. e extreme value of 11.49 A in the eccentricity moment is more than 40% higher than the current without eccentricity. Due to the uneven distribution of the air gap, the rotation flux fluctuates violently. Air gap eccentricity changes the frequency of magnetic flux and produces different depths of skin eddy current. It has a direct influence on current fluctuation. Figure 3(c) shows that the motorized spindle current presents a rapid upward trend with the increase in the speed difference. When the loading slip ratio reaches 0.1111, the current reaches 55.875 A. When the spindle is running under a larger load, the stator and rotor current of the spindle will increase significantly. is will cause a sharp increase in leakage of the spindle magnetic field. e iron core of stator and rotor tooth positions is seriously saturated. e iron loss at the rotor tooth position of the spindle will increase under load, which makes the motor rotor current not only increase in amplitude but also increase in its fluctuation frequency. e distribution diagram of rotor magnetic pole density is shown in Figure 4(a). It can be seen that the magnetic pole density of the ceramic motorized spindle rotor is smaller than the metal motorized spindle. e maximum magnetic pole of the metal motorized spindle rotor can reach 1.097 T. e maximum magnetic pole of the ceramic motorized spindle rotor is 50.1% lower than that of the metal motorized spindle. e ceramic shaft has reverse magnetic properties. Ceramic ZrO 2 molecules cannot be arranged in order. e magnetization phenomenon of the magnetic sensing line formed by the stator and rotor is weak. e material of the metal motorized spindle rotating shaft is 40Cr. e metal molecules are arranged in an orderly manner under the action of magnetic torque, and the magnetization phenomenon is obvious. Figure 4(b) shows the magnetic density distribution diagram about the middle position of the motorized spindle shell. Compared with Figures 4(b) and 4(a), the trend is consistent, though the magnetic density of the metal motorized spindle is higher than the ceramic motorized spindle. e magnetic density of the motorized spindle shell is much smaller than the rotor's magnetic pole.
e maximum magnetic density of the metal motorized spindle shell is 15.598 μT. e maximum magnetic density of the ceramic motorized spindle shell is 10.094 μT. e main reason is that the magnetic induction line is affected by the transmission distance in the transmission process. At the same time, the interference of the water channel, bearing housing, and other related parts is also the reason that the shell magnetic density is small. e distribution diagram of ceramic motorized spindle eccentric magnetic density is shown in Figure 4(c). Figure 4(c) shows that the maximum magnetic density of the air gap eccentricity of 0.1 mm is 0.749 T. Under the condition of 0.2 mm air gap eccentricity, the maximum magnetic density is 1.173 T, and the magnetic density alternating velocity is higher than the normal air gap magnetic density fluctuation. Due to the uneven distribution of air gap, the magnetic flux fluctuation is obvious. Especially in places with a small air gap, the alternating frequency of the magnetic field increases and the flux fluctuates dramatically, appearing like the phenomenon of magnetic density increases instantaneously.
Additionally, due to the ceramic reverse magnetic and the harmonic effect, the current and magnetic density of spindle is simulated through the magnetization model. In Figures 3 and 4, it is obvious that the magnetic and electric fields of ceramic spindle are smaller than those from metal spindle, which shows better adaptation for certain circumstances.

Loss Model of a Ceramic Motorized Spindle
e loss of a motorized spindle has a direct relationship with the thermal value of the spindle, which directly affects its thermal characteristics. e loss of a ceramic motorized spindle mainly includes mechanical loss and electromagnetic loss. Electromagnetic loss is directly related to alternating magnetization and rotational magnetization [28]. e establishment of a modified empirical coefficient model can improve the accuracy of predicting energy loss and efficiency under different working conditions.

Iron Consumption Model of a Ceramic Motorized Spindle.
Iron loss is divided into hysteresis loss, eddy current loss, and abnormal loss. Hysteresis loss is accompanied by external magnetic field change. e rotor is made by silicon steel sheets. In a high-frequency alternating magnetic field, the skin effect is the main factor affecting the loss.

Hysteresis Loss.
e fundamental reason for hysteresis loss is that the change in magnetic induction intensity lags behind the magnetic field intensity. e hysteresis of ceramic material is much smaller than the silicon steel material. e magnetic domain has a direct influence on the hysteresis loss mechanism. e magnetic domain distribution of ceramic and metal materials is shown in Figure 5. Figure 5 shows the arrangement of the motorized spindle molecules at a random moment under the working condition of H magnetization intensity. When the external magnetic field M is applied to the ceramic material, only a few magnetic domain direction inside the ceramic materials can be arranged in order, and other molecules are arranged in disorder.
e results show that only a few ceramic molecules in the magnetic domain change with the direction of magnetization. ere is a weak relationship between the internal ceramic molecular structure and the external magnetic field. Under the action of external magnetic field M, there are many tiny spontaneous magnetization regions inside ferromagnetic materials. e molecules of silicon steel can be arranged in an orderly way. With the change in magnetization H inside the motorized spindle, most silicon steel molecules can achieve the same direction as the magnetization intensity M. A small number of silicon steel molecules cannot be quickly aligned with the M direction, which is mainly due to the intrinsic properties of ferromagnetic materials, magnetic crystal anisotropy interactions, and high-frequency magnetic field variation. It indicates that ferromagnetic materials have a direct influence on the external magnetic field, resulting in hysteresis of ferromagnetic materials. Hysteresis loss of the ceramic spindle is shown in the following equation: where K c is the hysteresis coefficient of ceramic (K c < 0.5), l m is the effective length of the stator and rotor, and A j is the area of the j unit.

Loss of Eddy Current and Abnormal.
A silicon steel rotor can induce an electromotive force in a high-frequency magnetic field. e electromotive force causes an eddy current in the closed rotor material. e eddy current is a circular current line induced in the silicon steel conductor and is one of the ceramic spindle main loss. e eddy current loss is directly related to the changing mode of the magnetic field and the permeability of the conductor. Reverse magnetic materials have less eddy current loss when rotating in an uneven magnetic field than paramagnetic materials. e molecules of ceramic materials cannot be arranged in large numbers, and induced eddy current generation can be effectively reduced so that large eddy current loss is avoided. e influence of magnetic domain movement on current is shown in Figure 6. According to the microscopic level of Figure 7, the eddy current loss generated by moving magnetization is analyzed. e amplitude of magnetic density produced by domain wall movement in silicon steel is twice that of the ceramic shaft.
Factors such as saturation, slot effect, and winding arrangement produce a large number of harmonics. A harmonic magnetic field leads to eddy current loss on the surface of the stator and rotor. A large number of voltage harmonics and higher harmonic voltage output by the frequency converter will have a great impact on the iron consumption of the spindle. e magnetic density in the iron core of the motorized spindle is distorted. e reverse magnetic effect of ceramic makes the magnetic density distribution of rotor tooth and yoke change, which also makes the harmonic magnetic field and the phenomenon of rotating magnetization change. During the rotation process of the ceramic motorized spindle, the magnetic intensity vector of the trajectory coordinate changes with the rotation θ angle. Eddy current loss of the magnetic density presents an elliptic shape. e v harmonic elliptic magnetic density diagram is shown in Figure 7. e magnetization loss proportional coefficient of the ceramic motorized spindle is K a . e density ellipticity of the ceramic motorized spindle δ is linearly related to the highfrequency harmonics. e hysteresis phenomenon of silicon steel sheet material is obvious and seriously affects the loss. e abnormal losses are shown in the following equation: Influenced by the reverse magnetic effect of ceramic, the high-order loss coefficient of the harmonic eddy current K 2v is used to analyze the harmonic magnetic density of the rotor and the influence of rotor electrical density on spindle energy loss under different frequencies. e harmonic effect was considered. Bertotti's classical iron consumption model considered the eddy current loss and is shown in the following equation: where B n is the magnitude of the n harmonic of the magnetic density through Fourier transform decomposition; K e is the eddy current loss coefficient of silicon, K e � π 2 cd 2 /6ρ; d is the length of silicon; ρ is the density of silicon; c is conductivity; and β e is the v harmonic additional eddy current loss coefficient of magnetic density.

Mathematical Problems in Engineering
According to equations (12)- (14), the iron consumption of the ceramic motorized spindle can be obtained, as shown in the following equation:

Copper and Mechanical Loss of a Ceramic Motorized
Spindle. Copper loss is an important part of the motorized spindle electromagnetic loss. e copper loss of a motorized spindle not only reduces the motor efficiency but also increases the temperature of the stator winding, which affects the motor's operation performance. According to Joule's law, the copper loss is the product of the stator winding current squared times the resistance. Stator current is mainly composed of excitation current and load current. e excitation current of I cu1 increases with the increase in rotating speed, while the load current of I cu2 decreases with the increase in rotating speed. According to the energy conduction law of a motorized spindle, there is a linear relationship between the ratio of excitation current and load current with the change in load. erefore, the copper loss of the stator and rotor is shown in the following equation: where K v is the resistance increasing coefficient of v harmonics (K v < 0.5), s is the slip ratio, P em is the electromagnetic power, n 0 is the preload speed, and n 1 is the load speed. Mechanical loss is mainly composed of bearing friction loss and loss of ventilation. Bearing friction loss is directly related to bearing speed, the friction coefficient, the number of rolling bodies, and bearing load. Mechanical loss is shown in the following equation: where P f is the mechanical loss, F is the bearing load, d is the diameter of the ball, v is the rotation speed of the ball, and q is the air volume flow rate of the motorized spindle.

Loss Factors of a Ceramic Motorized
Spindle. e influence of complicated ceramic motorized spindle loss was studied under different working conditions. e simulation condition is that the speed of 170SD30 ceramic motorized spindle is 9000 rpm. Calculation of the spindle loss model was based on ceramic reverse magnetic characteristics. Energy loss parameters, such as iron loss, copper loss, and output power, were calculated under no-load and different slip ratio loading conditions. e iron loss diagram of the ceramic motorized spindle under different working conditions is shown in Figure 8. e copper loss diagram is shown in Figure 9.  Figure 6: e influence of magnetic domain movement on current: eddy current of (a) metal spindle and (b) ceramic spindle. Figure 7: e v harmonic elliptic magnetic density diagram.

Mathematical Problems in Engineering
It can be seen from Figure 8(a) that the ceramic motorized spindle loss is about 100 W under no eccentric load condition, which is less than the iron loss of the ceramic motorized spindle under eccentric working condition and only 2/3 of the iron loss of the 0.2 mm eccentric. e lower iron loss in the noneccentric working condition is due to the lower magnetization degree of the reverse magnetic material. e ceramic motorized spindle rotor's silicon steel is affected by the ceramic reverse magnetic property, which reduces the hysteresis iron loss and hysteresis degree of the rotor silicon steel material. Due to the eccentric air gap, the distribution of the air gap is uneven, and the smallest gap is only 0.1 mm. e magnetic flux alternating frequency increases. At the same time, the flux change has an obvious influence on the current skin effect. e iron loss also increases with the increase in the air gap unevenness. Figure 8(b) shows that the iron loss of the ceramic motorized spindle increases nonlinearly with the increase in the slip ratio. When the slip ratio reaches s � 0.1111, the iron consumption exceeds 700 W and is almost twice the slip ratio under the condition of s � 0.0278. is is mainly because the rotor current increases rapidly in the process of a large loading, in which the current exceeds 50 A in a short time. e magnetic domain wall experiences a jump rotation. e two domains appear as vibration through rubbing against each other, which leads to loss. e movement of the magnetic domain wall makes the Copper loss is mainly affected by Ohm's law. It can be seen from Figure 9 that the change rule of copper loss is similar to iron loss. With the increase in eccentricity and slip ratio, the copper loss and load will increase. is makes the spindle rotor current not only increase in amplitude but also increase in alternating current frequency due to the uneven width caused by eccentricity. e copper loss of the stator winding increases rapidly. By analyzing Figures 9(a) and  9(b), it can be seen that the maximum slip ratio s � 0.1111 reaches 240 W under loading conditions, which is much higher than under no-load and low-load conditions. e influence of the motorized spindle load on copper loss is greater than that of eccentricity. It also shows that the load has a more obvious influence on current.
According to equation (13), the friction and mechanical loss distribution diagram is obtained, as shown in Figure 10.
ere is a nonlinear relationship between the mechanical loss and the speed of the ceramic motorized spindle. In the process of the spindle rotating at medium and low speed, the mechanical loss has a relatively low impact on the loss of the spindle system. e loss of the bearing to the spindle during high-speed rotation cannot be ignored, and the mechanical loss exceeds 50 W under the condition that rotation speed is more than 2000 rpm. e friction between the ball and the outer ring channel generates heat. With the increase in bearing speed, the relation of generating heat square increases nonlinearly. It can also be seen from the figure that the mechanical loss of the metal motorized spindle is greater than the ceramic motorized spindle. is is mainly because of the self-lubricating properties, especially during highspeed rotation. Under the condition of constant oil supply, ceramic bearings have a more prominent lubrication effect and lower generation of heat. Ceramic bearings are characterized by high hardness and small deformation. e stiffness of the ball is not easy to soften in high-speed rotation. Compared to the contact area of the metal bearing channel, the ceramic is smaller. is is the reason that ceramic bearings generate less heat.

Electromagnetic and Loss Experiments of a Ceramic
Motorized Spindle. Oil and gas lubrication is adopted for spindle bearing lubrication under the test conditions. e experimental conditions are input oil of 5 mL/h, air pressure of 0.25 MPa, cooling water temperature of 289.15 K, and cooling water flow rate of 5 L/min. e speed is 9000 rpm, and the dragging loading test was carried out by another spindle. e loading range is from 9000 rpm to 8000 rpm, with speed increments of 200 rpm. A HIOKI FT3470-52 magnetic field tester and HIOKI MR8875-30 data acquisition instrument were used to conduct a magnetic density test of the spindle magnetic field for a 170SD30 ceramic spindle and 170HT30 metal spindle. A FLUCK NORMA 4000/5000 power analyzer was used to measure spindle current and harmonics. e measure of the magnetic field is shown in Figure 11, and the power analyzer measure of the ceramic motorized spindle is shown in Figure 12.
e basic parameters of spindle are shown in Table 1, and the structural parameters of spindle are shown in Table 2.

Experimental Results and Analysis.
Magnetic field density was measured for the ceramic motorized spindle and metal motorized spindle shell. Magnetic density data were collected by using the HIOKI FT3470-52 magnetic field tester and HIOKI MR8875-30 data acquisition instrument. e magnetic density of spindle shell was measured under different load conditions; the comparison figure is shown in Figure 13. As can be seen from Figure 13(a), the magnetic density of the ceramic motorized spindle shell is smaller than the metal motorized spindle. e maximum magnetic density of the ceramic motorized spindle is only 9.946 μT, while for the metal motorized spindle is 16.37 μT. Figure 13(a) is comparable with Figure 5(b). e magnetic density error of the metal motorized spindle is 2.9%, while the ceramic motorized spindle magnetic density error is 1.89%. e magnetization model can be precisely verified through the comparisons with different spindles. Figure 13(b) is the comparison diagram of magnetic density simulation and experiment under different loading conditions. It can be seen that different loads are loaded at the operating condition of 9000 rpm. Magnetic density and loading show a nonlinear increasing relationship. e maximum magnetic density error of the motorized spindle shell in simulation and testing is 2.14 μT.
e magnetic density distribution of the shell can also reflect the magnetic density distribution change in the rotor pole. e harmonic frequency domain, current, and power of the motorized spindle were obtained using the FLUCK NORMA 4000/5000 power analyzer. e harmonic frequency domain and current of the motorized spindle are shown in Figure 14. e relationship between spindle loss and torque load is shown in Figure 15. Figure 14(a) analyzes the harmonic spectrum of the current under the condition of a motorized spindle at 9000 rpm, for which the speed difference is 400 rpm. e energy of the spindle in the lowfrequency region is greater than that in the middle and highfrequency regions. It can be concluded that the reverse magnetic properties of ceramic can effectively reduce the influence of harmonics on current compared with the metal paramagnetic properties. Figure 14(b) shows that, in the simulation calculation considering harmonics, the average error of 1.91% is in good agreement with the current value measured in the experiment. e simulation current without considering harmonics is smaller than the experimental current, with an average error of 9.43%. e current error was as high as 16.7% under some large load conditions. It is shown that harmonics are an important factor, affecting the fluctuation and magnitude of current. e input and output power of the ceramic motorized spindle can be measured through a power analyzer. Figure 15(a) shows the relationship between input power, torque, and the speed difference of the motorized spindle. In the figure, it can be seen that the input power exhibits a    nonlinear increasing relationship with the increase in speed difference and torque. Under a low-load condition, torque increases rapidly, while under middle load conditions, torque and input power increase slowly. Under the condition that the rotation speed difference reaches 1000 rpm, the input power reaches the maximum value of 5505 W. Figure 15(b) shows the ceramic motorized spindle variation diagram of the total loss under different loading conditions. At the speed difference of 500 rpm, the maximum power difference is 32.47 W, the total simulation loss is 639.84 W, and the total loss of the experiment is 607.37 W. e existing simulation model for spindle loss is usually larger than reality. And the simulation accuracy can be obviously improved by considering the characteristics of ceramic reverse magnetism and the harmonic effect. e proposed spindle loss model simulating the maximum loss error is 5.07% and the average error is 2.1%. It is meaningful to develop a new electromagnetic spindle loss mechanism model to simulate the complex conditions.

Conclusions
An electromagnetic loss model of a ceramic motorized spindle considering the reverse magnetic effect is proposed.
is model can be used to analyze the influence of parameters, such as current, magnetic density, and loss under different working conditions. e influence of harmonics, reverse magnetic characteristics, and eccentricity on the ceramic motorized spindle magnetic field is analyzed. In addition, the influence of different loads on the ceramic motorized spindle copper and iron and the mechanical loss is studied. e conclusions are as follows:  (1) e electromagnetic physical model of a spindle system was studied. e magnetization eccentricity model of a ceramic motorized spindle was established by combining the Park formula and the Jiles-Atherton theory. e magnetic density of a ceramic motorized spindle at different conditions was simulated. Air gap eccentricity and magnetic density showed a nonlinear increasing relationship. e air gap magnetic density of a ceramic motorized spindle was obviously affected by the reverse magnetic effect. e air gap magnetic density of ceramic motorized spindle was only 64.7% of metal motorized spindle.
(2) Considering the magnetization model of the reverse magnetic effect, the ceramic motorized spindle's current and loss were smaller than the metal motorized spindle. Due to the hysteresis and eddy current, the loss mechanism of the ceramic motorized spindle was changed by the ceramic reverse magnetic property. (3) e loss model of the ceramic motorized spindle was calculated by Ohm's law. Combined with Bertotti's classical iron loss model, the influence of a ceramic bearing on the ceramic motorized spindle mechanical loss was considered. e distribution of loss under different load conditions was studied. e loss of the ceramic motorized spindle increased nonlinearly with the increase in load. Compared with the experimental results, the average loss error was 2.1%. e electromagnetic model of a ceramic motorized spindle has high accuracy.