Next Article in Journal
A New Composite Control Strategy for an Astronaut Virtual Operation Training System Based on Cable-Driven Technology
Next Article in Special Issue
Practical and Rapid Motor Sizing Technique Using Existing Electrical Motor
Previous Article in Journal
Investigation and Optimization of Energy Consumption for Hybrid Hydraulic Excavator with an Innovative Powertrain
Previous Article in Special Issue
Experimental Validation of a Permanent Magnets Magnetorheological Device under a Standardized Worldwide Harmonized Light-Duty Test Cycle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Design Method of Post-Assembly Magnetizing Device with Field–Circuit Coupling Analysis

1
College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China
2
Zhejiang Provincial Key Laboratory of Electrical Machine Systems, Hangzhou 310027, China
3
Zhejiang Windey Co., Ltd., Hangzhou 311106, China
4
Bio Electro and Mechanical Systems Department, Université Libre de Bruxelles, 1000 Brussels, Belgium
*
Author to whom correspondence should be addressed.
Actuators 2023, 12(10), 383; https://doi.org/10.3390/act12100383
Submission received: 1 September 2023 / Revised: 29 September 2023 / Accepted: 3 October 2023 / Published: 11 October 2023
(This article belongs to the Special Issue Electromagnetic Actuators)

Abstract

:
Post-assembly magnetization can significantly simplify the manufacturing of the rotor of permanent magnet (PM) electrical machines. The optimization of the post-assembly magnetizing device using finite element analysis (FEA) is time-consuming; therefore, a globally optimal solution aiming at achieving an adequate magnetizing level and minimal energy consumption is difficult to achieve. In this paper, a field–circuit coupling analysis (FCCA) model is proposed to optimize the auxiliary stator-type magnetizing device for interior permanent magnet synchronous machines (IPMSMs). A reasonable simplification of the highly saturated magnetic circuit is made based on FEA results so that the magnetic equivalent circuit (MEC) model can be established. On the other hand, the eddy currents in the PMs are equivalent to an eddy current short-circuit; thus, by converting the field calculation into a circuit calculation, the time cost can be reduced significantly, which greatly improves the speed of multi-objective optimization of the magnetizing device with multiple degrees of freedom. A V-type IPMSM is taken as a study case, and its post-assembly magnetizing device is optimized with the proposed method. FEA and experimental results show that the PMs are fully magnetized, while the required energy consumption is greatly reduced when compared with an existing magnetizing device.

1. Introduction

Interior permanent magnet synchronous machines (IPMSMs) are widely utilized thanks to their outstanding advantages of high torque density and high efficiency. In the conventional manufacturing process, the permanent magnets (PMs) are magnetized before being assembled in the rotor. Strong magnetic forces then occur, making the process complicated, time-consuming, and risky of rotor component damage and personnel injury [1,2].
Post-assembly magnetization provides a solution to these problems. Unmagnetized magnets are inserted in the rotor, and then the latter are magnetized integrally with a dedicated magnetizing device. Due to the simple process and high safety with the unmagnetized magnets, post-assembly magnetization attracts growing attention.
Numerous post-assembly magnetization-related literature works have been published. Ferrite magnets can be easily magnetized with the machine stator winding, with the magnetization completeness depending on the PM shape [3]. The advantage of magnetization with the normal stator and its winding is that no additional equipment is required. For such a method, the magnetizing current is limited by the machine stator structure and, in particular, the mechanical strength of the winding. Thus, when rare-earth PMs are adopted, a much higher field strength is required and the stator needs to be especially designed to withstand a sufficient current for magnetization. Different stator winding currents for magnetization are studied, and a magnetization field of 2.5 times the coercivity is found to be sufficient for the full magnetization of rare-earth PMs [4]; then, a rare-earth surface-mounted PM synchronous machine (SPMSM), which is post-assembly magnetized with its stator, is designed. In [5], to simplify the design process, the stator winding is simplified to a filamentary conductor with the same number of turns, which represents the center of the magnetomotive force (MMF), and the relative permeability of the highly saturated iron core is approximated to that of the air; then, the parameters of the machine winding are optimized to obtain a sufficient magnetizing field, while the influence of the core structure is disregarded. However, the stator winding needs to be stronger when the large magnetizing current is considered, which is in conflict with the machine volume, cost, and operation performance. Thus, several improved methods are proposed to satisfy the magnetization requirement of IPMSMs and large-power SPMSMs, whilst those methods should be chosen according to the machine structure and magnetic field strength required for the magnetization of the magnets. For example, integrated magnetizing coils are directly wound around the PM poles for the magnetization of large-power SPMSMs [6]. For line-start IPMSMs, the utilization and influence of the squirrel cage and the copper bars are investigated [7,8].
However, a more common approach for the post-assembly magnetization of rare-earth IPMSMs is to use a dedicated device composed of a magnetizer and a magnetizing fixture, which is also called an auxiliary stator. The eddy currents generated during the magnetizing process have attracted the attention of researchers. A three-dimensional finite element (FE) model is established to analyze the eddy currents inside the PMs and metal sleeves [9]; the results show that the eddy currents may reduce the magnetizing field, resulting in the failure of the PMs’ magnetization. The impact of a soft magnetic material type of fixture is also compared [10]; the fixture using solid steel as the soft magnetic material consumes much more energy to generate a sufficient magnetizing field due to the induced eddy currents. Thus, the pulse-width is increased by optimizing the coil structure, in order to weaken the eddy currents [11,12]. For a line-start IPMSM, rotor bar dimensions are further researched, considering the induced currents in the squirrel cage [13]. The design process is further refined, where a preliminary design is made by deriving the required magnetic field energy when the magnets are fully magnetized, and then the coil turns and the capacitance of the magnetizer are focused on for optimization [14]. A similar two-step approach is also taken for the magnetization of V-type IPMSMs [15]; the winding structure is optimized with FE parametric modeling starting from a rough design, considering the influence of magnetizer capacitance and eddy currents in the PMs. A design method that has a fixed auxiliary stator structure during the optimization process is further compared with a design method with a fixed wire gauge [16]. The optimized structures of the magnetizing fixture obtained with the two methods are completely different, which illustrates that the parameters are mutually restrictive.
For these design processes, since the FEA is time-consuming, a few parameters are selected as degrees of freedom, and a local optimal solution is obtained to generate enough magnetizing field. However, the calculation method of the global optimal solution has not been proposed in the aforecited paper, as the minimization of energy consumption of the magnetizing device has not been considered. Therefore, an optimization method without FE parametric calculations is worth investigating.
In this paper, an auxiliary stator-type post-assembly magnetization device for a 1.3 kW, 15,000 r/min 6-slot 4-pole V-type IPMSM is investigated as a study case. Multi-parameter optimization is required to meet the design specifications and reduce the energy consumption, and a design approach based on the field–circuit coupling analysis (FCCA) is proposed. The magnetic equivalent circuit (MEC) models are improved to consider the highly saturated iron cores, and the eddy currents in the PMs are considered via an equivalent circuit model. The proposed design approach is used to optimize a magnetizing device, and the FE method is used to confirm its performance. A prototype of the optimal design is made and compared with an existing magnetizing device that is being used in a factory.

