1. Introduction
Magnetic cores of electrical rotating motors are constructed using non-oriented steel sheets which should have isotropic properties; however, these sheets have certain anisotropy [
1,
2]. The anisotropy of the magnetic cores causes fluctuations of the magnetic rotating field in the rotating motors. This phenomenon can be especially noticed in two-pole induction motors, which are distinguished by a short path of the magnetic flux in the air gap. The saturation level of the teeth and yokes of the magnetic core has a significant influence on the amplitudes of these fluctuations [
3]. The fluctuations of the rotating field cause disturbances in phase currents and voltages induced in stator windings, which result in parasitic components in the electromagnetic torque. Currently, these phenomena are considered in the design of electrical drive control systems [
4,
5]; model studies using Finite Elements Analysis techniques can predict and explain these phenomena [
6,
7,
8,
9]. The fluctuation of the co-energy of the rotating magnetic field can be treated as a good indicator for the anisotropy of both the stator and the rotor. This conclusion was confirmed by field calculations performed for a medium-sized squirrel cage induction motor. In order to do that, a circuit model of this motor with an anisotropic core was prepared, assuming that the magnetic flux density along the air gap has a sinusoidal distribution. The circuit model is based on properties of the co-energy function of the magnetic field excited by winding currents. This function has a simple form when the equivalent magnetising current is used and magnetic leakage fields are described separately. This enabled mathematical formulas to be simplified in cases where nonlinear magnetisation characteristics should be considered. The unambiguous magnetisation characteristics of the non-oriented sheets were determined on the basis of magnetic hysteresis loops measured for the direction of both the easy and hard magnetisation [
10]. In the above-mentioned model, the common features of relations concerning flux linkages and currents are a consequence of using the prime co-energy function of the main field as the function of the magnetising current. These relations depend on the magnetomotive (MMF) axis position. They can model a non-zero angle between the MMF axis and the flux density field axis. In this approach, the model of the induction motor with the anisotropic core is similar to the model of a two-pole motor with the salient poles of the rotor or the stator [
11]. It was applied in the formulation of the proposed circuital model of the induction motor.
The anisotropic properties of the rotor core also affected the leakage fluxes of the stator windings, which refer to parts of leakage fluxes that close through the air gap and rotor teeth. This effect may be evident in starting states, especially when induction motors are supplied by voltage source inverters [
5]. The study of these effects may be useful in the diagnostics of motor core quality [
12,
13] and in the assessments of the internal motor imbalance [
14,
15].
2. Flux Linkages—Currents Relations in Motors with an Anisotropic Rotor
The co-energy of the magnetic field in the core and in the air gap of the two-pole induction motor with an anisotropic rotor core is a function of
N independent currents
in of the windings and the angular position
ϕ of the rotor. In the circuit model of the rotating motors the main circuit and leakage magnetic circuits are described separately. This is also acceptable when the nonlinear magnetising characteristic of the core is taken into account due to the specified magnetic circuits being saturated under different conditions. Therefore, the global field co-energy can be expressed as a sum of the main field co-energy and the co-energies of the leakage fields:
The main co-energy function can be written as a function of only three equivalent variables: the magnetizing current
iµ, which corresponds to the amplitude of the fundamental harmonic of the total MMF; the angular position
αµ of the total MMF axis in the air-gap with respect to the axis of the stator easy magnetization; and the angle of the rotor position
ϕ (
Figure 1 and
Figure 2). The co-energies of the leakage fields depend on the individual winding currents and the angle of the rotor position. The space vector of the rotating magnetizing current is specified in a complex domain by the following formula:
where
N denotes the number of independent stator and rotor currents, the angle
αn is an angular position of the
nth winding magnetic axis versus the stator axis or the rolling direction (RD) rotor axis, and
denotes the winding coefficient for the fundamental spatial harmonic of the flux density distribution in the air gap. Therefore, the absolute value and the argument of the magnetizing current can be written as follows:
An approximation using the complex Fourier series is proposed to write separately the independent angular variables:
The double series (4) express the main co-energy function in dependence on the position of the magnetising current vector versus the stator and rotor RD axis. For a three-phase squirrel cage induction motor with a uniform air gap and a linear isotropic magnetic core:
k =
l = 0. For a nonlinear and anisotropic core where
k, l = {0, ±2, ±4, …}, the co-energies of the stator leakage fields can be approximated using a similar complex Fourier series:
The independent winding flux linkages are defined as the partial derivatives of the co-energy function with respect to the independent currents:
The main flux linkage of each separate winding can be written as follows:
Hence, the result of the differentiation operations above can be written as follows:
Subsequently, the main flux linkage of each separate winding can be expressed using two components: component ‘
x’ in the direction of the magnetizing current spatial vector and component ‘
y’ in the perpendicular axis (
Figure 2). The amplitudes of these components are expressed by the following series:
The functional coefficients
and
of the series are the equivalent internal magnetization characteristics in the perpendicular
x and
y axes, respectively. They are associated with the elements of the co-energy Fourier series by the following formulas:
Based on Equations (7)–(9), the series that approximates the main flux linkage of each
nth independent circuit can be written in the following complex form:
where
Owing to the nonlinearity of the main magnetic circuit, the third harmonic of the flux density distribution along the air gap generates a significant component of the main flux linkages. Hence, the formula approximating the main flux linkages of the stator winding should be completed by the following component:
where the constant
denotes the winding coefficient of the third harmonic.
The leakage fluxes of the stator windings can be approximated by Series (15), considering both the core anisotropy in the linear magnetic condition and the effect of the rotor slots:
where
denotes the coefficients of the series approximating the leakage inductance versus rotor position.
Formulas (12), (14), and (15) enable the frequency components of the voltage Fourier spectra caused by the motor core anisotropy to be determined.
3. Predicted Spectra of Induced Stator Phase Voltages
The magnetising current defined by formula (3) has a constant value when the motor operates in a steady-state; however, the angular positions of both the magnetising current space vector (2) and the rotor change linearly with time:
where
is the fundamental pulsation and
s denotes the rotor slip.
Substituting variables
in Equation (13) and in Series (15) by the assumed time functions (16), the pulsation of the predicted components of the main flux linkages can be specified as follows:
where
ρ denotes multiples of the fundamental pulsation.
The fundamental component was expressed for k = l = 0 resulting in ρ = ±1. The higher significant components were obtained for abs(k) + abs(l) ≤ 4. When only the stator core was magnetically anisotropic, the flux linkage pulsation components were indicated for k ≠ 0 and l = 0, resulting in ρ = (±1, ±3, ±5). Meanwhile, when the rotor core was magnetically anisotropic, the flux linkage components were indicated for k = 0 and l ≠ 0, resulting in ρ = ±(1 ± 2s), ±(1 ± 4s),…. When both the stator and the rotor core were magnetically anisotropic, the additional components of the series shown in (12) and (15) obtained for k = l = ±2 and k = −l should be considered; therefore, ρ = ±(1 ± 2s), ±(3 ± 2s),….
The main phase voltages induced in the three-phase stator windings are time derivatives of the corresponding phase flux linkages according to formula:
The Fourier series components of the induced main voltage are the same as those for the main flux linkages. Relative amplitudes of the induced voltage components are equal to those of the flux linkage components multiplied by the corresponding pulsation. The specified series components enable the separation of the positive-, negative-, or zero-sequence symmetrical components of the phase voltage. They can be determined using Formula (19):
Pulsation components of the phase voltage that fulfil the condition sign(k – m + ls) ≠ sign(m) form the positive-sequence components; otherwise, negative components are formed. The components of the zero-sequence of the phase voltages that were induced by the main flux linkages were not generated by the magnetic core anisotropy; however, these components were induced by the nonlinearity of the core magnetization characteristic.
The stator leakage fields induced additional phase voltage components that can be approximated by the following series:
The presented analysis enables the predicted voltage frequencies to be listed, which can be generated in a symmetrical induction motor with a magnetically anisotropic core supplied by symmetrical mono-harmonic current sources. The final results are presented in
Table 1.
These predicted frequencies are caused by
the rotating field in the nonlinear anisotropic core,
the leakage stator fields in the anisotropic core,
the rotating field in the nonlinear isotropic core.
4. Field Calculations of Effects Generated by an Anisotropic Motor Core
A qualitative and quantitative identification of the magnetic field co-energy, electromagnetic torque, and flux linkages of the stator winding was performed for a 5.5 kW, 400 V, 50 Hz two-pole induction motor (six slots per pole and phase), omitting rotor cage bars.
Figure 3 presents a cross-section of the motor core. The magnetic field distribution was calculated assuming that the rotor core was anisotropic and the magnetic permeability in the RD axis was twice as high as in the TD axis, while the stator core and the shaft were isotropic. The axis of the MMF produced by stator currents was parallel to the
y axis.
The unambiguous magnetisation characteristics of the core material were determined for the non-oriented sheet M530-50A (
Figure 4).
Calculations were performed using MagNet 2D software (Infolytica Corp. Montreal, Canada). Three-phase stator winding was supplied using symmetrical three-phase current sources. The amplitudes of these currents were increased until the flux density in the stator teeth was equal to 1.8 T; the unwound rotor was stationary. The field equation was solved in 200 points in one current period, i.e., in one rotation of the MMF axis in the air gap. The waveforms of the magnetic field co-energy and energy (
Figure 5) as well as the flux linkages of the stator phase windings were treated as input values in further analyses. In the next step, the Fourier spectra of these waveforms were calculated. It is noteworthy that for isotropic cores with a sinusoidal distribution of the MMF, the co-energy of the magnetic main field exhibited a constant value. The anisotropic properties of the stator core resulted in even harmonics in the co-energy waveform, which yielded fluctuations in both the electromagnetic torque and flux linkages of the windings.
Similar results were obtained for the anisotropic stator and isotropic rotor. Analogous components occurred in numerical calculations when the same core had two-phase sinusoidal distributed windings. Detailed calculations were performed for the anisotropic rotor core, while the stator core was isotropic. The harmonics of the electromagnetic torque at different values of the magnetising current for the anisotropic rotor core were calculated using the virtual work principle.
The results obtained for relative values of 0.5, 1.0, and 1.5 of the magnetising current (with respect to the rated value 2.7 A) are presented in
Table 2.
As shown, the amplitudes of the torque harmonics did not exceed 0.1 Nm at the rated magnetizing current, i.e. 1% of the rated torque. It is noteworthy that increasing the magnetizing current from 100% to 150% yielded approximately a 4X increase in the amplitude.
The waveforms of the phase flux linkages, obtained during one turn of MMF axis were examined based on their Fourier spectra. The harmonic components of the flux linkages were used to calculate the symmetrical components individually for each harmonic. Hence, the relative amplitudes of harmonics of the induced phase voltages were calculated using the following formula:
where
U1 denotes the fundamental component amplitude.
The harmonics of the symmetrical components of the phase voltages are presented in
Table 3.
The estimated values of the voltage components were similar to the level of the field calculation accuracy. However, the results were consistent with both the components caused by the core anisotropy and those predicted in
Table 1. The negative-sequence component of the third harmonic reached the highest value, i.e., 58 dB at the rated magnetizing current. An increase in the magnetizing current from 100% to 150% yielded an absolute amplitude that was 5x higher (72 dB) than that obtained previously.
The zero-sequence component of the fundamental harmonic (40 dB), as the effect of the anisotropy, was dominated by the zero-component of the third harmonic (52 dB), which was caused by the nonlinearity of the core magnetization characteristic. The results presented in
Table 3 indicate the harmonic levels predicted in
Table 1 for the motor with an unwound stationary rotor (slip 1). Therefore, the harmonics were generated by the main magnetic field rotating in the motor, and the level of analogical harmonics can be deduced for the motor with a rotor slip of approximately zero.
Additionally, field calculations were performed for cases where the leakage fields dominated over the main field. This refers to the motor operation when the wound rotor rotated with a slip of approximately one. In this case, the model of the motor core was modified by eliminating the main field from the inside the rotor. Furthermore, the leakage fields of the stator can be closed by the anisotropic rotor teeth. Tested effects were obtained when the rotor rotated slowly with respect to the pulsating leakage fields of the stator.
Figure 6 shows the spectrum of the zero-sequence components of the phase voltages that were calculated for stator currents with frequency 5 Hz and for the current amplitude that was 3x higher than the rated current; the motor slip was 0.9.
The slot harmonic (27–28
s)
f0 dominated in the spectrum; however, the harmonic (3–4
s)
f0, that was predicted in
Table 1 for the anisotropic rotor was evident as well. For the isotropic rotor teeth, the harmonic (3–4
s)
f0 did not occur.
It can be expected that harmonics analogous to that indicated in the voltage positive- and negative-sequence component spectra caused by the anisotropy of the core can occur in the stator phase currents when the induction motor is supplied by symmetrical sinusoidal voltage sources. However, the zero-sequence components of the phase currents can occur only for delta-connected windings.
5. Experimental Results
Laboratory experiments were performed to validate conclusions from the analysis of both series that approximated functions of the flux linkages of stator windings and the field calculations. The tested motor, which was supplied directly from a 3 × 400 V grid, operated in idle running at a slip of 0.001. The voltages and phase currents of the stator were saved at a sampling frequency of 10 kS/s; meanwhile, the recording times provided a resolution of at least 0.1 Hz [
16]. The magnetic core had a relatively low saturation when the motor was supplied with the rated voltage.
Figure 7,
Figure 8,
Figure 9 and
Figure 10 show results of the motor measurements during idle running. The relative values of the higher harmonics of the currents caused by the anisotropy of the rotor core did not exceed 1% of the fundamental harmonic value. The zero-component of the phase voltages at 50 Hz, which indicates the anisotropy of the stator core, was approximately 43 dB (
Figure 7). Additionally, the zero-component of the phase voltages at a frequency of (3–2
s)50 Hz = 149.9 Hz and 46 dB (
Figure 8) and the positive-sequence components of the phase currents at a frequency of (1 ± 2
s)50 Hz = {49.9, 50.1} Hz and above 60 dB (
Figure 9) indicated that the rotor core was anisotropic.
Analogical harmonics at frequency (1 ± 2
s) 50 Hz = {49.9, 50.1}Hz and values of approximately 50 dB occurred in the spectrum of the negative-sequence component of the phase currents (
Figure 10). It is noteworthy that they were predicted as shown in
Table 1 for both the stator and rotor core anisotropy.
The results of the measurements allowed us to state that the spectrum of the positive-sequence component of voltages was contaminated, whereas the spectrum of the negative-sequence component had only a fundamental component bigger than 40 dB. It results in the spectrum of the negative-sequence component of phase currents being a more appropriate signal to assess the rotor anisotropy.
6. Conclusions
The anisotropy of the nonlinear magnetization characteristics of the core of the three-phase induction motors contributed to the oscillatory components of both the co-energy and energy of the rotating magnetic field. Moreover, this anisotropy caused a certain asymmetry of the magnetic leakage fluxes, which resulted in oscillations in the electromagnetic torque. Additionally, modulations of the perimeter of the stator phase current when the motor was supplied by voltage sources or those of the perimeter of the phase voltages when the motor was supplied by current sources were observed. The frequencies of these oscillations in steady states can be determined based on the presented Fourier series that approximate the functions of the winding flux linkages with respect to the position of the rotating field axis.
Similarly, the series that approximates independent phase currents relative to the phase flux linkages can be written based on the functions of the main field energy. These series were formulated with the assumption that the anisotropy can be modelled using the internal magnetization characteristics in both the magnetizing current vector and perpendicular directions. This enabled the dependence between the components of the amplitude spectrum of the voltages and currents and the rotor speed to be determined. The predicted oscillations were separable using field modelling despite their relative values of <1%. A high accuracy of the field calculations was achieved when the windings were supplied by current sources. In this case, only the field equation for the steady states at the rotating spatial vector of the magnetizing current was solved. The amplitude spectra of the induced phase voltages corresponded to the amplitude spectra of the flux linkages of the windings. It is noteworthy that the presentation of the results in the form of the spectra of symmetrical voltage components (field calculations) or in the form of the spectra of symmetrical phase current components (based on measurements) was the most representative for diagnostics purposes. This enabled an analysis of experimental results that were dependent on supply voltage deformations. The frequencies and levels of signals caused by the core anisotropy can be useful in diagnostics pertaining to the asymmetry of motor windings as it generates analogical harmonics of the motor currents.