Аэродинамика летающих крыльев

111

to GreenGo
ну а как по Вашему kirk30?

GreenGo
111:

ну а как по Вашему kirk30?

Толстоват по сравнению с PW106.
Насколько я понимаю, профиль это компромисс под определенные задачи. Может для медленно летающего прямого ЛК это то, что нужно? 😃

111

ну а как Cm по Alfa у Kirk30 по сравнению со 106, Ваше мнение?

GreenGo
111:

ну а как Cm по Alfa у Kirk30 по сравнению со 106

Все таки, нужно относительно какой то модели говорить. Для стреловидного ЛК, например, момент не так важен, а сопротивление оно везде мешает.

Хочу, в следующий раз, кратко сравнить несколько профилей для стреловидных ЛК (в этой теме), и не против квалифицированных замечаний.

GreenGo

Профили для стреловидных ЛК.

  1. Профили для стреловидных ЛК, с похожими характеристиками.

MH64, MH22 (delta400), PW51, DeltaN8 «быстрые» профили, с малым коэффициентом сопротивления, и небольшим отрицательным моментом. Одинаково не очень хорошо работают на малой скорости (PW51 и DeltaN8 немного лучше). Хорошо подходят для стреловидных ЛК, из за малого момента профиля и хорошего соотношения Cl/Cd на малых углах атаки.
На двухметровой модели стреловидного ЛК (в XFLR5) «побеждает» профиль с минимальным сопротивлением. Потому, что остальные характеристики мало отличаются.

  1. Профили для стреловидных ЛК, различающихся по скорости максимального качества.

PW106 лучше ведет себя на малой скорости, MH64 хорошо на большой и хуже на малой скорости.
DeltaN8 занимает среднее место, т.к. почти не отстает от PW106 на малой скорости и от MH64 на большой (даже слегка опережает).

В этом сравнением, я хотел показать, по каким характеристикам можно оценивать профили для стреловидных ЛК. Коэффициент момента Cm не так важен (график 4), но желательно, чтобы он не был большим (отрицательным). Качество Cl/Cd (график 5), имеет большое значение на малых углах атаки, поэтому симметричные профили не очень хороший выбор для стреловидных ЛК. Сопротивление Cd (график 6)на малых углах атаки так же имеет большое значение поэтому хорошо подходят тонкие профили 8-9%. Поведение на малой скорости и на больших углах атаки (пилотаже) зависит от критического угла атаки и максимального значения коэффициента подъемной силы Cl (график 2). По графику Cl от Cd (график 1)можно оценить работу профиля на малых и отрицательных углах атаки (в инверте).

DeltaN8 на практике не проверялся, поэтому координаты пока не выкладываю.

GreenGo

Для большей “наглядности” сравнение моделей стреловидных ЛК с профилями MH 64, PW51 и PW106 (размах 2000мм, вес 1300гр.):


График скорости снижения Vz в зависимости от скорости Vx. Скорость снижения меньше 0,4м/сек - это не хуже чем у классики, на таком же размахе.


График качества в зависимости от скорости.

Различие по максимальному качеству, намного заметнее чем разница в скорости снижения. MH 64 выглядит лучше, но его преимуществом будет сложно воспользоваться т.к. скорость на которой у него минимальное снижение очень близко к скорости срыва.

alix
Стрела_RoSa:

Выкладываю на русском полезную книжку про ЛК depositfiles.com/ru/files/mxloqhgn4?redirect книжка Капковского.Читаю - настроение поднимается и душа радуется!Рекомендую всем!

Ссылка на эту книгу есть на конференции. Если не ошибаюсь, даже в этой ветке. Однако там не забыли отметить, что Капковский в своей книге рассматривает СВОБОДНОЛЕТАЮЩИЕ крылья. И, к сожалению, начинающим моделистам не всегда ясна суть отличий. Так что титать её здорово, а вот советовать бы нужно с осторожностью 😉

GreenGo

Еще одна интересная на мой взгляд тема.

Анализ модели стреловидного ЛК в XFLR5.http://www.xflr5.sourceforge.net

Предисловие ​****Franchesco Meschina.
Многие моделисты знакомы с XFoil, известной программой Марка Дрелы и Гарольда Юнгена.
Эта программа выполняет симуляцию двухмерного воздушного потока вокруг изолированного профиля без учета многих 3D эффектов в том числе индуктивного сопротивления.

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

По причине своей большой важности не только для моделистов, но и для большой авиации эта «проблема» исследовалась уже давно, и для ее решения было предложено много методов. В том числе - VLM (vortex lattice method).

Не вдаваясь в подробности этого метода можно сказать, что в основе VLM, лежит разбиение несущих поверхностей ЛА на набор панелей. Каждая панель имеет присоединенный подковообразный вихрь. Используя несколько граничных условий можно рассчитать подъемную силу и сопротивление, создаваемое каждым присоединенным вихрем. Просуммировав вклады всех отдельных вихрей можно получить характеристики всей несущей поверхности.

По причине сложности расчетов в VLM, он не мог быть широко использован в «до цифровую эпоху». Только с появлением компьютеров после 60х годов прошлого столетия VLM показал себя как мощный инструмент и был использован для изучения и создания множества различных аэродинамических проектов.

