1 Introduction

Magnetic flux density harmonics in the airgap of electric machines create distorted voltage waveforms and acts as exciting forces that can cause unwanted vibrations. These problems are difficult to mitigate when a machine is already in operation and therefore a lot of effort is made during the design phase to eliminate them. Still, many in-operation machines experience problems related to harmonics and often the solution is to mechanically reinforce and change modal shapes which is expensive and inconvenient. Vibrations and noise created by the magnetic flux in the mechanical parts are defined as magnetic. The magnetic forces have two origins, one is due to magnetostriction and the other is explained by Maxwell’s stress tensor. It shows that the force exerted on a material is proportional to the magnetic flux density squared  [1,2,3,4]. The magnetic flux density in the airgap of an electric machine can be described as

$$\begin{aligned} B_\delta (\theta ,t) = \varLambda (\theta ,t)F(\theta ,t), \end{aligned}$$
(1)

where \(B_\delta \) is the airgap magnetic flux density, \(\varLambda \) is the permeance per unit area and F is the magnetomotive force  [5]. From (1) and Maxwell’s stress tensor, the forces acting on the structure of the electric machine can be found. The main parameters affecting the magnetic field in the airgap are listed by Traxler-Samek  [6] among others.

Analytical and experimental work on reducing noise and vibration for grid connected as well as frequency converter connected induction machines by injecting proper current harmonics into the stator has been done previously  [7,8,9,10,11,12]. A lot of research is currently being done on noise and vibration reduction in permanent magnet synchronous machines  [13, 14]. It has been shown that forces acting on the stator due to unbalanced magnetic pull can be compensated for by splitting the rotor winding and actively controlling the current to produce a resulting force vector  [15, 16]. Injection of current harmonics into the field winding of synchronous machines can be done to create an electromechanical filter, a machine that is connected in parallel to a non linear load and compensate for the harmonics produced by the load  [17,18,19]. Reduction of harmonics in the phase voltage waveform of a single-phase synchronous machine by injection of current harmonics into the field winding has been shown in  [20].

In this contribution a method to change the flux density harmonics in the airgap by injection of current harmonics in the field winding is presented. The focus is on affecting the phase voltage and the exciting magnetic forces. The appropriate currents are calculated through a minimization routine utilizing an in-house FEM model of an experimental generator. The results are then validated experimentally on a 12-pole synchronous machine in the lab. A simulation of a 33 MVA full-scale generator with pronounced Walker noise  [21] is done, to see what type of currents a full scale generator would require.

2 Theory

The idea is based on modulating the rotor field current with one or several AC components that targets specific magnetic flux density harmonics to either eliminate certain voltage harmonics or exciting force harmonics. The magnetic flux density in the airgap of an electric machine can be expressed as

$$\begin{aligned} B_{s} (\theta ,t) = \bigg [\sum _{n=1}^\infty A_n \sin \left[ p\big (n\theta + \omega _\mathrm{m} t\big ) \right] \bigg ]B_\mathrm{r}(t), \end{aligned}$$
(2)

where \(A_n\) is a scaling factor for the different space harmonics, \(\omega _\mathrm{m}\) is the mechanical angular frequency, p is the number of pole pairs in the machine, and \(B_\mathrm{r}\) is the amplitude of the radial magnetic flux density in the airgap produced by the rotor. The field current, \(I_\mathrm{f}\), in the machine is changed by adding harmonics

$$\begin{aligned} I_\mathrm{f}(t) = I_{\mathrm{f},0} + \sum _{k=1} \hat{I}_{\mathrm{f},k} \sin (\omega _{k} t+ \alpha _k), \end{aligned}$$
(3)

where \(\omega _k\) and \(\alpha _k\) are the angular frequency and the phase of the added current harmonic. The contribution from the rotor to the magnetic flux density in the airgap depend on the field current according to

