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

УДК 534.322.3

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

© 2006 Л.В. Родионов1, А Н. Крючков2, Е.В. Шахматов1, М.С. Гаспаров1

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

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

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

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

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

Гидравлическая компенсация торцевых зазоров с применением подвижного подпят-

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

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

2 3 7 4

Рис. 1. Конструктивная схема комбинированного насосного агрегата: 1- корпус; 2 - шнек; 3 - центробежное колесо; 4 - шестеренная ступень; 5 - привод; 6 - датчик

На рис. 2 представлен спектр вибрации на фланце подвода топлива при разных давлениях на выходе из шестеренной ступени.

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

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

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

Рис. 2. Спектры вибрации на фланце подвода топлива при различных давлениях

Рис. 3. Схема действия осевых сил на подпятники

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

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

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

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

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

Давление разбивается на три зоны (рис. 4):

1-я - зона полного давления нагнетания Рнагн (у=90°);

2-я - переходная зона от нагнетания к всасыванию, где давление падает от Рнагн до Рвс пропорционально углу (ф=180°);

3 - я зо на давления вс асывания Рвс (Х=90°);

Каждая зона, в свою очередь, разбивается на два участка:

1-й участок, ограниченный радиусами головок Rе и впадин Ri, т.е. участок, охваты-

Таблица 1. Основные геометрические и режимные параметры качающего шестеренного узла

№ п/п Геометрические параметры Значение

1 Число зубьев г 11

2 Модуль зацепления т, мм 6

3 Зазор по спинкам зубьев, мм 0,5

4 Ширина шестерни Ь, мм 21

5 Диаметр окружности головок De , мм 84

6 Диаметр начальной окружности D, мм 72

7 Диаметр делительной окружности DдЕЛ, мм 66

8 Диаметр основной окружности Do, мм 62

9 Диаметр окружности впадин Di, мм 57,9

10 Межцентровое расстояние А, мм 72

11 Угол зацепления а , град 30,527

12 Угол зацепления по вершинам зубьев ае, град 42,4

13 Угол дуги по начальной окружности фгеом, град 15,966

14 Высота зуба ^ мм 13

15 Шаг зацепления по основной окружности t0, мм 17,7

16 Толщина зуба у вершины, мм 1,9

17 Коэффициент перекрытия е 1,1338

18 Радиальный зазор в запертом объеме, мм 1,04

19 Минимальный радиус контакта г2, мм 20

20 Угол зоны полного давления нагнетания у, град 90

21 Угол переходной зоны от нагнетания к всасыванию ф, град 180

22 Угол зоны давления всасывания х, град 90

23 Давление нагнетания Рнагн, МПа 10

24 Давление всасывания Рвс, МПа 1

25 Частота вращения п, об/мин 4800

26 Круговая частота ю, 1/сек 502,6

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

2-й участок, ограниченный радиусами впадин Ш и минимальным радиусом контакта г2, на котором давление падает по логарифмическому закону.

Задаем P

Сила давления жидкости на первом участке 1-й зоны [1]:

нагн

2

(Яе2 - Щ2 ).

Сила давления жидкости на втором участке 1-й зоны:

F2 =

нагн

2

т1

т2 - г22

2  1п

т

\ г2 уу

.

1

Рис. 4. Схема разбиения шестерен на расчетные зоны

Значение суммарной силы для 3-й зоны определяется по формуле:

ПП = ^ + ^.

Полное значение силы определяется (рис. 5):

Рполн = Р1 + РП + ¥Ш.

Угол запертого объема рассчитывается по формуле:

2 ти

Рzap =--

z

Площадь одного рассматриваемого элемента:

Значение суммарной силы для 1-й зоны определяется по формуле:

Р1 = Р + Р2

Яе2- т2

Р30) = Рнагн -Ф)--

4

