Анализируйте данные о зубчатой передаче и извлеките спектральные функции Используя задачи 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 Функций задача, вычислите спектр мощности моторных текущих данных. Рассмотрите здоровые данные сначала.

[healthyMagnitudes, healthyFrequencies] = pspectrum(healthyData.MotorCurrent, healthyData.Time);

Используйте Извлечение Спектральный 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('sampleData.mat')
spectralMetrics_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.018807           77.106         0.044029        0.66383            90.842          2.2909        0.0077442           104.58         0.026416        0.001917           168.13        0.0076742        0.038105           182.6          0.13883        0.0034982           192.86         0.013928        0.026813           258.24         0.11464          0.192             270.7          0.68323         0.033756           284.07         0.12347         0.0026502           347.62          0.0099213        0.0048213           360.81          0.020783         0.0022192           375.64          0.0086095        0.018201            435.53          0.074231         0.041176            449.08           0.12665         0.018558            467.22          0.069264          0.96719            526.01           3.7322           0.91169            539.19           3.1072           0.29472            558.24           1.0803           1.0341             11.905           0.75471          4.2117             14.286           3.0372          0.076119            17.399          0.078405         0.037845            26.557          0.038047          0.75325            29.121           0.63976         0.010537            32.234          0.010985         0.0058578           55.311          0.006089         0.014194            58.425          0.011771         0.0033205           61.355          0.0034375         0.51023            8.7912           0.27615          16.529    

Эти метрики могут оказаться полезными в обнаружении отказов в настройке сервомотора. Существенные изменения в данных о спектре мощности часто указывают, что некоторый компонент изменяется или перестал работать. Если существует сдвиг в одной из пиковых частот или если пиковая амплитуда понижается значительно в зависимости от времени, то это может быть знаком отказа.

Чтобы исследовать этот сценарий, вычислите данные о спектре мощности для дефектного набора данных.

[faultyMagnitudes, faultyFrequencies] = pspectrum(faultyData.MotorCurrent, faultyData.Time);

Постройте спектр дефектных данных в задаче и настройте основные частоты и разделительные значения боковой полосы компонентов вала на основе измеренной выходной скорости настройки сервомотора.

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 teeth
shaft4Speed = 0.8072
shaft3Speed = shaft4Speed * 35 / 10 % 10 pinion teeth, 35 gear teeth
shaft3Speed = 2.8251
shaft2Speed = shaft3Speed * 50 / 10 % 10 pinion teeth, 50 gear teeth
shaft2Speed = 14.1254
shaft1Speed = shaft2Speed * 62 / 10 % 10 pinion teeth, 62 gear teeth
shaft1Speed = 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.0011035           75.641        0.0035389        0.032588           88.095         0.095866       0.00095266          101.47        0.0039494       0.00041734          158.97        0.0016878       0.0010561           174.73        0.0038409       0.0004178           190.11        0.0018162       0.0051797           249.27         0.017516        0.028642           261.54          0.1203        0.0051032           274.91         0.015495        0.0014915           336.08          0.0062556        0.0034171           350.37          0.015104         0.0014217           363.19          0.0049667        0.0022728           423.63          0.0092141        0.0058438           437.36          0.020561         0.0031571           450.37          0.010481          1.9389             511.36           6.3414           1.3398             523.81           5.1905           0.79257            539.38           2.9432           0.9969             11.538           0.82388           2.728             13.919           1.9338          0.030943            17.033          0.032899         0.019372            25.641          0.018553          0.27041            28.205           0.26875         0.0056339           31.136          0.0060698        0.0029558           54.029          0.0028722        0.024586            56.593          0.025098         0.0027373           58.974          0.0027412         0.74618            8.4249           0.40386          18.324    

Как альтернатива обновлению спектральных данных в задаче 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.018807           77.106         0.044029        0.66383            90.842          2.2909        0.0077442           104.58         0.026416        0.001917           168.13        0.0076742        0.038105           182.6          0.13883        0.0034982           192.86         0.013928        0.026813           258.24         0.11464          0.192             270.7          0.68323         0.033756           284.07         0.12347         0.0026502           347.62          0.0099213        0.0048213           360.81          0.020783         0.0022192           375.64          0.0086095        0.018201            435.53          0.074231         0.041176            449.08           0.12665         0.018558            467.22          0.069264          0.96719            526.01           3.7322           0.91169            539.19           3.1072           0.29472            558.24           1.0803           1.0341             11.905           0.75471          4.2117             14.286           3.0372          0.076119            17.399          0.078405         0.037845            26.557          0.038047          0.75325            29.121           0.63976         0.010537            32.234          0.010985         0.0058578           55.311          0.006089         0.014194            58.425          0.011771         0.0033205           61.355          0.0034375         0.51023            8.7912           0.27615          16.529    

Этот код может легко быть настроен для нового дефектного набора данных.

