Анализируйте данные о зубчатой передаче и извлеките спектральные функции Используя задачи Live Editor
В этом примере показано, как использовать Извлечение Спектральный Live Editor Функций задача анализировать данные из текущего сигнала, полученного из управления зубчатой передачей сервомотора класса хобби. Пример также показывает, как извлечь спектральные функции из данных, чтобы помочь в обнаружении отказа и идентификации.
Задачи Live Editor позволяют вам в интерактивном режиме выполнить итерации на параметрах и настройках при наблюдении их эффектов на результате расчета. Задачи автоматически генерируют код MATLAB®, который достигает отображенных результатов. Для получения дополнительной информации о задачах Live Editor обычно, смотрите, Добавляют Интерактивные Задачи к Live Script.
В частности, этот пример использует Извлечение Спектральный Live Editor Функций задача. Эта задача помогает с анализом и пониманием спектральных данных. Используя всесторонний интерфейс, можно добавить компоненты, чтобы представлять различные подшипники, сетки механизма или другие части настройки оборудования. Когда вы устанавливаете физические параметры этих компонентов, Извлечение Спектральный Live Editor Функций диапазоны частот отказа графиков задачи на характеристических частотах компонентов.
Можно наложить данные о спектре мощности по графику полосы отказа сопоставить различный peaks в данных с характеристическими частотами компонентов. Это сравнение может сделать обнаружение отказа и изоляцию отказа легче, когда можно легко проследить изменения в данных о спектре мощности назад к физическим компонентам, вызывающим их.
В дополнение к графику характеристических частот и данных о спектре мощности, задача может сгенерировать спектральные метрики данных в каждом характеристическом диапазоне частот. Выходная метрическая таблица, содержащая пиковую амплитуду, пиковую частоту и мощность полосы каждого лейкопластыри в охарактеризовании потенциальных механических отказов.
Обзор оборудования
В данном примере данные об электрическом токе были собраны от стандартного сервомотора хобби Futaba S3003, который был изменен для непрерывного вращения. Сервомоторы преобразуют высокую скорость внутреннего двигателя постоянного тока к высокому крутящему моменту в выходном сплайне. Чтобы достигнуть этого, сервомоторы состоят из двигателя постоянного тока, набора нейлона или металлических механизмов диска и схемы управления. Схема управления была демонтирована, чтобы позволить текущему сигналу двигателю постоянного тока быть непосредственно проверенным. Сигнал тахометра в выходном сплайне сервомотора был собран с помощью инфракрасного фотопрерывателя наряду с 35 мм в диаметре, колеса с шестнадцатью пазами. Эти шестнадцать пазов в колесе были равномерно распределены, и фотопрерыватель IR был помещен таким образом, что это испустило точно шестнадцать импульсов на вращение колеса со слотами. Сервомотор и фотопрерыватель считались на месте пользовательским, 3-D распечатанные, монтируется.
Двигатель постоянного тока управлялся на уровне постоянных 5 вольт, и с четырьмя парами механизмов, обеспечивающих 278:1 сокращение скорости, скорость вала в сплайне составляла приблизительно 19,5 об/мин. Второй сервопривод закорачивался и использовался в качестве загрузки для системы. Текущее потребление было вычислено с помощью закона Ома путем измерения падения напряжения через резистор на 0,5 Ома. Поскольку изменение в текущих значениях измерения слишком было мало для обнаружения, текущий сигнал был усилен с помощью усилителя интерфейса датчика одно предоставления AD22050. Усиленный текущий сигнал был затем отфильтрован с помощью пятого порядка сглаживания MAX7408 эллиптический фильтр lowpass, чтобы сглаживать его и устранить шум прежде, чем отправить его в Uno Arduino через аналого-цифровой конвертер (ADC).
Когда блок-схема показывает, текущий сигнал был сначала усилен и отфильтровал использование усилителя и сглаживающегося фильтра lowpass, соответственно. Uno Arduino произвел текущий сигнал через ADC на уровне 1,5 кГц и передал его потоком к компьютеру наряду с импульсами тахометра как последовательные данные на скорости в бодах 115 200 бит/с. Скрипт MATLAB выбрал последовательные данные от Uno Arduino, предварительно обработал его и записал его в.MAT файл. Извлечение Спектральный Live Editor Функций задача затем использовалось, чтобы извлечь спектральные метрики.
Зубчатая передача сервомотора
Сервомотор Futaba S3003 состоит из четырех пар нейлоновых механизмов, как проиллюстрировано в этом рисунке. Шестерня P1 на вале двигателя постоянного тока поймала в сети со ступенчатым механизмом G1. Шестерня P2 является прессованной частью ступенчатого механизма G1 и поймал в сети со ступенчатым механизмом G2. Шестерня P3, который является прессованной частью механизма G2, поймала в сети со ступенчатым механизмом G3. Свяжите P4, который прессуется с G3, сетками с итоговым механизмом G4, который присоединен к выходному сплайну. Ступенчатый механизм устанавливает G1 и P2, G2 и P3, и G3 и P4 являются свободными механизмами вращения то есть, они не присоединены к своим соответствующим валам. Набор механизмов диска обеспечивает 278:1 сокращение, идущее от частоты вращения двигателя 5 414,7 об/мин приблизительно к 19,5 об/мин в выходном сплайне, когда двигатель управляется на уровне 5 вольт. Следующая таблица обрисовывает в общих чертах зубное количество и теоретические значения выходной скорости, частот mesh механизма и совокупного сокращения механизма в каждой mesh механизма.
Предварительная Обработка Данных
Файл servoData.mat содержит два расписания, соответствующие данным из сервомотора. Одно расписание содержит здоровые данные, в то время как второе расписание содержит дефектные данные. Каждый набор данных содержит приблизительно 11 секунд данных, произведенных на уровне 1 500 Гц.
Загрузите данные.
load('servoData.mat', 'healthyData', 'faultyData') healthyData
healthyData=16384×2 timetable
Time MotorCurrent TachoPulse
______________ ____________ __________
0 sec 307.62 1
0.00066667 sec 301.27 1
0.0013333 sec 309.08 1
0.002 sec 315.92 1
0.0026667 sec 304.2 1
0.0033333 sec 311.04 1
0.004 sec 311.52 1
0.0046667 sec 305.18 1
0.0053333 sec 315.43 0
0.006 sec 310.06 0
0.0066667 sec 305.66 0
0.0073333 sec 310.55 0
0.008 sec 304.69 0
0.0086667 sec 310.55 0
0.0093333 sec 310.06 0
0.01 sec 299.8 0
⋮
faultyData
faultyData=16384×2 timetable
Time MotorCurrent TachoPulse
______________ ____________ __________
0 sec 313.48 0
0.00066667 sec 304.2 0
0.0013333 sec 303.22 0
0.002 sec 319.34 0
0.0026667 sec 304.2 0
0.0033333 sec 303.22 0
0.004 sec 319.82 0
0.0046667 sec 303.22 0
0.0053333 sec 306.64 0
0.006 sec 321.29 0
0.0066667 sec 303.71 0
0.0073333 sec 308.11 0
0.008 sec 319.34 0
0.0086667 sec 301.76 0
0.0093333 sec 309.08 0
0.01 sec 319.34 0
⋮
Каждое расписание содержит один столбец с моторным током и один столбец с tacho импульсом от настройки сервомотора. Для того, чтобы визуализировать данные в Извлечении Спектральный Live Editor Функций задача, вычислите спектр мощности моторных текущих данных. Рассмотрите здоровые данные сначала.
fs = 1500;
[healthyMagnitudes,healthyFrequencies] = pwelch(healthyData{:,1},16384,[],[],fs);Используйте Извлечение Спектральный Live Editor Функций задача отобразить данные о спектре мощности на графике. В задаче задайте healthyFrequencies для вектора частоты и healthyMagnitudes для величины спектра мощности.
Анализируйте Peaks спектра мощности с гармоническими диапазонами частот отказа
График спектра мощности текущих данных сервомотора содержит несколько значимых peaks. Можно сопоставить этот peaks с вращающимися валами в настройке сервомотора. Чтобы определить источник различного peaks, добавьте компоненты в Извлечение Спектральный Live Editor Функций задача.
Чтобы добавить компонент, чтобы представлять первый вал вращения в сервомоторе, введите имя для компонента, выберите его тип как Пользовательский, и нажмите Add. Используйте выходные скорости в приведенной выше таблице, чтобы выбрать частоту компонента вала. Выходные скорости были вычислены на основе измеренной скорости выходного вала и известных сокращений механизма настройки.
Частота первого вала составляет 90,24 Гц. После устанавливания значения частоты компонента вала обратите внимание, что диапазон частот отказа перекрывается с одним из peaks в данных о спектре мощности на уровне приблизительно 90 Гц. Можно поэтому сопоставить этот пик в значительной степени с первым валом. Добавление в еще нескольких гармониках основной частоты первого вала создает больше диапазонов частот отказа, которые перекрывают другой peaks в данных. Гармонические диапазоны частот сосредоточены вокруг целочисленных множителей основной частоты и могут все еще быть сопоставлены с тем же компонентом. Установите гармоники первого вала быть вектором [1 2 3 4 5 6] так, чтобы полосы отказа, распространенные по большей части частотного диапазона. Первый вал, вращается на самой высокой частоте, таким образом, peaks спектра мощности на более высоких частотах происходит из-за гармоник этого вала.
С учетом меньшего peaks спектра мощности, такого как те приблизительно 13 Гц или 29 Гц, добавьте в компоненте для второго вала вращения. Этот компонент является также собственным компонентом, и его основная частота составляет 14,56 Гц. Вы не должны добавлять как много гармоник для второго вала, поскольку большинство более высоких пиковых частот в основном составляется гармониками первого вала. Установите гармоники второго вала быть вектором [1 2 3 4]. Первые, вторые, и четвертые гармоники этой частоты вала выравниваются приятно с peaks в графике спектра мощности. Однако, поскольку третья гармоника является менее видной в данных, вы не должны включать эту гармонику. Измените гармоники второго вала, чтобы быть вектором [1 2 4].
Так же к второму валу, добавьте компонент для третьего вала вращения. Основная частота этого компонента составляет 2,91 Гц, как замечено в таблице. Начните с первых четырех гармоник определять, выравниваются ли они с каким-либо видным peaks в данных. Обратите внимание на то, что третья гармоника третьего вала совпадает со скачком спектра мощности приблизительно 8 Гц. Другие гармоники являются менее видными и могут быть удалены. Может иметь место, что разрешение частоты спектра мощности не позволяет отличать более низкие частоты. Установите гармоники третьего компонента вала быть только третьей гармоникой.
Поскольку выходные скорости остающихся валов являются также низкими частотами, которые не могут быть различимы от разрешения частоты спектра мощности, добавление, что компоненты для этих валов не необходимы, чтобы анализировать главный peaks моторных текущих данных.
Анализируйте Peaks боковой полосы
Путем изменения масштаба в график в задаче вы видите, что данные о спектре мощности содержат peaks стороны рядом с частью главного peaks. Например, меньший peaks стороны приблизительно 76 Гц и 104 Гц окружает пик на уровне 90 Гц. Этот peaks, вероятно, сопоставлен с боковыми полосами первого компонента вала. Боковые полосы вызываются вторым связанным источником частоты, влияющим на первичный гармонический источник частоты. С настройкой сервомотора это наблюдение приводит к предположению, что боковые полосы для каждого вала вызываются следующим валом в зубчатой передаче.
Отредактируйте первые два компонента вала, чтобы включать их первую боковую полосу. Для первого вала разделительное значение боковой полосы должно быть равно номинальной выходной частоте второго вала, 14,56 Гц.
Так же для второго вала разделительное значение боковой полосы должно быть равно номинальной частоте третьего вала, 2,91 Гц. Увеличивание масштаб снова на графике показывает, что многие новые боковые полосы перекрываются хорошо с peaks стороны в данных. Это легче видеть в низких частотах такой как между 0 Гц и 120 Гц.
Извлеките спектральные метрики, чтобы обнаружить отказы
Извлечение Спектральный Live Editor Функций задача генерирует различные спектральные метрики данных о спектре мощности в частотных диапазонах отказа. Для каждого диапазона частот отказа вычислите пиковую амплитуду, пиковую частоту и мощность полосы наряду с общей мощностью полосы всех диапазонов частот отказа.
load('dataSample.mat')
spectralMetrics_healthyspectralMetrics_healthy=1×85 table
PeakAmplitude1 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 PeakAmplitude16 PeakFrequency16 BandPower16 PeakAmplitude17 PeakFrequency17 BandPower17 PeakAmplitude18 PeakFrequency18 BandPower18 PeakAmplitude19 PeakFrequency19 BandPower19 PeakAmplitude20 PeakFrequency20 BandPower20 PeakAmplitude21 PeakFrequency21 BandPower21 PeakAmplitude22 PeakFrequency22 BandPower22 PeakAmplitude23 PeakFrequency23 BandPower23 PeakAmplitude24 PeakFrequency24 BandPower24 PeakAmplitude25 PeakFrequency25 BandPower25 PeakAmplitude26 PeakFrequency26 BandPower26 PeakAmplitude27 PeakFrequency27 BandPower27 PeakAmplitude28 PeakFrequency28 BandPower28 TotalBandPower
______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________
0.088268 77.179 0.031107 2.0259 91.37 1.4892 0.026872 104.74 0.017091 0.0040207 165.62 0.0040427 0.063531 178.99 0.088016 0.0079751 193.54 0.0093892 0.07461 254.61 0.080139 0.60862 270.72 0.42026 0.098992 284.09 0.08489 0.0068446 347.99 0.0072084 0.01066 362.09 0.013542 0.0059992 375.46 0.006807 0.055461 438.45 0.052725 0.071491 449.16 0.089713 0.0569 467.1 0.0561 2.0596 525.97 2.513 1.9851 539.52 2.3441 0.56849 554.81 0.64988 0.41738 11.627 0.083331 3.8511 14.557 0.67359 0.41719 17.487 0.060687 0.080245 26.184 0.03093 4.2605 29.114 0.70964 0.021784 32.043 0.0082336 0.018917 55.664 0.0047859 0.036751 58.594 0.01207 0.0087469 61.432 0.0022404 1.9068 8.6975 0.31922 9.862
Эти метрики могут оказаться полезными в обнаружении отказов в настройке сервомотора. Существенные изменения в данных о спектре мощности часто указывают, что некоторый компонент изменяется или перестал работать. Если существует сдвиг в одной из пиковых частот или если пиковая амплитуда понижается значительно в зависимости от времени, то это может быть знаком отказа.
Чтобы исследовать этот сценарий, вычислите данные о спектре мощности для дефектного набора данных.
[faultyMagnitudes,faultyFrequencies] = pwelch(faultyData{:,1},16384,[],[],fs);Постройте спектр дефектных данных в задаче и настройте основные частоты и разделительные значения боковой полосы компонентов вала на основе измеренной выходной скорости настройки сервомотора.
Fs = 1500; % 1500 Hz [outputSpeed,t] = tachorpm(faultyData.TachoPulse,Fs,'PulsesPerRev',16,'FitType','linear'); meanOutputSpeed = mean(outputSpeed)/60 % convert from rpm to Hz
meanOutputSpeed = 0.3150
shaft4Speed = meanOutputSpeed * 41 / 16 % 16 pinion teeth, 41 gear teethshaft4Speed = 0.8072
shaft3Speed = shaft4Speed * 35 / 10 % 10 pinion teeth, 35 gear teethshaft3Speed = 2.8251
shaft2Speed = shaft3Speed * 50 / 10 % 10 pinion teeth, 50 gear teethshaft2Speed = 14.1254
shaft1Speed = shaft2Speed * 62 / 10 % 10 pinion teeth, 62 gear teethshaft1Speed = 87.5772
Для первого компонента вала используйте shaft1Speed как основная частота и shaft2Speed как разделение боковой полосы. Для второго компонента вала используйте shaft2Speed как основная частота и shaft3Speed как разделение боковой полосы. Для третьего компонента вала используйте shaft3Speed как основная частота.
Как видно в визуализации дефектных данных о спектре мощности, несколько из peaks уменьшились в величине. Например, пик, выровненный в здоровом наборе данных со второй гармоникой первого вала приблизительно 180 Гц, почти незначителен в дефектном наборе данных. Поскольку было ранее определено, что этот пик, вероятно, сопоставлен с первым валом, это указывает на потенциальный отказ в первом вале. Дальнейшее исследование спектральной метрической таблицы может обеспечить более подробную информацию о пиковых частотах, пиковых амплитудах и мощностях полосы.
spectralMetrics_faulty
spectralMetrics_faulty=1×85 table
PeakAmplitude1 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 PeakAmplitude16 PeakFrequency16 BandPower16 PeakAmplitude17 PeakFrequency17 BandPower17 PeakAmplitude18 PeakFrequency18 BandPower18 PeakAmplitude19 PeakFrequency19 BandPower19 PeakAmplitude20 PeakFrequency20 BandPower20 PeakAmplitude21 PeakFrequency21 BandPower21 PeakAmplitude22 PeakFrequency22 BandPower22 PeakAmplitude23 PeakFrequency23 BandPower23 PeakAmplitude24 PeakFrequency24 BandPower24 PeakAmplitude25 PeakFrequency25 BandPower25 PeakAmplitude26 PeakFrequency26 BandPower26 PeakAmplitude27 PeakFrequency27 BandPower27 PeakAmplitude28 PeakFrequency28 BandPower28 TotalBandPower
______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________
0.0015834 71.411 0.0017557 0.082928 88.806 0.070886 0.0023192 101.62 0.0028621 0.00094302 161.77 0.0010918 0.0030192 174.68 0.0024061 0.00071009 190.98 0.0011196 0.009749 248.84 0.0121 0.071268 261.84 0.077265 0.022904 274.84 0.011766 0.0039476 334.9 0.0043915 0.008027 348.91 0.010401 0.0032387 362.82 0.0032159 0.0044798 422.24 0.0065658 0.010303 436.98 0.013423 0.0062224 451.45 0.008025 3.7506 511.23 3.6514 4.8463 526.89 4.1478 2.4573 540.07 1.8774 1.8378 11.353 0.33072 1.085 14.099 0.4306 0.1326 16.937 0.021304 0.051885 25.452 0.012906 1.6184 28.29 0.25581 0.018916 31.128 0.0052321 0.0043149 53.741 0.0012379 0.058117 56.58 0.02 0.0057658 59.052 0.0022057 3.0935 8.5144 0.49943 11.483
Как альтернатива обновлению спектральных данных в задаче Live Editor, можно также использовать автоматически сгенерированный код MATLAB для detemine спектральные метрики дефектных данных. Код ниже был автоматически сгенерирован, когда задача Live Editor использовалась, чтобы сгенерировать здоровые спектральные метрики. Выполните код.
% Generate the fault bands and information for each component [FB_Shaft1, info_Shaft1] = faultBands(90.24, 1:6, 14.56, 0:1); [FB_Shaft2, info_Shaft2] = faultBands(14.56, [1 2 4], 2.91, 0:1); [FB_Shaft3, info_Shaft3] = faultBands(2.91, 3); % Combine the fault bands of each component FB_healthy = [FB_Shaft1; ... FB_Shaft2; ... FB_Shaft3]; % Combine the information regarding the fault bands of each component info_healthy.Centers = [info_Shaft1.Centers, ... info_Shaft2.Centers, ... info_Shaft3.Centers]; info_healthy.Labels = [info_Shaft1.Labels, ... info_Shaft2.Labels, ... info_Shaft3.Labels]; info_healthy.FaultGroups = [info_Shaft1.HarmonicGroups, ... info_Shaft2.HarmonicGroups, ... info_Shaft3.HarmonicGroups]; % Clear temporary outputs from the workspace clear FB_Shaft1 info_Shaft1; clear FB_Shaft2 info_Shaft2; clear FB_Shaft3 info_Shaft3; % Compute fault band metrics of the power spectrum healthyMagnitudes spectralMetrics_healthy = faultBandMetrics(healthyMagnitudes, healthyFrequencies, FB_healthy)
spectralMetrics_healthy=1×85 table
PeakAmplitude1 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 PeakAmplitude16 PeakFrequency16 BandPower16 PeakAmplitude17 PeakFrequency17 BandPower17 PeakAmplitude18 PeakFrequency18 BandPower18 PeakAmplitude19 PeakFrequency19 BandPower19 PeakAmplitude20 PeakFrequency20 BandPower20 PeakAmplitude21 PeakFrequency21 BandPower21 PeakAmplitude22 PeakFrequency22 BandPower22 PeakAmplitude23 PeakFrequency23 BandPower23 PeakAmplitude24 PeakFrequency24 BandPower24 PeakAmplitude25 PeakFrequency25 BandPower25 PeakAmplitude26 PeakFrequency26 BandPower26 PeakAmplitude27 PeakFrequency27 BandPower27 PeakAmplitude28 PeakFrequency28 BandPower28 TotalBandPower
______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________
0.088268 77.179 0.031107 2.0259 91.37 1.4892 0.026872 104.74 0.017091 0.0040207 165.62 0.0040427 0.063531 178.99 0.088016 0.0079751 193.54 0.0093892 0.07461 254.61 0.080139 0.60862 270.72 0.42026 0.098992 284.09 0.08489 0.0068446 347.99 0.0072084 0.01066 362.09 0.013542 0.0059992 375.46 0.006807 0.055461 438.45 0.052725 0.071491 449.16 0.089713 0.0569 467.1 0.0561 2.0596 525.97 2.513 1.9851 539.52 2.3441 0.56849 554.81 0.64988 0.41738 11.627 0.083331 3.8511 14.557 0.67359 0.41719 17.487 0.060687 0.080245 26.184 0.03093 4.2605 29.114 0.70964 0.021784 32.043 0.0082336 0.018917 55.664 0.0047859 0.036751 58.594 0.01207 0.0087469 61.432 0.0022404 1.9068 8.6975 0.31922 9.862
Этот код может легко быть настроен для нового дефектного набора данных.
% Generate the fault bands and information for each component [FB_Shaft1, info_Shaft1] = faultBands(shaft1Speed, 1:6, shaft2Speed, 0:1); [FB_Shaft2, info_Shaft2] = faultBands(shaft2Speed, [1 2 4], shaft3Speed, 0:1); [FB_Shaft3, info_Shaft3] = faultBands(shaft3Speed, 3); % Combine the fault bands of each component FB_faulty = [FB_Shaft1; ... FB_Shaft2; ... FB_Shaft3]; % Combine the information regarding the fault bands of each component info_faulty.Centers = [info_Shaft1.Centers, ... info_Shaft2.Centers, ... info_Shaft3.Centers]; info_faulty.Labels = [info_Shaft1.Labels, ... info_Shaft2.Labels, ... info_Shaft3.Labels]; info_faulty.FaultGroups = [info_Shaft1.HarmonicGroups, ... info_Shaft2.HarmonicGroups, ... info_Shaft3.HarmonicGroups]; % Clear temporary outputs from the workspace clear FB_Shaft1 info_Shaft1; clear FB_Shaft2 info_Shaft2; clear FB_Shaft3 info_Shaft3; % Compute fault band metrics of the power spectrum healthyMagnitudes spectralMetrics_faulty = faultBandMetrics(faultyMagnitudes, faultyFrequencies, FB_faulty)
spectralMetrics_faulty=1×85 table
PeakAmplitude1 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 PeakAmplitude16 PeakFrequency16 BandPower16 PeakAmplitude17 PeakFrequency17 BandPower17 PeakAmplitude18 PeakFrequency18 BandPower18 PeakAmplitude19 PeakFrequency19 BandPower19 PeakAmplitude20 PeakFrequency20 BandPower20 PeakAmplitude21 PeakFrequency21 BandPower21 PeakAmplitude22 PeakFrequency22 BandPower22 PeakAmplitude23 PeakFrequency23 BandPower23 PeakAmplitude24 PeakFrequency24 BandPower24 PeakAmplitude25 PeakFrequency25 BandPower25 PeakAmplitude26 PeakFrequency26 BandPower26 PeakAmplitude27 PeakFrequency27 BandPower27 PeakAmplitude28 PeakFrequency28 BandPower28 TotalBandPower
______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________
0.0015834 71.411 0.0017557 0.082928 88.806 0.070886 0.0023192 101.62 0.0028621 0.00094302 161.77 0.0010918 0.0030192 174.68 0.0024061 0.00071009 190.98 0.0011196 0.009749 248.84 0.0121 0.071268 261.84 0.077265 0.022904 274.84 0.011766 0.0039476 334.9 0.0043915 0.008027 348.91 0.010401 0.0032387 362.82 0.0032159 0.0044798 422.24 0.0065658 0.010303 436.98 0.013423 0.0062224 451.45 0.008025 3.7506 511.23 3.6514 4.8463 526.89 4.1478 2.4573 540.07 1.8774 1.8378 11.353 0.33072 1.085 14.099 0.4306 0.1326 16.937 0.021304 0.051885 25.452 0.012906 1.6184 28.29 0.25581 0.018916 31.128 0.0052321 0.0043149 53.741 0.0012379 0.058117 56.58 0.02 0.0057658 59.052 0.0022057 3.0935 8.5144 0.49943 11.483
Объединение здоровых и дефектных спектральных метрических таблиц улучшает видимость в изменения в пиковых амплитудах и мощностях полосы спектральных данных.
spectralMetrics_total = [spectralMetrics_healthy; spectralMetrics_faulty]
spectralMetrics_total=2×85 table
PeakAmplitude1 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 PeakAmplitude16 PeakFrequency16 BandPower16 PeakAmplitude17 PeakFrequency17 BandPower17 PeakAmplitude18 PeakFrequency18 BandPower18 PeakAmplitude19 PeakFrequency19 BandPower19 PeakAmplitude20 PeakFrequency20 BandPower20 PeakAmplitude21 PeakFrequency21 BandPower21 PeakAmplitude22 PeakFrequency22 BandPower22 PeakAmplitude23 PeakFrequency23 BandPower23 PeakAmplitude24 PeakFrequency24 BandPower24 PeakAmplitude25 PeakFrequency25 BandPower25 PeakAmplitude26 PeakFrequency26 BandPower26 PeakAmplitude27 PeakFrequency27 BandPower27 PeakAmplitude28 PeakFrequency28 BandPower28 TotalBandPower
______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ ______________ ______________ __________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ _______________ _______________ ___________ ______________
0.088268 77.179 0.031107 2.0259 91.37 1.4892 0.026872 104.74 0.017091 0.0040207 165.62 0.0040427 0.063531 178.99 0.088016 0.0079751 193.54 0.0093892 0.07461 254.61 0.080139 0.60862 270.72 0.42026 0.098992 284.09 0.08489 0.0068446 347.99 0.0072084 0.01066 362.09 0.013542 0.0059992 375.46 0.006807 0.055461 438.45 0.052725 0.071491 449.16 0.089713 0.0569 467.1 0.0561 2.0596 525.97 2.513 1.9851 539.52 2.3441 0.56849 554.81 0.64988 0.41738 11.627 0.083331 3.8511 14.557 0.67359 0.41719 17.487 0.060687 0.080245 26.184 0.03093 4.2605 29.114 0.70964 0.021784 32.043 0.0082336 0.018917 55.664 0.0047859 0.036751 58.594 0.01207 0.0087469 61.432 0.0022404 1.9068 8.6975 0.31922 9.862
0.0015834 71.411 0.0017557 0.082928 88.806 0.070886 0.0023192 101.62 0.0028621 0.00094302 161.77 0.0010918 0.0030192 174.68 0.0024061 0.00071009 190.98 0.0011196 0.009749 248.84 0.0121 0.071268 261.84 0.077265 0.022904 274.84 0.011766 0.0039476 334.9 0.0043915 0.008027 348.91 0.010401 0.0032387 362.82 0.0032159 0.0044798 422.24 0.0065658 0.010303 436.98 0.013423 0.0062224 451.45 0.008025 3.7506 511.23 3.6514 4.8463 526.89 4.1478 2.4573 540.07 1.8774 1.8378 11.353 0.33072 1.085 14.099 0.4306 0.1326 16.937 0.021304 0.051885 25.452 0.012906 1.6184 28.29 0.25581 0.018916 31.128 0.0052321 0.0043149 53.741 0.0012379 0.058117 56.58 0.02 0.0057658 59.052 0.0022057 3.0935 8.5144 0.49943 11.483
Например, если вы смотрите на PeakAmplitude2 в таблице амплитуда пика спектра мощности понижается с 2,0259 до 0,0829. Используя PeakFrequency2 значение, вы знаете, что это отбрасывание происходит приблизительно 90 Гц. Постройте эти два спектра мощности на тех же осях, чтобы визуализировать отбрасывание за пределами задачи Live Editor.
plot(healthyFrequencies,10*log10(healthyMagnitudes),'b-'); % plot in decibels hold on; plot(faultyFrequencies,10*log10(faultyMagnitudes),'r-'); % plot in decibels legend('Healthy','Faulty') xlabel('Frequency (Hz)') ylabel('Power Spectrum (dB)') title('Comparison of Power Spectrums') hold off;
Когда метрическая таблица показала, пиковые отбрасывания на приблизительно 90 Гц значительно в амплитуде. Чтобы определить, какая частота компонента вызвала это, перепроверьте в предыдущем Извлечении Спектральный Live Editor Функции задача.
Полоса отказа приблизительно 90 Гц является первой гармоникой первого вала вращения. Таким образом вы знаете, что что-то изменяется в этом вале, и он может отклоняться к отказу.
Графический вывод здоровых и дефектных спектров мощности вместе может быть полезным способом подсветить изменения в пиковых амплитудах. В дополнение к пиковым приблизительно 90 Гц для первой гармоники первого вала другие значительные уменьшения в пиковых амплитудах могут рассматриваться такой как со второй гармоникой этого вала приблизительно 180 Гц. Этот пик чрезвычайно не существует в дефектном наборе данных.
Пиковые амплитуды из здоровых и дефектных данных могут также быть сравнены с помощью столбчатой диаграммы.
PeakFrequencies = spectralMetrics_total(:,2:3:end-1).Variables'; PeakAmplitudes = spectralMetrics_total(:,1:3:end-1).Variables'; bar(PeakFrequencies, PeakAmplitudes); legend('Healthy','Faulty') xlabel('Frequency (Hz)') ylabel('Peak Amplitude') title('Peak Amplitudes of Healthy and Faulty Power Spectrum Data')
Увеличьте масштаб, чтобы видеть изменение в пиковой амплитуде в первой гармонике первого вала вращения.
xlim([87 93]) ylim([0 0.7])
Использование извлечения спектральный Live Editor функций задача
Как показано в этом примере Извлечение Спектральный Live Editor Функций задача может оказаться полезной для нескольких различных приложений. С задачей Live Editor можно легко совпадать со спектральным peaks известным частотам компонента машины. Это помогает вам лучше изучить свои данные и механические компоненты, вызывающие различные функции в данных.
Другое приложение Извлечения Спектральный Live Editor Функций, задача состоит в том, чтобы сгенерировать метрики, которые характеризуют ваши спектральные данные в частотных диапазонах интереса. Задача производит выходную таблицу, содержащую пиковую амплитуду, пиковую частоту, и мощность полосы каждого диапазона частот отказа, а также общую мощность полосы всех полос отказа. Однако эти метрики характерны для ввода данных спектра мощности в задаче.
Чтобы расширить это использование так, чтобы можно было отслеживать метрики в зависимости от времени, как, новые наборы данных собраны, можно или обновить данные о спектре мощности в задаче или использовать автоматически сгенерированный код MATLAB, чтобы произвести метрическую таблицу. Копирование сгенерированного кода MATLAB является простым способом продолжить вычислять метрики полосы отказа для многих новых наборов данных.
Третье использование задачи комбинирует преимущества двух ранее обсужденного использования. Поскольку задача сопоставляет различные механические компоненты с peaks в спектральных данных, можно быстро определить, какие компоненты вызывают существенные изменения в спектральных данных и таким образом потенциальных отказах. Например, один общий индикатор аппаратного сбоя, выведенного из спектральных данных, является изменением в амплитуде спектрального peaks. В спектральной метрической таблице, если вы замечаете значительное понижение пиковой амплитуды или мощности полосы в зависимости от времени, можно проследить соответствующую пиковую частоту до графика видеть, какой диапазон частот отказа компонента выравнивается с тем пиком.
Смотрите также
faultBands | bearingFaultBands | gearMeshFaultBands | faultBandMetrics | Извлеките спектральные функции