Second Order Differential Equations

3.9 Application: RLC Electrical Circuits

In Section 2.5F, we explored first-order differential equations for electrical circuits consisting of a voltage source with either a resistor and inductor (RL) or a resistor and capacitor (RC). Now, equipped with the knowledge of solving second-order differential equations, we are ready to delve into the analysis of more complex RLC circuits, which incorporate a resistor, inductor, and capacitor.

Previously, we established that:

  • Ohm’s law dictates that the voltage drop [asciimath]E_R[/asciimath] ​ across a resistor is proportional to the current I flowing through it, expressed as  [asciimath]E_R=RI[/asciimath], where [asciimath]R[/asciimath]  is the resistance.
  • Faraday’s law, complemented by Lenz’s law, describes that the voltage drop [asciimath]E_L[/asciimath] ​across an inductor is proportional to the rate of change of current, given as [asciimath]E_L=L(dI)/(dt)[/asciimath], where [asciimath]L[/asciimath] is the inductance.
  • The voltage drop [asciimath]E_C[/asciimath] ​across a capacitor is proportional to the electric charge [asciimath]q[/asciimath] stored on it, represented as [asciimath]E_C=1/Cq[/asciimath], with [asciimath]C[/asciimath]  being the capacitance

Schematic of an RLC circuit consisting of a resistor, a capacitor, and an inductor in a series circuit

Figure 3.9.1 Schematic of an RLC series circuit

With these foundations, consider [asciimath]E(t)[/asciimath] as the external voltage supplied to the RLC series circuit in Fig. 3.9.1. By applying Kirchhoff’s voltage law, we have

 [asciimath]E_L+E_R+E_c=E(t)[/asciimath]

Substituting [asciimath]E_R=RI[/asciimath] ,[asciimath]E_L=L (dI)/dt[/asciimath] , and [asciimath]E_C=1/Cq[/asciimath] into this equation yields

 [asciimath]L (dI)/dt+RI+1/Cq=E(t)[/asciimath](3.9.1)

Differentiating this equation with respect to time and substituting [asciimath]I=(dq)/dt[/asciimath] transform it into a second-order differential equation.

 [asciimath]L(d^2I)/(dt^2)+R (dI)/dt+1/CI=(dE)/dt[/asciimath](3.9.2)

Alternatively, Equation 3.9.1 can be expressed in terms of charge [asciimath]q(t)[/asciimath].

 [asciimath]L (d^2q)/(dt^2)+R(dq)/(dt)+1/Cq=E(t)[/asciimath](3.9.3)

Given [asciimath]E(t)[/asciimath] and an initial condition, such as initial current[asciimath]I(0)[/asciimath] and initial charge [asciimath]q(0)[/asciimath], we can solve the equation for [asciimath]I(t)[/asciimath] using techniques discussed in previous sections, such as the method of undetermined coefficients. Once [asciimath]I(t)[/asciimath] is determined, the voltage across different components of the circuit can be calculated.

 

Example 3.9.1: RLC Series Circuit

Consider an RLC series circuit with a resistor of [asciimath]0.06\ Omega[/asciimath] and an inductor of [asciimath]0.01\ "H"[/asciimath], and a capacitor of [asciimath]50/89\ "F"[/asciimath] powered by a voltage source [asciimath]E(t)=0.1sin(10t)\ "V"[/asciimath]. Initially, the current and charge on the capacitor are zero. Determine the current in the circuit as a function of time.

Show/Hide Solution

 

Given information:

  • Resistor: [asciimath]R=0.06\ Omega[/asciimath]
  • Inductor: [asciimath]L=0.01 \ "H"[/asciimath]
  • Capacitor: [asciimath]C=50/89\ "F"[/asciimath]
  • Voltage source: [asciimath]E(t)=0.1sin(10t)\ "V"[/asciimath]
  • Initial current on capacitor: [asciimath]I(0)=0 \ "A"[/asciimath]
  • Initial charge on capacitor: [asciimath]q(0)=I'(0)=0 \ C[/asciimath]

 

The differential equation for an RLC series circuit is given by Equation 3.9.1.

 [asciimath]L(d^2I)/(dt^2)+R (dI)/dt+1/CI=(dE)/dt[/asciimath]

The initial value problem is then

  [asciimath]0.01(d^2I)/(dt^2)+0.06 (dI)/dt+89/50I=cos(10t),[/asciimath]   [asciimath]I(0)=0, \ I'(0)=0[/asciimath]

Multiplying the equation by 100, we get

 [asciimath](d^2I)/(dt^2)+6 (dI)/dt+178I=100cos(10t),[/asciimath]    [asciimath]I(0)=0, \ I'(0)=0[/asciimath]

