end

T(i+1,nR+1)=T(i+1,nR);

T(i+1,nR+2)=T(i+1,nR);

end

Изменим порядок расположения столбцов обработав массив в двойном цикле For :

for i=1:nt

for j=1:nR

TT(i,j)=T(i,nR-j+1);

end

end

Построим поверхность описывающую полученную функциональную зависимость T(t,r):

figure(1)

mesh(TT)

Подпишем координатные оси

xlabel('R, MM')

ylabel('t, cek')

zlabel('T C')

6 АНАЛИЗ МОДЕЛИРОВАНИЯ И РАСЧЁТОВ

В результате численного решения дифференциального уравнения с помощью составленной программы получены данные, хорошо согласующиеся с аналитическим решением дифференциального уравнения приведенным во второй главе данной пояснительной записки.

Результаты получаемые с помощью данной программы можно использовать для моделирований реальных технологических процессов связанных с охлаждением и нагреванием цилиндрических каналов.

СПИСОК ИСПОЛЬЗОВАНОЙ ЛИТЕРАТУРЫ

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

2. Карслоу Г., Егер Д. Теплопроводность твердых тел. М., «Наука», 1964. 487 с.

3. Кирпичев М. В., Михеев М. А. Моделирование тепловых устройств. М.,изд-во АН СССР, 1936. 255 с.

4. Тябин Н. В. и др. В кн.: Теплообмен. 1974. Советские исследования. М., «Наука», 1975, с. 195—198.

5. Торнер «Технология переработки пластмасс­», Москва, Московский политехи, ин-т, 1965, № 1, с. 138—143.

ПРИЛОЖЕНИЕ1

clear, clc, clf, clg

t0=0;

tk=120;

nt=120;

Tc=30;

Tpol=170;

nR=10;

R=0.01;

dr=R/(nR-1);

dt=(tk-t0)/(nt-1);

a=0.00000056;

for i=1:nt+1

T(i,1)=Tc;

end

for j=1:nR+2

T(1,j)=Tpol*exp(-2000*(R-dr*(j-1))^2);

end

for i=1:nt

for j=1:nR

r=R-dr*(j-1)+0.0001*dr;

T(i+1,j+1)=T(i,j+1)+(a*dt/dr)*(((T(i,j+2)-2*T(i,j+1)+T(i,j))/dr)+((1/r)*(T(i,j+2)-T(i,j+1))));

end

T(i+1,nR+1)=T(i+1,nR);

T(i+1,nR+2)=T(i+1,nR);

end

for i=1:nt

for j=1:nR

TT(i,j)=T(i,nR-j+1);

end

end

figure(1)

mesh(TT)

xlabel('R, MM')

ylabel('t, cek')

zlabel('T C')

ПРИЛОЖЕНИЕ2