1)
Cd0 = 0.0358;
k = 0.0404;
m = 2180;
rho = 1.225;
S = 15;
W = m*9.81;
Cl = sqrt(Cd0/k);
Cd = Cd0 + k*(Cl)^2;
x = Cl/Cd;
V = sqrt((2*w)/(rho*S*Cd));
Tr = (1/2)*rho*(V^2)*S*Cd;
Pr = Tr*V;
display(x)
display(V)
display(Pr)
2)
Cd0 = 0.0358;
k = 0.0405;
m = 2180;
rho = 1.225;
S = 15;
W = m*9.81;
v = 10:0.01:100;
v1 = v';
Cl = (2*W)/(rho*S.*(v1.^2));
Cdi = k.*Cl.^2;
Cd = Cd0+k.*Cl.^2;
A=rho*S*Cd0/2;
B=(2*k*W^2)/(rho*S);
Tre=A*v1.^2;
Tri=B./v1.^2;
Tr=Tre+Tri;
TA = W./(Cl./Cd);
figure;
subplot(2,1,1);plot(v1,Tr);
hold on;
plot(v1,TA);
legend("Tr","TA");
Pr = Tr.*v1;
PA = TA.*v1;
hold off;
subplot(2,1,2);plot(v1,Pr);
hold on;
plot(v1,PA);
● legend("Pr","PA")
2)
% Given values
S = 15;
m = 2180;
W = m * 9.81;
CDo = 0.0358;
k = 0.0405;
rho_inf = 1.225;
V_inf1 = 1:400;
% Calculate Thrust Required (TR) and Power Required (PR) as functions of velocity
TR = 0.5 * rho_inf * (V_inf1.^2) * S * CDo + (2 * k * W^2) ./ (rho_inf * (V_inf1.^2) * S);
PR = TR .* V_inf1;
% Thrust Available (TA) and Power Available (PA)
TA = W * ones(size(V_inf1));
PA = TA .* V_inf1;
% Find min thrust and power required
[min_TR, min_TR_idx] = min(TR);
[min_PR, min_PR_idx] = min(PR);
% Find stall velocity
stall_velocity = sqrt((2 * W) / (rho_inf * S * CDo));
% Create plots
figure;
% Plot Thrust Required and Thrust Available vs Velocity
subplot(2,1,1);
plot(V_inf1, TR, 'b', 'LineWidth', 1.5);
hold on;
plot(V_inf1, TA, 'r', 'LineWidth', 1.5);
plot(V_inf1(min_TR_idx), min_TR, 'bo'); % Mark min Thrust Required
plot(V_inf1(min_TR_idx), TA(min_TR_idx), 'ro'); % Mark min Thrust Required
plot(stall_velocity, TR(round(stall_velocity)), 'gx'); % Mark stall velocity
xlabel('Velocity (m/s)');
ylabel('Thrust (N)');
title('Thrust Required and Available vs Velocity');
grid on;
legend('Thrust Required', 'Thrust Available', 'Min Thrust Required', 'Min Thrust Available', 'Stall
Velocity');
% Plot Power Required and Power Available vs Velocity
subplot(2,1,2);
plot(V_inf1, PR, 'b', 'LineWidth', 1.5);
hold on;
plot(V_inf1, PA, 'r', 'LineWidth', 1.5);
plot(V_inf1(min_PR_idx), min_PR, 'bo'); % Mark min Power Required
plot(V_inf1(min_PR_idx), PA(min_PR_idx), 'ro'); % Mark min Power Required
plot(stall_velocity, PR(round(stall_velocity)), 'gx'); % Mark stall velocity
xlabel('Velocity (m/s)');
ylabel('Power (W)');
title('Power Required and Available vs Velocity');
grid on;
legend('Power Required', 'Power Available', 'Min Power Required', 'Min Power Available', 'Stall
Velocity');
% Mark maximum velocity
for i = 1:2
subplot(2,1,i);
plot(max(V_inf1), max(get(gca, 'ylim')), 'kx'); % Mark max velocity
end
3)
Cd0 = 0.0358;
k = 0.0405;
s = 15;
d = 1.225;
m = 2180;
w= m*9.81;
vi=10:0.0001:100;
v=vi';
A=d*s*Cd0/2;
B=(2*k*w^2)/(d*s);
Cde=A*v.^2;
% plot(v,Cd)
Cdi=B./v.^2;
CD = Cde + Cdi;
plot(v,Cde,v,Cdi,v,CD)
xlabel('velocity')
ylabel('drag')
legend('PARASITE DRAG','INDUCED DRAG','TOTALÂ DRAG')