2. Investigated Magnetizing Fixture and Design Guidelines

2.1. Structure of Investigated Post-Assembly Magnetizing Device

The post-assembly magnetizing device is composed of a magnetizing fixture and a magnetizer. For a fractional-slot machine such as the one investigated, its stator winding is unable to generate a magnetizing field corresponding to the rotor poles, which results in a lot of spatial harmonics and makes the magnetization process much more complicated [17]. Thus, a magnetizing fixture with the same number of poles as the investigated rotor is adopted as an auxiliary stator, which offers a convenient solution. For the ease of understanding, the auxiliary stator will be referred to as the magnetizing fixture for the remainder of this article. Unmagnetized PMs are assembled in the rotor core of the IPMSM and then magnetized as a whole. The parameters of the investigated rotor are shown in Table 1. A quarter cross-section of the auxiliary stator and the rotor is shown in Figure 1a. The coils of the auxiliary stator are connected to the magnetizer, whose power circuit structure is shown in Figure 1b and the current pulse is generated. When thyristor SCR1 is turned on, the capacitor is charged to a high voltage through the transformer and the diode rectifier, with SCR1 turning off automatically at the end of the charging. Then, when thyristor SCR2 is turned on, the capacitor is discharged, serving as a pulse power supply. The power diode in reverse parallel is used to ensure that the capacitor energy is completely drained without oscillation.

2.2. Design Guidelines

For Nd-Fe-B magnets, it is recognized that a peak flux density value (reached at some instant during the magnetization) of 3 T is sufficient [2,14]. In this paper, the region at the bottom of the V-type PMs is focused on due to the difficulty to magnetize. In this paper, the doubled time that the charging current takes to rise from zero to the peak is further referred to as the magnetization pulse width T M , which should be sufficiently long for the domain rotation of the PMs.
The permeability of the iron core is approximated as the one of air during the magnetization pulse in some papers. However, for higher accuracy, the B–H curve of the FeSi laminations needs to be considered. Due to the rapid change of the magnetizing field, eddy currents are induced in the PMs, weakening the magnetizing field. The eddy currents in the iron cores are neglected here because they are much smaller than those of magnet blocks. The electrical conductivity of the Nd-Fe-B material is assumed to be 0.625 MS / m .
The airgap between the auxiliary stator and the investigated rotor should be as small as possible to strengthen the magnetizing field, and 0.5 mm is selected as the airgap length. For the magnetizing winding, it is imperative to use thickly insulated copper wires to withstand the high currents. Additionally, employing epoxy resin curing for stabilizing the winding is crucial to prevent deformation resulting from the stress of the current, particularly at the winding end coils. These requirements limit the potential use of the machine stator winding for magnetizing PMs, considering both the manufacturing cost and their impact on the machine performance.
The number of turns per coil N is restricted by the wire gauge and the slot dimensions. To make full use of the slot space, the number of coil layers N layer is:
N layer = K s 8 l p l 0 D w
where the operator represents rounding down, K s is the coefficient to adjust the slot fill factor considering main insulation, l p is the pole length, l 0 is the length of the pole shoe, and D w is the wire diameter considering inter-turn insulation; see Figure 1a. Due to the non-parallel slot walls, the number of turns in each layer is generally different; thus, the total number of turns is given by:
N = n = 1 N layer K s π 1 2 p θ p 360 r si + l 0 + n 1 D w D w
with p the number of pole pairs, θ p the pole width angle, and r si the inner radius of the auxiliary stator.
The device is supposed to fully magnetize the PMs while minimizing the energy consumption W M . As the magnetizer is completely discharged after the magnetization, the energy consumption is equal to the energy initially stored in the capacitor:
W M = 1 2 C U C 2
with C the magnetizer capacitance, and U C the discharging voltage.
It has been verified via simulation that a specific peak value of flux density can be achieved with either a smaller capacitance but higher voltage or a larger capacitance but lower voltage. However, the former would lead to a lower pulse width T M and a lower copper loss during the rising of the current. Therefore, these magnetizers have advantages in W M compared with the latter. However, T M has a minimum value, as mentioned above. It is necessary to ensure a sufficiently high T M while minimizing W M .
Six parameters remain as degrees of freedom for the optimization: θ p , l p , yoke thickness t y , D w , C, and U C . The first four parameters define the structure of the auxiliary stator along with the inductance and resistance of its winding. Pulse width T M and the completeness of magnetization depend on all six parameters, and the best parameter combination is searched based on the FCCA method.

3. Equivalent Models and Field–Circuit Coupling Analysis

3.1. Magnetic Equivalent Circuit Model of Rotor