$$\begin{aligned} B_\mathrm{r}(t) = B_{\mathrm{r},0}(I_{\mathrm{f},0}) + \sum _{k=1} \hat{B}_{\mathrm{r},k}(I_{\mathrm{f},k}) \sin (\omega _{k} t+ \alpha _k). \end{aligned}$$
(4)

By inserting (4) into (2) the complete expression for the magnetic field in the airgap of the machine when the rotor field current is modulated can be obtained as

$$\begin{aligned} B_{s} (\theta ,t)&= \bigg [\sum _{n=1}^\infty A_n \sin \left[ p\big (n\theta + \omega _\mathrm{m} t\big ) \right] \bigg ]\bigg [B_{\mathrm{r},0}(I_{\mathrm{f},0}) \nonumber \\&\qquad {}+ \sum _{k=1} \hat{B}_{\mathrm{r},k}(I_{\mathrm{f},k}) \sin (\omega _{k} t+ \alpha _k)\bigg ]. \end{aligned}$$
(5)

In order to show the effect on stator output voltage (dependent on B directly) and exciting magnetic force (proportional to \(B^2\)), use \(\omega _k = kp\omega _\mathrm{m}\), and consider adding one current harmonic k to the DC-rotor current. Only the fundamental of the zeroth order field current, \(A_n = 1, n = 1\) was considered for clarity. The magnetic flux density for this simplified case can then, from (5), be expressed as

$$\begin{aligned} B_s(\theta ,t)&= B_{\mathrm{r},0}(I_{\mathrm{f},0})\sin (p\omega _\mathrm{m} t) +\frac{B_k(\hat{I}_{\mathrm{f},k})}{2} \Big (\cos \big [(k+1)\nonumber \\&\qquad {} p\omega _\mathrm{m} t - p\theta + \alpha _k\big ] + \cos \big [(k-1)p\omega _\mathrm{m} t + p\theta \nonumber \\&\qquad {} + \alpha _k\big ]\Big ). \end{aligned}$$
(6)

The stator output voltage of the machine can be obtained by integrating \(\theta \) for one coil span and the number of poles. This means that a “\(k+1\)”, or “\(k-1\)”, component added in the field current will affect the voltage component “k”. The method can easily be generalized to include higher harmonics.

Equation (6) was squared to get an expression that is directly proportional to the exciting force, since the Maxwell stress tensor depends partly on the square of the flux density. One point per pole on the stator side was of interest, \(\theta = 0^\circ \) was chosen for simplicity. The amplitude and order of the \(B^2\) harmonics can be extracted from

$$\begin{aligned} B_s^2(t)&= \Bigg [B_{\mathrm{r},0}(I_{\mathrm{f},0})\cos (p\omega _\mathrm{m} t) + \frac{B_k(\hat{I}_{\mathrm{f},k})}{2} \bigg (\cos \big [(k+1) \nonumber \\&\qquad {} p\omega _\mathrm{m} t)\big ] + \cos \big [(k-1)p\omega _\mathrm{m} t)\big ]\bigg )\Bigg ]^2 \nonumber \\ {}&= \frac{B_{0}^2(I_{\mathrm{f},0})}{2}\bigg [1 + \cos (2p\omega _\mathrm{m}t)\bigg ] + \frac{B_k(I_{\mathrm{f},k}) B_{0}(I_{\mathrm{f},0})}{2} \nonumber \\&\qquad {} \bigg [\cos \big [(k+2)p\omega _\mathrm{m}t\big ] + \cos \big [(k-2)p\omega _\mathrm{m}t\big ] + 2\cos \nonumber \\&\qquad {}\big [kp\omega _\mathrm{m}t\big ]\bigg ]+\frac{B_k^2(I_{\mathrm{f},k})}{8}\bigg [2+\cos \big [2(k+1)p\omega _\mathrm{m} t\big ] \nonumber \\&\qquad {}+ \cos \big [2(k-1)p\omega _\mathrm{m} t\big ] + 2\cos (2kp\omega _\mathrm{m} t) \nonumber \\&\qquad {}+ 2\cos (2p\omega _\mathrm{m} t)\bigg ], \end{aligned}$$
(7)

