Circuit Theory/Phasors/Examples/Example 8

From Wikibooks, open books for an open world
Jump to: navigation, search
Series RL circuit for example 8

Given that the current source is defined by  I_s(t) = 120 \sqrt{2} cos(377t+120^\circ), find all other voltages, currents and check power.

Label Loops Junctions[edit]

Series RL circuit marked up for analysis

The important point is that nothing changes even though the source is oscillating. + and - still have to capture the circuit topology and work their way into the equations.

Knowns, Unknowns and Equations[edit]

Knowns: I_s, R, L
Unknowns: V_s, v_R, v_L
Equations: v_R(t) = R * i(t), v_L = L * {d \over dt}i(t), v_R(t) + v_L(t) - V_s(t) = 0

Calculus Symbolic[edit]

Evaluate the terminal relations in this order:

v_R(t) = R * i(t)
v_L = L * {d \over dt}i(t)
V_s(t)= v_R(t) + v_L(t)

Calculus Numeric[edit]

phasor matlab solution of RL circuit .... m file
v_R(t) = 10 * 120 \sqrt{2} cos(377t+\frac{2\pi}{3}) = 1700.0*cos(377.0*t + 2.09)
v_L = .01 * {d \over dt}120 \sqrt{2} cos(377t+\frac{2\pi}{3})=-640.0*sin(377.0*t + 2.09)
V_s(t)= v_R(t) + v_L(t) = 1700.0*cos(377.0*t + 2.09) - 640.0*sin(377.0*t + 2.09)
phasor math avoids trig! m-file

Now the trouble is how to wrestle the above answer into a form where it can be compared to the other answers and our intuition built up. There are three possible ways:

Euler
Trig
Phasors (which is derived from Euler).

Here is the phasor method. Notice the math is done in the phasor domain ... with imaginary numbers.

1700cos(377t+2.09) \Leftrightarrow 1700cos(2.09) + j1700sin(2.09)
-640sin(377t +2.09) \Leftrightarrow -640*sin(2.09) + j640*cos(2.09)

Now the numbers are in the phasor world. The next step is to add them all up.

\mathbb{V}_s = 1814\angle2.45

Now put back into the time domain:

V_s = 1814cos(377t+2.45)

Phasors are also an alternative to the calculus above.

Because no integrals were evaluated in this calculus solution, there is no integration constant!

Phasor Symbolic[edit]

This problem could be done with Laplace transforms, but Calculus functions are less complex. It can also be done with phasors. The question is "Which mathmetical tool should be used?" Calculus above leaves the answer in a form that would require some trig to get back into the form:

V_s = V_m cos(\omega t + \phi)

The phasor solution enables us to stay close to the above solution form.

The phasor solution provides a uniform way (one tool) for doing all inductor and capacitor problems. We will use it through out the rest of the course.

time domain[edit]

v_R(t) = R * i(t)
v_L = L * {d \over dt}i(t)
v_R(t) + v_L(t) - V_s(t) = 0

phasor domain[edit]

i(t) \rightarrow \mathbb{I}
\mathbb{I} = I_m\angle\phi
\mathbb{V}_R = R*\mathbb{I} = R*I_m\angle\phi
{d \over dt}i(t)\rightarrow j\omega\mathbb{I}=\omega\mathbb{I}\angle\frac{\pi}{2}= \omega I_m \angle (\frac{\pi}{2} + \phi)
\mathbb{V}_L = jL\omega\mathbb{I} = L*\omega I_m \angle (\frac{\pi}{2} + \phi)
\mathbb{V}_S = jL\omega\mathbb{I} + R*\mathbb{I}=(R+jL\omega)\mathbb{I}=\sqrt{R^2 + (L\omega)^2}\angle arctan(\frac{L\omega}{R})*I_m\angle\phi
\mathbb{V}_S=I_m*\sqrt{R^2 + (L\omega)^2}\angle (arctan(\frac{L\omega}{R}) + \phi)

back to time domain[edit]

V_s(t) = \operatorname{Re}(\mathbb{V}_s e^{j\omega t})= I_m \sqrt{R^2 + (L\omega)^2} cos(\omega t + \phi + arctan(\frac{L\omega}{R}))
v_L(t) = \operatorname{Re}(\mathbb{V}_L e^{j\omega t}) = I_m*L *\omega*cos(\omega t + \phi+\frac{\pi}{2})
 v_R(t) = \operatorname{Re}(\mathbb{V}_R e^{j\omega t}) = R*I_m cos(\omega t + \phi)

Phasor Numeric[edit]

time domain[edit]

v_R(t) = 10*120\sqrt{2}cos(377t + \frac{2*\pi}{3}) =1700 cos(377t + 2.09)

phasor domain[edit]

Rather than substituting into the time domain as last time, showing phasor (imaginary math) by substituting into the phasor domain.

