В релизе Aero-PM 0.2.4 поправлена ошибка с минусом. Она хоть и не приводила к неверным результатам, сильно нарушала логику алгоритма. Теперь — все правильно. Рассматривается механизация — закрылок. Обычный отклоняемый закрылок теперь применяется редко, но его проще всего смоделировать. Два расчета: закрылок убран (FlappedFoil_1.txt), закрылок выпущен (FlappedFoil_2.txt). Профиль — без названия (несуществующий).
Механизация — закрылок (входные файлы)
Файл FlappedFoil_1.txt:
FlappedFoil_1 0.0000 0.0000 0.1000 0.0500 0.3000 0.0750 0.6000 0.0500 0.7000 0.0250 0.9000 0.0000
Файл FlappedFoil_2.txt:
FlappedFoil_2 0.0000 0.0000 0.1000 0.0500 0.3000 0.0750 0.6000 0.0500 0.7000 0.0250 0.8500 -0.1000
Запуск и результаты
Анализ результатов подтверждает увеличение подъемной силы при отклонении закрылка для неизменной скорости полета. Этим обусловлена необходимость закрылков самолета. Закрылки служат для увеличения подъемной силы на малых скоростях (взлетных или посадочных).
Закрылок (не щелевой) сдвигает кривую Cy(α) вверх на определенное расстояние (dCy). Относительный (%) прирост подъемной силы на малых углах атаки больше, т.к. величина сдвига кривой больше, чем Cy на малых углах. В решенной задаче угол атаки равен нулю. На рисунке ниже (условный пример) можно видеть, что при нулевом угле атаки расстояние от начала координат до точки А в несколько раз меньше расстояния до точки Б.
В вывод добавил ещё приращение к скорости набегающего потока. Пока только в наборе точек сверху и снизу от профиля. Каждая строка — две компоненты приращения по осям x и y. Сначала — верхние точки, а после пропуска строки — точки под профилем. По 11 точек.
>> Run(2, 0, 1.225, 'FlappedFoil_1.txt'); FlappedFoil_1 Angle of attack: 0 Lift, Moment, Moment1/4chord, LiftCoeff, PressureDiff: 1.2674 0.4935 0.2083 0.5748 5.2077 Velocities: 0.0756 0.2030 0.1531 0.2226 0.2340 0.1920 0.2997 0.1393 0.3523 0.0386 0.3178 -0.0764 0.2524 -0.1235 0.2070 -0.1475 0.1595 -0.1648 0.1138 -0.1643 0.0791 -0.1521 -0.0822 0.1367 -0.1254 0.1445 -0.1745 0.1276 -0.2126 0.0888 -0.2317 0.0381 -0.2293 -0.0130 -0.2122 -0.0547 -0.1876 -0.0866 -0.1574 -0.1092 -0.1250 -0.1201 -0.0959 -0.1205 >> Run(2, 0, 1.225, 'FlappedFoil_2.txt'); FlappedFoil_2 Angle of attack: 0 Lift, Moment, Moment1/4chord, LiftCoeff, PressureDiff: 4.9949 1.7981 0.7367 2.3985 27.4879 Velocities: 0.5164 0.5743 0.6796 0.5142 0.7988 0.3841 0.8691 0.2600 0.9517 0.0963 0.9216 -0.1278 0.8079 -0.2461 0.7461 -0.3128 0.6863 -0.4099 0.5790 -0.4982 0.4471 -0.5366 -0.4957 0.4148 -0.5834 0.3507 -0.6452 0.2663 -0.6832 0.1766 -0.7029 0.0856 -0.7045 -0.0022 -0.6977 -0.0805 -0.6941 -0.1619 -0.6779 -0.2677 -0.6106 -0.3846 -0.4926 -0.4593 >>