Matematyka i fizyka od podstaw
PLEN

Dynamika lotu › II. Samolot › 8 szkic · PDF do druku

Zanim zaczniesz Automatyka 2.2 Bieguny i zera pętli zamkniętej, linie pierwiastkoweAutomatyka 3.3 Nieliniowa rzeczywistość: nasycenie, windup, feedforward i gain scheduling7 Ruch boczny: przechylenie, spirala i holendrowanie

Nowe pojęcia w tej lekcji (1)

augmentacja stateczności

Dynamika lotu 8

Augmentacja stateczności: tłumik odchylenia i pochylenia

Gdy mody samolotu są za słabo tłumione, dokłada się tłumienie elektronicznie: żyroskop mierzy prędkość kątową, a mały serwomechanizm wychyla ster proporcjonalnie do niej. To najprostsza pętla sprzężenia zwrotnego w samolocie i ma jedną pułapkę: w zakręcie prędkość odchylania jest potrzebna, a tłumik próbuje ją zgasić. Ta lekcja buduje tłumik odchylenia z filtrem washout, tłumik pochylenia, programuje wzmocnienie w funkcji prędkości i kończy na samolocie niestatecznym statycznie, który lata tylko dzięki regulatorowi.

Czas
ok. 2 h teorii + ok. 2,5 h zadań
Zadania
14 zadań w treści + 5 do przemyślenia, z odpowiedziami
Źródła
E. H. J. Pallett, S. Coyle, Automatic Flight Control, rozdz. 2 (s. 79), 3 (s. 122–124), 7 (s. 222–224), 11; A. Tewari, rozdz. 3.6 i 3.9

0 Prolog: osobny podsystem

Pallett opisuje tłumik odchylenia (ang. yaw damper) jako osobny podsystem, niezależny od autopilota: pilot może go włączyć, nawet gdy pilotuje ręcznie, bo holendrowanie niektórych samolotów jest za słabo tłumione, żeby człowiek mógł nad nim panować. W samolotach z dwiema częściami steru kierunku, jak Boeing 747, każda ma własny tłumik. Typowy tłumik jest „szeregowy”: dodaje swoje wychylenie do wychylenia pilota, nie ruszając pedałów, więc pilot go nie czuje. Pallett odnotowuje też zamieszanie z nazwą: niektóre samoloty mają „tłumiki odchylenia”, które są zwykłymi amortyzatorami pedałów.

Taki podział, osobno augmentacja stateczności (SAS, ang. stability augmentation system), osobno autopilot, jest też u Tewariego: SAS to wewnętrzna pętla, która poprawia mody, a autopilot to zewnętrzna pętla, która wydaje polecenia. Tę lekcję poświęcamy pętli wewnętrznej, następną zewnętrznej.

Pytanie na start. Tłumik odchylenia wychyla ster kierunku proporcjonalnie do prędkości odchylania rr, żeby gasić holendrowanie. Co zrobi w ustalonym zakręcie, w którym rr jest stałe i niezerowe? Odpowiedź w części 2.

1 Tłumik odchylenia

W lekcji 7 Ruch boczny: przechylenie, spirala i holendrowanie holendrowanie samolotu lekkiego miało ζ=0,19\zeta = 0{,}19: znośnie, ale wyraźnie kołysze. Dodajemy sprzężenie od prędkości odchylania do steru kierunku:

Tłumik odchylenia

δr=K r⟹Nr  →  Nr+NδrK,Yr,Lr podobnie (mniejsze zmiany).\delta_r = K\,r\qquad\Longrightarrow\qquad N_r \;\to\; N_r + N_{\delta r}K,\qquad Y_r, L_r \text{ podobnie (mniejsze zmiany)}.

Przy Nδr<0N_{\delta r} < 0 (ster wychylony w lewo obraca nos w lewo) dodatnie KK zwiększa tłumienie odchylania. W macierzy zamkniętej pętli zmienia się tylko trzecia kolumna: Acl=A+br[0 0 K 0]A_{cl} = A + \mathbf b_r[0\ 0\ K\ 0].

W naszym modelu Nδr=−5,69N_{\delta r} = -5{,}69/s² na radian. Przy K=0,5K = 0{,}5 (pół radiana steru na radian na sekundę, czyli 0,5∘0{,}5^\circ steru na każdy 1∘1^\circ/s) efektywne NrN_r rośnie z −1,04-1{,}04 do −3,89-3{,}89/s, a holendrowanie dostaje ζ=0,65\zeta = 0{,}65 przy prawie niezmienionej częstości (3,03 rad/s). To tak, jakby statecznik pionowy miał prawie dwa razy dłuższe ramię (Cnr∝lV2C_{n_r} \propto l_V^2), ale tylko dla tłumienia, bez zmiany stateczności kierunkowej. Linie pierwiastkowe (lekcja 2.2 Bieguny i zera pętli zamkniętej, linie pierwiastkowe) pokazują, jak bieguny holendrowania przesuwają się w lewo przy rosnącym KK.

Symulacja 1 Tłumik odchylenia na samolocie z lekcji 7. Lewy panel: bieguny holendrowania i spirali przy wybranym KK (krzyżyki) i ścieżka przy KK od 0 do 1 (kropki). Prawy panel: odpowiedź na podmuch dający ślizg 5∘5^\circ i wychylenie steru kierunku. Przełącz filtr washout i programowanie wzmocnienia, zmieniaj prędkość lotu.

