Condiciones Iniciales en Circuitos Eléctricos
Condiciones Iniciales en Circuitos Eléctricos
CONDICIONES INICIALES
1
2 CONDICIONES INICIALES
II Determinar las corrientes en todos los inductores y los voltajes en todos los
capacitores, es decir, el estado energético, en t = 0− a partir de:
3.2.1. DEFINICIÓN
R R 2 6
R R
C 8
C1 vC1 4
L2
L1 L2 L1
C2 vC2 iL2 v(t) 1 7
v(t)
R iL1 L3
R 3 5
L3
iL3
(a) (b)
Figura 3.1 Circuito del Ejemplo 3.1 y su correspondiente gráfico orientado
Los parámetros del circuito de la figura 3.1(a), cuyo gráfico orientado se muestra
en la Figura 3.1(b), son los siguientes:
R1 = 3 Ω C1 = 300 µ C2 = 1, 500 µ
L1 = 6 · 0 mH L2 = 1 · 5 mH L3 = 3 · 0 mH
M12 = 1 · 0 mH M13 = 2 · 0 mH (3.4)
(3.5)
6 + 4 · 5 × 10−3 D 6 + ·0 × 10−3 D 3 1
i4
6 + 3 · 0 × 10 D 6 + 6 · 0 × 10 D
−3 −3
3 i7 = 1 (3.6)
3D 3D 6D + 4000 i8 D
40 80
vL1ss (t) = 40 Cos 103t vL2ss (t) = Cos 103t vL3ss (t) = Cos 103t
3 3
100 20
vC1ss (t) = Sen(103 t + 16,26◦ ) vC2ss (t) = Sen(103 t + 16,26◦ ) (3.8)
3 3
8 20
i4p (t) = iL2ss (t) ⇒ A Sen 103 t A
3 3
4 10
i7p (t) = iL1ss (t) ⇒ A Sen 103 t A
3 3
R R
C1 vC1
L1 L2
C2 vC2 iL2
V
R iL1 L3
iL3
Figura 3.2 Circuito Para hallar respuestas de estado estacionario del circuito de
la Figura 3.1(a) excitado con fuente de valor constante
En la solución del circuito resistivo resultante debe tenerse en cuenta que la corrien-
te y por lo tanto también el voltaje a través del resistor en serie con los capacitores
son nulos y que el terminal negativo de la fuente independiente de voltaje de valor
constante y el terminal común a los inductores y el resitor es el mismo nodo.
Vale la pena hacer notar, además, que la solución de este circuito resistivo permite
obtener únicamente la suma algebraica de los voltajes en capacitores conectados
en serie y de las corrientes a través de inductores en paralelo. Para determinar los
valores individuales se requieren ecuaciones adicionales que se obtienen aplicando
los siguientes criterios:
4
Nótese, sin embargo, que la energía almacenada en ellos es diferente de cero ya que existe una
corriente (6= 0) a través de los inductores y un voltaje (6= 0) en terminales de los capacitores
Como en este caso también se cumple que la suma algebraica de las corrientes a
través de los elementos del corte es nula se puede re-escribir la siguiente ecuación
para la correspondiente superficie:
Z t XZ t
ik (t) = 0 ⇒ ik (τ )dτ = ik (τ )dτ = qk (t) = Ck vk (t) = 0 ∀ t
X X X X
k −∞ k k −∞ k k
(3.10)
(3.10) establece que la siguiente LEY GENERAL: “cuando un corte está for-
mado exclusivamente de capacitores la suma algebraica de sus cargas
es igual a cero en todo instante”. En la aplicación de (3.10) debe tenerse en
cuenta que el signo de cada término lo define el sentido de referencia de cada vk (t).
La aplicación de (3.10) al caso particular de una secuencia de n capacitores co-
nectados en serie como los mostrados en la Figura 3.3 se resume en las ecuaciones
(3.11) y (3.12).
iab(t) C1 C2 Ck
a
v1 v2 vk vn Cn
vab(t)
b
Figura 3.3 Aplicación de (3.10) al caso de capacitores en serie
n
1 t qk (t)
Z
vab (t) = vk (t) vk (t) = iab (τ )dτ = (3.11)
X
k=1 Ck −∞ Ck
De (3.11) y dado que en una conexión serie la corriente es la misma para todos los
elementos se concluye que:
Z t
qk (t) = iab (τ )dτ = C1 v1 (t) = C2 v2 (t) = · · · = Cn vn (t)(3 − 14b) (3.12)
−∞
EJEMPLO 3.2
C1 C2
1 2 5 6
vC1 vC2
vC4 C4 C6 vC6
C3 V2
V1 3 4
vC5 vC3
C5 C7 vC7
0
Obtener los voltajes (con los sentidos de referencia indicados) para el circuito de
la Figura 3.4, cuyos parámetros son los de la ecuación (3.13)
V1 = 12 V V2 = 24 V C1 = 30 µF
C2 = 30 µF C3 = 270 µF C4 = 45 µF
C5 = 90 µF C6 = 45 µF C7 = 90 µF (3.13)
S7
1 1 2 S5
5 2 6
S1 S6
4 6 4
3 3 4
8 5 7 9
0
Z t XZ t
vk (t) = 0 ⇒ vk (τ )dτ = vk (τ )dτ = φk (t) = 0 ∀ t (3.17)
X X X
k −∞ k k −∞ k
~ =
donde las expresiones para los flujos magnéticos individuales se obtienen de φ(t)
[L]~i(t), donde [L] es la matriz simétrica de inductancias propias y mutuas, φ(t) ~
es el vector de flujos magnéticos asociados e ~i(t) el vector de corrientes en el caso
de inductores mutuamente acoplados [ecuación (1.33)] y para los no acoplados de
φk (t) = Lk ik (t) [definición (1.29)]
(3.17) establece la siguiente LEY GENERAL: “Cuando un anillo está forma-
do exclusivamente de inductores (algunos de los cuales podrían estar
inductivamente acoplados) la suma algebraica de los flujos magnéticos
asociados a ellos es igual a cero en todo instante”
vab(t) L1 L2 Lk Ln
iL1 iL2 iLk iLn
b
Figura 3.6 Conjunto de n inductores conectados en paralelo
1 t φk (t)
Z
iLk (t) = vab (τ ) dτ =
Lk −∞ Lk
Z t
⇒ vab (τ ) dτ = φk (t) = Lk iLk (t) k = 1, 2, · · · , n (3.18)
−∞
L1 iL1 = L2 iL2 = · · · = Lk iLk = · · · = Ln iLn
3.2.4. EXCEPCIONES
Las condiciones de que un circuito esté excitado por fuentes independientes de va-
lor constante y tenga un elemento resistivo no redundante, por sí mismas, no no
garantizan se puedan reemplazar los inductores por un corto-circuito y los capa-
citores por circuitos abiertos para obtener las respuestas del estado estacionario.
Es decir, ellade la misma frecuencia angulars son condiciones necesarias, pero no
suficientes.
Los inductores que hagan parte de trayectorias cerradas que contengan por lo menos
una fuente independiente de voltaje de valor constante y exclusivamente inducto-
res no se pueden reemplazar por un corto circuito. Sin embargo, los inductores
de una trayectoria cerrada formada exclusivamente por inductores se
pueden reemplazar por corto-circuide la misma frecuencia angulartos.
Similarmente, los capacitores que hagan parte de cortes que contengan por lo menos
una fuente independiente de corriente de valor constante y exclusivamente capacito-
res no se pueden reemplazar por un circuito abierto (ver ejercicio 3.3). Sin embargo,
los capacitores de un corte formado exclusivamente por capacitores se
pueden reemplazar por circuitos abiertos.
Se aclara que las situaciones descritas anteriormente son muy hipotéticas o teó-
ricas por cuanto es imposible construir en la práctica bobinas que no presenten
resistencia interna lo que se modela colocándo ésta en serie con el parámetro pre-
dominante inductancia. Así mismo todos los capacitores presentan fugas que se
modelan mediante una resistencia en paralelo con el parámetro capacitancia.
0 ∀t 6= 0
(
δ(t) =
−∞ δ(t) dt = 1
R∞
Z ∞
f (t) δ(t − t0 ) dt = f (t0 ) (3.20)
−∞
1 0+ 1 0+
Z Z
iL (0+) = i[ (0−) + vL (x)dx vC (0+) = vC (0−) + iC (x)dx
L 0− C 0−
1 Z 0+
1 0+ Z
= iL (0−) + KL δ(x)dx = vC (0−) + KC δ(x)dx
L 0− C 0−
KL KC
= iL (0−) + = vC (0−) + (3.21)
L C
1 t 1 t
Z Z
iL (t) = vL (x)dx = iL (0+) + vL (x)dx
L −∞ L 0+
Z t
~iL (t) = ~iL (0+) + [L]−1 ~vL (x)dx
0+
1 t 1 t
Z Z
vC (t) = iC (x)dx = vC (0+) + iC (x)dx (3.22)
C −∞ C 0+
respuestas y(t) (la corriente o voltaje en una cualquiera de las puertas), la cual
puede expresarse de cualquiera de las siguientes formas
π
y(t) = Am Sen(wt + α) = Am Cos(wt + β) − π ≤ α, β ≤ π β = α− (3.23)
2
donde tanto α como β son valores algebraicos, es decir, pueden ser valores mayores
o menores que cero.
La acción de los interruptores se produce cuando esta señal se encuentra en cual-
quiera de los siguientes casos:
(8) y(t) tiene una valor negativo dado n (< 0) y está creciendo.
Tabla 3.2
8 16
iL1 (0−) = − A iL2 (0−) = iL3 (0−) = − A
3 3
vC1 (0−) = −20 V vC2 (0−) = −4 V (3.24)
Reemplazando (3.24) en
x
− π2 − π4 0 π − π2 π 2π 3π x = ωt + α
−π − 3π
4 4
3π
4
5π
4
3π
2
7π
4
9π
4
5π
2
11π
4
−δ
x
x = ωt + α
π
−
7 δ= 2 φ=π
3π
2 ! !
|n| |n|
2π − ArcSen π + ArcCos
Am Am !
8 δ= ! φ=
|n| |n|
−π + ArcCos
−ArcSen
Am Am
Tabla 3.3 Expresión para las fuentes independientes sinusoidales que hacen parte
tanto del circuito previo como del conmutado
R R t=0
a
C1 vC1 L2
L1 iL2
C2 vC2
v(t) L3
R iL1
iL3
Reemplazando 103 tc = 53 · 13◦ − 106 · 26◦ = −53 · 13◦ de (3.25) en las expresiones en
función del tiempo para los voltajes en los capacitores [(3.8)] y en las mostradas para
las corrientes en los inductores de la 3.1 (tercera columna) se obtiene el siguiente
estado energético en t = 0−.
vC1 (0−) = −20 V vC2 (0−) = −4 V ∴ vC (0−) = vC1 (0−) + vC2 (0−) = −24 V
8 16 16
iL1 (0−) = − A iL2 (0−) = − A iL3 (0−) = − A (3.26)
3 3 3
CAMBIO DE REFERENCIA
Reemplazando 103 t = 103 tc = −53 · 13◦ en v(t) se obtiene el valor de vx (0) y del
d
signo de la derivada v(t) se sabe si en el instante de la conmuta-
dt 103 t=103 tc =−53·13◦
ción la fuente independiente está creciendo o decreciendo. En nuestro caso se tienen
los siguientes valores:
En la Figura 3.9 se muestran las gráficas de algunas del las respuestas del circuito
previo en el estado estacionario. Nótese que mientras la escla del eje horizontal es
la misma para todas ellas la del eje vertical para las corrientes se ha afectado en
factor de 5.
100
90
80
70
60
50
40
30
20
iL2(t) = 20 3
3 Sen(10 t)
10
ic(t) = 10Sen(103t + 106.26◦)
0 103t
− π2 − π4 0 π − π2 3π π 5π 3π 7π 2π 9π 5π 11π 6π 13π 7π 15π 4π
−10 4 4 4 2 4 4 2 4 2 4 2 4
vC1(t) = 80 3
3 Sen(10 t + 1
−20
−30
−60
−70
−80
−90
−100
Figura 3.9 Respuestas de estado estacionario del circuito previo del Ejemplo 3.4
para ilustrar el cambio de referencia.
1. (a) Para cada capacitor que no haga parte de trayectorias cerrada impropias
aplicar el principio de continuidad de los voltajes. Es decir,
(b) Para cada inductor que no pertenezca a ningún corte impropio y que
no esté acoplado a otros que hagan parte de cortes impropios aplicar el
principio de continuidad en las corrientes a través de ellos. Es decir,
2. Para determinar los voltajes en los capacitores que hagan parte de trayectorias
cerradas impropias en t = 0+ se procede de la siguiente manera:
qk (0+) = qk (0−)
X X
k k
k k
φ( 0+) = φ( 0+)
X X
k k
~ k = [L] ~iLk
φk = Lk iLk = φ ∀t (3.32)
vC3
R1 C
vC1 10
2 3 iL1 4 vC iL3 L3 8
M 9
C1 L1 iL2 im(t)
vx(t) C2 vC2 iL4
L2 C4 vC3
vy(t)
R3 5 L R2
7
1 11 iL 6
v
C in(t)
vC1 C3 1 3 3 9
2 3 9 2
C1 vC3
C2 vC2 C4 vC3 2 4
vx(t) 6
vC5
6
1
C5
6
1 5
(a) (b)
Figura 3.11 Circuito reducido para obtener los voltajes en los capacitores que
hacen partes de trayectorias impropias en el Ejemplo de la sección
3.8.1 y su correspondiente gráfico orientado
C2 vC2 (0+) − C1 vC1 (0+) − C3 vC3 (0+) = C2 vC2 (0−) − C1 vC1 (0−) − C3 vC3 (0−)
C3 vC3 (0+) + C4 vC4 (0+) = C3 vC3 (0−) + C4 vC4 (0−)
C5 vC5 (0+) − C3 vC3 (0+) = C5 vC5 (0−) − C3 vC3 (0−) (3.36)
2. (a) Como im (t), L1, L2, L3 y L4 forman parte de cortes impropios, reempla-
zando todo lo demás por corto-circuitos se obtiene el circuito reducido
de la Figura 3.12(a). (Nótese que todos los demás elementos se han
reemplazado por corto-circuitos).
(b) La figura 3.12(b) muestra el gráfico orientado del circuito reducido en el
que se ha elegido la fuente de corriente como enlace y se ha respetado los
sentidos de referencia especificados para describir el estado energético
inicial en t = 0−. Nótese que iLk = ik 1 ≤ k ≤ 4.
(c) i Aplicando la primera ley de Kirchhoff ik (t) = 0 a los cortes
P
k
independientes en el instante t = 0+ se obtiene:
4 8
M 5 8
12
iL2
L1 L2
L3 iL4 3
im(t) iL3 L4
iL1
1 2 3 5 6 7 9
4
1 2 3 5 6 7 9
(a) (b)
Figura 3.12 Circuito reducido para obtener las corrientes en los inductores que
hacen parte de cortes impropios del ejemplo de la subsección 3.8.1 y
su correspondiente gráfico orientado
(L1 − M)iL1 (0+) + (−M + L2 )iL2 (0+) = (L1 − M)iL1 (0−) + (−M + L2 )iL2 (0
L3 iL3 (0+) + L4 iL4 (0+) = L3 iL3 (0−) + L4 iL4 (0−) (3.
(a) considerar primero los elementos que hacen parte de la malla exterior
eligiendo primero los que también son enlaces y después los que también
las ramas:
i Si el elemento de la malla exterior es un enlace la ecuación toma la
forma ±ik = imx , donde k identifica el enlace en el gráfico orientado
y x la malla.
ii Si el elemento exterior es una rama distinguida en el gráfico orienta-
do con el número p se expresa la corriente a través de ella en función
de las de enlace (aplicando equilibrio de corrientes al corte que le
corresponde) en cuyo caso la ecuación toma la forma
EJEMPLO 3.6
R
5
im3 2
iL 2 C 1 3
1 3 4
L
iL1 vC 6
r
L1
6 R 6 4
3
im1 4 im2
v(t) R1 2
1
5 5
d2 i6 d2 i6 d2 i6
Hallar (0+), (0+), y (0+)
dt2 dt2 dt2
Si se describe el circuito de la Figura 3.13(a) mediante las corrientes de enlace
definidas en la Figura 3.13(b) se obtiene:
R 0+
0+ [i7 (τ ) + i8 (τ )]dτ = 0
i7 (0+) = IL0 i8 (0+) = IL10 v4 (0+) = vC (0+) = VC0 (3.48)
di
L dt7(0+) VC0 C dv4 (0+)
2 dt
1 3
r IL0
IL10 di
6
L1 dt8(0+) R
4
v(0+) R1
! !
dv di6 di6 di7 di8 di6
(t) = r +R − +R +
dt dt dt dt dt dt
2
d i7 i7 + i8
!
di6 di7
0 = L 2 + −R − ∀t≥0+
dt C dt dt
d2 i8 i7 + i8
!
di6 di8 di8
0 = L1 2 + −R + + R1 (3.50)
dt C dt dt dt
Evaluando (3.50) en t = 0+
" # " #
dv di6 di6 di7 di8 di6
(0+) = r (0+) + R (0+) − (0+) + R (0+) + (0+)
dt dt dt dt dt dt
d2 i7 i7 (0+) + i8 (0+)
" #
di6 di7
0 = L 2 (0+) + −R (0+) − (0+)
dt C dt dt
d2 i8 i7 (0+) + i8 (0+)
" #
di6 di8 di8
0 = L1 2 (0+) + −R (0+) + (0+) + R1 (0+)
dt C dt dt dt
(3.51)
d2 v d2 i6 d2 i6 d2 i7 d2 i8 d2 i6
" # " #
(0+) = r (0+) + R (0+) − (0+) + R (0+) + (0+)
dt2 dt2 dt2 dt2 dt2 dt2
(3.52)
d2 i6
de la cual se puede hallar (0+)
dt
De la inspección del circuito y del gráfico orientado de la Figura 3.13 se puede
obtener:
1 0 0
i6 i6 im1 im1 i6
i3 = i6 + i8 = im2 ⇒ im2 = 1 0 1 i7 (3.53)
i5 i6 − i7 −im3 im3 −1 1 0 i8
Observaciones:
(a) El problema se simplifica cuando los parámetros del circuito tienen valores
numéricos definidos.
di
Ri + L =V ∀ t ≥ 0+ (3.55)
dt
cuya solución es de la forma:
V R
i(t) = + Ke− L t ∀ t ≥ 0+ (3.56)
R
donde K se escoge de modo que (3.56) satisfaga la condición inicial (3.54). Es decir,
V
i(0+) = +K =0 (3.57)
R
Reemplazando K obtenida de (3.57) en (3.56) se obtienen las siguientes respuestas:
V R
i(t) = 1 − e− L t
R
di R
vL (t) = L = V e− L t ∀ t ≥ 0+ (3.58)
dt
di
Nótese que aunque i(0+) = 0 vL (0+) = L (0+) 6= 0. Además, en el estado
dt
estacionario las respuestas de régimen permanente son las siguientes:
V
iss (t) = lim i(t) = vLss (t) = lim vL (t) = 0 (3.59)
t→∞ R t→∞
V − 1 t 1
i(t) = e RC vC (t) = V 1 − e− RC t ∀ t ≥ 0+ (3.60)
R
vC
Nótese que aunque vC (0+) = 0, iC (0+) = C (0+) 6= 0. Obsérvese también que
dt
en el estado estacionario las respuestas de régimen permanente son las siguientes:
iss (t) = lim i(t) = 0 vCss (t) = lim vC (t) = V (3.61)
t→∞ t→∞
t t
fx (t) = F0 1 − e− T fy (t) = G0 e− T ∀ t ≥ 0+ (3.62)
t
0 1 2 3 4
T
fx (t) t
= 1 − e− T 0 · 0 0 · 63 0 · 86 0 · 95 0 · 982
F0
fy (t) t
= e− T 1 · 0 0 · 37 0 · 14 0 · 05 0 · 018
G0
y (t)
La Tabla 3.4 muestra los valores de fxF(t)
0
y fG 0
, de la ecuación (3.62), en función
de T . Se puede notar que “la constante de tiempo es aquel valor para el
t
régimen permanente en menos del 2 %. La Figura 3.16 muestra las curvas de las
ecuaciones (3.62).
f ( Tt )
fx(t) t
F0 = 1 − e
1 − T
0.75
0.632
0.5
0.367
0.25
fy
= e − Tt
G0
0
t
0 1 2 3 4 T
Vale la pena hacer notar que los circuitos discutidos contienen solamente un ele-
mento almacenador de energía y que existe una relación directa entre su constante
de tiempo y la única raíz de la ecuación característica de la ecuación diferencial.
Por lo tanto, es de esperarse que haya varias constantes de tiempo en circuitos con
más de un elemento almacenador de energía debido a la relación entre el número
de éstos y el orden de la ecuación diferencial.
R1 = 5 Ω I = 10 A L1 = 6 mH
L2 = M12 = 3 mH L3 = 16 mH L4 = 9 · 4 mH
M34 = 4 mH C = 500 µF
R = 3 Ω v(t) = 25 Cos(500t − 36 · 87◦ ) (3.63)
1. iL1 (t) ∀ t ≥ 0+
2. iL2 (t) ∀ t ≥ 0+
3. vC (t) ∀ t ≥ 0+
4. iL3 (t) ∀ t ≥ 0+
I R1
L1 L2
iL1 iL2
Figura 3.18 Circuito previo excitado con fuente de valor constante
6 −3
" # " #" #
φL1 iL1
= 10 −3
×
φL2 −3 3 iL2
9 iL1 − 6 iL2 = 0 (a) iL1 = 4 A = iL1 (0−)
)
φL1 − φL2 = 0 ∴ ⇒ (3.64)
iL1 + iL2 = 10 A (b) iL2 = 6 A = iL2 (0−)
16 4
" # " # " #
v3 d i3
= 10 −3
v4 9·4 dt i4
Rt
v1 = v(t) v5 = v5 (0) + 2000 0 −i3 (τ )dτ v2 = 3 i3 (3.65)
v3 − v1 + v2 − v5 (3.66)
di3 t
Z
16 × 10−3 + 3i3 − {v5 (0) − 2000 i3 (τ )dτ } = v(t) ∀ t ≥ 0 (3.67)
dt 0
d2 i3 di3 dv
16 × 10−3 2
+3 + 2000i3 = (3.68)
dt dt dt
cuya solución particular se obtiene suponiendo una respuesta de la forma su-
gerida en (3.69) y reemplazándola con sus derivadas temporales en la ecuación
diferencial lineal ordinaria (edlo) (3.68)
Comparando (3.73), (3.74), (3.75) y (3.76) se puede notar que todos los casos se
refieren al mismo instante de conmutación y, por lo tanto, el estado energético inicial
en t = 0− se obtiene reemplazando en las expresiones (3.72) para las corrientes en
todos los inductores y los voltajes en todos los capacitores 500t = 500tc = −53·13◦,
es decir,
v(t) = 25 Cos(−53 · 13◦ − 36 · 87◦ ) = vx (0) = 0
500t = −53 · 13◦
d dvx
◦ = −25 × 500Sen(−53 · 13 − 36 · 87 ) = (0) > 0
◦ ◦
v(t)
dt 500t = −53 · 13 dt
(3.78)
La expresión para una función sinusoidal que en vale t = 0 vale 0 y está creciendo,
de amplitud Am = 25 V y frecuencia angular w = 500 rad/seg, suponiendo una
expresión de la forma Am = Am Sen(wt + α) es:
6 −3 φL1 (0−) 6
" # " #" # " # " #
φL1 iL1
= 10−3 ∴ = × 10−3
φL2 −3 3 iL2 φL2 (0−) 6
16 4 φL3 (0−)
" # " #" # " # " #
φL3 iL3 −64
= 10 −3
∴ = × 10−3
φL4 4 9·4 iL4 φL4 (0−) −16
(3.81)
Aplicando la Primera ley de K. ik (t) = 0 en el nodo b evaluada en t = 0+, y
P
k
la ley de continuidad de flujos magnéticos en los anillos independientes formados
exclusivamente por inductores se obtiene el siguiente sistema lineal de ecuaciones:
cuya solución para las ecuaciones (a), (b), (c) y (d) en (3.82) es la siguiente:
L4 C R
L3 vx(t)
L1 L2
ECUACIONES PRIMITIVAS
6 −3 16 4
" # " # " # " # " # " #
v1 d i1 v3 d i3
= 10 −3
= 10 −3
v2 −3 3 dt i2 v4 4 9·4 dt i4
Z t
v5 = v5 (0+) − 2000 [i1 (τ ) + i2 (τ ) + i3 (τ )]dτ v6 = 3(i1 + i2 + i3 )
0
v7 = vx (t) (3.84)
v1 − v7 + v6 − v5 + v4 = 0
v2 − v7 + v6 − v5 + v4 = 0
v3 − v7 + v6 − v5 = 0 (3.85)
!
di1 di2 t
Z
10 −3
6 −3 − vx (t) + 3(i1 + i2 + i3 ) − v5 (0+) − 2000 [i1 (τ ) + i2 (τ ) + i3 (τ )]dτ +
dt dt 0
" !#
di3 di1 di2
+10 −3
4 +9·4 + =0
dt dt dt
!
di1 di2 t
Z
10 −3
−3 +3 − vx (t) + 3(i1 + i2 + i3 ) − v5 (0+) − 2000 [i1 (τ ) + i2 (τ ) + i3 (τ )]dτ +
dt dt 0
" !#
di3 di1 di2
+10 −3
4 +9·4 + =0 ∀ t ≥ 0+
dt dt dt
" !#
di3 di1 di2
10 −3
16 +4 + − vx (t) + 3(i1 + i2 + i3 )+
dt dt dt
Z t
− v5 (0+) − 2000 [i1 (τ ) + i2 (τ ) + i3 (τ )]dτ =0
0
(3.86)
16 × 10−3 D 2 + 3D + 2000 i3 D
Las ecuaciones diferenciales lineales ordinarias (3.89) son de grado 3 por lo que la
solución completa requiere la evaluación de cada respuesta y el sus derivadas de
primer y segundo orden en el instante de referencia t = 0+.
Evaluando el CEGIDLI (3.86) en t = 0+ y expresando en forma matricial se
obtiene:
di1
(0+)
15 · 4 × 10 6 · 4 × 10 4 × 10 dt
−3 −3 −3
−3 di2
6 · 4 × 10 12 · 4 × 10 4 × 10 =
(0+)
−3 −3
4 × 10−3 4 × 10−3 16 × 10−3 dt
di3
(0+)
dt
1 1 1
1 1 1
Reemplazando (3.79) [vx (0+)], (3.80) [v5 (0+)] y (3.83) [i1 (0+), i2 (0+) e i3 (0+)] en
(3.89) y resolviendo, se obtiene:
d2 i1 d2i2
! !
dvx di1 di2 di3
10 −3
6 −3 − +3 + + − {−2000[i1 + i2 + i3 ]} +
dt dt dt dt dt dt
2 2 2
" !#
d i3 d i1 d i2
+10−3 4 2 + 9 · 4 + 2 =0
dt dt2 dt
d2 i1 d2i2
! !
dvx di1 di2 di3
10−3 −3 2 + 3 2 − +3 + + − {−2000[i1 + i2 + i3 ]} +
dt dt dt dt dt dt
2 2 2
" !#
d i3 d i1 d i2
+10−3 4 2 + 9 · 4 + 2 =0 ∀ t ≥ 0+
dt dt2 dt
d2i3 d2 i1 d2 i2
" !# !
dvx di1 di2 di3
10 −3
16 2 + 4 + 2 − +3 + + +
dt dt2 dt dt dt dt dt
− {−2000[i1 + i2 + i3 ]} = 0
(3.91)
d2 i1
(0+)
dt2
15 · 4 × 10−3 6 · 4 × 10−3 4 × 10−3 1 "
d2 i2 di1
6 · 4 × 10 12 · 4 × 10 4 × 10 = −3 1 (0+)
(0+)
−3 −3 −3
dt2 dt
4 × 10−3 4 × 10−3 16 × 10−3
1
d2 i3
(0+)
dt2
1 1
#
di2 di3 dvx
+ (0+) + (0+) − 2000 1 [i1 (0+) + i2 (0+) + i3 (0+)] + 1 (0+)
dt dt dt
1 1
(3.92)
d2 i1
(0+) = 3 · 6873 × 105 a/seg 2
dt2
d2 i1
(0+) = 5 · 5313 × 105 a/seg 2
dt2
d2 i1
2
(0+) = 4 · 6091 × 105 a/seg 2 (3.93)
dt
Cuando se ensayan soluciones de la forma irp (t) = A Sen 500t + B Cos 500t, r =
1, 2, 3, el primer miembro de las ecuaciones diferenciales que se obtendrían de
(3.89) son de la misma forma para todas las variables. Es decir,
Los segundos miembros de las ecuaciones diferenciales para cada una de las res-
puestas tomarían las siguientes formas:
dvx d2 vx
vx (t) = 25 Sen 500t = 53 × 102 Cos 500t = 54 × 104 Sen 500t
dt dt2
−7 · 2 × 10−5 × 54 × 104 Sen 500t = −450 Sen 500t ⇒ i1p
−1 · 08 × 10−4 × 54 × 104 Sen 500t = −675 Sen 500t ⇒ i2p
−21 · 0 × 10−5 × 54 × 104 Sen 500t = −1312 · 5 Sen 500t ⇒ i3p
(3.95)
−450
A= = 2 · 2222 B = 0 i1p (t) = iL1p (t)
−202 · 5
−675
A= = 3 · 3333 B = 0 i2p (t) = iL2p (t)
−202 · 5
−562 · 5
A= = 2 · 7778 B = 0 i3p (t) = iL3p (t) (3.96)
−202 · 5
3
irh (t) = Ckr eDk t = C1r + eαt (C2r Sen βt + C3r Cos βt)
X
k=1
d
irh (t) = αeαt (C2r Sen βt + C3r Cos βt) + eαt β (C2r Cos βt − C3r Sen βt)
dt
d2
irh (t) = α2 eαt (C2r Sen βt + C3r Cos βt) + αeαt β (C2r Cos βt − C3r Sen βt) +
dt2
βαeαt (C2r Cos βt − C3r Sen βt) + eαt β 2 (−C2r Sen βt − C3r Cos βt) r = 1, 2, 3
(3.98)
Sumando a cada una de las ecuaciones (3.98) las correspondientes soluciones par-
ticulares (3.95), evaluando en t = 0+ tanto las respuestas como sus derivadas
temporales de primer y segundo orden, se obtiene:
1 0 1
C1r dir
dirp
(0+) (0+)
0 463 · 51 −187 · 5
C
2r =
−
r = 1, 2, 3
dt dt
5 5
0 −1 · 7382 × 10 −1 · 7969 × 10 C3r
d2 ir d2 irp
(0+) (0+)
dt dt
(3.100)
Reemplazando las condiciones iniciales para las respuestas totales (3.83), (3.90) y
(3.93) para construir el primer término del segundo miembro en (3.100) y, además,
los valores (3.96) de las soluciones particulares en (3.94) para calcular el segun-
do término del segundo miembro en (3.100) y expresando en forma matricial, se
obtienen los siguientes sistemas lineales de ecuaciones:
1 0 1 0 · 2667 0
C11
0 463 · 51 −187 · 5 C21 = 749 · 99 − 1111 · 1111
Es decir,
ir (t) = C1r + eαt (C2r Sen βt + C3r Cos βt) + irp (t)
(3.103)
i1 (t) = 1 · 1999 + e (−1 · 1567 Sen βt + 0 · 9332 Cos βt) + 2 · 2222 Sen ωt α = −187 · 5
αt
i2 (t) = 1 · 8001 + eαt (−1 · 7348 Sen βt − 1 · 4001 Cos βt) + 3 · 3333 Sen ωt β = 463 · 51
i3 (t) = −3 · 0002 + eαt (−1 · 4458 Sen βt − 1 · 1667 Cos βt) + 2 · 7778 Sen ωt ω = 500
EJERCICIOS
3.1 Demostrar que:
L
(a) y RC tienen dimensiones de tiempo
R
(b) ip (t) = ASenwt + BCoswt = C Sen(wt + α) = CCos(wt + β), donde
√ B
π
C = A2 + B 2 , α = ArcT an y β =α−
A 2
P 3.2 Cada uno de los circuitos mostrados en la Figura P.3.1 alcanzan el estado
estacionario con el interruptor S abierto, el cual se cierra en cierto instante
que se toma como referencia t = 0. Hallar i(t), vC (t) y vL (t) para cada caso.
Suponer que un instante antes de la conmutación vC (0−) = 0.
P 3.3 El estado energético inicial del circuito mostrado en la Figura P.3.2 es nulo.
El interruptor S se abre en un instante que se toma como referencia (t = 0+).
Hallar vC2 (t) y su valor en estado estacionario. (Nótese de los resultados que,
aunque el circuito contiene un elemento resistivo no redundante y la fuente
de corriente I es de valor constante, C2 no se comporta como un circuito
abierto en estado estacionario).
P 3.4 Después de que el circuito de la Figura P.3.3 ha alcanzado el estado estacio-
nario con el interruptor S abierto, éste se cierra instantáneamente. Tomando
Figura P.3.2
Figura P.3.3
Figura P.3.4
Figura P.3.5
di
(b) i(t) ∀t ≥ 0+ y evaluar (0+) y comparar con el resultado obtenido en
dt
(a).
R1 3
vC2 1
2
3
1 3 2 4
C2 iL
6 1
C1 5
2 7
r vC1 R 0
i(t)
0
Figura P.3.6
Determinar
d2 v10 d2 v20 d2 v30
(0+) (0+) (0+) (3.107)
dt2 dt2 dt2
P 3.8 Para el circuito y gráfico orientado mostrados en la Figura P.3.7 determinar:
v20 v2 iN +2 im2
~vn = ~vR = ~iL = ~im = (3.108)
.. .. .. ..
. . . .
vN 0 vN iB imL
vC
1 2 3 2 2 3
1 3
R C
im1 im2
v(t) L1 L2 1 4 5
iL1 M
0 iL2
0
Figura P.3.7
P 3.9 Determinar las condiciones iniciales necesarias para resolver las ecuaciones
diferenciales obtenidas en el Ejercicio P 2.9 suponiendo estado energético
inicial nulo y resolverlas para las siguientes excitaciones correspondientes:
FIGURA EXCITACIÓN
√ √
P.2.14(a) v1 (t) = √5 Cos 2t v2 (t) = 5 Sen 2t
P.2.14(b) v1 (t) = 5 Cos 2t √
P.2.14(c) √ v2 (t) = 5 Sen 2t
P.2.14(d) y P.2.14(e) v(t) = 5 Sen 2t
P.2.14(f) y P.2.14(g) i(t) = 0 · 5 Sen(104 t − 36 · 87◦ )
P 3.10 Sean ra (t), rb (t), rc (t), rd (t), re (t), rf (t) y rg (t) las soluciones obtenidas en el
ejercicio P 3.9. Establecer las siguientes comparaciones y verificar que:
Figura P.3.8
Figura P.3.9
Figura P.3.10
C1 = 1 F C2 = 2 F C3 = 1 F
L1 = 1 H L2 = 2 H M = 1H
R = 1000 Ω v(t) = 100 Sen 4t V i(t) = 2 Cos 3t A