Research on Vibration Characteristics of a Ceramic Spindle Based on the Reverse Magnetic Effect

/e electromagnetic radial force about a ceramic spindle affects the spindle dynamic, which determines the quality of processing. Using a Timoshenko beam unit to build the dynamic model for the ceramic spindle, the dynamic characteristic of an angular contact ball bearing is analyzed using a nonlinear bearing model. /e electromagnetic magnetization model was established based on Maxwell’s theory to calculate electromagnetic magnetic density and radial force. /e influence about the reverse magnetic field characteristic of the ceramic rotating shaft and dynamic stiffness of the contact ball bearing on the dynamic phenomena of the spindle is analyzed, which is verified by experiments. /e results show that the magnetic effect produced by the reverse magnetic of ceramic rotation shaft has a great influence on the electromagnetic radial force. Compared with the paramagnetic effect of the metal shaft, the dynamic characteristics of the spindle can be significantly improved. Considering the coupling relationship between the radial force of the magnetic field and the bearing contact force, dynamic stiffness, and other factors, the accuracy of the model simulation is highly consistent with the test results. In particular, the ceramic spindle model has been successful in predicting with high accuracy and is suitable for multiple extreme working conditions. /e parameters, such as initial eccentricity of the rotor, bearing preload, and rotating speed, can be adjusted to restrain the vibration of spindle. /e ceramic spindle model provides a theoretical basis for the dynamics development of a high-speed spindle.


Introduction
High-speed machine tools play an important role in today's manufacturing industries and can significantly improve product surface quality and machining accuracy [1][2][3].Due to the ceramic materials characteristics of high stiffness, small deformation, and reverse magnetic, in special extreme environments such as Aerospace area (ultrahigh temperature and ultralow temperature), high corrosion area, and high requirement for magnetic area (ultralow magnetic field), the performance of ceramic spindle is better than metal spindle.However, the vibration of the ceramic spindle is large in these extreme environments.It is of great significance to study the dynamic performance of the highspeed ceramic spindle.
Considerable research studies related to improve a more efficient machine tool, energy savings, and vibration avoidance of the spindle have been published.e purpose of modeling the spindle is to find the shaft deformation, bearing characteristics, and magnetic effect which have the greatest effect on the static and dynamic properties of the spindle systems.Lee and Hoang [4] studied the consideration of the continuous mass effects, eccentricity, and rotational inertia give a more realistic model that was established through the Timoshenko beam for spindle system and investigated the coupled lateral and torsional effects on the system vibration behavior.However, the model of An-Chen Lee neglects the radial force on the shaft.Zivkovic et al.'s model [5] is based on the mechanical model of the bearing and the numerical model FEM (Finite Element Method) of the spindle to predict bearing characteristics and the change in static stiffness of the spindle accuracy.Different speed responds are verified by experiments.But the dynamic stiffness is not mentioned in the model.Hu et al. [6] studied the effects of preload and preload method on the rotational performance of the spindle-bearing system under different speeds.However, the paper is not to solve the coupled problems of the spindle-bearing systems.Hong et al. [7] described a kinematic model of a five-axis machine tool under position-independent and positiondependent geometric errors associated with rotary axes through the experiments to verify the enlargement of a periodic radial error motion of C-axis with B-axis rotation.e model can not perform the vibration accurately in some extreme cases.Some papers [8][9][10] proposed speed-dependent FRF (Frequency Response Function) spindle model.
e coupled spindle model based on Timoshenko beam can be predicted.But the predicted models are not better adapted to the high stiffness of the ceramic spindle.
Vyroubal [11] presented a method focused on compensation of machine's thermal deformation in spindle axis direction based on decomposition analysis that introduces cheap and effective strategy how to compensate machine tool thermal deformation to reduce the residual thermal deformation greatly.Yang et al. [12] proposed least square support vector machine LS-SVM (Least Square Support Vector Machine) to establish axial elongation and radial thermal yaw about motorized spindle thermal error model which results demonstrated that the axial (in the X and Y directions) and radial accuracies were improved greatly through experiments.Both the papers did not consider the small deformation of the ceramic materials.Grossi [13] studied the specific cutting force coefficients in the entire speed range by a single milling test through aluminum coefficients traditional approach to validate the efficiency and effectiveness of the proposed technique, which is not adaptive in high-temperature environment.Werner [14] employed the electromagnetic model with an electromagnetic moment to calculate radial and angular electromagnetic stiffnesses, the mass moment of inertia, and the gyroscopic effect of the rotor, but magnetic radical force is not considered.Liu et al. [15] proposed a model including Hertz contact force and unbalanced electromagnetic force to study the dynamic properties of the spindle system.
e studies [16][17][18] developed Maxwell stress tensor nonlinear harmonic methods and evaluated the vibration influence caused by UMP (Unbalanced Magnetic Pull).However, the magnet characteristics of different materials were ignored.Some traditional metal dynamic models could not be better suited for the ceramic spindle.
However, Oliver et al., Mironov et al., and Khazdozian et al. [19][20][21] have analyzed some magnetic characteristics of ceramic materials and reverse magnetic effect on the ceramic spindle, and bearing has not been investigated.In this paper, we give a general method that considers reverse magnetic effect that can be used to predict the vibration characteristics of the ceramic spindle systems.Based on the application of Timoshenko beam unit, a dynamic model is established in this paper.e relation between the dynamic stiffness, gyroscopic moment, and centrifugal force change is analyzed using the bearing raceway control theory.
e relationship between the magnetic density and electromagnetic radial force of ceramic materials is studied based on Maxwell's theory.
e radial force of the reverse magnetic and paramagnetic effects and the contact force of bearing Hertz on the dynamics of spindle are discussed.e complex dynamic mechanism of a high-speed ceramic spindle system is revealed through a ceramic spindle model simulation and multimodal test.
e ceramic spindle model can predict the magnetic field and dynamic characteristics well in some special extreme situations.