Zadania · część 1

3 zadania · ok. 20 min

  1. Przy jakim KK efektywne NrN_r jest trzy razy większe co do modułu niż bez tłumika (Nr=−1,04N_r = -1{,}04/s, Nδr=−5,69N_{\delta r} = -5{,}69/s²)?

    Odpowiedź

    −1,04−5,69K=−3,12-1{,}04 - 5{,}69K = -3{,}12, więc K=0,37K = 0{,}37.

  2. Oszacuj z modelu dwustanowego holendrowania (lekcja 7), jakie KK daje ζ=0,5\zeta = 0{,}5: 2ζω≈−(Yβ/U+Nr+NδrK)2\zeta\omega \approx -(Y_\beta/U + N_r + N_{\delta r}K), ω≈2,8\omega \approx 2{,}8 rad/s, Yβ/U=−0,159Y_\beta/U = -0{,}159.

    Odpowiedź

    2⋅0,5⋅2,8=2,8=0,159+1,04+5,69K2\cdot 0{,}5\cdot 2{,}8 = 2{,}8 = 0{,}159 + 1{,}04 + 5{,}69K, więc K=0,28K = 0{,}28. Pełny model: K=0,34K = 0{,}34 (przykład 6.1). Model dwustanowy zawyża tłumienie, więc zaniża potrzebne wzmocnienie.

  3. Dlaczego tłumik odchylenia bez filtru stabilizuje też spiralę (w symulacji −0,019→−0,15-0{,}019 \to -0{,}15/s przy K=0,5K = 0{,}5)? Skorzystaj z warunku LβNr−LrNβ>0L_\beta N_r - L_rN_\beta > 0.

    Odpowiedź

    Tłumik zwiększa ∣Nr∣|N_r|, a Lβ<0L_\beta < 0, więc LβNrL_\beta N_r rośnie: warunek jest spełniony z większym zapasem. Ster kierunku zmienia też LrL_r (w naszym modelu z 2,65 do 3,81/s, bo Lδr>0L_{\delta r} > 0), co działa w przeciwną stronę, ale słabiej: LβNr−LrNβL_\beta N_r - L_rN_\beta rośnie i tak z 13 do 95. Ta „korzyść” to druga strona problemu z części 2: tłumik przeciwdziała każdej prędkości odchylania, także tej potrzebnej w zakręcie.

2 Filtr washout: tłumik, który nie walczy z zakrętem

W ustalonym zakręcie skoordynowanym prędkość odchylania wynosi r=(gtan⁡φ/U)cos⁡φr = (g\tan\varphi/U)\cos\varphi (lekcja 7 Ruch boczny: przechylenie, spirala i holendrowanie), przy 30∘30^\circ i 55 m/s ok. 5∘5^\circ/s. Tłumik bez filtru widzi ją jako zakłócenie i wychyla ster kierunku przeciwnie do zakrętu. Odpowiedź na pytanie z prologu: samolot leci z ciągłym ślizgiem, a pilot (albo autopilot) musi tę walkę przeważyć. W symulacji 2: po 40 s zakrętu tłumik z K=0,5K = 0{,}5 trzyma ster kierunku na 2,4∘2{,}4^\circ przeciw zakrętowi, ślizg wynosi 2,9∘2{,}9^\circ, a przechylenie nie dochodzi do zadanego.

Rozwiązanie: tłumik ma reagować na zmiany prędkości odchylania o częstości holendrowania, a nie na jej wartość stałą. Wystarczy filtr górnoprzepustowy pierwszego rzędu, nazywany w automatyce lotniczej washout („wymywanie”):

Tłumik odchylenia z filtrem washout

δr=K τsτs+1 rczylix˙w=r−xwτ,δr=K(r−xw).\delta_r = K\,\frac{\tau s}{\tau s + 1}\,r\qquad\text{czyli}\qquad \dot x_w = \frac{r - x_w}{\tau},\quad \delta_r = K(r - x_w) .

Przy stałym rr filtr daje zero; przy ω≫1/τ\omega \gg 1/\tau przepuszcza sygnał prawie bez zmian. Typowo τ=1\tau = 1–4 s.

Dobór τ\tau to kompromis. Przy częstości holendrowania 3 rad/s filtr z τ=2\tau = 2 s przepuszcza 99% amplitudy z przesunięciem fazy 9∘9^\circ, więc tłumienie się nie pogarsza (ζ=0,70\zeta = 0{,}70 przy K=0,5K = 0{,}5), choć filtr dokłada do pętli własny biegun i częstość holendrowania spada z 3,03 do 2,67 rad/s. Za krótkie τ\tau (0,5 s) „wymywa” już samo holendrowanie (83% amplitudy, 34∘34^\circ fazy): tłumienie spada do 0,39. Za długie τ\tau zbliża nas do tłumika bez filtru: wejście w zakręt trwa dłużej, bo filtr długo „pamięta” początkowe rr. Pallett opisuje to samo w urządzeniach: sygnał żyroskopu prędkościowego przechodzi przez filtr, który przepuszcza tylko krótkotrwałe zmiany, a do koordynacji zakrętu używa się osobnego sygnału o stałej prędkości odchylania.

