Next Article in Journal
Typhoon Trajectory Prediction by Three CNN+ Deep-Learning Approaches
Next Article in Special Issue
Automatic Classification of Rotating Rectifier Faults in Brushless Synchronous Machines
Previous Article in Journal
Research on Multiple AUVs Task Allocation with Energy Constraints in Underwater Search Environment
Previous Article in Special Issue
Bearing Faults Diagnosis by Current Envelope Analysis under Direct Torque Control Based on Neural Networks and Fuzzy Logic—A Comparative Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Start-Up and Steady-State Regimes Automatic Separation in Induction Motors by Means of Short-Time Statistics

by
Jonathan Cureño-Osornio
1,
Carlos A. Alvarez-Ugalde
1,
Israel Zamudio-Ramirez
1,2,
Roque A. Osornio-Rios
1,
Larisa Dunai
3,
Dinu Turcanu
4 and
Jose A. Antonino-Daviu
2,*
1
Engineering Faculty, San Juan del Río Campus, Autonomous University of Querétaro, Av. Río Moctezuma 249, San Juan del Río 76807, Querétaro, Mexico
2
Instituto Tecnológico de la Energía, Technical University of Valencia (UPV), Camino de Vera s/n, 46022 Valencia, Spain
3
Department of Graphic Engineering, Technical University of Valencia (UPV), Camino de Vera s/n, 46022 Valencia, Spain
4
Faculty of Electronics and Telecommunications, Technical University of Moldova (UTM), MD-2004 Chisinau, Moldova
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(19), 3850; https://doi.org/10.3390/electronics13193850
Submission received: 21 August 2024 / Revised: 14 September 2024 / Accepted: 26 September 2024 / Published: 28 September 2024

Abstract

:
Induction motors are widely used machines in a variety of applications as primary components for generating rotary motion. This is mainly due to their high efficiency, robustness, and ease of control. Despite their high robustness, these machines can experience failures throughout their lifespan due to various mechanical, electrical, and environmental factors. To prevent irreversible failures and all the implications and costs associated with breakdowns, various methodologies have been developed over the years. Many of these methodologies have focused on analyzing various physical quantities, either during start-up transients or during steady-state operations. This involves the use of specific techniques depending on the focus of the methodology (start-up transients or steady-state) to obtain optimal results. In this regard, it is of great importance to develop methods capable of separating and detecting the start-up transient of the motor from the steady state. This will enable the development of automatic diagnostic methodologies focused on the specific operating state of the motor. This paper proposes a methodology for the automatic detection of start-up transients in induction motors by using magnetic stray flux signals and processing by means of statistical indicators in time-sliding windows, the calculation of variances with a proposed method, and obtaining optimal values for the design parameters by using a Particle Swarm Optimization (PSO). The results obtained demonstrate the effectiveness of the proposed method for the start-up and steady-state regimes automatic separation, which is validated on a 0.746 kW induction motor supplied by a variable frequency drive (VFD).

1. Introduction