Dynamic Model of a High-Speed Ceramic Spindle
2.1.Spindle System Beam Element Model.Ceramic spindle bearings are installed on the ceramic shaft.A Timoshenko beam element is used to establish the dynamic model.e rotating shaft of the ceramic spindle is equivalent to stimulate vibration in the beam under the action of multiple factors, such as centrifugal force of the bearing and electromagnetic radial force [22].e model uses the Hamiltonian principle to calculate the kinetic energy, potential energy, and damping work of the beam element.
According to the damping and stiffness matrix, the influence of the electromagnetic radial force and bearing contact stress of the spindle system vibration is obtained.e discrete model of the ceramic spindle beam element is shown in Figure 1. e front bearing and rear bearing are equivalent to a spring and damping, which are K n1 , C n1 and K n2 , C n2 , respectively.
e discrete model assumes that the xyz axis is the inertial reference frame.Every axis unit is equivalent to 5 degrees of freedom: three directions of translational freedom are x, y, and z (u, v, w) and the degrees of freedom of rotation around the y and z axes (θ z and θ y ).Beam element nodes can be represented by point P in ball coordinates, φ is the point P which rotated about the axis of rotation, and d T is the distance between the point P and the center of the section circle.e spindle model of Timoshenko beam is shown in Figure 2. e displacement of point P is represented by q, q � (u, v, w, θ z , θ y ) T .e point P coordinate is Point P generates kinetic energy, strain energy, and external force work as the beam element changes between the instantaneous t 1 and t 2 moments.T and V are, respectively, shown in Figures 2(a) and 2(b).According to Hamilton's principle, formula (2) is derived similarly to the one reported in [23]: where ), θ z is the rotation angle of the x-y section, θ z c z + (z]/zx), and θ y is the rotation angle of the x-z section, θ y c y − (zω/zx).
Considering the formulas of motion of each axis, the dynamic model of the spindle is obtained: e contact stress between the ball and the inner ring is nonlinear.According to the Hertz contact theory, the ball contact point inner surface with the bearing inner and outer ring curvature center is positive and that with the outer where ρ i and ρ o are the total curvature of the bearing inner and outer rings with the contact point of the ball; ρ ir1 and ρ or1 are the normal curvature of the bearing inner and outer rings with the contact point of the ball, ρ ir1 2 cos θ/(d bc − d i ) and ρ or1 2 cos θ/(d o − d bc ); ρ ir2 and ρ or2 are the tangential curvature of the bearing inner and outer rings with the contact point of the ball, ρ ir2 −1/r ir , ρ or2 −1/r or ; and ρ b1 , ρ b2 are the normal and tangential curvature of bearing with the contact point of the ball, ρ b1 ρ b2 2/d b : Bearing friction is produced by centrifugal force acting on the ball at high speed.
e load and displacement of bearing under combined loading are shown in Figure 4. Figure 4(a) clearly shows that the gyroscopic moment can balance the friction generated between the ball and the contact of the inner ring and outer ring raceway surface.Figure 4(b) re ects the relative displacement relation of inner and outer rings under the conditions of some force and gyroscopic moment.According to the Jones raceway control theory, ball balance formula (8) is derived similarly to the one reported in [24].
e bearing inner ring under internal load is in a state of balance, and the bearing inner ring balance formula is where F cj is the centrifugal force of ball,

