Environmental Engineering Reference
In-Depth Information
% chord(i) - chord of blade element i
% twist(i) - twist of blade element i
r_hub = r_hub_in/r_tip;
delr = ( 1.0 - r_hub)/nbes;
for i = 1:nbes
if (i == 1)
rad(1) = r_hub + delr/2;
else
rad(i) = rad(i-1) + delr;
end
chord(i) = chd(rad(i));
twist(i) = t1(rad(i)) + pitch;
end
data_out =[rad' chord' twist'];
% Make the first line of the data contain nbes r_tip and r_hub
data_out =[nbes r_tip r_hub_in;data_out];
save('rad_ch_tw.dat','data_out','-ascii')
function out = t1(x)
% Twist distribution. Fit is 4th order poly to r/R = 0.7, then linear.
if (x <= 0.7)
out = 54.16632 - x*(307.42939 -x*(719.549614 - …
x*(785.971096-x*326.673372)));
else
out = 5.318999 - 7.059999*x;
end
return
function out = chd(x)
% Chord distribution - 4th order poly fit to data in Anderson et al.
out = 5*(0.16165732 - x*(0.5847727 - x*(1.0327255 -…
x*(0.8756711 -x*0.2844545))))/3;
return
The running and output of
tcdist.m
is shown in the box below. Note that the
double greater than signs ''[['' is the Matlab command prompt:
Search WWH ::
Custom Search