CN118190397B - Gear abrasion rapid measurement method for industrial cooling fan - Google Patents
Gear abrasion rapid measurement method for industrial cooling fan Download PDFInfo
- Publication number
- CN118190397B CN118190397B CN202410621416.9A CN202410621416A CN118190397B CN 118190397 B CN118190397 B CN 118190397B CN 202410621416 A CN202410621416 A CN 202410621416A CN 118190397 B CN118190397 B CN 118190397B
- Authority
- CN
- China
- Prior art keywords
- parameter
- data
- period
- difference
- analyzed
- 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
- 238000001816 cooling Methods 0.000 title claims abstract description 33
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 238000005299 abrasion Methods 0.000 title claims abstract 19
- 238000005070 sampling Methods 0.000 claims abstract description 46
- 230000002159 abnormal effect Effects 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims description 31
- 230000002596 correlated effect Effects 0.000 claims description 16
- 230000000875 corresponding effect Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 10
- 238000012545 processing Methods 0.000 abstract description 3
- 230000005856 abnormality Effects 0.000 description 22
- 238000004458 analytical method Methods 0.000 description 11
- 238000013507 mapping Methods 0.000 description 7
- 230000007423 decrease Effects 0.000 description 5
- 230000011218 segmentation Effects 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000314 lubricant Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000000007 visual 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/021—Gearings
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2123/00—Data types
- G06F2123/02—Data types in the time domain, e.g. time-series data
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及风机运行数据处理技术领域,具体涉及一种面向工业冷却风机的齿轮磨损快速测量方法。The invention relates to the technical field of fan operation data processing, and in particular to a gear wear rapid measurement method for an industrial cooling fan.
背景技术Background Art
工业冷却风机齿轮通常是指用于传动工业冷却系统中风机的齿轮。在风机运行使用过程中,风机的齿轮可能由于使用时间的增加导致磨损,而磨损的齿轮则会对风机运行的稳定性等带来一定的影响,因此需要对风机齿轮进行磨损检测,当其磨损程度较大时,则需要及时维护或者更换,以避免其磨损程度较大,对风机的正常工作带来影响。Industrial cooling fan gears usually refer to the gears used to drive the fan in the industrial cooling system. During the operation of the fan, the fan gears may wear out due to the increase in usage time, and the worn gears will have a certain impact on the stability of the fan operation. Therefore, it is necessary to detect the wear of the fan gears. When the degree of wear is large, it is necessary to maintain or replace it in time to avoid the degree of wear and tear affecting the normal operation of the fan.
现有齿轮磨损检测一般是通过视觉或者采集齿轮相关振动数据进行检测,而利用齿轮振动相关数据在进行齿轮磨损检测时,由于风机的运转速率则与实际冷却需求有关,即齿轮在不同时间段的转速存在差异,采集的振动数据也会发生变化,利用振动数据对齿轮磨损检测时,难以统一影响因素的变量,导致齿轮磨损检测结果不准确。Existing gear wear detection is generally carried out through visual means or by collecting gear-related vibration data. When using gear vibration-related data for gear wear detection, since the operating speed of the fan is related to the actual cooling demand, that is, the speed of the gear varies in different time periods, the collected vibration data will also change. When using vibration data for gear wear detection, it is difficult to unify the variables of the influencing factors, resulting in inaccurate gear wear detection results.
发明内容Summary of the invention
为了解决齿轮转速影响振动数据,导致齿轮磨损检测结果不准确的技术问题,本发明的目的在于提供一种面向工业冷却风机的齿轮磨损快速测量方法,所采用的技术方案具体如下:In order to solve the technical problem that the gear speed affects the vibration data and causes inaccurate gear wear detection results, the purpose of the present invention is to provide a rapid measurement method for gear wear of industrial cooling fans. The technical solution adopted is as follows:
一种面向工业冷却风机的齿轮磨损快速测量方法,所述方法包括:A method for quickly measuring gear wear of an industrial cooling fan, the method comprising:
以预设采样周期获取冷却风机齿轮的振动数据和转速数据;将首个采样周期作为基准周期,其他采样周期作为待分析数据周期;The vibration data and speed data of the cooling fan gear are obtained at a preset sampling period; the first sampling period is used as a reference period, and the other sampling periods are used as data periods to be analyzed;
根据每个待分析数据周期与所述基准周期之间所述振动数据的差异特征,结合每个待分析数据周期与所述基准周期之间所述转速数据的差异特征,获取每个待分析数据周期的振动异常参数;根据每个待分析数据周期与所述基准周期之间的采样时间间隔,结合对应的所述振动异常参数,获取每个待分析数据周期的齿轮磨损速率;According to the difference characteristics of the vibration data between each data cycle to be analyzed and the reference cycle, combined with the difference characteristics of the speed data between each data cycle to be analyzed and the reference cycle, the vibration abnormality parameter of each data cycle to be analyzed is obtained; according to the sampling time interval between each data cycle to be analyzed and the reference cycle, combined with the corresponding vibration abnormality parameter, the gear wear rate of each data cycle to be analyzed is obtained;
获取所述齿轮磨损速率的时序序列;在所述时序序列中,以每个所述齿轮磨损速率为中心,根据预设邻域参数获取局部数据段;根据每个局部数据段内所述齿轮磨损速率的变化特征,结合所述时序序列中所述齿轮磨损速率的变化特征,获得齿轮磨损程度。A time series of the gear wear rate is obtained; in the time series, a local data segment is obtained based on preset neighborhood parameters with each gear wear rate as the center; and the gear wear degree is obtained based on the change characteristics of the gear wear rate in each local data segment and in combination with the change characteristics of the gear wear rate in the time series.
进一步地,所述振动异常参数的获取方法包括:Furthermore, the method for obtaining the abnormal vibration parameters includes:
根据每个待分析数据周期与所述基准周期之间的所述振动数据中幅值的差异特征,获取每个待分析数据周期的第一差异参数;Acquire a first difference parameter for each period of data to be analyzed according to a difference characteristic of the amplitude in the vibration data between each period of data to be analyzed and the reference period;
根据每个待分析数据周期与所述基准周期之间的所述振动数据中频率的差异特征,获取每个待分析数据周期的第二差异参数;Acquire a second difference parameter for each period of data to be analyzed according to a difference characteristic of the frequency in the vibration data between each period of data to be analyzed and the reference period;
根据每个待分析数据周期与所述基准周期之间的所述振动数据拟合的数据曲线的差异特征,获取每个待分析数据周期的第三差异参数;Acquire a third difference parameter for each data period to be analyzed according to difference characteristics of the data curve fitted by the vibration data between each data period to be analyzed and the reference period;
根据每个待分析数据周期与所述基准周期之间的所述转速数据的差异特征,获取每个待分析数据周期的第四差异参数;Acquire a fourth difference parameter of each data period to be analyzed according to the difference characteristics of the rotation speed data between each data period to be analyzed and the reference period;
根据每个待分析数据周期与所述基准周期之间的所述第一差异参数、所述第二差异参数、所述第三差异参数和所述第四差异参数,获取每个待分析数据周期的所述振动异常参数;所述第一差异参数、所述第二差异参数和所述第三差异参数均与所述振动异常参数正相关;所述第四差异参数与所述振动异常参数负相关。The vibration abnormality parameter of each data period to be analyzed is obtained according to the first difference parameter, the second difference parameter, the third difference parameter and the fourth difference parameter between each data period to be analyzed and the reference period; the first difference parameter, the second difference parameter and the third difference parameter are all positively correlated with the vibration abnormality parameter; the fourth difference parameter is negatively correlated with the vibration abnormality parameter.
进一步地,所述第二差异参数的获取方法包括:Furthermore, the method for obtaining the second difference parameter includes:
根据每个待分析数据周期与所述基准周期之间的所述振动数据中极值点数量的差异特征,获取每个待分析数据周期的第二差异参数。According to the difference characteristics of the number of extreme value points in the vibration data between each data period to be analyzed and the reference period, a second difference parameter of each data period to be analyzed is obtained.
进一步地,所述齿轮磨损程度的获取方法包括:Furthermore, the method for obtaining the gear wear degree includes:
根据每个局部数据段内所述齿轮磨损速率的变化特征,获得参考数据段;According to the variation characteristics of the gear wear rate in each local data segment, a reference data segment is obtained;
根据所述参考数据段与其他局部数据段内所述齿轮磨损速率的差异特征,获得第一磨损参数;Obtaining a first wear parameter according to the difference characteristics of the gear wear rate in the reference data segment and other local data segments;
根据所述时序序列中所述齿轮磨损速率的统计特征,获得第二磨损参数;Obtaining a second wear parameter according to the statistical characteristics of the gear wear rate in the time series;
以所述参考数据段的中心对所述时序序列进行分割,根据分割后的左右两侧时序序列中所述齿轮磨损速率的差异特征,获取第三磨损参数;The time series is segmented at the center of the reference data segment, and a third wear parameter is obtained according to the difference characteristics of the gear wear rates in the time series on the left and right sides after the segmentation;
根据所述第一磨损参数、所述第二磨损参数和所述第三磨损参数,获取齿轮磨损程度;所述第一磨损参数、所述第二磨损参数和所述第三磨损参数均与所述齿轮磨损程度正相关。The gear wear degree is obtained according to the first wear parameter, the second wear parameter and the third wear parameter; the first wear parameter, the second wear parameter and the third wear parameter are all positively correlated with the gear wear degree.
进一步地,所述参考数据段的获取方法包括:Furthermore, the method for obtaining the reference data segment includes:
获取每个局部数据段内所述齿轮磨损速率的方差作为第一稳定参数;Obtaining the variance of the gear wear rate in each local data segment as a first stable parameter;
获取每个局部数据段内所述齿轮磨损速率的均值与所述时序序列的均值的差值绝对值,作为第二稳定参数;Obtaining an absolute value of a difference between a mean value of the gear wear rate in each local data segment and a mean value of the time series as a second stable parameter;
根据所述第一稳定参数和所述第二稳定参数,获取每个局部数据段的稳定程度参数;将所述稳定程度参数最大的局部数据段作为参考数据段;所述第一稳定参数和所述第二稳定参数均与所述稳定程度参数负相关。According to the first stability parameter and the second stability parameter, a stability parameter of each local data segment is obtained; the local data segment with the largest stability parameter is used as a reference data segment; the first stability parameter and the second stability parameter are both negatively correlated with the stability parameter.
进一步地,所述第一磨损参数的获取方法包括:Furthermore, the method for obtaining the first wear parameter includes:
将所述参考数据段内所述齿轮磨损速率的均值与其他局部数据段内所述齿轮磨损速率的均值的差值绝对值作为运行波动参数,将所有所述运行波动参数的平均值作为所述第一磨损参数。The absolute value of the difference between the mean value of the gear wear rate in the reference data segment and the mean value of the gear wear rate in other local data segments is taken as an operation fluctuation parameter, and the average value of all the operation fluctuation parameters is taken as the first wear parameter.
进一步地,所述第三磨损参数的获取方法包括:Furthermore, the method for obtaining the third wear parameter includes:
将分割后的右侧时序序列中所述齿轮磨损速率的均值与左侧时序序列中所述齿轮磨损速率的均值的比值作为第三磨损参数;其中右侧为采样时间最大的一侧。The ratio of the mean value of the gear wear rate in the right time series after segmentation to the mean value of the gear wear rate in the left time series is used as the third wear parameter; wherein the right side is the side with the largest sampling time.
进一步地,所述齿轮磨损程度的计算方法包括:Furthermore, the method for calculating the gear wear degree includes:
将所述第一磨损参数、所述第二磨损参数和所述第三磨损参数的乘积进行归一化后,归一化结果作为所述齿轮磨损程度。After the product of the first wear parameter, the second wear parameter and the third wear parameter is normalized, the normalized result is used as the gear wear degree.
进一步地,所述第三差异参数获取时使用的算法为DTW算法。Furthermore, the algorithm used when acquiring the third difference parameter is a DTW algorithm.
进一步地,所述齿轮磨损速率的获取方法包括:Furthermore, the method for obtaining the gear wear rate includes:
将每个待分析数据周期与所述基准周期的之间的采样时间间隔作为分母,将对应的所述振动异常参数作为分子,将分式作为每个待分析数据周期的齿轮磨损速率。The sampling time interval between each data cycle to be analyzed and the reference cycle is used as the denominator, the corresponding vibration abnormality parameter is used as the numerator, and the fraction is used as the gear wear rate of each data cycle to be analyzed.
本发明具有如下有益效果:The present invention has the following beneficial effects:
本发明首先以预设采样周期获取冷却风机齿轮的振动数据和转速数据,为后续分析提供分析基础;进一步将风机齿轮磨损程度最小的首个采样周期作为基准周期,其他采样周期作为待分析数据周期,更有利于后续了解齿轮磨损情况;进一步根据每个待分析数据周期与基准周期之间的振动数据和转速数据的差异特征,获取每个待分析数据周期的振动异常参数,消除齿轮转速对振动数据的影响,提升最终检测结果的准确性;进一步根据每个待分析数据周期与基准周期之间的采样时间间隔,结合对应的振动异常参数,获取每个待分析数据周期的齿轮磨损速率,为后续分析齿轮磨损的变化情况提供依据;进一步获取齿轮磨损速率的时序序列并获取局部数据段,便于后续从整体和局部两个角度分析齿轮磨损的变化情况;最后根据每个局部数据段内齿轮磨损速率的变化特征,结合时序序列中齿轮磨损速率的变化特征,获得齿轮磨损程度。本发明通过分析其他采样周期与基准周期内振动数据的差异,结合转速数据的差异,消除齿轮转速对振动数据的影响;然后根据齿轮磨损速率构成的序列的变化特征获取齿轮磨损程度,提升了检测方法的准确性。The present invention first obtains the vibration data and speed data of the cooling fan gear in a preset sampling period, so as to provide an analysis basis for subsequent analysis; further, the first sampling period in which the fan gear wear is the smallest is used as the reference period, and other sampling periods are used as the data period to be analyzed, which is more conducive to the subsequent understanding of the gear wear condition; further, according to the difference characteristics of the vibration data and speed data between each data period to be analyzed and the reference period, the vibration abnormality parameters of each data period to be analyzed are obtained, the influence of the gear speed on the vibration data is eliminated, and the accuracy of the final detection result is improved; further, according to the sampling time interval between each data period to be analyzed and the reference period, combined with the corresponding vibration abnormality parameters, the gear wear rate of each data period to be analyzed is obtained, so as to provide a basis for the subsequent analysis of the change of gear wear; further, the time series sequence of the gear wear rate is obtained and the local data segment is obtained, so as to facilitate the subsequent analysis of the change of gear wear from both the overall and local perspectives; finally, according to the change characteristics of the gear wear rate in each local data segment, combined with the change characteristics of the gear wear rate in the time series sequence, the gear wear degree is obtained. The present invention eliminates the influence of gear speed on vibration data by analyzing the difference between vibration data in other sampling periods and the reference period, combined with the difference in speed data; then, the gear wear degree is obtained according to the changing characteristics of the sequence composed of gear wear rate, thereby improving the accuracy of the detection method.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.
图1为本发明一个实施例所提供的一种面向工业冷却风机的齿轮磨损快速测量方法的流程图;FIG1 is a flow chart of a method for rapid measurement of gear wear for industrial cooling fans provided by an embodiment of the present invention;
图2为本发明一个实施例所提供的一种振动异常参数的获取方法的流程图;FIG2 is a flow chart of a method for obtaining abnormal vibration parameters provided by an embodiment of the present invention;
图3为本发明一个实施例所提供的一种齿轮磨损程度的获取方法的流程图。FIG. 3 is a flow chart of a method for obtaining the degree of gear wear provided by an embodiment of the present invention.
具体实施方式DETAILED DESCRIPTION
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种面向工业冷却风机的齿轮磨损快速测量方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。In order to further explain the technical means and effects adopted by the present invention to achieve the predetermined invention purpose, the following is a detailed description of the specific implementation method, structure, features and effects of a gear wear rapid measurement method for industrial cooling fans proposed by the present invention in combination with the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" does not necessarily refer to the same embodiment. In addition, specific features, structures or characteristics in one or more embodiments may be combined in any suitable form.
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
下面结合附图具体的说明本发明所提供的一种面向工业冷却风机的齿轮磨损快速测量方法的具体方案。The specific scheme of the method for rapid measurement of gear wear for industrial cooling fans provided by the present invention is described in detail below with reference to the accompanying drawings.
请参阅图1,其示出了本发明一个实施例提供的一种面向工业冷却风机的齿轮磨损快速测量方法的流程图,具体包括:Please refer to FIG1 , which shows a flow chart of a method for rapid measurement of gear wear for industrial cooling fans provided by an embodiment of the present invention, which specifically includes:
步骤S1:以预设采样周期获取冷却风机齿轮的振动数据和转速数据;将首个采样周期作为基准周期,其他采样周期作为待分析数据周期。Step S1: Acquire the vibration data and speed data of the cooling fan gear in a preset sampling period; take the first sampling period as a reference period, and the other sampling periods as data periods to be analyzed.
在本发明实施例中,考虑到不同时间的制冷需求不同,工业冷却风机的工作模式不同,所以以预设采样周期获取监测数据;又为了分析齿轮转速与振动数据之间的关系,所以以预设采样周期获取冷却风机齿轮的振动数据和转速数据,为后续分析提供分析基础。考虑到风机齿轮在使用过程中,随着使用时间的增加,其磨损程度也会增加,为了了解齿轮磨损变化的情况,将首个采样周期作为基准周期,其他采样周期作为待分析数据周期,在后续分析其他采样周期与基准周期中数据的差异。In the embodiment of the present invention, considering the different refrigeration demands at different times and the different working modes of the industrial cooling fan, the monitoring data is obtained in a preset sampling period; and in order to analyze the relationship between the gear speed and the vibration data, the vibration data and speed data of the cooling fan gear are obtained in a preset sampling period to provide an analysis basis for subsequent analysis. Considering that the degree of wear of the fan gear increases with the increase of the use time during use, in order to understand the change of gear wear, the first sampling period is used as the reference period, and the other sampling periods are used as the data period to be analyzed, and the difference between the data in other sampling periods and the reference period is analyzed later.
需要说明的是,在本发明一个实施例中,预设采样周期为1分钟,采样频率为1秒/次;在本发明其他实施例中,实施者可以自行调整采样周期和采样频率。It should be noted that, in one embodiment of the present invention, the preset sampling period is 1 minute and the sampling frequency is 1 second/time; in other embodiments of the present invention, implementers can adjust the sampling period and sampling frequency by themselves.
步骤S2:根据每个待分析数据周期与基准周期之间振动数据的差异特征,结合每个待分析数据周期与基准周期之间转速数据的差异特征,获取每个待分析数据周期的振动异常参数;根据每个待分析数据周期与基准周期之间的采样时间间隔,结合对应的振动异常参数,获取每个待分析数据周期的齿轮磨损速率。Step S2: According to the difference characteristics of the vibration data between each data cycle to be analyzed and the reference cycle, combined with the difference characteristics of the speed data between each data cycle to be analyzed and the reference cycle, the vibration abnormality parameters of each data cycle to be analyzed are obtained; according to the sampling time interval between each data cycle to be analyzed and the reference cycle, combined with the corresponding vibration abnormality parameters, the gear wear rate of each data cycle to be analyzed is obtained.
在本发明实施例中,考虑到待分析数据周期与基准周期之间的振动数据差异越大,说明振动差异越大,待分析数据周期表现出的异常特征越明显;而待分析数据周期与基准周期之间的转速数据差异越大,说明风机的运行模式差异越大,振动数据表现出的异常特征的可信度越低,所以根据每个待分析数据周期与基准周期之间振动数据的差异特征,结合每个待分析数据周期与基准周期之间转速数据的差异特征,获取每个待分析数据周期的振动异常参数,消除风机不同运行模式的不同转速对振动数据造成的影响,为后续分析齿轮磨损特征提供依据。In the embodiment of the present invention, it is considered that the greater the difference in vibration data between the data cycle to be analyzed and the reference cycle, the greater the vibration difference, and the more obvious the abnormal characteristics exhibited by the data cycle to be analyzed; and the greater the difference in speed data between the data cycle to be analyzed and the reference cycle, the greater the difference in the operating mode of the fan, and the lower the credibility of the abnormal characteristics exhibited by the vibration data. Therefore, according to the difference characteristics of vibration data between each data cycle to be analyzed and the reference cycle, combined with the difference characteristics of speed data between each data cycle to be analyzed and the reference cycle, the vibration abnormality parameters of each data cycle to be analyzed are obtained, and the influence of different speeds of different operating modes of the fan on the vibration data is eliminated, thereby providing a basis for subsequent analysis of gear wear characteristics.
优选地,在本发明一个实施例中,振动异常参数的获取方法包括:Preferably, in one embodiment of the present invention, the method for obtaining abnormal vibration parameters includes:
请参阅图2,其示出了本发明一个实施例所提供的一种振动异常参数的获取方法的流程图:Please refer to FIG. 2 , which shows a flow chart of a method for obtaining abnormal vibration parameters provided by an embodiment of the present invention:
步骤S101:根据每个待分析数据周期与基准周期之间的振动数据中幅值的差异特征,获取每个待分析数据周期的第一差异参数。Step S101: acquiring a first difference parameter of each data period to be analyzed according to a difference characteristic of the amplitude in the vibration data between each data period to be analyzed and a reference period.
考虑到振动数据的幅值是在振动分析时的重要指标,幅值可以反映出风机的运行状态和健康状况,所以根据每个待分析数据周期与基准周期之间的振动数据中幅值的差异特征,获取每个待分析数据周期的第一差异参数,第一差异参数与振动异常参数正相关。Considering that the amplitude of vibration data is an important indicator in vibration analysis, the amplitude can reflect the operating status and health of the fan. Therefore, according to the difference characteristics of the amplitude in the vibration data between each data period to be analyzed and the reference period, the first difference parameter of each data period to be analyzed is obtained, and the first difference parameter is positively correlated with the vibration abnormality parameter.
在本发明一个实施例中,将待分析数据周期内振动数据的平均幅值与基准周期内振动数据的平均幅值的差值绝对值作为第一差异参数。在本发明其他实施例中,也可以通过幅值的波动范围差异、幅值变化速率差异等其他幅值差异角度,获取第一差异参数。In one embodiment of the present invention, the absolute value of the difference between the average amplitude of the vibration data in the data period to be analyzed and the average amplitude of the vibration data in the reference period is used as the first difference parameter. In other embodiments of the present invention, the first difference parameter can also be obtained through other amplitude difference angles such as the difference in amplitude fluctuation range and the difference in amplitude change rate.
步骤S102:根据每个待分析数据周期与基准周期之间的振动数据中频率的差异特征,获取每个待分析数据周期的第二差异参数。Step S102: acquiring a second difference parameter of each data period to be analyzed according to a difference characteristic of the frequency in the vibration data between each data period to be analyzed and a reference period.
考虑到振动数据的频率也是振动分析的重要指标,当机械部件出现故障或磨损时,会导致振动信号中的频率特征发生变化,所以根据每个待分析数据周期与基准周期之间的振动数据中频率的差异特征,获取每个待分析数据周期的第二差异参数,第二差异参数与振动异常参数正相关。Considering that the frequency of vibration data is also an important indicator of vibration analysis, when mechanical parts fail or wear out, the frequency characteristics in the vibration signal will change. Therefore, according to the difference characteristics of the frequency in the vibration data between each data period to be analyzed and the reference period, the second difference parameter of each data period to be analyzed is obtained, and the second difference parameter is positively correlated with the vibration abnormality parameter.
在本发明一个实施例中,考虑到振动信号的频率越高,极值点数量越多,极值点数量从侧面反映出了振动数据的频率特征,所以根据每个待分析数据周期与基准周期之间的振动数据中极值点数量的差异特征,将待分析数据周期与基准周期之间的振动数据中极值点数量的差值绝对值,作为第二差异参数。In one embodiment of the present invention, considering that the higher the frequency of the vibration signal, the more extreme points there are, and the number of extreme points indirectly reflects the frequency characteristics of the vibration data, the absolute value of the difference in the number of extreme points in the vibration data between each data period to be analyzed and the reference period is used as the second difference parameter based on the difference characteristics of the number of extreme points in the vibration data between the data period to be analyzed and the reference period.
在本发明其他实施例中,实施者也可以借助频谱图获取其他现有技术,提取出待分析数据周期中的振动数据的主要频率以及基准周期中振动数据的主要频率,分析主要频率的差异特征,获取第二差异参数;或者分析频率的范围差异、分布差异获取第二差异参数。In other embodiments of the present invention, the implementer may also obtain other existing technologies with the help of the spectrum diagram, extract the main frequency of the vibration data in the data period to be analyzed and the main frequency of the vibration data in the reference period, analyze the difference characteristics of the main frequencies, and obtain the second difference parameter; or analyze the range difference and distribution difference of the frequency to obtain the second difference parameter.
步骤S103:根据每个待分析数据周期与基准周期之间的振动数据拟合的数据曲线的差异特征,获取每个待分析数据周期的第三差异参数。Step S103: acquiring a third difference parameter of each data period to be analyzed according to the difference characteristics of the data curve fitted by the vibration data between each data period to be analyzed and the reference period.
考虑到对振动数据进行拟合曲线,数据曲线的差异特征同样能够反映出待分析数据周期中振动数据的异常特征,所以根据每个待分析数据周期与基准周期之间的振动数据拟合的数据曲线的差异特征,获取每个待分析数据周期的第三差异参数,第三差异参数与振动异常参数正相关。Considering the fitting curve of the vibration data, the difference characteristics of the data curve can also reflect the abnormal characteristics of the vibration data in the data period to be analyzed. Therefore, the third difference parameter of each data period to be analyzed is obtained based on the difference characteristics of the data curve fitted by the vibration data between each data period to be analyzed and the reference period. The third difference parameter is positively correlated with the vibration abnormality parameter.
在本发明一个实施例中,利用DTW(Dynamic Time Warping)算法获取待分析数据周期与基准周期之间的振动数据拟合的数据曲线的DTW距离,DTW距离越大,说明两条数据曲线之间的相似度越低,差异程度越大,所以将DTW距离作为第三差异参数。In one embodiment of the present invention, a DTW (Dynamic Time Warping) algorithm is used to obtain the DTW distance of a data curve fitted by vibration data between a data period to be analyzed and a reference period. The larger the DTW distance, the lower the similarity between the two data curves and the greater the degree of difference. Therefore, the DTW distance is used as the third difference parameter.
在本发明其他实施例中,实施者也可以通过现有技术中的余弦相似度、皮尔逊相关系数等算法,间接获取第三差异参数,其均为本领域技术人员所熟知的技术手段,在此不再进行赘述。In other embodiments of the present invention, the implementer may also indirectly obtain the third difference parameter through algorithms such as cosine similarity and Pearson correlation coefficient in the prior art, which are technical means well known to those skilled in the art and will not be described in detail here.
步骤S104:根据每个待分析数据周期与基准周期之间的转速数据的差异特征,获取每个待分析数据周期的第四差异参数。Step S104: acquiring a fourth difference parameter of each data period to be analyzed according to the difference characteristics of the rotation speed data between each data period to be analyzed and the reference period.
考虑到转动速度代表了冷却风机的制冷需求,也反映了冷却风机的工作模式,而不同的转动速度,必然会引起振动特征的变化,所以根据每个待分析数据周期与基准周期之间的转速数据的差异特征,获取每个待分析数据周期的第四差异参数,第四差异参数越大,说明越有可能是工作模式引起的振动数据变化,借助第四差异参数进行修正,第四差异参数与振动异常参数负相关。Taking into account that the rotation speed represents the cooling demand of the cooling fan and also reflects the working mode of the cooling fan, and different rotation speeds will inevitably cause changes in the vibration characteristics, the fourth difference parameter of each data period to be analyzed is obtained according to the difference characteristics of the rotation speed data between each data period to be analyzed and the reference period. The larger the fourth difference parameter is, the more likely it is that the vibration data change is caused by the working mode. The fourth difference parameter is used for correction. The fourth difference parameter is negatively correlated with the abnormal vibration parameter.
在本发明一个实施例中,将待分析数据周期内转速数据的均值与基准周期内转速数据的均值的差值绝对值作为第四差异参数。在本发明其他实施例中,考虑到转速变化也会引起振动数据波动,所以还可以从转速变化角度,根据转速变化率和转速的分布范围获取第四差异参数。In one embodiment of the present invention, the absolute value of the difference between the mean of the rotational speed data in the data period to be analyzed and the mean of the rotational speed data in the reference period is used as the fourth difference parameter. In other embodiments of the present invention, considering that the rotational speed change will also cause the vibration data to fluctuate, the fourth difference parameter can also be obtained from the perspective of the rotational speed change, according to the rotational speed change rate and the distribution range of the rotational speed.
步骤S105:根据每个待分析数据周期与基准周期之间的第一差异参数、第二差异参数、第三差异参数和第四差异参数,获取每个待分析数据周期的振动异常参数。Step S105: Obtain the vibration abnormality parameter of each data period to be analyzed according to the first difference parameter, the second difference parameter, the third difference parameter and the fourth difference parameter between each data period to be analyzed and the reference period.
获得待分析数据周期与基准周期之间的差异参数后,就可以结合多个角度的差异特征,获取待分析数据周期的振动异常参数。After obtaining the difference parameters between the data period to be analyzed and the reference period, the vibration abnormality parameters of the data period to be analyzed can be obtained by combining the difference features of multiple angles.
在本发明一个实施例中,将第一差异参数、第二差异参数和第三差异参数的和值利用指数函数映射,将映射结果与第四差异参数的比值作为振动异常参数。振动异常参数的计算公式包括:In one embodiment of the present invention, the sum of the first difference parameter, the second difference parameter and the third difference parameter is mapped using an exponential function, and the ratio of the mapping result to the fourth difference parameter is used as the vibration abnormality parameter. The calculation formula of the vibration abnormality parameter includes:
; ;
其中,表示待分析数据周期的序号;表示第个待分析数据周期的振动异常参数;表示以自然常数为底数的指数函数;表示第个待分析数据周期的第一差异参数;表示第个待分析数据周期的第二差异参数;表示第个待分析数据周期的第三差异参数;表示第个待分析数据周期的第四差异参数,为除零参数,在本发明一个实施例中,。in, Indicates the serial number of the data period to be analyzed; Indicates The vibration anomaly parameters of the data period to be analyzed; Indicated by natural constant The exponential function with base ; Indicates The first difference parameter of the data period to be analyzed; Indicates A second difference parameter of the data period to be analyzed; Indicates The third difference parameter of the data period to be analyzed; Indicates The fourth difference parameter of the data period to be analyzed, is a division by zero parameter. In one embodiment of the present invention, .
需要说明的是,在本发明其他实施例中,也可选用其他基础数学运算或者函数映射实现相关映射,其均为本领域技术人员熟知的技术手段,在此不做赘述。It should be noted that in other embodiments of the present invention, other basic mathematical operations or function mappings may also be used to implement relevant mappings, which are technical means well known to those skilled in the art and will not be elaborated herein.
其中,正相关关系表示因变量会随着自变量的增大而增大,因变量会随着自变量的减小而减小,具体关系可以为相乘关系、相加关系、指数函数的幂等,由实际应用进行确定;负相关关系表示因变量会随着自变量的增大而减小,因变量会随着自变量的减小而增大,可以为相减关系、相除关系等,由实际应用进行确定。Among them, a positive correlation means that the dependent variable will increase as the independent variable increases, and the dependent variable will decrease as the independent variable decreases. The specific relationship can be a multiplication relationship, an addition relationship, or the power of an exponential function, which is determined by actual application; a negative correlation means that the dependent variable will decrease as the independent variable increases, and the dependent variable will increase as the independent variable decreases. It can be a subtraction relationship, a division relationship, etc., which is determined by actual application.
考虑到风机齿轮在使用过程中,随着使用时间的增加,其磨损程度也会增加,为了分析齿轮磨损的变化情况,根据每个待分析数据周期与基准周期之间的采样时间间隔,结合对应的振动异常参数,获取每个待分析数据周期的齿轮磨损速率,对齿轮磨损的变化快慢进行度量,同时便于在后续步骤中分析齿轮磨损速率的变化,分析齿轮磨损程度。Taking into account that the degree of wear of fan gears will increase with the increase of usage time during use, in order to analyze the changes in gear wear, the gear wear rate of each data cycle to be analyzed is obtained according to the sampling time interval between each data cycle to be analyzed and the reference cycle, combined with the corresponding vibration abnormality parameters, and the speed of change of gear wear is measured. At the same time, it is convenient to analyze the changes in gear wear rate and the degree of gear wear in subsequent steps.
优选地,在本发明一个实施例中,将每个待分析数据周期与基准周期之间的采样时间间隔作为分母,将对应的振动异常参数作为分子,将分式作为每个待分析数据周期的齿轮磨损速率。Preferably, in one embodiment of the present invention, the sampling time interval between each data cycle to be analyzed and the reference cycle is used as the denominator, the corresponding vibration abnormality parameter is used as the numerator, and the fraction is used as the gear wear rate of each data cycle to be analyzed.
需要说明的是,将不同采样周期的首个采样时刻之间时间间隔作为待分析数据周期与基准周期之间的采样时间间隔。It should be noted that the time interval between the first sampling moments of different sampling periods is used as the sampling time interval between the period of the data to be analyzed and the reference period.
步骤S3:获取齿轮磨损速率的时序序列;在时序序列中,以每个齿轮磨损速率为中心,根据预设邻域参数获取局部数据段;根据每个局部数据段内齿轮磨损速率的变化特征,结合时序序列中齿轮磨损速率的变化特征,获得齿轮磨损程度。Step S3: Obtain a time series of gear wear rates; in the time series, take each gear wear rate as the center and obtain a local data segment according to preset neighborhood parameters; obtain the gear wear degree according to the change characteristics of the gear wear rate in each local data segment and the change characteristics of the gear wear rate in the time series.
由经验可知,实际使用过程中的风机齿轮等机械设备,其磨损速度不会保持不变,而是随着机械系统的运行,材料部件出现老化、疲劳损伤不断累积,同时齿轮润滑剂、齿轮之间的啮合度也在不断降低,齿轮磨损速率会发生变化,同时风机不同的运行模式产生的损耗速度也不同,所以齿轮磨损速率的变化特征能够反映出齿轮的磨损情况,因此获取齿轮磨损速率的时序序列。Experience shows that the wear rate of mechanical equipment such as fan gears will not remain constant during actual use. Instead, as the mechanical system operates, aging and fatigue damage of material components continue to accumulate. At the same time, the meshing degree between gear lubricants and gears continues to decrease, and the gear wear rate will change. At the same time, different operating modes of the fan produce different loss rates. Therefore, the changing characteristics of the gear wear rate can reflect the wear condition of the gear, so the time series of the gear wear rate is obtained.
需要说明的是,在本发明实施例中,时序序列是由齿轮磨损速率构成,时序序列中的齿轮磨损速率根据所属的采样周期的首个采样时间,按照从小到大的顺序进行排序。It should be noted that, in the embodiment of the present invention, the time series is composed of gear wear rates, and the gear wear rates in the time series are sorted in ascending order according to the first sampling time of the sampling period to which they belong.
在本发明实施例中,考虑到风机的运行过程中,短时间内的运行模式的切换的频率越高,越容易对齿轮造成磨损,所以在时序序列中,以每个齿轮磨损速率为中心,根据预设邻域参数获取局部数据段,便于后续在分析局部数据段中分析。In the embodiment of the present invention, considering that during the operation of the fan, the higher the frequency of switching the operating mode in a short period of time, the more likely it is to cause wear on the gears, so in the time series sequence, with the wear rate of each gear as the center, local data segments are obtained according to preset neighborhood parameters, which is convenient for subsequent analysis in the analysis of the local data segments.
优选地,在本发明一个实施例中,预设邻域参数为5,即局部数据段的长度为5,在本发明其他实施例中,实施者可以设置其他的预设邻域参数;对于时序序列首尾处,可以利用多项式拟合等其他现有技术中的预测或拟合方法,补充数据值,便于获取对应的局部数据段。Preferably, in one embodiment of the present invention, the preset neighborhood parameter is 5, that is, the length of the local data segment is 5. In other embodiments of the present invention, the implementer may set other preset neighborhood parameters; for the beginning and end of the time series sequence, other prediction or fitting methods in the prior art such as polynomial fitting may be used to supplement the data values to facilitate the acquisition of the corresponding local data segment.
获取齿轮磨损速率的时序序列并获取局部数据段后,就根据每个局部数据段内齿轮磨损速率的变化特征,结合时序序列中齿轮磨损速率的变化特征,获得齿轮磨损程度。After obtaining the time series of the gear wear rate and the local data segments, the gear wear degree is obtained according to the change characteristics of the gear wear rate in each local data segment combined with the change characteristics of the gear wear rate in the time series.
优选地,在本发明一个实施例中,齿轮磨损程度的获取方法包括:Preferably, in one embodiment of the present invention, the method for obtaining the gear wear degree includes:
请参阅图3,其示出了本发明一个实施例所提供的一种齿轮磨损程度的获取方法的流程图:Please refer to FIG3 , which shows a flow chart of a method for obtaining the degree of gear wear provided by an embodiment of the present invention:
步骤S201:根据每个局部数据段内齿轮磨损速率的变化特征,获得参考数据段。Step S201: Obtain a reference data segment according to the variation characteristics of the gear wear rate in each local data segment.
考虑到所有的局部数据段中,存在相对平稳的局部数据段,将其作为参考数据段,便于对比出其他局部数据段的波动特征。Considering that there are relatively stable local data segments among all the local data segments, they are used as reference data segments to facilitate comparison of the fluctuation characteristics of other local data segments.
在本发明一个实施例中,考虑到方差越小,说明数据波动越稳定,所以获取每个局部数据段内齿轮磨损速率的方差作为第一稳定参数,第一稳定参数与稳定程度参数负相关;In one embodiment of the present invention, considering that the smaller the variance is, the more stable the data fluctuation is, the variance of the gear wear rate in each local data segment is obtained as the first stability parameter, and the first stability parameter is negatively correlated with the stability parameter;
考虑到局部数据段内齿轮磨损速率的均值与时序序列的均值的差异越小,说明局部数据段越能代表整体序列特征,所以获取每个局部数据段内齿轮磨损速率的均值与时序序列的均值的差值绝对值,作为第二稳定参数,第二稳定参数与稳定程度参数负相关;Considering that the smaller the difference between the mean of the gear wear rate in the local data segment and the mean of the time series is, the more the local data segment can represent the overall sequence characteristics, so the absolute value of the difference between the mean of the gear wear rate in each local data segment and the mean of the time series is obtained as the second stability parameter, and the second stability parameter is negatively correlated with the stability degree parameter;
根据第一稳定参数和第二稳定参数,获取每个局部数据段的稳定程度参数;将稳定程度参数最大的局部数据段作为参考数据段。According to the first stability parameter and the second stability parameter, a stability parameter of each local data segment is obtained; and the local data segment with the largest stability parameter is used as a reference data segment.
在本发明一个实施例中,将第一稳定参数和第二稳定参数的乘积负相关映射并归一化,将归一化结果作为稳定程度参数;稳定程度参数的计算公式包括:In one embodiment of the present invention, the product of the first stability parameter and the second stability parameter is negatively correlated and normalized, and the normalized result is used as the stability parameter; the calculation formula of the stability parameter includes:
; ;
其中,表示局部数据段的序号;表示第个局部数据段的稳定程度参数;表示以自然常数为底数的指数函数;表示第个局部数据段的第一稳定参数;表示第个局部数据段的第二稳定参数。in, Indicates the sequence number of the local data segment; Indicates The stability parameter of a local data segment; Indicated by natural constant The exponential function with base ; Indicates The first stable parameter of a local data segment; Indicates The second stable parameter of a local data segment.
需要说明的是,在本发明其他实施例中,也可以获取平均绝对偏差、变异系数等其他数据特征分析数据的稳定性,也可选用其他基础数学运算或者函数映射实现相关映射,其均为本领域技术人员熟知的技术手段,在此不做赘述。It should be noted that in other embodiments of the present invention, other data characteristics such as mean absolute deviation and coefficient of variation can also be obtained to analyze the stability of the data, and other basic mathematical operations or function mappings can also be used to achieve relevant mappings. These are technical means well known to those skilled in the art and will not be elaborated here.
需要说明的是,当存在多个稳定程度参数最大时,考虑到运行时间越短,风机齿轮的磨损程度越低,参考性更强,选取稳定程度参数最大的所有局部时间段中,采样时间最靠前的局部时间段作为参考时间段。It should be noted that when there are multiple maximum stability parameters, considering that the shorter the running time, the lower the degree of wear of the fan gear and the stronger the reference, the local time period with the earliest sampling time among all local time periods with the largest stability parameters is selected as the reference time period.
步骤S202:根据参考数据段与其他局部数据段内齿轮磨损速率的差异特征,获得第一磨损参数。Step S202: Obtain a first wear parameter according to the difference characteristics of the gear wear rates in the reference data segment and other local data segments.
考虑到参考数据段与其他局部数据段内齿轮磨损速率的差异越大,说明磨损速率的一致性越差,反映出齿轮运行模式发生变化的频率越高,齿轮磨损程度就越大,第一磨损参数与齿轮磨损程度正相关。Considering that the greater the difference between the gear wear rate in the reference data segment and other local data segments, the worse the consistency of the wear rate, it reflects that the higher the frequency of changes in the gear operation mode, the greater the gear wear degree, and the first wear parameter is positively correlated with the gear wear degree.
在本发明一个实施例中,考虑到局部数据段内齿轮磨损速率的均值代表了数据段的整体特征,所以将参考数据段内齿轮磨损速率的均值与其他局部数据段内齿轮磨损速率的均值的差值绝对值作为运行波动参数,将所有运行波动参数的平均值作为第一磨损参数。In one embodiment of the present invention, considering that the mean of the gear wear rate in a local data segment represents the overall characteristics of the data segment, the absolute value of the difference between the mean of the gear wear rate in the reference data segment and the mean of the gear wear rate in other local data segments is taken as an operating fluctuation parameter, and the average of all operating fluctuation parameters is taken as the first wear parameter.
在本发明其他实施例中,实施者也可以根据稳定程度参数的差异,将其他局部数据段的稳定程度参数与参考数据段的稳定程度参数的差值求平均后,作为第一磨损参数。In other embodiments of the present invention, the implementer may also average the difference between the stability parameters of other local data segments and the stability parameters of the reference data segment according to the difference in the stability parameters, and use the average as the first wear parameter.
步骤S203:根据时序序列中齿轮磨损速率的统计特征,获得第二磨损参数。Step S203: Obtain a second wear parameter according to the statistical characteristics of the gear wear rate in the time series.
考虑到整体时序序列的齿轮磨损速率越大,齿轮的磨损程度越大,所以根据时序序列中齿轮磨损速率的统计特征,获得第二磨损参数,第二磨损参数与齿轮磨损程度正相关。Considering that the greater the gear wear rate of the overall time series, the greater the gear wear degree, the second wear parameter is obtained according to the statistical characteristics of the gear wear rate in the time series, and the second wear parameter is positively correlated with the gear wear degree.
在本发明一个实施例中,将时序序列的所有齿轮磨损速率的均值作为第二磨损参数。在本发明其他实施例中,实施者也可以基于中程数、众数等获取第二磨损参数,例如将所有齿轮磨损速率的均值、众数和中程数的平均值作为第二磨损参数。In one embodiment of the present invention, the mean of the wear rates of all gears in the time series is used as the second wear parameter. In other embodiments of the present invention, the implementer may also obtain the second wear parameter based on the mid-range number, mode, etc., for example, the mean, mode, and average of the mid-range number of the wear rates of all gears are used as the second wear parameter.
步骤S204:以参考数据段的中心对时序序列进行分割,根据分割后的左右两侧时序序列中齿轮磨损速率的差异特征,获取第三磨损参数。Step S204: segmenting the time series at the center of the reference data segment, and obtaining a third wear parameter according to the difference characteristics of the gear wear rates in the time series on the left and right sides after segmentation.
考虑到齿轮随着工作时间的增加,相同运行条件下的磨损速率不断增大,时序序列整体呈现增大的趋势,并且增大的速度越快,说明齿轮的磨损程度越大;所以以参考数据段的中心对时序序列进行分割,根据分割后的左右两侧时序序列中齿轮磨损速率的差异特征,获取第三磨损参数。Considering that the wear rate of gears increases with the increase of working time under the same operating conditions, the time series shows an increasing trend as a whole, and the faster the increase, the greater the degree of gear wear; therefore, the time series is segmented at the center of the reference data segment, and the third wear parameter is obtained based on the difference characteristics of the gear wear rates in the time series on the left and right sides after segmentation.
在本发明一个实施例中,将分割后的右侧时序序列中齿轮磨损速率的均值与左侧时序序列中齿轮磨损速率的均值的比值作为第三磨损参数;其中右侧为采样时间最大的一侧。In one embodiment of the present invention, the ratio of the mean gear wear rate in the right time series after segmentation to the mean gear wear rate in the left time series is used as the third wear parameter; wherein the right side is the side with the largest sampling time.
右侧时序序列中齿轮磨损速率整体相较于左侧时序序列越大,说明磨损速率增大的越快,越符合实际风机运行中齿轮出现磨损的特征,第三磨损参数越大,齿轮磨损程度越大。The overall gear wear rate in the right timing sequence is larger than that in the left timing sequence, which means that the wear rate increases faster, which is more consistent with the characteristics of gear wear in actual wind turbine operation. The larger the third wear parameter, the greater the degree of gear wear.
在本发明其他实施例中,实施者也可以根据时序序列的中位数,将时序序列两等分,获取第三磨损参数;其中时序序列为奇数时,中位数处的齿轮磨损速率不参与计算。In other embodiments of the present invention, the implementer may also divide the time sequence into two equal parts according to the median of the time sequence to obtain the third wear parameter; wherein when the time sequence is an odd number, the gear wear rate at the median is not included in the calculation.
步骤S205:根据第一磨损参数、第二磨损参数和第三磨损参数,获取齿轮磨损程度。Step S205: Obtain the gear wear degree according to the first wear parameter, the second wear parameter and the third wear parameter.
在本发明一个实施例中,将第一磨损参数、第二磨损参数和第三磨损参数的乘积进行归一化后,归一化结果作为齿轮磨损程度。In one embodiment of the present invention, after the product of the first wear parameter, the second wear parameter and the third wear parameter is normalized, the normalized result is used as the gear wear degree.
在本发明其他实施例中,也可选用其他基础数学运算或者函数映射实现相关映射,其均为本领域技术人员熟知的技术手段,在此不做赘述。In other embodiments of the present invention, other basic mathematical operations or function mappings may be used to implement relevant mappings, which are technical means well known to those skilled in the art and will not be elaborated herein.
在本发明其他实施例中,实施者还可以结合工业冷却风机的运行时间分析齿轮磨损程度,由于运行时间越长,对齿轮的损耗越大,所运行时间与齿轮磨损程度正相关。In other embodiments of the present invention, the implementer may also analyze the degree of gear wear in combination with the operating time of the industrial cooling fan. Since the longer the operating time, the greater the loss to the gear, the operating time is positively correlated with the degree of gear wear.
需要说明的是,本发明实施例也适用于实时检测,只需要将实时获取的采样数据作为新的待分析数据周期中的数据,然后执行步骤S2,更新时序序列,然后执行步骤S3,即可输出实时的齿轮磨损程度。It should be noted that the embodiment of the present invention is also applicable to real-time detection. It only needs to use the sampling data acquired in real time as the data in the new data cycle to be analyzed, and then execute step S2, update the timing sequence, and then execute step S3 to output the real-time gear wear degree.
需要说明的是,在本发明一个实施例中,当齿轮磨损程度大于0.3时,认为冷却风机的齿轮工作状态异常,需要及时维护或更换。It should be noted that, in one embodiment of the present invention, when the gear wear degree is greater than 0.3, it is considered that the gear of the cooling fan is in an abnormal working state and needs to be maintained or replaced in time.
综上所述,本发明针对齿轮转速影响振动数据,导致齿轮磨损检测结果不准确的技术问题,提出了一种面向工业冷却风机的齿轮磨损快速测量方法。本发明首先以预设采样周期获取冷却风机齿轮的振动数据和转速数据;将首个采样周期作为基准周期,其他采样周期作为待分析数据周期;进一步根据每个待分析数据周期与基准周期之间的振动数据和转速数据的差异特征,获取每个待分析数据周期的振动异常参数;进一步根据每个待分析数据周期与基准周期之间的采样时间间隔,结合对应的振动异常参数,获取每个待分析数据周期的齿轮磨损速率;进一步获取齿轮磨损速率的时序序列;进一步在时序序列中,以每个齿轮磨损速率为中心,根据预设邻域参数获取局部数据段;最后根据每个局部数据段内齿轮磨损速率的变化特征,结合时序序列中齿轮磨损速率的变化特征,获得齿轮磨损程度。本发明通过分析其他采样周期与基准周期内振动数据的差异,结合转速数据的差异,消除齿轮转速对振动数据的影响;然后根据齿轮磨损速率构成的序列的变化特征获取齿轮磨损程度,提升了检测方法的准确性。In summary, the present invention aims at the technical problem that the gear speed affects the vibration data and causes the gear wear detection result to be inaccurate, and proposes a fast gear wear measurement method for industrial cooling fans. The present invention first obtains the vibration data and speed data of the cooling fan gear with a preset sampling period; takes the first sampling period as the reference period, and the other sampling periods as the data period to be analyzed; further, according to the difference characteristics of the vibration data and speed data between each data period to be analyzed and the reference period, obtains the vibration abnormality parameters of each data period to be analyzed; further, according to the sampling time interval between each data period to be analyzed and the reference period, combined with the corresponding vibration abnormality parameters, obtains the gear wear rate of each data period to be analyzed; further obtains the time series of the gear wear rate; further in the time series, with each gear wear rate as the center, obtains the local data segment according to the preset neighborhood parameters; finally, according to the change characteristics of the gear wear rate in each local data segment, combined with the change characteristics of the gear wear rate in the time series, obtains the gear wear degree. The present invention eliminates the influence of gear speed on vibration data by analyzing the difference between vibration data in other sampling periods and the reference period, combined with the difference in speed data; then, the gear wear degree is obtained according to the changing characteristics of the sequence composed of gear wear rate, thereby improving the accuracy of the detection method.
需要说明的是:上述本发明实施例先后顺序仅仅为了描述,不代表实施例的优劣。在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。It should be noted that the sequence of the above embodiments of the present invention is for description only and does not represent the advantages and disadvantages of the embodiments. The processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from other embodiments.
Claims (7)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410621416.9A CN118190397B (en) | 2024-05-20 | 2024-05-20 | Gear abrasion rapid measurement method for industrial cooling fan |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410621416.9A CN118190397B (en) | 2024-05-20 | 2024-05-20 | Gear abrasion rapid measurement method for industrial cooling fan |
Publications (2)
Publication Number | Publication Date |
---|---|
CN118190397A CN118190397A (en) | 2024-06-14 |
CN118190397B true CN118190397B (en) | 2024-09-10 |
Family
ID=91400449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410621416.9A Active CN118190397B (en) | 2024-05-20 | 2024-05-20 | Gear abrasion rapid measurement method for industrial cooling fan |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118190397B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118375711B (en) * | 2024-06-20 | 2024-09-27 | 摩多利智能传动(江苏)有限公司 | Wear detection method and alarm device for planetary gear speed reducer |
CN118410450B (en) * | 2024-07-02 | 2024-09-06 | 艾力斯特健康科技有限公司 | Massage chair remote health monitoring method based on cloud computing |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112304597A (en) * | 2020-10-13 | 2021-02-02 | 潍柴动力股份有限公司 | Gear detection method and device and vehicle |
CN112621382A (en) * | 2020-12-17 | 2021-04-09 | 北京工业大学 | Online cutter wear monitoring method based on cutting vibration damping characteristics |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0663916B2 (en) * | 1989-03-14 | 1994-08-22 | 富士電機株式会社 | Gear abnormality diagnosis device |
DE4023663C2 (en) * | 1989-07-31 | 1996-02-01 | Siemens Ag | Method for diagnosing the mechanical properties of a machine that has rotating components |
JP4848803B2 (en) * | 2006-03-13 | 2011-12-28 | 日本精工株式会社 | Abnormality diagnosis device, rotation device, and abnormality diagnosis method |
EP2498076A1 (en) * | 2011-03-11 | 2012-09-12 | Hexagon Technology Center GmbH | Wear-Monitoring of a Gearbox in a Power Station |
CN208012648U (en) * | 2018-01-11 | 2018-10-26 | 珠海市鑫世达测控技术有限公司 | Rotating machinery monitors system |
CN113468721B (en) * | 2021-06-07 | 2024-03-29 | 太原科技大学 | Method for predicting residual life of gears and bearings in gear reduction box |
CN113468801B (en) * | 2021-06-07 | 2024-03-26 | 太原科技大学 | Gear nuclear density estimation residual life prediction method |
CN114295367B (en) * | 2021-11-12 | 2024-09-06 | 华能新能源股份有限公司 | Wind turbine generator system gear box working condition on-line monitoring method |
CN114674551B (en) * | 2022-03-08 | 2023-10-13 | 广州机械科学研究院有限公司 | Gear abrasion energy monitoring method, device and system and electronic equipment |
CN115238829B (en) * | 2022-09-21 | 2022-12-20 | 南通佳晟宏机械配件有限公司 | A method for analyzing roll wear degree of hot rolling mill |
CN116258075B (en) * | 2023-02-16 | 2024-04-23 | 重庆工业职业技术学院 | Fault learning method for speed reducer of industrial robot |
CN116378910A (en) * | 2023-03-24 | 2023-07-04 | 华能广西清洁能源有限公司 | Real-time monitoring and fault diagnosis system for gearbox of wind driven generator |
CN117708748B (en) * | 2024-02-05 | 2024-05-03 | 苏州众志新环冷却设备有限公司 | Operation monitoring system and method for centrifugal fan |
CN117875191B (en) * | 2024-03-08 | 2024-05-28 | 东莞市星火齿轮有限公司 | Gearbox service life assessment method and system based on big data and storage medium |
-
2024
- 2024-05-20 CN CN202410621416.9A patent/CN118190397B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112304597A (en) * | 2020-10-13 | 2021-02-02 | 潍柴动力股份有限公司 | Gear detection method and device and vehicle |
CN112621382A (en) * | 2020-12-17 | 2021-04-09 | 北京工业大学 | Online cutter wear monitoring method based on cutting vibration damping characteristics |
Also Published As
Publication number | Publication date |
---|---|
CN118190397A (en) | 2024-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN118190397B (en) | Gear abrasion rapid measurement method for industrial cooling fan | |
JP5460160B2 (en) | Equipment diagnostic equipment | |
RU2484442C1 (en) | Vibration diagnostics and forecasting method of sudden failure of engine, and carrier | |
WO2020136823A1 (en) | Abnormality diagnosis device and abnormality diagnosis method | |
JP2002371989A (en) | Monitoring and controlling method, and monitoring device for compressor | |
CN117390972B (en) | Electronic cigarette atomizer service life prediction method | |
CN105114334A (en) | Method for monitoring abrasion loss of impeller wear ring of multi-stage centrifugal pump based on computational fluid dynamics theory | |
CN118312884B (en) | Performance evaluation method and system for electro-absorption modulator | |
CN118535945B (en) | A method and system for monitoring and diagnosing the operating status of a circulation pump | |
CN118391249B (en) | Method and device for processing running state monitoring data of large water pump | |
CN116952584B (en) | Motor bearing performance degradation assessment method based on vibration signal amplitude spectrum similarity | |
CN117235678B (en) | Method for detecting blade abrasion degree in wood cutting process | |
CN110398384A (en) | The machine health monitoring of rotating machinery | |
CN118520704B (en) | Method, system, equipment and readable storage medium for optimizing comprehensive performance of horizontal decanter centrifuge | |
Liu et al. | Fault diagnosis of ball bearing elements: A generic procedure based on time-frequency analysis | |
CN117459073A (en) | Intelligent management method for heat pump system operation data | |
CN118734229B (en) | Abnormal identification method for motor operation data | |
CN115393276A (en) | Intelligent detection and analysis system for surface morphology of gear | |
CN115438688A (en) | Fault detection method and device for key equipment, electronic equipment and storage medium | |
CN117848488A (en) | Online processing method for equipment vibration data | |
CN117560015B (en) | Efficient processing system for air conditioner installation engineering data | |
JP7383367B1 (en) | Vibration data analysis method and analysis system for rotating equipment | |
CN110146267B (en) | Method for detecting flutter of valve plate of totally-enclosed refrigeration compressor | |
CN106126840A (en) | Rotating machinery rotating speed based on orbit of shaft center calculates device and method | |
CN119362960A (en) | Centralized control method for variable frequency starting and speed regulation system of compressed air energy storage compressor |
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 |