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

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

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

Учитываем теперь, что в соответствии с формулой (5.69) имеют место равенства

Mσ

∫∫ L(gσ , h03 )ξl dΩ = Cmσ ∫∫ h03 ξm ξl dΩ,

Ω m=1 Ω

∫∫ L(sΛ , 3h02hσ )ξl dΩ = – 3∫∫ h02hσ sΛ ξl dΩ,

Ω

Ω

 

 

Mk

 

∫∫ L(gk , fr , h03 )ξl dΩ = Cm(k ) ∫∫ h03 ξm

ξl dΩ

Ω

m=1 Ω

 

∫∫ h03 fk ξl dΩ .

Ω

Нетрудно показать также, что

∫∫ (au)ξl dΩ = ∫∫ au ξl dΩ ,

Ω Ω

вследствие чего справедливы равенства

 

 

Λ∫∫

 

 

2h s

Λ

+ h g

σ

 

 

 

 

 

 

 

 

 

σ

0

u

ξl dΩ =

 

 

 

 

sΛ

 

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mσ

 

 

 

h0

 

 

= –2Λ∫∫ hσ sΛ u ξl dΩ Cmσ ∫∫

 

ξm u ξl dΩ,

 

 

 

 

Ω

 

 

 

 

 

m=1

 

Ω

 

sΛ

 

 

gk + fk

 

 

 

Mk

 

 

 

 

h0ξm

 

Λ∫∫

h0

u

ξl dΩ = Cm(k ) Λ∫∫

u ξl dΩ

 

 

Ω

 

sΛ

 

 

 

 

m=1

 

 

 

Ω

sΛ

Λ∫∫ h0sfk u ξl dΩ.

ΩΛ

161

Очевидны также равенства

 

 

σ ∫∫

2hσ sΛ + h0 gσ

ξl dΩ = σ ∫∫ 2hσ sΛ ξl dΩ +

 

 

 

 

 

 

 

 

 

 

 

 

 

Ω

 

 

 

sΛ

 

 

 

 

 

 

 

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mσ

 

 

 

 

 

hσ ξm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+Cm(σ )σ

∫∫

ξl dΩ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1

Ω

 

 

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk + fk

 

 

 

 

 

Mk

 

 

 

 

 

h0ξm

 

 

 

 

 

 

 

 

σ ∫∫ h0

ξl dΩ = Cm(k )σ ∫∫

ξl dΩ +

 

 

 

 

 

 

 

 

 

 

 

 

Ω

 

 

sΛ

 

 

 

 

 

 

m=1

 

 

Ω

 

 

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

∫∫

h0 fk

ξl dΩ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ω

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

 

Вводим обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

alm = ∫∫{(Λ

h0ξm

u h03 ξm ) ξl

σ h0ξm

ξl } dΩ,

 

 

 

 

 

 

 

 

 

 

 

Ω

 

 

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

 

sΛ

 

 

 

l

∫∫

{(

 

0 σ

Λ

 

 

 

σ

 

Λ

u

)

l

 

 

 

σ

Λ l }

(5.102)

b(σ ) =

 

 

3h

2h s

 

2Λh s

 

 

 

ξ + 2σ h s

ξ

dΩ,

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(k ) =

∫∫

{(h3 f

k

Λ

h0 fk

u) ξ + σ

h0 fk

ξ } dΩ.

 

 

 

 

 

 

l

 

0

 

 

 

 

 

sΛ

 

 

 

 

 

l

 

 

 

 

 

l

 

 

 

 

 

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sΛ

 

 

 

Системы уравнений (5.100)–(5.101) в принятых обозначениях записываются в виде

Mσ

almCm(σ ) = bl(σ ) , m=1

Mk

almCm(k ) = bl(k ) , m=1

l = 1, 2, ..., Mσ ;

l = 1, 2, ..., Mk

,

k = 1, 2, ..., N.

(5.103)

 

162

После решения систем уравнений (5.103) можно считать процесс установления приближенных значений функций gσ è gk по формулам (5.99) законченным. Безусловно, на точность расчетов влияет число функций Галеркина ξm, удерживаемых в рядах (5.99). Окончательный выбор значений этих чисел Mσ è Mk целесообразно провести после определения интегральных характеристик смазочного слоя (жесткости, демпфирования и т. п.), оценив влияние числа функций Галеркина в разложениях (5.99) на точность определения этих характеристик.

