Main Content

gearMeshFaultBands

Construct frequency bands around the characteristic fault frequencies of meshing gears for spectral feature extraction

所有的页面崩溃

Syntax

FB = gearMeshFaultBands(FR,Ni,No)
FB = gearMeshFaultBands(___,Name,Value)
[FB,info] = gearMeshFaultBands(___)
gearMeshFaultBands (___)

Description

example

FB= gearMeshFaultBands(FR,Ni,No)generates characteristic fault frequency bandsFBof gear mesh using the rotational speed of the input gearFRand the number of teeth on the inputNiand output gearNorespectively. The values inFBhave the same implicit units asFR

example

FB= gearMeshFaultBands(___,Name,Value)allows you to specify additional parameters using one or more name-value pair arguments.

example

[FB,info] = gearMeshFaultBands(___)also returns the structureinfocontaining information about the generated fault frequency bandsFB.

example

gearMeshFaultBands (___)with no output arguments plots a bar chart of the generated fault frequency bandsFB.

Examples

collapse all

Open Live Script

For this example, consider a simple gear set with an 8-toothed pinion on the input shaft meshing with a 42-toothed spur gear on the output shaft. Assume that the input shaft is spinning at 20 rpm. Construct the gear mesh frequency bands using the physical characteristics of the gear set.

Ni = 8; No = 42; FR = 20; [FB,info] = gearMeshFaultBands(FR,Ni,No)
FB =5×219.0000 21.0000 2.8095 4.8095 79.0000 81.0000 159.0000 161.0000 159.0000 161.0000
info =struct with fields:Centers: [20 3.8095 80 160 160] Labels: ["1Fi" "1Fo" "1Fa" "1Fm" "1Fm"] FaultGroups: [1 2 3 4 5]

FBis a 5x2 array which includes the primary frequencies1Fi,1Fo,1Faand1Fmrespectively. The structure info contains the center frequencies and labels of each frequency range inFB.

Open Live Script

For this example, consider a simple gear set with an 8-toothed pinion on the input shaft meshing with a 42-toothed spur gear on the output shaft. Assume that the input shaft is driven at 20 Hz. The datasetmotorSignal.mat包含vibration data for the gear mesh sampled at 1500 Hz.

First, construct the gear mesh frequency bands using the physical characteristics of the gear set. Construct the frequency bands with the first 3 sidebands and specify the 'Domain' as '订单'.

Ni = 8; No = 42; FR = 20; FB = gearMeshFaultBands(FR,Ni,No,'Sidebands',1:3,'Domain','order')
FB =15×20.9500 1.0500 0.1405 0.2405 3.9500 4.0500 4.9500 5.0500 5.9500 6.0500 6.9500 7.0500 8.9500 9.0500 9.9500 10.0500 10.9500 11.0500 7.3786 7.4786 ⋮

FBis a 15x2 array which includes the primary frequencies and their sidebands.

Load the vibration data and compute PSD and frequency grid usingpspectrum. Use a frequency resolution of 0.5.

load('motorSignal.mat','C'); fs = 1500; [psd,freqGrid] = pspectrum(C,fs,'FrequencyResolution',0.5);

Now, use the frequency bands and PSD data to compute the spectral metrics.

spectralMetrics = faultBandMetrics(psd,freqGrid,FB)
spectralMetrics=1×46 tablePeakAmplitude1 PeakFrequency1 BandPower1 PeakAmplitude2 PeakFrequency2 BandPower2 PeakAmplitude3 PeakFrequency3 BandPower3 PeakAmplitude4 PeakFrequency4 BandPower4 PeakAmplitude5 PeakFrequency5 BandPower5 PeakAmplitude6 PeakFrequency6 BandPower6 PeakAmplitude7 PeakFrequency7 BandPower7 PeakAmplitude8 PeakFrequency8 BandPower8 PeakAmplitude9 PeakFrequency9 BandPower9 PeakAmplitude10 PeakFrequency10 BandPower10 PeakAmplitude11 PeakFrequency11 BandPower11 PeakAmplitude12 PeakFrequency12 BandPower12 PeakAmplitude13 PeakFrequency13 BandPower13 PeakAmplitude14 PeakFrequency14 BandPower14 PeakAmplitude15 PeakFrequency15 BandPower15 TotalBandPower ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________ 0.82564 1 0.1542 0.057165 0.1875 0.011175 0.29169 4 0.055249 0.011486 5 0.0021583 0.070117 6 0.013877 1.0514 7 0.21675 0.0077621 9 0.001577 0.004752 10 0.0010282 0.012155 11 0.0025085 7.7318 7.4375 1.4057 4.2222 7.625 0.79678 0.92456 7.8125 0.1924 0.030489 8.1875 0.0060835 0.069138 8.375 0.012642 0.068649 8.5625 0.012578 2.8848

