Python-based Bode diagram data processing method and related product

Abstract

The invention discloses a Python-based Bode diagram data processing method and related products, which can be applied to the technical field of data processing, and the method comprises the following steps: loading rotating speed data and vibration data in the starting process of a rotating machine; analyzing the rotating speed data, and obtaining corresponding average rotating frequency and average rotating speed; analyzing the vibration data, extracting an amplitude and a phase of a frequency spectrum component corresponding to the average rotation frequency, and obtaining a vibration vector corresponding to the average rotation speed; and drawing a Bode diagram based on the average rotating speed and the vibration vector. Thus, on-site workers can process and draw the Bode diagram data in the starting and stopping process of the rotating machine only through the Python script without the help of professional software tools, and the flexibility of data processing is improved.

Description

Technical Field

[0001] The present application relates to the field of data processing technology, and in particular to a Python-based Bode diagram data processing method and related products. Background Art

[0002] The Bode diagram is a graphical representation method of linear discontinuous systems in frequency domain analysis. As a common tool for vibration analysis of rotating machinery, it displays the vibration data of the machine start-up or stop process. It can not only be used to determine the critical speed of the shaft system, the amplitude and phase at the critical speed, but also to determine the axial position and imbalance form of the unbalanced mass in the shaft system.

[0003] The existing Bode diagram data processing and drawing methods require the use of professional vibration analysis software, which has the problems of low versatility and inability to process third-party vibration data.

[0004] Therefore, how to improve the flexibility of Bode diagram data processing is an urgent problem to be solved by those skilled in the art. Summary of the invention

[0005] Based on the above problems, the present application provides a Python-based Bode diagram data processing method and related products, which can realize the processing and drawing of Bode diagram data during the start-up and shutdown of rotating machinery only through Python scripts, thereby improving the flexibility of data processing.

[0006] In a first aspect, an embodiment of the present application provides a Bode diagram data processing method based on Python, comprising:

[0007] Load the speed data and vibration data of the rotating machinery during startup;

[0008] Analyze the rotation speed data and obtain the corresponding average rotation frequency and average rotation speed;

[0009] Analyzing the vibration data, extracting the amplitude and phase of the frequency spectrum component corresponding to the average rotation frequency, and obtaining the vibration vector corresponding to the average rotation speed;

[0010] A Bode diagram is plotted based on the average rotation speed and the vibration vector.

[0011] Optionally, before loading the speed data and vibration data during the startup of the rotating machinery, the method further includes:

[0012] The speed data and vibration data of the rotating machinery during the startup process are acquired by sampling with sensors, and corresponding data groups are formed in chronological order and stored in a memory;

[0013] The sampling frequency and sampling time of the rotation speed data and the vibration data are consistent;

[0014] The data group includes at least one data.

[0015] Optionally, parsing the rotation speed data and obtaining a corresponding average rotation frequency includes:

[0016] Set detection threshold;

[0017] Based on the detection threshold, detecting a pulse position in the rotation speed data;

[0018] The average rotational frequency corresponding to the rotational speed data is determined by combining the time corresponding to each pulse position and the sampling frequency.

[0019] Optionally, the detecting the pulse position in the rotation speed data based on the detection threshold comprises:

[0020] Check whether the elements at index 1 and later in the data group corresponding to the rotation speed data are greater than or equal to the detection threshold, and construct a first Boolean array;

[0021] Check whether the index 0 to the second to last element in the data group corresponding to the rotation speed data is less than the detection threshold, and construct a second Boolean array;

[0022] Performing an AND operation on the first Boolean array and the second Boolean array to obtain a third Boolean array;

[0023] The index of the element in the third Boolean array whose value is true is determined and returned as the pulse position in the speed data.

[0024] Optionally, before analyzing the vibration data, the method further includes:

[0025] According to the pulse position in the speed data, the vibration data is subjected to full-cycle truncation processing;

[0026] Performing a de-averaging operation on the vibration data after the whole cycle truncation processing;

[0027] A windowing operation is performed on the vibration data after the mean removal operation, and a correction coefficient is assigned.

[0028] Optionally, the method further includes:

