УДК 621.6.01
ХАРАКТЕРИСТИКИ АКУСТИЧЕСКИХ ГАСИТЕЛЕЙ ПРИ НЕСТАЦИОНАРНОЙ МОДЕЛИ ГИДРАВЛИЧЕСКИХ ПОТЕРЬ
© 2004 А Н. Головин
Самарский государственный аэрокосмический университет
Приведена методика исследования динамических характеристик акустических гасителей колебаний при различных моделях движения рабочей среды.
Для исследования динамических процессов в трубопроводных линиях применяют три модели течения рабочей среды. Модель гидравлической линии без потерь. Модель с линейным трением, когда потери пропорциональны средней скорости. И модель, учитывающую эффекты вязкого трения и частотную диссипацию колебаний. Диссипа-тивная модель наиболее точно описывает протекание колебательных процессов и является в ряде случаев предпочтительнее других, особенно при изучении измерительных систем с распределенными параметрами. Основная сложность применения диссипа-тивной модели связана с использованием математических вычислений с комплексными функциями. Однако существующая необходимость уменьшения габаритов устройств, при увеличении их эффективности действия оправдывает детальное исследование динамических процессов, в том числе с эффектами вязкого трения.
Целью статьи является изучение собственных характеристик гасителей при введении в расчетные формулы параметров, учитывающих вязкость рабочей среды и сравнение этих характеристик с идеальными зависимостями, полученными при отсутствии потерь в элементах гасителей.
Акустические характеристики гасителей, в том числе при распределенных параметрах обычно исследуют без учета гидравлических потерь в реактивных элементах устройств, иногда при квазистационарной модели гидравлических потерь. В качестве реактивных элементов применяют отрезки трубопроводов различных протяженностей и
диаметров, которые определенным образом стыкуются друг с другом. Типовые гидравлические схемы таких гасителей приведены на рис. 1. У рассматриваемых гасителей цилиндрический проточный инерционный канал и упругая полость цилиндрической формы, имеют одинаковую протяженность и конструируются либо коаксиально, либо байпас-но. Переход к различным схемам устройств осуществляется за счет изменения сопротивлений дросселей, шунтирующих поток через расширительную полость. Так приняв
Zl = 0, Z2 = ю получаем схему реактивного акустического фильтра низких частот, и т.д.
Схемы гасителей и их характеристики, построенные на базе модели распространения колебаний без потерь, изложены в работе [1]. Не учет диссипации энергии колебаний упрощает формулы, описывающие исследуемые процессы. Это позволяет преобразовать зависимости характеристик до выражений, в которых соотношения между параметрами элементов устанавливаются из условия унификации всех исследуемых схем гасителей. Такие соотношения также определены в
работе [1]. Поэтому для обобщения результатов исследований характеристик гасителей при нестационарных процессах будем использовать уже установленные взаимосвязи между параметрами элементов этих устройств полученные при не учете явлений диссипации колебаний.
Обобщенная расчетная схема гасителей в виде соединения четырехполюсников, описывающих свойства элементов гасителей, приведена на рис. 2. На этой схеме динами-
Рис.1. Принципиальные гидравлические схемы базовых однозвенных гасителей и их электрические аналоги: 1 - расширительная полость; 2 - проточный канал; 3 - гидродроссели; Бк, Бцтр - площади поперечных сечений расширительной полости и проточной трубки; I - протяженность гасителя, Х!_ - инерционное сопротивление проточного канала; ХСПР - приведенное упругое сопротивление расширительной полости; R. - сопротивление гидродросселей; Т . - волновые сопротивления гасителей
ческие процессы в инерционной трубке и в расширительной полости представляются четырехполюсниками с коэффициентами передачи, полученными из решения уравнения движения вязкой жидкости по цилиндрическому трубопроводу [2].
Для сохранения адекватности зависимостей движение вязких жидкостей в цилиндрических трубопроводах и в реактивных элементах конструкции гасителя, примем у гасителя байбасное расположение расшири-
тельной полости и проточного канала.
Основные формулы, на базе которых сформированы параметры "телеграфных" уравнений следующие.
Коэффициент распространения колебаний является комплексным и имеет вид
у0Ю)=±(8 + je). (1)
Знак "+" в формуле (1) принимается в расчетах, когда исследуются колебания по направлению распространения потока. При изучении колебаний, распространяющихся
Рис. 2. Расчетная схема гасителя колебаний обобщенной структуры
против направления потока в формуле (1) принимается противоположный знак. Параметр § является коэффициентом затухания, а параметр в - коэффициентом фазы. Коэффициент затухания определяет уменьшение по длине трубопровода амплитуды давления или амплитуды скорости в волне возмущения, распространяющегося по трубопроводу с фазовой скоростью
ш
с = — в '
§ = в 0
Г ,- Л
к2 1 1 1
1 + —г -1
2кх [1 к2 ]
(3)
в = в 0
Г 1- Л
к2 1 1 1
1 + —г +1
2кх [1 к2 J
(4)
в 0 = ш
Р 0
В
тр
(5)
к =
(2л/ш - 1|4ш - 2л/ш +1) 4ш - л/ш '
4ш(2л/ш -1)
(6)
ко =
(7)
Коэффициенты к^ и ко зависят от безразмерной частоты колебаний
(2)
где ш - круговая частота колебаний.
Параметры § и в для гармонических колебаний потока рассчитываются по формулам:
— Rтр ш
ш =
8у
(8)
В формуле (5) р0 - плотность рабочей
среды; Втр - приведенная упругость рабочей
жидкости.
Коэффициент в 0 определяет фазовую
скорость колебаний в трубопроводе при не учете потерь. Соотношения для вычисления
коэффициентов к^ и ко следующие:
где ^Ктр - радиус внутреннего сечения трубопровода, по которому распространяются колебания; V - кинематическая вязкость рабочей среды.
При исследованиях не гармонических колебаний потока уравнение движения жидкости необходимо разложить в ряд Фурье и к каждой гармонике применить вычисления аналогичные изложенным для гармонических колебаний. Весь колебательный процесс получиться в результате суммирования параметров распространения составляющих гармоник.
Применительно к расчету характеристик гасителей уравнения (3)...(8) используются для вычисления коэффициентов распространения колебаний (1) как для инерционной трубки, так и для трубки, образующей расширительную полость гасителей. Для каждого расчета учитываются геометрические размеры соответствующего отрезка трубопровода. Соотношения между ними следующие
* = * цт =* к; S = вк/вцт. (9)
В формулах (9) £ цт, £ к, Sцт, Sk - соответственно, протяженности и площади поперечных сечений инерционной трубки и трубки, являющейся расширительной полостью гасителя.
Волновое сопротивление трубопровода с вязкой средой вычисляется по формуле
^ О® ) = Zeol
2.
Л
1 - Ь+
2.
Л
I- ь - \
(10)
где
2в0 =
л/р^в
тр
S
тр
(11)
трубопровода, образующего инерционный
канал гасителей
Упт О® ) • £\ 2вцт О®) • ^\учт О®) • £
^цт () ' £
Гчт( у®:» •£
2 (/®)
вцт ^ '
(15)
Для трубопровода, образующего расширительную камеру устройств
м=
(/®) • £] () • С/'®) • £ ^^[г.(*)•<]
(У®»
. (16)
Матрицы передачи, описывающие динамические свойства дросселей Х^ (1 = 1;2), шунтирующих с обеих сторон расширительную камеру, имеют следующий вид
- волновое сопротивление трубопровода с невязкой средой. Если формула (10) применяется для расчета нестационарного волнового сопротивления инерционной трубки гасителя, то обозначения параметров следующие:
2в 0®)= 2вцт 0®). (12)
Когда вычисляется нестационарное волновое сопротивление расширительной полости формула (10) преобразуется к виду
2в 0®)= 2вк 0®). (13)
При изучении характеристик гасителей необходимо определить зависимости волновых сопротивлений по формуле (10) соответственно для инерционной трубки и трубки, образующей расширительную полость устройств в параметрах этих элементов. В расчетах следует учесть, что волновые сопротивления исследуемых трубок при невязких средах связаны между собой соотношением
2вцт0 = ^2вк0. (14)
Из выражений для коэффициентов распространения и волновых сопротивлений
составим матрицы передачи в форме ||Л||, которые имеют следующий вид. Для отрезка
1 X 0 1
Xi•Z
вцт0
л/У+Т
01
(17)
В формуле (17) параметром Х1 обозначено сопротивление дросселя в безразмерном виде. Безразмерное сопротивление Х1 связывает в единый комплекс параметры реактивных элементов гасителей и сопротивление дросселя. В общем случае сопротивление Х1 комплексное. Например, если сопротивление Х1 инерционное, то имеем гаситель с резонансным контуром.
После формирования зависимостей коэффициентов передачи элементов гасителей, как четырехполюсников необходимо определить коэффициенты передачи в форме ||Л|| в
целом для схемы гасителя. Имея соотношения для коэффициентов передачи в форме
||л||, зависимости собственных характеристик гасителей рассчитываются по формулам.
]
1
Коэффициент собственного затухания гасителя
Кс О) = = V АО) •
. (18)
Волновое сопротивление гасителя со стороны входа
Zc1(jш) =
Л(]ш) • В(]ш)
(19)
С(]ш) • D(jш) и волновое сопротивление со стороны выхода
Zc2(jш) =
D(jш) • В^ш)
(20)
ССш) • Л(jш) •
Безразмерные зависимости для расчета волновых сопротивлений, полученные на базе соотношений (19), (20) имеют следующий вид:
Zci (ш) =
Zci(ш) •Л/8+Г
Zl
(21)
та гидравлических потерь, вычисления проводились при значении вязкости V = 0 . При квазистационарной модели течения колебательного потока коэффициенты к^ и к2 для соответствующих выражений определялись по формулам
2
к1 = к2 =
шR
тр
8v
(22)
вцт Сш)
Аналитические зависимости для характеристик (18)...(21) при значениях § = §н являются громоздкими, а их исследования трудоемкие. Поэтому расчеты характеристик гасителей были проведены с использованием программы "МаШС^", вычисляющей комплексные функции при каскадном и параллельном соединении четырехполюсников по определенному алгоритму. Расчеты собственных характеристик гасителей проводились при вариациях сопротивлений гидродросселей Ri и следующих параметрах устройств и рабочей среды. Радиус поперечного сечения инерционной трубки
Rцт = 0,5см . Протяженность гасителя
* = 1м. Такая длина гасителя выбрана для усиления эффекта влияния продольной распределенности параметров на исследуемые характеристики. Степень расширения гасителя S = 10 . Приведенный модуль упругости рабочей жидкости внутри элементов гасителя Втр = 1470МПа . Вязкость рабочей среды
V = 50 • 10-2 см2/с. В варианте исследования собственных характеристик гасителя без уче-
Результаты расчетов коэффициентов собственного затухания гасителей для модели колебательного движения жидкости без потерь, § = 0 , и для нестационаргой модели гидравлических потерь, § = § н, приведены на рис.За; б; в; г. Графики коэффициентов собственного затухания гасителей без потерь построены точками, а зависимости коэффициентов собственного затухания, рассчитанные с использованием нестационарной модели течения, построены сплошными линиями. Значения шунтирующих сопротивлений
Ri, при которых проводились расчеты, приведены на графиках.
Исходя из значений сопротивлений Ri, применена следующая классификация гасителей. На рис.За графики отражают зависимости Кс "классического" реактивного акустического фильтра низких частот -" АФНЧ". На рис. 3 г зависимости построены для симметричного гасителя с постоянными, активными волновыми сопротивлениями - " ПАС". Два промежуточных графика показывают этапы трансформации зависимостей Кс при изменении параметров Ri. Все графики являются периодическими функциями относительной длины гасителя
* = 1
(23)
где * - протяженность гасителя; % - длина волны подавляемых колебаний.
В окрестностях значений * = 1/4 + и/2
(п = 1;2...) характеристики Кс имеют особые точки. В этих точках реактивные составляющие модулей входных импедансов инерцион-
Рис. 3. Сравнение зависимостей коэффициентов собственного затухания гасителей при нестационарной модели гидравлических потерь (сплошные линии) и при отсутствии потерь в инерционной трубке и расширительной камере
ного канала и расширительном емкости достигают экстремальных значений. У "АФНЧ", при 5 = 0 и относительной протяженности
устройства £ = 1/4 максимальное значение модуля входного импеданса инерционной трубки равно " бесконечности", а минимальное значение модуля входного импеданса расширительной полости равно нулю. Запасаемая расширительной полостью потенциальная энергия равна количеству кинетической энергии, в которую преобразуется в инерционной трубке потенциальная энергия расширительной полости. Этому процессу соответствует волнообразная кривая Кс, проведенная точками на рис.За.
В реальной трубопроводной системе на работу гасителя оказывают влияние гидравлические потери в инерционной трубке и расширительной полости. В инерционной трубке нестационарные потери больше, чем нестационарные потери в расширительной по-
лости устройства. Существующие потери в реактивных элементах гасителей изменяют условия преобразования потенциальной энергии в кинетическую энергию и наоборот. В результате трансформируются собственные характеристики гасителя. При значении
5 = 5 н появляется местное "пиковое" изменение графика зависимости Кс в окрестности £ = 14 + п/2 . Эта зависимость Кс приведена на графике рис.За сплошной линией. При увеличении суммарных гидравлических потерь в расширительной полости, например за счет увеличения сопротивления дросселя
Rl "пик" на зависимости Кс трансформируется в " впадину". Графики Кс, в пределах интервалов частот 0 < £ < 1/2 и т. д. приобретают вид "т - образных" кривых, которые являются специфичными для гасителей с несимметричным шунтированием расшири-
тельной полости. Изменения другого сопротивления дросселя R2, а именно его уменьшение также изменяет характеристики Кс гасителя. При значениях сопротивлений дросселей R1 = R2 = 1, когда реализуется схема симметричного гасителя, зависимость Кс снова приобретает "синусоидальные" вид, имеющий максимумы при относительной частоте * = 1/4 + п/2 . Отличия этой
кривой от зависимостей Кс "АФНЧ" состоят в следующем. Отсутствуют "пиковые" трансформации характеристик. Максимальные значения Кс у симметричного гасителя примерно в два раза меньше чем у гасителя " АФНЧ" при значении § = 0 .
Общие закономерности изменения характеристик гасителей при моделях гидравлических потерь § = 0 и § = § н за исключением характеристик гасителя близкого по параметрам к схеме "АФНЧ" следующие.
Графики Кс отличаются незначительно. До частот колебаний, при которых зависимости Кс достигают максимумов текущие значения
Кс несколько больше при расчетной модели
§ = § н, а для частот, при которых графики Кс изменяются от максимальных значений до минимальных значений, наоборот, рис.3. Зависимости Кс при квазистационарной модели течения жидкости занимают промежуточное положение между зависимостями, построенными при значениях § = 0 и § = §
Характеристики волновых сопротивлений (19).(21), при различных моделях
течения жидкостей в трубопроводных цепях имеют ту же тенденцию изменения, что и
характеристики Кс . В рабочих диапазонах частот затухания характеристики волновых сопротивлений гасителей, построенные при
значении § = §„ незначительно отличаются
XI
от характеристик волновых сопротивлений, построенных при значении § = 0 .
Полученные результаты позволяют сделать следующие выводы.
Влияние нестационарности процессов в гидравлических системах с типовыми акустическими и гидромеханическими параметрами на собственные характеристики гасителей проявляются незначительно. Поэтому для проектных расчетов гасителей целесообразно использовать модель распространения колебаний без потерь. Эта модель существенно проще и позволяют отыскать аналитическим путем базовые оптимальные соотношения между характеристиками гасителей и параметрами его элементов. Сформированные базовые параметры гасителей в последующем могу быть скорректированы по результатам испытаний с учетом нестационарной модели динамических процессов.
СПИСОК ЛИТЕРАТУРЫ
1. Головин А.Н. Разработка гасителей колебаний жидкости для трубопроводных цепей двигателей и систем летательных аппаратов: Диссертация на соиск. учен. степ. канд. техн. наук. Куйбышев, 1983.
2. Попов Д. Н. Динамика и регулирование гидро - и пневмосистем. Учеб. для машиностроительных вузов. М..: Машиностроение, 1976.
ACOUSTIC DAMPER CHARACTERISTICS UNDER NONSTATIONARY
MODEL OF RESISTACE HEAD
© 2004 A.N. Golovin
Samara State Aerospase University
The research method of acoustic damper characteristics under nonstationary model of resistance head ispresented.