WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

 

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

Прокофьев Андрей Брониславович

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

Специальность 01.02.06 – «Динамика, прочность машин, приборов и аппаратуры»

Автореферат

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

доктора технических наук

Самара 2008 г.

Работа выполнена на кафедре автоматических систем энергетических установок ГОУ ВПО «Самарский государственный аэрокосмический университет»

Научный консультант:                доктор технических наук, профессор

                               Шахматов Евгений Владимирович

Официальные оппоненты:  доктор технических наук, профессор

Балякин Валерий Борисович,

доктор технических наук, профессор

Кулаков Геннадий Алексеевич,

доктор физико-математических наук, профессор

Радченко Владимир Павлович

Ведущая организация: ФГУП «ГНП РКЦ «ЦСКБ-Прогресс»

Защита состоится 28 ноября 2008 г. в 1200 часов на заседании диссертационного совета Д212.215.02 при Государственном образовательном учреждении высшего профессионального образования «Самарский государственный аэрокосмический университет им. акад. С.П. Королева» по адресу: 443086, г. Самара, Московское шоссе, 34.

С диссертацией можно ознакомиться в библиотеке Государственного образовательного учреждения высшего профессионального образования «Самарский государственный аэрокосмический университет им. акад. С.П. Королева».

Автореферат разослан «____» ______________ 2008 г.

Ученый секретарь

диссертационного совета Д212.215.02

доктор технических наук, профессор                                Д.Л. Скуратов

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

Актуальность темы. Современный гидропривод различных технических объектов характеризуется высокой энергонасыщенностью его узлов и элементов. Трубопроводы, являющиеся важным элементом гидропривода, работают в условиях воздействий как статического, так и динамического характера. При этом требования к надежности и работоспособности трубопроводов непрерывно растут. Одним из основных факторов, снижающих работоспособность трубопроводных систем, являются виброакустические нагрузки. К ним относятся вибрационные нагрузки, колебания давления и расхода рабочей жидкости, приводящие к разгерметизации трубопроводных соединений, появлению усталостных трещин. Виброакустические процессы в трубопроводных системах характеризуются широким спектром частот (20 Гц – 5 кГц) и амплитуд колебаний (размах колебаний давления рабочей жидкости – до 20 МПа, виброускорение – до 1000 м/с2). Причем именно колебания давления, амплитуда которых может достигать 50% от среднего значения, являются одной из причин, приводящих к вибрации, усталостным поломкам трубопроводов.

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

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

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

Основные задачи диссертационной работы состоят в следующем:

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

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

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

  1. Разработан метод комплексного исследования виброакустического состояния трубопроводных систем, позволяющий анализировать влияние гасителей колебаний рабочей жидкости как на пульсационную, так и на вибрационную нагруженность систем.
  2. Разработана математическая модель виброакустических процессов в трубопроводных системах сложной пространственной конфигурации, учитывающая одновременное воздействие на систему силового возбуждения со стороны рабочей жидкости и кинематического возбуждения от присоединенных опор и агрегатов и позволяющая оценивать прочностные характеристики таких трубопроводных систем.
  3. Разработан комплексный метод проектирования трубопроводных систем, обеспечивающий требуемые виброакустические и прочностные характеристики в условиях одновременного воздействия на трубопроводы статического нагружения давлением и динамического нагружения пульсациями рабочей жидкости, кинематического возбуждения со стороны опор и присоединенных агрегатов, температурных и внешних силовых нагрузок.
  4. На основе развития теоретического описания процессов виброакустического взаимодействия в гидромеханических системах разработана математическая модель динамики и прочности трубопроводных систем с гасителями колебаний.
  5. Разработаны конечно-элементные модели гасителей колебаний жидкости, позволяющие определять комплекс собственных характеристик гасителей с учетом распределенности параметров всех их элементов, а также их прочностные характеристики в условиях пульсирующего потока рабочей жидкости.
  6. Разработан алгоритм выбора схемы и конструктивных характеристик гасителей колебаний жидкости, мест их установки в систему, обеспечивающий требуемое снижение вибрационных и пульсационных нагрузок при выполнении условий технического задания по гидравлическому сопротивлению, габаритам и массе.
  7. Разработаны методы определения интегральных виброакустических характеристик пневмогидромеханических систем на основе вейвлет-анализа экспериментальных данных.

Практическая ценность и реализация результатов работы.

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

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

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

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

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

Результаты, полученные в работе, использованы:

- в ФГУП «ГНП РКЦ «ЦСКБ-Прогресс» при разработке мероприятий по исключению усталостных разрушений трубопроводов агрегата привода для наземных испытаний рулевых машин ракетоносителя;

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

- на Сызранской ТЭЦ при разработке мероприятий по исключению разрушений маслонасоса турбогенератора;

- в Институте акустики машин при проектировании гасителей колебаний рабочей жидкости;

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