[0029] If there is a problem with the rotor of the rotating machine, a slow rolling vector corresponding to the average speed being within a preset slow rolling speed range is calculated;

[0030] Compensating the vibration vector based on the slow rolling vector to obtain a compensated vibration vector;

[0031] The step of drawing a Bode diagram based on the average rotation speed and the vibration vector comprises:

[0032] A Bode diagram is plotted using the average rotation speed and the compensated vibration vector.

[0033] Optionally, compensating the vibration vector based on the slow rolling vector includes:

[0034] Among the vibration vectors, the vibration vector having a value higher than the slow rolling vector is subtracted from the slow rolling vector.

[0035] In a second aspect, an embodiment of the present application provides a Bode diagram data processing device based on Python, comprising:

[0036] A loading module is used to load the speed data and vibration data during the starting process of the rotating machinery;

[0037] A first analysis module, used for analyzing the rotation speed data and obtaining the corresponding average rotation frequency and average rotation speed;

[0038] A second analysis module, used to analyze the vibration data, extract the amplitude and phase of the frequency spectrum component corresponding to the average rotation frequency, and obtain the vibration vector corresponding to the average rotation speed;

[0039] A drawing module is used to draw a Bode diagram based on the average rotation speed and the vibration vector.

[0040] In a third aspect, an embodiment of the present application provides a Bode diagram data processing device based on Python, comprising:

[0041] Memory for storing computer programs;

[0042] A processor is used to implement the steps of the Bode diagram data processing method based on Python as described above when executing the computer program.

[0043] In a fourth aspect, an embodiment of the present application provides a readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the Python-based Bode diagram data processing method as described in claims 1 to 7 are implemented.

[0044] It can be seen from the above technical solutions that compared with the prior art, the present application has the following advantages:

[0045] This application first loads the speed data and vibration data during the startup of the rotating machinery. Then parse the speed data and obtain the corresponding average rotational frequency and average speed. Then parse the vibration data, extract the amplitude and phase of the spectral component corresponding to the average rotational frequency, and obtain the vibration vector corresponding to the average speed. Finally, draw the Bode diagram based on the average speed and vibration vector. In this way, on-site staff can process and draw the Bode diagram data during the start-up and shutdown of the rotating machinery through Python scripts without the help of professional software tools, thereby improving the flexibility of data processing. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 A flowchart of a Bode diagram data processing method based on Python provided in an embodiment of the present application;

[0047] Figure 2 An ideal rotor vibration Bode diagram provided in an embodiment of the present application;

[0048] Figure 3 A schematic diagram of a speed increase theoretical curve provided in an embodiment of the present application;

[0049] Figure 4 A schematic diagram of an ideal rotation speed data provided in an embodiment of the present application;

[0050] Figure 5 A flow chart of Bode diagram data processing provided in an embodiment of the present application;

[0051] Figure 6 A schematic diagram of the structure of a Bode diagram data processing device based on Python provided in an embodiment of the present application. DETAILED DESCRIPTION

[0052] As mentioned above, the existing Bode diagram data processing and drawing methods have the problem of insufficient flexibility. Specifically, the existing Bode diagram data processing and drawing methods require the use of specific instruments and professional vibration analysis software, which are not very versatile and cannot process third-party vibration data. In addition, they are inefficient when processing large amounts of data, so there is a problem of insufficient flexibility in Bode diagram data processing.

[0053] In order to solve the above problems, the embodiment of the present application provides a Bode diagram data processing method based on Python, including: first loading the speed data and vibration data during the startup process of the rotating machinery. Then parsing the speed data, and obtaining the corresponding average rotation frequency and average rotation speed. Then parsing the vibration data, extracting the amplitude and phase of the spectral component corresponding to the average rotation frequency, and obtaining the vibration vector corresponding to the average rotation speed. Finally, drawing the Bode diagram based on the average rotation speed and the vibration vector.

[0054] In this way, on-site staff can process and draw Bode diagram data during the start-up and shutdown of rotating machinery through Python scripts without the need for professional software tools, thereby improving the flexibility of data processing.