Несколько лет назад метод стал доступен для моделистов, когда Марк Дрела и его группа выпустили программу AVL. Но аналогично как и в случае с программой Profili2 (Stefano Duranti), VLM стал более доступен для неспециалистов, с появлением программы Andre Deperrois XFLR5.

Программа Andre Deperrois внесла новый вклад в анализ. В то время как классический VLM метод принимает поток воздуха как не имеющий вязкости и поэтому не совсем реалистичен для малых, модельных чисел Re. В XFLR5 принимается, что вклад вязкого и не вязкого обтекания воздушным потоком, линейно не зависимы и VLM метод может быть дополнен XFoil анализом, для получения более реалистичной математической модели.

Deperrois предупреждает, что эта «гипотеза»- о независимости вкладов вязкого и невязкого обтекания, не подтверждается теорией аэродинамики, поэтому к результатам VLM нужно относится с осторожностью и они требуют дополнительной проверки. Поэтому пользователи программы XFLR5 должны относиться с пониманием некоторых ограничений математических методов положенных в ее основу.

В XFLR5 должна быть «загружена» геометрическая модель несущей поверхности и набор поляр используемых профилей, для диапазона чисел Re достаточного, чтобы перекрыть все полетные режимы.
Другими словами, прежде чем начать использовать программу XFLR5, пользователь должен найти файлы координат профилей которые он использует и подготовить размеры трехмерной модели ЛА.

Так же нужно помнить, что XFLR5 это просто симулятор для набора физических переменных и не может дать ответы кроме тех на которые мы готовы задать вопросы. Это удобный инструмент, но задача научить нас хорошо проектировать ЛА выходит за рамки его применения.

На много лучше подходить к программе с практической точки зрения, зная ее возможности и ограничения, и то какие ей нужно задавать вопросы. Поэтому возьмем практический случай, чтобы на его примере найти ответы на некоторые общие вопросы: определение поляры модели, анализ устойчивости по тангажу, анализ распределения подъемной силы вдоль размаха крыла, определение балансировочных отклонений управляющих поверхностей.

Часть 1. Моделирование тестовой модели.
В качестве тестовой модели возьмем стреловидное ЛК, размахом 1800мм, весом 1200 гр., со следующими “проектными” характеристиками: рабочая скорость 18 м/сек, минимальная скорость 7 м/сек, максимальная скорость 28 м/сек.

Моделирование начинается с профиля, поэтому считаем, что с этим мы уже определились, это будут: MH64мод. и PW75.
Когда профиль найден и загружен на ваш компьютер, можно запустить XFLR5 и начать с импорта файла координат. Из меню File выбираем Load File и открываем файл типа .dat, (.plr, .wpa).
Теперь, когда у нас есть профили, можно используя модуль XFoil, рассчитать их поляры. Из меню выбираем Application -> XFoil Direct Analysis.
1.Выбираем тип анализа Type 1 (для фиксированной скорости и величины хорды, и переменной величины коэффициента подъемной силы).
Для каждого профиля мы должны рассчитать набор поляр, перекрывающий диапазон чисел Re, соответствующий желаемым полетным режимам.
2. Выбираем из меню Polars -> Run Batch Analysis -> Edit List, и добавляем в список числа от 60000 до 800000, как на рисунке.

3.Дальше нужно задать диапазон и шаг изменения угла атаки. В разделе Range, выбираем alpha, min -3 град. , max 16 град, increment 1 град.
4.Нажимаем Analize и ждем несколько минут в зависимости от мощности компьютера.
По окончании расчета, мы получим «семейство» поляр, и сможем рассмотреть их подробнее, нажав на цифровой клавиатуре число соответствующее графику от 1 до 5.
После того как семейство поляр рассчитано для всех необходимых профилей, пора приступить к созданию трехмерной модели нашего ЛК. К этому моменту мы должны знать геометрические размеры, стреловидность, величину крутки и величину поперечного V.
Выбираем из меню Application -> Wing Design
Wing/Plane -> Define a Plane. В появившемся окне нам нужно ввести имя модели, затем выбрать
Main Wing -> Define. В этом окне мы вводим данные по трехмерной модели. Форма крыла задается в виде одного или нескольких трапецевидных панелей, для каждого задается: размах, хорда, смещение от передней кромки, поперечное V, крутка и профиль.
Наша модель вначале будет состоять только из одной панели. Введем ее размеры:
Для корневого профиля:
Pos – 0,
Chord – 270мм,
Foil Name - MH64mod.
Для концевого профиля:
Pos-900мм,
Chord – 160мм,
Offset – 390мм,
Twist - -4град.
Foil Name – MH64mod.

При расчете используется разбиение крыла на определенное количество панелей - Total VLM Panels. Чем больше панелей тем расчет точнее. Но большое количество панелей отнимает у программы много ресурсов, поэтому перед закрытием окна желательно нажимать на кнопку Reset VLM Mesh, для оптимизации количества панелей.
После закрытия окна, когда все размеры введены правильно, мы можем увидеть нашу тестовую модель в 3****D View режиме программы. В изометрии Iso или в одной из трех проекций – X, Y, Z.

В следующей части рассмотрим расчет и анализ поляры модели.

wing_1800test.zip

modsley
111:

ну это вряд ли, куда же он денется

А вы замечали, насколько по-разному реагируют на РН пилотажки и тренеры?

GreenGo