where \(\alpha _k\) is set to \(0^\circ \). The expression (7) contains only harmonic waves of double order (e.g. for a 50Hz system the fundamental is 100 Hz). To affect a harmonic of order k, whose amplitude is \(B_k(I_{\mathrm{f},k}) B_{0}(I_{\mathrm{f},0})\), a current harmonic of the same order, k, should be added, and the phase adjusted accordingly. However, there is a drawback as adding a k-th component will also affect the amplitude of harmonics \(k\pm 2\).

3 Method

The method was first tested and verified in a finite element code and then comparative experiments were carried out on a smaller synchronous machine in the lab. A simulation was then made of a 33 MVA generator with known vibration problems, the specifications for each machine is shown in Table 1. The machine in the lab had a solid rotor rim where eddy currents develop if an AC-component is added to the field current, this limited the airgap magnetic flux density at higher frequencies.

Table 1 Machine data for the experimental generator and a full-scale generator with noise problems

Measurements were made at no-load operation at a lower speed in order to demonstrate the method. The frequency response from the applied field current to airgap magnetic flux was done to characterize the machine’s response to the injected field current. Another aspect of the added harmonics is worth to note. Since adding the harmonic affects the flux density, it will automatically affect everything that depends on the flux density. So if a minimization with regard to a special harmonic is done side-bands will develop, see (6) (7), that can increase instead. It is possible to add additional harmonics to reduce the effect of the developed side bands. As an example, assume that we want to affect the 5th harmonic in the phase voltage by adding a 4th and 6th harmonic in the field current. The 4th harmonic will affect the 3rd and the 5th harmonic in the phase voltage and the 4th harmonic in \(B^2\). The 6th harmonic will affect the 5th and the 7th harmonic in the phase voltage and the 6th harmonic in \(B^2\).

3.1 Finite element simulation

The method was implemented in an in-house finite element code. The code is adapted for synchronous machines and have been extensively verified against different machines, and operating conditions. The experimental and the full-scale generator were modelled using 5776 and 107,170 mesh nodes respectively, see Fig. 1. The time step was commensurate with the mesh in order to avoid numerical noise, and were 6.945\(\times 10^{-5}\) s, and 6.287\(\times 10^{-5}\) s for the two machines, respectively. For the experimental generator only one pole was modelled due to the integer slots per pole per phase compared to the full-scale generator where a large part of the geometry needed to be modelled due to fractional slots per pole per phase. Second order basis functions were used in order to better resolve higher order harmonics.

Fig. 1
figure 1

Simulated geometry with flux lines and mesh for the 200 kVA experimental generator

In order to find amplitudes and phase angles of the added harmonics a minimization procedure was used. The procedure sets the amplitude and phase of each added harmonic in the field current, then a time stepped simulation was performed. Thereafter the quantity that should be minimized was Fourier transformed and the harmonics were extracted. The minimization returned with updated values for the amplitude and phase angle and the procedure iterated until the desired harmonic component was below the minimization target. The minimization procedure used a Newton method and calculated the Jacobian numerically and made step-wise improvements based on the direction of the minima. Several numerical implementations and standard Fortran implementations were tested, and they all converged to the same result and eventually the most effective one was used.

3.2 Experimental setup

A test-rig, shown in Fig. 2, with a frequency converter controlled induction machine as the prime mover of the shaft was used. The test-rig had a 12 pole synchronous generator mounted on the shaft, and the pole shoes were laminated and the rotor rim was solid [22]. The machine was constructed to have a large airgap between stator and rotor resulting in low nominal airgap flux density compared to conventional machines. The airgap magnetic flux was measured by an in-house designed flux density sensor mounted on a stator tooth  [23]. The angular position of the shaft is required to be able to generate the appropriate current waveform, for this purpose an optical sensor was used that produced a pulse once every revolution. The angular position was approximated from this pulse and the rotational speed of the machine.

