How to use contourf to plot a mesh not generated in matlab, i.e. imported mesh coordinates from abaqus. The issue is that coordinates are not sequential.

Question

Subham Prasad 2024-7-7，19:20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2135241-how-to-use-contourf-to-plot-a-mesh-not-generated-in-matlab-i-e-imported-mesh-coordinates-from-abaq

评论： Subham Prasad 2024-7-9，5:50

I'm trying to use contourf with x and y coordinates being random or non sequential. I have z value for each X y coordinate. Normaly we would create a rectangular mesh and then fill up the z values. Here the issue is that the geometry is meshed in abaqus and the mesh nodes and element related information is fed to my code. Now the nodes are not sequential matrix.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2024-7-7，21:09

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2135241-how-to-use-contourf-to-plot-a-mesh-not-generated-in-matlab-i-e-imported-mesh-coordinates-from-abaq#answer_1482346

It would help tto have your data.

Since there are (x,y,z) coordinates (there can be duplicates), I would use the scatteredInterpolant function to create an interpolant, and then use either ndgrid or meshgrid tto create the relevant argument matrices using the ‘x’and ‘y’.vectors to define them.

There are several examples on MATLAB Answers using that approach, that create column vectors for ‘x’, ‘y’ and ‘z’ to use with scatteredInterpolant,.and then generate the ‘z’ matrix from the argument matrices. Then use surf or surfc to plot them.

.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2024-7-8，12:46

在 MATLAB Online 中打开

This sort of works, however it takes 605.3 seconds to complete and about as much time to render the surface plots. It fails, because it fills (interpolates) the voids areas. I do not see any way to correct for that.

imshow(imread('boundary conditions.jpg'))
LD = load('plotdata.mat');
plot = LD.plot;
x = plot(:,1);
y = plot(:,2);
z = plot(:,3);
figure
scatter3(x, y, z, 5, z, 'filled')
colormap(turbo)
% return
tic
xv = linspace(min(x), max(x), size(plot,1));
yv = linspace(min(y), max(y), size(plot,1));
F = scatteredInterpolant(x, y, z);
[X,Y] = ndgrid(xv, yv);
Z = F(X,Y);
figure
surfc(X, Y, Z, 'EdgeColor','none')
T = delaunay(X, Y);
figure
trisurf(T, X, Y, Z, 'EdgeColor','k', 'LineWidth',0.2)
toc

This also ended up crashing my computer.

.

Subham Prasad 2024-7-9，5:50

Thank you for the answer. The scatter plot (2D) will do for now. I think contourf kind of plot is not possible in this case due to the randomness of data. Had it been a sequential set of data points it would have worked. Do let me know if you find any way to get contour plot kind of image. I am attaching an example of how the plot should have looked and also one image to show how it looks with scatter plot (2D).

请先登录，再进行评论。