Розрахунок температурного поля необмеженої плоскої стінки

Стіна знаходиться під дією одностороннього охолодження (нагрівання) при незмінній температурі на другій поверхні стінки.

Розрахунок поля температур виконується для 6 моментів часу з 6 інтервалами по товщині стінки. Температурне поле задане диференціальним рівнянням теплопровідності:

(1.1).

де a — температуропроводність матеріалу, м2/c, Д2(T) — оператор Лапласа — сума других часткових похідних температури за координатами та початковими Т1k = 60 єК і граничними умовами:

Т0 = Т1k = 60 єК = const.

Т1N = Т1k S (I) 1.2.

Початкові умови описують поле параметрів процесу в деякий початковий момент часу, вибраний користувачем з тих чи інших міркувань. Це поле характеризує вихідний стан системи. Як правило, початкові умови являють собою функції розподілу в просторі показників, які входять до складу диференційного рівняння, тобто є функціями координат. Наприклад, при рішенні диференційного рівняння теплопровідності початкова умова може являти собою константу, рівняння або систему рівнянь, які дозволяють обчислити значення температури в усіх точках середовища в деякий початковий момент часу. У нашому випадку це константа:

Т0 (х, y, z)= const єК, (1.2).

Граничні умови показують, яким чином зовнішнє середовище впливає на показники процесу з часом. В загальному випадку вони являють собою функції часу і мають вигляд:

Т (х, y, z)= F (, x o, y o, z o), (1.3).

де x o, y o, z o — деякі значення координат на границях системи.

Наприклад, при рішенні диференційного рівняння теплопровідності для одномірного температурного поля (температура змінюється тільки вздовж однієї координати) гранична умова може мати вигляд:

T (z, x) = const єК, (1.4).

це свідчить, що на границі середовища (тіла) де координата Х=0, температура Т не змінюється з часом.

Граничні умови описують зміну температури на границях тіла з плином часу. Такою умовою може бути наведене вище рівняння (3), або деяке інше, яке зв’язує значення температури на поверхні і час:

T (z, x) = F (z, x, T1k) (1.5).

Сформульована таким чином задача являє собою математичну модель процесу.

Серед багатьох чисельних методів рішення диференційних рівнянь найбільше розповсюдження отримали кінцево-різничні. Вони засновані на апроксимації диференційних операторів відповідними різничними співвідношеннями, що призводить до створення системи алгебраїчних співвідношень, визначених в деяких точках розглядуваного тіла. Таким чином задача переростає із аналізу безперервного монотонного процесу в дискретну задачу, або, як кажуть, в задачу аналізу системи з параметрами згуртованими в вузлах. Такий підхід дозволяє ефективно застосовувати ітераційні методи обчислень.

Якщо задачу дещо спростити, вважаючи що температура змінюється тільки вздовж однієї з координат (випадок одномірної задачі), то рівняння (1) спроститься до вигляду:

(1.6).

Замінюючи в ньому диференціали різницями згідно з тлумаченням похідної одержимо:

(1.7).

(1.8).

Прийняті позначки в рівняннях:

і - індекс моментів часу;

k — індекс просторових вузлів;

h — крок часу, h = і+ - і-1.

— просторовий крок, = x k — x k-1.

Якщо співвідношення (7) та (8) підставити в рівняння (6) то одержимо його запис в кінцевих різницях:

= (1.9).

Вважаючи, що одна з величин h чи може бути задана користувачем в залежності від потрібних значень температури, або заданих моментів часу, та прийнявши, що:

одержимо рекурентну формулу для організації обчислень:

(1.11).

З рівняння (11) слідує, що температура в точці k тіла в кожний наступний i+1-й момент часу дорівнює її середньо арифметичному значенню в даний, i-й, момент часу в суміжних просторових вузлах — k-1 та k+1.

