
книги из ГПНТБ / Нестеров Ю.Ф. Теория и расчет судовой тепловой изоляции
.pdfРис. 15. Изоляционная конструкция (а) и ее прямая (б) и обрат ная (в) электрические модели
Коэффициенты теплопроводности материалов можно |
принимать |
по средней рабочей температуре, так как зависимость к |
^ / (/) яв |
ляется незначительной. |
|
Будем предполагать, что отдельные детали конструкции с различ ными коэффициентами і плотно прилегают друг к другу и благодаря хорошему контакту соприкасающиеся точки различных поверхностей имеют одну и ту же температуру.
Таким образом, при решении задач теплопроводности, по суще
ству, не будет учитываться лишь |
влияние |
крепежных деталей, мон |
|||
тажных дефектов, диффузии водяного пара |
|
и увлажнения |
изоляции, |
||
а также продувания ее воздухом. |
|
|
|
|
|
Явление теплопроводности подчиняется основному закону тепло |
|||||
проводности (закону |
Фурье): |
|
|
|
|
^ |
= - ^ ^ |
= - ^ |
Б |
dlT, |
(36) |
где dQ — элементарное количество тепла, проходящее через пло щадку dFT (в м2) на изотермической поверхности в направлении
нормали |
пт |
(dFT ~ BdlT); |
— = grad |
t — температурный |
гра |
|
диент, |
°С/м |
(пт — нормаль |
к изотермической |
поверхности, |
м)\ |
|
dlT — элементарная длина изотермической |
линии |
в плоскости |
хтуг |
|||
(рис. 15, а). |
|
|
|
|
|
Применение закона Фурье к случаю двухмерной задачи при ста ционарном режиме и независимости коэффициентов К от температуры
приводит к дифференциальному |
уравнению |
Лапласа |
- J i + |
^ r = 0, |
(37) |
дх; |
ду; |
|
где i=f |
(хт, ут)— гармоническая температурная функция координат. |
|
За |
начало координат принят левый нижний |
угол конструкции |
(см. рис. 35, а), ось хт направлена вправо, ось ут |
— вверх. |
|
Уравнение Лапласа описывает распределение температур в темпе |
ратурном поле. Оно представляет собой дифференциальное уравне ние в частных производных второго порядка эллиптического типа.
Для того чтобы из бесчисленного количества решений системы уравнений (36) и (37) выделить одно-единственное решение, описы вающее отдельный конкретный процесс теплопроводности, необхо димо к этой системе присоединить условия однозначности.
Полная совокупность условий однозначности для стационарных процессов теплопроводности распадается на:
1) геометрические условия, задающие конкретную форму и раз меры изоляционной конструкции;
2)физические условия, которые задают физические параметры материалов (Яи, л.д, л.с), существенные для процесса;
3)граничные условия, характеризующие особенности распределе ния температур на границах конструкции.
Вычислять коэффициент теплопередачи для всей изолированной поверхности судна, содержащей периодически повторяющиеся эле-
менты набора, неудобно. С целью упрощения расчетов всегда можно ограничиваться вычислением k лишь для периодически повторяю щегося участка с единичным профилем набора. Обычно шаг такого участка равен шпации s. При этом стенку набора высотой h распола
гают |
посередине |
шпации. |
|
В симметричных изоляционных конструкциях линии АТВТ |
и |
||
CT DT |
(рис. 15, а), |
выделяющие периодически повторяющийся участок |
и располагаемые посередине расстояния между стенками профилей набора, одновременно являются и граничными линиями тока тепла.
Или иначе, по боковым поверхностям |
АТВТ |
и CT DT |
участок |
адиа- |
||
батно изолирован и тепловой |
поток |
через |
боковые |
границы |
равен |
|
нулю. При этом на левой и правой боковых |
границах |
|
|
|||
|
|
= ( # - ) |
= 0 - |
|
О») |
|
\ |
ОХт / л е в |
\ ОХт / прав |
|
|
|
|
Для симметричных |
изоляционных |
конструкций |
= 0 |
не |
только на боковых границах, но и на оси симметрии. Поэтому для них достаточно рассматривать половину конструкции.
На боковых границах несимметричных |
|
конструкций |
|
( * \ |
= ( J L \ |
. |
( 3 9 ) |
\ ОХт / л е в |
\ ОХт / п р а в |
|
v ' |
Однако такие граничные условия реализовать сложно, так как пришлось бы моделировать стенку с большим количеством элементов набора. Поэтому для несимметричных конструкций строгие гранич ные условия (39) заменяют условиями (38). Такая замена приводит к незначительной положительной погрешности. Чем больше шпация s, тем меньше погрешность, возникающая вследствие замены гранич ных условий, так как при увеличении шпации практически осущест вляемые границы (прямые, перпендикулярные к наружным поверх ностям) приближаются к линиям тока у краев участка.
Вычисления, |
сделанные |
для |
участков, содержащих различное |
||
количество |
повторяющихся |
элементов набора |
[2], показали, что |
||
при s/h > |
2 погрешность не превосходит 1 %. Поэтому влиянием крае |
||||
вого эффекта у |
несимметричных |
конструкций |
практически можно |
пренебречь даже при малых шпациях. Обычно же s/h > (2,5—3,0); при этом погрешность, вызываемая заменой граничных условий, не превосходит 0,1%.
Из тепловых сеток (рис. 52, а—54, а; 58, а и др.) видно, что область влияния стального профиля на температурное поле ограничена рас стоянием s = (2-Г-4) h. При s > (2н-4) h линии тока у свободных краев модели становятся прямыми (неискаженными). Это означает, что при указанном условии соседние профили набора вообще не влияют на температурное поле исследуемой части конструкции. Сле довательно, вместо изучения влияния периодически повторяющегося набора можно исследовать влияние его единичного профиля.
Таким образом, во всех случаях граничные условия однознач ности задаются в следующем виде:
на |
наружной поверхности Л Т С Х |
|
|
t = |
ta |
= |
const; |
|
||
на |
внутренней поверхности |
BTDT |
|
|
t = |
tB |
= |
const; |
(40) |
|
|
|
|
|
|
|
|
|
|
|
|
на левой и правой боковых границах АТВТ |
и |
C T D T |
|
|
|
|
||||
При этом пусть |
для определенности |
ta |
>• |
|
|
|
|
|||
Следовательно, |
на наружной и |
внутренней поверхностях кон |
||||||||
струкции заданы граничные |
условия |
первого рода, а на ее боковых |
||||||||
нетеплопроводных границах — второго |
рода. |
|
|
|
|
Для определения количества тепла, проходящего через поверх ность Fr конечных размеров, нужно проинтегрировать уравне ние (36):
Длина конструкции В (вдоль оси zT) вынесена за знак интеграла, так как В = const для всех частей изоляционной конструкции.
Таким образом, чтобы найти тепловой поток Q, необходимо знать
подынтегральную |
функцию |
а следовательно, и |
распределение |
температур t — f |
(хт, ут). |
Знание поля температур |
является необ |
ходимым условием для решения всех без исключения задач тепло проводности. Для этого предварительно следует найти решение диф ференциального уравнения Лапласа (37).
С |
математической точки зрения решение уравнения (37) для судо |
||||
вой изоляции обычно сводится к плоской внутренней задаче |
Дирихле, |
||||
т. е. к отысканию гармонической температурной функции |
координат |
||||
t — f |
(хт, уг), |
удовлетворяющей |
уравнению Лапласа при заданных |
||
граничных |
условиях. |
|
|
||
Следовательно, чтобы определить Q и k, необходимо решить урав |
|||||
нения |
(37) |
и |
(41) совместно. |
|
|
|
|
|
§ |
20 |
|
|
|
|
Теория метода |
|
|
|
|
|
электротепловой аналогии |
|
Выявление аналогии. Выявим существование аналогии между тепловыми и электрическими явлениями.
Явление электропроводности управляется основным законом элек тропроводности (законом Ома):
dl |
1 |
dv dF, = — |
б |
ду |
Ш9; |
(42) |
|
р |
дпэ |
9 |
дп. |
|
|
здесь dl — сила Тока, о; р — удельное электрическое |
сопротивление |
||
материала, |
ом-м (1/р — удельная электропроводность, |
1/ом-м); |
|
J^- = grad |
v — градиент электрического потенциала |
*, |
в/м; v — |
электрический потенциал, в; пэ — нормаль к изопотенциальной по верхности, м\ dF3 — элементарная площадка, расположенная на изопотенциальной поверхности, м2; б — толщина электрической модели, м; dl3 — элементарная длина изопотенциальной линии в пло скости хэуэ (рис. 15, б), м.
Дифференциальное уравнение электропроводности, описывающее распределение потенциалов, в случае стационарного плоского поля при постоянных физических параметрах представляет собой уравне ние Лапласа:
Т ^ + ^ |
= 0 ' |
<4 3 > |
где хэ иэ у — координаты электрического поля; v = f (хэ, |
уэ) — гар |
моническая потенциальная функция координат, удовлетворяющая уравнению (43).
Математические выражения физических законов (36) и (42), кото рым подчиняются явления тепло- и электропроводности, имеют совер шенно одинаковую структуру. Они отличаются только тем, что вхо дящие в них физические переменные и параметры имеют различную размерность. Такое сходство имеет своим следствием то, что и описы вающие эти явления дифференциальные уравнения (37) и (43) (в основе которых лежат упомянутые законы) также имеют одинаковое строе ние. Это обстоятельство и позволяет провести формальную матема тическую аналогию между сравниваемыми явлениями. Сила тока / является аналогом теплового потока Q, электрический потенциал v — аналогом температуры t, а электропроводность 1/р — аналогом коэф фициента теплопроводности X.
Правила осуществления аналогии. Рассмотрим содержание тре бований, которые необходимо и достаточно выполнить для осущест вления аналогии.
Метод ЭТА базируется на теории подобия, так как аналогию физи ческих явлений можно рассматривать как наиболее общий случай подобия. Правила достижения аналогии принципиально не отли чаются от общих правил осуществления подобия (в узком смысле этого слова), которые даются третьей (основной) теоремой подобия. Согласно теореме, для того чтобы модель стала подобной конструкции, необходимо и достаточно соблюсти подобие условий однозначности (геометрических, физических и граничных).
Чтобы оказались подобными граничные условия однозначности, должны быть одинаковыми законы распределения потенциалов и температур на границах электрической модели и тепловой конструк ции.
* Индексом «э» будем отмечать величины, относящиеся к электрическому полю.
75
Подобие граничных условий осуществляют путем поддержания постоянных потенциалов v„ и vB вдоль границ модели АС и BD (рис. 15, б), соответствующих наружной и внутренней поверхностям конструкции, для чего эти границы должны соприкасаться во всех точках с шинами, на которых v = const. Свободные края модели АВ и CD, соответствующие боковым граничным линиям тока в конструк ции, должны быть электрически изолированы, вследствие чего на
свободных краях будет соблюдаться необходимое |
условие |
г-= 0. |
|||||
Таким образом, граничные условия для тонкой плоской электри |
|||||||
ческой модели записываются в следующем виде: |
|
|
|
||||
на |
наружной границе AC |
v = vH= |
const; |
|
|||
на |
внутренней |
границе |
BD |
v = vH |
= |
const; |
^ |
|
|
|
|
dv |
|
|
|
на |
свободных |
боковых |
краях модели АВ и CD |
— = |
0. |
|
|
|
|
|
|
охэ |
|
|
|
Математическое описание явлений формально одинаковыми диф ференциальными уравнениями и условиями однозначности является необходимой предпосылкой теории подобия.
Если аналогия достигнута, то числовые значения всех тепловых величин (температуры, теплового потока и др.) находятся в постоян ных отношениях с аналогичными электрическими величинами (по тенциалом, силой тока и пр.).
Введение масштабных преобразований. Аналогию явлений про анализируем методом теории подобия. В первой теореме теории подо бия утверждается, что нет необходимости изучать непосредственную -связь между отдельными размерными величинами, существенными для процесса. Значительно проще найти связь между безразмерными инвариантами подобия (симплексами и комплексами), составленными из этих размерных величин. Уравнения, выражающие связь между физическими величинами, представленные в безразмерном виде, ока зываются одинаковыми для всех подобных явлений. Поэтому для установления количественной связи между аналогичными физиче скими величинами приведем математические описания к безразмерной форме.
Для того чтобы привести к безразмерному виду геометрические условия однозначности, необходимо ввести вместо абсолютных отно сительные размеры конструкции и модели. В качестве определяющего размера тепловой конструкции возьмем высоту набора hT (рис. 15, а). Тогда относительные координаты температурного поля окажутся равными:
Х ^ Ь |
У* = Ъ - |
( 4 5 ) |
В качестве масштаба для линейных размеров электрической модели выберем сходственный размер h3 (рис. 15, б):
Хэ = ^ ; Ул = £ - . |
(46) |
Большими буквами будем обозначать значения всех величин, вы раженных в относительном масштабе.
Если геометрическое подобие достигнуто, то относительные раз меры и координаты сходственных точек конструкции и модели будут численно равными:
ХТ = ХЭ = Х = idem; Г т = Y3 = Y = idem. |
(47) |
При выборе размеров модели нет никаких ограничений. Приведем к безразмерному виду физические условия однознач
ности для температурного и электрического полей. Физическим пара метром, характеризующим свойства теплопроводной среды, является коэффициент к. В качестве масштаба для различных значений к возь
мем коэффициент теплопроводности |
изоляционного |
материала |
ки. |
|
Тогда получим следующие |
значения |
относительных |
коэффициентов |
|
теплопроводности: для изоляционного материала — Л и — kjkn |
— 1; |
|||
для дерева — Л д = кд/ки; |
для стали — Л с = кс/ки. |
|
|
Если бы все области плоской модели имели одинаковую толщину б, то физические свойства материала определялись бы только его электропроводностью 1/р. Модель из бумаги с областями различной электропроводности получают склеиванием внакладку специальным электропроводным клеем разных (по сопротивлению) сортов бумаги с неодинаковыми толщинами б или нескольких слоев одной и той же бумаги. При этом различные области плоской модели обладают раз
ной толщиной б в направлении, параллельном оси z3, тогда как |
все |
|||
детали изоляционной конструкции имеют одинаковую длину |
В |
|||
вдоль оси zT. В этом |
заключается |
существенное отличие |
модели |
от |
конструкции. По этой |
причине в |
качестве физического |
параметра, |
характеризующего свойства электропроводной бумаги, необходимо
брать |
не величину 1/р, |
а отношение б/р, которое назовем электриче |
|||
ской |
проводимостью бумаги g6 = б/р \1ом. Тогда электрическое со |
||||
противление бумаги R6 |
= |
р/б ом. При подборе бумаги удобнее поль |
|||
зоваться этой величиной, |
а не проводимостью |
g6. |
|||
В качестве масштаба для величины g6 |
возьмем электрическую про |
||||
водимость бумаги gfi. иг имитирующей |
область |
изоляционного мате |
|||
риала с коэффициентом |
ка. |
Тогда получим следующие значения отно |
сительных электрических проводимостей листов бумаги, воспроизво дящих различные области изоляционной конструкции: для изо
ляционного |
материала — G6. и |
= |
g6.Jg6, |
и = Re.JRe.* |
= |
^ |
Д л я |
дерева — Сб .д = £ б . д/#б .и = R6, JR6. |
д ; |
для |
стали — G6c= |
g6 |
Jg6.u |
= |
|
= #б. J Кб. с- |
|
|
|
|
|
|
|
Для того чтобы форма и расположение изолиний и линий тока в температурном и электрическом полях оказались подобными, необ ходимо подобрать листы бумаги таким образом, чтобы относительные
физические параметры-аналоги обоих полей |
получились равными: |
Л и = G6. и = idem; Л д == <Зб - д = idem; Л с |
= G6 .c = idem. (48) |
Таким образом, физическое подобие достигается подбором прово димостей или количества листов бумаги, имитирующих те или иные детали конструкции, по соотношениям (48).
Приведем к безразмерному виду граничные условия однознач ности для конструкции и модели.
Все температуры поля будем отсчитывать от наименьшей темпера туры на внутренней поверхности tB как от нуля. В качестве масштаба для избыточной температуры t — tB возьмем разность между наи большей и наименьшей температурами tH — tB наружной и внутрен ней поверхностей конструкции. Тогда для относительных разностей температур (представленных в долях общей разности температур)
получим |
выражение |
|
|
|
Т = |
. |
(49) |
|
|
' н < в |
|
Такая |
структура инварианта |
подобия |
(49) объясняется тем, что |
для явлений теплопроводности существенны только разности темпе ратур, а истинные температуры выпадают из рассмотрения, так как процессы распространения тепла происходят всецело под воздейст вием разностей температур.
Для электрических потенциалов примем за начало отсчета соот ветствующий аналог Ув , а в качестве масштаба для избыточных потен циалов v — vB выберем разность vH — vB. Тогда для относительных разностей электрических потенциалов (измеренных в долях от общей
разности потенциалов) получим |
выражение |
|
у = |
v ~ v * . |
(50) |
После приведения к безразмерному виду граничные условия одно значности (40) и (44) принимают следующий вид:
для тепловой конструкции
на поверхности |
АТСТ |
при t — tH |
= const |
Т = TU= \; |
|
на |
поверхности |
BTDT |
при t = tB |
= const |
Т = Тв = 0; |
|
|
|
|
|
дТ |
на |
поверхностях |
АТВТ |
и CTD.t |
|
= 0; |
для электрической модели
на |
границе |
АС |
при |
v = vH = const |
на |
границе |
BD |
при v = vB = const |
|
на |
свободных краях |
АВ и CD |
V = VW~ 1; V —- VB = 0;
-4^- = 0.
дХ
(51)
(52)
Следовательно, граничные условия оказались тождественно равными:
Тп |
= VH = |
idem; Тв |
= VB = idem; - J ^ - = - j ^ - = і dem. |
|
Из-за |
отсутствия ограничений масштабы для разностей |
t — tB |
||
и v — vB |
можно |
выбирать |
произвольно, необходимо только, |
чтобы |
величины |
ta — tB |
и УН — vB |
были аналогами, т. е. представляли со |
бой значения соответствующих разностей в сходственных точках. Такая автомодельность сравниваемых явлений позволяет произ вольно выбирать значения потенциалов в электрической модели, при
этом решение задачи теплопроводности будет пригодно для любых температурных напоров.
Введем масштабные преобразования (49) и (45) в дифференциаль ное уравнение (37):
Так как масштаб |
t„ — tB |
отличен от нуля, |
|
а размер |
hT |
конечен, то |
||||||
уравнение |
Лапласа (37) принимает |
безразмерный вид |
|
|||||||||
|
|
|
|
д 2 Т |
• |
д 2 Т |
0. |
|
|
|
|
(53) |
|
|
|
|
дХ2т |
' |
дҐ'т |
|
|
|
|
|
|
После |
подстановки |
масштабных |
преобразований |
(50) и (46) диф |
||||||||
ференциальное уравнение |
электропроводности |
(43) |
в |
безразмерной |
||||||||
форме принимает вид |
|
|
|
|
|
|
|
|
|
|
||
|
дЪ |
, |
дЧ |
_ |
v« — vB |
( &V |
, |
d*V \ _ Q |
|
|||
|
<э*2 |
^ |
^ |
|
A* |
^ ах2э |
1 |
ак2 |
э |
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дХ2э |
|
д¥2э |
|
|
|
|
|
Введение масштабных преобразований в дифференциальные урав нения и условия однозначности не дает никаких уравнений связи
между |
масштабами. Поэтому |
все сходственные |
масштабы могут быть |
|||||
выбраны произвольно. |
|
|
|
|
|
|
|
|
Если условиться рассматривать безразмерные переменные Т и V |
||||||||
только |
в сходственных точках, то при помощи |
равенства (47) урав |
||||||
нения |
(53) и (54) можно переписать |
в следующем виде: |
||||||
|
# |
Г |
+ |
^ |
= |
0 |
(55) |
|
|
<ЭХ2 |
1 |
дУ |
|
|
|
||
и |
|
|
|
|
|
|
|
|
|
^ - |
|
+ |
- |
^ |
= |
0. |
(56) |
После |
приведения к безразмерному |
виду дифференциальные урав |
нения (55) и (56) получаются тождественно равными и отличаются
только обозначениями |
безразмерных |
аналогов: |
|
д*Т . |
д*Т |
d2V |
, d2V |
+ |
-жтг = |
|
+ -svr = idem. |
Для того чтобы решения этих уравнений были тождественными, необходимо и достаточно, чтобы безразмерные условия однозначности были численно равными.
Если модель удовлетворяет всем условиям подобия, то относи тельные условия однозначности получаются одинаковыми. При этом решения уравнений (55) и (56) тождественно совпадают как по форме
функциональных зависимостей, так и по численным значениям пере менных:
Т = V |
=--f (X, Y) = idem. |
(57) |
Измерив распределение |
относительных разностей |
потенциалов V |
в модели и представив результаты измерений в виде графической совокупности изопотенциальных линий, можно рассматривать най
денное поле V как поле относительных разностей температур Т, |
а изо- |
|||||||||
потенциальные линии |
как |
изотермы. |
В любых сходственных |
точ |
||||||
ках конструкции и модели, имеющих |
относительные |
координаты |
X |
|||||||
и Y, безразмерные разности температур и потенциалов численно |
||||||||||
одинаковы: |
|
|
|
|
у |
|
|
|
|
|
гр |
t |
t& |
V |
Ув |
|
|
|
|
|
|
|
~~ tu — tB ~ Ун — Va |
~~ |
|
|
|
|
|
|||
С помощью этих равенств можно вычислять иоле |
температур |
t, |
||||||||
зная поле V (а также |
температуры |
tH |
и |
tB): t = |
tB |
+ |
V (tH — |
tB). |
||
Температурное поле, определяемое |
уравнением |
(57), |
является |
ре |
шением внутренней задачи Дирихле для уравнения Лапласа. Обычно температурное поле представляют в графической форме.
Таким образом, формальная аналогия между дифференциальными уравнениями, описывающими процессы различной физической при роды, приводит к формально одинаковым их решениям. Метод ЭТА позволяет изучать процесс теплопроводности, происходящий в кон струкции, на плоской модели, в которой протекает электрический ток, и заменять измерение тепловых величин электрическими.
Вычисление теплового потока. Рассмотрим электрический спо соб вычисления теплового потока, основанный на аналогии между
основными |
законами теплопроводности (36) и |
электропровод |
|||||
ности (42). |
|
|
|
|
|
|
|
Введем масштабные преобразования (49) и (45) в выражение за |
|||||||
кона Фурье |
(41) |
под знак |
интеграла: |
|
|||
Q = — X B \ |
{ І н h T |
. d |
L r |
= |
№(*„-tB) |
j ^ - d L T . |
|
Безразмерный |
интеграл со знаком минус |
|
|||||
|
|
|
|
|
Ч |
|
|
назовем критерием формы температурного поля. |
|
||||||
Следовательно, количество тепла, |
проходящее |
через изоляцион |
|||||
ную конструкцию в единицу |
времени, |
|
|||||
|
|
Q=KB |
(tH ~ |
ts) |
Ф т ккаліч. |
(58) |
Для определения общей силы тока, протекающего в модели, не
обходимо проинтегрировать |
уравнение |
закона |
Ома (42): |
||
I = |
і |
-=— aFb — |
—5— dL. |
||
р |
дп3 |
* |
р .1 |
дпэ |
|
|
|
|
|
|
э |
F3 |
/э |