Центробежные насосы нефтепереработки (fb2)

- Центробежные насосы нефтепереработки 2.89 Мб, 36с.  (читать) (читать постранично) (скачать fb2) (скачать исправленную) - Константин Владимирович Ефанов

Настройки текста:



Введение

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

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

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

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

__

Посвящается Богу Творцу Троице!

Выражаю благодарность моей маме, работавшей в нефтяном машиностроении.

Монтажная привязка насосных агрегатов на местах установки

На технологических установках в цехах распространена открытая установка насосов под эстакадами, навесами и этажерками с учетом в настоящее время недействующего ОСТ 26-1141-74 «Насосы. Основные требования к установке и эксплуатации вне помещений на химических, нефтехимических и нефтеперерабатывающих производствах» (могут быть заменены на легкие кожуха с отверстиями для вентиляции); для общезаводского хозяйства распространена закрытая установка насосов в насосных станция [2]. Для защиты от снега открытая площадка ограждается по периметру щитами.

Схема наружной установки насоса [5]:



Схема продувки вспомогательного трубопровода [5]:



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

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

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

Уровень площадки выше уровня прилегающей территории. Поверхность площадки выполняют с уклоном для стоков осадков и жидкостей. Реализуются условия по отсутствию образования наледи на полу, например, обогревом. На площадке делают канавы для отвода стоков и незамерзающий дренаж (с обогревающим трубопроводом-спутником) с емкостью.

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

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

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

Для соблюдения санитарных требований при наружной установке насосных агрегатов выполняют:

– автоматическое управление работой насосов и дистанционный контроль параметров процесса,

– контроль режима и условий работы обслуживающего и ремонтного персонала.

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

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

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

– шумы от работы подшипниковых узлов

– шумы при работе вентилятора системы охлаждения,

– колебаниями магнитопровода статора и ротора под действием магнитного поля.

Для обеспечения шумности электродвигателем 80 дБА, электродвигатель по всем трем указанным параметрам должен иметь соответствующие конструктивные решения.

На нефтеперерабатывающих заводах насосы делят на перекачивающие нефтепродукты, перекачивающие кислоты и щелочи, перекачивающие воду, очищенные стоки и химически неагрессивные вещества [3].

В промышленных условиях применяют в основном центробежные насосы [4].

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

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

Основные требования к обвязке насосов насосных станций:

– наличие резервного насоса на случай поломки основного насоса,

– удобство и безопасность демонтажа и обслуживания насоса,

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

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

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

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

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

Наиболее распространенные схемы обвязки насосов:

– схема обвязки для пропариваемого насоса:



– схема обвязки для продувки газом




– схема обвязки для пропариваемых и продуваемых насосов:



– схема обвязки насосов продуваемых и промываемых:



– схема обвязки для циркуляции через запасной насос (остальная часть обвязки условно не показана):




Варианты расположения трубопроводов:

– правильная трассировка труб для насосов типа НК:



– для насоса типа НДМ:



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

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

В проектном институте технологом проектируется общая технологическая схема установки, в которую интегрируется схема от насоса.

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

Схемы разрабатываются и рассчитываются в специализированных программах для технологических расчетов, по результатам расчетов выполняется графическая часть разработки схемы PI&D.

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

Приведем пример привязки насоса в технологии



Исходными данными для выбора насоса являются [2]:

– напор во всасывающей и нагнетательной линиях трубопровода,

– расход и температура перекачиваемого потока,

– вязкость и плотности жидкости.

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

Напор насоса, рассчитанный, по приведенной выше формуле Эйлера:



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

во всасывающей



и нагнетательной линиях



Технологические схемы обвязки насосов

Для комплектации отечественного насосного оборудования рекомендуется использовать торцовые уплотнения и системы их обвязки российских фирм-производителей Анод и Герметика. Отметим, что специалистами Герметики была выпущена монография «Уплотнения и уплотнительная техника».

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

Описание работ по проектированию блоков нефтяных аппаратов на современном уровне подробно изложено в работе [1].

Проектирование объектов нефтепереработки в целом приведено в работе [2].

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

Схемы обвязки современных торцовых уплотнений разрабатывают на основании приведенных данных в переведенном на русский язык стандарте ГОСТ 32600 (оригинальный стандарт API 610).

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

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

Одинарные торцовые уплотнения не комплектуются системой обвязки как таковой. В сальниковую камеру торцового уплотнения из нагнетательного патрубка насоса по трубопроводу подается жидкость для поддержания давления по плану 21, и с обеспечением циркуляции импеллером по схемам планов из ГОСТ 32600:





Защитная жидкость может подаваться на одинарное уплотнение с внешней линии:



– план 52



– план 53



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

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

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

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

Теплообменник проектируется по нормам на аппараты до 21МПа и поставляется с полным комплектом документации по ГОСТ. Чертеж на теплообменный аппарата должен быть выполнен так же, как принято выполнять и оформлять чертежи на теплообменные аппараты нефтепереработки с указанием главного вида, всех видов, разрезов, сечений, таблицы технической характеристики, таблицы штуцеров. К теплообменнику может быть приложен теплогидравличекий расчет, выполненный методом конечных элементов в специальном программном пакете или другим способом (вручную).

КИПиА измеряют отклонение параметров циркулирующей жидкости в непрерывном режиме, при необходимости выдается сигнализация и задержка отключения для перехода на резервный агрегат. Интерфейс – RS-485.

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