Electric motors are widely used across industries due to their robustness, efficiency, and ease of control. According to data from the ABB Group [1], approximately 300 million electric motors are installed worldwide, with an annual increase of 10%. These motors account for 28% of global electricity consumption. Moreover, 90% of these motors operate continuously at maximum speed. Keeping them in good working condition to ensure operational efficiency and to avoid costly downtime is an essential task. Although electric motors are considered robust machines, they are susceptible to failures due to various factors such as environmental conditions, operational stresses, and insufficient maintenance. Moreover, research by García et al. [2] indicates that the performance of these machines can be adversely affected by such malfunctions, leading to elevated energy consumption costs. If left unchecked, these failures can develop into irreversible damage, directly impacting the kinematic and production chains in which they operate. Therefore, it is necessary to develop techniques for monitoring and detecting faults in electric machines using different physical quantities inherent to the motors, such as current analysis, mechanical vibrations, acoustic signals, infrared imaging, and stray magnetic flux, among others [3]. In recent years, two types of methodologies have been developed for detecting faults in induction motors: on the one hand, there are those focused on the steady state, and on the other hand, there are those focused on the transient state. Within these methodologies, different techniques have been applied depending on the working regime of the analyzed machine (i.e., transient state or steady-state) [4]. Then, there arises the need to develop automatic diagnostic strategies focused on the specific operating state of the motor, as it is not feasible for a specialist technician constantly to monitor the machine for faults. In addition, this type of diagnostic practice allows for efficiency in spare parts, repair times, and monetary resources by obtaining a clear idea about the state of the motor in real-time and by knowing the evolution of any existing breakdown in it and predicting the level of severity [5].
Concerning the aforementioned, some methodologies have been developed in order to diagnose faults in electric machines. For example, the authors in [6] developed an online monitoring system for diagnosing faults, utilizing Motor Current Signature Analysis (MCSA) and the Extended Park’s Vector Approach (EPVA). By combining these two techniques, they successfully detected faults and diagnosed the condition of the stator and rotor in electric motors. Similarly, the authors in [7] developed a method to detect incipient mechanical faults in induction motors through vibration analysis using the Fast Fourier Transform (FFT) and artificial intelligence systems. This method includes a neural network algorithm that can identify common faults such as imbalance, misalignment, and bent shafts, and determine the severity of the problem based on the amplitude of vibration components. Likewise, the authors in [8] explored fault detection in 6 kV induction motors used for water pumping. Their analysis focused on motor current signals and magnetic stray flux, which proved to be a more sensitive alternative to many MCSA misdiagnoses. For instance, they highlighted the misinterpretation of faults in rotor air cooling ducts, where broken bars were detected within the motor. Additionally, in [9] the authors describe the use of MCSA and FFT to detect and analyze faults in induction motors caused by unsymmetrical supply voltages and broken rotor bars.
Although techniques based on steady-state signal analysis, mainly using the FFT algorithm, have proven highly effective in detecting diverse faults appearing in induction motors, the sensitivity of this technique poses challenges in identifying other types of faults [10]. It can lead to erroneous diagnostics when analyzing the machine’s steady regime, since the FFT-based diagnosis focuses solely on specific amplitudes and peaks, which can be masked by unrelated phenomena such as axial air ducts and oscillating loads [11]. In recent years, innovative techniques have been developed for analyzing various types of signals in various motor operating regimes, including the start-up transient. This is of great importance as several faults have been reported in the literature that cannot be detected in the steady-state or that tend to generate similar fault frequencies between one or more faults, which makes them difficult to distinguish. However, if the signals are analyzed during the start-up transient, there is a clear difference between the time-frequency pattern, thus allowing false diagnoses to be ruled out and showing a more robust behavior than in the steady state. Thus arises the need to develop an automatic method that allows for the separation of the start-up transient state—acting as the first step in an automatic methodology based on the diagnosis of transient state faults [12]. In this regard, Daviu [13] provides an overview of monitoring electric motors under transient conditions, noting that transient signals analysis has provided valuable diagnostic information, helping to eliminate the false interpretations that may arise with traditional methods. This research concludes that this new method measures and analyzes the corresponding variables using time-frequency transforms. The results of this analysis provide reliable patterns of failure, which are highly unlikely to be caused by other phenomena unrelated to the anomaly. Various investigations, such as [14], have been conducted for fault detection in induction motors using the analysis of different physical quantities during transient start-up. Among the analyzed techniques are Gabor distribution (Time-Frequency Distribution of Gabor), Morlet time-frequency scalogram (Time-Frequency Morlet Scalogram), multiple-signal classification (MUSIC), and Fast Fourier Transform [15]. Consequently, it can be observed that the diagnosis and interest in the analysis of faults in electric motors during transient states have increased in recent years, aiming to address false fault diagnoses and to detect one or more faults in the study of a single magnitude. An example of this kind of methodology is the paper presented in [16] where a study of several works is carried out. The authors base their research on the analysis of motor transient current signatures (MTCSA), demonstrating and highlighting its effectiveness in the diagnosis of faults in this regime. Similarly, in [17] an array of convolutional neural networks (CNN) and the analysis of the current signature during the transient state are used to achieve the automatic detection of broken rotor bars considering different severity levels and highlighting the importance of analyzing the transient regime for applications where the operating regime varies continuously, or in situations where it is required to diagnose the equipment before experiencing a prolonged uptime. Pointing out the diagnosis of faults in transient regimes is the work of Fernandez [18], who mentions how the detection of faults in induction motors operating in non-stationary regimes has become a necessity in the current industry due to the large number of motors fed by inverters. Then, to carry out this diagnosis successfully, it is necessary to use appropriate techniques. These mathematical techniques belong to the time-frequency decomposition (TFD) tools field and must be applied only in the transient state to obtain optimal results. In conclusion, as observed in the works above, it is necessary to apply specialized techniques depending on whether the machine is operating in a steady state or in a start-up transient state in order to have an automatic methodology. Therefore, it is necessary to have a methodology/technique that allows for isolating the signal obtained during the start-up transient from that obtained in the steady state in order to be able to apply the appropriate techniques.
Throughout this section we discussed several methodologies that have proposed using different processing tools in this field. Currently, a way to adapt these technologies to problems related to electric machines and to obtain optimal results has been sought, such as bio-inspired optimization algorithms. See, for example, the work reported in [19] where the authors propose a method for motor fault detection using the wavelet transform (WT), combining a back propagation neural network with a particle swarm optimization (PSO) algorithm analyzing motor current signals with bearing damage, winding damage, and bar breakage—obtaining results close to 97% of correct classification. Another work is shown in [20], using a model reference adaptive system and simulated annealing particle swarm optimization (MRAS-SAPSO) to identify the electrical parameters of the motor, such as stator winding resistance, transverse inductance, and magnetic linkage—obtaining simulation values very close to the real ones. Similarly, in [21], the authors employ an improved particle swarm optimization method (IPM) to facilitate the parameterization of active disturbance rejection controllers in a 2-motor control system, thus improving the response time and robustness of the system. As it was possible to observe in these works, this type of bio-inspired algorithm has been of great help in the diagnosis of faults in different electrical machines, allowing an increase in the percentages of effectiveness of diverse methodologies in development. Due to the interest in studying the transient state of an induction motor, discriminating between steady-state and transient start-up has become necessary. While the literature includes studies on the analysis of start-up transient state, this regime is often manual or relies on known motor start-up conditions due to the employed technologies. In this regard, Valtierra et al. [22] propose an expert system for early detection of broken bars in induction motors, including half-broken bars and one and two broken bars, by analyzing current signals and considering transient start-up and steady-state to enhance result reliability. In the study, different methodologies focused on diagnosing faults applied only in the transient state have been observed, solving a latent need for automatic diagnosis. This work applies an automatic detection and separation of transient start-up and steady-state, using the Hilbert transform and particle swarm optimization method to separate both current regimes. However, there is not yet a technique that allows for the automatic discrimination of the transient state to which to apply these methodologies, which would be the first step toward achieving a fully automated methodology. Unfortunately, while isolating the transient state of an electric motor is crucial for analysis, few works in the literature have been conducted to detect this regime automatically, and most need to specify the technology used for motor start-up. Since the transient regime fluctuates and is unique to each specific motor, developing a novel methodology for automatically detecting start-up transient would be convenient. The main contribution of this paper is a novel methodology for automatically detecting start-up transients in induction motors using magnetic stray flux signals and their processing through the calculation of statistical indicators in sliding time windows, the calculation of variances with a proposed method, and obtaining optimal values for the design parameters by Particle Swarm Optimization (PSO).

