4. Metoda operatorowa analizy stanów nieustalonych w obwodach elektrycznych

4.6. Obliczenia prądów i napięć w stanie nieustalonym metodą operatorową

Obliczenia prądów i napięć w stanie nieustalonym obwodu metodą operatorową sprowadzać się będą do wyznaczenia transformaty odpowiedniej wielkości a następnie obliczenia transformaty odwrotnej Laplace’a dla określenia zmiennej w dziedzinie czasu. Do obliczenia transformat prądów i napięć można stosować wszystkie poznane dotąd metody analizy obwodów, w tym metodę równań Kirchhoffa, oczkową, potencjałów węzłowych, Thevenina i Nortona operujące transformatami Laplace’a zamiast wartościami zespolonymi czy wartościami w dziedzinie czasu (dla obwodu rezystancyjnego).
    Podstawowymi zaletami metody operatorowej jest łatwość uwzględnienia niezerowych warunków początkowych (przez wprowadzenie źródeł napięciowych w modelu operatorowym) oraz sprowadzenie operacji różniczkowych do działań algebraicznych. 
W ogólności rozwiązując stan nieustalony w obwodzie metodą operatorową należy wyróżnić kilka etapów.

  1. Określenie warunków początkowych w obwodzie, poprzez wyznaczenie rozwiązania ustalonego obwodu przed przełączeniem i obliczenie wartości napięć na kondensatorach i prądów cewek w chwili \(t=0^-\), to jest \(i_L (0^-)\) oraz \(u_C (0^-)\)
  2. Określenie rozwiązania obwodu w stanie ustalonym po przełączeniu przy zastosowaniu metody symbolicznej z wykorzystaniem dowolnej metody analizy. Wynikiem jest postać czasowa rozwiązania ustalonego prądów cewek \(i_Lu (t)\) i napięć kondensatorów \(u_Cu (t)\). Przez założenie t=0 otrzymuje się wartości prądów i napięć w chwili początkowej, to jest \(i_Lu (0^+)\) oraz \(u_Cu (0^+)\).
  3. Określenie rozwiązania obwodu w stanie przejściowym po przełączeniu przy zastosowaniu metody operatorowej. Dla otrzymania takiego rozwiązania należy wykonać następujące etapy:
  • utworzenie schematu obwodu dla składowej przejściowej poprzez wyeliminowanie źródeł zewnętrznych wymuszających (zwarcie źródeł napięcia i rozwarcie źródeł prądu); obwód rzeczywisty dla składowej przejściowej w dziedzinie czasu nie zawiera żadnych źródeł wymuszających
  • określenie warunków początkowych dla składowej przejściowej przy wykorzystaniu praw komutacji, zgodnie z którymi \(x(0^-)=x_u (0^+)+x_p (0^+)\); z równania tego wynikają następujące wzory na warunki początkowe dla składowych przejściowych prądu cewki i napięcia kondensatora
\(i_Lp (0^+)=i_L (0^-)-i_Lu (0^+)\) (4.10)
\(u_Cp (0^+)=u_C (0^-)-u_Cu (0^+)\) (4.11)
  • utworzenie schematu operatorowego obwodu w stanie przejściowym poprzez zastąpienie elementów rzeczywistych obwodu ich modelami operatorowymi dla składowej przejściowej i rozwiązanie obwodu względem poszukiwanych prądów i napięć operatorowych
  • wyznaczenie transformaty odwrotnej Laplace’a dla poszukiwanych wielkości przejściowych określonych w punkcie poprzednim; w wyniku otrzymuje się i\(i_Lp (t) \) oraz \(u_Cp (t).\)
  1. Określenie rozwiązania obwodu w stanie przejściowym po przełączeniu przy zastosowaniu metody operatorowej. Dla otrzymania takiego rozwiązania należy wykonać następujące etapy:
\(i_Lp (0^+)=i_L (0^-)-i_Lu (0^+)\) (4.12)
\(u_Cp (0^+)=u_C (0^-)-u_Cu (0^+)\) (4.13)

