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

Время измерения вибрации нельзя выбирать только по принципу «чем дольше, тем устойчивее». Сначала определяют требуемое частотное разрешение спектра, которое задаёт минимальную длину одного временного блока, а затем рассчитывают, сколько таких блоков нужно обработать для усреднения. Для обычного FFT-анализа основная связь проста: чем меньше требуемый шаг по частоте, тем длиннее должен быть анализируемый отрезок сигнала.

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

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

Полезно разделять два понятия: время одного FFT-блока и полную продолжительность измерения. Они отвечают на разные вопросы.

• Длина FFT-блока определяет частотное разрешение одного спектра.
• Количество усреднений влияет на стабильность оценки при наличии случайной составляющей.
• Перекрытие блоков определяет, сколько новых данных нужно получить для следующего спектра.
• Время установления режима приходится учитывать отдельно, если после пуска или изменения нагрузки оборудование не сразу выходит на стационарный режим.
• Характер сигнала определяет, имеет ли смысл усреднение вообще: для стационарной вибрации и кратковременного удара подходы различаются.

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

Первый шаг: определить необходимое частотное разрешение

Для стандартного FFT шаг между соседними частотными линиями связан с длительностью анализируемого блока:

Δf = 1 / Tбл

где Δf — шаг спектральных линий в герцах, а Tбл — продолжительность одного временного блока в секундах.

Отсюда длина блока рассчитывается обратным выражением:

Tбл = 1 / Δf

Если требуется спектр с шагом 1 Гц, один блок должен охватывать 1 секунду сигнала. Для шага 0,5 Гц потребуется 2 секунды, для 0,2 Гц — 5 секунд, а для 0,1 Гц — 10 секунд.

Желаемый шаг спектра	Минимальная длина FFT-блока	Типичная задача
2 Гц	0,5 с	Обзор спектра и поиск выраженных составляющих
1 Гц	1 с	Общий анализ стационарной вибрации
0,5 Гц	2 с	Разделение относительно близких частот
0,2 Гц	5 с	Более детальное исследование низкочастотной области
0,1 Гц	10 с	Разделение очень близких спектральных составляющих

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

Шаг спектральных линий и способность разделить пики — не одно и то же

Частая ошибка — считать, что два источника на частотах 100 и 101 Гц обязательно будут хорошо различимы при Δf = 1 Гц. Реальная способность разделить соседние составляющие зависит также от применённого временного окна, отношения амплитуд двух сигналов, шума и стабильности их частот.

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

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

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

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

Tбл = N / Fs

а шаг частотной сетки:

Δf = Fs / N

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

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

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

Второй шаг: определить необходимое число усреднений

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

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

Усреднение особенно полезно, когда присутствуют:

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

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

Как понять, что усреднение уже достаточно

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

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

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

Расчёт полного времени с учётом усреднения

Без перекрытия временных блоков приблизительное время, необходимое для получения M отдельных спектров, равно:

Tизм = M × Tбл

Если один FFT-блок длится 2 секунды и требуется обработать 10 последовательных неперекрывающихся блоков, для самих данных понадобится около 20 секунд стационарной записи.

При перекрытии часть предыдущего блока используется повторно. Если доля перекрытия обозначена P, где 50 % соответствует P = 0,5, то для последовательности из M спектров продолжительность исходной записи можно оценить по выражению:

Tизм = Tбл × [1 + (M − 1) × (1 − P)]

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

Условный пример расчёта

Предположим, требуется Δf = 0,5 Гц. Тогда длительность одного FFT-блока составляет 2 секунды. Планируется получить 10 спектров с перекрытием 50 %.

1. Рассчитываем длину одного блока: Tбл = 1 / 0,5 = 2 с.
2. Доля новых данных для каждого следующего блока равна 1 − 0,5 = 0,5.
3. После первого блока необходимо получить ещё девять половинных блоков.
4. Продолжительность записи составит 2 × [1 + 9 × 0,5] = 11 секунд.

Без перекрытия для тех же десяти FFT понадобилось бы 20 секунд данных.

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

Зачем нужно перекрытие временных блоков

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

В результате перекрытие полезно для усреднённого анализа и увеличивает частоту обновления спектров. Но оно не улучшает базовое частотное разрешение: величину Δf по-прежнему определяет длина отдельного FFT-блока.

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