The IPMSM rotor with V-type magnets has a complex geometry. The distances from each part of the magnets to the outer surface of the rotor are different. Thus, it is difficult to establish the MEC model directly, and FEA is used to aid in the modeling.
The outer surface of the one-eighth rotor is divided into six segments as shown in Figure 2. The dashed line indicates the dividing line of the flux, which enters the rotor at each segment when the core is saturated. The flux through each segment is denoted by φ r 1 to φ r 6 . Taking the right radial boundary of the model as the zero-potential surface, the average magnetic potential drop is denoted by F r 1 to F r 6 . According to the rotor geometry, the rotor is simplified to several reluctances connected in parallel, which are denoted by R rotori , R PM , and R Barrier in Figure 3. As the permeability of the PMs is much smaller than the one of the unsaturated silicon steel, the reluctances of the PMs and the magnetic barriers are especially denoted to clarify the magnetic circuit under different saturation conditions. Each of the reluctances is defined by the cross-section area S R and length l R :
F R = f φ R / S R · l R
with F R the potential drop on the reluctance, φ R the flux through the reluctance, and f the H ( B ) function of the different materials. Then, the potential drop of each segment is calculated as:
F r 1 = F R r o t o r 1 + F R PM F r 5 = F R r o t o r 5 + F R PM F r 6 = F R r o t o r 6 F R PM = F R Barrier
The shape of the PMs and the magnetic barriers is regular so that their dimensions are readily substituted into the reluctance formula as S R and l R , while the irregular shape of the iron core complicates the parameter determination. FEA is used to calculate F r j and φ r j of each segment under different currents. F r j is obtained by integrating H from the zero-potential surface to the rotor surface, and φ r j is obtained by integrating B on the segment of the rotor surface. Then, the FEA results are obtained and parametric regression is adopted to calculate the undetermined parameters. In this way, the MEC model is established.
The influence of the auxiliary stator structure on the rotor MEC model is also considered. The middle point of the third surface segment is selected as the test point, see the red dot in Figure 2. Flux φ r 3 is replaced by the normal magnetic flux density B r 3 , and the B r 3 F r 3 curves with different pole lengths and pole widths are compared in Figure 4. Obviously, the pole width has a significant influence on the curve, whereas the effect of pole length can be neglected. Thus, it is needed to calculate S R and l R under different pole widths. Then, the curves of these two parameters as a function of pole width are obtained via interpolation.

3.2. Magnetic Equivalent Circuit Model of Stator

The structure of the auxiliary stator is relatively simple and is defined with several degrees of freedom. In Figure 5, the magnetic circuit of the auxiliary stator is saturated and divided into six branches, with the latter categorized into two types. The first type of branches closes through the airgap, stator pole, and yoke, while the second type closes through the airgap, slot, and yoke. The permeability of the coils in the slot is assumed the same as air so that the two parts are merged. Each branch is simplified to several reluctances and a source of MMF connected in series as shown in Figure 3. The formula of reluctance is the same as the rotor MEC model, and the source of magnetomotive force M M F = K N N I , with K N the proportion of the half coil enclosed by the branch, N the number of turns per half coil, and I the value of the current. For the first type of magnetic circuit, K N is equal to 1.
In the traditional MEC method and with a normal level of core saturation, the second type of branches is neglected. For the study at hand, with the peak flux density in the stator core exceeding 5 T during the magnetization, it is necessary to consider the flux branches through the air. When the core is highly saturated, the flux density distribution is significantly different from the unsaturated condition, especially in the stator yoke and slot. Thus, several coefficients are proposed to make the MEC model suitable for both situations.
In the stator yoke and the air outside the stator, the flux circumferentially passes the symmetrical boundary. To simplify the model, the magnetic field strength on the boundary H boundary is assumed to decay as an inversely proportional function with the radius of the flux lines:
H boundary ( x ) = r si + l p r si + l p + x H boundary 0
with x the distance from the inner radius of the yoke, and H boundary 0 the magnetic field strength on the inner radius of the yoke. Then, the flux φ yoke passing the boundary through the yoke reads:
φ yoke = l rotor 0 t y f sl H boundary ( x ) d x
with f sl the H ( B ) relation of the laminations, and l rotor the rotor stack length. The air outside the stator deserves attention as an external circuit when the core is saturated. A thick enough section of air circuit is used to calculate the flux passing through the air φ outair :
φ outair = l rotor t y t y + t air μ 0 H boundary ( x ) d x
where μ 0 is the permeability of vacuum, and t air the adopted empirical value of the air thickness, with the latter determined as one-eighth of the outer circumference of the stator:
t air = π r si + l p + t y / 2 p
The air outside the stator is equivalent to an additional yoke. The total equivalent thickness t y eq satisfies:
l rotor 0 t y eq f sl H boundary ( x ) d x = φ yoke + φ air
The sum of the flux flowing through each branch is equal to the flux passing through the boundary:
j = 1 6 φ s j = φ yoke + φ air
with φ s j the flux of the jth stator magnetic branch. Then, the yoke width occupied by each branch is calculated based on the distribution of H boundary , and the yoke reluctances lof the branches are obtained.
In the stator slots, the flux passes in parallel to the pole and then becomes curved with increasing saturation level. Figure 6a illustrates a magnetic branch when the core is saturated. The gray line is the cross-section of the branch at the interface of the yoke and slot, and the black dot is the center of the branch region in the slot. There is a deviation from the dot to the perpendicular bisector of the cross-section, which means that the proportion of the coil enclosed by the flux is affected by the bending of flux lines. The deviation is measured to guarantee the accuracy of MMF and flux linkage. The curved flux line is assumed to be a circle arc, as shown in Figure 6b. The tangential flux density of the jth branch at the midpoint of the cross-section B t a n j is:
B t a n j = k = 1 j φ s k φ s j k = 1 6 φ s k B yoke 0
with B yoke 0 the flux density on the inner radius of the yoke,
B yoke 0 = f sl H boundary 0
and the normal flux density of the j th branch at the midpoint of the cross-section B n o r j is:
B n o r j = φ s j / S n o r j
with S n o r j the the cross-section area. According to the trigonometric relationship, the angular offset of the j th branch θ o f f s e t j is:
θ o f f s e t j = B n o r j 2 + B t a n j 2 B n o r j B t a n j l p / 2 r si + l p / 2
Therefore, the proportion of the coil enclosed by the j th branch is:
K N j = 1 , θ o f f s e t j + θ j > θ s θ o f f s e t j + θ j / θ s , θ o f f s e t j + θ j < θ s
where θ j is the angle corresponding to the perpendicular bisector of the j th cross-section, and θ s is the angle corresponding to the slot width, which is:
θ s = 90 θ p / 2
With the proposed improvements, the calculations of the MMF and the reluctances are suitable for different saturation conditions. Further, the results under different numbers of equivalent circuit branches are compared. We observe little difference when more than three branches are considered. Thus, it is recommended to contain three branches in the MEC model to reduce the calculation time. At least one branch of each type is required, and the type of the remaining one is determined according to the pole width.

3.3. Equivalent Circuit of Eddy Currents in Permanent Magnet

In this section, the analytical method is taken to evaluate the reaction of PM eddy currents to the applied magnetizing field. Then, an RL series circuit is used as the equivalent eddy current circuit to approximate the analytical results. Based on Maxwell’s equations, the basic field equation in a PM is:
2 H PM = μ PM γ PM H PM t
where H PM is the magnetic field strength in PMs; μ PM and γ PM are the permeability and electrical conductivity of PMs. To simplify the analysis, several assumptions are made:
  • The applied magnetic field is uniform and perpendicular to the surface of the PMs.
  • μ PM and γ PM are constant.
