Вычисление лучше всего скользит количества

Открыть скрипт

В этом примере показано, как выполнить вычисления скольжения для Примера следующего Cessna 172 9.1 в Ссылке 1 использование программного обеспечения Aerospace Toolbox™.

Лучшие вычисления скольжения вводят значения (скорость и угол скольжения), которые минимизируют, перетаскивают и максимизируют, перетаскивают лифт отношение (также названный отношением скольжения).

Технические требования самолета

Параметры самолета объявляются можно следующим образом.

W = 2400; % weight, lbf
S = 174;  % wing reference area, ft^2;
A = 7.38; % wing aspect ratio
C_D0 = 0.037; % flaps up parasite drag coefficient
e = 0.72; % airplane efficiency factor

Условия

Установите текущие условия самолета. Угол крена (phi) является нулем для этого случая.

h = 4000; % altitude, ft
phi = 0; % bank angle, deg

Преобразуйте высоту в метры с помощью convlength. Атмосферные вычисления на следующем шаге требуют значений в метрических модулях.

h_m = convlength(h,'ft','m');

Вычислите атмосферные параметры на основе высоты с помощью atmoscoesa:

[T, a, P, rho] = atmoscoesa(h_m, 'Warning');

Преобразуйте плотность от метрики до английских модулей с помощью convdensity:

rho = convdensity(rho,'kg/m^3','slug/ft^3');

Лучшие данные о скольжении

Лучшая скорость скольжения вычисляется с помощью следующего уравнения. TAS (истинная скорость полета в ногах в секунду) является скоростью самолета относительно окружающей массы воздуха.

$TAS_{bg}=\sqrt{\frac{2W}{\rho S}}\times \lbrack\frac{1}{4{C_{D_0}}^2+C_{D_0}\pi e A\cos^2{\phi}}\rbrack^\frac{1}{4}$

TAS_bg = sqrt((2*W) / (rho*S))...
         *(1./(4*C_D0.^2 + C_D0.*pi*e*A*cos(phi)^2)).^(1/4); % TAS, fps

Преобразуйте скорость от кадр/с до kts использование convvel. KTAS является истинной скоростью полета в узлах.

KTAS_bg = convvel(TAS_bg,'ft/s','kts')';

Преобразуйте KTAS в KCAS использование correctairspeed. KCAS (калиброванная скорость полета в узлах) является скоростью, откорректированной для ошибки положения и инструментальной погрешности. Эта ошибка положения прибывает из погрешностей в статических измерениях давления в различных точках в конверте рейса.

KCAS_bg = correctairspeed(KTAS_bg,a,P,'TAS','CAS')';

Лучший угол скольжения вычисляется с помощью:

$\sin{\gamma_{bg}}=-\sqrt{\frac{4C_{D_0}}{\pi e A\cos^2{\phi}+4C_{D_0}}}$

Это - угол между углом тангажа и землей, которая обеспечивает самое высокое отношение L/D.

gamma_bg_rad = asin( -sqrt((4.*C_D0')./(pi*e*A*cos(phi)^2 + 4.*C_D0')) );

Преобразуйте угол скольжения от радианов до степеней с помощью convang:

gamma_bg = convang(gamma_bg_rad,'rad','deg');

Лучшее перетаскивание скольжения вычисляется с помощью:

$D_{min}=D_{bg}=\frac{1}{2}\rho (TAS^2_{bg})S(2C_{D_0})=-W\sin{\gamma_{bg}}$

D_bg = -W*sin(gamma_bg_rad);

Лучший лифт скольжения вычисляется с помощью:

$L_{bg}=L_{max}=W\cos{\gamma_{bg}}=\sqrt{W^2-D^2_{bg}}$

L_bg =  W*cos(gamma_bg_rad);

Вычислите динамическое давление с помощью dpressure:

