Методы математической физики. Куликова И.В - 7 стр.

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

Составители: 

Рубрика: 

7
0В данной задаче нет источников тепла, поэтому
. =
i
f
Ниже приведен один из вариантов кода программы, позволяющей
решить СЛАУ (3) и вывести результаты решения, которые приведены на
рисунке 1.
close all
clear all
% Решение уравнения теплопроводности
k=120; % удельная теплоемкость Si (Вт/м/К)
f=0; % плотность мощности источника тепла Вт/м^3
L=20e-2; % толщина стенки (м)
Tenv=300; % температура (К)при х=0
T0=310; %температура (К) при х=L
% Создание геометрии
N=20; % число точек по толщине;
dx=L/(I-1);
X=0:dx:L;
% Решение СЛАУ
a=zeros(I,I);
b=zeros(1,I);
for i=2:I-1
a(i,i) =-k*2/dx^2;
a(i,i+1)=k/dx^2;
a(i,i-1)=k/dx^2;
b(i) =f;
end
% граничные условия
i=1;
a(i,i) = 1;
b(i) = Tenv;
i=I;
a(i,i) = 1;
b(i) = T0;
T=b/a';
plot(X,T)
grid on
xlabel('m')
ylabel('K')
На рисунках 2-4 приведено распределение температуры при наличии
внутренних источников тепла и граничных условиях, приведенных выше, для
различного числа точек I = 30 – рисунок 2, I = 5 – рисунок 3 и I = 3 – рисунок
4.

Страницы