femreaddat.m

单元网格划分 一个最基本的有限元计算程序 二维传热问题

%FemReadDat

N=22;
str='dat.txt';
fid=fopen(str,'r');

M=fscanf(fid,'%f',[N N]);
fclose(fid);

nx=N-1;ny=N-1;
dx=1/nx;
dy=1/ny;
x=0:(1/nx):1;
y=0:(1/ny):1;
[X,Y]=meshgrid(x,y);
subplot(2,2,1)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('温度T-axis')
title('有限元法解二维热传导问题计算结果数据可视化 图 1')

subplot(2,2,2)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('温度T-axis')
title('图 2')
view(30,20)
subplot(2,2,3)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('温度T-axis')
title('图 3')
view(-30,45)
subplot(2,2,4)
surfl(X,Y,M)
view(-60,45)
shading interp
xlabel('X-axis'),ylabel('Y-axis'),zlabel('温度T-axis')

title('图 4')