CN114993670B - Gear speed determination method and system and computer readable storage medium - Google Patents
Gear speed determination method and system and computer readable storage medium Download PDFInfo
- Publication number
- CN114993670B CN114993670B CN202210583340.6A CN202210583340A CN114993670B CN 114993670 B CN114993670 B CN 114993670B CN 202210583340 A CN202210583340 A CN 202210583340A CN 114993670 B CN114993670 B CN 114993670B
- Authority
- CN
- China
- Prior art keywords
- frequency
- gear
- characteristic
- determining
- amplitude
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001228 spectrum Methods 0.000 claims abstract description 73
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本申请实施例提供一种齿轮转速确定方法及其系统及计算机可读存储介质，包括：采集待测时间内待测齿轮的齿轮振动信号；对齿轮振动信号进行傅里叶变换，得到频谱，频谱包括齿轮频率及其对应的幅值；根据频谱，确定齿轮频率对应的频谱能量；确定频谱能量中的最大值对应的齿轮频率为齿轮啮合频率；根据齿轮啮合频率和待测齿轮的齿数，确定待测齿轮在待测时间内的平均齿轮转速。本申请的齿轮转速确定方法计算简便快速，拥有较高的准确性且成本较低。
The embodiment of the present application provides a method for determining gear speed, a system thereof, and a computer-readable storage medium, including: collecting a gear vibration signal of a gear to be tested within a test time; performing Fourier transform on the gear vibration signal to obtain a spectrum, the spectrum including the gear frequency and its corresponding amplitude; determining the spectrum energy corresponding to the gear frequency according to the spectrum; determining the gear frequency corresponding to the maximum value in the spectrum energy as the gear meshing frequency; determining the average gear speed of the gear to be tested within the test time according to the gear meshing frequency and the number of teeth of the gear to be tested. The gear speed determination method of the present application is simple and fast to calculate, has high accuracy, and is low in cost.
Description
Technical Field
The embodiment of the application relates to the field of mechanical transmission, in particular to a gear rotating speed determining method and system and a computer readable storage medium.
Background
Along with the gradual exhaustion of energy sources such as coal, petroleum and the like, people pay more attention to the utilization of renewable energy sources. Wind energy is becoming increasingly important worldwide as a clean renewable energy source. With the continuous development of wind power technology, wind power generation sets are increasingly applied to power systems. Wind power generation sets are large-scale devices that convert wind energy into electrical energy, and are typically located in areas where wind energy resources are abundant.
The existing wind generating set is installed at the wind gap of mountain, wilderness, beach, island and the like, is affected by irregular turning and load-changing wind force and strong gust, is subjected to the influence of severe summer heat, severe cold and extreme temperature difference throughout the year, and various parts of the wind generating set are prone to faults. The gear box is an important component of the wind driven generator, a plurality of gear boxes are applied to the wind driven generator, the main functions of the gear boxes are to transmit power generated by the wind wheel under the action of wind power to the generator and enable the generator to obtain corresponding rotating speeds, and when the gear in the gear box breaks down, the gear is arranged in the main machine, so that the gear is difficult to find when the fault occurs, and serious consequences are extremely easy to cause. Most gear failure methods based on vibration signals need to determine the gear rotation speed first, and the existing gear rotation speed determination methods have low accuracy when the gears are failed.
Disclosure of Invention
The embodiment of the application aims to provide a gear rotating speed determining method and system and a computer readable storage medium, and the method is simple, convenient and quick to calculate and has higher accuracy.
One aspect of an embodiment of the present application provides a gear rotational speed determining method, including:
collecting gear vibration signals of a gear to be tested in the time to be tested;
Performing Fourier transform on the gear vibration signal to obtain a frequency spectrum, wherein the frequency spectrum comprises gear frequency and corresponding amplitude;
According to the frequency spectrum, determining the frequency spectrum energy corresponding to the gear frequency;
Determining the gear frequency corresponding to the maximum value in the frequency spectrum energy as gear meshing frequency;
and determining the average gear rotating speed of the gear to be tested in the time to be tested according to the gear meshing frequency and the tooth number of the gear to be tested.
Optionally, the determining, according to the spectrum, spectral energy corresponding to the gear frequency includes:
and obtaining the product of the amplitude corresponding to each gear frequency and the amplitude corresponding to the frequency doubling of the gear frequency according to the frequency spectrum, and determining the frequency spectrum energy corresponding to the gear frequency.
Optionally, the determining, according to the spectrum, spectral energy corresponding to the gear frequency includes:
Determining a first characteristic frequency and a second characteristic frequency corresponding to each of a plurality of gear frequencies in the frequency spectrum, wherein the first characteristic frequency and the second characteristic frequency are integer multiples of the corresponding gear frequency and are not equal to each other;
Determining a first characteristic amplitude from the frequency spectrum and the first characteristic frequency and a second characteristic amplitude from the frequency spectrum and the second characteristic frequency, and
And determining the product of the first characteristic amplitude and the second characteristic amplitude to obtain the frequency spectrum energy corresponding to the gear frequency.
Optionally, the determining the first characteristic amplitude according to the spectrum and the first characteristic frequency includes:
determining the maximum value of the amplitude corresponding to the frequency range including the first characteristic frequency in the frequency spectrum as the first characteristic amplitude, or determining the amplitude corresponding to the first characteristic frequency as the first characteristic amplitude, and/or
Said determining a second characteristic amplitude from said spectrum and said second characteristic frequency comprises:
And determining the maximum value of the amplitude corresponding to the frequency range including the second characteristic frequency in the frequency spectrum as the second characteristic amplitude, or determining the amplitude corresponding to the second characteristic frequency as the second characteristic amplitude.
Optionally, the difference between the upper limit value of the frequency range and the characteristic frequency included in the frequency range is [1,3], the difference between the characteristic frequency included in the frequency range and the lower limit value of the frequency range is [1,3], and the characteristic frequency is a first characteristic frequency or a second characteristic frequency.
Optionally, the determining the first characteristic frequency and the second characteristic frequency corresponding to each of the plurality of gear frequencies in the frequency spectrum includes:
determining a plurality of frequency multiplication of each gear frequency and frequency multiplication amplitudes corresponding to the plurality of frequency multiplication respectively, wherein the plurality of frequency multiplication is 1 to n frequency multiplication, and n is a positive integer greater than or equal to 2;
For each gear frequency, determining the frequency multiplication of the gear frequency corresponding to the maximum value in a plurality of frequency multiplication amplitudes as the first characteristic frequency; and determining the frequency multiplication of the gear frequency corresponding to a second maximum value in the multiple frequency multiplication amplitudes as the second characteristic frequency.
Optionally, the determining the first characteristic frequency and the second characteristic frequency corresponding to each of the plurality of gear frequencies in the frequency spectrum includes:
determining a gear meshing frequency range according to the rotating speed range of the gear and the number of teeth of the gear;
A plurality of the gear frequencies in the gear mesh frequency range is determined.
Optionally, said determining a plurality of said gear frequencies in said gear mesh frequency range comprises:
determining a start frequency of the gear engagement frequency range, and a plurality of frequencies spaced apart by a fixed frequency from the start frequency of the gear engagement frequency range as a start, as a plurality of the gear frequencies, or
And determining a plurality of frequencies corresponding to the amplitudes larger than an amplitude threshold value in the gear meshing frequency range on the frequency spectrum as a plurality of gear frequencies.
In another aspect, the embodiment of the application provides a gear speed determining system, which comprises one or more processors and is used for realizing the gear speed determining method.
Yet another aspect of an embodiment of the present application provides a computer-readable storage medium. The computer-readable storage medium has stored thereon a program which, when executed by a processor, implements the gear speed determination method as described above.
The gear rotational speed determining method processes the gear vibration signal to obtain the frequency spectrum energy of the gear, uses the frequency spectrum energy as a basis to obtain the gear engagement frequency, and then calculates the gear rotational speed according to the gear engagement frequency. The frequency spectrum energy of the gear can accurately reflect the state of the gear, so that the obtained rotating speed of the gear has higher accuracy, and the gear is simple and convenient to calculate and easy to use.
Drawings
FIG. 1 is a flow chart of a method of determining gear rotational speed according to one embodiment of the present application;
FIG. 2 is a schematic diagram of the frequency spectrum of the embodiment of FIG. 1;
FIG. 3 is a specific flowchart of step S3 in the embodiment shown in FIG. 1;
FIG. 4 is a schematic diagram of spectral energy of the embodiment of FIG. 1;
FIG. 5 is a schematic diagram of the error of the embodiment of FIG. 1;
FIG. 6 is a schematic block diagram of a gear speed determination system according to one embodiment of the application.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus consistent with aspects of the application as detailed in the accompanying claims.
The terminology used in the embodiments of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. Unless defined otherwise, technical or scientific terms used in the embodiments of the present application should be given the ordinary meaning as understood by one of ordinary skill in the art to which the present application belongs. The terms first, second and the like in the description and in the claims, are not used for any order, quantity or importance, but are used for distinguishing between different elements. Likewise, the terms "a" or "an" and the like do not denote a limitation of quantity, but rather denote the presence of at least one. "plurality" or "plurality" means two or more. Unless otherwise indicated, the terms "front," "rear," "lower," and/or "upper" and the like are merely for convenience of description and are not limited to one location or one spatial orientation. The word "comprising" or "comprises", and the like, means that elements or items appearing before "comprising" or "comprising" are encompassed by the element or item recited after "comprising" or "comprising" and equivalents thereof, and that other elements or items are not excluded. The terms "connected" or "connected," and the like, are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. As used in this specification and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.
Most gear fault diagnosis methods based on vibration signals need rotating speed information, and accurate gear rotating speed can be obtained through processing the gear vibration signals, so that cost is reduced. The application provides a gear rotational speed determining method, fig. 1 is a flowchart of a gear rotational speed determining method according to an embodiment of the application, please refer to fig. 1, the gear rotational speed determining method of the application includes steps S1-S5:
Because the gear rotational speed is further obtained by processing the gear vibration signal to obtain the more accurate gear meshing frequency, firstly, in step S1, the gear vibration signal of the gear to be measured in the time to be measured is collected. In some embodiments, the vibration signal of the gear includes a vibration acceleration. An acceleration sensor can be arranged on the outer surface of the gear box to monitor the vibration acceleration of the gear, vibration acceleration data of a time period are collected at fixed time intervals, for example, data are collected every 0.5 hour, the sampling frequency is 51200, and the single collection time is 2.56 seconds. In some embodiments, the captured gear vibration signal is primarily a gear mesh vibration signal and resembles a cosine waveform.
In step S2, fourier transform is performed on the gear vibration signal to obtain a frequency spectrum. The frequency is the abscissa of the gear vibration signal after fourier transform processing, the amplitude is the ordinate, fig. 2 is a schematic diagram of the frequency spectrum of the embodiment shown in fig. 1, please refer to fig. 2, and after fourier transform, the frequency spectrum includes the gear frequency and the amplitude corresponding to the gear frequency.
Since the gear meshing frequency and the amplitude corresponding to the frequency multiplication thereof are generally larger than other gear frequencies, some gear rotational speed determining methods select the frequency corresponding to the position with the largest amplitude in the frequency spectrum as the gear meshing frequency, however, when the gear fails, the frequency corresponding to the position with the largest amplitude may be far from the gear meshing frequency due to the failure. Compared with the method that the gear engagement frequency can be accurately determined through the frequency spectrum energy through the amplitude, the gear rotation speed determination method is high in accuracy, and the gear engagement frequency is determined according to the frequency spectrum energy and then the gear rotation speed is determined.
First, the spectral energy corresponding to each gear frequency needs to be determined, and in step S3, the spectral energy corresponding to the gear frequency is determined according to the spectrum. In some embodiments, the spectral energy corresponding to the gear frequency may be determined from the frequency spectrum by obtaining a product of the magnitude corresponding to each gear frequency and the magnitude corresponding to a frequency doubling thereof. When a fault occurs, a part of gear frequencies may be caused because of the fault, the corresponding amplitude of the gear frequencies is larger than the corresponding amplitude of the gear meshing frequency, but the fault frequency is not the correct gear meshing frequency, and the frequency spectrum energy is determined to be the gear meshing frequency after the corresponding amplitude of each gear frequency is multiplied by the corresponding amplitude of the frequency doubling of the gear frequency, because the corresponding amplitude of the frequency doubling of the gear frequency is smaller and is close to 0, and the corresponding amplitude of the frequency doubling of the gear meshing frequency is larger. The gear rotating speed determining method is high in accuracy.
In other embodiments, in order to further improve the accuracy of the gear rotational speed determining method of the present application, the method for determining spectral energy may be further refined, where step S3 includes steps S31-S35, and fig. 3 is a specific flowchart of step S3 in the embodiment shown in fig. 1.
In step S31, a gear mesh frequency range is determined based on the rotational speed range of the gear and the number of gear teeth. Taking a gear box of a wind driven generator as an example, an input gear of the wind driven generator is meshed with blades of the wind driven generator, the gear meshing frequency range can be determined through detection of the actual rotation speed of the blades of the wind driven generator and the gear ratio of gears meshed with each other in the gear box, in other embodiments, the rotation speed range of the gears can be determined through the power generation of the wind driven generator, and then the gear meshing frequency range is determined, so that the calculation is simpler and more convenient. When determining the spectral energy of the gear in step S2, it may be determined only in the gear engagement frequency range, thus reducing the calculation amount.
In step S32, a plurality of gear frequencies in the gear mesh frequency range are determined. These gear frequencies will be used to calculate the spectral energy. In some embodiments, a starting frequency of the gear mesh frequency range is determined, and a plurality of frequencies, spaced apart by a fixed frequency, starting from the starting frequency of the gear mesh frequency range, are the plurality of gear frequencies. Thus, the frequency in the gear meshing range can be completely calculated to judge the gear meshing frequency. In some embodiments, the accuracy of the speed determination may be further improved by reducing the fixed frequency of the interval. In other embodiments, to reduce the amount of computation, an amplitude threshold, for example, 0.25m/s 2, may be set, and only the frequencies corresponding to the amplitudes greater than the amplitude threshold in the gear mesh frequency range on the frequency spectrum are determined to be the plurality of gear frequencies.
In step S33, a first characteristic frequency and a second characteristic frequency corresponding to each of a plurality of gear frequencies in the frequency spectrum are determined, wherein the first characteristic frequency and the second characteristic frequency are integer multiples of the corresponding gear frequencies and are not equal to each other. In some embodiments, the first characteristic frequency is the frequency itself, the second characteristic frequency is a frequency that is a double of the frequency, and in other embodiments, the first characteristic frequency is the frequency itself, and the second characteristic frequency is a frequency that is a quadruple of the frequency. In still other embodiments, the first characteristic frequency is a frequency multiplication by three of the frequency and the second characteristic frequency is the frequency itself.
For each frequency in the gear meshing frequency range, the characteristic frequency selected by the gear meshing frequency range is not directly related to the multiple of the frequency, in some embodiments, in order to improve accuracy, the specific determination method comprises the steps of firstly determining multiple frequency multiplication and multiple frequency multiplication corresponding to each gear frequency, wherein the multiple frequency multiplication is 1 to n frequency multiplication, n is a positive integer greater than or equal to 2, then determining that the frequency multiplication of the gear frequency corresponding to the maximum value in the multiple frequency multiplication is a first characteristic frequency for each gear frequency, and determining that the frequency multiplication of the gear frequency corresponding to the second maximum value in the multiple frequency multiplication is a second characteristic frequency. When a fault occurs, it is known that a part of gear frequency may be caused because of the fault, the corresponding amplitude is larger than the corresponding amplitude of the gear engagement frequency, but because the fault frequency is not the actual gear engagement frequency, the corresponding amplitude of the frequency multiplication is smaller and is close to 0, the actual gear engagement frequency is the largest amplitude, and because only the magnitude of the frequency spectrum energy is qualitatively compared, two places with the largest amplitude of the frequency multiplication may not be one frequency multiplication or two frequency multiplication. In the illustrated embodiment, n=4, so that the range is reasonable, the calculation is simple, and in other embodiments, the calculation amount can be reduced to 3.
In step S34, a first characteristic amplitude is determined from the frequency spectrum and the first characteristic frequency, and a second characteristic amplitude is determined from the frequency spectrum and the second characteristic frequency. Since the frequency spectrum of the gear frequency has been established in step S2, it may be determined directly that the amplitude corresponding to the first characteristic frequency is the first characteristic amplitude, and that the amplitude corresponding to the second characteristic frequency is the second characteristic amplitude. The calculation is simple and convenient. In other embodiments, due to sampling accuracy or error, a certain error exists between the first characteristic frequency or the second characteristic frequency and the rotation speed, and step S34 may specifically include determining that the maximum value of the amplitude corresponding to the frequency range including the first characteristic frequency in the frequency spectrum is the first characteristic amplitude, and determining that the maximum value of the amplitude corresponding to the frequency range including the second characteristic frequency in the frequency spectrum is the second characteristic amplitude. Therefore, the influence of sampling precision and errors on gear meshing frequency determination can be reduced, and the accuracy of the gear meshing frequency determined by the gear rotating speed determination method is improved. Since the error tends to be small, the frequency range is not too large, and in some embodiments, the difference between the upper limit value of the frequency range and the characteristic frequency included in the frequency range is [1,3], and the difference between the characteristic frequency included in the frequency range and the lower limit value of the frequency range is [1,3]. In some embodiments, the frequency range may be positive or negative 2 of the characteristic frequency, i.e., the frequency range is [ characteristic frequency-2, characteristic frequency +2], and the maximum amplitude value within this range may be used as the characteristic amplitude value.
In other embodiments, the determining the respective multiple frequency multiplication amplitudes in step S33 may specifically be determining, for each gear frequency, that the maximum frequency in the multiple frequency multiplication frequency ranges is the corresponding frequency multiplication amplitude, where the obtained first characteristic amplitude and second characteristic amplitude have eliminated errors, and step S34 may also be omitted adaptively.
In step S35, a product of the first characteristic amplitude and the second characteristic amplitude is determined, so as to obtain spectral energy corresponding to the gear frequency. Fig. 4 is a schematic diagram of spectral energy of the embodiment shown in fig. 1, please refer to fig. 4, in which the first characteristic amplitude and the second characteristic amplitude are obtained and then multiplied to obtain spectral energy for determining a rotation speed.
In some embodiments, according to the gear engagement frequency range and the value of n in step S33, the value range of the spectrum frequency can be determined, please continue to refer to fig. 2, in this embodiment, n=4, and in addition, since the collected vibration signal is distorted at a low rotation speed, the degree of gear failure cannot be reflected, and the generated power of the wind turbine generator is very low or even zero at this time, the gear engagement frequency range can be set to be [100,800] hz, the value range of the abscissa frequency of the spectrum is [400,3200] hz, and the value range of the spectrum energy is determined according to the gear engagement frequency range, please refer to fig. 4, and the value range of the abscissa of the spectrum energy is [100,800] hz, so that the detection data amount can be reduced, thereby reducing the calculation complexity.
Since the failure frequency is not the true gear mesh frequency, its multiplied frequency corresponds to a smaller amplitude and is closer to 0, so its corresponding spectral energy is smaller, while the correct gear mesh frequency, its multiplied frequency corresponds to the largest spectral energy. In step S4, it is determined that the gear frequency corresponding to the maximum value in the spectrum energy is the gear mesh frequency. In some embodiments, the first characteristic frequency is a multiple of the first characteristic frequency, and the gear engagement frequency is the first characteristic frequency.
In step S5, according to the gear engagement frequency and the number of teeth of the gear to be measured, an average gear rotational speed of the gear to be measured in the time to be measured is determined. The formula isWherein speed is the average gear speed, f mes h is the gear mesh frequency obtained in step S4, and f teeth is the number of teeth.
According to the gear rotating speed determining method, the frequency spectrum energy of the gear is obtained by processing the gear vibration signal, the frequency spectrum energy is used as a basis to obtain the gear meshing frequency, and then the gear rotating speed is obtained by calculating according to the gear meshing frequency, so that the principle that the frequency spectrum energy of the gear can accurately reflect the state of the gear is skillfully utilized, and the obtained gear rotating speed has higher accuracy no matter when the gear normally operates or fails, and is simple and convenient to calculate and easy to use. FIG. 5 is a schematic diagram of the error of the embodiment of FIG. 1. Referring to FIG. 5, 150 gearbox vibration signals are verified using the method, the error being defined as The maximum error is only 0.0128.
The embodiment of the application also provides a gear rotating speed determining system 200, which can be applied to a wind driven generator. Fig. 6 is a schematic block diagram of a gear speed determination system 200 according to one embodiment of the application. As shown in fig. 6, the gear speed determination system 200 may include one or more processors 201 for implementing the gear speed determination method described in any of the embodiments above. In some embodiments, gear speed determination system 200 may include a computer readable storage medium 202, computer readable storage medium 202 may store a program that may be invoked by processor 201, and may include a non-volatile storage medium. In some embodiments, gear speed determination system 200 may include memory 203 and interface 204. In some embodiments, the gear speed determination system 200 of embodiments of the present application may also include other hardware depending on the application.
The gear speed determination system 200 according to the embodiment of the present application has similar advantages as the gear speed determination method described above, and therefore will not be described herein.
The embodiment of the application also provides a computer readable storage medium. A computer-readable storage medium has stored thereon a program which, when executed by a processor, implements the gear rotational speed determination method described in any of the above embodiments.
Embodiments of the application may take the form of a computer program product embodied on one or more storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having program code embodied therein. Computer-readable storage media include both non-transitory and non-transitory, removable and non-removable media, and information storage may be implemented by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer-readable storage media include, but are not limited to, new types of memory, such as phase change memory/resistive memory/magnetic memory/ferroelectric memory (PRAM/RRAM/MRAM/FeRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, read only memory (CD-ROM), digital Versatile Disks (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device.
The method for determining the rotational speed of the gear and the system thereof and the computer readable storage medium provided by the embodiment of the application are described in detail above. Specific examples are set forth herein to illustrate the method and system for determining gear rotational speed and the computer readable storage medium according to embodiments of the present application, and the description of the above embodiments is merely for aiding in understanding the core concept of the present application and is not intended to limit the present application. It should be noted that it will be apparent to those skilled in the art that various changes and modifications can be made herein without departing from the spirit and principles of the application, which should also fall within the scope of the appended claims.
Claims (9)
1. A gear rotational speed determination method, characterized by comprising:
collecting gear vibration signals of a gear to be tested in the time to be tested;
Performing Fourier transform on the gear vibration signal to obtain a frequency spectrum, wherein the frequency spectrum comprises gear frequency and corresponding amplitude;
According to the frequency spectrum, determining the frequency spectrum energy corresponding to the gear frequency; determining a first characteristic frequency and a second characteristic frequency corresponding to each of a plurality of gear frequencies in the frequency spectrum, wherein the first characteristic frequency and the second characteristic frequency are integer multiples of the corresponding gear frequencies and are not equal to each other; determining a first characteristic amplitude according to the frequency spectrum and the first characteristic frequency, and determining a second characteristic amplitude according to the frequency spectrum and the second characteristic frequency; determining the product of the first characteristic amplitude and the second characteristic amplitude to obtain the frequency spectrum energy corresponding to the gear frequency;
Determining the gear frequency corresponding to the maximum value in the frequency spectrum energy as gear meshing frequency;
and determining the average gear rotating speed of the gear to be tested in the time to be tested according to the gear meshing frequency and the tooth number of the gear to be tested.
2. The gear rotational speed determination method as set forth in claim 1, wherein the determining the spectral energy corresponding to the gear frequency from the frequency spectrum includes:
and obtaining the product of the amplitude corresponding to each gear frequency and the amplitude corresponding to the frequency doubling of the gear frequency according to the frequency spectrum, and determining the frequency spectrum energy corresponding to the gear frequency.
3. The gear rotational speed determination method as defined in claim 1, wherein said determining a first characteristic amplitude from said frequency spectrum and said first characteristic frequency comprises:
determining the maximum value of the amplitude corresponding to the frequency range including the first characteristic frequency in the frequency spectrum as the first characteristic amplitude, or determining the amplitude corresponding to the first characteristic frequency as the first characteristic amplitude, and/or
Said determining a second characteristic amplitude from said spectrum and said second characteristic frequency comprises:
And determining the maximum value of the amplitude corresponding to the frequency range including the second characteristic frequency in the frequency spectrum as the second characteristic amplitude, or determining the amplitude corresponding to the second characteristic frequency as the second characteristic amplitude.
4. The method of determining a rotational speed of a gear according to claim 3, wherein the range of values of the difference between the upper limit value of the frequency range and the characteristic frequency included in the frequency range is [1,3], the range of values of the difference between the characteristic frequency included in the frequency range and the lower limit value of the frequency range is [1,3], and the characteristic frequency is the first characteristic frequency or the second characteristic frequency.
5. The gear rotational speed determination method as defined in claim 1, wherein the determining the first and second characteristic frequencies corresponding to each of the plurality of gear frequencies in the frequency spectrum comprises:
determining a plurality of frequency multiplication of each gear frequency and frequency multiplication amplitudes corresponding to the plurality of frequency multiplication respectively, wherein the plurality of frequency multiplication is 1 to n frequency multiplication, and n is a positive integer greater than or equal to 2;
For each gear frequency, determining the frequency multiplication of the gear frequency corresponding to the maximum value in a plurality of frequency multiplication amplitudes as the first characteristic frequency; and determining the frequency multiplication of the gear frequency corresponding to a second maximum value in the multiple frequency multiplication amplitudes as the second characteristic frequency.
6. The gear rotational speed determination method as defined in claim 1, wherein the determining the first and second characteristic frequencies corresponding to each of the plurality of gear frequencies in the frequency spectrum comprises:
determining a gear meshing frequency range according to the rotating speed range of the gear and the number of teeth of the gear;
A plurality of the gear frequencies in the gear mesh frequency range is determined.
7. The gear rotational speed determination method as defined in claim 6, wherein said determining a plurality of said gear frequencies in said gear mesh frequency range comprises:
determining a start frequency of the gear engagement frequency range, and a plurality of frequencies spaced apart by a fixed frequency from the start frequency of the gear engagement frequency range as a start, as a plurality of the gear frequencies, or
And determining a plurality of frequencies corresponding to the amplitudes larger than an amplitude threshold value in the gear meshing frequency range on the frequency spectrum as a plurality of gear frequencies.
8. A gear speed determination system comprising one or more processors configured to implement the gear speed determination method of any one of claims 1-7.
9. A computer-readable storage medium, characterized in that a program is stored thereon, which when executed by a processor, implements the gear rotational speed determination method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210583340.6A CN114993670B (en) | 2022-05-25 | 2022-05-25 | Gear speed determination method and system and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210583340.6A CN114993670B (en) | 2022-05-25 | 2022-05-25 | Gear speed determination method and system and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114993670A CN114993670A (en) | 2022-09-02 |
| CN114993670B true CN114993670B (en) | 2025-01-24 |
Family
ID=83028637
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210583340.6A Active CN114993670B (en) | 2022-05-25 | 2022-05-25 | Gear speed determination method and system and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114993670B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110987404A (en) * | 2019-12-23 | 2020-04-10 | 安徽容知日新科技股份有限公司 | Rotating speed identification method and computing equipment |
| CN113029553A (en) * | 2019-12-24 | 2021-06-25 | 国家能源投资集团有限责任公司 | Method, system and device for extracting rotating speed information of gearbox shaft and storage medium |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103336140A (en) * | 2013-06-20 | 2013-10-02 | 苏州大学 | Rotational speed estimation method and device based on gear meshing vibration |
| CN104077474A (en) * | 2014-06-23 | 2014-10-01 | 华南理工大学 | Meshing frequency and spectrum correction technology based wind power gear box order tracking method |
| CN104865400B (en) * | 2015-04-14 | 2018-08-28 | 华北电力大学 | A kind of detection recognition method and system of Wind turbines rotating speed |
| CN105548595B (en) * | 2015-12-18 | 2019-04-30 | 国网河北能源技术服务有限公司 | A detection method for extracting the rotational speed of each shaft of a wind turbine gearbox |
| DE102017200761A1 (en) * | 2017-01-18 | 2018-07-19 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | DEVICE FOR MONITORING A ROTATING ELEMENT AND CORRESPONDING METHOD |
| CN109116048A (en) * | 2018-08-29 | 2019-01-01 | 西安科技大学 | A kind of rotating speed of automobile engine detection and control method |
| CN109268214B (en) * | 2018-10-29 | 2020-04-17 | 国电联合动力技术有限公司 | Intelligent monitoring method for centering state of wind driven generator coupler |
| CN111256993A (en) * | 2018-11-30 | 2020-06-09 | 中国电力科学研究院有限公司 | Method and system for diagnosing fault type of main bearing of wind turbine generator |
| CN112345238B (en) * | 2020-10-29 | 2023-07-25 | 上海电气风电集团股份有限公司 | Gearbox vibration monitoring method and system and computer readable storage medium thereof |
| CN114018573B (en) * | 2021-11-16 | 2024-11-08 | 西安因联信息科技有限公司 | A parameter-free method for estimating gearbox gear meshing frequency |
-
2022
- 2022-05-25 CN CN202210583340.6A patent/CN114993670B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110987404A (en) * | 2019-12-23 | 2020-04-10 | 安徽容知日新科技股份有限公司 | Rotating speed identification method and computing equipment |
| CN113029553A (en) * | 2019-12-24 | 2021-06-25 | 国家能源投资集团有限责任公司 | Method, system and device for extracting rotating speed information of gearbox shaft and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114993670A (en) | 2022-09-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102288286A (en) | Method for analyzing and evaluating measure point precision of gearbox in vibration acceleration sensor | |
| WO2015196735A1 (en) | Wind power gear box order tracking method based on meshing frequency and spectrum correction technology | |
| CN108362492B (en) | A vibration separation method suitable for fault diagnosis of planetary gear train at low speed | |
| CN115031959A (en) | Gear fault diagnosis method, system and computer readable storage medium | |
| CN103411774A (en) | On-line early warning method of wind turbine generating unit on fluctuation working condition | |
| CN106970264B (en) | An Improved Phase Difference Correction Method Considering Grid Frequency Change Rate | |
| CN110598269B (en) | Discrete spectrum parameter correction method in low sampling point | |
| Verma et al. | Modeling and prediction of gearbox faults with data-mining algorithms | |
| Peng et al. | Speed estimation in planetary gearboxes: A method for reducing impulsive noise | |
| CN113804354B (en) | Holographic synchronous dynamic balance method and system for flexible rotor | |
| CN108317052B (en) | Method and device for detecting damage factor of gears, wind turbine | |
| CN105604805B (en) | Method and device for measuring yaw angle of wind generating set | |
| CN110219816A (en) | Method and system for Fault Diagnosis of Fan | |
| CN102945318B (en) | A kind of ultra-short term wind speed dynamic prediction method based on cascade blower fan | |
| Sun et al. | Current-aided dynamic time warping for planetary gearbox fault detection at time-varying speeds | |
| CN114993670B (en) | Gear speed determination method and system and computer readable storage medium | |
| CN102901596B (en) | Method for testing photoelectric reflection type dynamic torque of equal-diameter rotary shaft | |
| CN109883692B (en) | Generalized differential filtering method based on built-in encoder information | |
| CN112345250B (en) | Bearing vibration monitoring method and system and computer readable storage medium | |
| WO2023206860A1 (en) | Method and apparatus for determining mechanical device fault | |
| CN113029553B (en) | Gearbox shaft speed information extraction method, system, device and storage medium | |
| CN112345238B (en) | Gearbox vibration monitoring method and system and computer readable storage medium thereof | |
| Houjoh et al. | Utilization of synchronous averaging for inspection of tooth surface undulations on gears (localization of nonmesh harmonic components to individual gear) | |
| CN108572074B (en) | Bearing fault detection method and device and wind generating set | |
| Huang et al. | Fault detection of servo motor bearing based on speed signal under variable-speed conditions |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |