Как рассчитать изменение вибрации при изменении массы ротора: методы, формулы и практические нюансы

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

Физическая суть: от массы к вибрации

Вибрация ротора на рабочей частоте вызвана центробежной силой небаланса. Для дискретной массы m, расположенной на радиусе r от оси вращения, амплитуда центробежной силы:

F = m · r · ω²

где ω = 2πn/60 — угловая скорость (рад/с), n — частота вращения (об/мин).

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

На практике используют понятие коэффициента влияния (influence coefficient) — комплексную величину, связывающую вектор небаланса с вектором вибрации:

V = α · U

где V — вектор вибрации (амплитуда и фаза), U = m·r — вектор небаланса (масса × плечо), α — коэффициент влияния (комплексный, учитывает амплитудное отношение и фазовый сдвиг).

Основные сценарии изменения массы

В инженерной практике встречаются три типичные ситуации:

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

Во всех случаях задача сводится к оценке вектора изменения небаланса ΔU и прогнозированию ΔV = α · ΔU.

Определение коэффициента влияния

Коэффициент влияния α можно получить тремя способами:

1. Пробный запуск (метод пробных масс) — устанавливают известную пробную массу mt на известном радиусе rt, измеряют изменение вектора вибрации ΔVt. Тогда α = ΔVt / (mt·rt). Это самый надёжный способ для конкретного ротора на конкретной частоте.
2. Расчёт по модели (FEM/модальный анализ) — используется на стадии проектирования или для сложных роторов. Требует точной модели жёсткостей, масс, демпфирования и граничных условий.
3. Справочные данные и эмпирические формулы — для типовых роторов (турбины, компрессоры, электродвигатели) существуют ориентировочные значения чувствительности к небалансу, но их разброс велик (±50% и более). Пригодны только для предварительной оценки.

Важно: коэффициент влияния сильно зависит от скорости вращения. Вблизи критических скоростей (резонансов) он возрастает на порядки, фаза меняется на ~180°. Все расчёты должны выполняться для конкретной рабочей скорости.

Расчёт исправительной массы (одиночная плоскость)

Для роторов, допускающих балансировку в одной плоскости (короткие роторы, диски, вентиляторы), вектор исправительной массы:

Ucorr = -Vmeas / α

где Vmeas — измеренный вектор вибрации до коррекции. Модуль даёт произведение массы на плечо m·r, аргумент — фазу установки (противоположную фазе небаланса).

Если исправительная масса устанавливается на радиусе rcorr, то требуемая масса:

mcorr = |Ucorr| / rcorr

Двухплоскостная балансировка и матрица влияния

Для длинных роторов (турбины, генераторы, компрессоры) небаланс распределён по длине. Коррекция производится в двух плоскостях. Связь записывается матрицей коэффициентов влияния 2×2:

[V₁] [α₁₁ α₁₂] [U₁]

[V₂] = [α₂₁ α₂₂] [U₂]

где V₁, V₂ — вибрация на двух измерительных плоскостях (обычно подшипники), U₁, U₂ — небалансы в корректирующих плоскостях.

Для расчёта исправительных масс решают систему:

[U₁] [α₁₁ α₁₂]⁻¹ [V₁]

[U₂] = [α₂₁ α₂₂]⁻¹ [V₂]

Матрица определяется двумя пробными запусками (поочерёдно в каждой плоскости) или одним запуском с двумя пробными массами одновременно (метод moindre carrés).

Прогноз вибрации при известном изменении массы

Если масса изменяется не как исправительная, а как данность (замена детали, набег отложений), прогноз делается прямо:

ΔV = α · ΔU

где ΔU = Δm · r · ejφ — вектор изменения небаланса. Зная α (от пробного запуска или модели), получают ожидаемый вектор вибрации после изменения.

Пример: при замене импеллера масса изменилась на +0.5 кг на радиусе 0.2 м под углом 45°. ΔU = 0.1∠45° кг·м. Если α = 0.8∠-30° мм/с на кг·м, то ΔV = 0.08∠15° мм/с. Вектор прибавляется к текущему вектору вибрации.

Учёт допусков: классы балансировки по ISO 21940-11

Международный стандарт ISO 21940-11 (ранее ISO 1940-1) задаёт допустимые остаточные небалансы через класс балансировки G (мм/с). Допустимый остаточный небаланс на плоскость:

Uper = (G · M) / (ω · 9549)

где M — масса ротора (кг), ω — угловая скорость (рад/с), 9549 — коэффициент перевода единиц. Результат в г·мм.

Типичные классы: G6.3 — общее машиностроение, G2.5 — турбины, компрессоры, шлифовальные шпиндели, G1.0 — гироскопы, высокоточные шпиндели, G0.4 — сверхточные приборы.

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

Практические нюансы и источники ошибок

• Нелинейность и зазоры — при больших небалансах или слабых опорах (зазоры в подшипниках, фрикционное демпфирование) коэффициент влияния зависит от амплитуды. Пробная масса должна давать измеримое, но не слишком большое изменение вибрации (обычно 10–30% от исходной).
• Термическая чувствительность — ротор нагревается, меняются зазоры, жёсткости, асимметрия. Коэффициенты влияния, измеренные на холодном роторе, могут не работать в горячем состоянии. Для критических машин проводят «горячую» балансировку или используют термокомпенсацию.
• Скольжение фазы — фазовый маркер (отражательная метка, датчик ключа) должен давать стабильный отсчёт. Ошибка фазы 10° даёт ошибку вектора небаланса ~17%.
• Взаимное влияние плоскостей — при двухплоскостной балансировке недиагональные члены матрицы (α₁₂, α₂₁) часто сравнимы с диагональными. Игнорирование их (метод «по одной плоскости») приводит к недоуравновешиванию.
• Изменение собственных частот — добавление массы (особенно большой, например, втулки) меняет массу ротора и его собственные частоты. Если рабочая скорость близка к критической, даже малая добавленная масса может сдвинуть резонанс и резко изменить вибрацию.

Пошаговый алгоритм расчёта для инженерной задачи

1. Соберите исходные данные: масса ротора M, рабочая скорость n, текущий вектор вибрации V0 (амплитуда и фаза на измерительных точках), радиусы корректирующих плоскостей r₁, r₂, требуемый класс балансировки G.
2. Определите коэффициенты влияния: проведите пробный запуск(и) с известными пробными массами. Зафиксируйте ΔV для каждого пробы. Вычислите α (или матрицу α). Проверьте адекватность: повторный пробный запуск должен давать тот же α в пределах 10–15%.
3. Рассчитайте требуемый вектор небаланса: Ureq = -V0 / α (одна плоскость) или решение матричной системы (две плоскости).
4. Переведите в массы: m₁ = |U₁|/r₁, m₂ = |U₂|/r₂. Фазы установки = arg(U₁), arg(U₂).
5. Проверьте по допуску: вычислите остаточный небаланс после установки рассчитанных масс (учитывая погрешность установки ±5–10° по фазе, ±5% по массе). Убедитесь, что он ниже Uper для класса G.
6. Выполните коррекцию и контрольный замер: после установки масс измерьте вибрацию. Если остаточная вибрация выше расчётной — анализируйте причины (нелинейность, ошибка фазы, изменение условий, второй гармоник).
7. Документируйте: зафиксируйте итоговые векторы вибрации, установленные массы, фазы, коэффициенты влияния. Это база для будущих балансировок того же ротора.

Типичные ошибки и как их избежать

Ошибка	Последствие	Правильный подход
Использование коэффициентов влияния с другой скорости или другого ротора	Грубая ошибка расчёта массы (кратная), ухудшение вибрации	Определять α на данном роторе на рабочей частоте пробным запуском
Игнорирование фазового сдвига (работа только с амплитудами)	Невозможно определить угол установки массы; метод «проб и ошибок»	Всегда измерять и использовать фазу вибрации относительно маркера
Пробная масса слишком мала (изменение вибрации < 5–10%)	Низкое отношение сигнал/шум, большая погрешность α	Подбирать пробную массу для изменения вибрации 15–30%
Пробная масса слишком велика (близко к нелинейности)	α зависит от амплитуды, расчёт неверен	Не превышать амплитуду вибрации 50–70% от допуска или пределов линейности
Балансировка в одной плоскости для длинного ротора	Остаточный динамический небаланс, вибрация на второй опоре	Оценка отношения длины к диаметру: L/D > 0.5–1 требует двухплоскостной балансировки
Неучёт изменения собственных частот при большой добавленной массе	Непредвиденный резонанс, рост вибрации в разы	Модальный анализ или расчёт критических скоростей после массогабаритных изменений

