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

Дыбрин, Александр Андреевич

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

кандидата технических наук: Дыбрин, Александр Андреевич
город: Ижевск
год: 2014
специальность ВАК РФ: 05.13.01

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

Автореферат диссертации по теме "Анализ и моделирование уровня вибраций и шума в трубах переменного сечения газопроводов"

На правах рукописи

ДЫБРИН Александр Андреевич

АНАЛИЗ И МОДЕЛИРОВАНИЕ УРОВНЯ ВИБРАЦИЙ И ШУМА В ТРУБАХ ПЕРЕМЕННОГО СЕЧЕНИЯ ГАЗОПРОВОДОВ

Специальности:

05.13.01 - Системный анализ, управление и обработка информации (в пауке и технике) 01.02.06 -Динамика, прочность машин, приборов и аппаратуры

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

1 и АПР 23М

005546948

Ижевск 2014

005546948

Работа выполнена в ФГБОУ ВПО «Ижевский государственный технический университет имени М.Т. Калашникова» (ИжГТУ им. М.Т. Калашникова)

Научный руководитель:

- Лялин Вадим Евгеньевич, заслуженный изобретатель РФ, доктор технических наук, доктор геолого-минералогических наук, профессор

Официальные оппоненты: -

Костарев Сергей Николаевич, доктор технических наук, доцент, ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет», профессор кафедры «Информационные технологии и автоматизированные системы»

Миловзоров Георгий Владимирович, доктор технических наук, профессор, Институт нефти и газа им. М.С. Гуцериева ФГБОУ ВПО «Удмуртский государственный университет», зав. кафедрой бурения нефтяных и газовых скважин

Ведущая организация:

- ФГБУН «Институт механики УрО РАН»

Защита состоится «24» апреля 2014 г. в 12:00 часов на заседании диссертационного совета Д 212.065.06 в ИжГТУ имени М.Т. Калашникова по адресу: 426069, г. Ижевск, ул. Студенческая, 7.

С диссертацией можно ознакомиться в библиотеке университета и на официальном сайте http://istu.ru. С авторефератом можно ознакомиться по адресу http://vak.ed.gov.ru/ш/с^еПаиоп.

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

Автореферат разослан «21» марта 2014 г.

Ученый секретарь диссертационного совета, кандидат технических наук, доцент

В.Н. Сяктерев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

приводит к ускоренному их износу.

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

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

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

Магистральные газопроводы (МГ) относятся к опасным производственным объектам. Часть оборудования, например газоперекачивающие агрегаты (ГПА), являются не только источником потенциальной опасности, но и источником'образования вредности. Следствием их работы является высокий уровень шума (90... 130 дБ) аэродинамического и механического характера, создаваемый истечением рабочей среды в газовоздушных трактах всасывания и выхлопа, которые, с точки зрения акустики, представляют собой своеобразные

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

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

Степень научной разработанности темы исследования. Исследования механикн пространственно-криволинейных трубопроводов и фланцевых соединений рассматриваются в работах Светлицкого В.А., Власова В.В., Гольденвей- " зера A.J1., Башты О.Т., Аксельрада Э.Л., Айнбивдера А.Б., Камерштейна А.Г., Герштейна М.С., Самарина A.A.

Теоретическими и экспериментальными исследованиями снижения уровня шума на компрессорных станциях занимались многие известные ученые: Кравчун П.А., Занченко В.И., Терехов А.Л., Рейнольде А.Ж. Леонтьев В.А., Григорян Ф.И., Юдин Е.А., Бэтчерлор Д., и др. Созданные ими теоретические положения и технические средства позволяют снижать уровень шума. Однако при многообразии конструкций оборудования компрессорных станций эти разработки не обеспечивают необходимую эффективность вследствие невозможности учета всех факторов, влияющих на шумообразование.

Область исследования. Диссертационная работа выполнена в соответствии с пунктами «6. Методы идентификации систем управления на основе ретроспективной, текущей и экспертной информации», «11. Методы и алгоритмы прогнозирования и оценки эффективности, качества и надежности сложных систем» паспорта специальности 05.13.01 - Системный анализ, управление и обработка информации (в науке и технике), и пунктами «1. Динамика машин, приборов, аппаратуры, систем и комплексов машин и приборов», «9. Математическое моделирование поведения технических объектов и их несущих элементов при статических, динамических, тепловых, коррозионных и других воздействиях» паспорта специальности 01.02.06 - «Динамика, прочность машин, приборов и аппаратуры».

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

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

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

Для достижения цели требуется решить следующие задачи:

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

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

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

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

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

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

Методы исследования. В диссертации использованы теоретические и экспериментальные методы исследования.

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

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

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

На защиту выносятся результаты разработки подходов, методик и математических моделей статики и динамики пространственно-криволинейных трубопроводов и снижения уровня шума в МГ, в том числе:

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

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

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

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

- сопоставление распределения импеданса по частоте и по длине трубы, найденное из эксперимента по предложенному методу, с теоретическим распределением импеданса, рассчитанным по уравнениям акустики;

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

Научная новизна результатов диссертационного исследования заключается в следующем:

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

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

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

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

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

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

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

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

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