[0055] It should be noted that the Python-based Bode diagram data processing method and related products provided in this application can be applied to the technical field of power terminal equipment. The above is only an example and does not limit the application field of the Python-based Bode diagram data processing method and related products provided in this application.

[0056] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0057] Figure 1 A flowchart of a Bode diagram data processing method based on Python provided in an embodiment of the present application. Figure 1 As shown, a Bode diagram data processing method based on Python provided in an embodiment of the present application may include:

[0058] S101: Loading speed data and vibration data during the startup of the rotating machinery.

[0059] In practical applications, Python, as a powerful and easy-to-learn programming language, is becoming more and more widely used in the fields of scientific computing and data analysis. It has a wealth of scientific computing libraries and data processing tools, such as NumPy, SciPy, and Matplotlib, which provide a convenient way for data processing and visualization of Bode diagrams, and also improve the flexibility of data processing. Rotating machinery refers to machinery that mainly relies on rotation to complete specific functions. The main function of rotating machinery is completed by rotating parts, generally rotors. Typical rotating machinery includes steam turbines, gas turbines, centrifugal and axial compressors, fans, etc. Figure 2 An ideal rotor vibration Bode diagram provided in an embodiment of the present application. Figure 2As shown, the Bode diagram consists of two rectangular coordinate systems, and the two coordinate diagrams will have vibration amplitude and phase respectively, and the horizontal axis of the rectangular coordinate diagram represents speed or frequency. It can be used for balance analysis, resonance identification and resonance analysis of rotating machinery. Specifically, in the data processing and drawing process of the Bode diagram in the embodiment of the present application, it is first necessary to load the collected speed data and vibration data of the start / stop process of the rotating machinery from the memory based on the Python program.

[0060] In addition, since the methods for acquiring rotation speed data and vibration data are different, the embodiment of the present application can illustrate a possible acquisition method.

[0061] In one embodiment, before loading the speed data and vibration data during the startup of the rotating machinery, the method further includes:

[0062] The speed data and vibration data of the rotating machinery during the startup process are acquired by sampling with sensors, and corresponding data groups are formed in chronological order and stored in a memory;

[0063] The sampling frequency and sampling time of the rotation speed data and the vibration data are consistent;

[0064] The data group includes at least one data.

[0065] In practical applications, the speed data and vibration data of rotating machinery are generally collected by eddy current sensors or photoelectric sensors, and then the collected data are stored in a memory. It should be noted that the speed data and vibration data in the embodiment of the present application must be collected by synchronous sampling, that is, the sampling frequency and sampling time of the two data must be consistent. Figure 3 A schematic diagram of a speed increase theoretical curve provided in an embodiment of the present application. Figure 3 As shown, taking the start-up process of a steam turbine as an example, its rated speed is 1500r / min, and the acceleration rate is 150r / min, then the frequency change per second is 150 / 60 / 60=0.042Hz, and the time required for the speed to change from 0 to 1500 is 10 minutes. If the sensor acquisition interval (sampling frequency) is set to 10s to save speed data and vibration data once, and the length of the saved data (sampling time) is 1s, a data group consisting of 60 speed data and a data group consisting of 60 vibration data will be generated in chronological order, and the sampling time of the elements at the same position in the two data groups is the same.

[0066] S102: Analyze the rotation speed data and obtain the corresponding average rotation frequency and average rotation speed.

[0067] Figure 4 A schematic diagram of an ideal speed data provided by an embodiment of the present application. Figure 4As shown, when the speed is 600 rpm, there are 10 complete amplitude values ​​from 0 to 1 per second. In the embodiment of the present application, since the speed data is obtained by interval sampling, the speed data loaded is generally multiple and discontinuous. In addition, since the loaded data is collected during the startup of the rotating machinery, there is a certain rate of increase, resulting in that the various speed data loaded are generally different. Based on the Python program, these speed data are further parsed to obtain the average frequency and average speed corresponding to each speed data. Among them, average speed = average frequency * 60.

[0068] In addition, since the methods for analyzing the rotation speed data are different, the embodiment of the present application can illustrate a possible analysis method.

[0069] In one embodiment, the analyzing the rotation speed data and obtaining the corresponding average rotation frequency includes:

