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

книги из ГПНТБ / Попов, В. М. Теплообмен через соединения на клеях

.pdf
Скачиваний:
47
Добавлен:
22.10.2023
Размер:
10.76 Mб
Скачать

основывается на рассмотрении потока тепла через одну контактную площадку и смежные твердые тела (см. рис. 1-7). Обычно полное термическое сопротивление контакта для такого элементарного канала в условиях вакуума выражается зависимостью вида (1-13).

Из соотношения (1-13) следует, что термическое со­ противление контакта обратно пропорционально тепло­ проводности материала у поверхности раздела. Таким образом, всякое снижение теплопроводности в области непосредственного контакта будет повышать контактное сопротивление для данной геометрии последнего. Следо­ вательно, образование на контактирующей площадке ма­ лотеплопроводной окисной пленки обусловливает повы­ шение контактного термического сопротивления.

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

При выборе тепловых моделей для двух сопряжен­ ных окисленных поверхностей наиболее целесообразным представляется в качестве упрощающего предположе­ ния рассматривать термические сопротивления отдельно в области окисленной пленки и основного металла для половины элементарного теплового канала (рис. 4-44,а), приравнивая при этом тепловые потоки и температуры во всех точках вдоль 0 < r < r 0, г= 60. Очевидно, что полное термическое сопротивление стягивания R cт в со­ ответствии с принятой моделью представляется как ре­ зультирующее для последовательно включенных сопро­ тивлений от стягивания линий теплового потока в обла­ сти основного металла Яст.м и окисной пленки R cт .о, т. е. R ct = R ct.м+^ст.о- Проанализируем составляющие правой

части этого выражения.

Сопротивление стягивания

13-745

193

в основном металле Rcт.м представим через сопротивле­ ние от стягивания линий теплового потока от сечения площадью яг20 при z » 60 до сечения площадью яг201 при z = 60, когда соблюдается условие постоянства теплового потока (рис. 4-44,в). Сопротивление стягивания в обла­ сти окисной пленки Rcт.0 выразим через эквивалентное сопротивление стягивания в цилиндре из материала, по сопротивлению идентичного основному металлу при

Рис. 4-44. К модельному представлению термического сопротив­ ления металлического контакта с окисной пленкой.

0-<2< бо (рис. 4-44,г). При этом стягивание линий те­ плового потока начинается от сечения площадью яг2<и

при

z = б0 и кончается в сечении с площадью

я а2 при

2 = 0

и постоянном тепловом потоке.

при этом

Рассмотрим сопротивление стягивания Рст.м,

математическая формулировка задачи сводится к реше­ нию уравнения Лапласа в цилиндрических координатах

(4-111) со следующими граничными условиями:

 

при z = 0 и 0 < У < г 01

Т = Т^ = const;

(4-130)

 

АТ

(4-131)

при 2 = 0 и r0l < r < r 0 -Я м1^ - = 0 ;

при г— оо для любого г

—ям1^ - д / ^ ; ]

(4-132)

194

flprt r =

r0

для любого z

—AM1 - f e — 0;

(4-133)

при

r = 0

 

 

- X Mid- L = 0.

(4-134)

Решение уравнения (4-111),

удовлетворяющее усло­

виям (4-132)

и

(4-133), представляется в виде

[Л. 122]

 

 

 

 

00

 

 

T =

A0+ - ? —2z +

Y

Ane-k*J0 (kr).

(4-135)

 

 

 

л^М1го

»

 

 

 

 

 

 

п=1

 

Здесь значения произвольной величины k в соответ­

ствии с условием

(4-134) являются корнями уравнения

 

 

 

J„(£ro) = 0,

(4-136)

где J п — функция Бесселя n-го порядка.

Смешанный характер условий (4-130), (4-131) накла­ дывает определенные трудности при определении вхо­

дящих в (4-135)

постоянных Ао и А п.

По этой причине

на области 2 = 0

и 0 < ;г < г О1 применим

граничное усло­

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

_______Q______

2nW oi ('щ — Г*)112

Следовательно, вместо условия (4-131) в первом при­ ближении для умеренных значений r0i/r0 имеем при z0= 0 и Го1< г < г 0:

дТ_

Q

(4-137)

А

2w0i Ош -г2)1/2

Ml дг

 

