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

МЕХАНИКА И МАШИНОСТРОЕНИЕ

УДК 621.64

МОДЕЛИРОВАНИЕ ВИБРОАКУСТИЧЕСКИХ ХАРАКТЕРИСТИК ТРУБОПРОВОДА С ИСПОЛЬЗОВАНИЕМ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ

© 2002 Г.М. Макарьянц1, А.Б. Прокофьев1, Е.В. Шахматов2

1 Самарский государственный аэрокосмический университет 2 Институт акустики машин

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

Вопрос изучения виброакустических характеристик трубопроводов заслуживает серьёзного внимания, так как он непосредственно связан с обспечением надёжности в условиях периодически действующих сил. На данный момент существует ряд теоретических работ посвящённых рассматриваемой тематике [1-6]. Одной из первых появилась работа [3]. В ней рассматривалась задача определения амплитуды параметрических колебаний прямолинейного участка трубопровода на податливых опорах с линейными упругими характеристиками с учётом демпфирования. Н.С. Кондрашовым было представлено аналитическое решение дифференциального уравнения, описывающего процесс установившихся параметрических колебаний. Полученное решение было ограничено случаем основного параметрического резонанса (юеын = 2асоб). Кроме того, в указанной математической модели пренебрегалось распределением параметров пульсаций давления и скорости рабочей жидкости по длине исследуемого трубопровода. Последнее допущение существенно ограничивало область применения модели и было корректно лишь для низких частот, когда длина исследуемого трубопровода не превышала 1/8 длины волны акустических колебаний в рабочей жидкости, а механическая подсистема находилась в дорезонансной области. В работе [2] решалась задача определения вибропараметров колебаний криволинейного трубопровода, вызван-

ных пульсациями давления протекающей жидкости. В частности были получены результаты по определению возбуждающей силы. Однако, полученная методика справедлива только для низкочастотных колебаний, так как в указанной модели не учтено изменение параметров пульсаций давления и скорости жидкости по длине трубопровода. Кроме того не учитывалось влияние присоединённой гидравлической цепи, вида опор и некоторых геометрических параметров системы. В работе [6] рассмотрены методы снижения вибрации трубопроводов. Была разработана математическая модель, позволяющая производить расчёт влияния пульсаций давления на динамические характеристики трубопровода. За основу были приняты уравнения Н.Е. Жуковского, которые в дальнейшем были приведены к виду Матье-Хилла с периодическими коэффициентами. Однако, полученные соотношения позволили произвести только расчёт зон динамической неустойчивости. Таким образом, математическая модель работы [6] позволяет оценить динамику системы только на резонансных режимах. В работе [1] рассмотрены вопросы динамики и устойчивости трубопровода и его элементов. При этом учитывается большое число характеристик исследуемых физических объектов и внешних воздействий, таких как старение материалов, воздействие потока жидкости, вязкоупругого основания и сжимающего (растягивающего) продольного уси-

лия. Однако, полученная автором аналитическая модель не позволяла производить расчёт параметров вибрации при вынужденных колебаниях от пульсаций рабочей жидкости, что представляет большой практический интерес. Сам автор отмечает, что решение подобных задач возможно только на базе численных методов расчёта. Развитием исследований Кондрашова Н.С. является работа [4], где представлено решение ранее полученного уравнения [3] для случая вынужденных колебаний трубопровода при возбуждении пульсациями рабочей среды произвольной частоты. При этом в учёт принималась и распределённость параметров пульсаций. Однако данное решение было справедливо только для однородного прямолинейного трубопровода. В работе [5] решена задача о поперечных колебаниях неоднородного прямолинейного трубопровода в условиях пульсаций давления протекающей жидкости. Предложена методика интегрирования полученной системы дифференциальных уравнений в частных производных с помощью приближённого метода Бубнова-Галеркина. Однако описанная модель была справедлива только для прямолинейного участка трубопровода. Кроме того, учёт сил трения производился с помощью подхода, изложенного в [3], который требует развития и уточнения.

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

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

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

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

