MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Mgr inż. Witold BUŻANTOWICZ DOI: 10.17814/mechanik.2015.7.217
Mgr inż. Jakub MIERNIK
Dr hab. inż. Jan PIETRASIECSKI, prof. WAT
Wojskowa Akademia Techniczna
MODELOWANIE I APROKSYMACJA
FUNKCJI PRZENOSZENIA MASZYNEK STEROWYCH
RAKIETY PRZECIWLOTNICZEJ
Streszczenie: W artykule przedstawiono procedurę numerycznego
wyznaczania charakterystyk dynamicznych maszynek sterowych,
stanowiących element wyposażenia stanowiska laboratoryjnego do badań
aparatury pokładowej rakiet przeciwlotniczych. Występujące
w rozpatrywanym silniku DC nieliniowości, spowodowane zastosowaniem
wewnętrznego korektora prędkości kątowej, sprawiają, że znalezienie
adekwatnej aproksymacji jego charakterystyk wymaga specyficznego
podejścia. Do ich wyznaczenia wykorzystano metody iteracyjne,
aproksymujące odpowiedzi badanego układu funkcjami sklejanymi.
Otrzymana w ten sposób funkcja przenoszenia układu wykonawczego
zaimplementowana zostanie w modelu numerycznym projektowanego
układu pilota automatycznego.
FIN ACTUATOR TRANSFER FUNCTION OF
THE ANTI-AIRCRAFT MISSILE:
MODELING AND APPROXIMATION
Abstract: This paper presents the procedure for numerical determining
dynamic characteristics describing the actuators of short-range anti-aircraft
missile, which are the elements of research laboratory stand. Internal
compensator of angular velocity, included in considered non-linear DC
servo, requires a specific approach to find an adequate solution. Therefore
an iterative method is used, which approximates the object s step and
impulse responses by spline functions. The transfer function obtained in this
way will be implemented in the numerical model of the proposed missile
autopilot system.
Słowa kluczowe: teoria sterowania, układy wykonawcze, rakieta
przeciwlotnicza
Keywords: control theory, actuators, anti-aircraft missile
1. WPROWADZENIE
Właściwy dobór współczynników transmitancji członów automatyki ma dla procesu
modelowania kluczowe znaczenie, stanowi bowiem podstawowy warunek zaprojektowania
poprawnego układu regulacji [1, 3]. Często zdarza się jednak, że parametry charakteryzujące
właściwości elementów wykonawczych nie są dostatecznie znane i wyznaczenie na ich
podstawie funkcji przenoszenia jest niemożliwe. Ponadto uzyskiwane dla niektórych
83
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
rodzajów obiektów dynamicznych wyższych rzędów rozwiązania analityczne mają
ograniczone znaczenie praktyczne, bowiem mogą być trudne do wykorzystania ze względu na
złożoną formę matematyczną. We wskazanych przypadkach alternatywne rozwiązanie
stanowi doświadczalne wyznaczenie i aproksymacja odpowiedzi badanych układów.
Celem przedstawionego w niniejszym referacie eksperymentu było wyznaczenie
charakterystyk czasowo-częstotliwościowych maszynek sterowych, przewidywanych do
zastosowania jako element stanowiska laboratoryjnego do badania aparatury pokładowej
rakiet przeciwlotniczych.
Najczęściej spotykany schemat układu stabilizacji rakiety przeciwlotniczej ilustruje rys. 1.
W rozpatrywanym przypadku obiektem badań jest maszynka sterowa, spełniająca funkcję
układu wykonawczego. Poszukiwana jest jej odpowiedz impulsowa.
Rys. 1. Schemat funkcjonalny pilota automatycznego rakiety przeciwlotniczej [2]
2. EKSPERYMENT
Stanowisko badawcze zestawiono zgodnie ze schematem przedstawionym na rys. 2.
Identyfikowany układ wykonawczy jest zbudowany w układzie komutatorowego silnika
prądu stałego, wyposażonego w wewnętrzny korektor prędkości kątowej. Silnik sterowany
jest cyfrowo, przy czym wymagane jest wstępne zaprogramowanie jego położeń. Maszynka
sterowa współpracuje ze sterownikiem mikroprocesorowym. Stanowisko tworzą ponadto:
zasilacz 24 V oraz komputer do wizualizacji i archiwizacji danych pomiarowych. Wchodzące
w skład zestawu oprogramowanie zapewnia możliwość rejestracji parametrów pracy silnika,
m.in. prędkości i położeń kątowych (zarówno zadanych, jak i uzyskanych). Aplikacja
komunikuje się z układem wykonawczym poprzez konwerter USB/CAN.
84
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Rys. 2. Schemat układu pomiarowego
Podczas badań rejestrowano odpowiedzi układu na wymuszenia skokowe dla różnych
wzmocnień regulatora wewnętrznego (rys. 3).
Rys. 3. Przykładowa odpowiedz badanego układu wykonawczego na wymuszenie skokowe
85
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
3. WYZNACZANIE CHARAKTERYSTYK BADANEGO UKAADU
Osobliwością rozpatrywanego elementu wykonawczego są nieliniowości wprowadzane przez
cyfrowy korektor wewnętrzny, szczególnie widoczne w końcowej fazie ustalania położenia.
Stąd znalezienie adekwatnej aproksymacji wymaga specyficznego podejścia.
Do określenia postaci transmitancji zastępczego, zlinearyzowanego modelu silnika
i wyznaczenia jej współczynników wykorzystano numeryczne metody obliczeniowe.
W pierwszej kolejności dla każdego z otrzymanych eksperymentalnie przebiegów określono
przedziały aproksymacji. Przeprowadzone badania symulacyjne wykazały, że
w rozpatrywanym przypadku wystarczające jest przyjęcie dwóch przedziałów, 1 = ( ", s)
oraz 2 = (s, + "). Z tak zdefiniowanymi 1 i 2 związano funkcje aproksymujące,
odpowiednio G1(s) dla 1 i G2(s) dla 2, oraz przyjęto punkt ich sklejenia P(ts, s), por. rys. 4.
Następnie ze względu na podobieństwo uzyskanych eksperymentalnie przebiegów do
charakterystyk elementów inercyjnych i przewidywaną konieczność uwzględniania opóznień
założono ogólne równania funkcji przenoszenia w przedziałach postaci:
y s
( ) k
G s = = e-sT
( )
N
(1)
u s
( )
s +1
)
(Tn
n=1
Skorzystano przy tym z faktu, że układy inercyjne wyższych rzędów mogą być otrzymywane
poprzez szeregowe połączenie członów inercyjnych rzędu pierwszego.
Rys. 4. Przyjęty sposób doboru przedziałów aproksymacji
86
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Dyskretny element zastępczy układu o inercyjności pierwszego rzędu został opisany za
pomocą równania różnicowego o zmiennych współczynnikach w postaci funkcji języka
skryptowego MathWorks"! MATLAB. Jako parametry wejściowe funkcji przyjęto: czas
próbkowania Tp, współczynnik wzmocnienia k, stałą czasową elementu T oraz sygnał
wymuszający u. W wyniku wykonania programu otrzymuje się sygnał odpowiedzi y
symulowanego obiektu inercyjnego w funkcji zadanych parametrów wejściowych.
Przekazanie sygnału wyjściowego y jako sygnału wymuszenia u w kolejnym wywołaniu
funkcji wraz z nowymi wartościami pozostałych parametrów pozwala uzyskać człon
inercyjny drugiego rzędu, kolejne wywołania rzędu trzeciego, czwartego itd.
Współczynniki funkcji aproksymującej poszukiwane są metodą iteracyjną. Dla zadanego
rzędu n elementu inercyjnego wyznacza w n pętlach przebiegi odpowiedzi na skok
jednostkowy u dla wszystkich wartości parametrów (stałych czasowych i wzmocnienia).
Wektor sygnałów odpowiedzi Y = { y1, y2, & , yI } układu porównywany jest następnie
z sygnałem odniesienia y0, tj. odpowiedzią układu rzeczywistego poszukuje się przebiegu yi,
dla którego błąd Ńi jest najmniejszy:
I L
2
1
Jmin = min Ł Ji = yi l - y0 l
( [ ] [ ] (2)
)
ż
UJi
L
i=1 l=1
gdzie: L liczba próbek sygnału.
W wyniku wykonania program zwraca następujące wartości: poszukiwany rząd układu n,
wzmocnienie k, wektor stałych czasowych T = { T1, T2, & , Tn } oraz wskaznik jakości
w postaci minimum błędu średniokwadratowego Ńmin aproksymacji.
Opisana procedura wykonywana jest dla dwóch funkcji aproksymujących, tj. G1(s)
w przedziale 1 oraz G2(s) w przedziale 2. Wartości czasów opóznień oraz położenie punktu
sklejania funkcji P są wyznaczane poprzez znalezienie minimum błędu średniokwadratowego
przebiegu jako całości. Spełnione być muszą przy tym dwa warunki nadrzędny, narzucający
brak przeregulowań, oraz dodatkowy, nakazujący osiągnięcie przez funkcję aproksymującą
wartości ust w czasie tr, równym trzykrotności stałej czasowej T zastępczego elementu
inercyjnego pierwszego rzędu:
fT
p
(3)
tr = fT = 3T T = Ł ymax = y Ą = dust
( )
p
3
gdzie: T stała czasowa zastępczego elementu inercyjnego I rzędu,
Tp czas próbkowania,
f liczba cykli obliczeniowych,
dust wartość sygnału y(t) w stanie ustalonym.
Dzięki tak sformułowanemu kryterium możliwe jest zarówno skuteczne odtworzenie
dynamiki odpowiedzi modelowanego układu w pierwszej fazie narastania sygnału, jak
i zapewnienie osiągnięcia zadanej wartości sygnału ust w założonym czasie i bez
przeregulowań.
87
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
W wyniku wykonania symulacji dla zbioru kryteriów (2)-(3) otrzymano krzywą sklejaną y(t),
związaną z sygnałem wejściowym u(t) transmitancją postaci:
k1
G s =
1
e-sT y t < ds
( ) ( )
1
4
s +1
)
(T1i
y s
( )
i=1
G s = =
( ) (4)
u s
( )
2
G2 s = k2 e-sT y t ł ds
( ) ( )
5
T2is +1
( )
i=1
gdzie: ds wyznaczona symulacyjnie wartość rzędnej punktu sklejenia funkcji P.
Uzyskaną na podstawie równania (4) aproksymację odpowiedzi badanego układu na skok
jednostkowy przedstawiono na rys. 5.
Wyznaczenie wartości współczynników zależności (4) pozwala znalezć postać odpowiedzi
impulsowej h(t), umożliwiającej z kolei określenie przebiegu wyjściowego y(t) poprzez całkę
splotową:
+Ą
y t = x t h t -t dt
( ) ( ) ( ) (5)
-Ą
W przypadku układów dyskretnych wyrażenie (5) przyjmuje postać:
Ą
y n = x k h n - k
[ ] [ ] [ ] (6)
k=0
88
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Rys. 5. Wyznaczona metodą iteracyjną aproksymacja rzeczywistego przebiegu
odpowiedzi badanego silnika na wymuszenie skokowe
Uwzględniając powyższe, na podstawie wyrażenia (4) wyznaczono numerycznie odpowiedzi
impulsowe h1[n] i h2[n] układów opisanych transmitancjami dyskretnymi:
k1
1
G1 z = Z G1 s = Z e-sT ż
( ) { ( )} (7a)
4
s +1
)
(T1i
i=1
k2
2
G2 z = Z G2 s = Z e-sT ż
( ) { ( )} (7b)
5
s +1
)
(T2i
i=1
i w oparciu o nie (por. rys. 6) złożono funkcję wynikową h[n] postaci:
h1 n nT < ts
[ ]
p
h n =
[ ] (8)
h n nT ł ts
[ ]
2 p
ts
gdzie: wyznaczona symulacyjnie wartość odciętej punktu sklejenia funkcji P.
89
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Dla oszacowania jakości aproksymacji wygenerowano testowy sygnał wyjściowy y[n] jako
splot funkcji skoku jednostkowego i odpowiedzi impulsowej h[n] określonej wyrażeniem (8):
Ą
y n =
[ ] (9)
1[k]h[n - k]
k =0
uzyskując wysoką zgodność z przebiegiem aproksymowanym zależnością (4), przy błędzie
średniokwadratowym z wartości równania (9) odniesionych do wartości wyrażenia (4)
równym Ń = 3,710-6. Obliczoną na podstawie zależności (8) charakterystykę amplitudowo-
-fazową zweryfikowano symulacyjnie, podając na wejście modelu sygnał sinusoidalny u[n]
o częstościach z przedziału od 0,01 do 35 rad/s. Po ustaleniu się odpowiedzi odczytywano
wzmocnienie oraz przesunięcie fazowe sygnału wyjściowego y[n], potwierdzając wyniki dane
na rys. 7.
Rys. 6. Odpowiedz impulsowa modelowanego układu
90
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
Rys. 7. Logarytmiczna charakterystyka amplitudowo-fazowa badanego silnika
91
MECHANIK 7/2015
XIX Międzynarodowa Szkoła Komputerowego Wspomagania Projektowania, Wytwarzania i Eksploatacji
4. PODSUMOWANIE
W artykule zaproponowano sposób aproksymacji charakterystyk nieliniowego silnika DC
funkcją sklejaną z liniowych funkcji przenoszenia, obejmujący etapy: (a) rejestracji
charakterystyki skokowej rozpatrywanego obiektu, (b) aproksymacji otrzymanych wyników
funkcjami sklejanymi założonej postaci, (c) symulacyjnego badania odpowiedzi układu na
wymuszenie impulsowe oraz (d) wyznaczania charakterystyki amplitudowo-fazowej.
Na podstawie przeprowadzonych symulacji można przedstawić następujące wnioski:
modelowanie nieliniowych elementów wykonawczych za pomocą członów
inercyjnych wyższych rzędów i funkcji sklejanych zapewnia odwzorowanie ich
dynamiki,
opracowany sposób umożliwia określanie wartości parametrów funkcji
aproksymujących badany układ inercyjny z uwzględnieniem stałych czasowych
opóznień.
Przedstawione w referacie oprogramowanie symulacyjne charakteryzuje się strukturą otwartą,
stąd możliwa jest jego dalsza rozbudowa w celu pozyskania efektywnego narzędzia do
identyfikacji obiektów o innych właściwościach dynamicznych.
LITERATURA
[1] Cannon R.H.: Dynamika układów fizycznych, WNT, Warszawa 1973.
[2] Siouris G.M.: Missile Guidance and Control Systems, Springer-Verlag, New York 2004.
[3] Żelazny M.: Podstawy automatyki, PWN, Warszawa 1976.
92
Wyszukiwarka
Podobne podstrony:
092 Sanders Glenda Doktor Hunk092 93Lesson Plan 083 Text087 092 zzhjfir6ptmapevkxbqtgyde3cqjmuka5gvkruaF G 083092 ADMM083 085092 093312[01] 0X 092 CZERWIEC 09 CZ 2v 04 083092 095F G 092091 092092 Kardynalny błądwięcej podobnych podstron