Анализ модели стреловидного ЛК в XFLR5.

Часть 2. 3D поляра модели.

Нам остался только один шаг до получения поляры модели.
Для расчета выбираем из меню: Polar -> Define Polar Analysis.

Мы должны определить несколько условий. Поэтому для начала мы должны решить какой вид анализа мы хотим выполнить:
для постоянной скорости - Fix Speed
для постоянной подъемной силы - Fix Lift
для постоянного угла атаки - Fix Alpha.
Для анализа планирующего полета, когда крыло создает столько подъемной силы сколько необходимо для поддержания горизонтального полета, хорошо подходит второй вариант - Fix Lift.
Так же мы должны указать вес модели – Plane Weight 1200гр.
Положение ЦТ относительно передней кромки крыла - Mom. Ref. location 215мм.
Далее выбираем метод расчета VLM или LLT. Метод LLT доступен только для изолированного крыла (без оперения).
Не забываем проверить «чекер» Viscous.

Подтверждаем внесение изменений нажимая на кнопку – Ok.

Теперь мы готовы произвести анализ. Но прежде чем нажать кнопку Analyze, нужно заполнить раздел Analysis в правой вертикальной панели, мы должны задать условия расчета:

Например: угол атаки от 0 до 11 градусов с шагом 0,5 градуса. Проверяем «чекер» Store op.points и запускаем анализ.

В процессе анализа будет выводится журнал в которым могут быть сообщения, что некоторые рабочие точки – **Op.**points выходят за пределы подготовленных поляр. С этим мы разберемся позже.

После завершения работы программы анализа, мы можем посмотреть на результаты ее работы в нескольких режимах просмотра:
Polar View –в виде графиков «поляр»,
3****D View – в виде трехмерного представления,
Operating Pints View – в виде графиков рабочих точек.

На данном этапе мы будем анализировать поляру модели в режиме Polar View.
Если не открыт нужный график (Vz от Vx), через контекстное меню выбираем Graph -> Variables (или на клавиатуре V) и задаем Vz для оси Y-axis и Xinf для оси X-axis.

Полученный график показывает зависимость скорости снижения от горизонтальной скорости модели.

Можно заметить, что минимальная скорость снижения 0,4 м/сек соответствует скорости модели около 10 м/сек - чуть больше чем скорость срыва.
Другая точка – точка максимального качества около 11,5м/сек.

Эти две точки соответствуют двум важным полетным режимам планера: полету на максимальную продолжительность и полету на максимальную дальность.

Другой удобный для анализа график – это зависимость качества Cl/Cd от скорости модели Vx.

На этом графике не только хорошо видна скорость максимального качества 11,5м/сек, но и величина максимального качества модели 24,5.

На этом этапе мы подошли к очень важному вопросу: на какой скорости должна лететь наша модель, с настройками «по умолчанию», заданным положением ЦТ и при нулевом отклонении элевонов?

Ответ может быть получен если посмотреть на график зависимости суммарного момента GCm от скорости Vx. Точка пересечения этого графика с осью X – это и есть скорость на которую настроена модель. В нашем случае это 18м/сек. Это скорость равномерного полета с данной центровкой и положением РВ.

График так же показывает, что более высокая скорость приведет к кабрированию (нос модели вверх) и потере скорости и более низкая скорость приведет к «пикированию» (нос модели вниз) и набору скорости.

То, что модель настроена на 18м/сек, а не на наивыгоднейший режим планирования (10 или 11м/сек), неудивительно так мы хотим получить не паритель, а быструю модель .

На скорости 18м/сек скорость снижения Vz будет выше – 0,8 м/сек, это означает, что для поддержания модели в воздухе в безмоторном полете нужен восходящий поток меньше 1 м/сек.

Часть 3. Анализ продолжительности полетного времени с использованием MotoCalc.

В желательных характеристиках модели, в части 1., была указана максимальная скорость 28м/сек. Предполагалась не скорость планирования, а скорость горизонтального моторного полета, при разумном расходе емкости аккумулятора.

Предположим, что в наличии есть такой моторсет: и 1200об/в 45А RimFire 3536 с винтом APC 9x6, регулятор Phoenix ICE 50 и аккумулятор 2200ma 4S. Попробуем определить продолжительность полетного времени нашей модели с этим моторсетом.

Можно конечно найти формулы, которые учитывают характеристики мотора, аккумулятора и эффективность воздушного винта, но мы сделаем проще – воспользуемся хорошим инструментом, предназначенным для этого, - программой MotoCalc. Программу и русификацию смотрте в дневнике у boroda_de.

Сначала нужно внести в программу характеристики мотора, аккумулятора и винта (характеристики модели использовать не будем).

Затем нажимаем кнопку «Рассчитать» и анализируем полученные данные. Когда мы не указываем параметры модели, MotoCalc рассчитывает только моторсет – получаемую тягу в зависимости от скорости и положения стика RC передатчика. Результаты представляются или в виде таблицы или в виде графика.

Исходными данными для анализа моторсета будут: скорость модели и сопротивление модели Drag(N) и соответствующая тяга мотора на этой скорости.
Определяем в программе XFLR5, по графику Drag(N) от Vx, значения сопротивления на скоростях 12м/сек, 18м/с и 28м/с . Значения увеличиваем на 15%, для учета сопротивления дополнительных конструктивных элементов на реальном ЛК.