Thus, a two-dimensional (2D) field on the x–y plane of local coordinate is used as the solution domain, as shown in Figure 7. The H PM has the component in the z-direction only, denoted by H PM , which satisfies:
2 H PM ( x , y , t ) x 2 + 2 H PM ( x , y , t ) y 2 = μ PM γ PM H PM ( x , y , t ) t
The boundary and initial conditions are:
H PM ( ± l PM / 2 , y , t ) = H 0 ( t ) H PM ( x , ± w PM / 2 , t ) = H 0 ( t ) H PM ( x , y , 0 ) = 0
with l PM the length of PMs, w PM the width of PMs, and H 0 the magnetic field strength of the applied magnetizing field. The method of separation of variables [18] is used to solve the equations. H PM is expanded into a Fourier series, and the unit impulse response function h PM is derived as:
h PM ( x , y , t ) = m = 1 , 3 , 5 , n = 1 , 3 , 5 , K m n τ m n cos m π x l PM cos n π y w PM e t / τ m n
with K m n the magnitude of the Fourier coefficient, τ m n the decay time constant:
K m n = 16 m n π 2 sin m π 2 sin n π 2
τ m n = μ PM γ PM π 2 m 2 l PM + n 2 w PM
The Fourier series in (21) is composed of sinusoidal plane waves of order m in the x direction and order n in the y direction. With the increase in the spatial order, the sinusoidal waves have smaller decay time constants.
A half-sinusoidal magnetic field with the pulse width of T and the amplitude of H mag is assumed to be applied to the PMs. By the Duhamel’s integral,
H PM ( x , y , t ) = 0 t H mag sin π τ T h PM ( t τ ) d τ
Therefore,
H PM ( x , y , t ) = H mag sin π t T H mag m = 1 , 3 , 5 , n = 1 , 3 , 5 , K m n C m n ( t ) cos m π x l PM cos n π y w PM
with C m n the time coefficient, which varies with the orders:
C m n ( t ) = π T τ m n · π T τ m n sin π t T + cos π t T e t / τ m n 1 + π T τ m n 2
The first term of H PM in (25) is the applied field, and the second term is the reaction field of eddy currents in the Fourier series. The first two terms of C m n are the steady state components, which are sine waves that lag behind the applied field. And the third term is the transient component, which decays exponentially.
Further, the flux counteracted by the reaction field of PMs φ re is:
φ re = φ mag PM sin π t T μ PM H PM ( x , y , t ) d x d y
where:
φ mag PM = μ PM l PM w PM H mag
Therefore,
φ re = 16 π 2 φ mag PM m = 1 , 3 , 5 , n = 1 , 3 , 5 , C m n ( t ) m 2 n 2
The total counteracted flux can be decomposed into the spatial harmonics of different orders, whose amplitude decrease rapidly with the increasing orders.
Equivalent eddy current circuit is used to simulate the reaction field of the PMs. A unit circuit consisting of a resistance R u and an inductance L u will generate an equivalent eddy current i u as well as the counteracted flux φ ^ ru when an induced electromotive force E u is applied:
L u d i u d t + R u i u = E u
φ ^ ru = L u i u
Therefore,
φ ^ ru = φ mag PM π T L u R u π T L u R u sin ω t + cos ω t e R u L u t 1 + π T L u R u 2
Comparing the counteracted flux obtained by the analytical method and unit-equivalent eddy current circuit, φ ^ ru has a similar formula structure with C m n . Thus, any spatial harmonic of the counteracted flux can be represented by a unit equivalent circuit with a resistance R u m n and an inductance L u m n , which satisfies:
L u m n R u m n = τ m n
Despite the unit-equivalent circuit accurately representing the counteracted flux of each spatial order, an equivalent circuit composed of m × n unit circuits is complicated and the computation will be time-consuming. An aggregate equivalent circuit with the same structure of the unit circuit is proposed as an alternative:
L e d i e d t + R e i e = E e
with L e the inductance of the aggregate equivalent circuit, i e the equivalent eddy current, R e the resistance which varies with the frequency of applied magnetic field, and E e the induced electromotive force. In the aggregate equivalent circuit, the eddy currents are equivalent to flowing along the edge of the PMs; thus, the E e and L e are:
E e = d d t φ mag PM sin π t T
L e = μ PM l PM w PM / t PM
with t PM the thickness of the PM.
Due to the multiple spatial harmonics contained in the total counteracted flux, (33) is not suitable for the aggregate equivalent eddy current circuit. Thus, R e is calculated with parametric regression to make the results obtained with the equivalent circuit close to the analytical results.
In Figure 8, φ re under different frequencies is analyzed with (29), and the φ ^ re calculated with the equivalent eddy current circuit and the corresponding R e is compared. When the frequency of the applied field is close to 100 kHz, the spatial harmonics of low orders, such as the 3rd and 5th, cannot be ignored and have waveforms different from the fundamental one. Thus, there is a deviation between the analytical method and equivalent circuit. The accuracy of the aggregate equivalent circuit can be improved by increasing the number of unit equivalent circuits at such a frequency. However, as the magnetization process requires an adequate pulse width of the order of a millisecond, the magnetizing field time-harmonics at the frequencies close to 100 kHz are insignificant and the aggregate equivalent circuit has a sufficient accuracy for such a scenario. The curve of R e with the frequency of the applied field is plotted in Figure 8e. The variation of R e below 100 kHz is relatively little, which means that R e can be regarded as a constant.

3.4. Field–Circuit Coupling Analysis