Given the characteristic equation has complex conjugates [asciimath]r_(1,2)=-3+-13i[/asciimath], the complementary solution is

 [asciimath]I_c(t)=e^(-3t)(c_1cos(13t)+c_2sin(13t))[/asciimath]

 

Finding the particular solution:

To find the particular solution, we use undetermined coefficients. Given the forcing cosine function, we guess the form of the particular solution to be

 [asciimath]I_p=Acos(10t)+Bsin(10t)[/asciimath]

The derivatives are

 [asciimath]I'_p=-10Asin(10t)+10Bcos(10t)[/asciimath]

 [asciimath]I''_p=-100Acos(10t)-100Bsin(10t)[/asciimath]

Substituting  [asciimath]I_p[/asciimath] and its derivatives into the differential equation yields

 [asciimath]-100Acos(10t)-100Bsin(10t) +6(-10Asin(10t)+10Bcos(10t) )[/asciimath] [asciimath]+178(Acos(10t)+Bsin(10t))=[/asciimath] [asciimath]100cos(10t)[/asciimath]

Simplifying gives

 [asciimath](78A+60B)cos(10t)+(-60A+78B)sin(10t)=100cos(10t)[/asciimath]

By matching coefficients of sine and cosine terms and solving the system of two equations in unknowns [asciimath]A[/asciimath] and [asciimath]B[/asciimath], we get

[asciimath]A=650/807,\ B=500/807[/asciimath]

Therefore, the particular solution is

 [asciimath]I_p=650/807cos(10t)+500/807sin(10t)[/asciimath]

 

Combining the particular and complementary solutions gives the general solution

 [asciimath]I(t)=650/807cos(10t)+500/807sin(10t) +[/asciimath] [asciimath]e^(-3t)(c_1cos(13t)+c_2sin(13t))[/asciimath]

Applying the initial conditions:

 [asciimath]I_0(0)=0 \ ->[/asciimath]   [asciimath]c_1=-650/807[/asciimath]

 [asciimath]I'_0(0)=0\ ->[/asciimath] [asciimath]c_2=-6950/10491[/asciimath]

The equation of the object’s displacement is then

 [asciimath]I(t)=650/807cos(10t)+500/807sin(10t) +[/asciimath] [asciimath]e^(-3t)(-650/807 cos(13t)-6950/10491 sin(13t))[/asciimath]

As with forced mechanical vibration scenarios, the current in an RLC circuit is composed of two distinct parts: the transient current, represented by the complementary solution that diminishes to zero as time progresses to infinity, and the steady-state current, described by the particular solution which is sinusoidal and persists over time.

 

Try an Example

 

Section 3.9 Exercises

  1. Consider an RLC circuit with a [asciimath]17/50\Omega[/asciimath] resistor, a [asciimath]1/100\H[/asciimath] inductor, and a [asciimath]100/93\F[/asciimath] capacitor driven by the voltage [asciimath]E(t) = 0.06sin(3t)V[/asciimath]. a) Write the differential equation associated with this circuit in terms of current [asciimath]I[/asciimath] . b) If the initial charge and initial current on the capacitor are both zero, find the current [asciimath]I[/asciimath] and the voltages across the resistor [asciimath]E_R[/asciimath] in terms of time [asciimath]t[/asciimath].
    Show/Hide Answer

    a) [asciimath]I''+34I'+93I = 18cos(3t)[/asciimath]

    b) [asciimath]I(t)=0.0205 e^(-31t)-0.1071 e^(-3t)+[/asciimath] [asciimath]0.1052sin(3t)+0.08660cos(3t)[/asciimath]

    c) [asciimath]E_R(t)=0.34(0.0205 e^(-31t)-0.1071 e^(-3t)+[/asciimath] [asciimath]0.1052sin(3t)+0.08660cos(3t) )[/asciimath]

  2. Consider an RLC circuit with a [asciimath]1/50\Omega[/asciimath] resistor, a [asciimath]1/100\H[/asciimath] inductor, and a [asciimath]50/61\F[/asciimath] capacitor driven by the voltage [asciimath]E(t) = 0.09t^2\ V[/asciimath]. a) Write the differential equation associated with this circuit in terms of current [asciimath]I[/asciimath]. b) If the initial charge and initial current on the capacitor are both zero, find the current [asciimath]I[/asciimath].
    Show/Hide Answer

    a) [asciimath]I''+2I'+122I = 18t[/asciimath]

    b) [asciimath]I(t)=e^(-t)(9/3721cos(11t)-540/40931sin(11t))+9/61t-9/3721[/asciimath]

License

Icon for the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License

Differential Equations Copyright © 2024 by Amir Tavangar is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License, except where otherwise noted.

Share This Book