3D surface plots
3 views (last 30 days)
Show older comments
ANOSH PHIROZ AMARIA
on 1 Oct 2018
Commented: Shafahat Ali
on 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 Comment
Shafahat Ali
on 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.
Accepted Answer
ANKUR KUMAR
on 1 Oct 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 Comments
ANKUR KUMAR
on 2 Oct 2018
Edited: ANKUR KUMAR
on 2 Oct 2018
Just use max command in the place of nanmean.
surf(xx,yy,max(P,[],3)','FaceAlpha',0.5,'EdgeColor','none')
More Answers (0)
See Also
Categories
Find more on Wind Power in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!