НАУКА и ОБРАЗОВАНИЕ

УДК 62-52.

МГТУ им. НЭ. Баумана

popov@bmstu.ru

sosn@bmstu.ru

m_s_@mail.ru

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

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

Приведенный в работе [2] анализ виброакустических характеристик центробежного насоса показывает, что вибрации могут быть вызваны кавитацией, которая, к тому же, ухудшает напорные характеристики насоса. По этим причинам при создании насосов всегда учитывается требуемый кавитационный запас, что позволяет в качестве основного фактора, влияющего на вибрации в насосном агрегате, рассматривать гидродинамическую нагрузку на ротор центробежного насоса. При этом следует иметь в виду, что гидродинамические процессы могут быть нестационарными [3].  Например, при «потере насосом привода» происходит самоторможение ротора насоса, сопровождающееся изменением расхода жидкости, которая течет в каналах рабочего колеса. Изменение расхода жидкости может происходить и при постоянной частоте вращения ротора насоса, если регулирование насоса выполняется с помощью задвижки. В обоих случаях гидродинамическая нагрузка на ротор насоса будет изменяться. Чтобы рассчитать в указанных условиях гидродинамическую нагрузку на ротор насоса, необходимо предварительно построить физическую и математическую модели исследуемого процесса.

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

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

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

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

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

,                                                                                                     (1)

где  – сила, вызывающая вибрацию ротора центробежного насоса, Н;

;

;

где  – подача насоса, м³/с;

        – плотность перекачиваемой жидкости, кг/ м³;

        – число лопастей у рабочего колеса насоса;

        – частота вращения ротора насоса, 1/с;

      – диаметр входа в насос, м;

       – углы соответственно;

      – параметр, определяемый конструкцией насоса, =2 – насос со средним расположением рабочего колеса; =1 – консольный насос.

Значения сил, вызывающих вибрацию ротора насосного агрегата, и подачи насоса типа НМ 5000-210 приведены в таблице 1.

Таблица 1

Размерные и безразмерные силы вибрации

 – относительное значение подачи насоса;

*= – относительное значение силы, вызывающей вибрацию ротора насоса, при  эта сила равна .

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

Обращаясь к приведенным в таблице результатам расчета, следует учитывать, что все силы были вычислены в предположении работы насоса с постоянной частотой вращения ротора. Регулирование подачи насоса в таком случае осуществляется путем изменения характеристик сети (нагрузки на насос). Необходимо также заметить, что для расчетов по формуле (1) принята квазистационарная модель течения жидкости в насосе. В такой модели не учитываются возможные изменения эпюр скоростей жидкости в насосе и волновые процессы в подключенных к насосу трубопроводах [3].

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

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

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

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

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

Для численного исследования влияния изменений частоты вращения ротора насоса и подачи насоса на силу, вызывающую вибрации, необходимо знать функции, входящие в уравнение (1) в виде коэффициентов. Их определение с помощью формул, приведенных вместе с уравнением, не возможно без подробных, чертежей насоса. Поэтому, исходя из подобия насосов типа НМ с одной номинальной подачей, требуемые коэффициенты были определены по приведенным в работе [1] графикам зависимостей силы, вызывающей вибрацию в насосе, от его относительной подачи. Такие графики в названной работе даны для всего ряда насосов типа НМ, в том числе для насоса НМ 7000-210. Идентификация коэффициентов была выполнена для насоса НМ 5000-210, так как для него ранее построенная математическая модель имела экспериментальное подтверждение [4]. Затем определялись параметры для нескольких значений относительной подачи насоса НМ 7000-210 при n=const(табл. 2).

Таблица 2

Результаты идентификации параметров математической модели для вычисления сил,
вызывающих вибрацию ротора насоса НМ 7000-210 при n=3000 об/мин

Здесь ; ; ; ; ;

            Математическая модель для расчета переходных процессов в насосном агрегате представлена в виде структурной схемы [4].

Рис. 1. Структурная схема насосного агрегата

Возмущающим воздействием в данном случае  является величина начальной мощности (N_н1)₀, подводимой к насосу до того, как произошло отключение его электродвигателя. Мгновенное отключение электродвигателя представлено ступенчатым (скачком) уменьшением этой мощности от начального значения до нуля. Регулируемой величиной служит Δp_вх, а основным управляющим воздействием - изменение расхода ΔQ_кл. Незначительное  увеличение этого воздействия создает расход  ΔQ_рег  жидкости (нефти), поступающей в систему сглаживания волн давления (ССВД).

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

,

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

            Переходные процессы при отключении электродвигателя насосного агрегата НМ 7000-210 представлены на рисунке 2.  Расчеты выполнены с учетом влияния частоты вращения  n ротора насоса на значение коэффициента B1. 1.94 м³/c (7000 м³/ч), = 0,88; = кВт;  k5=3.75. На том же рисунке приведен график зависимости  , который получен с помощью формулы (1) при постоянном значении частоты вращения  n ротора насоса НМ 7000-210 (см. таблицу 2).

а)                                                                               б)

в)                                                                               г)

Рис. 2. Переходные процессы и графики зависимости . Насос НМ7000-210:
а, б, г - переходные процессы при внезапном отключении электродвигателя насоса;
в1 - график зависимости , значения  соответствуют переходному процессу «г»;
в2 - график зависимости  - соответствуют таблице 2 при  n=const.

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

1. Рассматриваемые силы значительно увеличиваются в случае уменьшении нагруженности насоса по производительности (подаче) при постоянной частоте вращения ротора насоса. Для  насоса типа  НМ 7000– 210 указанное увеличение может быть четырехкратным в том случае, если относительная подача насоса снижается до 0,2  своего номинального значения. Вычисленные  изменения  нагруженности  насосов  могут иметь место, когда регулирование их подачи  осуществляются с помощью задвижек.
2. При переходных процессах, вызванных  отключением электродвигателя насоса, в том числе и при нештатных ситуациях с «потерей» насосом привода, происходят изменения во времени частоты вращения ротора насоса, а также его подачи.
   Компьютерное моделирование показывает, что во время рассмотренных переходных процессов сила, способная вызвать вибрацию ротора насоса,  уменьшается  по  сравнению  с  силой, действующей на ротор насоса в начале процесса.
3. При выборе конструкции насоса необходимо обратить внимание на то, что силы, вызывающие в насосе вибрации, неопределенно ориентированы по отношению к ротору насоса. Вследствие этого в значительной мере усложняется реализация эффективной системы демпфирования ротора насоса. Применение в  конструкции насоса вместо подшипников качения подшипников жидкостного трения  позволяет решить проблему. Расчет подшипников необходимо выполнять по наибольшему значению  (см. график 2 на рисунке 2 в).

СПИСОК ИСПОЛЬЗОВАННЫХ  ИСТОЧНИКОВ

1.     Перевощиков С.И. Разработка научных основ управления вибрацией гидродинамического происхождения в центробежных насосах магистральных трубопроводов: Дис.  …д-ра техн. наук, Тюмень. 2004. 347 с.

2.     Покровский Б.В, Рубинов В.Я. О виброакустическом проявлении кавитации в центробежном насосе// Труды ВНИИ ГИДРОМАШ. Исследование, расчет и технология изготовления гидромашин. М.: Энергия, 1977.

3.     Попов Д.Н. Нестационарные гидромеханические процессы. М.: Машиностроение, 1982. 240 с.

4.     Попов Д.Н., Сосновский Н.Г. Структурный метод моделирования на ЭВМ нестационарных процессов в системах с лопастными насосами// Докл. на конференции 4-го Международного форума PCVEXPO’2005 «Насосы. Компрессоры. Арматура».