- •Додаток1. Математична модель просторового руху літака
- •Рівняння сил
- •Рівняння моментів
- •Геометричні рівняння
- •Кінематичні рівняння
- •Рівняння траєкторного руху
- •Додаток 2. Вивід лінеарізованих рівнянь на основі повної системи рівнянь ізольованого поздовжнього руху
- •Додаток 3. Методика аналізу статичних та динамічних характеристик контурів автоматичного управління
Додаток1. Математична модель просторового руху літака
При отриманні рівнянь просторового руху приймаються такі припущення:
-літак являє собою абсолютне жорстке тіло з постійними інерційно-масовими характеристиками;
- осі зв'язаної системи координат збігаються з головними осями інерції, тобто, відцентрові моменти інерції дорівнюють нулю;
- вектор тяги прикладений до центра мас і не створює додаткових моментів;
- Земля плоска, не обертається і не переміщується в інерціальному просторі, тобто, не враховуються переносні та коріолісові сили й моменти;
- повітряне середовище нерухоме.
Рівняння руху літака, як твердого тіла, можуть бути одержані з законів збереження кількості та момента кількості руху, які у векторній формі мають вигляд:
, (Д.1)
де
- вектор повітряної швидкості;
- вектор зовнішніх сил, діючих
на літак;
- вектор момента кількості
руху;
- аеродинамічний момент;
m - маса літака.
З векторного аналізу відомо,
що похідна вектора
у нерухомій системі координат -
може
бути подана як сума похідної вектора в
рухомій системі координат
та векторного добутку кутової швидкості
обертання рухомої системи координат
відносно нерухомої на вектор
.
. (Д.2)
Саме в такому вигляді векторні рівняння (Д.1) проектують на осі рухомих систем координат.
Складові векторного добутку
(
)
визначаютьсчя як алгебраїчне доповнення
першої строки детермінанта
і у проекціях на осі зв`язаної системи координат мають вигляд:
на вісь OX
;
на вісь OY
; (Д.3)
на вісь OZ
.
Тут x, y, z - проекції вектора абсолютної кутової швидкості на осі зв`язаної системи координат;
Ax, Ay, Az - проекції вектора A на осі зв`язаної системи координат.
Рівняння сил
Використовуючи формулу повної похідної вектора (Д.2), (Д.3) можна записати перше векторне рівнянння (Д.1) - закон збереження кількості руху, як рівняння сил у проекціях на осі зв`заної системи координат:
(Д.4)
де Fx, Fy, Fz - проекції вектора зовнішніх сил на осі зв`язаної системи координат;
Vx, Vy, Vz - проекції вектора швидкості руху центра мас ЛА на осі зв`язаної системи координат.
Але в траєкторній системі координат рівняння сил мають більш простий вигляд, оскільки вісь ОХк траєкторної системи координат збігається з вектором швидкості, а проекції вектора швидкості на осі ОYк , ОZк дорівнюють нулю:
Крім того при записуванні в траєкторній системі координат у рівняннях сил добре проглядається динаміка зміни параметрів польоту під впливом аеродинамічних сил. Отже, рівняння сил, що записані в траєкторній системі координат, мають вигляд:
(Д.5)
де Fxк, Fyк, Fzк,yк, zк - проекції вектора зовнішніх сил і вектора абсолютної кутової швидкості обертання ЛА на осі траєкторної системи координат.
Розкриваючи праві частини рівнянь, будемо враховувати силу тяги двигуна Р, силу ваги G та повну аеродинамічну силу R, яка складається з підйомної сили Ya, сили лобового опору Хa, бокової сили Za, тобто
Fxк Pxк + Gxк + Rxк ; Fyк Pyк + Gyк + Ryк ; Fzк Pzк + Gzк + Rzк
Аеродинамічні сили визначаються в швидкісній системі координат і спрямовані за відповіднимі осямі цієї системи, причому напрямок сили лобового опору Xa збігається з негативним напрямком осі OXa , яка у свою чергу збігається з напрямком осі OXк траєкторної системи координат. Швидкісна система координат нахилена відносно траєкторної на швидкісний кут крену a, тому проекції повної аеродинамічної сили R на вісі траєкторної системи координат (див. мал. Д.1) визначаються як:
(Д.6)
Вектор сили тяги P
збігається з напрямком поздовжньої осі
ЛА, тому його проекції на вісі траєкторної
системи координат можна отримати за
допомогою матриці
переходу від зв`язаної до траєкторної
системи координат:
(Д.7)
Ці співвідношення можна також отримати, проектуючи вектор сили тяги P на осі траєкторній системи координат:
(Д.8)
Напрямок вектора сили ваги G збігається з негативним напрямком осі Yg нормальної системи координат, яка нахилена відносно площини горизонту на кут (див. мал. Д.2), тому проектуючи вектор сили ваги на вісі траєкторної системи координат можна отримати:
(Д.9)
Використовуючи співвідношення (Д.6), (Д.8), (Д.8), можна записати рівняння сил в траєкторній системі координат у вигляді:
(Д.10)
