CN112417774B - Identification method and device for rotating stall of compressor and computer equipment - Google Patents

Identification method and device for rotating stall of compressor and computer equipment Download PDF

Info

Publication number
CN112417774B
CN112417774B CN202011095061.2A CN202011095061A CN112417774B CN 112417774 B CN112417774 B CN 112417774B CN 202011095061 A CN202011095061 A CN 202011095061A CN 112417774 B CN112417774 B CN 112417774B
Authority
CN
China
Prior art keywords
frequency
compressor
rotating stall
working condition
rotating
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
Application number
CN202011095061.2A
Other languages
Chinese (zh)
Other versions
CN112417774A (en
Inventor
王睿
王海伦
李宏坤
孟继纲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Turbo Machinery Co Ltd
Original Assignee
Shenyang Turbo Machinery Co Ltd
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 Shenyang Turbo Machinery Co Ltd filed Critical Shenyang Turbo Machinery Co Ltd
Priority to CN202011095061.2A priority Critical patent/CN112417774B/en
Publication of CN112417774A publication Critical patent/CN112417774A/en
Application granted granted Critical
Publication of CN112417774B publication Critical patent/CN112417774B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04BPOSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
    • F04B51/00Testing machines, pumps, or pumping installations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Positive-Displacement Air Blowers (AREA)

Abstract

The application discloses a method, a device and computer equipment for identifying rotating stall of a compressor, relates to the field of compressor control, and can solve the problem that the rotating stall of the compressor cannot be accurately identified and accurately predicted. The method comprises the following steps: acquiring each rotating speed frequency of the operation of the compressor impeller; performing constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve; determining a rotating stall datum point corresponding to each rotating speed frequency through unsteady CFD analysis and Fourier transform analysis based on the flow-pressure ratio characteristic curve; and constructing an early warning datum line by using the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line. The method is suitable for identifying and controlling the rotating stall of the compressor.

Description

Identification method and device for rotating stall of compressor and computer equipment
Technical Field
The present disclosure relates to the field of compressor control, and in particular, to a method and apparatus for identifying rotating stall of a compressor, and a computer device.
Background
Rotating stall is one of the most common instabilities of compressors. Rotating stall can deteriorate the flow conditions in the compressor, the pressure ratio drops, and the flow and pressure fluctuate over time. At a certain rotational speed, when the inlet flow is reduced to a certain value, the unit will generate a strong rotating stall, which will further cause an unstable aerodynamic phenomenon, i.e. surge, with a greater risk of the whole compressor system. In addition, when the rotating stall is carried out, the compressor blade is subjected to a periodic exciting force, and if the frequency of the rotating stall is matched with the natural frequency of the blade, strong vibration is caused, so that the blade is fatigued and damaged to cause accidents.
At present, the rotating stall of the compressor cannot be effectively and accurately identified, so that how to combine production practice and effectively early warn the rotating stall in advance in the actual running process of the compressor, thereby avoiding the occurrence of the rotating stall phenomenon, and becoming the current difficult problem to be solved urgently.
Disclosure of Invention
In view of this, the application provides a method, a device and a computer device for identifying rotating stall of a compressor, which mainly solve the problem that the rotating stall of the compressor cannot be accurately and effectively identified.
According to one aspect of the present application, there is provided a method of identifying rotating stall of a compressor, the method comprising:
acquiring each rotating speed frequency of the operation of the compressor impeller;
performing constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve;
determining a rotating stall datum point corresponding to each rotating speed frequency through unsteady CFD analysis and Fourier transform analysis based on the flow-pressure ratio characteristic curve;
and constructing an early warning datum line by using the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line.
According to another aspect of the present application, there is provided an apparatus for identifying rotating stall of a compressor, the apparatus comprising:
The acquisition module is used for acquiring each rotating speed frequency of the operation of the compressor impeller;
the analysis module is used for respectively carrying out constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve;
the determining module is used for determining the rotating stall datum points corresponding to the rotating speed frequencies based on the flow-pressure ratio characteristic curves through unsteady CFD analysis and Fourier transform analysis;
and the construction module is used for constructing an early warning datum line by utilizing the rotating stall datum point and identifying the rotating stall of the compressor based on the early warning datum line.
According to yet another aspect of the present application, there is provided a non-transitory readable storage medium having stored thereon a computer program which when executed by a processor implements the above-described method of identifying rotating stall of a compressor.
According to yet another aspect of the present application, there is provided a computer device comprising a non-volatile readable storage medium, a processor and a computer program stored on the non-volatile readable storage medium and executable on the processor, the processor implementing the above method for identifying rotating stall of a compressor when executing the program.
By means of the technical scheme, the identification method, the identification device and the computer equipment for the rotating stall of the compressor relate to the field of compressor control, and can firstly obtain each rotating speed frequency of a compressor impeller in a working state, and further obtain a characteristic curve of the compressor pressure ratio changing along with flow by respectively carrying out constant CFD analysis on each rotating speed frequency. Based on the flow-pressure ratio characteristic curve, a preset number of working condition points which are likely to generate rotating stall can be extracted, and the rotating stall datum point corresponding to the maximum flow of the rotating stall phenomenon is extracted through unsteady CFD analysis and Fourier transform analysis on the extracted working condition points. After the rotating stall datum points corresponding to the rotating speed frequencies are extracted, an early warning datum line can be constructed by utilizing the rotating stall datum points, and the early warning datum line is preset in an alarm control system, so that when the early warning datum line is triggered or is about to be triggered in actual operation of the compressor, an alarm is automatically started, personnel in a control room are prompted to take corresponding measures, and therefore rotating stall of the compressor is restrained. Through the technical scheme in this application, can apply the CFD analog method of computer technology, the early warning datum point that each rotational speed frequency corresponds is accurately discerned to establish the early warning datum line through the early warning datum point, in order to realize when different rotational speed frequencies, carry out effective early warning through the early warning datum line rotatory stall phenomenon, thereby avoid the harm that rotatory stall phenomenon caused to the compressor.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute an undue limitation to the present application. In the drawings:
FIG. 1 is a flow chart illustrating a method for identifying rotating stall of a compressor according to an embodiment of the present application;
FIG. 2 is a flow chart illustrating another method for identifying rotating stall of a compressor according to an embodiment of the present application;
FIG. 3 illustrates an example flow diagram of a compressor rotating stall identification provided by embodiments of the present application;
FIG. 4 is a schematic diagram of a device for identifying rotating stall of a compressor according to an embodiment of the present disclosure;
FIG. 5 shows a schematic structural diagram of another compressor rotating stall identification apparatus provided by an embodiment of the present application;
fig. 6 shows a schematic structural diagram of a computer device according to an embodiment of the present invention.
Detailed Description
The present application will be described in detail hereinafter with reference to the accompanying drawings in conjunction with embodiments. It should be noted that, without conflict, the embodiments and features of the embodiments in the present application may be combined with each other.
Aiming at the problem that the rotating stall of the compressor cannot be accurately and effectively identified at present, the embodiment of the application provides an identification method of the rotating stall of the compressor, as shown in fig. 1, the method comprises the following steps:
101. and acquiring each rotating speed frequency of the operation of the compressor impeller.
In a specific application scenario, according to actual application requirements of clients, the same impeller can correspond to a plurality of rotating speed frequencies in the working process of the compressor. After each rotational speed frequency of the operation of the compressor impeller is determined, an early warning reference line can be constructed for each determined rotational speed frequency through the following steps of embodiments 102-104, so that effective early warning and recognition of the phenomenon of rotating stall can be realized through the early warning reference line in the use process of the compressor by a user.
It should be noted that, according to the principle of centrifugal compressor, keemunda, definition of compressor stage in the mechanical industry press: "an impeller and all the fixed elements associated therewith constitute a stage". And definition of compressor stage in the university of western traffic press, edited by axial compressor principle and pneumatic design Li Chaojun Yu Wenlong: "a row of moving blades and a row of stationary blades immediately following the row constitute one stage of an axial compressor", the compressor to which the present application refers is not limited to a centrifugal compressor, but may include an axial compressor and the like.
For the identification system of the rotating stall of the compressor, the execution main body of the method can acquire each rotating speed frequency of the operation of the compressor impeller in the early warning system, and each flow-pressure ratio characteristic curve is obtained by respectively carrying out constant CFD analysis on each rotating speed frequency; based on the flow-pressure ratio characteristic curve, determining a rotating stall datum point corresponding to each rotating speed frequency through unsteady CFD analysis and Fourier transform analysis; and constructing an early warning reference line by utilizing the rotating stall reference point, and storing the early warning reference line in an early warning module so as to effectively monitor and early warn the rotating stall of the compressor based on the early warning reference line in the actual running process of the compressor.
102. And (3) respectively carrying out constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve.
In the CFD calculation method, steady CFD analysis is steady-state analysis corresponding to the condition of large span of working condition points, and the calculation process does not involve time variable. Therefore, for the present embodiment, by performing the constant CFD analysis on each rotational frequency, a flow-pressure ratio characteristic curve with no time-dependent parameter change can be obtained.
103. Based on the flow-pressure ratio characteristic curve, the rotating stall datum point corresponding to each rotating speed frequency is determined through unsteady CFD analysis and Fourier transformation analysis.
For the embodiment, in a specific application scenario, after a flow-pressure ratio characteristic curve is constructed, a preset number of possible stall operating points may be initially extracted from the flow-pressure ratio characteristic curve, and a rotating stall reference point may be further determined from the extracted operating points by performing unsteady CFD analysis and fourier transform analysis on the extracted operating points, where the rotating stall reference point is a reference point when the rotating stall phenomenon corresponds to the maximum flow.
104. And constructing an early warning datum line by utilizing the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line.
In a specific application scenario, when the compressor is judged to generate the rotating stall phenomenon according to the early warning datum line, early warning prompt information of the rotating stall can be output, wherein the prompt information can comprise one or more of text prompt information, picture prompt information, audio prompt information, video prompt information, lamplight prompt information, vibration prompt information and the like.
Through the identification method of the rotating stall of the compressor in the embodiment, each rotating speed frequency of the compressor impeller in the working state can be obtained first, and further, the characteristic curve of the compressor pressure ratio changing along with the flow is obtained through respectively carrying out constant CFD analysis on each rotating speed frequency. Based on the flow-pressure ratio characteristic curve, a preset number of working condition points which are likely to generate rotating stall can be extracted, and the rotating stall datum point corresponding to the maximum flow of the rotating stall phenomenon is extracted through unsteady CFD analysis and Fourier transform analysis on the extracted working condition points. After the rotating stall datum points corresponding to the rotating speed frequencies are extracted, an early warning datum line can be constructed by utilizing the rotating stall datum points, and the early warning datum line is preset in an alarm control system, so that when the early warning datum line is triggered or is about to be triggered in actual operation of the compressor, an alarm is automatically started, personnel in a control room are prompted to take corresponding measures, and therefore rotating stall of the compressor is restrained. Through the technical scheme in this application, can apply the CFD analog method of computer technology, the early warning datum point that each rotational speed frequency corresponds is accurately discerned to establish the early warning datum line through the early warning datum point, in order to realize when different rotational speed frequencies, carry out effective early warning through the early warning datum line rotatory stall phenomenon, thereby avoid the harm that rotatory stall phenomenon caused to the compressor.
Further, as a refinement and expansion of the foregoing embodiment, for a complete description of the implementation procedure in this embodiment, another method for identifying rotating stall of a compressor is provided, as shown in fig. 2, and the method includes:
201. and acquiring each rotating speed frequency of the operation of the compressor impeller.
The generation mechanism of the rotating stall is known, and the rotating stall mainly generates stall separation groups in impeller flow channels by means of unsteady airflow pulsation in the impeller machinery, and the stall separation groups sequentially occur in each flow channel of the impeller along the opposite direction of rotation of the impeller at a certain propagation speed. Thus, to detect rotating stall, each rotational frequency at which the compressor wheel operates may first be obtained.
202. And (3) respectively carrying out constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve.
In a specific application scenario, each rotational speed frequency can correspondingly generate a flow-pressure ratio characteristic curve, wherein the abscissa of the flow-pressure ratio characteristic curve corresponds to a flow value and the ordinate corresponds to a pressure ratio value.
203. And extracting a preset number of first working condition points meeting preset conditions based on the flow-pressure ratio characteristic curve.
Wherein, according to the experience coefficient, the preset number may be set to be 7 or more, and in this embodiment, in order to reduce the statistical workload, the preset number may be preferentially set to be 7.
For the present embodiment, in a specific application scenario, the embodiment step 203 may specifically include: extracting a first target working condition point corresponding to the maximum pressure ratio and each second target working condition point which is separated from the first target working condition point by a preset flow interval from the first target working condition point in a flow-pressure ratio characteristic curve; and determining the first target operating point and the second target operating point as the first operating point.
The specific number of the preset flow intervals can be set according to the preset extraction number of the first working condition points, the value of each preset flow interval can be set according to the preset flow difference value between the adjacent first working condition points, for example, when the preset extraction number of the first working condition points corresponds to 7, in order to ensure equal flow difference between the adjacent first working condition points, if the preset flow difference value between the first working condition points is set to be 1%, the preset flow intervals between the first working condition points can be determined to comprise +/-1%, +/-2% and +/-3%, so that after the first target working condition point corresponding to the maximum pressure ratio is extracted according to the parabola corresponding to the flow-pressure ratio characteristic curve, the flow value corresponding to the first target working condition point can be taken as 100%, the mass flow of the inlet is m, and the second target working condition point can be extracted according to the preset flow interval +/-1% firstly when the preset flow difference value between the first working condition points is set to be equal to m: 99% m,101% m; extracting a second target working condition point according to the preset flow interval of +/-2 percent: 98% m,102% m; extracting a second target working condition point according to the preset flow interval of +/-3 percent: 97% m,103% m; and uniformly determining the first target working point and the second target working point as first working points, and further obtaining seven first working points with the flow difference of 1% between adjacent working points, namely 97% m,98% m,99% m,101% m,102% m and 103% m.
204. A rotating stall reference point is determined from the first operating point by performing an unsteady CFD analysis and a Fourier transform analysis on the first operating point.
In a specific application scenario, the embodiment step 204 may specifically include:
s1, performing unsteady CFD analysis and Fourier transform analysis on a first working condition point to obtain a rotating stall analysis result.
For an embodiment, the rotating stall reference point when the rotating stall phenomenon corresponds to the maximum flow rate may be determined from the first operating points by sequentially performing an unsteady CFD analysis and a fourier transform analysis on each of the first operating points. In a specific application scene, when the compressor rotates and stalls, a relatively obvious stall characteristic frequency appears at a low frequency position within the blade pass frequency in the frequency spectrum of the pressure pulsation signal, and modulation phenomena appear at two sides of the blade pass frequency. Therefore, in the detection of the rotating stall of the compressor, whether the rotating stall of the compressor occurs can be further determined by detecting whether the modulation phenomenon occurs on both sides of the blade pass frequency in the frequency spectrum of the pressure pulse signal and comparing the specific characteristic frequency with the amplitude of the She Tong frequency.
Correspondingly, the embodiment step S1 may specifically include: sequentially performing unsteady CFD simulation on the first working condition points according to the sequence from the large flow value to the small flow value to obtain pressure signals of the first working condition points corresponding to preset time steps; performing Fourier transform analysis on the pressure signal to obtain a frequency spectrum of the pressure pulse signal; according to She Tong frequency in the frequency spectrum of the pressure pulse signal, carrying out band-pass filtering and envelope demodulation analysis, and detecting whether modulation occurs at two sides of She Tong frequency; if the modulation phenomenon appears on the two sides of the leaf-pass frequency, comparing the specific characteristic frequency with the amplitude of She Tong frequency in the frequency spectrum of the pressure pulsation signal; if the amplitude of the specific characteristic frequency is determined to be larger than the amplitude of She Tong frequency of the preset multiple, judging that the compressor is in rotating stall at the first working point; and if the modulation phenomenon does not occur on the two sides of the blade pass frequency, judging that the rotating stall of the compressor does not occur at the first working condition point.
In a specific application scenario, performing unsteady CFD simulation on the first working condition point in sequence according to the sequence of the flow values from large to small to obtain a pressure signal of the first working condition point corresponding to a preset time step, which specifically may include: performing iterative convergence calculation on the first working condition point according to preset phase parameters of unsteady CFD simulation; and if the unsteady CFD simulation is judged to reach a convergence state, extracting a pressure value of the first working condition point corresponding to a preset time step as a pressure signal. According to the empirical coefficient, the preset time step may be set to 8192 time steps, and the preset phase parameter may correspond to the phase setting of the unsteady CFD simulation, where 1 ° in the circumferential direction is used as 1 phase position, that is, at least 360 phase positions are set for one rotation of the impeller, and each phase position is set to at least 20 steps to perform iterative convergence calculation.
Correspondingly, the most obvious characteristic frequency can be extracted from the frequency spectrum of the pressure pulsation signal as stall characteristic frequency through bandpass filtering and envelope demodulation, and in general, the most obvious characteristic frequency in the frequency spectrum of the pressure pulsation signal is the frequency component modulated at both sides of She Tong frequency. Analysis shows that before the compressor stalls, no frequency component with the amplitude being greater than 1/3 leaf passband amplitude exists in the low frequency range of 5 Hz-1/2 rotational speed frequency in the frequency spectrum of the pressure pulsation signal; in the stall range from the initial stall to the surge, a frequency component with an amplitude greater than 1/3 of the lobe passband amplitude appears in the low frequency range of 5 Hz-1/2 of the rotational speed frequency in the frequency spectrum of the pressure pulsation signal, and the frequency component is a specific characteristic frequency. Specifically, for the embodiment in step S1, according to She Tong frequencies in the frequency spectrum of the pressure pulse signal, band-pass filtering and envelope demodulation analysis are performed, and detecting whether modulation occurs at two sides of She Tong frequencies may specifically include: extracting specific characteristic frequency in the frequency range of 5 Hz-1/2 of the frequency conversion range of the pressure pulse signal frequency spectrum; detecting whether side bands with difference values of one characteristic frequency in specific characteristic frequencies appear on two sides of She Tong frequency; if a sideband occurs, it is determined that modulation is occurring on both sides of the She Tong frequency. Further, if it is determined that modulation occurs on both sides of the foldback frequency, comparing magnitudes of the specific characteristic frequency and the She Tong frequency in the frequency spectrum of the pressure pulsation signal, and if it is determined that the magnitude of the specific characteristic frequency is greater than the magnitude of She Tong frequency of a preset multiple, it is determined that rotating stall occurs in the compressor at the first operating point, where the corresponding preset multiple may refer to the magnitude of the foldback frequency that is greater than 1/3 times the magnitude of the specific characteristic frequency.
S2, determining a first working condition point with the maximum flow value corresponding to the occurrence of the rotating stall phenomenon as a rotating stall reference point under the current rotating speed frequency according to the rotating stall analysis result.
For the embodiment, in a specific application scenario, since the rotating stall reference point is the working point when the rotating stall phenomenon corresponds to the maximum flow, in order to save operation resources, when the first working point is subjected to unsteady CFD analysis and fourier transform analysis to obtain the rotating stall analysis result, the first working point can be sequentially subjected to unsteady CFD simulation according to the sequence from the large flow value to the small flow value, and after the first working point with the rotating stall is determined, the first working point can be determined as the rotating stall reference point, so that the rotating stall analysis of the first working point with the subsequent smaller flow value is not required to be performed again, the operation workload can be effectively reduced, and the extraction efficiency of the rotating stall reference point is improved.
For example, based on the example of embodiment step 203, seven first operating points are extracted, respectively: 97%m,98%m,99%m, m,101%m,102%m and 103%m, so that when the unsteady CFD analysis and the Fourier transform analysis are carried out, the analysis can be sequentially carried out according to the order from the large flow value to the small flow value corresponding to the first working point, namely, the first working point 103%m is firstly analyzed, if the first working point 103%m is judged to be the rotating stall of the compressor, the first working point 103%m can be directly determined as the rotating stall reference point under the current rotating speed frequency, and the rotating stall analysis of the first working point with other smaller flow values is not needed to be continuously carried out; in contrast, if it is determined that the compressor does not stall at the first operating point 103%m, further analysis of the rotating stall at the first operating point 102%m is required, and if it is determined that the compressor does not stall at the first operating point 102%m as well, analysis of the rotating stall condition … … at the first operating point 101%m is required, and after the rotating stall analysis is performed on the first operating points one by one in the order of the flow values from the high to the low, the first operating point where the rotating stall phenomenon occurs is determined as a rotating stall reference point at the current rotational frequency.
205. And constructing an early warning datum line by utilizing the rotating stall datum points corresponding to the rotating speed frequencies.
For this embodiment, after each rotational speed frequency passes through the embodiment steps 202-204, the respective rotating stall reference point may be obtained. After the rotating stall datum points corresponding to the rotating speed frequencies are obtained, all the datum points can be connected in a spline curve mode in sequence to further form an alarm datum line, and the alarm datum line is preset in a control system so as to early warn and control the rotating stall phenomenon of the compressor based on the rotating speed frequency and the pressure ratio information in the alarm datum line. It should be noted that, the spline curve can correspond to multiple curve types, and specifically can be selected in a replacement manner according to an actual application scene, and in the scheme, the spline curve can be preferably a B-spline curve.
206. And acquiring operation information of the compressor in real time, and judging that the compressor has rotating stall if the rotation speed frequency and the pressure ratio value in the operation information are determined to trigger an early warning datum line.
For the embodiment, in a specific application scenario, the rotational speed frequency and the pressure ratio of each rotating stall reference point are marked in the alarm reference line, after the current operation information of the compressor is collected, the rotational speed frequency and the pressure ratio of the current working condition point in the operation information are further compared with the marked data of each rotating stall reference point, if the existence of the rotating stall reference point overlapped with the current working condition point is determined, the triggering early warning reference line is determined, namely, the existence of the rotating stall of the compressor is identified, and further early warning information can be output.
In a specific application scenario, an example flow chart of the identification of the rotating stall of the compressor is shown in fig. 3, and the flow chart of the identification of the rotating stall of the compressor can be firstly analyzed by a constant CFD corresponding to each rotating speed frequency of the compressor stage, so as to obtain a flow-pressure ratio characteristic curve corresponding to each rotating speed frequency, a possible stall point with the maximum pressure ratio is extracted through the flow-pressure ratio characteristic curve, 3 working condition points with the flow difference of 1% are respectively taken before and after the pressure ratio value according to the pressure ratio epsilon and the mass flow m of the possible stall point, namely, 97% m,98% m,99% m,101% m,102% m and 103% m. And sequentially performing unsteady CFD simulation on the seven points according to the sequence of the mass flow from high to low, namely firstly performing unsteady CFD simulation on 103% m working points to obtain pressure signals of 8192 time steps corresponding to the working points. Further fast fourier transformation is performed on the pressure signal to obtain a frequency domain result of the pressure pulsation, and She Tong frequency (Blade Passing Frequency) of the impeller operation can be identified. Detecting whether modulation occurs at two sides of She Tong frequency by carrying out band-pass filtering and envelope demodulation analysis on She Tong frequency in the frequency spectrum of the pressure pulsation signal; if yes, comparing the specific characteristic frequency in the low frequency range with the modulation frequency; if the amplitude of the specific characteristic frequency is larger than the She Tong frequency amplitude of the preset multiple, determining that the compressor is in rotating stall, further extracting a rotating stall reference point corresponding to each rotating speed frequency, and realizing identification and control of the rotating stall of the compressor by constructing an alarm reference line.
By the identification method of the rotating stall of the compressor, each rotating speed frequency of the compressor impeller in the working state can be obtained first, and further, the characteristic curve of the compressor pressure ratio changing along with the flow is obtained by respectively carrying out constant CFD analysis on each rotating speed frequency. Based on the flow-pressure ratio characteristic curve, a preset number of working condition points which are likely to generate rotating stall can be extracted, and the rotating stall datum point corresponding to the maximum flow of the rotating stall phenomenon is extracted through unsteady CFD analysis and Fourier transform analysis on the extracted working condition points. After the rotating stall datum points corresponding to the rotating speed frequencies are extracted, an early warning datum line can be constructed by utilizing the rotating stall datum points, and the early warning datum line is preset in an alarm control system, so that when the early warning datum line is triggered or is about to be triggered in actual operation of the compressor, an alarm is automatically started, personnel in a control room are prompted to take corresponding measures, and therefore rotating stall of the compressor is restrained. In the technical scheme, the CFD simulation method applied to the computer technology can accurately identify the early warning datum points corresponding to each rotating speed frequency, so that the early warning datum points are created, the effect of effectively early warning the rotating stall phenomenon through the early warning datum points when the rotating stall phenomenon is caused at different rotating speed frequencies is realized, and the damage of the rotating stall phenomenon to the compressor is avoided.
Further, as a specific implementation of the method shown in fig. 1 and fig. 2, an embodiment of the present application provides a device for identifying rotating stall of a compressor, as shown in fig. 4, where the device includes: the device comprises an acquisition module 31, an analysis module 32, a determination module 33 and an identification module 34;
an acquisition module 31, configured to acquire each rotational frequency at which the compressor wheel operates;
the analysis module 32 is configured to perform a constant CFD analysis on each rotational frequency to obtain each flow-pressure ratio characteristic curve;
a determining module 33, configured to determine a rotating stall reference point corresponding to each rotational speed frequency based on the flow-pressure ratio characteristic curve and through unsteady CFD analysis and fourier transform analysis;
the identification module 34 may be configured to construct an early warning baseline using the rotating stall reference point, and identify the rotating stall of the compressor based on the early warning baseline.
In a specific application scenario, in order to determine the rotating stall reference point corresponding to each rotational speed frequency, as shown in fig. 5, the determining module 33 may specifically include: an extraction unit 331, a determination unit 332;
the extracting unit 331 is configured to extract a preset number of first operating points meeting a preset condition based on the flow-pressure ratio characteristic curve;
A determining unit 332, configured to determine a rotating stall reference point from the first operating point by performing unsteady CFD analysis and fourier transform analysis on the first operating point;
correspondingly, in order to extract a preset number of first working condition points meeting preset conditions, the extracting unit 331 may be specifically configured to extract, in the flow-pressure ratio characteristic curve, a first target working condition point corresponding to a maximum pressure ratio value and each second target working condition point spaced from the first target working condition point by a preset flow interval; and determining the first target operating point and the second target operating point as the first operating point.
In a specific application scenario, in order to determine the rotating stall reference point from the first working point, the determining unit 332 may be specifically configured to perform unsteady CFD analysis and fourier transform analysis on the first working point, to obtain a rotating stall analysis result; and according to the rotating stall analysis result, determining a first working condition point with the maximum flow value, corresponding to the occurrence of the rotating stall phenomenon, as a rotating stall datum point under the current rotating speed frequency.
Correspondingly, in order to obtain the rotating stall analysis result by performing unsteady CFD analysis and fourier transform analysis on the first operating point, the determining unit 332 may be specifically configured to sequentially perform unsteady CFD simulation on the first operating point according to the order from the large flow value to the small flow value, so as to obtain a pressure signal corresponding to a preset time step of the first operating point; performing Fourier transform analysis on the pressure signal to obtain a frequency spectrum of the pressure pulse signal; according to She Tong frequency in the frequency spectrum of the pressure pulse signal, carrying out band-pass filtering and envelope demodulation analysis, and detecting whether modulation occurs at two sides of She Tong frequency; if the modulation phenomenon appears on the two sides of the leaf-pass frequency, comparing the specific characteristic frequency with the amplitude of She Tong frequency in the frequency spectrum of the pressure pulsation signal; if the amplitude of the specific characteristic frequency is determined to be larger than the amplitude of She Tong frequency of the preset multiple, judging that the compressor is in rotating stall at the first working point; and if the modulation phenomenon does not occur on the two sides of the blade pass frequency, judging that the rotating stall of the compressor does not occur at the first working condition point.
In a specific application scenario, in order to sequentially perform unsteady CFD simulation on the first working condition points according to the order of the flow values from large to small, to obtain a pressure signal of the first working condition point corresponding to a preset time step, the determining unit 332 may be specifically configured to perform iterative convergence calculation on the first working condition point according to a preset phase parameter of the unsteady CFD simulation; and if the unsteady CFD simulation is judged to reach a convergence state, extracting a pressure value of the first working condition point corresponding to a preset time step as a pressure signal.
Correspondingly, in order to perform bandpass filtering and envelope demodulation analysis according to the She Tong frequency in the pressure pulse signal spectrum, detect whether modulation occurs on two sides of She Tong frequency, and determine the unit 332, which can be specifically used to extract a specific characteristic frequency in the frequency range of 5 Hz-1/2 of the pressure pulse signal spectrum; detecting whether side bands with difference values of one characteristic frequency in specific characteristic frequencies appear on two sides of She Tong frequency; if a sideband occurs, it is determined that modulation is occurring on both sides of the She Tong frequency.
In a specific application scenario, in order to construct an early warning reference line by using the rotating stall reference point to realize the identification of the rotating stall of the compressor, as shown in fig. 5, the identification module 34 may specifically include: a construction unit 341 and a determination unit 342;
A construction unit 341, configured to construct an early warning reference line by using the rotating stall reference points corresponding to the rotational speed frequencies;
the construction unit 342 may be configured to collect operation information of the compressor in real time, and determine that the compressor has rotating stall if it is determined that the rotation speed frequency and the pressure ratio in the operation information trigger the early warning reference line.
It should be noted that, other corresponding descriptions of each functional unit related to the identification device for rotating stall of a compressor provided by the embodiment of the present invention may refer to corresponding descriptions in fig. 1 and fig. 2, and are not repeated herein.
Based on the above-described methods as shown in fig. 1 and 2, correspondingly, the embodiment of the present invention further provides a computer readable storage medium having a computer program stored thereon, the program when executed by a processor implementing the steps of: acquiring each rotating speed frequency of the operation of the compressor impeller; performing constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve; based on the flow-pressure ratio characteristic curve, determining a rotating stall datum point corresponding to each rotating speed frequency through unsteady CFD analysis and Fourier transform analysis; and constructing an early warning datum line by utilizing the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line.
Based on the above embodiment of the method shown in fig. 1 and the apparatus shown in fig. 4, the embodiment of the present invention further provides a computer device, as shown in fig. 6, including a processor (processor) 41, a communication interface (Communications Interface) 42, a memory (memory) 43, and a communication bus 44. Wherein: processor 41, communication interface 42, and memory 43 communicate with each other via a communication bus 44. A communication interface 44 for communicating with network elements of other devices, such as clients or other servers. The processor 41 is configured to execute a program, and may specifically perform the relevant steps in the above-described embodiment of the method for identifying rotating stall of a compressor. In particular, the program may include program code including computer-operating instructions. The processor 41 may be a central processing unit CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present invention.
The one or more processors included in the terminal may be the same type of processor, such as one or more CPUs; but may also be different types of processors such as one or more CPUs and one or more ASICs. A memory 43 for storing programs. The memory 43 may comprise a high-speed RAM memory or may further comprise a non-volatile memory (non-volatile memory), such as at least one disk memory. The program may be specifically for causing the processor 41 to: acquiring each rotating speed frequency of the operation of the compressor impeller; performing constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve; based on the flow-pressure ratio characteristic curve, determining a rotating stall datum point corresponding to each rotating speed frequency through unsteady CFD analysis and Fourier transform analysis; and constructing an early warning datum line by utilizing the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line.
The identification method, the identification device and the computer equipment for the rotating stall of the compressor can firstly acquire each rotating speed frequency of the compressor impeller in the working state, and further acquire the characteristic curve of the compressor pressure ratio changing along with the flow by respectively carrying out constant CFD analysis on each rotating speed frequency. Based on the flow-pressure ratio characteristic curve, a preset number of working condition points which are likely to generate rotating stall can be extracted, and the rotating stall datum point corresponding to the maximum flow of the rotating stall phenomenon is extracted through unsteady CFD analysis and Fourier transform analysis on the extracted working condition points. After the rotating stall datum points corresponding to the rotating speed frequencies are extracted, an early warning datum line can be constructed by utilizing the rotating stall datum points, and the early warning datum line is preset in an alarm control system, so that when the early warning datum line is triggered or is about to be triggered in actual operation of the compressor, an alarm is automatically started, personnel in a control room are prompted to take corresponding measures, and therefore rotating stall of the compressor is restrained. Through the technical scheme in this application, can apply the CFD analog method of computer technology, the early warning datum point that each rotational speed frequency corresponds is accurately discerned to establish the early warning datum line through the early warning datum point, in order to realize when different rotational speed frequencies, carry out effective early warning through the early warning datum line rotatory stall phenomenon, thereby avoid the harm that rotatory stall phenomenon caused to the compressor.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and for parts of one embodiment that are not described in detail, reference may be made to related descriptions of other embodiments.
It will be appreciated that the relevant features of the methods and apparatus described above may be referenced to one another. In addition, the "first", "second", and the like in the above embodiments are for distinguishing the embodiments, and do not represent the merits and merits of the embodiments.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.
The algorithms and displays presented herein are not inherently related to any particular computer, virtual system, or other apparatus. Various general-purpose systems may also be used with the teachings herein. The required structure for a construction of such a system is apparent from the description above. In addition, the present invention is not directed to any particular programming language. It will be appreciated that the teachings of the present invention described herein may be implemented in a variety of programming languages, and the above description of specific languages is provided for disclosure of enablement and best mode of the present invention.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the invention, various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. However, the disclosed method should not be construed as reflecting the intention that: i.e., the claimed invention requires more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Those skilled in the art will appreciate that the modules in the apparatus of the embodiments may be adaptively changed and disposed in one or more apparatuses different from the embodiments. The modules or units or components of the embodiments may be combined into one module or unit or component and, furthermore, they may be divided into a plurality of sub-modules or sub-units or sub-components. Any combination of all features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or units of any method or apparatus so disclosed, may be used in combination, except insofar as at least some of such features and/or processes or units are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings), may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
Furthermore, those skilled in the art will appreciate that while some embodiments herein include some features but not others included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments can be used in any combination.
Various component embodiments of the invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that some or all of the functions of some or all of the components in a subject seating chart generation apparatus in accordance with an embodiment of the present invention may be implemented in practice using a microprocessor or Digital Signal Processor (DSP). The present invention can also be implemented as an apparatus or device program (e.g., a computer program and a computer program product) for performing a portion or all of the methods described herein. Such a program embodying the present invention may be stored on a computer readable medium, or may have the form of one or more signals. Such signals may be downloaded from an internet website, provided on a carrier signal, or provided in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the words first, second, third, etc. do not denote any order. These words may be interpreted as names.