Апробация работы. Основные результаты работы докладывались на Международной научно-технической конференции «Проблемы и перспективы развития двигателестроения» (Самара, 1997, 2001, 2006 гг.), 139-ой конференции Акустического общества Америки (Атланта, США, 2000 г.), IV Всероссийской научно-практической конференции «Решетневские чтения» (Красноярск, 2000 г.), 140-ой конференции Акустического общества Америки (Ньюпорт Бич, США, 2000 г.), Международной научной конференции «XXVIII Гагаринские чтения» (Москва, 2002 г.), III конференции пользователей программного обеспечения CAD-FEM GmbH (Москва, 2003 г.), Международной научно-технической конференции «Современное состояние и перспективы развития гидромашиностроения в XXI веке» (Санкт-Петербург, 2003 г.), IV конференции пользователей программного обеспечения CAD-FEM GmbH (Москва, 2004 г.), VII Международном симпозиуме «Transport Noise and Vibration» (Санкт-Петербург, 2004 г.), IV Школе-семинаре молодых ученых и специалистов «Проблемы тепломассообмена и гидродинамики в энергомашиностроении» (Казань, 2004 г.), V конференции пользователей программного обеспечения CAD-FEM GmbH (Москва, 2005 г.), XIV Международном научно-техническом семинаре «Современные технологии в задачах управления, автоматики и обработки информации» (Алушта, Украина, 2005 г.), VI Международной научно-практической конференции «Людина i Космос» (Днепропетровск, Украина, 2005 г.), III Международной научно-технической конференции «Гидравлические машины, гидроприводы и гидропневмоавтоматика» (Санкт-Петербург, 2005 г.), 5-ой международной научно-технической конференции «Обеспечение и повышение качества машин на этапах их жизненного цикла» (Брянск, 2005 г.), XIII международном конгрессе «Sound and Vibration (ICSV13)» (Вена, Австрия, 2006 г.), VII Международной научно-практической конференции «Людина i Космос» (Днепропетровск, Украина, 2006 г.), 5-ой международной конференции «Авиация и космонавтика – 2006» (Москва, 2006 г.), X международной научной конференции «Решетневские Чтения» (Красноярск, 2006 г.), VIII Международной научно-практической конференции «Людина i Космос» (Днепропетровск, Украина, 2007 г.), Международной конференции «Power Transmission and Motion Control 2007» (Бас, Великобритания, 2007 г.).

Публикации. По теме диссертации опубликовано 72 научные работы, в том числе 2 монографии, 1 патент на изобретение, 1 патент на полезную модель.

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

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

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

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

Проведен анализ методов снижения виброакустических нагрузок. Показано, что одним из методов снижения интенсивности динамических процессов является подавление колебаний в источнике. Этому направлению посвящены работы Башты Т.М., Орлова Ю.М., Шахматова Е.В., Юдина Е.М., Шорина В.П., Задонцева В.А., Овсянникова Б.В., Пилипенко В.В., Крючкова А.Н. и др. Однако устранение колебаний в источнике часто связано с необходимостью значительных изменений в конструкции отдельных агрегатов, а в некоторых случаях сам характер организации рабочего процесса предопределяет генерирование колебаний. Кроме того, эффективность мероприятий, связанных с уменьшением генерируемых источником колебаний, зачастую имеет противоречивый характер. Например, увеличение зазора между рабочим колесом и статором центробежного насоса способствует уменьшению амплитуды колебаний давления, но приводит к снижению производительности и КПД насоса. Развитие рассматриваемого пути снижения виброакустической нагруженности сдерживает отсутствие адекватных методов моделирования и научно-обоснованных методов проектирования устройств, рабочий процесс которых сопровождается генерированием колебаний рабочей среды. Оценивая возможности устранения колебаний в источнике, следует отметить, что задача далека от разрешения.

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

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

Проведен анализ методов математического моделирования виброакустических характеристик элементов гидромеханических систем. Вопросы математического описания виброакустических процессов рассмотрены в работах Колесникова К.С., Светлицкого В.А., Попова Д.Н., Шорина В.П., Шахматова Е.В., Крючкова А.Н., Гликмана Б.Ф., Кондрашова Н.С. Однако до настоящего времени задачи моделирования процессов виброакустического взаимодействия гидравлической и механической подсистем решены не полностью. Связано это с тем, что традиционно различные научные школы делали акцент на детальном математическом описании какой-либо одной из подсистем, принимая достаточно серьезные упрощения при моделировании другой и упуская некоторые особенности их взаимодействия.

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

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

Проведен анализ динамических характеристик гасителей колебаний давления. В работах Шорина В.П., Головина А.Н. показано, что акустические свойства гасителя однозначно описываются матрицей передачи:

определяющей взаимосвязь комплексных амплитуд давления Pi и расхода Qi на входе (i=1) и выходе (i=2) устройства и имеющей комплексные частотнозависимые коэффициенты. Коэффициенты матрицы передачи определяют комплекс собственных характеристик гасителя: волновые сопротивления со стороны входа Zc1 и выхода Zc2 и коэффициент собственного затухания Kc. Это – внутренние параметры. Внешними параметрами гасителя, характеризующими его эффективность по снижению пульсаций давления являются коэффициент вносимого затухания Kвн и коэффициент бегущей волны Kб. Они определяются как комплексом собственных характеристик гасителя, так и параметрами присоединенной системы.

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