Подбираем в MotoCalc значение Throtle (отклонение стика) для получения нужной тяги на этих скоростях и получаем ток и полетное время:


Полученное время выглядит слишком оптимистично, т.к. это время для идеального случая равномерного горизонтального полета, без внешних и управляющих воздействий. Если принять, что расход емкости аккумулятора при маневрировании на 30-50% больше, то значение полетного времени становится более реалистичным .

GreenGo

Часть 4. Распределение подъемной силы.

В предидущей части мы анализировали графики характеризующие поведение модели от переменных изменяющихся в полете, в частности скорости. Есть еще один способ анализа модели по результатам расчета XFLR5.

Каждая точка отображаемая на графиках поляр модели содержит не одно числовое значение, а набор данных включающий все рассчитанные программой величины. Этот набор называется Operating point.

Б олее того, в Operating point хранятся не только скалярные величины отдельных характеристик модели, а набор данных распределенных по всем единичным панелям разбиения VLM метода, по всей площади крыла.
XFLR5 позволяет нам использовать эти данные для интересного и довольно «фундаментального» исследования по аэродинамике.

Например, если мы посмотрим на набор кривых описывающих распределение коэффициента подъемной силы вдоль размаха крыла, мы сможем определить какой участок крыла «нагружен» сильнее и где начнется срыв потока.

Этот вид анализа выполняется в режиме программы Operating Point View. Если у нас только одна модель в программе то нам достаточно выключить Show Current Opp Only и включить все Operating Point режимом Show All Operating Point.

Если отображение графика Cl не включено, выберем его через контекстное меню. По изменению графика с ростом угла атаки видно увеличение и перераспределение коэффициента подъемной силы ближе к концу крыла. В случае отсутствия крутки этот эффект еще заметнее и срыв при этом происходит намного раньше.

На начало срыва потока XFLR5 указывает когда в логе расчета поляры модели, на больших углах атаки в появляется сообщение подобное такому:
Span pos = -608.83 mm, Re = 92 552, Cl = 1.00 could not be interpolated
Это означает, что программа не может рассчитать очередную Operating point, косвенно показывая, что приближается срыв потока. Из этой строчки мы можем определить место на крыле где начинается срыв и при каком числе Re это происходит. Это может дать нам подсказку как исправить ситуацию и отодвинуть срыв на большие углы атаки.

Еще один набор графиков которые полезно рассмотреть – это графики распределения локальной подъемной силы, Lift.

Для сравнения, полученного на нашей модели распределения, с «идеальным» эллиптическим распределением включим через контекстное меню его отображение - Show Elliptic Loading. Как видим распределение у нашей модели не идеальное - есть возможность увеличения качества модели за счет получения лучшего распределения.

Часть 5. Анализ устойчивости по тангажу.

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

Для получения графика характеризующего устойчивость, вернемся в режим Polar View.
Выберем из контекстного меню Graph -> Variables, переменную GCm для оси Y и Alpha для оси X.

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

Кривая проходит через ось X в точке 2,6град. Alpha, что говорит нам об устойчивом положении модели с таким углом атаки.

Если кривая будет наклонена с возрастанием момента при увеличении угла атаки, модель будет неустойчивой – увеличение Alpha будет приводить к увеличению кабрирующего момента.

Горизонтальная кривая говорит о нейтральной устойчивости модели.
Кривые на Рис. Сделаны при различном положении ЦТ. Когда мы сдвигаем ЦТ назад устойчивость модели уменьшается.

В соответствии с этими графиками, в нашей тестовой модели центр давления находится в точке расположенной в 232 мм от передней кромки крыла. ЦТ находится в 16мм впереди ЦД, что создает величину запаса устойчивости 8% САХ. Оптимальным для стреловидных ЛК, считается запас устойчивости от 1 до 5 процентов средней аэродинамической хорды.

Результаты расчета в WinLaengs4 незначительно отличаются (различие в расчете расположения ЦД, два миллиметра).

GreenGo

Часть 6. Управляющие поверхности.

В XFLR5 мы можем добавить нашей модели управляющие поверхности. Для модели обычной схемы анализ закрылков или флапернов может быть интересен для определения характеристик модели с измененной кривизной профиля, при отклонении флаперонов вниз или вверх.

В нашем случае мы рассмотрим необходимые отклонения элевонов в качестве РВ для балансировки, при настройке на различные полетные режимы и изменение характеристик ЛК с отклоненными элевонами.

Сначала необходимо добавить управляющие поверхности в профиль. Возвращаемся к началу и выбираем из меню Application -> Foil Direct Design.
Из контекстного меню выбираем Set Flap и выбираем: **T.**E. Flap, -1,5град., 75%.

Выполняем расчет семейства поляр для нового профиля, как мы делали это для основного профиля.

По окончании расчета мы увидим «мешанину» поляр. Их отображение можно убрать: View -> Polars -> Hide all polars.

Теперь включим отображение только двух поляр, для сравнения, выбрав нужный профиль и нужное число Re из списка вверху и поставив «чекер» Show Curve в панели справа внизу.
На графиках видна существенная разница по сравнению с профилем без отклонения управляющей поверхности.

