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

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

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

Ищем, как и ранее, функцию П в виде

 

Π = Π0 (z) + εΠ1 (ϕ, z) + O(ε 2 ).

(4.71)

Функция П0 совпадает с найденной ранее (4.31)–(4.32), так как отвечает случаю ε = 0, при котором вращение вала не влияет на распределение давлений в смазочном слое модельного подшипника.

Подставляем выражение (4.71) в уравнение (4.69) и граничные условия (4.70), учитываем при этом формулу (4.43) и равенство

h3 = 1 3ε cos ϕ + O(ε 2 ).

Для определения функции П1 получаем уравнение

 

 

 

2Π

1

+

2

Π

1

= 2

Λ(sin ϕ Π

 

+

 

1

 

∂Π

1

)

(4.72)

 

 

 

∂ϕ 2

z2

 

0

 

Π0

 

∂ϕ

 

 

 

 

 

 

 

 

 

2

 

 

 

и краевые условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Π1 (ϕ, 0) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Π1 (ϕ, 2λ1 ) =

Re(Te

iϕ

),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.73)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Π

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|z= λ

= 0 (для двухрядного подшипника).

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим сначала область 0 < z < 2λ1. В этой области функцию

Ï1 èùåì â âèäå

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Π1 = Re[ y(z)eiϕ ].

 

 

 

 

 

(4.74)

Подставляем это выражение для П1 в уравнение (4.72) и первые два условия (4.73) и учитываем, что sin φ = Re(ie). Для определения функции получаем краевую задачу

y′′ − (1

Λi

) y

= 2Λi

Π0

 

 

 

,

 

Π0

 

 

 

 

 

 

(4.75)

y(0)

= 0, y(2λ1 ) = T.

 

 

 

 

 

 

111

Эта краевая задача не может быт решена точно аналитически, так как П0 есть функция от z. Для численного или приближенного аналитического решения задачи удобно функцию y представить в виде

y = T (

z

+ w1 ) + w2 ,

(4.76)

 

 

2λ1

 

где функции w1 è w2, в силу произвольности параметра Т, являются решениями таких краевых задач

w′′ − w (1

Λ

) = g

 

,

w (0) = w (2λ

) = 0, k = 1, 2,

(4.77)

 

k

k

k

Π0

 

 

 

k

k 1

 

 

 

 

 

 

 

 

 

 

 

 

ãäå

g1

= (1 +

Λi

 

 

)

z

,

g2 = 2Λi

Π0 .

 

Π0

 

 

 

 

 

 

 

 

 

2λ1

 

 

 

Ограничимся приближенным решением уравнений (4.77) по методу Бубнова–Галеркина с использованием только одной базисной функции. Полагаем

wk

= akζ (1 − ζ ), ζ =

z

, k = 1, 2.

(4.78)

 

 

2λ1

 

Коэффициенты ak находим следующим образом. Подставляем выражения (4.78) в уравнения (4.77), затем умножаем обе части этих уравнений на ζ(1 – ζ) и проводим интегрирование по ζ в пределах от 0 до 1. Получаем

1

(2λ1 )2 g1ζ (1 − ζ )dζ

a1 = − 1

0

.

 

[2 + (2λ1 )2 (1

Λi

)ζ (1 − ζ )]ζ (1 − ζ )dζ

 

 

0

Π0

Проведя вычисления, можем записать

2

 

1

 

 

 

 

a1 = −(2λ1 )

(

 

− ΛiM 21 )b,

 

12

(4.79)

 

 

 

 

 

 

 

 

2

 

 

 

a2 = −2(2λ1 )

ΛiN11b,

 

 

 

 

112

ãäå

b = [

1

+

 

1

 

(2λ1 )2

(2λ1 )2 ΛiM 22 ]1 ,

 

 

 

 

 

 

 

3

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

ζ k (1 − ζ )m dζ

 

 

 

 

 

 

 

 

M km

=

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

( p 2

 

p 2 )ζ + p 2

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 0

 

1

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N11 = ζ (1 − ζ ) ( pm2 0 p12 )ζ + p12 dζ.

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обращаемся теперь к межрядной зоне (в случае двухрядного

подшипника) и учитываем, что в этой зоне

 

 

Π0

= pm 0 . Äëÿ îïðå-

деления функции у получаем краевую задачу

 

 

 

 

 

 

y′′ − (1

 

Λi

) y

= 2Λipm 0

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

pm 0

 

 

 

 

 

(4.80)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(2λ1 ) = T , y(λ ) = 0.

 

 

 

 

 

 

 

Эта краевая задача имеет точное решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch[(λ − z) 1

Λi

 

]

 

 

 

 

 