Рис. 4. Конструктивные решения элементов гасителей колебаний.

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

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

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

а)

б)

в)

Рис. 5.  Варианты конструктивных соединений фланцев с корпусом гасителя колебаний жидкости: а) сварное; б) болтовое; в) резьбовое.

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

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

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

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

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

Разработана методика конечно-элементного моделирования виброакустических характеристик трубопроводных систем с пульсирующим потоком рабочей жидкости в программном комплексе ANSYS. Основные допущения к модели: отсутствие вязкости жидкости, малый средний расход жидкости, постоянство среднего уровня давления и плотности, гармонический закон изменения давления во времени, адиабатичность волновых процессов. Как показали расчетные и экспериментальные исследования, неучет вязких свойств жидкости не влияет на точность определения резонансных областей на АЧХ трубопровода, а также на значения виброакустических параметров в межрезонансной области. При этом наблюдается некоторое завышение расчетных значений вибрации и пульсаций давления на частотах, близких к собственным. Обоснован выбор типов конечных элементов, позволяющих проводить процесс моделирования виброакустических характеристик в трубопроводных системах с учетом принятых допущений. Даны рекомендации по выполнению процедуры генерации конечно-элементной сетки. Показано, что наибольшая точность и наименьшая трудоемкость расчетов обеспечивается при использовании регулярной сетки конечных элементов (рис. 7). Разработан способ задания граничных условий механической и гидравлической подсистем. С помощью разработанной методики численного моделирования проведено исследование виброакустических характеристик элементов трубопроводных систем при их силовом нагружении пульсациями давления рабочей жидкости и кинематическом возбуждении. Проведен анализ сходимости результатов моделирования по разработанной и известным методикам. Показана высокая степень сходимости в области низких частот (до частот в 1.5 – 2 раза ниже частот первого резонанса по гидравлической или механической подсистеме трубопровода – рис. 8). В области высоких частот адекватность разработанной конечно-элементной модели доказана экспериментально (рис. 9).

К сложностям использования такого метода моделирования следует отнести его высокую трудоёмкость. Для адекватного описания поведения трубопровода необходимо использовать регулярную сетку конечных элементов.  При этом в поперечном сечении (см. рис. 7) видно 32 конечных элемента, описывающих стенку трубопровода, и 48 конечных элементов, моделирующих рабочую жидкость. Таким образом, полная конечно-элементная модель трубопровода содержит сотни тысяч конечных элементов и является довольно громоздкой при проведении расчётов. Расчёт амплитудной частотной характеристики вибропараметров трубопровода на персональном компьютере Pentium IV занимает несколько часов. Большую сложность представляет также создание сетки конечных элементов в местах разветвления трубопроводов. В ряде случаев это приводит к ситуациям, когда размерность системы уравнений решаемой задачи превышает возможности программного комплекса и используемых персональных компьютеров. Особенно актуально это для сложных и разветвлённых трубопроводных систем.

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

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

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

,

(1)

,  ,  ,  , ,  ,

где - вектор перемещения точек осевой линии трубопровода; - вектор внутренних сил; w – относительная скорость рабочей жидкости; p – давление жидкости; - вектор сил взаимодействия между трубопроводом и рабочей жидкостью; τ - время; ε - координата, отсчитываемая вдоль осевой линии; n – безразмерный коэффициент, равный отношению погонной массы трубопровода в каком-либо сечении к погонной массе трубопровода в начальном сечении (для однородного трубопровода n=1); - единичные векторы, направленные по касательной (i=1), нормали (i=2) и бинормали (i=3) к осевой линии трубопровода.

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

Система (1) в проекциях на координатные оси представляет собой систему 22 уравнений. Однако при расчете используется только 20 уравнений. В связи с принятой идеализацией движения рабочей жидкости (гипотеза одномерности движения) проекции уравнения движения жидкости на оси нормальной и бинормальной кривизны не учитываются. Разработана методика задания краевых условий. Для решения рассматриваемой задачи необходимо 14 краевых условий: 12 – для механической подсистемы (по 6 на каждом конце трубопровода) и два – для гидравлической (акустической) подсистемы.

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

,

,                                                        (2)

,  , ,

,

где - символы Леви – Чивита, w0 и Δw – постоянная и пульсационная составляющая скорости жидкости соответственно; p и Δp – постоянная и переменная составляющая давления соответственно; , , и т.д., где , , - динамические составляющие соответствующих векторов; векторы , , имеют компоненты в базисе , равные статическим значениям; φ1, φ2, φ3 – углы поворота подвижных координатных осей, τтрн – нестационарное касательное напряжение на стенке трубы. При выводе системы уравнений (2) учитывалось, что пульсационные составляющие давления и скорости рабочей жидкости могут быть одного порядка или даже больше, чем соответствующие статические составляющие. Особенность системы (2) заключается в зависимости ее соотношений от членов и , которые определяются из условий равновесия трубопровода, заполненного стационарным потоком жидкости, с действующими на него статическими силами , , . Для жестких трубопроводов приближенно можно считать, что статические силы и стационарный поток практически не изменяют форму осевой линии, т.е. вектор кривизны для нагруженного трубопровода равен вектору для трубопровода в естественном состоянии.

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

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

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

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

