Circuit Theory/Phasors/Examples/Example 10

Parallel RL circuit for example 10

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

Label Loops Junctions

Parallel RL circuit marked up for analysis

In a parallel circuit, all devices have the exact same voltage across them.

Knowns, Unknowns and Equations

Knowns: $I_s,R,L$

Unknowns: $V_s,i_R,i_L$

Equations: $V_s(t)=R*i_R(t),V_s=L*\frac{d}{dt}{i}_L(t),i_R(t)+i_L(t)-I_s(t)=0$

Phasor Symbolic

mupad attempt to do some of phasor algebra ... need to find a way to solve symbolically for the real and imaginary parts in rectangular form and polar form, plus find a way to numerically evaluate .. m-file

time domain

Evaluate the terminal relations in this order:

$i_R(t)=\frac{V_s}{R}$

$i_L(t)=\int \frac{V_s}{L}dt$

Substitute into this equation:

$I_s(t)=i_R(t)+i_L(t)$

$I_s(t)=\frac{V_s}{R}+\int \frac{V_s}{L}dt$

Need to formulate as a differential equation (to keep mathematicians happy) so differentiate all sides (don't let them see you doing this):

$\frac{d{I}_s(t)}{dt}=\frac{1}{R}\frac{d{V}_s}{dt}+\frac{V_s}{L}$

So have a differential equation that can be solved for $V_s$.

Phasor domain

now transform both operations into the phasor domain ..

$I_s(t)\to {𝕀}_s\to I_m\mathrm{\angle }\varphi$

$V_s(t)\to {𝕍}_s$

$\int V_s(t)dt\to \frac{1}{j\omega }{𝕍}_s$

So:

${𝕀}_s=\frac{{𝕍}_s}{R}+\frac{{𝕍}_s}{j\omega L}$ or $j\omega {𝕀}_s=j\omega \frac{{𝕍}_s}{R}+\frac{{𝕍}_s}{L}$

Solving for ${𝕍}_s$

${𝕍}_s=\frac{{𝕀}_{s}j\omega LR}{R+j\omega L}$

Putting in Rectangular form:

${𝕍}_s={𝕀}_s*\frac{{\omega }^2{L}^{2}R+j\omega L{R}^2}{R^2+{\omega }^2{L}^2}$

Putting in Polar form:

${𝕍}_s={𝕀}_s\frac{\sqrt{({\omega }^2{L}^{2}R{)}^2+(\omega L{R}^2{)}^2}}{R^2+{\omega }^2{L}^2}\mathrm{\angle }\arctan (\frac{R}{\omega L})$

back to time domain

$V_s(t)=I_m\frac{\sqrt{({\omega }^2{L}^{2}R{)}^2+(\omega L{R}^2{)}^2}}{R^2+{\omega }^2{L}^2}\cos (\varphi +\arctan (\frac{R}{\omega L}))$

There will be an integration constant, but it is impossible to calculate now.

This is from a differential equation.

Will calculate the integration constant after adding the homogeneous solution to the above particular solution.

Phasor Numeric

phasor numeric solution combines all phasors in rectangular form then goes straight to the polar answer form. Angle is in the third quadrant... m-file

time domain

Evaluate the terminal relations in this order:

$i_R(t)=\frac{V_s}{10}$

$i_L(t)=\int \frac{V_s}{.01}dt$

Substitute into this equation:

$I_s(t)=i_R(t)+i_L(t)$

So have a differential equation that can be solved for $V_s$:

$I_s(t)=\frac{V_s}{10}+\int \frac{V_s}{.01}dt$

Phasor domain

plugging into the symbolic time domain solution created from phasor analysis ... same angle in the third quadrant ... positive this time m-file

Choosing to do calculation in time domain

back to time domain

$V_s(t)=120*\sqrt{2}\frac{\sqrt{(377^2*.01^2*10{)}^2+(377*.01*10^2{)}^2}}{10^2+377^2*.01^2}\cos (377t+\frac{2*\pi }{3}+\arctan (\frac{10}{377*.01}))$

$V_s(t)=599\cos (377t+3.30)$

$i_R(t)=59.9\cos (377t+3.30)$

$i_L(t)=159\sin (377t+3.30)=159\cos (377t+1.74)$

Laplace Symbolic

time domain

Evaluate the terminal relations in this order:

$i_R(t)=\frac{V_s}{R}$

$i_L(t)=\int \frac{V_s}{L}dt$

Substitute into this equation:

$I_s(t)=i_R(t)+i_L(t)$

So have a differential equation that can be solved for $V_s$:

$I_s(t)=\frac{V_s}{R}+\int \frac{V_s}{L}dt$

laplace domain

now transform both into the laplace domain ..

$I_s(t)\to ℒ\{I_s(t)\}$

$V_s(t)\to ℒ\{V_s(t)\}$

$\int V_s(t)dt\to \frac{1}{s}ℒ\{V_s(t)\}$

So:

$ℒ\{I_s(t)\}=\frac{ℒ\{V_s(t)\}}{R}+\frac{ℒ\{V_s(t)\}}{sL}$

Solving for $ℒ\{V_s(t)\}$

$ℒ\{V_s(t)\}=$$\frac{ℒ\{I_s(t)\}LsR}{R+sL}$

At this point the Laplace symbolic solution has to stop. The next steps depend upon the form of $ℒ\{I_s(t)\}$.

back to time domain

This can not be done without $ℒ\{I_s(t)\}$. The form of the function $ℒ\{I_s(t)\}$ determines the inverse mapping.

Laplace Numeric

laplace solution has the same problems as before ... constants ... sin term with wrong sign m-file

Time Domain

$I_s(t)=120\sqrt{2}cos(377t+2.09)$

$R=10$

$L=.01$

Laplace Domain

$I_s(t)\to ℒ\{I_s(t)\}=ℒ\{120\sqrt{2}cos(377t+120^{\circ })\}=-\frac{(170.0*(0.5*s+326.0))}{(s^2+142129)}$

$ℒ\{V_s(t)\}=\frac{\frac{-17.0*s*(0.5*s+326.0)}{(s^2+142129)}}{10+.01*s}$

Time Domain

$V_s(t)={ℒ}^{-1}\{i(t)\}=(-2.58)*e^{-\frac{t}{.001}}+97.2sin(377t)-591cos(377*t)$

Comparison of answers

$V_s(t)=599\cos (377t-2.98)$ ... from numeric solution out of phasor domain

$V_s(t)=599\cos (377t+3.30)$ ... solution from substituting into symbolic time domain

Phase angles are exactly the same place in the third quadrant.

Simulation

$i_r$ is in phase with $V_s$ as it should. $V_s$ is leading $i_L$ through the inductor by $\frac{\pi }{2}$ as it should. Everything appears centered around 0 volts, so no constants are in play.

period check

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.7ms$

magnitude check

The magnitude of V_s is above 500 volts .. although perhaps not 600 volts.

transient response check

Everything starts out at zero. This means that only the middle to right side of the simulation windows are going to display steady state information that we can compare to calculated values.

phase check

$I_s$ brown, $i_R$ blue (with same phase angle as source voltage), and $i_L$ (orange). $i_R{V}_s$ (blue) leads $i_L$ (orange) by $\frac{\pi }{2}$. The brown current (the starting point) is the starting point at 2.09 radians. The blue current should lead the brown by 3.30-2.09=1.21 radians which is 69.3 degrees which is about 3.21 ms or 1.5 squares above. This can not be seen at the beginning where the simulation software is dealing with the initial energy embalance, but by the middle of the simulation, the blue is leading the brown by about 1.5 squares. So all is well.

Power Analysis

The constants are going to add some DC power or real power to this analysis. Without knowing what they are, we can not compute their impact. So for now we stick with phasor domain power analysis:

${𝕀}_s=120\sqrt{2}\mathrm{\angle }2.09$

${𝕍}_s=599\mathrm{\angle }3.30$

${𝕀}_s^{*}=120\sqrt{2}\mathrm{\angle }-2.09$

if :$𝕍=M_v\mathrm{\angle }{\varphi }_v$, and $𝕀=M_i\mathrm{\angle }{\varphi }_i$ then

$𝕊=𝕍{𝕀}^{*}=\frac{M_v{M}_i}{2}\mathrm{\angle }({\varphi }_v-{\varphi }_i)=\frac{120\sqrt{2}*599}{2}\mathrm{\angle }(3.30-2.09)=10,164\mathrm{\angle }1.21=3,586+9,510j$

and

$cos(1.21)=.357$

This is a bad power factor. Power factors below .9 will cause you a visit from the utility company or the transformer on the power pole might blow. The utility's apparent power is much larger than what the customer is willing to pay.

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

Intuition

Should be way to do phasor math in symbol form in mupad ... probably with matrices.