How can i create efficiency isolines in a compressor map?

2018 9 12
2 个回答
回答已采纳
更新时间：2020 10 19
23 次查看（30 天）

0 个投票

I have mass flow, pressure ratio and efficiency values of about 30 points for 3 different rotation speed. I need to plot a compressor map as in the link by using these datas: https://www.uniplot.de/_images/metafiles-compressormap-80.png
I used interp1 with spline as method to create constant rotation speed lines (black lines in link above). And now i need to create efficiency isolines. I tried with different commands(for example; contour, griddata) but it did not work so far.
How can i create efficiency isolines?

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

 采纳的回答

jonas
jonas 2018-9-12

0 个投票

This is what I'm getting from the data you uploaded.
%%Load data
rpm{1}=xlsread('data_compressor_map.xlsx','A3:D10');
rpm{2}=xlsread('data_compressor_map.xlsx','A14:D25');
rpm{3}=xlsread('data_compressor_map.xlsx','A29:D33');
data=vertcat(rpm{:});
x=data(:,2)
y=data(:,4)
z=data(:,3)
%%Interpolate on grid
[X,Y]=meshgrid(min(x):.01:max(x),min(y):.01:max(y))
Z=griddata(x,y,z,X,Y)
%%Plot contour
contourf(X,Y,Z);hold on
%%Plot lines
plot([x(1:8);NaN;x(9:20);NaN;x(21:end)],[y(1:8);NaN;y(9:20);NaN;y(21:end)],'-k','linewidth',2)
cb=colorbar

更多回答（1 个）

chapuisat
chapuisat 2018-9-16

0 个投票

Your answer helped me a lot! Thank you!
I think the map in the link below looks visually better with regards to efficiency area.
How can i change efficiency area as in the link below?
https://www.google.de/imgres?imgurl=x-raw-image%3A%2F%2F%2F405abd0c340c8e0e2f65d63f06a2eb7b6cc5c02a70d88a6f246ae45c93063160&imgrefurl=http%3A%2F%2Fpublications.lib.chalmers.se%2Frecords%2Ffulltext%2F196400%2F196400.pdf&docid=CpI22Qr1H9koZM&tbnid=Do4MNAy9WkBVRM%3A&vet=12ahUKEwiEmYOb6b_dAhVGkCwKHancCrU4ZBAzKBYwFnoECAEQFw..i&w=727&h=625&bih=754&biw=1536&q=compressor%20map%20contour&ved=2ahUKEwiEmYOb6b_dAhVGkCwKHancCrU4ZBAzKBYwFnoECAEQFw&iact=mrc&uact=8

3 个评论
显示 1更早的评论 隐藏 1更早的评论

jonas
jonas 2018-9-16
编辑：jonas 2018-9-27
The link brought me to a master thesis pdf but I think I know which fig you are referring to. Basically its a surface plot with some added features.
There are a million things you can change to make the surface look more visually appealing. Here's an example of some things you could try:
%%Load data
rpm{1}=xlsread('data_compressor_map.xlsx','A3:D10');
rpm{2}=xlsread('data_compressor_map.xlsx','A14:D25');
rpm{3}=xlsread('data_compressor_map.xlsx','A29:D33');
data=vertcat(rpm{:});
x=data(:,2)
y=data(:,4)
z=data(:,3)
%%Interpolate on grid
[X,Y]=meshgrid(min(x):.002:max(x),min(y):.01:max(y));
Z=griddata(x,y,z,X,Y)
%%Tighten surface to avoid extrap
b = boundary(x,y,1)
in=inpolygon(X,Y,x(b),y(b));
Z(~in)=NaN;
%%Plot contour
h=pcolor(X,Y,Z);hold on
set(h,'linestyle','none')
%%Plot lines
plot([x(1:8);NaN;x(9:20);NaN;x(21:end)],[y(1:8);NaN;y(9:20);NaN;y(21:end)],'o-k','linewidth',1.5,'markerfacecolor','k')
cb=colorbar
%%Plot boundary
plot(x(b),y(b),'k','linewidth',1)
xlim([.2 1.8])
ylim([1 4.5])
set(gca,'layer','top')
chapuisat
chapuisat 2018-9-27
How did you define x(b) and y(b)?
jonas
jonas 2018-9-27
编辑：jonas 2018-9-27
If I remember correctly I used.
b = boundary(x,y,1);
That line got lost somehow, sorry.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by