Electronics/RCL time domain

From testwiki
Revision as of 05:53, 22 July 2023 by imported>Minorax
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
Figure 1: RCL circuit
Figure 1: RCL circuit

Template:- When the switch is closed, a voltage step is applied to the RCL circuit. Take the time the switch was closed to be 0s such that the voltage before the switch was closed was 0 volts and the voltage after the switch was closed is a voltage V. This is a step function given by Vu(t) where V is the magnitude of the step and u(t)=1 for t0 and zero otherwise.

To analyse the circuit response using transient analysis, a differential equation which describes the system is formulated. The voltage around the loop is given by:

Vu(t)=vc(t)+di(t)dtL+Ri(t) (1)

where vc(t) is the voltage across the capacitor, di(t)dtL is the voltage across the inductor and Ri(t) the voltage across the resistor.

Substituting i(t)=dvc(t)dt into equation 1:

Vu(t)=vc(t)+d2vc(t)dt2LC+Rdvc(t)dtC

d2vc(t)dt2+RLdvc(t)dt+1LCvc(t)=Vu(t)LC (2)

The voltage vc(t) has two components, a natural response vn(t) and a forced response vf(t) such that:

vc(t)=vf(t)+vn(t) (3)

substituting equation 3 into equation 2.

[d2vn(t)dt2+RLdvn(t)dt+1LCvn(t)]+[d2vf(t)dt2+RLdvf(t)dt+1LCvf(t)]=0+Vu(t)LC

when t>0s then u(t)=1:

[d2vn(t)dt2+RLdvn(t)dt+1LCvn(t)]=0 (4)

[d2vf(t)dt2+RLdvf(t)dt+1LCvf(t)]=VLC (5)

The natural response and forced solution are solved separately.

Solve for vf(t):

Since VLC is a polynomial of degree 0, the solution vf(t) must be a constant such that:

vf(t)=K

dvf(t)dt=0

d2vf(t)dt=0

Substituting into equation 5:

1LCK=VLC

K=V

vf=V (6)

Solve for vn(t):

Let:

RL=2α

1LC=ωn2

vn(t)=Aest

Substituting into equation 4 gives:

d2Aestdt2+2αdAestdt+ωn2Aest=0

s2Aest+2αAest+ω22Aest=0

s2+2αs+ωn2=0

s=2α±4α24ωn22=α±α2ωn2

Therefore vn(t) has two solutions Aes1t and Aes2t

where s1 and s2 are given by:

s1=α+α2ωn2

s2=αα2ωn2

The general solution is then given by:

vn(t)=A1es1t+A2es2t

Depending on the values of the Resistor, inductor or capacitor the solution has three posibilies.

1. If α>ωn the system is said to be overdamped

2. If α=ωn the system is said to be critically damped

3. If α<ωn the system is said to be underdamped

Example

Given the general solution

R L C V
0.5H 1kΩ 100nF 1V

α=R2L=1000

ωn=1LC4472

s1=10004359j

s2=1000+4359j

vn(t)=A1e(10004359j)t+A2e(1000+4359j)t

Thus by Euler's formula (ejϕ=cosϕ+jsinϕ):

vn(t)=e1000[A1(cos(4359t)+jsin(4359t))+A2(cos(4359t)+jsin(4359t))]

vn(t)=e1000t[(A1+A2)cos(4359t)+j(A1+A2)sin(4359t)]

Let B1=A1+A2 and B2=j(A1+A2)

vn(t)=e1000t[B1cos(4359t)+B2sin(4359t)]

Solve for B1 and B2:

From equation \ref{eq:vf}, vf=1 for a unit step of magnitude 1V. Therefore substitution of vf and vn(t) into equation \ref{eq:nonhomogeneous} gives:

vc(t)=1+e1000t[B1cos(4359t)+B2sin(4359t)]

for t=0 the voltage across the capacitor is zero, vc(t)=0

0=1+B1cos(0)+B2sin(0)

B1=1 (7)

for t=0, the current in the inductor must be zero, i(0)=0

i(t)=dvc(t)dtC

i(0)=100109[e1000t(4359B1sin(4359t)+4359B2cos(4359t))1000e1000t(B1cos(4359t)+B2sin(4359t))]

0=100109[4359B21000B1]

substituting B1 from equation \ref{eq:B1} gives

B20.229

For t>0, vc(t) is given by:

vc(t)=1e1000t[cos(4359t)+0.229sin(4359t)]

vout is given by:

vout=Vinvc(t)

vout=Vu(t)vc(t)

For t>0, vout is given by:

vout=e1000t[cos(4359t)+0.229sin(4359t)]

Template:BookCat