Геометрические размеры исследуемого трубопровода следующие: длина прямолинейного входного участка: 11 = 200 мм; радиус колена Я = 200 мм; длина прямолинейного выходного участка 12 = 200 мм; наружный и внутренний диаметры = 22мм и = 20 мм.

В качестве опор принимается жёсткая

Рис.1. Основные геометрические размеры исследуемого трубопровода

заделка по обоим концам трубопровода.

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

Традиционно, при проведении расчётов параметров пульсаций рабочей жидкости в трубопроводных системах, в качестве динамической нагрузки используется входной импеданс присоединённой гидравлической цепи [9]. Особенностью программного комплекса АКБУБ является то, что акустическую нагрузку здесь можно задать только с помощью коэффициента звукопоглощения а. Импеданс и коэффициент звукопоглощения связаны следующим соотношением [10]:

ячеи волны;

х

рс

+

У

рс

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

^ - площадь поперечного сечения трубопровода;

] = 4-1 - мнимая единица;

2Ы

У

"N2 +1)-(2 -1) cos8 ’

(n2 -1) sin 8 (N2 + l)-(N2 -1) cos28 ;

N

a

1 +У1-

1 -V 1 -a

- коэффициент сто-

8 =

x-

ЯЛ

4

2п

—— - фаза коэффициента отра-Л

жения;

a - коэффициент звукопоглощения;

х0 - координата ближайшего к концу трубы

минимума давления.

Таким образом, по известному импедансу можно рассчитать коэффициент звукопоглощения. Для случая бегущей волны импеданс присоединённой гидравлической цепи равен волновому сопротивлению трубопровода Z=Z При этом амплитуда волны давления постоянна по длине трубопровода (при условии отсутствия трения), то есть Хо=0, тогда a =1. Для случая стоячей волны Z=<^>. Известно [8], что для рассматриваемой стоячей волны минимум пульсаций расположен в сечении, удалённом на четверть длины волны от акустически закрытого конца, то есть хо=Л/4, что по умолчанию реализовано в программном комплексе ANSYS.

Задача решается в трёхмерной постановке. Для моделирования жидкости используется восьмиузловой пространственный элемент FLUID30. Этот элемент предназначен для описания жидкой среды, а также её взаимодействия с твёрдым телом. Обычно применяется для исследования звуковых волн и различных динамических процессов. Узлы используемого элемента имеют 4 степени свободы: перемещения по координатам X, Y, Z и давление. Математическое описание элементов данного типа базируется на волновом уравнении:

~г-df-V2 P = 0,

с dt

где Р - давление в жидкости; t - время;

V = {L}T

дх dy dz

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

Mp Jp }+ [kp k}+ Po [r. j u. )={0}

/Э э эЛ

Рис.2. Краевые условия исследуемой системы где |мр ]=4 jln }n }Td (voi) - матрица

vol

s

s

x

max

масс жидкости;

[кр ] = JW№ (vol) - матрица жёсткости

уо!

жидкости;

ро &]=РоN7 (я) - матрица

масс виброакустического взаимодействия; {Ре} - вектор узловых давлений;

{/е} - вектор узловых перемещений;

{V} - функция формы элемента давления; {/V'} - функция формы элемента перемеще-

ния;

М={Ф };

{п} - единичный вектор нормали к поверхности 5.

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

[M, ] [0]

] MPA

[C, ] [0Г

[0] C,P J

[к, ] к J

[0] [кр]

+

+

где [к, ] - матрица масс твёрдого тела;

к"\=р. [к, ]т;

, ] - вектор узловых перемещений твёрдого тела;

[С, ] - матрица демпфирования твёрдого тела; матрица жёсткости твёрдого тела;

[к,]

[к" ]=-[R, ];

Р} - вектор внешних сил.

Нагрузка (коэффициент звукопоглощения) задаётся операцией ГМРБ. Математически это выражается добавлением слагаемого

СР } в левую часть волнового уравнения

в матричной форме.

