Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теоретические основы работы газостатических опор. Монография

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
912 Кб
Скачать

5. РАСЧЕТ СЕГМЕНТНЫХ ГАЗОСТАТИЧЕСКИХ ОПОР

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

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

Наддув газа в смазочный слой сегментных подшипников обыч- но осуществляется через один или несколько питателей типа кольцевой или простой диафрагмы. Если смазываемые поверхности в сегментном газостатическом подшипнике остаются параллельными во время работы машины, а зазор между ними остается постоянным, то в случаях малых и умеренных чисел сжимаемости Λ вращение вала незначительно влияет на распределение давлений в смазочном слое опоры и можно принять, что квадрат давлений в слое является гармонической функцией. В этих случаях для расчета поля давлений в смазочном слое каждого сегмента можно использовать методы источников и конформных отображений, изложенные в разделе 3. Основная трудность при этом связана с построением соответствующего комплексного потенциала для прямоугольной области, а также в нахождении комплексной функции, реализующей конформное отображение сегмента на прямоугольник. Для

121

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

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

5.1. Комплексный потенциал для прямоугольной газостатической опоры с дискретным поддувом

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

Если зазор в такой опоре постоянный, а эффект относительного перемещения смазываемых поверхностей незначителен и его можно не учитывать, то квадрат давлений р2 в смазочном слое удовлетворяет уравнению Лапласа (3.20), что позволяет ввести в

рассмотрение комплексный потенциал потока U = p2 + iΨ. Воспользуемся декартовой системой координат ξ è η, îñè Oξ è

Oη которой проходят вдоль границ прямоугольника. Пусть центры питателей, диаметры выходных кромок которых dj расположены в точках с комплексными координатами ςj = ξj + iηj (j = 1, 2, …, N), L1 è L2 – размеры опоры вдоль осей Oξ è Oη соответственно, N – число питателей в сегменте.

Комплексный потенциал потока U(ς) = p2 + iΨ в смазочном слое прямоугольной газостатической опоры должен удовлетворять условиям

 

Re (U (ς ))= pa2 ïðè ξ = 0, ξ = 0, η = 0, η = L2,

(5.1)

 

 

1

 

 

 

 

 

2

 

2

 

 

Re

 

 

 

 

 

U (ς ) | dς |

=

pdk

>

pa

, k = 1, 2, …, N.

(5.2)

 

 

 

 

π

 

dk

 

ς −ς k

 

= dk / 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем в рассмотрение комплексный потенциал Uj(ς) потока, порожденного только j-м питателем, помещенным в точку A (ςA = ςj) опоры.

122

Поток газа в смазочном слое опоры от питателя, помещенного в точку A, удовлетворяющий условиям (5.1)–(5.2), можно рассматривать в плоскости ς как результат наложения потоков, порожденных четырьмя семействами дорожек источников и стоков, параллельных оси , имеющих шаг 2L1 и период по η, равный 2L2. При этом в базовые точки этих семейств A è B (ςB = –ςj) попадает источник мощности