С помощью условий (4-136), (4-137), а также орто­ гональных свойств функций Бесселя и интегральной за­ висимости

гJo (kr)

dr —

sin {krоi)

(4-138)

I

 

k

 

 

 

 

13*

195

находим из (4-i35), что

 

 

 

 

 

 

 

 

 

 

^ 2- (" rJ0(^nO rfr +

A Jl„ j

r Jg (V ) & =

 

 

 

7t\,

о

 

 

 

 

 

0

 

 

 

 

 

 

10

 

 

rox

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

C

 

r j 0

( k n r)

,

 

 

 

 

 

-

2яХМ1г0,

J

 

(r^ — r2)1^2

 

 

после чего

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

sin ( У ,,)

 

(4-139)

 

 

2

О

 

 

 

2тгА.М1Гоj

kn

 

 

 

 

 

 

 

 

 

 

 

Подставляя значение Ап из (4-139) в (4-135), полу­

чаем:

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р __ д j_____ Q

g _i

Q

 

 

^

 

kn Z

Sin (fen^oi) Jo {&nr)

 

0

" V ' o

 

" V ^ oi

 

 

 

 

( V o) 2 Jo ( V o)

 

 

 

 

 

 

 

 

 

 

 

 

 

(4-140)

 

Очевидно,

что при z ^ c o ,

T —►^40-j—

 

при лю-

 

 

 

 

 

 

 

 

 

 

”V lr0

 

 

бом г.

 

 

 

 

 

 

 

 

 

 

 

 

 

Принятые

здесь граничные

условия

не

позволяют

определить Л0. В то же время из

(4-140)

и (4-136) вид­

но,

что А 0 представляет собой среднюю температуру при

z— >-оо и для малых значений r0i/r0 можем написать:

Г /j*

j

д

I

2Qr0

 

sin (^wr01) J; ( V

,)

 

J

Z”°

 

°

 

 

 

 

2j

( V

o ) 3 J o

( V

o)

 

0

 

 

 

 

 

 

 

n = 1

 

 

 

(4-141)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соотношение (4-141) можно представить в виде

 

 

ГМ=

ЛЙ

 

 

 

 

 

 

 

 

(4-142)

где

 

 

м

0 ^4АМ1г01 Ф»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф

 

 

V I s i n ( V o i ) Ji ( V o i )

 

/ 4 1 4 3 )

 

 

 

 

n = 1

( V o ) 3 J ;? (V o )

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

— коэффициент

стягивания

линий

теплового

потока

к мнимой площадке на границе раздела

металл — окис-

ная пленка (рис. 4-44,в ) .

 

 

 

 

 

 

 

 

 

196

Половина полного температурного перепада (рйС. 4-44,6) в зоне раздела мнимого контакта двух металли­ ческих тел находится согласно соотношениям (4-140) и (4-141) в виде

4АМНо 1 откуда полный температурный перепад представляется как

 

ь т = 2 ь т 1

Q

Ь,

(4-144)

 

2^мГо 1

где

приведенный

коэффициент

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

условно контактирующих тел с теплопроводностями материалов A,Mi и %к2-

Наконец,

полное контактное

У

 

 

 

 

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

единичного

 

 

 

 

 

теплового канала представля­ 0,8

 

 

 

 

ется в виде

 

 

 

 

 

 

 

 

 

=

£ =

2Кго,

(4-145)

 

 

 

 

 

 

у

 

 

¥

 

 

 

 

Измерение

коэффициента

 

 

 

 

стягивания

ipj

в

диапазоне

 

 

 

 

 

изменения отношения r0i/ro отО

¥

 

 

 

 

до 1 для вакуума приводится

 

 

 

 

 

на рис.

4-45.

 

 

 

 

0,2

0,9

0,6

0,9

1,0

Перейдем к анализу сопро­

,

roih

 

 

тивления стягивания в области

 

 

окисной пленки./?ст.о- Модель­

Рис. 4-45. Зависимость

ная схема зоны, имитирующей

коэффициента

стягивания

окисную

пленку

теплопровод­

в основном металле от от­

ностью Хо, предполагает фигу­

ношения /"oi/Го.

 

 

 

ру в форме

 

диска

диаметром

 

 

 

 

 