Рівняння вигляду (11) можна записати для i-1 моментів часу і для k-2 вузлів просторової сітки. В крайніх вузлах температура визначається за допомогою граничних умов, а первинний розподіл температури в товщі середовища визначаються початковою умовою.

Для рішення системи рівнянь (2), (4), (5) та (11) застосовують методи чисельного рішення систем лінійних рівнянь (Гауса, прогонки і т.д.).

Блок — схема алгоритму рішення систем рівнянь методом прогонки наведена нижче.

Рис 1.1. Блок-схема алгоритму для рішення рівняння теплопровідності методом прогонки

Текст програми для розрахунку температурного поля плоскої ізотропної стінки:

• 10 CLS
• 20 PRINT «************ ПРОГРАММА POLE-T ***********»
• 50 PRINT «Програма призначена для розрахунку температурного поля «
• 60 PRINT «плоскої стінки методом кінцевих різниць.»
• 70 INPUT «Вкажіть товщину стінки, М» ;DD
• 80 INPUT «Вкажіть число кроків просторової координати «;N
• 90 INPUT «Задайте коефіцієнт температуропровідності матеріалу м2/с» ;A
• 100 INPUT «Задайте число часових інтервалів «;M
• 110 DK=DD/N:DH=DK2/(2*A):DH=DH/3600:P=M+1:L=N+1
• 120 DIM T (P, L), T1(L), S (P), D (L)
• 130 REM’ПОЧАТКОВА УМОВА ЗАВДАЄТЬСЯ В РЯДКУ 490
• 140 REM' ГРАНИЧНІ УМОВИ ЗАВДАЮТЬСЯ В РЯДКУ 520
• 150 GOSUB 480
• 170 I=0
• 180 FOR K=0 TO N: X=K:D (K)=X*DK:T (I, K)=T1(K):NEXT K
• 190 FOR I=1 TO M: S (I)=I*DH
• 195 GOSUB 510
• 200 FOR K=0 TO N: IF K=0 THEN T (I, K)=T0:GOTO 240
• 210 IF K=N THEN T (I, K)=TN:GOTO 240
• 220 X=K:D (K)=X*DK
• 230 T (I, K)=(T (I-1,K-1)+T (I-1,K+1))/2
• 240 NEXT K
• 250 NEXT I
• 260 CLS: PRINT «РОЗРАХУНОК ТЕМПЕРАТУРНОГО ПОЛЯ ЗАКІНЧЕНО.»
• 270 PRINT «»
• 280 PRINT «РЕЗУЛЬТАТИ РОЗРАХУНКУ»
• 290 PRINT «————————————————————————————————— «
• 300 PRINT «Координати сітки ! Значення температур «
• 310 PRINT «——————————————————-! в вузлах сітки, «
• 320 PRINT «Час, годин ! Товщина, м ! T (i, k), K «
• 330 PRINT «——————————————————-!——————————————- «
• 340 FOR I=0 TO M
• 350 PRINT TAB (4);:PRINT USING «####.###» ;S (I);
• 360 FOR K=0 TO N: PRINT TAB (18);:PRINT USING «##.##» ;D (K);:PRINT TAB (33);:PRINT USING «####.###» ;T (I, K) 370 NEXT K
• 375 STOP
• 377 NEXT I
• 380 END
• 480 FOR K=0 TO N
• 490 T1(K)=273 + 60 NEXT K
• 500 RETURN
• 510 FOR K=0 TO N
• 520 T0= T1(K):TN= T1(K) S (I)1.2 NEXT K
• 530 RETURN

За програмою 1_POLE-T. BAS, яка призначена для розрахунку температурного поля одновимірних суцільних масивів методом кінцевих різниць ми отримуємо дані по розподілу температури по товщини стінки на підставі яких можемо побудувати графік. Вихідними даними для розрахунку були: товщина огорожі, число кроків просторових координат, коефіцієнт температуропроводності матеріалу огорожі, а також число інтервалів часу.

Рис. 1.2.