[ср ]=а ]{# л# (я) - матрица дем-

с я

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

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

Для моделирования стенки трубопровода используется десятиузловой тетраэдральный (с промежуточными узлами на рёбрах) пространственный элемент с неплоскими гранями 80ЬГО92. Этот элемент используется для трёхмерного моделирования твёрдых тел. Его узлы имеют 3 степени свободы: перемещения по осям X, У, Z. Математическое описание данного типа элемента для рассматриваемого случая акустического анализа определяется уравнением:

[М м+[с ]{и }+[к{ }={ра }

[к ] - матрица масс твёрдого тела;

[С ] - матрица демпфирования твёрдого тела; [К ] - матрица жёсткости твёрдого тела;

{Ра} - вектор внешней нагрузки.

Задаются следующие свойства материала трубопровода: модуль упругости, коэффициент Пуассона, плотность. При этом учёт физических свойств материала трубопровода позволяет смоделировать известное свойство уменьшения скорости звука в жидкости при наличии податливости стенки трубопровода [7]:

с =

1 +

Рвн E ж

8 E

ст тр

где с0 - скорость распространения волн дав-

с

0

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

Ввн - внутренний диаметр трубопровода;

8ст - толщина стенки трубопровода;

Еж - модуль упругости рабочей жидкости;

Етр - модуль упругости материала трубопровода.

Величина скорости звука в жидкости при условии абсолютно жёстких стенок трубопровода принимается равной с=1300 м/с. Аналитический расчёт для случая податливых стенок по вышеприведенной зависимости показал следующее значение скорости звука с=1214 м/с. При численном моделировании скорость звука получилась равной с=1202 м/с. Расхождение последних двух величин связана с погрешностью численных методов и не влияет на физическую адекватность получаемых результатов.

Геометрическое моделирование трубопровода и жидкости имеет ряд особенностей. Сначала производится построение поперечного сечения стенки трубопровода и заключённой в нём рабочей жидкости. Затем выполняется вытягивание построенных поверхностей вдоль заранее заданного контура при помощи операции УОЯАО. Далее, при помощи операции УОШЕ, производится склеивание объёмов жидкости по наружным поверхностям с соответствующими объёмами стенок трубопровода. При этом полученные объёмы стенки и рабочей жидкости имеют единую общую граничную поверхность. Это необходимо для дальнейшего задания взаимодействия пристеночных элементов со стенкой трубопровода. После этого на геометрические объекты наносится сетка конечных элементов. Следует отметить, что нанесение сетки конечных элементов на этапе построения сечения с дальнейшим вытягиванием вдоль оси трубопровода нецелесообразно, так как полученные при этом объекты не будут иметь общей единой поверхности, и моделирование взаимодействия стенки и жидкости станет невозможным.

Расчёты для случая акустически закрытого конца на частоте колебаний 150 Гц по-

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

А = —

А„

Распределение пульсаций давления по длине (рис.3) является характерным для стоячей волны: на акустически закрытом конце реализуются максимальные амплитуды. Моделируемая частота колебаний на первом этапе была специально выбрана достаточно низкой - ниже первой резонансной частоты как по гидравлической, так и по механической подсистеме. В этом случае амплитуды пульсаций давления монотонно уменьшаются по длине от выходного к входному сечению. Для оценки адекватности результатов численного моделирования пульсационного состояния были проведены аналитические расчёты по известным моделям [7, 8]. В соответствии с [8] отношение амплитуды пульсаций на акустически закрытом конце к амплитуде в текущем сечении определяется по соотношению:

2

1 + 008

I А=1,00...1 1 А=1,02...1 ! А=1,04...1 I А=1,06...1 I А=1,08...1 I А=1Д0...1 I I А=1Д2...1

I А=1Д4...1 М А=1,16...1

Рис.З. Распределение пульсаций подлине трубопровода для случая акустически закрытого конца при частоте 150Г ц

р

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

£ - коэффициент фазы; е = —;

с

— - круговая частота колебаний.

Различие в результатах численного и аналитического моделирования не превышало 0,4% и объясняется погрешностями численных методов. Следует также отметить, что расчёт с использованием программного ком -плекса АКБУБ показал софазный характер колебаний давления в трубопроводе с акустически закрытым концом на дорезонансных частотах, что также согласуется с известными аналитическими моделями [7, 8].

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

2п • / • х

роскоростеи V = -

м2 • с

у=(0...1,01)10'9

У= (1,01. ..2,02)10°

V— (2,02...3,03)10 9

У=(3,03...4,03)10'9

у= (4,ОЗ...5,05) 10 9

(5,05...б,05)10 9

(б,05...7,07) 10 9

\г= (7,07...8,09) 10 9

ш (8,09...9,09) 10 9

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

В ходе моделирования также было получено распределение вибропараметров по длине трубопровода для обоих исследованных случаев (рис.4, 5).

Результаты представлены в виде облас-теи равных значении для относительных виб-

Рис.4. Распределение виброскорости подлине трубопровода для случая акустически закрытого конца при частоте 150Гц

носительнои виброскорости для случая реализации бегущей волны составляет 7,46-10-9, в то время как для нагрузки на акустически закрытый конец - 9,09-10-9.

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

Авх кг

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

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

\ї=(0...8,29)10 10

її=(8,29...1б,б0)10'1

У= (1,66...2,49) 10"9

у= (2,49...3,32) 10'9

У= (3,32...4,15) 10"9

У=(4,15...4,98)10"9

ш ч= (4,98...5,81)109

ш у=(5,81...б,64)10"9

ш У=(б,б4...7,4б)10 9

Рис.5. Распределение виброскорости по длине трубопровода для случая нагрузки на неотражающее сопротивление при частоте 150 Гц

с

V

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

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

СПИСОК ЛИТЕРАТУРЫ

1. АкимовМ.Ю., Вельмисов П.А., Милушева С.Д. О днамических и статических деформациях трубопровода // Proceedings of the XXII Summer School "Application of Mathematics in Engineering". Technical University of Sofia. Bulgaria, Sozopl.

2. Гладких А.Г., Хачатурян С.А. Вибрации в трубопроводах и методы их устранения. М.: Машгиз, 1959.

3. Кондратов Н.С. О параметрических колебаниях трубопроводов // Вибрационная прочность и надёжность авиационных двигателей. Вып. XIX. Куйбышев, 1965.

4. ЛеньшинВ.В., ШахматовЕ.В. Виброаку-стическая модель трубопровода // Техно-

парки, инновации, конверсия. Самара: СГАУ, 1995.

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

6. Сапожников В.М., Лагосюк Г.С. Прочность и испытания трубопроводов гидросистем самолётов и вертолётов. М.: Машиностроение, 1973.

7. Попов Д.Н. Динамика и регулирование пневмо- и гидросистем. М.: Машиностроение, 1977.

8. Прокофьев А.Б., ШестаковГ.В. Исследование частотных характеристик отрезка трубопроводной линии при несогласованной нагрузке в виде акустически закрытого трубопровода // Ракетно-космическая техника. Расчет, проектирование, конструирование и испытания космических систем. Серия XII. Вып.1. Самара: ВКБ РКК Энергия, 1998.

9. Шорин В. П. Устранение колебаний в авиационных трубопроводах. М.: Машиностроение, 1980.

10.ЛепендинЛ.Ф. Акустика. М.: Высш. школа, 1978.

MODELING OF PIPELINE VIBROACOUSTICAL CHARACTERISTICS BY THE FINITE ELEMENT METHOD

© 2002 G.M. Makaryarnts1, A.B. Prokofjev1, E.V. Shakhmatov2

1 Samara State Aerospace University

2 Institute of Machines Acoustics, Samara

Method of modeling of pipeline vibroacoustical characteristics with loads of working liquid pulse flow is proposed. Method uses modern CAE system ANSYS based on FEA. Calculations of elbow pipe vibration for different cases of connected hydraulic circuit dynamic loads are made.