Апробация работы. Результаты работы докладывались на: IV международной научной конференции «Актуальные вопросы современной науки» (Санкт-Петербург, 2012); VII Международной научно-практической конференции «Актуальные вопросы науки» (Москва, 2012); XVI молодежной международной научно-практической конференции «Интеллектуальный потенциал 21 века: ступени познания» (Новосибирск, 2012); международной научно-практической конференции «Техника и технологии: роль в развитии современного общества» (Москва, 2013).

Публикации. Основные полученные автором научные результаты отражены в 14 научных публикациях общим объемом 3,9 п.л., авторский вклад -3,5 п.л. Автор имеет 6 научных трудов в изданиях, выпускаемых в РФ и рекомендуемых ВАК для публикации основных результатов диссертаций.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы из 145 наименований. Работа изложена на 164 страницах, содержит 43 рисунка и 2 таблицы.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Во второй главе исследованы динамика фланцевых соединений трубопроводов (рис. 2) и вопросы определения взаимных перемещений уплотняемых поверхностей и деформаций элементов соединения. Изолирующее фланцевое соединение состоит из двух основных фланге. 2. Фланцевое соединение цев 12 и 13, приваренных к концам

трубопроводов с плоской прокладкой:

г н газопровода, и третьего специально-

го фланца 11 толщиной 16-20 мм.

Предлагается подход, согласно которому трубопровод с фланцевым со-

I — газосборные ссги: 2 — промысловый пункт сбора газа: 3 — головные сооружении: 4—компрессор! шя станция: 5 — газораспределительная станция; в — подомные хранилища; 7— магистральный трубопровод: в—ответвление: 9—лниейнал арматура: 10—дяуипгточный иереход через воднут преграду Рис. I. Схема магистрального газопровода

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

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

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

Уравнения продольных колебаний участков 1 и 2 трубопровода имеют вид: а2и, /сЯ2 = а2 д\ /дх;; д2и2 ¡д(г = а\ д\¡дх\ , (1)

Где а, = 4Е\!Р\ ^ «2=7^7 Pj, - соответственно плотность и модуль упругости материала трубы у - го участка трубопровода, у = 1,2.

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

м

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

E^dujdx,)^ = -m{d2uJdt2)^ -с[(«,)Х1Ц +(«2¡u]; (2) E2F2(du2/dx2)u__,2 =-m{d\/dt2)xi i +(",)„-,,]'

где с - жесткость фланцевого соединения; 2т - масса фланцевого соединения; Fj, lj - соответственно площадь поперечного сечения и длина у-го участка трубопровода, у = 1,2.

Собственные колебания динамической системы находились в виде: ы, =-Vl(^l)(/4lcosp/ + 51sinpí); и2 = X2(x2)(A2cospt + B2s\npt), (3)

где Xj{x¡) - собственная форма колебаний ./-го участка трубопровода; Ajy Bj - произвольные постоянные, определяемые по начальным условиям, у = 1,2; р - частота собственных колебаний системы.

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

= eos(pxt/о,) + D, sin(рх,/а,); Х2 (х2) = С2 cos(/w2/а2) + D2 sin(рх2¡a2) ;(4) где С/( D( - постоянные, зависящие от граничных условий, у = 1,2. Для удовлетворения первых двух граничных условий (2) необходимо в выражении (4)

принять С, =С2 =0. Из двух последних граничных условий (2) с учетом (3) и (4) после несложных преобразований получаем

[£¡f¡ р/о, сск(р/,/йг,) +csin(p/,/a, )-пр1 sin [pija, )]d,H, (/) +csin(plJa2)D2H1(t) =0; (5) csm{plJa^)[\Hx {t)+{E1F1 p/a2ax(p!2/a2)+cs\n(pl2/a2)-mp2sm(pl2/a2)]jD1H2(t) =0.

Выводим частное уравнение:

{^Фi «»(p'iM) +csin(p/1/<al) -»p2 sin(p/,/a,)) x , (6) pfa cosl^/flj) +c-sin(p/,/a,)-»p2 svn^pljc^)) -c2 sin( Д /ú'1)sin(p/2/a2)=0

В случае изотропной системы, что наиболее часто встречается в реальных условиях эксплуатации трубопроводных систем, £, = Е2 = Е; Ft-F2=F; р{ - р2 = р \ а, = а, = а. Тогда частотное уравнение (6) можно представить в следующем виде: тр2{тр2 -2c)/g(p/,/fl)#(p/2/a) +EF/a(c-mp2){tg(plja) + tg{plja)) +E2F2/a2 р = 0.

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

И, (*,.') = Y Sin-^(/4,,' eos P,t + Bj sin p,t); м Я

* O JC

U2 (X2 >0 = ZS¡n C0S A' + ^2,' SÍn P¡t) ,

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

В частном случае, если система симметрична, т.е. /, = /2 = / = l¡2 (рис. 3), где F - межопорное расстояние, зависимости (5) после несложных преобразований можно привести к следующему виду:

a = ptgp-,a = {p-k//3)tgp, (7)

где а = pFl/m; р = pl/a; к = 2ср1г/(тЕ).

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

Поэтому первое частотное уравнение (7), соответствующее кососимметричным колебаниям, может быть получено на основе отдельной динамической модели (рис. 4а). Второе частотное уравнение (7) соответствует симметричным колебаниям и может быть получено на основе динамической модели, изображенной на рис. 46.

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

1

2с

W

Рис. 4. Расчетная схема трубопровода с фланцевым соединением; а) для кососимметричных продольных колебаний; б) для симметричных продольных колебаний

Численные расчеты частного уравнения (6) проведены на компьютере в широком диапазоне вариаций конструктивно-технологических параметров рассматриваемой динамической системы. Расчеты проводились для однопролетного трубопровода, выполненного из трубы, предназначенной для рабочих давлений до 20 МПа, параметры которой, как и параметры фланцевого соединения, соответствовали соответствующим отраслевым стандартам. Некоторые из этих расчетов, полученные для различных диаметров условного прохода - Ду с жестким закреплением концов трубопровода-на опорах и межопорном расстоянии 1 = 2м представлены графически на рисунках 4-6.

Р> с

7500

66001

620С

58°Ъ 50 ДуГмм

Рис. 5. Зависимость частоты колебаний от диаметра условного прохода трубопровода и фланцевого соединения; продольная жесткость фланцевого соединения С = 6 • Ю4 Н/М ; 1) фланцевое соединение у опоры; 2) фланцевое соединение в середине пролета

р. «

7500

7000

6500

/ / \

in« /-> :г/.<

Рис. 6. Зависимость частоты основного тона колебаний трубопровода с фланцевым соединением от его жесткости, фланцевое соединение у опоры; диаметры условного подхода: 1) Ду 32мм; 2) Ду 50мм; 3) Ду 80мм; 4) Ду 100мм

