УДК 621.6.01
ЭФФЕКТИВНОСТЬ ГАСИТЕЛЕЙ КОЛЕБАНИЙ
© 2004 А Н. Головин
Самарский государственный аэрокосмический университет
Одним из методов устранения колебаний рабочей среды в трубопроводных системах машин является применение акустических гасителей. Гасители изменяют динамические свойства систем и структуры волновых процессов. В статье показано, что для эффективного применения гасителей необходимо учитывать характеристики источников колебаний и нагрузок или применять такие устройства, которые обеспечивают гарантированное ослабление колебаний в трубопроводных цепях с любыми динамическими свойствами.
Распространяющиеся в трубопроводах колебания давления и расхода, образуют совокупность бегущих и стоячих волн. Введение в трубопроводную систему гасителя изменяет ее динамические свойства. Происходит трансформация волнового процесса. В результате перераспределяются амплитуды колебаний давления и расхода на участке системы до гасителя и на участке системы после гасителя. Очевидно, что на входном участке амплитуды колебаний зависят от входного импеданса гасителя. На выходном участке структура волнового процесса более сложная. Она определяется отражениями в сечениях стыка до гасителя и за гасителем и его аккумулирующими и демпфирующими свойствами. Причем характеристики гасителя, системы и все изменения в волновом процессе взаимосвязанные.
Задача влияния гасителей на колебания в выходном участке системы известна как задача о вносимом затухании и наиболее пол-
но изучена в работе [1] при сосредоточенных параметрах гасителя. Допущение о сосредоточенности параметров позволило упростить формулу коэффициента вносимого затухания КВН , но одновременно с этим была потеряна оценка влияния волнового процесса на вносимое затухание. Для устранения недостатков существующей модели КВН исследуем соответствующие зависимости при распределенных параметрах гасителя и системы.
Рассмотрим это подробно.
Установим в произвольную трубопроводную систему, рис.1а, имеющую независимый источник колебаний с импедансом ZИ вместо отрезка трубопровода длиной I с активным волновым сопротивлением ZВ, гаситель колебаний произвольной структуры. После введения гасителя амплитуды колебаний давления и переменного расхода рабочей жидкости в системе изменятся. Расчетная схема для определения их изменения изображена на рис. 1б.
3
б)
Рис. 1. Гидравлическая и расчетная схемы трубопроводной цепи: На расчетной схеме параметром ||А|| без дополнительных обозначений и с подстрочными индексами I и L обозначены, соответственно, матрицы передачи: гасителя, "удаляемого" отрезка трубопровода и трубопровода, соединяющего источник колебаний с гасителем
2
2
На схеме:
- часть трубопроводной системы, расположенная слева от гасителя (входной участок), характеризуется входным импедансом
г1=г1(®);
- часть трубопроводной системы, расположенная справа от гасителя (выходной участок) характеризуется входным импедансом 22=22(ю ).
Используя результаты работ [1], запишем выражения для амплитуд колебаний давления, которые были в сечении 2-2 (рис. 1б) до установки гасителя и стали в этом сечении трубопроводной системы после его установки (соответственно для Р 2 и Р 2* ):
Р2 = Р X
2 • г 2 • г
(1)
Р2 = Р X
Р
К - 2
ЛВН р** .
(3)
После соответствующих преобразований получим:
КВН --
Кс 7 71 + 7с1 7в 7 + 7 с 2
К стр V 7 с1 71 + 7в 7с 2 7 2 + 7е
- (
1 -
К
7 - 7
с1
-7 ^
2 с 2
7 + 7
V1 ^ ^с1 у
V 7 2 + 7с 2 У
- ] ^ • Утр (
1 -
К
2
" стр
(4)
7 - 7 ^в
^ - 7 ^
V 71 + 7в у
7 + 7
V 2 ^ ^в У
где:
,Гтр
- К
стр >
,Уг
-К
К — К
стр стр *
х _2 72 • 2 '7с1_
^ ^ + ^ ) ■ (7 2 + 7с2 )-е- Уг • (7с2 - 72 ) •(7с1 - 71)]
(2)
где: i — 1 или i — 2; — (ш) - волновое
сопротивление гасителя; У об — У об (ш) - постоянная передачи трубопровода, в который устанавливается гаситель; у а — у а (ш) - постоянная передачи гасителя.
Подставим выражения (1) и (2) в формулу для Квн (3).
,(ю); Кс — Кс (ю) -
коэффициенты собственного затухания удаляемого отрезка трубопровода и гасителя;
УГ — уГ(ю); (рТР — (рТР(ю) -характеристические фазовые сдвиги гасителя и удаляемого отрезка трубопровода;
Уг (ю) — (ю) ] — Л(ю)в(ю) + л/ В{ю)С (ю) ] ш • £
УТР ; а - скорость звука; j= 1.
В общем случае, при известных характеристиках системы, выражение (4) позволяет:
- определить эффективность действия гасителей конкретных схем или
- определить характеристики гасителя, необходимые для обеспечения требуемой величины КВН и синтезировать структуру гасителя.
Для решения этих задач необходимо знать зависимости входных импедансов участков системы . Однако их определение на практике затруднено и ограничено в большинстве случаев лишь приближенными оценками. Поэтому проведем анализ выражения (4) с целью его упрощения и отыскания способов выбора типа гасителя и его характеристик при ограниченной или неполной информации о динамических свойствах и характеристиках системы. Это также позволит получить обобщенные рекомендации проектирования гасителей для систем с известными характеристиками.
Введем дополнительные обозначения и,
учитывая, что КСТР « 1, формулу (4) запишем в виде
X
г
е
х
Квн к Кс •©! •©2. (5)
В соотношение (5) входят параметры, которые имеют следующие зависимости:
©, =
7В 7 -г + 7а
i 7а 7,г + 7В
(6)
М =
1 -
- j'2'Фг
К
(
7 - 7
c1
7 + 7
V 1 ^
\ (
Л
7 - 7
2 c 2 V 7 2 + 7c 2 J
- (7)
1 - q-1■ 2■ ртр
f7i - 7е Л ( 7 - 7 ^
V 71 + 7e j
7 + 7
V^ 2 ^ ^е J
- (8)
М=
Т Т ^
1 _ 1 2 ^
1Р
K
(9)
В формуле (9): р = р1 + р2 - 2 ■ рг;
Т =
1 - m.,. ■q1
1 + m , ■Q1
Рг = arg
1 - m , ■ Q
'■■ РсЙ Л
v1 + mcii-q1 рс" j
; mcii =
Z
ci
Z
Z
ci
функции взаимодействия гасителя и системы в сечении входа i = 1 и в сечении выхода i = 2;
функция, характеризующая затухание, обусловленное несовпадением входных импедан-сов участков системы и соответствующих волновых сопротивлений гасителя;
функция, характеризующая затухание, обусловленное несовпадением входных сопротивлений участков системы и волнового сопротивления удаляемого отрезка трубопровода или трубопровода, в который устанавливается гаситель. В формулах (5), (8) знак приближенного равенства стоит из-за пренебрежения активными потерями по длине удаляемого отрезка трубопровода, т.е. принято, что К ~ 1
Л СТР ~ 1 •
Значения функций "м" и "К" зависят: " м" - от характеристик гасителя и его места установки в системе;
" К" - от предполагаемого места размещения гасителя в системе и относительной
длины 1 = 1/X удаляемого отрезка трубопровода. Параметр X - длина волны подавляемых колебаний.
Рассмотрим выражение (7). Перепишем его в следующем виде
Феи = 7
V Zl У
Из соотношения (9) следует, что на величину функции М влияют параметры Т, Ф1
и коэффициент Кс . Значения Т, ф1, в свою очередь, определяются параметрами тсц и Фс11.
В общем случае параметр тсц может принимать любые значения, определяемые неравенством 0 < тсц < ^.
Диапазон изменения аргумента фсц зависит от характера (инерционный, упругий или активный) импедансов Zci и 21. Для определенности будем считать:
i i к
- если |фси| < —, то импедансы 2С1 и
21 одного характера;
К | |
- если — < |фсц| < К, то импедансы 2С1
и 21 разного характера.
Запишем выражения для параметров Т и ф1 в следующем виде:
Ti =
1 - 2 ■ mcii ■ cosФей + mcii
1 + 2 ■ mcii ■ cosФей + mai
(10)
Фi = ± arceos
1 - m2
д/1 + ■ cos2■ Фей + m^i " (11)
В формуле (11) знак "+" следует выбирать, когда sin< 0, а знак "минус" при
sin Фей ^ 0
Анализ выражения (10) показывает,
что величина коэффициента Т с учетом пределов изменения параметров шсц и фсИ , может изменяться:
при однохарактерном сочетании характеристик - в диапазоне
1 - ш,
XII
1 + ш
< т < 1
(12)
С11
или в пределах
0 < Т1 < 1;
при разнохарактерном сочетании волновых сопротивлений гасителя и соответствующих входных сопротивлений участков системы, т.е. когда — < |фсц| < —, коэффициент
Т1 может принимать значения в диапазоне
1 + ш,
1 < Т1 <
С11
|1 - ш
(13)
С11
или
1 < т < го.
Аргумент ф 1 при изменении величины
шсц в диапазоне 0 < шсц < го в обоих вариантах сочетания характеристик может принимать любые значения в пределах
0 < |ф¿| < го.
На рис. 2а, б приведены графики, иллюстрирующие изменение предельных значе-
ний функций Т и | ф11 в зависимости от параметра шсц при различных значениях аргумента ф сц . Графики симметричны относительно ординаты и абсциссы со значениями
Т — 1
— 1, — . При условиях
ш
с11
шсц ^ 0 и шсц ^ го диапазон изменения
параметров Т и | ф 11 сужается, а сами величины стремятся к своим предельным значениям: Т ^ 1, |ф11 ^ 0 или ф ^ — |.
Для исследования возможного диапазона изменения аргумента фГ воспользуемся известным соотношением, справедливым для пассивных четырехполюсников
А(ш) • D(ш) - В(ш) • С(ш) — 1 (14) и запишем выражение для аргумента фГ в виде
.(15)
У Г —
— ащ[^е[л(ю)п(ю)] + у \т[л(ю)о(ю)] +
+ Л]^е[л(ю)п(ю)]-1} + у \ш[ л(ю)й (ю)] ]
Из анализа формулы (15) следует, что параметр уГ может изменяться в пределах
0 < уГ < п . При вычислениях аргумента уГ для конкретной схемы гасителя следует учитывать, что извлечение квадратного корня из комплексного числа действие неоднозначное.
Рис. 2. Зависимости предельных значений параметров Т и Ц.| от соотношения волновых сопротивлений гасителя и соответствующих им входных импедансов участков системы
Поэтому фактическую величину фГ можно
определить только после дополнительных исследований. Например, построением векторных диаграмм изменения амплитуд давления и расхода в сечениях трубопровода на входе и выходе гасителя.
Проведенный анализ показывает, что
параметры |ф| и Т в формуле (9), соответственно, могут принимать значения в пределах 0 < |ф| < 4к и 0 < Т < ^ . Причем, если на диапазон изменения Т оказывает влияние характер фсц , то на пределы изменения
Ф данное влияние менее существенно.
На основе полученных результатов определим диапазон и характер изменения функции М. После алгебраических преобразований соотношение (9) примет вид
м =
1 „ т т
1 - 2—— - cos фн—
2
к2
к4
(16)
где Т=Т1«Т2.
Допуская реализации предельных значений функции cos(ф), запишем в виде неравенства границы возможных изменений функции М в виде
1 -
Т
К2
< М <
1 +
Т
К2
(17)
Проведем анализ соотношения (17) при однохарактерных, |фсц| < —, т < 1, и разно-
характерных, к < |ф с111 < К; т > 1,
сочетани-
ях, соответствующих характеристик гасителя и примыкающих к нему частей системы.
В первом случае, с ростом значения коэффициента Кс, границы возможных значений функции М непрерывно сужаются, а ее величина стремится к единице. Во втором случае диапазон изменения функции М с возрастанием коэффициента Кс вначале расши-
ряется, достигая при условии Т = К2С значений [0; 2], а затем, как и в первом случае, сужается. Функция М стоит в числителе выражения для Квн (5). Следовательно, для получения больших значений коэффициента Квн, необходимо иметь максимальные гарантированные значения функции М. При одноименных характерах аргументов взаимодействующих характеристик гасителя и системы это будет всегда при условии Кс >> 1.
При различных характерах взаимодействующих характеристик гасителя и системы большие значения функция М принимает при значениях коэффициента Кс много меньших и
много больших значений ^/Т. Иначе, даже при увеличенных значениях коэффициента собственного затухания, Кс, возможна реализация минимальных значений функции М. На рис. 3 приведены графики изменения предельных значений функции М в зависимости от изменения коэффициента Кс и при
различных параметрах комплекса Т.
Выражение для функции N подобно выражению для функции М. Ее величина зависит от относительной длины гасителя 1 и
места размещения гасителя в трубопроводной системе. Так как гаситель устанавлива-
Рис. 3. Зависимости предельных значений функции М при различных значениях параметра Т
ется вместо отрезка трубопровода, волновое сопротивление которого с достаточной для практики точностью принимается активным,
всегда
arg
^ Z ^
7
V^ в
Ж
- . Поэтому максимально
возможный диапазон изменения коэффициента N равен [0,2], а его границы в зависимости от значения произведения параметров
N • N2: Ni -
Z - Z 1 В - Z 2 - Zb
Zi + Zb Z 2 + Zb
опре-
©i-
1 Zci 2 + IN 2 + 2- Zci IN •cos^ci -Фzi)
Zci 1 + 2- Zi •cos Фа 1- |2 , + Zi
(19)
где
Zci
Z
ci
Z,
Zi
Zi
Z,
Фс1 - arg(Zci); 9zi - arg(Zi) - соответственно, модули и аргументы волновых со-
противлении гасителя и входных сопротивлений участков системы.
Примем в формуле (19) возможные предельные значения косинусов и запишем выражения для наименьших и наибольших значений 0.:
inf © -
1 N1 2 + Zci 2
Zci 1 + 2- N1 + Zi 2
(20)
деляются следующим неравенством
1 - N -N2 < N < 1 + N1•N2. (18)
Функция N стоит в знаменателе выражения (5), поэтому для получения наибольшего КВН необходимо, чтобы значение N было минимальным или, по крайней мере, N<1. Это возможно, например, при установке гасителя, продольные размеры которого много меньше длины волны колебаний в пучности или узле давления. Однако, если частоты колебаний переменны и, следовательно, переменны коэффициенты N N что характерно, например, для авиационных трубопроводных систем, то расчет гасителя целесообразно вести исходя из худшего случая, чему соответствует максимальное значение N=2. При этом реальная эффективность гасителя будет всегда не ниже расчетной. Данное упрощение является справедливым также на этапе выбора типа гасителей, так как в формулу для функции N не входят параметры, определяющие характеристики гасителя.
Рассмотрим функции взаимодействия
01. Приведем выражение (6) к виду
sup ©1 -
_ 2 _ _ _
1 Zi + 2- Zci Zi + Zci
Zci
1 +
Zi
(21)
Графики функции inf ©i и sup ©i приведены на рис. 4.
Из анализа графиков следует, что зависимости имеют сложную пространственную конфигурацию с экстремумами и симметрич-
ны относительно значении
Zi
- 0
и
Zci
- 0. Для расчетов и исследовании фун-
кции inf ©i и sup ©i следует принимать зна-
чения
Zi
> 0
и
Zci
> 0 . В этой области изменения параметров функции взаимодействия 01 могут принимать любые значения из диапазона
0 < 01 < ю. (22)
Разделим этот диапазон на два интервала. В первом интервале функция 01 имеет значения в пределах
0 <©1 <[©i ]
(23)
а во втором интервале значения функции 01 определяются неравенством
[01 ]<01 <ю . (24)
В соотношениях (23) и (24) параметр [01 ] устанавливает допустимую величину функции взаимодействия, которая рассчитывает-
2
2
НАИМЕНЬШИЕ НАИБОЛЬШИЕ
Рис. 4. Зависимости наибольших и наименьших значений функций взаимодействия от изменения волновых сопротивлений гасителя и входных импедансов присоединенных к гасителю частей
ся из соотношения (5) при требуемом коэффициенте вносимого затухания, Квн.
Условие (24) определяет гарантированное уменьшение амплитуд колебаний за гасителем. Если значения функции 01 находятся в интервале (23), то амплитуды колебаний вместо уменьшения могут увеличиться и возникнут условия для резонансных режимов в системе.
Исследуем выражение (19) в зависимости от характеристик гасителя, а именно в
зависимости от изменения модуля, 2с1 , и
аргумента, фс1, волнового сопротивления гасителя. В результате получим, что функция (19) при выполнении равенства
Zci
Zi
(25)
min 0 i =
2-
Zci
•f1 + C°s(9 ci -Ф zi )1
1 + 2-
Zci -cosФzi +
^ |2 . (26) Zci
противления гасителя и от аргумента входного импеданса системы. Дальнейшие исследования на экстремумы показывает, что функция выражения (26) при условии
Zci
=1
(27)
имеет минимум. Формула, описывающая траектории минимальных значений 01, имеет следующий вид
Параметры функции min 0 j зависят от величин модуля и аргумента волнового со-
принимает наибольшие значения, которые рассчитываются по формуле
sup[minQ;]= I'+ Г(ФС'~ФZ1} . (28) \ 1 + cos ф z'
Функция sup[min 0' ] зависит от
аргументов фС' и ф z', которые могут варьироваться в пределах
о <|ф С'| , (29)
0 <|ф zi| • (30)
Для этих изменений аргументов интервал значений функции sup[min 0' ] максимальный и находится в следующих пределах
0 < sup[min 0' ]< V2 • (31)
Анализ формулы (28) показывает, что интервал значений функции sup[min 0' ], при уменьшении интервала изменения аргумента ФС', будет сужаться и при значении аргумента
Фс1 = 0 (32)
функция
sup[min 01 ]=1, (33)
то есть становится независимой от акустических характеристик системы.
Условия (27), (32) определяют требования к волновому сопротивлению гасителя, при котором реализуется гарантированное
значение sup[min 01 ]. Для этого волновое сопротивление гасителя должно быть активным и равным волновому сопротивлению трубопровода, в который устройство будет монтироваться. При таких параметрах гаситель будет эффективно работать в максимально широком диапазоне частот колебаний.
Определим требования к гасителям, когда они проектируются для систем, имеющих специфические динамические характеристики.
Так для гасителя, устанавливаемого в систему с предельными характеристиками,
когда Ъ1 = 0 или Ъ1 = ю, выражения для функций взаимодействия 01 следующие:
если = 0, то 01 =
^С1
если = ю, то
01 =
Ъс1
(34)
(35)
Из соотношений (34) и (35) следует, что если гаситель проектируется для систем с предельными характеристиками (или близкими к ним), то для получения максимальных
значений функций 01 необходимо:
при условии
Ъ1
^ 0
применять гаси-
тель с возможно большей величиной
Ъс1
при условии
Ъ1
^ ю применять гаси-
тель с возможно меньшей величиной
Ъс1
На основании проведенных исследований примем М» 1 и N=2 и приведем выражение (5) к виду
КВН >КС-0-0,
'1
(36)
и рассмотрим три важнейших случая.
1. Если в выражениях (6) и (19) принять
Ъс1 = ЪВ (ъС1 = 1), то значение функции
01 = 1, и соотношение (36) примет вид
Квн> Кс 02
(37)
Из формулы (37) следует, что минимальная величина вносимого затухания не зависит от импеданса источника колебаний. Гаситель можно проектировать без учета характеристик источника колебаний. Устройство, спроектированное для конкретной гидравлической системы, будет эффективно работать при любых источниках колебаний.
2. Полагаем в выражениях (6) и (19)
2с2 = 2В (ъс2 = 1). Тогда 02 = 1 и из (36) получим
Квн> ^ 0:
(38)
Теперь минимальная величина вносимого затухания не зависит от импедансных характеристик нагрузки или присоединенной цепи. Гаситель, спроектированный для конкретного источника колебаний, будет эффективно работать в системе с любыми характеристиками.
3. При значениях волновых сопротивлений 2Л = 2 с2= 2В функции 01 = 02 = 1 и условие эффективности гасителя на выходном участке системы запишется следующим образом
Квн > К-
(39)
В этом случае минимальная величина вносимого затухания определяется только коэффициентом собственного затухания, т.е. схемой и конструкцией самого гасителя. Такой гаситель будет эффективен в любой системе с любым источником колебаний.
Из анализа выражений (36)...(39), с учетом возможного диапазона изменения
функции 01 следует, что гасители, у которых одно волновое сопротивление является активным и равным по величине волновому сопро-
1
тивлению "удаляемого" отрезка трубопровода, могут быть эффективнее симметричного активного гасителя. Для этого у несимметричного гасителя необходимо выбрать волновое сопротивление с противоположной стороны так, чтобы обеспечивалось наибольшее значение функции 01 или, по крайней мере, 01>1.
Выводы:
1. Для достижения наибольшей эффективности действия гасителя в трубопроводах его волновые сопротивления и соответствующие входные акустические сопротивления участков системы по обе стороны от гасителя должны иметь одинаковый характер. Если входной импеданс участка трубопроводной цепи трудноопределим или изменяется в широких пределах, то необходимо чтобы обращенное в сторону этого участка цепи волновое сопротивление гасителя было активным. Величина активного волнового сопротивления гасителя должна быть равна величине волнового сопротивления "удаляемого" отрезка трубопровода, вместо которого будет установлен гаситель.
2. Показано, что для управления волновыми процессами в трубопроводных системах с произвольными или изменяющимися в широких пределах динамическими характеристиками необходимо применять симметричные гасители с активными волновыми сопротивлениями. Такие гасители являются универсальными. Однако фактическая эффективность действия активных симметричных гасителей определяется величиной коэффициента собственного затухания Кс и, как правило, невелика. В трубопроводных системах, в которых входные импедансы могут быть измерены или рассчитаны, целесообразно использовать несимметричные гасители. Эффективность действия таких гасителей, за счет дополнительных отражений, выше эффективности действия, которую могут обеспечить симметричные активные гасители.
СПИСОК ЛИТЕРАТУРЫ
1. Шорин В.П. Устранение колебаний в авиационных трубопроводах. М.: Машиностроение, 1980.
EFFICIENCY DAMPERS OF FLUCTUATIONS
© 2004 A.N. Golovin
Samara State Aerospace University
One of methods of elimination of fluctuations of a working environment in pipeline systems of machines is application acoustic dampers. Dampers change dynamic properties of systems and structures of wave processes. In article it is shown, that for effective application dampers it is necessary to take into account characteristics of sources of fluctuations and loadings or to apply such devices which provide guaranteed easing fluctuations in pipeline circuits with any dynamic properties.