Основное изменение в графиках это балансировочный положительный момент Cm, который появляется при отклонении элевона вверх, на графике Cm от Alpha.
Но кроме полезного балансировочного изменения момента произошли и нежелательные изменения в характеристиках профиля. Уменьшился коэффициент подъемной силы Cl и увеличился коэффициент сопротивления Cd, в результате уменьшилось соотношение Cl/Cd.

Посмотрим как это отразится на поведении модели, переходим в Application -> Wing Design.
Создадим копию нашей модели: Wing/Plane -> Current Wing/Plane -> Duplicate.
Назовем ее: **wing_1800****MH64_**up, и внесем в нее изменения как на рисунке.

Создадим для новой модели поляру для анализа: Polars -> Define Polar Analisis, с такими же параметрами как и для основной модели.
Запустим расчет поляры кнопкой – Analize.

На полученных графиках GCm от Alpha, видим, что модель теперь настроена стабильный полет при скорости 11 м/сек, приблизительно с углом атаки 6 град.

Но изменения затронули не только балансировку. На графике качества видно, что максимальное качество уменьшилось. Минимальная скорость снижения Vz почти не изменилась. Увеличилось сопротивление модели Drag на малых углах атаки.

Можно сделать вывод, что использовать отклонение элевонов для настройки ЛК на основную «рабочую» скорость нежелательно, такая настройка должна быть сделана при проектировании выбором профиля с положительным моментом и выбором необходимой крутки.

Для примера график с использованием профиля PW75. Без отклонения элевонов модель сбалансирована на скорость 12м/сек.

Часть 7. Винглеты.
Проведем с помощью XFLR5 анализ использования винглетов на стреловидном ЛК, чтобы определить какой вклад они вносят в сопротивление модели.
Размеры винглетов выбираем минимальными для получения максимального качества модели.

После добавления винглетов в Wing/Plane -> Current Wing/Plane -> Edit Fin ->Double Fin, произведем расчет поляры измененной модели.
На графике качества Cl/Cd от Vx видно небольшое увеличение качества на больших углах атаки и уменьшение качества на малых углах атаки.
Действие винглет, аналогично небольшому увеличению размаха и площади крыла. Догрузка небольшим количеством баланса, в нашем случае 70 грамм, и модель летит так же как и без винглет.

GreenGo
GreenGo:

Догрузка небольшим количеством баланса, в нашем случае 70 грамм

Догрузка конечно балластом.

Adekamer

Приветствую, а не подскажите как считать центровку вот такого крыла ?

GreenGo
Adekamer:

как считать центровку вот такого крыла ?

Нужны размеры для подстановки в WinLaengs4.
Померяйте корневую, концевую хорды, и смещение концевой (назад) от передней кромки (носа). Размах 1400мм плюс вставка 200мм?

Adekamer
GreenGo:

Размах 1400мм плюс вставка 200мм?

все верно

GreenGo
Adekamer:

все верно

Остальные размеры нужны.

Чем бы продолжить тему топика - аэродинамика ЛК?

boroda_de
GreenGo:

Чем бы продолжить тему

Можно углубиться в разбор таких животрепещущих тем, как напр. голландский шаг или перевести статьи с aerodesign.de

GreenGo
boroda_de:

напр. голландский шаг или перевести статьи с aerodesign.de

Голландский шаг был в первой статье, посмотрю, что там еще можно добавить.
На aerodesign.de нет английского варианта. Перевести не смогу:(

Вот такие темы будут в продолжении:

  1. управляющие плоскости и механизация,
  2. динамическая устойчивость (по тангажу и курсу),
  3. флатер,
  4. профили для ЛК.
GreenGo

На форуме уже довольно много тем посвящено флаттеру, достатточно набрать “флаттер” в поиске.
Поэтому очень кратко только для галочки:)

Флаттер – автоколебания возникающие при взаимодействии аэродинамических, упругих и инерционных сил, в результате несовпадения центра жесткости с центром давления и недостаточной жесткости конструкци крыла

  1. Изгибно-крутильный флаттер крыла, возникает когда центр тяжести сечений крыла находится позади центра жесткости.
  2. Изгибно-элеронный флаттер крыла, возникает когда несбалансированный элерон имеет недостаточную жесткость проводки тяг или люфты.

Для предотвращения флаттера модель должна иметь легкую и прочную (жесткую) конструкцию. В конструкции консолей должно быть предусмотрено расположение линии центров жесткости позади линии центров тяжести.
Проводка тяг управления должна быть достаточно жесткой и не иметь люфтов.
В качестве дополнительных мер, возможно использование противофлаттерных грузов и весовой компенсации элеронов.

kvs-vm.narod.ru/uchob/l7/7.html
aeroconstruction.ru/flatter/1
В теории все достаточно просто. 😃

На практике могут быть различные особенности в зависимости от конструкции. Особенно при разборной конструкции - конструктивный элемент обеспечивающий жесткость должен быть сплошным и узел стыковки не должен иметь люфта.

Полезные темы о флаттере на нашем форуме.
1.Флаттер
2.Флаттер
3.Флаттер

Прежде чем продолжить об управляющих поверхностях, вот интересное достижение сделанное в сентябре 2010 года:
Новости.

GreenGo
Adekamer:

снял размеры

Не хватило одного размера - 320мм? Если с фотки определил его правильно то ЦТ 170мм-8%, 180мм-4% от передней кромки.

Посмотрел внимательнее:) похоже смещение 340мм (а не 320), тогда ЦТ должен быть 177-189мм. Для начала я бы ставил 8% центровку, а потом подбирал по поведению модели.
А 220мм для модели без вставки это тогда похоже нейтральный пункт - центр давления.

GreenGo

Martin Hepperle.
Автор серии профилей MH_, программы для разработки профилей JavaFoil._
www.mh-aerotools.de/airfoils/index.htm

Профили для ЛК.

В принципе , модели бесхвосток и ЛК могут иметь любой профиль. Но если качество модели имеет значение, как в случае модели для соревнований, профиль должен быть тщательно подобран.

На ЛК возможно использование профилей для F3B моделей, если момент профиля компенсирован соответствующей комбинацией стреловидности и крутки. Однако излишняя стреловидность и крутка обычно уменьшают эффективность (качество) модели. Потери в эффективности могут быть минимизированы выбором профиля с коэффициентом Cm близким к нулю.

В начале 1980 годов, когда команда LOGO присоединилась ко мне в поиске новых профилей, нами были использованы различные профили на бесхвостых моделях планеров.
При хорошем качестве, модели имели проблемы с затяжкой леером и не очень хорошую управляемость. Это было вызвано малым Clmax и срывным характером профилей. Поскольку высота старта на леере очень важна, бесхвостые модели сильно отставали, из за этих двух характеристик использованных профилей.
Типичный асимметричный срыв (с одной консоли) приводил к полной потере управления с быстрым спиральным вращением модели на леере, которое заканчивалось с «катастрофическими» последствиями для модели.

Основными целями при создании серии профилей MH60 MH40 были:

  •      Малое сопротивление по сравнению с F3B профилями;  
    
  •      Малый коэффициент момента Cm;  
    
  •      Улучшение Clmax, по сравнению с другими профилями имеющими малый Cm.
    

Эти профили были опубликованы в 1988 году и широко использовались:
MH60 10%
MH61 10,28%
MH62 9,3%
MH64 8,61%

Следующие профили не планировались к публикации, но как то вышли за пределы моей лаборатории. Все они довольно хорошо работают, даже на Re ниже 200000.
MH44 9,66%
MH45 9.85%
MH46 11,39%
MH49 10,5%

Создание и выбор профиля.
Типы моделей ЛК и коэффициент Cm.

Вместе с коэффициентами подъемной силы и сопротивления, Cl и Cd, коэффициент момента Cm, так же имеет большое значение, особенно для ЛК, так как существенно влияет на продольную устойчивость модели.
ЛА обычной схемы могут компенсировать момент создаваемый крылом при помощи горизонтального оперения. На бесхвостых ЛА, по понятной причине, этого сделать нельзя.
Бесхвостые ЛА можно разделить на три группы в зависимости от способа получения продольной устойчивости.

  •      ЛК без стреловидности – прямые крылья. Продольная устойчивость достигается передним расположением ЦТ и балансировочным моментом профиля.  
    
    •      Стреловидные ЛК. На стреловидных ЛК возможно использовать любой профиль, т.к. стабильность может быть достигнута соответствующей комбинацией стреловидности и крутки. Но для получения хороших результатов подходят профили с малым Cm.  Такие профили требуют меньшей величины крутки, что увеличивает скоростной диапазон без дополнительных балансировочных потерь.  
      
    •      ЛК с низким расположением ЦТ (параплан). Cm менее важен, возможно использование традиционных профилей с отрицательным Cm. Обычно используют профили со средней величиной Cm т.к. они имеют меньшее сопротивление и  более широкий диапазон скоростей.
      

Коэффициент Cm и форма профиля.
Форму профиля можно представить как комбинацию формы средней линии и распределение толщины вдоль средней линии профиля. Требования для перечисленных трех классов ЛК, можно перевести в соответствующую форму профиля и крутку крыла.
Единственный способ получить положительный коэффициент момента Cm и требуемую величину подъемной силы – использовать S-образную форму средней линии профиля.

  •      ЛК без стреловидности. Требуемый положительный момент  Cm, приводит к необходимости наличия S- образности средней линии у подходящих, для данного класса ЛК, профилей.  
    
    •      Стреловидные ЛК. Малый Cm и малая требуемая крутка, могут быть получены при использовании профилей с малой кривизной (вогнутостью) средней линии и нейтральной или слегка S-образной формой.  
      
    •      Параплан. Нет строгого ограничения на форму профиля.
      

Профили с большой S-образностью обычно не используются на ЛА обычной схемы – они «уникальны» в использованиии только на бесхвостках.
С аэродинамической точки зрения S-образные профили сложны в дизайне и модификации, в частности, и потому, что они очень чувствительны по отношению к числам Re.

S- образность и коэффициент Cm.
Мы уже знаем, что Cm и форма средней линии профиля тесно связаны. Если мы посмотрим на S- образную среднюю линию более внимательно, мы увидим, что форма задней части средней линии очень сильно влияет на коэффициент Cm.
Фактически - отклоняя заднюю кромку профиля, мы можем получить практически любой желаемый коэффициент Cm.
Рис. Внизу показывает, как можно управлять значением Cm: используя отклонение задней части профиля. Задняя часть профиля отклоняется вверх на 5 и 10 градусов и в результате Cm следует за этим изменением.