Magnetization Model of High-Speed Ceramic Spindle
Model. e stator of spindle is hot-pressed on the shell.e rotor is driven by a variable frequency magnetic eld to rotate at high speed.e magnetic potential distribution of spindle stator and rotor is shown in Figure 5. e two points α, β produce the z point magnetic potential.
where k pm is the pitch factor of m harmonic k pm sin m(y/τ)90 Air gap tangential magnetic eld density B θ is given by It can be seen on the surface of the stationary rotor that the spindle magnetic eld is orthogonal to the stator and  Shock and Vibration rotor surface.Even if eccentricity occurs, its value is ignored.e spindle coil is placed on the stator, and the air gap density of the radial magnetic eld is approximated by the Fourier transform: where χ m is the medium of magnetic susceptibility; χ c is the ceramic magnetic susceptibility 10 −5 cm 3 /mol; χ Fe is the steel magnetic susceptibility 14000 cm 3 /mol; μ r is the relative permeability; and μ 0 is the air permeability.e stator and rotor magnetic media exist from Ampere's loop theorem: where H 1 is the air magnetic eld intensity vector; H 2 is the ceramic magnetic eld intensity vector; and S is the stator and rotor contact area, S 2πrΔl; t → is the tangential unit vector; and Δl is the length of AB, AB CD. e tangential components of the air gap magnetic density on both sides are continuous.erefore, the magnetic eld density after magnetization B 1 has a linear relation: e initial static eccentricity is caused by installation error and stator deformation of the spindle structure.Electromagnetic vibration is one of the main factors of spindle vibration.e actual eccentricity is mainly composed of initial static eccentricity and vibration eccentricity [25,26].e rotor is eccentric relative to the stator and an eccentricity diagram of the spindle is shown in Figure 6. e formula ( 15) is derived similarly to the one reported in [27].
e air gap δ varies with the angle c.As the air gap is not uniform, the air gap magnetic permeability will change periodically with rotor rotation to generate additional harmonic magnetic elds.Harmonic magnetic elds interact to produce low-order electromagnetic radial forces [28].According to the Maxwell principle, the e ect of an air gap magnetic eld on a spindle stator inner surface produces a radial electromagnetic force, which is linear with the square of air gap magnetic eld density.e radial electromagnetic force of the spindle F ra (θ, t) is derived similarly to one reported in [29].
e spindle air gap magnetic density is where B(θ, t) is the air gap magnetic density, λ(θ, t) is the air gap magnetic permeance, and f(θ, t) is the air gap magnetic potential.

Synthetic Potential of Stator and Rotor f(θ, t).
When spindle is running at high speed, the gap magnetic potential is directly related to the spindle current, pole logarithm, and number of stator and rotor slots.e synthetic potential of the stator and rotor is where f 0 (θ, t) is the fundamental wave magnetic potential, f v (θ, t) is the stator harmonic magnetic potential, and f u (θ, t) is the rotor harmonic magnetic potential.
Considering the bearing damping e ect of the rotor system, the fundamental wave synthesis magnetic potential, stator winding harmonic magnetic potential, and rotor winding harmonic magnetic potential are obtained.
where I s is the stator current, I r is the rotor current, and I 0 is the synthesis current, I 0 I s + I r .

Air Gap Permeance λ(θ, t).
e spindle stator and rotor are provided with di erent number of tooth slots and tooth forms.
e air gap permeance is equivalent to where Λ 0 is the xed part of air gap permeance, k 1 λ k 1 is the air gap harmonic magnetic permeability caused by stator slotting, k 2 λ k 2 is the air gap harmonic magnetic permeability caused by rotor slotting, k 1 k 2 λ k 1 k 2 is the air gap harmonic magnetic permeability caused by both stator and rotor slotting interaction.e air gap permeability is a ected by the shape, number of tooth slots, concentricity, and symmetry of the stator and rotor.
e gap permeability is expanded by Fourier transform, keeping the rst two items in the decomposition item: e stator and rotor gap magnetic densities ignore factors such as a full magnetic circuit, and the slot a ects their interaction with each other.Substituting formulas (18) and (20) in formula (17) obtains where ε is the relative eccentricity, ε d rq /R s − r r .Substituting formula (22) in formula (16) gives the radial force of spindle:

