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

книги из ГПНТБ / Бровкин Л.А. Температурные поля тел при нагреве и плавлении в промышленных печах учеб. пособие

.pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
11.54 Mб
Скачать

быотроло определения из решетя (11-II') средней по наосе.температуря пластины и температуры адиабатически изолированной плоскости по фо}~

цулам

 

 

V

f i J

- T

. -

С Ю Р Ю Ю )

. . ,25-Ш

 

 

 

 

 

 

n= e

 

 

 

 

 

т

(e,Fo)= % + T k b iJ L W * > ) .

. . . ^

 

 

 

 

 

n=e

 

 

n=e

 

 

Помимо номограмм для четных значений

в расчетах м<жно воспели-

еоваться также формулами:

 

 

 

 

 

 

' Ср -

J - txpf-efFo)

. m

с Z СГД

 

'no

 

п с Z г-

 

)

!

/

 

W

 

 

Z b p F o

 

0,4

S ~ h o

 

 

 

ОО

 

 

 

 

 

 

ОО

 

O - L - F

 

 

СЮ ^ r n - L ^ o \ j m t ^ i F o )

>F

-

'

a

4

 

Ч.ч

и

3

 

z

f

 

 

 

 

 

 

 

 

^ Lop

<op

Co

 

Fo *r 2

ехэ

/ .- ^ ^—

c p

 

OJCpl-^t ro

г.г

 

г

ч а

>

4

Z

 

3

15

t4

<5.

 

3,4

'гА

 

 

 

 

 

 

 

 

 

 

 

 

2.

2»_Номограммы и формулы рля учета начальных

 

 

 

_уолош й_нагрвЕа_

 

Значение FJoJjLfFl в форму

найдем ив (4-Н):

 

оо

3 X

формулах С15—0’) и (16-Ш)

ЭTx(f, Fb)_

1.Л г,

 

 

 

 

 

 

3 X

 

 

 

 

В частит случаях при 1 (* ,0 )-Т о и T f b F o ) ^

имеем;

Э Х =Ф

- T C t exP ( - C F o ) C 4 - i ) V j Fo ) - - ^

- 150

При начальном усл овии

 

1 2 4 -1 ) получим;

 

 

 

 

 

 

 

 

 

 

 

й7р д > М

 

Канальное условие

(25-1) дает для пластины

 

 

 

 

 

 

 

 

 

 

 

ъ

%

 

щ

. . . 30-Ш

где

 

ОО

 

 

 

 

 

 

 

 

Ср '

 

^ P ' ^ p f - s f F o ) ,

 

 

 

 

 

 

 

 

 

$Ункши

; Q 01

 

и ^

j для пластины показаны на рис. (10-Щ).

 

6._Прн:8£ рвсчета нагрева тела_диснретнш

 

 

 

уцовдетвореi-vov. границьк условий

 

 

 

Задача:

 

 

 

 

 

 

 

 

 

 

 

Найти температурное поле стальной пластины (

= 0,1

м; С =

=1000 вдж/меград;

Л =

10 вт/м град), нагреваемой от исходной постоян­

ной теьтературы

Т0 = ЗОСРа в среде,

температура

которой растет по

задаиноьу закону

 

Тг =[l000 + X ] °ii,

где Т

в секундах.

Параметры

теплообмена среда-тало:

с / =

100 вт/м^град;

6 = 4 10"® вт/м^градЧ

Время нагрева принять

Г£К -

1000 секунд.

 

 

 

 

Решение:

 

 

 

 

 

 

 

Fo ~

 

Расчеты будем вести в безразмерном времени

. Время

нагрет

запишется как

Fn =—15-1225—_ = 1

а закон изменения тен-

1

 

 

 

1ик

1000000

О, lz

 

.

 

..

пературы среды в виде

 

=

1000 ь 1000 Fo • Ограничим тоодость, рас­

чета температурного поля значением

= 3.

и

запишем искомое поле

в форме

(11-Е) с

тремя неизвестными коэффициентами /1( ,

 

/дя [йхоудеггия коэфрипиентов репим систему Цб-Р5)

ив трех уравнений.

- 151 -

- 152 -

Первое уравнение

( 0=0)

отвечает

времени

&-°1>.ое

Го =0рС~ t = l ) -

отвечает времени

Fo-

 

=0,5 и

третье

( О =2)-еначе-

нию Fc = FoK=1. Предварительно по рио.(

7 -Ш най­

дем

 

% (с),0 ,443,

 

0;

QOjol0; ^ /0 ,5 ) ^ 4 0 2 ) }

