- •Реферат
- •Ключевые слова
- •Аннотация
- •Оглавление
- •Расчёт и анализ аэродинамических характеристик ла
- •Введение
- •Лобовое сопротивление ла
- •Общая характеристика лобового сопротивления
- •Коэффициент сопротивления трения ла
- •Коэффициент донного сопротивления ла
- •Коэффициент волнового сопротивления ла
- •Коэффициент лобового сопротивления ла
- •2.3 Производная от коэффициента подъёмной силы ла по углу атаки
- •Интерференция корпуса и несущих поверхностей
- •Скос потока
- •Производная коэффициента подъёмной силы ла по углу атаки.
- •Фокус ла
- •Нелинейность аэродинамических коэффициентов
- •Зависимость коэффициента подъёмной силы ла от угла атаки
- •Поляра ла
- •2.4.5. Максимальное аэродинамическое качество
- •Расчёт и анализ траектории наведения ла
- •Общие сведения о траектории наведения
- •Определение исходных данных для перехватчика и цели
- •Формирование файла аэродинамических характеристик
- •Результаты расчёта траектории перехвата на эвм
- •Анализ параметров траектории
- •Численное интегрирование уравнений движения центра масс ла для участка траектории
- •Определение угла атаки и вычисление левых частей уравнений
- •Общие сведения о численном интегрировании обыкновенного дифференциального уравнения
- •Метод Эйлера
- •Модифицированный метод Эйлера
- •Метод Рунге-Кутта
- •Дополнительный участок интегрирования
- •Определение значения сил в заданный момент времени
- •Определение угла атаки и угла отклонения рулевых поверхностей
- •Выбор основных параметров для расчёта динамических коэффициентов
- •Определение параметров, связанных с отклонением управляющих поверхностей
- •Определение параметров, связанных с демпфирующими моментами
- •Уточнение угла атаки и определение угла отклонения управляющих поверхностей
Дополнительный участок интегрирования
Из-за простоты метода наведения и, следовательно, простоты вычислений при итерациях, а также для закрепления методов численного интегрирования дополнительно проинтегрируем участок введения в режим наведения.
Таблица 4.17. – Метод средней ординаты
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
t,s |
V,m/s |
X,m |
H,m |
Xc,m |
Rc,m |
β° |
q° |
e° |
4 |
123 |
164 |
160 |
|
|
|
19 |
|
4,5 |
139 |
227 |
176 |
11200 |
10973,26 |
-0,39683 |
9,6 |
-9,99683 |
5 |
155,7127 |
301,0534 |
183,1809 |
|
|
|
-0,16707 |
|
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
cos(q) |
sin(q) |
V*cos(q) |
V*sin(q) |
m,kg |
P(H) |
a(H) |
g(H) |
mg |
0,985996 |
0,166769 |
137,0534 |
23,18086 |
38,675 |
99640 |
339,753 |
9,806 |
379,2471 |
19 |
20 |
21 |
22 |
23 |
24 |
25 |
26 |
27 |
M |
qS |
Cxtr |
Cxд |
Cxво |
Cx0 |
Cyalp |
Alp1 |
Cx1 |
0,409121 |
922,278 |
0,04173 |
0 |
0 |
0,04173 |
2,9061 |
-16 |
0,28977 |
28 |
29 |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
Cxa |
Xa |
Nxa |
Vdot |
DCY/AR |
Cya |
Ya |
Nya |
TetaDot |
0,3315 |
305,7352 |
3,502756 |
32,71269 |
0,8055 |
-1,0364 |
-955,847 |
-3,75594 |
-19,1671 |
Таблица 4.18. – Определение угла атаки в точке t=4.5с
Teta(dot) |
Nya |
A |
Alpha0 |
Dcya(alp) |
Cya(alp) |
Ya(alp) |
B |
Alpha1 |
|
-0,00981 |
-3,74845 |
-1421,59 |
0 |
0 |
2,9061 |
2680,232 |
4380,232 |
-0,32455 |
-18,5952 |
-0,334 |
-3,75 |
-1421,59 |
-0,32455 |
0,8055 |
3,7116 |
3423,127 |
5123,127 |
-0,27748 |
-15,8987 |
-0,33 |
-3,75 |
-1421,59 |
-0,27748 |
0,803 |
3,7091 |
3420,821 |
5120,821 |
-0,27761 |
-15,9059 |
При вычислении принимаем значение α(t=4.5c) известным и равным -16
Применив все три метода численного интегрирования, мы убедились, что наиболее точным является метод Рунге-Кутта. Хотя данный метод является наиболее трудоёмким. Менее трудоёмким является модифицированный метод Эйлера. Наименее точным является метод Эйлера.