С помощью этого трюка, проблема кажется решенной. Мы просто отгибаем заднюю кромку профиля, пока не получим нужный балансировочный момент у нашего ЛК и все…
Но мы обычно предпочитаем ЛА, который не только сбалансирован (настроен на нужную скорость полета), но так же имеет высокое качество и малое сопротивление, и вот здесь то и начинаются проблемы.

S-образность и коэффициенты подъемной силы и сопротивления Cl и Cd.
Форма и расположение кривой Cl от Cd, это ключ к эффективности (качеству) ЛА.
Рис. Показывает, как Cl – Cd кривая изменяется в зависимости от формы средней линии профиля (степени S-образности).

Теперь у нас появляется новая проблема: с добавлением S-образности средней линии профиля для получения положительного коэффициента Cm, мы сдвигаем Cl – Cd «поляру» вниз. Это означает, что мы уменьшаем подъемную силу и что еще хуже мы уменьшаем Clmax. А это ведет к увеличению скорости срыва, что нам совершенно не нужно.
Конечно, профессионалы в аэродинамике знают, что делать, в данном случае, при создании профиля. Необходимо увеличение кривизны средней линии. Это увеличит подъемную силу и немного уменьшит коэффициент момента Cm.

Положение точки максимальной кривизны средней линии и коэффициент Cm.
К счастью, у нас еще есть один параметр профиля для компенсации дестабилизирующего эффекта увеличенной кривизны средней линии – место расположения максимальной кривизны.
Место расположения максимальной кривизны средней линии имеет небольшое влияние на Cl –Cd, но очень сильное влияние на коэффициент Cm.
На Рис. Показано семейство профилей с разным положением максимальной кривизны среденей линии Xc/c и графики коэффициента Cm.

Как мы видим, сдвиг Xc/c назад. Так же уменьшает Cm (сдвигает график вниз), к отрицательным значениям.
Поэтому разумным будет иметь координату Xc/c в первой четверти хорды 0-25%, если мы хотим компенсировать потерю подъемной силы от использования S-образности средней линии.

Сложность проектирования – выбора профиля с малым коэффициентом Cm.
Мы уже узнали, что несколько параметров, профиля связаны вместе, и влияют на дизайн-выбор профиля с малым коэффициентом Cm и большим коэффициентом подъемной силы.
Но мы еще не поговорили о дополнительных проблемах, связанных с поведением пограничного слоя.