Установив вид функций gσ è gk , можно перейти к определению значений коэффициентов *Gk и величин динамических составляющих давлений в «карманах» и на выходных кромках питателей xˆσj è xσj, для чего служат последние три системы условий на питателях, входящих в уравнения (5.95), содержащие величины xˆσj è xσj.

В указанные условия необходимо подставить значение функции sσ, определяемое по формуле (5.96), и учесть, что на контурах питателей Гj функции gσ è gk обращаются в нуль. Получаем следующие уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ jk Gk = xσ j ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J jk Gk = 2πγ j q*(xΛj )xσ j =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xσ j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= – 2πσ

K + h

K

 

 

+ h

K

 

xΛj

V

,

 

 

 

 

Sj )

 

 

 

Sj

 

 

 

 

 

 

( Vj

0 j

 

 

xΛj

σ j

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

= 2πγ j h0 j q* (xΛj

 

 

) xΛj

 

 

 

 

 

 

 

 

 

 

 

(5.104)

J jk Gk

xΛj

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xΛj )

xΛj x

 

x

Λj

xσ j

 

xσ j

 

h

 

 

 

 

 

 

 

 

 

σ j

 

 

 

 

+

+

σ j

Vj

.

q*(xΛj xΛj ) q* (xΛj

 

 

xΛj

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2xΛj

 

h0 j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

163

Здесь

 

1

 

 

 

fk dΓ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ jk =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πd

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

Γ j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk + fk

 

 

 

 

 

gk

 

fk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

J jk = σ ∫∫ h0

 

 

 

dΩ h0

 

 

 

+

 

 

 

dl

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sΛ

 

 

 

 

 

v

 

 

v

 

 

 

 

 

Ω j

 

 

 

 

 

 

 

l j

 

 

 

 

 

 

 

 

 

h0

g

k

 

+ f

k

uvdl,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.105)

 

 

 

 

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

 

l j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2hσ sΛ + h0 gσ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vj = σ ∫∫

dΩ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l j

 

 

 

 

 

sΛ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

s

 

 

 

 

2h s

 

+ h g

 

 

 

 

h02

 

 

 

 

+ Λ

Λ

σ

 

 

 

h0

 

 

 

σ

+ 3hσ

 

 

Λ dl

 

σ

 

 

0

uvdl.

 

 

 

v

 

 

 

 

 

 

sΛ

 

 

l j

 

 

 

 

 

v

 

 

l j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Напомним, что если не учитывается эффект двойного дросселирования в питателях с «карманами», то последние N уравнений полученной системы должны быть отброшены. Необходимо также

принять xˆΛj = xΛj, xˆσj = xσj.

В этом случае исключением величин xσj из первых 2N уравнений системы (5.102) получаем систему линейных алгебраических уравнений для определения величин Gk:

 

 

N

I j Φ jk )Gk

 

 

 

 

 

 

 

 

(J jk

= Bj ,

(5.106)

 

 

k =1

 

 

 

 

 

 

 

 

I = 2πγ qxΛj

(K + h K

 

)

πσ

,

 

 

 

 

 

j

j * (

) Vj

0 j

Sj

 

 

 

 

ãäå

 

xΛj

(5.107)

Bj = –2πσ hσ j KSj

xΛj Vj .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определив значения величин Gk из этих уравнений, из уравнений (5.104) находим

N

 

xσ j = Φ jk Gk .

(5.108)

k =1

164

Если же рассматриваются питатели типа «кольцевая диафрагма», то необходимо решать первые N уравнений совместно с последними N уравнениями системы (5.104), положив xˆΛj = 1, xˆσj = 0.

При этом исключением величин xσj, легко получаем систему линейных алгебраических уравнений для определения величин Gk:

 

 

 

N

 

I *Φ

 

G

= B* ,

 

 

 

 

J

 

(5.109)

 

 

(

jk

 

j

 

 

jk ) k

0

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

I

*

= 2πγ

р

 

 

q′ x

,

 

 

 

 

j

 

 

j

0 j

* ( Λj )

 

 

ãäå

B* = 2πγ h

 

q′

x

V .

(5.110)

 

 

j

 

σ j

 

*

(

Λj )

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае питателей с карманами для учета эффекта двойного дросселирования необходимо решать всю систему 3N уравнений (5.104). При этом из первых 2N уравнений находим

 

N

 

 

 

 

 

xσ j = Φ jk Gk ,

 

 

 

 

k =1

 

 

 

 

 

 

1

 

N

 

(5.111)

xσ j =

 

 

 

 

J jk Gk

Bj

.

 

I

 

 

 

 

j

 

k =1

 

 

 

 

 

 

 

 

 

После подстановки этих значений величин xσj è xˆσj в последние N уравнений системы (5.104) получаем систему линейных алгебраических уравнений для определения величин Gk:

N

 

 

 

 

*

 

 

πγ j h0 j q* (xΛj

 

 

 

*

 

 

 

 

 

 

 

I j xΛj

 

 

xΛj )

 

 

J jk

1+

 

 

 

 

 

 

 

 

 

I j

Φ jk Gk

=

I j xΛj

 

I j xΛj

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.112)

 

 

