Get a low cost access to TINACloud to edit the examples or create your own circuits
Thévenin’s Theorem for AC circuits with sinusoidal sources is very similar to the theorem we have learned for DC circuits. The only difference is that we must consider impedance instead of resistance. Concisely stated, Thévenin’s Theorem for AC circuits says:
Any two terminal linear circuit can be replaced by an equivalent circuit consisting of a voltage source (VTh) and a series impedance (ZTh).
In other words, Thévenin’s Theorem allows one to replace a complicated circuit with a simple equivalent circuit containing only a voltage source and a series connected impedance. The theorem is very important from both theoretical and practical viewpoints.
It is important to note that the Thévenin equivalent circuit provides equivalence at the terminals only. Obviously, the internal structure of the original circuit and the Thévenin equivalent may be quite different. And for AC circuits, where impedance is frequency dependent, the equivalence is valid at one frequency only.
Using Thévenin’s Theorem is especially advantageous when:
· we want to concentrate on a specific portion of a circuit. The rest of the circuit can be replaced by a simple Thévenin equivalent.
· we have to study the circuit with different load values at the terminals. Using the Thévenin equivalent we can avoid having to analyze the complex original circuit each time.
We can calculate the Thévenin equivalent circuit in two steps:
1. Calculate ZTh. Set all sources to zero (replace voltage sources by short circuits and current sources by open circuits) and then find the total impedance between the two terminals.
2. Calculate VTh. Find the open circuit voltage between the terminals.
Norton’s Theorem, already presented for DC circuits, can also be used in AC circuits. Norton’s Theorem applied to AC circuits states that the network can be replaced by a current source in parallel with an impedance.
We can calculate the Norton equivalent circuit in two steps:
1. Calculate ZTh. Set all sources to zero (replace voltage sources by short circuits and current sources by open circuits) and then find the total impedance between the two terminals.
2. Calculate ITh. Find the short circuit current between the terminals.
Now let’s see some simple examples.
Example 1
Find the Thévenin equivalent of the network for the points A and B at a frequency: f = 1 kHz, vS(t) = 10 cosw×t V.
The first step is to find the open circuit voltage between points A and B:
The open circuit voltage using voltage division:
= -0.065 – j2.462 = 2.463 e-j91.5º V
Checking with TINA:
The second step is to replace the voltage source by a short circuit and to find the impedance between points A and B:
Here is the Thévenin equivalent circuit, valid only at a frequency of 1kHz. We must first, however, solve for CT’s capacitance. Using the relationship 1/wCT = 304 ohm, we find CT = 0.524 uF
Now we have the solution: RT = 301 ohm and CT = 0.524 m F:
Next, we can use TINA’s interpreter to check our calculations of the Thévenin equivalent circuit:
VM:=10;
f:=1000;
om:=2*pi*f;
Z1:=R1+j*om*L;
Z2:=R2/(1+j*om*C*R2);
VT:=VM*Z2/(Z1+Z2);
VT=[-64.0391m-2.462*j]
abs(VT)=[2.4629]
abs(VT)/sqrt(2)=[1.7415]
radtodeg(arc(VT))=[-91.49]
ZT:=Replus((R1+j*om*L),replus(R2,(1/j/om/C)));
ZT=[301.7035-303.4914*j]
Abs(ZT)=[427.9393]
radtodeg(arc(ZT))=[-45.1693]
Ct:=-1/im(ZT)/om;
Ct=[524.4134n]
import math as m
import cmath as c
#Lets simplify the print of complex
#numbers for greater transparency:
cp= lambda Z : “{:.4f}”.format(Z)
#Define replus using lambda:
Replus= lambda R1, R2 : R1*R2/(R1+R2)
VM=10
f=1000
om=2*c.pi*f
Z1=complex(R1,om*L)
Z2=R2/complex(1,om*C*R2)
VT=VM*Z2/(Z1+Z2)
print(“VT=”,cp(VT))
print(“abs(VT)= %.4f”%abs(VT))
print(“abs(VT)/sqrt(VT)= %.4f”%(abs(VT)/m.sqrt(2)))
print(“degrees(arc(VT))= %.4f”%m.degrees(c.phase(VT)))
ZT=Replus(complex(R1,om*L),Replus(R2,1/1j/om/C))
print(“ZT=”,cp(ZT))
print(“abs(ZT)= %.4f”%abs(ZT))
print(“degrees(arc(ZT))= %.4f”%m.degrees(c.phase(ZT)))
Ct=-1/ZT.imag/om
print(“Ct=”,Ct)
Note that in the listing above we used a function “replus.’ Replus solves for the parallel equivalent of two impedances; i.e., it finds the product over the sum of the two parallel impedances.
Example 2
Find the Norton equivalent of the circuit in Example 1.
f = 1 kHz, vS(t) = 10 cosw×t V.
The equivalent impedance is the same:
ZN=(0.301-j0.304) kW
Next, find the short-circuit current:
IN = (3.97-j4.16) mA
And we can check our hand calculations against TINA’s results. First the open circuit impedance:
Then the short-circuit current:
And finally the Norton equivalent:
Next, we can use TINA’s interpreter to find the Norton equivalent circuit components:
VM:=10;
f:=1000;
om:=2*pi*f;
Z1:=R1+j*om*L;
Z2:=R2/(1+j*om*C*R2);
IN:=VM/Z1;
IN=[3.9746m-4.1622m*j]
abs(IN)=[5.7552m]
abs(IN)/sqrt(2)=[4.0695m]
radtodeg(arc(IN))=[-46.3207]
ZN:=Replus((R1+j*om*L),replus(R2,(1/j/om/C)));
ZN=[301.7035-303.4914*j]
Abs(ZN)=[427.9393]
radtodeg(arc(ZN))=[-45.1693]
CN:=-1/im(ZN)/om;
CN=[524.4134n]
import math as m
import cmath as c
#Lets simplify the print of complex
#numbers for greater transparency:
cp= lambda Z : “{:.4f}”.format(Z)
#Define replus using lambda:
Replus= lambda R1, R2 : R1*R2/(R1+R2)
VM=10
f=1000
om=2*c.pi*f
Z1=complex(R1,om*L)
Z2=R2/complex(1,om*C*R2)
IN=VM/Z1
print(“IN=”,cp(IN))
print(“abs(IN)= %.4f”%abs(IN))
print(“degrees(arc(IN))= %.4f”%m.degrees(c.phase(IN)))
print(“abs(IN)/sqrt(2)= %.4f”%(abs(IN)/m.sqrt(2)))
ZN=Replus(complex(R1,om*L),Replus(R2,1/1j/om/C))
print(“ZN=”,cp(ZN))
print(“abs(ZN)= %.4f”%abs(ZN))
print(“degrees(arc(ZN))= %.4f”%m.degrees(c.phase(ZN)))
CN=-1/ZN.imag/om
print(“CN=”,CN)
Example 3
In this circuit, the load is the series-connected RL and CL. These load components are not part of the circuit whose equivalent we are seeking. Find the current in the load using the Norton equivalent of the circuit.
v1(t) = 10 cos wt V; v2(t) = 20 cos (wt+30°) V; v3(t) = 30 cos (wt+70°) V;
v4(t) = 15 cos (wt+45°) V; v5(t) = 25 cos (wt+50°) V; f = 1 kHz.
First find the open circuit equivalent impedance Zeq by hand (without the load).
Numerically
Below we see TINA’s solution. Note that we replaced all the voltage sources with short circuits before we used the meter.
Now the short-circuit current:
The calculation of the short-circuit current is quite complicated. Hint: this would be a good time to use Superposition. An approach would be to find the load current (in rectangular form) for each voltage source taken one at a time. Then sum the five partial results to get the total.
We will just use the value provided by TINA:
iN(t) = 2.77 cos (w×t-118.27°) A
Putting it all together (replacing the network with its Norton equivalent, reconnecting the load components to the output, and inserting an ammeter in the load), we have the solution for the load current that we sought:
By hand calculation, we could find the load current using current division:
Finally
I = (- 0.544 – j 1.41) A
and the time function
i(t) = 1.51 cos (w×t – 111.1°) A{The shortcircuited current by mesh current method}
om:=2000*pi;
V1:=10;
V2:=20*exp(j*pi/6);
V3:=30*exp(j*pi/18*7);
V4:=15*exp(j*pi/4);
V5:=25*exp(j*pi/18*5);
Sys J1,J2,J3,J4
J1*(R-j*2/om/C)+V1+J2*j/om/C+J3*j/om/C=0
J1*j/om/C+J2*(j*om*L-j/om/C)+V4-V2=0
J1*j/om/C+J3*(R+j*om*L-j/om/C)-J4*j*om*L+V3+V5-V4=0
-J3*j*om*L+J4*(R+j*om*L)-V3=0
end;
J3=[-1.3109E0-2.4375E0*j]
{The impedance of the ‘killed’ network}
ZLC:=j*om*L/(1-sqr(om)*L*C);
ZRL:=j*om*L*R/(R+j*om*L);
ZN:=(R+ZLC)/(1+j*om*C*(R+ZLC))+R+ZRL;
ZN=[1.3923E1-5.8456E0*j]
I:=J3*ZN/(ZN+RL-j/om/C);
I=[-5.4381E-1-1.4121E0*j]
import math as m
import cmath as c
#Lets simplify the print of complex
#numbers for greater transparency:
cp= lambda Z : “{:.4f}”.format(Z)
om=2000*c.pi
V1=10
V2=20*c.exp(1j*c.pi/6)
V3=30*c.exp(1j*c.pi/18*7)
V4=15*c.exp(1j*c.pi/4)
V5=25*c.exp(1j*c.pi/18*5)
#We have a linear system of equations
#that we want to solve for J1,J2,J3,J4:
#J1*(R-j*2/om/C)+V1+J2*j/om/C+J3*j/om/C=0
#J1*j/om/C+J2*(j*om*L-j/om/C)+V4-V2=0
#J1*j/om/C+J3*(R+j*om*L-j/om/C)-J4*j*om*L+V3+V5-V4=0
#-J3*j*om*L+J4*(R+j*om*L)-V3=0
import numpy as n
#Write up the matrix of the coefficients:
A=n.array([[complex(R,-2/om/C),1j/om/C,1j/om/C,0],
[1j/om/C,1j*om*L-1j/om/C,0,0],
[1j/om/C,0,R+1j*om*L-1j/om/C,-1j*om*L],
[0,0,-1j*om*L,R+1j*om*L]])
b=n.array([-V1,V2-V4,V4-V3-V5,V3])
J1,J2,J3,J4=n.linalg.solve(A,b)
print(“J3=”,cp(J3))
#The impedance of the ‘killed’ network
ZLC=1j*om*L/(1-om**2*L*C)
ZRL=1j*om*L*R/(R+1j*om*L)
ZN=(R+ZLC)/(1+1j*om*C*(R+ZLC))+R+ZRL
print(“ZN=”,cp(ZN))
I=J3*ZN/(ZN+RL-1j/om/C)
print(“I=”,cp(I))