The magnetization process is actually the discharge of a magnetizing circuit as shown in Figure 9, where i ( t ) is the instantaneous magnetizing current, u C ( t ) is the instantaneous capacitor voltage, with initial value U C , and R coil is the resistance of the coil. The latter is related to the number of turns N and the wire gauge as follows:
R coil = 4 p N l m / γ copper S wire
with γ copper the electrical conductivity of copper, l m the stack length of the magnetizer, and S wire the cross-section area of the wire.
The MEC models of the rotor and auxiliary stator are connected to calculate the flux linkage. The stator circuit segmentation is decided by pole width, and then rotor segmentation is adjusted accordingly. The parameters of the adjusted rotor MEC model are obtained via parameter smoothing based on the rotor MEC model studied in Section 2.
Then, the flux φ s j in each branch considering different magnetizing currents I is calculated with the Newton–Raphson method. The flux linkage Ψ is:
Ψ ( I ) = 2 p i K N j N φ s j ( I )
from which, the resulting current-dependent inductance L ( I ) = Ψ ( I ) / I is obtained. The eddy currents are generated by the varying flux through the PMs φ R PM ; this is similar to the induced current in the short-circuited secondary winding of a transformer. The inductance corresponding to φ R PM is denoted by L 1 :
L 1 ( I ) = 4 p K N PM N φ R PM ( I ) I
with K N PM the average proportion of coils surrounded by the magnetic circuits through the PMs. Then, the inductance corresponds to the remainder of the part L 2 ( I ) = L ( I ) L 1 ( I ) .
R e , L e , and i e are the reduction values of the equivalent eddy resistance, inductance, and current on the coil side, which are converted from the equivalent eddy current circuit:
R e = 4 p K N PM N 2 R e
L e = 4 p K N PM N 2 L e
Based on the magnetizing circuit, the equations of the magnetization process after turning on the thyristor SCR2 are divided into two stages. When the capacitor is discharging:
d L 1 i i e · i i e d t + d L 2 i · i d t = u C R coil i d L 1 i i e · i i e d t d L e i e · i e d t = R e i e C d u C d t = i
After the capacitor is completely discharged, the power diode in reverse parallel turns on and thyristor SCR2 turns off:
d L 1 i i e · i i e d t + d L 2 i · i d t = R coil i d L 1 i i e · i i e d t d L e i e · i e d t = R e i e u C = 0
Since the inductance varies with the current, the equations are discretized and solved iteratively. The results obtained are verified with the FE method. Figure 10 shows the operating parameters calculated with FEA and FCCA. The accuracy and efficiency of the FCCA method are demonstrated by how closely the FCCA results match those of the FEA.

4. Optimization Design of Auxiliary Stator-Type Magnetizing Device

The stages of the design method proposed in this paper are illustrated in Figure 11. The PM equivalent eddy current circuit and rotor MEC model are simplified from the rotor to be magnetized. The stator MEC model is established according to the auxiliary stator structure, which is decided by multiple degrees of freedom. And the transient analysis is carried out according to the FCCA model with the software MATLAB R2021a. The FCCA greatly reduces the time consumption for magnetizing simulation, providing the foundation for the multi-objective optimization of the auxiliary stator-type magnetizing device with multiple degrees of freedom.
The six degrees of freedom mentioned in Section 2 are considered for optimizing the device. To fully magnetize the PMs, the peak flux density B PM _ min at the bottom region of the V-type PM is estimated with the rotor MEC model and is required to be above 3 T. Meanwhile, a pulse width of the order of millisecond is also essential. Provided that the magnetizing ability is sufficient, the lower the volume, cost, and energy consumption, the better. Thus, the objective function f obj and the constraints are:
m i n f obj = b 1 r so + b 2 N S wire + b 3 W M s . t . B PM _ min 3 T T M 1 ms
where b 1 , b 2 , and b 3 are the weights corresponding to the indicators, r so is the indicator of device volume, N S wire is the indicator of copper quantity of the fixture which is a significant part of the device cost, and W M is the energy consumption, which can be calculated by (3).
A genetic algorithm (GA) is used for the global searching of the optimal solution. A population of 40 potential solutions, also called individuals, is randomly generated as the first generation. The performance of each individual is calculated with the FCCA method, and a penalty term f pen is introduced to transform the optimization problem into an unconstrained one:
f pen = b 4 3 T min B PM _ min , 3 T 1 + b 5 1 ms min T M , 1 ms 1
where b 4 and b 5 are the weights corresponding to the indicators. Therefore, the fitness function S is obtained:
S = S 0 f obj r so , N , S wire , W M f pen B PM _ min , T M
where S 0 is the fitness baseline, which is used to adjust the exploration and convergence speed of the GA. The tournament selection is adopted to form a new generation, where the best individual always makes it into the next generation and the remaining individuals undergo the roulette-wheel selection. Crossover and mutation steps are then performed with the probability of 0.4 and 0.01, respectively.
The convergence of the optimization process is shown in Figure 12, where the fitness gradually increases and stabilizes at a value with the 44th generation. The individuals of the last 10 generations, the Pareto frontier, and the selected solution are shown in Figure 13 with the three objectives involved. Then, an individual of the Pareto frontier with high fitness value, suitable wire gauge, and appropriate magnetizer capacitance is selected as the final solution. The parameters of the obtained auxiliary stator-type magnetizing device are shown in Table 2.
The optimized design is compared to an empirical design that is being used in a factory. The empirical design has the same outer diameter r so and the number of turns N as the optimized design, but differs in the shape of the auxiliary stator and the arrangement of the coils. The winding resistance and inductance of the two device designs are similar. The magnetizer capacitance adopted by the empirical design is also consistent with the optimized design.
Ansys Maxwell is used to simulate the magnetization of both designs, and the discharging voltage is set to 700 V. In Figure 14a, the flux density of the PMs is above 3 T in most areas, which indicates the PMs are close to saturation. Compared with Figure 14b, the optimized design has been better magnetized. The reason is that the pole width of the empirically designed stator is equal to the angle of the top of V-type PM, which means that slots must be deeper to accommodate the coils. Consequently, the yoke is relatively thin and becomes more saturated during the magnetization. Figure 15 compares the flux density distribution at the PMs’ middle line, which is drawn with black dotted lines in Figure 14. The narrower pole adopted in the optimized stator helps to concentrate the magnetic field into the middle, which has a better magnetizing effect in the bottom region.

5. Experiment