[0070] Set detection threshold;

[0071] Based on the detection threshold, detecting a pulse position in the rotation speed data;

[0072] The average rotational frequency corresponding to the rotational speed data is determined by combining the time corresponding to each pulse position and the sampling frequency.

[0073] In practical applications, combined with the above speed data, it can be known that the amplitude of the speed data varies between 0 and 1, so the detection threshold can be set to 0.5, of course, it can also be set to other reasonable values. Then, based on the Python program and the detection threshold, all the loaded speed data are analyzed to determine all the pulse positions that reach the detection. Among them, the moment corresponding to the pulse position is the moment when one revolution begins, recorded as K = [K0, K1, ..., K n ], assuming that the speed data shows n revolutions. The time difference of the pulse position is the time it takes for the rotor to rotate once, so the average speed frequency f in a set of speed data can be calculated as shown in the following formula: Where f s is the sampling frequency.

[0074] In addition, since the methods for detecting the pulse position are different, the embodiment of the present application can illustrate a possible detection method.

[0075] In one embodiment, the detecting the pulse position in the rotation speed data based on the detection threshold comprises:

[0076] Check whether the elements at index 1 and later in the data group corresponding to the rotation speed data are greater than or equal to the detection threshold, and construct a first Boolean array;

[0077] Check whether the index 0 to the second to last element in the data group corresponding to the rotation speed data is less than the detection threshold, and construct a second Boolean array;

[0078] Performing an AND operation on the first Boolean array and the second Boolean array to obtain a third Boolean array;

[0079] The index of the element in the third Boolean array whose value is true is determined and returned as the pulse position in the speed data.

[0080] In practical applications, continuing with the above example, if the loaded speed data is a data set consisting of sixty elements, then during the parsing process, the Python program is used to check whether the 2nd to 60th elements (elements with index 1 and later) in the data set are greater than or equal to the detection threshold (0.5), and the check results are saved in a Boolean array (the first Boolean array); the 1st to 59th elements (index 0 to the second to last element) in the data set are checked to be less than the detection threshold, and the check results are saved in another Boolean array (the second Boolean array). A Boolean array is an array whose element type is a Boolean type, that is, each element in the array can only be true or false. Then, based on the Python program, the first Boolean array and the second Boolean array above are ANDed to obtain a Boolean array (the third Boolean array) indicating which positions meet the above two conditions at the same time. Finally, the index of the true element in the third Boolean array is returned, and an array is formed to represent all pulse positions that have reached the detection threshold in the speed data.

[0081] S103: Analyze the vibration data, extract the amplitude and phase of the frequency spectrum component corresponding to the average rotation frequency, and obtain the vibration vector corresponding to the average rotation speed.

[0082] In practical applications, the loaded vibration data can be subjected to fast Fourier transform based on the Python program to extract the amplitude and phase of the spectral components corresponding to the average rotation frequency, and then obtain the vibration vector at the average rotation speed.

[0083] In addition, since the preprocessing methods of vibration data are not the same, the embodiment of the present application can illustrate a possible preprocessing method.

[0084] In one embodiment, before analyzing the vibration data, the method further includes:

[0085] According to the pulse position in the speed data, the vibration data is subjected to full-cycle truncation processing;

[0086] Performing a de-averaging operation on the vibration data after the whole cycle truncation processing;

[0087] A windowing operation is performed on the vibration data after the mean removal operation, and a correction coefficient is assigned.

[0088] In practical applications, in order to reduce the error introduced by the algorithm during fast Fourier transform. The embodiment of the present application first pre-processes the loaded vibration data based on the Python program. Specifically, first based on the Python program, the whole cycle of the head and tail of the vibration data is truncated according to the pulse position of the speed data. Then the vibration data after the stage is de-meaned. Finally, the vibration data after the de-meaning operation is windowed. The Python program can provide window functions such as rectangular window, Hamming window, Hanning window, Blackman window, flat-top window, triangular window, etc., and also provides correction coefficients corresponding to the window functions. The default window function is the Hanning window, and the correction coefficient is 2.

