Студопедия КАТЕГОРИИ: АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция |
Рівняння просторового руху літакаСтр 1 из 16Следующая ⇒
Системи координат, їх взаємозв'язок При розв'язанні задачі управління польотом літака рекомендовано використовувати нерухомі та рухомі праві прямокутні декартові системи координат. Нерухомі системи координат, у яких осі та початок координат фіксовані відносно Землі, використовуються як системи відліку швидкостей, прискорень, переміщень ЛА. Як нерухома частіше всього застосовується нормальна система координат O0XgYgZg , вісь O0Yg якої спрямована уверх за місцевою вертикаллю, а напрями осей O0Xg i O0Zg обираються відповідно із задачею. Рухомі системи координат, початок яких розташовується в центрі мас ЛА (звідси іх друга назва - ЛАцентричні системи координат), як правило, використовуються для запису рівнянь руху ЛА. Як рухомі частіше всього використовуються: нормальна, зв'язана, швидкісна та траєкторна системи координат. Нормальна система координат OXgYgZg - це рухома система координат, вісь OYg якої спрямована уверх за місцевою вертикаллю, а напрям осей OXg i OZg обирається відповідно із задачею. Зв'язана система координат OXYZ - це рухома система координат, осі якої фіксовані відносно ЛА. Поздовжня вісь OX розташовується в площині симетрії ЛА та спрямована від хвостової до носової його частини. Нормальна вісь OY розташовується в площині симетрії ЛА та спрямована до верхньої його частини. Поперечна вісь OZ перпендикулярна площині симетрії ЛА та спрямована у правий його бік. Осі зв'язаної системи координат, як правило, збігаються з головними осями інерції ЛА. Швидкісна система координат OXaYaZa - це рухома система координат, швидкісна вісь OXa якої збігається з напрямком вектора повітряної швидкості (швидкість ЛА відносно повітряного середовища), вісь піднімальної сили OYa розташовується в площині симетрії ЛА та спрямована до верхньої його частини, а бічна вісь OZa доповнює осі OXa i OYa до правої системи координат. Траєкторна система координат OXкYкZк - це рухома система координат, вісь OXк якої збігається з напрямком вектора земної швидкості (швидкість ЛА відносно земної поверхні), а вісь OYк спрямована уверх від поверхні Землі й лежить у вертикальній площині, яка проходить крізь вісь OXк. При відсутності вітру вісь OXктраєкторної системи кооринат сполучається з віссю OXa швидкісної системи координат. Рівняння просторового руху літака Рух літака у просторі визначається з урахуванням поступального руху центра мас й обертального руху відносно центра мас. При записуванні диференціальних рівнянь руху в проекціях на осі нерухомої системи координат виникають труднощі проектування сил і моментів, що діють на літак, на осі нерухомої системи координат. До того ж сильно ускладнюється подальший аналіз таких рівнянь. Тому в динаміці польоту віддають перевагу рівнянням руху, що записані в рухомих системах координат. При отриманні рівнянь просторового руху приймаються такі припущення: - літак являє собою абсолютне жорстке тіло з постійними інерційно-масовими характеристиками; - осі зв'язаної системи координат збігаються з головними осями інерції, тобто, відцентрові моменти інерції дорівнюють нулю; - вектор тяги прикладений до центра мас і не створює додаткових моментів; - Земля плоска, не обертається і не переміщується в інерціальному просторі, тобто не враховуються переносні та коріолісові сили й моменти; - повітряне середовище нерухоме. Рівняння руху літака як твердого тіла можуть бути одержані із законів збереження кількості та моменту кількості руху, які у векторній формі мають вигляд:
, (1.10) де m - маса літака; - вектор повітряної швидкості; - вектор зовнішніх сил, діючих на літак; - вектор моменту кількості руху; - аеродинамічний момент. Із векторного аналізу відомо, що похідна вектора у нерухомій системі координат може бути подана як сума похідної вектора в рухомій системі координат та векторного добутку кутової швидкості обертання рухомої системи координат відносно нерухомої на вектор : . Саме в такому вигляді векторні рівняння (1.10) проектують на осі рухомих систем координат. Із урахуванням прийнятих припущень рух літака як твердого тіла з шістьма ступенями свободи описується шістьма рівняннями Ейлера з доповненням цих рівнянь сил і моментів кінематичними та геометричними співвідношеннями, а також рівняннями, які описують траєкторний рух центра мас літака. Додаткові рівняння дозволяють замкнути систему рівнянь, що описує просторовий рух літака:
Докладний вивід системи рівнянь, що описує просторовий рух літака, наведений в дод. 1. У системі рівнянь просторового руху літака (1.11) перші три рівняння сил розкривають закон збереження кількості руху. Рівняння записані в траєкторній системі координат і визначають зміни швидкості польоту V , кута нахилу траєкторії Q і кута шляху Y під впливом тяги Р, сили ваги G і аеродинамічних сил: піднімальної сили Ya, сили лобового опору Хa, бічної сили Za . При записуванні в траєкторній системі координат у рівняннях сил добре проглядається динаміка зміни параметрів польоту під дією аеродинамічних сил і, крім того, рівняння сил мають простий вигляд, оскільки вісь ОХк траєкторної системи координат збігається з вектором швидкості, а проекції вектора швидкості на осі ОYк, ОZк дорівнюють нулю. Наступні три рівняння системи розкривають закон збереження моменту кількості руху. Ці рівняння записані у зв'язаній системі координат і описують зміну швидкості обертання літака під впливом аеродинамічних моментів Мx, My, Mz . Осі зв'язаної системі координат збігаються з головними осями інерції, тому рівняння моментів спрощуються, оскільки відцентрові моменти інерції Ixy , Iyz , Izx дорівнюють нулю, і проекції моменту кількості руху на осі зв'язаної системи координат мають простий вигляд:
де Ix, Iy, Iz - осьові моменти інерції; wx, wy, wz - проекції вектора кутової швидкості обертання літака на осі зв'язаної системи координат. Наступні три рівняння - це геометричні співвідношення для обчислення кута атаки a , кута ковзання b та швидкісного кута крену ga , які необхідні для визначення сил, що діють на літак, у проекціях на осі траєкторної системи координат. Геометричні співвідношення отримані шляхом прирівнювання матриці напрямляючих косинусів прямого переходу від зв'язаної системи координат до швидкісної до добутку матриць послідовного переходу від зв'язаної системи координат до нормальної і далі від нормальної до швидкісної. Така процедура може бути записана у вигляді матричного рівняння й ілюструється графом (рис. 1.7). Кінематичні співвідношення доповнюють систему рівнянь до повної й описують зміну кутів крену g, тангажа J і рискання y. Кінематичні рівняння отримані шляхом запису проекцій швидкостей зміни кутів крену , тангажа і рискання ( проекцій кутової швидкості обертання зв`язаної системи координат навколо нормальної ) на осі зв'язаної системи координат у вигляді: з подальшим розв`язанням отриманої системи рівнянь. Аналогічно для рівняння сил отримані проекції кутової швидкості обертання траєкторної системи координат навколо нормальної: Останні три рівняння системи (1.11) необхідні для дослідження траєкторного руху літака. Крім цього, перше з останніх трьох рівнянь необхідне для виявлення функціональних залежностей аеродинамічних сил і моментів, а також тяги двигуна від густини повітря, яка змінюється з висотою польоту. Щоб система рівнянь (1.11) була замкненою необхідно описати залежності аеродинамічних сил і моментів, а також тяги двигуна від параметрів польоту. Загальна характеристика сил і моментів, що діють на літак Розглядаючи залежність аеродинамічних сил і моментів від параметрів польоту, вважають: (1.12) , де - коефіцієнт аеродинамічної піднімальної сили, який залежить від кута атаки a, швидкості польоту V та кута відхилення руля висоти dв; - коефіцієнт сили лобового опору; - коефіцієнт аеродинамічної бічної сили, який залежить від кута ковзання b та кута відхилення руля направлення dн; r = r(H)- густина повітря; S - площа крила; Залежність Р = Р(V, Н) має вигляд, проілюстрований рис.1.8, і описує висотно-швидкісні характеристики авіаційного двигуна. Зі збільшенням висоти польоту витрата повітря Gпов крізь двигун зменшується. Це спричиняє зменшення тяги Р. У діапазоні висот від 0 до 11 км внаслідок зменшення температури тяга падає менш інтенсивно (зменшуються витрати енергії на роботу стиснення в компресорі). Починаючи з 11 км, температура повітря залишається постійною, тому інтенсивність падіння тяги двигуна зростає. Характер зміни тяги Р = Рпит Gпов від швидкості польоту V визначається відповідними змінами питомої тяги Рпит і витрати повітря Gпов. Зі збільшенням швидкості польоту Рпит » Сс _ V зменшується, оскільки швидкість виходу газу з реактивного сопла Ссчерез зниження ступеню підігріву в камері згорання збільшується повільніше, ніж зростає швидкість польоту V. При наближенні до деякої крайньої швидкості польоту Рпит прагне до нуля. Витрата повітря Gпов при збільшенні швидкості польоту збільшується. Але на малих швидкостях, коли стисливість повітря виявляється незначно, збільшення Gпов невелике, тому тяга Р = Рпит Gпов через зменшення Рпит може знижуватися. На великих швидкостях, коли витрата повітря починає інтенсивно збільшуватися, тяга зростає. Однак після певної швидкості зменшення Рпит починає перевищувати збільшення Gпов і тяга швидко падає. Залежність Р = Р(dс.г ) описує закон зміни тяги двигуна від положення сектора газу.
Найбільш істотні залежності коефіцієнтів аеродинамічних сил і моментів такі: Залежність характеризує зміну аеродинамічної піднімальної сили від кута атаки, яка залежить від розподілу коефіцієнтів тиску вздовж поверхні крила. При безвідривному обтіканні профілю крила
,
При виникненні зриву потоку з верхньої поверхні крила Залежність характеризує зміну коефіцієнта аеродинамічної піднімальної сили від стисливості повітря. На рис. 1.10 зображена залежність похідної від числа Маха – числа M польоту. На малих швидкостях стисливість повітря практично не впливає на характер розподілу тиску вздовж профілю крила, тобто не змінюється. При Надзвукова зона на нижній поверхні крила швидко досягає задньої крайки крила (нижня поверхня більш плоска ніж верхня), коефіцієнт знову починає зростати. При М > 1профіль крила вже обтікається надзвуковим потоком без місцевих стрибків згущення повітря. Збільшення числа М надзвукового потоку призводить до зменшення коефіцієнта через втрату повного тиску на головному стрибку згущення. Залежність характеризує зміну коефіцієнта лобового опору
при зміні кута атаки та швидкості польоту. Тут - коефіцієнт сили лобового опору при нульовій піднімальній силі характеризує профільний опір (опір сили тертя), який не залежить від швидкості польоту, і хвильовий опір стрибків згущення повітря, які виникають при М > Мкр ( рис. 1.11) і збільшують коефіцієнт . При
Залежності , характеризують додаткові зміни піднімальної та бічної сили літака з урахуванням сил вертикального та горизонтального оперення при відхиленнях органів управління. Залежність mz(a) характеризує поздовжню статичну стійкість літака з перевантаження. Коефіцієнт аеродинамічного моменту тангажа (без урахування демпфірування) , (1.13) де - коефіцієнт аеродинамічного моменту тангажа при = 0; xт - координата центра мас літака; xF - координата аеродинамічного фокуса F . Узявши похідну від mz за a, отримаємо: . Таким чином, mz(a) визначається, в основному, відстанню між центром мас літака й аеродинамічним фокусом F, в якому прикладається приріст піднімальної сили DYa (рис.1.12). При задньому розташуванні фокуса літак є статично стійким з перевантаження, тобто, при зміні кута атаки на величину Da приріст піднімальної сили DYa спричиняє момент Mz(a), спрямований на усунення відхилення кута атаки. На рис. 1.13 показаний характер зміни коефіцієнта аеродинамічного моменту тангажа mz від a . Видно, що при a = aн раніше стійкий літак стає статично нейтральний з перевантаження, оскільки приріст кута атаки не викликає приріст моменту Mz, а при a > aн літак стає статично нестійкий: позитивний приріст кута атаки спричиняє позитивний приріст моменту Mz. Тут можливий підхват літака з перевантаження, тобто, несвавільне збільшення перевантаження при фіксованому положенні руля висоти. Залежності mz(wz), mx(wx), my(wy) характеризують власні аеродинамічні демпфіруючі моменти тангажа, крену та рискання, які гасять кутові коливання літака. При обертанні літака його частини обтікаються додатковим повітряним потоком. Цей потік змінює місцеві кути атаки та ковзання, що викликає додаткові аеродинамічні сили, які відносно центра мас створюють демпфіруючі моменти. Демпфіруючі моменти діють тільки при наявності кутової швидкості й спрямовані у бік, протилежний обертанню. У літаків нормальної схеми демпфіруючі моменти створюються головним чином горизонтальним і вертикальним опереннями (г.о і в.о). Рис.1.14,а ілюструє створення демпфіруючого моменту тангажа Mz(wz). Тут , , – відповідно прирости піднімальних сил горизонтального оперення, передньої та задньої частин фюзеляжу, які створюються за рахунок зміни місцевих кутів атаки. Місцеві кути атаки елементів конструкції літака виникають за рахунок приросту лінійних швидкостей цих елементів при обертанні літака навколо поперечної осі ОZ. При обертанні літака навколо поздовжньої осі ОХ (рис. 1.14,б) демпфіруючий момент створюється крилом. У крила, що опускається, піднімальна сила за рахунок збільшення місцевих кутів атаки вздовж розмаху крила зростає, у крила, що підіймається, піднімальна сила зменшується. Це викликає аеродинамічний момент крену Mx(wx), спрямований проти кутової швидкості wx. Залежність характеризує демпфіруючий момент тангажа, що обумовлюється запізненням скосу потоку. При створенні піднімальної сили крило відкидає повітряний потік униз зі швидкістю Vy, внаслідок чого вектор швидкості повертається на кут скосу потоку e . Піднімальна сила горизонтального оперення, яке обтікається цим потоком (рис.1.15), визначається як .
Залежність mz(V) характеризує поздовжню статичну стійкість літака з швидкості. Аналізуючи формулу (1.13), відмітимо: . Залежність коефіцієнта аеродинамічного моменту тангажа від швидкості польоту виявляється, в основному, на надзвукових швидкостях і визначається зміною положення аеродинамічного фокуса при М > 1, а також залежностями і . Залежність mz = f(М) проілюстрована на рис. 1.16. У статично стійкого зі швидкості літака збільшення швидкості викликає приріст піднімальної сили без порушення балансування за моментами, які створюються крилом і стабілізатором. Приріст піднімальної сили викликає вертикальну швидкість, і траєкторія польоту буде викривлюватися вверх. Почне збільшуватись кут нахилу траєкторії з одночасним зменшенням кута атаки, що в свою чергу викликає у статично стійкого з перевантаження літака кабрируючий момент, який збільшує кут тангажа й припиняє зменшення кута атаки. Літак переходить у кабрирування, при цьому з’являється складова сили ваги, яка діє проти сили тяги й відновлює початкову швидкість. У статично нестійкого зі швидкості літака на числах М > 1 при збільшенні швидкості аеродинамічний фокус зміщується назад і приріст піднімальної сили викликає пікіруючий момент. Літак затягується в пікірування і починає зменшувати кут нахилу траєкторії, а складова сили ваги у цьому випадку додатково збільшує швидкість польоту. Залежність my(b) характеризує шляхову (флюгерну) статичну стійкість літака і визначається в основному характеристиками вертикального оперення та його відстанню до центра мас літака. Момент шляхової стійкості My(b) викликається поперечною силою Z(b) ( рис. 1.17), яка виникає при ковзанні літака і точка прикладання якої у статично стійкого у шляховому відношенні літака розташовується поза центром мас. У цьому випадку момент шляхової стійкості My(b) утримує літак наче флюгер у повітряному струмі. Залежність mx(b) характеризує поперечну статичну стійкість літака й визначається ( рис.1.18 ) затінком півкрила фюзеляжем, а також силою Zво, яка виникає при несиметричному обтіканні вертикального оперення, і точка прикладання якої розміщується вище центра мас літака. Крім того mx(b) залежить від поперечного ”V” та стрілоподібності крила, оскільки за рахунок ковзання змінюється фактична стрілоподібність півкрила. Залежність mx(wy) характеризує спіральний момент крену. При обертанні літака навколо нормальної осі ОY ( рис. 1.19) змінюються місцеві швидкості вздовж розмаху крила. Одне півкрило, що рухається назустріч повітряному потоку, придбає, а друге втрачає прирости лінійної швидкості DVx. Це супроводжується зміною піднімальної сили півкрил на величину DYa і появою моменту Мx(wy), який нахиляє літак у бік розвороту. Вертикальне оперення також бере участь у створенні спірального моменту, тому що при обертанні літака змінюється місцевий кут ковзання вертикального оперення й виникає сила Zво, точка прикладання якої розташовується вище поздовжньої осі. Залежність my(wx) характеризує спіральний момент рискання. Визначається двома факторами. З одного боку, при обертанні літака з кутовою швидкістю wx ( рис. 1.20) вздовж розмаху крила змінюються кути атаки й піднімальні сили півкрил (у крила, що опускається, вона збільшується), як наслідок змінюються поздовжні сили півкрил Х » Хa - – Yasin a (у крила, що опускається, вона зменшується), з'являється момент, який розвертає літак у бік крила, що підіймається. З іншого боку вертикальне оперення при обертанні літака з кутовою швидкістю wx обтікається додатковим повітряним потоком, що в свою чергу викликає поперечну силу Zв.о і момент Мy, який розвертає літак у бік крила, що опускається. Співвідношення цих двох моментів обумовлює знак спірального моменту рискання. На великих кутах атаки вертикальне оперення затінюється і головну роль відіграє момент від зміни поздовжніх сил півкрил, а на малих кутах атаки цей момент незначний і визначну роль відіграє момент вертикального оперення. Залежності mz(dв), mx(dе), my(dн) характеризують ефективність руля висоти, елеронів і руля направлення відповідно. Залежності mx(dн), my(dе) характеризують перехресний вплив органів управління бічним рухом. |
|||||||||||||||
Последнее изменение этой страницы: 2018-05-31; просмотров: 1620. stydopedya.ru не претендует на авторское право материалов, которые вылажены, но предоставляет бесплатный доступ к ним. В случае нарушения авторского права или персональных данных напишите сюда... |