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

книги из ГПНТБ / Жаров Г.Г. Судовые высокотемпературные газотурбинные установки

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.85 Mб
Скачать

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

Характер

зависимости

термического

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

контакта

от силы сжатия почти не зависит

от свойств окружающей

среды.

іІаг-ч-кІкіж

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і/а.-101,м'-ч/кдж

 

 

 

 

0,20

 

 

 

 

 

 

 

 

 

 

 

1,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,15

 

 

 

 

 

 

 

 

 

І,О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

0,10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

W0

 

 

200

р,Мн/м

300

 

 

 

323

 

373

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.

101.

Зависимость

термического

со­

 

Рис.

102.

Зависимость

терми­

противления

от

давления

и

материала.

 

ческого сопротивления

от

тем­

I —

непосредственный

 

контакт;

2

про­

 

пературы.

 

 

 

 

 

 

кладка из

меди;

3

прокладка

 

из

латуни.

 

1

сталь

гладкая

 

поверхно­

 

 

 

 

 

 

 

 

 

 

 

 

сть;

2 сталь

 

шероховатая

 

 

 

 

 

 

 

 

 

 

 

 

поверхность.

 

 

 

 

 

 

Это говорит о том, что большая часть

теплового

потока

проходит

через

газовые

прослойки.

 

 

 

 

 

 

 

 

 

 

 

 

 

В работе [77] проводились исследования по контактному сопро­

тивлению

некоторых

пар

металлов. Для

всех пар

характерно, что

с увеличением давления и чистоты обработки

поверхности

терми­

на,-/ff^-v-tf/fte

 

 

 

 

 

 

1/а-Ю!г/кдж

 

 

 

 

 

 

 

0.15

 

 

 

 

 

 

 

 

 

0,015

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,050

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,025

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

30

 

. .40

 

О

3

100

 

 

200

 

 

300

 

 

 

 

 

 

 

 

 

 

 

р,Мн/мг

 

 

 

 

 

 

р,Мн/м2

 

 

 

 

 

 

 

 

 

 

Рис.

103.

Зависимость

 

термиче­

 

Рис.

104.

 

Термическое

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

ского

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

от

газовой

 

для

различных

материалов.

 

 

 

прослойки

и давления.

 

 

 

 

 

1 — ЭИ2И;

2

ЭЖ2;

3

сталь

 

30.

 

 

 

 

 

 

 

 

 

 

 

 

 

1 — С02;

2 — воздух; 3

Я 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

ческое сопротивление уменьшается. В наибольшей степени терми­ ческое сопротивление уменьшается для пар сталь 30—сталь 30 и медь—медь, в наименьшей — для таких пар, как латунь—латунь и дюралюминий—дюралюминий.

Исследованиями установлено, что скачок температур на стыке двух поверхностей имеет существенное значение. Так, при удельных потоках порядка 580 000 вт/(м*-К) для стали 30 при обработке

поверхности

V V 6 температура изменяется от 673 до 373 К, а для

пары сталь

30медь при тех же исходных данных — от 443 до

277 К. Поэтому при расчетах температурных полей и температурных напряжений необходимо учитывать контактное сопротивление и пре­ небрегать им только в том случае, если тепловой поток менее 23 200 вт/(м2-К).

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

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

В общем можно сделать вывод о том, что термическое сопротивле­ ние контакта обусловливается:

— значительно меньшей теплопроводностью среды (меди, газа

ит. д.) по отношению к металлу деталей;

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

стей;

уменьшением площади проводника тепла.

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

§43. Теплообмен в цилиндрических капиллярных каналах и хвостовиках рабочих лопаток

В охлаждаемых газовых турбинах встречаются каналы диаметром менее 1 мм. В таких каналах теплообмен не­ сколько отличен от общих закономерностей теплообмена в трубах и каналах общего диаметра. Чаще всего каналы диаметром 1 мм расположены в местах стыковки отдельных деталей. Особое значение эти каналы имеют в зоне хвостовиков рабочих лопаток. Они служат протоком охлаждающего воздуха, который отводит от лопатки опре­ деленное количество тепла с целью поддержания заданной темпера­ туры хвоста. Поэтому знание характера теплообмена в этой области для расчета температурных полей является важным и необходимым. Передача тепла в зоне хвостовиков рабочих лопаток идет двумя путями: первый — через охлаждающий воздух, протекающий через

капиллярные каналы, и второй — через контактирующиеся

поверх­

ности хвостовика лопатки и ротора.

 

 

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

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

