gearMeshFaultBands

Создайте диапазоны частот вокруг характеристических частот отказа запутывающих механизмов для спектрального извлечения признаков

свернуть все на странице

Синтаксис

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

Описание

пример

FB = gearMeshFaultBands(FR,Ni,No) генерирует характеристические диапазоны частот отказа FB из mesh механизма с помощью скорости вращения входного механизма FR и количество зубов на входе Ni и выходной механизм No соответственно. Значения в FB имейте те же неявные модули как FR

пример

FB = gearMeshFaultBands(___,Name,Value) позволяет вам задавать дополнительные параметры с помощью одного или нескольких аргументов пары "имя-значение".

пример

[FB,info] = gearMeshFaultBands(___) также возвращает структуру info содержа информацию о сгенерированных диапазонах частот отказа FB.

пример

gearMeshFaultBands(___) без выходных аргументов строит столбчатую диаграмму сгенерированных диапазонов частот отказа FB.

Примеры

свернуть все

Скрипт Open Live Script

В данном примере рассмотрите простой набор механизма с 8-зубой шестерней на входном вале, сцепляющемся с 42-зубым цилиндрическим механизмом на выходном вале. Примите, что входной вал вращается на уровне 20 об/мин. Создайте диапазоны частот mesh механизма с помощью физических характеристик набора механизма.

Ni = 8;
No = 42;
FR = 20;
[FB,info] = gearMeshFaultBands(FR,Ni,No)
FB = 5×2

   19.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]

FB 5x2 массив, который включает первичные частоты 1Fi, 1Fo, 1Fa и 1Fm соответственно. Информация о структуре содержит центральные частоты и метки каждого частотного диапазона в FB.

Скрипт Open Live Script

В данном примере рассмотрите простой набор механизма с 8-зубой шестерней на входном вале, сцепляющемся с 42-зубым цилиндрическим механизмом на выходном вале. Примите, что входной вал управляется на уровне 20 Гц. Набор данных motorSignal.mat содержит данные о вибрации для mesh механизма, произведенной на уровне 1 500 Гц.

Во-первых, создайте диапазоны частот mesh механизма с помощью физических характеристик набора механизма. Создайте диапазоны частот с первыми 3 боковыми полосами и задайте 'Domain'как 'order'.

Ni = 8;
No = 42;
FR = 20;
FB = gearMeshFaultBands(FR,Ni,No,'Sidebands',1:3,'Domain','order')
FB = 15×2

    0.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
      ⋮

FB 15x2 массив, который включает первичные частоты и их боковые полосы.

Загрузите данные о вибрации и вычислите PSD и сетку частоты с помощью pspectrum. Используйте разрешение частоты 0,5.

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

Теперь используйте диапазоны частот и данные PSD, чтобы вычислить спектральные метрики. 

spectralMetrics = faultBandMetrics(psd,freqGrid,FB)
spectralMetrics=1×46 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    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

spectralMetrics 1x46 таблица с пиковой амплитудой, пиковой частотой и мощностью полосы, вычисленной для каждого частотного диапазона в FB. Последний столбец в spectralMetrics общая мощность полосы, вычисленная через все 15 частот в FB.

Скрипт Open Live Script

В данном примере рассмотрите простой набор шестерни и механизма с входной скоростью вала 1 800 об/мин. Полагание, что шестерня на входном вале имеет 6 зубов и механизм на выходном вале, имеет 8 зубов, визуализируйте диапазоны частот для mesh механизма.

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

Из графика наблюдайте следующее:

  • Выведите частоту дефекта вала, 1Fo на уровне 1 350 Гц

  • Введите частоту дефекта вала, 1Fi на уровне 1 800 Гц

  • Частота дефекта фазы Assembly, 1Fa на уровне 5 400 Гц

  • Частота дефекта mesh механизма, 1Fm на уровне 10 800 Гц

Входные параметры

свернуть все

Скорость вращения входного механизма в виде положительной скалярной величины. FR основная частота вокруг который gearMeshFaultBands генерирует диапазоны частот отказа. Задайте FR или в Герц или в оборотах в минуту.

Количество зубов на входном механизме в виде положительного целого числа.

Количество зубов на выходном механизме в виде положительного целого числа.

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: ...,'Harmonics',[1,3,5]

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

Боковые полосы вокруг основной частоты и ее гармоник, которые будут включены в виде разделенной запятой пары, состоящей из 'Sidebands'и вектор неотрицательных целых чисел. Значение по умолчанию 0. Задайте 'Sidebands'когда это необходимо, чтобы создать диапазоны частот с боковыми полосами вокруг основной частоты и ее гармоник.

Ширина диапазонов частот, сосредоточенных на номинальных частотах отказа в виде разделенной запятой пары, состоящей из 'Width'и положительная скалярная величина. Значение по умолчанию составляет 10 процентов основной частоты. Постарайтесь не задавать 'Width'с большим значением так, чтобы полосы отказа не перекрывались.

Единицы частот полосы отказа в виде разделенной запятой пары, состоящей из 'Domain'и любой 'frequency' или 'order'. Выбор:

  • 'frequency' если вы хотите FB быть возвращенным в тех же модулях как FR.

  • 'order' если вы хотите FB быть возвращенным как количество вращений относительно FR.

Выходные аргументы

свернуть все

Диапазоны частот отказа, возвращенные как Nx2 массив, где N количество частот отказа. FB возвращен в тех же модулях как FR, или в Герц или в порядках. Используйте сгенерированные диапазоны частот отказа, чтобы извлечь спектральные метрики с помощью faultBandMetrics. Сгенерированные полосы отказа, [F−W2, F+W2], сосредоточены на характеристических дефектных частотах и их гармониках и боковых полосах для:

  • Введите частоту дефекта вала, Fi

  • Выведите частоту дефекта вала, Fo

  • Частота дефекта mesh механизма, Fm

  • Фаза Assembly передает дефектную частоту, Fa

Когда вы задаете боковые полосы, gearMeshFaultBands вычисляет боковые полосы относительно частот дефекта вала ввода и вывода:

  • Диапазоны частот отказа для входного механизма дезертируют с его гармониками и первой боковой полосой в Fi

  • Диапазоны частот отказа для выходного механизма дезертируют с его гармониками и первой боковой полосой в Fo

gearMeshFaultBands обрезает отрицательные диапазоны частот отказа автоматически и генерирует предупреждающее сообщение.

Значение W ширина диапазонов частот, которые можно задать использование 'Width'пара "имя-значение".

Информация о диапазонах частот отказа в FB, возвращенный как структура со следующими полями:

  • Centers — Центральные частоты отказа

  • Labels — Метки, описывающие каждую частоту

  • FaultGroups — Числа группы отказа равняются количеству частот

Алгоритмы

gearMeshFaultBands вычисляет различные характеристические частоты отказа можно следующим образом:

  • Введите частоту дефекта вала, Fi=FR

  • Выведите частоту дефекта вала, Fo=NiNoFR

  • Частота дефекта mesh механизма, Fm=NiFR=NoFo

  • Фаза Assembly передает дефектную частоту,Fa=Fmgcd(Ni,No)

Ссылки

[1] Ленг, Джордж Фокс. “Геометрия S&V 101”. Звук и Вибрация 33 (1999): 16-26.

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

| | |

Темы

Введенный в R2019b

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