Введение
@@В доступных статьях по авиамоделизму можно найти расчет крыла, и много всяческой полезной информации по компоновкам и постройке моделей. Всем ясно, что в лобовое сопротивление модели вносят вклад все ее детали. Но в доступной литературе практически нет информации по этой части лобового сопротивления. Возникает вопрос: почему? Причина простая: даже лобовое сопротивление (и его остальные параметры) крыла авиамодели можно только оценить. Причины: нельзя с достаточной степенью точности изготовить профиль крыла; нет (или почти нет) достоверных данных о характеристиках профилей на модельных скоростях и размерах. Вот и получается, что учитывать вклад "мелочей" вроде бы бессмысленно. На мой взгляд, считать вклад мелочей действительно бесполезно, а вот его оценка может быть полезна. Польза заключается в том, что позволит не "полировать" или усложнять что-то до бесконечности с одной стороны, и не допускать ощутимой потери летных свойств из-за "мелочей", которые можно легко сделать лучше.
Как учитываются аэродинамические "мелочи"
@@В книгах по аэродинамике предлагается следующий подход.
@@Аэродинамическое качество самолета равно Cy/Cx. Где Cy - коэффициент подъемной силы крыла, а Cx - коэффициент лобового сопротивления самолета.
@@То есть, Cx отражает сумму сопротивлений всех элементов деталей самолета. Но подъемная сила крыла равна S*P*Cy*V*V/2, где S - площадь крыла. Очевидно, что Р (плотность воздуха) и V (скорость) для всего самолета одна и та же (если "забыть", что некоторые детали самолета могут обдуваться струей винта).
@@А вот площади у всех деталей самолета разные. Так что, если мы определим Сх для какой-то детали (а именно это есть возможность сделать), то прямо в сумму Сх самолета мы не можем ее подставить. Выход простой: полученный Сх детали нужно умножить на ее характерную площадь и поделить на площадь крыла. Пример:
@@Сх крыла Схк=0.017; Сх стабилизатора и киля (которые рассматривались как крыло) Схс=0.012; Сx фюзеляжа=0.0037; площадь крыла Sк=30дм**2 ("**" у меня будет обозначать степень), площадь стабилизатора и киля Sс=5дм**2. Поверхность фюзеляжа (рассматривается планер с фюзеляжем обтекаемой формы) Sф=10дм**2. Определим Cх планера в целом Cxпл:
@@Схпл=Схк+Схс*Sс/Sк+Схф*Sф/Sк=0.017+0.012*5/30+0.0037*10/30=0.0202
@@А дальше мы можем уже анализировать результат и принимать (или не принимать) какие-то меры.
@@Допустим, на этот планер мы решили поставить видеокамеру в виде кубика 5*5*5см. считаем, что "кубик" поставлен удачно и дополнительного сопротивления от влияния на другие детали планера у него нет. Попробуем посчитать. Cx кубика = 0.8; Sкубика=0.25дм**2 (берется площадь "лба" куба). Посчитаем добавку от него в коэффициент лобового сопротивления планера =0.8*0.25/30=0.0067. Легко заметить, что видеокамера в виде куба испортит аэродинамическое качество модели существенно сильнее, чем его "испортило" наличие фюзеляжа и хвостового оперения. Придется для камеры делать бокс обтекаемой формы.
Пограничный слой
@@Воздух имеет вязкость, и имеются силы трения о поверхность.
@@В аэродинамике есть "магическое число" - число Рейнольца, обычно пишется Re. Оно зависит от вязкости и плотности воздуха, размера тела (для которого считается Re), и его скорости. В авиамоделизме, чтобы не усложнять себе жизнь, обычно Re вычисляется очень просто: Re=69*X*V, где X - размер вдоль потока (для крыла - хорда) в миллиметрах, V - Скорость в метрах в секунду.
@@Перейдем теперь непосредственно к пограничному слою. Им называется слой воздуха вдоль поверхности, который из-за сил трения движется медленнее остального потока воздуха. Выглядит он приблизительно следующим образом: непосредственно вблизи поверхности воздух неподвижен, чуть дальше от нее начинает относительно поверхности ползти потихонечку, еще выше слой "ползет" по этому слою и так далее. На каком-то расстоянии от поверхности воздух движется уже со своей "нормальной" скоростью. Вот такой пограничный слой (когда все скользит по слоям) называется ламинарным. Если длина поверхности больше какой-то, и скорость воздуха больше какой-то, (попросту - Re больше некоторой величины), то вместо скольжения внутри пограничного слоя получаются маленькие вихри. Такой пограничный слой называется турбулентным. Для обоих видов пограничного слоя есть формулы для определения его толщины и силы трения.
@@Толщина слоя для ламинарного слоя: d=0.7*((X/V)**0.5), где d - толщина пограничного слоя в миллиметрах, X - расстояние от передней кромки в миллиметрах, V - скорость в метрах в секунду. "**0.5" читается как корень квадратный из отношения X к V.
@@Аналогично, толщина турбулентного пограничного слоя d=0.16*(X**0.8)/(V**0.2), все размерности величин те же.
@@Коэффициенты трения:
@@Для ламинарного слоя Cx=1.33/((Re)**0.5)
@@Для турбулентного слоя Cx=0.074/((Re)**0.2)
@@Возникает вопрос: "А зачем нужно знать толщину пограничного слоя?". Ответ довольно простой - если Вы подставите какие-нибудь конкретные числа для своей модели (скорость, хорду крыла), то увидите, что крыло находится в довольно толстой "шубе" из пограничного слоя (если Вы подставите эти же данные для какого-нибудь настоящего самолета или планера, то увидите, что, относительно их размеров, эта "шуба" намного тоньше).
@@Внимание! Эти формулы дают полную толщину пограничного слоя, а реальная толщина, в которой воздух движется существенно медленнее, меньше толщины "шубы" раз в десять.
@@Следствие. Чем ближе в передней кромке (или к носу модели), тем тщательнее надо обрабатывать поверхность и соблюдать форму, так как там "шуба" тоньше.
@@Закономерность. Коэффициент трения при турбулентном пограничном слое выше, чем при ламинарном (при прочих равных условиях).
@@Внимание! Если у Вас получился коэффициент трения турбулентного слоя МЕНЬШЕ, чем ламинарного, то пользуйтесь ламинарным - на этих Re формула для турбулентного слоя просто не работает, так как она выведена для "большой" авиации.
@@Немного о точке перехода ламинарного пограничного слоя в турбулентный. Для модельных Re в случае плоской гладкой пластинки все обтекание было бы ламинарным.
@@Но.
@@Шероховатость поверхности приводит к ранней турбулизации пограничного слоя.
@@Наличие неровностей и уступов приводит к тому же.
@@Все, что обтекается струей винта, имеет турбулентный пограничный слой.
"Мелочи" и не только они
@@Все дальнейшие данные я взял из учебников для ВУЗов. Написаны они для "больших самолетов, так что, применительно к моделям, могут врать. Так как я не умею рисовать, то буду давать словесное описание "изделия" и, вместо графиков будут таблицы.
Рисунок 1
@@"Изделие 1".
@@Диск или пластина (линейка) поперек потока . То есть она поставлена так, чтобы создавать наибольшее сопротивление (ветер перпендикулярен плоскости).
@@В качестве "площади для расчета" сопротивления модели берется площадь пластинки.
@@Коэффициент сопротивления от Re не зависит. Для круглой пластинки Сх=1.12, причем, избыток давления спереди дает Сх=0.82, разрежение за пластинкой дает Сх=0.3 (в сумме получается Сх=1.12).
@@Для "линейки" коэффициент сопротивления увеличивается с ростом отношения длины к ширине.
Таблица для "линейки"
Отношение длины к ширине
|
Сх
|
Примечание
|
1
|
1,12
|
Квадрат, круг
|
3
|
1,15
|
.
|
5
|
1,18
|
.
|
10
|
1,25
|
.
|
Очень много
|
2
|
.
|
Квадратную пластинку (воздушный змей) постепенно поворачивают
Угол, градусы
|
Сх
|
Примечание
|
0
|
0
|
Пластинка вдоль потока, реальное сопротивление рассчитывается из сил трения. Не забывайте, что у пластинки две стороны!
|
5
|
0,2
|
..
|
10
|
0,35
|
..
|
20
|
0,6
|
..
|
30
|
0,8
|
.
|
45
|
0,9
|
.
|
60
|
1
|
.
|
90
|
1,12
|
Пластинка перпендикулярна ветру
|
@@Пластинку превращают в кирпич. График дан для цилиндра, вдоль оси которого дует ветер (фюзеляж типа "полено", нос и хвост обрублен).
Таблица для "полена". За размер площади принята площадь торца полена (круга)
Отношение длины цилиндра к его диаметру
|
Сх
|
Примечание
|
0
|
1,12
|
диск
|
1
|
0,8
|
"кубик"
|
2
|
0,73
|
..
|
3
|
0,73
|
..
|
5
|
0,8
|
.
|
8
|
0,9
|
длинное "полено"
|
@@"Изделие 2".
@@Цилиндр поперек потока (круглый подкос крыла, стойка шасси из проволоки). Re считается исходя из диаметра. Площадь - произведение диаметра на длину.
Таблица для очень длинной проволоки
Re
|
Сх
|
Примечание
|
10
|
2,7
|
Диаметр 0.14мм, скорость 1м/с
|
30
|
1,8
|
.
|
100
|
1,35
|
D=0.7мм, V=2м/с
|
300
|
1,2
|
.
|
1000
|
1
|
D=1.5мм, V=10м/с
|
3000
|
0,9
|
.
|
10000
|
1,1
|
D=15мм, V=10м/с
|
30000
|
1,2
|
.
|
100000
|
1,2
|
.
|
@@При уменьшении отношения длины цилиндра к его ширине сопротивление слегка уменьшается.
@@На этот коэффициент надо умножить полученный выше, если у Вас короткий "чурбачок" (если он торчит из чего-то, надо взять его удвоенную длину для определения коэффициента, за площадь, естественно, берется произведение его диаметра на длину).
Таблица для коротких "чурбачков", поставленных поперек потока
Отношение длины цилиндра к диаметру
|
Поправочный коэффициент
|
1
|
0,4
|
2
|
0,5
|
5
|
0,6
|
10
|
0,8
|
20
|
0,9
|
@@"Изделие 3".
@@Стойка шасси или подкос каплевидной формы. Отношение длины капельки к толщине капельки 3-4. Площадь считается умножением толщины капельки на длину подкоса или стойки шасси.
@@К сожалению, минимальное Re, для которого есть информация - 50000. Для меньших Re лучше применить эллипс (тупой носик и хвостик) или "чечевицу" (острый носик и хвостик)
Re
|
Cx
|
200000
|
0,08
|
150000
|
0,09
|
100000
|
0,1
|
75000
|
0,12
|
50000
|
0,15
|
@@@"Изделие 4".
@@Шар. На числах небольших Re (<200000) поток отрывается, образуется много вихрей, в результате Cx=0.5 На больших числах Re обтекание шара становится турбулентным (турбулентный поток лучше "прилипает" к кривым поверхностям), вихрей становится меньше и Сх=0.12. Поэтому в "большой авиации" раньше были популярны носы фюзеляжа в форме полусферы (большие Re).
@@В моделях "полусфера" будет давать большое сопротивление, поэтому носы надо "вытягивать". Число Re рассчитывается по диаметру шара.
@@"Изделие 5".
@@Удобообтекаемые тела (фюзеляж планера).
@@Считается трение. Если фюзеляж действительно удобообтекаемый, то поверхность до самого толстого места считается с ламинарным пограничным слоем, а после него - с турбулентным. Я бы брал просто полусумму этих двух Сх. Трение хвостовой балки вообще бы не учитывал (Re считается из той длины фюзеляжа, которая участвует в расчете площади поверхности).
Общие замечания
@@Стабилизатор типа "плоская пластина". Он довольно распространен. Часто бывает "вырубленным" без закругления кромок (передней, заднюю закруглять нет смысла, эффект уменьшения сопротивления даст только сведение ее на нет клином с отношением длины к толщине не меньше 5).
"Обрубленная" передняя кромка дает два эффекта:
-
тот самый коэффициент Сх=0.8 для линейки (то есть, надо посчитать площадь передней кромки и учесть с коэффициентом Сх=0.8; впрочем, я считаю, что в этом случае коэффициент надо в пару раз уменьшить);
-
весь стабилизатор имеет турбулентный пограничный слой; считается, как плоская пластинка (не забудьте учесть, что у стабилизатора две стороны);
@@Обрубленная задняя кромка дает по аналогии с передней кромкой сопротивление с Сх=0.3 (я считаю, что его тоже надо уменьшить в пару раз). Это же касается и "обрубленной" задней кромки крыла (см. Рис. 2).
@@Закругления передних кромок. Популярны просто полукруглые. На мой взгляд, надо закруглять по эллипсу с отношением длины к толщине 2/1 (носик "острее" полукруглого).
@@В случае хорошо закругленного носика стабилизатор имеет ламинарный пограничный слой.
@@Пример оценки дополнительного сопротивления:
Рисунок 2
@@О моделях - полукопиях "лобастых" самолетов (И16, Зеро, Ла5, Як18). У них середина носа имеет вид диска (реально там воздухозаборник системы охлаждения), а края закруглены. Если у модели нет действующего воздухозаборника, то надо считать, что 2/3 "диска" имеют Сх=0.8 (см. Рис. 2).
О здравом смысле
@@Когда Вы начнете оценивать "мелочи", увидите довольно много неожиданного. Я в свое время "неожиданностей" увидел много. Не забывайте, что все это - только оценки (хотя и считается точно, но сами исходные формулы выведены для "большой" авиации). И совсем немного о "мелочах" которые "мелочами" не являются.
@@Первая "мелочь" касается верхней поверхности крыла и всего, с чего воздух может на нее попасть (учтите при этом, что перед крылом встречный поток движется под углом вверх).
@@Из воспоминаний. Есть такой планер - "БРОшка". Очень простенький, на нем школьников учили летать с лебедки. У этого планера одно колесо. Первое упражнение - "пробежки" по земле без взлета. Тянет лебедка, а задача пилота - не дать планеру завалиться на крыло. Управление: руль направления (педали) и элероны (ручка). С одной стороны, планер должен двигаться быстро, чтобы рули работали, а с другой - ему нельзя взлетать (при падении даже с высоты двух метров можно основательно покалечиться). Задачу решили просто - на середину крыла, сверху, (хорда около одного метра, размах - около десяти), приблизительно на половину размаха, были поперек привязаны деревянные бруски сечением 5см*5см. В результате этот планер принципиально не мог взлететь на скоростях больше взлетной!!!
@@Так что - "вылизывайте" все, с чего воздух может попасть на верхнюю поверхность крыла! Естественно, это касается и "мелочей", которые Вы хотите на верхнюю поверхность крыла пристроить. Можно не только схлопотать лишнее сопротивление (причем, в разы больше ожидаемого), но и "потерять" кусок крыла в создании подъемной силы.
@@О подвеске рулей. Если Вы их подвешиваете на шарнирах, то будьте очень внимательны к размеру щелей и смещений по вертикали подвешенного руля относительно плоскости, к которой он подвешен. Конечно, в "большой" авиации есть щелевые закрылки, хитро (со щелями) подвешенные рули.
НО ТАМ ПРИМЕНЯЮТСЯ СПЕЦИАЛЬНО СПРОЕКТИРОВАННЫЕ ЩЕЛИ И СМЕЩЕНИЯ!!!
@@"Неправильная" щель или смещение может "убить" эффективность руля.
@@Тут тоже у меня есть пример из жизни. Приятель сделал планер (довольно простенький). И подвесил руль высоты на шарнирах так, что при отклонении вверх руль изгибался, возникала приличная щель и смещение между ним и стабилизатором. В результате модель почти не реагировала на руль высоты. После подвески руля просто полоской скотча руль оказался даже избыточно эффективным.
Заключение
@@Здесь я хочу еще раз напомнить об очевидной вещи: авиамодель нельзя посчитать, ее можно только оценить. Если, конечно, у Вас нет доступа ко всяческим "закрытым" программам, аэродинамическим трубам и оборудованию, позволяющему изготовить многие детали (например, профиль крыла) с точностью в сотую миллиметра.
@@В "открытом доступе" есть довольно много расчетных программ, которые, как кажется на первый взгляд, могут довольно точно посчитать модель. К сожалению, это только кажется.
@@Например, при сравнеии поляр профиля на небольших Re (60000), вычисленных неплохой программой, с реально измеренными полярами, взятыми с сайта, получилось большое расхождение (в некоторых точках поляры расхождение по Cx получилось разница в пару раз).
@@Тем не менее, глупо не проводить ОЦЕНКУ, если есть такая возможность. Производить оценку вручную утомительно (и глупо, когда под руками есть PC). Пока у меня на РС был DOS, я в качестве "калькулятора" использовал Турбо Бейсик. Написав за 5минут - 1 час простенькую программку, я мог потом легко повторять расчеты, варьируя исходные данные, при желании мог получать графики нужных мне в данный момент зависимостей. К сожалению, языками программирования под Windows я не владею - просто времени не хватает освоить их. "Калькулятор" я сразу отмел - достаточно ошибиться в наборе одной цифры - и можно часами разбираться в полученных результатах. Стал искать замену Бэйсику. И нашел - в состав офиса входит Exel, который удалось в достаточной степени освоить за несколько часов. Существенный плюс Exel (и программы, написанной на каком-нибудь языке) в том, что видны ОДНОВРЕМЕННО и исходные данные и результат (не говоря уже о том, что изменение входных данных автоматически изменяет результат). Это позволяет легко найти ошибку типа "опечатка". Для примера привожу два простеньких "калькулятора".
@@Каждый из них был сделан менее, чем за час (не считая "раскраски", которую я обычно не делаю).
Скачать первый калькулятор (zip-архив, 7 кБ)
Скачать второй калькулятор (zip-архив, 5 кБ)
@@Как что считается, Вы можете сами посмотреть, просто щелкнув мышкой по любой желтой ячейке (только ничего не вводите в нее, так как там окажется то, что Вы введете). RyCx прямо указывает на координаты ячейки (Ry - номер строки, Cx - номер столбца) откуда берутся данные; R[y]C[x] - координаты относительно текущей ячейки.
Авторам опубликованных статей предоставляются скидки в нашем магазине