Как выбрать частоту дискретизации для сбора вибрационных данных

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

Теорема Найквиста задаёт только теоретический нижний предел: для сигнала с максимальной частотой Fmax частота дискретизации должна быть выше 2 × Fmax. Для реальной вибродиагностики этого условия самого по себе недостаточно. Между полезной полосой и половиной частоты дискретизации обычно требуется запас, позволяющий антиалиасинговому фильтру подавить более высокочастотные составляющие до преобразования сигнала в цифровую форму.

Сначала определяют не частоту дискретизации, а полезную полосу

Главный исходный параметр — верхняя частота, до которой в сигнале находится диагностически значимая информация. Она зависит от конструкции машины и от того, какую неисправность требуется обнаружить.

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

Поэтому при подготовке измерения полезно заранее определить:

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

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

Почему правила «дискретизировать в два раза быстрее» недостаточно

Для идеально ограниченного по полосе сигнала теорема отсчётов требует, чтобы частота дискретизации превышала удвоенную максимальную частоту сигнала. Частота fs/2 называется частотой Найквиста.

Например, если интерес представляет диапазон до 1 кГц, математический минимум располагается выше 2 кГц. Однако настройка системы непосредственно возле этого предела обычно оставляет слишком мало места для реального фильтра.

У физического антиалиасингового фильтра нет вертикальной границы между пропускаемой и подавляемой полосами. Его коэффициент передачи уменьшается постепенно. Нужна переходная область, в которой высокочастотные компоненты успевают достаточно ослабнуть до частоты, после которой они могут отразиться внутрь анализируемого диапазона.

Поэтому практический порядок другой:

1. Определить максимальную полезную частоту Fmax.
2. Выяснить характеристики антиалиасингового фильтра измерительной системы.
3. Выбрать режим дискретизации, при котором вся полезная полоса находится внутри гарантированной рабочей полосы прибора.
4. Проверить полосу датчика, усилителя и остальных элементов тракта.
5. После этого рассчитать длительность записи и параметры FFT.

Нельзя назначить один универсальный коэффициент между Fmax и частотой дискретизации для всех приборов. В одной системе полезная полоса может занимать значительную часть интервала до fs/2, в другой производитель устанавливает более существенный запас из-за особенностей фильтрации и преобразования.

Что происходит при алиасинге

Алиасинг возникает, когда в аналого-цифровой преобразователь поступают частоты, которые система с выбранной частотой дискретизации не способна однозначно представить. Высокочастотные компоненты при этом могут появиться в цифровом спектре как ложные более низкие частоты.

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

Отсюда следует принципиальное требование: защита от алиасинга должна действовать до того, как нежелательная составляющая необратимо попадёт в оцифрованный диапазон.

Обычный цифровой фильтр, применённый к уже полученной записи, не может восстановить первоначальную частоту компонента, который успел наложиться на другую область спектра. В современных измерительных системах антиалиасинговая обработка может быть частью сложного аналого-цифрового тракта с передискретизацией и внутренней цифровой фильтрацией, поэтому следует ориентироваться на заявленную производителем полезную полосу конкретного режима измерения.

Частота дискретизации и частотное разрешение FFT — не одно и то же

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

Для блока из N отсчётов при частоте дискретизации fs номинальный шаг между линиями FFT составляет:

Δf = fs / N

Поскольку длительность блока T = N / fs, ту же зависимость можно записать как:

Δf = 1 / T

Следовательно, если необходимо лучше разделять близкие частоты, обычно нужно увеличивать длительность анализируемого фрагмента, а не бесконечно повышать fs.

Условный пример: два измерения могут иметь одинаковую длительность, но разные частоты дискретизации. Более высокая fs даст больше отсчётов и расширит доступную высокочастотную область, однако номинальный шаг FFT при одинаковой длительности записи останется тем же.

На реальную способность различать составляющие дополнительно влияет выбранное спектральное окно. Оконная функция уменьшает утечку спектра, но одновременно изменяет форму спектральных линий, поэтому шаг FFT нельзя считать единственным показателем фактического разрешения.

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

При анализе медленно вращающихся машин часто возникает противоположная проблема: частота дискретизации достаточно велика, но запись слишком короткая.

Если интересующая составляющая имеет очень низкую частоту, за короткий интервал наблюдения система увидит лишь небольшую часть её циклов. В результате ухудшается устойчивость оценки амплитуды и частоты, а близкие спектральные компоненты оказываются трудноразличимыми.

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