__

Выполняется 3D-компоновка трубной обвязки насоса в соответствии с материалом работ [22], [23]. Затем должна выполняется выборка трубопроводных линий, как это принято в строительстве технологических установок. На каждую линию выпускается отдельный чертеж с припусками по трем осям координат для обеспечения сборки. В строительстве принято выпускать изометрические проекции трубопроводов по ВСН. В машиностроении следует выпускать документацию по требованиям ЕСКД в виде ортогональных проекций с указанием размеров и конфигурации. В строительстве изометрические схемы получают автоматически из 3D-программ таких как PDMS. В нефтяном машиностроении применяют 3D-программы Компас 3D, SolidWorks, и др. в которых из 3D трубопроводная линия может быть открыта отдельно и разработан в 2D чертеж в ортогональных проекциях с глубокой конструкторской и технологической проработкой.

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

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

Система КИПиА насоса

Безопасность и корректность работы насосного агрегата обеспечивается установкой средства КИПиА.

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

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

Насосный агрегат в общем случае комплектуется средствами КИПиА:

– сигнализатора уровня в верхней части насоса для остановки агрегата при снижении уровня ниже допустимого,

– при недопустимости установки сигнализатора уровня производят визуальный контроль со смотровым фонарем перед пуском, регулятор уровня устанавливают на приемном резервуаре или емкостном аппарате,

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

– термометр устанавливают на подшипники электродвигателя,

– электродвигатели комплектуются термисторной защитой,

– может использоваться электроконтактный манометр, установленный к полости статора,

– на напорном трубопроводе может быть установлен электроконтактный маномер для отключения насоса при снижении давления Рвс+(0,8)Рном.

– отбор давления блокировки на линии отбора жидкости после фильтра и направляемой обратно в зазор между ротором и статором.

Расчет характеристик насосов

Сравнительное рассмотрение теорий методов расчета проточной части и построение характеристик методом конечных элементов приведены в работе [6]. На современном уровне развития техники проводить экспериментальное определение точек и выполнять построение характеристик необязательно. Достаточно точные результаты можно получить расчетом в программном пакете ANSYS методом конечных объемов. Методом конечных объемов решаются задачи вычислительной гидродинамики.

В гидродинамике описание движения потока жидкости производится описанием вектора скорости и двух термодинамических величин (давление и плотность) в зависимости от координат и времени. Система уравнений должна содержать 5 уравнений, в том числе уравнение неразрывности, уравнение Навье-Стокса (уравнение Эйлера для идеальной жидкости), уравнение переноса тепла (уравнение сохранения энтропии для идеальной жидкости).

– уравнение неразрывности потока:



– уравнение Навье-Стокса для сжимаемой жидкости (уравнение движения вязкой среды):



для несжимаемой жидкости при :



– уравнение переноса тепла



Закон сохранения для идеальной жидкости (при отсутствии вязкости и теплопроводности правая часть уравнения становится равной нулю и получится уравнение сохранения энтропии):



Для описания турбулентного течения потока используются четыре подхода:

– прямое численное решение уравнений Навье-Стокса,

– применение аналитических теорий турбулентности,

– применение моделей переноса турбулентности,

– применение моделей замыкания движений мелкого масштаба.

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

Конечно-разностный метод расчета сравнивается со спектральным по пяти параметрам [6]:

– скорость сходимости,

– эффективность (затраты на расчет для заданной погрешности результата),

– граничные условия (точность конечно-разностных методов нарушается около границ за счет необходимости расчёта точек вне области течения, поэтому сетка корректируется вдоль границ и усложняется),

– разрывы (сглаживание разрывов при локальных ошибках),

– априорная оценка точности (для конечно-разностных методов точность сравнивается на сетках с разным числом конечных элементов).

Аналитические теории турбулентности строятся на статическом подходе к описанию турбулентности. Динамические параметры в этих теориях являются средними характеристиками течения потока.

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

Метод прямого численного моделирования DNS – Direct Numerical Simulation предложен в работе Orszag, S. A., and Patterson, G. S. в 1972 г.

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

По методу DNS решаются уравнения Навье-Стокса напрямую непосредственно без применения моделей турбулентности (например, модели «k-ε») в отличие от других методов расчета.

При решении уравнений Навье-Стокса находят для любой точки в потоке скорость течения и давление. Результатом расчета по методу DNS является нахождение этих параметров потока.

По методу DNS возможно выполнение расчета течения для различных значений числа Re.

Существует модель однородной изотропной турбулентности, но с помощью её нельзя провести описание реального потока [6]. Существует модель локально изотропной турбулентности. Согласно этой модели турбулентные пульсации для мелких масштабов с большим числом Рейнольдса можно рассматривать как однородные изотропные. Колмогоров ввел гипотезу о том, что статический режим для мелких масштабов зависит от коэффициента вязкости k и скорости (средней) диссипации энергии ε.

Масштаб вихрей, на который влияет вязкость получается из этой гипотезы Колмогорова с учетом соображений размерности [6]:



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

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

__

Методики расчета проточной части по струйной теории Эйлера являются устаревшими и расчет необходимо выполнять методом конечных объемов в программном пакете.

Горизонтальный агрегат погружного насоса