Claims (7)

1. A method of identifying rotating stall of a compressor, comprising:
acquiring each rotating speed frequency of the operation of the compressor impeller;
performing constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve;
Extracting a first target working condition point corresponding to the maximum pressure ratio and each second target working condition point which is separated from the first target working condition point by a preset flow interval from the first target working condition point in the flow-pressure ratio characteristic curve; determining the first target operating point and the second target operating point as first operating points; sequentially performing unsteady CFD simulation on the first working condition point according to the sequence from the high flow value to the low flow value to obtain a pressure signal of the first working condition point corresponding to a preset time step; performing Fourier transform analysis on the pressure signal to obtain a pressure pulse signal frequency spectrum; according to She Tong frequency in the frequency spectrum of the pressure pulse signal, carrying out band-pass filtering and envelope demodulation analysis, and detecting whether modulation phenomenon occurs at two sides of the She Tong frequency; if the modulation phenomenon is judged to occur on two sides of the She Tong frequency, comparing the specific characteristic frequency with the amplitude of the She Tong frequency in the frequency spectrum of the pressure pulse signal; if the amplitude of the specific characteristic frequency is determined to be larger than the amplitude of the She Tong frequency of the preset multiple, determining that rotating stall of the compressor occurs at the first working condition point; if the modulation phenomenon is not judged to occur on the two sides of the She Tong frequency, judging that the rotating stall of the compressor does not occur at the first working condition point; according to the rotating stall analysis result, determining a first working condition point with the maximum flow value corresponding to the occurrence of the rotating stall phenomenon as a rotating stall datum point under the current rotating speed frequency;
And constructing an early warning datum line by using the rotating stall datum point, and identifying the rotating stall of the compressor based on the early warning datum line.
2. The method of claim 1, wherein the step of sequentially performing unsteady CFD simulation on the first operating point in order of flow values from high to low to obtain a pressure signal of the first operating point corresponding to a preset time step specifically includes:
performing iterative convergence calculation on the first working condition point according to preset phase parameters of unsteady CFD simulation;
and if the unsteady CFD simulation is judged to reach a convergence state, extracting a pressure value of the first working condition point corresponding to a preset time step as a pressure signal.
3. The method according to claim 2, wherein the detecting whether modulation occurs on both sides of the She Tong frequency by performing bandpass filtering and envelope demodulation analysis according to She Tong frequencies in the frequency spectrum of the pressure pulse signal comprises:
extracting specific characteristic frequency within a frequency range of 5 Hz-1/2 of the frequency conversion range of the pressure pulse signal frequency spectrum;
detecting whether side bands with difference values of one characteristic frequency in the specific characteristic frequencies appear on two sides of the She Tong frequency;
If the side band appears, the modulation phenomenon appears on two sides of the She Tong frequency.
4. The method of claim 3, wherein said constructing an early warning reference line using said rotating stall reference point, identifying rotating stall of said compressor based on said early warning reference line, comprises:
constructing an early warning datum line by utilizing the rotating stall datum points corresponding to the rotating speed frequencies;
and acquiring the operation information of the compressor in real time, and judging that the compressor is rotating stall if the rotation speed frequency and the pressure ratio value in the operation information are determined to trigger the early warning datum line.
5. An apparatus for identifying rotating stall of a compressor, comprising:
the acquisition module is used for acquiring each rotating speed frequency of the operation of the compressor impeller;
the analysis module is used for respectively carrying out constant CFD analysis on each rotating speed frequency to obtain each flow-pressure ratio characteristic curve;
the determining module is used for extracting a first target working condition point corresponding to the maximum pressure ratio and each second target working condition point which is separated from the first target working condition point by a preset flow interval from the first target working condition point in the flow-pressure ratio characteristic curve; determining the first target operating point and the second target operating point as first operating points; sequentially performing unsteady CFD simulation on the first working condition point according to the sequence from the high flow value to the low flow value to obtain a pressure signal of the first working condition point corresponding to a preset time step; performing Fourier transform analysis on the pressure signal to obtain a pressure pulse signal frequency spectrum; according to She Tong frequency in the frequency spectrum of the pressure pulse signal, carrying out band-pass filtering and envelope demodulation analysis, and detecting whether modulation phenomenon occurs at two sides of the She Tong frequency; if the modulation phenomenon is judged to occur on two sides of the She Tong frequency, comparing the specific characteristic frequency with the amplitude of the She Tong frequency in the frequency spectrum of the pressure pulse signal; if the amplitude of the specific characteristic frequency is determined to be larger than the amplitude of the She Tong frequency of the preset multiple, determining that rotating stall of the compressor occurs at the first working condition point; if the modulation phenomenon is not judged to occur on the two sides of the She Tong frequency, judging that the rotating stall of the compressor does not occur at the first working condition point; according to the rotating stall analysis result, determining a first working condition point with the maximum flow value corresponding to the occurrence of the rotating stall phenomenon as a rotating stall datum point under the current rotating speed frequency;
And the identification module is used for constructing an early warning datum line by utilizing the rotating stall datum point and identifying the rotating stall of the compressor based on the early warning datum line.
6. A non-transitory readable storage medium having stored thereon a computer program, wherein the program when executed by a processor implements the method of identifying rotating stall of a compressor of any of claims 1 to 4.
7. A computer device comprising a non-volatile readable storage medium, a processor and a computer program stored on the non-volatile readable storage medium and executable on the processor, characterized in that the processor implements the method of identifying rotating stall of a compressor according to any of claims 1 to 4 when executing the program.
CN202011095061.2A 2020-10-14 2020-10-14 Identification method and device for rotating stall of compressor and computer equipment Active CN112417774B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011095061.2A CN112417774B (en) 2020-10-14 2020-10-14 Identification method and device for rotating stall of compressor and computer equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011095061.2A CN112417774B (en) 2020-10-14 2020-10-14 Identification method and device for rotating stall of compressor and computer equipment

