Scribd
Upload
12 views4 pages

Flight Dynamics

Codee

Uploaded by

attom997.799
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
12 views4 pages

Flight Dynamics

Codee

Uploaded by

attom997.799
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

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')

You might also like