2. Materials and Methods

2.1. Magnetic Stray Flux

The researchers in [23] describe stray flux as the magnetic flux emitted outside a machine’s boundaries, produced by the currents in both the stator and the rotor. These currents, and thus the stray flux, change when there is a fault in the electric motor. As noted in [24], the stray flux can be dissected into two magnetic elements: axial and radial. The axial stray flux originates from the currents in either the end windings of the stator or in the cage of the end ring, occurring along the axis of the machine. Conversely, the radial stray flux exists in a plane orthogonal to the axis. According to the findings in [25], it is feasible to detect stray flux components—individually or combined—by positioning sensors appropriately on the motor’s frame. Illustrated in Figure 1, the presumed pathways of the axial and radial field magnetic lines are shown along with the sensor positions for their detection. The sensor at position 1 is for the axial field magnetic, at position 3 for the radial field magnetic, and at position 2 for capturing both types of stray flux.

2.2. Statistical Indicators

The primary objective of employing statistical indicators is to extract relevant information on the behavior of signals in the time domain. When a system operates under various conditions, its signals manifest through diverse statistical parameters. It is anticipated that a state of failure will be evident in these characteristics or indicators [26]. A primary goal of using statistical indicators is to identify failure patterns by measuring parameters of any physical magnitude related to the system, such as vibrations, electric current and voltage, magnetic stray flux, etc. [27]. Statistical indices in the time domain have been used to analyze the steady-state regime of motors, including the mean or average, RMS value, standard deviation, and variance. Table 1 presents various statistical indicators and their mathematical equations. These indicators are recommended for their ability to provide insights into the general trend of system behavior. It is worth noting that these indicators can be easily implemented in signal processing devices due to their low computational demand and minimal memory resources required to compute them. In the equations presented in Table 1, x k represents the data series in the time domain, and N denotes the number of data points per series. These indicators are straightforward to compute, enhancing their applicability in real-time systems where computational efficiency is critical.

2.3. Particle Swarm Optimization (PSO)

The PSO algorithm is one of the variants of evolutionary optimization algorithms, which is inspired by the behavior of swarms of animals, such as flocks or shoal, which exhibit intelligent collective behavior. This algorithm has several advantages in that it is able to deal with complex, high-dimensional problems. The general operation of the PSO is shown in Figure 2, by means of a flowchart, which consists of a group of particles (known as a swarm) navigating a multi-dimensional search space (called design space), where the position of each particle (PBest) depends on the design variables and each one represents a potential solution to the problem. During each iteration, the particles adjust their velocity and position by considering both their current and previous positions (GBest). Finally, the position of the optimal GBest is obtained by performing this process until the convergence condition given by a maximum number of iterations is met, carrying out a comparison between the maximum number of iterations (T) and the current iteration (t). This process allows each particle to be influenced by the entire set of particles in the swarm, with the goal of finding the best global position that optimizes the objective function, seeking both its maximum and minimum [28].

2.4. Short-Time Statistics and Short-Time Variance

Short-time statistics are used to calculate the indicators shown in Table 1. In addition, Figure 3 shows an example of how to obtain the indicators proposed here. As shown there, it is a signal captured in the time domain, over which square sliding windows are generated with or without overlapping occurring between each one of them. The size of the windows is selected by the proposed optimization algorithm (PSO) in order to avoid as much as possible the processing of disturbances in the captured signal. Then, for each window created, the proposed indicators are calculated using the expressions (1) to (9). Finally, a new raw signal corresponding to the statistical analysis in a short time is obtained. This new raw signal will correspond to the statistical amplitude along the acquired signal in the full-time domain, including the start-up transient and the steady state. Likewise, for the calculation of the short-time variance, the same process is carried out concerning the segmentation of the windows and obtaining the number of windows of the particle optimization algorithm (PSO), with the difference that once the signal has been windowed only the variance in each of them is calculated. The input signal at this stage is the signal resulting from the process above (Shor-time statistics) and results in a new signal based on the variance.

2.5. Proposed Variance Method