(3)

Распределение параметров пульсаций давления и скорости в пространстве и во времени при этом рассчитывается отдельно по известным формулам. Разработана конечно-разностная схема решения системы (3). Показано, что разработанная схема является условно устойчивой. Разработаны рекомендации по выбору шагов дискретизации в пространстве и во времени, обеспечивающие устойчивость схемы. Результаты расчета временных реализаций безразмерного нормального вибросмещения для рассмотренного модельного случая трубопровода, вибрирующего под действием пульсирующего потока жидкости из положения покоя, представлены на рис. 12. Рис. 12,а соответствует системе без трения, в которой наряду с установившимися вынужденными колебаниями реализуются незатухающие свободные колебания. Введение в систему трения приводит к быстрому затуханию свободных колебаний (рис. 12,б). Частота вынужденных колебаний на временных реализациях вибросмещения соответствует частоте пульсаций жидкости. На рис. 13 представлена амплитудная частотная характеристика виброускорения трубопровода, колеблющегося под действием пульсирующего потока жидкости с амплитудой 105 Па во входном сечении, полученная в результате расчетов по разработанной модели. Для проверки адекватности расчетов  по  разработанной  математической  модели  было  проведено моделирование рассматриваемого трубопровода в программном комплексе ANSYS на базе указанной выше методики. Результаты расчетов в программном комплексе ANSYS и по разработанной математической модели совпали полностью. Однако вычислительные затраты при моделировании виброакустических характеристик трубопровода в программном комплексе ANSYS были на три порядка выше, чем по разработанной математической модели, что позволяет судить об ее высокой эффективности.

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

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

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

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

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

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

Рис. 15. Частотные зависимости комплекса собственных характеристик гасителя:

1 – аналитическая модель в сосредоточенных параметрах; 2 – аналитическая модель, учитывающая распределенность параметров центрального канала; 3 – численная модель.

Сравнение результатов расчета по аналитическим моделям в сосредоточенных параметрах и с учетом распределенности проточного канала показывает их некоторое количественное различие. Однако качественно почти все характеристики этих двух моделей схожи друг с другом. Если же сравнивать результаты расчетов по аналитической и численной моделям, то можно сделать вывод, что при <2 различие результатов по всем моделям незначительное. Однако при >2 появляется их качественное расхождение. Такое различие результатов объясняется отсутствием учета в аналитических моделях распределенности параметров расширительной полости. Проведенные расчеты позволяют сделать вывод, что в области низких частот (или ) наиболее целесообразно использование аналитической модели, позволяющей достаточно легко  анализировать  зависимости  собственных  характеристик  гасителей  от свойств входящих в их структуру элементов, реализовать процедуру оптимизации структуры. Однако в высокочастотной области трудности учета распределенности параметров в аналитической модели приводят к некоторому завышению расчетной эффективности гасителя колебаний. Поэтому для расчета собственных характеристик гасителя и выбора оптимальной его конструкции здесь более целесообразно применение конечно-элементной модели.

Разработана методика конечно-элементного моделирования динамических характеристик емкостей с податливыми стенками. Целью расчетов является определение коэффициента k, представляющего собой отношение приведенного и геометрического объемов полости. На рис. 16 представлены расчетные частотные характеристики этого коэффициента для ряда сферических емкостей с внутренним радиусом Rвн и толщиной стенки δ. Анализ результатов моделирования показывает, что во всех случаях k>1, т.е. приведенный объем полости больше, чем геометрический. При этом в исследованном диапазоне частот до 400 Гц для различных емкостей значения этого коэффициента составляют от 1.1 до 1.25. Большие значения коэффициента k соответствуют емкостям с большими величинами отношения и меньшими толщинами стенки. Коэффициент k является частотнозависимым, однако эта зависимость выражена слабо. Разработанная методика моделирования динамических свойств емкостей и полученные на основе ее использования результаты позволяют сделать вывод, что аналитическое соотношение, описывающее динамические свойства емкостей:

,

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

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

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

               (4)

где -  безразмерная частота колебаний, ; ; а и b – размеры поперечного сечения; и -  параметры, определяющие длину волны колебаний в направлениях x и y соответственно; ; h – толщина стенки канала. Система (4) может быть решена методом простых итераций. В качестве начального приближения для первой итерации целесообразно брать результаты, полученные для случая с абсолютно жесткими стенками (рис. 17):

,

где с – скорость звука в среде; с0 -  скорость звука в материале ограничивающей поверхности; n и m – число узловых линий в направлениях y и x соответственно; . Решение быстро сходится. Получены аналитические выражения для собственных форм колебаний . Полученные результаты обобщены для трехмерного случая.

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

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

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

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

