6th Order RLC

A DC voltage source excites an RLC circuit, as shown. At time t = 0 t = 0 , the inductors and capacitors are de-energized. The quantity I S I_S is the current flowing out of the DC source. Let I S 0 I_{S0} be the value of I S I_S at time t = 0 t = 0 . Let I S I_{S \infty} be the limiting value of I S I_S as the elapsed time approaches infinity. Let I S m i n I_{S min} and I S m a x I_{S max} be the smallest and largest values of I S I_S over all time.

Determine the value of Q Q :

Q = I S 0 I S m i n ( I S m a x I S ) \large{ Q = I_{S0} \,I_{Smin} \, (I_{Smax} - I_{S \infty})}

Details and Assumptions:
1) All R L C RLC values are 1 1 . V S = 10 V_S = 10
2) All four current values are positive
3) The circuit shown is a modified version of a circuit from this web page


The answer is 20.415.

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

Steven Chase
Mar 27, 2020

I will post my solution as well, since there are minor stylistic differences relative to the solution given by @Karan Chatrath . The state space representation is:

V L 1 = V S R 1 I L 1 V C 3 = L 1 I ˙ L 1 V L 2 = V C 3 + V C 1 R 3 I L 2 = L 2 I ˙ L 2 V L 3 = V C 3 = L 3 I ˙ L 3 I R 2 = V S ( V C 3 + V C 1 ) R 2 I C 1 = I R 2 I L 2 = C 1 V ˙ C 1 I R 4 = V C 3 V C 2 R 4 I C 2 = I R 4 = C 2 V ˙ C 2 I C 3 = I L 1 + I C 1 I R 4 I L 3 = C 3 V ˙ C 3 V_{L1} = V_S - R_1 I_{L1} - V_{C3} = L_1 \dot{I}_{L1} \\ V_{L2} = V_{C3} + V_{C1} - R_3 I_{L2} = L_2 \dot{I}_{L2} \\ V_{L3} = V_{C3} = L_3 \dot{I}_{L3} \\ I_{R2} = \frac{V_S - (V_{C3} + V_{C1})}{R_2} \\ I_{C1} = I_{R2} - I_{L2} = C_1 \dot{V}_{C1} \\ I_{R4} = \frac{V_{C3} - V_{C2}}{R_4} \\ I_{C2} = I_{R4} = C_2 \dot{V}_{C2} \\ I_{C3} = I_{L1} + I_{C1} - I_{R4} - I_{L3} = C_3 \dot{V}_{C3}

Isolate the time derivatives of state variables ( I L 1 , I L 2 , I L 3 , V C 1 , V C 2 , V C 3 ) (I_{L1}, I_{L2}, I_{L3}, V_{C1}, V_{C2}, V_{C3}) and numerically integrate. Then the source current is:

I S = I L 1 + I R 2 I_S = I_{L1} + I_{R2}

@Steven Chase Sir how are you? Are you alright? I have seen the cases in America are increasing very fastly.

A Former Brilliant Member - 1 year, 2 months ago

I'm doing well, thanks. Hopefully you are too.

Steven Chase - 1 year, 2 months ago

Log in to reply

Yes, I am also doing well. Thanks

A Former Brilliant Member - 1 year, 2 months ago
Karan Chatrath
Mar 27, 2020

Currents and charges on the capacitor are as shown in the diagram. As charge accumulates in the capacitors, one obtains:

Q ˙ 1 = I 1 ; Q ˙ 2 = I 2 ; Q ˙ 3 = I 3 \dot{Q}_1 = I_1 \ ; \ \dot{Q}_2 = I_2 \ ; \ \dot{Q}_3= I_3

From Kirchoff's current law:

I = I 6 + I 7 I = I_6 + I_7 I 7 = I 1 + I 5 I_7 = I_1 + I_5 I 1 + I 6 = I 2 + I 8 I_1 + I_6 = I_2 + I_8 I 8 = I 3 + I 4 I_8 = I_3 + I_4

From Kirchoff's voltage law, ignoring all passive element notation, as they all have a value of unity, one gets:

V S + I 7 + I 5 + I ˙ 5 = 0 -V_S + I_7 + I_5 + \dot{I}_5 = 0 I ˙ 6 + I 6 = I 7 + Q 1 \dot{I}_6 + I_6 = I_7 + Q_1 I ˙ 4 = I 2 + Q 2 \dot{I}_4 = I_2 + Q_2 I ˙ 4 = Q 3 \dot{I}_4 = Q_3 Q 1 + I 2 + Q 2 = I 5 + I ˙ 5 Q_1 + I_2 + Q_2 = I_5 + \dot{I}_5

This system of equations can be re-arranged to obtain a sixth-order state-space form. However, this re-arrangement has been avoided. These differential-algebraic equations are numerically integrated to obtain the required result which is as follows:

The required answer turns out to be:

Q 20.4146 \boxed{Q \approx 20.4146}

@Karan Chatrath What's the difference between code and numerical integration. I am very excited to do numerical integration.

A Former Brilliant Member - 1 year, 2 months ago

@Karan Chatrath How you do numerical integration??

A Former Brilliant Member - 1 year, 2 months ago

Log in to reply

Numerical integration comprises of techniques used to solve integrals and differential equations, especially those which do not have exact solutions, on in cases when solving for the exact solution is time-intensive, such as this problem. I perform numerical integration by writing a few lines of code. So while solving double/triple integrals or the above circuit, I numerically integrate by writing code. In this context, code and numerical integration are terms I use interchangeably.

As for how I perform it, you can refer to any of the solutions posted earlier. The theory of numerical integration can be easily found all over the internet.

Karan Chatrath - 1 year, 2 months ago

Log in to reply

@Karan Chatrath Sir which rule did you use in numerical integration like Simmons/3rd rule, trapezoid rule???

A Former Brilliant Member - 1 year, 2 months ago

Log in to reply

@A Former Brilliant Member Both Trapezoidal and Simpson's rule are useful methods. I frequently use the Explicit Euler integration method.

Karan Chatrath - 1 year, 2 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...