Fig. 2
figure 2

Experimental test rig used for the experiments

The field winding of the test generator was connected through slip rings to two power electronic legs on a Powerex PP200T120-ND IGBT stack with DC-link voltage \(V_\mathrm{{DC}}=60\) V and switching frequency \(f_s=10~\mathrm{kHz}\). The shaft angular position sensor and current measurement was used as input to a National Instruments sbRIO-9606 based field current controller. It controlled the field current utilizing a fixed frequency bang-bang control scheme; the polarity of the field winding voltage was set each iteration depending on if the measured current was above or below a set reference. The set-point value was based on the shaft angular position measurement and the targeted harmonics. Current was measured with a LEM LA55-p hall effect based current sensor sampled at 40 kHz with a negligible amount of noise. In Fig. 3 a simplified block diagram of the system setup is shown.

Fig. 3
figure 3

Simplified block diagram showing the field winding excitation circuit. The connection between the H-bridge and the field winding is done through slip rings

3.3 Full-scale generator, 33 MVA, with noise problems due to electromagnetic Walker harmonics

The method was tested by simulating a full-scale generator that emitted audible noise at \(f_w=700~\mathrm{Hz}\), it had been confirmed that the noise was due to Walker noise [24] that excited a mechanical resonance. The Walker noise arises due to interaction of permeance harmonics and rotor mmf. The Walker noise frequency is

$$\begin{aligned} f_w=2f_{el}nint(qm), \end{aligned}$$
(8)

where \(f_{el}\) is the stator electric frequency of the machine, q the number of slots per pole and phase, m the number of phases in the machine, and nint() denotes the nearest integer. By injecting a 14th order field current harmonic the noise can potentially be eliminated by reducing the exciting force at the resonance frequency.

4 Results

4.1 Airgap magnetic flux density waveform at constant field current

In Fig. 4 one period of the airgap magnetic flux density at a point with rated constant field current is shown. It is used as reference to compare both simulation and experimental results. It is seen that the current ripple is quite small. This is due to the large time constant of the field winding.

Fig. 4
figure 4

Measured current and magnetic flux density waveforms in the 200 kVA generator for one period measured at one point in the airgap with constant \(I_\mathrm{f} = 12.48~A\), with \(w_\mathrm{m}=2\pi f_\mathrm{m} = 2\pi \frac{1}{3}\)

4.2 Magnitude and phase response from applied field current to magnetic flux density

The FEM simulations were done in 2D and do not take into account the 3D effects that occur, e.g. eddy currents. To get those characteristics a magnitude and phase response from applied field current (\(I_\mathrm{f}\)) to magnetic flux density (B) on one point at the stator was done, see Fig. 5. A considerable reduction in measured magnetic flux density is observed as the frequency of the applied field current is increased. This can largely be attributed to the solid rotor rim where eddy currents develop.

Fig. 5
figure 5

Magnitude and phase response in the 200 kVA generator from applied field current (\(I_\mathrm{f}\)) to magnetic flux density (B). The data is normalized with the simulated flux density for easy comparison

4.3 Harmonic minimization in the 200 kVA experimental generator

In Table 2 the parameters for the added field current harmonics as described in (3) is listed, for all cases \(I_{\mathrm{f},DC} = 12.48~A\). The harmonic field current amplitudes injected in the experiments were higher than the ones from the minimization to compensate for the eddy current effects in the solid rotor rim.

Table 2 Amplitude and phase of the injected field current for different minimization targets applied to the experimental generator

4.3.1 Minimization of harmonics in the output phase voltage

Minimization of the 3rd harmonic in one of the phase voltages is done by adding 2nd and 4th order harmonics to the applied field current. See Fig. 6 for the actual measured waveforms over one period of magnetic flux density.

