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

книги из ГПНТБ / Макаров Г.В. Уплотнительные устройства

.pdf
Скачиваний:
67
Добавлен:
23.10.2023
Размер:
8.39 Mб
Скачать

37.ВЯЗКОСТЬ НЕКОТОРЫХ ВЕЩЕСТВ

ВСЖИЖЕННОМ И ГАЗООБРАЗНОМ СОСТОЯНИИ

Для изучения 'Процесса истечения среды при наличии воз­ можности фазовых превращений существенное значение имеют

изменения вязкости

при разных

состояниях.

 

В работе [463 приведены данные, характеризующие изменение

динамической вязкости

для некоторых веществ, например N 2 , Н 2 ,

0 2 , Не,

воздух, при разных фазовых состояниях (газ,

жидкость).

Как видно из графиков, вязкость указанных веществ в жидком

состоянии превышает вязкость газа в 5—30 раз.

 

При

газообразном

состоянии

азота

при низких

давлениях

= 10-f-20 кгс/см2 )

с повышением

температуры вязкость увели­

чивается

в 2—3 раза.

 

 

 

 

 

При

давлениях р ^

150 кгс/см2

с

повышением температуры

до 200 К вязкость резко уменьшается, а при дальнейшем повыше­

нии до Т = 283 К остается

малой и

практически постоянной,

хотя ее значение будет выше,

чем при

р = 1-М0 кгс/см2 .

Отмеченные зависимости являются существенными для опреде­ ления расхода среды через уплотнение.

Зависимость вязкости газов от температуры приближенно вы­

ражается

формулой Сезерленда

 

з

 

 

 

 

 

 

 

 

 

 

 

где при 0° С н 1 атм для азота С =

114, цй

=

167 мкП;

для кисло­

рода

С =

131, г)о = 191 мкП.

Зависимость

вязкости

сжиженных

газов

от

температуры

(исключая

гелий)

выражается

формулой

 

 

 

 

г) =

 

с

 

 

 

 

 

 

 

At

т ,

 

 

 

где А

и С — константы.

 

 

 

 

 

 

38,

ОПРЕДЕЛЕНИЕ РАСХОДА

СЖИЖЕННОГО

 

 

 

ГАЗА

ЧЕРЕЗ

УПЛОТНЕНИЕ

 

С

целью

выявления

характера

истечения сжиженных газов

при возможности фазовых превращений рассмотрим величины расхода сжиженного и сжатого газа через щелевое уплотнение при неподвижном вале (п — 0 об/мин), получаемые при расчетах по различным методикам.

Определение расхода сжатого газа при изоэнтропийном истечении

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

186

Для

двухатомных газов при

/г =

1,40

 

 

 

 

G = 0,685Й2

YlT

к г

/ с '

 

где рх

давление

подводимого

к соплу

газа;

vx—удельный

объем подводимого

газа; й 2 — п о п е р е ч н о е

сечение

сопла.

Признаком критического режима

является неравенство р ^

^р'В, где р' — давление среды, в которую происходит истечение,

В= 0,258 для двухатомных газов и В = 0,484 для одноатомных

газов.

Критическое

давление р к р

= рр х .

 

Если р к р < ; р',

то имеет-место

некритическое истечение, если

Ркр >

Р ' критическое истечение. Расход зависит от

проходного

сечения Q 2 , увеличивается прямо

пропорционально

увеличению

зазора s в первой степени. Расход не зависит от длины уплотнителыюго кольца I.

В термодинамике при определении расхода сжатого газа через сопло считают, что имеет места истечение идеального газа, не об­ ладающего вязкостью и теплопроводностью.

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

В этом случае при радиальных зазорах s << 26 (где 6 — тол­ щина пограничного слоя) существенное влияние на расход жид­ костей оказывает трение вязкой среды и перенос тепла в процессе истечения.

Определение расхода через микрозазоры при ламинарном истечении

Определение расхода запираемой среды при ламинарном исте­ чении исходя из потерь напора на трение. Газ можно считать несжимаемым, если скорость движения его мала по сравнению со скоростью звука [14] . "

Практически внутри любой области газо-жидкостной системы движение среды определяется обычными уравнениями гидро­ динамики [111].

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

Если бы весь перепад давлений Ар =

р — ра расходовался

на преодоление сил трения в зазоре, тогда

 

Ар

=

Ар/,

 

 

где р а — давление в полости,

в

которую

происходит истечение.

В этом случае, согласно исследованиям Т. М. Башты,

ApfKdssk3

 

 

 

< ? У =