We have proposed the variance method to detect the point at which the steady state of magnetic stray flux signals begins. This algorithm performs a first-degree polynomial fit between several signal points (i.e., window size). It must be ensured that the number of samples used for the fit is manageable in order to avoid a significant reduction of the total sample size, which could lead to a loss of information. Another critical aspect of this methodology is the proposal of a threshold, which corresponds to the average value of the variances of the last 10% of the signal (this portion is considered to be in a steady state), as shown at the bottom of the flowchart in Figure 4 as the “threshold of last samples”. This percentage can differ for different IMs depending on different aspects, such as the connected load and power supply technology. To this end, the methodology employs a Particle Swarm Optimization (PSO) algorithm, which allows the optimal estimate of the window size to perform the calculation of statistical and non-statistical indicators (Nwins) in the first stage of the algorithm, the window size for the variance calculation (NwinsVariance), and a second threshold to discard a part of the signal noise that could result in false positives as shown in the second block of the flowchart in Figure 4. This approach allows for an optimal threshold that adjusts to the values of each analyzed signal, as seen in the last block of the flowchart in Figure 4 labeled “Optimized threshold”. Once the stationary value of the signal has been determined, each of the selected variance values is compared with the first threshold calculated in the previous step. The transition from the transient state to the steady state is identified by looking for the most significant difference between two consecutive indices that has coincided with the value obtained at the threshold mentioned above, which will correspond to the time that the steady state ends. This process is illustrated graphically in Figure 4. As an additional step, the relative error of each sample was calculated using Equation (10). The real value corresponds to the time value of the start-up ramp programmed in the frequency inverter to carry out the tests, while the approximate value is the value obtained by the proposed methodology.
R e l a t i v e   e r r o r = R e a l   v a l u e A p p r o x i m a t e   v a l u e R e a l   v a l u e     100
This relative error value is subsequently used to calculate the root mean square error (RMSE) of the total of the case studies, which is calculated by means of Equation (11). The RMSE value allows for the percentage of assertiveness of the proposed method.
R M S E = i = 1 n R e a l   v a l u e A p p r o x i m a t e   v a l u e 2 R e a l   v a l u e
The signal received by the proposed variance method as input is the signal directly acquired by the magnetic flux sensor, and upon completion of the above process it delivers a new signal based on the variance as well as the optimized threshold.

3. Proposed Methodology

This part describes the proposed methodology, which is shown by a flow chart in Figure 5. In the first block, a normalization process is carried to the signals generated and stored by the data acquisition system (DAS): axial stray flux ( 1 ), radial stray flux ( 2 ), and axial+radial stray flux ( 3 ). This normalization is performed by scaling the amplitude of the input signal to a range [−1, 1]. Then, a Savitzky-Golay type filtering is applied to this signal with a data window size of 49 points and order 5 for the filter’s polynomial approximation. These parameters are chosen due to the nature of the filter, which requires an odd window size and which is larger than the order to be used for the approximation. Once this process has been carried out, a pre-processed signal is obtained SG ( i ), this signal will be the input of the next block, which is the proposed variance method (described in detail in the previous chapter). As explained above, the output of this block consists of a signal and an optimized threshold, which are used in the last block of automatic start-up detection explained below. The fourth block of the methodology corresponds to the detection of the transient state. For this purpose, a threshold (stable threshold) is established, which is obtained with the absolute value of the average of the last 10% of the total signal, since it is considered that at this point the signal is already in a stable state and will be used as a reference. In addition, the optimized threshold obtained in the previous block is used to overcome the noise peaks that may appear in the stable zone, as they can produce false positives and affect the start-up time detected by the methodology. Once the stable threshold has been obtained, this threshold is compared with each of the amplitudes of the signal resulting from the variance method. For those indices with a value lower than the threshold, a dashed line is drawn as a graphical reference only. Subsequently, once the entire signal has been compared, the most significant difference between two consecutive indices that has fulfilled the condition is sought, which will correspond to the time that the transient state lasts—allowing to find the point where the transient state ends and the steady state begins.

4. Experimental Setup

The experimental setup involves a WEG 00136APE48T induction motor (0.74 kW, 2.9 A, 220 V AC at 60 Hz) connected to a standard alternator, loading it at 25% of its nominal capacity. The IM is supplied via a variable frequency drive (VFD) that manipulates the start-up profiles through various ramps and target frequencies. Magnetic stray flux components: radial, axial, and radial + axial are captured using a triaxial Hall-effect sensor, with data captured at a rate of 1 kHz by a proprietary Data Acquisition system (DAS). This sensor is placed at a distance of 10 cm to avoid saturation and to obtain an accurate measurement of the magnetic flux. The sensor operates with a direct current power supply of 1.7 V to 3.6 V, has a measurement range of 1200 µT, and a sensitivity of 0.042 µT. The setup, illustrated in Figure 6, aims to validate the accuracy of the slope method (described in Section 2.3) for predicting start-up times, comparing these with controlled measurements according to the testbench matrix shown in Table 2. Each configuration is tested under seven conditions to ensure robust statistical analysis and reliability of the results. The matrix includes tests conducted at three different frequencies: 45 Hz, 50 Hz, and 60 Hz, with four different start-up times: 4 s, 6 s, 8 s, and 10 s. In addition to the existing case studies, a further analysis involving a Direct-On-Line (DOL) induction motor was conducted. This allowed for a broader evaluation of the methodology, as the motor’s kinematic chain differed from those used in the other cases.

5. Results and Discussion

This section presents the results obtained by applying the proposed methodology to the samples collected, as detailed in the test matrix in Table 2. The matrix includes tests conducted at three different frequencies: 45 Hz, 50 Hz, and 60 Hz, with four different start-up times: 4 s, 6 s, 8 s, and 10 s, with a ramp-up profile. After acquiring the samples of the start-up profiles, the three components of each test are processed to obtain the start-up times for each of the statistics listed in Table 1.
The signal to be shown as an example in the following section corresponds to a magnetic flux component in the radial direction for a power supply frequency at 50 Hz. Furthermore, the following graphs represent the results obtained from the step-by-step methodology proposed in Section 3, showing the acquired signal without any processing, the normalization and filtering stage, the calculation of statistical and non-statistical indicators, the proposed variance method, the automatic detection block for the transient state, and the percentage error results of each case study.
In this regard, Figure 7 corresponds to the radial magnetic flux signal mentioned above, which does not have any type of processing other than the segmentation of the 30 s with which we worked throughout the proposed methodology.
Following the steps of the proposed methodology, a normalization stage is conducted, allowing scaled signals to the same range of values, which was between −1 to 1. In addition to implementing a Savitzky-Golay filter, this filter allows the numerical data to be smoothed and eliminates noise while preserving the main characteristics of the data. This graph is shown in Figure 8.
The next step corresponds to the results obtained from the calculation of the statistical and non-statistical short-time indicators mentioned in Table 1. As an example, Figure 9 shows the results obtained for the calculation of the arithmetic mean of the normalized and filtered signal, corresponding to the radial stray flux with a power supply frequency of 50 Hz, as shown above. As seen in Figure 9, in this step it is already possible to distinguish when the IM is operating in a steady state, however there is still a certain amount of noise, which could limit the accuracy of the proposed methodology. For this, the proposed variance method is applied.
Next, Figure 10 shows the graph obtained at the end of the signal processing stage, where the short-time variance is applied to the signal previously shown from the calculation of the statistical and non-statistical indicators. This signal clearly shows the peaks of change that exist in the transient state. This is due to the variance calculation because as the amplitude of the variance tends to zero, it can be inferred that this corresponds to the steady state of the signal, while the segment with a more significant variance belongs to the transient state. Once these limits are defined, several vertical lines (shown in red in Figure 11) are drawn for observation and clarification purposes, indicating the points on the graph where there is no considerable change in variance and, therefore, it is considered a steady state. This allows the detection algorithm to focus on finding the part with a more negligible variance difference, which it selects as the definitive point at which the steady-state begins.
The following graph, shown in Figure 11, corresponds to the result of the detection stage for which the signal calculated by the variance is taken, and two horizontal lines are drawn whose values are given by the threshold calculated by the bio-inspired algorithm. This algorithm used ten initial particles and a maximum of five iterations. For each iteration, these particles are moved along the search space, looking for the best position obtained as well as the best global position, which would be defined by the error obtained in each case of study and modifying the three variables to optimize, namely: the window size for the calculation of indicators, the window size for the calculation of variance, and the size of the threshold to avoid noise in the signal. The window sizes are bounded between 1 and 100 with integers, and from 0 to 2 with floats for the threshold, looking for optimal values according to the PSO algorithm.
The following section presents an extra case to those proposed in the test matrix, which consists of the induction motor used in the previous tests but connected to a different kinematic chain which represents 20% load. In this case, the motor was fed at direct voltage. First, in Figure 12 the flux signal acquired can be observed, already normalized and filtered, and then the result obtained by applying the method proposed in this article is shown. It can be seen that the optimized threshold discards the noisy signal that is already in steady state and positions the end of the transient state in an approximate time of 2.39 s when the real steady state time corresponds to 2.5 s. It is important to note that these results were obtained by applying the parameters obtained by the methodology proposed for the previous cases, without the need to perform a new training and allowing for the corroboration of the efficiency of the method.
As final results, Figure 13 shows a box and whiskers plot, resulting from the calculation of the relative error obtained by completing the proposed methodology, calculating the error for each repetition of the four cases of study (4 s, 6 s, 8 s, and 10 s of start-up time), and for the power frequency of 45 Hz. This graph shows the analysis of variance (ANOVA) resulting from the eight repetitions for each start-up ramp used in the experimentation. It is noticeable the relatively low error percentages obtained, having an average of 2.7% error for the case of 4 s, 4.8% error for 6 s, 6.5% error for the case of 8 s, and 6.7% relative error for the 10 s start = up transient.
The next power supply frequency tested was 50 Hz, for which the same four motor starting times and the same number of repetitions are used. Similar to the previous figure, Figure 14 shows the graph obtained for the radial magnetic stray flux with the error percentages for these study cases, with a 1.8% error for the 4 s case, 0.5% error for the 6 s, 2.8% error for 8 s, and finally 2.7% error for the 10 s start ramp.
To conclude this section, Figure 15 shows the error percentages of the radial magnetic stray flux for a power frequency of 60 Hz for the cases of the study mentioned above. The results show a 9% error for the 4 s start-up ramp, being the highest error of all the cases, 5.5% error for the 6 s ramp, 1.25% error for the 8 s start-up, and 1.20% error for the 10 s start-up.

6. Conclusions

As could be observed throughout this work, the proposed methodology is capable of automatically detecting the onset of the steady state for an induction motor, with a high level of assertiveness considering the large number of case studies available, since by applying this methodology for all starting times (4 s, 6 s, 8 s, and 10 s) and for the three power supply frequencies (45 Hz, 50 Hz, and 60 Hz) a total mean square error of 11. 95% is obtained. These results are largely due to the application of the particle optimization algorithm that allowed the calculation of the optimum values for the three variables to be controlled (statistical and non-statistical window size, variance window size, and threshold value)—obtaining an optimal value of 26 samples for the statistical and non-statistical window size, 15 samples for the variance window size, and 1.970302843426175 for the threshold value. This represents a contribution in the field of automatic fault detection, as the calculations necessary for the development of this methodology do not represent large computational resources such as memory, so they can be easily implemented in a microcontroller or in a field-programmable gate array (FPGA). This characteristic allows the use of this methodology to be adapted to automatic diagnostic methods that require applying appropriate techniques depending on whether the IM is in a transient state or in the steady state. It should be noted that this methodology promises to be useful to achieve the automatic separation of the transient state, highlighted for its high percentage of assertiveness in general and the low percentages of error obtained in most of the particular cases: for the frequency of 45 Hz an average error of 5.175% was obtained, for the frequency of 50 Hz 1.95% was obtained, and finally for 60 Hz 4.23% was obtained. These percentages were obtained by averaging the results obtained for the four proposed start-up times.

