МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ MACHINE BUILDING AND MACHINE SCIENCE
Ф
и Ö О
тз
л
'S
s (U
ü £ л
Ц) Check for updates
Научная статья
УДК 62-19
https://doi.org/10.23947/2687-1653-2022-22-2-130-141
Определение динамических характеристик шестерённого насоса по методу вариации нагрузок с использованием специальных стендовых систем
В. И. Санчугов , П. Д. Рекадае Ш
Самарский национальный исследовательский университет имени академика С. П. Королева, Российская Федерация, г. Самара, ул. Московское шоссе, 34 И Rekadze 1993@gmail.com
Аннотация
Введение. Пульсации рабочей среды в гидросистемах станков, системах топливоподачи авиационных двигателей и подачи жидкостного ракетного двигателя являются одним из главных факторов возникновения вибрации агрегатов, приводящей к усталостному разрушению корпусов элементов, трубопроводов и выходу из строя элементов насосов. Целью работы являлись реализация метода определения динамических характеристик объёмного насоса с использованием специальных стендовых систем и сравнение результатов расчёта с экспериментальными данными. Описаны этапы расчёта динамических характеристик объёмного насоса с учётом разработанных ранее специальных стендовых систем на примере шестерённого насоса внешнего зацепления производительностью 14 см3/об. Показана реализация метода вариации нагрузок В. П. Шорина с использованием разработанных авторами специальных стендовых систем с заранее определёнными динамическими характеристиками. Излагаются основные этапы методики определения динамических характеристик шестеренного насоса.
Материалы и методы. В работе применены методы спектрального анализа пульсирующего давления. Пульсации расхода жидкости на выходе из насоса определялись с помощью импедансного метода, метода вариации нагрузок и специальных стендовых систем.
Результаты исследования. В работе реализована методика определения динамических характеристик шестеренного насоса в диапазоне оборотов приводного вала 500-2500 об/мин для четырёх гармонических составляющих спектра колебаний в широком диапазоне динамических нагрузок (от инерционного до емкостного характера). Представлен анализ стендовых систем, приводящих к расчёту динамических характеристик насоса с минимальной погрешностью на основании условия согласования динамической нагрузки и источника колебаний. Проведена проверка разработанного подхода к оценке динамических характеристик насоса посредством сравнения расчётных и экспериментальных данных пульсаций давления в стендовых системах с дросселем, полостью и протяжённым трубопроводом на выходе насоса. Обсуждение и заключения. Реализован метод определения динамических характеристик объёмного насоса с использованием разработанных авторами специальных стендовых систем. Результаты исследований показывают, что исследуемый шестерённый насос может рассматриваться как независимый источник колебаний расхода, для которого отклонение собственных динамических характеристик от средних значений не превышает 10 % по первой гармонической составляющей.
Ключевые слова: объёмный насос, метод вариации нагрузок, стендовые системы, источник колебаний расхода, динамические характеристики насоса, независимый источник колебаний, пульсации расхода, импеданс. Финансирование: работа второго автора выполнена по теме № FSSS-2020-0015 в рамках выполнения госзадания Минобрнауки России.
© Санчугов В. И., Рекадзе П. Д., 2022
Для цитирования: Санчугов, В. И. Определение динамических характеристик шестерённого насоса по методу вариации нагрузок с использованием специальных стендовых систем / В. И. Санчугов, П. Д. Рекадзе // Advanced Engineering Research. — 2022. — Т. 22, № 2. — С. 130-141. https://doi.org/10.23947/2687-1653-2022-22-2-130-141
Original article
Determination of the dynamic characteristics of a gear pump by the load variation method using special bench systems
Valeriy I. Sanchugov , Pavel D. Rekadze ED
Samara National Research University, 34, Moskovskoye Shosse,Samara, Russian Federation И Rekadze 1993@gmail.com
Abstract
Introduction. Some of the major factors in the occurrence of vibration of units, causing fatigue failure of the housings of elements, pipelines, and failure of pump elements, are pulsations of the working medium in the hydraulic systems of machine tools, fuel feed systems of aircraft engines and liquid-propellant engine supply. This study aimed at the implementation of a method for determining the dynamic characteristics of a volumetric pump using special bench systems, and the comparison of the calculation results to the experimental data. The stages of calculating the dynamic characteristics of a volumetric pump were described, taking into account the pre-developed special bench systems on the example of an external gear pump with a capacity of 14 cm3/rev. The implementation of V. P. Shorin's load variation method using special bench systems developed by the authors with predetermined dynamic characteristics was shown. The main stages of the methodology for determining the dynamic characteristics of a gear pump were described. Materials and Methods. Methods of spectral analysis of pulsating pressure were applied in the work. Pulsations of fluid flow at the pump outlet were determined using the impedance method, the method of load variation, and special bench systems.
Results. The paper implemented a technique for determining the dynamic characteristics of a gear pump in the drive shaft speed range of 500-2500 rpm for four harmonic components of the vibration spectrum in a wide range of dynamic loads (from inertial to capacitive nature). The bench systems yielding the calculation of the dynamic characteristics of the pump with a minimum error based on the condition of matching the dynamic load and the source of vibrations were analyzed. The developed approach to the evaluation of the dynamic characteristics of the pump was verified through comparing the calculated and experimental data of pressure pulsations in the bench systems with choke, cavity and an extended pipeline at the pump outlet.
Discussion and Conclusions. The method for determining the dynamic characteristics of a volumetric pump was implemented using special bench systems developed by the authors. The research results show that the gear pump under study can be considered as an independent source of flow fluctuations, for which the deviation of its own dynamic characteristics from the average values does not exceed 10% for the first harmonic component.
Keywords: volumetric pump, load variation method, bench systems, source of flow fluctuations, pump dynamic § characteristics, independent source of fluctuations, flow pulsations, impedance.
<u
Funding information: the research of the second author was carried out on theme no. FSSS-2020-0015 within the g
frame of the government task of the Ministry of Education and Science of the Russian Federation. §
For citation: V. I. Sanchugov, P. D. Rekadze. Determination of the dynamic characteristics of a gear pump by the load Э
variation method using special bench systems. Advanced Engineering Research, 2022, vol. 22, no. 2, pp. 130-141. S
https://doi.org/10.23947/2687-1653-2022-22-2-130-141 ^
К
Введение. Основные недостатки насосов заключаются в сопровождаемом их работу акустическом шуме, вибрации корпуса, присоединенных трубопроводов и возникающих вследствие этого разрушениях конструкции. ср
Традиционный путь совершенствования конструкций насосов, используемый в машиностроении и g заключающийся в сглаживании неравномерности течения жидкости, затруднен отсутствием средств прямого 55 измерения расхода жидкости на выходе насоса. Это ограничивает любую экспериментальную оценку Э эффективности конструктивно-доводочных работ по насосам. Кроме того, это ограничивает аналитическое описание насоса как источника колебаний, необходимое разработчикам гидросистем для обеспечения требуемого динамического качества самих систем.
По этой причине исследования, связанные с моделированием динамических процессов в гидромеханических системах [1, 2] и применением новых подходов к решению задач динамики гидравлических систем [3, 4] на основе физических принципов, имеют актуальность. 131
Материалы и методы. Представлена реализация метода вариации нагрузок В. П. Шорина [4], заключающаяся в том, что определяются пульсации давления за насосом, взаимодействующим с несколькими стендовыми системами, динамические характеристики которых рассчитаны заранее. Далее обработанные амплитудно-частотные характеристики давления с использованием импеданса стендовых систем пересчитываются в собственные динамические характеристики насоса — пульсационную производительность и импеданс. Развитие подхода состоит в том, что предложены специальные стендовые системы, учитывающие соединительные штуцеры, переходники и внутренние каналы насоса, что позволяет точнее определять динамические характеристики насосов, а также осуществлять проверку применимости используемой модели насоса в виде эквивалентного источника колебаний.
Результаты исследования
Методика определения динамических характеристик шестеренного насоса
В Самарском университете разработана расчетно-экспериментальная методика определения динамических характеристик шестеренного насоса (пульсационной производительности и импеданса), которая базируется на исследованиях Л. Тевенена, Э. Нортона и В. П. Шорина [5, 6]. Разработанная методика включает в себя следующие основные этапы: создание моделей и расчет динамических характеристик специальных стендовых систем; экспериментальное определение пульсаций давления на выходе насоса; расчет спектра возбуждаемых колебаний в различных стендовых системах; определение аппроксимирующих зависимостей для отдельных гармонических составляющих; расчёт собственных динамических характеристик насоса (пульсационной производительности и импеданса); оценка независимости (стабильности) динамических характеристик насоса от характеристик стендовых систем; проверка применимости используемой модели насоса в виде эквивалентного источника колебаний расхода (или давления). Далее приведен пример реализации предложенной методики с использованием шестерённого насоса внешнего зацепления удельной производительностью 14 см3/об.
Модели и расчет динамических характеристик специальных стендовых систем
Рассмотрены базовые варианты стендовых систем на выходе жидкости при испытаниях насосов. В сосредоточенных параметрах стендовая система может реализовывать активную расходно-зависимую, инерционную и упругую нагрузки. В распределённых параметрах стендовая система может быть реализована в виде протяжённого цилиндрического трубопровода, реализующего активную частотно-независимую неотражающую нагрузку.
Модели насоса в виде эквивалентных источников колебаний расхода (ЭИКР) или давления (ЭИКД) определяются с использованием источника постоянного расхода О, идеального источника колебаний давления р (или расхода д) и внутреннего импеданса [4]. На рис. 1 представлены модели насоса в виде эквивалентного источника колебаний давления (и расхода) с присоединенной системой.
и Й О
тз
Л 'й
и
£ л
а) б)
Рис. 1. Модели насоса в виде эквивалентного источника колебаний: а — давления; б — расхода с присоединенной системой: О,Р0 — средний расход и давление; рИ— собственная пульсационная производительность насоса по
давлению и расходу, 2п, 2ст с — импеданс насоса и стендовой системы (адаптировано из [4])
При этом задачей выбора стендовой системы является обеспечение таких параметров системы, при которых реализуются наилучшие условия высокой точности измерений. Задачу выбора параметров присоединенной системы необходимо проводить при изменении импеданса в диапазоне от -да до +да. Принципиальные схемы стендовых систем представлены на рис. 2.
Рис. 2. Принципиальные схемы стендовых систем для определения собственных динамических характеристик насосов: а — с активной дроссельной нагрузкой; б — с активной неотражающей нагрузкой; в — емкостной нагрузкой: 1 — испытуемый насос; 2 — нагрузочный дроссель; 3 — расходомер; 4 — датчик давления; 5, 6 — соединительные штуцеры; 7 — магистраль; 8 — протяжённый трубопровод; 9 — полость
Расчётные модели импедансов специальных стендовых систем: - с активной дроссельной нагрузкой (дроссель):
Тст.с.к (00 ' /)
(1др (00) + КеТтр (00 ))2 +
\2
Рк
(1)
£
Ъдр (0о )+Яе2тр (0о )'
тр
- с активной распределенной нагрузкой (протяжённый трубопровод):
ра
- с емкостной нагрузкой (полость):
Тст.с.к (0О ' /к ) =
1
]атс,%
(КеТам (00)+ Яе2тУрб (0о) + Ке1мс (00) + Тр (0О))
(2ж/к )2 1С-1_
(2я/к )21С -1 2ж/кС
(2)
(3)
2ж/кс(кетаам + кеттурб (/к)+яеги.с.(/к)+г^ (0о )) е
В формулах обозначены: ^ (0О) — импеданс дросселя; (0О ) — действительная часть импеданса
трубопровода, включающего соединительную магистраль, штуцеры, внутренние каналы в насосе и дросселе; ЯеТ^ (0О ) (0О ) ,Яе2шс (0О ) — действительная часть импеданса штуцеров и внутреннего канала в
насосе и полости при ламинарном, турбулентном течении жидкости за насосом и местных сопротивлениях соответственно; / — частота к - й гармонической составляющей спектра колебаний за насосом; р — плотность рабочей среды; /Е — длина трубопровода, включающая длину магистрали 6, соединительных штуцеров 5 и 7 и внутренних каналов насоса и дросселя (не показаны на рис. 7); £ — площадь поперечного
сечения трубопровода; Ь — «инерционность» соединительных штуцеров и внутреннего канала в насосе и полости; С — упругость полости; к — номер гармонической составляющей спектра колебаний; ] — мнимая
единица (] = >/-1).
и К X
<и «
и ю о X
к 3
и
к
X
<и О
Н О
о X X
3
133
х
2
С использованием стендовой системы с дросселем на выходе насоса реализованы три инерционно-активные динамические нагрузки, в которых изменялось активное сопротивление изменением площади перекрытия дросселя Б0:
- для стендовой системы №1 площадь перекрытия F01=9,290x 10-6 м2,
- для стендовой системы №2 площадь перекрытия F02=7,299x 10-6 м2,
- для стендовой системы №3 площадь перекрытия F03=2,247x 10-5 м2.
С использованием стендовой системы с полостью на выходе насоса объёмом 2 л реализована инерционно-упругая динамическая нагрузка.
С использованием стендовой системы с протяжённым трубопроводом длиной 106 м реализована активная частотно-независимая нагрузка.
По моделям стендовых систем по формулам (1)-(3) рассчитаны зависимости импеданса (модуля и фазы) от частоты. В таблице 1 для 1 -й гармоники приведены функции регрессии модуля импеданса стендовой системы и фаз Р^ для анализируемых стендовых систем.
Таблица 1
Функции регрессии динамических характеристик стендовых систем по 1-й гармонике
№ Модуль импеданса
стендовой Тип стендовой стендовой системы Фаза импеданса стендовой системы
системы по рисунку 2 системы кг 2 (1), - 1 4 4 2 м с в да, рад
а - №1 С дросселем на выходе насоса 1,41210^ + 2,201 • 109 7,066 10-^ + 0,199
а - №2 1,61010'^ + 4,258 109 6,906 10-^ +0,09039
а - №3 4,980 10'^ + 1,503 • 109 2,723^0-^ + 0,251
б - №4 С полостью на выходе насоса 3,308 10¥ - 1,41610^ + 1,873 • 109 ,4,234 -106 Г + 1,495 -10% атс1я\ ) 1,43 0-105 Г + 2,930-106
в - №5 С протяжённым трубопроводом на выходе насоса 1,897 1010 0
и Й О
тз
л 'й
и
£ Л
Экспериментальное определение пульсации давления на выходе насоса
Исходной информацией для определения динамических характеристик насоса являлось измерение давления на выходе насоса. Пример осциллограммы и спектра пульсаций давления на выходе из насоса в стендовой системе с активной распределенной нагрузкой представлены на рис. 3. Режим работы насоса при этом: среднее давление Р = 23,8 МПа, частота вращения приводного вала п = 1000об / мин .
24,30 -\ =-:-■=:-;:■":-
н
а
0,03 0,04 0,05 0,06 0,07 0,08 Время поворот ротора, с
60 120 180 240
Угол поворота ротора, град а)
300
360
0
н
S3
й <
0,42 0,35 0,28 0,21 0,14 0,07 0,00 1
2
4 I5 6 7 8
167 334 501 668 835 1002 1169 1336 1503 1670 1837 Частота составляющих спектра колебаний, Гц б)
Рис. 3. Давление на выходе насоса: а — осциллограмма; б — спектрограмма
Амплитудно-частотная характеристика пульсаций давления на выходе насоса в диапазоне оборотов приводного вала насоса n=500-1250 об/мин представлена на рис. 4. Сплошной линией нанесены функции регрессии, построенные в среде Microsoft Excel при помощи инструмента «Линия тренда».
^ 0,6
е
л
лва
д
0,5 0,4 0,3 0,2
1 гармоника
2 гармоника
3 гармоника
.4 гармоника
Г 0,1
<
83 166 249 332 415 498 581
Частота составляющих спектра колебаний 13, Гц
I-1-1-1
0 500 1000 1500
Частота вращения приводного вала п, об/мин Т"
664 747 830
X
0 7,3 14,6 21,6
Средний расход жидкости насоса <3, л/мин
Рис. 4. Амплитудно-частотная характеристика пульсаций давления на выходе насоса по отдельным гармоническим составляющим (точками обозначены данные на: О — первой, О — второй, О — третьей, • — четвёртой гармониках,
линиями обозначены функции регрессии)
Анализ осциллограммы, спектрограммы и амплитудно-частотной характеристики пульсаций давления показал следующее: тип процесса полигармонический установившийся; число регистрируемых гармоник — 8, анализируемых1 — 4; диапазон амплитуд колебаний давления 0,28-0,52 МПа; частотный диапазон приводного вала 500-1250 об/мин, анализируемых гармонических составляющих спектра 83-833 Гц.
Для гармонического анализа периодических сигналов использовалось разложение в ряд Фурье [7]. Для разложения в ряд Фурье использованы участки осциллограммы давления, измеряемые относительно опорного сигнала на роторной частоте, т.е. за полный оборот ротора извлекался один участок осциллограммы. Опорный сигнал регистрировался при фиксированном положении шестерни, выставляемом при начале записи осциллограммы давления для пяти стендовых систем. Таким образом определялись амплитудно-частотные
(АЧХ) А (®) и фазо-частотные (ФЧХ) ^ (ю) характеристики.
1
<и К X <и
<и и о X
к 3
eö
<и К X <и О
а н о о X X
3
CÖ
1 В работе принято, что крайняя высшая анализируемая гармоника по амплитуде не должна быть более 5 % от амплитуды первой гармоники.
0
Расчет спектра возбуждаемых колебаний в различных стендовых системах и определение аппроксимирующих зависимостей для отдельных гармонических составляющих колебаний давления
На основе экспериментального определения давления в пяти стендовых системах: с дросселем (№1-3), с полостью (№4) и протяжённым трубопроводом (№5) на выходе насоса рассчитаны спектры возбуждаемых колебаний. В таблице 2 приведены функции регрессии амплитуд Арф и начальных фаз пульсаций давления для 1-й гармоники.
Таблица 2
Функции регрессии амплитуд и начальных фаз пульсаций давления
№ стендовой системы по рисунку 2 Тип стендовой системы Амплитуда пульсаций давления кг Api (f), — мс Начальная фаза пульсаций давления ф рад Частотный диапазон
а - №1 С дросселем на выходе насоса 1,003 10-3f -0,543 83-417
а - №2 1,243-10"Jf -0,280
а - №3 6,129- 10-4f -5,941
б - №4 С полостью на выходе насоса 1,343 10-6f2 - 4,298 10-4f + 3,91310-2 -5,373
в - №5 С протяжённым трубопроводом на выходе насоса 2,694 10-3f -0,658
Расчёт динамических характеристик насоса
Для расчета динамических характеристик насоса (импеданса и пульсационной производительности источника) необходимо провести эксперимент как минимум с двумя динамическими нагрузками2. При этом импеданс насоса по уточнённой авторами модели эквивалентного источника колебаний расхода с использованием аналитических зависимостей р. (/) и 7,. (/) из таблиц 1 и 2 рассчитывается по формуле:
1 - А
7 - 21-ртк ■ <4)
Р2 71
где 7а — импеданс насоса (7и — ) — модуль импеданс насоса (|7И| Яе7И2 + 1ш7И2 ),Яе7и,1т7и —
действительная и мнимая часть импеданс насоса, рш — аргумент импеданса насоса, 7Х, 72 — входной импеданс динамических нагрузок № 1-2 (71 = \7\е Д, 72 =|72|еД), Д, Д — аргументы импедансов динамических нагрузок № 1-2, р1, р2 — пульсации давления за насосом в стендовой системе № 1-2 ( р = Лр1еР , р = Лр2еуР), Лр1, А 2 — амплитуды пульсаций расхода за насосом в стендовой системе № 1-2, рх ,р2 — начальные фазы пульсаций давления динамических нагрузок № 1-2.
Затем с использованием аналитических зависимостей из таблиц 1, 2 определяется пульсационная производительность насоса по уточнённой авторами модели эквивалентного источника колебаний (переменная составляющая объёмного расхода да и давления ра):
7 7
1 - ^ 1 - ^
р а — р1-— , р —-, (5)
| Чи 721_71 , (5)
й р1 р2 71
тз где а, ри — переменная составляющая объёмного расхода (ди — А еР ) и давления (ри — А^еРр) насоса, •д Аи ,Аи — амплитуда пульсаций расхода и давления насоса, рч ,рр — аргумент переменной составляющей объёмного расхода и давления насоса.
Из
рассчитанных характеристик по модели ЭИКР (аи (/)) и ЭИКД (ри (/)) по формуле (5) только одна
будет отличаться стабильностью, т.к. эквивалентный источник колебаний имеет независимую динамическую характеристику только по одному из параметров (пульсирующий расход или давление).
^^ 2 Для достаточного объёма выборки необходимо использование хотя бы четырёх динамических нагрузок 7^ , формирующих шесть неповторяющихся (в том числе / Ф ]') комбинаций: 712 , 713, 714, 723, 724 , 734 .
Расчет импеданса 1 и переменной составляющей объёмного расхода д насоса по формулам (4)-(5) проведён для пяти динамических нагрузок, т.е. с использованием десяти различных комбинаций этих нагрузок: «1_2», «1_3», «2_3», «1_4», «2_4», «3_4», «1_5», «2_5», «3_5», «4_5» (см. таблицу 1).
Для дальнейшего расчёта и анализа из десяти кривых амплитуд пульсаций расхода получена характеристика амплитуд пульсаций расхода насоса А . Значения амплитуд А (на каждой частоте) рассчитывались по формуле из источника [8]:
А = А + ДА ,
ди ди ди'
где Ади — средняя амплитуда пульсирующего расхода, ДА — доверительный интервал амплитуды пульсаций расхода насоса (ДАди = ), а — среднеквадратическое отклонение амплитуды пульсирующего расхода от
среднего значения, ш — число комбинаций динамических нагрузок (ш = 10), г — коэффициент Стъюдента ( г = 2,262, рассчитан для доверительной вероятности 0,95).
Доверительный интервал амплитуд ДА формируется из погрешностей (инструментальной и методической) и степени зависимости динамической характеристики насоса Ад (/) от стендовой системы. Инструментальная
погрешность состоит из погрешности определения давления и фаз, методическая погрешность состоит из погрешности неучтённого влияния динамических свойств приводной механической системы.
Фазы пульсаций расхода р и их отклонения рассчитывались аналогичным образом. При расчёте фаз комбинации «1_2» и «1_5» исключены по причине вырождения переменной составляющей объёмного расхода насоса д при использовании близких значений характеристик на одной частоте: р1хр1иР1& Р2.
Результаты расчёта собственной динамической характеристики насоса в виде пульсационной производительности и импеданса представлены на рис. 5. Фазы для удобства представлена в градусах.
«
к .у
7,0Е-05 6,0Е-05 5,0Е-05
Е? 8 4,0Е-05
а
й
£ д. 3,0Е-05
I 8 2,0Е-05 К О
§ Й 1,0Е-05 э л
0,0Е+00
й «
о
8 й Л
о М-
Й 3 & &
е Й
& Я и1 °
и к
Е?
й п й
О
83 264 445 626 807 988 1169 1350 1531 Частота гармонических составляющих спектра колебаний, Гц а)
90 60 30 0 -30 -60 -90
83 264 445 626 807 988 1169 1350 1531 Частота гармонических составляющих спектра колебаний, Гц б)
и К X
<и «
и и о X
к
3
й
и К X <и О
Н О
о X
к
3
й
и К
1,6Е+10 1,4Е+10 1,2Е+10 1,0Е+10 8,0Е+09 6,0Е+09 4,0Е+09 2,0Е+09 1,0Е+00
83 264 445 626 807 988 1169 1350 1531 Частота гармонических составляющих спектра колебаний, Гц в)
св о о о св
к
св
° _г
св Й о 1-н
Св
е
90 60 30 0 -30 -60 -90
и Й О
тз
л 'й
и
ср £ Л
83 264 445 626 807 988 1169 1350 1531 Частота гармонических составляющих спектра колебаний, Гц
г)
Рис. 5. Динамические характеристики насоса: а — амплитуда пульсаций расхода; б — фаза пульсаций расхода; в — модуль импеданса; г — фаза импеданса (сплошными линиями обозначены данные по: — первой;
— второй; — третьей; — четвёртой гармониках)
Амплитуда пульсаций расхода на 1-й гармонике имеет монотонно возрастающий характер. С увеличением номера гармоники наклон кривой амплитуд пульсаций расхода снижается до практически постоянного значения к 4-й гармонике. Значения амплитуд Л на высших гармониках на порядок меньше,
чем на 1 -й гармонике.
Модуль импеданса насоса на 1-й гармонике имеет монотонно возрастающий характер.
Для оценки реализованных динамических нагрузок использован критерий согласования модуля импеданса |7и| источника колебаний и входного модуля импеданса динамической нагрузки |7стс| [99, 10]. Модули
импедансов должны быть сравнимы (|7и| = \7тс |), чтобы количественные значения динамической характеристики насоса определялись с минимальной погрешностью. Для сравнения использовано отношение
= 1 динамическая нагрузка считается согласованной с источником.
Ввиду различия форм кривых модулей импеданса насоса и динамических нагрузок в диапазоне частот 83-
7„
7 7
и . При и
7 ст.с. 7 ст.с.
417 Гц рационально использовать среднее отношение
7,,
которое позволит учитывать вклад большей
части частотного диапазона. Т.е. отношение
7„
7 „
по 1 -й гармонике составит:
- для активной дроссельной нагрузки №1 — 1,16;
- для активной дроссельной нагрузки №2 — 0,8;
- для активной дроссельной нагрузки №3 — 2,65;
- для емкостной нагрузки №4 — 4,64;
- для активной распределенной нагрузки №5 — 0,35.
Т. е. динамические нагрузки стендовых систем №1, 2 являются наиболее согласованными с источником колебаний. Таким образом, динамическая характеристика насоса, рассчитанная с использованием динамических нагрузок стендовых систем №1, 2 («1_3», «2_3», «1_4», «2_4», «1_5», «2_5») определена с минимальной погрешностью.
Оценка стабильности динамических характеристик насоса от характеристик стендовых систем Стабильность динамических характеристик А и ^ (и <ри) насоса оценивалась с использованием инструментов математической статистики [8]. Так стабильность характеристики Л оценена через расчёт
отношения доверительного интервала ДА на частоте / к среднему значению характеристики 1АдиМ1 на этой
частоте ' Полученное относительное отклонение Д не должно превышать заданное отклонение Д
дощ '
\ =
ДА
ди
А
диМе
< Ддопд,
где \АсиМ1е — среднее медианное значение амплитуды пульсаций расхода насоса; Дйоид — задаваемое относительное отклонение амплитуды пульсаций расхода насоса от среднего значения.
При относительном отклонении Д ^ 0 (когда |ДА^| ^ 0) стабильность характеристики А?и максимальна. Исходя из опыта других исследователей, предпочтительно выбирать величину Д^ в пределах 5-30 % [11-17].
Для одночисловой оценки динамических характеристик насоса рационально использовать величину относительного отклонения. Расчёт такой величины возможен по оценке среднего медианного значения
ДА.
дМе
относительных отклонений в диапазоне частот каждой гармоники (/тМ.../тах): Д Ме = —_
А
. Величину
ДА
Ааи
дМе рационально использовать, если -=-
Ади
< 0,3 , т.е. когда доверительный интервал ДА^ не превышает
/.ах
30 % от средней амплитуды пульсирующего А .
Результаты оценки стабильности динамических характеристик насоса от характеристик стендовых систем4
составляют: по 1-й гармонике — 2 %, по 2-й — 8 %, по 3-й — 15 %, по 4-й — 48 %. Результаты расчёта
динамических характеристик насоса показали, что по четырём гармоникам в диапазоне частот 500-2500 об/мин
исследуемый шестерённый насос целесообразно рассматривать по модели эквивалентного источника
колебаний расхода. «и
к
Проверка применимости используемой модели насоса в виде эквивалентного источника колебаний расхода Е
Для проверки разработанного подхода оценки динамических характеристик насоса было выполнено <5 .. „и
расчётное определение пульсаций давления в стендовых системах с дросселем, полостью и протяжённым о
трубопроводом на выходе насоса. С использованием найденных динамических характеристик насоса К
рассчитаны амплитуды пульсаций давления в анализируемых стендовых системах, которые сравнивались с ^
экспериментальными. Ниже представлены результаты расчёта амплитуд Ар (/) в стендовой системе с ^
протяжённым трубопроводом. Амплитуды пульсаций давления р рассчитаны в соответствии с моделью ЭИКР к
К
по формуле: §
а • 2 • 2 (6)
Ч.и и i гл
Р< =
2, + 2.
обозначены амплитуды, полученные экспериментально, а сплошной линией — рассчитанные по формуле (6).
Результаты расчёта амплитуд пульсаций давления за насосом приведены на рис. 6. Пунктирной линией
й
3 Стабильность других собственных характеристик насоса оценивается аналогично.
4 Необходимо учитывать методическую погрешность подхода и инструментальную погрешность определения давления.
£ 1,0
0,8
0,6
0,4
к ч с
8 <
0,2
0,0
83
283
483
683
883
1083
1283 1483
1683
Частота гармонических составляющих спектра колебаний, Гц
и Й О
тз
л 'й
и
£ -Й
Рис. 6. Амплитудно-частотная характеристика давления для стендовой системы с протяжённым трубопроводом длиной /=106 м (сплошными линиями обозначены расчётные данные по: — первой; — второй;
— третьей; — четвёртой гармониках; пунктирными линиями обозначены экспериментальные данные)
Оценка сходимости результатов проверочного расчёта проведена по формуле отклонения разницы экспериментальных и расчётных значений от экспериментальных данных. Результаты проверочного расчёта для стендовых систем с дросселем, полостью и протяжённым трубопроводом на выходе насоса показали, что по 1-й гармонике относительное отклонение от экспериментальных данных для рассмотренных стендовых систем не превышает 10 %, минимальное относительное отклонение характерно стендовой системе с протяжённым трубопроводом на выходе насоса — 2%. Амплитуды пульсаций давления 2-й-4-й гармоник находятся в пределах
коридора, рассчитанного по формуле А = А - А (где А — средняя амплитуда пульсирующего давления).
Обсуждение и заключения. Реализована методика расчёта динамических характеристик насоса в виде эквивалентного источника колебаний расхода жидкости на примере шестерённого насоса производительностью 14 см3/об с использованием разработанных авторами специальных стендовых систем и уточнённой модели импеданса и пульсационной производительности по расходу. Показано, что исследуемый шестерённый насос внешнего зацепления рационально рассматривать по модели эквивалентного источника колебаний расхода. Проведена проверка предлагаемого подхода в виде расчёта амплитуд пульсаций давления в стендовых системах и сравнения с экспериментальными значениями. Отклонение расчётных амплитуд пульсаций давления от экспериментальных значений не превысило 10 %.
Стоит отметить, что за рамками обзора осталось влияние на динамические характеристики шестерённого насоса приводной механической системы и износа деталей шестерённого насоса. На основе проведённой работы можно сформулировать направления развития предлагаемого подхода:
- оценка динамических характеристик шестерённого насоса в течение всего срока службы;
- разработка четырёхполюсной модели шестерённого насоса как эквивалентного источника колебаний, учитывающей приводную механическую и выходную гидравлическую системы;
- разработка модели шестерённого насоса как эквивалентного источника колебаний, учитывающей его конструктивные особенности, процесс вытеснения, обратный гидроудар в межзубовую камеру, работу запертого объёма и перетечки из области нагнетания в область всасывания.
При этом для подтверждения работоспособности предлагаемого подхода на весь класс машин объёмного вытеснения необходимо расширение номенклатуры испытуемых объёмных насосов.
Список литературы
1. Рыбак, А. Т. Моделирование и экспериментальные исследования гидромеханической системы со знакопеременной нагрузкой / А. Т. Рыбак, В. П. Жаров, Р. А. Фридрих // Вестник Донского государственного технического университета. — 2006. — Т. 6, № 1. — С. 17-25.
2. Рыбак, А. Т. Объемная жесткость и ее влияние на динамику гидромеханической системы / А. Т. Рыбак // Вестник Донского государственного технического университета. — 2006. — Т. 6, № 3. — С. 200-207.
3. Озерский, А. И. Применение подхода Лагранжа к решению задач динамики гидравлических систем гидроприводных и теплоэнергетических установок / А. И. Озерский // Вестник Донского государственного технического университета. — 2010. — Т. 10, № 6. — С. 914-925.
4. Формирование динамических свойств трубопроводных цепей / Н. Д. Быстров, А. Г. Гимадиев, А. Н. Головин [и др.]; под ред. В. П. Шорина. — Самара : Изд-во Самар. ун-та, 2020. — 328 с.
5. The equivalent circuits Thévenin and Norton / M. Staniloiu, H. Popescu, G. Rezmeri^a, M. Iordache // The Scientific Bulletin of Electrical Engineering Faculty. — 2021. — Vol. 21. — P. 40-48.
6. Thévenin, L. Sur un nouveaux théorème d'électricité dynamique / L. Thévenin // Comp. Rendus hebdomadaires des Séances de l'Académie des Sci. — 1883. — Vol. 97. — P. 159-161.
7. Generalized Functions and Fourier Analysis / M. Oberguggenberger, J. Toft, J. Vindas, P. Wahlberg (ed.) — Springer International Publishing AG, 2017. — 276 p.
8. Handbook of Mathematics, 5th ed. / I.N. Bronshtein, K.A. Semendyayev, G. Musiol, H. Muhlig. — Berlin: Springer, 2015. — 1207 p.
9. Артюхов, А. В. Особенности частотных испытаний гидравлических насосов / А. В. Артюхов, Л. И. Брудков // Вибрационная прочность и надежность двигателей и систем летательных аппаратов: межвуз. сб. науч. тр. — Куйбышев : КуАИ, 1982. — Вып. 9. — С. 12-17.
10. Sundararajan, D. Introductory circuit theory / D. Sundararajan. — Switzerland : Springer, 2020. — 297 p. https://doi.org/10.1007/978-3-030-31985-4
11. Ichiyanagi, T. The temperature characteristics of flow ripple and source impedance in an external gear pump / T. Ichiyanagi, T. Nishiumi, Sh. Nakagawa // In: Proc. 10th JFPS Int. Symposium on Fluid Power, 2017. — 5 p.
12. Артюхов, А. В. Методика определения динамических характеристик гидравлических насосов / А. В. Артюхов, В. П. Шорин // Динамические процессы в силовых и энергетических установках летательных аппаратов аппаратов : межвуз. сб. науч. тр. — Куйбышев : КуАИ, 1988. — С. 70-77.
13. Jinghao Liu. Source flow ripple and source impedance measurement for different hydraulic pumps / Jinghao Liu // In: Proc. 47th Int. Congress and Exposition on Noise Control Engineering (InterNoise 2018). — Chicago, Illinois, 2018. — 7 p.
14. Liselott Ericson. On Fluid Power Pump and Motor Design - Tools for Noise Reduction. PhD Thesis / Liselott Ericson. — Linkoping University, Sweden, 2011. — 130 p.
15. Bramley, C. Comparison of methods for measuring pump flow ripple and impedance / C. Bramley, D. N. Johnston // In: Proc. ASME/BATH 2017 Symposium on Fluid Power and Motion Control. — 2017. — 11 p. http://dx.doi.org/10.1115/FPMC2017-4223
16. Jinghao Liu. Hydraulic Fluid-Borne Noise Measurement and Simulation for Off-Highway Equipment / Jinghao Liu, Sanghoon Suh, Yuzhen Yang // In: Proc. 2017 SAE/NOISE-CON Joint Conf. (NoiseCon 2017), Grand Rapids, MI. —2017. — Vol. 9. — P. 435-443.
17. Theoretical and experimental studies of a switched inertance hydraulic system including switching transition dynamics, non-linearity and leakage / Min Pan, Nigel Johnston, Andrew Plummer [et al.] // Proc. Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering. — 2014. — Vol. 228. — P. 802-815. https://doi.org/10.1177/0959651814548299
Поступила в редакцию 27.04.2022 ш
К
Поступила после рецензирования 20.05.2022 К
Принята к публикации 20.05.2022
ю
Об авторах: щ
К
Санчугов Валерий Иванович, профессор кафедры «Автоматические системы и энергетических установок» д
Самарского национального исследовательского университета имени академика С. П. Королева (443086, РФ, ^
г. Самара, ул. Московское шоссе, 34), доктор технических наук, ОЯСЮ, Sanchugovv.i@vandex.ru К
<о К
Рекадзе Павел Дмитриевич, инженер-конструктор НИИ-201, Самарского национального щ
исследовательского университета имени академика С. П. Королева (443086, РФ, г. Самара, ул. Московское о^
шоссе, 34), Я^еаг^егГО, ScopusID, ОЯСЮ, Яек^е 1993@gmail.com £
о
Заявленный вклад соавторов: §
В. И. Санчугов — научное руководство, формирование основной концепции, цели и задачи исследования; Э П. Д. Рекадзе — формирование основной концепции, проведение расчетов, подготовка текста, формирование выводов, анализ результатов исследований, доработка текста, корректировка выводов.
Конфликт интересов.
Авторы заявляют об отсутствии конфликта интересов.
Все авторы прочитали и одобрили окончательный вариант рукописи