Symulacja 2 Wejście w zakręt prostym regulatorem przechylenia (lotki) i trzy wersje tłumika odchylenia. Górny wykres: przechylenie (linia przerywana: zadane). Dolny: ślizg, wychylenie steru kierunku i prędkość odchylania. Porównaj ster i ślizg po 40 s.

Zadania · część 2

3 zadania · ok. 25 min

  1. Policz wzmocnienie i przesunięcie fazy filtru τsτs+1\frac{\tau s}{\tau s + 1} przy ω=3\omega = 3 rad/s dla τ=0,5\tau = 0{,}5 s i τ=2\tau = 2 s.

    Odpowiedź

    ∣H∣=ωτ1+ω2τ2|H| = \frac{\omega\tau}{\sqrt{1 + \omega^2\tau^2}}, faza 90∘−arctan⁡ωτ90^\circ - \arctan\omega\tau. τ=0,5\tau = 0{,}5: 1,5/3,25=0,831{,}5/\sqrt{3{,}25} = 0{,}83, 34∘34^\circ. τ=2\tau = 2: 6/37=0,996/\sqrt{37} = 0{,}99, 9,5∘9{,}5^\circ.

  2. Tłumik bez filtru, K=0,5K = 0{,}5, ustalony zakręt z r=5∘r = 5^\circ/s. Jakie wychylenie steru „przeciw zakrętowi” daje tłumik? Porównaj z wynikiem symulacji (2,4∘2{,}4^\circ przy r=4,8∘r = 4{,}8^\circ/s).

    Odpowiedź

    δr=0,5⋅5∘=2,5∘\delta_r = 0{,}5\cdot 5^\circ = 2{,}5^\circ. W symulacji samolot ustala się przy trochę mniejszym rr (4,8∘4{,}8^\circ/s), bo ślizg i mniejsze przechylenie zmieniają równowagę: 0,5⋅4,8=2,4∘0{,}5\cdot 4{,}8 = 2{,}4^\circ ✓.

  3. Skojarz filtr washout z parą S+T=1S + T = 1 z lekcji 2.5 Zakłócenia, szum i wrażliwość: S, T i kształtowanie pętli i z filtrem komplementarnym z lekcji 4.5 automatyki. Jakie dwa filtry sumują się do jedności w układzie Palletta, w którym tłumik bierze rr przez filtr górnoprzepustowy, a koordynacja zakrętu stałe rr przez dolnoprzepustowy.

    Odpowiedź

    τsτs+1+1τs+1=1\frac{\tau s}{\tau s + 1} + \frac{1}{\tau s + 1} = 1, jak w filtrze komplementarnym z lekcji 4.5 automatyki. Część szybka (holendrowanie) trafia do tłumika, część wolna (ustalony zakręt) do koordynacji zakrętu, która wychyla ster w stronę zakrętu. Każde pasmo częstości ma swoje zadanie.

3 Wzmocnienie zależne od prędkości

Skuteczność steru rośnie z ciśnieniem dynamicznym, czyli jak V2V^2. Pallett pisze, że wzmocnienie wzmacniacza serwa jest programowane w funkcji prędkości: przy większej prędkości potrzeba mniejszego wychylenia. To jest programowanie wzmocnienia z lekcji 3.3 Nieliniowa rzeczywistość: nasycenie, windup, feedforward i gain scheduling. Ale co dokładnie ma być stałe?

Jeśli chcemy stałego tłumienia holendrowania, potrzebny przyrost 2ζω2\zeta\omega rośnie razem z częstością holendrowania, a ta rośnie mniej więcej jak VV. Przyrost ∣Nr∣|N_r| od tłumika to ∣Nδr∣K∝V2K|N_{\delta r}|K \propto V^2K. Żeby był proporcjonalny do VV, potrzeba K∝1/VK \propto 1/V. W modelu:

prędkośćbez tłumikaK=0,5K = 0{,}5 stałeK=0,5⋅55/VK = 0{,}5\cdot 55/V
35 m/sζ=0,17\zeta = 0{,}17ζ=0,44\zeta = 0{,}44ζ=0,61\zeta = 0{,}61 (K=0,79K = 0{,}79)
55 m/sζ=0,19\zeta = 0{,}19ζ=0,65\zeta = 0{,}65ζ=0,65\zeta = 0{,}65
75 m/sζ=0,19\zeta = 0{,}19ζ=0,84\zeta = 0{,}84ζ=0,66\zeta = 0{,}66 (K=0,37K = 0{,}37)

W tabeli pochodne zależne od CLC_L (CnpC_{n_p}, ClrC_{l_r}) wzięto z lotu przy 55 m/s; symulacja 1 przelicza CLC_L dla każdej prędkości, więc pokazuje nieco inne liczby (np. przy 35 m/s bez tłumika ζ=0,21\zeta = 0{,}21 i ω=2,28\omega = 2{,}28 zamiast 0,170{,}17 i 2,032{,}03). Wniosek o programowaniu K∝1/VK \propto 1/V pozostaje ten sam.