в каналах.

В работе [107], где проведено

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

0,28—0,622 мм, предполагается,

что ламинарный

режим

течения

прекращается при числе Рейнольдса, близком к 2000. При этом в ка­ налах меньших диаметров турбулизация потока начинается при

меньших значениях чисел Рейнольдса. Все эксперименты проводи­ лись при изотермическом течении жидкости и, соответственно, без турбулнзации потока конвективными силами. В работах, проведен­ ных Институтом теплоэнергетики АН УССР, исследовались капил­ лярные щелевые каналы сложного профиля. Опыты проводились с воздухом при изменении М на выходе из канала до 0,15. Воздух подавался давлением до 1 Мнім", предварительно очищенный. Испытывались два типа профилей монтажных зазоров елочных хвосто­ виков рабочих лопаток: с острой входной кромкой и закрученной входной кромкой. Диаметры каналов не превышали 0,25 мм. Вход­ ная закрученная кромка имела радиус 3 мм. Для определения влия­ ния формы поперечного сечения на теплоотдачу одновременно были испытаны цилиндрические трубки близких гидравлических диа­ метров (0,39—2,0 мм), входная кромка которых была исполнена острой.

Весь исследуемый диапазон изменений чисел Рейнольдса можно разделить в зависимости от вида течения: ламинарный, переходный и турбулентный.

Для области ламинарного течения воздуха критериальная зави­

симость имеет вид

 

Nu = l,025Re°.2 .

(123)

В переходной области для различных форм каналов эта зависи­

мость

различная:

 

 

канал с острой входной кромкой

 

 

 

Nu =

0,00045Re1.2;

(124)

— канал

с закрученной входной кромкой

 

 

 

Nu =

0,00095Re1.1;

(125)

канал

цилиндрический

капиллярный

 

 

 

Nu = 0,00002Re4

(126)

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

Nu = 0,018Re0.8,

(127)

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

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

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

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

где AT— разность температур

между узлами;

q — тепловой поток.

 

В работе [77] произведено

определение действительного терми­

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

где h — расстояние, на котором измерялась разность температур AT; % — коэффициент теплопроводности материала лопаток.

Зная действительное и расчетное сопротивления, можно опре­ делить дополнительное термическое сопротивление как разность между ними

Дополнительное термическое сопротивление определяется вели­

чиной монтажных

зазоров в соединениях. Для сравнительной

 

 

 

Таблица 34

 

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

 

Т е р м и ч е с к о е с о п р о т и в л е н и е ,

м3-К/в'п

Тип хвостовика

 

 

 

 

д е й с т в и т е л ь н о е

д о п о л н и т е л ь н о е

р а с ч е т н о е

Т-образный

0,00322

0,00205

0,00117

Гребенчатый

0,00283

0,00043

0,0024

Елочный

0,00204

0,00068

0,00136

Лаваля

0,00150

0,00125

0,00025

оценки в табл. 34 приведены некоторые величины рассмотренных термических сопротивлений для четырех типов хвостовых соеди­ нений.

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

На термическое сопротивление оказывает влияние и окисляемость поверхностей замков. Проведенными исследованиями на авиацион­ ных турбинах типа PD-45 было установлено, что после работы в те­ чение 150—200 ч толщина окисной пленки на опорных поверхно­ стях хвостовиков составляет 2—3 мк, что соответствует условному увеличению толщины металла лопаток на 0,5 мм. Это не оказывает значительного влияния на увеличение термического сопротивления. Обычно термическое сопротивление рабочих лопаток оценивается методами аналогии.

Глава VI

МЕТОДЫ ОПРЕДЕЛЕНИЯ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В ОХЛАЖДАЕМЫХ УЗЛАХ ТУРБИНЫ

§ 44. Уравнение теплопроводности

Для нахождения температурного поля охла­ ждаемой детали необходимо решить уравнение теплопроводности. Под дифференциальным уравнением теплопроводности понимают математическую зависимость между физическими величинами, кото­ рые являются функциями пространства и времени. Дифференциаль­ ное уравнение теплопроводности описывает зависимость между температурой, временем и координатами элементарного объема. Такое уравнение характеризует изменение температуры в любой точке рассматриваемого тела и в любой момент времени. Обычно дифферен­ циальное уравнение твердого тела выводится упрощенным методом на основании теплового баланса элементарного параллелепипеда. В данной книге уравнение теплопроводности выведено исходя из закона сохранения и превращения энергии для тепловых процессов. Полученное таким образом уравнение теплопроводности будет при­ годно и при расчете температурного поля в любой вещественной среде [42].