Author Contributions

Conceptualization, I.Z.-R. and J.A.A.-D.; methodology, C.A.A.-U. and J.C.-O.; software, J.C.-O. and C.A.A.-U.; validation, I.Z.-R., D.T. and L.D.; formal analysis, C.A.A.-U.; investigation, J.C.-O., L.D. and D.T.; resources, R.A.O.-R. and J.A.A.-D.; data curation, J.A.A.-D.; writing—original draft preparation, J.C.-O.; writing—review and editing, I.Z.-R. and R.A.O.-R.; visualization, J.C.-O., J.A.A.-D., R.A.O.-R. and D.T.; supervision, I.Z.-R., R.A.O.-R. and J.A.A.-D.; project administration, R.A.O.-R., J.A.A.-D. and L.D.; funding acquisition, J.A.A.-D. and L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish “Ministerio de Ciencia e Innovación”, Agencia Estatal de Investigación, and FEDER program in the framework of the “Proyectos de Generación de Cono-cimiento 2021” of the “Programa Estatal para Impulsar la Investigación Científico-Técnica y su Transferencia”, belonging to the “Plan Estatal de Investigación Científica, Técnica y de Innovación 2021–2023”. (Ref: PID2021-122343OB-I00).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Dataset available on request from the authors.

Acknowledgments