An empirically designed device and a prototype of the optimized design are tested with the 2 mF source shown in Figure 16, which is composed of eight 250 μ F capacitors. The magnetizing winding of the fixture is stabilized with epoxy resin curing to prevent the deformation of the coil end turns.
An increasing initial voltage is set to the magnetizer to find the minimal energy for full magnetization. Since the repeated unsaturated magnetization of the same PM will increase the remanence, new unmagnetized PMs produced in the same batch are used for each magnetization in the experiments. After each completion of magnetization, all eight PMs are taken out from the rotor, placed in the center of a Helmholtz coil, and quickly pulled out. A voltage is induced in the Helmholtz coil as the flux linkage changes. Thus, the magnetic flux is obtained by integrating the induced voltage E H ( t ) , and the magnetic moment M of a PM is calculated as follows [19]:
μ 0 M = α · E H ( t ) d t
with α the coil constant of the Helmholtz coil.
As illustrated in Figure 17, the magnetic moment M of the PMs magnetized by the optimized design reaches its saturation at 600 V, which is slightly less than the 700 V estimated by the proposed circuit–field coupling analysis. To evaluate the magnetization completeness with the two post-assembly magnetizing devices, PMs that are fully magnetized in the factory are used as a reference; these are also further referred to as pre-magnetized PMs. The black dotted line indicates the magnetic moment M on the pre-magnetized PMs, and the coincidence indicates that the post-assembly magnetization achieves the same effect of the pre-magnetization. The results of the empirical design are also tested where the PMs are unsaturated until 1000 V. Comparatively, the voltage required by the optimized design is 60% of that by the empirical design. According to (3), the energy consumption W M of the optimal design is 360 J, while that of the empirical design is 1000 J, which means that 64% of the energy is saved. Additionally, the aging of the insulation paper and the wires is slowed down thanks to the decrease in electro-thermal stress, thereby preventing insulation breakdown and increasing the device’s lifespan [20,21,22].
Considering the uneven distribution of the magnetizing field throughout the PMs, the surface flux density of the PMs is sampled with a Gauss meter. Four positions on the front and back sides of the top and bottom regions of the V-shaped PM are selected as sampling points, as shown in Figure 18a. PMs magnetized with both devices are tested and the surface flux density is compared in Figure 18b. The surface flux density at points B and D is lower than that at points A and C, indicating that the bottom region of the PMs is more difficult to magnetize when the discharge voltage is insufficient; several values measured at the bottom region are negative due to the field generated by the top region being closed through to the bottom region. Therefore, special attention should be paid to the bottom region for the magnetization of V-type IPMSMs. Comparing the result of the two designs, the empirical design has a similar saturation voltage to the optimized design at points A and C, while the difference in saturation voltage at points B and D is much larger. Consequently, as validated in the optimal design, the narrow pole is advised for the magnetizing devices of V-type IPMSM, whilst the multi-parameter optimization design should be performed.
Moreover, the flux density on the surface of the magnetized V-type PM rotor is measured. A rotor with the pre-magnetized PMs is used as the reference, and the rotors with unmagnetized PMs are magnetized with both prototypes at 2 mF and 600 V for comparison. The surface magnetic flux density of the three rotors is shown in Figure 19. The rotor surface flux density of the optimized prototype almost coincides with that of the pre-magnetized PMs, indicating that the PMs achieve saturation with the proposed design. In contrast, the flux density of the empirical design is slightly lower. The results of the rotor surface flux density are consistent with the experiments above.

6. Conclusions

A design method of an auxiliary stator-type magnetizing device based on the FCCA is proposed. The highly saturated iron core and the eddy currents inside the PMs are considered with equivalent models. It is verified that the results of the FCCA have an excellent agreement with the FEA and experimental results. The fast computing ability of the FCCA model allows for rapid global searching optimization of the magnetizing device. A multi-objective optimized solution is successfully obtained with a GA, resulting in a 40% reduction of the capacitor voltage and a 64% reduction of the energy consumption when compared with the empirical design. Therefore, this design method may become a valuable tool to accelerate the discovery of optimal solutions, which is not limited to V-type IPMSMs, and which can help the designer to clarify the intricate interactions between the parameters of the rotor and the magnetizing device.
In this paper, the optimization method of a V-type IPMSM is proposed, with a focus on the magnetic field. In further research, the coil stress and temperature rise caused by the large currents will be studied. To this end, a multi-physics model combining magnetic, heat transfer, and mechanical aspects will be established, and more objectives and constraints will be considered in the optimization process.

Author Contributions

Conceptualization, J.-X.S. and Z.-A.Z.; methodology, Z.-A.Z. and X.-F.Q.; software, Z.-A.Z. and L.Y.; validation, Z.-A.Z. and Y.-C.W.; formal analysis, Z.-A.Z.; investigation, Z.-A.Z.; resources, J.-X.S. and Y.-C.W.; data curation, Z.-A.Z.; writing—original draft preparation, Z.-A.Z.; writing—review and editing, J.G. and Z.-A.Z.; visualization, Z.-A.Z.; supervision, J.-X.S.; project administration, J.-X.S.; funding acquisition, J.-X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China grant number 2022YFB3403100.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PMPermanent Magnet
FEAFinite Element Analysis
FCCAField–Circuit Coupling Analysis
IPMSMInterior Permanent Magnet Synchronous Motor
MECMagnetic Equivalent Circuit
SPMSMSurface-Mounted Permanent Magnet Synchronous Motor
MMFMagnetomotive Force
FEFinite Element
GAGenetic algorithm

