WO2020007375A1 - Continuous monitoring-based method and device for identifying 1p signal of wind turbine, terminal, and computer readable storage medium - Google Patents

Continuous monitoring-based method and device for identifying 1p signal of wind turbine, terminal, and computer readable storage medium Download PDF

Info

Publication number
WO2020007375A1
WO2020007375A1 PCT/CN2019/098094 CN2019098094W WO2020007375A1 WO 2020007375 A1 WO2020007375 A1 WO 2020007375A1 CN 2019098094 W CN2019098094 W CN 2019098094W WO 2020007375 A1 WO2020007375 A1 WO 2020007375A1
Authority
WO
WIPO (PCT)
Prior art keywords
frequency
fan
wind turbine
signal
continuous monitoring
Prior art date
Application number
PCT/CN2019/098094
Other languages
French (fr)
Chinese (zh)
Inventor
胡卫华
滕军
张笑
Original Assignee
哈尔滨工业大学(深圳)
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 哈尔滨工业大学(深圳) filed Critical 哈尔滨工业大学(深圳)
Publication of WO2020007375A1 publication Critical patent/WO2020007375A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones

Definitions

  • the present application belongs to the technical field of fan modal analysis, and specifically relates to a method, a device, a terminal, and a computer-readable storage medium for fan 1p signal recognition based on continuous monitoring.
  • this application provides a method, device, terminal, and computer-readable storage medium for fan 1p signal recognition based on continuous monitoring, which can effectively identify 1p signals and extract 1p frequency, providing effective for wind turbine structure monitoring. And reliable means of automation.
  • a method for identifying a 1p signal of a fan based on continuous monitoring includes:
  • the eigenmode function whose average instantaneous frequency is located in the low frequency section is taken as the effective eigenmode function, and the sum of the effective eigenmode functions is used as the target eigenmode function, the low frequency section is Frequency range below 1Hz;
  • the frequency corresponding to the target eigenmode function is identified and a steady state map is established, and a 1p frequency is extracted according to the steady state map.
  • the "identifying the frequency corresponding to the target eigenmode function and establishing a steady state map” includes:
  • a steady state diagram is established according to the frequency and modal shape of the fan.
  • the transform is a Hilbert transform.
  • the structural dynamic response of the wind turbine is an acceleration, speed or displacement response.
  • obtaining the structural dynamic response of the wind turbine includes:
  • a fan 1p signal identification device based on continuous monitoring includes:
  • An acquisition module configured to acquire a structural dynamic response of the wind turbine
  • a decomposition module configured to decompose the structural dynamic response of the wind turbine into several eigenmode functions
  • a transformation module configured to transform each eigenmode function to obtain a corresponding instantaneous frequency, and calculate an average value of the instantaneous frequency of each eigenmode function
  • a screening module configured to use the eigenmode function with the instantaneous frequency mean value located in a low frequency section as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function, so
  • the low frequency range is a frequency range below 1 Hz;
  • the identification module is configured to identify a frequency corresponding to the target eigenmode function and establish a steady state map, and extract a 1p frequency according to the steady state map.
  • the transform is a Hilbert transform.
  • the identification module includes:
  • a modeling sub-module configured to establish a discrete state space equation
  • An identification sub-module configured to calculate a frequency, a damping ratio, and a modal shape of the fan according to the discrete state space equation and a structural dynamic response of the fan;
  • a steady-state graph sub-module configured to establish a steady-state graph according to the frequency and modal shape of the fan
  • An extraction submodule configured to extract a 1p frequency according to the steady-state map.
  • a terminal includes a memory and a processor.
  • the memory is configured to store a computer program, and the processor executes the computer program to enable the terminal to implement the continuous monitoring-based fan 1p signal identification method according to any one of the above. .
  • a computer-readable storage medium stores the computer program executed by the terminal.
  • the corresponding Hilbert spectrum is obtained by decomposing and transforming the structural dynamic response of the wind turbine, and the eigenmode function corresponding to the highest concentration in the low frequency range is selected as the target eigenmode function.
  • the eigenmode function is used to identify the corresponding frequency and establish a steady-state map.
  • the 1p frequency is extracted according to the steady-state map to realize the automatic recognition of the 1p signal, which has significant effectiveness and recognition accuracy.
  • FIG. 1 is a flowchart of a method for identifying a 1p signal of a fan based on continuous monitoring according to Embodiment 1 of the present application;
  • FIG. 2 is a flowchart of step A of a method for identifying a fan 1p signal based on continuous monitoring according to Embodiment 1 of the present application;
  • FIG. 3 is a schematic diagram of an axonometric structure of sensors involved in a method for identifying a fan 1p signal based on continuous monitoring according to Embodiment 1 of the present application;
  • FIG. 4 is a time-domain distribution diagram of an acceleration response based on a continuous monitoring of a 1p signal recognition method of a fan provided in Embodiment 1 of the present application;
  • FIG. 5 is a time-domain distribution diagram of an eigenmode function based on a continuous monitoring of a fan 1p signal recognition method provided in Embodiment 1 of the present application;
  • step E is a flowchart of step E of a method for identifying a fan 1p signal based on continuous monitoring provided in Embodiment 1 of the present application;
  • FIG. 7 is a steady state diagram obtained by the method for identifying a 1p signal of a fan based on continuous monitoring provided in Embodiment 1 of the present application;
  • FIG. 9 is a Campbell diagram obtained by the method for identifying a 1p signal of a wind turbine based on continuous monitoring provided in Embodiment 1 of the present application;
  • FIG. 10 is a schematic structural diagram of a fan 1p signal recognition device based on continuous monitoring provided in Embodiment 2 of the present application;
  • FIG. 11 is a schematic structural diagram of an acquisition module of a fan 1p signal identification device based on continuous monitoring provided in Embodiment 2 of the present application;
  • FIG. 12 is a schematic structural diagram of an identification module of a fan 1p signal identification device based on continuous monitoring provided in Embodiment 2 of the present application;
  • FIG. 13 is a schematic structural diagram of a terminal provided in Embodiment 3 of the present application.
  • 100-continuous monitoring-based fan 1p signal identification device 110-acquisition module, 111-cycle submodule, 112-acquisition submodule, 120-decomposition module, 130-transformation module, 140-screening module, 150-identification module, 151 -Modeling sub-module, 152-recognition sub-module, 153-steady-state graph sub-module, 154-extraction sub-module, 200-terminal, 210-memory, 220-processor, 230-input unit, 240-display unit, 300 -Fan, 310-tower, 320-cabin, 400-sensor.
  • this embodiment provides a method for identifying a 1p signal of a fan based on continuous monitoring.
  • the method includes steps A to E:
  • Step A Obtain the structural dynamic response of the fan 300.
  • the structural dynamic response refers to the speed, acceleration, and displacement generated by the fan 300 under the dynamic load.
  • the structural dynamic response to be obtained may be an acceleration signal collected by an acceleration sensor, or a speed signal or a displacement signal measured by a displacement sensor.
  • step A may include steps A1 to A2:
  • Step A1 Determine the time constant period of the fan 300.
  • the time constant period is determined according to the rotation speed of the fan 300, the pitch angle and the turning angle of the nacelle 320.
  • the fan 300 needs to adjust the relevant parameters such as the rotation speed, the pitch angle, and the cabin 320 in real time according to the wind direction and wind speed to ensure better energy conversion efficiency, which belongs to a time-varying structure that is difficult to accurately measure.
  • a time-invariant cycle needs to be determined.
  • the fan 300 conforms to the time-invariant assumption, that is, its structural characteristics can be considered to not change with time.
  • Step A2 Obtain the structural dynamic response of the fan 300 in the time-invariant period.
  • the tower 310 of the wind turbine 300 is divided into three parts along the height, and the acquisition units are respectively arranged at four points of the three parts. Among them, the bottom end of the tower 310 is 4Q.
  • the sensor 400 may be an acceleration sensor or a displacement sensor. Additionally, FIG. 4 shows the acceleration response signal collected by the sensor 400.
  • Step B Decompose the structural dynamic response of the wind turbine 300 into several eigenmode functions.
  • the eigenmode function, Intrinsic Mode Function (imf) should satisfy the following conditions: (1) the number of extreme points of the function and the zero crossings are equal or at most one difference over the entire signal length; (2) at any time, The average value of the upper envelope defined by the maximum value of the function and the lower envelope defined by the minimum value is zero, that is, the upper and lower envelopes are symmetrical to the time axis.
  • FIG. 5 shows time-domain distribution diagrams (imf1 to 6) of each eigenmode function obtained by decomposition.
  • Step C Transform each eigenmode function to obtain a corresponding instantaneous frequency, and calculate the mean instantaneous frequency value of each eigenmode function.
  • the transform is a Hilbert transform. It should be understood that after each eigenmode function undergoes Hilbert transformation, a series of instantaneous frequencies will be obtained. The mean value of the instantaneous frequency of the eigenmode function is obtained by averaging the instantaneous frequencies corresponding to the same eigenmode function.
  • Step D Use the eigenmode function whose mean frequency is in the low frequency region as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function.
  • the band is a frequency range below 1 Hz.
  • the eigenmode functions that satisfy the following conditions are all valid eigenmode functions: the mean instantaneous frequency is located in the low frequency section. Then add all the effective eigenmode functions, and the sum is the target eigenmode function.
  • Step E Identify the frequency corresponding to the target eigenmode function and establish a steady-state map, and extract a 1p frequency according to the steady-state map.
  • the steady state diagram is used to indicate the position of the poles of the system. Since the poles are the global characteristics of the system (represented by the frequency of the wind turbine 300 here), as the model order increases, the system poles extracted by the mathematical model of the order increase will appear repeatedly and be characterized on the same graph in order to Find the physical poles of the structure by observing the pole distribution. Furthermore, by performing feature extraction on the steady-state graph, a frequency of 1p can be obtained. It should be understood that the 1p frequency is less than the fundamental frequency of the fan 300.
  • step E includes steps E1 to E4:
  • Step E1 Establish a discrete state space equation.
  • the discrete state space equation is established based on the discretized structural dynamic response, and a mathematical relationship is established between the structural dynamic response of the wind turbine 300 and the frequency, damping ratio, and modal shape.
  • the specific form of the discrete state space equation is:
  • x k is a discrete state vector
  • y k is a discrete output vector
  • w k and ⁇ k are white noise terms
  • is one of the acceleration, speed, and displacement of the fan 300
  • a ⁇ is a discrete state matrix
  • C ⁇ is a discrete output matrix.
  • Step E2 Calculate the frequency, damping ratio and modal shape of the fan 300 according to the discrete state space equation and the structural dynamic response of the fan 300.
  • Step E3 Establish a steady-state map according to the frequency and modal shapes of the fan 300.
  • FIG. 7 shows an actual steady state diagram obtained by performing step E3.
  • Step E4 Extract the 1p frequency according to the steady-state map.
  • the fundamental frequency of the fan 300 can be accurately extracted to be 0.6 Hz, and the 1p frequency is 0.41 Hz.
  • FIG. 8 shows a steady-state map obtained based on the conventional method.
  • the steady-state map lacks an obvious extreme point distribution in a low frequency region, and it is difficult to extract a 1p frequency.
  • the steady-state map ( Figure 7) obtained based on the continuous monitoring of the 1p signal recognition method of the fan provided in this embodiment has a significant and concentrated extreme point distribution in the low-frequency region, which provides reliable observation and extraction of the 1p frequency.
  • the basics The basics.
  • the 1p frequency is caused by slight quality differences between the leaves.
  • the slight mass difference will generate an excitation force with a fixed frequency to the wind tower, and the fixed frequency is an excitation frequency with a double speed.
  • the average speed of the fan 300 during the time-invariant period is calculated, and a required number of frequency-average speed numbers are obtained, and a Campbell diagram is established according to the frequency-average speed numbers.
  • a preset number of time-invariant periods are continuously monitored, and the frequency and average speed of the fan 300 corresponding to each time-invariant period are obtained, so as to obtain a frequency-average speed number pair.
  • a Campbell diagram based on continuous monitoring can be established to analyze the structural mode of the fan 300.
  • FIG. 9 shows a Campbell diagram obtained based on continuous monitoring of acceleration. According to the Campbell diagram, the 1p frequency (0.41 Hz) of the fan 300 can be obtained.
  • this embodiment provides a device 1p signal identification device 100 based on continuous monitoring.
  • the device includes:
  • An obtaining module 110 configured to obtain a structural dynamic response of the wind turbine 300
  • a decomposition module 120 configured to decompose the structural dynamic response of the wind turbine 300 into several eigenmode functions
  • the transformation module 130 is configured to transform each eigenmode function to obtain an instantaneous frequency corresponding thereto, and calculate an average instantaneous frequency value of each eigenmode function;
  • the screening module 140 is configured to use the eigenmode function with the instantaneous frequency mean value located in the low-frequency section as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function,
  • the low frequency section is a frequency range below 1 Hz;
  • the identification module 150 is configured to identify a frequency corresponding to the target eigenmode function and establish a steady-state map, and extract a 1p frequency according to the steady-state map.
  • the obtaining module 110 includes:
  • a period sub-module 111 configured to determine a time-invariant period of the fan 300
  • the obtaining sub-module 112 is configured to obtain a structural dynamic response of the wind turbine 300 in the time-invariant period.
  • the identification module 150 includes:
  • a modeling sub-module 151 configured to establish a discrete state space equation
  • An identification submodule 152 configured to calculate a frequency, a damping ratio, and a modal shape of the fan 300 according to the discrete state space equation and the structural dynamic response of the fan 300;
  • a steady-state map sub-module 153 configured to establish a steady-state map according to the frequency and modal shapes of the fan 300;
  • the extraction sub-module 154 is configured to extract a 1p frequency according to the steady-state map.
  • this embodiment provides a terminal 200.
  • the terminal 200 includes a memory 210 and a processor 220.
  • the memory 210 is configured to store a computer program, and the processor 220 executes the computer program to enable the terminal 200 to implement the continuous-based Identification method of the monitored 1p signal of the fan.
  • the terminal 200 includes terminal devices (such as computers and servers) that do not have mobile communication capabilities, and also includes mobile terminals (such as smart phones, tablet computers, on-board computers, and smart wearable devices).
  • terminal devices such as computers and servers
  • mobile terminals such as smart phones, tablet computers, on-board computers, and smart wearable devices.
  • the memory 210 may include a storage program area and a storage data area.
  • the storage program area can store an operating system and applications required for at least one function (such as a sound playback function and an image playback function), and the like; the storage data area can store data (such as audio data and Backup files, etc.).
  • the memory 210 may include a high-speed random access memory, and may further include a non-volatile memory (for example, at least one magnetic disk storage device, a flash memory device, or other volatile solid-state storage device).
  • the terminal 200 further includes an input unit 230 and a display unit 240.
  • the input unit 230 is configured to receive various instructions or parameters (including a preset scrolling method, a preset time interval, and a preset scrolling number) input by a user, including a mouse, a keyboard, a touch panel, and other input devices.
  • the display unit 240 is configured to display various output information (including a webpage page and a parameter configuration interface, etc.) of the terminal 200, including a display panel.
  • a computer-readable storage medium is also provided, which stores the computer program executed by the terminal 200.
  • each block in the flowchart or block diagram may represent a module, a program segment, or a part of code, which contains one or more components for implementing a specified logical function Executable instructions.
  • each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart can be a dedicated hardware-based system that performs a specified function or action To achieve, or can be implemented by a combination of dedicated hardware and computer instructions.
  • the functional modules or units in the various embodiments of the present application may be integrated together to form an independent part, or each of the modules may exist separately, or two or more modules may be integrated to form an independent part.
  • the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium.
  • the technical solution of the present application is essentially a part that contributes to the existing technology or a part of the technical solution can be embodied in the form of a software product.
  • the computer software product is stored in a storage medium, including Several instructions are used to cause a computer device (which may be a smart phone, a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in the embodiments of the present application.
  • the foregoing storage media include: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disks or optical disks and other media that can store program codes .
  • any specific value should be construed as exemplary only and not as a limitation, so other examples of the exemplary embodiments may have different values.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Wind Motors (AREA)
  • Control Of Positive-Displacement Air Blowers (AREA)

Abstract

A continuous monitoring-based method for identifying a 1p signal of a wind turbine, comprising: acquiring the structural dynamics response of a wind turbine (A); decomposing the structural dynamics response of the wind turbine into several intrinsic mode functions (B); transforming each intrinsic mode function so as to obtain the instantaneous frequency corresponding to said function, and calculating the instantaneous frequency average for each intrinsic mode function (C); treating intrinsic mode functions of which the instantaneous frequency average is in a low frequency segment as valid intrinsic mode functions, and using the sum of said valid intrinsic mode functions as the target intrinsic mode function, said low frequency segment being a frequency range below 1Hz (D); identifying the frequency corresponding to the target intrinsic mode function and establishing a stability graph, and extracting the 1p frequency on the basis of the stability graph (E). Also disclosed is a continuous monitoring-based device for identifying a 1p signal of a wind turbine, comprising an acquisition module (110), a decomposition module (120), a transformation module (130), a screening module (140), and an identification module (150). Also disclosed are a continuous monitoring-based terminal and computer readable storage medium for identifying a 1p signal of a wind turbine. The present5 invention can effectively identify a 1p signal and extract the 1p frequency, thereby providing effective and reliable automated means for wind turbine structure monitoring.

Description

基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质Method, device, terminal and computer-readable storage medium for fan 1p signal recognition based on continuous monitoring
相关申请的交叉引用Cross-reference to related applications
本申请要求于2018年07月06日提交中国专利局的申请号为CN201810738808.8、名称为“基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。This application requires a Chinese patent application filed with the Chinese Patent Office on July 06, 2018 under the application number CN201810738808.8 and titled "Continuous Monitoring-based Fan 1p Signal Identification Method, Device, Terminal, and Computer-readable Storage Medium" Priority, the entire contents of which are incorporated herein by reference.
技术领域Technical field
本申请属于风机模态分析技术领域,具体地来说,是一种基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质。The present application belongs to the technical field of fan modal analysis, and specifically relates to a method, a device, a terminal, and a computer-readable storage medium for fan 1p signal recognition based on continuous monitoring.
背景技术Background technique
随着人们环保意识的提高,清洁能源的需要日益增长,使风电行业迎来了良好的发展机遇。风机多建设于荒野无人区域,风速较大且环境恶劣,难以到达现场进行维护。尤其是海上风机,环境更为恶劣且更难到达,使测量维护面临巨大挑战。With the increase of people's awareness of environmental protection, the demand for clean energy is increasing, which brings good development opportunities to the wind power industry. The wind turbines are mostly built in unmanned areas, where the wind speed is high and the environment is harsh, making it difficult to reach the site for maintenance. Especially for offshore wind turbines, the environment is more severe and harder to reach, making measurement and maintenance face huge challenges.
为了实现风机的连续监测,一些在线监测系统应运而生。目前,在线监测系统多运用于监测风机的3p(三倍转速激励频率)以及6p(六倍转速激励频率)等信号,对于1p(一倍转速激励频率)信号难以有效地识别观察,具有明显的缺陷。In order to achieve continuous monitoring of wind turbines, some online monitoring systems have emerged. At present, the online monitoring system is mostly used to monitor the 3p (three-speed excitation frequency) and 6p (six-speed excitation frequency) and other signals. defect.
发明内容Summary of the invention
为了克服现有技术的不足,本申请提供了一种基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质,可有效识别1p信号而提取1p频率,为风机结构监测提供有效而可靠的自动化手段。In order to overcome the shortcomings of the prior art, this application provides a method, device, terminal, and computer-readable storage medium for fan 1p signal recognition based on continuous monitoring, which can effectively identify 1p signals and extract 1p frequency, providing effective for wind turbine structure monitoring. And reliable means of automation.
本申请的目的通过以下技术方案来实现:The purpose of this application is achieved by the following technical solutions:
一种基于连续监测的风机1p信号识别方法,包括:A method for identifying a 1p signal of a fan based on continuous monitoring includes:
获取风机的结构动力响应;Obtain the structural dynamic response of the fan;
将所述风机的结构动力响应分解为若干本征模态函数;Decomposing the structural dynamic response of the wind turbine into several eigenmode functions;
对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值;Transform each eigenmode function to obtain the instantaneous frequency corresponding to it, and calculate the mean instantaneous frequency value of each eigenmode function;
以所述瞬时频率均值位于低频区段的本征模态函数作为有效本征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围;The eigenmode function whose average instantaneous frequency is located in the low frequency section is taken as the effective eigenmode function, and the sum of the effective eigenmode functions is used as the target eigenmode function, the low frequency section is Frequency range below 1Hz;
识别所述目标本征模态函数对应的频率并建立稳态图,根据所述稳态图提取1p频率。The frequency corresponding to the target eigenmode function is identified and a steady state map is established, and a 1p frequency is extracted according to the steady state map.
作为上述技术方案的改进,所述“识别所述目标本征模态函数对应的频率并建立稳态图”包括:As an improvement of the above technical solution, the "identifying the frequency corresponding to the target eigenmode function and establishing a steady state map" includes:
建立离散状态空间方程;Establish discrete state space equations;
根据所述离散状态空间方程与所述风机的结构动力响应计算风机的频率、阻尼比与模态振型;Calculating the fan frequency, damping ratio and modal shape according to the discrete state space equation and the structural dynamic response of the fan;
根据所述风机的频率与模态振型建立稳态图。A steady state diagram is established according to the frequency and modal shape of the fan.
作为上述技术方案的进一步改进,所述变换为希尔伯特变换。As a further improvement of the above technical solution, the transform is a Hilbert transform.
作为上述技术方案的进一步改进,所述风机的结构动力响应为加速度、速度或位移响应。As a further improvement of the above technical solution, the structural dynamic response of the wind turbine is an acceleration, speed or displacement response.
作为上述技术方案的进一步改进,“获取风机的结构动力响应”包括:As a further improvement of the above technical solution, "obtaining the structural dynamic response of the wind turbine" includes:
确定风机的时不变周期;Determine the time constant period of the fan;
获取所述风机于所述时不变周期内的结构动力响应。Obtain the structural dynamic response of the wind turbine in the time-invariant period.
一种基于连续监测的风机1p信号识别装置,包括:A fan 1p signal identification device based on continuous monitoring includes:
获取模块,配置成获取风机的结构动力响应;An acquisition module configured to acquire a structural dynamic response of the wind turbine;
分解模块,配置成将所述风机的结构动力响应分解为若干本征模态函数;A decomposition module configured to decompose the structural dynamic response of the wind turbine into several eigenmode functions;
变换模块,配置成对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值;A transformation module configured to transform each eigenmode function to obtain a corresponding instantaneous frequency, and calculate an average value of the instantaneous frequency of each eigenmode function;
筛选模块,配置成以所述瞬时频率均值位于低频区段的本征模态函数作为有效本征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围;A screening module configured to use the eigenmode function with the instantaneous frequency mean value located in a low frequency section as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function, so The low frequency range is a frequency range below 1 Hz;
识别模块,配置成识别所述目标本征模态函数对应的频率并建立稳态图,根据稳态图提取1p频率。The identification module is configured to identify a frequency corresponding to the target eigenmode function and establish a steady state map, and extract a 1p frequency according to the steady state map.
作为上述技术方案的改进,所述变换为希尔伯特变换。As an improvement of the above technical solution, the transform is a Hilbert transform.
作为上述技术方案的进一步改进,所述识别模块包括:As a further improvement to the above technical solution, the identification module includes:
建模子模块,配置成建立离散状态空间方程;A modeling sub-module configured to establish a discrete state space equation;
识别子模块,配置成根据所述离散状态空间方程与所述风机的结构动力响应计算风机的频率、阻尼比与模态振型;An identification sub-module configured to calculate a frequency, a damping ratio, and a modal shape of the fan according to the discrete state space equation and a structural dynamic response of the fan;
稳态图子模块,配置成根据所述风机的频率与模态振型建立稳态图;A steady-state graph sub-module configured to establish a steady-state graph according to the frequency and modal shape of the fan;
提取子模块,配置成根据所述稳态图提取1p频率。An extraction submodule configured to extract a 1p frequency according to the steady-state map.
一种终端,包括存储器以及处理器,所述存储器配置成存储计算机程序,所述处理器执行所述计算机程序以使所述终端实现以上任一项所述的基于连续监测的风机1p信号识别方法。A terminal includes a memory and a processor. The memory is configured to store a computer program, and the processor executes the computer program to enable the terminal to implement the continuous monitoring-based fan 1p signal identification method according to any one of the above. .
一种计算机可读存储介质,其存储有所述终端所执行的所述计算机程序。A computer-readable storage medium stores the computer program executed by the terminal.
本申请的有益效果是:The beneficial effects of this application are:
通过对风机的结构动力响应进行分解与变换而得到对应的希尔伯特谱,并选择其中与 低频范围的集中程度最高者对应的本征模态函数作为目标本征模态函数,对目标本征模态函数进行识别而得到对应频率并建立稳态图,根据稳态图提取1p频率,实现1p信号的自动识别,具有显著的有效性与识别精度。The corresponding Hilbert spectrum is obtained by decomposing and transforming the structural dynamic response of the wind turbine, and the eigenmode function corresponding to the highest concentration in the low frequency range is selected as the target eigenmode function. The eigenmode function is used to identify the corresponding frequency and establish a steady-state map. The 1p frequency is extracted according to the steady-state map to realize the automatic recognition of the 1p signal, which has significant effectiveness and recognition accuracy.
为使本申请的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。In order to make the above-mentioned objects, features, and advantages of this application more comprehensible, preferred embodiments are described below in conjunction with the accompanying drawings and described in detail below.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本申请的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。In order to explain the technical solutions of the embodiments of the present application more clearly, the drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present application, so It should be regarded as a limitation on the scope. For those of ordinary skill in the art, other related drawings can be obtained according to these drawings without paying creative work.
图1是本申请实施例1提供的基于连续监测的风机1p信号识别方法的流程图;1 is a flowchart of a method for identifying a 1p signal of a fan based on continuous monitoring according to Embodiment 1 of the present application;
图2是本申请实施例1提供的基于连续监测的风机1p信号识别方法的步骤A的流程图;FIG. 2 is a flowchart of step A of a method for identifying a fan 1p signal based on continuous monitoring according to Embodiment 1 of the present application;
图3是本申请实施例1提供的基于连续监测的风机1p信号识别方法涉及的传感器的布置轴测结构示意图;3 is a schematic diagram of an axonometric structure of sensors involved in a method for identifying a fan 1p signal based on continuous monitoring according to Embodiment 1 of the present application;
图4是本申请实施例1提供的基于连续监测的风机1p信号识别方法的加速度响应的时域分布图;4 is a time-domain distribution diagram of an acceleration response based on a continuous monitoring of a 1p signal recognition method of a fan provided in Embodiment 1 of the present application;
图5是本申请实施例1提供的基于连续监测的风机1p信号识别方法的本征模态函数的时域分布图;FIG. 5 is a time-domain distribution diagram of an eigenmode function based on a continuous monitoring of a fan 1p signal recognition method provided in Embodiment 1 of the present application; FIG.
图6是本申请实施例1提供的基于连续监测的风机1p信号识别方法的步骤E的流程图;6 is a flowchart of step E of a method for identifying a fan 1p signal based on continuous monitoring provided in Embodiment 1 of the present application;
图7是本申请实施例1提供的基于连续监测的风机1p信号识别方法得到的稳态图;7 is a steady state diagram obtained by the method for identifying a 1p signal of a fan based on continuous monitoring provided in Embodiment 1 of the present application;
图8是基于传统方法得到的稳态图;8 is a steady state diagram obtained based on a conventional method;
图9是本申请实施例1提供的基于连续监测的风机1p信号识别方法得到的坎贝尔图;9 is a Campbell diagram obtained by the method for identifying a 1p signal of a wind turbine based on continuous monitoring provided in Embodiment 1 of the present application;
图10是本申请实施例2提供的基于连续监测的风机1p信号识别装置的结构示意图;FIG. 10 is a schematic structural diagram of a fan 1p signal recognition device based on continuous monitoring provided in Embodiment 2 of the present application;
图11是本申请实施例2提供的基于连续监测的风机1p信号识别装置的获取模块的结构示意图;11 is a schematic structural diagram of an acquisition module of a fan 1p signal identification device based on continuous monitoring provided in Embodiment 2 of the present application;
图12是本申请实施例2提供的基于连续监测的风机1p信号识别装置的识别模块的结构示意图;12 is a schematic structural diagram of an identification module of a fan 1p signal identification device based on continuous monitoring provided in Embodiment 2 of the present application;
图13是本申请实施例3提供的终端的结构示意图。FIG. 13 is a schematic structural diagram of a terminal provided in Embodiment 3 of the present application.
主要元件符号说明:Explanation of main component symbols:
100-基于连续监测的风机1p信号识别装置,110-获取模块,111-周期子模块,112-获取子模块,120-分解模块,130-变换模块,140-筛选模块,150-识别模块,151-建模子模块,152-识别子模块,153-稳态图子模块,154-提取子模块,200-终端,210-存储器,220-处理 器,230-输入单元,240-显示单元,300-风机,310-塔架,320-机舱,400-传感器。100-continuous monitoring-based fan 1p signal identification device, 110-acquisition module, 111-cycle submodule, 112-acquisition submodule, 120-decomposition module, 130-transformation module, 140-screening module, 150-identification module, 151 -Modeling sub-module, 152-recognition sub-module, 153-steady-state graph sub-module, 154-extraction sub-module, 200-terminal, 210-memory, 220-processor, 230-input unit, 240-display unit, 300 -Fan, 310-tower, 320-cabin, 400-sensor.
具体实施方式detailed description
为了便于理解本申请,下面将参照相关附图对基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质进行更全面的描述。附图中给出了基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质的优选实施例。但是,基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质可以通过许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质的公开内容更加透彻全面。In order to facilitate the understanding of the present application, a method, a device, a terminal, and a computer-readable storage medium for identifying a fan 1p signal based on continuous monitoring will be described more fully below with reference to related drawings. Preferred embodiments of a method, a device, a terminal, and a computer-readable storage medium for identifying a fan 1p signal based on continuous monitoring are shown in the drawings. However, the method, device, terminal, and computer-readable storage medium for fan 1p signal recognition based on continuous monitoring can be implemented in many different forms and are not limited to the embodiments described herein. On the contrary, the purpose of providing these embodiments is to make the disclosure of the method, device, terminal and computer-readable storage medium of the fan 1p signal recognition based on continuous monitoring more thorough and comprehensive.
需要说明的是,当元件被称为“固定于”另一个元件时,它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件时,它可以是直接连接到另一个元件或者可能同时存在居中元件。相反,当元件被称作“直接在”另一元件“上”时,不存在中间元件。本文所使用的术语“垂直的”、“水平的”、“左”和“右”以及类似的表述只是为了说明的目的。It should be noted that when an element is referred to as being "fixed to" another element, it may be directly on the other element or there may be a centered element. When an element is considered to be "connected" to another element, it can be directly connected to the other element or intervening elements may be present concurrently. In contrast, when an element is referred to as being "directly on" another element, there are no intervening elements present. The terms "vertical", "horizontal", "left" and "right" and similar expressions used herein are for illustrative purposes only.
除非另有定义,本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在基于连续监测的风机1p信号识别方法、装置、终端与计算机可读存储介质的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在限制本申请。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。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 application belongs. The terminology used in the description of the method, device, terminal, and computer-readable storage medium for fan 1p signal recognition based on continuous monitoring herein is only for the purpose of describing specific embodiments and is not intended to limit the application. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.
实施例1Example 1
请参阅图1,本实施例提供一种基于连续监测的风机1p信号识别方法,该方法包括步骤A~E:Referring to FIG. 1, this embodiment provides a method for identifying a 1p signal of a fan based on continuous monitoring. The method includes steps A to E:
步骤A:获取风机300的结构动力响应。其中,结构动力响应是指风机300在动力荷载作用下产生的速度、加速度和位移等。实际应用中,所需获取的结构动力响应,可以是由加速度传感器采集的加速度信号,也可以是由位移传感器测量的速度信号或位移信号。Step A: Obtain the structural dynamic response of the fan 300. Among them, the structural dynamic response refers to the speed, acceleration, and displacement generated by the fan 300 under the dynamic load. In practical applications, the structural dynamic response to be obtained may be an acceleration signal collected by an acceleration sensor, or a speed signal or a displacement signal measured by a displacement sensor.
示范性地,在本实施例中,风机300的结构动力响应为加速度响应。请参阅图2,在连续监测条件下,步骤A可包括步骤A1~A2:Exemplarily, in this embodiment, the structural dynamic response of the fan 300 is an acceleration response. Referring to FIG. 2, under continuous monitoring conditions, step A may include steps A1 to A2:
步骤A1:确定风机300的时不变周期。其中,时不变周期根据风机300的转速、桨距角与机舱320转角确定。风机300需要根据风向与风速,实时调整转速、桨距角与机舱320等相关参数,以保证较佳的能量转换效率,属于难以准确测量的时变结构。为了进行有效的精确监测,需要确定时不变周期。在时不变周期内,风机300符合时不变假定,即其结构特性可认为不随时间变化。Step A1: Determine the time constant period of the fan 300. Among them, the time constant period is determined according to the rotation speed of the fan 300, the pitch angle and the turning angle of the nacelle 320. The fan 300 needs to adjust the relevant parameters such as the rotation speed, the pitch angle, and the cabin 320 in real time according to the wind direction and wind speed to ensure better energy conversion efficiency, which belongs to a time-varying structure that is difficult to accurately measure. For effective and accurate monitoring, a time-invariant cycle needs to be determined. During the time-invariant period, the fan 300 conforms to the time-invariant assumption, that is, its structural characteristics can be considered to not change with time.
步骤A2:获取风机300于所述时不变周期内的结构动力响应。Step A2: Obtain the structural dynamic response of the fan 300 in the time-invariant period.
请参阅图3,在一个实际监测例中,风机300的塔架310沿高度三等分,并于三等分的四个点处分别布置采集单元。其中,塔架310的底端为4Q,在一个至简结构中,仅需于4Q处的外周布置一个传感器400,即可实现结构动力响应的采集。示范性地,传感器400可以是加速度传感器或位移传感器等类型。补充说明,图4示出了传感器400采集到的加速度响应信号。Please refer to FIG. 3. In an actual monitoring example, the tower 310 of the wind turbine 300 is divided into three parts along the height, and the acquisition units are respectively arranged at four points of the three parts. Among them, the bottom end of the tower 310 is 4Q. In a simple structure, only one sensor 400 needs to be arranged on the periphery of the 4Q to realize the structural dynamic response collection. Exemplarily, the sensor 400 may be an acceleration sensor or a displacement sensor. Additionally, FIG. 4 shows the acceleration response signal collected by the sensor 400.
步骤B:将所述风机300的结构动力响应分解为若干本征模态函数。本征模态函数即Intrinsic Mode Function(imf),应当满足以下条件:(1)在整个信号长度上,函数的极值点与过零点数目相等或至多相差一个;(2)在任意时刻,由函数极大值定义的上包络线与由极小值定义的下包络线的平均值为零,即上下包络线对称于时间轴。补充说明,图5示出了分解得到的各本征模态函数的时域分布图(imf1~6)。Step B: Decompose the structural dynamic response of the wind turbine 300 into several eigenmode functions. The eigenmode function, Intrinsic Mode Function (imf), should satisfy the following conditions: (1) the number of extreme points of the function and the zero crossings are equal or at most one difference over the entire signal length; (2) at any time, The average value of the upper envelope defined by the maximum value of the function and the lower envelope defined by the minimum value is zero, that is, the upper and lower envelopes are symmetrical to the time axis. In addition, FIG. 5 shows time-domain distribution diagrams (imf1 to 6) of each eigenmode function obtained by decomposition.
步骤C:对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值。示范性地,所述变换为希尔伯特变换。应当理解,每一本征模态函数经过希尔伯特变换后,将得到一系列的瞬时频率。将对应于同一本征模态函数的瞬时频率进行均值求取,即得到该本征模态函数的瞬时频率均值。Step C: Transform each eigenmode function to obtain a corresponding instantaneous frequency, and calculate the mean instantaneous frequency value of each eigenmode function. Exemplarily, the transform is a Hilbert transform. It should be understood that after each eigenmode function undergoes Hilbert transformation, a series of instantaneous frequencies will be obtained. The mean value of the instantaneous frequency of the eigenmode function is obtained by averaging the instantaneous frequencies corresponding to the same eigenmode function.
步骤D:以所述瞬时频率均值位于低频区段的本征模态函数作为有效本征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围。Step D: Use the eigenmode function whose mean frequency is in the low frequency region as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function. The band is a frequency range below 1 Hz.
换言之,满足以下条件的本征模态函数均为有效本征模态函数:瞬时频率均值位于低频区段。进而将所有有效本征模态函数进行相加,其和即为目标本征模态函数。In other words, the eigenmode functions that satisfy the following conditions are all valid eigenmode functions: the mean instantaneous frequency is located in the low frequency section. Then add all the effective eigenmode functions, and the sum is the target eigenmode function.
步骤E:识别所述目标本征模态函数对应的频率并建立稳态图,根据稳态图提取1p频率。其中,稳态图用于指明系统极点的位置。由于极点是系统的全局特征(在此表现为风机300的频率),随着模型阶数的增加,由阶数增加的数学模型提取到的系统极点将重复出现,并于同一图上表征,以便通过观察极点分布而找出结构的物理极点。进而,对稳态图进行特征提取,即可得到1p频率。应当理解,1p频率小于风机300的基频。Step E: Identify the frequency corresponding to the target eigenmode function and establish a steady-state map, and extract a 1p frequency according to the steady-state map. Among them, the steady state diagram is used to indicate the position of the poles of the system. Since the poles are the global characteristics of the system (represented by the frequency of the wind turbine 300 here), as the model order increases, the system poles extracted by the mathematical model of the order increase will appear repeatedly and be characterized on the same graph in order to Find the physical poles of the structure by observing the pole distribution. Furthermore, by performing feature extraction on the steady-state graph, a frequency of 1p can be obtained. It should be understood that the 1p frequency is less than the fundamental frequency of the fan 300.
请参阅图6,示范性地,步骤E包括步骤E1~E4:Referring to FIG. 6, for example, step E includes steps E1 to E4:
步骤E1:建立离散状态空间方程。离散状态空间方程基于离散化的结构动力响应而建立,于风机300的结构动力响应与频率、阻尼比、模态振型之间建立数学联系。在本实施例中,离散状态空间方程的具体形式为:Step E1: Establish a discrete state space equation. The discrete state space equation is established based on the discretized structural dynamic response, and a mathematical relationship is established between the structural dynamic response of the wind turbine 300 and the frequency, damping ratio, and modal shape. In this embodiment, the specific form of the discrete state space equation is:
Figure PCTCN2019098094-appb-000001
Figure PCTCN2019098094-appb-000001
式中,x k为离散状态向量,y k为离散输出向量,w k、ν k为白噪声项,ε为风机300 的加速度、速度与位移中的一者,A ε为离散状态矩阵,C ε为离散输出矩阵。 Where x k is a discrete state vector, y k is a discrete output vector, w k and ν k are white noise terms, ε is one of the acceleration, speed, and displacement of the fan 300, A ε is a discrete state matrix, C ε is a discrete output matrix.
步骤E2:根据所述离散状态空间方程与风机300的结构动力响应计算风机300的频率、阻尼比与模态振型。Step E2: Calculate the frequency, damping ratio and modal shape of the fan 300 according to the discrete state space equation and the structural dynamic response of the fan 300.
步骤E3:根据风机300的频率与模态振型建立稳态图。图7示出了执行步骤E3得到的一个实际稳态图。Step E3: Establish a steady-state map according to the frequency and modal shapes of the fan 300. FIG. 7 shows an actual steady state diagram obtained by performing step E3.
步骤E4:根据稳态图提取1p频率。根据图7,可以精确提取得到风机300的基频为0.6Hz,1p频率为0.41Hz。图8示出了基于传统方法得到的稳态图,该稳态图于低频区域缺乏明显的极值点分布,难以提取得到1p频率。而基于本实施例提供的基于连续监测的风机1p信号识别方法得到的稳态图(图7),在低频区域内具有显著而集中的极值点分布,为1p频率的观察与提取提供了可靠的基础。Step E4: Extract the 1p frequency according to the steady-state map. According to FIG. 7, the fundamental frequency of the fan 300 can be accurately extracted to be 0.6 Hz, and the 1p frequency is 0.41 Hz. FIG. 8 shows a steady-state map obtained based on the conventional method. The steady-state map lacks an obvious extreme point distribution in a low frequency region, and it is difficult to extract a 1p frequency. The steady-state map (Figure 7) obtained based on the continuous monitoring of the 1p signal recognition method of the fan provided in this embodiment has a significant and concentrated extreme point distribution in the low-frequency region, which provides reliable observation and extraction of the 1p frequency. The basics.
补充说明,1p频率由叶片之间的微小质量差异引起。当叶片高速旋转时,该微小质量差异会对风塔产生具有固定频率的激励力,该固定频率即为一倍转速的激励频率。It is added that the 1p frequency is caused by slight quality differences between the leaves. When the blade rotates at a high speed, the slight mass difference will generate an excitation force with a fixed frequency to the wind tower, and the fixed frequency is an excitation frequency with a double speed.
根据前述识别结果,可进一步进行相关分析。例如,计算风机300于所述时不变周期内的平均转速,并获取所需数量的频率-平均转速数对,根据所述频率-平均转速数对建立坎贝尔图。具体地,连续监测预设数量的时不变周期,获取每一时不变周期对应的风机300频率与平均转速,从而得到频率-平均转速数对。根据复数个频率-平均转速数对,即可建立基于连续监测的坎贝尔图,从而分析风机300的结构模态。Based on the aforementioned recognition results, further correlation analysis can be performed. For example, the average speed of the fan 300 during the time-invariant period is calculated, and a required number of frequency-average speed numbers are obtained, and a Campbell diagram is established according to the frequency-average speed numbers. Specifically, a preset number of time-invariant periods are continuously monitored, and the frequency and average speed of the fan 300 corresponding to each time-invariant period are obtained, so as to obtain a frequency-average speed number pair. According to a plurality of frequency-average speed number pairs, a Campbell diagram based on continuous monitoring can be established to analyze the structural mode of the fan 300.
例如,基于长达两年的加速度连续监测,得以建立长期监测得到的坎贝尔图。图9示出了一个根据加速度连续监测得到的坎贝尔图,根据该坎贝尔图可以得到所监测的风机300的1p频率(0.41Hz)。For example, based on continuous acceleration monitoring for up to two years, a Campbell diagram obtained from long-term monitoring can be established. FIG. 9 shows a Campbell diagram obtained based on continuous monitoring of acceleration. According to the Campbell diagram, the 1p frequency (0.41 Hz) of the fan 300 can be obtained.
实施例2Example 2
请参阅图10,本实施例提供一种基于连续监测的风机1p信号识别装置100,该装置包括:Referring to FIG. 10, this embodiment provides a device 1p signal identification device 100 based on continuous monitoring. The device includes:
获取模块110,配置成获取风机300的结构动力响应;An obtaining module 110 configured to obtain a structural dynamic response of the wind turbine 300;
分解模块120,配置成将所述风机300的结构动力响应分解为若干本征模态函数;A decomposition module 120 configured to decompose the structural dynamic response of the wind turbine 300 into several eigenmode functions;
变换模块130,配置成对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值;The transformation module 130 is configured to transform each eigenmode function to obtain an instantaneous frequency corresponding thereto, and calculate an average instantaneous frequency value of each eigenmode function;
筛选模块140,配置成以所述瞬时频率均值位于低频区段的本征模态函数作为有效本征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围;The screening module 140 is configured to use the eigenmode function with the instantaneous frequency mean value located in the low-frequency section as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function, The low frequency section is a frequency range below 1 Hz;
识别模块150,配置成识别所述目标本征模态函数对应的频率并建立稳态图,根据稳态图提取1p频率。The identification module 150 is configured to identify a frequency corresponding to the target eigenmode function and establish a steady-state map, and extract a 1p frequency according to the steady-state map.
请参阅图11,示范性地,所述获取模块110包括:Referring to FIG. 11, for example, the obtaining module 110 includes:
周期子模块111,配置成确定风机300的时不变周期;A period sub-module 111 configured to determine a time-invariant period of the fan 300;
获取子模块112,配置成获取风机300于所述时不变周期内的结构动力响应。The obtaining sub-module 112 is configured to obtain a structural dynamic response of the wind turbine 300 in the time-invariant period.
请参阅图12,示范性地,所述识别模块150包括:Referring to FIG. 12, for example, the identification module 150 includes:
建模子模块151,配置成建立离散状态空间方程;A modeling sub-module 151 configured to establish a discrete state space equation;
识别子模块152,配置成根据所述离散状态空间方程与所述风机300的结构动力响应计算风机300的频率、阻尼比与模态振型;An identification submodule 152 configured to calculate a frequency, a damping ratio, and a modal shape of the fan 300 according to the discrete state space equation and the structural dynamic response of the fan 300;
稳态图子模块153,配置成根据所述风机300的频率与模态振型建立稳态图;A steady-state map sub-module 153 configured to establish a steady-state map according to the frequency and modal shapes of the fan 300;
提取子模块154,配置成根据所述稳态图提取1p频率。The extraction sub-module 154 is configured to extract a 1p frequency according to the steady-state map.
实施例3Example 3
请参阅图13,本实施例提供一种终端200,该终端200包括存储器210以及处理器220,存储器210配置成存储计算机程序,处理器220执行计算机程序以使终端200实现以上所述的基于连续监测的风机1p信号识别方法。13, this embodiment provides a terminal 200. The terminal 200 includes a memory 210 and a processor 220. The memory 210 is configured to store a computer program, and the processor 220 executes the computer program to enable the terminal 200 to implement the continuous-based Identification method of the monitored 1p signal of the fan.
其中,终端200包括不具备移动通信能力的终端设备(比如计算机和服务器等),亦包括移动终端(比如智能电话、平板电脑、车载电脑和智能穿戴设备等)。The terminal 200 includes terminal devices (such as computers and servers) that do not have mobile communication capabilities, and also includes mobile terminals (such as smart phones, tablet computers, on-board computers, and smart wearable devices).
存储器210可包括存储程序区和存储数据区。其中,存储程序区可存储操作系统以及至少一个功能所需的应用程序(比如声音播放功能和图像播放功能等)等;存储数据区可存储根据终端200的使用所创建的数据(比如音频数据和备份文件等)等。此外,存储器210可以包括高速随机存取存储器,还可以包括非易失性存储器(例如至少一个磁盘存储器件、闪存器件或其他易失性固态存储器件)。The memory 210 may include a storage program area and a storage data area. The storage program area can store an operating system and applications required for at least one function (such as a sound playback function and an image playback function), and the like; the storage data area can store data (such as audio data and Backup files, etc.). In addition, the memory 210 may include a high-speed random access memory, and may further include a non-volatile memory (for example, at least one magnetic disk storage device, a flash memory device, or other volatile solid-state storage device).
优选地,终端200还包括输入单元230与显示单元240。其中,输入单元230配置成接收用户输入的各项指令或参数(包括预设滚动方式、预设时间间隔与预设滚动次数),包括鼠标、键盘、触控面板及其他输入设备。显示单元240配置成显示终端200的各种输出信息(包括网页页面和参数配置界面等),包括显示面板。Preferably, the terminal 200 further includes an input unit 230 and a display unit 240. The input unit 230 is configured to receive various instructions or parameters (including a preset scrolling method, a preset time interval, and a preset scrolling number) input by a user, including a mouse, a keyboard, a touch panel, and other input devices. The display unit 240 is configured to display various output information (including a webpage page and a parameter configuration interface, etc.) of the terminal 200, including a display panel.
在此一并提供一种计算机可读存储介质,其存储有终端200所执行的所述计算机程序。A computer-readable storage medium is also provided, which stores the computer program executed by the terminal 200.
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和结构图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may also be implemented in other ways. The device embodiments described above are merely schematic. For example, the flowcharts and structure diagrams in the accompanying drawings show the possible architectures and functions of the devices, methods, and computer program products according to various embodiments of the present application. And operation. In this regard, each block in the flowchart or block diagram may represent a module, a program segment, or a part of code, which contains one or more components for implementing a specified logical function Executable instructions.
也应当注意,在作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中 所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。It should also be noted that in alternative implementations, the functions labeled in the blocks may also occur in a different order than those labeled in the figures. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending on the functions involved.
也要注意的是,结构图和/或流程图中的每个方框、以及结构图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。It should also be noted that each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart can be a dedicated hardware-based system that performs a specified function or action To achieve, or can be implemented by a combination of dedicated hardware and computer instructions.
另外,在本申请各个实施例中的各功能模块或单元可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或更多个模块集成形成一个独立的部分。In addition, the functional modules or units in the various embodiments of the present application may be integrated together to form an independent part, or each of the modules may exist separately, or two or more modules may be integrated to form an independent part.
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是智能手机、个人计算机、服务器或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially a part that contributes to the existing technology or a part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including Several instructions are used to cause a computer device (which may be a smart phone, a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in the embodiments of the present application. The foregoing storage media include: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disks or optical disks and other media that can store program codes .
在这里示出和描述的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制,因此,示例性实施例的其他示例可以具有不同的值。In all the examples shown and described herein, any specific value should be construed as exemplary only and not as a limitation, so other examples of the exemplary embodiments may have different values.
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。It should be noted that similar reference numerals and letters indicate similar items in the following drawings, so once an item is defined in one drawing, it need not be further defined and explained in subsequent drawings.
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation manners of the present application, and the descriptions thereof are more specific and detailed, but cannot be understood as limiting the scope of the present application. It should be noted that, for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of this application shall be subject to the appended claims.

Claims (10)

  1. 一种基于连续监测的风机1p信号识别方法,其特征在于,包括:A method for identifying a 1p signal of a fan based on continuous monitoring is characterized in that it includes:
    获取风机的结构动力响应;Obtain the structural dynamic response of the fan;
    将所述风机的结构动力响应分解为若干本征模态函数;Decomposing the structural dynamic response of the wind turbine into several eigenmode functions;
    对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值;Transform each eigenmode function to obtain the instantaneous frequency corresponding to it, and calculate the mean instantaneous frequency value of each eigenmode function;
    以所述瞬时频率均值位于低频区段的本征模态函数作为有效本征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围;The eigenmode function with the instantaneous frequency mean value located in the low frequency section is taken as the effective eigenmode function, and the sum of the effective eigenmode functions is used as the target eigenmode function. Frequency range below 1Hz;
    识别所述目标本征模态函数对应的频率并建立稳态图,根据所述稳态图提取1p频率。The frequency corresponding to the target eigenmode function is identified and a steady state map is established, and a 1p frequency is extracted according to the steady state map.
  2. 根据权利要求1所述的基于连续监测的风机1p信号识别方法,其特征在于,“识别所述目标本征模态函数对应的频率并建立稳态图”包括:The method for identifying a 1p signal of a fan based on continuous monitoring according to claim 1, characterized in that "identifying the frequency corresponding to the target eigenmode function and establishing a steady state map" includes:
    建立离散状态空间方程;Establish discrete state space equations;
    根据所述离散状态空间方程与所述风机的结构动力响应计算风机的频率、阻尼比与模态振型;Calculating the fan frequency, damping ratio and modal shape according to the discrete state space equation and the structural dynamic response of the fan;
    根据所述风机的频率与模态振型建立稳态图。A steady state diagram is established according to the frequency and modal shape of the fan.
  3. 根据权利要求1所述的基于连续监测的风机1p信号识别方法,其特征在于,所述变换为希尔伯特变换。The method for identifying a fan 1p signal based on continuous monitoring according to claim 1, wherein the transform is a Hilbert transform.
  4. 根据权利要求1所述的基于连续监测的风机1p信号识别方法,其特征在于,所述风机的结构动力响应为加速度、速度或位移响应。The method for identifying a 1p signal of a fan based on continuous monitoring according to claim 1, wherein the structural dynamic response of the fan is an acceleration, speed, or displacement response.
  5. 根据权利要求4所述的基于连续监测的风机1p信号识别方法,其特征在于,“获取风机的结构动力响应”包括:The method for identifying a 1p signal of a wind turbine based on continuous monitoring according to claim 4, wherein "obtaining a structural dynamic response of the wind turbine" includes:
    确定风机的时不变周期;Determine the time constant period of the fan;
    获取所述风机于所述时不变周期内的结构动力响应。Obtain the structural dynamic response of the wind turbine in the time-invariant period.
  6. 一种基于连续监测的风机1p信号识别装置,其特征在于,包括:A fan 1p signal identification device based on continuous monitoring is characterized in that it includes:
    获取模块,配置成获取风机的结构动力响应;An acquisition module configured to acquire a structural dynamic response of the wind turbine;
    分解模块,配置成将所述风机的结构动力响应分解为若干本征模态函数;A decomposition module configured to decompose the structural dynamic response of the wind turbine into several eigenmode functions;
    变换模块,配置成对每一本征模态函数进行变换以得到与之对应的瞬时频率,并计算所述每一本征模态函数的瞬时频率均值;A transformation module configured to transform each eigenmode function to obtain a corresponding instantaneous frequency, and calculate an average value of the instantaneous frequency of each eigenmode function;
    筛选模块,配置成以所述瞬时频率均值位于低频区段的本征模态函数作为有效本 征模态函数,并以所述有效本征模态函数之和作为目标本征模态函数,所述低频区段为低于1Hz的频率范围;A screening module configured to use the eigenmode function with the instantaneous frequency mean value located in a low frequency section as the effective eigenmode function, and use the sum of the effective eigenmode functions as the target eigenmode function The low frequency range is a frequency range below 1 Hz;
    识别模块,配置成识别所述目标本征模态函数对应的频率并建立稳态图,根据稳态图提取1p频率。The identification module is configured to identify a frequency corresponding to the target eigenmode function and establish a steady state map, and extract a 1p frequency according to the steady state map.
  7. 根据权利要求6所述的基于连续监测的风机1p信号识别装置,其特征在于,所述变换为希尔伯特变换。The fan 1p signal identification device based on continuous monitoring according to claim 6, wherein the transformation is a Hilbert transformation.
  8. 根据权利要求7所述的基于连续监测的风机1p信号识别装置,其特征在于,所述识别模块包括:The fan 1p signal identification device based on continuous monitoring according to claim 7, wherein the identification module comprises:
    建模子模块,配置成建立离散状态空间方程;A modeling sub-module configured to establish a discrete state space equation;
    识别子模块,配置成根据所述离散状态空间方程与所述风机的结构动力响应计算风机的频率、阻尼比与模态振型;An identification sub-module configured to calculate a frequency, a damping ratio, and a modal shape of the fan according to the discrete state space equation and a structural dynamic response of the fan;
    稳态图子模块,配置成根据所述风机的频率与模态振型建立稳态图;A steady-state graph sub-module configured to establish a steady-state graph according to the frequency and modal shape of the fan;
    提取子模块,配置成根据所述稳态图提取1p频率。An extraction submodule configured to extract a 1p frequency according to the steady-state map.
  9. 一种终端,其特征在于,包括存储器以及处理器,所述存储器配置成存储计算机程序,所述处理器执行所述计算机程序以使所述终端实现权利要求1~5中任一项所述的基于连续监测的风机1p信号识别方法。A terminal, comprising a memory and a processor, wherein the memory is configured to store a computer program, and the processor executes the computer program to enable the terminal to implement the device according to any one of claims 1 to 5. Method for identifying 1p signal of wind turbine based on continuous monitoring.
  10. 一种计算机可读存储介质,其特征在于,其存储有权利要求9所述的终端所执行的所述计算机程序。A computer-readable storage medium, characterized in that it stores the computer program executed by the terminal according to claim 9.
PCT/CN2019/098094 2018-07-06 2019-07-29 Continuous monitoring-based method and device for identifying 1p signal of wind turbine, terminal, and computer readable storage medium WO2020007375A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810738808.8A CN109100103B (en) 2018-07-06 2018-07-06 Fan 1p signal identification method, device, terminal and computer readable storage medium based on continuous monitoring
CN201810738808.8 2018-07-06

Publications (1)

Publication Number Publication Date
WO2020007375A1 true WO2020007375A1 (en) 2020-01-09

Family

ID=64845754

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/098094 WO2020007375A1 (en) 2018-07-06 2019-07-29 Continuous monitoring-based method and device for identifying 1p signal of wind turbine, terminal, and computer readable storage medium

Country Status (2)

Country Link
CN (1) CN109100103B (en)
WO (1) WO2020007375A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109100103B (en) * 2018-07-06 2020-04-14 哈尔滨工业大学(深圳) Fan 1p signal identification method, device, terminal and computer readable storage medium based on continuous monitoring

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102467654A (en) * 2010-11-03 2012-05-23 西安电子科技大学 Structural modal parameter identification method
US20120253694A1 (en) * 2011-03-28 2012-10-04 Hong-Tsu Young Method and apparatus for judging status of mechanical system
CN106650218A (en) * 2016-10-19 2017-05-10 上海电机学院 Harmonic analysis method based on complementary ensemble empirical mode decomposition (CEEMD) algorithm and Hilbert transform
CN107192524A (en) * 2017-04-05 2017-09-22 天津大学 A kind of wind-powered electricity generation structure operational modal parameter recognition methods for considering strong harmonic wave interference
CN108229382A (en) * 2017-12-29 2018-06-29 广州供电局有限公司 Vibration signal characteristics extracting method, device, storage medium and computer equipment
CN109100103A (en) * 2018-07-06 2018-12-28 哈尔滨工业大学(深圳) Based on blower 1p signal recognition method, device, terminal and the computer readable storage medium continuously monitored

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102467654A (en) * 2010-11-03 2012-05-23 西安电子科技大学 Structural modal parameter identification method
US20120253694A1 (en) * 2011-03-28 2012-10-04 Hong-Tsu Young Method and apparatus for judging status of mechanical system
CN106650218A (en) * 2016-10-19 2017-05-10 上海电机学院 Harmonic analysis method based on complementary ensemble empirical mode decomposition (CEEMD) algorithm and Hilbert transform
CN107192524A (en) * 2017-04-05 2017-09-22 天津大学 A kind of wind-powered electricity generation structure operational modal parameter recognition methods for considering strong harmonic wave interference
CN108229382A (en) * 2017-12-29 2018-06-29 广州供电局有限公司 Vibration signal characteristics extracting method, device, storage medium and computer equipment
CN109100103A (en) * 2018-07-06 2018-12-28 哈尔滨工业大学(深圳) Based on blower 1p signal recognition method, device, terminal and the computer readable storage medium continuously monitored

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
LIU XIN: "Automated operational modal parameter extraction with stochastic subspace identification", SCIENCE-ENGINEERING (B), CHINA MASTER'S THESES FULL-TEXT DATABASE, no. 6, 15 June 2018 (2018-06-15), pages 18 - 19 *

Also Published As

Publication number Publication date
CN109100103B (en) 2020-04-14
CN109100103A (en) 2018-12-28

Similar Documents

Publication Publication Date Title
Devriendt et al. Structural health monitoring of offshore wind turbines using automated operational modal analysis
Bódai et al. Annual variability in a conceptual climate model: Snapshot attractors, hysteresis in extreme events, and climate sensitivity
WO2021227468A1 (en) State preview method and system for key equipment of nuclear power plant, and device and storage medium
Dong et al. Vibration source identification of offshore wind turbine structure based on optimized spectral kurtosis and ensemble empirical mode decomposition
CN105005581A (en) Method and system for processing wind resource data of wind power plant
WO2020007375A1 (en) Continuous monitoring-based method and device for identifying 1p signal of wind turbine, terminal, and computer readable storage medium
Meddage et al. Predicting transient wind loads on tall buildings in three-dimensional spatial coordinates using machine learning
Hu et al. Resonance monitoring of a horizontal wind turbine by strain-based automated operational modal analysis
Zhu et al. Wind turbine blade defect detection based on acoustic features and small sample size
Zhang et al. Influence of the heights of low-level jets on power and aerodynamic loads of a horizontal axis wind turbine rotor
Castellani et al. Wind turbine loads induced by terrain and wakes: an experimental study through vibration analysis and computational fluid dynamics
CN112651622A (en) Electric energy quality evaluation method and system
Li et al. Weak signal detection method based on the coupled Lorenz system and its application in rolling bearing fault diagnosis
CN109100102B (en) Fan modal analysis method, device, terminal and computer readable storage medium based on strain continuous monitoring
Zhang et al. The method for determining optimal analysis length of vibration data based on improved multiscale permutation entropy
TW201245997A (en) Method and apparatus for evaluating efficiency of wind generator
Hassanzadeh Best location of the small wind turbine around a single high-rise building
Wang et al. A method of modal parameter identification for wind turbine blade based on binocular dynamic photogrammetry
CN114841077A (en) Wind power prediction method, device and medium
Aldereguía Sánchez et al. Experimental Simulation of Thunderstorm Profiles in an Atmospheric Boundary Layer Wind Tunnel
Wang et al. Mechanical Fault Prognosis through Spectral Analysis of Vibration Signals
Zhao et al. Refined composite hierarchical multiscale Lempel-Ziv complexity: A quantitative diagnostic method of multi-feature fusion for rotating energy devices
Zhou et al. Data Inspecting and Denoising Method for Data‐Driven Stochastic Subspace Identification
Abu Hasan et al. Automated harmonic signal removal technique using stochastic subspace-based image feature extraction
Li et al. Experimental investigation of wind loads on wind turbine blade under various turbulent flows

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19831408

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19831408

Country of ref document: EP

Kind code of ref document: A1