The present invention relates to the field of electrical equipment monitoring and diagnosis, especially to the field of monitoring operation status of an electric motor.
With the development of modern industry, more and more electric motors have been used in the manufacturing system. These motors are important transmission elements and executive elements of equipment in the system. Its nominal operation is of great significance to ensure the safe, high-efficiency, high-quality and low-consumption operation of the manufacturing process. An error or fault in the motor will not only damage the motor itself, but also affect the entire production. The error may even endanger safety of a person operating the motor.
Therefore, it is necessary to monitor the operation status of the electric motor, so as to reduce motor failure rate, maintenance costs, avoid dangers, to improve the stability of the system, and also to provide data for improving the performance or design of the electric motor.
In the art, it is common to use electric motor drive application signals to detect abnormal or special behavior. In particular, it is often to read the signals, such as internal drive specific features including, flux, stator frequency, estimated torque as well as external drive-specific features from bus voltage, currents, vibration sensor, position sensor, and then to store these signals for later post-processing and data analysis. The solution in the art can monitor and analyze the operation status of the electric motor. However, such a solution requires high memory and high processing power processor, and the results cannot be obtained in real time.
The present invention aims to improve these drawbacks, especially to effectively and timely monitor the operation status of an electric motor.
In this regard, according to one aspect of the invention, it is provided a method for monitoring operation status of an electric motor in real time, the method comprising:
Preferably, different frequencies can be selected and a spectrum of these frequencies is obtained;
With such a method, it is possible to monitor the operation status of an electric motor in real time. In addition, the spectrum analysis on a wide frequency range can be done using parallel calculation for faster detection. Consequently, the detection time is strongly reduced compared to methods that store data and analyse them after. Moreover, the spectrum analysis can be also performed on different drive signals simultaneously for a more comprehensive monitoring.
In addition, since it is not necessary to store the signals for later post-processing and data analysis, only the computation results of the analyzed spectrum of the signals are stored for comparison and fault detection, if necessary. Therefore, there is no need for a high memory, neither high processing power processor, such as CPU.
Optionally, the step of analyzing the observed frequency to obtain a spectrum thereof is carried out by means of Fourier Transforms, especially by discrete Fourier Transforms (DFT). In particular, the step of analyzing the observed frequency to obtain a spectrum thereof may comprise obtaining discrete time Fourier Transforms component for the frequency of the drive application signal.
In an embodiment, obtaining the discrete Fourier Transforms (DFT) component comprises calculating the following integrators:
where x(t) is the input signal, T is an integration time, f is the frequency to be observed, wherein the integration time is inversely proportional to the corresponding frequency value and it, for example, varies from a few milliseconds to a few seconds, N is the number of iterations for the computation for the DFT component and is determined by the following formular:
where k is a number of periods at the observed frequency, Ts is sample time;
and
calculating the following modulus
|Xf|=2×√{square root over (ReX(f)2+ImX(f)2)}
Optionally, discrete time Fourier Transforms can be run on different frequency ranges simultaneously or in sequence. Once the modulus is obtained, it is possible to compare the calculated modulus of the observed frequencies, so as to detect normal or abnormal operation status of the electric motor.
Preferably, in the abovementioned equations, TS=250 μs, and other values can also be selected, as long as the observed frequency is lower than
Preferably, in the abovementioned equations, T is in a range between 0.5 ms to 100 s, preferably 15 ms and 30 s.
Preferably, in the abovementioned equations, k is in a range between 1 and 100, preferably 30.
Preferably, in the abovementioned equations, the frequency f to be observed is in a range between 0.1 Hz and 2 kHz.
Optionally, the step of detecting the operation status of the electric motor on the basis of the spectrum analysis comprises:
Moreover, the method according to the present invention may further comprise:
Preferably, when an abnormal operation status is detected, the method may further comprise warning a user of the issue, displaying recommendations, modifying motor parameters to mitigate risks, and/or shutting down the motor.
Optionally, according to the method of the present invention, only the analysed spectrum is stored in a memory, such as a RAM, and the drive application signals are discarded. Therefore, there is no need for a high memory, neither high processing power processor.
According to another aspect of the invention, it is also provided a computer program product comprising instructions to be executed by a processor, the instructions being adapted to performance of a method as mentioned above when executed by the processor.
According to a third aspect of the invention, it is further provided a system for monitoring operation status of an electric motor in real time, comprising:
In this regard, the present invention intends to catch the signal frequency signature in real time to detect the operation status of the electric motor, for monitoring, analyzing or diagnosing the electric motor, especially for the early detection of abnormal behavior and the prevention of damage to the motor, drive, or other components in the system, without additional hardware, as well as for monitoring tool for unexpected failures.
Other features and advantages of the present invention will appear in the description hereinafter, in reference to the appended drawings, where:
In one embodiment of the system according to the present invention, a F28M36x microcontroller can be used. It has 2 cores, i.e an ARM M3 core and a DSP C28 core. The algorithm according to the invention (along with the control law) runs on the C28 core.
In step S1, signal X, such as bus voltage, flux, stator speed, estimated torque, and/or currents of the electric motor, is read, and a first frequency to be observed in the signal is selected,
Then, in step S2, the observed frequency is analyzed so as to obtain its corresponding component of a spectrum and the spectrum thereof by means of Fourier Transform, especially by increment Fourier transform integrators, until a preset integration time T is reached, so as to obtain the component of discrete Fourier Transforms (DFT).
In Step S3, a modulus is calculated based on the obtained integrators of discrete Fourier Transforms (DFT).
In Step S4, a new frequency to be observed in the signal is selected, and steps S2 and S3 are repeated.
When all the frequencies to be observed in a frequency range, such as between 0 Hz to 2 kHz, are selected and their corresponding modulus are calculated wherein a spectrum is obtained and analyzed, a result is then output and displayed in Step S5.
Alternatively, in one exemplary implementation, the 0-2 kHz range can be split into several blocks, such as 6, and steps S2 to S4 can run in 6 parallel instances. Accordingly, the number of CPU cycles needed to increment the 12 integrators is 1675 per iteration. When one of the moduli is computed, 2000 cycles are required.
The spectrum analysis of the signal according to the present invention is further detailed in the below algorithm:
Let x be the signal we want to know the spectrum of.
Each DFT component is computed over a time T given as a number k of periods of the observed frequency f. The integrators are updated every sample time TS=250 μs. The total number N of samples used is then:
TS: sampling time, time interval between consecutive iterations of the algorithm (TS=0.25 ms in our implementation)
f: observed frequency, component of the DFT that is being computed, varies from 1 Hz to
(for example, from 1 to 2000 Hz)
k: number of periods used for the DFT computation, integer greater than 1 (for example, from 1 to 100, preferably k=30). If k is high, the accuracy of the computation is high, but it takes more time
total duration of the computation for the component f of the DFT (for example, from 30 s to 15 ms)
total number of iterations needed for the computation for the component f of the DFT (for example, from 120 000 to 60)
Xf is called the component at frequency f of the DFT of x. It is a complex number so we will only consider its module |Xf|. The module is computed through its real and imaginary parts ReX(f) and ImX(f).
Before the computation, a window is needed so that the computed components of the DFT do not interfere with each other. For instance, the window is a Hann window, which consists in multiplying the input signal x by a squared sine function of period T. However, other windows can also be used.
The algorithm then computes the component f of the DFT as follows:
In the above equations, the term sin2
is a windowing term, and the rest is the classic DFT.
After N samples, output the modulus:
|X|=2×√{square root over (ReX(f)2+ImX(f)2)}
In the above equation, the factor “2” compensates the ½ factor introduced by the windowing term.
Then the algorithm stores or displays the value of |Xf|, it changes the value of f and runs again.
In one embodiment, the method according to the present invention can be used to detect a normal operation status of the electric motor, i.e. nominal operating conditions.
For a properly working electric motor, different frequencies of the spectrum are dispatched in several blocks, each block containing between 1 and M frequencies. The frequencies contained in a block are computed by the method according to the present invention one after the other. As an example,
In another embodiment, the method according to the present invention can be used to detect a fault of the electric motor.
When a perturbation or an issue appears, the spectrum obtained by the present invention changes. The detection time is equal to the integration time plus a jitter shorter than the integration time. In this regard, the detection time is lower than twice the integration time. As shown in
It should be noted that in the method according to the invention, for example, during the comparison between components of normal and abnormal signals, measurement tolerance, computational tolerance, and/or system tolerance shall be considered.
After detecting the fault, several reactions are possible to take, for instance to warn the user of the fault, display recommendations, modify drive parameters to mitigate risks; and/or shut down the electric motor.
A further exemplary embodiment according to the present invention is discussed below to show how to monitor operation status of a motor running at 50 Hz and driving a pulley with a belt, especially to demonstrate how to detect the belt of the motor is tight or loose belt by means of present invention.
In this embodiment, the spectrum analysis covering 0 Hz to 2 kHz for squared voltages sum is acquired after 30 s using parallel computation based on the method of the present invention. In other words, the whole spectrum covering 0 Hz to 2 kHz is computed.
After the calculation according to the present invention, it is found that, among all the frequencies, the amplitude of the 100 Hz component of the spectrum of the squared voltages sum is different between a motor driving a pulley with a tight belt and a motor driving a pulley with a loose belt. The results are shown in table 1 below, wherein for the frequency with the same amplitude, only 50 Hz component, is shown, to avoid redundancy.
Therefore, only the computation on the 50 Hz and 100 Hz component is now discussed as follows:
let s be the observed signal:
s(t)=v1(t)2+v2(t)2+v3(t)2
We can rewrite s as follows:
s(t)=A×cos(2π50t)+B(t)
where A is the amplitude of the 50 Hz component of s, and B is the rest of the signal.
The algorithm computes the spectrum of s. Let's focus on the 50 Hz component. The integration time is chosen to be 30 periods long, so T=30/50=0.6 s. The integrators are calculated every 250 μs, so the total number of points is N=2400. With t incrementing by 250 μs at each step, we have:
We replace s by its expression. After some trigonometry we obtain the following equations:
Where C is the part
C will be close to 0 since it does not contain the 50 Hz component. We will not focus on the imaginary part, it is easy to show that it is close to 0 but it takes a lot of trigonometry formulas. Mathematically it should be 0 but the algorithm works with integers, there are rounding errors.
The sum of the term
will be close to 0 as well (which can be proven through trigonometry).
The sum of a cos( ) function over an integer number of periods is 0, so what remains is:
The ½ factor comes from the windowing (the sin2 term).
The modulus then becomes:
|S
50 Hz|=2×√{square root over (Re(50 Hz)2+Im(50 Hz)2)}
|S50 Hz|=A
In the example above, we have A=2 with either the tight or loose belt.
Similarly, let s be the observed signal:
s(t)=v1(t)2+v2(t)2+v3(t)2
We can rewrite s as follows:
s(t)=A×cos(2π100t)+B(t)
where A is the amplitude of the 100 Hz component of s, and B is the rest of the signal.
The integration time is chosen to be 30 periods long, so T=30/100=0.3 s. The integrators are calculated every 250 μs, so the total number of points is N=1200. With t incrementing by 250 μs at each step, we have:
After 0.3 s, the integration is complete. We can use the values of the real and imaginary part to compute the modulus of the 100 Hz component of the DFT:
A=2×√{square root over (Re(100 Hz)2+Im(100 Hz)2)}
In the example above, we have A=8 with the tight belt.
If the belt suddenly loses tension, the value of A will change. In the present embodiment, the value of A for 100 Hz changes to 14 for the motor with loose belt.
Furthermore, in the present embodiment, it may take at most 30 s to computer all the components. However, since T=0.3 s, the new value of A could be obtained after 0.3 second, meaning a real time monitoring is achieved.
Alternatively, another exemplary embodiment according to the present invention is to monitor and detect an unbalanced motor running at 50 Hz. In this embodiment, as an example, the current spectrum components of a motor (balanced or unbalanced) is analyzed, and it is found that among 25 Hz, 50 Hz and 300 Hz components, the 25 Hz component of the spectrum differs from the balanced motor one, as shown in Table 2 below.
Moreover, it is known to those skilled in the art, the aforementioned exemplary embodiments according to the present invention can be implemented in many ways, such as program instructions for execution by a processor, as software modules, microcode, as computer program product on computer readable media, as logic circuits, as application specific integrated circuits, as firmware, etc. The embodiments of the invention can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment containing both hardware and software elements. In a preferred embodiment, the invention is implemented in software, which includes but is not limited to firmware, resident software, microcode, etc.
Furthermore, the embodiments of the invention can take the form of a computer program product accessible from a computer-usable or computer-readable medium providing program code for use by or in connection with a computer, processing device, or any instruction execution system. For the purposes of this description, a computer-usable or computer readable medium can be any apparatus that can contain, store, communicate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The medium can be electronic, magnetic, optical, or a semiconductor system (or apparatus or device). Examples of a computer-readable medium include, but are not limited to, a semiconductor or solid state memory, magnetic tape, a removable computer diskette, a RAM, a read-only memory (ROM), a rigid magnetic disk, an optical disk, etc. Current examples of optical disks include compact disk-read-only memory (CD-ROM), compact disk-read/write (CD-RAN) and DVD.
The embodiments described hereinabove are illustrations of this invention. Various modifications can be made to them without leaving the scope of the invention which stems from the annexed claims.
Number | Date | Country | Kind |
---|---|---|---|
21306906.5 | Dec 2021 | EP | regional |