1) соотношения (3), описывающие виброакустические характеристики входного и выходного участков;

2) соотношения, описывающие динамику механической подсистемы гасителя:

, ,  ,

,,

где .

3) условия стыковки рассматриваемых участков трубопровода с гасителем колебаний:

,

где ;

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

, , ;

; ;;

; ; ; ; ; ; ; ; ; ; ; ,

где - постоянная распространения волнового процесса; Aг, Bг, Сг, Dг, - частотнозависимые коэффициенты матрицы передачи  гасителя; Zв – волновое сопротивление трубопровода; n – число энергетически значимых гармоник в спектре пульсационного процесса.

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

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

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

- критерии эффективности снижения вибрационных нагрузок -  коэффициент изменения вибрации;

- критерии повышения прочностных характеристик - коэффициент увеличения прочности.

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

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

;

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

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

В качестве критерия, оценивающего изменение прочностных характеристик трубопроводной системы после установки в нее гасителя колебаний, предлагается коэффициент увеличения прочности Kуп,  представляющий собой отношение минимального запаса  прочности n2min  в каком-либо сечении после установки гасителя к минимальному запасу прочности в каком-либо (возможно другом) сечении до установки n1min. Минимальный запас прочности трубопроводной системы после установки гасителя определяется как наименьшее значение запасов прочности собственно для трубопровода и для самого гасителя колебаний. Следует отметить, что введение гасителя, снижающего уровень и пульсационных, и вибрационных нагрузок, может сопровождаться значениями Kуп<1. Это может быть обусловлено либо низкой прочностью элементов самого гасителя, либо внесением дополнительных концентратов напряжений в местах соединений гасителя и трубопроводов. Разработана методика расчета коэффициента увеличения прочности, основанная на использовании математических моделей и расчетных соотношений глав 2 и 3.

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

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

Исходными данными являются свойства материала трубопровода и рабочей жидкости, геометрия трубопроводной системы, динамические характеристики гидромеханической системы и ограничения на параметры проектируемого устройства. Алгоритм (см. рис. 19) состоит из трех основных блоков: блока предварительного расчета виброакустических характеристик системы, в результате работы которого определяются потребные уровни снижения виброакустической нагруженности; блока выбора места установки, схемы, расчета оптимального соотношения акустических характеристик, конструктивных параметров гасителя; блока проверочных расчетов виброакустических и прочностных характеристик.

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

Разработаны методики анализа данных исследований быстропеременных, в том числе неустановившихся, процессов в элементах гидромеханических систем, основанные на применении цифрового (дискретного) вейвлет-анализа. В частности, разработана экспериментальная методика определения собственных частот гидромеханической системы на базе вейвлет-анализа данных теста с вариацией скорости. По сравнению с традиционными процедурами оконного преобразования Фурье она позволяет решать задачу определения собственных частот с большей точностью в частотной области и достоверностью, при высокой степени автоматизации вычислительных алгоритмов и незначительном увеличении вычислительных затрат. Расчетные соотношения разработанной методики построены на базе использования вейвлета Морлета. Однако после некоторой доработки методику можно адаптировать и к другим типам вейвлетов. Разработаны рекомендации к выбору параметров вейвлет-преобразования. Методика предусматривает выполнение следующих вычислительных процедур: расчет коэффициентов вейвлет-преобразования; выделение гребней на поверхности амплитуд вейвлет-коэффициентов; построение скелетных кривых; определение зависимостей амплитуд вейвлет-коэффициентов на выделенных скелетных кривых от мгновенного значения соответствующей частоты колебательной компоненты; выделение из этих зависимостей собственных частот. Зависимость амплитуды вейвлет-коэффициентов на одной из выделенных скелетных кривых от мгновенного значения частоты пульсаций давления в топливной системе авиационного двигателя представлена на рис. 20. На рисунке наложено две кривые, одна из которых соответствует разгону, другая – снижению частоты вращения привода насоса. Анализ данной зависимости позволяет выделить собственные частоты – 470, 873 и 1680 Гц.

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

Рис. 21. Вейвлет-преобразование собственных затухающих колебаний трубопровода: а) проекция поверхности амплитуды вейвлет-коэффициентов;
б) сечение поверхности амплитуды вейвлет-коэффициентов плоскостью ; 1 - аппроксимация сечения при b>b0.

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

На рис. 22 представлен график полуэмпирической зависимости модуля объемной упругости металлического рукава от частоты колебаний, экспериментальные данные, полученные лично автором (♦) и другими исследователями (•). На рис. 23,а изображена полученная экспериментально амплитудная частотная характеристика металлического рукава, нагруженного на акустически закрытый конец. Максимальные зарегистрированные величины относительных амплитуд на резонансах не превышают 3…4. В то же время их максимальные значения, полученные автором  при  исследовании  АЧХ  металлических  трубопроводов  аналогичной длины, на резонансах доходят до 20…30. Это свидетельствует о существенном демпфировании колебаний рабочей жидкости по длине металлического рукава, которое может быть объяснено потерями энергии на трение в оплетке. Получены полуэмпирические зависимости для амплитудо-фазовых частотных характеристик пульсаций жидкости. В частности, для случая нагрузки на акустически закрытый конец можно записать следующие выражения для амплитудной и фазовой частотных характеристик:

, ,

где ε - коэффициент фазы; l – длина металлического рукава.

Для рукава, амплитудная и фазовая частотные характеристики которого представлены на рис. 23 для и k получены следующие эмпирические формулы:

,  ,

где δ - коэффициент затухания для металлического трубопровода таких же размеров; f – частота колебаний.

а)

б)

Рис. 23. Амплитудная (а) и фазовая (б) частотные характеристики тупикового металлического рукава:

- эксперимент;

- расчет.

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

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

Проведены экспериментальные исследования динамических процессов в гидросистеме пресса Erfurt. При работе пресса происходило разрушение трубопроводов сливной магистрали (рис. 26), приводящее к сбою работы конвейера сборки автомобилей и потере нескольких тонн гидравлического масла. Было установлено, что причиной поломок являются гидроударные процессы в трубопроводной системе и вызываемая ими вибрация трубопроводов (рис. 27). На основе разработанного метода анализа виброакустических характеристик трубопроводных систем был предложен комплекс мероприятий по уменьшению динамической нагруженности - введение в систему гидроаккумуляторов, снижающих интенсивность гидроударных процессов, и установка дополнительных опор трубопроводов, исключающих их резонансы. Реализация указанного комплекса мероприятий позволила снизить пиковые значения вибрации трубопроводов в 2.5 - 3.5 раз (рис. 28).

Рис. 27. Временные реализации давления и виброускорения. Коллектор сливной магистрали.

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

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

а)

б)

Рис. 29. Амплитуды гармоник пульсаций давления (а) и виброскорости (б) до и после установки гасителя колебаний:

- до установки; - после установки гасителя.

Для эффективного снижения интенсивности всех колебательных компонент в спектре давления на выходе из насоса был установлен гаситель колебаний жидкости. Его применение снизило среднеквадратичное значение пульсаций более чем в 30 раз (рис. 29). Замеры вибрации на трубопроводе в месте его разрушения показали снижение виброскорости более чем в 6 раз. Введение в гидросистему гасителя колебаний позволило обеспечить работоспособность и заданный ресурс испытательного стенда.

РЕЗУЛЬТАТЫ  И  ВЫВОДЫ

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

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

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

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

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

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

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

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

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

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

10. Разработанные методы анализа и моделирования, средства коррекции динамических характеристик гидромеханических систем и гасители колебаний жидкости внедрены в ФГУП ГНП РКЦ «ЦСКБ-Прогрес», в ОАО «Самарское конструкторское бюро машиностроения», на Сызранской ТЭЦ. Они позволили снизить амплитуды пульсаций давления рабочей среды до 30 раз, амплитуды вибрации – до 6 раз.

Основные публикации по теме диссертации

Монографии:

1. Гимадиев А.Г., Крючков А.Н., Леньшин В.В., Прокофьев А.Б., Шахматов Е.В., Шестаков Г.В., Шорин В.П. Снижение виброакустических нагрузок в гидромеханических системах. – Самара, СГАУ, 1998. – 270 с.

2. Иголкин А.А., Крючков А.Н., Макарьянц Г.М., Прокофьев А.Б., Прохоров С.П., Шахматов Е.В., Шорин В.П. Снижение колебаний и шума в пневмогидромеханических системах. – Самара: Изд-во СНЦ РАН, 2005. – 314 с.

Статьи в изданиях, рекомендуемых ВАК РФ:

3. Прокофьев А.Б. Расчет собственных частот и форм колебаний трубопроводов с помощью программного комплекса // Изв. СНЦ РАН. – г. Самара, СНЦ РАН, 1999, №2. – С. 335-342.

4. Прокофьев А.Б., Шестаков Г.В. Оптимизация профиля инерционного элемента гасителей колебаний // Вестник СГАУ. Серия: Проблемы и перспективы развития двигателестроения. Вып. 3, Ч.2. – Самара, СГАУ, 1999. – С. 60-68.

5. Прокофьев А.Б., Шахматов Е.В. Виброакустическая модель прямолинейного неоднородного трубопровода при его силовом возбуждении пульсациями рабочей жидкости // Изв. СНЦ РАН. – Самара, 2000, Т.2, №1. – С.135-140.

6. Прокофьев А.Б., Шестаков Г.В. Оценка погрешности расчета коэффициента вносимого затухания при неучете распределенности параметров гасителя колебаний // Вестник СГАУ. Серия: Проблемы и перспективы развития двигателестроения. – Самара, СГАУ, 2000, Вып. 4, Ч. 2. – С. 86-90.