References

  1. Dorrell, D.G.; Hsieh, M.F.; Hsu, Y.C. Post Assembly Magnetization Patterns in Rare-Earth Permanent-Magnet Motors. IEEE Trans. Magn. 2007, 43, 2489–2491. [Google Scholar] [CrossRef]
  2. Li, L.; Huang, X.; Wu, A.; Thompson, P.; Ranze, R.; Yu, H.; Laskaris, T.; Amm, K. Development of a Large Bore Long Pulse Magnet at GE Global Research Center. IEEE Trans. Appl. Supercond. 2006, 16, 1672–1675. [Google Scholar] [CrossRef]
  3. Kim, K.S.; Park, M.R.; Kim, H.J.; Chai, S.H.; Hong, J.P. Estimation of Rotor Type Using Ferrite Magnet Considering the Magnetization Process. IEEE Trans. Magn. 2016, 52, 8101804. [Google Scholar] [CrossRef]
  4. Hsieh, M.F.; Hsu, Y.C.; Chen, P.T. Analysis and Experimental Study of Permanent Magnet Machines with In-Situ Magnetization. IEEE Trans. Magn. 2013, 49, 2351–2354. [Google Scholar] [CrossRef]
  5. Hsieh, M.F.; Hsu, Y.C.; Dorrell, D.G. Design of Large-Power Surface-Mounted Permanent-Magnet Motors Using Postassembly Magnetization. IEEE Trans. Ind. Electron. 2010, 57, 3376–3384. [Google Scholar] [CrossRef]
  6. Wang, Q.; Ding, H.; Zhang, H.; Zhou, J.; Wang, G.; Lv, Y.; Li, L. Research of Post-Assembly Magnetization of Large-Power Surface-Mounted Rare-Earth Permanent Magnet Machines With Integrated Magnetizing Winding. IEEE Trans. Appl. Supercond. 2020, 30, 5204505. [Google Scholar] [CrossRef]
  7. Ho, S.L.; Li, H.L.; Fu, W.N. A Post-Assembly Magnetization Method of Direct-Start Interior Permanent Magnet Synchronous Motors and Its Finite-Element Analysis of Transient Magnetic Field. IEEE Trans. Magn. 2012, 48, 3238–3241. [Google Scholar] [CrossRef]
  8. Fu, W.N.; Chen, Y. A Post-Assembly Magnetization Method for a Line-Start Permanent-Magnet Motor. IEEE Trans. Appl. Supercond. 2016, 26, 0602604. [Google Scholar] [CrossRef]
  9. Lv, Y.; Xia, D.; Wang, Q.; Li, L. Post-assembly magnetization of high speed PM rotor with metallic sleeve. In Proceedings of the 2017 20th International Conference on Electrical Machines and Systems (ICEMS), Sydney, Australia, 11–14 August 2017; pp. 1–5. [Google Scholar] [CrossRef]
  10. Angara, R.C.S.B.; Hsu, K.W.; Villar, P.J.; Sheth, N.K. Influence of Soft Magnetic Material Type in Fixture Components on the Magnetization of Bonded Neo Magnet and Motor Performance. IEEE Trans. Magn. 2018, 54, 2101305. [Google Scholar] [CrossRef]
  11. Lv, Y.; Yang, Y.; Xia, D.; Xu, B.; Peng, T.; Xiao, H.; Li, L. Saddle-Shaped Post-Assembly Magnetization Coil for a 300 kW 2-Pole High-Speed Permanent Magnet Rotor. IEEE Trans. Appl. Supercond. 2020, 30, 5206705. [Google Scholar] [CrossRef]
  12. Kwon, S.J.; Lee, B.H.; Kim, K.S.; Jung, J.W. Design Process of Post-Assembly 3-Times Magnetizer for 10-Poles of Flux Concentrated Rotor Considering Eddy Current Effect. IEEE Access 2023, 11, 34476–34485. [Google Scholar] [CrossRef]
  13. Lee, C.; Kwon, B.; Kim, B.T.; Woo, K.; Han, M. Analysis of magnetization of magnet in the rotor of line start permanent magnet motor. IEEE Trans. Magn. 2003, 39, 1499–1502. [Google Scholar] [CrossRef]
  14. Lee, C.K.; Kwon, B.I. Design of post-assembly magnetization system of line start permanent-magnet motors using FEM. IEEE Trans. Magn. 2005, 41, 1928–1931. [Google Scholar] [CrossRef]
  15. Wang, Q.; Ding, H.; Zhang, H.; Lv, Y.; Guo, H.; Li, L. Study of a Post-Assembly Magnetization Method of a V-Type Rotor of Interior Permanent Magnet Synchronous Motor for Electric Vehicle. IEEE Trans. Appl. Supercond. 2020, 30, 5206205. [Google Scholar] [CrossRef]
  16. Negahdari, A.; Toliyat, H.A. Post-assembly magnetization of rare-earth permanent magnet materials in permanent magnet assisted synchronous reluctance motors. In Proceedings of the 2017 IEEE International Electric Machines and Drives Conference (IEMDC), Miami, FL, USA, 21–24 May 2017; pp. 1–6. [Google Scholar] [CrossRef]
  17. Hsieh, M.F.; Lien, Y.M.; Dorrell, D.G. Post-Assembly Magnetization of Rare-Earth Fractional-Slot Surface Permanent-Magnet Machines Using a Two-Shot Method. IEEE Trans. Ind. Appl. 2011, 47, 2478–2486. [Google Scholar] [CrossRef]
  18. Bland, D.R. The Method of Separation of Variables. In Solutions of Laplace’s Equation; Bland, D.R., Ed.; Springer: Dordrecht, The Netherlands, 1961; pp. 15–28. [Google Scholar] [CrossRef]
  19. Trout, S. Use of Helmholtz coils for magnetic measurements. IEEE Trans. Magn. 1988, 24, 2108–2111. [Google Scholar] [CrossRef]
  20. Ramu, T.S. On the Estimation of Life of Power Apparatus Insulation Under Combined Electrical and Thermal Stress. IEEE Trans. Electr. Insul. 1985, EI-20, 70–78. [Google Scholar] [CrossRef]
  21. Sun, P.; Sima, W.; Yang, M.; Wu, J. Influence of thermal aging on the breakdown characteristics of transformer oil impregnated paper. IEEE Trans. Dielectr. Electr. Insul. 2016, 23, 3373–3381. [Google Scholar] [CrossRef]
  22. Polyakov, D.A.; Nikitin, K.I.; Novozhilov, A.N.; Kletsel, M.Y. Research of Electro-Thermal Aging Process of Cross-Linked Polyethylene. J. Phys. Conf. Ser. 2018, 1050, 012061. [Google Scholar] [CrossRef]
