- •Реферат
- •Ключевые слова
- •Аннотация
- •Оглавление
- •Расчёт и анализ аэродинамических характеристик ла
- •Введение
- •Лобовое сопротивление ла
- •Общая характеристика лобового сопротивления
- •Коэффициент сопротивления трения ла
- •Коэффициент донного сопротивления ла
- •Коэффициент волнового сопротивления ла
- •Коэффициент лобового сопротивления ла
- •2.3 Производная от коэффициента подъёмной силы ла по углу атаки
- •Интерференция корпуса и несущих поверхностей
- •Скос потока
- •Производная коэффициента подъёмной силы ла по углу атаки.
- •Фокус ла
- •Нелинейность аэродинамических коэффициентов
- •Зависимость коэффициента подъёмной силы ла от угла атаки
- •Поляра ла
- •2.4.5. Максимальное аэродинамическое качество
- •Расчёт и анализ траектории наведения ла
- •Общие сведения о траектории наведения
- •Определение исходных данных для перехватчика и цели
- •Формирование файла аэродинамических характеристик
- •Результаты расчёта траектории перехвата на эвм
- •Анализ параметров траектории
- •Численное интегрирование уравнений движения центра масс ла для участка траектории
- •Определение угла атаки и вычисление левых частей уравнений
- •Общие сведения о численном интегрировании обыкновенного дифференциального уравнения
- •Метод Эйлера
- •Модифицированный метод Эйлера
- •Метод Рунге-Кутта
- •Дополнительный участок интегрирования
- •Определение значения сил в заданный момент времени
- •Определение угла атаки и угла отклонения рулевых поверхностей
- •Выбор основных параметров для расчёта динамических коэффициентов
- •Определение параметров, связанных с отклонением управляющих поверхностей
- •Определение параметров, связанных с демпфирующими моментами
- •Уточнение угла атаки и определение угла отклонения управляющих поверхностей
Определение параметров, связанных с демпфирующими моментами
Для первых несущих поверхностей:
где
,
и
.
Таблица 5.4. –
Определение коэффициента
,
Cyak1(alp) |
(Xfk1-Xц.м.)^2 |
-mzk1(omeg) |
0,56228 |
0,23863225 |
0,134178142 |
Cyak1(f)(alp) |
(Xfk1(f)-Xц.м.)^2 |
-mzk1(f)(omeg) |
0,22042 |
0,23687689 |
0,052212404 |
Cyaf(k1)(alp) |
(Xff(k1)-Xц.м.)^2 |
-mzf(k1)(omeg) |
0,37788 |
0,2304 |
0,087063552 |
|
|
-mz1(omeg) |
|
|
0,273454098 |
.
Для вторых несущих поверхностей:
;
где
,
и
Таблица
5.5 – Определение коэффициента
Cyak2(alp) |
(Xfk2-Xц.м.)^2 |
-mzk2(omeg) |
0,87952 |
3,65689129 |
3,216309027 |
Cyak2(f)(alp) |
(Xfk2(f)-Xц.м.)^2 |
-mzk2(f)(omeg) |
0,27535 |
3,64619025 |
1,003978485 |
Cyaf(k2)(alp) |
(Xff(k2)-Xц.м.)^2 |
-mzf(k2)(omeg) |
0,46301 |
3,6864 |
1,706840064 |
|
|
-mz2(omeg) |
|
|
5,927127577 |
.
Значение производной демпфирующего момента от запаздывания скоса потока в области вторых несущих поверхностей:
Значение производной демпфирующего момента для корпуса:
Значение производной демпфирующего момента для всего ЛА от его вращения относительно оси OZ:
Уточнение угла атаки и определение угла отклонения управляющих поверхностей
В задачах первой
группы динамики ЛА моделировался
материальной точкой. При расчётах
использовались аэродинамические силы
ЛА, которые были получены в стационарных
условиях, без учёта угловой скорости
вращения ЛА и отклонения управляющих
поверхностей. Для нахождения угла атаки
и определения угла отклонения рулевых
поверхностей можно считать, что траектории
ЛА, посчитанные при решении первой
группы задач динамики полёта, будут
отличаться от аналогичных траекторий
для задач второй группы. Если попробовать
сохранить хотя бы форму траектории, то
можно предположить, что для обеих
траекторий в исследуемой точке должны
быть одинаковыми
,
а, следовательно, и
.
При исследовании продольного движения ЛА:
Если предположить, что последние три слагаемых продольного момента линейно зависят от своих аргументов, то можно получить следующую итерационную формулу для определения угла атаки:
Здесь
;
;
,
и
.
При определении
производных
можно воспользоваться следующими
формулами:
;
;
;
;
;
;
;
;
;
Таблица 5.6. – определение угла атаки
Определение |
угла атаки |
|
|
|
|
|
A |
alp0 |
deltaCya/alp |
Cya(alp)+ deltaCya/alp |
(Cya(alp)+ deltaCya/alp) *B |
C |
Alp(i+1) |
0,027321721 |
0 |
0 |
3,2802 |
3,75440714 |
3,913479097 |
0,400007 |
|
0,400007067 |
0,00159 |
3,28179 |
3,756227 |
3,915298958 |
0,399821 |
Определяем угол отклонения управляющей поверхности:
;