Analysis of Dynamic Influencing Factors of the Spindle System
Bearing Hertz contact force and electromagnetic radial force will be produced during high-speed rotation of the spindle.
To study the e ects of the reverse magnetic e ect and bearing Hertz contact force on spindle dynamics, we combine every shaft segment of the Timoshenko beam unit with the bearing motion formulas and magnetic eld formulas of spindle.Formulas ( 9) and ( 23) are substituted in formula ( 4) to obtain the dynamic formula (24) under the 6 Shock and Vibration action of the bearing contact force and electromagnetic radial force.e dynamic formula of the spindle is solved by an iterative method.Formulas ( 9), (23), and (24) of the ceramic spindle model show that there is a direct relationship between the vibration and speed of the spindle, centrifugal force of the bearing, pretension force, gyroscopic moment, air gap, electromagnetic radial force, and magnetic field density.e structural parameters of spindle and the length of the rotating shaft discrete system, respectively, are shown in Tables 1 and 2.
where M d is the element rotational mass matrix, M is the element translational mass matrix, n is the number of shafts, K d is the element dynamic stiffness matrix, and F(t) is the spindle distribution force vector that considers electromagnetic radial force and bearing contact force.

Influence of Electromagnetic Radial Force on the Spindle.
Consider that the rotor rotates at a high speed in the stator, according to the electromagnetic model of the spindle.e stator and rotor are symmetrically distributed, so 1/4 is taken to analyze the electromagnetic field of the spindle.e conditions of the simulations are that the shell of the ceramic spindle is 170 mm, the air gap of stator and rotor is 0.3 mm, the motor pole logarithm is 4, the inner diameter of the stator is 80 mm, the outer diameter of the rotor is 79.4 mm, and bearing preload force is 400 N. e speed range, for which each increment is 1000 rpm, is from 9000 rpm to 18000 rpm.e magnetic density of the ceramic and metal spindle is shown in Figure 7. e radial force of the ceramic and metal spindle is shown in Figure 8. Figure 7 shows that the magnetic density of the spindle fluctuates with time.e main reason for this is that the spindle cuts the magnetic induction line under the highspeed frequency converter to form an eddy current.e magnetic force line is bent to form a closed magnetic circuit.e magnetic force line is limited by the air gap magnetic density.
e magnetic pole reaches the maximum saturation state.e magnetic pole is affected by the saturation state to show the skin effect and adjacent effect.
ere is a phenomenon of demagnetization that the current movement in the stator tends to the radial direction of the outer circle in the stator.
e magnetic density fluctuation degree of the ceramic spindle is far less than the metal spindle.Comparing the ZrO 2 ceramic spindle with the 40Cr metal spindle, the ceramic materials are the reverse of the magnetic materials.A ceramic material's permeability magnetic is only 10 −12 , which is much smaller than that of a metal material.e ceramic magnetic permeability is too low, and there is no magnetic force line passing through the ceramic spindle rotor.e magnetic force line inside the stator and rotor decreases so that the diffusion flux is suppressed.e air gap magnetic potential produces demagnetization effect to reduce eddy current loss.However, due to the high permeability of metal 40Cr rotor, under the magnetization effect, the internal molecules in the magnetic medium can be arranged in an orderly manner and the magnetic force line closes quickly, and the eddy current effect of the stator and rotor is obvious.e magnetic density fluctuation frequency of the ceramic reverse magnetic material is less than that of the metal paramagnetic material.
e magnetic density distribution of the ceramic spindle is less than that of the metal spindle.
Figure 7(a) shows the maximum magnetic density of the ceramic spindle stator pole gets to 7.5 μT at 0.047 s, which is less than the maximum magnetic density that the metal spindle gets to (17.5 μT at 0.033 s). e maximum magnetic density of the ceramic spindle is lower than that of metal spindle.
e variable frequency current of the spindle changes the flux direction of the rotor flux chain.Figure 7(b) shows the magnetic density distribution trend of spindle shell which the ceramic spindle shell can get to 0.011 μT at 0.0232 s and the metal spindle gets to (0.157 μT at 0.0218 s) are similar to the stator pole.
e magnetic density of shell is much smaller than stator pole due to the increasing distance, which results in loss of magnetic lines.
According to Lenz's law, the molecular magnetic moment in the ceramic magnetic medium is stored in the opposite direction of the external magnetic field.e reverse induction electromotive force is generated and an additional antimagnetic field phenomenon occurs.e change speed magnetic density of the stator magnetic pole slows down.
e magnetic density amplitude at the magnetic pole decreases; the metal spindle has the variable frequency current of the spindle that is consistent with the rotor flux chain.e magnetic field density of the magnetic pole increases under the paramagnetic action.

