3D surface plots
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
ANOSH PHIROZ AMARIA
am 1 Okt. 2018
Kommentiert: Shafahat Ali
am 11 Jan. 2022
I want to plot a 3D surface plot for Power (P) of a wind turbine while its gear ratio (GR) and Diameter (D) are varying over a range of 14.6:2:20.6 amd 16.5:1:20.5 respectively. I am using a nested for loop to get the power at each of the GR and D. How do I get the surface plot of GR vs D vs P. Code is as below.
clear all close all clc %WIND TURBINE SYSTEM PARAMETERS % D = 16.5:0.1:20.5; %Rotor diameter in m %GR= 14.6:0.2:20.6; %Gear Ratio of the gearbox connecting rotor to generator beta=0; % Angle of attack for maximum efficiency of the turbine %AIR PROPERTIES rho = 1.27; % Density of air kg/m^3 Vw = [0.001:0.01:35]; Wgen = 1800; %Synchronous Speed of the Generator in rpm for D=16.5:0.1:20.5 for GR=14.6:0.2:20.6 Vt= ((Wgen.*pi/30)./GR).*(D/2); Lambda = Vt./Vw; %Heier Approach to Calculate the Coefficient of Performance (Cp) %Matrix of Coefficient C=[0.5 116 0.4 0 0 6 21 0.08 0.035]; Lambda_i = ((1./(Lambda+(C(8).*beta)))-(C(9)./((beta.*beta.*beta)+1))).^(-1); Cp=0.5*((C(2)./Lambda_i)-(C(3).*beta)-C(6)).*exp(-C(7)./Lambda_i); %Calculation of Power P=0.5*rho*((pi/4)*(D*D/2).*(Vw.^3)).*(Cp); end end
1 Kommentar
Shafahat Ali
am 11 Jan. 2022
i want to make 3d plots with input cutting speed and feed rate of 3 levels. Output of three levels can be anything. How can i make it please tell me so it will come in grid shape.
Akzeptierte Antwort
ANKUR KUMAR
am 1 Okt. 2018
I have made few modifications. In your program, you are not saving the output of P for each diameter and gear ratio.
clear all
close all
clc
%WIND TURBINE SYSTEM PARAMETERS
% D = 16.5:0.1:20.5; %Rotor diameter in m
%GR= 14.6:0.2:20.6; %Gear Ratio of the gearbox connecting rotor to generator
beta=0; % Angle of attack for maximum efficiency of the turbine
%AIR PROPERTIES
rho = 1.27; % Density of air kg/m^3
Vw = [0.001:0.01:35];
Wgen = 1800; %Synchronous Speed of the Generator in rpm
i=0;
for D=16.5:0.1:20.5
i=i+1;
j=0;
for GR=14.6:0.2:20.6
j=j+1;
Vt= ((Wgen.*pi/30)./GR).*(D/2);
Lambda = Vt./Vw;
%Heier Approach to Calculate the Coefficient of Performance (Cp)
%Matrix of Coefficient
C=[0.5 116 0.4 0 0 6 21 0.08 0.035];
Lambda_i = ((1./(Lambda+(C(8).*beta)))-(C(9)./((beta.*beta.*beta)+1))).^(-1);
Cp=0.5*((C(2)./Lambda_i)-(C(3).*beta)-C(6)).*exp(-C(7)./Lambda_i);
%Calculation of Power
P(i,j,:)=0.5*rho*((pi/4)*(D*D/2).*(Vw.^3)).*(Cp);
end
end
[xx,yy]=meshgrid(16.5:0.1:20.5,14.6:0.2:20.6);
surf(xx,yy,nanmean(P,3)','FaceAlpha',0.5,'EdgeColor','none')
xlabel('Rotor diameter')
ylabel('Gear Ratio')
colorbar
3 Kommentare
ANKUR KUMAR
am 2 Okt. 2018
Bearbeitet: ANKUR KUMAR
am 2 Okt. 2018
Just use max command in the place of nanmean.
surf(xx,yy,max(P,[],3)','FaceAlpha',0.5,'EdgeColor','none')
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Wind Power finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!