Ze stałym KK samolot jest przy małej prędkości słabo tłumiony, a przy dużej przetłumiony i ospały. Z K∝1/VK \propto 1/V tłumienie jest prawie stałe. Program K∝1/qˉK \propto 1/\bar q (stałe wychylenie momentu) dałby przy małej prędkości za dużo, przy dużej za mało. Wniosek ogólny: zanim zaprogramujesz wzmocnienie, ustal, która wielkość zamkniętej pętli ma być stała.

Zadania · część 3

3 zadania · ok. 20 min

  1. Nδr=−5,69N_{\delta r} = -5{,}69/s² przy 55 m/s. Ile wynosi przy 35 i 75 m/s?

    Odpowiedź

    Skaluje się jak V2V^2: −5,69⋅(35/55)2=−2,30-5{,}69\cdot(35/55)^2 = -2{,}30 i −5,69⋅(75/55)2=−10,6-5{,}69\cdot(75/55)^2 = -10{,}6/s².

  2. Pokaż, że przy K∝1/VK \propto 1/V iloczyn NδrKN_{\delta r}K rośnie jak VV, i wyjaśnij, dlaczego to daje stałe ζ\zeta, jeśli ω∝V\omega \propto V.

    Odpowiedź

    NδrK∝V2⋅V−1=VN_{\delta r}K \propto V^2\cdot V^{-1} = V. Przyrost 2ζω2\zeta\omega od tłumika jest równy ∣NδrK∣|N_{\delta r}K| (w modelu dwustanowym), więc Δζ∝V/ω=\Delta\zeta \propto V/\omega = const. Własne tłumienie samolotu (Nr∝VN_r \propto V, Yβ/U∝VY_\beta/U \propto V) skaluje się tak samo, więc całe ζ\zeta jest prawie stałe.

  3. Skąd tłumik ma wiedzieć o prędkości? Co się stanie, jeśli czujnik prędkości ulegnie awarii i pokaże zero?

    Odpowiedź

    Z centrali aerometrycznej (ciśnienie całkowite i statyczne, lekcja Dział 10: Płyny: ciśnienie, przepływ, lepkość i siły aerodynamiczne, przykład 10.2). Przy wskazaniu zero program dałby wzmocnienie dążące do nieskończoności: tłumik nasyciłby ster i mógł wzbudzić oscylacje. Dlatego wzmocnienie ogranicza się od góry, prędkość bierze się z kilku źródeł, a przy rozbieżności tłumik się wyłącza (Pallett: logika monitorująca i lampka awarii).

4 Tłumik pochylenia

Ten sam pomysł w ruchu podłużnym: δe=Kqq\delta_e = K_qq. Zmienia efektywne MqM_q o MδKqM_\delta K_q i trochę ZqZ_q. W modelu z lekcji 5 (przykład 3.3 z naszym Mδ=−25M_\delta = -25/s²):

  • Kq=0,05K_q = 0{,}05 s: krótkookresowa ζ=0,69→0,79\zeta = 0{,}69 \to 0{,}79, ω=3,54→3,87\omega = 3{,}54 \to 3{,}87 rad/s;
  • Kq=0,1K_q = 0{,}1 s: ζ=0,89\zeta = 0{,}89, ω=4,17\omega = 4{,}17 rad/s; przy 0,20{,}2 s krótkookresowa jest już aperiodyczna.

Samolot lekki tego nie potrzebuje. Inaczej jest wysoko: przy tej samej prędkości przyrządowej (stałe qˉ\bar q) prędkość rzeczywista rośnie, a pochodne tłumiące maleją jak 1/U1/U (zadanie 8.1 z lekcji 5). Ten sam samolot przeliczony rachunkowo na 9000 m miałby krótkookresową z ζ=0,47\zeta = 0{,}47, a z Kq=0,1K_q = 0{,}1 s znowu 0,76. Samoloty odrzutowe latające wysoko mają więc tłumik pochylenia, zwykle zintegrowany z układem sterowania.

Tłumik pochylenia też walczy z manewrem: w ustalonym zakręcie i w wyrwaniu q≠0q \ne 0 przez długi czas. Rozwiązania są te same: washout albo sprzężenie od qq w pętli wewnętrznej regulatora, który za polecenie bierze zadaną prędkość pochylania lub przeciążenie (lekcja 9 Autopilot: pętle wewnętrzne i zewnętrzne).

Zadania · część 4