qbar = dpressure([TAS_bg' zeros(size(TAS_bg,2),2)], rho);

Вычислите перетаскивают и снимают содействующее использование:

$C_{D_{bg}}=\frac{D_{bg}}{\bar{q}S}$

$C_{L_{bg}}=\frac{L_{bg}}{\bar{q}S}$

C_D_bg = D_bg./(qbar*S);
C_L_bg = L_bg./(qbar*S);

Сводные данные лучших значений скольжения

Вот лучшие значения скольжения:

$KCAS_{bg}=71.9 \hspace{1pt} KCAS$

$\gamma_{bg}=-5.38\deg$

$C_{D_{bg}}=0.074$

$C_{L_{bg}}=0.7859$

$D_{bg}=224.9 \hspace{1pt} lbf$

$L_{bg}=2389.4 \hspace{1pt} lbf$

Верификация

Эти графики показывают, перетаскивают и перетаскивают лифт графики отношения для самолета в зависимости от KCAS. Графики используются, чтобы проверить лучшие вычисления скольжения.

Установите область значений скоростей полета и преобразуйте в KCAS использование convvel и correctairspeed:

TAS = (70:200)'; % true airspeed, fps
KTAS = convvel(TAS,'ft/s','kts')'; % true airspeed, kts
KCAS = correctairspeed(KTAS,a,P,'TAS','CAS')'; % corrected airspeed, kts

Вычислите динамическое давление для новых скоростей полета с помощью dpressure:

qbar = dpressure([TAS zeros(size(TAS,1),2)], rho);

Вычислите паразитное использование перетаскивания:

$D_p=\frac{1}{2}\rho S C_{D_0}(TAS^2)$

Dp = qbar*S.*C_D0;

Вычислите вызванный, перетаскивают использование:

$D_i=\frac{2W^2}{\rho S\pi e A}\frac{1}{(TAS^2)}$

Di = (2*W^2)/(rho*S*pi*e*A).*(TAS.^-2);

Вычислите общее использование перетаскивания:

$$D=D_p+D_i$$

D = Dp + Di;

Аппроксимированный лифт как вес (принимающий маленький угол скольжения и маленький угол нападения). На этой скорости, принимая

$C_L=2\pi\alpha$

и использование

$C_{L_{bg}}$

сверху, угол нападения является приблизительно 7 градусами. Добавление угла угла тангажа (i.e. лучший угол скольжения), сверху показывает тангаж фюзеляжа (угловая theta отношения), чтобы быть приблизительно 2 градусами.

L = W;

Постройте L/D по сравнению с KCAS

Как ожидалось максимальный L/D происходит приблизительно при лучшей скорости скольжения, вычисленной выше.

h1 = figure;
plot(KCAS,L./D);
title('L/D vs. KCAS');
xlabel('KCAS'); ylabel('L/D');
hold on
plot(KCAS_bg,L_bg/D_bg,'Marker','o','MarkerFaceColor','black',...
    'MarkerEdgeColor','black','Color','white');
hold off
legend('L/D','L_{bg}/D_{bg}','Location','Best');
annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','KCAS_{bg}');

Постройте паразита, вызванные, и общие кривые перетаскивания

Заметьте минимальное общее перетаскивание (i.e. D_bg), происходит приблизительно при той же лучшей скорости скольжения, вычисленной выше.

h2 = figure;
plot(KCAS,Dp,KCAS,Di,KCAS,D);
title('Parasite, induced, and total drag curves');
xlabel('KCAS'); ylabel('Drag, lbf');
hold on
plot(KCAS_bg,D_bg,'Marker','o','MarkerFaceColor','black',...
    'MarkerEdgeColor','black','Color','white');
hold off
legend('Parasite, D_p','Induced, D_i','Total, D','D_{bg}','Location','Best');
annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','KCAS_{bg}');

close(h1,h2);

Ссылка

[1] Lowry, J. T., "Performance of Light Aircraft", AIAA(R) Education Series,
    Washington, DC, 1999.

Документация