Горизонтальные агрегаты погружных (скважинных) насосов, изготовленные в соответствии с требованиями API 610 (ГОСТ 32601), например, как указано в работе [13], могут использоваться в насосном парке нефтеперерабатывающих заводов.

ГОСТ 31840 определяет требования безопасности к погружным и полупогружным насосам. Эти требования аналогичны требованиям для полупогружных насосов, проектируемых по API 610 (ГОСТ 32601) в том числе по выбору систем торцовых уплотнений. Какие-либо препятствия для применения погружных насосов в нефтепереработке отсутствуют.

Такие агрегаты могут выглядят возможно конкурентоспособнее высоконапорных насосов типа НДМ за счет простоты конструкции и широкого диапазона напорной характеристики.

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

В нефтедобыче погружные насосы используются в составе скважинного агрегата установки электроцентробежного насоса [11], [12].

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

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

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

__

Вертикальные полупогружные насосы типа ХП.

Полупогружные насосы в нефтепереработке устанавливают на люках крышек емкостных аппаратов или резервуаров.

Функцию крышки аппарата выполняет плита насоса, но при этом для плиты необходим прочностной расчет по нормам на сосуды и аппараты для обеспечения прочности оболочки аппарата.

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

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

Компоновка по типу погружного насоса очевидно выигрывает у компоновки с корпусом от насоса НК как конструктивно, так и по технологически параметрам насоса.

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

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

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

Расчет на резонанс валов насосов

Конструкции нефтяных насосов подробно приведены и описаны в классической литературе [7], [8], [9], [10].

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

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

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

Для выполнения связанного междисциплинарного расчета рекомендуется пакет ANSYS. Для сквозного проектирования по результат расчета рекомендуется российский пакет КОМПАС 3D.

Расчет валов на резонанс по теории колебаний

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



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

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

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

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

Изменение жесткости вала связано с изменением статической деформации, которая связана со свободной частотой по формуле:



На резонансной частоте амплитуда вынужденных колебаний неограниченно возрастает при отсутствии внешних сопротивлений:



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

Для значений частот, близких к резонансной возникают биения вала. Для случая вала мешалки при отсутствии сопротивлений биению, колебания имеют вид:



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



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

Упругие колебания системы с одной степенью свободы в общем случае (вторые два члена формулы относятся к вынужденным колебаниям):



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

– собственные колебания без учета сопротивлений (f = 0, q = 0)



– собственные затухающие колебания (вынуждающая сила W = 0, )



– вынужденные колебания без учета сопротивлений (, , в формуле получается, что первый член является вынужденными колебаниями, остальные два члена свободными колебаниями)



Формула вынужденных колебаний получается из вторых двух членов уравнения упругих колебания после отбрасывания свободных колебаний и замены в формуле



Т.е. вынужденные колебания являются гармоническими (так же как и собственные)



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



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

При m<p амплитуда С приближается к статической деформации вала.

При m=p амплитуда С достигает больших величин, наступает явление резонанса вала.

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

При m>p амплитуда С стремиться к нулевому значению, колебания отсутствуют.

Приведем график амплитуд колебаний:



Как видно из рисунка, при резонансной частоте происходит разрыв кривой прогиба вала и разрушение вала.

При расчете вала необходимо не допускать наличия расчетных частот в пределах биения, то есть в пределах близких к резонансной частоте для недопущения разрушения вала. Запас может превышать критическую частоту на 20%. Такой запас, например, установлен для валов центробежных нефтяных насосов в ГОСТ 32601.

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



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

Обобщенной силой является сила, которая полностью определяет действующую систему сил.

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

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

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





Для кинетическая энергия системы



находится производная по обобщенным координате и скорости и после преобразований:



Уравнение движения запишется в виде



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



По теореме Кастильяно [17,с.319] прогиб точки приложения сосредоточенной силы (P) равен частной производной потенциальной энергии деформации по этой силе, а производная потенциальной энергии деформации по обобщенной силе равна обобщенному перемещению:



В результате получается уравнение движения Лагранжа:



__

Равновесное положение системы вала принимается за начало обобщенных координат, т.е.



Кинетическая и потенциальная энергии системы:





-

коэффициенты инерции,



– коэффициенты жесткости.

Существует форма записи обобщенного закона Гука [5,с.314], связывающая все силы и перемещения:



В условиях равновесия:





С учетом этого, уравнение Лагранжа можно записать в виде системы линейных однородных дифференциальных уравнений второго порядка с постоянными коэффициентами:



Частными решениями уравнений системы будут уравнения:



В частных решениях (j = 0, 1,2,3…s):

Частным решениям соответсвуют резонансные частоты колебаний.

Для неизвестных получают систему линейных однородных уравнений подстановкой полученного частного решения в приведенную систему уравнений (основные уравнения система малых колебаний с s степенями свободы):





Полученная система уравнений имеет решение, отличное от нуля в случае равенства нулю определителя этой системы.

На этом основании записывается вековое уравнение (уравнение частот). Вековое уравнение является уравнением s-степени относительно :





Искомые частота колебаний р и амплитуды μ, возникающие при этой частоте (k = 1,2,3…n), находятся из:

– основных уравнений системымалых колебаний с s степенями свободы,

– векового уравнения.

Вековое уравнение является уравнением s степени относительно k2. И из этого уравнения находятся все частоты свободных колебаний k системы.