The magnitude of the harmonics in output phase voltage and the magnetic flux density squared at one point on a stator tooth is presented in Fig. 7. The values are normalized to the fundamental component of each quantity. Minimizing the 3rd harmonic in one phase voltage will also affect the surrounding harmonics, the sidebands. Thus, it will also affect the exciting forces. Experimental verification was done for elimination of the 3rd voltage harmonic, but in principle an arbitrary amount of harmonics can be minimized. The minimization of harmonic 3 and 5 was also simulated, the result is shown in Fig. 8. Adding a relatively high 6th order harmonic affect the 7th harmonic in the phase voltage. Due to the minimization in one phase voltage the other two phase voltages will be affected.

In Fig. 9 the stator voltages as a function of time obtained from simulations and experiments before and after minimization of the 3rd harmonic in \(U_a\) are shown. In Fig. 10 the Fourier transform of the measured voltages in the experiment is presented. An increase in the other two phase voltages is observed. In Table 3 the total harmonic distortion of the phase voltages for constant field current compared to minimization of 3rd harmonic in \(U_a\) is presented.

4.3.2 Minimization of harmonics in \(B^2\)

Minimization of the 4th harmonic in the magnetic flux density squared is done by adding a 4th order harmonic to the applied field current. See Fig. 11 for the actual measured waveforms over one period of magnetic flux density. The magnitude of the harmonics in output phase voltage and the magnetic flux density squared at one point on a stator tooth is presented in Fig. 12. The values are normalized to the fundamental component of each quantity. Experimental verification was done for elimination of the 4th harmonic but, as mentioned previously, in principle an arbitrary amount of harmonics can be minimized. In Fig. 13 the simulated result for minimization of both 4th and 8th harmonic of the magnetic flux density squared is shown. It can be observed that adding additional harmonics in \(B^2\) reduces the overall harmonic content of the exciting forces.

Fig. 6
figure 6

Measured current and magnetic flux density waveforms in the 200 kVA generator for one period measured at one point in the airgap with \(I_\mathrm{f}\) modulated to eliminate the 3rd harmonic in the phase voltage, with \(w_\mathrm{m}=2\pi \frac{1}{3}\)

Fig. 7
figure 7

Magnitude of harmonics at constant and modulated field current in the 200 kVA generator from simulation and experiment, the 3rd harmonic in the phase voltage is minimized. Top part: Output phase voltage. Bottom part: \(B^2\) at one point on the stator. Harmonic order 1 equals 2 Hz

Fig. 8
figure 8

Magnitude of harmonics at constant and modulated field current in the 200 kVA generator from simulation, the 3rd and 5th harmonic in the phase voltage is minimized. Top part: Output phase voltage. Bottom part: \(B^2\) at one point on the stator

Fig. 9
figure 9

Stator output phase voltages from experiments and simulation. Top part: Phase voltages with constant \(I_\mathrm{f} = 12.48~A\). Bottom part: Phase voltages with minimization of 3rd harmonic in \(U_a\). The high frequency noise on the waveforms from the experiments is measurement noise

Fig. 10
figure 10

Magnitude of harmonics in the measured phase voltages in the 200 kVA generator with constant \(I_\mathrm{f} = 12.48~A\), with minimization of 3rd harmonic in \(U_a\)

4.3.3 Minimization of electromagnetic walker harmonics in the 33 MVA full-scale generator

The minimization routine was done in FEM for the electric machine and the resulting field current is presented in Table 4, where it can be seen that a harmonic of 4.2% of the DC-component has to be added to remove the desired exciting force. In Fig. 14 the results of the minimization is shown, the 14th harmonic in \(B^2\) is eliminated.

5 Discussion

The time constant of the field winding, available DC-voltage for the rotor inverter, and possible eddy currents on the magnetic circuit limits the maximum field current amplitudes and frequencies that can be modulated.