2Λip

m 0

 

 

 

p 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = −

 

 

 

 

+ C

 

 

 

 

 

 

 

 

 

 

m 0

 

 

.

(4.81)

 

1

 

Λi

 

 

 

 

 

 

 

 

 

 

 

Λi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pm 0

 

 

 

 

ch[(λ − 2λ1 ) 1 pm2 0

]

 

 

 

 

 

 

 

 

 

 

 

 

 

Из первого краевого условия находим постоянную интегри-

рования С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C = T +

2λ1 pm 0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.82)

 

 

 

 

 

 

 

 

 

 

 

1

 

Λi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pm 0

Второе краевое условие выполняется тождественно.

Итак, в соответствии с формулами (4.76) и (4.81) функцию y окончательно можно представить в виде

y = TA (z) + B (z),

(4.83)

113

 

 

 

+ w1 (ζ ), ζ =

 

z

 

 

 

 

0 z 2λ1 ,

 

ζ

 

 

 

,

 

 

 

 

 

2λ1

ãäå

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(z) =

ch(λ − z)

β

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

2λ1 z ≤ λ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch(λ − 2λ1 )β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

w2

(ζ ), ζ

=

 

 

 

 

,

 

 

 

 

0 z 2λ1

 

 

2λ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B(z) =

 

2Λipm 0

 

 

 

 

 

ch(λ − z)β

 

 

 

 

 

 

[1

 

 

],

2λ ≤ z ≤ λ

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

β

 

 

 

 

ch(λ − 2λ1 )β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β =

 

1

Λi

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pm 0

 

 

 

 

Расходы газа Q(1)j

è Q(j

2 ) , порождаемые участком (φ –

/N, φ + /N)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

j

линии поддува, вычисляются по формулам (4.36) и (4.38). Учитываем, что согласно полученному решению

∂Π

|

 

 

 

=

1

{p 2

p 2

+ ε Re[(T (1 + w(1)) + w(1))eiϕ ]} + O(ε 2 ),

 

z= 2

λ1

0

 

z

 

 

m 0

1

1

2

 

 

 

 

2λ1

 

 

 

 

∂Π

|z= 2 λ +0 = −ε Re{th[(λ − 2λ1 )β ](βT +

 

2Λipm 0

)eiϕ ]} + O(ε 2 ).

 

 

 

 

 

 

 

z

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q(1)

=

 

κπ p2 H

3

 

{p 2

p 2

+ ε

N

 

sin(

π

) Re[(T (1 + w(1)) +

 

 

 

 

 

 

S

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

24λ1μaS2 N

 

m 0

1

 

π

 

 

N

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ w(1)

3( p 2

p 2 ))eiϕ j ]} + O(ε 2 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

m 0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.84)

 

 

 

 

2

H

3

 

 

N

 

 

π

 

 

 

 

 

 

 

 

 

 

2Λipm 0

 

 

Q( 2 ) =

κπ pS

0

{ε

 

sin(

) Re{th[(λ

2λ

)β

](βT +

)eiϕ

]} +

 

 

 

 

 

 

 

 

j 12μaS2 N

 

π

 

N

 

 

 

 

 

1

 

 

 

 

 

β

 

 

+O(ε 2 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для расхода газа, протекающего через питатель, справедлива формула

114

 

 

 

2

 

 

κ +1

 

κ p

A(ν )

 

 

M j =

α (

)κ −1

 

 

 

 

 

 

 

S

 

0

 

 

κ + 1

 

 

aS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+T

N

 

 

π

 

E

 

 

q( x

)

 

sin(

 

)

 

 

 

 

 

*

 

0

 

 

π

N

nNλ

1

q ( x

)

 

 

 

 

 

 

 

 

 

 

*

 

0

 

 

q* ( x0 ){1 + ε Re[(−ν +

(4.85)

)eiϕ ]} + O(ε 2 ).

Выполняя условие баланса расходов (4.40), получаем уравнение для определения Т, аналогичное уравнению (4.41), из которого находим

 

 

 

 

 

 

 

 

 

 

 

 

 

thλ β

 

 

 

 

 

 

 

π

 

 

 

( p2

p2 ) + w(1) 4

 

 

 

 

3

− ν

 

 

 

 

 

1

Λλ p

i

 

 

 

π

 

 

β

 

 

 

N sin(

 

)

m0

1

2

 

 

1 m0

 

 

 

 

 

 

 

 

 

 

 

(4.86)

 

 

 

 

 

 

 

 

 

T =

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

2)q

 

 

 

 

1+ w