° Р ц }0,ъи 0 ,ш '}

^?5(0' 5)=

° '3;5*у 9 ] , Ш = Ч з05 1

9 ^ г (1)

= 0,465;

< ^ а(1)

.

0,590.

 

 

 

Ив форцулы (17-ffi) имеем:

 

 

 

 

 

 

 

 

_

Q

 

 

 

 

 

—1

^? =0,443’

 

аоос^-зоо4)

 

аооо-зоо) *1238

2-10

 

 

 

2‘ 10

 

 

 

 

 

 

 

 

 

 

 

Формула (19-ffi) 8апишетоя в виде;

 

 

 

Т

 

=

300 + 1238-1^ +

/1*1

+ Д -l"*

 

 

 

 

 

 

Z .

 

ь

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

4 - Х - 1 5 3 8 - / 1 3 .

 

 

 

 

 

Уравнение (20-Ш)

при подстановке числовых вначений дает

1238О0,305+( 7^-1538-

0,465-Ю,59оД3 «гЛСГ^гООО4- 7^) +

+ 0,5 ( 2000 - Тк )•

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

 

 

 

 

Л =

36300 - 7,72 Т -16* 10"10 . Т,4

 

 

 

 

'

 

 

 

к

 

 

' к

 

 

. и,

следовательно,

 

 

 

 

 

 

 

 

/}

=*

-37830

+ 8,72Х< + 16-10~20 J^4

 

 

Найдем предварительно:

 

 

 

 

 

а

t ^

 

j - т . + л , г ^

 

i

4 f f j - 4 Г & г -