Как выбор окна влияет на необходимое время

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

Оконные функции уменьшают этот эффект, но обладают разной шириной главного лепестка и разной точностью определения амплитуды. Поэтому математического условия Δf = 1 / Tбл недостаточно для выбора реального разрешения.

В практическом анализе обычно учитывают следующие особенности:

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

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

Низкие частоты требуют более длинного наблюдения

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

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

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

Почему слишком длинный FFT-блок тоже может быть ошибкой

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

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

Получается противоречие:

• длинный блок улучшает частотное разрешение;
• короткий блок лучше локализует изменение процесса во времени.

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

Что делать при изменяющейся частоте вращения

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

Признаки проблемы можно заметить по самому спектру:

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

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

Усреднение не исправляет нестационарный режим

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

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

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

Если изменение режима является предметом исследования, его не нужно искусственно «успокаивать» усреднением. В таком случае требуется метод представления спектра во времени или относительно другого изменяющегося параметра.

Как выбрать время измерения на практике

Рабочую методику удобно строить от диагностической задачи, а не от доступного размера FFT.

1. Задайте диапазон интересующих частот. Определите, какие механические процессы должны быть видны и насколько близко могут находиться диагностически значимые компоненты.
2. Выберите требуемый шаг Δf. Он должен быть достаточно мал для разделения нужных составляющих с учётом используемой оконной функции.
3. Рассчитайте длину блока. Используйте Tбл = 1 / Δf или соответствующий параметр анализатора.
4. Проверьте стационарность на этом интервале. Убедитесь, что обороты и нагрузка не меняются настолько, чтобы пик заметно перемещался внутри блока.
5. Выберите тип окна. Он должен соответствовать характеру временного сигнала и цели анализа.
6. Добавьте усреднение. Для случайной составляющей накапливайте последовательные спектры и наблюдайте за сходимостью ключевых показателей.
7. Назначьте перекрытие при необходимости. Оно позволяет эффективнее использовать запись и повышает частоту формирования спектров, но не заменяет реальное время наблюдения.
8. Проверьте повторяемость. Повторите измерение при том же рабочем режиме и сравните диагностически значимые параметры.

Пример выбора параметров от диагностической задачи

Предположим, нужно различать две составляющие вибрации, расположенные примерно в нескольких герцах друг от друга. Начинать с очень мелкого шага вроде 0,01 Гц без необходимости невыгодно: один FFT-блок занял бы 100 секунд, а за это время рабочий режим машины может измениться.

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

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

Как оценивать устойчивость спектра без субъективной оценки графика

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

Например, можно наблюдать:

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

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

Распространённые ошибки при расчёте времени

Выбирать длительность записи только по верхней частоте

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

Принимать число FFT-линий за реальное физическое разрешение

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

Увеличивать число усреднений при плавающих оборотах

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

Считать перекрывающиеся блоки полностью независимыми

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

Использовать одно время измерения для всех режимов машины

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

Когда обычного усреднённого FFT недостаточно

Расчёт по формуле Tбл = 1 / Δf наиболее понятен для относительно стационарного процесса. Другой метод анализа стоит рассмотреть, если диагностическая информация существенно меняется во времени.

Характер измерения	Что ограничивает время	Подход к анализу
Стабильные обороты и нагрузка	Требуемое Δf и статистическая устойчивость	FFT с подходящим усреднением
Медленно меняющиеся обороты	Смещение частот внутри FFT-блока	Короткие последовательные спектры или порядковый анализ
Разгон или выбег	Сильная нестационарность	Порядковый или время-частотный анализ
Одиночный переходный процесс	Необходимость сохранить событие во времени	Триггерная запись и анализ переходного сигнала
Случайная стационарная вибрация	Разброс спектральной оценки	FFT с усреднением нескольких временных блоков

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

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

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

Для воспроизводимых измерений полезно фиксировать вместе с результатом как минимум диапазон анализа, частоту дискретизации, длину FFT-блока или Δf, оконную функцию, тип и количество усреднений, перекрытие и рабочий режим оборудования. Тогда два спектра можно сравнивать как результаты одной методики, а не просто как похожие графики.