Publications (2)

Publication Number Publication Date
CN112417774A CN112417774A (en) 2021-02-26
CN112417774B true CN112417774B (en) 2024-02-23

Family

ID=74854455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011095061.2A Active CN112417774B (en) 2020-10-14 2020-10-14 Identification method and device for rotating stall of compressor and computer equipment

Country Status (1)

Country Link
CN (1) CN112417774B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117688696B (en) * 2024-02-02 2024-05-24 陕西空天信息技术有限公司 Stall assessment method for radial flow centrifugal impeller

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106321465A (en) * 2016-08-22 2017-01-11 江苏大学 Kinetics characterization method for rotating stall of centrifugal pump impeller
CN107165850A (en) * 2017-06-27 2017-09-15 西北工业大学 A kind of rotating stall of axial flow compressor method for early warning recognized based on frequency domain hump
CN109030849A (en) * 2018-05-22 2018-12-18 大连透平机械技术发展有限公司 The monitoring method and device of compressor rotating stall
CN110953178A (en) * 2019-12-06 2020-04-03 华北电力科学研究院有限责任公司 Stall early warning method and system for parallel running fans
CN110987404A (en) * 2019-12-23 2020-04-10 安徽容知日新科技股份有限公司 Rotating speed identification method and computing equipment

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106321465A (en) * 2016-08-22 2017-01-11 江苏大学 Kinetics characterization method for rotating stall of centrifugal pump impeller
CN107165850A (en) * 2017-06-27 2017-09-15 西北工业大学 A kind of rotating stall of axial flow compressor method for early warning recognized based on frequency domain hump
CN109030849A (en) * 2018-05-22 2018-12-18 大连透平机械技术发展有限公司 The monitoring method and device of compressor rotating stall
CN110953178A (en) * 2019-12-06 2020-04-03 华北电力科学研究院有限责任公司 Stall early warning method and system for parallel running fans
CN110987404A (en) * 2019-12-23 2020-04-10 安徽容知日新科技股份有限公司 Rotating speed identification method and computing equipment

