Прогрев тел под действием лучистого тепла

ИЗВЕСТИЯ

ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО ИНСТИТУТА имени С. М. КИРОВА

Том 89

1957

ПРОГРЕВ ТЕЛ ПОД ДЕЙСТВИЕМ ЛУЧИСТОГО ТЕПЛА

Г.П. БОЙКОВ

В работе рассматривается прогрев тел классической формы при граничном условии типа 2-го рода, когда теплообмен происходит по закону Стефана-Больцмана и лучистый поток симметричен.

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

Даются расчетные зависимости и соотношения, а также методика выбора расчетного интервала времени.

Известно, что при больших температурах источника тепла (приблизительно выше 800°С) передача тепла к нагреваемому телу происходит в основном лучистой энергией по закону Стефана-Больцмана.

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

Получить решение непосредственно из системы (1) - (4) затруднительно, так как пока еще не найдено пути, который позволил бы удовлетворительно согласовать решение дифференциального уравнения (1) с граничным условием вида (3). В связи с тем, что определение температурного поля, описанного условиями (1) - (4), имеет большой практический и принципиальный интерес, приходится прибегать к искусственным приемам, которые в некотором приближении дают возможность рассчитать прогрев тел под действием лучистого тепла.

Одним из наиболее распространенных приемов решения рассматриваемого вопроса является способ сведения граничных условий типа второго рода (3) к граничным условиям типа третьего рода путем введения условного коэффициента теплоотдачи излучением. Граничное условие, выраженное законом

\

дх

ЗЛ {Тс — Т,7^) ) ,

(30

позволяет получить решение системы (1), (2), (3'), (4) в виде:

которое в некотором приближении дает возможность вести практически расчеты.

Предлагаемая методика исследования позволяет иначе подойти к рассматриваемому вопросу и дает возможность несколько более приблизить расчетные данные к практическим. Согласно этой методике, температурное поле в пространстве рассматривается в виде непрерывной функциональной зависимости, дающей плавное изменение функции в зависимости от аргумента. Расчет же температурного поля во времени производится зонально с использованием принципа конечных разностей (см. [9], [10], [12]). Исходными данными при получении распределения температуры служат решения А. В. Лыкова для пластины, цилиндра и шара при постоянном лучистом потоке (см. [6]). Решение, например, для пластины им дано в виде:

Взяв выражение (6) за исходное, полагаем, что постоянный лучистый поток ga определяется по закону:

[7 Т, / Тп

ё С, = С о

I 100

100

и действует не на протяжении всего процесса прогрева, а лишь на протяжении малого конечного отрезка времени, равного т1 Подставив т1 л выражение (6) вместо т, получим распределение температуры по истечении первого отрезка времени от начала прогрева. Здесь же, заменив х на Я; получим значение температуры на поверхности пластины, Т1 (Я).

Имея эти данные, получаем распределение температуры по истечении второго такого же интервала времени т1 путем решения системы:

и подставкой в полученное решение:

вместо тзначения т1 (см. [1], [2], [6], [12]).

Продолжая аналогичные рассуждения, можно убедиться, что распределение температуры в пластине по истечении т-го момента времени от начала прогрева имеет вид:

или в сокращенной записи:

Последнее в критериях подобия1' запишется в форме

Для случая "тонких"2' тел (^<0,12) вести расчет по выражениям (14) и (15) будет нерационально. В этом случае удобнее пользоваться понятием средней температуры всей массы тела, которая, например, для пластины определится из соотношения:

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

Выражения [16] и [17] получены как среднее интегральное от распределений соответственно (14) и (15).

Как видно из зависимостей (14) и (16), при расчете распределения температуры требуется задаваться расчетным интервалом времени т1. Несомненно, что увеличение численного значения расчетного интервала времени т1 может значительно сократить объем вычислений. С другой стороны, увеличение расчетного интервала времени влечет за собой увеличение погрешности, превышение которой сверх допустимой величины снижает ценность расчета.

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

РЛ.ТС.Я

ёс • а

где Р - постоянное число, различное для пластины, цилиндра и шара

Второй шаг заключается в проверке полученного из формулы (18) расчетного интервала времени путем оценки приближения с использованием неравенства

где фт - действительная погрешность расчета, выраженная через лучистые потоки

Дт) - некоторая функция, обладающая следующими особенностями:

1)см. [3], [4], [11], [12].

2) О "тонких" телах см. [7], [8].

- во-первых, эта функция до некоторой степени отражает действительную погрешность расчета, рт приближаясь к ней по своей численной величине;

- во-вторых, действительная погрешность расчета рт всегда меньше или равна этой функции. Функция /(т1) находится из анализа изменения температур и лучистых потоков на поверхности тела.

Найдем эту функцию и докажем, что действительная погрешность рт< /(т1).

Пусть имеем некоторую ломаную кривую у1 = / (х) (см. рис. 1), которая обладает следующими свойствами:

1) Перелом кривой происходит через одинаковые интервалы аргумента, имеющие какое-то значение Дх1".

2) При уменьшении численной величины интервала Дх1" до Дх1' и далее в сторону нуля и, таким образом, при увеличении числа этих интервалов в сторону бесконечности, кривая у1 =/ (х) через ряд промежуточных ломаных кривых у2=/2(х); у3=/3(х) и т. д. стремится к некоторому пределу, который является "ис-тиным" и представляет уже плавную кривую зависимости у=/(х).