% 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.0011035           75.641        0.0035389        0.032588           88.095         0.095866       0.00095266          101.47        0.0039494       0.00041734          158.97        0.0016878       0.0010561           174.73        0.0038409       0.0004178           190.11        0.0018162       0.0051797           249.27         0.017516        0.028642           261.54          0.1203        0.0051032           274.91         0.015495        0.0014915           336.08          0.0062556        0.0034171           350.37          0.015104         0.0014217           363.19          0.0049667        0.0022728           423.63          0.0092141        0.0058438           437.36          0.020561         0.0031571           450.37          0.010481          1.9389             511.36           6.3414           1.3398             523.81           5.1905           0.79257            539.38           2.9432           0.9969             11.538           0.82388           2.728             13.919           1.9338          0.030943            17.033          0.032899         0.019372            25.641          0.018553          0.27041            28.205           0.26875         0.0056339           31.136          0.0060698        0.0029558           54.029          0.0028722        0.024586            56.593          0.025098         0.0027373           58.974          0.0027412         0.74618            8.4249           0.40386          18.324    

Объединение здоровых и дефектных спектральных метрических таблиц улучшает видимость в изменения в пиковых амплитудах и мощностях полосы спектральных данных.

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.018807           77.106         0.044029         0.66383           90.842           2.2909        0.0077442          104.58         0.026416         0.001917          168.13        0.0076742        0.038105            182.6          0.13883       0.0034982           192.86         0.013928        0.026813           258.24          0.11464           0.192            270.7         0.68323         0.033756           284.07          0.12347        0.0026502           347.62          0.0099213        0.0048213           360.81          0.020783         0.0022192           375.64          0.0086095         0.018201           435.53           0.074231         0.041176           449.08           0.12665          0.018558           467.22          0.069264          0.96719            526.01           3.7322           0.91169            539.19           3.1072           0.29472            558.24           1.0803           1.0341             11.905           0.75471          4.2117             14.286           3.0372          0.076119            17.399          0.078405         0.037845            26.557          0.038047          0.75325            29.121           0.63976          0.010537           32.234           0.010985        0.0058578           55.311           0.006089        0.014194            58.425          0.011771         0.0033205           61.355          0.0034375         0.51023            8.7912           0.27615          16.529    
      0.0011035           75.641        0.0035389        0.032588           88.095         0.095866       0.00095266          101.47        0.0039494       0.00041734          158.97        0.0016878       0.0010561           174.73        0.0038409       0.0004178           190.11        0.0018162       0.0051797           249.27         0.017516        0.028642           261.54          0.1203        0.0051032           274.91         0.015495        0.0014915           336.08          0.0062556        0.0034171           350.37          0.015104         0.0014217           363.19          0.0049667        0.0022728           423.63          0.0092141        0.0058438           437.36          0.020561         0.0031571           450.37          0.010481           1.9389            511.36           6.3414            1.3398            523.81           5.1905           0.79257            539.38           2.9432           0.9969             11.538           0.82388           2.728             13.919           1.9338          0.030943            17.033          0.032899         0.019372            25.641          0.018553          0.27041            28.205           0.26875         0.0056339           31.136          0.0060698        0.0029558           54.029          0.0028722        0.024586            56.593          0.025098         0.0027373           58.974          0.0027412         0.74618            8.4249           0.40386          18.324    

Например, если вы смотрите на PeakAmplitude2 в таблице, амплитуда пика спектра мощности понижается с 0,6638 до 0,0326. Используя значение 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;

Figure contains an axes. The axes with title Comparison of Power Spectrums contains 2 objects of type line. These objects represent Healthy, Faulty.

Когда метрическая таблица показала, пиковые отбрасывания на приблизительно 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])

Figure contains an axes. The axes with title Peak Amplitudes of Healthy and Faulty Power Spectrum Data contains 2 objects of type bar. These objects represent Healthy, Faulty.

Использование извлечения спектральный Live Editor функций задача

Как показано в этом примере Извлечение Спектральный Live Editor Функций задача может оказаться полезной для нескольких различных приложений. С задачей Live Editor можно легко совпадать со спектральным peaks известным частотам компонента машины. Это помогает вам лучше изучить свои данные и механические компоненты, вызывающие различные функции в данных.

Другое приложение Извлечения Спектральный Live Editor Функций, задача состоит в том, чтобы сгенерировать метрики, которые характеризуют ваши спектральные данные в частотных диапазонах интереса. Задача производит выходную таблицу, содержащую пиковую амплитуду, пиковую частоту, и мощность полосы каждого диапазона частот отказа, а также общую мощность полосы всех полос отказа. Однако эти метрики характерны для ввода данных спектра мощности в задаче.

Чтобы расширить это использование так, чтобы можно было отслеживать метрики в зависимости от времени, как, новые наборы данных собраны, можно или обновить данные о спектре мощности в задаче или использовать автоматически сгенерированный код MATLAB, чтобы произвести метрическую таблицу. Копирование сгенерированного кода MATLAB является простым способом продолжить вычислять метрики полосы отказа для многих новых наборов данных.

Третье использование задачи комбинирует преимущества двух ранее обсужденного использования. Поскольку задача сопоставляет различные механические компоненты с peaks в спектральных данных, можно быстро определить, какие компоненты вызывают существенные изменения в спектральных данных и таким образом потенциальных отказах. Например, один общий индикатор аппаратного сбоя, выведенного из спектральных данных, является изменением в амплитуде спектрального peaks. В спектральной метрической таблице, если вы замечаете значительное понижение пиковой амплитуды или мощности полосы в зависимости от времени, можно проследить соответствующую пиковую частоту до графика видеть, какой диапазон частот отказа компонента выравнивается с тем пиком.

Смотрите также

| | | |

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте