Finding General Solution for Rlc Circuit
6.3: The RLC Circuit
-
- Last updated
- Save as PDF
- Page ID
- 9427
-
- Andrew G. Cowles Distinguished Professor Emeritus (Mathematics) at Trinity University
In this section we consider the \(RLC\) circuit, shown schematically in Figure 6.3.1 . As we'll see, the \(RLC\) circuit is an electrical analog of a spring-mass system with damping.
Nothing happens while the switch is open (dashed line). When the switch is closed (solid line) we say that the circuit is closed . Differences in electrical potential in a closed circuit cause current to flow in the circuit. The battery or generator in Figure 6.3.1 creates a difference in electrical potential \(E=E(t)\) between its two terminals, which we've marked arbitrarily as positive and negative. (We could just as well interchange the markings.) We'll say that \(E(t)>0\) if the potential at the positive terminal is greater than the potential at the negative terminal, \(E(t)<0\) if the potential at the positive terminal is less than the potential at the negative terminal, and \(E(t)=0\) if the potential is the same at the two terminals. We call \(E\) the impressed voltage .
At any time \(t\), the same current flows in all points of the circuit. We denote current by \(I=I(t)\). We say that \(I(t)>0\) if the direction of flow is around the circuit from the positive terminal of the battery or generator back to the negative terminal, as indicated by the arrows in Figure 6.3.1 \(I(t)<0\) if the flow is in the opposite direction, and \(I(t)=0\) if no current flows at time \(t\).
Differences in potential occur at the resistor, induction coil, and capacitor in Figure 6.3.1 . Note that the two sides of each of these components are also identified as positive and negative. The voltage drop across each component is defined to be the potential on the positive side of the component minus the potential on the negative side. This terminology is somewhat misleading, since "drop" suggests a decrease even though changes in potential are signed quantities and therefore may be increases. Nevertheless, we'll go along with tradition and call them voltage drops. The voltage drop across the resistor in Figure 6.3.1 is given by
\[\label{eq:6.3.1} V_R=IR,\]
where \(I\) is current and \(R\) is a positive constant, the resistance of the resistor. The voltage drop across the induction coil is given by
\[\label{eq:6.3.2} V_I=L{dI\over dt}=LI',\]
where \(L\) is a positive constant, the inductance of the coil.
A capacitor stores electrical charge \(Q=Q(t)\), which is related to the current in the circuit by the equation
\[\label{eq:6.3.3} Q(t)=Q_0+\int_0^tI(\tau)\,d\tau,\]
where \(Q_0\) is the charge on the capacitor at \(t=0\). The voltage drop across a capacitor is given by
\[\label{eq:6.3.4} V_C={Q\over C},\]
where \(C\) is a positive constant, the capacitance of the capacitor.
Table 6.3.1 names the units for the quantities that we've discussed. The units are defined so that
\[\begin{aligned} 1\mbox{volt}&= 1 \text{ampere} \cdot1 \text{ohm}\\ &=1 \text{henry}\cdot1\,\text{ampere}/\text{second}\\ &= 1\text{coulomb}/\text{farad}\end{aligned} \nonumber \]
and
\[\begin{aligned} 1 \text{ampere}&=1\text{coulomb}/\text{second}.\end{aligned} \nonumber \]
| Symbol | Name | Unit |
|---|---|---|
| \(E\) | Impressed Voltage | volt |
| \(I\) | Current | ampere |
| \(Q\) | Charge | coulomb |
| \(R\) | Resistance | ohm |
| \(L\) | Inductance | henry |
| \(C\) | Capacitance | farad |
According to Kirchoff's law, the sum of the voltage drops in a closed \(RLC\) circuit equals the impressed voltage. Therefore, from Equation \ref{eq:6.3.1}, Equation \ref{eq:6.3.2}, and Equation \ref{eq:6.3.4},
\[\label{eq:6.3.5} LI'+RI+{1\over C}Q=E(t).\]
This equation contains two unknowns, the current \(I\) in the circuit and the charge \(Q\) on the capacitor. However, Equation \ref{eq:6.3.3} implies that \(Q'=I\), so Equation \ref{eq:6.3.5} can be converted into the second order equation
\[\label{eq:6.3.6} LQ''+RQ'+{1\over C}Q=E(t)\]
in \(Q\). To find the current flowing in an \(RLC\) circuit, we solve Equation \ref{eq:6.3.6} for \(Q\) and then differentiate the solution to obtain \(I\).
In Sections 6.1 and 6.2 we encountered the equation
\[\label{eq:6.3.7} my''+cy'+ky=F(t)\]
in connection with spring-mass systems. Except for notation this equation is the same as Equation \ref{eq:6.3.6}. The correspondence between electrical and mechanical quantities connected with Equation \ref{eq:6.3.6} and Equation \ref{eq:6.3.7} is shown in Table 6.3.2 .
| Electrical | Mechanical |
|---|---|
| charge \(Q\) | displacement \(y\) |
| current \(I\) | velocity \(y'\) |
| impressed voltage \(E(t)\) | external force \(F(t)\) |
| inductance \(L\) | mass \(m\) |
| resistance \(R\) | damping \(c\) |
| 1/capacitance \(1/C\) | spring constant \(k\) |
The equivalence between Equation \ref{eq:6.3.6} and Equation \ref{eq:6.3.7} is an example of how mathematics unifies fundamental similarities in diverse physical phenomena. Since we've already studied the properties of solutions of Equation \ref{eq:6.3.7} in Sections 6.1 and 6.2, we can obtain results concerning solutions of Equation \ref{eq:6.3.6} by simply changing notation, according to Table 6.3.1 .
Free Oscillations
We say that an \(RLC\) circuit is in free oscillation if \(E(t)=0\) for \(t>0\), so that Equation \ref{eq:6.3.6} becomes
\[\label{eq:6.3.8} LQ''+RQ'+{1\over C}Q=0.\]
The characteristic equation of Equation \ref{eq:6.3.8} is
\[Lr^2+Rr+{1\over C}=0,\nonumber\]
with roots
\[\label{eq:6.3.9} r_1={-R-\sqrt{R^2-4L/C}\over2L}\quad \text{and} \quad r_2= {-R+\sqrt{R^2-4L/C}\over2L}.\]
There are three cases to consider, all analogous to the cases considered in Section 6.2 for free vibrations of a damped spring-mass system.
Case 1: Underdamped Motion
The oscillation is underdamped if \(R<\sqrt{4L/C}\). In this case, \(r_1\) and \(r_2\) in Equation \ref{eq:6.3.9} are complex conjugates, which we write as
\[r_1=-{R\over2L}+i\omega_1\quad \text{and} \quad r_2=-{R\over2L}-i\omega_1,\nonumber\]
where
\[\omega_1={\sqrt{4L/C-R^2}\over2L}.\nonumber\]
The general solution of Equation \ref{eq:6.3.8} is
\[Q=e^{-Rt/2L}(c_1\cos\omega_1 t+c_2\sin\omega_1 t),\nonumber\]
which we can write as
\[\label{eq:6.3.10} Q=Ae^{-Rt/2L}\cos(\omega_1 t-\phi),\]
where
\[A=\sqrt{c_1^2+c_2^2},\quad A\cos\phi=c_1,\quad \text{and} \quad A\sin\phi=c_2.\nonumber\]
In the idealized case where \(R=0\), the solution Equation \ref{eq:6.3.10} reduces to
\[Q=A\cos\left({t\over\sqrt{LC}}-\phi\right),\nonumber\]
which is analogous to the simple harmonic motion of an undamped spring-mass system in free vibration.
Actual \(RLC\) circuits are usually underdamped, so the case we've just considered is the most important. However, for completeness we'll consider the other two possibilities.
Case 2: Overdamped Motion
The oscillation is overdamped if \(R>\sqrt{4L/C}\). In this case, the zeros \(r_1\) and \(r_2\) of the characteristic polynomial are real, with \(r_1 < r_2 <0\) (see \ref{eq:6.3.9}), and the general solution of \ref{eq:6.3.8} is
\[\label{eq:6.3.11} Q=c_1e^{r_1t}+c_2e^{r_2t}.\]
Case 3: Critically Damped
The oscillation is critically damped if \(R=\sqrt{4L/C}\). In this case, \(r_1=r_2=-R/2L\) and the general solution of Equation \ref{eq:6.3.8} is
\[\label{eq:6.3.12} Q=e^{-Rt/2L}(c_1+c_2t).\]
If \(R\ne0\), the exponentials in Equation \ref{eq:6.3.10}, Equation \ref{eq:6.3.11}, and Equation \ref{eq:6.3.12} are negative, so the solution of any homogeneous initial value problem
\[LQ''+RQ'+{1\over C}Q=0,\quad Q(0)=Q_0,\quad Q'(0)=I_0,\nonumber\]
approaches zero exponentially as \(t\to\infty\). Thus, all such solutions are transient , in the sense defined Section 6.2 in the discussion of forced vibrations of a spring-mass system with damping.
Example 6.3.1
At \(t=0\) a current of 2 amperes flows in an \(RLC\) circuit with resistance \(R=40\) ohms, inductance \(L=.2\) henrys, and capacitance \(C=10^{-5}\) farads. Find the current flowing in the circuit at \(t>0\) if the initial charge on the capacitor is 1 coulomb. Assume that \(E(t)=0\) for \(t>0\).
Solution
The equation for the charge \(Q\) is
\[{1\over5}Q''+40Q'+10000Q=0, \nonumber \]
or
\[\label{eq:6.3.13} Q''+200Q'+50000Q=0.\]
Therefore we must solve the initial value problem
\[\label{eq:6.3.14} Q''+200Q'+50000Q=0,\quad Q(0)=1,\quad Q'(0)=2.\]
The desired current is the derivative of the solution of this initial value problem.
The characteristic equation of Equation \ref{eq:6.3.13} is
\[r^2+200r+50000=0, \nonumber \]
which has complex zeros \(r=-100\pm200i\). Therefore the general solution of Equation \ref{eq:6.3.13} is
\[\label{eq:6.3.15} Q=e^{-100t}(c_1\cos200t+c_2\sin200t).\]
Differentiating this and collecting like terms yields
\[\label{eq:6.3.16} Q'=-e^{-100t}\left[(100c_1-200c_2)\cos200t+ (100c_2+200c_1)\sin200t\right].\]
To find the solution of the initial value problem Equation \ref{eq:6.3.14}, we set \(t=0\) in Equation \ref{eq:6.3.15} and Equation \ref{eq:6.3.16} to obtain
\[c_1=Q(0)=1\quad \text{and} \quad -100c_1+200c_2=Q'(0)=2;\nonumber\]
therefore, \(c_1=1\) and \(c_2=51/100\), so
\[Q=e^{-100t}\left(\cos200t+{51\over100}\sin200t\right)\nonumber\]
is the solution of Equation \ref{eq:6.3.14}. Differentiating this yields
\[I=e^{-100t}(2\cos200t-251\sin200t).\nonumber\]
Forced Oscillations With Damping
An initial value problem for Equation \ref{eq:6.3.6} has the form
\[\label{eq:6.3.17} LQ''+RQ'+{1\over C}Q=E(t),\quad Q(0)=Q_0,\quad Q'(0)=I_0,\]
where \(Q_0\) is the initial charge on the capacitor and \(I_0\) is the initial current in the circuit. We've already seen that if \(E\equiv0\) then all solutions of Equation \ref{eq:6.3.17} are transient. If \(E\not\equiv0\), we know that the solution of Equation \ref{eq:6.3.17} has the form \(Q=Q_c+Q_p\), where \(Q_c\) satisfies the complementary equation, and approaches zero exponentially as \(t\to\infty\) for any initial conditions, while \(Q_p\) depends only on \(E\) and is independent of the initial conditions. As in the case of forced oscillations of a spring-mass system with damping, we call \(Q_p\) the steady state charge on the capacitor of the \(RLC\) circuit. Since \(I=Q'=Q_c'+Q_p'\) and \(Q_c'\) also tends to zero exponentially as \(t\to\infty\), we say that \(I_c=Q'_c\) is the transient current and \(I_p=Q_p'\) is the steady state current. In most applications we are interested only in the steady state charge and current.
Example 6.3.2
Find the amplitude-phase form of the steady state current in the \(RLC\) circuit in Figure 6.3.1 if the impressed voltage, provided by an alternating current generator, is \(E(t)=E_0\cos\omega t\).
Solution
We'll first find the steady state charge on the capacitor as a particular solution of
\[LQ''+RQ'+{1\over C}Q=E_0\cos\omega t.\nonumber\]
To do, this we'll simply reinterpret a result obtained in Section 6.2, where we found that the steady state solution of
\[my''+cy'+ky=F_0\cos\omega t \nonumber\]
is
\[y_p={F_0\over\sqrt{(k-m\omega^2)^2+c^2\omega^2}} \cos(\omega t-\phi), \nonumber\]
where
\[\cos\phi={k-m\omega^2\over\sqrt {(k-m\omega^2)^2+c^2\omega^2}}\quad \text{and} \quad \sin\phi={c\omega\over\sqrt{(k-m\omega^2)^2+c^2\omega^2}}. \nonumber\]
(see Equations \ref{eq:6.3.14} and Equation \ref{eq:6.3.15}.) By making the appropriate changes in the symbols (according to Table 6.3.2 ) yields the steady state charge
\[Q_p={E_0\over\sqrt{(1/C-L\omega^2)^2+R^2\omega^2}}\cos(\omega t-\phi), \nonumber\]
where
\[\cos\phi={1/C-L\omega^2\over\sqrt{(1/C-L\omega^2)^2+R^2\omega^2}} \quad \text{and} \quad \sin\phi={R\omega\over\sqrt{(1/C-L\omega^2)^2+R^2\omega^2}}. \nonumber\]
Therefore the steady state current in the circuit is
\[I_p=Q_p'= -{\omega E_0\over\sqrt{(1/C-L\omega^2)^2+R^2\omega^2}}\sin(\omega t-\phi). \nonumber\]
Finding General Solution for Rlc Circuit
Source: https://math.libretexts.org/Bookshelves/Differential_Equations/Book%3A_Elementary_Differential_Equations_with_Boundary_Value_Problems_(Trench)/06%3A_Applications_of_Linear_Second_Order_Equations/6.03%3A_The_RLC_Circuit