12 tiL

М / С >

 

187

и весовой расход Gy = Qyyp =

Qyyp = P l кгс/с, где Ар/ в кгс/м2 ;

d, s п L в м, 11 в кгс-с/м2 , 7 в

кгс/м3 .

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

Определение расхода при ламинарном истечении с учетом мест­ ных и скоростных потерь напора. Потери давления при проте­

кании среды через

зазор

 

 

 

 

 

 

 

Ap = p —P a = Apf + Ари + Ар в х + Лр в ы х ,

 

 

где Apf — перепад

давлений на длине

втулки,

расходуемый на

преодоление сил трения;

Арц

перепад давлений, расходуемый

на создание скоростного напора; Л р в х , Д р в

ы х

перепад давлений,

расходуемый на

преодоление местных

сопротивлений

на

входе

в зазор и выходе из зазора.

 

 

 

 

 

 

На основе этой зависимости имеем

 

 

 

 

 

~2~ (1 +

Sbx + W ) " 2 +

J f j J

L -

(Р -

Ра) =

О-

(163)

Решая квадратное уравнение, определим среднюю скорость

течения в зазоре

и.

 

 

 

 

 

 

 

Секундный расход

G =

uFyp.

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение

расхода

запираемой среды

 

 

 

при турбулентном истечении

 

 

 

Определение

расхода

исходя

из потерь напора на трение.

При условии Ар =

р — ра = Ар/, расход

жидкости

 

 

 

 

 

 

Х1уР

 

 

 

 

где К определяется

по формуле Блазиуса

 

 

 

 

* 0,3164 n

us

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

Определение расхода с учетом местных и скоростных потерь напора. Ранее было получено

+$)+hf

+1* =

Y - a 2

2g

где и — средняя

скорость потока в зазоре;

и0 — скорость перед

входом в зазор.

1

 

На основании проведенных расчетов можно сделать следу­

ющие выводы.

 

 

1. Кривые изменения вязкости сжиженного газа с изменением

температуры и давления имеют различный

характер при низких

иповышенных давлениях. В соответствии с ними будут изменяться

изначения расхода.

2.При движении сжиженного газа по длине зазора происхо­ дит одновременное повышение температуры и падение давления. Благодаря этому в диапазоне малых зазоров можно вести расчет

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

39.УРАВНЕНИЯ ПЕРЕПАДОВ ДАВЛЕНИЯ

ПО ДЛИНЕ ЗАЗОРА

И

ИЗМЕНЕНИЯ СКОРОСТИ

ПО ВЫСОТЕ ЗАЗОРА

ПРИ

ЛАМИНАРНОМ ИСТЕЧЕНИИ

Рассмотрим неизотермическое течение в микрозазорах несжи­ маемой среды с переменной вязкостью. Условно считаем, что

уплотнительное кольцо и вал абсолютно_жесткие и идеально гладкие. Температура ра­

бочей среды изменяется по

 

 

 

 

 

 

пути течения,

и

при

этом

 

'//////////////////////А

 

возможен

переход

части

 

 

 

т-<£т а>

 

 

или всей движущейся сре­

 

 

 

ды из

 

жидкой фазы

в га­

—»-

о -—

^

-P-dp

— "о

 

ШМ

 

зообразную.

Радиальные

 

т

 

 

 

 

зазоры

s <

28,

где б

0

 

 

 

У*

X

 

 

 

 

 

толщина

 

пограничного

 

 

 

 

 

 

слоя.

Рассматриваем

слу­

 

 

 

 

 

 

чай течения среды в малом

 

 

 

 

 

 

зазоре, при котором

зазор

Рис.

106. Схема

к

расчету перепадов дав­

s мал по сравнению с дру­

гими

размерами

щели

L

 

лений

 

и d.

Плотность

рабочей

среды принимаем

постоянной и

равной

плотности при давлении р перед входом

в

 

зазор.

 

Как показали проведенные эксперименты, в газ превращается

только

4%

всего

протекающего

в щели жидкого газа.

 

С изменением температуры рабочей среды на 60° С плотность изменяется на 0,5—4%, поэтому погрешность расчета давления

при у = const будет мала

[14, 36] .

 

 

Определение перепада

давлений по пути движения

среды

 

в узких щелях, обусловленного потерями энергии на

трение

-

слоев среды, без учета инерции жидкости произведено Т. М. Баш-

той [ 6 ] .

 

 

 

189

Условие равновесия элементарного объема жидкости, нахо­ дящейся в нижней части зазора (рис. 106),