Figure 1. Auxiliary stator-type post-assembly magnetizing device: (a) Quarter cross-section of the rotor to be magnetized and the auxiliary stator. (b) Power circuit structure of magnetizer.
Figure 1. Auxiliary stator-type post-assembly magnetizing device: (a) Quarter cross-section of the rotor to be magnetized and the auxiliary stator. (b) Power circuit structure of magnetizer.
Actuators 12 00383 g001
Figure 2. One-eighth model of rotor flux distribution with saturated iron core.
Figure 2. One-eighth model of rotor flux distribution with saturated iron core.
Actuators 12 00383 g002
Figure 3. MEC model of rotor and auxiliary stator.
Figure 3. MEC model of rotor and auxiliary stator.
Actuators 12 00383 g003
Figure 4. Influence of auxiliary stator structure on the rotor MEC model: (a) B r 3 F r 3 curve with different pole widths and 4.4 mm pole length. (b) B r 3 F r 3 curve with different pole lengths and 23.5° pole width.
Figure 4. Influence of auxiliary stator structure on the rotor MEC model: (a) B r 3 F r 3 curve with different pole widths and 4.4 mm pole length. (b) B r 3 F r 3 curve with different pole lengths and 23.5° pole width.
Actuators 12 00383 g004
Figure 5. One-eighth model of auxiliary stator flux distribution when iron core is being saturated.
Figure 5. One-eighth model of auxiliary stator flux distribution when iron core is being saturated.
Actuators 12 00383 g005
Figure 6. Distribution of flux lines around the auxiliary stator slot: (a) One of the magnetic branches with saturated iron core. (b) Geometric relationship of θ o f f s e t j , B t a n j , and B n o r j .
Figure 6. Distribution of flux lines around the auxiliary stator slot: (a) One of the magnetic branches with saturated iron core. (b) Geometric relationship of θ o f f s e t j , B t a n j , and B n o r j .
Actuators 12 00383 g006
Figure 7. The 2D magnetic field within the x–y plane of PMs.
Figure 7. The 2D magnetic field within the x–y plane of PMs.
Actuators 12 00383 g007
Figure 8. Comparison between φ re and φ ^ re at different frequencies: (a) Reaction flux with f = 1 kHz and R e = 35.8 m Ω . (b) Reaction flux with f = 10 kHz and R e = 37.2 m Ω . (c) Reaction flux with f = 100 kHz and R e = 42.5 m Ω . (d) Reaction flux with f = 1 MHz and R e = 90.0 m Ω . (e) Variation of R e with frequency.
Figure 8. Comparison between φ re and φ ^ re at different frequencies: (a) Reaction flux with f = 1 kHz and R e = 35.8 m Ω . (b) Reaction flux with f = 10 kHz and R e = 37.2 m Ω . (c) Reaction flux with f = 100 kHz and R e = 42.5 m Ω . (d) Reaction flux with f = 1 MHz and R e = 90.0 m Ω . (e) Variation of R e with frequency.
Actuators 12 00383 g008
Figure 9. Equivalent electric circuit model considering the eddy currents.
Figure 9. Equivalent electric circuit model considering the eddy currents.
Actuators 12 00383 g009
Figure 10. Dynamic results obtained with FCCA and FEA: (a) Current with time. (b) Flux linkage with time. (c) Eddy currents losses with time. (d) Capacitor voltage with time.
Figure 10. Dynamic results obtained with FCCA and FEA: (a) Current with time. (b) Flux linkage with time. (c) Eddy currents losses with time. (d) Capacitor voltage with time.
Actuators 12 00383 g010
Figure 11. Stages and processes of the design method.
Figure 11. Stages and processes of the design method.
Actuators 12 00383 g011
Figure 12. Convergence of fitness with generation.
Figure 12. Convergence of fitness with generation.
Actuators 12 00383 g012
Figure 13. Individuals of last 10 generations, Pareto frontier, and selected solution.
Figure 13. Individuals of last 10 generations, Pareto frontier, and selected solution.
Actuators 12 00383 g013
Figure 14. Comparison of the two device designs when the capacitance is 2 mF and discharge voltage is 700 V. (a) Flux density map of empirical design. (b) Flux density map of optimized design.
Figure 14. Comparison of the two device designs when the capacitance is 2 mF and discharge voltage is 700 V. (a) Flux density map of empirical design. (b) Flux density map of optimized design.
Actuators 12 00383 g014
Figure 15. Peak value distribution of flux density in PMs along width direction.
Figure 15. Peak value distribution of flux density in PMs along width direction.
Actuators 12 00383 g015
Figure 16. Magnetizer and prototypes.
Figure 16. Magnetizer and prototypes.
Actuators 12 00383 g016
Figure 17. Flux generated by post-magnetized PMs with discharge voltage for empirical design and optimized design.
Figure 17. Flux generated by post-magnetized PMs with discharge voltage for empirical design and optimized design.
Actuators 12 00383 g017
Figure 18. Flux density of post-assembly PMs at different sampling points. (a) Positions of sampling points of a PM. (b) Comparison of the flux density at different sampling points under empirical design and optimal design.
Figure 18. Flux density of post-assembly PMs at different sampling points. (a) Positions of sampling points of a PM. (b) Comparison of the flux density at different sampling points under empirical design and optimal design.
Actuators 12 00383 g018
Figure 19. Measured flux density distribution at the rotor surface with pre-magnetized PMs, and with the PMs post-magnetized with the optimized and empirically designed devices.
Figure 19. Measured flux density distribution at the rotor surface with pre-magnetized PMs, and with the PMs post-magnetized with the optimized and empirically designed devices.
Actuators 12 00383 g019
Table 1. Parameters and ratings of the investigated IPMSM.
Table 1. Parameters and ratings of the investigated IPMSM.
Parameter NameParameter Value
Power1300 W
Speed15,000 r/min
Number of poles4
Rotor outer diameter42 mm
Rotor stack length46.2 mm
Magnet thickness2.8 mm
Magnet width11.2 mm
Iron core materialFeSi laminations
PM materialNd-Fe-B (N42SH)
Table 2. Parameters of the optimized auxiliary stator-type magnetizing device.
Table 2. Parameters of the optimized auxiliary stator-type magnetizing device.
Parameter NameParameter Value
Number of poles4
Air gap length0.5 mm
Auxiliary stator outer diameter72.4 mm
Auxiliary stator inner diameter43 mm
Auxiliary stator length46.2 mm
Pole angle43 degree
Pole length4.4 mm
Yoke thickness6.6 mm
Number of coil layers2
Number of turns per coil10
Number of parallel branches1
Wire gauge1.8 mm
Bare wire diameter1.5 mm
Magnetizer capacitance2 mF
Discharging voltage700 V
Auxiliary stator materialSilicon Laminations
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, Z.-A.; Wang, Y.-C.; Qin, X.-F.; Yao, L.; Gyselinck, J.; Shen, J.-X. Optimal Design Method of Post-Assembly Magnetizing Device with Field–Circuit Coupling Analysis. Actuators 2023, 12, 383. https://doi.org/10.3390/act12100383

AMA Style

Zhu Z-A, Wang Y-C, Qin X-F, Yao L, Gyselinck J, Shen J-X. Optimal Design Method of Post-Assembly Magnetizing Device with Field–Circuit Coupling Analysis. Actuators. 2023; 12(10):383. https://doi.org/10.3390/act12100383

Chicago/Turabian Style

Zhu, Zi-Ang, Yun-Chong Wang, Xue-Fei Qin, Lei Yao, Johan Gyselinck, and Jian-Xin Shen. 2023. "Optimal Design Method of Post-Assembly Magnetizing Device with Field–Circuit Coupling Analysis" Actuators 12, no. 10: 383. https://doi.org/10.3390/act12100383

APA Style

Zhu, Z. -A., Wang, Y. -C., Qin, X. -F., Yao, L., Gyselinck, J., & Shen, J. -X. (2023). Optimal Design Method of Post-Assembly Magnetizing Device with Field–Circuit Coupling Analysis. Actuators, 12(10), 383. https://doi.org/10.3390/act12100383

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop