Next Article in Journal
District Heating Challenges for the UK
Next Article in Special Issue
An Analysis of Non-Isolated DC-DC Converter Topologies with Energy Transfer Media
Previous Article in Journal
Comprehensive Hydraulic Improvement and Parametric Analysis of a Francis Turbine Runner
Previous Article in Special Issue
Robust Design Optimization with Penalty Function for Electric Oil Pumps with BLDC Motors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Analysis of a Novel Brushless Hybrid Excited Adjustable Speed Eddy Current Coupling

1
School of Electrical Engineering, Southeast University, Nanjing 210096, China
2
Jiangsu Magnet Valley Technologies Co., Ltd., Zhenjiang 212009, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(2), 308; https://doi.org/10.3390/en12020308
Submission received: 19 December 2018 / Revised: 13 January 2019 / Accepted: 17 January 2019 / Published: 19 January 2019

Abstract

:
A novel brushless hybrid excited adjustable speed eddy current coupling is proposed for saving energy in the drive systems of pumps and fans. The topology and operation principle of the coupling are presented. Based on the real flux paths, the fluxes excited by permanent magnet (PM) and field current are analyzed separately. A magnetic circuit equivalent (MEC) model is established to efficiently compute the no-load magnetic field of the coupling. The eddy current and torque are calculated based on the proposed MEC model, Faraday’s law, and Ampere’s law. The resultant magnetic fields, eddy currents, and torques versus slip speeds under different field currents are studied by the MEC-based analytical method and verified by finite element analysis (FEA). The copper loss, core loss, and efficiency were investigated by FEA. The analytically predicted results agree well with the FEA, and the analysis results illustrate that a good speed regulation performance can be achieved by the proposed hybrid excited control.

1. Introduction

The eddy current coupling forms a basic variable speed drive, which is widely applied to regulate the flows from large pumps and fans for a remarkable energy-saving effect [1]. With the development of permanent magnet (PM) material, a high torque density adjustable speed PM eddy current coupling (AS-PMECC) (Oregon State University, Corvallis, OR, USA) has been developed and applied in the field of variable speed drive [2]. Compared with other adjustable speed devices, such as gearboxes [3] and variable frequency drives [4], AS-PMECC has many advantages, such as lower sensitivity to environmental conditions, more reliable overload protection, and better energy-saving performance [5,6].
Based on the magnetized directions of PMs, AS-PMECC can be generally classified as axial flux [7] and radial flux couplings [8]. It is usually installed between prime motor and load, which commonly consists of a PM rotor (PMR), a conductor rotor (CR), and an additional mechanical manipulator (AMM) [6]. When a relative rotation happens between the two rotors, eddy currents can be induced in the copper sheet (CS) mounted on the iron core of the CR, which yields an electromagnetic torque from the interaction with the primary PM magnetic field [9,10,11,12]. For a given load torque, the speed of AS-PMECC can be adjusted by changing either the air-gap length or coupling area between the two rotors by the AMM [6,10].
In essence, the mechanical-based solutions usually realize speed regulation by adjusting the air-gap flux between the PMR and the CR. However, it is relatively hard to regulate the axial displacement between the two rotors, which requires a relative complicated design of AMM. Further, the complicated AMM is not reliable enough and takes more axial space, which is urgent to be simplified and improved for the energy-saving reconstruction of aged pump and fan systems. In addition, the local overheat problem produced by shifting a rotor is very serious in the conventional radial flux AS-PMECC [8,13]. Therefore, this paper proposes a novel brushless hybrid excited adjustable speed eddy current coupling (HE-ASECC) (Southeast University, Nanjing, China ) based on a hybrid excited concept to solve the above problem [14,15].
The novelty of the proposed coupling lies on the brushless hybrid excited (HE) control by an additional field excitation stator (AFES), which can provide a flexible air-gap field adjustment. The AFES just occupies some radial space rather than axial space, which is very suitable for the energy-saving reconstruction of aged pump and fan systems. Due to the absence of brush, slip ring, and complicated AMM, the magnetic field control of the HE-ASECC becomes highly simple and reliable. For illustrating the advantages of the HE-ASECC more clearly, it is compared with valve and baffling vane (VB), AS-PMECC, and variable frequency drive (VFD) in Appendix A. Due to the special geometry of the coupling, three-dimensional (3D) finite element analysis (FEA) is usually required to accurately compute its electromagnetic characteristics, but it is very time-consuming. For improving computational efficiency, an analytical model based on the magnetic equivalent circuit (MEC) method is proposed to analyze the electromagnetic characteristics of the coupling. The eddy current and torque are calculated and analyzed based on Ampere’s laws under an asymmetric magnetic field. Finally, the analytically predicted results are presented and compared with the FEA.

2. Structure and Principle

Figure 1a shows the structure of the studied HE-ASECC. It consists of an AFES, a PMR, and a CR, all of which are coaxial. Different from the AMM used in the conventional AS-PMECC, the AFES embedded with a toroidal field winding is located outside the PMR. Usually, the PMR and the CR are connected with the load mover and the prime motor, respectively. The PMR is characterized by two suits of axially parallel consequent-poles mounted on the inner surface of the PMR core that consists of two annular iron cores fixed together with an axial distance. The PMs mounted on the two annular iron cores are radial magnetized in opposite directions, respectively. The copper sheet (CS) is tightly mounted on the surface of iron core of the CR to provide paths for the induced eddy currents.
The field control principle can be explained by the magnetic fluxes passing through the PM and iron poles as shown in Figure 1b–d under different field currents. The PMs serve as constant magnetomotive force (MMF) sources, while the field winding acts as a changeable MMF source to perform a flexible air-gap flux adjustment. As a result, the fluxes produced by the field current can enhance or weaken the air-gap fluxes produced by PMs alone based on the directions and amplitudes of the applied field currents. Consequently, the slip speed between the PMR and the CR under a given load torque can be adjusted with the aid from HE control.

3. Analytical Model

In the design process of an electromagnetic device, a reliable and effective theoretical model is desirable for evaluating its electromagnetic characteristics. Generally, the theoretical models based on the magnetic vector potential [9,16] and the MEC method [17,18,19,20,21] are very popular for the design and analysis of electromagnetic devices. Considering the special 3D structure of the proposed HE-ASECC, the latter is adopted in this paper.

3.1. No-Load Magnetic Field Calculation