Składowa przejściowa zanika z czasem do zera i pozostaje jedynie składowa ustalona określająca przebieg wielkości w stanie ustalonym. Taka metodyka rozwiązania stanów nieustalonych przy zastosowaniu transformacji Laplace’a nosi nazwę metody superpozycji stanów, gdyż rozdziela w sposób jawny stan ustalony od stanu przejściowego. Jest szczególnie zalecana przy wymuszeniach sinusoidalnych, choć obowiązuje również dla obwodów prądu stałego. Zaletą takiego podejścia jest jej uniwersalność i stosowalność do każdego obwodu liniowego RLC niezależnie od rodzaju wymuszenia (wymuszenia stałe lub sinusoidalne mają jedynie wpływ na stan ustalony i są wyeliminowane przy rozwiązywaniu stanu przejściowego). 
 Należy podkreślić, że rozbicie stanu nieustalonego na ustalony i przejściowy jest zalecane jedynie przy istnieniu wymuszeń sinusoidalnych w obwodzie po przełączeniu. Jeśli źródła takie nie występują schemat operatorowy może dotyczyć obwodu całkowitego, bez rozbijania go na schemat dla składowej ustalonej i przejściowej. W takim przypadku pozostawia się zewnętrzne źródła wymuszające w obwodzie przyjmując ich model operatorowy, czyli zastępując postać czasową źródła (wartość stała A przy wymuszeniu stałym) przez funkcję \(\frac{A}{s}\). Warunki początkowe również nie podlegają modyfikacji, co oznacza, że \(i_L (0^+)=i_L (0^-)\) oraz \(u_C (0^+)=u_C (0^-)\).

Wyznaczyć przebieg czasowy napięcia na kondensatorze w stanie nieustalonym w obwodzie z rys. 4.4 po przełączeniu. Dane liczbowe parametrów obwodu są następujące: \(R_1=R_2=1Ω\), \(L=1H\), \(C=1F\), \(e_2 (t)=10V\). Źródło wymuszające sinusoidalne dane jest w następującej postaci \(e_1 (t)=5√2  sin⁡( t+π/4) V\).
 

Uzupelnij opis obrazka

Rys. 4.4. Schemat obwodu do przykładu 4.10.


Rozwiązanie
W rozwiązaniu problemu wyznaczymy najpierw warunki początkowe w obwodzie rozwiązując stan ustalony przed przełączeniem. Ponieważ przed przełączeniem w obwodzie występowały dwa źródła: stałe i sinusoidalne w obliczeniu warunków początkowych (stan ustalony przed przełączeniem) należy zastosować metodę superpozycji źródeł. 

Uzupelnij opis obrazka Uzupelnij opis obrazka
Uzupelnij opis obrazka Uzupelnij opis obrazka

Rys. 4.5 Schematy obwodu: a) w stanie ustalonym przed przełączeniem (źródło sinusoidalne), b) w stanie ustalonym przed przełączeniem (źródło stałe), c) w stanie ustalonym po przełączeniu,
d) schemat operatorowy dla składowej przejściowej

Schemat obwodu w stanie ustalonym przed przełączeniem przy wymuszeniu sinusoidalnym przedstawiony jest na rys. 4.5a. Wobec rezonansu równoległego w gałęzi LC prąd wydawany przez źródło jest równy zeru a napięcie na tej gałęzi jest równe napięciu źródła. Stąd 

\(u_{Cu}^{(1)}(t)=5\sqrt2sin{(}t+\pi/4)\)

Prąd cewki (wartość skuteczna zespolona) dany jest wzorem

\(I_{Lu}^{(1)}=\frac{5e^{j\pi/4}}{j1}=5e^{-j\pi/4}\)

co odpowiada postaci czasowej 

\(i_{Lu}^{(1)}(t)=5\sqrt2sin{(}t-\pi/4)\)

Uwzględniając źródło stałe e2(t) uzyskuje się znaczne uproszczenie obwodu (cewka dla prądu stałego w stanie ustalonym stanowi zwarcie a kondensator przerwę) jak to przedstawiono na rys. 4.5b. Rozwiązanie na prąd cewki i napięcie kondensatora ma więc postać: 