*300+1238*0,5^+(-37830+8,727>16* 10"10 •7^)0,5 +

+(36300- 7,727^ -16.10*10 Тк4) О»6 *'

=-4940+1,63 7^ +2,45.10“ 10 J 4.

Расчетное уравнение относительно

на основе уравнешя

( Ь =1) запишется в виде:

- 153 - •

•1238 0,4 С 2 + (-3 7 8 3 0 + 8 ,72 7^ +16 10~10 Т ? ) 0 ,3 6 в +

+(36300-7,72 Тк -16,10-10 Т4) 0,315 = 1Ж 1ДА (

 

к

2 10

1

(1000 + 1000 0,5)4 -

(-4940 + 1,63 7^ + 2,45

10“10 Т*4)4 +

(ЮОО + 1000

0,5 + 4S40 -

1,637К -

2,45

Ю"10'/-4)

или

34,87 - 8,8(-1к.) - 1,195(-Ik.)4 _ pt,94 + 1,63—Zk +

1000 ЮОО L ' ЮОО

Уравнение (31-ПО решем подбором. Цусть Тк = 1900°К-. Тогда невязка в уравнении (31-Ш) составит +0,75. Бели принять

Тка 1900. то невязка составит -0,32. С точностью до градуса корень уравнения (31-Ш) найдется равным Тк = ШЗ°К. Искомые значения коэффициентов Д^и /Д окажутся равными соответственно

т37830 + 8,72 1893 + 16 10"10 1©34 = -125

Ьъ - 36300 * 7,72 1893 - 16 1СГ10 1©34 = 480.

Искомое температурное поле построено для Х= 1 и Х= 0 на рио.11-Ш.

£*-^5с291 нагрева терыически_тонкого тела

Граничное условие нагрева термически тонкого тела в любой момент времени Т : определяется через скорость нагрева:

с 1 Т )К ^ Ж - Ф г ^ Ы ( Т г- . ^ ) .

. .32-Ш

- 154 -

Рис.(IM l) Температурное поле пластины, нагреваемой излучением и конвекцией в среде о переменной температурой по рао^ чету методом дискретного удовлетворения

(З.з ).

- 153 -

Искомый аакон нагрева 7 ^ '? ) можно поучить прямым интегри-

ровашем (32-Ш) и метод дискретного удовлетворения в этом слу­ чае следует рассматривать как метод приближенного интегрирова­ ния. Подстановка в (32-Ш) выражения

7

= 7 +

/L

( n~j J

,.33-Ш

7

0

 

t

 

дает оиотеыу Cj нелинейных алгебраических уравнений относитель­ но коэффициентов Д)7.

10. _05ратная_аадача для тершчески_тонкого_тела

Представляет интерес

обратная задача,

когда

по ваданноцу

еакону

]~{Т)

можно определить етчения

~fj

и

 

и иопольэовать уравнеше

(32-Ш) для нахождения таких неиевест-

ных параметров "внешней"

теплоотдачи, как 7 * cL, <о j?4j.

Если

неизвестными явятся все

3 параметра, то, очевидно, уравнение

(32-Ш),

записанное для ГС.~

образует систему, решение ко-

 

 

 

J

<>

 

времени

с *

торой даст эффективные, усредненные для отреэка

параметры 7 , 6

, 7 .

 

 

 

 

 

11._Срраничения в применении метода дискретного удовлетворения

Вваключеше отметим, что расчет методом дискретного удов­

летворения может окаваться неустойчивым и привести к большим

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

не соблкщается прин-

- 156 -

цип равномерного на отрывке FoK выбора точек удовлетворения

 

 

и если

принимаемое число

Q

не обеспечено соответствующей точ­

ностью

определения функций

ср

 

 

 

Foh

 

 

Значение

определяет

отрезки

времени Fo ~Fo' ~

 

'

n-F

При достаточно малых

Fo

,

сое

 

^

J'*

*

г

>

F -

дние еначения Функций

9 Р ( Fb:J

« FVJ F o j . , ) мало

 

П- j

 

 

 

 

 

 

 

d

отучаются друг от друга и, если это

отли­

 

 

Hi,о

 

 

 

 

 

 

 

 

 

 

 

 

чие сопоставимо с погрешностью определения самих санкций 0Рт п

,

расчеты вести

бессмысленно.

3 этих условиях можно получить аб­

 

сурдные ревулътатн,

а увеличение

Cj не ведет к повышению точнос­

ти и надежности расчета.

 

 

 

 

 

 

 

 

По номограммам,

приведенным выла,

можно оценить функции

 

 

с точностью максимум третьего 8нака. Такая точность позволяет

 

принимать q

ив неравенства ;

 

 

 

 

 

 

 

Fo,

 

0,2

v 0,3 .

 

 

...34-Ш

 

 

1

 

 

 

 

 

 

 

 

 

 

Для расчетов

о малыми отрееками

Fo*

 

долж-

у - 1

функции

'*>"

 

 

 

 

 

 

 

■*— г'"

 

 

 

ны приниматься о точностью до четырех-шести знаков по таблицам,

составленным А.К.Соколовым и приведенным в Приложении.

§ б. ЩСЛШНЫВ МЕТОДЫ. МЕТОД КОНЕЧШХ РАЗНОСТЕЙ

деленные методы повволяют решать любые аедачи теплопровод­ ности и их иопольвование непрерывно раопиряетоя о внедрением

в инженерию практик^ вычислительных цифровых и аналоговых ишин. 1._Ожовные_понятия_в_мет^де конечных цаеноотей

Идея метода 8включается в еамене бесконечно малых приращений

- 157 -

Э Т ,Э г - ,Э Х , входящих в дифференциальное уравнение теплопро­ водности и краевые условия,на малые, но конечные приращения ДТ , А Т , А Х . Температурное поле тела при втом,естественно,

приходится рассматривать дискретно только в точках X ^ i -ЛХ

и только в моменты времени rC j - J ' A rt .

Метод конечных разностей удобно связывать с некоторой рас­

четной сеткой, имеющей шаги в пространстве А X и времени

Д'Г (рие.12-111). Декретное температурное поле представляется

температурами i J соответствующих узлов

сетки ( L , J

при-

V

7 7 X

&г).

 

В- расчетной сетке различают рады температур во времени год

номерами

0 ,1 ,2 ,... J - 1, J , j +1,.. .,$ -1, £/.

Нулевой ряд отража­

ет дискретно начальное температурное поле тела.

Расчет заключается в последовательном определении воех темпе­ ратур последующего ряда. ^ +1 во времени по известным темпера-

. турам в предцвущем ряду j , начиная с нулевого рада и кончая

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

Расчетными формулами являются записанные в коночных раёностях

для внутренних уадов сетки

( 6X0 и I V J ) само дифференциальное

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

а для узлов,отвечающих границам

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

- 158 -

 

 

 

-------

 

 

3

С

г г

*

пH Z F

 

 

г

 

 

п

 

 

 

 

 

‘J*2

 

 

 

 

!1

i

 

 

 

 

I i

 

 

 

 

i

t 1

 

 

 

 

 

!

 

 

 

 

 

1

 

 

 

 

1

 

 

 

I.

 

j

«

 

 

с .

 

 

 

г—

 

 

 

 

1

t-

 

 

 

'

~ А > Г

<

 

 

 

 

3 : i

й

 

Гн

0

 

 

 

Рио.(12-Ш) Расчетная сетка в одномерной задаче теплопроводности.

- 159 -

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