实验6平板稳态温度的计算.docx
- 文档编号:12574721
- 上传时间:2023-04-20
- 格式:DOCX
- 页数:9
- 大小:92.28KB
实验6平板稳态温度的计算.docx
《实验6平板稳态温度的计算.docx》由会员分享,可在线阅读,更多相关《实验6平板稳态温度的计算.docx(9页珍藏版)》请在冰豆网上搜索。
实验6平板稳态温度的计算
实验6-平板稳态温度的计算
实验6平板稳态温度的计算
6.1实验指导
本实验在简单的平板热传导模型下,利用线性方程组计算平板稳态温度的分布。
当选取节点较多时,MATLAB软件的应用就更显重要。
最后利用MATLAB的绘图命令绘制平板的温度分布,形象地表示出计算结果。
表6.1给出了与本实验相关的MATLAB命令或函数。
表6.1与本实验相关的MATLAB命令或函数
命令
功能说明
位置
mod(i,N)
计算i除以N的余
例6.2
ifmod(i,N)==0
判断i是否为N的整数倍
例6.2
mesh(T)
T为一个矩阵,以T的行标为x轴值,以T的列标为y轴值,以T的元素值为z轴值,绘制三维图形。
例6.2
6.2实验内容
例6.1在钢板热传导的研究中,常常用节点温度来描述钢板温度的分布。
假设图6.1中钢板已经达到稳态温度分布,上下、左右四个边界的温度值如图所示,而
表示钢板内部四个节点的温度。
若忽略垂直于该截面方向的热交换,那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值,如
。
请计算该钢板的温度分布。
图6.1钢板的节点分布(4个节点)
解:
根据已知条件可以得到以下线性方程组:
化简为标准的矩阵形式如下:
在MATLAB命令窗口输入:
A=[4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4];
b=[30;50;60;80];
U=rref([A,b])
结果为:
U=
1.000000021.2500
01.00000026.2500
001.0000028.7500
0001.000033.7500
得到方程组的解为:
℃,
℃,
℃,
℃。
例6.2在例6.1中,把钢板内部分成了2×2个节点,本例把钢板内部分为5×5个节点,如图6.2所示。
求钢板的稳态温度分布,并绘制温度分布图形。
图6.2钢板的节点分布(25个节点)
解:
根据例6.1中的讨论,知:
5×5个节点就构成了一个具有25个方程的线性方程组。
在MATLAB的M文件编辑器中编写la19.m文件:
%计算钢板的稳态温度分布
clear
closeall
N=input('N=');%输入节点数,共有N×N个节点
t_u=input('temperatureup:
');%输入四个边界的温度值
t_d=input('temperaturedown:
');
t_l=input('temperatureleft:
');
t_r=input('temperatureright:
');
A=zeros(N*N);b=zeros(N*N,1);%构造N^2×N^2零矩阵A;构造N^2维零向量b
fori=1:
N*N%矩阵A的主对角线元素都是4
A(i,i)=4;
end
fori=1:
N*N;%给矩阵A和向量b赋值
ifi<=N%给向量b中和上边界节点对应的分量赋值
b(i)=t_u;
end
ifmod(i,N)==0%给向量b中和右边界节点对应的分量赋值
b(i)=b(i)+t_r;
end
ifmod(i,N)==1%给向量b中和左边界节点对应的分量赋值
b(i)=b(i)+t_l;
end
ifi>N*(N-1)%给向量b中和下边界节点对应的分量赋值
b(i)=b(i)+t_d;
end
ifi>N%给矩阵A中和上边界无关的节点所对应的元素赋值
A(i,i-N)=-1;
end
ifmod(i,N)~=1%给矩阵A中和左边界无关的节点所对应的元素赋值
A(i,i-1)=-1;
end
ifmod(i,N)~=0%给矩阵A中和右边界无关的节点所对应的元素赋值
A(i,i+1)=-1;
end
ifi<=N*(N-1)%给矩阵A中和下边界无关的节点所对应的元素赋值
A(i,i+N)=-1;
end
end
U=rref([A,b]);%对增广矩阵进行行初等变换,化为最简行阶梯矩阵
fori=1:
N%把矩阵U的最后一列按矩阵形式赋给B
forj=1:
N
B(i,j)=U(N*(i-1)+j,N*N+1);
end
end
T(2:
N+1,2:
N+1)=B;%把钢板内部温度和和四周温度值放入矩阵T中
T(1,2:
N+1)=t_u;
T(N+2,2:
N+1)=t_d;
T(2:
N+1,1)=t_l;
T(2:
N+1,N+2)=t_r;
T([1,N+2],[1,N+2])=NaN;%矩阵四个角没有温度值,故把非数NaN放入
T%显示计算结果
mesh(T)%对钢板温度值绘制成曲面图形
在MATLAB的命令窗口中输入:
la19
人机对话及运算结果为:
N=5
temperatureup:
20
temperaturedown:
50
temperatureleft:
10
temperatureright:
30
T=
NaN20.000020.000020.000020.000020.0000NaN
10.000016.565219.866721.840023.341525.312530.0000
10.000016.395821.060624.153826.212127.911130.0000
10.000017.958323.826127.500029.442630.119430.0000
10.000021.608728.787932.576933.939433.123330.0000
10.000029.687537.139540.083340.614038.434830.0000
NaN50.000050.000050.000050.000050.0000NaN
钢板的温度分布如图6.3所示。
其中x、y坐标分别表示钢板横、纵方向的节点数,高度表示节点的温度值,该三维图形形象地反映了钢板的温度分布。
图6.3钢板的温度分布
6.3实验习题
1、如图6.4所示,假设钢板已经达到稳态温度分布。
分别用
来表示钢板的上、下、左、右四个边界的温度值。
请编写一个通用Matlab程序,求出钢板内部六个节点的温度值,要求四个边界的温度值由用户输入。
(1)
,单位℃;
(2)
,单位℃;
(3)
单位℃。
图6.4钢板的节点分布(6个节点)
解:
先列写对此6点的方程:
写成矩阵形式
程序为
A=[4,-1,0,-1,0,0;-1,4,-1,0,-1,0;0,-1,4,0,0,-1;-1,0,0,4,-1,0;0,-1,0,-1,4,-1;0,0,-1,0,-1,4],T_u=input('T_u='),T_d=input('T_d='),T_l=input('T_u='),T_r=input('T_r='),
B=[T_l+T_u;T_u;T_r+T_u;T_l+T_d;T_d;T_d+T_r],
X=A\B
(1)X=
19.4203
20.4348
22.7536
27.2464
29.5652
30.5797
(2)X=
38.8406
40.8696
45.5072
54.4928
59.1304
61.1594
其余按比例增加。
(此题用了近一小时)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 平板 稳态 温度 计算