+ 2λ βth((λ − 2λ )β) (x p

(x ) / q (x )

 

 

1

1

 

 

 

 

1

0

1

*

0

*

0

 

Воздействие на шип со стороны смазочного слоя определяется по формулам (1.91):

Fax + iFay = ∫∫ (p + iτ rϕ )dS,

Lax = ∫∫ r0τ rϕ dS.

(4.87)

S

S

 

В первом интеграле касательными напряжениями по сравнению с давлениями можно пренебречь, так как они имеют по сравнению с последними порядок ΛH0 / D.

Согласно формуле (4.43), давления в смазочном слое распределены по закону

p = pS ( Π0

+ ε

yeiϕ + yeiϕ

) + O(ε 2 ).

 

 

 

4 Π0

В связи с этим первую формулу (4.87) можно представить в

âèäå

λ

2π

y + ye

2iϕ

 

Fx + iFy = −2 pS r02 dz ( Π0 eiϕ + ε

 

)dϕ.

4 Π0

0

0

 

115

 

2π

Òàê êàê

eikϕ dϕ = 0, k = ± 1, ± 2, ... ,

 

0

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

 

 

 

 

 

λ

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fx

 

+ iFy = −

pS r02

 

 

dz.

 

 

 

 

 

 

 

(4.88)

 

 

 

 

 

 

 

 

 

2

 

 

Π0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитываем, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

y

 

 

 

2 λ1

 

A(z)

 

 

 

λ

 

A(z)

 

 

 

 

 

 

2 λ1

B(z)

 

 

 

λ

B(z)

 

 

 

 

dz = T[

 

dz +

 

dz] +

 

dz +

dz. (4.89)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pm 0

0

 

Π0

 

0

 

 

 

 

Π0

2 λ1

 

 

pm 0

 

 

 

 

0

 

Π0

 

 

 

 

 

2 λ1

 

 

Здесь при условии, что функции wk определены в виде

(4.78),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 λ1

A(z)

 

 

 

 

1

 

 

 

 

 

 

ζ dζ

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

ζ (1 − ζ )dζ

 

 

 

dz = 2λ1 [

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a1

 

 

 

 

],

 

 

 

( p 2

p 2 )ζ + p 2

 

 

 

( p 2

 

 

 

 

 

0

 

Π0

 

 

 

0

 

 

 

 

 

 

 

0

 

0

p 2 )ζ + p 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 0

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

m

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

λ

A(z)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

dz] =

 

 

 

 

th((λ − 2λ1 )β ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 λ1

 

p

m 0

 

 

 

β p

m 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 λ1

B(z)

 

 

 

 

 

1

 

 

 

 

 

ζ (1 − ζ )dζ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz = 2λ1a2

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p 2

 

 

p 2 )ζ + p

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Π0

 

 

 

 

 

0

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

λ

B(z)

 

 

 

 

 

 

2λi

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

dz =