Also Published As

Publication number Publication date
CN112417774A (en) 2021-02-26

Similar Documents

Publication Publication Date Title
CN106769033B (en) Variable speed rolling bearing fault recognition methods based on order envelope time-frequency energy spectrum
EP4123284A1 (en) Fault diagnosis method and apparatus therefor
CN112417774B (en) Identification method and device for rotating stall of compressor and computer equipment
CN109765055B (en) Rolling bearing fault detection method and system based on EWT, spectrum effective value and KNN
KR20220116051A (en) Method and apparatus for inspecting wind turbine blades, and apparatus and storage medium therefor
CN104483602A (en) Local discharge signal identification method and device
US20220283556A1 (en) Sensor based data set method and apparatus
CN112131907A (en) Method and device for training classification model
CN109030849B (en) Method and device for monitoring rotating stall of compressor
CN103375421A (en) System and method for monitoring the health of stator vanes
US20190332102A1 (en) Machine health monitoring of rotating machinery
CN112782421B (en) Audio-based rotating speed identification method
CN115788936A (en) Aerodynamic force analysis method of vaneless diffuser
CN112881044A (en) Abnormality diagnosis method for variable speed device
CN112903296B (en) Rolling bearing fault detection method and system
CN112462179A (en) Narrow band-pass filtering method, device, equipment and storage medium
Xie et al. Fan bearing fault diagnosis based on continuous wavelet transform and autocorrelation
CN113542200A (en) Risk control method, risk control device and storage medium
CN205449521U (en) Range hood failure diagnosis's system based on EMD falls and makes an uproar
CN112733351B (en) Method and device for processing axial force data of impeller of centrifugal compressor
CN116242612B (en) Fault diagnosis method, device, medium and equipment
Lin et al. Experimental investigation of unsteady rotor tip flows in a high speed compressor throttled to stall
US20230235743A1 (en) Incipient compressor surge detection using artificial intelligence
CN113177514B (en) Unmanned aerial vehicle signal detection method and device and computer readable storage medium
CN109695579B (en) Early warning method and device for wind driven generator

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