i(t) \rightarrow {\mathbb{I}} = 120\sqrt{2}\angle\frac{2*\pi}{3}
R + jL\omega = 10+.01*377j = 10.6870\angle 0.3605
{\mathbb{V}_s} = \mathbb{I}*(R + jL\omega)= 120\sqrt{2}*10.6870 \angle(\frac{2*\pi}{3}+0.3605) = 1814\angle 2.45
\mathbb{V}_L = I_m*\omega*L\angle(\phi+\frac{\pi}{2}) = 120\sqrt{2}*377*.01\angle(\frac{2*\pi}{3}+\frac{\pi}{2})=640\angle3.67
doing problem using phasors rather than calculus m-filec

back to time domain[edit]

V_s(t) =  1811 cos(377t + 2.45)
v_L(t) = 640 cos(377t + 3.67)
v_R(t) = 1700 cos(377t + 2.09)

comparison[edit]

Only V_s was calculated with both the phasor 1811 and calculus/phasor 1814 method. The 1814 number can be found in the phasor math of both the calculus method and the phasor method. But the phasor math in both sections had an intermediate calculation of R + jLw in the phasor domain. A matlab script was written for the phasor solution that merely plugged into the time domain symbolic solution. This did not allow errors to build up in matlab caused by reusing numbers .. and came up with 1811.

The goal is to type numbers into the calculator once, or into matLab once and calculate the answers directly from the given numbers with no intermediary calculations.

Simulation[edit]

Example8f.png Vs,VL and I of the simulation graphed above.

period check[edit]

The period above looks to be between 16ms and 17ms, closer to 17ms. This agrees with the formula:

\omega = 2\pi*f
f=\frac{\omega}{2\pi}
T =\frac{1}{f} = \frac{2\pi}{\omega} = \frac{2\pi}{377} = 16.7 ms

magnitude check[edit]

The magnitude of Vs above appears to be close to 2000 volts. This is close to 1811/1814 of the math.

The magnitude of i(t) appears to be 120*\sqrt{2}.

The magnitude of VL appears to be above 500 volts which could be the 640 volts from the math.

Vr can not be graphed because of the choice of ground. Notice how Vr and VL share the same -. They can be measured with an oscilloscope or the above simulator. Vr can not.

transient response check[edit]

The transient response is different. The straight lines and starting above 0 seems odd. Since this is a steady state analysis, going to save digging into this until later. Otherwise, everything is the same as when driven by a voltage source.

phase check[edit]

Voltage will always lead the current through an inductor (think of the terminal relationship or ELI) by 90^{\circ}, \frac{\pi}{2} or \frac{1}{4} of a period, whether a current source or voltage source is driving the circuit.

Voltage Check[edit]

Could do the phasor math again to make sure they check, but this is the math that produced the values in the first place. Here the goal is to do a quick spot check. This check is not good, but it is quick. It can detect many mistakes. At the least it builds our confidence in the answer. Need to plug into:

v_L = L * {d \over dt}i(t), v_R(t) + v_L(t) - V_s(t) = 0

And make sure that they add to zero.

V_s(t) =  1811 cos(377t + 2.45)
v_L(t) = 640 cos(377t + 3.67)
v_R(t) = 1700 cos(377t + 2.09)

Pick t = 1:

 V_s(t) = 1811 cos(377+2.45)= -1410
v_L(t) = 640 cos(377 + 3.67)=-560
v_R(t) = 1700 cos(377 + 2.09)= -857

The equation above is true. The actual numbers above would need more accuracy (decimal places) to get closer to zero.

Power Analysis[edit]

Power Analysis is rooted in the phasor domain!

\mathbb{V}_s =  1814\angle 2.45
\mathbb{I} = 120 \sqrt{2}\angle 2.09
\mathbb{I}^* = 120 \sqrt{2}\angle -2.09

if :\mathbb{V} = M_v\angle\phi_v\quad, and \quad\mathbb{I} = M_i\angle\phi_i then

\mathbb{S} = \mathbb{V}\mathbb{I}^* = \frac{M_vM_i}{2}\angle(\phi_v - \phi_i) = \frac{120 \sqrt{2} * 1814}{2}\angle(2.45 - 2.09)= 15400\angle 0.36=14400 + 5420j

and

cos(.36) = .936

The power is much larger ... by a factor of 100 because of the resistor, but the phase angle is still the same. This means the power factor is still the same.

Value Units Description
15400 volt-ampere va apparent power what utility companies manage: peak power they design for, peak power they have to deliver
.936 unitless power factor, ratio of real power to apparent power, ideally 1
14400 watt W real, average, active power ... what consumers want to pay for (watt-hours)
5420 volt-amp-reactive var reactive power ... why not all outlets in a room are on the same circuit breaker

Intuition[edit]

ELI ... voltage leads the current through an inductor

Derivatives cause a lag ... a delay in time ... which is a positive angle in sinusoidal.