CN114485914B - Vibration order tracking analysis method - Google Patents
Vibration order tracking analysis method Download PDFInfo
- Publication number
- CN114485914B CN114485914B CN202111587611.7A CN202111587611A CN114485914B CN 114485914 B CN114485914 B CN 114485914B CN 202111587611 A CN202111587611 A CN 202111587611A CN 114485914 B CN114485914 B CN 114485914B
- Authority
- CN
- China
- Prior art keywords
- vibration
- sequence
- time
- sampling
- rotating speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 238000005070 sampling Methods 0.000 claims abstract description 79
- 238000012952 Resampling Methods 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000000819 phase cycle Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000001360 synchronised effect Effects 0.000 claims abstract description 5
- 238000012544 monitoring process Methods 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000010586 diagram Methods 0.000 description 6
- 239000002699 waste material Substances 0.000 description 5
- 238000007405 data analysis Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The embodiment of the application provides a vibration order tracking analysis method, which comprises the following steps: under the pseudo-synchronous sampling condition, real-time vibration and rotation speed signal acquisition are carried out; adopting a high-frequency counter to perform timing compensation on the rotating speed pulse to acquire an instantaneous rotating speed and an instantaneous phase sequence synchronous with rotating speed sampling; obtaining a sampling time list sequence corresponding to the phase to be interpolated by using an interpolation resampling method according to the rotating speed sampling frequency and the compensated instantaneous phase sequence; and according to the interpolated phase sequence and the corresponding sampling time list, combining the vibration sampling frequency and the time sequence generated by corresponding vibration, and obtaining a new interpolated vibration sequence at the corresponding moment of the interpolated phase sequence by using an interpolation method again. The invention can greatly improve the calculation efficiency of the traditional tracking algorithm.
Description
Technical Field
The embodiment of the application belongs to the technical field of vibration data analysis of rotary machines, and particularly relates to a vibration analysis method under a working condition of variable rotation speed of the rotary machines.
Background
The variable rotation speed working condition contains rich dynamic characteristic information of mechanical products, and is an important content for fault diagnosis and performance evaluation. In order to overcome the problem of vibration data spectrum blurring caused by fluctuation or severe change of rotating speed of a rotating machine, an angle domain resampling algorithm is proposed and solves the problem to a certain extent, and related algorithms are widely applied to vibration data analysis of the rotating machine, and the main stream acquisition and analysis strategies at present are as follows: and determining the sampling frequency according to the highest analysis frequency of the vibration data, synchronously collecting and storing the vibration signal and the rotating speed signal at the same sampling frequency, and finally, reading the data segment by means of analysis software to perform angle domain resampling on the vibration data to obtain a vibration sequence with equal angle intervals.
Thus, the current traditional acquisition and analysis has major drawbacks: the sampling frequency of data acquisition is denoted by Fs, with F r The rotation frequency of the rotation shaft of the tracked shaft is represented, the number of pulses generated by the rotation speed sensor for one rotation of the rotation speed monitoring shaft is represented by N, the transmission ratio between the tracked shaft and the rotation speed monitoring shaft is represented by R (when the tracked shaft and the rotation speed monitoring shaft are coaxial, the R value is 1), and the acquired rotation speed signal per second pulse number M is represented as follows: m=r·f r ·N。
The existing scheme mainly has the following defects:
1. when part of equipment is designed in structural design or scheme, the rotation speed sensor is arranged on a high-speed shaft, the tracked shaft is a low-speed shaft, and at the moment, the rotation speed of the rotation speed monitoring shaft and the tracked shaft are asynchronous, so that a new challenge is brought to the traditional order tracking, and if the equipment transmission ratio R is larger and the rotation frequency F is higher r When the number N of pulses per revolution of the rotating speed sensor is higher or the number N of pulses per second of the rotating speed signal under the corresponding working condition is larger, the number M of pulses per second of the rotating speed signal under the corresponding working condition is synchronously increased, when the sampling frequency Fs is smaller than M, the current sampling strategy can not effectively capture the rotating speed pulses, the resampling of the vibration data angle domain can not be carried out, the actual sampling frequency must be improved for meeting the requirement of the acquisition of the rotating speed signal, the requirements on data acquisition, analysis and storage are higher, the hardware cost is increased, and the resource waste is caused by the excessive sampling frequency of the vibration data.
2. The phase information extracted from the rotating speed signal has a time error of a sampling interval in the time dimension, and for the working condition of high rotating speed, the introduced phase error is larger, and a new acquisition strategy and a phase compensation strategy are required to be provided to improve the phase precision, so that the effective performance of an angle domain resampling algorithm can be ensured.
3. The angle domain resampling is carried out on the vibration data acquired before and after through analysis software, the data segments after the angle domain resampling cannot be spliced smoothly due to the fact that data truncation exists when the phase is acquired, unnecessary frequency components are introduced, and the existing algorithm is not suitable for real-time analysis of the data segments stored in a segmented mode.
Disclosure of Invention
An object of the embodiment of the application is to provide a vibration order tracking analysis method, which adjusts a collection strategy of vibration and rotation speed signals according to actual working condition parameters, avoids the problem of resource waste caused by synchronous increase of sampling frequency on the premise of guaranteeing collection of effective signals, collects the rotation speed signals through a matching board card high-frequency counter, records the value of the counter at each sampling point, and obtains rotation speed trigger information with higher precision through a compensation algorithm, thereby solving the problem in the background technology.
In order to solve the above technical problems, the technical solution of the vibration order tracking analysis method provided in the embodiment of the present application is specifically as follows:
the embodiment of the application discloses a vibration order tracking analysis method, which comprises the following steps:
step 1: under the pseudo-synchronous sampling condition, real-time vibration and rotation speed signal acquisition are carried out;
step 2: adopting a high-frequency counter to perform timing compensation on the rotating speed pulse to acquire an instantaneous rotating speed and an instantaneous phase sequence synchronous with rotating speed sampling;
step 3: obtaining a sampling time list sequence corresponding to the phase to be interpolated by using an interpolation resampling method according to the rotating speed sampling frequency and the compensated instantaneous phase sequence;
step 4: according to the interpolated phase sequence and the corresponding sampling time list, combining the vibration sampling frequency and the time list generated by corresponding vibration, and obtaining a new interpolated vibration sequence at the corresponding moment of the interpolated phase sequence by using an interpolation method again;
step 5: repeating the steps 1 to 5 to complete the continuous order tracking of the signal.
In a preferred embodiment of any of the foregoing solutions, the sampling rates of the real-time vibration and the rotational speed sampling may be inconsistent, and the following relationship needs to be satisfied:
when the device pulses m=r·f per second r N, wherein F r The rotation frequency of the tracking shaft is set, N is the pulse number of each rotation of the rotation speed sensor, and R is the transmission ratio between the tracked shaft and the rotation speed monitoring shaft; vibration sampling frequency F s ＝D·F m (D≥2)，F m Monitoring the highest analysis frequency of the system; sampling frequency F of rotating speed channel signal s ′＝D 1 ·F s ；D 1 Is an integer.
In a preferred embodiment of any of the foregoing solutions, the method for implementing counting by using a high-frequency counter to perform timing compensation on a rotational speed signal in the rotational speed channel includes:
the acquisition board counts the rotating speed signals at high frequency, the counting start time is the rotating speed pulse triggering time, the counting stop timing time is the rotating speed sampling frequency reaching time, and the steps are repeated in turn.
In a preferred embodiment of any of the foregoing aspects, the timing compensation for the rotation speed signal using a high frequency counter includes the steps of:
taking the abscissa as each sampling time point t k 1 The ordinate is the count value m of the high-frequency counter corresponding to each sampling time point k 1 Coordinate point
The corresponding minimum count value and the corresponding time point, which represents the sampling interval, the rotation speed signal passes through the real rotation speed triggering position, the counter data is reset to zero and then the time is re-counted, and the counter is stopped>
For sampling time +.>
The corresponding count value is used for obtaining the actual rotation speed trigger moment +.>
Wherein F is h For the main frequency of the board high-frequency counter, the error range of the compensated rotating speed trigger time is from 1/F s ' reduced to 1/F h For the rotation speed monitoring shaft, every time N pulse numbers pass, the rotation angle of the rotation speed monitoring shaft is 2 pi, and according to the transmission ratio relation, the rotation angle of the tracked shaft is R.2pi.
In a preferred embodiment of any of the foregoing aspects, the method further comprises:
according to the obtained time sequence after the compensation of each rotating speed;
using the trigger phase sequence R.cndot.02 pi.4 pi.6 pi]And corresponding time sequence [ t ] 0 t 1 t 2 t 3 ]Performing primary spline interpolation;
obtain equal angular interval
Corresponding interpolation time sequence t 0 ′ t 1 ′ t 2 ′ … t n′ ′]Wherein the equal angle value ∈ ->
Where order represents the number of resampling points per turn of the setting.
In a preferred embodiment of any of the foregoing aspects, the method further comprises:
according to the vibration sampling frequency F s Acquiring a time sequence [0 delta t 2 delta t … n delta t ]]And vibration sequence [ vib ] 0 vib 1 vib 2 …vib n ]Wherein the time interval deltat in the time sequence is the sampling frequency F s Is the reciprocal of (2);
for time sequence [ 0.DELTA.t2DELTA.t … nDELTA.t]And vibration sequence [ vib ] 0 vib 1 vib 2 … vib n ]Performing secondary cubic spline interpolation calculation;
calculating a time sequence [ t ] obtained by first interpolation 0 ′ t 1 ′ t 2 ′ … t n′ ′]Corresponding vibration signal [ vib ] 0 ′ vib 1 ′ vib 2 ′ … vib n′ ′]Obtaining the resampled angle domainVibration sequences, the angle intervals between the vibration sequences are all
In a preferred embodiment of any of the foregoing aspects, the method further comprises:
taking the first derivative of the calculated front and rear data segment splice points as the prior condition of cubic spline interpolation comprises the following steps:
in the data section
In the method, the same data point is continuous in first order conduction between the front section and the rear section, and the first section interpolation sequence selects data section +.>
First data point [ x ] at default acquisition 0 ,y 0 ]The interpolation is not participated, and the first order derivative of the first point of the first section interpolation sequence is calculated;
the first-end two-point first-order derivative calculation formula of the first-segment interpolation sequence is as follows:
for the second segment of interpolation sequence
The first-order derivative calculation formula of the first and the last points is as follows:
data point x 4 ,y 4 ]As the splicing point of the front and rear data, the first derivative value of the two interpolation operations at the point is equal.
Compared with the prior art, the vibration order tracking analysis method adjusts the acquisition strategy of vibration and rotating speed signals according to actual working condition parameters, avoids the problem of resource waste caused by synchronous increase of sampling frequency on the premise of guaranteeing acquisition of effective signals, acquires the rotating speed signals through the matched plate card high-frequency counter, records the value of the counter at each sampling point, and acquires key phase information with higher precision through a compensation algorithm.
According to the method, the prior condition of the traditional cubic spline interpolation is optimized, interpolation operation is carried out by selecting the first derivative of the first point and the second point as the prior condition, the first derivative value of the splicing point of each data segment is ensured to be equal, the splicing effect after the segmented interpolation of the data segment is consistent with the direct interpolation result of the whole segment, and the practical requirements of segmented acquisition, analysis and storage of the data are met.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application 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 application. Some specific embodiments of the present application will be described in detail hereinafter by way of example and not by way of limitation with reference to the accompanying drawings. The same reference numbers in the drawings identify identical or similar components or portions thereof, and it should be understood by those skilled in the art that the drawings are not necessarily drawn to scale and that:
fig. 1 is a flow chart of a vibration order tracking analysis method according to an embodiment of the present application.
Fig. 2 is a schematic diagram of phase compensation of a vibration order tracking analysis method according to an embodiment of the present application.
Fig. 3 is a schematic diagram of cubic spline prior condition acquisition of the vibration order tracking analysis method according to the embodiment of the present application.
Detailed Description
In order to enable those skilled in the art to better understand the present application, the following description will make clear and complete descriptions of the technical solutions in the embodiments of the present application with reference to the accompanying drawings in the embodiments of the present application. It is to be understood that the described embodiments are merely partial embodiments of the present application and not full embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments herein without making any inventive effort, shall fall within the scope of the present application.
The following embodiments of the present application will describe the embodiments of the present application in detail by taking the example of the vibration order tracking analysis method having a front wheel and a rear wheel, but the embodiments should not limit the scope of the present application.
The invention provides a high-precision vibration order tracking analysis method under an asynchronous rotating speed, which aims to solve the problem of resource waste caused by excessive vibration sampling frequency, optimize vibration and rotating speed acquisition modes under an asynchronous rotating speed working condition, compensate key phase time by matching with an algorithm, improve phase precision, optimize priori conditions of cubic spline interpolation, realize smooth splicing of data segments before and after angle domain resampling, and improve data quality.
According to the invention, different sampling frequencies are adopted to collect vibration and rotating speed signals according to actual analysis requirements, a traditional mode of recording high-level signals and low-level signals is changed, a high-frequency counter is used for storing counting values of rotating speed pulse signals, a time compensation algorithm is used for obtaining more accurate key phase time, further improving phase precision, meanwhile, a traditional angle domain resampling priori condition is optimized, a first derivative value of a splicing point of a data segment before transmission and a first derivative value of a data segment after transmission is used as a priori condition of cubic spline interpolation, smooth splicing of vibration data of the angle domain of the data segment before transmission and the data segment after transmission is realized, the problem of discontinuous front data segment after data segment caused by original data interception is solved, and an order spectrum with higher signal to noise ratio is further obtained for subsequent analysis.
Principle elucidation: the angular domain resampling involves two interpolation operations, the most widely used at present is cubic spline interpolation, and the cubic spline interpolation specific algorithm is not the main content of this patent and is not specifically described herein. The angular domain resampling algorithm is input as vibration and rotating speed signals, and the invention provides a method for acquiring vibration and rotating speed signals by utilizing different sampling frequencies to meet the requirement of angular domain resampling data analysis under the high rotating speed working condition, and acquiring the rotating speed signals by utilizing higher sampling frequencies to acquire each pulse moment. The invention provides a method for counting rotating speed pulses by adopting a high-frequency counter arranged in an acquisition board, when the pulse signal zero crossing point of a rotating speed sensor is used for setting the counter value to zero, re-counting, recording the current counter value at each data sampling point, calculating each key phase moment by matching with an algorithm, and when the rotating speed signal is acquired by utilizing the acquisition strategy, the higher the counting frequency of the counter is, the more accurate the obtained phase information is, and the common high-frequency counting frequency of the current acquisition board is 100 megameters.
In addition, the invention provides a priori condition of cubic spline interpolation by using the first order derivative of the splicing points of the front data segment and the rear data segment, so as to realize smooth splicing after resampling processing of the angle domain of the front data segment and the rear data segment.
Examples
The embodiment of the application provides a vibration order tracking analysis method, data acquisition and analysis are realized through a board, input is vibration and rotation speed signals, output is vibration signals in an equal angle domain, different sampling frequencies are adopted for the vibration and rotation speed signals to acquire, on one hand, the relatively higher sampling frequency is adopted for acquiring the rotation speed signals to ensure effective acquisition of rotation speed pulses, on the other hand, the lower sampling frequency is adopted for acquiring the vibration signals on the premise of meeting analysis requirements, waste of acquisition, analysis and storage resources is avoided, a mode of recording the rotation speed pulses by adopting a high-frequency counter built in the board is adopted, phase compensation is carried out by collocation algorithm, and phase precision is improved. And taking the first order derivative of the splicing points of the front data segment and the rear data segment as the priori condition of cubic spline interpolation, and improving the problem of unconnected splicing of the front data segment and the rear data segment before and after the angle domain resampling.
In order to make the content of the present invention more clear, the following will further describe the present invention with reference to the accompanying drawings, the strategy of the acquisition and data analysis algorithm is shown in fig. 1, the schematic diagram of phase compensation by using high-frequency count values is shown in fig. 2, and the schematic diagram of obtaining the cubic spline interpolation first order derivative priori condition is shown in fig. 3.
As shown in fig. 1, an embodiment of the present application provides a vibration order tracking analysis method, including the following steps:
step 1: according to the tracked axis rotation frequency F r The pulse number N of each revolution of the revolution speed sensor, the transmission ratio R between the tracked shaft and the revolution speed monitoring shaft and the most monitoring systemHigh analysis frequency F m And confirming the acquisition strategy by the system parameters.
The specific mode is as follows: let the actual sampling frequency F s ＝D·F m According to the sampling theorem, the coefficient D must be greater than 2, the actual industrial application coefficient D generally takes a value of 2.5-10, and D is an integer and represents the proportionality coefficient between the sampling frequency and the highest analysis frequency. When the number of pulses per second is m=r·f r N satisfies D 0 M<F s Time (D) 0 For a reserved safety factor, preferably taking a value greater than 1.5), the current sampling frequency F s The vibration and rotating speed sampling requirements can be met simultaneously, and at the moment, the vibration and rotating speed signals adopt the same sampling frequency. When D is 0 M≥F s When the sampling frequency is low, the current sampling frequency can not collect effective rotating speed pulse, so that phase information is lost, and a sampling strategy needs to be adjusted: vibration sampling frequency holding F s The sampling frequency of the rotating speed channel signal is unchanged and is adjusted to F s ′＝D 1 ·F s Preferred coefficient D 1 Get to satisfy D 0 M<F s ' minimum integer, and coefficient D according to the FPGA downsampling rate hierarchical relationship 1 Further limiting the value of (c). Here coefficient D 1 The rounding is to ensure that the vibration sampling moments have one-to-one corresponding rotating speed pulse count values, so that the interception of data segments in the two interpolation processes is facilitated, and the calculation error introduced by interpolation is reduced.
Step 2: in order to acquire phase information with higher precision, a high-frequency counter is used for counting the rotating speed signal, and the counting rule is as follows: the acquisition board card carries out high-frequency judgment on the rotating speed signal, the numerical value of a high-frequency counter is recorded at each sampling moment, when the zero crossing point of the rotating speed signal is judged, the numerical value of the counter is set to zero and counted again, fig. 2 is a stored rotating speed signal time domain diagram, the abscissa is the time point of each sampling, the ordinate is the count value of the high-frequency counter corresponding to each sampling time point, the value of the counter is increased continuously from small to small every time when one pulse is experienced, and when the next pulse arrives, the value of the counter is set to zero and enters the next counting period.
According to the counting principle, the accurate key phase time is found to be at the maximum and minimum of the count valueBetween time periods corresponding to the values, taking fig. 2 as an example, coordinate points
Corresponding to the minimum value of the count value and the corresponding moment, wherein the moment point represents that the rotating speed signal passes through the real key phase position in the sampling interval, and the counter is re-timed after the data of the counter is set to zero, and the counter is used for counting the time of the counter>
For sampling time +.>
The corresponding count value can be used to obtain the actual key phase time +.>
Wherein F is h For the counting frequency of the board card high-frequency counter, the compensated key phase time error range is from 1/F s ' reduced to 1/F h The phase accuracy for angular domain resampling is greatly improved. For the rotation speed monitoring shaft, the rotation angle of the rotation speed monitoring shaft is 2 pi after every N pulse numbers, and the rotation angle of the tracked shaft is R.2pi according to the transmission ratio relation.
Step 3: obtaining a time sequence after each key phase signal compensation according to the step 2, and utilizing the key phase sequence R [02 pi 4 pi 6 pi ]]And corresponding time sequence [ t ] 0 t 1 t 2 t 3 ]Performing a first cubic spline interpolation (taking three rotations of the rotation speed monitoring shaft as an example) to obtain an equal angle interval
Corresponding interpolation time sequence t 0 ′ t 1 ′ t 2 ′ … t n′ ′]Δφ, which represents the angular interval, determined by the final analysis of the highest order and the rotation rate scaling factor; n 'represents the sequence length, n' is related to the number of rotational speed sequences calculated in step 2, there are 5 n, d=5, here the equiangular value +.>
According to the formula
And calculating, wherein order represents the number of resampling points of each circle, and generally taking a value according to 2.5-10 times of the highest analysis order. In order to reduce calculation errors introduced by phase interpolation, vibration data segments corresponding to integer rotation periods are selected by key phase signals for analysis (corresponding to three circles in the example) when the segments are calculated.
Step 4: according to the vibration sampling frequency F s Acquiring a time sequence [0 delta t 2 delta t … n delta t ]]And vibration sequence [ vib ] 0 vib 1 vib 2 … vib n ]Wherein the time interval deltat in the time sequence is the sampling frequency F s Is the inverse of (c). For time sequence [ 0.DELTA.t2DELTA.t … nDELTA.t]And vibration sequence [ vib ] 0 vib 1 vib 2 … vib n ]Performing secondary cubic spline interpolation calculation, and calculating a time sequence [ t ] obtained by the primary interpolation in the step (3) 0 ′ t 1 ′ t 2 ′ … t n′ ′]Corresponding vibration signal [ vib ] 0 ′ vib 1 ′ vib 2 ′ … vib n′ ′]Obtaining vibration sequences after angle domain resampling, wherein the angle intervals between the vibration sequences are all
Step 5: step 3 and step 4 both involve cubic spline interpolation, and for angle domain resampling of a single data segment, most of the conventional algorithms currently take zero first order of the first and last data points as prior conditions of cubic spline interpolation. However, for continuously collected data, because the buffer space is limited, the analyzed data segments must be collected, analyzed and stored in a segmented manner, so that new requirements are put forward for splicing the data segments after resampling of different data segments, the traditional prior condition that the first order leads of the head and tail points of the data segments are zero can cause the first order leads of the data splicing points to be discontinuous, and the data segments are spliced after interpolation to have mutation points, so that the data quality is poor.
In order to solve the problem, the invention provides a corresponding solution, and the following is a schematic diagram of fig. 3, which illustrates the acquisition mode of prior conditions when performing cubic spline interpolation operation on the front and rear data segments:
to ensure that the two data sets can be smoothly spliced after cubic spline interpolation, the first order conduction continuity of the same data point between the front section and the rear section is ensured, and the data section of FIG. 3
For example, the first interpolation sequence selects data segment +.>
First data point [ x ] at default acquisition 0 ,y 0 ]The interpolation is not participated, the first order guide of the first point of the first section interpolation sequence is calculated, and the first order guide calculation formula of the first point and the last point of the first section interpolation sequence is as follows:
whereas for the second segment interpolation sequence +.>
The first-order derivative calculation formula of the first and the last points is as follows: />
Data points [ x ] can be seen 4 ,y 4 ]As the splicing point of the front section data and the rear section data, the first derivative value of the two interpolation operations at the point is equal, and the smooth transition of the curve after the two interpolation operations at the point is ensured. The calculation strategy is suitable for the step 3 and the step 4, and the data segments obtained by the front interpolation operation and the rear interpolation operation can be smoothly spliced by ensuring that the first order leads of the first point and the tail point of the data segments adopt the same calculation mode in each interpolation link.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or the technical features of the components or the whole components can be replaced equivalently; such modifications and substitutions do not depart from the spirit of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present application.
Claims (3)
1. A vibration order tracking analysis method, the method comprising the steps of:
step 1: under the pseudo-synchronous sampling condition, real-time vibration and rotation speed signal acquisition are carried out;
step 2: adopting a high-frequency counter to perform timing compensation on the rotating speed pulse to acquire an instantaneous rotating speed and an instantaneous phase sequence synchronous with rotating speed sampling;
step 3: obtaining a sampling time list sequence corresponding to the phase to be interpolated by using an interpolation resampling method according to the rotating speed sampling frequency and the compensated instantaneous phase sequence;
step 4: according to the interpolated phase sequence and the corresponding sampling time list, combining the vibration sampling frequency and the time list generated by corresponding vibration, and obtaining a new interpolated vibration sequence at the corresponding moment of the interpolated phase sequence by using an interpolation method again;
step 5: repeating the steps 1 to 5 to complete continuous order tracking of the signals;
the sampling rates of the real-time vibration and the rotation speed sampling can be inconsistent, and the following relation needs to be satisfied:
when the device pulses m=r·f per second r N, wherein F r The rotation frequency of the tracking shaft is set, N is the pulse number of each rotation of the rotation speed sensor, and R is the transmission ratio between the tracked shaft and the rotation speed monitoring shaft; vibration sampling frequency F s ＝D·F m ，D≥2，F m Monitoring the highest analysis frequency of the system; sampling frequency F of rotating speed channel signal s ′＝D 1 ·F s ，D 1 Is an integer;
the rotating speed channel adopts a high-frequency counter to carry out timing compensation on rotating speed signals, and the counting implementation method comprises the following steps:
the acquisition board counts the rotating speed signal at high frequency, the counting start time is the rotating speed pulse triggering time, the counting stop timing time is the rotating speed sampling frequency reaching time, and the steps are repeated in turn;
the method adopts a high-frequency counter to carry out timing compensation on the rotating speed signal, and the compensation comprises the following steps:
taking the abscissa as each sampling time point t k 1 The ordinate is the count value m of the high-frequency counter corresponding to each sampling time point k 1 Coordinate point
The corresponding minimum count value and the corresponding time point, which represents the sampling interval, the rotation speed signal passes through the real rotation speed triggering position, the counter data is reset to zero and then the time is re-counted, and the counter is stopped>
For sampling time +.>
The corresponding count value is used for obtaining the actual rotation speed trigger moment +.>
Wherein F is h For the main frequency of the board high-frequency counter, the error range of the compensated rotating speed trigger time is from 1/F s ' reduced to 1/F h For the rotation speed monitoring shaft, the rotation angle of the rotation speed monitoring shaft is 2 pi after N pulse numbers pass, and according to the transmission ratio relation, the rotation angle of the tracked shaft is R.2pi;
the method further comprises the steps of:
according to the obtained time sequence after the compensation of each rotating speed;
using the trigger phase sequence R.cndot.02 pi.4 pi.6 pi]And corresponding time sequence [ t ] 0 t 1 t 2 t 3 ]Performing primary spline interpolation;
obtain equal angular interval
Corresponding interpolation time sequence t 0 ′ t 1 ′ t 2 ′ … t n′ ′]Wherein the equal angle value ∈ ->
Where order represents the number of resampling points per turn of the setting.
2. The vibration order tracking analysis method according to claim 1, characterized in that the method further comprises:
according to the vibration sampling frequency F s Acquiring a time sequence [0 delta t 2 delta t … n delta t ]]And vibration sequence [ vib ] 0 vib 1 vib 2 … vib n ]Wherein the time interval deltat in the time sequence is the sampling frequency F s Is the reciprocal of (2);
for time sequence [ 0.DELTA.t2DELTA.t … nDELTA.t]And vibration sequence [ vib ] 0 vib 1 vib 2 … vib n ]Performing secondary cubic spline interpolation calculation;
calculating a time sequence [ t ] obtained by first interpolation 0 ′ t 1 ′ t 2 ′ … t′ n′ ]Corresponding vibration signal [ vib ] 0 ′ vib 1 ′ vib 2 ′ … vib′ n′ ]Obtaining vibration sequences after angle domain resampling, wherein the angle intervals between the vibration sequences are all
3. The vibration order tracking analysis method according to claim 2, characterized in that the method further comprises:
taking the first derivative of the calculated front and rear data segment splice points as the prior condition of cubic spline interpolation comprises the following steps:
in the data section
In the method, the same data point is continuous in first order conduction between the front section and the rear section, and the first section interpolation sequence selects data section +.>
First data point [ x ] at default acquisition 0 ,y 0 ]The interpolation is not participated, and the first order derivative of the first point of the first section interpolation sequence is calculated;
the first-end two-point first-order derivative calculation formula of the first-segment interpolation sequence is as follows:
for the second segment of interpolation sequence
The first-order derivative calculation formula of the first and the last points is as follows:
data point x 4 ,y 4 ]As the splicing point of the front and rear data, the first derivative value of the two interpolation operations at the point is equal. />
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111587611.7A CN114485914B (en) | 2021-12-23 | 2021-12-23 | Vibration order tracking analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111587611.7A CN114485914B (en) | 2021-12-23 | 2021-12-23 | Vibration order tracking analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114485914A CN114485914A (en) | 2022-05-13 |
| CN114485914B true CN114485914B (en) | 2023-05-30 |
Family
ID=81493882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111587611.7A Active CN114485914B (en) | 2021-12-23 | 2021-12-23 | Vibration order tracking analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114485914B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6351714B1 (en) * | 1998-03-03 | 2002-02-26 | Entek Ird International Corporation | Order tracking signal sampling process |
| CN103196547A (en) * | 2013-03-11 | 2013-07-10 | 安徽新力电业科技咨询有限责任公司 | Method achieving rotary machinery vibration signal synchronization order ratio tracing analysis |
| CN203132688U (en) * | 2013-03-11 | 2013-08-14 | 安徽新力电业科技咨询有限责任公司 | Device for realizing synchronous order-domain trace analysis on rotating mechanical vibration signals |
| CN110887663A (en) * | 2019-10-30 | 2020-03-17 | 中国石油化工股份有限公司 | Bearing fault diagnosis method combining variable working condition calculation order tracking and spectral kurtosis |
| CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
| CN112665712A (en) * | 2020-12-22 | 2021-04-16 | 浙江省交通投资集团有限公司智慧交通研究分公司 | Wide-area order tracking method and system for monitoring train running gear |
-
2021
- 2021-12-23 CN CN202111587611.7A patent/CN114485914B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6351714B1 (en) * | 1998-03-03 | 2002-02-26 | Entek Ird International Corporation | Order tracking signal sampling process |
| CN103196547A (en) * | 2013-03-11 | 2013-07-10 | 安徽新力电业科技咨询有限责任公司 | Method achieving rotary machinery vibration signal synchronization order ratio tracing analysis |
| CN203132688U (en) * | 2013-03-11 | 2013-08-14 | 安徽新力电业科技咨询有限责任公司 | Device for realizing synchronous order-domain trace analysis on rotating mechanical vibration signals |
| CN110887663A (en) * | 2019-10-30 | 2020-03-17 | 中国石油化工股份有限公司 | Bearing fault diagnosis method combining variable working condition calculation order tracking and spectral kurtosis |
| CN111307460A (en) * | 2020-03-14 | 2020-06-19 | 中国石化销售股份有限公司华南分公司 | Rolling bearing fault diagnosis method based on calculation order tracking and spectral kurtosis |
| CN112665712A (en) * | 2020-12-22 | 2021-04-16 | 浙江省交通投资集团有限公司智慧交通研究分公司 | Wide-area order tracking method and system for monitoring train running gear |
Non-Patent Citations (1)
| Title |
|---|
| 基于自适应计算阶次跟踪的旋转机械阶次分析系统;魏玉果;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第5期);第10-35页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114485914A (en) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102270961B (en) | Signal subdividing method and subdividing device for sine-cosine encoder, and motor control device | |
| CN103196547B (en) | One realizes the synchronous Computed order tracking analytical approach of characteristic of rotating machines vibration signal | |
| CN109945819B (en) | Method for measuring position of rotor of permanent magnet synchronous motor | |
| CN107560724B (en) | Vibration signal analysis method | |
| CN104006962A (en) | Gear fault feature extraction method and system | |
| CN108181482A (en) | Real-time low-speed detection device based on virtual sine wave | |
| CN205991836U (en) | A kind of characteristic of rotating machines vibration signal acquisition analysis system | |
| CN109682597A (en) | A kind of gear-box vibration signal processing and analysis method | |
| CN102967270A (en) | Method and system for measuring engine tip clearance | |
| CN106645780B (en) | DSP-based rotating speed detection method and system | |
| CN110750484B (en) | Synchronous acquisition system and acquisition method for data of rotating speed and multiple vibration channels | |
| CN114509159B (en) | Order tracking analysis method, system and computer readable storage medium | |
| CN105245203A (en) | System and method for duty ratio detection employing high precision and low speed clock | |
| CN111953252A (en) | Synchronous sampling method for three-stage brushless synchronous motor variable-frequency alternating-current power generation system | |
| CN114485914B (en) | Vibration order tracking analysis method | |
| CN109188077A (en) | A kind of signal frequency measuring method, system, device and readable storage medium storing program for executing | |
| CN203132688U (en) | Device for realizing synchronous order-domain trace analysis on rotating mechanical vibration signals | |
| CN109959407B (en) | Method and device for quickly detecting linearity of steering engine position sensor | |
| CN104677486B (en) | The aero-engine vibration signal Method for Phase Difference Measurement reconstructed based on tacho-pulse | |
| CN109558041A (en) | Tip clearance signal acquisition, processing and the transmission method accelerated based on GPU | |
| CN114858429B (en) | Equal-angle sampling method and order tracking analysis method for vibration signals of rotary machine | |
| CN1502167A (en) | Apparatus and method for adjusting filter frequency in relation to sampling frequency | |
| CN103278332B (en) | The engine startup detection method to cylinder misfire | |
| CN102832911B (en) | A kind of digital signal recovery method and device | |
| CN115270853A (en) | Bearing fault diagnosis self-adaptive input method and system based on deep learning |
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 |