SYMULACJA ZAKŁÓCEŃ W UKŁADACH AUTOMATYKI UTWORZONYCH ZA POMOCĄ OBWODÓW ELEKTRYCZNYCH W PROGRAMACH MATHCAD I PSPICE

(1)

__________________________________________

* Uniwersytet Szczeciński.

Piotr FRĄCZAK*

SYMULACJA ZAKŁÓCEŃ W UKŁADACH AUTOMATYKI UTWORZONYCH ZA POMOCĄ OBWODÓW ELEKTRYCZNYCH W PROGRAMACH MATHCAD

I PSPICE

Praca zawiera opis analityczny metodą zmiennych stanu (zapis macierzowy) obwodu elektrycznego. Jako zmienne stanu przyjęto prądy w cewkach i napięcie na kondensatorze.

Zmienne stanu opisują określony układ automatyki. Wybrane układy automatyki wielowymiarowe poddane działaniu sygnałów zakłócających ujęto analitycznie układami równań różniczkowych pierwszego rzędu (zapis macierzowy). Sygnały zakłócające utworzono z funkcji trygonometrycznych i procedury rnd(x) programu Mathcad generującej liczby pseudolosowe o rozkładzie naturalnym. Obliczenia symulacyjne sygnałów wyjściowe wybranych układów automatyki z zakłóceniami wyznaczono w oparciu o ich opisy analityczne (układ równań różniczkowych pierwszego rzędu). Układ równań różniczkowych pierwszego rzędu rozwiązano numerycznie metodą równań różnicowych w programie Mathcad. Uzyskane wyniki obliczeń symulacyjnych sygnałów wyjściowych wybranego układu automatyki zweryfikowano w programie PSpice. Weryfikacja w programie PSpice polegała na stworzeniu schematu zastępczego danego obwodu elektrycznego stanowiącego układ automatyki i zadeklarowania odpowiedniej analizy.

1. WSTĘP

Obwody elektryczne opisuje się między innymi układami n równań różniczkowych pierwszego rzędu (metoda zmiennych stanu) z czasem ciągłym.

Układ n równań różniczkowych w zapisie wektorowo - macierzowym ma postać:

1 11 12 1n 1 11 12 1m 1

2 21 22 2n 2 21 22 2m 2

n n1 n2 nn n n1 n2 nm m

x (t) a a a x (t) b b b u (t)

d dt

x (t) a a a x (t) b b b u (t)

         

         

       

         

         

 

          

 

(1)

Stąd w ujęciu uproszczonym

d (t) (t) (t)

dtx A x B u (2)

(2)

układ automatyki), x jest wektorem stanu (wektor stanu – układ automatyki).

Równanie (2) ujmujące opis obwodu elektrycznego w postaci macierzowej odpowiada określonej strukturze schematu blokowego (rys.1). Schemat blokowy zamieszczony na rys.1 stanowi układ automatyki. W układach automatyki x(t) oznacz pochodną wektora stanu.

Rys. 1. Struktura układu opisanego równaniem (2)

2. OPIS ANALITYCZNY WYBRANYCH UKŁADÓW AUTOMATYKI Z SYGNAŁAMI ZAKŁÓCAJĄCYMI

Klasyczne badanie układów automatyki w relacji: wejście sterujące – wyjście lub wejście zakłócające – wyjście nie ujmują bezpośrednio w całości obrazu dynamiki danego układu. Aby uzyskać całościowe ujęcie dynamiki danego układu, opisuje się go metodą zmiennych stanu. Metoda ta pozwala opisać stan układu automatyki w każdej chwili

[4]

.Poniżej zamieszczono wybrane układy automatyki z uwzględnieniem sygnałów zakłócających i odpowiednio je opisano.

2.1. Działanie sygnału zakłócającego na wejście układy automatyki

Dynamikę układu liniowego ciągłego stacjonarnego wielowymiarowego (rys.

2), w którym występuje oddziaływanie sygnałów zakłócenia na wejście układu ujmują równanie stanu (3) i równanie wyjścia układu (4).

(t) =  (t) +  (t) +  (t)

x A x B u E z (3)

(t) =  (t)

y C x (4)

przy czym: (t)x  ⁿ- wektor stanu; u(t) ^m- wektor wymuszeń; y(t)  ^l- wektor odpowiedzi; z(t)  ^m-wektor sygnałów zakłócających; A - macierz stanu o wymiarze nn; B - macierz wejść o wymiarze nm;C - macierz wyjść o wymiarze ln ; E - macierz wejść sygnałów zakłócających o wymiarze nm.

(3)

Rys. 2. Schemat zastępczy układu automatyki z sygnałami zakłócającymi na wejściu