3 zadania · ok. 20 min

  1. O ile zmienia się efektywne Mq/JyM_q/J_y przy Kq=0,1K_q = 0{,}1 s i Mδ/Jy=−25M_\delta/J_y = -25/s²?

    Odpowiedź

    MδKq=−2,5M_\delta K_q = -2{,}5/s: z −2-2 do −4,5-4{,}5/s, ponad dwa razy więcej.

  2. Z modelu dwustanowego (lekcja 6) policz ωs\omega_s i ζs\zeta_s z tłumikiem Kq=0,1K_q = 0{,}1 s (pomiń zmianę ZqZ_q).

    Odpowiedź

    Mq/Jy=−4,5M_q/J_y = -4{,}5: ωs2=1,964⋅4,5+8,6=17,44\omega_s^2 = 1{,}964\cdot 4{,}5 + 8{,}6 = 17{,}44, ωs=4,18\omega_s = 4{,}18; 2ζω=1,964+4,5+0,9=7,362\zeta\omega = 1{,}964 + 4{,}5 + 0{,}9 = 7{,}36, ζ=0,88\zeta = 0{,}88. Pełny model: 4,17 i 0,885 ✓. Tłumik pochylenia zwiększa też częstość, bo MqM_q wchodzi do ωs2\omega_s^2 przez człon ZαMqZ_\alpha M_q.

  3. Dlaczego tłumienie krótkookresowej maleje z wysokością przy stałej prędkości przyrządowej, a częstość prawie nie?

    Odpowiedź

    Częstość zależy głównie od Mα/Jy∝qˉM_\alpha/J_y \propto \bar q, a ono jest stałe. Tłumienie zależy od MqM_q, Mα˙M_{\dot\alpha}, Zα/(mU)Z_\alpha/(mU), a każde z nich ma UU w mianowniku: przy 9000 m UU jest 1,6 razy większe. W modelu ωs\omega_s: 3,54→3,183{,}54 \to 3{,}18 rad/s, ζs\zeta_s: 0,69→0,470{,}69 \to 0{,}47.

5 Fly-by-wire i niestateczność statyczna

Skoro regulator może dodać tłumienie, może dodać też sztywność. Sprzężenie od kąta natarcia, δe=Kαα\delta_e = K_\alpha\alpha, zmienia efektywne CmαC_{m\alpha} o CmδKαC_{m\delta}K_\alpha: działa jak przesunięcie punktu neutralnego do tyłu. Samolot może więc mieć środek ciężkości za punktem neutralnym geometrycznym i nadal latać stabilnie. Zyski: mniejsze usterzenie, mniejsza siła w dół na usterzeniu, czyli mniejszy opór trymu, lepsza zwrotność. Takie samoloty mają sterowanie elektryczne (ang. fly-by-wire, FBW): drążek nie rusza steru, tylko wysyła polecenie do komputera, który wychyla ster tak, żeby samolot zachowywał się jak stateczny. Pallett w rozdziale 11 opisuje FBW razem z redundancją: kilka niezależnych kanałów, czujników i źródeł zasilania, bo awaria oznacza utratę samolotu.

W modelu z lekcji 4 przy h=0,55h = 0{,}55 (4% cięciwy za punktem neutralnym) samolot ma pierwiastek +0,157+0{,}157/s, czyli podwojenie w 4,4 s. Sprzężenie Kα=0,4K_\alpha = 0{,}4 przesuwa wszystkie pierwiastki na lewą półpłaszczyznę; efektywne Cmα=0,217+(−0,995)(0,4)=−0,18C_{m\alpha} = 0{,}217 + (-0{,}995)(0{,}4) = -0{,}18, czyli równoważny zapas 3,4% cięciwy. Przy samym Kα=0,4K_\alpha = 0{,}4 pierwiastki krótkookresowe są już rzeczywiste (−7,5-7{,}5 i −2,4-2{,}4/s); dodatkowe Kq=0,1K_q = 0{,}1 s przesuwa szybszy z nich do −9,8-9{,}8/s, czyli głównie przyspiesza odpowiedź.

Symulacja 3 Samolot z lekcji 4 z regulatorem δ=Kαα+Kqq\delta = K_\alpha\alpha + K_qq. Lewy panel: bieguny bez regulatora (kółka) i z regulatorem (krzyżyki). Prawy panel: kąt natarcia i wychylenie steru po podmuchu +2∘+2^\circ. Przesuń środek ciężkości za punkt neutralny i znajdź najmniejsze KαK_\alpha, przy którym samolot jest stateczny.

Zadania · część 5

2 zadania · ok. 15 min

  1. Przy h=0,55h = 0{,}55: Cmα=+0,217C_{m\alpha} = +0{,}217, Cmδ=−0,995C_{m\delta} = -0{,}995. Jakie KαK_\alpha daje równoważny zapas 5% cięciwy (a=5,39a = 5{,}39)?

    Odpowiedź

    Potrzebne Cmα,ef=−5,39⋅0,05=−0,27C_{m\alpha,\text{ef}} = -5{,}39\cdot 0{,}05 = -0{,}27. 0,217−0,995Kα=−0,270{,}217 - 0{,}995K_\alpha = -0{,}27, więc Kα=0,49K_\alpha = 0{,}49 rad steru na rad kąta natarcia.

  2. Czujnik kąta natarcia zamarza i pokazuje stałą wartość. Co dzieje się z samolotem z części 5? Jak zabezpieczają się przed tym prawdziwe systemy?

    Odpowiedź

    Sprzężenie od α\alpha znika (stały sygnał to tylko stałe wychylenie steru), więc wraca niestateczność statyczna z czasem podwojenia kilku sekund, a do tego stałe wychylenie wytrąca samolot z trymu. Zabezpieczenia: kilka czujników kąta natarcia i porównywanie ich wskazań, estymacja kąta natarcia z innych sygnałów (przyspieszenie normalne, prędkość, pochylenie: filtr z lekcji 4.4 automatyki), przejście na prostsze prawo sterowania przy niezgodności.

6 Metoda i przykłady rozwiązane