Shock and Vibration
Based on the magnetic density values of the metal and ceramic spindles, the radial forces of the metal and ceramic spindles can be calculated by substituting it in formula (24).Figure 8 shows that the maximum unbalanced radial force of the metal spindle can reach 139.96N at 0.048 s. e maximum unbalanced radial force of the ceramic spindle at 0.025 s can reach 25.6 N. e radial force distribution of the whole ceramic spindle is smaller than that of metal spindle.
e radial force uctuation of the ceramic spindle is less than that of metal spindle.e magnetization rate of the selenium metal shaft is di erent with the 40Cr metal rotor.e unevenness in magnetization is the main reason that the radial force of the metal spindle uctuates so greatly.
e eccentric radial force of the ceramic spindle is shown in Figure 9. e radial force of the spindle points to the minimum air gap between the stator and rotor, and the air gap directly a ects the radial force of the spindle.e radial force of the air gap eccentricity and the magnetic eld density of the stator and rotor increase with the increase in air gap eccentricity.e normal air gaps of a 170HT30 metal spindle and 170SD30 ceramic spindle are 0.3 mm. e radial force of the uneccentric air gap of the ceramic spindle is 24.4256N, and the radial force of 0.05 mm is 31.3097N. e maximum eccentric air gap reaches 0.1 mm. e minimum air gap between the stator and rotor is only 0.2 mm, at which the maximum radial force at the minimum air gap is 35.9246N. e radial force of electromagnetic imbalance is 50% higher than that of a noneccentric electromagnetic force.Eccentricity causes uneven air gap distribution and leads to a rapid increase in magnetic eld density at the   Shock and Vibration magnetic poles.A small air gap leads to a rapid increase in bypass flux density, rapidly forming closed magnetic lines of force.Diffusion of the magnetic field density is inhibited, and the peak magnetic pole is rapidly increased.e radial force has a great impact on the vibration of spindle.e eddy current phenomenon inside the stator and rotor is accelerated, the air gap magnetic density is inhibited, and the skin effect gradient is increased.Compared with a normal air gap, the magnetic pole amplitude of an eccentric air gap decreases significantly to the surrounding area, which is the main factor why the fluctuation degree of eccentric radial force is greater than that of normal operation.