Based on the real fluxes shown in Figure 1b–d and the reluctance elements of MEC, the fluxes created by the PMs and the field current can be analyzed separately. In addition, it should be noted that the CR core and the stator core are communal for the MEC branches excited by the PMs and the field current. Figure 2a,b give the MEC models of the coupling excited by both the PMs and the field current when the slip speed between the two rotors is 0 rpm, respectively.
As shown in Figure 2a, the PM fluxes can be divided into two parts. One passes through the adjacent PM poles located on the two annular iron cores, the inner and outer air-gaps and the CS, and loops with the CR core and the stator core. The other begins with the PM poles and passes through the adjacent iron poles located on the same annular iron core, the inner air-gap and the CS, and loops with the PMR core and the CR core. As shown in Figure 2b, the fluxes created by the field winding pass through the outer air-gap, the PMR cores, the iron pole, and the PM-air region, including the PM poles and the air region between PM and iron poles located on the same annular iron core, the inner air-gap and the CS, and loop with the CR core. All elements of the MEC model are calculated in the cylindrical coordinate to permit higher accuracy. In addition, the symbols describing the structure parameters of the coupling are presented in its axial section and axial cross section shown in Figure 3a,b, respectively.
Assuming the iron parts of the coupling are unsaturated, the reluctances of the stator core, the PMR core, the iron poles, and the CR core can be ignored for simplifying calculation. The MMF sources are constituted by the PM poles and the field winding located in the AFES, which can be respectively expressed as
F P M = H c h P M
F f = I f N f
where If, Nf, Hc, and hPM represent the field current, the turns of the field winding, the PM coercive force, and the thicknesses of the PM poles, respectively. Applying Kirchhoff’s voltage (KVL) law to Loops 1–6 of the MEC model shown in Figure 2 yields [17,18]
[ R 11 R 12 0 R 21 R 22 R 23 0 R 32 R 33 ] × [ ϕ 1 ϕ 2 ϕ 3 ] = [ F P M 2 F P M F P M ]
[ R 44 R 45 0 R 54 R 55 R 56 0 R 65 R 66 ] × [ ϕ 4 ϕ 5 ϕ 6 ] = [ 0 F f 0 ]
where ϕ1 to ϕ6 are the fluxes flowing along Loops 1–6 in the MEC model.
{ R 11 = R 33 = 2 ( R P M + 2 R C S , g + R I r o n ) R 22 = 4 ( R P M + R C S , g + R g o u t ) R 12 = R 21 = R 23 = R 32 = 2 ( R P M + R C S , g )
{ R 44 = R 66 = 2 ( k α m R P M + 2 k α m R C S , g + R I r o n ) R 55 = 4 ( k α m R P M + k α m R C S , g + R g o u t ) R 45 = R 54 = R 56 = R 65 = 2 ( k α m R P M + k α m R C S , g )
where RPM, RIron, RCS,g, Rgout, and kαm denote the reluctances of the PM pole, the iron pole, the copper sheet and the inner air-gap, the outer air-gap, and the coefficient to modify RPM and RCS,g under the field current excitation alone, respectively. Assuming the permeabilities of the CS and PM are equal to μ0, the above parameters can be calculated as [17,18]
{ R P M = ln [ 1 + h P M / ( r C S o + l g i n ) ] / ( μ 0 α m τ p l P M a ) R C S = ln ( 1 + h C S / r C S i ) / ( μ 0 α m τ p l P M a ) R g i n = ln ( 1 + l g i n / r C S o ) / ( μ 0 α m τ p l P M a ) R g o u t = 2 ln ( 1 + l g o u t / r P M R c o ) / [ μ 0 α m τ p ( 2 l P M a + w P M a w s l o t ) ] R C S , g = R C S + R g i n k α m = ( 2 α m ) / α m
where rCSi, rCSo, rPMRco, αm, τp, lPMa, and μ0 is the inner and outer radii of the CS, the outer radius of the PMR core, the pole-arc coefficients of PM and iron poles, the pole pitch in the circumferential direction, the axial length of PM poles, and the vacuum permeability, respectively.
Based on the analysis of the flux paths and reluctance parameters of the coupling, the flux density distributions can be obtained by a superposition of the fields produced by the PMs and the field current. Based on the MEC model in Figure 2 and the flux density model in Figure 4, the no-load air-gap flux density of the coupling can be expressed as [17,18,20]
B P M ( θ ) = { k t B P M _ P M ( θ θ 2 ) / ( θ 2 θ 1 ) + B P M _ I f ;   θ 1 θ < θ 2 k t B P M _ P M + B P M _ I f ;           θ 2 θ θ 3 k t B P M _ P M ( θ θ 4 ) / ( θ 3 θ 4 ) + B P M _ I f ;   θ 3 < θ θ 4
B I r o n ( θ ) = { k t B I r o n _ P M ( θ θ 2 ) / ( θ 2 θ 1 ) + B P M _ I f ;   θ 1 θ < θ 2 k t B I r o n _ P M + B I r o n _ I f ;           θ 2 θ θ 3 k t B I r o n _ P M ( θ θ 4 ) / ( θ 3 θ 4 ) + B P M _ I f ;   θ 3 < θ θ 4
where kt denotes the coefficient to modify the trapezoid flux density waveforms shown in Figure 4a,b. θ1 to θ4 and kt can be obtained as follows:
{ θ 1 = θ 4 = τ p / 2 θ 2 = θ 3 = τ p α m / 2 k t = 2 α m / ( 1 + α m )
In addition, BPM_PM, BIron_PM, BPM_If, and BIron_If denote the average flux densities facing the PM and iron poles produced by the PM and the field current, respectively, which can be derived as
{ B P M _ P M = 2 ( ϕ 2 ϕ 1 ) / [ τ p α m ( r C S o + l g / 2 ) l P M ] B I r o n _ P M = 2 ϕ 1 / [ τ p α m ( r C S o + l g / 2 ) l P M ] B P M _ I f = 2 ( ϕ 5 ϕ 4 ) / [ τ p ( 2 α m ) ( r C S o + l g / 2 ) l P M ] B I r o n _ I f = 2 ϕ 4 / [ τ p α m ( r C S o + l g / 2 ) l P M ]

3.2. Eddy Current Field Calculation

Thanks to FEA, the distributions of eddy currents under different field currents can be easily calculated and shown in Figure 5 when the slip speed is 150 rpm. It can be seen that the directions of eddy currents do not change, while the direction of the field current is reversed. In addition, the eddy current densities rise with the increase in field current.
As a relative rotation between the PMR and the CR happens, the eddy current will be induced in the CS. Since the air-gap flux densities facing the PM and iron poles are different, the eddy current densities can be calculated separately, which can be respectively given as [17,18,20]
J P M ( r , θ ) = r σ ω B r _ P M ( θ )
J I r o n ( r , θ ) = r σ ω B r _ I r o n ( θ )
where σ and ω denote the conductivity of the CS and relative angular velocity, respectively. Br_PM (θ) and Br_Iron (θ) are the radial components of the resultant flux density in the inner air-gap facing the PM and iron poles, respectively. In fact, the axial components of the eddy current densities facing the PM and iron poles in the same loop are the same, i.e., the average eddy current densities of JPM and JIron are equal, which can be expressed as
J a v ( r , θ ) = ( J P M ( r , θ ) J I r o n ( r , θ ) ) / 2 = r σ ω ( B r _ P M ( θ ) B r _ I r o n ( θ ) ) / 2 .
In fact, the PM magnetic field is interacted by the magnetic field produced by the eddy current in the CS. Thus, an iterative process should be taken in the air-gap field calculation of the coupling. The air-gap flux density is first determined using the no-load MEC model. The eddy current density in the CS is then calculated by means of Faraday’s law under the different slip speeds between the two rotors. Afterward, the impact of the magnetic field yielded by the eddy currents on the original air-gap flux density is taken into account through Ampere’s law. The radial components of the resultant flux densities in the air-gap can be expressed as [17,18,20]
B r _ P M ( n ) ( θ ) = B P M ( θ ) + B C S _ P M ( n 1 ) ( θ )
B r _ I r o n ( n ) ( θ ) = B I r o n ( θ ) + B C S _ I r o n ( n 1 ) ( θ )
where the positive integer n represents the iteration number, and BCS_PM (θ) and BCS_Iron (θ) denote the flux densities produced by the eddy current in the CS facing the PM and iron poles, respectively. The eddy current field can be calculated based on Ampere’s laws expressed as [17,18,20]
C P M H d l = θ C S 1 θ C S 2 r C S i r C S o J P M ( r , θ )   r d r d θ
C I r o n H d l = θ C S 3 θ C S 4 r C S i r C S o J I r o n ( r , θ )   r d r d θ
where the right terms are the total currents enclosed in flux paths CPM and CIron. θCS1 to θCS4 are the edges of the eddy current field flux paths shown in Figure 3b.
Based on Equations (17) and (18), the eddy currents under the iron and PM poles are integrated along the different lengths of the flux paths (CPM and CIron) shown in Figure 3b. For simplifying calculation, the magnetic reluctances of the iron materials are ignored. Thus, the average lengths of the flux paths facing the PM and iron poles can be expressed as le_PM = π(τp/4) and le_Iron = 2(lgin + hCS), respectively [20]. Based on Equations (14), (17), and (18) and the selections of paths CPM and CIron in Figure 3b, the flux densities BCS_Iron(θ) and BCS_PM(θ) can be respectively obtained as [17,18,20]
B C S _ I r o n = k l e B C S _ I r o n
B C S _ P M ( θ ) = { B C S 1 = k 1 e m θ ;           θ 1 θ < θ 2 B C S 2 = k 2 e m θ B m / ( 1 + k l e ) ; θ 2 θ θ 3 B C S 3 = k 3 e m θ ;           θ 3 < θ θ 4
where
B m = ( B P M _ P M + B P M _ I f B I r o n _ P M B I r o n _ I f ) / 2 m = m P M ( 1 + k l e ) / 2 k l e = l e _ P M / l e _ I r o n m P M = μ 0 σ ω ( r C S o 3 r C S i 3 ) / ( 3 l e _ P M ) .
The following boundary conditions are taken as
B C S 2 ( θ 0 ) = 0
B C S 1 ( θ = θ 2 ) = B C S 2 ( θ = θ 2 )
B C S 2 ( θ = θ 3 ) = B C S 3 ( θ = θ 3 )
where Equation (22) is the main boundary condition of the problem referring to a particular point where the total currents enclosed in the intervals [θ1, θ0] and [θ0, θ4] are equal, and Equations (23) and (24) express the continuity of BCS(θ). θ0 can be determined by
θ 1 θ 0 r C S i r C S o J ( r , θ ) r d r d θ = θ 0 θ 4 r C S i r C S o J ( r , θ ) r d r d θ .
Thus,
θ 0 = 1 m ln { 1 + ( 1 e α m ) cosh [ ( 1 α m ) m τ p / 2 ] ( 1 e α m ) cosh ( m τ p / 2 ) + cosh ( m α m τ p / 2 ) } .
The coefficients k1, k2 and k3 in (20) can then be expressed as [17,18]
{ k 1 = B m ( e m θ 0 e α m m θ 1 ) / ( 1 + k l e ) k 2 = B m e m θ 0 / ( 1 + k l e ) k 3 = B m ( e m θ 0 e α m m θ 4 ) / ( 1 + k l e ) .
Thus, the radial components of the air-gap flux densities can be obtained by substituting Equations (19) and (20) into Equations (15) and (16), respectively.

3.3. Torque Calculation

The torque of the coupling can be finally determined by using the total ohmic losses dissipated in the CS [17,18,19,20,21]:
T = K s p c u _ c s ω = 2 K s l P M a σ ω C S | J ( r , θ ) | 2 r d r d θ = 2 K s l P M a π ( r C S o 2 r C S i 2 ) J a v 2 / σ ω
where Ks, pcu_CS, and Jav denote the 3D coefficient, the copper loss in the CS, and the average eddy current density in the same loop. Ks can be derived as [22]
K s = 1 tanh [ p l P M a / ( 2 r a v ) ] / [ p l P M a / ( 2 r a v ) ] 1 + tanh [ p l P M a / ( 2 r a v ) ] tanh [ p l o / ( 2 r a v ) ]
where rav = (rCSo + rCSi)/2, and lo and p denote the overhang-length of the CS and the pole pair numbers of the PMR, respectively.

3.4. Efficiency

The output and input powers of the coupling can be respectively expressed as
P o u t = T o u t Ω o u t
P i n = P o u t + p c u _ C S + p c u _ f w + p o t h e r
where Tout, Ωout, and pcu_fw denote the output torque, the output angular speed, and the copper loss of the field winding, respectively. pother denotes the other losses including core loss, stray loss, mechanical loss, and so on. The coupling efficiency can then be obtained as
η C = P o u t / P i n × 100 % .
Ignoring other losses, the input torque is almost equal to the output torque. Thus, the input power of the coupling can be estimated as [20]
P i n = ( 1 s ) 2 P N
where PN and s denote the rated power of pump and fan loads, and the slip between the PMR and the CR, respectively. The efficiency of the coupling can be simplified as
η C = ( 1 s ) × 100 % .
With the increase in slip, Pin declines much faster than ηc. Based on [6], the energy-saving effect of the system can be still more pronounced when compared to the method through throttling valves and baffles.

4. Results and Discussions

The major design parameters of the HE-ASECC are tabulated in Table 1. The air-gap magnetic fields, eddy currents, and torque-slip characteristics are calculated by the proposed analytical method and verified by FEA. The loss and efficiency are subsequently investigated by FEA.

4.1. Air-Gap Magnetic Field

Figure 6 shows the radial components of the no-load inner air-gap flux densities between the PMR and the CR under different field currents. When If increases from −10 to 10 A, the amplitude of BPM decreases from about 0.7 to 0.5 T (about 30%), while BIron decreases from 0.2 to −0.3 T. The variation of BPM is smaller that of BIron, which can be explained by the fact that the reluctance of PM is much larger than that of iron poles. In addition, the difference between analytical and FEA methods is found, which is caused by the flux leakage between adjacent PM and iron poles, which is difficult to be considered in the analytical modeling. The above analysis illustrates a fairly good effect of the proposed approach on air-gap magnetic field regulation.
Figure 7a,b show the average values of Br_PM and Br_Iron as functions of slip speed under different field currents. It is obvious that the average value of Br_PM declines from about 0.62 to 0.41 T, while that of Br_Iron increases from about −0.3 to 0.15 T when the slip speed is equal to 50 rpm. With the increase in slip speed, the average values of Br_PM and Br_Iron decline and rise gradually, respectively. Based on Equations (14) and (17), the eddy current field is enhanced with the increase in slip speed. Although the direction of the flux passing through iron pole is changed by the opposite field currents (e.g., If = −10 and 10 A), the direction of the eddy current field remains unchanged. The reason is that the direction of eddy current loop is mainly determined by JPM, which is much larger than JIron.

4.2. Eddy Current Density

As expected, the eddy current densities facing the PM and iron poles increase gradually with the increase in field current as shown in Figure 8. In addition, the eddy current versus rotor position distributions under the PM and iron poles are not uniform, which can be explained by Equations (12)–(14) and the air-gap flux density distributions shown in Figure 6. It should be noted that the directions of eddy currents facing the iron pole are unchanged when the field current changes from −10 to 10 A. The reason is that JPM and JIron are respectively proportional to Br_PM and Br_Iron, while Br_PM is much larger than Br_Iron, which determines the flowing directions of eddy currents.
Figure 9 shows the variations in average eddy current densities facing the PM and iron poles of the proposed coupling under different slip speeds. With the increase in field current, both JPM and JIron rise gradually, and they also increase with the increase in slip speed under a given field current. It should be noted that the average eddy current densities calculated by analytical and FEA methods agree well, which illustrates that the eddy current calculation model shown in Equation (14) is sufficiently reliable.

4.3. Torque Characteristics

Figure 10 shows the torque-slip characteristics of the proposed coupling under different field currents. In the lower slip speed range, the torque presents an ascending trend, the reason of which is that the eddy current density increases proportionally with slip speed. Due to the skin effect, the torque increases more and more slowly. Moreover, it can be seen that the torque experiences a gradually reduction with the decrease in If from 10 to −10 A at a given slip speed. This implies that a wide torque regulation range can be achieved by field control. The above analysis confirms the validity of HE control for the load speed regulation by the proposed coupling.

4.4. Losses and Efficiency

Figure 11 shows the core losses of the coupling versus slip speed under different field currents. With the increase in slip speed, the core losses of the coupling gradually increase, and are proportional to the frequencies of the alternating magnetic fluxes passing through the cores. In fact, the magnetic fluxes passing through the PMR core and the stator core are almost not alternating, which results in low core losses in the PMR and stator cores. Thus, the core loss is mainly produced in the CR core. In addition, the core loss also increases with the increase in field current, which enhances the alternating magnetic fluxes passing through the CR core.
Figure 12 shows the copper losses of the CS versus slip speed under different field currents. With the increase in slip speed and field current, the copper losses rise gradually. The copper loss of field winding is proportional to the field current, which can obtain the maximum value as 97.5 W when the magnitude of the field current is set to 10 A. It takes up a small percentage of the total copper loss of the coupling, which illustrates that the HE control is acceptable for the speed regulation application.
Figure 13 shows the efficiencies of the coupling under different field currents when the speed of the prime motor is 1500 rpm. The efficiency of the coupling reaches the highest when If is zero, which is attributed to the fact that the copper loss of the field winding is the lowest. It decreases from 96.7 to 56.9% when the slip speed increases from 50 to 600 rpm. When If decreases from 10 to 0 A, the efficiency of the coupling slightly declines, while its reduction is enlarged with the further decrease in If from 0 to −10 A. This can be attributed to the fact that Pout decreases with the decrease in the field current. In addition, pcu_fw declines when the field current decreases from 10 to 0 A, and it then rebounds to a higher value with the further decrease in the field current. In a lower slip speed range, the variations in efficiencies are obvious under different field currents. It is responsible for this that the copper loss of field winding occupies a large percentage of the relatively small output power when the torque of the coupling is small. Although the efficiency of the HE-ASECC declines with the increase in slip speed linearly, for a pump system, an energy-saving effect can still be achieved based on Equation (33) compared with the valves and baffles controls.

5. Conclusions

This paper proposes a novel brushless HE-ASECC with a simple AFES, which makes the online speed regulation reliably and effectively for the energy-saving of the aged pump and fan systems. The MEC-based analytical method is presented to calculate the magnetic field and torque of the proposed coupling efficiently. In addition, the analytical method provides an approach to solve the eddy current field problem under an asymmetric magnetic field of the PMR. Finally, the magnetic field and torque adjustment ability of the coupling illustrates that a wide speed regulation range can be achieved by the HE control. It also reveals that that the proposed speed regulation method can be applied in the field of variable speed drive for energy-saving reconstruction.

Author Contributions

Conceptualization: Y.L., H.L., and H.Y.; methodology: Y.L., Q.T., and H.L; software: Y.L. and Q.T.; validation: H.Y. and H.L.; formal analysis: Y.L. and H.Y.; investigation: Y.L., H.Y., and H.H.; resources: H.L.; data curation: H.L.; writing—original draft preparation: Y.L.; writing—review and editing: H.L., H.Y., and S.F; visualization: Y.L.; supervision: H.L. and H.H.; project administration: H.L.; funding acquisition: H.L.

Funding

This work was supported by the National Natural Science Foundation of China (51577026).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1 shows a comparison of HE-ASECC with VB, AS-PMECC, and VFD in energy-saving effect, space, efficiency, power electronic converter capacity, reliability, and cost.
Table A1. Full comparison among different energy saving devices.
Table A1. Full comparison among different energy saving devices.
ItemVBAS-PMECCVFDHE-ASECC
Energy-saving effect★★★★★★★★★★★
Space occupation★★★★★★★★★★★★★★
Efficiency★★★★★★★★★★★
Power electronic converter capacityNoneNoneLargeNone
Reliability★★★★★★★★★★★★★★
Cost★★★★★★★★★★★★★★★

References

  1. Bloxham, D.A.; Wright, M.T. Eddy-current coupling as an industrial variable-speed drive. Proc. IEE 1972, 119, 1149–1154. [Google Scholar] [CrossRef]
  2. Wallace, A.; Wohlgemuth, C.; Lamb, K. A high efficiency, alignment and vibration tolerant, coupler using high energy-product permanent magnet. In Proceedings of the Seventh International Conference on Electrical Machines and Drives, Durham, UK, 11–13 September 1995; pp. 232–236. [Google Scholar]
  3. Abdelaziz, E.A.; Saidur, R.; Mekhilef, S. A review on energy-saving strategies in industrial sector. Renew. Sustain. Energy Rev. 2011, 15, 150–168. [Google Scholar] [CrossRef]
  4. Khalid, N. Efficient energy management: Is variable frequency drives the solution. Procedia Soc. Behav. Sci. 2014, 145, 371–376. [Google Scholar] [CrossRef]
  5. Saidur, R.; Mekhilef, S.; Ali, M.B.; Safari, A.; Mohammed, H.A. Applications of variable speed drive (VSD) in electrical motors energy-savings. Renew. Sustain. Energy Rev. 2012, 16, 543–550. [Google Scholar] [CrossRef]
  6. Wallace, A.; Jouanne, A. Industrial speed control: Are PM couplings an alternative to VFDs? IEEE Trans. Ind. Appl. 2001, 7, 57–63. [Google Scholar] [CrossRef]
  7. Canova, A.; Vusini, B. Design of axial eddy-current couplers. IEEE Trans. Ind. Appl. 2003, 39, 725–733. [Google Scholar] [CrossRef]
  8. Canova, A.; Vusini, B. Mathematic modeling of rotating eddy-current couplers. IEEE Trans. Magn. 2005, 41, 24–35. [Google Scholar] [CrossRef]
  9. Wang, J.; Lin, H.; Fang, S.; Huang, Y. A general analytical model of permanent magnet eddy current couplings. IEEE Trans. Magn. 2014, 50, 1–9. [Google Scholar] [CrossRef]
  10. Lubin, T.; Rezzoug, A. Steady-state and transient performance of axial-field eddy-current coupling. IEEE Trans. Ind. Electron. 2015, 4, 2287–2296. [Google Scholar] [CrossRef]
  11. Wang, J.; Lin, H.; Fang, S. Analytical prediction of torque characteristics of eddy current couplings having a quasi-Halbach magnet structure. IEEE Trans. Magn. 2016, 52, 1–9. [Google Scholar] [CrossRef]
  12. Dai, X.; Cao, J.Y.; Long, Y.J.; Liang, Q.H.; Mo, J.Q.; Wang, S.G. Analytical modeling of an eddy-current adjustable-speed coupling system with a 3-segment Halbach magnet array. Electr. Power Comp. Syst. 2015, 43, 1891–1901. [Google Scholar] [CrossRef]
  13. Zheng, D.; Wang, D.; Li, S.; Zhang, H.; Yu, L.; Li, Z. Electromagnetic-thermal model for improved axial-flux eddy current couplings with combine rectangle-shaped magnets. IEEE Access 2018, 6, 26383–26390. [Google Scholar] [CrossRef]
  14. Yang, H.; Zhu, Z.Q.; Lin, H.; Chu, W.Q. Flux adjustable permanent magnet machines: A technology status review. Proc. CJEE 2016, 2, 14–30. [Google Scholar]
  15. Aydin, M.; Huang, S.; Lipo, T.A. Design, analysis, and control of a hybrid field-controlled axial-flux permanent-magnet motor. IEEE Trans. Ind. Electron. 2010, 57, 78–87. [Google Scholar] [CrossRef]
  16. Lubin, T.; Rezzoug, A. 3-D analytical model for axial-flux adjust speed eddy-durrent couplings and brakes under steady-state conditions. IEEE Trans. Magn. 2015, 51, 1–12. [Google Scholar] [CrossRef]
  17. Mohammadi, S.; Mirsalim, M.; Vaez-Zadeh, S. Nonlinear modeling of eddy-current couplers. IEEE Trans. Energy Convers. 2014, 29, 224–231. [Google Scholar] [CrossRef]
  18. Mohammadi, S.; Mirsalim, M. Double-sided permanent-magnet radial-flux adjust speed eddy-current couplings: Three-dimensional analytical modelling, static and transient study, and sensitivity analysis. IET Electr. Power Appl. 2013, 7, 665–679. [Google Scholar] [CrossRef]
  19. Li, Z.; Wang, D.; Zheng, D.; Yu, L. Analytical modeling and analysis of magnetic field and torque for novel axial flux eddy current couplers with PM excitation. AIP Adv. 2017. [Google Scholar] [CrossRef]
  20. Li, Y.; Lin, H.; Yang, H.; Fang, S.; Wang, H. Analytical analysis of a novel flux adjustable permanent magnet eddy-current coupling with a movable stator ring. IEEE Trans. Magn. 2018, 3, 1–4. [Google Scholar] [CrossRef]
  21. Wang, J.; Zhu, J. A simple method for performance prediction of permanent magnet eddy current couplings using a new magnetic equivalent circuit model. IEEE Trans. Ind. Electron. 2018, 65, 2487–2495. [Google Scholar] [CrossRef]
  22. Russell, R.L.; Norsworthy, K.H. Eddy currents and wall losses in screened-rotor induction motors. Proc. IEE Part A Power Eng. 1958, 105, 163–175. [Google Scholar] [CrossRef]
Figure 1. The magnetic field control principle of the coupling. (a) The topology. (b) If = −10 A. (c) If = 0 A. (d) If = 10 A.
Figure 1. The magnetic field control principle of the coupling. (a) The topology. (b) If = −10 A. (c) If = 0 A. (d) If = 10 A.
Energies 12 00308 g001
Figure 2. The magnetic circuit equivalent (MEC) model of the hybrid excited adjustable speed eddy current coupling (HE-ASECC). (a) permanent magnet (PM) excited. (b) If excited.
Figure 2. The magnetic circuit equivalent (MEC) model of the hybrid excited adjustable speed eddy current coupling (HE-ASECC). (a) permanent magnet (PM) excited. (b) If excited.
Energies 12 00308 g002
Figure 3. The 2D view of the HE-ASECC. (a) Axial section. (b) Axial cross section.
Figure 3. The 2D view of the HE-ASECC. (a) Axial section. (b) Axial cross section.
Energies 12 00308 g003
Figure 4. The model of air-gap flux density. (a) BPM excited by PMs. (b) BIron excited by PMs. (c) BPM excited by If. (d) BIron excited by If.
Figure 4. The model of air-gap flux density. (a) BPM excited by PMs. (b) BIron excited by PMs. (c) BPM excited by If. (d) BIron excited by If.
Energies 12 00308 g004
Figure 5. The distributions of eddy currents. (a) If = −10 A. (b) If = 0 A. (c) If = 10 A.
Figure 5. The distributions of eddy currents. (a) If = −10 A. (b) If = 0 A. (c) If = 10 A.
Energies 12 00308 g005
Figure 6. No-load air-gap magnetic flux density distributions.
Figure 6. No-load air-gap magnetic flux density distributions.
Energies 12 00308 g006
Figure 7. The average values of flux densities under different field currents. (a) The average values of Br_PM. (b) The average values of Br_Iron.
Figure 7. The average values of flux densities under different field currents. (a) The average values of Br_PM. (b) The average values of Br_Iron.
Energies 12 00308 g007
Figure 8. The eddy current densities versus the rotor position.
Figure 8. The eddy current densities versus the rotor position.
Energies 12 00308 g008
Figure 9. The average eddy current densities under different slip speeds.
Figure 9. The average eddy current densities under different slip speeds.
Energies 12 00308 g009
Figure 10. The torque–slip characteristics under different field currents.
Figure 10. The torque–slip characteristics under different field currents.
Energies 12 00308 g010
Figure 11. The core losses under different field currents.
Figure 11. The core losses under different field currents.
Energies 12 00308 g011
Figure 12. The copper losses of the copper sheet (CS) under different field currents.
Figure 12. The copper losses of the copper sheet (CS) under different field currents.
Energies 12 00308 g012
Figure 13. The efficiencies of the coupling under different field currents.
Figure 13. The efficiencies of the coupling under different field currents.
Energies 12 00308 g013
Table 1. Major design parameters of the coupling.
Table 1. Major design parameters of the coupling.
ItemSymbolValueUnit
Inner radical of CSrCSi46mm
CS thicknesshCS3mm
CS conductivityσ5.77 × 107S/m
Over hang lengthlo5mm
CR yoke thicknesshCRc10mm
Air-gap lengthlgout/lgin1mm
PM inner radicalrmi50mm
Magnetic pole heighthPM6mm
PM axial lengthlPMa30mm
PM coercive forceHc−890kA/m
Pole-arc coefficientαm0.8
Pole pair numbers of the PMRp16
Axial distance between the two annular iron coreswPMa10mm
PMR yoke thicknesshPMRc10mm
Stator yoke thicknesshsc10mm
Stator slot widthwslot40mm
Stator slot heighthslot10mm
Field coil turns numberNf300turns

Share and Cite

MDPI and ACS Style

Li, Y.; Lin, H.; Huang, H.; Yang, H.; Tao, Q.; Fang, S. Analytical Analysis of a Novel Brushless Hybrid Excited Adjustable Speed Eddy Current Coupling. Energies 2019, 12, 308. https://doi.org/10.3390/en12020308

AMA Style

Li Y, Lin H, Huang H, Yang H, Tao Q, Fang S. Analytical Analysis of a Novel Brushless Hybrid Excited Adjustable Speed Eddy Current Coupling. Energies. 2019; 12(2):308. https://doi.org/10.3390/en12020308

Chicago/Turabian Style

Li, Yibo, Heyun Lin, Hai Huang, Hui Yang, Qiancheng Tao, and Shuhua Fang. 2019. "Analytical Analysis of a Novel Brushless Hybrid Excited Adjustable Speed Eddy Current Coupling" Energies 12, no. 2: 308. https://doi.org/10.3390/en12020308

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