УДК 621.65
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДИНАМИКИ РАБОТЫ ЗАПЕРТОГО ОБЪЕМА ШЕСЕРЕННОГО НАСОСА С УЧЕТОМ ДВИЖЕНИЯ ПОДПЯТНИКА
© 2010 Г.О. Белов1, А.Н. Крючков2, Л.В. Родионов1, Е.В. Шахматов1
1 Самарский государственный аэрокосмический университет 2 Институт акустики машин, г. Самара
Поступила в редакцию 7.12.2009
В статье проводится теоретическое исследование влияния геометрии разгрузочных канавок на импульсы давления в запертом объеме и на торцевые утечки. Предложен профиль разгрузочных канавок, обеспечивающий снижение импульсов давления в запертом объеме и торцевых утечек. Ключевые слова: математическая модель, запертый объем, разгрузочные канавки, подпятники, шестеренный насос, давление нагнетания, износ пары трения, диагностика износа, виброакустическое взаимодействие.
Работе шестеренного насоса сопутствует образование запертого объема между двумя парами контактирующих зубьев (рис. 1). Согласно [2, 4, 5, 6, 7], образование и деформации запертого объема приводит к резким импульсным скачкам и провалам давления в этой полости, что значительно сокращает ресурс службы насоса в целом. Зачастую импульсные процессы провала давления приводят к появлению кавитационных каверн, схло-пывание которых способствует выкрашиванию зубьев и увеличению уровня виброакустических колебаний всей системы. В связи с этим недостатком часто ресурс больших сложных машин ограничен ресурсом входящих в их состав шестеренных насосов, однако из-за сравнительно простой конструкции и, как следствие, малой стоимости данный тип насосов получил самое широкое распространение. Для того чтобы снизить пульсации давления в запертом объеме в настоящее время используют два основных метода:
соединение запертого объема с полостью нагнетания и всасывания посредствам специальных каналов в шестернях насоса или их опорах (подпятниках);
разгрузка запертого объема за счет включения в конструкцию разгрузочных полостей.
Несмотря на то, что эти методы исследуются долгое время, единых рекомендаций не суще-Белов Глеб Олегович, аспирант кафедры "Автоматические системы энергетических установок". E-mail: glebbelov@mail.ru
Крючков Александр Николаевич, доктор технических наук, профессор, заместитель директора. E-mail: kan@ssau.ru
Родионов Леонид Валерьевич, кандидат технических наук, ассистент кафедры "Автоматические системы энергетических установок". E-mail: rl63@bk.ru Шахматов Евгений Владимирович, доктор технических наук, профессор, проректор по науке и инновациям. E-mail: iam@ssau.ru
ствует до сих пор. Это происходит вследствие сложности процессов протекающих в запертом объеме. Многие исследователи создавали свои модели [1, 2, 3, 5, 6, 7] для исследования аналогичных систем, однако до сих пор нет общепринятого четкого понимания всех процессов, протекающих при работе шестеренного насоса. Стоит отметить, что для тихоходных насосов существуют, доказавшие свою эффективность, четкие рекомендации по разработке мероприятия для увеличения ресурса за счет разгрузки запертого объема, которые, однако, не пригодны для быстроходных насосов. Для них в условиях производства чаще всего подбирают конфигурацию разгрузочных канавок опытным путем, проводя огромное количество опытов.
В данной же работе для определения эффективных мероприятий использовалась, созданная авторами математическая модель. Эффективность предлагаемых мероприятий оценивалась по размаху импульса давления в запертом объеме. Для того чтобы определить давление в запертом объеме необходимо было определить его величину запертого объема и площади контакта запертого объема с областями нагнетания и всасывания через разгрузочные канавки. Значения площадей контакта разгрузочных канавок определялись численно.
Рис. 1. Схема образования запертого объема
Стоит отметить, что на динамику работы влияли не только разгрузочные канавки, но и передвижения подпятников из за повышения давления в запертом объеме.
Численно была найдена функция запертого объема от угла поворота шестерни (рис. 2). Здесь и далее все зависимости от угла поворота шестерни представлены в диапазоне от 0 (что соответствует моменту образования запертого объема) до 8,5 градусов (что соответствует моменту раскрытия запертого объема).
Давление в запертом объеме можно найти из следующей системы уравнений [3]:
Озо = Осж + Он + Об ; ЖУп
Осж
йг '
Он =±МН Рн.
2|Л^н|.
Р
Об = ±МбКб
2|Лрв|.
Р
1 ЖУп
К
Уп &р
где, Озо - геометрический расход запертого объема; р - давление в запертом объеме; Он -расход в область нагнетания; 0Б - расход в область всасывания; УП - приведенный запертый объем; ¡лн , ¡лБ - коэффициенты расхода площадей контакта запертого объема с областями нагнетания и всасывания соответственно; К , К -
Н Б
площадь контакта запертого объема с полостями нагнетания и всасывания соответственно; Лрн , ЛрБ - разность между давлением в запертом объеме и давлением нагнетания и всасыва-
ния соответственно; р - плотность рабочей жидкости; К - объемный модуль упругости.
В модели приняты следующие допущения: коэффициенты расхода в область нагнетания ¡лн и в область всасывания ¡лБ из запертого объема приняты постоянными.
В представленной модели есть возможность конфигурирования разгрузочных канавок путем изменения следующих пяти параметров (рис. 3): а - угол наклона перемычки; 8НГ - смещение перемычки в сторону области нагнетания; 8ВС - смещение перемычки в сторону области всасывания; ННГ - глубина разгрузочной канавки в области нагнетания; НВС - глубина разгрузочной канавки в области всасывания.
На основе предложенной модели, стало возможно подобрать конфигурацию разгрузочных канавок, обеспечивающую требуемое давление в запертом объеме.
Стоит также отметить, что изменение запертого объема приводит к передвижению подпятника. Схема к расчету перемещения подпятника показана на рис. 4. Из основных уравнений динамики можно определить ускорение подвижного подпятника. Таким образом, в любой момент времени ускорение подпятника может быть определено из следующего выражения:
Уп С т —т- = р- Сх, Ж2 Ъ
где т - масса подпятника; х - перемещение подпятника; С - жесткость системы.
Скорость подпятника вычислялась следующим образом:
йг
Жх1-1 . & х1-1
йг
+
йг2
• Лг, при 0 < х^ < х
'МАХ ~>
- к&х=1
Жг
при
х = 0;
х = X
МАХ'
где Лг - дискретизация по времени; хМАХ - мак-
1
Рис. 3. Параметры, определяющие геометрию разгрузочных канавок
симально допустимое системой перемещение подпятника; к - коэффициент потерь при соударении.
Перемещение подпятника определятся следующим образом:
Лх., , Ох, Аг „
ох, , О'х,, Аг2
!—А +-2т--
Л Л! 2
ОХ , . сРх,, Аг2
х , +—— А +---, при0<х , +—— А! +-'—■--<х,,,.:
0* = (Рн - Рв ):
х3 (вн-вв ) г
где вн и вв - углы соответствующие полости нагнетания и всасывания; г, и г2 - радиусы внутренней и внешней границы торцевого зазора.
На основе предложенной модели была проанализирована эффективность применения большого ряда вариантов разгрузочных канавок,
Ох , , Л2х,, А2 Хмх при 0 < х—, +Л- А ■—.
Согласно [6], образование такого зазора влечет к появлению торцевых утечек. В процессе работы действительные торцовые зазоры являются переменными. Изменения величины зазоров вызываются биением торцов шестерен, пульсацией давления в нагнетательном трубопроводе, отклонениями от плоскостности сопряженных торцов роторов и уплотняющих деталей, неровностями на торцовых поверхностях деталей, компрессией жидкости во впадинах зубьев и упругой деформацией поверхностей скольжения в зоне контактирования. Поэтому расчетная зависимость для потерь через торцовые зазоры должна включать их величины в качестве переменных функций. При всех возможных изменениях торцового зазора его минимальная величина не должна быть меньше значений, гарантирующих сохранение масляной пленки между трущимися поверхностями. Это одно из основных требований, предъявляемых к конструкции шестеренных насосов любой разновидности. В расчетах использовалась следующая формула для определения величины утечек через торцевые неплотности:
Рис. 4. Расчетная схема
для определения смещения подпятника: 1- подшипники; 2 - корпус; 3 - неподвижные подпятники; 4 - зубчатые шестерни; 5 - упругий элемент(прокладка); 6 - шпонка; 7- подвижный подпятник; 8 - валы
подобрана конфигурация, обеспечивающая максимальное снижение импульсов давления в запертом объеме. Геометрия исходного и предложенного профиля разгрузочной канавки представлены на рисунках 5а и 5б. На рис. 6 и 7 показаны пульсации давления в запертом объеме на одном и том же режиме при исходном и предложенном профиле.
По полученным функциям давления от угла поворота шестерни были определены перемещения подпятника. Эти зависимости показаны на рис. 8 и 9.
По зависимостям смещения подпятника от угла поворота шестерни (рис. 8 и 9), а также с учетом постоянной составляющей зазора (для этого насоса 10 мкм), были получены импульсы расхода через торцевые неплотности. Они представлены на рис. 10 и 11.
Разработанная математическая модель позволила выполнить канавки, обеспечивающие снижением амплитуды пульсаций давления в запертом объеме в три раза. В настоящее время как за рубежом [1, 2], так и у нас ведуться разработки активних и адаптивних систем по снижению виброакустической нагруженности ше-стеренних насосов, хотя стоит отметить, что пока ни одна такая система не внедрена в серийное производство.
Еще одним важным результатом модели является получение зависимости торцевых утечек от угла поворота ведущей шестерни. В работе было показано, что порядок торцевых утечек зависит от эффективности работ разгрузочных канавок. Стоит отметить, что большинство специалистов полагают, что торцевые утечки являются преобладающими в шестеренных насосах, а в их уменьшении есть некий запас по увеличению КПД.
х
г
а) б)
Рис. 5. Исходный (а) и предлагаемый (б) профиль разгрузочных канавок
! * 3 I
1 А 111111 -1.ПО -П| ---Ри 1 1
1 ч_ V
пг'т Т П |1 П Г11 Г|Т 1 ТП ГГ' Н 11ГГ-ГТ^'|»ГТ туит
Рис. 6. Рассчитанное давление в запертом объеме при исходном профиле
Рис. 7. Рассчитанное давление в запертом объеме при предложенном профиле
В настоящее время за рубежом [8] активно осваивают методы определения пульсаций утечек по пульсациям давления для последующей диагностики состояния качающего узла. Для этого сперва определяют пульсационную производительность качающего узла в лабораторных условиях, затем в рабочей системе, зная пульса-ционную производительность качающего узла и пульсации давления на выходе из насоса, определяют импеданс системы. В процессе работы, зная импеданс и измеряя пульсации давления,
вычисляют пульсационную производительность. По изменению пульсаций потока определяют пульсации утечек диагностируют состояние насоса. Однако стоит отметить, что импеданс системы после определенного времени работы будет меняться, к тому же пульсации потока, определенные в лабораторных условиях могут отличаться от пульсационной производительность в системе, таким образом этот метод диагностики насосов вряд ли имеет перспективы вследствие, вероятно, малой точности.
V
/ > ч
У / ч \
/
/
/
/
/ /
3 4 5
Посопг шестерни, град
Рис. 8. Смещение подпятника при исходном профиле разгрузочных канавок
ф I
/
/
l1
/
/
J
J / \
/ \
у \
/ \
1
Поверг ш«т»рнн. ФЛД
Рис. 9. Смещение подпятника при предложенном профиле разгрузочных канавок
Несмотря на то, что этот метод малопригоден для диагностики износа всего насоса с помощью него можно диагностировать насколько эффективно работают разгрузочные канавки и насколько изношены пары трения шестерня подпятник. Современные насосы позволяют снижать радиальные утечки и утечки через неплотности контакта зубьев практически до нуля. Торцевые же утечки составляют основную долю утечек насоса. Если считать импеданс насоса не зависящем от конфигурации разгрузочных канавок, то описанным выше методом можно определять изменение утечек через торцевой зазор. Этот метод будет не точен, зато он даст порядок числа, а как было показано выше эффективное снижение заброса давления в запертом объеме в несколько раз дает изменение пульсационных
торцевых утечек на несколько порядков. Таким образом по изменению утечек можно судить об эффективности работы разгрузочных канавок и об износе пары трения шестерня подпятник.
СПИСОК ЛИТЕРАТУРЫ
1. Edge K.A., Lipscombe B.R. The reduction of gear pump pressure ripple // Proc. IMechE. 1987 Vol. 201. No B2.
2. Casoli P., Vacca A., Franzioni G. A numerical model for the simulation of external gear pumps // The six JFPS international symposium on fluid power, Tsukuba, Japan 2005.
3. Шейпак А.А. Гидравлика и гидропневмопривод: Учебное пособие. Ч.1. Основы механики жидкости и газа. М.: МГИУ, 2005.
4. Башта Т.М. Объёмные насосы и гидравлические двигатели гидросистем. М.: Машиностроение, 1974.
5. Чиняев И.А. Роторные насосы. Справочное пособие. Л.: Машиностроение, 1969.
— f /
/
/ /
0 12345678
Угол поворота ведущей ujç(.rç|jMH. грод
Рис. 10. Торцевые утечки при исходном профиле разгрузочных канавок
Г Vi
Рис. 11. Торцевые утечки при предложенном профиле разгрузочных канавок
6. Рыбкин Е.А., Усов А.А. Шестеренные насосы для металлорежущих станков. М.: Машиностроение, 1960.
7. Юдин ЕМ. Шестеренные насосы. М.: Машиностроение, 1964.
8. Yang M, Edge K.A., Johnston D.N. Condition monitoring and fault diagnosis for vane pumps using flow ripple measurement. Fluid Power and Motion Control. 2008.
NUMERICAL MODEL OF THE LOCKED VOLUME WORKING OF GEAR PUMP WITH PROVISION FOR MOTION OF BRIDGES
© 2010 G.O. Belov1, A.N. Kruchkov2, L.V. Rodionov1, E.V.Shakhmatov1
1 Samara State Aerospace University 2 Institute of Acoustic of Machines, Samara
Paper considers the numerical model of the locked volume of gear pump. That model gave us research of influence of geometry of grooves on pressure ripples in locked volume. The profile of grooves provided decrease of pressure pulsations in locked volume. It was designed method of monitoring work of grooves in gear pump. Key words: mathematical model, locked volume, discharge gouges, gear pump, pressure outlet, wear-out pairs friction, diagnostics of the wear-out, hydro-acoustic interaction
Gleb Belov, Graduate Student. E-mail: glebbelov@mail.ru. Alexander Kruchkov, Deputy Director, Doctor of Technics, Associate Professor. E-mail: kan@ssau.ru. Leonid Rodionov, Candidate of Technics, Assistant Lecturer. E-mail: rl63@bk.ru.
Eugene Shakhmatov, Vice Rector for Science and Innovation, Doctor of Technics, Professor. E-mail: iam@ssau.ru