2.2. Działanie sygnału zakłócającego na wyjście układy automatyki

Układ automatyki liniowy ciągły stacjonarny wielowymiarowy (rys. 3), w którym występuje oddziaływanie sygnałów zakłócenia na wejście układu ujmują równanie stanu (5) i równanie wyjścia układu (6).

(t) =  (t) +  (t)

x A x B u (5)

(t) =  (t) +  (t)

y C x H z (6)

przy czym:x(t),u(t), y(t),z(t)- wektory określone wzorami (3) i (4); A,B,C- macierze określone równaniami (3) i (4); H - macierz wyjść sygnałów zakłócających o wymiarze lm.

Rys. 3. Schemat zastępczy układu automatyki z sygnałami zakłócającymi na wyjściu

2.3. Działanie sygnałów zakłócających na wejście i wyjście układu automatyki Dynamikę układu liniowego ciągłego stacjonarnego wielowymiarowego (rys. 4), w którym występuje oddziaływanie sygnałów zakłócenia na wejście i wyjście układu ujmują równanie stanu (7) i równanie wyjścia układu (8).

(t) =  (t) +  (t) +  (t)

x A x B u E z (7)

(t) =  (t) +  (t)

y C x H z (8)

(4)

Rys. 4. Schemat zastępczy (blokowy) układu automatyki z sygnałami zakłócającymi na wejściu i wyjściu

3. OBLICZENIA SYMULACYJNE UKŁADÓW AUTOMATYKI

3.1. Symulacja sygnałów wyjściowych układu automatyki utworzonego za pomocą obwodu elektrycznego w programie Mathcad

Do obliczeń symulacyjnych sygnałów wyjściowych w układach automatyki utworzonych za pomocą obwodów elektrycznych wykorzystano obwód elektryczny, którego schemat zastępczy przedstawiono w programie PSpice i zamieszczono na rys. 5.

Zamieszczony na rys. 5 obwód elektryczny opisano za pomocą praw Kirchhoffa stosownie do metody zmiennych stanu. Prawa Kirchhoffa dla oczek i węzłów mają postać:

1

1 3 2 1 1 3 1 2

2 1

3 2 2 1 2 1 , 1

       

           

c

c c

di (t)

E R i (t) R i (t) L u (t) 0 , i (t)= i (t)+ i (t) dt

du (t)

di (t) di (t) 1

R i (t) L u (t) L R i (t) 0 i (t)

dt dt dt C

(9)

Jako zmienne stanu przyjęto prądy (i (t) ,₁ i (t) ) w cewkach indukcyjnych i ₂ napięcieu (t) na kondensatorze. W celu ujednolicenia oznaczeń podczas obliczeń _c numerycznych obwodów elektrycznych i układów automatyki zmienne stanu oznaczono następująco: x1(t) i (t) ,  ₁ x2(t) i (t) oraz ₂ x3(t) u (t) . Wprowadzając  _c wymienione oznaczenia prądów i napięcia w rozpatrywanym obwodzie elektrycznym otrzymano układ trzech równań różniczkowych pierwszego rzędu w postaci macierzowej (10).

(5)

1 1 1

1

2 2 2

, ,

1 2 1

1 3

1

(R + R ) R 1

dx1(t)

1 0

L L L

dt x1(t)

L

E R (R + R ) 1

dx2(t)

x2(t) 0

0 L E

dt L L

x3(t)

dx3(t) 1 0 0

0 0

dt C

 

     

       

   

 

   

           

         

       

   

 

   

A B A B

(10)

Rys. 5. Schemat zastępczy obwodu elektrycznego w programie PSpice

Przyjmującdx(t) x^{n 1} xⁿ

dt t

 

 oraz dokonując dyskretyzacji czasu z krokiem t ( t T_end/N; n0 1, ,..., ;N t_n   n t ) otrzymano równania różnicowe (11).

,

n 1 n n 1 n n

dx1(t)

dt x1 x1 x1 x1 x1

dx2(t) 1 E

x2 x2 x2 x2 x2 t t

dt t E

x3 x3 x3 x3 x3

dx3(t) dt

 



 

 

 

        

 

         

             

            

         

 

 

A B (11)

Obliczenia symulacyjne oparte o równania różnicowe (11) przeprowadzono w programie Mathcad.

(6)

Rys. 6. Symulacja sygnałów wyjściowych układu automatyki utworzonego za pomocą obwodu elektrycznego

3.2. Weryfikacja obliczeń zmiennych stanu obwodu elektrycznego

stanowiącego sygnały wyjściowe układu automatyki w programie PSpice W celu sprawdzenia wyników obliczeń numerycznych zmiennych stanu (x1, x2, x3) obwodu elektrycznego (rys. 5) w programie Mathcad przeprowadzono dodatkowe obliczenia symulacyjne w programie PSpice. Obliczenia symulacyjne

