Теоретические основы работы газостатических опор. Монография
.pdf
Учитываем теперь, что в соответствии с формулой (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πγ q′ xΛ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+∫i∞ sσ eiσ dσ.
c–i∞
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
На рис. 5.3 представлены результаты расчета несущей способности коэффициентов жесткости и демпфирования упорного шестисегментного подшипника.
Как показывают представленные графики, коэффициент жесткости подшипника возрастает с ростом частоты осевых колебаний подвижного элемента, а коэффициент демпфирования снижается, что означает, что для демпфирования осевых колебаний в газостатической опоре необходимо предусматривать внешние демпферы.
