Induction Motor Supplied with Voltage Containing Symmetrical Subharmonics and Interharmonics

: Sinusoidal voltage ﬂuctuations can be considered a speciﬁc result of the occurrence of voltage subharmonics and interharmonics, which are components of low frequency or not being an integer multiple of the frequency of the fundamental voltage harmonic. These components— symmetrical subharmonics and interharmonics—are of the same magnitude, while their frequencies are symmetrical with respect to the fundamental frequency. Depending on their phase angles, various kinds of voltage ﬂuctuations can be distinguished: amplitude modulation, phase modulation and intermediate modulation. In this study, the effect of phase angles on noxious phenomena in induction motors was analyzed. Additionally, torque pulsations and vibrations of an induction motor under sinusoidal voltage ﬂuctuation and a single voltage subharmonic or interharmonic were compared. The investigations were performed with the ﬁnite element method and an experimental method. Among other ﬁndings, it was found that for some phase angles torque pulsations could be about ten times higher than for other angles, roughly corresponding to the amplitude modulation. interharmonic voltage components; m , m ih , m sh are auxiliary variables depending on the amplitude of speed ﬂuctuations, the modulation frequency, the fundamental frequency and the rotational speed.


Introduction
A widely used industrial prime mover is a cage induction motor. Two of its most important advantages are reliability and durability. However, adverse working conditions, such as for example the occurrence of significant power quality disturbances, may result in its premature destruction [1][2][3]. A common power quality disturbance is a voltage waveform distortion connected with the presence of undesirable components in the supply voltage-usually harmonics. Other harmful components that may exist in the voltage waveform are interharmonics and subharmonics. The former are defined as components of the frequency that are not an integer multiple of the frequency of the fundamental harmonics. The latter are components of the frequency less than the fundamental harmonics. Subharmonics are often regarded as a kind of interharmonics (subsynchronous interharmonics), although in some works (e.g., [4][5][6][7]) they are considered separately because of their specific impact on electrical equipment [8]. In this study, subsynchronous components are referred to as subharmonics, and supersynchronous ones as interharmonics.
A specific case of the simultaneous presence of SaIs is sinusoidal voltage fluctuations [8,[22][23][24][25]. The voltage waveform contains two additional components of the same magnitude, while their frequencies are equal to the frequency of the fundamental harmonic (f 1 ) ± the frequency of voltage fluctuations (f f ) [8,[22][23][24][25]. For instance, for the voltage fluctuations of a 10 Hz frequency in a 50 Hz system, the voltage waveform contains SaIs of 40 Hz and 60 Hz frequencies. Considering [22] (where they were called symmetrical interharmonics), in this study such couples of SaIs are referred to as symmetrical subharmonics and interharmonics (SSaIs). It should be noted that, depending on their phase angles [8,22], various kinds of voltage fluctuations can occur: amplitude modulation (AM), phase modulation (PhM) and intermediate modulation.
In spite of the extraordinary harmfulness of SaIs, their permissible values generally have not been introduced into power quality standards and regulations. According to the standard EN 50,160 voltage characteristics of electricity supplied by public distribution systems [37], "Levels are under consideration, pending more experience". Further, in the standard IEEE-519: IEEE Recommended Practices and Requirements for Harmonic Control in Electric Power Systems [38], only limits based on lamp flicker are included in an informative annex. As per the standard [38], due consideration should be given to the effect of interharmonics on various types of electrical equipment, including motors.
The previous works concerning the effect of voltage SaIs on induction motors were generally restricted to either the perfect AM [23,[27][28][29][30] or a single subharmonic/interharmonic component occurring in the voltage waveform [5,26,[31][32][33][34][35]. The exceptions are [24,25], in which currents and speed fluctuations were analyzed for both the above cases. The appropriate calculations were performed by employing the dq transformation. The presented results of the investigations showed that the local maxima of current subharmonics, current interharmonics and speed fluctuations could be significantly lower for the AM than for the injection of single voltage subharmonics.
It should be stressed that phase angles and possible interaction of current SaIs under voltage fluctuations were not taken into consideration in the previous works. At the same time, these issues could be especially important during the analysis of an induction motor. This is so because symmetrical voltage subharmonics and voltage interharmonics result in current components and torque pulsations of the same frequencies. For example, let us consider voltage SSaIs of frequencies 40 Hz and 60 Hz, occurring in a 50 Hz system. The voltage subharmonic of frequency f sh = 40 Hz causes the flow of SaIs currents of frequencies 40 Hz, 60 Hz and torque pulsations of frequency f p = 10 Hz, exactly as voltage interharmonic of frequency f ih = 60 Hz. Consequently, for some phase angles torque pulsations and SaIs currents (caused separately by voltage subharmonic and voltage interharmonic) are expected to supercharge, and for other angles, to weaken. It is also worth mentioning that in the previous papers the effects of the AM, PhM and a single voltage subharmonic/interharmonic on torque pulsations and vibration of induction motors were not compared.
The main purpose of this work is to point out that phase angles significantly affect undesirable phenomena occurring in the induction motor under voltage SSaIs, and in some cases their values are almost as important as the values of voltage subharmonics and interharmonics. Additionally, a comparison of torque pulsations and vibrations caused by sinusoidal voltage fluctuations and a single voltage subharmonic/interharmonic is presented. The investigations were carried out with the finite element method and an experimental method.
The structure of this article is as follows. Sections 2 and 3 describe voltage fluctuations and an analytical model of an induction motor. Section 4 presents the methodology. Sections 5-7 provide the results of the research on currents, torque pulsations and vibrations. Section 8 concludes this study.

Voltage Fluctuations
Voltage fluctuations occur in various power systems, both land and marine [8,39] and are considered as one of the most common and critical power quality issues (based on [8,23,40,41]). The fluctuations are usually defined as fast changes in the rms voltage value [8,40,42]. However, under voltage fluctuations the amplitude as well as the phase can be modulated [8].
In the conditions of voltage fluctuations, the instantaneous voltage value can be described as (based on [8,43]) where V 1A and f 1 are the amplitude and the frequency of the fundamental voltage component, m v (t) and ϕ v (t) are functions modulating the amplitude and the phase. Let us assume that phase angle ϕ 1 of the fundamental component V 1A cos(2πf 1 t + ϕ 1 ) is null and that the frequency of voltage fluctuations is less than f 1 . For the perfect amplitude or phase modulation, (1) takes the forms (based on [8,23,40,42]): A specific case of voltage fluctuations is a sinusoidal modulation of the amplitude and the phase, which can be described as (based on [43]): where k v and k ph are the factors describing the depth of modulation, f m is the modulation frequency, φ v and φ ph are phase angles of the modulating function. For the assumption that the sine of a small angle is equal to its value (expressed in radians), the cosine is equal to unity and after omitting some negligible frequency components, it can be pointed out (based on [8]) that the sinusoidal modulation (4) may be considered as the superposition of the fundamental voltage component and SSaIs [8,22,23]: where a is the value of voltage SSaIs, related to the fundamental voltage harmonic, for the AM a is equal to k v /2, and for the PhM-k ph /2, φ sh and φ ih are phase angles, f sh , f ih are frequencies of subharmonic and interharmonic components, equal to: Note that for the assumed phase angle φ 1 = 0 • , angles φ sh and φ ih are evaluated with reference to the cosines of frequencies f sh , f ih and of null phase angles at t = 0 [22]. Additionally, for the assumption that the frequency of voltage fluctuations is lower than f 1 , both the subharmonic and interharmonic components are of the positive sequence (based on [15]). For the voltage waveform described with (5), the AM occurs if φ sh + φ ih = 0 • , while the PhM occurs if φ sh + φ ih = 180 • [22]. Exemplary waveforms of SSaIs and their superposition with the fundamental component are presented in Figures 1 and 2. The former corresponds to phase angles equal to φ sh = 0 • , φ ih = 0 • (AM) and the latter-to Note that for the assumed phase angle ϕ1 = 0°, angles ϕsh and ϕih are evaluated reference to the cosines of frequencies fsh, fih and of null phase angles at t = 0 [22]. A tionally, for the assumption that the frequency of voltage fluctuations is lower than f1, the subharmonic and interharmonic components are of the positive sequence (base [15]). For the voltage waveform described with (5), the AM occurs if ϕsh + ϕih = 0°, while PhM occurs if ϕsh + ϕih = 180° [22]. Exemplary waveforms of SSaIs and their superpos with the fundamental component are presented in Figures 1 and 2. The former co sponds to phase angles equal to ϕsh = 0°, ϕih = 0° (AM) and the latter-to ϕsh = 0°, ϕih = (PhM).    Note that for the assumed phase angle ϕ1 = 0°, angles ϕsh and ϕih are evaluated with reference to the cosines of frequencies fsh, fih and of null phase angles at t = 0 [22]. Additionally, for the assumption that the frequency of voltage fluctuations is lower than f1, both the subharmonic and interharmonic components are of the positive sequence (based on [15]). For the voltage waveform described with (5), the AM occurs if ϕsh + ϕih = 0°, while the PhM occurs if ϕsh + ϕih = 180° [22]. Exemplary waveforms of SSaIs and their superposition with the fundamental component are presented in Figures 1 and 2. The former corresponds to phase angles equal to ϕsh = 0°, ϕih = 0° (AM) and the latter-to ϕsh = 0°, ϕih = 180° (PhM).    It is worth adding that even the perfect PhM results in little amplitude volatility, the maximal per-unit value of which is [22]: The main reason for voltage fluctuations is current variation [8]. For the current modulation described with (9) [8]: where I 0 is the rms value of the non-fluctuating current component, m i (t) and ϕ i (t) are functions modulating the amplitude and the phase, ψ is the phase angle, the function modulating the voltage amplitude and the phase can be expressed as [8]: where R and X are per-unit fundamental-frequency source impedances, L is per-unit source inductance, P 0 and Q 0 are per-unit non-fluctuating active and reactive power. For the slow current modulation: (in the case of the sinusoidal modulation), (10) and (11) take the form [8]: The above equations show that even the perfect amplitude or phase modulation in current causes both the AM and the PhM in voltage [8]. At the same time, voltage fluctuations are usually considered as the AM (based on [42]). As the PhM is not expected to result in excessive flickers (based on [22]), this approach is validated for the assessment of their severity. In the case of the induction motor, the assumption that voltage fluctuations are limited to the perfect AM, de facto results in the motor analysis for phase angles φ sh + φ ih = 0 • , while their real values may considerably differ.

Analytical Model of Induction Motor under Voltage Fluctuations
Various methods can be used for the analysis of induction motors; for example, they can be based on a transformer-type equivalent circuit, the dq transformation and a dedicated equivalent circuit.
The simplest method employs a transformer-type equivalent circuit [34]. The method is characterized by simplicity of computations and comparatively high accuracy (based on the authors' experience, for example [34]). Its main drawback is the omission of speed fluctuations, which consequently causes that the use of this method to be limited. In the case of positive sequence subharmonics [15], it can be applied only for motors driving a load with a high inertia moment [34]. For an exemplary 3 kW motor investigated in [34], the load's moment should be at least three times greater than the motor's moment. For the lower load moment of inertia, the application of this method may result in an unacceptable computational error.
Speed fluctuations are included in the method based on the dq transformation, which is widely used to examine electrical machinery. For modeling the induction motor under SaIs, the synchronously rotating reference frame is usually applied (for example [24,25,34]). A detailed description of this method is provided in [24,25]. Its disadvantage is the omission of non-linear phenomena appearing in the induction motor.
Ghaseminezhad et al. [27] elaborated and experimentally verified another method that makes use of a dedicated equivalent circuit ( Figure 3). The circuit consists of three parts, constituting coupled circuits for the fundamental, interharmonic and subharmonic components (respectively, the upper, middle and lower part of Figure 3). The appropriate couplings are possible due to the presence of additional electromotive sources depending on voltages, which are denoted in Figure 3 as V x , V xih , V xsh . The circuits together with couplings correspond to the following equations (based on [27]): where L M is mutual inductance, Ω R0 and Ω Rm are the average angular speed and the amplitude of speed fluctuations due to voltage fluctuations, subscripts S, R refer to a stator and a rotor, respectively, and subscripts 1, sh and ih to the fundamental, subharmonic and interharmonic components. The above equations show how voltage subharmonics induce current interharmonics due to speed fluctuations and vice versa. Of note, if the speed fluctuations are neglected (for example, for high load moment of inertia), the circuits become three separate transformertype equivalent schemes for the fundamental, subharmonic and interharmonic components.
In the above methods, the presence of stator and rotor slots is neglected. In the next method, developed by Ghaseminezhad et al. [29], this is taken into account during the magnetic analysis. The calculation results show a good correlation with the experimental and numerical data. According to [29], this method can be applied for optimization, the assessment of the life expectancy of mechanical components and studies on vibrations.
Apart from the analytical methods, field methods are commonly used to investigate induction motors under voltage fluctuations as well, to analyze the generation of SaIs by electrical machinery driving a fluctuating load [9]. The applied field model is presented in the next section.

Methodology
Currents and torque pulsations under voltage SaIs were investigated with th element method (FEM), and vibration with an experimental method.
Numerical computations were carried out for cage induction motors, type SL ML2B (200 kW) and TSg100L-4B (3 kW), referred to in this study as motor1 and m Their name plate parameters are presented in Table 1. During the 2D FEM analy ANSYS Electronics Desktop environment (ANSYS Electromagnetics Suite 18.0.0) MAXWELL-ANSYS environment were employed. The appropriate calculations w formed with a transient-type solver for tau-type meshes. For motor1 this consisted triangle elements, and for motor2, of 22,094 elements. The maximal side dimension finite elements are specified in [32,33]. Of note, ANSYS Electromagnetics Suite 18.

Methodology
Currents and torque pulsations under voltage SaIs were investigated with the finite element method (FEM), and vibration with an experimental method.
Numerical computations were carried out for cage induction motors, type SLgm 315 ML2B (200 kW) and TSg100L-4B (3 kW), referred to in this study as motor1 and motor2. Their name plate parameters are presented in Table 1 8082 triangle elements, and for motor2, of 22,094 elements. The maximal side dimensions of the finite elements are specified in [32,33]. Of note, ANSYS Electromagnetics Suite 18.0.0 was also used to calculate current SaI for motor1, while for motor2 a separate application was employed, generally in compliance with the standard [44]. The parameters of the supply voltage, including phase angles of SSaIs, were set in the used environments. The model parameters were determined on the basis of the construction data and the results of empirical tests [32,33]. It is worth mentioning that one of the co-authors is an employee of the manufacturer of motor1 (Zakład Maszyn Elektrycznych EMIT S.A. Cantoni Group) and took part in the design process of the motor. A comprehensive description of the applied FEM models is provided in [31][32][33][34], together with the experimental verification for motor 2 (motor1 was not tested experimentally under SaIs because of its high rated power).
For the experimental research, a laboratory setup was used, composed of two 3 kW induction motors, a vibration measurement system, a computer power quality analyzer and a programmable power source. Two investigated motors, type 3SIE100L4B, (see Table 1) were attached to a rigid frame. One of them was uncoupled and the other coupled with the DC generator ( Figure 4). In this study, they are referred to as motor3 Case A and motor3 Case B, respectively. For the vibration measurements the Bruel&Kjear (B&K) system was applied. It consisted of a standalone four-channel data acquisition module (B&K type: 3676-B-040), a three-axis accelerometer (B&K type: 4529-B of the following parameters: frequency range: 0.3-12,800 Hz, weight: 14.5 g, sensitivity: 10 mV/ms −2 , maximum shock level peak: 5100 g) with a calibrator (B&K type: 4294) and a computer equipped with the BK Connect software. The accelerometer was mounted to a steel stand (Figure 4) screwed to a motor casing made of cast aluminum. It is worth adding that the accelerometer was calibrated before each measurement and that the experimental investigations were carried out generally in accordance with the main provisions of the standards [44,45]. In particular, the vibration measurements were taken simultaneously in three axes. For each measurement point, the accelerometer indications were recorded three times. Based on statistical metrics such as the mean value, standard deviation, etc., the acceleration waveforms were evaluated to select the most appropriate one. The chosen waveform was passed through a low-pass filter and then converted into vibration velocity. Of note, the broad-band velocity [45] is required by the standards [45,46] for the assessment of vibration severity. A flowchart of the measurement procedure is presented in Figure 5.
The motors under consideration were supplied with an AC programmable power source type Chroma 61512+A615103, of the rated power 36 kVA. It enabled injection of a single subharmonic or interharmonic component into the output voltage and setting the voltage AM. Because of the limitations of the power source, the sinusoidal-shape function modulating voltage could not be set directly, but it was approximated with line segments. Exemplary spectra of the testing voltage are shown in Figure 6. They were determined with two methods: the computational (theoretical) method, and on the basis of the measurements of the voltage produced by the AC power source loaded with the investigated motor. Both the spectra show sufficient consistency.
The simplified diagram of the measurement stand is presented in Figure 7 (based on [31]).   The motors under consideration were supplied with an AC programmable power source type Chroma 61512+A615103, of the rated power 36 kVA. It enabled injection of a single subharmonic or interharmonic component into the output voltage and setting the voltage AM. Because of the limitations of the power source, the sinusoidal-shape function modulating voltage could not be set directly, but it was approximated with line segments. Exemplary spectra of the testing voltage are shown in Figure 6. They were determined with two methods: the computational (theoretical) method, and on the basis of the measurements of the voltage produced by the AC power source loaded with the investigated motor. Both the spectra show sufficient consistency.
The simplified diagram of the measurement stand is presented in Figure 7 (based on [31]).   The motors under consideration were supplied with an AC programmable power source type Chroma 61512+A615103, of the rated power 36 kVA. It enabled injection of a single subharmonic or interharmonic component into the output voltage and setting the voltage AM. Because of the limitations of the power source, the sinusoidal-shape function modulating voltage could not be set directly, but it was approximated with line segments. Exemplary spectra of the testing voltage are shown in Figure 6. They were determined with two methods: the computational (theoretical) method, and on the basis of the measurements of the voltage produced by the AC power source loaded with the investigated motor. Both the spectra show sufficient consistency.
The simplified diagram of the measurement stand is presented in Figure 7 (based on [31]).

Effect of SSaIs on Current
This section presents the effects of voltage fluctuations on current SaIs. As the load moment of inertia has a significant impact on undesirable phenomena occurring in an induction motor under voltage SaIs [31,32,34], the FEM computations were performed for two extreme cases: a load with a negligible moment of inertia and the work with a constant rotational speed. In practice, the constant rotational speed corresponds to the work with the load moment of inertia much greater [34] than the motor's moment. The numerical experiments were carried out for the load torque and the fundamental voltage components equal to the rated values, the angle ϕsh = 0° and voltage SaIs equal to 1%. It is worth adding that, in all the presented diagrams, the rms values of current SaI refer to the rated currents of the investigated motors.
In order to present the properties of the induction motor under voltage SaIs, Figure  8 shows the waveforms of electromagnetic torque, rotational speed and stator current for motor2, fm = 20 Hz, ϕih = 0°, a negligible load moment of inertia and a step change of load torque from 0 to 100%. The fluctuations in the electromagnetic torque, speed and current result from the SaIs influence and non-linear phenomena occurring in the induction motor [23,32,34,[47][48][49]. Transients of induction motors under SaIs will be the subject a separate study.
In the next diagram (Figure 9), the current SaIs versus the angle ϕih are presented for the modulation frequency fm corresponding to rigid-body resonance [9,23,32]. If the frequency of torque pulsations fp (based on [25]): is close to the natural frequency of the rigid-body mode of rotating mass (fNrb) [9,23,32], the speed fluctuations caused by torque pulsations boost current SaIs, which increases

Effect of SSaIs on Current
This section presents the effects of voltage fluctuations on current SaIs. As the load moment of inertia has a significant impact on undesirable phenomena occurring in an induction motor under voltage SaIs [31,32,34], the FEM computations were performed for two extreme cases: a load with a negligible moment of inertia and the work with a constant rotational speed. In practice, the constant rotational speed corresponds to the work with the load moment of inertia much greater [34] than the motor's moment. The numerical experiments were carried out for the load torque and the fundamental voltage components equal to the rated values, the angle φ sh = 0 • and voltage SaIs equal to 1%. It is worth adding that, in all the presented diagrams, the rms values of current SaI refer to the rated currents of the investigated motors.
In order to present the properties of the induction motor under voltage SaIs, Figure 8 shows the waveforms of electromagnetic torque, rotational speed and stator current for motor2, f m = 20 Hz, φ ih = 0 • , a negligible load moment of inertia and a step change of load torque from 0 to 100%. The fluctuations in the electromagnetic torque, speed and current result from the SaIs influence and non-linear phenomena occurring in the induction motor [23,32,34,[47][48][49]. Transients of induction motors under SaIs will be the subject a separate study.
In the next diagram (Figure 9), the current SaIs versus the angle φ ih are presented for the modulation frequency f m corresponding to rigid-body resonance [9,23,32]. If the frequency of torque pulsations f p (based on [25]): is close to the natural frequency of the rigid-body mode of rotating mass (f Nrb ) [9,23,32], the speed fluctuations caused by torque pulsations boost current SaIs, which increases torque pulsations and speed fluctuations even more. Consequently, at the rigid-body resonance, extraordinarily high speed fluctuations and current SaIs may occur. For motor1 and motor2, the frequencies f m corresponding to the rigid-body resonance (determined on the grounds of field computations) are 7 Hz and 30 Hz, respectively. Of note, the rigid-body resonance occurs only for motors driving loads with a comparatively low moment of inertia [34] and for this reason it is not the case for the constant rotational speed.
are of the opposite phase. For this reason, they cancel each other out. On the contrary, for ϕih ≈ 180° the CIs caused independently by voltage subharmonics and voltage interharmonics are of the same phase and add up algebraically. Consequently, their maximal value in Figure 9 is about twice as high as for a single-voltage subharmonic or interharmonic. Further, for motor2 (Figure 9), current SaIs are less affected by the angle ϕih than for motor1. Current SaIs reach the minima for the angles ϕih approximately equal to 100° and 50°, and the maxima for 280° and 230°. The next diagram ( Figure 10) presents the analogous characteristics as in Figure 9, but determined for the constant rotational speed. For both the motors, the angle ϕih has rather a moderate effect on current SaIs. Of note, in Figures 9 and 10 current subharmonics take considerably higher values for motor2 than for motor1.   The characteristics provided in Figure 9 show that current SaIs significantly depend on the angle φ ih . For motor1, the maximal and minimal values of current subharmonics are 22.1% of I rat and 3.3% of I rat for the angles φ ih ≈ 120 • and φ ih ≈ 300 • . Further, current interharmonics (CIs) reach the maximum and the minimum for φ ih ≈ 180 • and φ ih ≈ 0 • , respectively, and the corresponding values are 19.7% of I rat and 1.6% of I rat . It should be noted that, for this case, the CIs caused independently by voltage subharmonics and voltage interharmonics are of similar values (10.6% and 9.1% of I rat ). Additionally, for φ ih ≈ 0 • they are of the opposite phase. For this reason, they cancel each other out. On the contrary, for φ ih ≈ 180 • the CIs caused independently by voltage subharmonics and voltage interharmonics are of the same phase and add up algebraically. Consequently, their maximal value in Figure 9 is about twice as high as for a single-voltage subharmonic or interharmonic. Further, for motor2 ( Figure 9), current SaIs are less affected by the angle φ ih than for motor1. Current SaIs reach the minima for the angles φ ih approximately equal to 100 • and 50 • , and the maxima for 280 • and 230 • . The next diagram ( Figure 10) presents the analogous characteristics as in Figure 9, but determined for the constant rotational speed. For both the motors, the angle φ ih has rather a moderate effect on current SaIs. Of note, in Figures 9 and 10 current subharmonics take considerably higher values for motor2 than for motor1.  The calculated current SaIs versus the frequency of the voltage modulation f m are given in Figures 11-14. Figures 11 and 12 correspond to the load with a negligible moment of inertia, and Figures 13 and 14 to the constant rotational speed. For motor1 and the load with a negligible moment of inertia, current SaIs under the PhM (Figure 12) are up to above ten times higher than for the AM (Figure 11), while for the constant rotational speed (Figures 13 and 14) the differences are less significant. It should be noted that medium and high power induction motors are characterized by relatively low values of winding resistance and the frequency f Nrb [32], which explains the shape of the considered characteristics. Peaks of current result from the presence of the above-mentioned rigid body resonance. Similar peaks were observed in [5,24,25,27,31,32,34]. For example, in [27] the AM of depth k v = 5% and the frequency f m close to 30 Hz caused an increase in the rms value of the stator current from about 2.5 A to c. 3.2 A. It is also worth mentioning that for motor1 the PhM (Figures 12 and 14) results in significantly greater current subharmonics (up to about 55%) than a single voltage subharmonic (see [32]). Further, for motor2, current subharmonics are of comparable values for both the cases of modulation, while current interharmonics are generally considerably higher for the PhM than the AM (up to c. 50%).
( Figures 13 and 14) the differences are less significant. It should be noted that medium and high power induction motors are characterized by relatively low values of winding resistance and the frequency fNrb [32], which explains the shape of the considered characteristics. Peaks of current result from the presence of the above-mentioned rigid body resonance. Similar peaks were observed in [5,24,25,27,31,32,34]. For example, in [27] the AM of depth kv = 5% and the frequency fm close to 30 Hz caused an increase in the rms value of the stator current from about 2.5 A to c. 3.2 A. It is also worth mentioning that for motor1 the PhM (Figures 12 and 14) results in significantly greater current subharmonics (up to about 55%) than a single voltage subharmonic (see [32]). Further, for motor2, current subharmonics are of comparable values for both the cases of modulation, while current interharmonics are generally considerably higher for the PhM than the AM (up to c. 50%).

Effect of SSaIs on Torque Pulsations
The most harmful effects of voltage SaIs on induction motors, such as vibration torsional vibrations [5,20,23,31,32], are connected with torque pulsations [2,3,20]. Th In summary, the effect of sinusoidal voltage fluctuations on currents depend on the phase angles of SSaIs, especially at rigid-body resonance. For some phase angles, the impact of voltage subharmonics and voltage interharmonics adds up, which might lead to almost the duplication of current SaIs in comparison with their values caused by a single voltage subharmonic. Phase angles, for which the highest and the lowest current SaIs occur, differ for various motors. The detailed explanation of the observed phenomena will be the subject of future investigations.

Effect of SSaIs on Torque Pulsations
The most harmful effects of voltage SaIs on induction motors, such as vibration and torsional vibrations [5,20,23,31,32], are connected with torque pulsations [2,3,20]. The results of the investigations on torque pulsations caused by SSaIs are shown below. All the presented characteristics concern the pulsating torque component (∆T p ) of the frequency described with (21), (22) and refer to the rated torques of the motors under investigation. The numerical experiments were performed for the same assumptions as in Section 5.
The characteristics of the pulsating torque component ∆T p versus the angle φ ih are provided in Figures 15 and 16 for frequency f m = 7 Hz (motor1) and 30 Hz (motor1). As was mentioned in Section 4, these frequencies correspond to the rigid-body resonance. For the load with a negligible moment of inertia ( Figure 15) and motor1, the pulsating torque component T p reaches the maximum 44.3% of the rated torque (T rat ) for the angle φ ih ≈ 150 • . The minimal value is T p = 5.0% of T rat for φ ih ≈ 330 • . For motor2, the pulsating torque component T p is significantly greater than for motor1. The maximal and minimal values of the component are 75.1% and 41.5% of T rat and they occur for the angles φ ih ≈ 265 • and φ ih ≈ 90 • , respectively. It is also worth adding that for both the motors the extrema of T p appear for slightly bigger angles φ ih than the extrema of current subharmonics and slightly lesser angles φ ih than the extrema of current interharmonics (see Section 5).
As was pointed out in Section 5, for the constant rotational speed the phase angles have rather a moderate effect on current SaIs. However, the torque component ∆T p shows considerable dependency on the angle φ ih also for the constant rotational speed (Figure 16). The presence of the extrema in Figures 15 and 16 results from the fact that for some phase angles the torque component ∆T p induced by voltage interharmonics adds to the component ∆T p caused by voltage subharmonics, and for other phase angles it subtracts. It is worth adding that for the constant rotational speed the pulsating torque component ∆T p reaches the extrema for the same angles φ ih as for the load a with negligible moment of inertia ( Figure 15).  In the next diagrams (Figures 17 and 18) the characteristics of the p component ΔTp are given versus the frequency fm for the AM (ϕih = 0°) an 180°). Additionally, the torque component ΔTp due voltage subharmonic ure 18 for informative purposes (the analogous characteristics for voltage are presented in [31]). The peaks of the pulsating torque component occur sidered diagrams are caused by the rigid-body resonance (see the previou ilar peaks were observed in [27,31,32]. For an exemplary 1.1 kW motor u depth kv = 5% and frequency fm of about 30 Hz, the torque fluctuations re comparable with the rated torque [27]. In the next diagrams (Figures 17 and 18) the characteristics of the pulsating torque component ∆T p are given versus the frequency f m for the AM (φ ih = 0 • ) and the PhM (φ ih = 180 • ). Additionally, the torque component ∆T p due voltage subharmonic is shown in Figure 18 for informative purposes (the analogous characteristics for voltage interharmonics are presented in [31]). The peaks of the pulsating torque component occurring in the considered diagrams are caused by the rigid-body resonance (see the previous section). Similar peaks were observed in [27,31,32]. For an exemplary 1.1 kW motor under the AM of depth k v = 5% and frequency f m of about 30 Hz, the torque fluctuations reached the level comparable with the rated torque [27].
The results of computations presented in Figures 17 and 18 show that for motor1 the PhM leads to a considerably greater pulsating torque component ∆T p than the AM. For both the considered cases-a load with a negligible moment of inertia and the constant rotational speed-the torque component ∆T p under the PhM (Figure 18) is up to c. four times greater than for the AM (Figure 17). Further, for motor2 the pulsating torque component ∆T p is of a similar value to the AM (Figure 17), PhM ( Figure 18) as well as a single voltage subharmonic ( Figure 19).
In summary, from the point of view of torque pulsation, phase angles are almost as important as the value of voltage SSaIs. For some phase angles, the pulsating torque component ∆T p might be of an even greater order of magnitude than for the others.
ponent ΔTp is of a similar value to the AM (Figure 17), PhM (Figure 18) as well voltage subharmonic (Figure 19).
In summary, from the point of view of torque pulsation, phase angles are important as the value of voltage SSaIs. For some phase angles, the pulsating to ponent ΔTp might be of an even greater order of magnitude than for the others   In summary, from the point of view of torque pulsation, phase angles a important as the value of voltage SSaIs. For some phase angles, the pulsating ponent ΔTp might be of an even greater order of magnitude than for the othe

Effect of SSaIs on Vibration
According to the authors' experience [5,26,31], voltage SaIs cause the tion for motors under no load. At the same time, some motors idle for mos ational life, e.g., under the standard work S6 15% [50]. Therefore, the resul urements presented in this section concern no-load operation. The appro mental investigations were carried out for motor3 and two cases: Case A, motor and Case B, the motor coupled with the unloaded DC generator (se should be noted that all the tests were carried out for voltage fluctuations o equal to 2% or voltage SaIs of 1%. Because of the limitations of the used p power source (see Section 4), the measurements for the PhM and variable were not performed.
General guidelines concerning measurements and evaluation of vib cluded in the standards [45,46]. As the current standard ISO20816-1 Mecha -Measurement and evaluation of machine vibration-Part 1: General g does not contain the univocally specified recommendation concerning the vibration severity, for the purpose of this paper the recommendations incl terior version [46] were adopted.
The broad-band vibration velocity under the AM is provided in Figure  Case

Effect of SSaIs on Vibration
According to the authors' experience [5,26,31], voltage SaIs cause the highest vibration for motors under no load. At the same time, some motors idle for most of their operational life, e.g., under the standard work S6 15% [50]. Therefore, the results of the measurements presented in this section concern no-load operation. The appropriate experimental investigations were carried out for motor3 and two cases: Case A, the uncoupled motor and Case B, the motor coupled with the unloaded DC generator (see Section 4). It should be noted that all the tests were carried out for voltage fluctuations of an amplitude equal to 2% or voltage SaIs of 1%. Because of the limitations of the used programmable power source (see Section 4), the measurements for the PhM and variable phase angles were not performed.
General guidelines concerning measurements and evaluation of vibration are included in the standards [45,46]. As the current standard ISO20816-1 Mechanical vibration-Measurement and evaluation of machine vibration-Part 1: General guidelines [45] does not contain the univocally specified recommendation concerning the assessment of vibration severity, for the purpose of this paper the recommendations included in its anterior version [46] were adopted. Hz. It should be stressed that these vibration levels significantly exceed the threshold value of Zone D (4.5 mm/s [46]), defined in [45,46] for the evaluation of vibration severity. According to [45], the vibration within Zone D is normally considered to be of sufficient severity to cause damage to the machine [45]. For other frequencies f m the vibration velocity is generally considerably lower, but it virtually exceeds the lower boundaries of another evaluation zone, Zone C. Vibrations within this zone (1.8 mm/s-4.5 mm/s [46]) are considered unacceptable for long-term operation [45,46]. Of note for Case A, the vibration velocity ( Figure 20) roughly corresponds to current SaIs caused by voltage fluctuations (Figure 22), while for Case B (Figures 21 and 23) the relevant characteristics differ significantly. It should be stressed that the vibration level is significantly affected by the behavior of the mechanical structure [2,3] and excessive motor electromagnetic vibration is very often a result of a resonance condition on the structure of an entire unit or on the motor components, such as a stator core or frame [2]. Further, the peaks appearing in Figures 20-24 can be explained by the occurrence of rigid body resonance (see the previous sections). It is also worth mentioning that the maximal vibration velocity under the AM approximately corresponds to the vibration velocity caused by a single subharmonic injection (Figures 24 and 25 that the vibration level is significantly affected by the behavior of the mechanical structure [2,3] and excessive motor electromagnetic vibration is very often a result of a resonance condition on the structure of an entire unit or on the motor components, such as a stator core or frame [2]. Further, the peaks appearing in Figures 20-24 can be explained by the occurrence of rigid body resonance (see the previous sections). It is also worth mentioning that the maximal vibration velocity under the AM approximately corresponds to the vibration velocity caused by a single subharmonic injection (Figures 24 and 25). For Case A, the maximal vibration velocity is 6.41 mm/s, while for Case B, 7.49 mm/s.  that the vibration level is significantly affected by the behavior of the mechanical structure [2,3] and excessive motor electromagnetic vibration is very often a result of a resonance condition on the structure of an entire unit or on the motor components, such as a stator core or frame [2]. Further, the peaks appearing in Figures 20-24 can be explained by the occurrence of rigid body resonance (see the previous sections). It is also worth mentioning that the maximal vibration velocity under the AM approximately corresponds to the vibration velocity caused by a single subharmonic injection (Figures 24 and 25). For Case A, the maximal vibration velocity is 6.41 mm/s, while for Case B, 7.49 mm/s.  To sum up, the AM may cause significant vibration of velocity considerably exceeding the boundaries of evaluation Zone D [45,46]. For the investigated motor, both the AM and a single subharmonic injection result in a vibration of similar velocity.
Energies 2022, 15, x FOR PEER REVIEW To sum up, the AM may cause significant vibration of velocity considerably e ing the boundaries of evaluation Zone D [45,46]. For the investigated motor, both t and a single subharmonic injection result in a vibration of similar velocity.   To sum up, the AM may cause significant vibration of velocity considera ing the boundaries of evaluation Zone D [45,46]. For the investigated motor, b and a single subharmonic injection result in a vibration of similar velocity.

Conclusions
Periodical voltage fluctuations can be regarded as the superposition of symmetrical voltage subharmonics and interharmonics (SSaIs). Depending on their phase angles, the

Conclusions
Periodical voltage fluctuations can be regarded as the superposition of symmetrical voltage subharmonics and interharmonics (SSaIs). Depending on their phase angles, the following can be distinguished: amplitude modulation (AM), phase modulation (PhM) and intermediate cases. At the same time, voltage fluctuations are usually considered as the perfect amplitude modulation (based on [42]). This approach apparently seems justified, as the phase modulation has rather a little effect on light flickers (based on [22]), which are believed to be a particularly serious problem caused by voltage fluctuations. In fact, rotating machinery is especially susceptible to voltage fluctuations (or voltage subharmonics and interharmonics). Voltage fluctuations (e.g., in the output voltage in an inverter) may even lead to a mechanical failure of a power train [12,20].
In this study, noxious phenomena occurring in induction motors under SSaIs were analyzed taking into account their phase angles. The presented results of the research show that in some cases phase angles are almost as important as the value of voltage SSaIs. For some phase angles, the torque pulsations could be even about ten times greater than for other phase angles. It should be stressed that for some motors the most harmful effect of voltage SSaIs occurs for the phase angles approximately corresponding to the PhM, and the least detrimental for the phase angles corresponding to the AM. Consequently, from the point of view of induction motors, voltage fluctuations cannot be a priori considered as the perfect AM. This approach may lead to a significant underestimation of the effect of voltage fluctuations on induction motors.
Further, the most significant effect of phase angles on undesirable phenomena due to SSaIs was observed for motor1 (of the rated power 200 kW), while for motor2 (3 kW) it was less considerable. It was also found that voltage fluctuations may result in vibration, the velocity of which is normally considered to be of sufficient severity to cause damage to the machine [45]. Under the AM, the measured vibration velocity was comparable with its value caused by a single subharmonic injection.
The presented results of the investigations lead to two practical conclusions. Firstly, during the analysis of an induction motor under SSaIs of unknown phase angles, a safe approach is to assume such phase angles for which the torque pulsations are caused separately by voltage subharmonics and voltage interharmonics add up algebraically. Secondly, the power quality standards (PQSs) should not contain only the limits for single voltage subharmonics and interharmonics. They should also take into account the cumulative effect of voltage SSaIs; in practice, separate limits for single-voltage subharmonics, interharmonics and SSaIs could be included in the standards. The elaboration of a proposal for PQSs modification will be the subject of future research work.