3) Каждая последующая ломаная кривая у2=/2(х); у3=/3(х) и т. д. является более близкой к "истинной" кривой у=/(х) по сравнению с предыдущей и имеет более плавный ход.

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

Учитывая все сказанное в пунктах 1 - 4, рассмотрим отдельно участок рис. 1, обозначаемый точками:

Из рис. 1а видно, что кривая 0-1-2 является более близкой к истинной кривой, чем кривая 0—1'—2', так как а1 стоит по своей величине значительно ближе к а, чем а'1, о чем говорит более плавный ход кривой 0-1-2.

Очевидно, что степень несовпадения кривой 0-1-2 и истинной кривой может быть оценена из соотношения

Чем меньшее значение имеет это отношение, тем кривая 0-1-2 ближе подходит к истинной и наоборот. Если ломаная кривая состоит из отдельных плавных кривых (см. рис. 16), функциональная зависимость которых известна в первых двух интервалах [у = ф(х); у = ф(х)], то выражение (V) может быть представлено через угловые коэффициенты:

Значит, для нашего случая при расчетах изменения температуры, например, на поверхности пластины, зависимость, аналогичная (ДО7), примет вид:

Здесь, согласно выражений (6) и (12),

поэтому

Сравнивая теперь

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

потому что

(так как Т0 - начальная температура тела, а Тис1 - температура на поверхности по истечении некоторого отрезка времени от начала прогрева).

Следует отметить, что ряд в выражениях (14) и (15) (третий член в фигурных скобках) быстро сходится, его значение с увеличением времени становится очень малым по сравнению с двумя первыми слагаемыми в фигурных скобках. Поэтому, начиная с некоторого момента времени, им можно пренебречь. Тогда дальнейший расчет, связанный с учетом лишь первых двух слагаемых в фигурных скобках, значительно упрощается.

Удельный расход тепла для пластины может быть определен по формуле:

¿=т

(21

Выражения, аналогичные (14) - (21), для цилиндра имеют вид:

Гт

ёс-К (9 ёс.

. 2.

> с I Я3

<I ___

ёс 4 ' ёс

I"

п — 1

1—тЛ

_2_ \

е

»с/ К, . ;. ёс

, е

(т-

I = т

в„г (/?) = во -Г Л',- [ 1Р(>1 V д,. - А <2Ст (1 - 2/?2)-

га 1 ¿-от—1

Ря • Л> Ол)

Л)

Ус. .е

- /и /^о,

(<3е, ~ <3с{1 + и)  6

—(яг—г) ¡).п ./"о,

(15'

ТУ* _ 2

Ге й-

1\

г • —"г(>) . ¿г

I 2 8 с  К

тп V- г.

Я2 ' ^ '

(16'

I = /тг

9т

2 ^ ёс, )

я = 1

,2

]

IX

£ = 2

Подобные соотношения могут быть написаны и для шара.

На рис. 2 показана кривая безразмерной температуры центра цилиндра (пунктирная линия), построенная по формуле (15').

Для сравнения сплошной линией показана безразмерная расчетная кривая, построенная по графикам Шака для тех же условий.

ЛИТЕРАТУРА

1. Дидкин В.А. и Кузнецов Н.И. Справочник по операционному исчислению, М.—Л., 1951.

2. Карслоу Х.С. и Эгер Д.Т. Операционные методы в прикладной математике, ГИТТЛ, М., 1948.

3. Иванцов Г.П. Анализ подобия нагрева металла в муфельной печи, Сборник "Промышленные печи", Металлургиздат, 1953.

4. Кирпичев М.В., Михеев М.А., Эйгенсон Л.С. Теплопередача, ГЭИ, М., 1940.

5. Лыков А.В. Теплопроводность нестационарных процессов, ГИТТЛ, М.-Л., 1948.

6. Лыков А.В. Теория теплопроводности, ГИТТЛ, М., 1952.

7. Немчинский А.Л. Тепловые расчеты термической обработки, Судпромгиз, 1953.

8. Иванцов Г.П. Нагрев металлов, Металлургиздат, М., 1950.

9. Тайц. Н.Ю. Технология нагрева стали, Металлургиздат, М., 1950.

10. Шваб В.А. Нестационарные температурные поля в твердых телах при изменяющихся граничных условиях, Вестник инженеров и техников, 3, 1935.

11. Эйгенсон Л.С. Моделирование, ГИ "Советская наука", М., 1952.

12. Бойков Г.П. Прогрев тел под действием лучистого тепла (диссертация), Томск, 1955.

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ

Т — переменная температура тела, °К. Тс — температура источника тепла, °К. Т0 — начальная температура тела, °К.

gc — лучистый поток от источника тепла в пустоту при степени черноты системы е„,

gcL — лучистый поток на поверхности тела в итый момент времени, направленный внутрь тела,

т — время, час.

х, r — текущие координаты, м. R — половина толщины (радиус) тела, м. X — коэффициент теплопроводности, шал0.

еп - приведенная степень черноты системы. С0 - коэффициент излучения черного тела,

ккал

2 U

м час K

аизл - коэффициент теплоотдачи излучением, ал .

м час

¡¡п - корни характеристических уравнений: sin ¡ = 0; J1 (¡) = 0. a - коэффициент температуропроводности, м1/час.

gc - R

K¡ = Х- Tc - критерии Кирпичева. axl

F01 = R2 - критерии Фурье.

01 R T

T - критерий безразмерной температуры.

©=T -

Оа = -критерий безразмерного лучистого потока. X = - критерии безразмерной координаты.