spectralMetricsis a 1x46 table with peak amplitude, peak frequency and band power calculated for each frequency range inFB. The last column inspectralMetricsis the total band power, computed across all 15 frequencies inFB.

Open Live Script

For this example, consider a simple pinion and gear set with an input shaft speed of 1800 rpm. Considering that the pinion on the input shaft has 6 teeth and the gear on the output shaft has 8 teeth, visualize the frequency bands for the gear mesh.

FR = 1800; Ni = 6; No = 8; gearMeshFaultBands(FR,Ni,No)

Figure contains an axes. The axes with title Fault Frequency Bands contains 10 objects of type patch, text.

From the plot, observe the following:

  • Output shaft defect frequency,1Foat 1350 Hz

  • Input shaft defect frequency,1Fiat 1800 Hz

  • Assembly phase defect frequency,1Faat 5400 Hz

  • Gear mesh defect frequency,1Fmat 10800 Hz

Input Arguments

collapse all

Rotational speed of the input gear, specified as a positive scalar.FRis the fundamental frequency around whichgearMeshFaultBandsgenerates the fault frequency bands. SpecifyFReither in Hertz or revolutions per minute.

Number of teeth on the input gear, specified as a positive integer.

Number of teeth on the output gear, specified as a positive integer.

Name-Value Pair Arguments

Specify optional comma-separated pairs ofName,Valuearguments.Nameis the argument name andValueis the corresponding value.Namemust appear inside quotes. You can specify several name and value pair arguments in any order asName1,Value1,...,NameN,ValueN.

Example:...,'Harmonics',[1,3,5]

Harmonics of the fundamental frequency to be included, specified as the comma-separated pair consisting of 'Harmonics' and a vector of positive integers. The default value is 1. Specify 'Harmonics' when you want to construct the frequency bands with more harmonics of the fundamental frequency.

显然在基频及其周围harmonics to be included, specified as the comma-separated pair consisting of 'Sidebands' and a vector of nonnegative integers. The default value is 0. Specify 'Sidebands' when you want to construct the frequency bands with sidebands around the fundamental frequency and its harmonics.

Width of the frequency bands centered at the nominal fault frequencies, specified as the comma-separated pair consisting of 'Width' and a positive scalar. The default value is 10 percent of the fundamental frequency. Avoid specifying 'Width' with a large value so that the fault bands do not overlap.

Units of the fault band frequencies, specified as the comma-separated pair consisting of 'Domain' and either'frequency'or'order'. Select:

  • 'frequency'if you wantFBto be returned in the same units asFR.

  • 'order'if you wantFB返回,number of rotations relative toFR.

Output Arguments

collapse all

Fault frequency bands, returned as anNx2array, whereNis the number of fault frequencies.FBis returned in the same units asFR, in either Hertz or orders. Use the generated fault frequency bands to extract spectral metrics usingfaultBandMetrics. The generated fault bands, [ F − W 2 , F + W 2 ] , are centered at the characteristic defect frequencies and their harmonics and sidebands for:

  • Input shaft defect frequency,Fi

  • Output shaft defect frequency,Fo

  • Gear mesh defect frequency,Fm

  • Assembly phase pass defect frequency,Fa

When you specify the sidebands,gearMeshFaultBandscomputes the sidebands with respect to the input and output shaft defect frequencies:

  • Fault frequency bands for input gear defects with its harmonics and the first sideband atFi

  • Fault frequency bands for output gear defects with its harmonics and the first sideband atFo

gearMeshFaultBandstruncates negative fault frequency bands automatically and generates a warning message.

The valueWis the width of the frequency bands, which you can specify using the 'Width' name-value pair.

Information about the fault frequency bands inFB, returned as a structure with the following fields:

  • Centers——中心故障频率

  • Labels— Labels describing each frequency

  • FaultGroups— Fault group numbers equal to the number of frequencies

Algorithms

gearMeshFaultBandscomputes the different characteristic fault frequencies as follows:

  • Input shaft defect frequency, F i = F R

  • Output shaft defect frequency, F o = N i N o F R

  • Gear mesh defect frequency, F m = N i F R = N o F o

  • Assembly phase pass defect frequency, F a = F m gcd ( N i , N o )

References

[1] Lang, George Fox. “S&V geometry 101.”Sound and Vibration33 (1999): 16-26.

See Also

|||

Topics

Introduced in R2019b

Predictive Maintenance Toolbox Documentation

万博1manbetx

Predictive Maintenance with MATLAB

Predictive Maintenance with MATLAB

Download ebook