Table 3 Left: THD of phase voltage with constant field current. Right: THD of phase voltage with minimization of 3rd harmonic in \(U_a\) is presented
Fig. 11
figure 11

Measured current and magnetic flux density waveforms in the 200 kVA generator for one period measured at one point in the airgap with \(I_\mathrm{f}\) modulated to eliminate the 4th harmonic in magnetic flux density squared, with \(w_\mathrm{m}= 2\pi \frac{1}{9}\)

The voltage needed to modulate the 4th harmonic field current, \(\hat{I}_{\mathrm{f},4th} = 1.262~A\), in the experimental test generator at nominal speed is

$$\begin{aligned} U_{\mathrm{f},4\text {th}}&=L\frac{dI_{\mathrm{f},4\text {th}}}{dt} = -L\hat{I}_{\mathrm{f},4\text {th}}4p\omega _\mathrm{m} \sin (4p\omega _\mathrm{m} t + \alpha _4) \nonumber \\ {}&= -3847\sin (4p\omega _\mathrm{m} t + \alpha _4), \end{aligned}$$
(9)

which has a maximal value of 3847 V. For comparison the nominal field voltage is \(U_\mathrm{f} = 36.8\,\mathrm{V}\). Worth noting is the unusually high field winding inductance in the test generator \(L_\mathrm{f} = 1.95\,\mathrm{H}\), due to the large number of turns in the field winding. This estimate does not include the eddy current effect which also limits the bandwidth. For the method to be feasible the need for such high field voltage must be reduced. The eddy currents were not included in the FEM model but makes a big impact on the resulting magnetic flux density from a change in field current. To remedy these problems the following modifications can be done to the electric machine.

  • Increase the field winding supply voltage.

  • Split the field winding into parallel circuits to increase the supply voltage for each pole.

  • Laminate the complete rotor.

Fig. 12
figure 12

Magnitude of harmonics at constant and modulated field current in the 200 kVA generator from simulation and experiment, the 4th harmonic in the magnetic flux density squared is minimized. Top part: Output phase voltage. Bottom part: \(B^2\) at one point on the stator. Harmonic order 1 equals \(\frac{2}{3}\) Hz

Fig. 13
figure 13

Magnitude of harmonics at constant and modulated field current in the 200 kVA generator from simulation, the 4th and 8th harmonic in the magnetic flux density squared is minimized. Top part: Output phase voltage. Bottom part: \(B^2\) at one point on the stator

Table 4 Amplitude and phase of the injected field current for the 33 MVA full-scale generator
Fig. 14
figure 14

Magnitude of harmonics at constant and modulated field current from simulation of the 33 MVA generator, the 14th harmonic in the magnetic flux density squared is minimized. Top part: Output phase voltage. Bottom part: \(B^2\) at one point on the stator

A natural question arises as to the applicability of the suggested method. The need for a very high rate of change in the field current makes it very hard to implement in a normal synchronous machines at rated frequency. The minimization of the exciting force is done on one point per pole. This means that for adjacent tangential positions the exciting force will not be minimized. However, normally the mechanical vibrations form standing waves of different modal shapes on the stator which means that if the force minimization is done in a position where the amplitude peaks (not the standing wave nodes) the method can still be useful to reduce vibration amplitudes. When minimizing for vibration it is important to find a good balance between the amplitude reduction in flux density versus the additional voltage harmonics induced in the stator windings. The voltage harmonic method can be useful for single phase synchronous machines. It is also possible to increase a chosen harmonic using a similar approach to what has been presented here, maybe for so called energy harvesting. The FEM models and minimization shows good correspondence with the experiments. Injection of alternating current in the field winding create additional losses in damper bars and the rotor iron but the effect is small, at the levels used here.

6 Conclusion

In this paper a method to minimize a specific harmonic in one output phase voltage or magnetic flux density harmonic is shown. It was done by injecting harmonics into the field current to change the airgap magnetic flux density. The method has been demonstrated during no-load operation but would work equally well during loaded operation.