6th-Order LC Circuit

The circuit below consists of three capacitors ( C 1 = 2 , C 2 = 2 , C 3 = 3 ) (C_1 = 2, C_2 = 2, C_3 = 3) and three inductors ( L 1 = 2 , L 2 = 3 , L 3 = 1 ) (L_1 = 2, L_2 = 3, L_3 = 1) . At time t = 0 t = 0 the inductor currents are zero and the capacitor voltages are ( V C 1 = + 5 , V C 2 = + 7 , V C 3 = 12 ) (V_{C1} = +5, V_{C2} = +7, V_{C3} = -12) .

At time t = 4 t = 4 , what fraction of the total system energy is contained within the three capacitors combined?

Inspiration

Note: In the drawing, the plus and minus signs indicate voltage polarities, and the arrows indicate current polarities.


The answer is 0.15.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

2 solutions

Karan Chatrath
Jun 6, 2020

Slightly different approach from that of @Hosam Hajjir

The circuit equation by applying Kirchoff's current and voltage laws read:

Q ˙ 3 = I C 3 \dot{Q}_3 = I_{C3} Q ˙ 2 = I C 2 \dot{Q}_2 = I_{C2} Q ˙ 1 = I C 1 \dot{Q}_1 = I_{C1}

I C 3 = I C 2 + I L 3 I_{C3} = I_{C2} + I_{L3} I C 1 = I C 3 + I L 1 I_{C1} = I_{C3} + I_{L1} I C 2 = I C 1 + I L 2 I_{C2} = I_{C1} + I_{L2}

Q 3 C 3 + L 3 I ˙ L 3 = L 1 I ˙ L 1 \frac{Q_3}{C_3} + L_3 \dot{I}_{L3} = L_1 \dot{I}_{L1} Q 2 C 2 + L 2 I ˙ L 2 = L 3 I ˙ L 3 \frac{Q_2}{C_2} + L_2 \dot{I}_{L2} = L_3 \dot{I}_{L3} Q 1 C 1 + L 1 I ˙ L 1 = L 2 I ˙ L 2 \frac{Q_1}{C_1} + L_1 \dot{I}_{L1} = L_2 \dot{I}_{L2}

Taking Laplace transform on both sides of all equation leads to the following linear equations in the s-domain:

s Q 3 ( s ) C 3 V C 3 = I C 3 ( s ) sQ_3(s) - C_3V_{C3} = I_{C3}(s) s Q 2 ( s ) C 2 V C 2 = I C 2 ( s ) sQ_2(s) - C_2V_{C2} = I_{C2}(s) s Q 1 ( s ) C 1 V C 1 = I C 1 ( s ) sQ_1(s) - C_1V_{C1} = I_{C1}(s)

I C 3 ( s ) = I C 2 ( s ) + I L 3 ( s ) I_{C3}(s) = I_{C2}(s) + I_{L3}(s) I C 1 ( s ) = I C 3 ( s ) + I L 1 ( s ) I_{C1}(s) = I_{C3}(s) + I_{L1}(s) I C 2 ( s ) = I C 1 ( s ) + I L 2 ( s ) I_{C2}(s) = I_{C1}(s) + I_{L2}(s)

Q 3 ( s ) C 3 + s L 3 I L 3 ( s ) = s L 1 I L 1 ( s ) \frac{Q_3(s)}{C_3}+ sL_3 I_{L3}(s)=sL_1 I_{L1}(s) Q 2 ( s ) C 2 + s L 2 I L 2 ( s ) = s L 3 I L 3 ( s ) \frac{Q_2(s)}{C_2}+ sL_2 I_{L2}(s)=sL_3 I_{L3}(s) Q 1 ( s ) C 1 + s L 1 I L 1 ( s ) = s L 2 I L 2 ( s ) \frac{Q_1(s)}{C_1}+ sL_1 I_{L1}(s)=sL_2 I_{L2}(s)

Plugging in all constants and solving for Q 1 ( s ) Q_1(s) , Q 2 ( s ) Q_2(s) and Q 3 ( s ) Q_3(s) gives:

Q 1 ( s ) = 160 s 16 s 2 + 1 Q_1(s) = \frac{160\,s}{16\,s^2+1} Q 2 ( s ) = 4 s ( 616 s 2 + 31 ) ( 11 s 2 + 1 ) ( 16 s 2 + 1 ) Q_2(s) = \frac{4\,s\,\left(616\,s^2+31\right)}{\left(11\,s^2+1\right)\,\left(16\,s^2+1\right)} Q 3 ( s ) = 6 s ( 1056 s 2 + 71 ) ( 11 s 2 + 1 ) ( 16 s 2 + 1 ) Q_3(s) = -\frac{6\,s\,\left(1056\,s^2+71\right)}{\left(11\,s^2+1\right)\,\left(16\,s^2+1\right)}