2% и толщиной 60 (рис. 4-44,г). Подвод тепла произво­ дится равномерно к торцу диска при г= Д 0 через пло­ щадь jtr2oi. На другом торце диска, т. е. при г —0, тепло­ вой поток максимально стягивается к контактному пят­ ну площадью я а2.

Решение задачи сводится к интегрированию уравне­

ния Лапласа в цилиндрических координатах

(4-111) при

следующих граничных условиях:

 

 

при z =

0 и 0 < г < а

 

}

при г =

"

) (4-146)

0 и f l < r < r 01 _ я 0— =

0;

 

197

П ри 2 = &о Для

л ю б о го Г —

Я0 ^ - = - 4 - 5

(4-147)

 

 

>1

 

при г = 0 —Яо1^ -= 0;

(4-148)

при г — г01

для любого г

—Яо^ - = 0 .

(4-149)

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

Из

условия

(4-146)

—Я0 дт_

 

~ у ;

г9 ° ’ а]

сле­

дует,

что

при

 

 

 

дг z=0

 

 

 

z — 0 начинается диск с радиусом основа­

ния а.

При г > а

и 2=

0

форма

диска нарушается. Это

видно из условия (4-146)

 

 

 

 

 

 

 

 

 

&Г_

 

 

г & о , а);

 

 

 

 

— Яо dz г = о

0;

г^(а , г01).

 

 

 

 

 

 

 

 

При 2=

80

задается

диск с радиусом

основания

го1:

 

 

 

 

-Я ^

Z=5_

Q h rt

 

 

 

 

 

 

Л°

Oz

 

 

 

 

Из всех условий задачи получается форма теплового канала в области окисной пленки, представленная на рис. 4-46, для которой из условий (4-148) и (4-149) сле-

Рис. 4-46. Модель теплового канала по толщине окисной пленки.

198

дует:

dr

 

ОТ

= 0.

г—Г01

1 дг

 

 

r= 0

Перейдем к решению уравнения (4-111). Допустим, что функция Т (г, z ) представляет собой произведение двух функций, зависящих только от г и только от г:

T(r,

z) = G(z)B(r).

(4-150)

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

новки их в уравнение (4-111)

получаем:

 

п d2B ,

G dB

I

d2G

B = 0.

 

G 4 F + -

dr

'

dz2

 

 

Поскольку задача

имеет

смысл

при 'условии

в В ф О,

то после деления последнего уравнения на GB, разде­ ления переменных и приравнивания полученного выра­ жения числу 'Р2 имеем:

1

d 2B

|

1

1

dB

 

 

В

dr2

'

г

В

dr

 

 

 

 

 

1

d 2G __

 

 

 

 

 

G

J z 2

 

 

 

После преобразований имеем:

 

 

 

d2B'

 

у

d- ? z W B = 0;

(4-151)

 

dr2

 

dr

 

 

 

 

d 2G

 

•p3G =

0.

(4-152)

 

dz 2 '

 

 

 

 

Решение уравнений

(4-151)

и (4-152)

произведем

с помощью метода степенных рядов. Решение уравнения (4-151) в общем виде представляется как

B ( r ) = /4 iJ 0 (|3г) +i42Yo(|3r),

где Jo(pr)— функция Бесселя 1-го рода нулевого поряд­ ка; Y0(pr) — функция Бесселя 2-го рода нулевого поряд­ ка, обладающая свойством при г>-0 величина Yo—>-°о.

Но при этом условии В ( г ) —>-оо и T = G B *оо. Од­ нако условие Тг^о— *~оо невозможно, так как по условию

задачи Т — конечная величина, т. е. Л2=

0 и

В (г) = А До(рг).

(4-153)

