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

Методическое пособие 699

.pdf
Скачиваний:
7
Добавлен:
30.04.2022
Размер:
4.89 Mб
Скачать

 

 

kAe ak 2 sin k

 

 

 

Ae ak 2 cos k ,

(3.13)

 

 

 

 

 

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

 

 

 

 

ctg k

 

 

k

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

= Bi. Если обозначить k

=

 

 

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

можно записать следующим образом:

 

 

 

 

ctg

 

 

 

/ Bi .

(3.14)

 

y

y1 = ctg

y1

y1

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2 =

/Bi

0

 

2

 

3

 

4

1

2

3

4

Рис. 3.3 К решению уравнения (3.14)

Из анализа этого тригонометрического уравнения следует, что при каждом значении Bi существует бесконечное множество решений. Наиболее просто уравнение (3.14) можно решить графическим способом.

Обозначим левую часть уравнения (3.14) через y1 = ctg , а правую - через y1 = /Bi. Пересечение котангенсоиды y1 с прямой y2 дает нам значение корней характеристического

198

уравнения, т. е. (рис. 3.3).

Из рис. 3.3 следует, что мы имеем бесконечное множест-

во значении величины

n причем каждое последующее больше

предыдущего:

 

 

 

1

2

3 ...

n ...

Важно отметить, что каждому значению числа Bi отвечает своя совокупность корней уравнения (3.14).

198

Таблица 3.1

 

 

 

Значения n

для пластины

 

 

 

Bi

1

2

3

4

 

Bi

1

 

2

3

4

0

0,0000

3,1416

6,2832

9,4248

 

1,0

0,8603

 

3,4256

6,4373

9,5293

0,001

0,0316

3,1419

6,2833

9,4249

 

1,5

0,9882

 

3,5422

6,5097

9,5801

0,002

0,0447

3,1422

6,2835

9,4250

 

2,0

1,0769

 

3,6436

6,5783

9,6296

0,004

0,0632

3,1429

6,2838

9,4252

 

3,0

1,1925

 

3,8088

6,7040

9,7240

0,006

0,0774

3,1435

6,2841

9,4254

 

4,0

1,2646

 

3,9352

6,8140

9,8119

0,008

0,0893

3,1441

6,2845

9,4256

 

5,0

1,3138

 

4,0336

6,9096

9,8928

0,01

0,0998

3,1448

6,2848

9,4258

 

6,0

1,3496

 

4,1116

6,9924

9,9667

0,02

0,1410

3,1469

6,2864

9,4269

 

7,0

1,3866

 

4,1746

7,0640

10,0339

0,04

0,1987

3,1543

6,2895

9,4290

 

8,0

1,3978

 

4,2264

7,1263

10,0949

0,06

0,2425

3,1606

6,2927

9,4311

 

9,0

1,4149

 

4,2694

7,1806

10,1502

0,08

0,2791

3,1668

6,2959

9,4333

 

10,0

1,4289

 

4,3058

7,2281

10,2003

0,1

0,3111

3,1731

6,2991

9,4354

 

15,0

1,4729

 

4,4255

7,3959

10,3898

0,2

0,4328

3,2039

6,3148

9,4459

 

20,0

1,4961

 

4,4915

7,4954

10,5117

0,3

0,5218

3,2341

6,3305

9,4565

 

30,0

1,5202

 

4,5615

7,3058

10,6543

0,4

0,5932

3,2636

6,3461

9,4670

 

40,0

1,2325

 

4,5979

7,6647

10,6334

0,5

0,6533

3,2923

6,3616

9,4775

 

50,0

1,5400

 

4,6202

7,7012

10,7832

0,6

0,7051

3,3204

6,3770

9,4879

 

60,0

1,5451

 

4,6353

7,7259

10,8172

0,7

0,7561

3,3477

6,3923

9,4983

 

80,0

1,5514

 

4,6543

7,7573

10,8606

0,8

0,7910

3,3744

6,4074

9,5087

 

100

1,5552

 

4,6658

7,7764

10,8871

0,9

0,8274

3,4003

6,4224

9,5190

 

 

1,5708

 

4,7124

7,8540

10,9956

Первые четыре корня уравнения (3.14) 1, 2, 3, и 4 приведены в табл. 3.1 для различных значений числа Bi (от 0 до ).

