实验6-平板稳态温度的计算.docVIP

实验6 平板稳态温度的计算 6.1实验指导 本实验在简单的平板热传导模型下，利用线性方程组计算平板稳态温度的分布。当选取节点较多时，MATLAB软件的应用就更显重要。最后利用MATLAB的绘图命令绘制平板的温度分布，形象地表示出计算结果。表6.1给出了与本实验相关的MATLAB命令或函数。 表6.1 与本实验相关的MATLAB命令或函数 命令 功能说明 位置 mod(i,N) 计算i除以N的余 例6.2 if mod(i,N)==0 判断i是否为N的整数倍 例6.2 mesh(T) T为一个矩阵，以T的行标为x轴，以T的列标为y轴，以T的元素值为z轴，绘制三维图形。 例6.2 表示钢板内部四个节点的温度。若忽略垂直于该截面方向的热交换，那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值，如。请计算该钢板的温度分布。 图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.0000 0 0 0 21.2500 0 1.0000 0 0 26.2500 0 0 1.0000 0 28.7500 0 0 0 1.0000 33.7500 得到方程组的解为：℃，℃，℃，℃。 例6.2 在例6.1中，把钢板内部分成了2×2个节点，本例把钢板内部分为5×5个节点，如图6.2所示。求钢板的稳态温度分布，并绘制温度分布图形。 图6.2 钢板的节点分布（25个节点） 解：根据例6.1中的讨论，知：5×5个节点就构成了一个具有25个方程的线性方程组。在MATLAB的M文件编辑器中编写la19.m文件： % 计算钢板的稳态温度分布 clear close all N=input(N=); % 输入节点数，共有N×N个节点 t_u=input(temperature up:); % 输入四个边界的温度值 t_d=input(temperature down:); t_l=input(temperature left:); t_r=input(temperature right:); A=zeros(N*N);b=zeros(N*N,1); % 构造N^2×N^2零矩阵A；构造N^2维零向量b for i=1:N*N % 矩阵A的主对角线元素都是4 A(i,i)=4; end for i=1:N*N; % 给矩阵A和向量b赋值 if i=N % 给向量b中和上边界节点对应的分量赋值 b(i)=t_u; end if mod(i,N)==0 % 给向量b中和右边界节点对应的分量赋值 b(i)=b(i)+t_r; end if mod(i,N)==1 % 给向量b中和左边界节点对应的分量赋值 b(i)=b(i)+t_l; end if iN*(N-1) % 给向量b中和下边界节点对应的分量赋值 b(i)=b(i)+t_d; end if iN % 给矩阵A中和上边界无关的节点所对应的元素赋值 A(i,i-N)=-1; end if mod(i,N)~=1 % 给矩阵A中和左边界无关的节点所对应的元素赋值 A(i,i-1)=-1; end if mod(i,N)~=0 % 给矩阵A中和右边界无关的节点所对应的元素赋值 A(i,i+1)=-1; end if i=N*(N-1) % 给矩阵A中和下边界无关的节点所对应的元素赋值 A(i,i+N)=-1; end end U=rref([A,b]); %