Можно показать, что собственные значения J0(^/") являются корнями уравнения J,(finrQl)=0, Из предположе199

ния Т(г, z) = G(z)-B(r) находим выражение -^-— G-^-,

которое после подстановки в (4-148) дает:

 

йВ

 

 

dr Г=Г100.

С другой стороны,

B = AiJ0(lr), откуда

dB

 

(r=rol)J,o(Р0 = о .

dr

r=r0i

 

Используем рекуррентное соотношение для Бесселе­ вых функций в виде J i(p r)= —Л'оОЗг), которое после подстановки в последнее уравнение дает:

J i( ( V o i) = 0 .

В общем случае функция В для п корней будет иметь вид:

В»(г)=/-,Ле(р*г). (4-154)

Решение уравнения '(4-152) представляется в виде

G =D ich|3z—D2sh |3z.

(4-155)

Подставляя решения (4-153) и (4-155) в выражение

(4-150), находим:

Т = А А 0{ $г ) [ D i ch ( pz) D z sh (|3z) ] .

(4-156)

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

ОО

 

 

 

т='2i AnJ»( М ch IP" <So — г)1

(4-157)

n=l

 

 

 

Рассмотрим Т = В ( г ) .

Из

уравнения (4-111)

имеем:

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

В =

Л01пг + А0.

 

С учетом условий (4-148) и (4-149) последнее урав­ нение представляется как

В = А 0.

(4-158)

Рассмотрим T2= G { z ). Из уравнения (4-111) полу­ чаем выражение d2G/dz2= 0, которое после интегрирова­

л о

ния приводится к виду

 

G = N o ( z - 6 o ) .

(4-159)

Из начального условия

(4-147) имеем:

дТ_

Q

или Nn =

дг г=Ьп

nVoi

 

"Voi

подстановка которого в (4-159) дает:

 

Л

QT

( K ~ z ) .

(4-160)

 

*Voi

 

Решение задачи с учетом (4-158),

(4-160) и (4-157)

представляется в виде

 

 

 

 

СО

 

Т = А0+ - Л - (8о - z) + У Ап ch фп (80 - г)] J 0ф„г).

A/oi

 

bU

 

 

 

П=1

(4-161)

 

 

 

Граничное условие (4-146) можно представить через разложение в ряд Фурье по функции Бесселя в виде

_Я Ш

= Л

2Q

 

 

J | ( М

Jo (РпН

 

 

(4-162)

o U ; * = .

 

 

 

Л=1 (РпОн) r oiJo (?nr oi)

 

 

Взаимным

преобразованием (4-161)

и (4-162)

получаем

выражение

 

 

 

 

 

 

 

 

 

 

д

___

____________ J» (Рда )_________

 

 

/ 4 j 03")

"

a

»a (PnCo,)2J g ( ^ o i ) s h ( M o ) ’

 

 

 

Подстановкой (4-163) преобразуем (4-161)

к

виду

 

Т = А0 +

кк„Г,

 

■<8^

2) + £ х

 

 

 

 

 

 

■о'01

 

 

 

 

 

 

 

У^Г\

ch^(S0г) J, ( М ) Jo

 

 

(4-164)

 

U Sh(Mo)

(?«'■..)*Jo (P*r*i)

 

 

 

 

 

 

 

п = 1

 

 

 

 

 

 

 

 

Выше отмечалось,

что

собственные значения

(4-161)

являются корнями уравнения

Ji(P«r0i) = 0,

что

с уче­

том (4-164)

позволяет представить А0 как среднюю тем­

пературу при z = б0.

 

 

 

 

 

 

 

 

Контактное сопротивление стягивания в области окйсной пленки в соответствии с модельной схемой можно представить в виде

R

ст.о

(4-165)

nko'oi

 

 

где Д Г = Г 1Т2 температурный перепад в слое окисла толщиной боСредние температуры соответственно вдоль площадок яа2 при z — 0 и яг2<ц при г — б0 находятся из выражений

T, = £ j r , =,rJr,

и

2% г 01

т =

1 2

лг01 О

Решая выражение (4-165) с учетом (4-164) и значе­ ний средних температур площадок, получаем конечное значение контактного термического сопротивления в об­ ласти окисной пленки:

^ ст-°

пАоa 'l'2’

(4-166)

 

где

J? м

(4-167)

th ( Ш (pnr01)3Jq(P-n^oi)

n=l

— коэффициент стягивания линий теплового потока в окисной пленке к площадке контакта.

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

£c t . o = 2£'cT.0=-=U<!>2,

пК0а

2\ \

где Я0= -т—sj-Y------приведенный коэффициент Ло1~т~Ло2

водности окисных пленок контактирующих тел.

(4-168)

теплопро-

Анализ уравнения (4-167) показывает, что коэффи­ циент ф2 находится в функциональной зависимости от отношений а/г01 и 60/а- Графическая интерпретация вы­ ражения для ф2 в функции а/г01 для различных значений

202

Соседние файлы в папке книги из ГПНТБ