Using an online inverse Laplace transform calculator:

Q 1 ( t ) = 10 cos ( t 4 ) Q_1(t) = 10\,\cos\left(\frac{t}{4}\right) Q 2 ( t ) = 20 cos ( 11 t 11 ) 6 cos ( t 4 ) Q_2(t) = 20\,\cos\left(\frac{\sqrt{11}\,t}{11}\right)-6\,\cos\left(\frac{t}{4}\right) Q 3 ( t ) = 6 cos ( t 4 ) 30 cos ( 11 t 11 ) Q_3(t) = -6\,\cos\left(\frac{t}{4}\right)-30\,\cos\left(\frac{\sqrt{11}\,t}{11}\right)

The energy stored in the capacitors at any instant is:

E C ( t ) = Q 1 2 2 C 1 + Q 2 2 2 C 2 + Q 3 2 2 C 3 E_C(t) = \frac{Q_1^2}{2C_1} + \frac{Q_2^2}{2C_2} + \frac{Q_3^2}{2C_3}

At time t = 0 t=0 there is no current flowing through inductors, this means that the total energy of the system at that instant is stored in the capacitors. Therefore, the required answer is:

E C ( 4 ) E C ( 0 ) = 0.15 \boxed{\frac{E_C(4)}{E_C(0)} = 0.15}

