Engineering Finding Coefficients and Ratio of Over-damped Circuit Equation for RLC Circuit

[ResonantW]

Homework Statement

For an overdamped series RLC circuit, determine the coefficients \beta_1 and \beta_2 in the equation V=Ae^{-\beta_1t}+Be^{-\beta_2t} for the case where R=600 \Omega, L=100 \mu H, and C=.01 \mu F. Also determine the ratio of B to A.

Homework Equations

For a series RLC circuit, the general differential equation is \large \frac{d^2V}{dt^2} + \frac{R}{L} \frac{dV}{dt} + \frac{1}{LC} V = 0.

The Attempt at a Solution

I think I can just assume that V=Ae^{-\beta t}, take the necessary derivatives, and then solve the quadratic equation that results for \beta_1 , \beta_2.

After canceling the exponential terms, I get \beta^2 - \beta \frac{R}{L} + \frac{1}{LC} = 0.

The quadratic equation is \large \frac{\frac{R}{L} \pm \sqrt{\frac{R^2}{L^2} - \frac{4}{LC}}}{2}.

My numbers end up being \beta_1 = 5.828 * 10^6 sec, \beta_2 = 1.715*10^5 sec. Is this the correct answer? I don't know why they would be an entire order of magnitude separate...

Was it ok to assume the solution was not the sum of exponentials, but a single exponential which I take both roots of? I don't have much experience with differential equations.

Then for the A to B ratio, I know that \frac{dV}{dt} = 0, so I can use the \beta 's and I end up getting \frac{A}{B} = \frac{-\beta_2}{\beta_1}... is that right as well? I end up getting a huge ratio, like -33.98.

Thanks for the help!

 

[rude man]

Having problems here:
1. What is V?
2. I suspect you mean a parallel, not a series, RLC circuit, for which the equation is
d²V/dt² + (1/RC)dV/dt + V/LC= 0.
Got a diagram of your circuit handy?

 

[rude man]

Assuming L/R = RC for some reason having to do with your definition of "overdamped", what you did to get the two betas looks right, except the units of beta is 1/sec, not sec.

Anyway, without a diagram including indication of how the initial conditions (capacitor voltage and inductor current at t = 0+) are determined I don't know what else to say.
I don't know where you got dV/dt = 0 at t = 0+ or any other time.
]

 

[ResonantW]

V refers to the voltage across the capacitor. It really is a series circuit. The differential equation I wrote is exactly correct for a series RLC (undriven) circuit. The circuit is given in Purcell's textbook, Figure 8.4a on page 303. It is essentially a series RLC with charged capacitor and switch closed at time t=0. So there is no need to say L/R=RC, because the RC comes from the parallel differential equation with which are not concerned here.

So the voltage is instantaneously equal to the initial voltage of the capacitor at the t=0, and at t=0 dV/dt (proportional to current) is 0 because the inductor prevents any flow of current at that instant.

 

[rude man]

p;pp

V refers to the voltage across the capacitor. It really is a series circuit. The differential equation I wrote is exactly correct for a series RLC (undriven) circuit. The circuit is given in Purcell's textbook, Figure 8.4a on page 303. It is essentially a series RLC with charged capacitor and switch closed at time t=0. So there is no need to say L/R=RC, because the RC comes from the parallel differential equation with which are not concerned here.

So the voltage is instantaneously equal to the initial voltage of the capacitor at the t=0, and at t=0 dV/dt (proportional to current) is 0 because the inductor prevents any flow of current at that instant.

My gleanings from your post: C is pre-charged, then the switch throws C across the series L-R network. The current thru L zero at time t = 0. If this is the case the results you got for A/B agree with mine.

I can't verify your values of β₁ and β₂. The reason, which probably means nothing to you but might to some of my Laplace-savvy friends, is that the Laplace of v(t) is V(s) = v(0)(s + R/L)/[s+a)(s+b)] where (s+a)(s+b) = s² + (R/L)s + 1/LC. In other words, I can't vouch for Aexp(-β₁t) + Bexp(-β₂t) as being the correct solution without inverse-transforming V(s) which I am not in a position to do now. However, if you were given the correct form of v(t) then what you did looks all right, ignoring the possibility of numerical errors. Your equation for finding β₁ and β₂ (the quadrature in β) would be correct.

 

[aralbrec]

My numbers end up being \beta_1 = 5.828 * 10^6 sec, \beta_2 = 1.715*10^5 sec. Is this the correct answer? I don't know why they would be an entire order of magnitude separate...

The roots of a second order equation like this will lie on a semicircle of radius 1/√(LC) in the left half of the complex plane when the roots are complex. Keeping LC constant but manipulating other components (R) will cause the roots to move, following a fixed path. Just as the roots leave the semicircle and become real, they will both be located at -1/√(LC) = 1e6 /s on the -ve real axis in this case. This happens when the stuff under the square root in your β equation equals zero. After this, one root will move toward zero on the real axis (but never get there) and the other root will move toward -∞. You can see you found one root to the right of 1e6 and one root to the left. The root that moves toward -∞ is dominated by the root that moves toward zero. This is because the exponential term for the root toward -∞ settles more quickly.

Was it ok to assume the solution was not the sum of exponentials, but a single exponential which I take both roots of? I don't have much experience with differential equations.

What you did was right. You assumed the solution was exponential and found two exponentials that will satisfy the equation. The differential equation is second order so there must be two linearly independent functions in the solution.

Then for the A to B ratio, I know that \frac{dV}{dt} = 0, so I can use the \beta 's and I end up getting \frac{A}{B} = \frac{-\beta_2}{\beta_1}... is that right as well? I end up getting a huge ratio, like -33.98.

That's right too. You will still need one more initial condition to find values for A and B.

 

[rude man]

For what it's worth, I finished my computations. I agreed with your characteristic equation. But now for A/B I got
A = E(R/L-β₁₎)/(β₂-β₁)
B = E(β₂-R/L)/(β₂-β₁)
so that A/B = (R/L-β₁)/(β₂-R/L).
Note that this is the complete solution to the problem assuming initial L current = 0 ad initial C voltage = E.