(7)

w programie PSpice polegały na dołączeniu w chwili t = 0 (Sw tClose - wyłącznnik) napięcia stałego E = 10 V do rozpatrywanego obwodu elektrycznego oraz zadeklarowania odpowiedniej analizy. Wskutek dołączenia napięcia do obwodu elektrycznego w czasie t = 0 wystąpił w nim stan nieustalony. Stany nieustalone prądów i napięcia w obwodzie elektrycznym odpowiadają zmiennym stanu (x1, x2, x3). Stan nieustalony obwodu elektrycznego w programie PSpice ujmuje procedura Transient Analysis oraz parametry wyłącznika (Sw tClose).Transient Analysis - Print Step = 7 us; Final Time = 0.7 us; No – Print Deley = 0; Step Celling = 7 us. Sw tClose - tClose = 0; ttran = 0.01ns; Rclosed = 0.0001; Ropen = 1Meg. Obliczenia symulacyjne w programie PSpice prądów i napięcia będących odpowiednio zmiennymi stanu zamieszczono na rys. 7. W obu programach numerycznych Mathcad (rys.6) i PSpice (rys.7) uzyskano takie same wartości zmiennych stanu.

Rys. 7. Symulacja sygnałów wyjściowych układu automatyki utworzonego za pomocą obwodu elektrycznego w programie PSpice

3.3. Symulacja obliczeń sygnałów wyjściowych układu automatyki z oddziaływaniem sygnałów zakłócających

W celu dokonania obliczeń symulacyjnych sygnałów wyjściowych układu automatyki, który został utworzony z obwodu elektrycznego i podany działaniu sygnałów zakłócających na wejściu i wyjściu wykorzystano odpowiednio równania (7) i (8). Utworzono sygnały zakłócające z funkcji trygonometrycznych i liczb pseudolosowych o rozkładzie naturalnym. Równania (7) i (8) w zapisie numerycznym zaimplementowano w środowisko programu Mathcad.

(8)

Rys. 8. Symulacja sygnałów wyjściowych układu automatyki utworzonego za pomocą obwodu elektrycznego z oddziaływaniem sygnałów zakłócających

na jego wejściu i wyjściu w programie Mathcad

(9)

PODSUMOWANIE

 Za pomocą obwodów elektrycznych można modelować układy automatyki.

 Układy automatyki opisane równaniami różniczkowym pierwszego rzędu można analizować metodą równań różnicowych w ujęciu procedur obliczeniowych programu Mathcad.

 Za pomocą procedur obliczeniowych programu Mathcad można dokonywać wizualizacji dynamiki układów automatyki z uwzględnieniem sygnałów zakłócających.

LITERATURA

[1] Frączak P.: Obliczenia numeryczne obwodów elektrycznych i układów cyfrowych, Szczecin, Wydawnictwo PPH ,,Zapol’’ Dmochowski, Sobczyk Sp.j. 2012, 173 s., ISBN 978-83-7518-432-7.

[2] Król A., Moczko J.: PSpice Symulacja i optymalizacja układów elektronicznych, Wydawnictwo Nakom, Poznań 1999, ISBN 83 - 86969-28-8.

[3] Palczewski W.: Mathcad 12,11, 2001i, 2000 w algorytmach, Akademicka Oficyna Wydawnicza EXIT, Warszawa 2005, ISBN 83-87674-81-8.

[4] Pełczewski W.: Teoria sterowania, Wydawnictwo Naukowo-Techniczne, Warszawa 1980, ISBN 83 - 204-0187-9.

SIMULATION OF INTERFERENCE AUTOMATION SYSTEM CREATED BY ELECTRICAL CIRCUIT PROGRAM MATHCAD I PSPICE

The paper contains a description of the analytical method for the state variables (matrix notation) of circuit. As the state variables currents in the coils and the voltage across the capacitor has been adopted. State variables describe the specific automation system. Some automation systems which were subjected to multi-dimensional interference signals have been recognized analytically by systems of first-order differential equations (matrix notation). Spurious signals have been created with trigonometric functions and procedures rnd (x) in Mathcad program. These functions were generating pseudo-random numbers on a natural distribution. Output signals simulations of selected automation systems with the interference has been based on its analytical descriptions (first-order differential equations system). System of first-order differential equations has been solved numerically using the differential equations method in the Mathcad program. The results of simulation calculations of the output signals from the selected automation system have been verified in PSpice program. Verification in the PSpice program was to create equivalent circuit diagram constituting automation system and declare the appropriate analysis.