Как учитывать частотный диапазон акселерометра

Даже правильно выбранная частота дискретизации бесполезна, если необходимая полоса находится за пределами достоверной частотной характеристики датчика.

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

При выборе тракта следует согласовать:

• нижнюю и верхнюю рабочие частоты датчика;
• требуемый диапазон ускорений;
• частотный диапазон усилителя или IEPE-кондиционера;
• частоту среза встроенных и внешних фильтров;
• полезную полосу АЦП;
• частоту дискретизации;
• способ крепления акселерометра.

Последний пункт особенно важен для высокочастотных измерений. Нежёсткое крепление датчика способно ухудшить передачу высокочастотной вибрации и создать дополнительные резонансные эффекты. Увеличение частоты дискретизации такой недостаток механической установки не исправляет.

Как задача диагностики меняет требования к частоте

Задача	Что определяет верхнюю полосу	На что обратить внимание
Общий уровень вибрации	Полоса, в которой должен рассчитываться контролируемый параметр	Фильтрация и единообразие полосы при сравнении измерений
Обычный FFT-анализ машины	Наиболее высокая диагностически значимая гармоника или механическая составляющая	Антиалиасинг, длительность записи, спектральное окно
Редукторы	Зубцовые частоты, их гармоники и боковые полосы	Частоты вращения всех валов и количество зубьев
Подшипники качения	Выбранный метод диагностики и высокочастотная область отклика	Не путать полосу обычного FFT с полосой, используемой для анализа огибающей
Удары и переходные процессы	Скорость изменения временного сигнала и высокочастотный состав импульса	Полоса датчика, перегрузка, запись формы импульса
Машины с переменной скоростью	Максимальная анализируемая кратность частоты вращения при максимальных оборотах	Тахометрический канал, синхронизация и порядок пересчёта в угловую область

Особенности анализа подшипников и огибающей

При диагностике подшипников выбор частоты дискретизации требует особого внимания. Характерные частоты дефекта, вычисленные по геометрии подшипника и скорости вращения, сами по себе могут быть сравнительно низкими. Однако локальный дефект создаёт короткие воздействия, способные возбуждать значительно более высокочастотные резонансы конструкции и датчика.

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

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

Что меняется при ударах и коротких переходных процессах

Для стационарной гармонической вибрации можно сосредоточиться на её спектральном диапазоне. Для удара важна ещё и форма сигнала во времени.

Чем короче фронт импульса, тем шире его спектр. Система с недостаточной полосой сглаживает быстрые изменения: измеренный пик может стать ниже и шире реального, а временная форма — заметно искажённой.

При таких измерениях выбор тракта нужно начинать не только с предполагаемой основной частоты, но и с необходимой временной точности. Проверяют полосу акселерометра, собственный резонанс, диапазон амплитуд, перегрузочную способность, фильтры и скорость АЦП.

Просто назначить очень высокую частоту дискретизации тоже недостаточно: если аналоговый датчик или усилитель ограничивает полосу значительно раньше, дополнительные отсчёты не восстановят отсутствующую информацию.

Переменная скорость и порядок вращения

Для машин с изменяющейся скоростью обычный FFT во временной области может давать размазанные спектральные линии. Причина заключается не обязательно в неправильной частоте дискретизации: частота вращения меняется непосредственно во время записи.

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

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

Если нужен фазовый анализ, орбиты валов или сравнение нескольких каналов, появляется ещё одно требование: каналы должны иметь известную временную связь. Большая частота дискретизации сама по себе не исправит неизвестную задержку между каналами последовательно опрашиваемого мультиплексированного АЦП.

Нужно ли брать частоту дискретизации с большим запасом

Запас полезен, но его следует выбирать осмысленно. Передискретизация может упростить фильтрацию, дать более подробное представление временной формы и оставить резерв на неизвестные высокочастотные компоненты. Однако после определённого уровня увеличение fs перестаёт приносить пользу конкретной диагностической задаче.

Чрезмерная частота дискретизации имеет свои последствия:

• растёт объём записываемых данных;
• увеличиваются требования к пропускной способности интерфейса и накопителя;
• обработка FFT и архивов требует больше ресурсов;
• при многоканальном мониторинге быстрее заполняется хранилище;
• может записываться высокочастотный шум, который не несёт диагностической информации;
• при неизменном количестве точек FFT увеличение fs ухудшает шаг по частоте.

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