7. Журавлев О.А., Комаров С.Ю., Попов К.Н., Прокофьев А.Б. Разработка автоматизированного метода исследования вибрационных характеристик энергоустановок // Компьютерная оптика. – Самара, ИСОИ РАН, №21, 2001. – С. 7-11.

8. Комаров С.Ю., Прокофьев А.Б., Шапошников Ю.Н., Щеглов Ю.Д. Исследование колебаний трубопровода методом цифровой спекл-интерферометрии // Известия СНЦ РАН. – Самара, СНЦ РАН, 2002, Т.4, №1(7) – С. 87-90.

9. Иголкин А.А., Крючков А.Н., Прокофьев А.Б., Шахматов Е.В. Исследование влияния подачи воздуха во всасывающую магистраль центробежного насоса на его виброакустические характеристики // Вестник СГАУ. – Самара, СГАУ, 2002, №1. – С.78-83.

10. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Моделирование виброакустических характеристик трубопровода с использованием метода конечных элементов // Изв. СНЦ РАН. – Самара, СНЦ РАН, 2002, Т.4, №2(8) – С. 327-323.

11. Леньшин В.В., Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Выбор конструктивно-режимных критериев, наибольшим образом влияющих на работоспособность насосов ГТД // Вестник СГАУ. Спец. выпуск. Ч.2. – Самара, СГАУ, 2003. – С. 66-75.

12. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Анализ методик расчета виброакустических характеристик криволинейного трубопровода при его силовом возбуждении // Вестник СГАУ. Спец. выпуск. Ч.2. – Самара, СГАУ, 2003. – С. 342-345.

13. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Влияние динамических характеристик присоединенных гидравлических цепей на вибрационные характеристики трубопроводов // Вестник СГАУ. – Самара, СГАУ, 2004, №1. – С.96-101.

14. Прокофьев А.Б., Шахматов Е.В. Использование вейвлет-анализа свободных колебаний трубопроводной системы при оценке характеристик демпфирования// Изв. СНЦ РАН. – Самара, СНЦ РАН, 2004, Т.6, №2(12) – С. 354-363.

15. Прокофьев А.Б. Определение собственных частот системы на базе вейвлет-анализа данных теста с вариацией скорости // Изв. СНЦ РАН. – Самара, СНЦ РАН, 2006, Т.7, №2 – С. 405-414.

16. Миронова Т.Б., Прокофьев А.Б., Шахматов Е.В. Возбуждение пульсаций давления в рабочей жидкости при вибрации трубопровода // Вестник СГАУ. – Самара: СГАУ, №2(10), Ч.2, 2006. – С. 161-164.

17. Макарьянц Г.М., Прокофьев А.Б. Разработка принципиальной схемы экспериментального комплекса исследования виброакустических характеристик амортизаторов транспортных средств // Изв. СНЦ РАН. Спец. выпуск «Безопасность. Технологии. Управление». Т.2. – Самара: Изд-во СНЦ РАН, 2007. – С. 135-138.

18. Астафьев В.И., Прокофьев А.Б., Шахматов Е.В. Анализ собственных частот и форм колебаний жидкости в прямоугольной области в двухмерной задаче // Изв. СНЦ РАН. – Самара, СНЦ РАН, 2007, Т.9, №3 – С. 657-663.

Статьи и материалы конференций:

19. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Некоторые результаты экспериментальных исследований виброакустических характеристик элементов гидромеханических систем // Научно-технический сборник Ракетно-космическая техника. Расчет, проектирование, конструирование и испытания космических систем. Серия XII. Вып. 1. – Самара, ВКБ РКК Энергия, 1999. – С. 191-201.

20. Kruchkov A.N., Lenchine V.V., Prokofiev A.B., Shakhmatov E.V. The investigation of vibroacoustical activity of the combined fuel pump // The Journal of the Acoustical Society of America. 139th Meeting Acoustical Soc. of America. – Vol. 107, №5, Pt2 – Atlanta, USA, 2000 – P.2877.

21. Lenchine V.V., Prokofiev A.B., Shakhmatov E.V. Diagnostics of vibration excitation mechanisms for complex hydromechanical systems // The Journal of the Acoustical Society of America. 140th Meeting Acoustical Soc. of America. – Vol. 108, №5 – Newport Beach, USA, 2000 – P.2500.

22. Прокофьев А.Б., Шахматов Е.В. Виброакустическая модель прямолинейного участка трубопроводной системы с гасителем колебаний в условиях силового возбуждения пульсациями рабочей жидкости // Научно-технический сборник Ракетно-космическая техника. Расчет, проектирование, конструирование и испытания космических систем. Серия XII. Вып. 1. – Самара, ВКБ РКК Энергия, 2000. – С. 120-131.

23. Прокофьев А.Б. Моделирование акустических процессов в трубопроводных системах с использованием программного комплекса ANSYS // Сборник докладов международной научно-технической конференции, посвященной памяти генерального конструктора аэрокосмической техники Н.Д. Кузнецова. – Самара, СГАУ, 2001, Ч.2. – С. 89-96.

24. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Использование численных методов при моделировании собственных колебаний трубопроводных систем // Сборник трудов третьей конференции пользователей программного обеспечения CAD-FEM GmbH. – М.: Полигон-пресс, 2003. – С. 401-408.

25. Иголкин А.А., Крючков А.Н., Прокофьев А.Б. Исследование влияния особенностей пусковых операций на работоспособность и виброакустические характеристики насосного агрегата турбогенератора // Труды международной научно-технической конференции «Современное состояние и перспективы развития гидромашиностроения в XXI веке» – СПбг: СПбГПУ, 2003. – С. 238-246.

26. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В., Шестаков Г.В. Исследование виброакустических характеристик трубопровода при его силовом нагружении с использованием программного комплекса ANSYS // Сборник трудов четвертой конференции пользователей программного обеспечения CAD-FEM GmbH. – М.: Полигон-пресс, 2004. – С. 280-287.

27. Gasparov M.S., Igolkin A.A., Kruchkov A.N., Nazarov O.V., Prokofiev A.B., Shakhmatov E.V. Application of soundproof structures on the basis of Z-gofer panels to reduce transport noise. Proceedings of the 7-th International Symposium "Transport Noise and Vibration", St.Petersburg, 2004, s2-4.

28. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Исследование пульсаций рабочей жидкости в резиновом шланге с металлическими оплетками // Проблемы тепломассообмена и гидродинамики в энергомашиностроении: Труды IV Школы-семинара молодых ученых и специалистов под руководством академика РАН В.Е. Алемасова. – Казань, КГУ, 2004. – С. 480-489.

29. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Численное исследование динамических характеристик участка трубопроводной системы с четвертьволновым резонатором // Сборник трудов пятой конференции пользователей программного обеспечения CAD-FEM GmbH. – М.: Полигон-пресс, 2005. – С. 411-416.

30. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Обеспечение работоспособности трубопроводной обвязки стенда для наземных испытаний рулевых машин // Труды III Международной научно-технической конференции «Гидравлические машины, гидроприводы и гидропневмоавтоматика». – СПб.: Изд-во Политехн. ун-та, 2005. – С. 268-272.

31. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Особенности применения четвертьволнового резонатора в качестве средства снижения виброакустических нагрузок в трубопроводных системах // Обеспечение и повышение качества машин на этапах их жизненного цикла. Материалы 5-ой междунар. науч.-техн. конф. – Брянск: БГТУ, 2005. – С. 19-21.

32. Миронова Т.Б., Прокофьев А.Б. Использование вейвлет-анализа при экспериментальном исследовании собственных частот колебаний гидравлической системы //VIII Мiжнародна молодiжна науково-практична конференцiя «Людина i Космос»: Збiрник тез. – Днiпропетровськ, Украина, НЦАОМУ, 2006. – С. 18.

33. Prokofiev A., Lenchine V., Shakhmatov E. Calculation of logarithmic decrement by morlet wavelet of a decay curve // CD-ROM Proceedings of the Thirtheenth International Congress on Sound and Vibration (ICSV13), July 2-6, 2006, Vienna, Austria, Eds.: Eberhardsteiner, J.; Mang, H.A.; Waubke, H., Publisher: Vienna University of Technology, Austria, ISBN: 3-9501554-5-7.

34. Миронова Т.Б., Прокофьев А.Б Использование вейвлет-преобразования для определения характеристик демпфирования элементов гидромеханической системы // Тезисы докладов X международной научной конференции «Решетневские Чтения», Красноярск, 2006. – c.22-23.

35. Kruchkov A., Prokofiev A., Sverbilov V., Shakhmatov E. Vibroacoustic load reduction in hydro mechanical systems by use of flow oscillation dampers // Power Transmission and Motion Control 2007. – Bath, UK: Hadleys Ltd, 2007. – 351-365 pp.

36. Миронова Т.Б., Прокофьев А.Б. Математическая модель динамики криволинейного трубопровода с пульсирующим потоком рабочей жидкости //VIII Мiжнародна молодiжна науково-практична конференцiя «Людина i Космос»: Збiрник тез. – Днiпропетровськ, Украина, НЦАОМУ, 2007. – С. 15.

37. Макарьянц Г.М., Прокофьев А.Б. Математическая модель динамики системы с дозирующим предохранительным клапаном //VIII Мiжнародна молодiжна науково-практична конференцiя «Людина i Космос»: Збiрник тез. – Днiпропетровськ, Украина, НЦАОМУ, 2007. – С. 19.

Патенты:

38. Макарьянц Г.М., Прокофьев А.Б. Прибор для анализа результатов испытаний амортизаторов в подвеске транспортного средства. Патент на полезную модель №59247. Приоритет полезной модели 21.06.2006 г.

39. Макарьянц Г.М., Прокофьев А.Б., Шахматов Е.В. Двухчастотный настраиваемый гаситель колебаний давления. Патент на изобретение №2293246. Приоритет изобретения 20.04.2005 г.






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

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