The author Jonathan Cureño-Osornio acknowledges the financial support provided by the Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT) scholarships registered under CVU number 1184202.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Grupo-ABB. Powering the World Economy: Is There a Better Way to Use Electricity? Grupo-ABB: Zurich, Switzerland, 2015. [Google Scholar]
  2. Garcia, M.; Panagiotou, P.A.; Antonino-Daviu, J.A.; Gyftakis, K.N. Efficiency assessment of induction motors operating under different faulty conditions. IEEE Trans. Ind. Electron. 2018, 66, 8072–8081. [Google Scholar] [CrossRef]
  3. Nandi, S.; Toliyat, H.A.; Li, X. Condition monitoring and fault diagnosis of electrical motors—A review. IEEE Trans. Energy Convers. 2005, 20, 719–729. [Google Scholar] [CrossRef]
  4. Lee, S.B.; Hyun, D.; Kang, T.J.; Yang, C.; Shin, S.; Kim, H.; Kim, H.D. Identification of False Rotor Fault Indications Produced by Online MCSA for Medium-Voltage Induction Machines. IEEE Trans. Ind. Appl. 2016, 52, 729–739. [Google Scholar] [CrossRef]
  5. Kumar, R.R.; Andriollo, M.; Cirrincione, G.; Cirrincione, M.; Tortella, A. A Comprehensive Review of Conventional and Intelligence-Based Approaches for the Fault Diagnosis and Condition Monitoring of Induction Motors. Energies 2022, 15, 8938. [Google Scholar] [CrossRef]
  6. Acosta, G.G.; Verucchi, C.; Gelso, E. A current monitoring system for diagnosing electrical failures in induction motors. Mech. Syst. Signal Process. 2006, 20, 953–965. [Google Scholar] [CrossRef]
  7. Saucedo-Dorantes, J.J.; Zamudio-Ramirez, I.; Cureno-Osornio, J.; Osornio-Rios, R.A.; Antonino-Daviu, J.A. Condition monitoring 252 method for the detection of fault graduality in outer race bearing based on vibration-current fusion, statistical features and neural 253 network. Appl. Sci. 2021, 11, 8033. [Google Scholar] [CrossRef]
  8. Gyftakis, K.N.; Panagiotou, P.A.; Spyrakis, D. Detection of simultaneous mechanical faults in 6-kV pumping induction motors using combined MCSA and stray flux methods. IET Electr. Power Appl. 2021, 15, 643–652. [Google Scholar] [CrossRef]
  9. Gaeid, K.S.; Ping, H.W.; Khalid, M.; Salih, A.L. Fault diagnosis of induction motor using MCSA and FFT. Electr. Electron. Eng. 2011, 1, 85–92. [Google Scholar]
  10. Lee, S.B.; Stone, G.; Antonino-Daviu, J.; Gyftakis, K.; Strangas, E.; Maussion, P.; Platero, C. Recent challenges in condition monitoring of industrial electric machines. IEEE Ind. Electron. Mag. 2020, 14, 158–167. [Google Scholar] [CrossRef]
  11. Kim, H.; Lee, S.B.; Park, S.; Kia, S.H.; Capolino, G.-A. Reliable Detection of Rotor Faults Under the Influence of Low-Frequency Load Torque Oscillations for Applications with Speed Reduction Couplings. IEEE Trans. Ind. Appl. 2016, 52, 1460–1468. [Google Scholar] [CrossRef]
  12. Park, Y.; Choi, H.; Shin, J.; Park, J.; Lee, S.B.; Jo, H. Airgap Flux Based Detection and Classification of Induction Motor Rotor and Load Defects During the Starting Transient. IEEE Trans. Ind. Electron. 2020, 67, 10075–10084. [Google Scholar] [CrossRef]
  13. Antonino-Daviu, J. Electrical monitoring under transient conditions: A new paradigm in electric motors predictive maintenance. Appl. Sci. 2020, 10, 6137. [Google Scholar] [CrossRef]
  14. Delgado-Arredondo, P.A.; Garcia-Perez, A.; Morinigo-Sotelo, D.; Osornio-Rios, R.A.; Avina-Cervantes, J.G.; Rostro-Gonzalez, H.; Romero-Troncoso, R.d.J. Comparative study of time-frequency decomposition techniques for fault detection in induction motors using vibration analysis during startup transient. Shock Vib. 2015, 2015, 708034. [Google Scholar] [CrossRef]
  15. Jeon, H.; Jung, Y.; Lee, S.; Jung, Y. Area-Efficient Short-Time Fourier Transform Processor for Time–Frequency Analysis of Non-Stationary Signals. Appl. Sci. 2020, 10, 7208. [Google Scholar] [CrossRef]
  16. Yakhni, M.F.; Cauet, S.; Sakout, A.; Assoum, H.; Etien, E.; Rambault, L.; El-Gohary, M. Variable speed induction motors’ fault detection based on transient motor current signatures analysis: A review. Mech. Syst. Signal Process. 2023, 184, 109737. [Google Scholar] [CrossRef]
  17. Valtierra-Rodriguez, M.; Rivera-Guillen, J.R.; Basurto-Hurtado, J.A.; De-Santiago-Perez, J.J.; Granados-Lieberman, D.; Amezquita-Sanchez, J.P. Convolutional Neural Network and Motor Current Signature Analysis during the Transient State for Detection of Broken Rotor Bars in Induction Motors. Sensors 2020, 20, 3721. [Google Scholar] [CrossRef]
  18. Fernandez-Cavero, V.; Morinigo-Sotelo, D.; Duque-Perez, O.; Pons-Llinares, J. A Comparison of Techniques for Fault Detection in Inverter-Fed Induction Motors in Transient Regime. IEEE Access 2017, 5, 8048–8063. [Google Scholar] [CrossRef]
  19. Lee, C.-Y.; Cheng, Y.-H. Motor Fault Detection Using Wavelet Transform and Improved PSO-BP Neural Network. Processes 2020, 8, 1322. [Google Scholar] [CrossRef]
  20. Su, G.; Wang, P.; Guo, Y.; Cheng, G.; Wang, S.; Zhao, D. Multiparameter Identification of Permanent Magnet Synchronous Motor Based on Model Reference Adaptive System—Simulated Annealing Particle Swarm Optimization Algorithm. Electronics 2022, 11, 159. [Google Scholar] [CrossRef]
  21. Deng, Y.; Zhu, J.; Liu, H. The Improved Particle Swarm Optimization Method: An Efficient Parameter Tuning Method with the Tuning Parameters of a Dual-Motor Active Disturbance Rejection Controller. Sensors 2023, 23, 8605. [Google Scholar] [CrossRef]
  22. Valtierra-Rodriguez, M.; Rivera-Guillen, J.R.; De Santiago-Perez, J.J.; Perez-Soto, G.I.; Amezquita-Sanchez, J.P. Expert System Based on Autoencoders for Detection of Broken Rotor Bars in Induction Motors Employing Start-Up and Steady-State Regimes. Machines 2023, 11, 156. [Google Scholar] [CrossRef]
  23. Bellini, A.; Concari, C.; Franceschini, G.; Tassoni, C.; Toscani, A. Vibrations, currents and stray flux signals to asses induction motors rotor conditions. In Proceedings of the IECON 2006-32nd Annual Conference on IEEE Industrial Electronics, Paris, France, 7–10 November 2006; pp. 4963–4968. [Google Scholar]
  24. Romary, R.; Roger, D.; Brudny, J.F. Analytical computation of an AC machine external magnetic field. Eur. Phys. J.-Appl. Phys. 2009, 47, 31102. [Google Scholar] [CrossRef]
  25. Romary, R.; Pusca, R.; Lecointe, J.; Brudny, J. Electrical machines fault diagnosis by stray flux analysis. In Proceedings of the 2013, IEEE Workshop on Electrical Machines Design, Control and Diagnosis (WEMDCD), Paris, France, 11–12 March 2013; pp. 247–256. [Google Scholar]
  26. Caesarendra, W.; Tjahjowidodo, T. A review of feature extraction methods in vibration-based condition monitoring and its application for degradation trend estimation of low-speed slew bearing. Machines 2017, 5, 21. [Google Scholar] [CrossRef]
  27. Huerta-Rosales, J.R.; Granados-Lieberman, D.; Garcia-Perez, A.; Camarena-Martinez, D.; Amezquita-Sanchez, J.P.; Valtierra Rodriguez, M. Short-circuited turn fault diagnosis in transformers by using vibration signals, statistical time features, and support vector machines on FPGA. Sensors 2021, 21, 3598. [Google Scholar] [CrossRef]
  28. Ramírez-Ochoa, D.-D.; Pérez-Domínguez, L.A.; Martínez-Gómez, E.-A.; Luviano-Cruz, D. PSO, a Swarm Intelligence-Based Evolutionary Algorithm as a Decision-Making Strategy: A Review. Symmetry 2022, 14, 455. [Google Scholar] [CrossRef]