*

xΛj

 

πγ j h0 j q* (xΛj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

I j

xΛj )

B

 

+ B*.

 

 

 

 

 

 

 

 

 

 

 

 

 

I j

xΛj

I

j

xΛj

 

 

 

j

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры Ijk, Tj, Tj* , Bj è B*j , входящие в эти уравнения, определяются по формулам (5.105), (5.107) и (5.110).

165

После определения величин Gk из уравнений (5.112) по формулам (5.111) вычисляются величины xσj è xˆσj.

Определив теперь по формуле (5.96) функцию sσ, можно найти ее оригинал по формуле

st = c+isσ eiσ dσ.

ci

5.8.Динамические характеристики смазочного слоя упорных сегментных подшипников

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

образом:

.

 

 

 

 

W = W0 G H D H,

(5.113)

ãäå H – малое смещение подвижного элемента,

.

H – скорость этого

смещения, G – коэффициент динамической жесткости смазочного слоя, D – коэффициент демпфирования.

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

H = H0 (1 – εsin ωt) = Hh,

(5.114)

ãäå h = 1 – εsin ωt.

 

В этом случае

 

.

 

H = H0 εsin ωt, H = – H0

εωcos ωt.

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

W = W0 + Wsεsin ωt + Weεcos ωt.

(5.115)

Из сопоставления формул (5.114)–(5.115) следует, что

G =

Wε

,

D =

We

.

(5.116)

 

 

 

H0

 

H0ω

 

166

Определение динамических коэффициентов жесткости G и демпфирования D упорного сегментного газостатического подшипника в зависимости от частоты сжатия смазочного слоя ω* позволяет исследовать поведение системы «смазочный слой–подвижный элемент» и условия, при которых она устойчива.

Исследуем собственные колебания системы «смазочный слой–подвижный элемент». Уравнение собственных колебаний подвижного элемента имеют вид

m H + D ω

H + G(ω

) H = 0.

(5.117)

( * )

*

 

Ищем решение этого уравнения в виде

H = Ce(b+iω* )t .

На основе уравнения (5.117) получаем уравнение частот

m b + iω t 2

+ D(ω

) b + iω t

)

+ G(ω

) = 0,

(5.118)

(

* )

*

(

*

*

 

из которого путем разделения вещественной и мнимой частей получается система уравнений для определения b è ω*:

m(b2 ω 2t2 ) + D(ω

)b + G(ω

) = 0,

 

 

*

*

*

 

 

(5.119)

2mb + D(ω

) = 0.

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

Собственные колебания будут затухающими, а система устой- чивой, если b < 0, ò. å. ïðè D > 0, как это следует из последнего уравнения.

Ïðè b < 0 колебания будут развивающимися (система неустой- чива), при этом D*) < 0. Следовательно, граница области устойчи- вости (b = 0) может быть определена из условия D*) = 0.

Пусть ω* = ωкр – положительный корень этого уравнения (если он существует). Так как в этом случае b = 0, òî ïðè ω* = ωкр из первого уравнения системы (5.119) можно найти критическую массу подвижного элемента, определяющее границу области устойчивости:

mкр = G (ωкр ). ωкр

167

В области устойчивости D*) > 0. Если существует ωкр > 0, то, как показывают расчеты, в области устойчивости ω* > ωкр.

Из уравнений системы (5.119) получаем

 

 

 

 

2

 

 

dG

1

 

dD

 

 

 

 

 

dm

 

 

2m ω*

m

 

2

D(ω* )

 

 

 

 

 

 

 

=

dω*

dω*

.

(5.120)

 

 

 

 

 

 

 

G(ω* ) – 2mω*2

 

 

 

 

dω*

 

 

 

 