Spindle-Bearing Dynamic Stiffness
Effect.According to formulas ( 8) and ( 9) for raceway control theory of a spindlebearing system, the bearing correction coefficient is adjusted for the characteristics of a ceramic Si 3 N 4 bearing to establish a ceramic bearing model and a metal bearing model, respectively.e conditions of simulation are radial dynamic stiffness, gyroscopic moment, and inner and outer ring contact angles of the bearing calculated under the working conditions of a preload force of 300-500 N and speed of 3000-30000 rpm.
Figure 10 shows the influence of preload change on the radial dynamic stiffness of metal and ceramic bearings.e radial stiffness of the rotating speed bearing is approximately decreasing with an increase in speed and increasing with an increase in pretension.At the same time, the preload force can hinder the change of radial stiffness with velocity.e radial stiffness of the ceramic bearing is greater than that of the metal bearing.Figure 10(b) shows that the radial dynamic stiffness of the ceramic bearing reaches a maximum value of 2.252 × 10 8 under the working condition of a speed of 3000 rpm and preload of 500 N. Since the speed increase is affected by the centrifugal force of the rolling body, there is a nonmonotone decreasing law and the bearing stiffness has a "softening effect."e trend of metal bearings is the same as that of ceramic bearings at low speed and medium speed.
Under the condition of low speed and a large preload, Figure 10(a) shows that the maximum value is 1.676 × 10 8 .
e turning point occurs when the speed exceeds 18000 rpm under the action of the centrifugal force of the bearing.When the metal bearing speed is above 18000 rpm, the centrifugal force of the rolling body increases.e radial dynamic stiffness of metal bearings is greatly affected by the contact expansion deformation between the outer ring raceway and ball, and the contact angle with the external load decreases to lead to the bearing operation mechanism changing.e decrease in inner ring stiffness is less than the increase in outer ring stiffness.For the ceramic bearing, the characteristics of large material stiffness lead to the deformation of outer ring channel being small.e change degree of the decrease in the stiffness of the inner ring is greater than that of the increase in the stiffness of the outer ring.
e radial stiffness of the ceramic bearing tends to decrease, and the bearing stiffness softening effect is obvious.
Figure 11 shows the comparison of different parameters of metal and ceramic bearings.e bearing operating speed range is 3000-30000 rpm, and the preload force is 400 N. Figure 11(a) shows that as the friction torque and centrifugal force between the ball and the bearing raceway with the ceramic bearing and metal bearing speed increase, the gyroscopic moment increases significantly and shows a rapid upward trend.e high stiffness and low density of the ceramic bearing are the main reasons that the gyroscopic torque of the ceramic bearing is less than that of metal bearing.Figure 11(b) shows that as the outer ring contact angle of the metal and ceramic bearings increases with the speed, the inner ring contact angle decreases.e inner and outer rings of the high-speed rotating bearing are affected by the centrifugal force of the ball and the gyroscopic moment, which is changing the contact mechanism.e contact stress and amount of deformation between the bearing ball and the raceway of inner and outer ring are changed.
e metal bearing outer ring contact angle is increased from 15.44 °to 22.82 °. e ceramic bearing outer ring contact angle is increased from 15.04 °to 19.77 °. e contact angle of the metal bearing changes more with the increase of speed than that of the ceramic bearing.e ceramic bearing inner ring contact angle reduces to 4.63 °and is 52% of the reduction in the contact angle of the inner ring of the metal bearing.ere is an approximate linear relationship between the centrifugal force of the bearing ball, gyroscopic moment, and the contact angle of the bearing outer ring raceway.e ball mass of the ceramic bearing being smaller than that of metal bearing is the main factor that means the centrifugal force and gyroscopic moment of the ceramic bearing are less than those of metal bearing.At the same time, the high stiffness of the ceramic bearing reduces the inner and outer ring raceway deformation and improves the dynamic performance of the bearing.

Test Magnetic Field and Vibration of Spindle.
Oil and gas lubrication is adopted for spindle-bearing lubrication under the test conditions.e experimental conditions are that the input oil is 5 mL/h, air pressure is 0.25 MPa, cooling water temperature is 18 °, the cooling water flow rate is 5 L/min, and bearing preload force is 400 N. e speed range, for which each increment is 1000 rpm, is from 9000 rpm to 18000 rpm.A Japan HIOKI FT3470-52 magnetic field tester was used to conduct a magnetic density test of the spindle magnetic field for a 170SD30 ceramic spindle and 170HT30 metal spindle.Figure 12(a) shows the magnetic field density test of a ceramic spindle shell measured on two positions by an electromagnetic field tester.A and B points are, respectively, stator position and bearing position on the spindle shell surface which are shown in Figure 12(a).A German Polytec OFV-505 contactless laser vibration measurement system was used to measure the spindle end vibration based on the laser Doppler effect.Figure 12(b) shows the spindle end vibration velocity test of the ceramic spindle.
e basic parameters of the 170SD30 ceramic spindle and 170HT30 metal spindle are shown in Table 3.    10

Magnetic field tester
Shock and Vibration