@Karan Chatrath Y E A H ! \textcolor{#20A900}{YEAH! } Everyone is a legend until the real legend comes.

I don't know why why why. But sometimes your solution seems me as an ideal solution for a problem like this one .

Log in to reply

Thank you! Truth is that this solution is not ideal. Solving a 9 equation system is very tedious and time-consuming. I have left out those steps. Deriving a state-space form and numerical integration is a less tedious method to solve this problem.

Karan Chatrath - 1 year ago
Hosam Hajjir
Jun 5, 2020

Writing the loop equations, and using the voltages across the capacitors and the currents through the inductors as our variables, we get,

V C 1 L 2 i L 2 d t + L 1 d i L 1 d t = 0 ( 1 ) V_{C1} - L_2 \dfrac{i_{L2}}{dt} + L_1 \dfrac{di_{L1}}{dt} = 0 \hspace{12pt} (1)

V C 2 L 3 i L 3 d t + L 2 d i L 2 d t = 0 ( 2 ) V_{C2} - L_3 \dfrac{i_{L3}}{dt} + L_2 \dfrac{di_{L2}}{dt} = 0 \hspace{12pt} (2)

V C 3 L 1 i L 1 d t + L 3 d i L 3 d t = 0 ( 3 ) V_{C3} - L_1 \dfrac{i_{L1}}{dt} + L_3 \dfrac{di_{L3}}{dt} = 0 \hspace{12pt} (3)

We now note that

V C 3 = V C 1 V C 2 ( 4 ) V_{C3} = -V_{C1} - V_{C2} \hspace{12pt} (4)

and, with this last equation, if we add, equations (1) through (3), we get an identity, which implies that one of the equations is linearly dependent on the

other two. So we can eliminate the third equation and the variable V C 3 V_{C3} as it is linearly dependent on V C 1 V_{C1} and V C 2 V_{C2} .

further,

i L 3 = i L 1 i L 2 ( 5 ) i_{L3} = -i_{L1} - i_{L2} \hspace{12pt} (5)

Substituting this into equation (2),

V C 2 + L 3 ( i L 1 d t + i L 2 d t ) + L 2 d i L 2 d t = 0 ( 6 ) V_{C2} + L_3 (\dfrac{i_{L1}}{dt} + \dfrac{i_{L2}}{dt}) + L_2 \dfrac{di_{L2}}{dt} = 0 \hspace{12pt} (6)

equations (1) and (6) can be solved for d i L 1 d t \dfrac{di_{L1}}{dt} and d i L 2 d t \dfrac{di_{L2}}{dt} in terms of V C 1 V_{C1} and V C 2 V_{C2} .

Define the matrix,

L = [ L 1 L 2 L 3 L 2 L 3 ] L = \begin{bmatrix} -L_1 && L_2 \\ -L_3 && -L2 - L3 \end {bmatrix}

then

[ d i L 1 d t d i L 2 d t ] = L 1 [ V C 1 V C 2 ] \displaystyle \begin{bmatrix} \dfrac{di_{L1}}{dt} \\ \dfrac{di_{L2}}{dt} \end{bmatrix} = L^{-1} \begin{bmatrix} V_{C1} \\ V_{C2} \end{bmatrix}

Next, we have the following equations,

i C 1 = C 1 d V C 1 d t = i L 1 + i C 3 ( 7 ) i_{C1} = C_1 \dfrac{dV_{C1}}{dt} = i_{L1} + i_{C3} \hspace{12pt} (7)

i C 2 = C 2 d V C 2 d t = i C 1 + i L 2 ( 8 ) i_{C2} = C_2 \dfrac{dV_{C2}}{dt} = i_{C1} + i_{L2} \hspace{12pt} (8)

i C 3 = C 3 d V C 3 d t = i L 3 + i C 2 ( 9 ) i_{C3} = C_3 \dfrac{dV_{C3}}{dt} = i_{L3} + i_{C2} \hspace{12pt} (9)

Using equation (5), which is,

i L 3 = i L 1 i L 2 ( 5 ) i_{L3} = -i_{L1} - i_{L2} \hspace{12pt} (5)

and again, we can eliminate i L 3 i_{L3} from our variables, as it is linearly dependent on i L 1 i_{L1} and i L 2 i_{L2} . So, we can re-write equation (7) as follows,

C 1 d V C 1 d t = i L 1 C 3 ( d V C 1 d t + d V C 2 d t ) ( 10 ) C_1 \dfrac{dV_{C1}}{dt} = i_{L1} - C_3 (\dfrac{dV_{C1}}{dt} + \dfrac{dV_{C2}}{dt}) \hspace{12pt} (10)

also, we can write equation (8) as follows,

C 2 d V C 2 d t = C 1 d V C 1 d t + i L 2 ( 11 ) C_2 \dfrac{dV_{C2}}{dt} = C1 \dfrac{dV_{C1}}{dt} + i_{L2} \hspace{12pt} (11)

The last two equations can be solved for d V C 1 d t \dfrac{dV_{C1}}{dt} and d V C 2 d t \dfrac{dV_{C2}}{dt} in terms of i L 1 i_{L1} and i L 2 i_{L2}

Define the matrix C C as follows

C = [ C 1 + C 3 C 3 C 1 C 2 ] C = \begin{bmatrix} C_1 +C_3 && C_3 \\ -C_1 && C_2 \end{bmatrix}

then it follows that

[ d V C 1 d t d V C 2 d t ] = C 1 [ i L 1 i L 2 ] \displaystyle \begin{bmatrix} \dfrac{dV_{C1}}{dt} \\ \dfrac{dV_{C2}}{dt} \end{bmatrix} = C^{-1} \begin{bmatrix} i_{L1} \\ i_{L2} \end{bmatrix}

Now we have the system fully specified, so if we define the vector x x as follows,

x = [ i L 1 i L 2 V C 1 V C 2 ] x = \begin{bmatrix} i_{L1} \\ i_{L2} \\ V_{C1} \\ V_{C2} \end{bmatrix}

Then, the dynamics of the system are given by,

x ˙ = A x \dot{x} = A x , where the matrix A A is given by,

A = [ 0 2 × 2 L 1 C 1 0 2 × 2 ] A = \begin{bmatrix} 0_{2\times2} && L^{-1} \\ C^{-1} && 0_{2 \times 2} \end{bmatrix}

This coupled with x ( 0 ) = [ 0 , 0 , 5 , 7 ] T x(0) = [0, 0, 5, 7]^T

This system can be solved using a numerical integration method such as the 4th-order Runge-Kutta method (a widely known method), from t = 0 t = 0 to t = 4 t = 4 . The required fraction is computed based on the fact that the energy stored in a capacitor is given by 1 2 C V 2 \frac{1}{2} C V^2 .

Therefore,

Fraction of Energy = 1 2 C 1 V C 1 2 ( 4 ) + 1 2 C 2 V C 2 2 ( 4 ) + 1 2 C 3 ( V C 1 ( 4 ) + V C 2 ( 4 ) ) 2 1 2 C 1 ( 5 ) 2 + 1 2 C 2 ( 7 ) 2 + 1 2 C 3 ( 12 ) 2 \text{Fraction of Energy} = \dfrac{ \frac{1}{2} C_1 V_{C1}^2(4) + \frac{1}{2} C_2 V_{C2}^2(4) + \frac{1}{2} C_3( V_{C1}(4)+V_{C2}(4))^2 }{ \frac{1}{2} C_1 (5)^2 + \frac{1}{2} C_2 (7)^2+\frac{1}{2} C_3 (12)^2 }

This fraction comes out to be about 0.15 0.15 . In summary the dynamics of this circuit is described by a 4-th order differential equation which is implemented using the state space approach by a first order differential equation (with constant coefficients) of the state vector with has dimension 4.

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...