function oned
% 根据网上一维算例改了下边界条件，边界条件为绝热和恒定热流
% by hxg
clear all;clc;
%%%%%%%%%%%%%
%需要输入的物性参数
Lambda=10;%导热系数
cp=440;%热容
rou=7800;%密度
qw=500000;%壁面热流
%%%%%%%%%%%%%
a=Lambda/rou/cp;%定义中间系数
c=qw/Lambda;%定义中间系数
xspan=[0 0.012];%轴向坐标起止位置
tspan=[0 10];%仿真时间起止
ngrid=[1000 20];%空间网格数和时间网格数
%%%%%%%%%%%%%
%调用子函数
[T,x,t]=rechuandao(a,c,xspan,tspan,ngrid);
%画图
[x,t]=meshgrid(x,t);
figure(1)
mesh(x,t,T);
xlabel('x')
ylabel('t')
zlabel('T')
function [U,x,t]=rechuandao(a,c,xspan,tspan,ngrid)
% 热传导方程：
% Ut(x,t)=c^2*Uxx(x,t) a0.5
error('为了保证算法的收敛，请增大步长 h 或减小步长 k!')
end
%两个中间系数
s=1-2*r;
s2=c*h;
U=zeros(ngrid);%初步赋初值
% 初值条件，给了一个均温场
U(1,:)=20;
% 差分计算
for j=2:n
for i=2:m-1
U(j,i)=s*U(j-1,i)+r*(U(j-1,i-1)+U(j-1,i+1));
end
%%%%%%%下面是边界条件
U(j,1)=s2+U(j-1,2);%定热流边界条件
U(j,m)=U(j,m-1);%绝热边界条件
end
zxm
附注
j时间
zxm
附注
i为距离
zxm
附注
-k(dT/dx)=q
zxm
附注
r=(a x delta t)/(delta x)^2
本文档为【一维非稳态导热方程-热流绝热边界(附Matlab程序)】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483