Так как определитель Δk2 = 0, одно из уравнений системы при μ = 1 является следствием других уравнений системы. Последовательно подставляя в уравнения системы все полученные значения k2 получается система уравнений:



Находятся значения коэффициентов μ:





– определитель матрицы, получаемый вычеркиванием из определителя



первых столбца и строки.



– минор элемента первой строки и

j

–го столбца со знаком (-1) основного

определителя



– коэффициенты распределения равные 1.


В результате частные решения первой системы уравнений:



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

k

1

и начальной фазой β

1

.



– второе главное колебание с частотой

k

2

>

k

1

и начальной фазой β

2

.



– третье главное колебание с частотой

k

3

>

k

2

и начальной фазой β

3

.

…..



Коэффициенты определяют форму главных колебаний:



– форму первого главного колебания,



– форму второго главного колебания,



– форму третьего главного колебания,

и т.д.

Общее решение первой системы уравнений можно получить суммированием частных решений:



2s неизвестные постоянных определяются по 2s и по начальным обобщенным скоростям и координатам:



На основании приведенного выше, алгоритм полного исследования свободных колебаний системы с s степенями свободы состоит из следующих действий:

а) нахождение частот свободных колебаний k1, k2ks из векового уравнения,

б) нахождение коэффициентов распределения

в) нахождения амплитуд и начальных фаз

Применение программы MathCAD

Яблонский отмечает [15,с.143] если число степеней свободы превышает 4, то для полного решения задачи потребуется громадная вычислительная работы.

Однако, в настоящее время возможно применение математических пакетов таких как MathCAD.

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

MathCAD позволяет выполнять с матрицами символьные вычисления.

Для решения матричного уравнения типа:



необходимо записать матрицу



вставить определитель



, вызвать команду «→».

В результате получается запись многочлена из определителя. Многочлен копируется в отдельное место. Выделяют переменную «Х» в многочлене и в панели инструментов выбирают полиноминальный коэффициент. В результате этого получится матрица с коэффициентами из полученного многочлена:



Затем вызывается или записывается вручную команда polyroots, в которую добавляется полученная матрица в виде:



М1 и М2 –являются корнями матричного уравнения.

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

__

Рассмотрим пример построения эпюры свободных колебаний



Находим значение кинетической и потенциальной энергии:



Находим коэффициенты инерции и жесткости системы:



Для системы с 2 степенями свободы, уравнения частот записываются в виде:



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



Корни уравнения частот



и



определяют частоты свободных колебаний

k

1

и

k

2

(частоты главных колебаний системы).

Частота k1 (k1 < k2) является основной частотой колебаний.

Значения коэффициентов инерции и жесткости подставляются в полученное уравнение частот:



После преобразований:



В условии примера



Корни:



Значения частот k1 и k2 по результатам сопроматского расчета (см. работу Беляева [5]):



С учетом этого значения корней:



Коэффициенты распределения:





Эпюра главных колебаний:



__

Форма эпюр подчиняется теореме об узлах собственных форм колебаний [4,с.120]. По этой теореме амплитуды для разных частот колебаний не имеют одинакового знака. То есть, если амплитуда первой формы положительная, то амплитуда остальных форм должна иметь минимально одну перемену знака. Число перемен знака или число узлов собственной формы колебаний m-го порядка равно m-1.

Бабаков [16,с.124] для балки с 3 точечными нагрузками приводит три возможные формы колебаний:



__

Решение приближенным методом Релея

По методу Релея допускается:

– масса системы не изменяет типа колебаний

– перемещение системы при колебании имеют ту же форму, что и при статической деформации (сходство формы не означает равенство величин деформации).

Ошибка по методу Релея не превышает 1,5% [14,с.60].

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

__

Рассмотрим по методу Релея колебания консольной балки (вала) с защемленным концом [14,с.73].




р – круговая частота собственных колебаний в этом примере и ниже.


Обобщенное перемещение:



Кинетическая энергия груза:



в этом уравнении квадрат скорости

Кинетическая энергия элемента балки dc:



Уравнение упругой линии:



Минуя выкладки, полная кинетическая энергия системы:



Потенциальная энергия системы:



Уравнение Лагранжа:



В этом уравнении круговая р0 частота:



Статический прогиб на консоли балки:



И



Решение уравнения :

– период колебания



– частота



– круговая частота




__

Рассмотрим по методу Релея колебания двухопорной однопролетной балки (вала), нагруженной сосредоточенной силой посередине [14,с.65].




Обобщенное перемещение:



Кинетическая энергия груза:



Уравнение упругой линии:



Интегрируя последовательно:



Прогиб:



Прогиб посередине пролета:



Следовательно,



Как видно, прогибы x и xc являются динамическими прогибами, а не статическими, и имеют переменное значение, зависящее от времени.

Так, формула прогиба имеет переменное от времени значение так как сила Р, состоящая из веса груза и сил инерции зависит от времени.

Кинетическая энергия стержня:



Полная кинетическая энергия системы:



Потенциальная энергия системы:



Уравнение Лагранжа:



Эта формула аналогична формуле движения груза, подвешенного на пружине, имеющий общий интеграл .

Используя этот интеграл находим:

– период:



– частоту



– круговая частота



Если собственную массу балки не учитывать:



Т.е. к массе мешалки необходимо прибавить от веса вала.

