Рассмотрим кольцевую газостатическую опору с дискретным поддувом. Внешний радиус опоры равен RH, а внутренний - RB. Полагаем, что зазор между смазываемыми поверхностями постоянен и равен h. Газ в смазочный слой опоры попадает под давлением ps через N произвольно расположенных на одной из смазываемых поверхностей питателей типа «простая диафрагма». При этом диаметры питателей могут быть разные: d0j, j = 1, 2, ..., N - диаметры подводящих каналов, а dj, j = 1, 2, ..., N - диаметры «карманов». Поток газа в смазочном слое принимается неодномерным, установившимся, ламинарным, изотермическим, а в питателях - одномерным, установившимся, адиабатическим, подчиняющимся законам динамики идеального газа.
Любому потоку вязкого газа между неподвижными параллельными плоскостями соответствует определенная аналитическая функция w(z) комплексного переменного z = x + iy - комплексный потенциал:
(1)
Комплексный потенциал на границах опоры удовлетворяет условию P = pa, а на границах питателей
Γj - условиям:
(2)
Из представления комплексного потенциала (1) получаем формулу для определения поля давлений в смазочном слое:
(3)
Расчет кольцевых газостатических опор порой очень затруднителен. В таких случаях для нахождения комплексного потенциала удобно воспользоваться методом конформных отображений и перейти к расчету полосовой опоры. Для этого свяжем с плоскостью опоры систему координат r, φ и введем новую комплексную координату . Конформное отображение дает функция . Тогда кольцевой опоре в фиктивном потоке z соответствует безграничная полоса шириной , на которой расположены N дорожек точечных источников. Диаметры питателей в полосовой опоре равны для подводящих каналов и для «карманов». Базовые координаты питателей будут . В - шаг, с которым расположены питатели в каждой дорожке. Полагаем параметры газа и толщину смазочного слоя в фиктивном потоке такими же, как и для реального газа. Следовательно, давления в фиктивном потоке в соответствующих точках будут такими, как и в реальном.
Комплексный потенциал фиктивного потока в смазочном слое полосовой опоры легко строится методом источников и стоков. В результате получаем:
(4)
Функции f(z), определяющие комплексные потенциалы потоков в полосе z, которые породили дорожки точечных источников, полагаем известными. μ - динамический коэффициент вязкости, κ - показатель адиабаты Пуассона, pa - давление окружающей среды, as - скорость звука в газе.
Условие (2) и равенство расходов газа через смазочный слой и через питатели Qj = Mj, дают систему нелинейных уравнений для определения давлений pdj на кромках питателей и расхода газа Qj. Но нужно помнить, что мы рассматриваем питатели типа «простая диафрагма», поэтому целесообразно применять схему «двойного дросселирования». Это означает, что на входе в «карман» полагаем площадь минимального сечения и эмпирический поправочный коэффициент , а на выходе из «кармана» - . Тогда для определения давлений на кромках питателей составляем две системы нелинейных уравнений:
(5.1)
(5.2)
Здесь
,
q(x) - газодинамическая функция, р1 - отношение давлений на входе в питатель и на выходе.
В этом случае количество газа, поступающего в карман в единицу времени через подводящий канал с диаметром d0j < dj, вычисляется по формуле:
(6.1)
А количество газа, вытекающего в единицу времени из «кармана» через кольцевую диафрагму под действием давления находим из равенства:
(6.2)
Зная расходы газа, строим комплексный потенциал и определяем поле давлений в смазочном слое. Далее можно исследовать влияние неравномерного поддува на интегральные характеристики опоры. При этом расчет может быть произведен для фиктивного потока без возвращения к реальной плоскости, но с использованием якобиана перехода
.
Так, например, несущая способность газостатической кольцевой опоры вычисляется по формуле:
(7)
Можно в формуле для определения поля давлений вернуться к полярным координатам r, φ (для простоты вычислений) и найти центр давлений как центр параллельных сил:
. (8)
ПРИМЕР. Рассмотрим кольцевую газостатическую опору с внешним радиусом RH = 0,5 м и внутренним - RB = 0,005 м. Возьмем N = 5. Диаметры подводящих каналов равны d0 = [0,002; 0,008; 0,004; 0,01; 0,006] м, а диаметры «карманов» пусть будут в два раза больше, т.е. dj = 2⋅d0j, j = 1, 2, ..., N. Координаты питателей и толщину смазочного слоя задают массивы r = [0,1; 0,46; 0,18; 0,34; 0,22] м, φ = [π/12, π/2, 5π/6, 5π/4, 11π/6], h = [0,01; 0,02; 0,03,; 0,04; 0,05] м. Давление на входе в питатели равняется двум атмосферным давлениям, т.е. ps = 2pa. Определим несущую способность и центр давлений этой опоры при данных величинах смазочного слоя.
С помощью метода конформных отображений переходим от кольцевой опоры к безграничной полосе ширины L, на которой расположены пять дорожек питателей. Определяем по схеме «двойного дросселирования» поля давлений на кромках питателей и расход газа. Строим комплексный потенциал для течения газа в смазочном слое, находим поле давлений.
Зависимость несущей способности (ось ординат, Н) от толщины газового слоя (ось абсцисс, м) показана на графике:
На расположение центра давлений толщина смазочного слоя оказывает очень слабое влияние. На исследуемой опоре центр давлений располагается практически в начале координат (xц ≈ -0,003; yц ≈ 10-6).
Список литературы
- Снопов А.И. Методические указания к курсу «Динамика вязкой жидкости и газа» часть II - Ростов н/Д.: УПЛ РГУ, 1990. - С. 4-7, 16-26, 28-31.
- Снопов А.И., Ларикова Н.А., Миронова Е.В. Моделирование газостатических опор с неравномерным дискретным поддувом // Современные проблемы науки и образования: Материалы конференции. - 2009. - №6. - С. 34-36.