1. Introduction
With the rapid development of space technology, the requirements for the working performance of spacecraft are increasing, especially for high-resolution earth observation, satellite docking, rapid response and other tasks. The magnetically suspended control moment gyroscope (MSCMG) realizes attitude adjustment through the momentum exchange output torque of the high-speed rotor. It has the characteristics of high control accuracy and high reliability. In recent years, it has been extensively studied and has become the attitude control actuator of many spacecraft [
1,
2]. The performance of the motor that provides the angular momentum for MSCMG affects the accuracy of attitude control. The hollow-cup brushless permanent magnet motor has the advantages of being lightweight and having a high power density and no cogging torque pulsation. Therefore, in order to improve the accuracy and reduce power consumption, MSCMG mainly uses hollow-cup motors. At the same time, it has been widely used in other aerospace fields [
3,
4,
5] and civilian fields [
6,
7]. At present, the hollow-cup motor with trapezoidal characteristics of air-gap magnetic field waveforms has the problem that the waveform does not match the ideal trapezoidal wave due to the problem of its own structure, which leads to torque ripples, which limit its high-precision rate servo application [
8].
In order to solve the above problems, improving the air-gap magnetic field waveform is an important way to reduce torque ripple [
9,
10], and some studies have optimized the air-gap magnetic field of the motor [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24]. In [
11,
12], the shape of PM was designed, and the equivalent surface current method was used to obtain the analytical calculation formula of the air-gap magnetic field when the outer arc of the tile-shaped PM is eccentric. In [
13], two kinds of eccentric PMs were designed to optimize the trapezoidal wave waveform of the permanent magnet synchronous motor; at the same time, an improved method for the air-gap magnetic field model was proposed. In [
14], PMs are optimized for eccentric design, and obtaining the eccentricity under the optimal solution by FEA proved that the design can effectively reduce the harmonic components of the air-gap magnetic field. In [
15], through the chamfering design of PMs, the waveform of the air-gap magnetic field was optimized. In [
16], a double-pole PM was proposed to optimize the air-gap magnetic field waveform and sawtooth torque.
In addition to optimizing the shape of PMs, interior motors mainly optimize the air-gap magnetic field waveform by optimizing the shape of the outer surface of the rotor. In [
16], four different rotor shapes were designed, and the feasibility of improving the performance of the motor by changing the shape of the rotor was proved through experiments and finite element analysis. In [
18], the outer diameter of the rotor was eccentrically designed, and the analytical solution of the air-gap magnetic field was obtained through conformal transformation. It was proved that the design could reduce the harmonic component and optimize the air-gap magnetic field waveform. In [
19,
20,
21], the rotor surface was slotted design. This design can effectively reduce the harmonic components of the air-gap magnetic field and optimize the trapezoidal wave waveform and the cogging pulsation. In [
22], a rotor shape of a permanent magnet synchronous motor was designed as an octagon; it was proved that the air-gap magnetic field waveform was improved and the harmonic content and ripple torque were reduced. In [
23,
24], through genetic algorithms and other algorithms, the motor structure was improved and the air-gap magnetic field waveform was optimized.
In summary, since the hollow-cup motor adopts the PMs surface mount structure, at present, there is only one method for optimizing the air-gap magnetic field waveform, that is, to change the PM shape. However, this involves problems of magnetization as well as manufacturing difficulties, resulting in an air-gap magnetic field waveform that cannot achieve the ideal waveform. In this paper, on the premise of not changing the shape of PMs, according to the characteristics that the hollow-cup motor is composed of inner and outer rotors and the relative position of the inner and outer rotors remains unchanged during operation, we propose a novel structure design of the inner rotor for the hollow-cup motor. Compared with changing the shape of PMs, this design avoids the magnetization and manufacturing problems of PMs. Firstly, the influencing factors of the air-gap magnetic field are analyzed through a mathematical model, and the model of inner rotor eccentricity is proposed and established. Then, the finite element model of the motor rotor is established, the trapezoidal wave quality under different eccentricity values is evaluated, and the optimal eccentricity value is obtained by analysis to make it as close as possible to the ideal waveform. Finally, a test prototype is made and through experiments and the effectiveness of the novel design is verified.
2. Mathematical Model of Motor Air-Gap Magnetic Field
A sectional view and a top view of the hollow-cup motor are shown in
Figure 1. It is mainly composed of an inner rotor core, an outer rotor core, permanent magnets (PMs) and a hollow-cup stator. The PM is attached to the inner side of the outer rotor. Taking the hollow-cup motor in MSCMG as an example, the inner and outer rotors are fixed by the rotor fixing ring and connected to the flywheel. The cup-shaped stator is wound with windings and fixed on the gyro room, and the winding part is placed in the gap between the PMs and the inner rotor.
The air-gap magnetic field of the motor is excited by the PMs, and the surface current is used to replace the PMs by the equivalent surface current method; the integral method is then used to obtain the magnetic field generated by the equivalent surface current [
25,
26,
27]. This paper mainly analyzes the magnetic field generated by the parallel magnetized magnets in the air gap.
The basic law of a constant magnetic field can be expressed by the following basic equation:
Since the divergence of the magnetic induction intensity B is always zero, a vector magnetic potential
is introduced into the curl magnetic field, such that:
For a constant magnetic field, it is usually specified that:
In the motor, the magnetic field can usually be transformed into a two-dimensional parallel plane magnetic field for processing, and the motor is a centrosymmetric rotating structure; the current direction is thus the Z direction in the polar coordinate system. Combined with the above formula, we obtain:
From Equation (2),
where
is the radial air-gap magnetic field and
is the tangential air-gap magnetic field.
The general solution of Equation (4) is:
Using the separation of variables method to solve this equation, let:
Bring Equation (7) into Equation (6), and separate the variables:
where
m is any constant, when
m = 0:
Due to the centrosymmetric structure of the motor, the magnetic field in the motor is a periodic function about the
; m is thus a positive integer. Therefore,
where
are constant, determined according to boundary conditions.
Then, according to the boundary conditions of the specific mathematical model, the constants in the air-gap magnetic field of the hollow-cup motor are determined. For the magnetic field produced by a single-ended wire at the center of the polar coordinate system, according to Ampere’s law:
In
Figure 2a, the inner radius of the outer rotor is
and the outer radius of the inner rotor is
. Place a single long straight conductor at
in the air gap of the motor. The vector magnetic potential
generated by the conductor passing through the ferromagnetic boundary at any point in the air gap is:
At the interface between the air gap and the inner rotor, the tangential component of the air-gap magnetic field is 0:
At the interface between the air gap and the outer rotor, according to the tangential component of the air-gap magnetic field, the boundary conditions are:
By solving the boundary conditions:
Substitute Equations (14) and (18) into Equation (11), the vector magnetic potential of a single long straight wire at any point in the air-gap magnetic field can be obtained:
where:
is the coordinate of point
relative to point
.
Figure 2.
Current carrying line in the air gap of a hollow-cup motor. (a) A current-carrying coil in the air gap; (b) a pair of current-carrying coils in the air gap.
Figure 2.
Current carrying line in the air gap of a hollow-cup motor. (a) A current-carrying coil in the air gap; (b) a pair of current-carrying coils in the air gap.
The position of a pair of current-carrying coil element sides A and B in the air gap is
in the space plane polar coordinate system, as shown in
Figure 2b. Treat the current-carrying coil as two oppositely directed currents and add the resulting magnetic potentials, and a mathematical model of the magnetic field produced by a pair of current-carrying coils can be obtained:
Considering the centrifugal force, the magnetic steel is usually attached to the outer rotor, so that
; then:
Considering the centrifugal force, the magnetic steel is attached to the outer rotor:
Then, the radial flux density of a pair of coils in the air gap is:
The superposition calculation of the coils in the multi-stage coil can obtain the expression of the magnetic flux density produced by the interaction of the current-carrying coil in the air gap by 2
p:
The equivalent surface current of the PMs is shown in
Figure 3:
According to the equivalent surface current method, the surface currents of AB and CD along the surface of the magnet are equal in opposite directions; thus,
The surface currents along the AD side and BC side of the PM surface are in opposite directions and equal in magnitude, which is:
where
is the angle of the PM,
is the angle of the AD or BC side and
is the coercivity.
Since the AB and CD sides are exactly the same,
The BC and AD edge-to-surface flow micro-elements are:
The equivalent coil current micro-element is , the radial length micro-element of the AB and CD sides of the PM is and the angle micro-element is .
It can be obtained from the above formula that the PM thickness , the PM inner diameter b, the outer rotor inner diameter , the inner rotor outer diameter and the number of pole pairs have an influence on the rotor air-gap magnetic field. When the number of pole pairs and the structure size of the motor is determined, the change in the air-gap magnetic field at any point in the air gap can be achieved through PM eccentricity or rotor eccentricity. However, it is difficult to magnetize and mechanize the PMs after eccentricity; the eccentric model of the inner rotor is thus designed in this paper.
Figure 4 shows the structure design of the inner rotor based on the eccentricity of the inner rotor. The eccentric structure design of the inner rotor is an eccentric design of the inner rotor outer diameter and the corresponding part of the outer rotor PMs. Move the center O of the outer diameter of the inner rotor to the position O’ in the opposite direction along the center line of the PM, keep the original end position of the outer diameter of the inner rotor unchanged and draw an arc with O’ as the center of the circle, which is the new outer diameter.
From the geometric relationship, the position A(
) of any point on the outer diameter of the inner rotor can be obtained, where:
In the above formula, is the eccentricity value and is half of the opening angle of PMs.
Combining
with Formulas (23) and (27)–(29), we can obtain the equivalent currents of p-pair PMs on the AB and CD sides, the BC sides, and the AD sides when the outer diameter of the inner rotor is eccentric. The magnetic flux density generated at the air gap (
) is
Then, integrate and sum the above formula to obtain the air-gap magnetic density produced by the PMs in the air gap:
From the above formula, when the size of the motor structure is determined, the air-gap magnetic field is only related to the eccentric value.
For trapezoidal waveforms, there is currently no effective evaluation index. In this paper, the trapezoidal wave waveform is evaluated based on the proportion of the top part τ, which is defined as:
where
is the part above 99% of the maximum value of the air-gap magnetic field
in one cycle of the air-gap waveform, and α is the half-cycle width of the air-gap waveform, as shown in
Figure 5.
3. Finite Element Simulation Analysis
This paper discusses the simulation of the motor model used in the laboratory. The model parameters are shown in
Table 1.
The motor model was constructed in the Ansoft Maxwell, and the above-mentioned motor model was simulated to analyze the influence of the inter rotor eccentricity on the air-gap magnetic field waveform. The maximum length of the elements was set to 30 mm and the maximum number of elements was set to 30,000. The mesh is shown in
Figure 6 and the magnetic map and magnetic line distributions are shown in
Figure 7.
Comparing the analytical results and the simulation results, the difference was about 1%, which proves the correctness of the model, as shown in
Figure 8.
A waveform analysis was carried out for different eccentricity values. The eccentricity range was 0~18 mm and the interval was 2 mm. The magnetic field waveform at the midpoint of the air gap under different eccentricities is shown in
Figure 9.
It can be concluded from
Figure 9 that, with the increase in the eccentricity value, the flat top part of the air-gap magnetic field waveform gradually widened. When the eccentricity value increased to 14 mm, the air-gap magnetic field waveform was distorted.
According to Formula (1), we calculated the proportion of flat top
under different eccentric values. When the eccentricity value was 0 mm, the
was 58.8° and
was 32.66%. The relationship between
and eccentricity is shown in
Figure 10.
increased significantly with the increase in eccentricity within a certain range. The highest increase was from 32.66% in the 0 mm eccentric state to 70.34% in the 12 mm eccentric state, an increase of 115.37%. As the eccentricity value continued to increase, there was a decrease. The reason is that with the increase in the eccentricity value, the width of the air gap increases and the relative magnetic resistance increases, resulting in a decrease in the amplitude of the air-gap magnetic field, making the peak value of the waveform smaller and thus closer to the flat top. When the eccentricity value is too large, the reduction in the amplitude of the magnetic density of the air gap caused by the large magnetic resistance will be too large, and the phenomenon of depression will appear. Therefore, it can be concluded that within a certain range, the flat top width of the air-gap waveform increases with the increase in the eccentricity value of the outer diameter of the inner rotor, and the air-gap magnetic field waveform is optimized; the eccentricity value of 12 mm is the optimal solution in the model mentioned in this paper.
A comparison of the air-gap magnetic field waveform before and after optimization under the eccentric value of 12 mm is shown in
Figure 11.
In summary, when the inner rotor was 12 mm eccentric, the half-cycle of the air-gap magnetic field was significantly flatter and the proportion of the flat top part was increased from 32.66% to 70.34%. This shows that the eccentric design of the rotor in the motor can effectively improve the waveform of the air-gap magnetic field: the proportion of the flat top part makes the air-gap magnetic field waveform closer to the trapezoidal wave, thereby improving the performance of the motor.
4. Prototype and Experiment
The prototype parameters are shown in
Table 1, and the prototype is shown in
Figure 12. It consists of inner/outer rotors and permanent magnets. The clamp fixes the inner and outer rotor cores and keeps the distance between them fixed. The permanent magnets are attached to the inner side of the outer rotor core. The prototype of the inner rotor under the optimal eccentricity value (12 mm) in the simulation and the prototype before optimal design were manufactured.
The experimental platform was mainly composed of fixtures, rotating platforms, magnetometers, brackets and control measurement devices, as shown in
Figure 13. The test process was to install the prototype, adjust the position of the magnetometer probe, set the measurement parameters, perform the test and record the test data. The detailed testing process is as follows:
- (1)
The clamping mechanism consists of a three-jaw centering chuck to fix the rotor and stator of the hollow-cup motor;
- (2)
The rotating motor is fixedly connected with the three-jaw chuck to drive the three-jaw chuck to rotate;
- (3)
The magnetometer uses the Hall effect to measure the magnitude and direction of the magnetic field, with a measurement accuracy of 0.001 T and a sampling of 36,000 points for each rotation of the prototype;
- (4)
The manipulator clamps the magnetometer and adjusts the magnetometer probe to the middle of the air gap;
- (5)
The measurement data is output as a two-dimensional array of angle points and magnetic field strength values.
Figure 13.
Air-gap magnetic field test experimental platform.
Figure 13.
Air-gap magnetic field test experimental platform.
The experimental test was carried out on the motor with an inner rotor eccentricity value of 12 mm, and the abnormal period data caused by the uneven magnetization of the PMs were eliminated. A comparison between the experimental results and the simulation results under 12 mm eccentricity is shown in
Figure 14:
When the eccentricity value was 12 mm, the was 123.68°, calculated by Formula (1); was 68.71% from Formula (1), which is different from the simulation value of 70.34% by −2.32%. The difference in was small, which verifies the correctness of the finite element simulation results.
The before the optimized design of the inner rotor was 66.99°; when the inner rotor eccentricity value was 12 mm, τ increased from 37.22% to 68.71%, which proves that the hollow-cup motor model based on the eccentricity of the inner rotor can effectively improve the shape of the air-gap magnetic field waveform and make it better adopt a trapezoidal shape. In order to achieve a more ideal waveform, the requirements for the machining accuracy of the eccentric inner rotor are high, but the current machining accuracy can meet the requirements. At the same time, compared with the design of PMs eccentric, it is more convenient for magnetization and manufacturing; the structure proposed in this paper therefore has a certain engineering application value.