__

Рассмотрим по методу Релея колебания двухопорной однопролетной балки (вала), нагруженной сосредоточенной силой в произвольном положении [14,с.70].



Обобщенное перемещение:



Кинетическая энергия груза:



Кинетическая энергия элемента балки dc:



Уравнение изогнутой оси балки (вала):



В точке приложения груза:







При формула имеет вид, как для предыдущего примера:



Потенциальная энергия системы:



Уравнение Лагранжа:



Для статического удлинения k необходим груз:



Находим:

– период



– частоту



– круговая частота



__

Рассмотрим по методу Релея колебания двухопорного однопролетного вала, нагруженной двумя произвольно приложенными сосредоточенными силами [14,с.76].




Ограничения метода Релея приводят систему к системе с 1 степенью свободы. При точном рассмотрении системы, она имеет множество степеней свободы.

Перемещение каждого груза:



Наибольшие перемещения грузов являются амплитудой для, для

Скорости грузов:



Максимальная скорость при



Максимальная скорость соответсвует переходу точки через статическое равновесие, т.к. фаза pt равна 0° или 180° при положении точки с на оси балки.

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

Сумма энергий постоянна и является полной энергией системы при рассмотрении идеального случая без потерь:



Для какого-либо конкретного положения системы:



При нахождении точки на оси абсцисс (оси вала), потенциальная энергия равна нулю, кинетическая максимальная:



Т.е. вся полная энергия системы является максимальной кинетической энергией.

Для фазы pt равной 90° или 270° кинетическая энергия равна нулю, а потенциальная энергия максимальная:



Т.е. вся полная энергия системы является потенциальной энергией.

Можно записать:



Для случая рассматриваемого груза:



Из этой формулы находится круговая частота:



Период колебаний:



___

Для трех грузов на валу, круговая частота запишется по формуле:



__

Для n грузов круговая частота запишется по формуле:



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

Для одного груза круговая частота запишется по формуле:



__

Рассмотрим по методу Релея колебания двухопорного однопролетного вала, нагруженной распределенной нагрузкой [14,с.81].

Мешалки являются сосредоточенной нагрузкой на валу и пример приводится для сведения.




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

Колебания системы с распределенной нагрузкой находятся по приведенной выше формуле:



Точность решения зависит от числа n участков.

Прогибы находят по уравнению упругой линии с равномерно распределенной нагрузкой:



Для 8 участков (8 прогибов):



С учетом этого, уравнение упругой линии:



С учетом того, что



__

Рассмотрим по методу Релея колебания балки на нескольких опорах [14,с.87].

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

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





Форма прогиба такая же как у статического прогиба под действием сил, применяя принцип Даламбера (приводя динамическое нагружение к статическому приложению сил).

Силы инерции вызывают дополнительный прогиб х1 и х2. Их уравновешивают дополнительные силы упругости, возникшие из-за этого прогиба.



k1 – прогиб в сечении I от силы равной 1 и приложенной в сечении I,

k2 – прогиб в сечении I от силы равной 1 и приложенной в сечении II,

k3– прогиб в сечении II от силы равной 1 и приложенной в сечении I,

k4 – прогиб в сечении II от силы равной 1 и приложенной в сечении II,

Сила инерции в сечении I:



Сила инерции в сечении II:



Сила равная 1 приложенная в сечении I вызывает прогиб k1, а сила инерции в этом же сечении вызывает прогиб:



Прогиб в этом же сечении от силы инерции, приложенной в сечении II:



Полный прогиб в сечении I:



Полный прогиб в сечении II:



Полученные уравнения для х1 и х2 являются дифференциальными уравнениями движения для рассматриваемого случая трехопорного вала.

Коэффициенты в уравнениях находятся по принципу сложения сил, по которому прогиб в любой точке вала под действием сосредоточенных сил получается в виде суммы прогибов от каждой из силы по отдельности (для прогиба в сечении I находятся и суммируются прогибы от сил Q1, Q2, RC).

Уравнение упругой линии для левой части вала (с – расстояние между правой опорой и точкой приложением силы):



Прогиб в месте приложения груза:



Находится неизвестная реакция опоры RC для статически неопределимого трехопорного вала (балки). Для нахождения реакции RC принципом сложения сил отбрасывается средняя опора вала и заменяется направленной снизу вверх реакцией RC. Так получается статически определимая система, нагруженная 3 силами: известными Q1 и Q2 и неизвестной реакцией RC. Сумма прогибов от каждой силы в точке с равна нулю так как в этой точке находится опора. И из условия равенства нулю прогибов находится реакция RC.



Прогиб от силы Q1 в точке с:



Прогиб от силы Q2 в точке с:



Прогиб от силы RC в точке с:



Вместо прогибов в формулу подставляются их значения:



Из этоф формулы находится Rc



Находится прогиб в сечении I по известной RC. Прогиб равен сумме прогибов от сил Q1, Q2, RC



Прогиб в сечении I от силы Q1 (c = la1)



Прогиб в сечении I от силы RC (c = l2 и y = a1)



Подставляя значение RC



Прогиб в сечении I от силы Q2 (c = a2 и y = la2)



Суммарный прогиб в сечении




Формула прогиба в сечении I зависит от силы Q1 и силы Q2. Группируются члены, содержащие силу Q1 c получением формулы прогиба в сечении от силы равной Q1, приложенной в сечении I:



Если в эту формулу вести Q1 = 1, то формула покажет прогиб в сечении I от единичной силы, приложенной в сечении I:



Если в полученном уравнении Q2 = 1



если в эту формулу вести Q2 = 1,



Прогиб в сечении II от силы Q1



Прогиб в сечении II от силы RC



Прогиб в сечении II от силы Q2



Полный прогиб в сечении II



Группируя члены для сил Q1 и Q2 и принимая эти силы равными 1:



Теперь решаются уравнения прогибов х1 и х2. Коэффициент k3 заменяется на равный k2.



Вал совершает гармонические колебания:



Производные этих последних уравнений по времени:





Теперь в полученные ранее формулы для х1 и х2 подставляются вторые производные:



После преобразований:



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



После группировки членов, содержащих р2 и р4:



Полученная формула решается для нахождения р2:







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



Критические скорости вала:



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

__


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

Рассмотрим однопролетный вал с силой, приложенной посередине [14,с.97].

Вал жесткий:



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

При вращении возникает центробежная сила:



Внутренняя сила упругости:



Уравнение прогиба по условию равновесия:



После решения относительно х:



Вводится обозначение:



(р – круговая частота собственных колебаний)

Получается:



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

Скорость вала, равная частоте собственных поперечных колебаний, является критической скоростью.

Критическое число оборотов вала:





Нахождение критического числа оборотов вала состоит в задаче нахождения частоты собственных поперечных колебаний.

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

Гибкий вал:




В этом случае формула изменится на формулу:



т.е. между х и e поменяется знак с «+» на «-».

Из этой формулы:



Из формулы видно, что с ростом скорости за пределом критической частоты прогиб вала стремится выпрямится. В пределе при x = e вал имеет прямую ось.

Лунц указывает [2,с.99] о доказательстве этого положения в работе Фепля и в работе Зоммерфельда.

__

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

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

Из формулы собственной круговой частоты



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

Здесь под гибким валом не понимается вал со свободно перемещающимся сечением и осью с двоякой кривизной [14,с.100].

Для изменения жесткости вала изменяют его длину, размеры сечения (инерциальные характеристики).

__

Приведем несколько отличающееся описание выкладок расчета критических оборотов вала в работе Тимошенко [18].

Тимошенко указывает [18,с.256] о возникновении критических колебаний вследствие эксцентриситета масс, возникших при изготовлении вала (биение поверхности).

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



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

Центробежная сила:



Упругая сила:



Приравнивая:



На невысокой угловой скорости с эксцентриситетом близким к нулю, прогиб незначителен. С увеличением ω прогиб увеличивается и при становится.

В этом случае угловая скорость является критической скоростью:



При превышении критической скорости формула равновесия:



(изменился знак между y и e с «+» на «+»).



Формула показывает, что с увеличением частоты, прогиб уменьшается.

После этого Тимошенко [18,с.258] принимает для анализа вала модель, в которой сам вал вращается вокруг своей оси (изогнутой оси) с частотой ω, и плоскость вала вращается вокруг прямой оси с такой же частотой ω.

В этом случае на вал будет действовать сила



Работа центробежной силы:



Из этой формулы получается такая же формула для критической частоты.

Оценивается влияние массы вала на значение критической частоты. Используется метод Релея. Задается вид кривой изгиба вала. Этим система вала преобразуется в систему с одной степенью свободы. Для вала с одним импеллером (η – прогиб):



Для нескольких импеллеров на валу:



Второй член левой части формулы относится к работе центробежной силы.

Тимошенко [18,с.260] рассматривает вал с 4 дисками:



Горизонтальные силы уравновешиваются, вертикальные силы приводятся к паре сил и силе в плоскости xy. Пара сил:



Все пары приводятся к паре (θ – момент инерции мешалки относительно оси z).

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

Формула для определения критической частоты:





Тимошенко называет приведенную формулу общим решением о разыскании критической угловой скорости [18,с.260].

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

Совместное действие поперечных и крутильных колебаний на вал

Тимошенко С.П. в работе [19,с.427] подробно рассмотрел проблему совместного действия изгибных и крутильных колебаний на балку. Для рассматриваемого им случая изгибные колебания проходили не в плоскости симметрии стержня, в результате чего возникают крутильные колебания. В нашем случае крутильные колебания возникают при вращении вала с мешалками. Однако, выводы полученные Тимошенко могут быть применены для анализа совместного действия поперечных и крутильных колебаний вала с мешалками.

Для вертикальной нагрузки кривая прогиба:



(w – интенсивность распределения поперечной нагрузки, за положительное направление принимается верх)

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

Крутящий момент:



R – крутильная жесткость, R1 – жесткость стесненного кручения.

Дифференцируя получается:



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

Интенсивность поперечных сил инерции

Интенсивность моментов инерции

Iп – центральный полярный момент инерции сечения вала.

Формулы для совместных изгибных и крутильных колебаний:




Вал колеблется в одной из собственных форм колебаний.



р – круговая частота колебаний,

Х, Х1 – нормальные функции, решения которых отыскиваются для удовлетворения граничным условиям.

После подстановки:



Тимошенко приводит пример стержня со свободно опертыми концами:



Функции Х и Х1 в этом случае:



Ci и Di – произвольные постоянные.

Вводятся обозначения:



После подстановки получается:



Решения для Ci и Di находятся в случае, если определитель уравнений равен нулю.

В этом случае частотное уравнение:



Из этой формулы:



Для случая совпадения центра тяжести с центром сдвига, то есть с = 0 и λ =0:



Из формулы получаются две системы значений частот:



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

Связанные изгибно-крутильные колебания можно найти методом Релея-Ритца [19,с.430].

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

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

Расчет изгибно-крутильных колебаний вала с мешалками по данным [20].

Рассмотрим шарнирно опертый стержень [20,с.200]. Система уравнений распадется на две независимые системы. Уравнение, описывающее только изгибные колебания в плоскости симметрии:



Уравнения, описывающие изгибно-крутильные колебания:



Граничные условия при x = 0 и x = l:



Граничные условия удовлетворяются при:



Собственные частоты определяются из формулы:



Частоты изгибных и крутильных колебаний :





Собственные частоты колебаний:



При a3 = 0 центр тяжести и центр изгиба совпадают,



Как видно, формулы Тимошенко и по справочнику [19] для определения поперечных и изгибных колебаний почти полностью совпадают.

Однако, Тимошенко указывает о независимости от и необходимости применения метода Релея-Ритца.

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

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

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

Тимошенко указывает о необходимости численного выполнения расчетов колебаний в работе [18]. То есть в том числе маститый специалист признает превосходство численных методов над ручными расчетами.

Расчет валов методом конечных элементов

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

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

Метода Даламбера удобно применять для систем с небольшим числом степеней свободы [21,с.486], к которым относятся валы с мешалками. Но вариационный подход Гамильтона является обобщением методов. Поэтому расчет вала с мешалками методом конечных элементов приведем на основе вариационного подхода Гамильтона.

Принцип Гамильтона записывается в форме [21]:



(Т и П – кинетическая и потенциальная энергии, Wne – силы демпфирования).

Функционал Лагранжа [20]:





Функционал Лагранжа по принципу Гамильтона при возможных перемещениях удовлетворяет условиям совместности и граничным условиям на контуре в течении времени от t1 до t2 и имеет стационарное значение.

Начальное положение для вариационной формулировки МКЭ следует при Т = 0 и Wne = 0:



Введем зависимости для Т, П и Wne от обобщенных перемещений, скоростей и сил [20]:



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



Для конечного элемента объема V [20]

– кинетическая энергия в матричной форме:



– потенциальная энергия (складывающаяся из внутренней энергии деформации, потенциальной энергии внешних объемных и внешних поверхностных сил):



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



Скорость связана с обобщенной скоростью:



Силы демпфирования пропорциональны скоростям (являются неконсервативными):



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



Формулы для кинетической и потенциальной энергии можно записать после преобразований в виде:





После подстановки записанных формул в первую формулу вариационной формулировки, получается матричная формулировка конечного элемента [20]:



m – матрица масс, c – матрица демпфирования элемента, k – матрица жесткости, Qe – вектор обобщенных сил в узлах конечного элемента.

В результате составляется уравнение движения системы конечных элементов на основе уравнений движения одного (каждого) конечного элемента [20]:



М – матрица масс, С – матрица демпфирования, K – матрица жесткости, Q – вектор обобщённых сил.

__

Собственные колебания вала находят решением последней записанной системы дифференциальных уравнений. Для колебаний без затухания, система запишется в виде [21,с.500]:



Матричное уравнение запишется в виде т.к.:



Уравнение имеет решение при равном нулю детерминанте системы:

Матрица массы конечного элемента записывается формулой:



Для плоского линейного элемента перемещения описываются полиномами Гермита [21,с.491], матрица жесткости запишется:



После преобразований [20]:





Для конечного элемента, показанного на рисунке выше, с нагрузкой вдоль оси и с узлами на концах, с применением линейных интерполяционных функций, матрица масс записывается в виде [21,с.492]:



Запишем формулу для матрицы жесткости.

На рисунке показан стержневой элемент под действием изгиба [21,с.69]:



Вектор параметров перемещений в узлах элемента имеет два перемещения и два вращения:



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



Угол

Перемещения и вращения на концах стержня:



Матрица С [20,с.70]:








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



Делитация связана с перемещением:



Для вектора деформации:



(составляющие деформации в зависимости от составляющих перемещений находятся применением матрицы оператора над матрицей интерполяционных функций).







L – матрица-оператор, для плоских задач

Ar – матрица интерполяции

Для матрицы интерполяции могут быть приняты функции вида:



По уравнению :









Пропуская математические выкладки, получается:



Для конечного элемента так как перемещения на концах равны нулю, матрица жесткости записывается в виде [20,с.505]:





Теперь, подставив в уравнение матрицы получится:



Вводится обозначение:



Характеристическое уравнение:



В виде многочлена (см. о решении уравнений в программе MathCAD):







Для случая б), т.е. для второй части на рисунке выше, перемещение в узле 1 и вращение в узле 2 равны 0. С учетом этого матрицы k и m уменьшаются:



Характеристическое уравнение:



В виде многочлена:







Эпюра собственных колебаний вала:



__

Итак, в разделе показаны теоретические основы расчета методом конечных элементов валов на свободные колебания.

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

Расчет по методу МКЭ является более обоснованным теоретически и позволяет выполнять расчет валов с опорными узлами любой конфигурации.

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