[0089] S104: Drawing a Bode diagram based on the average rotation speed and the vibration vector.

[0090] In practical applications, a Bode diagram can be drawn based on the average velocity and vibration vector obtained by analysis, based on the drawing function in the matplotlib package of the Python program.

[0091] In addition, since the data processing methods corresponding to different states of the rotating machinery rotor are not exactly the same, the embodiment of the present application can illustrate a possible processing method.

[0092] In one embodiment, the method further comprises:

[0093] If there is a problem with the rotor of the rotating machine, a slow rolling vector corresponding to the average speed being within a preset slow rolling speed range is calculated;

[0094] Compensating the vibration vector based on the slow rolling vector to obtain a compensated vibration vector;

[0095] The step of drawing a Bode diagram based on the average rotation speed and the vibration vector comprises:

[0096] A Bode diagram is plotted using the average rotation speed and the compensated vibration vector.

[0097] In practical applications, the rotor of the rotating machinery may also have problems, such as bending or connection problems that cause mechanical or electrical oscillation, causing the rotor to generate a slow rolling vector at low speed. The slow rolling vector will distort the dynamic response data of the rotating machinery, which is generally non-dynamic data observed by the sensor. When analyzing rotor vibration, the dynamic response of the rotor is generally analyzed. Therefore, when processing and drawing the Bode diagram data of this type of rotating machinery, slow rolling vector compensation is also required. Specifically, the slow rolling range is generally 10% of the first-order resonance frequency of the rotating machinery. Continuing with the example of the steam turbine with a rated speed of 1500r / min, the slow rolling speed is generally selected at 120r / min. Therefore, based on the analytical results of the above vibration data, the first-order vibration vector of the vibration data with an average speed within the range of 120±5r / min can be selected, and the average is taken as the slow rolling vector. Then, the vibration vector is compensated based on the slow rolling vector, and the compensated vibration vector is obtained. Finally, based on the drawing function in the matplotlib package of the Python program, the Bode diagram can be drawn according to the average speed obtained by analysis and the compensated vibration vector.

[0098] In addition, since the vector compensation methods are not the same, the embodiment of the present application can illustrate a possible compensation method.

[0099] In one embodiment, compensating the vibration vector based on the slow rolling vector includes:

[0100] Among the vibration vectors, the vibration vector having a value higher than the slow rolling vector is subtracted from the slow rolling vector.

[0101] In practical applications, the vibration vector can be compensated by subtracting the slow rolling vector from the vibration vector having a higher speed than the slow rolling vector, thereby obtaining the vibration vector compensated by the slow rolling vector.

[0102] Figure 5 A flow chart of Bode diagram data processing provided in an embodiment of the present application. Figure 5 As shown, the embodiment of the present application is implemented based on Python program. First, the rotation speed data and vibration data are accelerated, and then the rotation speed data and vibration data are analyzed to obtain the average rotation frequency, average rotation speed, the amplitude and phase of the spectrum components corresponding to the average rotation frequency, and the vibration vector corresponding to the average rotation speed. Then, it is judged whether to perform slow rolling vector compensation. If compensation is required, the average rotation speed corresponding to the slow rolling vector is determined according to the rated speed, and then the corresponding slow rolling vector is selected, the slow rolling vector is taken out from the vibration data, and the Bode diagram is drawn; if compensation is not required, the Bode diagram is directly drawn according to the average rotation speed and the vibration vector.

[0103] In summary, the present application first loads the speed data and vibration data during the startup of the rotating machinery. Then the speed data is parsed, and the corresponding average rotational frequency and average rotational speed are obtained. Then the vibration data is parsed, the amplitude and phase of the spectral components corresponding to the average rotational frequency are extracted, and the vibration vector corresponding to the average rotational speed is obtained. Finally, the Bode diagram is drawn based on the average rotational speed and the vibration vector. In this way, on-site staff can process and draw the Bode diagram data during the start-up and shutdown of the rotating machinery through Python scripts without the help of professional software tools, thereby improving the flexibility of data processing.