Figure 1. Magnetic stray flux components: (a) radial stray flux; (b) axial stray flux.
Figure 1. Magnetic stray flux components: (a) radial stray flux; (b) axial stray flux.
Electronics 13 03850 g001
Figure 2. Flowchart of the PSO algorithm.
Figure 2. Flowchart of the PSO algorithm.
Electronics 13 03850 g002
Figure 3. Example of the calculation of Short-time Statistics in the signal obtained.
Figure 3. Example of the calculation of Short-time Statistics in the signal obtained.
Electronics 13 03850 g003
Figure 4. Flowchart of proposed slope method technique.
Figure 4. Flowchart of proposed slope method technique.
Electronics 13 03850 g004
Figure 5. Proposed method flowchart.
Figure 5. Proposed method flowchart.
Electronics 13 03850 g005
Figure 6. Experimental setup: (a) kinematic chain and IM starters; (b) side view of the installed sensor; (c) front view of kinematic chain.
Figure 6. Experimental setup: (a) kinematic chain and IM starters; (b) side view of the installed sensor; (c) front view of kinematic chain.
Electronics 13 03850 g006
Figure 7. Radial magnetic flux signal for 50 Hz power supply.
Figure 7. Radial magnetic flux signal for 50 Hz power supply.
Electronics 13 03850 g007
Figure 8. Radial magnetic flux signal for 50 Hz power supply, normalized and filtered.
Figure 8. Radial magnetic flux signal for 50 Hz power supply, normalized and filtered.
Electronics 13 03850 g008
Figure 9. Signal obtained from the calculation of the arithmetic mean following the short-time statistics method for the radial flux and a 50 Hz power supply.
Figure 9. Signal obtained from the calculation of the arithmetic mean following the short-time statistics method for the radial flux and a 50 Hz power supply.
Electronics 13 03850 g009
Figure 10. Signal obtained from the calculation of the variance over the arithmetic mean for a power supply frequency of 50 Hz.
Figure 10. Signal obtained from the calculation of the variance over the arithmetic mean for a power supply frequency of 50 Hz.
Electronics 13 03850 g010
Figure 11. Signal used for steady state detection at a power supply frequency of 50 Hz.
Figure 11. Signal used for steady state detection at a power supply frequency of 50 Hz.
Electronics 13 03850 g011
Figure 12. (a) Acquired magnetic flux signal for the direct online test (DOL). (b) Result obtained for the detection of the transient state.
Figure 12. (a) Acquired magnetic flux signal for the direct online test (DOL). (b) Result obtained for the detection of the transient state.
Electronics 13 03850 g012
Figure 13. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 45 Hz power supply cases.
Figure 13. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 45 Hz power supply cases.
Electronics 13 03850 g013
Figure 14. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 50 Hz power supply cases.
Figure 14. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 50 Hz power supply cases.
Electronics 13 03850 g014
Figure 15. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 60 Hz power supply cases.
Figure 15. Results obtained from the ANOVA for the variance analysis of the radial magnetic flux in the 60 Hz power supply cases.
Electronics 13 03850 g015
Table 1. Statistical time-domain features.
Table 1. Statistical time-domain features.
Time-Domain FeaturesMathematical Equations
Arithmetic mean T 1 = 1 N × i = 1 N | X i | (1)
Maximum value T 2 = max ( x ) (2)
RMS T 3 = 1 N × i = 1 N ( x i ) 2 (3)
SRM T 4 = 1 N × i = 1 N | x i | 2 (4)
Standard Deviation T 5 = 1 N × i = 1 N ( x i T 1 ) 2 (5)
Variance T 6 = 1 N × i = 1 N ( x i T 1 ) 2 (6)
Skewness T 7 = ( x i T 1 ) 3 T 5 3 (7)
5th Moment T 8 = ( x i T 1 ) 5 T 5 5 (8)
6th Moment T 9 = ( x i T 1 ) 6 T 5 6 (9)
Table 2. Test Matrix for Database Generation.
Table 2. Test Matrix for Database Generation.
Number of RepetitionsTime of Start-Up Transient (s)
46810
Power supply frequency45 Hz8888
50 Hz8888
60 Hz8888
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cureño-Osornio, J.; Alvarez-Ugalde, C.A.; Zamudio-Ramirez, I.; Osornio-Rios, R.A.; Dunai, L.; Turcanu, D.; Antonino-Daviu, J.A. Start-Up and Steady-State Regimes Automatic Separation in Induction Motors by Means of Short-Time Statistics. Electronics 2024, 13, 3850. https://doi.org/10.3390/electronics13193850

AMA Style

Cureño-Osornio J, Alvarez-Ugalde CA, Zamudio-Ramirez I, Osornio-Rios RA, Dunai L, Turcanu D, Antonino-Daviu JA. Start-Up and Steady-State Regimes Automatic Separation in Induction Motors by Means of Short-Time Statistics. Electronics. 2024; 13(19):3850. https://doi.org/10.3390/electronics13193850

Chicago/Turabian Style

Cureño-Osornio, Jonathan, Carlos A. Alvarez-Ugalde, Israel Zamudio-Ramirez, Roque A. Osornio-Rios, Larisa Dunai, Dinu Turcanu, and Jose A. Antonino-Daviu. 2024. "Start-Up and Steady-State Regimes Automatic Separation in Induction Motors by Means of Short-Time Statistics" Electronics 13, no. 19: 3850. https://doi.org/10.3390/electronics13193850

APA Style

Cureño-Osornio, J., Alvarez-Ugalde, C. A., Zamudio-Ramirez, I., Osornio-Rios, R. A., Dunai, L., Turcanu, D., & Antonino-Daviu, J. A. (2024). Start-Up and Steady-State Regimes Automatic Separation in Induction Motors by Means of Short-Time Statistics. Electronics, 13(19), 3850. https://doi.org/10.3390/electronics13193850

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