dp dy — (т — dx) dx — т dx.

Откуда

dp

_

dx

dx

 

dy

Имея

в виду, что

т =

 

(r\vx),

при

1] постоянном

по

высоте

зазора

у,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvjc

dx

 

 

 

d2vx

 

 

 

 

 

 

 

 

dy

dy

 

 

1

dy 2

 

 

Выражение

потерь

на жидкостное трение

 

 

 

 

 

 

 

 

dp

_

d?vx

 

 

 

 

 

 

 

 

 

 

dx

1

dy*

"

 

 

 

Имея

в виду, что s мало по сравнению с L и d, принимаем

 

 

 

 

 

 

 

- £ - = <>•

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

Ввиду

малости

зазора

s считаем

л

и

 

не зависящими

от у.

Дважды интегрируя по у выражение

находим

изменение

скорости vx

по

высоте

зазора s.

 

 

 

 

 

 

Исходное

выражение

 

 

 

 

 

 

 

После

двукратного

интегрирования

имеем

 

 

 

 

 

 

 

dx

2

 

 

 

 

 

 

 

Определим

коэффициенты С х

и С 2

из граничных условий:

 

1)

при

у

=

0 uA . = 0;

 

 

 

 

 

 

 

2)

при

у

=

s vx = 0.

 

 

 

 

 

 

 

Из

первого

условия

получаем

С 2

= 0.

 

 

Второе

условие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp

s2 _

 

р

_

 

 

 

ф S

Подставляя коэффициенты С\ и С 2 , имеем

dp о2

. dp s

190

Откуда

и

Расход среды через данное сечение кольцевой щели S S

Qy = j vxnd dy = nd J vx dy,

о0

где

 

vxdy

=

— -^--^ry(y

s)dy.

 

 

Расход на

единицу длины окружности кольцевой

щели

 

 

 

-

 

f v

du -

 

— L .

s3

 

 

 

 

 

 

о

 

 

 

 

 

 

 

Следовательно, расход среды через данное сечение кольцевой

 

щели

 

 

 

 

 

 

 

 

 

 

 

 

п — n d

d p

"3 G y

 

 

 

 

V

y

12г|

dx

 

у

 

 

 

и перепад

давлений

 

 

 

 

 

 

 

 

 

 

 

J f L

=

Q

^

nds3

=

А 1 | П _ ,

(165)

 

 

 

 

 

^ у

 

 

ynds3

'

v

7

где у = pg; Gy — весовой расход в единицу времени. Часть сжи­ женных газов при течении в зазоре превращается в пар, в этом случае в качестве условия неразрывности потока должно обеспе­ чиваться сохранение постоянства весового расхода в отдельных сечениях

°у = yPQy = Yp"cP ^ = const,

где ур — плотность среды при данном давлении (для пара ур =

^ р ^ а т и ^ ) *

Имея в виду условие сплошности потока, из уравнения (165) следует

• ( - 5 r ) i = c o n s t -

< 1 6 6 >

При течении сжиженных газов в микрозазоре и повышении их температуры часть сжиженных газов может превращаться в пар, при этом существенно изменяется коэффициент динамической вяз-

191

кости г). Поэтому, чтобы выявить влияние отдельных факторов на расход, необходимо провести исследование движения жидкости при переменных значениях вязкости, считая г) = f (Т, р).

40.УРАВНЕНИЕ ПРИТОКА ТЕПЛА

ПРИ НЕИЗОТЁРМИЧЕСКОМ ТЕЧЕНИИ СРЕДЫ

ВКОЛЬЦЕВОМ МИКРОЗАЗОРЕ

Уравнение изменения энергии в фиксированном элементар­ ном объеме (не изменяющемся во времени) вязкой жидкости запи­ сывается в следующем виде:

/ дЕ , дЕ . дЕ . дЕ \ „ .

+ т Н ^ ( * £ ) + £ ( т ) + ^ ( ^ ) М

