LabVIEW
TM
Order Analysis ToolsetUser Manual
LabVIEW Order Analysis Toolset User Manual
August 2003 Edition
Part Number 322879B-01
Страница 1: ...LabVIEW TM Order Analysis Toolset User Manual LabVIEW Order Analysis Toolset User Manual August 2003 Edition Part Number 322879B 01...
Страница 2: ...rmany 49 0 89 741 31 30 Greece 30 2 10 42 96 427 India 91 80 51190000 Israel 972 0 3 6393737 Italy 39 02 413091 Japan 81 3 5472 2970 Korea 82 02 3451 3400 Malaysia 603 9131 0918 Mexico 001 800 010 079...
Страница 3: ...ccident actions of third parties or other events outside reasonable control Copyright Under the copyright laws this publication may not be reproduced or transmitted in any form electronic or mechanica...
Страница 4: ...sis of Rotating Machinery 1 2 System Requirements 1 3 Installation 1 3 Example VIs 1 4 Acquiring Data for Example VIs 1 4 Configuring DAQ Hardware Used with Examples 1 4 Acquire Data Analog Tach VI 1...
Страница 5: ...8 Calculating Waveform Magnitude 3 10 Chapter 4 Resampling Based Order Analysis LabVIEW Order Analysis Toolset Resampling Method 4 1 Determining the Time Instance for Resampling 4 2 Resampling Vibrat...
Страница 6: ...learn about how to use the LabVIEW Order Analysis Toolset Refer to Chapter 1 Introduction to the LabVIEW Order Analysis Toolset for information about the example VIs If you want to learn more about t...
Страница 7: ...ters that you should enter from the keyboard sections of code programming examples and syntax examples This font is also used for the proper names of disk drives paths directories programs subprograms...
Страница 8: ...ther time or rotational speed Extraction of the order components from the original noise or vibration signal Measurement of the magnitude and phase of any order component as a function of rotational s...
Страница 9: ...s for the Analysis of Rotating Machinery Order analysis is a powerful tool for analyzing rotating machinery when the rotational speed might change over time However to successfully use the LabVIEW Ord...
Страница 10: ...isable any automatic virus detection programs before you install After installation check your hard disk for viruses and enable any virus detection programs you disabled Windows 2000 NT XP Complete th...
Страница 11: ...have a Boolean control named Data Source The Data Source control has two choices Example and DAQ When you choose Example the VI uses prerecorded data generated during a fan run up as the data source W...
Страница 12: ...out the Acquire Data with PXI 4472 and TIO VI Acquire Data Analog Tach VI In analog tachometer examples setting Data Source to DAQ and clicking the Run button opens the Acquire Data Analog Tach VI The...
Страница 13: ...Figure 1 2 Figure 1 2 Acquire Data Analog Tach VI Acquisition Tab The Acquisition tab shown in Figure 1 2 allows you to acquire and observe data Click the Acquire button to acquire data Continue to co...
Страница 14: ...e the DAQ devices Use one of the counters on a TIO device to receive TTL compatible tachometer pulses Use the controls in the TIO Board Setting Digital Tach Signal section of the Configuration tab to...
Страница 15: ...th PXI 4472 and TIO VI Acquisition Tab The Acquisition tab shown in Figure 1 4 allows you to acquire and observe data Click the Acquire button to acquire data Continue to configure the data acquisitio...
Страница 16: ...diagnosing heart problems You can use order analysis to study design and monitor rotating machinery By measuring and analyzing sound or vibration signals generated by a system with rotational compone...
Страница 17: ...ysis for rotating machinery Figure 2 1 shows the relationship between frequency and order spectra Figure 2 1 Order and Frequency Domain Display of a Shaft Rotating at 3 000 rpm The top graph in Figure...
Страница 18: ...lts in a spectral peak at the first order Misalignment or bending of the shaft generates a large second order Oil whirl might lead to strong fractional orders Gears belts and blades might enhance high...
Страница 19: ...ysis is effective only when the rotational speed remains constant Effect of Rotational Speed on Order Identification The ability to make a reliable identification of individual orders from the convent...
Страница 20: ...d on the fast Fourier transform FFT The upper right plot in Figure 2 3 shows the frequency time spectral map computed from the short time Fourier transform STFT with a 1 024 point Hanning window Becau...
Страница 21: ...change in the frequency and amplitude of the signal The faster the frequency changes the wider the overall frequency bandwidth becomes as measured from the power spectrum Figure 2 4 illustrates a con...
Страница 22: ...e overall frequency bandwidth of the signal whose frequency increases with time in Figure 2 5 is much wider than that of the signal whose frequency is constant in Figure 2 4 When frequency or amplitud...
Страница 23: ...onal speed increases both the fundamental frequency bandwidth and the frequency bandwidths of related orders widen causing orders to overlap Whether you can separate nearby orders in a power spectrum...
Страница 24: ...annot control the rotational speed of the system well or the system is not allowed to run at the critical rotational speed for a sufficient length of time Order Analysis An important goal of order ana...
Страница 25: ...bor expansion you cannot reconstruct the general STFT using an inverse Fourier transform You can compute the Gabor transform by either STFT or windowed Fourier transform However to ensure reconstructi...
Страница 26: ...r expansion the reconstructed data is longer than the original data To avoid the reconstructed data being longer than the original data consider the signal as periodic and use wrap padding shown in Fi...
Страница 27: ...out the Gabor transform and the Gabor expansion Refer to Chapter 3 Gabor Transform Based Order Tracking for information about how the LabVIEW Order Analysis Toolset uses a method based on the Gabor tr...
Страница 28: ...r or a multiplied tachometer signal to trigger the analog to digital conversion which ensures a sampling process spaced equally in the rotation angle The hardware solution requires additional hardware...
Страница 29: ...he spectrum of angle samples Figure 2 10 Spectrum of Angle Samples In Figure 2 10 the Angle Samples plot shows the signal sampled at a constant angle interval The Order Spectrum plot shows the order s...
Страница 30: ...onents change as the rotational speed changes you can consider the rotational speed and the frequency of order components to remain constant in a relatively short time interval The adaptive filter met...
Страница 31: ...ise or vibration signal Reconstruct the desired order components in the time domain Overview of Gabor Order Analysis The Gabor transform can give the power distribution of the original signal as the f...
Страница 32: ...hite curves have magnitudes larger than the magnitudes in local neighborhoods off the curves The white curves indicate the order components and are referred to as order curves From the frequency time...
Страница 33: ...performing Gabor expansion f Calculate the waveform magnitude and phase Refer to the Important Considerations for the Analysis of Rotating Machinery section of Chapter 1 Introduction to the LabVIEW O...
Страница 34: ...some coefficients in the neighborhood You must include the coefficients in the neighborhood when selecting the order component Constant frequency bandwidth and constant order bandwidth are two ways t...
Страница 35: ...er curves select a subset of the initial Gabor coefficient array The subset of the initial Gabor coefficient array contains only the coefficients in the neighborhood of the selected order curves You c...
Страница 36: ...number of elements set to TRUE determines the bandwidth of the passband Extracting Orders Figure 3 4 illustrates the order extraction and signal reconstruction process Figure 3 4 Gabor Coefficients Fi...
Страница 37: ...3 5 Original Signal and Extracted Order Component Unlike the reverse discrete Fourier transform the Gabor expansion in general is not a one to one mapping A Gabor coefficient is the subspace of a two...
Страница 38: ...outside the passband the magnitude of the reconstructed signal is no longer zero Instead the magnitude does have some certain value However the value of the magnitude quickly decreases as the frequenc...
Страница 39: ...ar to the rpm time function shown in Figure 3 7 b In Figure 3 7 a notice the horizontal white lines The horizontal white lines indicate the large power around the resonance frequencies The physical ch...
Страница 40: ...ral map and the order rpm spectral map you can clearly and efficiently observe all the order components in the entire time and rpm ranges Calculating Waveform Magnitude The reconstructed time waveform...
Страница 41: ...to complete the following tasks Resample even time space samples to even angle spaced samples Obtain magnitude and phase information for each order Figure 4 1 illustrates the resampling based order a...
Страница 42: ...r magnitudes and phases if necessary Refer to the Important Considerations for the Analysis of Rotating Machinery section of Chapter 1 Introduction to the LabVIEW Order Analysis Toolset for informatio...
Страница 43: ...me sequence for a smaller angle The transfer function of the CIC filter is given by the following equation where L is the interpolation factor and n is the order The CIC filter has the advantage of on...
Страница 44: ...le spaced samples you must be able to calculate the value of the vibration signal at any time instance According to the Nyquist sampling theorem you can exactly reconstruct the signal for all time ins...
Страница 45: ...s lower that the previous sampling rate The LabVIEW Order Analysis Toolset uses a digital adaptive interpolation filter to complete the entire resample process The bandwidth of the adaptive interpolat...
Страница 46: ...nsated signal Figure 4 4 Slow Roll Compensation The LabVIEW Order Analysis Toolset can perform the vector based slow roll compensation on order magnitudes and phases The toolset performs slow roll com...
Страница 47: ...al speed as a function of time during the data acquisition process The LabVIEW Order Analysis Toolset uses a digital differentiator to calculate the rotational speed The tachometer processing VIs buff...
Страница 48: ...wing equation you can present rotational speed as the first derivative of k While the time interval between pulses varies with the rotational speed the cumulative angle increment remains constant Calc...
Страница 49: ...ate the rotational speed also increases The LabVIEW Order Analysis Toolset can automatically adjust the order of the differentiator according to the number of pulses in the buffer Automatically adjust...
Страница 50: ...In Figure 5 2 the dotted line represents the results of calculating the rotational speed with no averaging The solid line in Figure 5 2 represents the results of calculating the rotational speed avera...
Страница 51: ...mplementation was an open research topic until Bastiaans discovered the relationship between the Gabor expansion and the short time Fourier transform STFT in the early 1980s Over the years many differ...
Страница 52: ...ansform in Equation A 2 contains redundancy from a mathematical point of view However the redundancy in Equation A 2 provides freedom for the selection of better window functions h k and k Notice that...
Страница 53: ...or LMSE For the LMSE solution the Euclidean distance between the dual functions is minimum as given by the following equation where denotes the matrix form of Equation A 4 When the error is small and...
Страница 54: ...t exist For example the following equation represents a modified two dimensional function where m n denotes a binary mask function whose elements are either 0 or 1 Applying the Gabor expansion to the...
Страница 55: ...the Gabor expansion from Equation A 1 and the Gabor transform from Equation A 2 in matrix form as shown in the following equations where H denotes the analysis matrix and G denotes the synthesis matri...
Страница 56: ...est in terms of the LMSE to the masked Gabor coefficients C The masked Gabor coefficients are the desirable Gabor coefficients The case of k h k usually implies considerable over sampling which result...
Страница 57: ...er Order Tracking Methods SAE Noise and Vibration Conference and Exposition Traverse City MI April 30 through May 3 2001 Gade S H Herlufsen H Konstantin Hansen and H Vold Characteristics of the Vold K...
Страница 58: ...ted Support Options Contact NI engineers and other measurement and automation professionals by visiting ni com support Our online system helps you define your question and connects you to the experts...
Страница 59: ...ency less than half the sampling frequency See also anti aliasing filter and Nyquist frequency anti aliasing filter Analog lowpass filters used before analog to digital conversion to filter out the fr...
Страница 60: ...int time frequency transforms H harmonic Frequencies that are integer or fractional multiples of a fundamental frequency L LMSE Least mean square error N Nyquist frequency Half the sampling frequency...
Страница 61: ...k over sampling Occurs in a Gabor transform when the window length is greater than the window shift step In over sampling the number of Gabor coefficients cm n is more than the number of original dat...
Страница 62: ...1 4 guidelines for 1 4 diagnostic resources C 1 displaying spectral maps 3 8 documentation conventions used in the manual vii getting started vii how to use this manual vii online library C 1 related...
Страница 63: ...Up Kit 1 1 Order Analysis Toolset using in analysis of rotating machinery 1 2 orders definition of 2 2 over sampling definition of A 2 P phone technical support C 1 processing tachometer signals 5 1 c...