\(u_{Cu}^{(2)}(t)=0\)

\(i_{Lu}^{(2)}(t)=\frac{10}{1}=10\)

Dokonując superpozycji obu rozwiązań otrzymuje się

\(i_{Lu}(t)=i_{Lu}^{(1)}(t)+i_{Lu}^{(2)}(t)=10+5\sqrt2sin{(}t-\pi/4)\)

\(u_{Cu}(t)=u_{Cu}^{(1)}(t)+u_{Cu}^{(2)}(t)=5\sqrt2sin{(}t+\pi/4)\)

Stąd warunki początkowe są następujące: \( u_C(0^-)=5\), \(i_L(0^-)=5\).
Po przełączeniu w obwodzie pozostaje jedynie źródło sinusoidalne e2(t). Schemat obwodu dla tego wymuszenia pokazany jest na rys. 4.5c. Z analizy tego obwodu wynika następująca procedura rozwiązania. Wobec rezonansu równoległego w gałęzi LC prąd wydawany przez źródło jest równy zeru a napięcie na tej gałęzi jest równe napięciu źródła. Stąd 

\(u_{Cu}(t)=5\sqrt2sin{(}t+\pi/4)\)

Prąd cewki (wartość skuteczna zespolona) dany jest wzorem

\(I_{Lu}=\frac{5e^{j45^\circ}}{j1}=5e^{-j\pi/4}\)

co odpowiada postaci czasowej 

\(i_{Lu}(t)=5\sqrt2sin{(}t-\pi/4)\)

Stan początkowy dla składowej ustalonej prądu cewki i napięcia kondensatora przyjmuje więc następujące wartości:

\(u_{Cu}(0^+)=5\)

oraz

\(i_{Lu}(0^+)=-5\)

Warunki początkowe dla składowej przejściowej prądu i napięcia są zatem równe:

\(u_{Cp}(0^+)=u_C(0^-)-u_{Cu}(0^+)=0\)

\(i_{Lp}(0^+)=i_L(0^-)-i_{Lu}(0^+)=10\)

Schemat operatorowy obwodu przedstawiono na rys. 4.5d (źródło wewnętrzne przy kondensatorze nie występuje, bo \(u_{Cp}(0^+)=0\). Zastosowanie metody potencjałów węzłowych do wyznaczenia postaci operatorowej rozwiązania prowadzi do wyniku 

\(U_{Cp}(s)=\frac{-\frac{10}{s}}{s+\frac{1}{s}+0,5}=\frac{-10}{s^2+0,5s+1}\)

Wobec zespolonych wartości własnych (pierwiastków mianownika transformaty napięcia) w wyznaczaniu oryginału zastosujemy metodę wykorzystującą tablice transformat. W związku z powyższym transformatę przedstawimy w postaci przekształconej

\(U_{Cp}(s)=\frac{-10}{s^2+0,5s+1}=-\frac{\left(\sqrt{\frac{15}{16}}\right)\cdot10\cdot\sqrt{\frac{16}{15}}}{\left(s+0,25\right)^2+\left(\sqrt{\frac{15}{16}}\right)^2}=-10,33\frac{\left(\sqrt{\frac{15}{16}}\right)}{\left(s+0,25\right)^2+\left(\sqrt{\frac{15}{16}}\right)^2}\)

Powyższej funkcji operatorowej można przyporządkować następującą postać czasową (patrz wiersz szósty tablicy 12.1)

\(u_{Cp}(t)=-10,33e^{-0,25t}sin{\left(\sqrt{\frac{15}{16}}t\right)}\)

Rozwiązanie całkowite określające napięcie kondensatora jest sumą składowej ustalonej i przejściowej

\(u_C(t)=u_{Cu}(t)+u_{Cp}(t)=5\sqrt2sin{(}t+\pi/4)-10,33e^{-0,25t}sin{\left(\sqrt{\frac{15}{16}}t\right)}\)

Składowa przejściowa zanika z biegiem czasu ze stałą czasową \(\tau=4\) i po około 5 stałych czasowych pozostaje jedynie składowa ustalona sinusoidalna.