Test Results and Analysis.
e magnetic eld density of the electric spindle is mainly concentrated in the magnetic pole of the stator and the radial shell region.e distance between the magnetic eld density and magnetic poles is reversely proportional.
e magnetic eld density of the outer stator shell is less than that of the magnetic pole.e spindle shell is tested by an electromagnetic eld tester that can re ect the variation law of the electromagnetic eld density between xed rotor magnetic poles.
Magnetic density comparisons of the ceramic and metal spindles in A and B positions are, respectively, shown in Figures 13(a) and 13(b).Under test conditions at 18000 rpm to measure the magnetic density, A point of the metal and ceramic spindles, respectively, are 0.09475 μT and 0.0089 μT; meanwhile, the B point are, respectively, 0.08231 μT and 0.00823 μT.It can be clearly seen that the magnetic eld density of the metal spindle shell is higher than that of ceramic spindle shell.e magnetic density of the B point is a little smaller than A point that is attributed to the loss of the magnetic induction line.e magnetic density error between the simulation and experiment of the ceramic and metal spindles in A point, respectively, is 0.0018 μT and 0.01095 μT and B point, respectively, is 0.00101 μT and 0.00799 μT. e prediction of the magnetization model can be high, precisely veri ed through the comparisons with di erent positions.
e magnetic density of the spindle uctuates under the action of an eddy current and reverse magnetism is shown in Figure 14.It shows the magnetic density time domain diagram of the ceramic spindle at 18000 rpm.It is consistent with the trend of magnetic density simulation diagram in Figure 7, which proves the accuracy of the electromagnetic model.
Spindle vibration is mainly composed of mechanical vibration and electromagnetic vibration.e vibration velocity frequency domain diagram for the shaft end of ceramic and metal spindles at 18000 rpm is shown in Figure 15.
e amplitude of the metal spindle in the low-frequency region and high-frequency region is larger than that of the ceramic spindle in the frequency domain diagram.is shows that the radial force of electromagnetic imbalance, friction force between the bearing and rotating shaft, centrifugal force, and other factors in the spindle system are the components that in uence the frequency distribution.One time the frequency of rotation mainly includes electromagnetic magnetic tension.Two times the frequency of rotation mainly includes the bearing centrifugal force and friction force excitation frequency.e energy of 1 time the frequency in the low-frequency region is greater than that in the middle-and high-frequency regions.
In Figure 15(a), the peak value of 1 time the ceramic spindle frequency is lower than that of 2 times the frequency.
is shows that the electromagnetic radial force of the ceramic spindle has less in uence on the vibration of the spindle than the Hertz force of the bearing.In Figure 15(b), the peak value of 1 time the frequency of the metal spindle is greater than that of 2 times the frequency of the metal spindle.It can be seen that the electromagnetic radial force of the metal spindle has a great impact on the vibration of the spindle.In Figures 15(a) and 15(b), the amplitude of 1 and 2 times the frequency of the metal spindle frequency domain with vibration shock energy is larger than that of the ceramic spindle frequency.
e in uence of vibration velocity on  e in uence of electromagnetic force and bearing force on the spindle is analyzed in the frequency domain.Under the working condition of 18000 rpm, the frequency domain is approximately 600 Hz, so a frequency range within 1000 Hz is selected.
Figure 16 shows the comparison with the diagram of spindle end vibration velocity test.It can be seen that the vibration speed of the shaft end presents an increasing trend with the increase in rotating speed.e vibration speed of the shaft end of the metal spindle is higher than that of the ceramic spindle. is indicates that the Hertz contact force and centrifugal force caused by the change of rotation speed parameters are the main factor a ecting the vibration.
As seen in Figure 16(a), without considering the dynamic characteristics of electromagnetic radial force, the simulation vibration velocity is roughly the same as the variation trend of spindle vibration.e vibration velocity of the shaft end calculated by the simulation is less than that of the test result, and there is a large error compared with the real value.
e amplitude di erence of vibration velocity between the metal spindle simulation and experiment is higher than that of the ceramic spindle.
From Figure 16(b), considering the dynamic characteristics of the traditional electromagnetic radial force, the simulation vibration speed of the shaft end of the metal spindle is slightly higher than the test result.However, the traditional model does not take into account the reverse magnetic e ect of ceramics.
e spindle end vibration amplitude calculated by theory is much larger than the test results, and the error even reaches 20% under the working condition of 14000 rpm, which cannot e ectively re ect the vibration characteristics of the ceramic spindle.
From Figure 16(c), considering the magnetization e ect, the simulation speed of electromagnetic radial force dynamics is in good agreement with the test results.e electromagnetic radial force calculated after the magnetization of the ceramic reverse magnetic e ect and metal paramagnetic e ect has a high accuracy in obtaining the vibration simulation speed.In particular, the calculation accuracy of the ceramic spindle is greatly improved.Under the test condition of 18000 rpm, the vibration speed at the shaft end of the metal spindle was measured as 23.4 mm/s, and that at the shaft end of the ceramic spindle was measured as 18.65 mm/s.e simulation accuracy of the vibration velocity of the metal and ceramic spindle ends is reduced from the traditional dynamic magnetic model error values of 2.4% and 16.1% to 1.4% and 1.17%, respectively.It shows that the experimental results are consistent with the theoretical results, and the magnetization theoretical model can e ectively simulate the vibration speed under di erent working conditions.