Можно ли сначала записать с высокой частотой, а потом уменьшить её

Да, такой подход полезен, если исходный широкополосный сигнал действительно нужен для нескольких видов анализа. Можно получить одну высокочастотную запись, а затем сформировать из неё потоки с меньшей частотой дискретизации.

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

Многоскоростной подход особенно удобен в системах мониторинга: широкополосные временные записи можно сохранять по событию, тогда как для долговременных трендов использовать существенно меньший объём обработанных параметров.

Как подобрать параметры для нового объекта

Если частотный состав оборудования заранее неизвестен, настройку лучше проводить последовательно.

1. Опишите диагностическую цель. Укажите, нужно ли контролировать дисбаланс, несоосность, зубчатую передачу, подшипники, ударные процессы или несколько типов дефектов.
2. Составьте карту механических частот. Зафиксируйте скорости валов, передаточные отношения, число зубьев и лопаток, а также другие источники периодического возбуждения.
3. Определите верхнюю полезную частоту. Учтите не только основные компоненты, но и необходимые гармоники или высокочастотный диапазон выбранного метода диагностики.
4. Проверьте датчик. Необходимая область должна находиться в его пригодной для измерений частотной полосе, а ожидаемые амплитуды — внутри измерительного диапазона.
5. Проверьте весь измерительный тракт. Ограничивающим элементом может оказаться не акселерометр, а кондиционер, фильтр или входной модуль.
6. Изучите антиалиасинговую характеристику системы. Ориентируйтесь на гарантированную полезную полосу выбранного режима, а не только на значение fs в настройках.
7. Выберите частоту дискретизации. Она должна обеспечивать требуемую полосу с необходимой переходной областью для фильтра.
8. Назначьте длительность записи. Рассчитайте её исходя из требуемого шага FFT и самых низких анализируемых частот.
9. Выполните пробное широкополосное измерение. Проверьте, есть ли значимые компоненты вблизи верхней границы выбранной полосы.
10. Оптимизируйте режим. Если верхняя область содержит только нерелевантный шум, полосу и объём данных можно уменьшить без потери необходимой информации.

Как проверить, что выбранной частоты хватает

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

После этого сравнивают:

• амплитуды основных спектральных компонентов;
• наличие пиков возле верхней границы диапазона;
• форму временного сигнала;
• результаты расчёта диагностических показателей;
• поведение спектра после увеличения рабочей полосы.

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

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

Типичные ошибки при настройке сбора вибрации

Ориентироваться только на обороты двигателя

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

Выбирать fs ровно вдвое выше интересующей частоты

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

Считать, что высокая fs автоматически даёт хороший FFT

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

Игнорировать характеристики датчика

АЦП способен записывать частоты, на которых установленный акселерометр уже не обеспечивает нужной точности. Верхняя граница системы определяется самым слабым звеном измерительного тракта.

Фильтровать алиасинг после записи

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

Использовать один режим для всех диагностических задач

Измерение общего уровня вибрации, анализ зубчатой передачи, поиск дефектов подшипников и запись ударов предъявляют разные требования к полосе. Универсальный режим возможен, но нередко приводит либо к недостатку информации, либо к неоправданно большому объёму данных.

Частота дискретизации — только одна часть измерительной системы

Достоверность вибрационных данных определяется всей цепочкой от механической поверхности до файла с отсчётами. Даже идеально рассчитанная fs не компенсирует плохое крепление акселерометра, насыщение входа, недостаточный динамический диапазон, неподходящий датчик или нарушение синхронизации каналов.

Перед постоянным мониторингом полезно проверить весь путь сигнала: объект — крепление — датчик — кабель — питание или кондиционирование — аналоговая фильтрация — АЦП — цифровая обработка — сохранение данных.

Особое внимание требуется к перегрузке. Если амплитуда превышает диапазон датчика либо входного каскада, возникает клиппирование или другое нелинейное искажение. Оно само создаёт дополнительные спектральные составляющие. Повышение частоты дискретизации в такой ситуации не делает запись достоверной.

Практический принцип выбора

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

Не стоит использовать коэффициент 2, 2,5, 5 или 10 как универсальную норму без привязки к аппаратуре. Удвоенная частота — теоретическая граница для ограниченного по полосе сигнала, тогда как реальная допустимая полоса определяется устройством АЦП и его фильтрами. Большой запас может быть полезен для временных форм и неизвестного спектра, но также увеличивает поток данных.

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