CN116203419A - Motor energy efficiency detecting system based on multi-force field coupling action algorithm - Google Patents
Motor energy efficiency detecting system based on multi-force field coupling action algorithm Download PDFInfo
- Publication number
- CN116203419A CN116203419A CN202310496561.4A CN202310496561A CN116203419A CN 116203419 A CN116203419 A CN 116203419A CN 202310496561 A CN202310496561 A CN 202310496561A CN 116203419 A CN116203419 A CN 116203419A
- Authority
- CN
- China
- Prior art keywords
- motor
- frequency
- representing
- follows
- energy efficiency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 30
- 230000001808 coupling effect Effects 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 47
- 230000008878 coupling Effects 0.000 claims abstract description 39
- 238000010168 coupling process Methods 0.000 claims abstract description 39
- 238000005859 coupling reaction Methods 0.000 claims abstract description 39
- 238000001514 detection method Methods 0.000 claims abstract description 39
- 230000005684 electric field Effects 0.000 claims abstract description 31
- 230000007613 environmental effect Effects 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims abstract description 8
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 49
- 238000001228 spectrum Methods 0.000 claims description 35
- 238000012546 transfer Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000011161 development Methods 0.000 claims description 13
- 230000018109 developmental process Effects 0.000 claims description 13
- 230000004907 flux Effects 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 8
- 238000009423 ventilation Methods 0.000 claims description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 6
- 238000012952 Resampling Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 230000033772 system development Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 238000004804 winding Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 4
- 238000005265 energy consumption Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 239000004973 liquid crystal related substance Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 101100391182 Dictyostelium discoideum forI gene Proteins 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000002195 synergetic effect Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
- G01R31/343—Testing dynamo-electric machines in operation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
Abstract
A motor energy efficiency detection system based on a multi-force field coupling action algorithm comprises a multi-force field information coupling module, a motor parameter identification module and a loss detection module; the method comprises the steps of carrying out feature identification and extraction on an electric field and a magnetic field where a motor is located, coupling multidimensional force field environmental factors, carrying out on-line identification through a sensorless method to obtain motor rotating speed and torque parameters, carrying out signal spectrum analysis on the obtained motor rotating speed and torque parameters, calculating the power of the motor, and subtracting stray loss and wind friction loss in motor operation to obtain motor energy efficiency.
Description
Technical Field
The invention relates to the field of motor energy efficiency detection, in particular to a motor energy efficiency detection system based on a multi-force field coupling action algorithm.
Background
At present, the energy problem is increasingly stressed, the climate problem is more and more important, and the energy problem is the focus of attention of various countries in the world. Energy conservation and emission reduction are important directions, the motor is widely applied to various devices, the electricity consumption accounts for most of the total electricity consumption, but the operation efficiency of the motor is low at present, so that the improvement of the operation efficiency of the motor is important in energy conservation and emission reduction. The problem of low motor operation efficiency is larger, but more main use problems are also included, such as low load rate and long service life of the motor. To achieve this objective, it is first necessary to accurately detect the actual operating efficiency of the motor without affecting the normal operation of the motor. The traditional detection method is difficult to finish on site and is more difficult to detect for a small motor, so that the detection method with simple operation is needed.
Disclosure of Invention
The invention aims to provide a motor energy efficiency detection system based on a multi-force field coupling action algorithm so as to solve the problems in the background technology.
In order to achieve the above purpose, a motor energy efficiency detection system based on a multi-force field coupling action algorithm is provided, which comprises a multi-force field information coupling module, a motor parameter identification module and a loss detection module; firstly, performing feature identification and extraction on an electric field and a magnetic field where a motor is positioned, and coupling multidimensional force field environmental factors; secondly, inputting a control signal to the motor, carrying out on-line identification through a sensorless method, acquiring the rotating speed and torque parameters of the motor, carrying out signal spectrum analysis on the acquired rotating speed and torque parameters of the motor, and calculating the power of the motor; and finally, calculating wind friction loss caused by bearing friction and ventilation when the motor rotates and load stray loss generated by the influence of relative motion of stator tooth grooves and rotor tooth grooves and higher harmonic components in a magnetic field in the stator and rotor iron cores, and obtaining energy efficiency after the stray loss and the wind friction loss are removed.
Further, the multi-force field information coupling module performs feature extraction on an electric field environment where the motor is located, performs electric field feature identification, and comprises the following detailed processes:
the electric field in the stator slot of the motor is concentrated at the conductor edge, and the maximum value of the electric field intensity is at the upper chamfer of the motor slot.
Further, the multi-force field information coupling module performs feature extraction on the magnetic field environment where the motor is located, performs magnetic field feature identification, and the detailed process is as follows
The magnetic field in the motor can be divided into an air gap magnetic field, a leakage magnetic field between motor magnetic poles, a leakage magnetic field in a slot, an electromagnetic field at the end of a winding, a magnetic field in a laminated iron core and a magnetic field in a solid rotor; the air gap magnetic field is the coupling of stator and rotor magnetic field, when the motor is in operation, under the action of air gap magnetic field, the armature winding will produce induced electromotive force, the motor magnetic field is symmetrically distributed when the rotor is unloaded under the action of electromagnetic force, the distribution of air gap magnetic field under one pole is a flat top wave magnetic, and the magnetic flux density at the geometric positions of two poles is zero.
Further, the multi-force field information coupling module couples the electric field and the magnetic field environmental factors, and the detailed process is as follows:
the invention provides a two-system coupling model, which comprises coordination and development between systems, wherein the coordination degree between the two systems is represented by a deviation difference coefficient, and the formula is as follows:
for the coefficient of deviation difference between systems X and Y, X represents the electric field system, Y represents the magnetic field system, further:
wherein the method comprises the steps of
M is the definition of two system cooperation schedules,
the value range is between 0 and 1;
deriving a system development model to make
As a state of the art of electric field systems,
and
respectively representing specific indexes and corresponding weights of the electric field system;
indicating the level of development of the magnetic field system,
and
respectively representing indexes and corresponding weights contained in the magnetic field system; weights corresponding to two systems
，
To express, the comprehensive development model of the two systems is as follows:
in summary, the coupling is a combination of two aspects of system coordination and development, and the derived coupling model is:
e is the coupling degree of the two systems;
further, the motor parameter identification module performs low-pass filtering on the acquired signals, and the detailed process is as follows
The butterworth low-pass filter has the largest flat amplitude-frequency characteristic in the pass band, and the transfer function is as follows:
represents the butterworth transfer function,
in order to be a frequency of the light,
for the filter cut-off frequency, N represents the order and i represents the frequency coefficient; converting transfer function into Laplace domain analysis to make
The formula translates into:
after the processing, the Laplace domain plane obtains 2N symmetrical poles, and all poles are distributed on a unit circle taking an origin as a center;
the update transfer function formula is as follows:
representing transfer function coefficients, designing a Butterworth low-pass filter, and establishing a constraint condition formula aiming at a cut-off frequency:
representing the passband cut-off frequency,
representing the stop-band cut-off frequency,
representing the maximum attenuation of the pass-band,
represents the maximum attenuation of the stop band by aboutThe beam condition, the theoretical order, is calculated as follows:
representing the number of theoretical steps in the order,
representing an adjustment coefficient; the theoretical order is rounded up to obtain the final theoretical order
The cut-off frequency is calculated by the following formula
：
Carry-in transfer function
Obtaining a digital filter transfer function z under the sampling frequency:
representing the sampling frequency, followed by simulation analysis by matlab to build a butterworth filter.
Further, the motor parameter identification module uses stator current spectrum analysis based on rotor frequency, and the detailed process is as follows:
performing fast Fourier transform on the filtered signals to obtain a spectrogram, and converting the signals from a time domain to a frequency domain; performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm;
and performing complex modulation frequency shift, and positioning the frequency domain starting point at the coordinate zero point position. Discrete signal x 0 The discrete fourier transform of (n) is:
（k=0，1,2,…,N-1) X (k) represents the discrete spectrum, X (n) represents the complex modulated frequency shifted signal, I is the number of fast fourier transform points, for a refined frequency band
Center frequency
The method comprises the following steps:
for x 0 (n) performing complex modulation to obtain a complex modulation frequency shift signal x (n), wherein the formula is as follows:
in the form of a fourier series,
is the sequence number of the center frequency in the spectrum.
Discrete spectra X (k) and X of X (n) 0 (n) discrete frequency spectrum
Satisfies the following formula:
resampling the original discrete signal to reduce the sampling frequency to
H is a thinning multiple;
resampling to obtain new discrete signals;
ensuring that I is unchanged through zero padding treatment, and performing FFT to obtain a discrete frequency spectrum;
then frequency shift is performed, and the formula is as follows:
g is a resolution calculation function, a ratio method is adopted to carry out spectrum correction, the spectrum function of a window function is f (x), and the window function is symmetrical relative to a Y axis and takes two points on the function
And (3) making:
by passing through
Solving for
I.e. spectral correction
The construction function is as follows:
w represents the ratio of the ordinate with an abscissa interval of 1, and the solution inverse function is:
h is
Will be the inverse function of
The carry-in function is obtained:
solving to obtain
。
Further, the motor parameter identification module uses a rotation speed identification algorithm based on a stator current frequency spectrum to perform online identification by a sensorless method, and the detailed process is as follows:
analyzing current frequency spectrum, extracting rotor frequency, analyzing mechanical characteristic curve of motor, wherein the minimum value of rotor frequency is defined by rotor frequency at rated rotation speed, maximum value is rotation magnetic field frequency, frequency search region is between minimum value and maximum value, and frequency corresponding to maximum amplitude in search region is rotor frequency
： 
Where p is the pole pair number of the motor,
is the rated rotating speed of the motor,
represents the rotor frequency and the rotating magnetic field frequency at the rated rotation speed,
representing the lowest frequency; rotor frequency
Interposed between
Between, according to the rotor frequency
Calculating the rotation speed of the motor
。
。
Further, the motor parameter identification module obtains motor torque by an air gap torque method, and the detailed process is as follows:
according to the fact that the instantaneous input power of the motor is equal to the sum of products of corresponding voltages and currents, the formula is as follows:
p represents the instantaneous input power and,
the type of voltage is indicated and,
representing the current type, the voltage formula is:
representing the component of the stator flux linkage on each phase, r representing the stator resistance, t representing time; the motor stator flux linkage expression is:
carry over into the instantaneous input power equation:
transforming the coordinates, transferring the target coordinate system from the stator to the rotor, performing Clark transformation, and transforming the original coordinate system to
The coordinate system, the formula is as follows:
representative of
The current column vector is obtained under the coordinates,
the conversion coefficient is represented by a number of coefficients,
in order to transform the matrix,
representing the current column vector under the original coordinate system; the conversion formula is brought into Park conversion and is converted into a rotation coordinate system d-q-0 to obtain the following formula:
the synchronous electrical angle of the motor is indicated,
and (3) expressing a current column vector under a rotating coordinate system to obtain a Park equation of the voltage:
wherein, the liquid crystal display device comprises a liquid crystal display device,
the components of the voltage, the current and the stator flux linkage value are respectively on the d axis, the q axis and the 0 axis; improved output electric power
The method comprises the following steps:
the air gap torque expression, representing the output coefficient, is:
which represents the electromagnetic torque and which is used to control the electromagnetic torque,
the system parameters representing the electromagnetic torque, the expression for the input power is:
。
further, the loss detection module detects that additional loss generated in the stator and rotor cores becomes load stray loss and mechanical loss caused by bearing friction and ventilation is called wind friction, and the detailed process is as follows:
the motor stray loss is obtained by adopting a recommended value method, the stray loss of the motor is measured by testing a series of motors, a typical motor stray loss database is formed, and the calculation formula of the motor stray loss is summarized:
as a result of the stray losses,
is rated power; the wind friction calculation formula obtained through the no-load test is as follows:
represents the wind-powered electricity generation,
indicating the rated output power.
Further, the loss detection module detects that additional loss generated in the stator and rotor cores becomes load stray loss and mechanical loss caused by bearing friction and ventilation is called wind friction, and the detailed process is as follows:
according to the input power, output power, stray loss and wind friction of the motor, the specific actual energy efficiency of the motor is calculated, and the calculation formula is as follows:
indicating the motor efficiency.
The invention has the beneficial effects that:
the invention provides a motor energy efficiency detection system based on a multi-force field coupling action algorithm, which comprises a multi-force field information coupling module, a motor parameter identification module and a loss detection module. The invention adopts a coupling model to perform characteristic identification and extraction on an electric field and a magnetic field where a motor is positioned, couples environmental factors of a multidimensional force field, analyzes coordination and development influence degree of two systems of the electric field and the magnetic field on the energy efficiency of the motor, designs a Butterworth low-pass filter, filters an input signal, and uses Laplace domain analysis to obtain sampling frequency. Performing fast Fourier transform on the filtered signals to obtain a spectrogram,designing different window functions to reduce spectrum leakage; and aiming at the dense signal area, performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm. After the current spectrum is obtained, the rotor frequency is extracted, the mechanical characteristic curve of the motor is analyzed, and after the rotor frequency is obtained, the rotating speed of the motor is calculated. The electromagnetic torque is calculated by an air gap torque method, the input voltage and the current quantity are obtained through the two ends of the motor, and then the original coordinate system is converted into the Clark through the Clark change
The coordinate system is converted into a rotating coordinate system d-q-0 by Park transformation. A calculation formula of the stray loss and the wind friction of the motor is obtained through a large number of test methods, and the stray loss and the wind friction of the motor are calculated through the input power and the rated output power of the motor. And subtracting the stray loss and the wind friction from the rated energy efficiency to obtain the actual motor energy efficiency. The invention provides a motor energy efficiency detection system based on a multi-force field coupling action algorithm, which can quickly and simply calculate the actual motor energy efficiency of a motor under the on-site condition, changes the problems of complexity and difficult operability in the traditional motor energy efficiency detection method, and can comprehensively consider the interaction and synergistic effect among all parts in the motor and the influence of external environment, thereby realizing high-precision detection of motor energy consumption, accurately evaluating the motor energy efficiency and providing guidance significance for the improvement of the motor efficiency. The system data acquisition module can acquire and transmit data in real time, and can better know information such as load distribution, heat distribution, abrasion condition and the like of each part of the motor under the condition that normal operation of the motor is not affected, and perform data interaction with equipment such as a computer, so that the system data acquisition module has good expansibility and flexibility, can be customized and adjusted according to actual requirements, and can accurately identify potential hazards and potential problems possibly existing in the motor through analysis of motor operation parameters, thereby effectively maintaining the motor and prolonging the service life of the motor. Has important application value in the field of motor energy consumption detection.
Drawings
The invention will be further described with reference to the accompanying drawings, in which embodiments do not constitute any limitation on the invention, and other drawings can be obtained by one of ordinary skill in the art without undue effort from the following drawings.
Fig. 1 is a schematic diagram of the structure of the present invention.
Detailed Description
The invention is further described in connection with the following examples.
Referring to fig. 1, the present invention aims to provide a motor energy efficiency detection system based on a multi-force field coupling action algorithm, so as to solve the problems set forth in the background art.
In order to achieve the above purpose, a motor energy efficiency detection system based on a multi-force field coupling action algorithm is provided, which comprises a multi-force field information coupling module, a motor parameter identification module and a loss detection module; firstly, performing feature identification and extraction on an electric field and a magnetic field where a motor is positioned, and coupling multidimensional force field environmental factors; secondly, inputting a control signal to the motor, carrying out on-line identification through a sensorless method, acquiring the rotating speed and torque parameters of the motor, carrying out signal spectrum analysis on the acquired rotating speed and torque parameters of the motor, and calculating the power of the motor; finally, the wind friction loss caused by bearing friction and ventilation when the motor rotates and the load stray loss generated by the influence of the relative motion of stator tooth grooves and rotor tooth grooves and the higher harmonic component in the magnetic field in the stator and rotor iron cores are calculated, and the energy efficiency is obtained after the stray loss and the wind friction loss are removed, wherein the processes of each module are as follows:
in the multi-force field information coupling module: the electric field in the stator slot of the motor is mainly concentrated at the edge of a conductor, the maximum value of the electric field intensity is positioned at the chamfer of the upper part of the motor slot, and the maximum value of the electric field intensity is positioned at the bottom of the slot when the winding sequence of the wire group is changed; in the stator winding end insulation, the electric field is concentrated mainly at the slots, and as the voltage amplitude increases, the end electric field starts to increase.
The magnetic field in the motor can be divided into an air gap magnetic field, a leakage magnetic field between motor magnetic poles, a leakage magnetic field in a slot, an electromagnetic field at the end of a winding, a magnetic field in a laminated iron core and a magnetic field in a solid rotor; the air gap field is the coupling of stator and rotor fields, and when the motor is running, the armature winding will generate induced electromotive force under the action of the air gap field. The motor magnetic field is symmetrically distributed in no-load state, the distribution of the air gap magnetic field under one pole is a flat top wave magnetic, and the magnetic flux density at the geometric positions of the two poles is zero.
The invention provides a two-system coupling model, which comprises coordination and development between systems, wherein the coordination degree between the two systems is represented by a deviation difference coefficient, and the formula is as follows:
for the coefficient of deviation difference between systems X and Y, X represents the electric field system, Y represents the magnetic field system, further:
wherein the method comprises the steps of
M is the definition of two system cooperation schedules,
the value range is between 0 and 1;
deriving a system development model, and for an electric field system and a magnetic field system, making
As a state of the art of electric field systems,
and
respectively representing specific indexes and corresponding weights of the electric field system;
indicating the level of development of the magnetic field system,
and
respectively representing indexes and corresponding weights contained in the magnetic field system; weights corresponding to two systems
，
To express, the comprehensive development model of the two systems is as follows:
representing the system development level, the derived coupling model is:
e is the coupling degree of the two systems;
the motor parameter identification module is used for: and (3) carrying out low-pass filtering on the obtained signal, designing a Butterworth low-pass filter with the maximum flat amplitude-frequency characteristic in a passband, and the transfer function is as follows:
represents the butterworth transfer function,
in order to be a frequency of the light,
for the filter cut-off frequency, N represents the order and i represents the frequency coefficient; conversion of transfer function into Laplce domain analysis, order
The formula translates into:
after the processing, the Laplace domain plane obtains 2N symmetrical poles, and all pole distribution takes the origin as the centerIs a unit circle of (2);
updating a transfer function formula:
representing transfer function coefficients, designing a Butterworth low-pass filter, and establishing a constraint condition formula aiming at a cut-off frequency:
representing the passband cut-off frequency,
representing the stop-band cut-off frequency,
representing the maximum attenuation of the pass-band,
representing the maximum attenuation of the stop band, calculating the theoretical order by constraint conditions, and the formula is as follows:
representing the number of theoretical steps in the order,
representing an adjustment coefficient; the theoretical order is rounded up to obtain the final theoretical order
The cut-off frequency is calculated by the following formula
：
Carry-in transfer function
Obtaining a digital filter transfer function z under the sampling frequency:
representing the sampling frequency, followed by simulation analysis by matlab to build a butterworth filter.
Based on stator current spectrum analysis of rotor frequency, performing fast Fourier transform on the filtered signal to obtain a spectrogram, and converting the signal from a time domain to a frequency domain; performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm;
and performing complex modulation frequency shift, and positioning the frequency domain starting point at the coordinate zero point position. Discrete signal x 0 The discrete fourier transform of (n) is:
(k=0, 1,2, …, N-1) X (k) represents discrete spectrum, X (N) represents complex modulated frequency shift signal, I is the number of fast fourier transform points, and the frequency band is required to be refined
Center frequency
The method comprises the following steps:
for x 0 (n) performing complex modulation to obtain a complex modulation frequency shift signal x (n), wherein the formula is as follows:
in the form of a fourier series,
representing the sequence number of the center frequency in the spectrum.
Discrete spectrum of x (n)X (k) and X 0 (n) discrete frequency spectrum
Satisfies the following formula:
then resampling the original discrete signal to reduce the sampling frequency to
H is a thinning multiple;
resampling to obtain new discrete signals;
the zero padding treatment ensures that I is unchanged, FFT is carried out, discrete frequency spectrum is obtained, frequency shift is carried out, and the formula is as follows:
g is a resolution calculation function, a ratio method is adopted to carry out spectrum correction, the spectrum function of a window function is f (x), and the window function is symmetrical relative to a Y axis and takes two points on the function
And (3) making:
by passing through
Solving for
I.e. spectral correction
The construction function is as follows:
w represents the ratio of the ordinate with an abscissa interval of 1, and the solution inverse function is:
h is
Will be the inverse function of
The carry-in function is obtained:
solving to obtain
。
Analyzing current frequency spectrum, extracting rotor frequency, analyzing mechanical characteristic curve of motor, the minimum value of rotor frequency is defined by rotor frequency at rated rotating speed, maximum value is rotating magnetic field frequency, frequency search region is between minimum value and maximum value, and the frequency correspondent to maximum amplitude in search region is rotor frequency
： 
Where p is the pole pair number of the motor,
is the rated rotating speed of the motor,
represents the rotor frequency and the rotating magnetic field frequency at the rated rotation speed,
representing the lowest frequency; rotor frequency
Interposed between
Between, obtain the rotor frequency
Then, the rotation speed of the motor is calculated
。
。
According to the motor, the instantaneous input power is equal to the sum of products of corresponding voltage and current, and the formula is as follows:
p represents the instantaneous input power and,
the type of voltage is indicated and,
representing the current type, the voltage formula is:
representing the component of the stator flux linkage on each phase, r representing the stator resistance, t representing time; the motor stator flux linkage expression is:
carry over into the instantaneous input power equation:
then coordinate transformation is carried out, the target coordinate system is transferred from the stator to the rotor, clark change is firstly carried out, and the original coordinate system is converted into
The coordinate system, the formula is as follows:
representative of
The current column vector is obtained under the coordinates,
the conversion coefficient is represented by a number of coefficients,
in order to transform the matrix,
representing the current column vector under the original coordinate system; the conversion formula is brought into Park conversion and is converted into a rotation coordinate system d-q-0 to obtain the following formula:
the synchronous electrical angle of the motor is indicated,
a Park equation representing the current column vector in the rotating coordinate system and the voltage is obtained:
wherein, the liquid crystal display device comprises a liquid crystal display device,
the components of the voltage, the current and the stator flux linkage value are respectively on the d axis, the q axis and the 0 axis; improved output electric power
The method comprises the following steps:
the output coefficient is represented, and the air gap torque expression is:
which represents the electromagnetic torque and which is used to control the electromagnetic torque,
a system parameter indicative of the electromagnetic torque,the expression for the input power is therefore:
。
at the loss detection module: the motor stray loss is obtained by adopting a recommended value method, the stray loss of the motor is measured by testing a series of motors, a typical motor stray loss database is formed, and then the calculation formula of the motor stray loss is summarized:
as a result of the stray losses,
is rated power; the wind friction calculation formula obtained through the no-load test is as follows:
represents the wind-powered electricity generation,
indicating the rated output power.
After the input power, the output power, the stray loss and the wind friction of the motor are obtained, the specific actual energy efficiency of the motor can be calculated, and the calculation formula is as follows:
indicating the motor efficiency.
The invention has the beneficial effects that:
the invention provides a motor energy efficiency detection system based on a multi-force field coupling action algorithm, which comprises a multi-force field information coupling module, a motor parameter identification module and a loss detection module. The invention adopts a coupling model to conduct characteristic distinction on the electric field and the magnetic field of the motorThe method comprises the steps of identifying and extracting, coupling environmental factors of a multidimensional force field, analyzing the coordination and development influence degree of an electric field system and a magnetic field system on the energy efficiency of a motor, designing a Butterworth low-pass filter, filtering an input signal, and analyzing by using a Laplace domain to obtain sampling frequency. Performing fast Fourier transform on the filtered signals to obtain a spectrogram, and designing different window functions to reduce frequency spectrum leakage; and aiming at the dense signal area, performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm. After the current spectrum is obtained, the rotor frequency is extracted, the mechanical characteristic curve of the motor is analyzed, and after the rotor frequency is obtained, the rotating speed of the motor is calculated. The electromagnetic torque is calculated by an air gap torque method, the input voltage and the current quantity are obtained through the two ends of the motor, and then the original coordinate system is converted into the Clark through the Clark change
The coordinate system is converted into a rotating coordinate system d-q-0 by Park transformation. A calculation formula of the stray loss and the wind friction of the motor is obtained through a large number of test methods, and the stray loss and the wind friction of the motor are calculated through the input power and the rated output power of the motor. And subtracting the stray loss and the wind friction from the rated energy efficiency to obtain the actual motor energy efficiency. The invention provides a motor energy efficiency detection system based on a multi-force field coupling action algorithm, which can quickly and simply calculate the actual motor energy efficiency of a motor under the field condition, solves the problems of complexity and difficult operability in the traditional motor energy efficiency detection method, and can comprehensively consider the interaction and synergistic effect among all parts in the motor and the influence of external environment, thereby realizing high-precision detection of motor energy consumption, accurately evaluating motor energy efficiency and providing guidance significance for the improvement of motor efficiency. The system data acquisition module can acquire and transmit data in real time, better understand information such as load distribution, heat distribution, abrasion condition and the like of each part of the motor under the condition of not influencing the normal operation of the motor, perform data interaction with equipment such as a computer and the like, has good expansibility and flexibility, can customize and adjust according to actual requirements, and can accurately identify possible existence of the motor through analysis of motor operation parametersHidden danger and potential problem to carry out effectual maintenance to the motor, extension motor life. Has important application value in the field of motor energy consumption detection.
The present invention also provides a computer readable storage medium having stored therein at least one instruction that is loaded and executed by a processor to implement the above-described method. The computer readable storage medium may be, among other things, ROM, random access memory, CD-ROM, magnetic tape, floppy disk, optical data storage device, etc. The instructions stored therein may be loaded by a processor in the terminal and perform the methods described above.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (10)
1. The motor energy efficiency detection system based on the multi-force field coupling action algorithm is characterized by comprising a multi-force field information coupling module, a motor parameter identification module and a loss detection module;
the workflow of the system is as follows: firstly, performing feature identification and extraction on an electric field and a magnetic field where a motor is positioned, and coupling multidimensional force field environmental factors; secondly, inputting a control signal into a motor, filtering the input signal through a designed Butterworth low-pass filter, performing fast Fourier transform on the filtered signal, performing signal spectrum analysis on the acquired motor rotating speed and torque parameters, and performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm; calculating the power of the motor to calculate the power of the motor; and finally, calculating wind friction loss caused by bearing friction and ventilation when the motor rotates and load stray loss generated by the influence of relative motion of stator tooth grooves and rotor tooth grooves and higher harmonic components in a magnetic field in the stator and rotor iron cores, and obtaining energy efficiency after the stray loss and the wind friction loss are removed.
2. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 1, wherein the multi-force field information coupling module is used for extracting characteristics of an electric field environment where the motor is located and identifying electric field characteristics.
3. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 2, wherein the multi-force field information coupling module performs feature extraction on a magnetic field environment where the motor is located, performs magnetic field feature identification, and comprises the following detailed procedures:
the motor magnetic field is symmetrically distributed in no-load state, the distribution of the air gap magnetic field under one pole is a flat top wave magnetic, and the magnetic flux density at the geometric positions of the two poles is zero.
4. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 3, wherein the multi-force field information coupling module couples the electric field and the magnetic field environmental factors, and the detailed process is as follows:
the two-system coupling model comprises coordination and development between systems, wherein the degree of coordination between the two systems is represented by a deviation difference coefficient, and the formula is as follows:
for the coefficient of deviation difference between systems X and Y, X represents the electric field system, Y represents the magnetic field system, further: />
Wherein->
M is the definition of two system cooperation schedules,
m is in the range of 0-1;
deriving a system development model to make
For the electric field system development level, +.>
And->
Respectively representing specific indexes and corresponding weights of the electric field system; />
Indicating the level of development of the magnetic field system,/->
And->
Respectively representing indexes and corresponding weights contained in the magnetic field system; />
、/>
The weight corresponding to the system is represented, and the comprehensive development model of the two systems is as follows:
representing the comprehensive development level of two systems, deriving a coupling model as follows: />
E is the coupling degree of the two systems.
5. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 1, wherein the motor parameter identification module performs low-pass filtering on the acquired signals, and the detailed process is as follows:
the butterworth low pass filter transfer function is as follows:
representing the Butterworth transfer function, < >>
For frequency +.>
For the filter cut-off frequency, N represents the order and i represents the frequency coefficient; transfer function is converted into Laplace domain analysis, let ∈>
The Laplace domain plane obtains 2N symmetrical poles, all poles are distributed on a unit circle with an origin as a center, and the update transfer function formula is as follows: />
Representing transfer function coefficients, establishing a constraint condition formula for the cut-off frequency:
represents passband cut-off frequency, < >>
Represents stop band cut-off frequency, < >>
Represents the passband maximum attenuation, < >>
Representing the maximum attenuation of the stop band, calculating the theoretical order by constraint conditions, and the formula is as follows:
representing theoretical order, +.>
Representing an adjustment coefficient; the theoretical order is rounded up to obtain the final theoretical order +.>
Then calculate the cut-off frequency +.>
：/>
Bringing the transfer function to obtain a digital filter transfer function z at the sampling frequency: />
Representing the sampling frequency, performing simulation analysis by matlab, and establishing a Butterworth filter.
6. The motor energy efficiency detection system based on the multi-force field coupling action algorithm according to claim 1, wherein the motor parameter identification module uses a stator current spectrum analysis based on the following detailed procedures: performing fast Fourier transform on the filtered signals to obtain a spectrogram, and performing spectrum refinement by adopting a complex modulation ZoomFFT algorithm; positioning the frequency domain starting point at the coordinate zero point position; discrete signal x 0 The discrete fourier transform of (n) is:
(k=0, 1,2, …, N-1) X (k) represents a discrete spectrum, and X (N) represents complex modulationFrequency shift signal, I is the number of fast Fourier transform points, for refined frequency band
Center frequency +.>
The method comprises the following steps: />
For x 0 (n) performing complex modulation to obtain a complex modulation frequency shift signal x (n), wherein the formula is as follows: />
For Fourier series +.>
Representing the sequence number of the center frequency in the frequency spectrum; discrete spectra X (k) and X of X (n) 0 Discrete spectrum +.>
Satisfies the following formula: />
Resampling the original discrete signal, reducing the sampling frequency to +.>
H is a thinning multiple; resampling to obtain new discrete signals; the zero padding treatment ensures that I is unchanged, FFT is carried out, discrete frequency spectrum is obtained for frequency shift, and the formula is as follows:
g is a resolution calculation function, a ratio method is adopted to carry out spectrum correction, the spectrum function of a window function is f (x), and the window function is symmetrical relative to a Y axis, and two points on the function are taken to be +.>
And (3) making:
by->
Solving->
I.e. spectral correction +.>
The construction function is as follows: />
W represents the ratio of the ordinate with an abscissa interval of 1, and the solution inverse function is: />
h is->
Will->
The carry-in function is obtained: />
Relieve->
。
7. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 6, wherein the motor parameter identification module uses a rotational speed identification algorithm based on a stator current spectrum to perform online identification by a sensorless method, and the detailed process is as follows:
analysis of the current spectrumExtracting rotor frequency, analyzing mechanical characteristic curve of motor, determining minimum value of rotor frequency from rotor frequency at rated rotation speed, maximum value of rotor frequency is rotation magnetic field frequency, frequency search area is between minimum value and maximum value, frequency corresponding to maximum amplitude in search area is rotor frequency
： />
Wherein->
For pole pair number of motor, < >>
For rated rotational speed of the motor, ">
Represents the rotor frequency and the rotating magnetic field frequency at the rated rotational speed, < >>
Representing the lowest frequency; rotor frequency->
Interposed between
Between, calculate motor speed +.>
。
8. The motor energy efficiency detection system based on the multi-force field coupling algorithm of claim 7, wherein the motor parameter identification module obtains motor torque by an air gap torque method, and the detailed process is as follows:
the instantaneous input power formula of the motor is calculated as follows:
p represents instantaneous input power, < >>
Indicating the voltage type>
Representing the current type, the voltage formula is: />
Representing the component of the stator flux linkage on each phase, r representing the stator resistance, t representing time; the motor stator flux linkage expression is: />
Carry over into the instantaneous input power equation:
transforming the coordinates, transferring the target coordinate system from the stator to the rotor, performing Clark transformation, and transforming the original coordinate system to +.>
The coordinate system, the formula is as follows: />
Represents->
Coordinate current column vector, ">
Representing conversion coefficient->
For the conversion matrix +.>
Representing the current column vector under the original coordinate system; the conversion formula is brought into Park conversion and is converted into a rotation coordinate system d-q-0 to obtain the following formula: />
Indicating the synchronous electrical angle of the motor, < >>
And (3) expressing a current column vector under a rotating coordinate system to obtain a Park equation of the voltage:
wherein (1)>
The components of the voltage, the current and the stator flux linkage value are respectively on the d axis, the q axis and the 0 axis; improved output electric power->
The method comprises the following steps: />
The air gap torque expression, representing the output coefficient, is: />
Representing electromagnetic torque +.>
The system parameters representing the electromagnetic torque, the expression for the input power is: />
。
9. The motor energy efficiency detection system based on the multi-force field coupling algorithm according to claim 1, wherein the loss detection module detects that additional loss generated in stator and rotor cores becomes load stray loss and mechanical loss caused by bearing friction and ventilation is called wind friction, and the detailed process is as follows: the motor stray loss is obtained by adopting a recommended value method, the stray loss of the motor is measured by testing a series of motors, a typical motor stray loss database is formed, and the calculation formula of the motor stray loss is summarized:
for stray losses +.>
For rated power, the wind friction calculation formula obtained by the no-load test is as follows: />
Indicating wind consumption->
Indicating the rated output power.
10. A system for detecting motor energy efficiency based on a multi-force field coupling action algorithm according to claim 1, wherein said lossThe detection module detects that additional loss generated in the stator and rotor iron cores becomes load stray loss and mechanical loss caused by bearing friction and ventilation is called wind friction, and the detailed process is as follows: the specific actual energy efficiency of the motor is calculated as follows:
indicating motor efficiency, +.>
Representing the output electric power. />
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310496561.4A CN116203419A (en) | 2023-05-05 | 2023-05-05 | Motor energy efficiency detecting system based on multi-force field coupling action algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310496561.4A CN116203419A (en) | 2023-05-05 | 2023-05-05 | Motor energy efficiency detecting system based on multi-force field coupling action algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116203419A true CN116203419A (en) | 2023-06-02 |
Family
ID=86513372
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310496561.4A Pending CN116203419A (en) | 2023-05-05 | 2023-05-05 | Motor energy efficiency detecting system based on multi-force field coupling action algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116203419A (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761117A (en) * | 2018-07-27 | 2018-11-06 | 江苏省特种设备安全监督检验研究院 | Portable type current detects rotational speed detector |
-
2023
- 2023-05-05 CN CN202310496561.4A patent/CN116203419A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761117A (en) * | 2018-07-27 | 2018-11-06 | 江苏省特种设备安全监督检验研究院 | Portable type current detects rotational speed detector |
Non-Patent Citations (4)
| Title |
|---|
| 宋庆启: "基于电磁力机理的电机能效现场检测方法研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, no. 07, pages 2 - 5 * |
| 廖重斌 等: "环境与经济协调发展的定量评判及其分类体系——以珠江三角洲城市群为例", 《热带地理》, vol. 19, no. 2, pages 171 - 177 * |
| 杨辉: "巴特沃斯低通滤波器的设计与仿真分析", 《电子产品世界》, no. 12, pages 1 - 4 * |
| 梁雯 等: "长三角地区制造业、物流业与新型城镇化发展耦合效应研究", 《武汉商学院学报》, vol. 36, no. 6, pages 30 - 35 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ramadhan et al. | An effect sensitivity harmonics of rotor induction motors based on fuzzy logic | |
| Prince et al. | Modeling, parameter measurement, and control of PMSG-based grid-connected wind energy conversion system | |
| CN102510263B (en) | Method for identifying practical parameters of synchronous generator on basis of load rejection test and numerical difference | |
| Lu et al. | A survey of efficiency-estimation methods for in-service induction motors | |
| CN104991161B (en) | Double-fed aerogenerator stator winding interturn short-circuit models and feature extracting method | |
| CN102998591B (en) | Positioning method of generator rotor winding dynamic inter-turn short circuit fault | |
| CN108761117B (en) | Portable current detection rotating speed tester | |
| CN108183651A (en) | A kind of control method of permanent magnet linear synchronous motor position-sensor-free | |
| Fang et al. | Improved ZSVC‐based fault detection technique for incipient stage inter‐turn fault in PMSM | |
| Sarma et al. | DFIG stator flux‐oriented control scheme execution for test facilities utilising commercial converters | |
| CN107656203B (en) | A kind of loss of excitation fault detection method of magneto alternator | |
| Smolin et al. | The energy method for monitoring the instantaneous state and the formation of a synchronous motor control variables | |
| Sarma et al. | Modeling, analysis, and validation of controller signal interharmonic effects in DFIG drives | |
| CN106610473B (en) | Electrical parameters measure method based on synchronous generator rest frequency response test | |
| Ito et al. | Simulation for design purposes of magnetic fields in turbine-driven generator end region | |
| Mahdi et al. | Parameter identification of a PMSG using a PSO algorithm based on experimental tests | |
| CN116203419A (en) | Motor energy efficiency detecting system based on multi-force field coupling action algorithm | |
| CN116865610A (en) | Permanent magnet synchronous motor control method and system based on supercoiled sliding mode observer | |
| He et al. | Stator current identification in generator among single and composite faults composed by static air‐gap eccentricity and rotor inter‐turn short circuit | |
| Yan et al. | Iron Losses Model for Induction Machines Considering the Influence of Rotational Iron Losses | |
| CN109918862A (en) | A kind of generator method for dynamic estimation filtered based on robust without mark H infinity | |
| CN109471026A (en) | A kind of more damping circuit Identifying Dynamical Parameters methods of solid rotor large size phase modifier | |
| CN106610474B (en) | Parameter test method based on synchronous generator asynchronous rotation frequency response test | |
| Babel et al. | Parametric sensitivity in the analysis and control of permanent magnet synchronous machines | |
| Xia et al. | Inter-turn short circuit fault detection of PMSM based on wavelet packet energy spectrum and CEEMDAN-HT |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230602 |
|
| RJ01 | Rejection of invention patent application after publication |