Metoda: projekt tłumika

  1. Diagnoza: który mod jest za słabo tłumiony (wartości własne, wymagania jakości pilotażu z lekcji 6).
  2. Sprzężenie: prędkość kątowa tej osi do steru tej osi (r→δrr \to \delta_r, q→δeq \to \delta_e); znak tak, żeby zwiększyć tłumienie.
  3. Linie pierwiastkowe względem KK; wybór KK z zapasem.
  4. Washout z 1/τ1/\tau kilka razy mniejszym niż częstość modu; sprawdzenie manewru ustalonego.
  5. Program wzmocnienia w funkcji prędkości (i wysokości); sprawdzenie w rogach obwiedni lotu.

Przykład 6.1 Wzmocnienie tłumika dla ζ=0,5\zeta = 0{,}5

Dobierz KK tłumika odchylenia samolotu z lekcji 7, tak żeby holendrowanie miało ζ=0,5\zeta = 0{,}5, bez filtru i z washout τ=2\tau = 2 s.

  1. Szacunek z modelu dwustanowego: K≈0,28K \approx 0{,}28 (zadanie 1.2).
  2. Pełny model, bez filtru: ζ(K)\zeta(K) rośnie monotonicznie; ζ=0,5\zeta = 0{,}5 przy K=0,34K = 0{,}34 (bisekcja).
  3. Z washout τ=2\tau = 2 s: K=0,33K = 0{,}33; filtr przy 3 rad/s przepuszcza prawie wszystko.

Wynik: K≈0,33K \approx 0{,}33–0,340{,}34, czyli ok. 0,34∘0{,}34^\circ steru na 1∘1^\circ/s prędkości odchylania.

Sprawdzenie: przy holendrowaniu z amplitudą r=5∘r = 5^\circ/s ster wychyla się o ok. 1,7∘1{,}7^\circ, daleko od ograniczeń ✓.

Przykład 6.2 Zakręt z tłumikiem i bez washout

Regulator przechylenia wprowadza samolot w zakręt 30∘30^\circ. Porównaj stan po 40 s bez tłumika, z tłumikiem K=0,5K = 0{,}5 bez filtru i z filtrem τ=2\tau = 2 s.

  1. Bez tłumika: φ=29,8∘\varphi = 29{,}8^\circ, β=0,77∘\beta = 0{,}77^\circ, r=5,2∘r = 5{,}2^\circ/s, ster kierunku 0; największy ślizg w czasie wejścia 7,7∘7{,}7^\circ.
  2. Tłumik bez filtru: φ=28,4∘\varphi = 28{,}4^\circ, β=2,9∘\beta = 2{,}9^\circ, r=4,8∘r = 4{,}8^\circ/s, ster kierunku +2,4∘+2{,}4^\circ (przeciw zakrętowi).
  3. Tłumik z washout: stan ustalony jak bez tłumika (ster wraca do zera), największy ślizg 5,6∘5{,}6^\circ.

Wynik: washout zachowuje korzyść tłumika w czasie wejścia w zakręt i usuwa szkodę w zakręcie ustalonym.

Sprawdzenie: żaden wariant nie daje zerowego ślizgu w ustalonym zakręcie: do tego potrzebna jest koordynacja zakrętu (sygnał wychylający ster w stronę zakrętu albo pętla zerująca ślizg) ✓.

Przykład 6.3 Sztuczna stateczność statyczna

Samolot z lekcji 4 ze środkiem ciężkości w 55% cięciwy. Jakie sprzężenie od kąta natarcia czyni go statecznym?

  1. Bez regulatora: pierwiastki −8,93-8{,}93; −0,81-0{,}81; −0,30-0{,}30; +0,157+0{,}157/s.
  2. Kα=0,2K_\alpha = 0{,}2: największy pierwiastek +0,026+0{,}026, wciąż niestateczny (Cmα,ef=0,217−0,199=+0,018C_{m\alpha,\text{ef}} = 0{,}217 - 0{,}199 = +0{,}018).
  3. Kα=0,4K_\alpha = 0{,}4: pierwiastki −7,47-7{,}47; −2,44-2{,}44; −0,022±0,12i-0{,}022 \pm 0{,}12i: stateczny, z wolną parą jak fugoida.

Wynik: granica leży przy Kα≈0,22K_\alpha \approx 0{,}22, gdzie efektywne CmαC_{m\alpha} przechodzi przez zero; z zapasem bierzemy 0,40{,}4–0,50{,}5.

Sprawdzenie: granica stateczności dynamicznej pokrywa się z Cmα,ef=0C_{m\alpha,\text{ef}} = 0, tak jak w lekcji 4 pokrywała się z punktem neutralnym ✓.

7 Błędne intuicje, których trzeba się oduczyć

