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

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

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

 

12αμasπdH

 

2

 

κ +1

γ =

(

)

2 (κ −1)

E.

 

 

 

π H 3 ps

κ + 1

 

 

 

Параметры γ, γ è γ называют конструктивными параметрами питателей.

Åñëè x = xd – решение уравнения (3.36) или уравнений (3.38) (в зависимости от типа питателя), то давление на кромке питателя

определяется по формуле pd

= ps

xd . При этом распределение дав-

лений в смазочном слое находится на основе формулы (3.31):

 

 

p 2

p 2

 

R

 

p =

 

d

 

a

ln

o

+ ps 2 .

(3.39)

 

 

E

 

 

 

 

 

 

 

r

 

Несущая способность рассматриваемой дисковой опоры вычис-

ляется с помощью ЭВМ по формуле

 

 

 

 

R

 

 

 

 

 

 

F = 2π O

( p pa )rdr.

 

 

 

 

0

 

 

 

 

 

Коэффициент при lnr в формуле (3.39) определяет коэффициент A в формуле (3.31):

A = − pd 2 pa 2 .

E

Поэтому расход газа для дискового подпятника можно определить по формуле

Q =

κπ H 3

 

p 2

p

2

.

 

 

 

d

a

 

2

 

 

 

 

12μas

 

 

E

 

 

Расчет полосовой опоры с щелевым дросселем

Рассмотрим (рис. 3.3) прямоугольную полосовую опору протяженности T и ширины L, причем L << T. Пусть газ в смазочный слой подается через продольную капиллярную щель, имеющую ширину b и такую же протяженность, что и опора. Давление на входе в щелевой дроссель принимаем равным ps, а давление окружающей среды – pa. Принимаем толщины смазочного слоя H и дросселирующей щели c

61

постоянными. Для упрощения исследования пренебрегаем концевыми эффектами в опоре и на входе в щелевой дроссель.

 

 

Х

У

а

б

 

Ps

Ðèñ. 3.3. Схема щелевой опоры

Åñëè îñü Ox направить поперек питающей щели, то комплексные потенциалы потоков в зазоре определятся, очевидно, по формуле (3.28), из которой следует, что для смазочного слоя в опоре справедливы равенства

 

 

p 2

 

p 2

 

X + pd 2 , l1 < X < 0,

 

 

 

 

 

d

 

 

a

 

 

 

 

 

 

 

l1

 

 

p

2

 

 

 

 

 

 

 

 

 

 

(3.40)

 

=

 

 

pd

2

pa

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

X + pd

, l2 > X > 0,

 

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2)Y , l < X < 0,

 

 

 

 

 

 