Р4(() =

нагн

ф

4

Я12 -

Я12 - г22

2  1п

гя^ v г2 уу

т(() = ¥3(() + ¥4(().

Р5 =

вс

2

Х (Яе2 - Я12).

Рвс X

2

Ш2 -

Яi2 - г22

V г2 У

2 1п

5 =

Сила давления жидкости на первом участке 2-й зоны:

где угол одного межзубового объема Р?дк = 360

.

11

Сила давления жидкости на втором участке 2-й зоны:

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

Для первого зуба (рис. 6):

Значение суммарной силы для 2-й зоны определяется по формуле:

РО) =

(РНАГН РВС ) ' 5

-Р(1) + Рвс- 5

(рвс-5- РР) • 2 (РВС - 5 -РР ) '2 . .

Рг

Рг

Сила давления жидкости на первом участке 3-й зоны:

Сила давления жидкости на втором участке 3-й зоны:

0.002 0.004 0.006 0.008 0.010 0.012

Время, с

Рис. 5. График полной силы

<Р

Р...„. - 5

РВС '5

х 10

1.9

0.9

0

6

Ч

и

Ч

и

0.002 0.004 0.006 0.008

Время, с

Рис. 6. Осевая сила первого зуба, действующая на неподвижный подпятник

k  Р  •

п нагн ° :

где Fp

предполагаемый скачок осевой силы в момент образования запертого объема; к=4. Для второго зуба (рис. 7):

А 2(0 = р(0 + А.

где межзубовый угол А

360

f2(t) =

^РНАГН РВС )  S

9

р  я

НАГН

(рвс ' s - f|¡) 2

р20) + рвс ' s

Рилр гвс 'S

(РНАГН — РВС )  S

9

«>' 2 (рвс  s — fй ) 2

^--р2(t) + Рр — ^-рр—

р20) + рвс ' s —

Рялг

(РНАГН — РВС S

'

Время, с

Рис. 7. Осевая сила второго зуба, действующая на неподвижный подпятник

11

Fsum(t) = 2Ъ(*) .

I=1

Аналогично для остальных зубьев смещение на Д.

Суммарная осевая сила, действующая на неподвижный подпятник (рис. 8):

Для подвижного подпятника вводим учет временной задержки Дt=0,000284с

t1

Рs

Осевая сила, действующая на подвижный подпятник, описывается следующей зависимостью:

Для первого зуба (рис. 9):

РФ) =

((РНАГН РВС )' S

+Рвс' S

Р 2ЛР

Рвс '

Р,

К

(Й ч к о

0.002 0.004

0.006 0.008 0.010 0.012

Время, с

Рис. 8. Суммарная осевая сила, действующая на неподвижный подпятник

К

ев Ч К О

0 0.002 0.004 0.006 0.008 0.010 0.012

Время, с

Рис. 9. Осевая сила первого зуба, действующая на подвижный подпятник

0

0

0

0.002

0.004

0.006

0.008

0.010

0.012

0

0.010

0.012

9

9

Р,..... •

X 10

3.8

3.6

3

2.8

2.6

0

0

К

ей Ч К О

Время, с

Рис. 10. Осевая сила второго зуба, действующая на подвижный подпятник

Для второго зуба (рис. 10): Р2s(t) = Рs(t) + А,

f2s(t) =

^РНАГН РВС ) ' Я '

9

р  я

1 НАГН °

(рвс ' я — Рр )' 2 '

Р1ЛР РВС ' я

(РНАГН — РВС )' я '

. 9

р2s(t) + рвс ' я

р>' 2

р2s(t) + рвс  я —

Рхлр

(РНАГН — РВС ) ' Я

[9 + у + х)

Fnр-Fл

к

Й н к о

0 0.002 0.004 0.006 0.008 0.010 0.012

Время, с

Рис. 12. Результирующая осевая сила

Рпр-Рл

3000

x

ев

И 2500

СО

Л

1- ев 2000

я

«

ев 1500

и

о

<Ц

(Г 1000

8

2

ев

я 500

я

«

А_

О 1000 2000 3000 *

Частота, Гц Рис. 13. Спектр результирующей осевой силы

К

ей Ч К

о

/

0.002 0.004

0.006 0.008

Время, с

Рис. 11. Суммарная осевая сила, действующая на подвижный подпятник

Аналогично для остальных зубьев смещение на Д.

Суммарная осевая сила, действующая на подвижный подпятник (рис. 11):

11

Fsums

I =1

Результирующая осевая сила, действующая в шестеренной ступени насосного агрегата равна

На рис. 13 представлен спектр результирующей осевой силы.

В ходе измерений установлено, что основной вклад в спектре осевой вибрации на фланце подвода топлива и фланце соединения ШЦС и ШС вносят составляющие с час-

1000 2000 3000 4000

Частота, Гц Рис. 14. Амплитудный спектр осевой вибрации на фланце соединения ШЦС и ШС при угле РУД = 400

X 10

3.8

3.6

3.4

3.2

2.8

2.6

2.4

2

0.002

0.004

0.006

0.008

0.010

0.012

0. 010

0. 012

9

X 10

1.5

0.5

0

-0.5

-1

тотой свыше 3 кГц (рис. 14).

Полученные графики свидетельствуют о том, что данная математическая модель может применятся в диапазоне частот до 4 кГц.

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

1. Юдин Е. М. Шестеренные насосы. М.: Машиностроение, 1964.

DEVELOPMENT OF MATHEMATICAL MODEL OF EXCITATION OF HIGH-FREQUENCY AXIAL VIBRATIONS OF THE CASE OF COMBINED WITH THE SWINGING UNIT GEAR

© 2006 L.V. Rodionov1, A.N. Kruchkov2, E.V.Shakhmatov1, M.S.Gasparov1

1 Samara State Aerospace University 2 Institute of Acoustic of Machines

In the article the reasons increased axial vibroloads the combined pump including pumping up screw-centrifugal and basic gear a step are researched. The developed mathematical model is presented. It is carried out calculation of axial forces with the purpose of research of influence of axial forces on bearing unit of the pump, operating in gear step.