„Tłumik odchylenia to autopilot kierunku.”
Tłumik tylko gasi szybkie wahania; nie trzyma kursu i działa także przy pilotażu ręcznym. Autopilot to pętla zewnętrzna (lekcja 9).
„Większe wzmocnienie tłumika jest zawsze lepsze.”
Za duże daje mod aperiodyczny i ospałą reakcję, nasyca ster, a z opóźnieniami serwa może wzbudzić drgania. Poza tym bez washout walczy z manewrem.
„Washout osłabia tłumik.”
Przy dobrze dobranym τ\tau prawie nie zmienia tłumienia holendrowania (przepuszcza 99% sygnału), a usuwa walkę z zakrętem.
„Samolot niestateczny statycznie nie może latać.”
Może, jeśli regulator działa szybciej niż niestateczność rośnie i nigdy nie zawodzi. Ceną jest pełna zależność od elektroniki.
„Wzmocnienie programuje się odwrotnie do ciśnienia dynamicznego.”
To zależy od tego, co ma być stałe. Dla stałego tłumienia holendrowania w naszym modelu potrzeba K∝1/VK \propto 1/V, nie 1/qˉ1/\bar q.

8 Zadania do przemyślenia

Te zadania wymagają więcej myślenia niż liczenia.

Zadania do przemyślenia

5 zadań · ok. 1,5 godziny

  1. Opóźnienie serwa. Serwo steru kierunku ma transmitancję 1/(0,1s+1)1/(0{,}1s + 1). Jak to zmienia projekt tłumika i dlaczego zbyt duże KK może wtedy doprowadzić do drgań?

    Rozwiązanie
    1. Przy częstości holendrowania 3 rad/s serwo daje opóźnienie fazy arctan⁡0,3=17∘\arctan 0{,}3 = 17^\circ: część wychylenia przychodzi „za późno” i mniej tłumi. Przy większych częstościach (np. drgania sprężyste kadłuba) opóźnienie rośnie do 90∘90^\circ.
    2. Z dużym KK pętla zaczyna oddziaływać na szybsze mody (serwo, drgania konstrukcji), gdzie przesunięcie fazy może dać dodatnie sprzężenie: zapas fazy (lekcja 2.4 automatyki) maleje do zera. Dlatego tłumiki mają też filtry dolnoprzepustowe i ograniczone wzmocnienie.
  2. Tłumik szeregowy i równoległy. Pallett opisuje tłumik szeregowy: jego wychylenie dodaje się do wychylenia pilota bez ruchu pedałów. Jaka jest zaleta takiego rozwiązania i dlaczego jego autorytet (zakres wychylenia) ogranicza się do kilku stopni?

    Rozwiązanie
    1. Pilot nie czuje pracy tłumika i zachowuje pełną kontrolę nad sterem; awaria tłumika (np. wychylenie maksymalne) daje tylko ograniczone zaburzenie, które pilot łatwo pokona.
    2. Do gaszenia holendrowania wystarczy kilka stopni (przykład 6.1: 1,7∘1{,}7^\circ przy r=5∘r = 5^\circ/s), a ograniczony autorytet czyni awarię nieszkodliwą. Pełny autorytet mają dopiero układy FBW, ale wtedy z redundancją.
  3. Sprzężenie od ślizgu. Zamiast rr można wprowadzić do steru kierunku ślizg β\beta (z czujnika albo z akcelerometru poprzecznego, lekcja 7 Ruch boczny: przechylenie, spirala i holendrowanie, zadanie 8.4). Jakie dwa efekty daje δr=−Kββ\delta_r = -K_\beta\beta?

    Rozwiązanie
    1. Zwiększa efektywne NβN_\beta: sztuczna stateczność kierunkowa, wyższa częstość holendrowania. Tłumienia sam nie dodaje (trzeba rr albo β˙\dot\beta).
    2. Zmniejsza ślizg w ustalonym manewrze (a z członem całkującym go zeruje), czyli koordynuje zakręt automatycznie: tego nie zrobi tłumik od rr. Dlatego nowoczesne układy łączą oba sygnały.
  4. Wysokość jako druga zmienna programu. W części 4 tłumienie krótkookresowej spadało z wysokością przy stałej prędkości przyrządowej. Jak powinno zależeć KqK_q od wysokości i prędkości, żeby ζs\zeta_s było stałe?

    Rozwiązanie
    1. Przyrost 2ζω2\zeta\omega od tłumika to ∣Mδ∣Kq|M_\delta|K_q, a Mδ∝qˉM_\delta \propto \bar q. Częstość ωs\omega_s zależy głównie od qˉ\sqrt{\bar q}. Brakujące tłumienie rośnie z wysokością (bo własne Mq∝qˉ/UM_q \propto \bar q/U maleje).
    2. Z warunku 2ζsωs=c qˉ/U+k qˉ Kq2\zeta_s\omega_s = c\,\bar q/U + k\,\bar q\,K_q (własne tłumienie plus tłumik, cc i kk stałe samolotu): Kq=2ζsωs−c qˉ/Uk qˉK_q = \dfrac{2\zeta_s\omega_s - c\,\bar q/U}{k\,\bar q}. Przy stałym qˉ\bar q (stała prędkość przyrządowa) i rosnącej wysokości rośnie prędkość rzeczywista UU, więc KqK_q też musi rosnąć: program zależy od dwóch zmiennych, qˉ\bar q i UU (albo wysokości).
    3. W praktyce wzmocnienie programuje się dwuwymiarową tablicą od qˉ\bar q i liczby Macha (albo wysokości), wyznaczoną z linearyzacji w wielu punktach obwiedni, jak w lekcji 3.3 automatyki.
  5. Ile niestateczności? Gdyby środek ciężkości samolotu z części 5 przesunąć do 60% cięciwy (podwojenie w 2 s), a serwo miało opóźnienie 0,1 s, czy regulator z sprzężeniem od α\alpha i qq nadal by sobie poradził? Jak to sprawdzić?

    Rozwiązanie
    1. Nie z tymi samymi wzmocnieniami: przy h=0,60h = 0{,}60 efektywne CmαC_{m\alpha} rośnie do ok. +0,49+0{,}49, a Kα=0,4K_\alpha = 0{,}4, Kq=0,1K_q = 0{,}1 zostawia pierwiastek +0,066+0{,}066/s. Granica stateczności przesuwa się do Kα≈0,5K_\alpha \approx 0{,}5; z zapasem trzeba ok. 0,7. Po przestrojeniu niestateczny pierwiastek (+0,34+0{,}34/s) jest wciąż wolny w porównaniu z serwem (10/s), więc samo opóźnienie 0,1 s nie jest przeszkodą. Kłopot pojawia się, gdy pasmo pętli (potrzebne, żeby szybko przeciwdziałać niestateczności) zbliża się do pasma serwa i drgań konstrukcji.
    2. Sprawdzić zapasy fazy i amplitudy pętli z serwem (lekcja 2.4 automatyki); dla obiektów niestatecznych istnieją też dolne granice zapasu wzmocnienia (zbyt małe wzmocnienie też destabilizuje). Na koniec symulacja nieliniowa z ograniczeniami steru.