__

Стандартом по умолчанию является программа ANSYS, описанная в работах [24], [25], [26]. Может быть использован пакет [27].

Заключение

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

Предложена горизонтальная установка погружного насоса для анализа применения в нефтепереработке.

Подробно изложена теория расчета валов насосов на резонанс, так как в классической литературе такие сведения отсутствуют.

Библиография

1. Ефанов К.В. Блоки нефтяных аппаратов. – М.: Литрес, 2020. – 27 с.

2. Капустин В.М., Рудин М.Г., Кудинов А.М. Основы проектирования нефтеперерабатывающих и нефтехимических заводов. М.: Химия. 2012. 440 с.

3. Капустин В.М., Рудин М.Г., Химия и технология переработки нефти. – М.: Химия, 2013. – 496 с.

4. Капустин В.М. Технология переработки нефти. В 4-х частях. Ч.1. Первичная переработка нефти. – М.: Колос, 2012. – 456 с.

5. ОСТ 26-1141-74 «Насосы. Основные требования к установке и эксплуатации вне помещений на химических, нефтехимических и нефтеперерабатывающих производствах». Минхиммаш СССР, 10.02.1975.

6 Ефанов К.В. Теория расчета нефтяных центробежных насосов. – М.: Литрес, 2020. – 28 с.

7. Айзенштейн М.Д. Центробежные насосы для нефтяной промышленности. М.: Гостоптехиздат. 1957. 363 с.

8. Михайлов А.К., Малюшенко В.В. Конструкции и расчет центробежных насосов высокого давления. М.: Машиностроение. 1971. 304 с.

9. Михайлов А.К., Малюшенко В.В. Лопастные насосы. Теория, расчет и конструирование. М.: Машиностроение. 1977. 288 с.

10. Малюшенко В.В. Динамические насосы. – М.: Машиностроение, 1984. -84 с.

11. Ивановский В.Н., Пекин С.С., Сабиров А.Л. Установки погружных центробежных насосов для добычи нефти. – М.: Нефть и газ, 2002. – 256 с.

12. Богданов А.А. Погружные центробежные электронасосы для добычи нефти (Расчет и конструкция). – М.: Недра, 1968. – 272 с.

13. Мясников А.Н. Горизонтальные насосные системы Шлюмберже в России // Инженерная практика, 2015 – N06-07.

14. Лунц Е.Б. Упругие колебания. – М.: Изд-во МАИ, 1935. – 182 с.

15. Яблонский А.А., Норейко С.С. Курс теории колебаний. 3-е изд. – М.: Высш. шк, 1975. – 248 с.

16. Бабаков И.М. Теория колебаний. 3-е изд. – М.: Наука, 1968. – 560 с.

17. Беляев Н.М. Сопротивление материалов. – М.: Наука, 1956. – 854 с.

18. Тимошенко С.П. Прочность и колебания элементов конструкции. – М.: Наука, 1975. – 704 с.

19. Тимошенко С.П., Янг Д.Х., Уивер У. Колебания в инженерном деле. – М.: машиностроение, 1985. – 472 с.

20. Вибрации в технике: Справочник. В 6-ти т. / Ред. совет: В. Н. Челомей (пред.). – М.: Машиностроение, 1978 – Т. 1. Колебания линейных систем/Под ред. В. В. Болотина. 1978. – 352 c.

21. Секулович М. Метод конечных элементов. – М.: Стройиздат, 1993. – 664 с.

22. Магалиф В.Я., Иткина Д.М., Корельштейн Л.Б. Монтажное проектирование химических, нефтехимических и нефтеперабатывающих производств. – Москва, 2010. – 346 с.

23. Магалиф В.Я. Теоретические основы конструирования трубопроводов (курс лекций. – Москва, 2011. – 35 с.

24. Бруяка В.А., Фокин В.Г., Солдоусова Е.А., Глазунова Н.А., Адеянов И.Е. Инженерный анализ в ANSYS Workbench.: Учеб. пособ. – Самара: Самар. гос. техн. ун-т., 2010. – 271 с.

25. Бруяка В.А., Фокин В.Г., Кураева Я.В. Инженерный анализ в ANSYS Workbench.: Учеб. пособ. – Самара: Самар. гос. техн. ун-т., 2013. – 149 с.

26. Пересыпкин В.П., Пересыпкин К.В., Иванова Е.А. Проектирование силовых конструкций ракет-носителей с применением метода конечных элементов. Самара: Самарский гос. аэрокосм. ун-т, 2012. – 95 с.

27. Алямовский А.А. SolidWorks Simulation. Как решать практические задачи. – СПб.: БХВ-Петербург, 2012. – 448 с.


Оглавление

  • Введение
  • Монтажная привязка насосных агрегатов на местах установки
  • Привязка насоса в технологической документации
  • Технологические схемы обвязки насосов
  • Система КИПиА насоса
  • Расчет характеристик насосов
  • Горизонтальный агрегат погружного насоса
  • Вертикальные полупогружные насосы типа ХП.
  • Расчет на резонанс валов насосов
  •   Расчет валов на резонанс по теории колебаний
  •   Совместное действие поперечных и крутильных колебаний на вал
  •   Расчет валов методом конечных элементов
  • Заключение
  • Библиография




  • MyBook - читай и слушай по одной подписке