Qj, а в базовые точки C (ςC = ςj) è D (ςD = –ς j) попадают стоки (источники мощности (–Qj). Обильности источника Qj надо подобрать

таким образом, чтобы выполнялось условие (3.54) при k = j. Базовый источник в m-й дорожке источников семейства A

имеет координату

ςm = ςA + 2miL2, m = 0, ±1, ±2, ...

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

Переменная часть комплексного потенциала потока, наводимого m-й дорожкой источников семейства A, в соответствии с формулой (3.49), равна

U

(m)

=

12μas2Qj

ln sin

π

(ς ς

 

).

j

κπ H 3

2L

m

 

 

 

 

 

 

 

 

 

 

1

 

 

 

Совокупность всех дорожек источников и стоков семейств A è C порождает комплексный потенциал, аналогичный потенциалу (3.50), но без слагаемого, отвечающего поступательному потоку. Вся совокупность дорожек источников и стоков семейств A, B, C è D порождает комплексный потенциал

 

12μa2Q

j

 

 

U j =

s

 

uj (ς ),

(5.3)

κπ H

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

ãäå

uj (ς ) =

V (ς , ς j

, ς j

) ln

sin

 

 

(ς − ς j

)

,

 

2L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

π

(ς − ς j

) sin

π

 

(ς + ς j )

 

 

 

 

 

 

V (ς , ς j , ς j ) = ln

 

 

2L1

 

 

 

 

2L1

 

 

+

 

 

 

 

sin

π

 

(ς + ς j

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2L1

 

 

 

 

 

 

 

123

 

 

 

 

 

+

 

ln

 

k =

1

 

 

 

 

 

 

 

π

 

 

 

 

 

2π

kL

 

 

π

 

 

 

 

2π

kL

 

1cos

(ς −

 

 

)

ch

2

1cos

(ς +

 

 

)

ch

2

 

ς

j

ς

j

 

L1

 

 

 

 

 

 

L1

+ ln

 

L1

 

 

 

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

2π

kL

 

π

 

 

 

 

2π

kL

.

1cos

(ς − ς j )

ch

2

1cos

(ς + ς j )

ch

2

 

L

L

L

L

 

 

1

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

j

= ξ j iη j.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ς

 

 

 

 

 

 

 

 

Комплексная функция uj(ς) обладает теми свойствами, что ее реальная часть fj (ξ, η) = Re (uj (ς)) обращается в нуль на всех сторонах рассматриваемого прямоугольника и что

f j

d = −2π δ kj ,

∂ν

Γ

k

 

ãäå Γj – контур питателя.

Первое свойство вытекает из способа построения функции uj(ς). Второе свойство следует из формулы (5.2), в которой надо учесть, что

 

2

 

p2j

=

12μas

Qj f j (ξ, η) + const.

3

 

 

κπ H

Комплексный потенциал потока, порожденного всеми N питателями опоры, находится суммированием всех потенциалов Uj(ς) j =

1, 2, ..., N è pa2 . В результате имеем

2

N

 

U (ς ) =

12μas

Qjuj (ς ) + pa2 .

(5.4)

3

 

κπ H

j=1

 

Мощности источников Qj, входящие в эту формулу, могут быть определены по значениям давлений pdj на выходных кромках питателей из системы уравнений (5.2), после подстановки в них выражения (5.4).

При выполнении условий (5.2) надо вычислять контурный

интеграл u(ς ) | dς | . Здесь удобно поступить следующим образом.

Γ k

124

Отдельно вычислить интеграл по контуру k-го питателя от той части функции U(ς), которая имеет логарифмическую особенность в точке ς = ςk, а значения интегралов от остальных слагаемых, входящих в функцию U(ς) и остающихся аналитическими в точке ς = ςk, могут быть вычислены по теореме о среднем, утверждающей, что среднее значение аналитической функции на окружности равно значению этой функции в центре окружности.

Ïðè

π | ς − ς k

|

<< 1 , что обычно имеет место, как показано в

2L1

 

 

 

 

разделе 3, справедливо соотношение

 

1

 

ln

 

π

 

 

 

π d

k

 

Re

 

 

sin

 

(ς − ς k

) | dς |

ln

 

.

 

 

2L1

4L1

 

π

dk |ς −ς k |= dk / 2

 

 

 

 

 

 

С учетом этой формулы можем записать, что

1

 

 

 

 

 

 

 

 

 

 

4L1

 

 

 

 

π

 

 

 

 

 

 

 

 

f j

d

= Re V (ς k , ς j , ς j ) + δ kj

ln

 

ln

δ kj

+ sin

 

(ς k

− ς j

)

,

(5.5)

 

 

 

π dk

2L1

 

π

dk Ck

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а условия (5.2) принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12μ

2

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aS

 

 

ΦkjQj = pkj2

pa2 ,

k = 1,

2, ,

N ,

 

 

 

(5.6)

 

 

 

 

 

κ

 

π

H

3

 

 

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

ãäå Φkj

=

 

1

 

 

 

f j

 

 

d .

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k C j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если давления pdk заданы, то из уравнений (5.5) могут быть определены мощности источников Qj, порождающих рассматриваемый поток в газостатической прямоугольной опоре, и наоборот, если мощности источников Qj заданы, то эти уравнения служат для определения давлений на выходных кромках питателей.

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

125

тому, как это делалось для упорных и радиальных полноохватных подшипников в предыдущих разделах.

Например, если не учитывается вторичное дросселирование в питателях, то условие баланса расходов записывается в виде

Qj = Mj, j = 1, 2, …, N,

(5.7)

ãäå Mj вычисляется по формуле ().

Полагая x j = pdj2 pS2 и исключая Qj из уравнений (5.6) и (5.5),

получаем систему уравнений для определения давлений на кромках питателей

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

xk2 p12

= Φkjγ j q ( x j ), k = 1, 2, , N ,

(5.8)

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

12αμ aS Aj

 

2

κ +1

 

 

ãäå γ

 

=

2 (κ −1)

конструктивные параметры пита-

 

 

 

 

 

 

 

 

j

 

π pS H 3

 

 

κ + 1

 

 

 

 

 

 

 

телей.

 

 

 

 

 

 

 

 

 

 

Определив из системы нелинейных уравнений (3.7) величины xk, по формуле (3.13) находим значения величин Mj, а следовательно, и Qj, что и завершает определение функции U(ς).

5.2. Комплексный потенциал для кольцевого сегмента упорного газостатического подшипника с дискретным поддувом

Рассмотрим один сегмент кольцевого упорного газостатического подшипника. Для удобства направим оси цилиндрической системы координат r, φ, z таким образом, чтобы φ = 0 отвечало передней кромке сегмента а φ = θ – задней кромке сегмента.

Пусть R = Ri и R = RÂ соответственно уравнения внутренней и внешней границ сегмента.

Ïðè H = const è Λ = 0 квадрат давления в смазочном слое является гармонической функцией и можно ввести в рассмотрение комплексный потенциал

126

W(z) = p2(r, φ) + iΨ(r, φ),

являющийся аналитической функцией комплексной координаты z = re.

Принимаем, что питатели имеют диаметры выходных кромок dj << Rj и расположены так, что их центрам отвечают координаты

z j = rj eiϕ j , j = 1, 2, …, N.

(5.9)

Контуры питателей при этом описываются уравнением

z = z j

+

1

d j eiθ .

 

 

2

 

С помощью конформного преобразования

ς = ln

z

(5.10)

RВ

 

 

сегмент отображается на прямоугольник (рис. 5.1) в плоскости фиктивного потока ς = η + iφ, границы которого ограничены линиями

η = 0 è η = λ = ln RН ,

RВ

а также линиями φ = 0 è φ = θ.

Y

 

 

 

R

H

 

 

 

R

B

 

φv

 

 

 

 

rj

d

 

 

 

φj

 

 

X

 

 

 

η

 

 

d/rj

j

 

η

 

2

 

L

 

0

ξ

ξj

 

L1

Ðèñ. 5.1. Схема конформного отображения кольцевого сегмента

на прямоугольную область

127

Центры питателей в фиктивном потоке оказываются расположенными в точках

ς j

= ln

rj

+ iϕ j ,

 

 

 

RВ

а контуры выходных кромок питателей описываются в фиктивном потоке уравнениями

 

 

1 d j

 

i (θ −ϕ j )

d 2j

 

 

ς = ς j

+

 

 

 

e

 

+ O

 

 

,

2 rj

 

2

 

 

 

 

rj

 

 

свидетельствующими о том, что с точностью до величин поряд-

êà o (d 2j rj2 ) кромкам питателей в плоскости фиктивного потока отвечают окружности диаметров dj/Rj, а сам фиктивный поток представляет собой поток в прямоугольной области, порожденный источниками мощности Qj, расположенными в точках ςj è èìåþ-

щими диаметры выходных кромок d j = d j rj . Потенциал такого потока представлен формулой (5.4)

2

N

 

U (ς ) =

12μas

Qjuj (ς ) + pa2 ,

(5.11)

3

 

κπ H

j=1

 

ãäå uj(ς) даются формулой (5.3), в которой надо положить L1 = λ,

L2 = θ.

Мощности источников Qj связаны с давлениями на кромках питателей pdk формулами (5.6), где

 

1

 

 

 

 

 

 

 

 

 

4λr

 

Φkj =

 

f j | dς |= Re

V (ς k , ς j , ς j )

+ δ kj ln

 

 

 

k

 

 

 

 

 

 

π dk |ς −ς k |= dk /( 2 Γ k )

 

 

 

 

 

 

 

π

 

 

 

 

(5.12)

 

 

 

 

 

 

 

 

 

 

 

πdk

 

 

 

 

ln

 

δ kj

+ sin

 

(ς k

− ς j

)

 

.

 

 

 

 

2λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При вычислении величин V(ς , ς

, ς

) по формуле (5.3) необхо-

 

 

 

k

 

 

j

 

j

 

 

 

 

 

 

 

 

 

 

димо принять L1 = λ è L2 = θ.

При выполнении условий баланса расходов газа на питателях в формулах (3.15)–(3.16), определяющих расходы газа через

128

питатели Mj è M j , должны использоваться диаметры питателей в реальном потоке.

Например, если не учитывается вторичное дросселирование, то условия баланса расходов записываются в виде

Qj = Mj, j = 1, 2, …, N,

(5.13)

где при вычислении Mj надо использовать диаметры реальных питателей, а при вычислении Qj использовать диаметры питателей в фиктивном потоке.

Положив x j = pdj2 ps2 и исключив Qj из уравнений (5.6) и (5.12), получаем систему уравнений для определения давлений на кромках питателей:

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

xk2

p12

= Φkjγ j q ( x j ), , j = 1, 2, …, N,

(5.14)

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

12αμ as

Aj

2

κ +1

 

 

ãäå γ

 

=

2 (κ −1)

конструктивные параметры пита-

 

 

 

 

 

 

 

 

 

j

 

π ps H

3 κ + 1

 

 

 

 

 

 

 

телей.

Здесь, в отличие от прямоугольной опоры, коэффициенты Фkj вычисляются по параметрам (5.11), отвечающим фиктивному потоку, а γj – по параметрам, отвечающим реальному потоку.

Для расчета же прямоугольной опоры понятие фиктивного потока нами не использовалось, и Фkj вычислялись по параметрам реального потока.

Определив из системы нелинейных уравнений (5.13) величи- ны xk, по формуле (3.13) найдем величины Mj, а следовательно, и Qj, что и завершит определение комплексного потенциала U(ς) для сегментной газостатической опоры. При необходимости знать комплексный потенциал W(z) для сегмента в координатах реального потока достаточно положить в формуле (5.10)

ς = ln z ,

RB A

129

тогда

W (z) = 12μa

2

Qjuj (ln(z / RB )).

3s

 

 

 

N

 

κπ H

 

j=1

5.3. Общая постановка задач о расчете полей давлений в смазочных слоях газовых опор

с дискретным поддувом

Рассмотрим случай произвольной в плане газостатической опоры с дискретным поддувом, обе смазываемые поверхности которой подвижны.

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

H = H*h, s = p2/pS2, t' = σ*t,

Vaτ = V v'aτ , V= V v', u = v'aτ + v'bτ ,

x j = pdj2

pS2 ,

 

 

 

2

 

pS2 ,

 

 

 

 

(5.15)

x j =

pdj

p1 = pa

pS

,

Ω =

2Ω ',

=

', L =

1

L ', =

1

',

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L =

1

L ', σ

=

 

p

H 2

 

,

Λ =

6μ

V

 

.

 

 

S

 

 

 

 

 

2

 

 

 

2

 

 

2

 

 

 

 

 

12μ

 

 

 

 

pS H

 

 

Âэтих формулах H* – характерная толщина смазочной пленки,

– характерная протяженность смазочной пленки, V* – характерная продольная скорость смазываемых поверхностей, σ* – характерная частота сдавливания слоя, L – оператор Рейнольдса, pdj – давление на выходной кромке питателя, dj – давление в кармане, если питатель типа «простая диафрагма».

Уравнение Рейнольдса (1.21) преобразуется к виду

L '(s, h

3

) = 2

(h s )

+ 2Λ

 

(5.16)

 

 

 

t

' (h s u).

 

 

 

 

 

 

 

130

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]