9 Sprawdź się

Odpowiedz w myślach, zanim rozwiniesz odpowiedź.

1 Jak działa tłumik odchylenia?

δr=Kr\delta_r = Kr: zwiększa efektywne ∣Nr∣|N_r|, więc tłumienie holendrowania.

2 Po co filtr washout?

Żeby tłumik nie reagował na stałą prędkość odchylania w ustalonym zakręcie; przepuszcza tylko zmiany szybsze niż 1/τ1/\tau.

3 Jak dobrać τ\tau filtru washout?

1/τ1/\tau kilka razy mniejsze niż częstość holendrowania (typowo τ=1\tau = 1–4 s); za krótkie wymywa holendrowanie, za długie walczy z zakrętem.

4 Jak programować wzmocnienie tłumika z prędkością?

Tak, żeby stała była wielkość docelowa; dla stałego ζ\zeta holendrowania w naszym modelu K∝1/VK \propto 1/V.

5 Co robi sprzężenie od kąta natarcia do steru wysokości?

Zmienia efektywne CmαC_{m\alpha} o CmδKαC_{m\delta}K_\alpha: sztuczna stateczność statyczna, przesunięcie punktu neutralnego do tyłu.

6 Czym SAS różni się od autopilota?

SAS (tłumiki, sztuczna stateczność) to pętla wewnętrzna poprawiająca mody; autopilot to pętla zewnętrzna wydająca polecenia (wysokość, kurs, prędkość).

10 Karta wzorów i co dalej

Tłumiki

δr=Kτsτs+1r\delta_r = K\frac{\tau s}{\tau s + 1}r: Nr→Nr+NδrKN_r \to N_r + N_{\delta r}K (w paśmie holendrowania)

δe=Kqq\delta_e = K_qq: Mq→Mq+MδKqM_q \to M_q + M_\delta K_q

washout: ∣H(iω)∣=ωτ1+ω2τ2|H(i\omega)| = \frac{\omega\tau}{\sqrt{1 + \omega^2\tau^2}}, faza 90∘−arctan⁡ωτ90^\circ - \arctan\omega\tau

stałe ζ\zeta: K∝ω/∣Nδr∣∝1/VK \propto \omega/|N_{\delta r}| \propto 1/V

Sztuczna stateczność

δe=Kαα\delta_e = K_\alpha\alpha: Cmα,ef=Cmα+CmδKαC_{m\alpha,\text{ef}} = C_{m\alpha} + C_{m\delta}K_\alpha

równoważny zapas Kn,ef=−Cmα,ef/aK_{n,\text{ef}} = -C_{m\alpha,\text{ef}}/a

nasz model: h=0,55h = 0{,}55, Kα=0,4K_\alpha = 0{,}4: zapas 3,4%3{,}4\%

Do przeczytania

  • Obowiązkowo: E. H. J. Pallett, S. Coyle, Automatic Flight Control, rozdz. 7 (s. 222–224, układ tłumika odchylenia, rys. 7.13) i rozdz. 3 (s. 122–124, koordynacja zakrętu i programowanie wzmocnienia).
  • Dla chętnych: Pallett, rozdz. 11 (fly-by-wire); A. Tewari, rozdz. 3.6 i 3.9 (SAS i autopilot w ujęciu przestrzeni stanów).

Pętla wewnętrzna jest gotowa: samolot ma dobrze tłumione mody, także wtedy, gdy natura dała mu za mało. Teraz możemy mu coś kazać. Następna lekcja buduje autopilota tak, jak opisuje go Pallett: wokół pętli trzymania pochylenia i przechylenia, z pętlami zewnętrznymi wysokości, prędkości pionowej, kursu i prędkości, z synchronizacją przy włączaniu i automatycznym trymem.

Korzystają z tej lekcji 9 Autopilot: pętle wewnętrzne i zewnętrzne