Сценарии: как действовать в типовых ситуациях

Сценарий 1: Плановая балансировка ротора в сборе

Есть доступ к корректирующим плоскостям, есть измерительная система. Действуете по алгоритму выше. Проводите пробный запуск, считаете массы, устанавливаете, замеряете. Обычно 1–2 итерации достаточно для попадания в допуск G6.3–G2.5.

Сценарий 2: Замена импеллера/колеса на турбине/компрессоре

Новый элемент имеет иную массу и небаланс. Если известен вектор остаточного небаланса нового элемента (по протоколу заводской балансировки) и коэффициенты влияния ротора — прогнозируете вибрацию: Vnew = Vold + α·ΔUelement. Если прогноз выше допуска — требуется балансировка в сборе после установки.

Сценарий 3: Ремонт с изменением массы ротора (наварка, наплавка, замена втулки)

Масса ротора меняется существенно. Пересчитываете критические скорости (МКЭ или аналитически). Если рабочая скорость уходит в зону резонанса — меняют жёсткость опор или массу. После механических работ проводите полную балансировку с новыми пробными запусками — старые α недействительны.

Сценарий 4: Накопление отложений / износ в процессе эксплуатации

Небаланс нарастает постепенно. Тренд вибрации позволяет оценить скорость нарастания небаланса. Плановая балансировка проводится по достижении пороговых значений (по ISO 10816 / ISO 20816 или корпоративным нормам). Важно: если нарастание неравномерное (сброс отложений, трещина) — вибрация меняется скачкообразно, требуется диагностика, а не просто балансировка.

Инструменты и программное обеспечение

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

• Встроенное ПО портативных балансировщиков (Bruel & Kjaer, Schenck, CSI/Emerson, Vibro-Meter, Baltech, «Вибромера», «Этрон» и др.) — автоматизируют пробные запуски, решают матрицу, выдают массы и фазы.
• Системы онлайн-мониторинга (Bently Nevada, Metrix, Pruftechnik, «Вибротех», «СПМ») — хранят историю векторов, позволяют отслеживать тренды и сравнивать с базой.
• CAE-пакеты (ANSYS, NX Nastran, SAMCEF, DyRoBeS, XLRotor) — для расчёта коэффициентов влияния на стадии проектирования и анализа изменений конструкции.

Выбор инструмента зависит от частоты задач, сложности роторов и требований к отчётности.

Контроль качества результата

После балансировки проверяют не только амплитуду на рабочей частоте (1×), но и:

• Гармоники 2×, 3× — их рост может указывать на невыравнивание, неравномерную жёсткость, трещины, слабость опор, а не на небаланс.
• Фазовую стабильность — при нескольких замерах подряд фаза не должна «гулять» более 5–10°.
• Вибрацию в вертикальной/горизонтальной/осевой направлениях — небаланс даёт преимущественно радиальную вибрацию; осевая компонента говорит о других дефектах.
• Повторяемость при перезапуске — идеально проверить вибрацию после остановки и повторного разгона (термические перепады, седление).

Резюме: от расчёта к решению

Расчёт изменения вибрации при изменении массы ротора строится на линейной модели ΔV = α · ΔU. Ключевые условия применимости: линейная область работы, постоянные параметры ротора и опор, известная рабочая скорость, достоверные коэффициенты влияния.

Практический алгоритм:

1. Получите достоверные α пробным запуском на рабочей частоте.
2. Определите вектор изменения небаланса ΔU (известная масса × плечо × фаза).
3. Рассчитайте прогноз ΔV = α · ΔU.
4. Для балансировки найдите Ucorr = -V/α (или решите матричную систему).
5. Проверьте результат по допуску класса G (ISO 21940-11).
6. Всегда делайте контрольный замер и анализируйте спектр, а не только амплитуду 1×.

Если масса меняется 크게 (ремонт, замена крупных элементов) — пересчитывайте критические скорости и определяйте α заново. Не переносите коэффициенты с других машин, других скоростей или прошлых лет без верификации. Качественная балансировка — это не подбор массы «на глаз», а измерение, расчёт, установка и подтверждение результата.

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