前言
在有限元求解完成后,需要对结果进行可视化,比如位移云图、应力云图和反力云图等。本文章使用MATLAB对二维四边形单元进行位移云图的可视化。
准备数据
首先要准备节点坐标(x,y),单元节点编号(逆时针编号)和节点位移(水平方向位移和竖直方向位移),下面提供一个实例数据。
% 坐标
XYZ = [ 0 0;
1 0;
2 0;
0 1;
1 1;
2 1;
0 2;
1 2;
2 2;];
% 单元
EL = [ 1 2 5 4;
2 3 6 5;
4 5 8 7;
5 6 9 8;];
% 水平方向位移
DX = [ 0;
0;
0;
0.0697349;
0.0604598;
0.0746214;
0.127212;
0.124422;
0.132233;];
% 竖直方向位移
DY = [ 0;
0;
0;
-0.00419673;
-0.0434221;
-0.0837584;
-0.0370525;
-0.0877267;
-0.138599;];
计算位移的模
% 位移的模
D = sqrt(DY.^2 + DX.^2);
与Abaqus一致的配色方案
% 颜色栏
Colormap = [0 0 1;
0 0.364705882000000 1;
0 0.725490196000000 1;
0 1 0.909803922000000;
0 1 0.545098039000000;
0 1 0.180392157000000;
0.180392157000000 1 0;
0.545098039000000 1 0;
0.909803922000000 1 0;
1 0.725490196000000 0;
1 0.364705882000000 0;
1 0 0;];
绘图
基本云图
利用patch函数即可完成基本位移云图的绘制,patch函数会根据节点位移自动进行双线性插值,其中XYZ为点坐标,EL为单元节点,D为节点位移。
patch('Faces',EL,'Vertices',XYZ,'FaceVertexCData',D,'FaceColor','interp','EdgeColor','k');
包含变形的云图
绘制变形后的位移云图只需要更新一下坐标即可。
XYZ = XYZ + [DX,DY]; % 更新坐标
patch('Faces',EL,'Vertices',XYZ,'FaceVertexCData',D,'FaceColor','interp','EdgeColor','k');
无框线(无网格)的云图
绘制无框线的云图只需要将EdgeColor设置为none即可。
XYZ = XYZ + [DX,DY]; % 更新坐标
patch('Faces',EL,'Vertices',XYZ,'FaceVertexCData',D,'FaceColor','interp','EdgeColor','none');
细节调整
云图铺满 & 关闭坐标轴
设置坐标轴紧凑格式,位移云图铺满整个绘图区。
axis tight % 坐标轴紧凑格式
取消坐标轴
axis off % 关闭坐标轴显示
云图左移
在设置颜色栏之前,首先要将云图向左移动,为颜色栏绘制腾出空间。
% 将图向左移动, 为绘制颜色栏做准备
ax = gca;
PositionNew = get(ax,'Position');
PositionNew(1) = 0.08; % 更改图形距离左边界的距离
set(gca,'Position',PositionNew);
绘制颜色栏
使用colormap函数即可绘制颜色栏,配色方案采用自定义配色方案、设置字体字号和设置标题等。
% 绘制颜色栏
Co = colorbar('eastoutside'); % 设置颜色栏位移图像外侧
colormap(Colormap) % 设置颜色栏
set(gca, 'FontName','Times New Roman', 'FontSize',14) % 设置字体字号
Co.Label.String = 'U, Magnitude'; % 设置颜色栏标题
设置颜色栏标签数
绘制了颜色栏后,可以看到标签并没有对齐,需要进行修改使之对齐。
% 设置颜色栏显示标签个数
NTicks = 13; % 标签个数
[cmin, cmax] = caxis;
Co.Ticks = linspace(cmin,cmax,NTicks)';
颜色栏标签格式
设置了标签数后,标签的格式还没有统一,模仿Abaqus标签采用科学计算法,设置方法如下。
% 设置颜色栏的刻度标签为科学计数法
yt = Co.Ticks; % 获取刻度位置
yt_labels = arrayfun(@(x) sprintf('%.3e', x), yt, 'UniformOutput', false); % 将刻度位置转换为科学计数法格式的字符串
Co.YTickLabel = yt_labels; % 设置刻度标签
保存云图
可以将位移云图保存成png图片格式或者svg矢量图格式。
% 保存图片
print(f,'-dpng','位移云图','-r600')%保存图像
print(f,'-dsvg','位移云图','-r600')%保存图像
MATLAB可视化结果
Abaqus可视化结果
可以看到,MATLAB可视化结果与Abaqus可视化结果一致。
总代码
可视化的总代码如下,后续还会更新MATLAB应力云图(高斯点外插、节点平均)、位移动画云图可视化的教程,制作不易,别忘了关注和点赞喔。
clc;clear;close all
set(0,'defaultfigurecolor','w');
%% 准备数据
% 坐标
XYZ = [ 0 0;
1 0;
2 0;
0 1;
1 1;
2 1;
0 2;
1 2;
2 2;];
% 单元
EL = [ 1 2 5 4;
2 3 6 5;
4 5 8 7;
5 6 9 8;];
% 水平方向位移
DX = [ 0;
0;
0;
0.0697349;
0.0604598;
0.0746214;
0.127212;
0.124422;
0.132233;];
% 竖直方向位移
DY = [ 0;
0;
0;
-0.00419673;
-0.0434221;
-0.0837584;
-0.0370525;
-0.0877267;
-0.138599;];
% 位移的模
D = sqrt(DY.^2 + DX.^2);
% 颜色栏
Colormap = [0 0 1;
0 0.364705882000000 1;
0 0.725490196000000 1;
0 1 0.909803922000000;
0 1 0.545098039000000;
0 1 0.180392157000000;
0.180392157000000 1 0;
0.545098039000000 1 0;
0.909803922000000 1 0;
1 0.725490196000000 0;
1 0.364705882000000 0;
1 0 0;];
%% 绘图
f = figure('Name','位移云图');
hold on
width = 600; %宽度,像素数
height = 400; %高度
left = 200; %距屏幕左下角水平距离
bottem = 200; %距屏幕左下角垂直距离
f.Position = [left,bottem,width,height];
XYZ = XYZ + [DX,DY]; % 更新坐标
patch('Faces',EL,'Vertices',XYZ,'FaceVertexCData',D,'FaceColor','interp','EdgeColor','k');
axis tight % 坐标轴紧凑格式
axis off % 关闭坐标轴显示
% 将图向左移动, 为绘制颜色栏做准备
ax = gca;
PositionNew = get(ax,'Position');
PositionNew(1) = 0.08; % 更改图形距离左边界的距离
set(gca,'Position',PositionNew);
% 绘制颜色栏
Co = colorbar('eastoutside'); % 设置颜色栏位移图像外侧
colormap(Colormap) % 设置颜色栏
set(gca, 'FontName','Times New Roman', 'FontSize',14) % 设置字体字号
Co.Label.String = 'U, Magnitude'; % 设置颜色栏标题
% 设置颜色栏显示标签个数
NTicks = 13; % 标签个数
[cmin, cmax] = caxis;
Co.Ticks = linspace(cmin,cmax,NTicks)';
% 设置颜色栏的刻度标签为科学计数法
yt = Co.Ticks; % 获取刻度位置
yt_labels = arrayfun(@(x) sprintf('%.3e', x), yt, 'UniformOutput', false); % 将刻度位置转换为科学计数法格式的字符串
Co.YTickLabel = yt_labels; % 设置刻度标签
% 保存图片
print(f,'-dpng','位移云图','-r600')%保存图像
print(f,'-dsvg','位移云图','-r600')%保存图像
总结
这只是一个基础的示例,实际中还会有更具体的、更细致的要求,这就需要再做额外调整;另外本人也仍在学习中,这只是个人的学习笔记,可能还有一些不足之处,欢迎指正。