How to arrainge data to use streamline?

6 次查看（过去 30 天）
David Winthrop
David Winthrop 2020-2-13
评论： MS 2020-5-11
Hello,
I have data exported from a finite element software that I want to use in MATLAB to plot streamlines. The format is
xc,yc,zc,vx,vy,vz
where xc, yc and zc are (M-by-1) vectors of the coordinates of the nodes. vx, vy and vz are the components of the velocity at those coordinate locations.
The help for streamline says the X, Y and Z inputs must be monotonic, but my vectors are not sorted. I think if I sort them to satisfy the monotonic condition then they will not match the vx, vy and vz sorting, and since it is very unlikely they will all sort the same way, there is no way to keep them meatched to the vx, vy and vz.
How can I get my data in a form that streamline can use?
Thanks

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

回答（3 个）

KSSV
KSSV 2020-2-13
I guess what you have is unstructured FEM mesh data. YOu need to convert it to grided data. You can do inteprolation using gridddata. Refer this link:
https://in.mathworks.com/help/matlab/ref/griddata.html
David Winthrop
David Winthrop 2020-2-13
编辑：David Winthrop 2020-2-13
I tried to use griddata but there are some results that MATLAB can't seem to find, such as the re-circulation zones. Here is the streamline plot from the FEA software. I have annotated to show the coordinate system used.
Here is the code I used. It reproduces some of the stream lines but not the re-circulation lines. Any idea if I am doing something wrong?
sf = xlsread('stream_dat2.csv');
f = min( sf(:,1:2), [], 1 ); %// XY grid starting points
t = max( sf(:,1:2), [], 1 ); %// XY endpoints
[X,Y] = meshgrid( linspace( f(1), t(1), 50 ), linspace( f(2), t(2), 50 ) ); %//grid
U = griddata( sf(:,1), sf(:,2), sf(:,3), X, Y );
V = griddata( sf(:,1), sf(:,2), sf(:,4), X, Y );
startx = [0,-0.001,-0.0014]; %// define streamline starting points
starty = [-0.002,0,-0.0024]; %//
S = streamline( X, Y, U, V, startx, starty );
set(S,'color','r','linewidth',4)
hold on
plot(sf(:,1),sf(:,2),'*')
axis equal
set(gca,'ydirection','reverse')
view([-90,90])
xlabel('X')
ylabel('Y')
I have attached the CSV file as well.
  1 个评论
KSSV
KSSV 2020-2-15
Increase the number in linspace. You have taken only 50. Keep increasing it and check.

请先登录，再进行评论。

David Winthrop
David Winthrop 2020-2-14
Bumping in the hopes that someone can help.
  1 个评论
MS
MS 2020-5-11
Hi David,
I am facing the same issue as you. My streamline plots also not picking up the recirculation zone. Would you able to solve your problem in the streamline plot. Please let me know how you solved the issue.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Vector Volume Data 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by