При Bi прямая y2 = /Bi совпадает с осью абсцисс и корни уравнения будут равны:

 

 

;

 

3

;

 

5

… ,

 

2n 1

 

.

1

 

2

 

3

 

n

 

 

 

 

 

2

 

2

 

2

 

2

 

 

 

 

 

 

 

 

 

При Bi

0 прямая y2 =

/Bi совпадает с осью ординат и

тангенс угла наклона прямой стремится к бесконечности, при этом корни уравнения (3.14) равны:

1

0 ;

2

;

3

2 … ;

n

n 1 ,

 

 

 

 

где n = 1, 2, 3 ...

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

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

198

 

 

 

 

x

2 a

 

 

 

 

 

 

 

 

e 1

 

 

 

 

 

1

A1 cos

1

 

 

2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

2 a

 

 

 

 

 

 

 

2

 

 

 

 

 

A 2 cos

 

 

 

e

2 ;

 

 

 

 

 

 

2

2

 

 

 

 

 

 

(3.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. . . . . . . . .

 

 

 

 

 

 

 

 

 

 

 

x

 

2 a

 

 

 

 

 

 

 

 

n

 

 

 

 

 

A n cos

 

 

 

 

e

 

2 .

 

 

 

 

 

 

 

n

n

 

 

 

 

 

 

 

Полученные частные решения (3.15) будут удовлетворять дифференциальному уравнению при любых значениях постоянных A1, A2, ... , An, но ни одно из этих решений не будет соответствовать действительному распределению температуры в начальный момент времени. Однако путем наложения бесконечного числа таких распределений при соответствующем выборе величин An можно воспроизвести любую действительную температурную зависимость в начальный момент времени.

На основании сказанного общее решение можно представить суммой бесконечного ряда:

n

x

 

2

a

 

 

A n cos n

e

n

2 .

(3.16)

 

 

n 1

 

 

 

 

 

 

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

Постоянная An в уравнении (3.16) найдется из начальных условий. Подчинив уравнение (3.16) начальному условию, получим:

 

n

x

 

 

0 F x

A n cos n

.

(3.17)

 

 

n 1

 

 

 

Уравнение (3.17) есть разложение четной функции в ряд Фурье с заданными параметрами n, определяемыми характе-

198

ристическим уравнением (3.14). Для этой последовательности чисел n справедлива формула

cos

 

x

cos

 

x

dx

0

при n

m,

 

 

 

 

 

n

m

0

при n

m,

 

 

 

 

 

 

 

 

 

 

с помощью которой можно определить все коэффициенты An в уравнении (3.17). Для этого умножим обе части уравнения (3.17) на cos ( nx/ ) dx и затем проинтегрируем полученное соотношение по толщине пластины. Тогда

F x cos

 

x

dx A

 

cos2

 

x

dx ,

(3.18)

n

 

n

n

 

 

 

 

 

 

ибо все остальные слагаемые в правой части, для которых n m, обращаются в нуль. Интеграл в правой части соотношения (3.18) равен

 

 

1

1

 

 

sin 2 n .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

n

 

 

 

Тогда

 

 

 

 

 

 

 

 

An

 

n

 

 

 

F x cos

 

x

dx . (3.19)

 

sin n cos

 

 

 

n

 

n

n

 

 

 

 

 

 

Из уравнения (3.19) следует, что An является функцией только корня характеристического уравнения и начального распределения температуры.

Подставив полученное выражение для постоянной An в уравнение (3.16), получим окончательное выражение для температурного поля при охлаждении однородной пластины:

n

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

n

F x cos

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

sin n cos n

n

 

 

 

 

 

 

 

 

n 1

 

 

n

 

 

 

(3.20)

 

 

 

 

 

a

 

 

 

 

 

 

x

 

2

 

 

 

 

cos

 

e

n

2 .

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (3.20) позволяет получить значение температуры в любой точке пластины для любого момента времени

198

при любом начальном распределении температуры

0.

 

Если в начальный момент времени (

= 0) температура в

пластине распределена равномерно (рис. 3.2), т. е. t0 . tж =

0 =

const, то интеграл в уравнении (3.19) равен ( 02 /

n)sin

n. С

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

принимает

вид:

 

 

 

 

 

 

 

 

A n

 

 

2 sin

n

 

.

(3.21)

0

 

sin

n cos

 

 

 

n

n

 

 

Подставляя значение An, полученное для случая равномерного распределения температуры в пластине в начальный момент времени, в уравнение (3.20), получаем:

n

0 2sin n

 

 

 

x

 

2 a

 

 

 

 

 

cos

n

 

e

n

 

.

(3.22)

 

 

 

 

 

 

 

n sin n cos

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

n

 

 

 

 

 

 

 

Уравнению

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

(3.22)

целесообразно

придать безразмерную форму. Для этого разделим правую и

левую части уравнения (3.22) на

0. При этом обозначим:

 

 

 

 

 

 

 

 

 

 

 

 

 

A n

D n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После этих преобразований получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

2sin

 

 

 

 

 

x

 

2 a

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

cos

 

e

n

 

2

.

(3.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

n cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

n 1

n

n

 

 

 

 

 

 

 

 

 

Входящие в уравнение температурного поля (3.23) вели-

чины

 

 

; Dn;

n;

 

a

;

x

являются безразмерными и имеют

 

 

 

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следующий смысл:

0 = безразмерная температура; х/ = Х . безразмерная координата;

a / 2 = Fo - число Фурье, представляющее собой безразмерное время; Dn - безразмерный коэффициент.

С учетом последних обозначений уравнение (3.23) запишется:

129

n

2 sin

 

 

 

 

 

 

 

 

 

n

 

cos

 

X exp

2 Fo

. (3.24)

 

 

 

 

n

n 1 n sin

n cos

 

 

 

n

 

 

n

 

 

 

 

 

Анализ полученного решения. Так как

1, 2,

3, …,

n,

представляет собой ряд возрастающих чисел,

то чем больше

,

тем меньше роль последующего члена ряда по сравнению с предыдущим. Кроме того, чем больше число Fo, тем члены ряда будут убывать быстрее с увеличением номера n.

Многочисленные исследования показали, что уже при Fo 0,3 ряд (3.24) становится настолько быстросходящимся, что распределение температуры достаточно точно описывается

первым членом ряда:

 

2 sin

1

 

cos

 

X exp

2 Fo .

(3.25)

 

 

 

 

1

 

1 sin

1 cos

 

 

 

1

 

 

1

 

 

 

 

 

Ранее обозначено 2sin

1/(

1 + sin 1 cos 1) = D1. С уче-

том этого обозначения уравнение (3.25) можно записать в следующем виде:

D1 cos 1X exp

12 Fo .

(3.25 )

Величина D1 является только функцией числа Bi и заранее может быть рассчитана и табулирована. Кроме того, если рассматривать температуру для определенного значения Х = х/ , то и cos ( 1Х) является функцией Bi. Конкретно для оси

пластины Х = х/

= 0 и cos ( 1 0) = 1, а для поверхности

Х = х/ = 1 и cos (

1 1) = cos 1.

Для оси пластины произведение D1cos(0) обозначим как некоторую функцию N(Bi). Тогда уравнение (3.25) можно за-

писать в следующем виде:

 

X 0 N Bi exp 12 Fo .

(3.26)

Для поверхности 'пластины произведение D1cos

1 обо-

значим как некоторую функцию P(Bi) и уравнение (3.25 ) запишется так:

X 1 P Bi exp

12 Fo .

(3.27)

Функции N(Bi) и P(Bi) в уравнениях (3.26) и (3.27) табулированы и для расчета могут быть взяты из справочников.

198

Кроме того, из уравнений (3.26) и (3.27) следует, что при заданной координате безразмерная температура является только функцией двух безразмерных параметров Bi и Fo:

X 0 f1 Bi, Foи X 1 f2 Bi, Fo .

Логарифмируя уравнение (3.26), получаем:

ln

X 0

ln N Bi

2 Fo .

(3.28)

 

 

1

 

Аналогичное уравнение может быть получено после логарифмирования уравнения (3.27).

Из уравнения (3.28) следует, что при заданном значении координаты и при заданном Bi натуральный логарифм безразмерной температуры линейно зависит от времени. Последнее обстоятельство дает возможность представить для уравнений (3.26) и (3.27) графическое решение (рис. 3.4 и 3.5).

198

Рис. 3.4 Зависимость = f1 (Fo, Bi) для середины пластины

198

Рис. 3.5 Зависимость = f2 (Fo, Bi) для поверхности пластины

198