(1

где D — дпссипатнвная функция, представляющая собой механи­ ческую работу сил вязкости, выделившуюся необратимо в виде тепла в единице объема жидкости за единицу времени,

+ ( * + - & ) ' +

( * + - £ ) ' ] - т * " =

8 — дивергенция вектора скорости при переходе от данной фикси­ рованной точки пространства к другой

а _ ^>х_

,

dvy_ ,

dvz_.

дх

~*

ду

dz ' -

Е — внутренняя энергия единицы массы; X — коэффициент тепло­ проводности жидкости; е — приток тепла за единицу времени в единице объема, происходящий по причинам, отличным от

теплопроводности;

А =

термический эквивалент

работы

в ккал/кгс-м; Т — температура

фиксированной частицы

с по­

стоянной массой в

К-

 

 

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

192

После ряда преобразований и упрощений [36] уравнение притока тепла для неизотермического течения в гидродинамиче­ ской теории смазки записывается в следующем виде:

 

( дс£_ ^_ п

дс£_

L п

дс-0Т

А

 

 

 

 

1

д

 

 

 

где р — давление жидкости; t—• время.

Уравнение притока тепла (168) при его осреднении по толщине

слоя для

неизотермического течения с учетом закона охлаждения

Ньютона

и интегрирования по у в пределах толщины слоя при­

нимает следующий вид:

 

 

 

 

[goJp\Oxdy^J+-jjr[gc0Tp

J

vzdy^j + gcvTpw

+ а(Т-Те)

+

+

* 4 " (goDTP)=-Ap

} в ^ +

Л т , | [ ( ^ . ) а +

( ^ - ) 2

] ^ ,

оо

где w

скорость перемещения рабочих поверхностей в направле­

нии оси оу.

В нашем случае w =

0 и /г =

s. Примем

vz = 0 (при

п — 0

об/мин).

 

 

 

 

 

При установившемся движении уравнение притока тепла запи­

шется

в

следующем

виде:

 

 

 

д

gcvTp

\ vxdy\

 

+ а (Т - Тс)

= -Ap\Qdy

+ Ar]\

) 2 dy,

 

 

дх \

 

о

)

о

 

о

 

 

 

 

 

 

 

 

(169)

где Тс — температура поверхностей, разделяемых слоем жидкости; а—коэффициент теплообмена между этими поверхностями и средой.

Подготовим необходимые зависимости для использования их

s

в уравнении (169). Определим выражение | G dy из рассмотрения

уравнения сплошности и выражения 6.

о

 

 

 

Уравнение

сплошности

 

 

 

 

 

 

dp

I

dJPVx)

1 д

(pvy) I

d (рог ) _

n

 

 

dt

 

dx

*

dy

dz

' '

При

= 0

vz =

0

н - g - =

0;

 

 

 

 

 

 

 

д (P%)

i

д (pvy)

Q_

 

 

 

 

 

dx

'

du

'

 

 

 

 

 

 

 

dy

 

 

193

дх

~

v*

дх

_ | "

^

5л:

'

д(/

~

иУ

ду

^

^

ду

'

Ох-^- + рв = 0;

Qdy = —

-j^rvxdy\

Принимая жидкость несжимаемой, имеем

=

0 и | 0 dy — 0.

Ранее было получено

 

о

 

 

5

 

 

J vxuy 1 2 т ) ^ ь — n d

p g n d

,

где p = - у - плотность среды, соответствующая давлению на входе в зазор.

°* = - - ^ - £ - 0 ( 5 - 0 ) ;

~2т] dx- l s

Ранее было получено

dp =

, l2Qyr\

dx

я as3

194

Имея в виду, что Qy

получаем

12G-;

тр у -

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

д (

cvTG

 

 

 

 

 

дх

 

 

 

 

 

 

или

 

 

 

 

 

 

иУС"

dx

= nd (p^r)

^

-a(T—Tc)nd

 

2

s 3

 

 

или

 

 

 

 

 

 

GycvdT--

 

ltd (pg)2 s 3

 

and

(T — Tc) dx.

(170)

 

 

 

 

 

 

Количество тепла, сообщаемое протекающей среде в единицу времени

 

 

dQ = Gycv dT.

(171)

Окончательное

уравнение притока тепла будет

 

 

dQ--

- ^ - a n d ( T - T c ) dx.

(172)

Принятые размерности Q в ккал/ч; л. в кг-ч/м2 ;- Gy

в кг/ч; -у в кг/м3 ;

s, х и d в м; 7

К; С0

в ккал/кг-град.

 

41 . ОСОБЕННОСТИ НЕИЗОТЕРМИЧЕСКОГО ТЕЧЕНИЯ СЖИЖЕННЫХ ГАЗОВ В КОЛЬЦЕВЫХ МИКРОЗАЗОРАХ ПРИ ВРАЩЕНИИ ВАЛА

Изменение давления среды при наличии поступательного и вращательного движения ее частиц

Л х + ж й г -

Согласно ранее приведенным выводам, при наличии малой величины у по сравнению с х, и, принимая г\ не изменяющимся в направлении оси оу, было получено

- # - = 0.

ду

195