Закон сохранения и превращения энергии для единицы объема

движущейся среды можно записать в

виде

 

 

 

 

qvdx + ALvdx

= pg (d3 n

+

,

(128)

где qv

количество тепла, втекающего за единицу времени в еди­

 

 

ницу объема;

 

 

 

 

A — тепловой эквивалент работы;

 

 

 

 

 

работа, совершаемая над единицей объема среды за еди­

 

 

ницу времени внешними силами;

 

 

T

время;

 

 

 

 

Pg — плотность среды;

 

 

 

 

э в

-

внутренняя энергия;

 

 

 

w — скорость движения

среды;

 

 

 

g —

ускорение свободного падения.

 

 

Как

известно, внутренняя

энергия среды

 

 

йЭ — di — Ad (рv)

пли

d3

= dl — A vdp — A pdv,

 

где v — удельный

объем;

p — давление;

 

 

/ — энтальпия.

 

Для совершенного

газа

 

 

dl — ср dt,

где ср—удельная

теплоемкость среды при постоянном давлении.

Тогда правая часть в скобках уравнения (128) может быть пред­ ставлена

 

 

d3B + М-^-

= Cpdt Avdp — Apdv -\-Ad~.

 

(129)

Выделим в теле некоторый объем V, который ограничен поверхно­

стью F. Уравнение

теплового баланса этого объема запишем:

 

 

 

\qvdV +

\qvdF

= \qVadV,

 

 

 

(130)

 

 

 

1'

 

F

 

 

v

 

 

 

 

 

где

q — вектор

теплового

потока,

направленный

по

нормали

 

 

к площадке в сторону, обратную направлению градиента

 

 

температур;

 

 

 

 

 

 

 

 

 

| qvdF

количество тепла,

перетекшее через

поверхность

F;

f

 

 

 

 

 

 

 

 

 

 

 

 

 

\ 4vBdV

количество тепла, выделенное

внутренними

источниками

і'

 

или поглощенное

объемом.

 

 

 

 

 

 

 

 

 

 

 

 

На

основании преобразования

Остроградского—Гаусса,

которое

справедливо, если в объеме нет разрывов функции,

выражение (130)

представим

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\qvdV+\diVqdV

 

= \

qVBdV.

 

 

 

(131)

 

 

 

V

v

 

 

v

 

 

 

 

 

Приняв во внимание, что

 

 

 

 

 

 

 

 

 

 

 

 

 

q =

—A, grad t,

 

 

 

 

 

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

выражение (131), получим'

 

 

 

 

 

 

 

 

<7„=div(a.grad0 + ?yB

 

 

 

(132)

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. /

дЧ .

d"t

,

дЧ

\ .

дХ

dt

,

 

 

дХ

dt

дХ

dt

.

/tQ9\

Подставив выражения (129) и (133) в (128), получим

і (JUL

_i_

J^L

J

JZ!L \

J

V dx2

^

б!/2

ґ

az2

J

r

_EL

_ 1

d* ' а* + ~dy'

dy +

, ЗЯ,

3/

,

і л I T

і

dp

,

сік \

(134)

Для общего решения теплообмена в вещественной среде к полу­ ченному уравнению необходимо присоединить еще три уравнения: движения, сплошности и состояния.

Поскольку полные производные от давления, энтальпии, квадрата скорости и удельного объема являются функциями координат и вре­ мени, то полный дифференциал их можно выразить:

 

 

 

dp

 

dp

 

JE..

dx

.

dp

dy

 

і

dp

 

dz

 

 

 

 

 

 

dx

 

dx

 

dx

 

~dT~^

 

dy

dx

 

 

dz

 

dx

'

 

 

 

 

 

dv

~

dv

 

dv

 

dx

 

.

dv

dy

 

,

dv

 

dz

'

 

 

 

 

 

dx

dx

 

dx

 

~dT

+

dy

dx

 

+

" a T -

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

_

d

t

+

dt

 

dx

 

,

dt

dy .

dt

 

dz

 

 

 

 

 

 

dx

~

dx

dx

 

dx

 

 

dy

dx

 

 

dz

 

dx

'

 

 

 

w2

 

 

d w2

 

 

 

a w2

 

dx

 

д w°-

 

dy

і

 

 

Ш2

dz

 

 

2g

 

 

~ *

 

 

1

 

 

 

 

l_

 

 

 

 

 

 

 

dx

 

 

 

 

dx

 

 

dx

 

dy

 

dx

 

 

 

dz

dx

 

 

 

 

dx

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, r

 

 

 

 

dx

 

du

\

 

dz

 

есть проекции вектора скорости

Учитывая, ч т о —

^

 

 

 

течения среды на соответствующие оси wx,

wy,

ш2,

выражение

(134)

запишем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- (

дЧ_

,

_ а ^ .

дЧ_\

, _dX_

dt_

 

dl_

_dt_,J&_

 

Jt_

 

_

\

dx2

+

dy2 +

az2

)

дх

'

дх

+

dy

'

ду

 

dz

'

dz +

q V e

~

 

 

 

 

 

(

 

dt

,

 

dt

 

.

аг

 

,

 

аг

\ .

 

 

 

 

 

 

/

и)2

 

 

 

 

ш3

 

 

^

Ш2

 

 

 

^ Jul2

 

 

 

+ т

Л

_ Ж - +

ч

, _ 1 1

+

^ _ і і _

+

.

4 .

 

-

 

 

 

 

 

 

 

dx

 

 

1

*

ax

1

 

" at/

 

'

'

 

az

 

 

 

 

 

 

— A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

dv

і

 

au

 

 

ay

 

 

 

аи

\ i

 

 

 

o r i

 

 

 

+ w ( ^ + » * - a r + B ' * i F + ^ - e r ) J -

 

 

( 1 3 5 )

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

12*

187

умеренных скоростях течения жидкости, уравнение (135) принимает вид:

V дх2 "т" ду" ~^ dz-

)

 

a.v * дх

 

5//

ді/ +

dz

' 3z +

~~

(

dt

,

dt

,

dt

,

dt

\

 

Е С Л И коэффициент

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

принять

постоянным,

то

 

 

 

 

 

.

/

дч .

з 2 /

.

а2 / \ .

 

 

 

 

 

 

 

/

dt

 

,

З/

,

 

3/ .

dt \

 

 

Для

твердого

тела,

которое

можно

принять

неподвижной

средой,

w = 0,

и тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л /

 

дЧ

,

дЧ .

дЧ

\

,

dt

 

 

И Л И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

T ^

+

9 , D = ^ .

 

(136)

При

отсутствии

источников тепла

уравнение (136) принимает

вид

 

 

 

 

 

 

 

- ^ V 2

i =

# - -

 

 

(137)

 

 

 

 

 

 

 

су

 

 

Зт

 

v

'

При

стационарном

процессе

 

=

0,

и тогда

 

 

 

 

 

 

 

 

 

 

V4 = О

 

 

 

 

 

 

 

 

 

 

 

су

 

 

 

 

 

 

или

 

 

 

 

 

 

 

аУ21 = О,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А.

 

 

 

 

 

 

 

 

 

 

 

 

где

а =

коэффициент

температуропроводности.

 

 

Коэффициент температуропроводности представляет собой отно­

шение коэффициента теплопроводности среды к ее объемной тепло­

емкости. Он пропорционален скорости распространения изотерми­

ческой поверхности. Величина, обратная коэффициенту температу­

ропроводности, характеризует инерционные свойства тела в отно­

шении распределений температур. Коэффициент температуропро­

водности зависит от температуры, а для

пористых тел —• от плот­

ности и влажности тела. Для материалов

охлаждаемых узлов физи­

ческие характеристики — коэффициенты теплопроводности и темпе­ ратуропроводности — нельзя считать постоянными, так как в неко­ торых случаях они значительно изменяются.

Так как коэффициент температуропроводности всегда положи­

тельная величина, то уравнение

(137) можно свести к виду

V а /

= 0.

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

Временные краевые условия сводятся к заданию скалярной функ­ ции, которая определяет распределение температуры в теле в неко­ торый момент времени

t = f (х; у; г).

Пространственные (граничные) краевые условия могут быть за­ даны четырьмя способами:

— температурой поверхности тела в любой момент времени (гра­ ничные условия первого рода);

•— тепловыми потоками на поверхности тела в любой момент времени (граничные условия второго рода);

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

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

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

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

§45. Решение дифференциального уравнения теплопроводности классическими методами

Выведенное в § 44 уравнение теплопровод­ ности (137) является однородным линейным дифференциальным уравнением второго порядка с частными производными. Такое диф­ ференциальное уравнение имеет бесчисленное множество частных решений. Классический метод решения состоит в том, что находятся