传热学数值计算.docxVIP

传热学数值计算作业 数值解程序：tw1=40 %三边温度 tw2=100 %一边温度正弦变化幅度l1=40 %板长L1:40厘米l2=20 %板宽L2:20厘米m=41 %分划成40*20的网格n=21 k=2dx=l1/(m-1)c=ones(n,m)for i=1:m a2(i)=tw1+tw2*sin(pi*dx*(i-1)/l1)c(1,i)=tw1 ,c(n,i)=a2(i)endfor j=1:nc(j,1)=tw1c(j,m)=tw1endwhile (abs(c(j,i)-k)0.0001)k=c(j,i)for i=2:m-1for j=2:n-1c(j,i)=0.25*(c(j,i-1)+c(j,i+1)+c(j-1,i)+c(j+1,i))endendend数值解中各网格点的温度值：数值二维温度分布图像：解析解程序：tw1=40tw2=100l1=40l2=20p=40q=20x(1)=0for i=1:px(i+1)=x(i)+1endy(1)=0for j=1:qy(j+1)=y(j)+1endfor i=1:p+1for j=1:q+1n(j,i)=tw1+tw2*sinh(pi*y(j)/l1)*sin(pi*x(i)/l1)/sinh(