[−λ + 2λ1 +

 

th((λ − 2λ1 )β ].

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2 λ1

p

m 0

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

Используем обозначения (4.79) для интегралов, входящих в правые части записанных выше равенств. Для вычисления несущей способности подшипника получаем формулу

116

F + iF = −επ p r2

{T[2λ

(M

 

 

+ a M

 

) +

 

1

th((λ − 2λ

)β ] +

10

11

 

 

x

 

y

 

 

S

0

 

1

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β pm 0

(4.90)

 

 

 

 

 

2Λi

 

 

 

 

 

 

 

1

 

 

 

 

 

 

+2λ

a

M

 

+

[−λ +

2λ

 

+

 

th((λ − 2λ

)β ]}.

 

 

β 2

 

 

 

1

2

 

11

 

 

 

 

 

1

 

 

β

 

 

1

 

 

 

Момент сопротивления вращения валу, создаваемый силами трения смазочного слоя, вычисляется по формуле

L / 2

2π

Lz = 2 dzτ rϕ r02dϕ.

0 0

Расход газа, как нетрудно показать, с точностью до членов порядка ε2, определяется по формуле (4.13) и не зависит от скорости вращения вала.

На рис. 4.3–4.5 представлено сравнение теоретических и экспериментальных нагрузочных характристик одно- и двухрядных радиальных газостатических подшипников.

3

 

 

 

 

2,5

 

 

 

 

2

 

 

 

1,1

 

 

 

 

cw1,5

 

 

 

1,1

 

 

 

 

1

 

 

 

0

0,5

 

 

 

0

0

 

 

 

 

0

0,2

0,4

0,6

0,8

Штриховой линией показаны экспериментальные данные

Ðèñ. 4.3. График зависимости безразмерной нагрузки Cw от относительного эксцентриситета ε при относительном давлении питания и числах сжимаемости Λ = 0., 1.1 (двухрядный подшипник)

117

 

50

 

 

 

 

 

 

 

 

1,7

 

45

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,7

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,3

 

35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

2,3

 

 

 

 

 

 

 

 

 

 

Fi,

25

 

 

 

 

 

 

 

 

3,6

град

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

3,6

 

15

 

 

 

 

 

 

 

 

5

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

ΛШтриховой линией показаны экспериментальные данные

Ðèñ. 4.4. График зависимости угла нагрузки от числа сжимаемости при ε = 0.1 и относительных давлениях питания pS = 1.7, 2.3, 3.6, 6 (двухрядный подшипник)

0,3

 

 

 

 

 

0,25

 

 

 

 

 

0,2

 

 

 

 

 

cw 0,15

 

 

 

 

0

 

 

 

 

 

0,1

 

 

 

 

0

 

 

 

 

2,2

 

 

 

 

 

0,05

 

 

 

 

2,2

 

 

 

 

 

0

 

 

 

 

 

0

0,1

0,2

0,3

0,4

0,5

 

 

 

ε

Штриховой линией показаны

 

 

 

экспериментальные данные

Ðèñ. 4.5. График зависимости безразмерной нагрузки от относительного эксцентриситета при относительном давлении питания pS = 4.6 è Λ = 0., 2.2 (однорядный подшипник). При расчете однорядного подшипника

использовались следующие значения параметрoв

L

= 2

,

D

l =

L

 

 

 

, D = 0,0508 м, d0 = 0,559 мм, h0 = 0,0157 м, N = 8

 

2

 

 

 

 

 

 

118

4.4. Вырожденное решение

Особого исследования требует задача расчета газостатических подшипников, когда Λ >> 1, так как в этом случае в уравнении Рейнольдса (4.69) при старших производных появляется малый параметр Λ–1 (для этого надо предварительно все члены уравнения (4.69) разделить на Λ). Ограничимся анализом вырожденного решения, отвечающего случаю Λ → ∞. Предполагая ограниченность давлений и их производных при всех Λ, из уравнения Рейнольдса (4.69) находим, что вне окрестностей границ (внешних и контуров питателей) давления должны быть при Λ → ∞ распределены так, что

( pH ) = 0.

∂ϕ

Откуда следует, что при Λ → ∞ вне рядов питателей давления в смазочном слое распределены по закону

p =

C(z)

,

(4.91)

 

 

H

 

в котором вид функции C(z) может быть установлен на основе следующих рассуждений.

Ïðè Λ → ∞ на линии поддува устанавливается постоянное давление pS, равное давлению в камере нагнетания. Объяснить это можно следующим образом. Учитывая вид параметра Λ, случай Λ → ∞ можно интерпретировать как случай, когда H0 → 0. Íî ïðè H0 → 0 давления на кромках питателей стремятся к pS, так как при этом расход газа Mj через j-тый питатель уменьшается и в пределе, согласно формуле (4.4), при H0 → 0

limq*(xj) = 0, j = 1, 2, …, nN.

Последнее равенство возможно лишь при pdj pS.

Òàê êàê ïðè Λ → ∞ на линиях поддува устанавливается давление pS, то газостатический подшипник при Λ → ∞ работает как газодинамический, у которого торцевые участки находятся в условиях осевого перепада давлений pS è pa, а межрядная зона работает как газодинамический подшипник длины L–2l, на торцах которого

119

поддерживается давление pS. Решение этих задач представлено в разделе 2, откуда следует, что на торцевых участках устанавливается давление

p = p

 

H 0

1 + 1, 5ε 2 (1

p 2 )

z

+ p 2

,

S H

 

 

 

1

2λ1

1

 

а в межрядной зоне

p = pS

H 0

1 + 1, 5ε 2 .

 

 

H

Вычисляя при этом несущую способность подшипника по формуле (4.87), находим

Fx = −π D2 pS

1 1 − ε

2

1 + 1, 5ε 2

[

4 (1 +

p1

 

)λ1 + λ − 2λ1 ],

 

 

 

 

 

 

 

 

 

2

 

 

ε 1 − ε 2

 

 

 

 

 

3

 

1 + p1

 

 

 

 

 

 

 

 

 

 

 

 

(4.92)

 

Fy = −

πμω D2 L( 1

 

1 − ε 2 )

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2H0

ε

 

1 − ε

2

 

 

 

 

 

 

 

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

Расход газа через подшипник при Λ → ∞ равен

 

Q =

κπ DH03 ( pS2

pa2 )

(1

3

ε 2 ).

(4.93)

12μlaS2

 

 

 

2

 

 

120

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