( p 2

 

 

ψ =

 

d

 

 

a

l

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

( p 2p 2) Y , l > X > 0.

 

 

 

 

 

 

 

d

 

 

 

a l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

Масса газа, протекающего в единицу времени из участка питающей щели шириной dY в область смазочного слоя X < 0, согласно формулам (3.24) и (3.38), равна

dQ =

κ H 3

 

( p 2

p 2 )

dY

,

(3.41)

24μa

 

 

1

2

d

a

l

 

 

s

 

 

1

 

 

62

а в область X > 0

dQ

 

=

κ H 3

 

( p

2 p 2 )

dY

.

(3.42)

 

24μas

 

 

 

2

 

2

d

a

l2

 

Так как поток газа в дросселирующей щели полностью аналогичен потоку газа в смазочном слое, но происходит под действием давлений ps è pd, то нетрудно записать, что через участок щели протяженности dX в единицу времени протекает масса газа, равная

dM =

κ H 3

 

( p 2

p 2 )

dY

.

(3.43)

24μas

 

 

 

2

s

d

b

 

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

dQ1 + dQ2 = dM,

из которого, учитывая формулы (3.39)–(3.41), следует уравнение для определения давления на выходе из щели

(

1

+

1

)H 3 ( pd

2

pa

2 ) =

c3

( ps

2 pd

2 ).

(3.44)

 

 

 

b

 

l1 l2

 

 

 

 

 

 

 

Легко находим

 

 

 

 

 

 

pd = ps

γ + p

2

,

 

 

 

 

 

 

1

 

 

 

 

 

 

 

γ + 1

 

 

 

 

 

 

 

 

 

 

ãäå γ =

 

l l

c3

, p1

= pa / ps .

 

 

 

 

1

2

 

 

 

 

(l

+ l

)bH 3

 

 

 

1

 

2

 

 

 

 

 

 

Несущая способность части полосовой опоры, лежащей в области X < 0, вычисляется по формуле

0

0

2

2

 

 

F1 = T pdX = T

pd

pa

X + pd

2 dX .

 

 

l1

l1

 

l1

 

Интеграл легко вычисляeтся. Получаем

F1 =

2

Tl1

( pd

+

pa

2

).

3

 

 

 

 

 

 

pd + pa

63

Аналогично запишется формула для несущей способности части полосовой опоры, лежащей в области X > 0:

F2 =

2

Tl2

( pd

+

pa

2

).

3

 

 

 

 

 

 

pd + pa

Очевидно, что полная несущая способность полосовой опоры с щелевым поддувом равна

F =

2

TL( pd

+

pa

2

).

(3.45)

 

 

 

 

3

 

 

pd + pa

 

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

Q = M =

κ c3

( ps 2

pd

2 )

T.

(3.46)

24μas

2b

 

 

 

 

 

3.4. Метод источников и стоков

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

64

Комплексный потенциал дорожки источников в плоскости

Используем декартовы координаты ξ, η. Îñü Oξ направляем вдоль границы полосы (рис. 3.4). Пусть ζ = ξ + . Рассматриваем случай однорядной опоры ширины L, когда все питатели одинаковые, имеют диаметры d и расположены равномерно с шагом B на линии η = l, так что их центры имеют координаты

ζk = ζ0 + Bk, k = 0, ±1, ±2,.. ,

ãäå ζ0 – координата центра питателя, выбираемого в качестве базового ζ0 = ξ0 + il.

η

2 /n

λ 2

L

 

1

 

λ

O

ξ

Ðèñ. 3.4. Схема полосовой однорядной опоры

Построим комплексный потенциал, порожденный рассматриваемой дорожкой питателей в безграничном слое, имитируя каждый питатель точечным источником мощности Q′. Для источника, помещенного в точку ς = ςk, имеем (см. формулы (3.32) и (3.48))

w = −

12μa

2

Qln(ζ − ζ

 

) + c

 

=

 

 

 

 

s

 

 

 

 

 

κπ H 3

 

 

 

 

k

 

 

 

 

 

 

 

k

 

 

k

 

 

 

= −

12μa 2

Qln(ζ − ζ

 

Bk) + c

 

=

 

 

 

s

 

 

κπ H 3

0

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

12μa

2

Qln(1

ζ − ζ

k

) + c

, (k 0).

 

 

s

 

 

 

κπ H

 

 

Bk

 

 

3

 

 

 

 

 

 

 

k

 

 

 

 

 

65

Найдем комплексный потенциал всей системы источников как сумму комплексных потенциалов каждого из них:

 

 

12μas

 

 

 

 

 

 

 

 

 

ζ

− ζ0

 

ζ − ζ0

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w = −

 

 

Q {[ln(ζ − ζ

0

) +

∑ln[(1

 

 

 

)(1 +

 

)]} +c =

κπ H 3

 

 

Bk

Bk

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

ζ − ζ0

 

 

 

 

 

 

= −

12μas

Qln{(ζ − ζ

 

)(1 (

)2 ]} + c.

 

 

 

0

 

 

 

 

κπ H 3

 

k =1

 

 

 

 

Bk

 

 

 

 

 

 

Так как имеет место равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

2

 

 

 

 

 

 

 

 

 

 

 

β∏(1

 

 

 

 

) = sin β,

 

 

 

 

 

 

 

 

k 2π 2

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

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

w = −

12μa

2

Qln sin

π (ζ − ζ

)

+ D,

s

 

0

 

κπ H 3

B

 

 

 

 

 

ãäå D = D0 + iD* – комплексная постоянная.

Комплексный потенциал дорожки источников в полосе

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

имеют координаты

ζ m = ζ0 + 2mLi, (m = ± 1, ± 2, ...), а базовые сто-

êè – ζ m=

 

 

(m′ = ± 1, ± 2, ...;

 

 

ζ0 + 2mLi,

ζ0 = ξ0 iη0 ) , и некоторого

поступательного потока, нормального к границам полосы. Здесь при построении искомого комплексного потенциала нельзя непосредственно суммировать все бесконечное множество комплексных потенциалов, так как такая сумма будет расходящейся. Поэтому рассмотрим сначала случай, когда в плоскости расположено ограниченное число (2K + 1) дорожек источников и столько же дорожек стоков.

Суммарный потенциал, порождаемый дорожками источников, равен

66

 

12μas

π (ζ − ζ0 )

 

π (ζ − ζ k )

 

π (ζ − ζ k )

 

 

2

 

 

K

 

 

 

 

u = −

 

Q {[ln sin

 

+ ln[sin

 

sin

 

]}.

κπ H 3

 

 

 

1

 

B

k =1

B

 

B

 

Используя формулу для произведения синусов, можно записать с точностью до постоянной

 

 

12μas

2

 

π (ζ − ζ0 )

K

cos

2π (ζ − ζ0 )

 

 

 

 

B

 

u

= −

 

 

Q {[ln sin

 

+ ln[1

 

 

 

 

]}.

κπ H 3

B

 

 

4kL

1

 

 

k =1

 

ch

 

 

 

 

 

 

 

 

 

 

 

 

B

Аналогично, суммарный потенциал, поставляемый дорожками стоков, равен

 

 

 

12μas

2

 

π (ζ −

 

 

 

cos

2π (ζ − ζ 0 )

 

 

 

 

ζ0 )

K

 

B

 

u

 

=

 

 

 

Q {[ln sin

 

 

 

+ ln[1

 

 

 

 

]}.

2

 

κπ H 3

B

 

ch

4kL

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

Введем еще в рассмотрение однородный поступательный поток, направленный поперек полосы, комплексный потенциал которого для удобства представляем в виде

u

 

=

12μa 2

Qb i

ζ

+ e ,

(3.47)

 

 

s

 

 

 

κπ H 3

L

 

3

 

 

k

k

 

ãäå bk è ek – постоянные, подлежащие определению. Рассмотрим теперь поток с суммарным потенциалом

 

 

 

 

 

 

 

 

 

 

uK = u1 + u2 + u3 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π (ζ −

 

 

 

 

 

4kL

 

2π (ζ −

 

 

 

 

 

 

12μas

2

 

 

 

ς

 

sin

ζ0 )

 

K

ch

cos

ζ0 )

 

 

 

 

 

 

 

B

 

 

 

B

 

B

 

 

 

 

 

=

 

 

 

Q {b i

 

+ ln

 

 

 

 

 

+ ln[

 

 

 

 

 

 

 

 

]} + e

 

.

 

κπ H 3

 

 

π (ζ − ζ0 )

 

 

4kL

 

2π (ζ − ζ

0 )

 

 

 

 

K

 

L

 

sin

 

k =1

ch

cos

 

K

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

B

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определим постоянные bK è eK, выполняя следующие условия:

 

 

 

 

 

 

Re[uK (ξ0 )] = pa

2 ,

 

Re[uK (ξ0

+ iL)] = pa

2 .

 

 

 

 

 

 

67

Используя формулы sin (ia) = isha è cos (ia) = cha, находим

 

sh

π

[(2k + 1)L + l]sh

π

[(2k 1)L l]

 

sh

π

[L + l]

K

B

B

 

B

uk = ln

 

 

 

 

 

 

 

+ ln

 

 

=

 

π

 

 

 

 

π

 

 

 

π

 

k =1

sh

[(2k + 1)L l]sh

[(2k 1)L

+ l]

 

sh

[L l]

 

 

 

 

B

B

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

π

[(2K + 1)L + l]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ln

 

B

 

 

.

 

 

 

 

 

 

 

 

sh

π

[(2K + 1)L l]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

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

 

 

 

 

 

lim bK

=

2πl

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w = lim u

 

=

12μa

2

Qf (ζ ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κπ H 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π (ζ −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

ζ

0 )

 

 

2πl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

2

 

ãäå

f (ζ ) = v(ζ ,ζ0 ,ζ0 ) + ln

 

 

 

 

 

 

 

 

 

+

iζ + pa

,

sin

 

 

π (ζ − ζ

0 )

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 cos

2π (ζ −

 

 

 

4kL

 

 

 

 

 

 

ζ0 )

/ ch

 

 

 

 

 

 

 

 

 

 

 

 

 

v(ζ ,ζ0 ,

ζ0 ) = ln[

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

B

]}.

 

 

 

 

 

 

 

 

 

 

2π (ζ − ζ0 )

 

 

4kL

 

 

 

 

 

k =1

1 cos

/ ch

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

B

 

 

 

Отметим два свойства аналитической функции важных для дальнейших исследований.

(3.48)

v(ζ ,ζ0 ,ζ0 ) ,

68

1. На границе полосы η = 0 реальная часть функции v равна нулю, что можно записать в виде

Re[v(ζ ,ζ0 ,

 

 

 

ζ0 )] = 0.

(3.49)

В справедливости этого свойства легко убедиться простой проверкой, учитывая равенство

1

Re[v(ξ,ζ0 ,ζ0 )] = [v(ξ,ζ0 ,ζ0 ) + v(ξ,ζ0 ,ζ0 )]. 2

Из вида функции v (3.48) следует, что

v(ξ,ζ0 ,ζ0 ) = −v(ξ,ζ0 ,ζ0 ),

что и доказывает справедливость утверждения (3.49). 2. Имеет место равенство

B / 2o

 

Re v(ζ ,ζ0 ,ζ0 ) η= const dζ = 0.

(3.50)

B / 2o

Для доказательства этого свойства рассмотрим интеграл

Коши от функции v(ξ, ς 0 , ς0 ) по прямоугольному контуру ABCD, ограниченному прямыми η = 0, (AB); ξ = B/2 + ξ0, (BC); η = const, (CD); ξ = – B / 2 + ξ0, (DA). Реальная часть этого интеграла равна нулю в силу аналитичности неоднозначной логарифмической подинтегральной функции

Re v(ζ ,ζ0 ,ζ0 )dζ = 0.

С

Можем записать

Re[ v(ς , ς 0 , ς0 )dς + v(ς , ς 0 , ς0 )dς + v(ς , ς 0 , ς0 )dς + v(ς , ς 0 , ς0 )dς ] = 0.

AB BC CD DA

В силу периодичности функции v ïî ξ имеем

v(ς , ς 0 , ς0 )dς + v(ς , ς 0 , ς0 )dς = 0.

BC DA

69

Âсилу первого свойства подынтегральная функция в интеграле

v(ς , ς 0 , ς0 )dς равна нулю, поэтому и сам этот интеграл равен нулю.

AB

Следовательно, имеет место равенство

Re v(ς , ς 0 , ς0 )dς = 0.

CD

η = const на отрезке CD, а концы этого отрезка совпадают с пределами интегрирования в формуле (3.50), справедливость которой таким образом доказана.

3.5.Расчет полосовых однорядных газостатических опор

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

Рассмотрим сперва однорядную полосовую опору протяженности T >> L, вдув газа в смазочный слой которой производится через n одинаковых питателей диафрагменного типа, расположенных равномерно по всей длине опоры в один ряд, отстоящий от боковой границы полосы на расстоянии L (ñì. ðèñ. 3.4).

При расчете такой опоры будем пренебрегать концевыми эффектами и рассматривать ее как часть бесконечно длинной полосовой опоры с одним рядом питателей, имеющем шаг B = T/n. В этом случае при H = const для комплексного потенциала потока в смазочном слое применима формула (3.48), а среднеквадратичое давление на кромке питателя может быть вычислено по формуле

p d 2 =

1

Re

 

w(ζ )

 

dζ

 

.

(3.51)

 

 

 

 

πd

 

ζ −ζO

 

= d / 2

 

 

 

 

 

 

 

 

 

 

 

Здесь контурные интегралы от аналитических в круге фун-

кций, входящих в w(ζ), вычисляем по теореме о среднем, а от функции, имеющей логарифмическую особенность в центре круга, – непосредственно. При вычислениях учитываем, что на кромке питателя

ζ − ζ0

=

1

deiθ ,

 

dζ

 

=

1

 

 

dθ

 

,

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

,

ln sin π (ζ ζ0 ) ln πd + iθ.

B 2B

70

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