Representation of wall shear stress in 3D

Question

Jacopo Rossi 2022-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688079-representation-of-wall-shear-stress-in-3d

回答： VINAYAK LUHA 2023-11-24

I have a set of data points (in 3D coordinates, so a nx3 matrix) that represents the boundary of a region, then I have a nx1 matrix that represents the magnitude of wall shear stress in the corresponding point. I want to represent these data as in the figure I attached. I also attach a set of coordinates and the corresponding set of magnitude values.

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

请先登录，再进行评论。

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

Answer 1

VINAYAK LUHA 2023-11-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688079-representation-of-wall-shear-stress-in-3d#answer_1359192

在 MATLAB Online 中打开

Hi Jacopo,

I understand that you have a set of coordinate points representing the boundary of a region and an array indicating the shear stress values at these points. Further, you wish to visualize this data in MATLAB using color mapping to represent the magnitude at the respective points.

Here is how you can create a 3D surface plot with color mapping based on magnitude values-

Create a 3D delaunay triangulation based on the input coordinates using the "delaunay" function.
Normalize the shear stress magnitude values to the range [0, 1] for color mapping.
Next, Create a 3D surface plot with color mapping based on the normalized magnitude values using the "trisurf" function.
Finally, use a colormap and include a colorbar to represent the magnitude values.

The following code demonstrates the aforementioned method-

load('coordinates.mat');  
load('wssdata.mat');  
tri = delaunay(coordinates(:,1), coordinates(:,2), coordinates(:,3));
normalized_magnitude = (wssdata - min(wssdata)) / (max(wssdata) - min(wssdata));
figure;
trisurf(tri, coordinates(:,1), coordinates(:,2), coordinates(:,3), normalized_magnitude, 'FaceColor', 'interp', 'EdgeColor', 'none');
colorbar('Ticks', linspace(min(wssdata), max(wssdata), 5), 'TickLabels', arrayfun(@num2str, linspace(min(wssdata), max(wssdata), 5), 'UniformOutput', false));
colormap('jet');  
xlabel('X-axis');
ylabel('Y-axis');
zlabel('Z-axis');
title('3D Surface Plot with Color-Mapped Magnitude Values');