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

UptoLike

ВУЗ: 

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

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

Рубрика: 

18
close all
clear all
% Решение уравнения теплопроводности
k=120; % коэффициент теплопроводности Si (Вт/м/К)
f=-1e5; % плотность источника тепла Вт/м^3
Lx=20e-2;
Ly=30e-2;
Lz=50e-2;
Tenv=300; % температура (К)
% Создание сетки
I=15; % число точек по x;
J=17; % число точек по y;
S=17; % число точек по t;
dx=Lx/(I-1);
dy=Ly/(J-1);
dz=Lz/(S-1);
X=0:dx:Lx;
Y=0:dy:Ly;
Z=0:dz:Lz;
%Заполнение коэффициентов матрицы А и В
N=I*J*S;
a=zeros(N,N);
b=zeros(1,N);
for s=2:S-1
for i=2:I-1
for j=2:J-1
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) =k*2/dx^2+k*2/dy^2+k*2/dz^2;
a(n,n+1)=-k/dx^2;
a(n,n-1)=-k/dx^2;
a(n,n+I)=-k/dy^2;
a(n,n-I)=-k/dy^2;
a(n,n-I*J)=-k/dz^2;
a(n,n+I*J)=-k/dz^2;
b(n)= -f;
end
end
end
% граничные условия
for i=1:I
for j=1:J
s=1;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
s=S;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
end
end
for i=2:I-1
for s=2:S-1

Страницы