[0104] Figure 6 A schematic diagram of the structure of a Bode diagram data processing device based on Python provided in an embodiment of the present application. Figure 6 As shown, the Bode diagram data processing device 600 based on Python provided in the embodiment of the present application includes:

[0105] The loading module 601 is used to load the speed data and vibration data during the starting process of the rotating machinery;

[0106] A first analysis module 602 is used to analyze the rotation speed data and obtain the corresponding average rotation frequency and average rotation speed;

[0107] A second analysis module 603 is used to analyze the vibration data, extract the amplitude and phase of the spectrum component corresponding to the average rotation frequency, and obtain the vibration vector corresponding to the average rotation speed;

[0108] The drawing module 604 is configured to draw a Bode diagram based on the average rotation speed and the vibration vector.

[0109] As an implementation method, with respect to how to obtain vibration data and vibration data, the above-mentioned Python-based Bode diagram data processing device 600 includes: an acquisition module;

[0110] An acquisition module is used to acquire the speed data and vibration data of the rotating machinery during the startup process by using sensor sampling, and form corresponding data groups in chronological order and store them in a memory;

[0111] The sampling frequency and sampling time of the rotation speed data and the vibration data are consistent;

[0112] The data group includes at least one data.

[0113] As an implementation manner, regarding how to parse the rotation speed data, the first parsing module 602 includes: a setting module, a detection module and a determination module;

[0114] A setting module, used to set the detection threshold;

[0115] A detection module, configured to detect a pulse position in the rotation speed data based on the detection threshold;

[0116] The determination module is used to determine the average rotational frequency corresponding to the rotational speed data in combination with the time corresponding to each pulse position and the sampling frequency.

[0117] As an implementation method, regarding how to detect the pulse position in the speed data, the above detection module is specifically used to:

[0118] Check whether the elements at index 1 and later in the data group corresponding to the rotation speed data are greater than or equal to the detection threshold, and construct a first Boolean array;

[0119] Check whether the index 0 to the second to last element in the data group corresponding to the rotation speed data is less than the detection threshold, and construct a second Boolean array;

[0120] Performing an AND operation on the first Boolean array and the second Boolean array to obtain a third Boolean array;

[0121] The index of the element in the third Boolean array whose value is true is determined and returned as the pulse position in the speed data.

[0122] As an implementation method, with respect to how to pre-process the vibration data, the above-mentioned Python-based Bode diagram data processing device 600 further includes: a pre-processing module;

[0123] A preprocessing module, used for performing whole-cycle truncation processing on the vibration data according to the pulse position in the speed data;

[0124] Performing a de-averaging operation on the vibration data after the whole cycle truncation processing;

[0125] A windowing operation is performed on the vibration data after the mean removal operation, and a correction coefficient is assigned.

[0126] As an implementation method, regarding how to perform slow rolling vector compensation, the above-mentioned Python-based Bode diagram data processing device 600 further includes: a calculation module and a compensation module;

[0127] A calculation module, for calculating a slow rolling vector corresponding to an average rotation speed within a preset slow rolling speed range if there is a problem with the rotor of the rotating machine;

[0128] The compensation module is used to compensate the vibration vector based on the slow rolling vector and obtain a compensated vibration vector.

[0129] As an implementation manner, regarding how to compensate the vibration vector, the compensation module is specifically used for:

[0130] Among the vibration vectors, the vibration vector having a value higher than the slow rolling vector is subtracted from the slow rolling vector.

[0131] As an implementation manner, with respect to how to draw a Bode diagram, the drawing module 604 is further specifically used for:

[0132] A Bode diagram is plotted using the average rotation speed and the compensated vibration vector.

[0133] In summary, the present application first loads the speed data and vibration data during the startup of the rotating machinery. Then the speed data is parsed, and the corresponding average rotational frequency and average rotational speed are obtained. Then the vibration data is parsed, the amplitude and phase of the spectral components corresponding to the average rotational frequency are extracted, and the vibration vector corresponding to the average rotational speed is obtained. Finally, the Bode diagram is drawn based on the average rotational speed and the vibration vector. In this way, on-site staff can process and draw the Bode diagram data during the start-up and shutdown of the rotating machinery through Python scripts without the help of professional software tools, thereby improving the flexibility of data processing.