Рис. И список показывают связь наиболее важных параметров.
- Смещение максимальной кривизны средней линии назад.
Плюс - более положительный Cm
Минус- большее влияние на пограничный слой у носка профиля – пик давления.
- Увеличение кривизны (вогнутости) средней линии.
Плюс- большее Cl (график Cl – Cd сдвигается в сторону положительных значений.
Плюс – Clmax больше.
Минус- более отрицательный коэффициент Cm.
Минус- большее влияние на пограничный слой.
- Увеличение S-образности.
Плюс- более положительный коэффициент Cm.
Минус- меньше подъемной силы.
Минус- Clmax уменьшается.

Распределение скорости обтекания профиля.
Чтобы «пролить» больше света на проблему влияния S-образности профиля на пограничный слой, посмотрим на распределение скорости потока обтекании профиля.
Рис. Показывает распределение четырех разных профилей с разными комбинациями S-образности и вогнутости средней линии профиля.

В общем распределение по верхней поверхности профиля имеет максимальную скорость в передней трети профиля, и постепенное уменьшение к задней кромке.

В зависимости от Re, слишком большой наклон линии скорости (которая по уравнению Бернулли показывает подъем давления) и слишком быстрый подъем давления не желателен. Когда подъем давления слишком резок, поток отделяется, формируя «ламинарный пузырь», что резко ухудшает эффективность профиля.

Распределение по нижней поверхности создает меньше проблем, с отделением потока вблизи задней кромки профиля.
Уникальной особенностью профилей с S-образной средней линией является пересечение графика распределения скорости на второй половине хорды. Увеличение рефлекса – S-образности, увеличивает локальную скорость потока на нижней поверхности и уменьшает скорость на верхней поверхности. Результатом наличия замкнутой области на графике распределения скорости является положительный коэффициент момента Cm.

Увеличение вогнутости в передней части профиля увеличивает локальную скорость на верхней поверхности и уменьшает скорость на нижней поверхности. Слишком большая кривизна средней линии в передней части профиля может приводить к ухудшению эффективности профиля, т.к. должна компенсироваться большим отклонением задней кромки вверх. Все вместе это создает сильный стресс для пограничного слоя и приводит к раннему срыву, что может быть опасно на взлете и посадке.

Выводы
В зависимости от типа басхвостого ЛА, обеспечение устойчивости ведет к различным критериям выбора профиля. Для большинства бесхвостых ЛА, профиль с малым коэффициентом Cm позволяет достичь лучших результатов.
Малый коэффициент Cm и большой коэффициент подъемной силы Cl могут быть получены использованием профиля с S- образностью, но получаемое распределение скорости (давления)на малых Re может в результате вызывать проблемы со срывным характером профиля.

Лучшим компромиссом может быть профиль со средней величиной S-образности, в комбинации с вогнутостью сдвинутой в переднюю часть профиля и тупым носком профиля.
Относительно срывных характеристик, профиль играет важную роль, но только в комплексе с распределением подъемной силы по размаху, формой в плане и распределением крутки.

Стреловидные ЛК с сужением имеют распределение с сдвигом к концам крыла, что создает склонность к концевому срыву, если не используется дополнительная крутка крыла.

GreenGo

Martin Hepperle.
Автор серии профилей MH, программы для разработки профилей JavaFoil.
www.mh-aerotools.de/airfoils/index.htm

Винглеты
Большинство современных высокоэффективных ЛК, с близким к эллиптическому распределением подъемной силы имеют винглеты.
ЛК с колоколообразным распределением Хортена, обычно не нуждаются в винглетах.
Винглеты иногда используют на ЛА обычной схемы т.к. позволяют уменьшить индуктивное сопротивление, что аналогично увеличению размаха. В сравнении с увеличением размаха, винглеты создают меньший изгибающий момент в лонжероне, что делает их полезными для улучшения существующих ЛА («увеличение размаха» без необходимости увеличения прочности крыла).
Цель этой статьи не в поиске оптимальной формы винглетов и ее влиянии на индуктивное сопротивление, а в том, чтобы понять, что происходит в области где винглет соединяется с крылом и почему уменьшение индуктивного сопротивления, при неудачном варианте соединения, может быть ухудшено дополнительным сопротивлением в этой области. Поэтому мы сконцентрируемся на эффектах пограничного слоя.
В журнале «Competitor Achmer News» 1992года, я нашел детальный эскиз трех ЛК с винглетами. На этом эскизе были показаны углы образованные соединением крыла и винглета. Этот эскиз подтолкнул меня к исследованию этой области аэродинамики ЛК.

Объект исследования

Для исследования было выбрано, простое ЛК, в конфигурациях с:

  • закругленной законцовкой без винглета;
  • винглетом с закругленным плавным соединением;
  • винглет с соединением углом;
  • винглет сдвинутый назад на половину хорды крыла.

Слишком большой стресс.

Пограничный слой очень чувствителен к подъему давления, которое происходит позади точки максимальной толщины профиля. Если давление нарастает слишком быстро (что соответствует быстрому спаду скорости потока), возможен ранний отрыв пограничного слоя с образованием «пузыря». Ранний отрыв потока создает большое дополнительное сопротивление.
При сравнении потока обтекания при отсутствии и при наличии винглета, можно увидеть большое различие в изменении распределения давления.

Там где свободное крыло «чувствует» подъем давления в одном направлении, при наличии винглета подъем давления происходит в двух направлениях одновременно и суммируясь производит увеличенный стресс на пограничный слой в месте соединения крыла и винглета.

Это подтверждается графиком распределения скорости потока в интересующей нас области. Так же можно заметить, что влияние формы законцовок распространяется на небольшую часть размаха (около 10% размаха).

Неровная дорога.

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

Начнем с черной линии, она относится к крылу без винглетов. Мы видим, что остальные графики проходят выше – винглеты увеличивают скорость потока в месте соединения, но по разному.

  • Винглет с загругленным соединением (красный цвет) сдвигает скорость потока вверх, равномерно без искажений формы графика, и не создает большого стресса для пограничного слоя Единственный недостаток в этом случае – сложность при строительстве модели.
    • Винглет с угловым соединением (зеленый цвет) создает искаженную форму распределения с пиком на носке профиля. Далее почти до 60% хорды идет плоское рапределение, которое предпочтительно для пограничного слоя, но затем следует быстрое уменьшение скорости, что приводит к отделению пограничного слоя.
    • Винглет со сдвигом назад (синий цвет) имеет распределение почти не вносящее изменений в сравнении с крылом без винглетов, до 50% хорды. Далее идет область до 80% с более плоским и менее подверженным отрыву распределением. Общий результат будет аналогичен соединению с закругленным плавным переходом крыла к винглету.

      Индуктивное сопротивление

Основная идея использования винглетов на ЛА обычной схемы, состоит в уменьшении индуктивного сопротивления, поэтому вопрос был бы не полно раскрыт без изучения индуктивного сопротивления винглетов.
С другой стороны использование винглетов на стреловидном ЛК заключается в получении устойчивости по курсу и уменьшение индуктивного сопротивления - это небольшой дополнительный бонус.

Поляры показывают, что все три варианта с винглетами дают очень близкий результат в уменьшении индуктивного сопротивления. Есть небольшое преимущество последнего варианта (со сдвигом), но оно почти в пределах погрешности расчетов (очень небольшое).
Все конфигурации с винглетами превосходят крыло без винглетов, и уменьшение индуктивного сопротивления становится больше на больших Cl (с ростом угла атаки).

Этот график, не включает профильного сопротивления, которое уменьшает выигрыш от винглетов, из за увеличенной площади поверхности.
С учетом профильного сопротивления, на малых углах атаки преимущество будет у крыла без винглетов, и наоборот на больших углах атаки преимущество у крыла с винглетами.
Кроме того наличие винглетов делает центральный киль (который так же создает дополнительное сопротивление) ненужным.

Большой брат

Если мы посмотрим на решения используемые в большой авиации, то увидим, что большинство используемых винглетов устанавливается со сдвигом назад и при возможности с плавным переходом в месте соединения с крылом.