Содержание
Введение
1. Подготовка исходных данных
2. Расчёт и построение зависимостей cya(α) для различных режимов полёта
2.1 Расчёт и построение зависимости критического числа Маха от коэффициента подъёмной силы Мкр(суа)
2.2 Расчёт и построение вспомогательной зависимости cya(α)
2.3 Расчёт и построение взлётных кривых cya(α)
2.4 Расчёт и построение посадочных кривых cya(α)
2.5 Расчёт и построение крейсерских зависимостей cya(α)
3. Расчёт и построение поляр самолёта
3.1 Расчёт и построение вспомогательной поляры
3.2 Расчёт и построение взлётных поляр
3.3 Расчёт и построение посадочных поляр
3.4 Расчёт и построение крейсерских поляр
Библиографический список
Введение
В данной работе рассматривается лёгкий спортивный самолет «T-30 Katana», представляющий собой одноместный одномоторный свободнонесущий среднеплан с закрытой кабиной и неубирающимся шасси. Вычисляются основные геометрические и аэродинамические параметры этого самолёта, на основании которых строятся теоретические зависимости коэффициента подъёмной силы от угла атаки и от коэффициента сопротивления (поляры) для взлетного, крейсерского и посадочного режимов полёта.
1. Подготовка исходных данных
Аэродинамические характеристики самолёта зависят от его геометрических параметров. Поэтому сначала по чертежу летательного аппарата (рисунок 1) узнаём необходимые размеры из заносим их в таблицу 1. В эту же таблицу заносим основные лётно-технические характеристики самолёта. Далее на основании имеющихся данных вычисляются прочие необходимые геометрические характеристики и также включаются в таблицу вместе с формулами, по которым они были вычислены.
Таблица 1.
Элемент самолета, параметр
Размерность
Обозначение, формула
Значение
1
2
3
4
1. Крыло:
1.1 Размах/ размах его консолей
мм
l / lk= l — Dф
7,70/6,77
1.2 Площадь
м2
S
10,60
1.3 Хорда средняя
мм
B= S/ l
1,38
1.4 Хорда центральная
мм
b
1,82
1.5 Хорда концевая
мм
bк
0,89
1.6 Сужение в плане
ηb= b0 / bк
2,04
1.7 Относительная толщина профиля центрального
/>
0,20
1.8 Относительная толщина профиля концевого
/>
0,12
1.9 Средняя относительная толщина профиля
/>= (/>∙ ηb+ />) / (ηb+ 1)
0,17
1.10 Относительная координата максимальной толщины
/>= />/ b
0,23
1.11 Стреловидность по линии
max-х толщин
град.
/>
-1--PAGE_BREAK--
--PAGE_BREAK----PAGE_BREAK----PAGE_BREAK--
10. Общие данные
10.1 Взлётная масса самолёта
кг
m
880
10.2 Расчетная скорость полета
км/ч
V
365
10.3 Расчетная высота полета
км
H
2,5
10.4 Тип и количество двигателей
n
1 проп. дв.
10.5 Стартовая тяга (мощность) одного двигателя при V=0, H=0
даН
(кВт)
Р0i
(N0i)
220
(300 )
10.6 Среднее за полет аэродина- мическое качество рассматриваемого самолета
К
~12,8
10.7 Относительная масса топлива
/>т= mт/ m
0,2
2. Расчёт и построение зависимостей cya(α)для различных режимов полёта
2.1 Расчёт и построение зависимости критического числа Маха от коэффициента подъёмной силы Мкр(суа)
Эта зависимость задаётся формулой:
/>
/>
Рисунок 2 — Зависимость критического числа Маха от режима полёта
2.2 Расчёт и построение вспомогательной зависимости суа(α)
Эта зависимость строится для полёта на нулевой высоте при отсутствии экранного эффекта с убранными средствами механизации крыла на минимальной скорости полёта, которая находится по следующей формуле:
/>
Ей соответствует число Маха:
/>
Удлинение крыла данного самолёта достаточно велико (λ>4),и поэтому для нахождения теоретического наибольшего значения коэффициента подъёмной силы можно применить формулу:
/>
/>
/>
/>
/>
/>
/>
Определяем три точки для построения графика суа(α):
/>
/>
/>
И строим по этим трём точкам график зависимости коэффициента подъёмной силы от угла атаки, аппроксимируя её параболой в области больших углов атаки (рисунок 3 и рисунок 4, кривая 1).
/>
Рисунок 3 — Вспомогательная зависимость коэффициента подъёмной силы от угла атаки.
2.3 Расчёт и построение взлётных кривых суа(α)
1) Во взлётном режиме закрылки выпущены под углом:
/>
Зная данную величину, а также относительную хорду закрылков />, по справочным данным определяем приращение угла атаки нулевой подъёмной силы из-за выпущенных закрылков:
/>
То есть на взлёте этот угол равняется:
/>
Теперь можно найти приращение максимального значения коэффициента подъёмной силы из-за выпущенных закрылков по формуле:
/> продолжение
--PAGE_BREAK--
где /> — величина определяемая типом механизации крыла. Данный самолёт оснащён простым безщелевым отклоняемым закрылком, для которого />. Для учёта влияния обдувки крыла винтом на подъёмную силу найдём сначала коэффициент нагрузки винта по тяге во взлётном режиме:
/>
Зная данную величину, а также относительную площадь крыла, обдуваемую винтом, />, по справочным данным определяем изменение максимального значения коэффициента подъёмной силы за счёт обдувки крыла винтом: />.
Теперь можно вычислить максимальное значение коэффициента подъёмной силы во взлётной конфигурации без учёта экранного эффекта земли:
/>
Исходя из найденных значений />и неизменившейся величины />, строим по аналогии со вспомогательной зависимостью суа(α)из пункта 2.2 взлётную кривую суа(α)без учёта экрана земли (рисунок 4, кривая 2).
2) Учтём теперь влияние экрана земли, которое вносит такую поправку в коэффициент подъёмной силы:
/>
А максимальное значение коэффициента подъёмной силы во взлётной конфигурации с учётом экранного эффекта земли равно:
/>
Находим фиктивное удлинение крыла, учитывающее влияние земли:
/>
Тогда производная коэффициента подъёмной силы по углу атаки с учётом экранного эффекта равна:
/>
Используя найденные значения />, строим взлётную кривую суа(α)с учётом влияния земли (рисунок 4, кривая 3).
2.4 Расчёт и построение посадочных кривых суа(α)
1) Во время посадки закрылки выпущены под углом:
/>
Зная данную величину, а также относительную хорду закрылков />, по справочным данным определяем приращение угла атаки нулевой подъёмной силы из-за выпущенных закрылков:
/>
То есть на посадке этот угол равняется:
/>
Теперь можно найти приращение максимального значения коэффициента подъёмной силы из-за выпущенных закрылков по формуле:
/>.
На посадке винт не влияет на подъёмную силу крыла. Тогда максимальное значение коэффициента подъёмной силы во взлётной конфигурации без учёта экранного эффекта земли равно:
/>
Исходя из найденных значений />и неизменившейся величины />, строим посадочную кривую суа(α)без учёта экрана земли (рисунок 4, кривая 4).2) Учтём теперь влияние экрана земли, которое вносит такую поправку в коэффициент подъёмной силы:
/>
А максимальное значение коэффициента подъёмной силы в посадочной конфигурации с учётом экранного эффекта земли равно:
/>
Используя найденные значения />, строим посадочную кривую суа(α)с учётом влияния земли (рисунок 4, кривая 5).
2.5 Расчёт и построение крейсерских зависимостей cya(α)
Высота полёта расчётная Н=2500 м, скорость звука на этой высоте равна аН=330,6 м/с. Скорость полёта расчётная V=101,4 м/с, при этом число Маха равно:
Мрасч=V/аН=0,31.
Самолёт находится в полётной конфигурации, то есть закрылки убраны. При этом для различных чисел Маха зависимость cya(α)задаётся формулой:
/>
Сводим в таблицу 2 параметры этой зависимости для нескольких чисел Маха.
Таблица 2.
М
0,00
Мрасч=0,31
0,40
0,50
0,60
/>
0,078
0,082 продолжение
--PAGE_BREAK--
0,085
0,090
0,097
/>
0,603
0,634
0,658
0,696
0,754
И по этим данным строим крейсерские зависимости cya(α)(рисунок 5).
3. Расчёт и построение поляр самолёта
3.1 Расчёт и построение вспомогательной поляры
1) При построении данной поляры принимают, что закрылки убраны, высота полёта нулевая, экранный эффект отсутствует, скорость полёта минимальна (М=Мmin).
2) Для нахождения профильного сопротивления фюзеляжа, сначала вычислим его число Рейносльдса:
/>
Поскольку воздушный винт находится спереди, то весь фюзеляж обдувается турбулентным потоком, т.е. />. Коэффициент сопротивления одной стороны плоской пластины в таком потоке при заданном числе Рейнольдса Re равен:
/>
Коэффициенты, учитывающие вклад сил давления и эффекта сжимаемости в профильное сопротивление, вычисляются так:
/>
/>/>
Профильное сопротивления фюзеляжа как тела вращения определяют по формуле:
/>
Далее учитываются конструктивные особенности фюзеляжа путём определения приращения коэффициента профильного сопротивления: а) из-за сужения кормовой части:
/>;
б) из-за её скошенности:
/>;
в) под влиянием фонаря кабины:
/>;
г) от установленного в носовой части ПД воздушного охлаждения:
/>
Итак, коэффициент профильного сопротивления фюзеляжа равен:
/>.
Вычислив его для
М=Мmin и Н=0,
получаем: />.
3) Для расчёта профильного сопротивления крыла найдём сначала его число Рейнольдса:
/>.
Коэффициенты, учитывающие вклад сил давления и эффекта сжимаемости в профильное сопротивление, вычисляются так:
/>
/>
Коэффициент профильного сопротивления крыла находится так:
/>
и равен />для М=Мmin и Н=0.
4) Для расчёта профильного сопротивления стабилизатора найдём сначала его число Рейнольдса:
/>.
Коэффициенты, учитывающие вклад сил давления и эффекта сжимаемости в профильное сопротивление, вычисляются так:
/>
/>
Коэффициент профильного сопротивления стабилизатора находится так:
/>
и равен />для М=Мmin и Н=0.
5) Для расчёта профильного сопротивления киля найдём сначала его число Рейнольдса:
/>
Коэффициенты, учитывающие вклад сил давления и эффекта сжимаемости в профильное сопротивление, вычисляются так:
/>
/>
Коэффициент профильного сопротивления киля находится так:
/>
и равен />для М=Мmin и Н=0.
6) Рассматриваемый самолёт является среднепланом. Коэффициент интерференции для крыла и фюзеляжа среднеплана равен />. Тогда приращение коэффициента профильного сопротивления крыла за счёт его взаимодействия с фюзеляжем равно:
/>
для М=Мmin и Н=0
7) Стабилизатор установлен вверху кормовой части фюзеляжа, а значит, коэффициент их интерференции равен />. Тогда приращение коэффициента профильного сопротивления стабилизатора из-за его взаимодействия с фюзеляжем равно:
/>
для М=Мminи Н=0 .8) Учтя вклад всех местных источников сопротивления: антенны, выхлопных патрубков, стыков между листами обшивки, щелей между крылом (оперением) и управляющими поверхностями — получаем следующее приращение к коэффициенту сопротивления, не зависящее от М, Н, и α:
/>
9) Вспомогательная поляра строится при отсутствии тяги двигателя. При этом коэффициенты торможения потока для крыла, стабилизатора и киля равны:
/>10) продолжение
--PAGE_BREAK--
Итак, коэффициент профильного сопротивления равен:
/>
Коэффициент сопротивления при нулевой подъёмной силе равен:
/>
11) Рассчитаем теперь индуктивное сопротивление самолёта, для чего найдём коэффициент отвала поляры:
/>
/>
/>
/>
При М=Мmin: />.
Коэффициент индуктивного сопротивления зависит от коэффициента подъёмной силы следующим образом: />
12) Приращение коэффициента сопротивления с увеличением угла атаки (подъёмной силы) оценивается следующим выражением:
/>
Итак, теперь можно найти коэффициент лобового сопротивления:
/>
Шасси данного самолёта является неубирающимся, поэтому его воздушное сопротивление (по статистике />) учитывается во всех режимах полёта. Вычислим />для нескольких значений угла атаки от />до />и занесём результаты в таблицу 3.1.
По данным этой таблицы строятся график вспомогательной зависимости />и вспомогательная поляра />с разметкой углов атаки на ней (рисунок 3.1 и рисунок 4, где кривая 6 — вспомогательная поляра).
Таблица 3.1.
/>
-2,77
-2
2
4
6
8
10
12
14
16,73
/>
0,060
0,215
0,370
0,525
0,680
0,836
0,991
1,145
1,259
1,315
/>
0,000
0,000
0,000
0,001
0,001
0,003
0,005
0,009
0,016
0,040
/>
0,000
0,003
0,009
0,019
0,031
0,047
0,067
0,089
0,107
0,117
/>
0,041
0,041
0,044
0,050
0,060
0,074
0,091
0,112
0,139
0,164
0,198
/>
Рисунок 3.1 — построение вспомогательной поляры
3.2 Расчёт и построение взлётных поляр
1) При построении данных поляр принимают, что закрылки выпущены под углом 20О, высота полёта нулевая, скорость полёта минимальна (М=Мmin). 2 — 8) Все промежуточные величины, вычисленные в пунктах 2 — 8 при построении вспомогательной поляры (раздел 3.1), вычисляются по тем же формулам и имеют то же численное значение и во взлётном режиме.9) Для учёта вклада обдувки самолёта винтом в лобовое сопротивление найдём коэффициент нагрузки винта по тяге во взлётном режиме:
/>
При этом коэффициенты торможения потока для крыла, киля и стабилизатора равны:
/>
/>
Коэффициент дополнительного сопротивления из-за обдува части крыла винтом определяется по формуле:
/>
/>
/>для М=Мminи Н=0 .10)
Итак, коэффициент профильного сопротивления равен:
/> продолжение
--PAGE_BREAK--
Коэффициент сопротивления при нулевой подъёмной силе равен:
/>
11) Без учёта экрана земли коэффициент отвала поляры рассчитывается по тем же формулам и имеет то же численное значение, что и при расчёте вспомогательной поляры (раздел 3.1, пункт 11): />. С учётом экранного эффекта коэффициент отвала поляры ищется по формулам:
/>
/>
/>
Итак, коэффициент индуктивного сопротивления без учёта и с учётом экранного эффекта ищется по формулам:
/>
12) Приращение коэффициента сопротивления с увеличением угла атаки (подъёмной силы) оценивается следующим выражениями (без учёта и с учётом экранного эффекта):
/>
/>
13) Выпущенные на 20Опри взлёте закрылки увеличивают коэффициент сопротивления на величину:
/>
где />для данного угла отклонения закрылков определяется по справочным данным на основе относительной хорды закрылков.Теперь находим коэффициент лобового сопротивления во взлётном режиме (без учёта и с учётом экранного эффекта):
/>
/>
Вычислим />для нескольких значений угла атаки от />до />и занесём результаты в таблицу 3.2.1. По данным этой таблицы строятся взлётные поляры без учёта экрана земли (рисунок 4, кривая 7).
Таблица 3.2.1
/>
-9,67
-6
-3
3
6
8
10
12
14
15,07
/>
0,285
0,518
0,750
0,983
1,218
1,371
0,991
1,515
1,662
1,669
/>
0,000
0,000
0,001
0,002
0,005
0,007
0,005
0,012
0,027
0,039
/>
0,005
0,018
0,038
0,065
0,100
0,127
0,067
0,155
0,187
0,188
/>
0,111
0,116
0,129
0,150
0,179
0,215
0,245
0,112
0,278
0,325
0,339
Вычислим />для нескольких значений угла атаки от />до />и занесём результаты в таблицу 3.2.2. По данным этой таблицы строятся взлётные поляры с учётом экрана земли (рисунок 4, кривая 8).
Таблица 3.2.2
/>
-9,67 продолжение
--PAGE_BREAK----PAGE_BREAK----PAGE_BREAK--
схi
сха
схi
сха
0,0
0,037
0,042
0,041
0,041
0,041
0,1
0,001
0.038
0,001
0.043
0,001
0.042
0,001
0.042
0,001
0.042
0,2
0,003
0,040
0,003
0,045
0,003
0,044
0,003
0,044
0,003
0,044
0,3
0,006
0.043
0,006
0.048
0,006
0.047
0,006
0.047
0,006
0.047
0,4
0.011
0.048
0.011
0.053
0.011
0.052
0.011
0.052
0.011
0.052
0,5
0.017
0.054
0.017
0.059
0.017
0.058
0.017
0.058
0.017
0.058
0,6
0.024
0.062
0.024
0.067
0.024
0.066
0.024
0.066
0.024
0.066
0,7
0.033
0.072
0.033
0.077
0.033
0.076
0.033
0.076
0.033
0.076
0,8
0.043
0.083
0.043
0.088
0.043
0.087
0.043
0.087
0.043
0.087
0,9
0.054
0,095
0.055
0,100
0.055
0,099
0.055
0,099
0.055
0,099
1,0
0.067
0.110
0.068
0.115
0.068
0.114
0.068
0.114
0.068
0.114
1,1
0.081
0.126
0.082
0.131
0.082
0.130
0.082
0.130
0.082
0.130
1,2
0.096
0.146
0.097
0.151
0.097
0.150
0.097
0.150
0.097
0.150
1,315
0.116
0.187
0.117
0.192
0.117
0.191
0.117
0.191
0.117
0.191
/>
Рисунок 5 — Крейсерские поляры и зависимости суа (α).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Остовлавский И.В. Аэродинамика самолёта. — М.: Оборонгиз, 1957.
2. Меньшиков В.И. Аэродинамические характеристики самолётов: Учебное пособие. — Харьков: Харьк. Авиац. Ин-т, 1984.
3. Мхитарян А.М. Аэродинамика. — М.: Машиностроение, 1976.