[0134] In addition, the present application also provides a Bode diagram data processing device based on Python, including:

[0135] Memory for storing computer programs;

[0136] A processor is used to implement the steps of the Bode diagram data processing method based on Python as described above when executing the computer program.

[0137] In addition, the present application also provides a readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the Python-based Bode diagram data processing method as described above are implemented.

[0138] The above description of the disclosed embodiments enables those skilled in the art to implement or use the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A Bode diagram data processing method based on Python, characterized in that: The method comprises: Load the speed data and vibration data of the rotating machinery during startup; Analyze the rotation speed data and obtain the corresponding average rotation frequency and average rotation speed; Analyzing the vibration data, extracting the amplitude and phase of the frequency spectrum component corresponding to the average rotation frequency, and obtaining the vibration vector corresponding to the average rotation speed; A Bode diagram is plotted based on the average rotation speed and the vibration vector.

2. The method according to claim 1, characterized in that Before loading the speed data and vibration data during the starting process of the rotating machinery, the method further includes: The speed data and vibration data of the rotating machinery during the startup process are acquired by sampling with sensors, and corresponding data groups are formed in chronological order and stored in a memory; The sampling frequency and sampling time of the rotation speed data and the vibration data are consistent; The data group includes at least one data.

3. The method according to claim 2, characterized in that The analyzing the rotation speed data and obtaining the corresponding average rotation frequency includes: Set detection threshold; Based on the detection threshold, detecting a pulse position in the rotation speed data; The average rotational frequency corresponding to the rotational speed data is determined by combining the time corresponding to each pulse position and the sampling frequency.

4. The method according to claim 3, characterized in that The step of detecting the pulse position in the rotation speed data based on the detection threshold comprises: Check whether the elements at index 1 and later in the data group corresponding to the rotation speed data are greater than or equal to the detection threshold, and construct a first Boolean array; Check whether the index 0 to the second to last element in the data group corresponding to the rotation speed data is less than the detection threshold, and construct a second Boolean array; Performing an AND operation on the first Boolean array and the second Boolean array to obtain a third Boolean array; The index of the element in the third Boolean array whose value is true is determined and returned as the pulse position in the speed data.

5. The method according to claim 3, characterized in that: Before analyzing the vibration data, the method further includes: According to the pulse position in the speed data, the vibration data is subjected to full-cycle truncation processing; Performing a de-averaging operation on the vibration data after the whole cycle truncation processing; A windowing operation is performed on the vibration data after the mean removal operation, and a correction coefficient is assigned.

6. The method according to claim 1, characterized in that The method further comprises: If there is a problem with the rotor of the rotating machine, a slow rolling vector corresponding to the average speed being within a preset slow rolling speed range is calculated; Compensating the vibration vector based on the slow rolling vector to obtain a compensated vibration vector; The step of drawing a Bode diagram based on the average rotation speed and the vibration vector comprises: A Bode diagram is plotted using the average rotation speed and the compensated vibration vector.

7. The method according to claim 6, characterized in that The compensating the vibration vector based on the slow rolling vector includes: Among the vibration vectors, the vibration vector having a value higher than the slow rolling vector is subtracted from the slow rolling vector.

8. A Bode diagram data processing device based on Python, characterized in that: include: A loading module is used to load the speed data and vibration data during the starting process of the rotating machinery; A first analysis module, used for analyzing the rotation speed data and obtaining the corresponding average rotation frequency and average rotation speed; A second analysis module, used to analyze the vibration data, extract the amplitude and phase of the frequency spectrum component corresponding to the average rotation frequency, and obtain the vibration vector corresponding to the average rotation speed; A drawing module is used to draw a Bode diagram based on the average rotation speed and the vibration vector.

9. A Bode diagram data processing device based on Python, characterized in that: include: Memory for storing computer programs; A processor, configured to implement the steps of the Bode diagram data processing method based on Python as claimed in any one of claims 1 to 7 when executing the computer program.

10. A readable storage medium, characterized in that: The readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the Python-based Bode diagram data processing method according to any one of claims 1 to 7 are implemented.

Images