Conclusion
e electromagnetic coupling dynamics model considered reverse magnetic e ect which can be used to predict magnetic eld and vibration on the ceramic spindle is presented.
e magnetic density of the ceramic spindle a ected by reverse magnetic e ect, Hertz contact force of the bearing,  Shock and Vibration and Maxwell stress of the rotor is analyzed.Furthermore, the e ects of the rotational speed, bearing sti ness, and the magnetic radial force on the dynamic characteristics of highspeed spindle are studied in detail.e following conclusions are drawn: (1) e actual working conditions of a high-speed rotating eccentric electromagnetic unbalanced radial force and bearing contact force excitation of a spindle system are studied.e coupling dynamics model considered reverse magnetic e ect of a ceramic spindle was established.e simulation results are basically consistent with the experimental results.
A more realistic simulation of speed and bearing contact force on spindle vibration has veri ed that the prediction of magnetization dynamic magnetic coupling model has high accuracy.
(2) Based on the Maxwell principle, the electromagnetic model of a ceramic spindle is established for the reverse magnetic characteristics of ceramic materials.e characteristics of ceramic reverse magnetic materials, electromagnetic magnetic density, and radial force generated by various working conditions are studied.e electromagnetic radial force of the ceramic spindle is smaller than that of the metal spindle.Compared with the experimental results, the simulation results of magnetization dynamic magnetic coupling model have smaller errors.
(3) e basic characteristics of a bearing, such as centrifugal force, gyroscopic moment, and contact angle of the inner and outer ring, are coupled together to establish a ceramic high-speed angular contact ball bearing model.It is calculated that the dynamic

Shock and Vibration 13
stiffness of a metal bearing is less than that of a ceramic bearing.With the increase in speed, the stiffness of the ceramic bearing appears to "soften."When the rotating speed reaches more than 18000 rpm, the centrifugal force of the metal bearing, gyroscopic moment, and ball contact angle of the inner ring increase significantly.e coupling effect changes the bearing operating mechanism to lead the stiffness to rise.

Nomenclature
T: Kinetic energy of spindle segment W: Damping work of spindle segment V: Deformation energy of spindle segment S: Cross-sectional area θ z : Rotation angle of the x-y section θ y : Rotation angle of the x-z section a: Inner diameter of rotating shaft cross section b: Outer diameter of rotating shaft cross section A: Magnetic

Figure 5 :
Figure 5: Magnetic potential distribution of spindle stator and rotor.

Figure 4 :
Figure 4: Load and displacement of bearing under combined loading: (a) ball load; (b) the relative displacement relation of inner and outer ring.

Figure 7 :Figure 8 :Figure 9 :
Figure 7: Magnetic density of the metal and ceramic spindles at 18000 rpm: (a) the stator pole of spindle; (b) the shell of spindle.

Figure 10 :Figure 11 :
Figure 10: Comparison of preload force of radial dynamic sti ness of (a) metal bearing and (b) ceramic bearing.

Figure 13 :
Figure 13: Magnetic density comparison of the ceramic and metal spindles.(a) A point: stator position on the spindle shell surface.(b) B point: bearing position on the spindle shell surface.

Figure 15 :Figure 14 :
Figure 15: Vibration velocity frequency domain diagram at the spindle shaft end at 18000 rpm: (a) ceramic spindle; (b) metal spindle.

Figure 16 :
Figure 16: Vibration velocity in the time domain between ceramic and metal spindles under 400 N bearing preload.(a) Without considering the dynamic characteristics of electromagnetic radial force comparison.(b) Considering the dynamic characteristics of traditional electromagnetic radial force comparison.(c) Considering the dynamic characteristics of the magnetization e ect.

Table 2 :
e length of the rotating shaft discrete system.

Table 1 :
e structural parameters of spindle.

Table 3 :
Basic parameters of spindle.