6400

\ ,1

) \\

\

<

0,1

0,2

0,3

0,4 LJL

ог о.з 6

Рис. 7. Зависимость частоты собственных продольных колебаний трубопровода от месторасположения фланцевого соединения относительно опор; диаметр условного прохода Ду 100мм; продольная жесткость фланцевого соединения: Т) 5 • 1010 Н/м; 2) 1 • I о10 Я/л«; 3) 4■10' И/м; 4) 2,5•10" Н/м; 5) 2 • 10" Н/м; а) для первого тона колебаний; б) для третьего тона колебаний

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

Уравнение свободных поперечных колебаний участков 1 и 2 трубопровода без учета поглощения энергии имеют вид:

E,l, tfyjdx] + PiF. d'yj/dt2 = 0; j = 1,2, (8)

где Ej - модуль упругости у-ого участка газопровода; /(. - момент инерции поперечного сечения у-ого участка трубопровода; р; - плотность материала трубопровода; F. - площадь поперечного сечения.

Граничные условия на концах участков 1 и 2 трубопровода для шарнирного закрепления на опорах (рис. 8) имеют вид:

(9)

-Ы,^} (Ю)

где с,, с, - поперечная и изгибная жесткость фланцевого соединения; М/ - сосредоточенная масса фланцевого соединения, приведенная к j-му участку трубопровода; у = 1,2.

Решение уравнений (1) ищем в следующем виде;,.

(id

В соответствии с выражениями (8) и (11) для данного случая уравнения движения трубопровода с фланцевым соединением можно представить в следующем виде:

d*hj (0/Л1 + р2Я,. (0 = 0; d*X\ (*, )/dx] - к'Х, {Xj) = 0; j = 1,2, (12) где р - частота собственных колебаний; к - частотный параметр

k = ijp1pF/(EI). (13)

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

xJ (*> 4ñs{kxt)+ 4Пт{кх,)+ /Í<%(A*,)+ = 1,2, (14)

где - произвольные постоянные, / = 1,2,3,4; S,T,U,V - функции А. Н. Крылова.

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

yí{xj't) = ÍXÁxJ-P,.)HlÁ^P,lY / = 1,2, « = 1,2,3..., (15)

«-i

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

XÁ*J'P.) = ^(¿.х^^'У^л,); j = 1,2, (16)

где / = 2,4 - производные постоянные /;-ой собственной формы.

Функции Hjn{t,pn)', j- 1,2 согласно первым двум уравнениям системы (12) подчиняются зависимостям

HJ„^PJ = B{¿)smpJ + B¡Í)cosp„n j —1,2,

где B\'J, - произвольные постоянные /г -ой формы собственных колебаний, определяемые по начальным условиям.

Таким образом, общее решение (15) уравнений (8) можно представить в следующем виде:

Ух М = ¿[^(М,) + 5иу (кл )](/?,<:> sin PJ + /?_!'„) eospj\

У г M = l[32„T(k„x2) + S}iiV(k„x2)] (¿e sin PJ + eos PA

Аналогично случаю шарнирного закрепления трубопровода на опорах получены частотное уравнение и общее решение уравнения (8) в случае жесткого закрепления.

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

Рассмотрена в общем виде задача о прохождении звука через трубу переменного сечения длиной d, соединяющую две полубесконечные трубы с площадями переменного сечения 5, и 5, (рис. 9).

Положим, что труба переменного сечения 2 расположена между точками

Гс, .2

'sJS^/

ЦлаХ

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

х = аих = Ь. Ее длина <1 = Ъ-а. Закон изменения площади поперечного сечения описывается функцией Б(х). Во всех трубах находится среда, характеризуемая

плотностью р и скоростью распространения звука с.

Для простоты рассмотрены гармонические колебания с круговой частотой со. Одновременные установившиеся колебания среды в трубе переменного сечения описываются уравнением Вебстера

л2 ад ад РЛ '

где р2(х) - звуковое давление;

к-(о 1с - волновое число; $'(*) = (х)/с£х - производная . Общее решение (20) можно записать в виде линейной комбинации двух функций <р,(дг) и <р2 (х), удовлетворяющих (20),

5

Труба! 0 а Труба II Ь Труба III д,

Рис. 9. Схема перехода между трубами разного поперечного сечения

р2{х)=См{х) + Сг(рг{х), (21)

где С, и С2 - постоянные, которые определяются из граничных условий.

Решение (21) образует фундаментальную систему тогда и только тогда, когда их определитель Вронского не равен нулю:

<Р,(кх) <р2{кх)

™[(рх{х),(рг{х)~] =

(22)

<р1(Ь) <р'(кх)'

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

Г,3=ЗД|ДЙ,/Д|\ (23)

где Д0и и Д определяются по формулам:

Д = [>, (ка)<р'2 (кЬ) + <р\ (ка)<р, (кЬ)-(р[{ка)<р2 (кЬ) - <р2 (ка)<р1(кЬ)] +

{ка)(р[(кЬ) + ср2 (ка)(рх (кЬ)- <р, (ка)<р2 (кЬ) - <р[(ка)<р'2 (№)];

МкЬ)

<р[{кЬ) (р2 (кЬ)

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

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

Л0 =2

= 2 {(р,<р2-<рг<р\) = 2"[>1 (**)]•

, где k0 = р/2, а при

q = А2 / 2рс, А - амплитуда давления бегущей плоской волне, не меняет коэффициента звукопрохождения по энергии г (соотношение подобия). Установлено, что не имеет значения, с какой трубой переменного сечения мы имеем дело - сужающейся или расширяющейся. Если только у них длина, отношения площадей поперечного сечения на концах и закон изменения S(x) будут одинаковы, то и коэффициенты звукопрохождения по энергии г будут равны.

В качестве примера применения полученных выражений рассмотрена труба переменного сечения, у которой площадь поперечного сечения меняется по экспоненциальному закону, т.е. 5(х) = S^' (экспоненциальный диффузор). При длине диффузора d 53 = S\epd, откуда р = 1п (53 'S, )/d. Получено выражение для коэффициентов звукопрохождения по энергии г = ly^l+A:2/^2 -/c2)sh2^/r(,

и звукоизоляции при к>к0 R = lOlg^l л-^Цк1 -£02 )sin 2{jk2 -k2dj

*->*„ /? = 101g[l + ^2] = 101g[l + (ln(5,/5l))Y4].

У экспоненциального диффузора существует частота <u0 = к0с = рс/2, называемая критической. Ниже нее диффузор ведет себя как фильтр низких частот, т.е. плохо пропускает шум на низких частотах. Выше а>0 звукоизоляция R —» 0 с увеличением частоты.

На практике более технологичным является конический диффузор, поэтому в работе выполнены расчеты его звукоизоляции и проведено сравнение ее со звукоизоляцией экспоненциального. У конического диффузора S(x) = Ах2, где А - некоторая постоянная, характеризующая степень раскрытия диффузора. Координаты а тлЬ входа и выхода диффузора связаны с длиной d и отношением площадей выходного 5, и входного 5, сечений a -S}/St соотношениями:

Определен коэффициент звукопрохождения по энергии

г = Ьг1аг |Д0/Д|2 = (1 + ((iW)2 + (1 + Ak-ab) sin2 kd - 2/Wsin tóeos tó) /(4k2a2b2))"'

и величина звукоизоляции R =-l01gr. (24)

Анализ выражения (24) показывает, что на частотах / <c/2d звукоизоляция меняется плавно. Затем появляются максимумы и минимумы в значениях R. На низких частотах, когда f < c/^lnd), (24) после преобразований можно записать как

/г*1018[|+5,(ЗД-1)7(45,)]

Следует отметить, что на низких частотах, когда по длине диффузора нет резонансов, звукоизоляция определяется только соотношением а = Si/Si и расчет R можно вести по графику рис. 10. Только тогда, когда / > cjljtd, для расчетов нужно пользоваться формулой (24).

R.4B 6,0 -

« = S,/S,

5,0 4.0 3,0 2,0 1,0

1 2 3 4 5 6 7 8 9 10 11 12 1314 15 16

■—«___

N

V5

N

s ■с4 \

\ \ \

3>

\Л

10"'

10"

Рис. II. Зависимость звукоизоляции от безразмерной частоты

Pua 10. Зависимость звукоизоляции от отношения площадей поперечного сечения труб

Проведены численные расчеты звукоизоляции различных типов диффузоров для пяти отношений а = = 2, 4, 6, 8, 16 (кривые 1-5 соответственно). Результаты их приведены на рис. 11. Здесь по оси абсцисс отложена безразмерная частота kd, по оси ординат - звукоизоляция R в дБ. Как видно из графиков, R тем больше, чем выше отношение S}/St. На частотах, которые выше критической о\, звукоизоляция R быстро падает до нуля.

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

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

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

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

t.z +1 Hu(xi,z)+H2J(x3,z)-h

(25)

/.С (1Нп {х3,г)/с1хз + йН2Ъ (дг3,7)/еЬс3-Ип (г)' где Л,2(г) = [/''(;с2'0]/[^'(л:|'0] " передаточная функция от пульсаций давления в точке х, к пульсациям давления в точке х2;

H»{*)=sh ffu{z) = sh

z(x2-x,)/cyJl + [/(t.z)

Фз -X, )/cy]l + l/(t.z)

/

sh

z{x2-xx)j'cj\ + \j{t.z) z(x2-xt)/cJl + l/(Kz)

z = a + iû> - комплексная переменная; со - циклическая частота; а - ре-

альная часть переменной z; i -мнимая единица;

f. = Gl/[2sAPg) - характерное время сил вязкостного сопротивления, где

G - средний весовой расход;.? - площадь поперечного сечения; g - ускорение силы тяжести; Ар - перепад давления на длине L.

Найдена зависимость между импедансом на мнимой оси z = /со при (о = й>0 и соотношением между амплитудами :

/(«ч ) = (1 + ту. )/(ту.) • (E{ico0) +1 )/{E{ico0) -1), (26)

где

Е (<Ч) = ( А ♦ /А-) ехР [<Ч {9* ~ 9. - 1к0х)] (27)

или зависимость, обратную (27), (28)

E(iw0) = {[ту. / (iw0)/'(1 + ico„t, )] + l}/{[icoatJ (ícü0 )/(1 + ia)0t.)] -1} (28)

Вводя обозначение íúj0/./(I + íoy,) /(í'ü)0) = re(&>0) + iim(&>0), получим расчетную формулу для A А,

AJA_= |)| = ^[Яе(ш0)+1]Ч[1т(а)0)]7¿Re(«0)-1]2 н-^шЦ,)]2 (29)

Нетрудно убедиться в том, что формула (29) характеризует связь импеданса со структурой волны. Так, импеданс стоячей волны, которая возникает при равенстве амплитуд Af и А , может быть только мнимым, если потери на вязкостное сопротивление несущественны (ау. » 1), и, наоборот, импеданс бегущей волны в приближении (&у. »1) может быть только действительным числом, равным +/ либо -/, в зависимости от направления распространения волны.

Такие же выводы следуют и из формулы (25), по которой импеданс получается мнимой величиной в случае стоячей волны, характеризующейся действительной передаточной функцией й12, и получается равным ±/ в случае бегущей волны, для которой передаточная функция hn описывает лишь транспортное запаздывание hn = exp[±ia)(x2 - xt )/с].

В общем случае неравенства амплитуд = А_ + Д (для определенности Д >0) часть волны, бегущей в положительном направлении с амплитудой, равной амплитуде А_, накладываясь на волну, бегущую в отрицательном направлении, дает стоячую волну, а часть волны с амплитудой Д остается бегущей акустической волной. В этом случае величина (&/2А)2 = (|£|)2/4 описывает распределение энергии акустической волны между бегущей и стоячей компонентами, а знак величины ([£■(-1) характеризует направление распространения

бегущей компоненты (положительное направление при (|£j-l)>0 и отрицательное при (|£|-l)<0. Из формулы (29) видно, что направление распространения бегущей компоненты однозначно определяется знаком комплекса («„/. Re/ —Im /) или знаком реальной части импеданса в случае слабого влияния сил вязкостно'гЪ сопротивления на акустическую волну (ау. »1).

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

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

Задача распознавания этих случаев решается установкой двух датчиков пульсаций давления на трубопроводе, примыкающем к источнику колебаний (для определенности считается, что источник при этом находится слева). Определяя из эксперимента передаточную функцию между пульсациями давления и подставляя ее в (25), рассчитаем акустический импеданс, а затем по формуле (29) - соотношение между амплитудой волны ( А4), излучаемой источником вправо, и амплитудой волны (А_), отраженной акустической системой, расположенной правее датчиков давления. При этом возможны три исхода:

1) » 1 - акустическая система слабо отражает звуковые колебания,

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

2) А+/А_ = 1 - акустическая система полностью отражает звуковые колебания и поэтому образуется стоячая волна, амплитуда которой определяется не только мощностью источника, но и близостью системы к резонансу.

3) А./А_>\ - промежуточный случай, когда акустическая система частично отражает звуковые колебания, излучаемые источником, и результирующее поле акустических колебаний представляет собой сумму бегущей и стоячей волны. Амплитуда результирующего колебания определяется в этом случае как мощностью источника, так и близостью системы к резонансу.

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

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

Звуковое давление, создаваемое излучателем, обозначим давление звукового поля, отраженного от оболочки, - р2, а давление внутри оболочки -рг Если считать стенки абсолютно жесткими, то г- составляющие скорости колебаний на них должны равняться 0. Таким граничным условиям удовлетворяют функции со%к„х, где кт =«;;г// и т = 0,1,2,... - целое число.

\чЧ\\\\ч

1

Рис. 12. Оси координат и обозначе- Рис. 13. Обозначения углов при переносе

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

Источник создает звуковое поле (рис. 12, 13), давление которого в точке (r,<p,z) записывается в виде

Р] (r,V, z) = xz pLJSKsV"* cos A:.2, (30)

= 4^1 ~ k2m ; kM = ткЦ, hi = 0,1,2,..., kh = (o/ch - волновое число в воздухе. Рассеянное и прошедшее сквозь оболочку поле запишем в виде

= (31)

ияО П--Г: Х- 'С

где - функция Ганксля первого рода порядка к, Jn - функция Бесселя. Звукоизоляцию оболочки Я определим как

/? = 101ё[р1(г,<р,2)/р3(г,<р,г)У. После преобразований получено

/г = ioig

2 phco

где Z - импеданс оболочки

I Z PlJAK„r)ehl«coskmzn

IS

(32)

Z„„, = -'«»'о x

1--

r \2

it: +

•4 у

я»

I 2 2

Km + 2 K0

"l ■ J

(33)

где к1 - о>/с,; кп=а>/сй\ к п = (о/сП ; Аи = ^¡со2т0/ В - волновые числа изгибных колебаний; В = £/гл/12(1-сг2) - изгибная жесткость материала, из которого изго-

товлена оболочка; Е и а - модуль Юнга и коэффициент Пуассона материала; с,,ся,с0 - скорости распространения сдвиговых и продольных в стержне и пластине волн соответственно.

Выражение (32) является сложным для практического применения. Однако оно значительно упрощается, если рассмотреть звукоизоляцию на оси оболочки, когда все 7„(0) = 0, за исключением Jll(0)=\. Тогда исчезает суммирование по п, т.е. зависимость от угла <р, и (32) можно записать в более простом виде:

Я = 10^

Х-

т=0

1 +

2 рьа

Ао^,2_____________-...........

Чктак)Н\\кгтак)

(34)

р(Я.<г,7)

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

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

В качестве источника выбран участок полуцилиндра длиной / и радиусом а0, на котором задано произвольное распределение радиальных смещений поверхности IV = (У0(г,<р). На остальной части полуцилиндра IV = 0. Кожух радиусом ак имеет бесконечную протяженность и расположен коаксиально с источником. Для полуцилиндров выбрана ось г цилиндрической системы координат (). Начало координат расположено в середине источника шума.

Радиальные смещения источника представлены в виде

Рис. 14. К расчету звукоизоляции полуцилиндра

(35)

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

ос *

-«,(1-0

Звуковое поле снаружи кожуха р](г,г,ер), внутри -р2(г,г,<р) представлены в виде разложений

х

мы?)=х I слтч'ф^-к^^^чк-,

*

р2(г.2,<р) = £ ¡[АЛкЖ'^к; - к2г)+ Вг,(к)Н)2](фг -к2г)]епГ"> ёк\

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

х

и(г,<р) = £ _[ ип ¿/*;

(37)

Здесь и, у и н> - осевое, азимутальное и радиальное смещения, =о>/с, с- скорость распространения звука в воздухе, и //,',2)- функции Ганкеля первого и второго рода порядка п . Неизвестные амплитуды ап, вп, сп и И'* определяются из уравнения движения кожуха и граничных условий, выражающих равенство радиальных смещений частиц воздуха и источника при г = а() и частиц воздуха и кожуха при г = ак.

Уравнения движения кожуха записаны в виде:

1 + сг д и сг о и'

ди 1 - ад2 +

кдг2 2а]д(р

со

О

2 ак дгд(р ак &

= 0,

1+<т дги (1-ст д2 1 д2 со*

-+ _ -Т+-Г—Г+-

2ак дгд(р 1 2 дг2 а; д<р2 с{

о ди 1 ди + ак дг а] д<р

со

■у

с:

/г 12

' д2 15

О /

1 дп п

и---— = 0,

Я- дер

(38)

Й 2 2 -> 2 г ак оср

р"с0

Здесь с0 скорость распространения продольных волн в

пластине; £, ст и р„, - модуль Юнга, коэффициента Пуассона и плотность материала, из которого изготовлен кожух; и, и и и'- составляющие вектора смещения нейтральной поверхности кожуха по осям координат.

Граничное условие на излучателе

(39)

где р - плотность среды (воздуха).

На тонком кожухе радиальные смещения частиц среды внутри и снаружи

равны w кожуха:

Решение этой системы можно значительно упростить, если воспользоваться выражением для импеданса тонкой цилиндрической оболочки, который определяется как отношение амплитуды силы, приложенной к оболочке, к амплитуде радиальнои скорости и; для бегущей волны е'1^"""1 (33). Окончательно звукоизоляция кожуха

Л-2 01g|1 + ^COs2Q^fa>x 4 р0соак Hw(aj

(41)

Для удобства проведения расчетов по формуле (41) можно воспользоваться показательной формой записи функции Ганкеля:

н]и(х) = 1С',Хх)в^\ н]\х) = -«С'„(х)в(42) Здесь точка над функциями Ганкеля обозначает производную по аргументу. С учетом этих выражений (41) записывается в виде

Шк1ак со^е^^^ду.

R„ - 20 lg

I--

2 рссо

(43)

где

кксаг cos©

2я к a, cos©

У

Я .г 1^/2 2\ л 71 ACOS0\2»/V' „г«\

При больших значениях аргумента

2

С:(*Ас os©)«

тт^я, COS©

ЗАКЛЮЧЕНИЕ

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

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

3. Выведенные аналитические зависимости и численные значения пара-

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

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

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

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

НАУЧНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

/. Статьи в изданиях, рекомендованных для публикации ВАК РФ

1. Дыбрин A.A. Моделирование течения двухфазной смеси в криволинейных каналах газопровода ii Интеллектуальные системы в производстве -Ижевск: Изд-во ИжГТУ им. М.Т. Калашникова - № 2. -2013. - С. 32-43.

2. Дыбрин A.A., Тененев В.А., Абдуллаев Р.В. Методы обеспечения без-гидратного режима работы аппаратов воздушного охлаждения на компрессорных станциях магистральных газопроводов // Нефтепромысловое дело - М.: ОАО «ВНИИОЭНГ». -№1. - 2013. - С. 42-47.

3. Дыбрин A.A., Лялин В.Е. Решение задачи о звукоизоляции цилиндрических оболочек на основе исследования их нормальных импедансов // Оборудование и технологии для нефтегазового комплекса - М.: ОАО «ВНИИОЭНГ». -№4.-2013.-С. 62-67.

4. Дыбрин A.A., Лялин В.Е. Снижение уровня шума при прохождении газа через трубы переменного сечения газопроводов // Защита окружающей среды в нефтегазовом комплексе - М.: ОАО «ВНИИОЭНГ». - №4. - 2013. - С. 20-26.

5. Дыбрин A.A., Лялин В.Е. Экспериментально-расчётный метод определения характеристик акустического поля в трубах энергетических установок ii Защита окружающей среды в нефтегазовом комплексе - М.: ОАО «ВНИИОЭНГ». - №12. - 2012. - С. 20-26.

6. Дыбрин A.A. Динамика участков газопровода с фланцевым соединением ii Интеллектуальные системы в производстве - Ижевск: Издательство ИжГТУ. - №2(20) - 2012 - С. 163-167.

24 6/

II. Материалы международных конференций

7. Дыбрин A.A. Исследование поперечных колебаний однопролетного трубопровода с фланцевым соединением // «Актуальные вопросы современной науки»: Материалы IV международной научной конференции, г. Санкт-Петербург. - Петрозаводск: Петропресс, 2012. - 109 с. - В надзаг.: Науч.-изд.

центр «Открытие», С. 47-52.

8. Дыбрин A.A. Принцип взаимности акустических полей для труб переменного сечения в газопроводах // «Актуальные вопросы современной науки»: Материалы IV международной научной конференции, г. Санкт-Петербург. - Петрозаводск: Петропресс, 2012. - 109 с. - В надзаг.: Науч.-изд. центр «Открытие», С. 52-56.

9. Дыбрин A.A. Решение уравнений течения газа с дисперсными частицами в криволинейных газопроводах // «Актуальные вопросы науки»: Материалы VII Международной научно-практической конференции. - М.: Издательство

«Спутник +», 2012. - С.49-56.

10. Дыбрин A.A. Определение зон гидратообразования в трубах газопровода // «Актуальные вопросы науки»: Материалы VII Международной научно-практической конференции. - М.: Издательство «Спутник +», 2012. - С.42-49.

11. Дыбрин A.A. Расчет коэффициента звукопрохождения в трубах переменного сечения газопроводов // Интеллектуальный потенциал 21 века: ступени познания: Материалы XVI Молодежной международной научно-практической конференции, Новосибирск, 2012, С. 90-93.

12. Дыбрин A.A. Анализ акустических колебаний больших амплитуд в энергетических установках газотранспортных систем // Техника и технологии: роль в развитии современного общества: Материалы Международной научно-практической конференции - Краснодар, 2013. - С. 69-74.

III. Статьи в других изданиях

13. Дыбрин A.A. Звукоизоляция цилиндрической оболочки от внешнего источника шума в ограниченном пространстве // НПЖ Приволжский научный вестник -

Ижевск, 2014-№2(30) - С. 34-39.

14. Дыбрин A.A. Звукоизоляция полуцилиндрическим кожухом при ограниченном источнике // НПЖ Приволжский научный вестник - Ижевск, 2014 -№2(30) - С. 40-46.

В авторской редакции

Подписано в печать 17. 03.14. Усл. печ. л. 1,4. Заказ №91. Тираж 100 экз. Издательство Ижевского государственного технического университета имени М. Т. Калашникова Отпечатано в типографии Издательства ИжГТУ. 426069, Ижевск, Студенческая, 7

Текст работы Дыбрин, Александр Андреевич, диссертация по теме Системный анализ, управление и обработка информации (по отраслям)

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ ФГБОУ ВПО «ИЖЕВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМЕНИ М.Т. КАЛАШНИКОВА»

АНАЛИЗ И МОДЕЛИРОВАНИЕ УРОВНЯ ВИБРАЦИЙ И ШУМА В ТРУБАХ ПЕРЕМЕННОГО СЕЧЕНИЯ ГАЗОПРОВОДОВ

Специальности:

05.13.01 - Системный анализ, управление и обработка информации

(в науке и технике)

01.02.06 - Динамика, прочность машин, приборов и аппаратуры

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель:

заслуженный изобретатель РФ, доктор технических наук, доктор геолого-минералогических наук, профессор Лялин В.Е.

Ижевск 2014

На правах рукописи

04201457664

ДЫБРИН Александр Андреевич

ОГЛАВЛЕНИЕ

Стр.

Введение.............................................................................................................5

Глава 1. Анализ причин возникновении вибраций и шумов в магистральных газопроводах.....................................................................................13

1.1. Состав сооружений магистральных газопроводов..........................13

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

1.3. Методы устранения вибраций трубопроводов................................22

1.4. Гидратообразование в газопроводах................................................25

1.4.1. Условия образования гидратов..................................................26

1.4.2. Определение зон гидратообразовання.................................28

1.4.3. Предупреждение образования гидратных пробок..............30

1.5. Моделирование процессов гидратообразовання при течении влажного газа в газопроводах..................................................................31

1.5.1. Уравнения для расчета осесимметричного течения газожидкостной смеси в дисперсно-кольцевом режиме.....................3 1

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

1.5.3. Определение местоположения начала образования гидратов при двумерной и одномерной зависимости влагосодержа-ния от длины трубы........................................................................39

1.5.4. Решение уравнений течения смеси природного газа и дисперсных частиц в пространственно-криволинейных трубопроводах..........................................................................................41

1.6. Снижения шума в газопроводах.......................................................48

1.6.1. Источники шума в газоперекачивающих агрегатах...........48

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

и газопроводах................................................................................50

Глава 2. Математические модели динамики однопролетного трубопровода с фланцевым соединением...............................................................55

2.1. Задача статики пространственно-криволинейных трубопроводов, заполненных стационарным потоком идеальной несжимаемой жидкости..........................................................................................................55

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

2.3. Динамика фланцевых соединений трубопроводов высокого давления .........................................................................................................65

2.4. Исследование поперечных колебаний однопролетного трубопровода с фланцевым соединением..............................................................74

2.5. Полученные результаты и выводы...................................................82

Глава 3. Методы снижения уровня шума в магистральных газопроводах..84

3.1. Снижение уровня шума при прохождении потока газа через трубопроводы переменного сечения............................................................84

3.1.1. Звукоизоляция трубы переменного сечения.......................85

3.1.2. Соотношение взаимности для труб переменного сечения . 87

3.1.3. Экспоненциальный диффузор (конфузор)..........................90

3.1.4. Конический диффузор (конфузор).......................................93

3.1.5. Параболический диффузор (конфузор)...............................97

3.1.6. Степенные диффузоры (конфузоры)...................................98

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

3.2.1. Уравнение Вебстера для вынужденных колебаний............100

3.2.2. Теорема взаимности для трубы переменного сечения.......102

3.2.3. О коэффициентах прохождения по энергии через диффузор и конфузор................................................................................103

3.3. Акустический импеданс бесконечной цилиндрической оболочки трубы.........................................................................................................106

3.4. Экспериментально-расчётный метод определения характеристик акустического поля..................................................................................115

3.5. Полученные результаты и выводы...................................................129

Глава 4. Методы снижения шумоизлучения трубопроводов вибропог-лощающими и звукоизолирующими конструкциями................................132

4.1. Введение.............................................................................................133

4.2. Звукоизоляция цилиндрической оболочки в ограниченном пространстве от внешнего источника шума..................................................89

4.3. Звукоизоляция полуцилиндрическим кожухом при ограниченном источнике.........................................................................140

4.4. Полученные результаты и выводы...................................................148

Заключение.......................................................................................................151

Литература........................................................................................................154

ВВЕДЕНИЕ

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

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

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

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

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

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

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

Магистральные газопроводы (МГ) относятся к опасным производственным объектам. Часть оборудования, например газоперекачивающие агрегаты (ГПА), являются не только источником потенциальной опасности, но и источником образования вредности. Следствием их работы является высокий уровень шума (90...130 дБ) аэродинамического и механического характера, создаваемый истечением рабочей среды в газовоздушных трактах всасывания и выхлопа, которые, с точки зрения акустики, представляют собой своеобразные волноводы, способные практически беспрепятственно транспортировать акустическую энергию. При создании и модернизации систем шумоглушения требуется проводить экспериментальные исследования, которые вследствие больших габаритов ГПА сложны, дорогостоящи и продолжительны по времени.

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

Степень научной разработанное in темы исследования. Исследования механики пространственно-криволинейных трубопроводов и фланцевых соединений рассматриваются в работах Светлицкого В.А., Власова В.В., Гольденвейзера А.Л., Башты О.Т., Акссльрада Э.Л., Айнбивдера А.Б., Камерштейна А.Г., Гершгейна М.С., Самарина A.A.

Теоретическими и экспериментальными исследованиями снижения уров-

пя шума на компрессорных станциях занимались многие известные ученые: Кравчун П.А., Занченко В.И., Терехов А.Л., Рейнольде А.Ж. Леонтьев В.А., Григорян Ф.И., Юдин Е.А., Бэтчерлор Д., и др. Созданные ими теоретические положения и технические средства позволяют снижать уровень шума. Однако при многообразии конструкций оборудования компрессорных станций эти разработки не обеспечивают необходимую эффективность вследствие невозможности учета всех факторов, влияющих на шумообразованис.

Область исследования. Диссертационная работа выполнена в соответствии с пунктами «6. Методы идентификации систем управления на основе ретроспективной, текущей и экспертной информации», «11. Методы и алгоритмы прогнозирования и оценки эффективности, качества и надежности сложных систем» паспорта специальности 05.13.01 - Системный анализ, управление и обработка информации (в науке и технике), и пунктами «1. Динамика машин, приборов, аппаратуры, систем и комплексов машин и приборов», «9. Математическое моделирование поведения технических объектов и их несущих элементов при статических, динамических, тепловых, коррозионных и других воздействиях» паспорта специальности 01.02.06 - «Динамика, прочность машин, приборов и аппаратуры».

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

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

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

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

Для достижения цели требуется решить следующие задачи:

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

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

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

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