Íà

границе

области

устойчивости ω* = ωкр,

Dкр) = 0,

G ω

= m ω 2 .

 

 

 

 

 

 

 

 

 

 

 

(

кр )

кр кр

 

 

 

 

 

 

 

 

 

 

 

Тогда из уравнения (5.120) следует, что в области устойчивости m < mêð, òàê êàê ïðè ýòîì

 

 

 

 

dm =

2mω*

 

dG

 

 

 

 

 

 

 

ω*2

dω*

< 0.

 

 

 

 

 

 

dω*

 

 

 

 

 

10–6, Н/м

 

 

 

 

10–2, Н

180

10–3, Н . c/м

 

 

360

 

 

 

 

36

 

 

 

 

300

W0

 

 

30

 

150

 

 

 

5000

 

 

 

 

 

 

 

 

 

 

 

 

w* = 50

 

 

240

 

 

 

 

24

 

120

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

18

 

90

500

 

 

180

 

 

 

2500

 

 

 

 

 

 

 

 

 

 

 

 

 

 

120

 

 

 

 

12

 

D60

2500

 

 

 

 

 

 

 

 

 

 

500

 

 

 

 

 

 

 

 

 

 

60

w*

= 50

 

 

6

 

30

 

 

 

 

 

 

 

 

 

5000

 

 

 

 

G0

 

 

 

 

 

 

 

0

20

Н0

40

60

80, мкм 0

20 Н0

40

60, мкм

Рис. 5.3. Зависимость от величины зазора Ho несущей способности Wo и динамических коэффициентов жесткости G и демпфирования D упорного сегментного газостатического подшипника: RH = 0.113, RB = 0.044, d = 0.0005, r1 = r2 =

= r3 = r4 = 0.08 ì, φv = 60°, φ1 = 13,5°, φ2 = 24,5°, φ3 = 35,5°, φ1 = 46,5°, pS = 4,19 . 105 Па, μ =18,15 . 10–6 Па с, aS = 341 м/ с, κ = 1,4, α = 0,8, n = 6

168

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

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

ЛИТЕРАТУРА

1.Руководящий технический материал. Расчет радиальных газостатических подшипников турбомашин атомной энергетики: РТМ 108.129.101-76. Издание официальное / Данильченко В. Ф., Лу- чин Г. А., Мордвинкин В. А., Снопов А. И. и др. – Л.: РИО НПО ЦКТИ, 1977.

2.Руководящий технический материал. Расчет и проектирование кольцевых газостатических подшипников турбомашин атомных энергетических установок: РТМ 108.129.104-78. Издание официальное / Данильченко В. Ф., Лучин Г. А., Мордвинкин В. А., Снопов А. И.

èäð. – Ë.: ÐÈÎ ÍÏÎ ÖÊÒÈ, 1980.

3.Руководящий технический материал. Уплотнения газовые плавающего типа валов центробежных и осевых компрессорных машин. Расчет и проектирование: РТМ 108.129.15-81. Издание официальное / Даниль- ченко В. Ф., Лучин Г. А., Снопов А. И. и др. – Л.: РИО НПО ЦКТИ, 1982.

4.Руководящий технический материал. Подшипники упорные сегментные газостатические валов центробежных и осевых компрессоров. Расчет и проектирование: РТМ 108.129.105.84 Издание официальное / Данильченко В. Ф., Иванов А. Н., Лучин Г. А., Снопов А. И.

èäð. – Ë.: ÐÈÎ ÍÏÎ ÖÊÒÈ, 1985.

5.Дидов В. В. Динамика роторов судовых турбомашин на подшипниках с газовой смазкой: монография. – Владивосток: Изд-во ДВГТУ, 2005.

6.Константинеску В. Н. Газовая смазка: пер. с рум. – М.: Машиностроение, 1968.

7.Лучин Г. А., Пешти Ю. В., Снопов А. И. Газовые опоры турбомашин. – М.: Машиностроение, 1989.

8. Пешти Ю. В. Газовая смазка. – М.: Изд-во МГТУ им. Н. Э. Баумана, 1993.

9.Подшипники с газовой смазкой / под ред. Н. С. Грессема и Дж. У. Пауэлла: пер. с англ. под ред. С. А. Харламова. – М.: Мир, 1966.

10.Пинегин С. В., Поспелов Г. А., Пешти Ю. В. Опоры с газовой смазкой. – М.: Наука, 1977.

170

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