14/4/25 9:07 PM C:\Users\mayan\Downloads\SPVWES_EXP2.
m 1 of 1
clc
clear all
close all
%defining constants
cl(1)=0.95;
cl(2)=1.2;
cl(3)=0.75;
cl(4)=0.46;
cd(1)=0.0105;
cd(2)=0.0143;
cd(3)=0.0092;
cd(4)=0.0078;
%defining parameters
v=6;
diam=9;
w=100;
blade=4;
tsr=5.23;
ro=1.25;
chord=0.45;
beta=5;
axis=0.5;
aerofoil=1;
for i=1:4
I(i)=atan(v/(2*pi*i*(w/60))); %angle of attack
inci=I(i)-beta; %angle of incidence
dp(i)=3*0.5*ro*chord*v^3*cot(I(i))*(1+(cot(I(i)))^2)*(cl(i)*sin(I(i))-cd(i)*cos(I
(i))); %power generated in i(th) segment
df(i)=3*0.5*ro*chord*v^2*(1+(cot(I(i)))^2)*(cl(i)*cos(I(i))+cd(i)*sin(I(i))); %
thrust generated in i(th) segment
end
dp
dP=sum(dp)
df
dF=sum(df)