Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Ejercicio Resuelto 7. Análisis Del Estado Transitorio de Circuitos RC, RL y RLC

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 10

Problema 15

En la figura 15 se muestra un circuito eléctrico con dos inductores, para el cual la magnitud de la
fuente de voltaje, vF , es:
5, t<0
vF = 
10 cos (30 t), t ≥ 0
a) Demuestre que la ecuación diferencial que determina la corriente del inductor superior para
t > 0 es:
d 2 i2 d i2
d t2
+ 70 dt
+ 600 i2 = -600 sin (30 t)

b) Verifique que la respuesta completa del voltaje del resistor de la derecha para t > 0 es:
v(t) = 0.6 sin (30 t) + 4.2 cos (30 t) - 2.4 e-10 t - 1.8 e-60 t
c) Obtenga con Mathematica la gráfica del voltaje obtenido en el inciso anterior y compruebe con
Proteus ISIS dicha gráfica.

Figura 15 Circuito de carga de las celdas fotovoltaicas de un satélite espacial.

a) Demostración de la ecuación diferencial que determina la


corriente del inductor superior para t > 0
Se tomará como base la ley de voltajes de Kirchhoff y las corrientes indicadas en el diagrama que
se muestra a continuación:

Las ecuaciones de malla serán:


R1 iF + vL1 = vF
2 Serie 8 Problema 15.nb

vL2 + R2 i2 = vL1
la intensidad de corriente que circula por el inductor L1 quedará:
i F = i1 + i2
Las expresiones para los voltajes en los inductores serán:
vL1 = L1 ddit1
d i2
vL2 = L2 dt
Luego de sustituir estas útlimas expresiones en las ecuaciones de malla se obtendrá:
R1 (i1 + i2 ) + L1 ddit1 = vF
d i2 d i1
L2 dt
+ R2 i2 = L1 dt

Se sumará miembro a miembro las dos ecuaciones anteriores, con lo que se tendrá:
R1 (i1 + i2 ) + L1 ddit1 + L2 ddit2 + R2 i2 = vF + L1 ddit1
de donde:
d i2
R1 i1 + R1 i2 + L2 dt
+ R2 i2 = vF
de esta última expresión se despejará i1 :
R1 i1 = vF - R1 i2 - L2 ddit2 - R2 i2
1 L2 d i2 R2
i1 = R1
vF - i2 - R1 d t
- R1
i2
1
i1 = R1
vF - (1 + RR21 ) i2 - L2
R1
d i2
dt
1
i1 = R1
vF - ( R1R+1R2 ) i2 - L2
R1
d i2
dt
y se le sustituirá en la ecuación de la malla derecha, con lo que se obtendrá:
L2 ddit2 + R2 i2 = L1 ddt ( R11 vF - ( R1R+1R2 ) i2 - RL21 ddit2 )
d i2 L1 d vF L1 (R1 + R2 ) d i2 L 1 L 2 d 2 i2
L2 dt
+ R 2 i2 = R1 d t
- R1 dt
- R1 d t2
luego de reescribir la ecuación de manera que todas las expresiones de i2 queden en el miembro
izquierdo, la expresión quedará:
L 1 L 2 d 2 i2
R1 d t2
+ L2 ddit2 + L1 (RR1 1+ R2 ) ddit2 + R2 i2 = RL11 ddvtF
L 1 L 2 d 2 i2 L2 R1 + L1 (R1 + R2 ) d i2 L1 d vF
R1 d t2
+ R1 dt
+ R 2 i2 =
R1 d t
Finalmente, se divide toda la expresión por L1R1L2 , con lo que la ecuación diferencial será:
d 2 i2
d t2
+ R1 [L2 R1R+1 LL11 L(R2 1 + R2 )] ddit2 + RL11 RL22 i2 = RR1 1L1L1L2 ddvtF
d 2 i2
d t2
+ L2 R1 +LL11 L(R2 1 + R2 ) ddit2 + RL11 RL22 i2 = L12 ddvtF
Ahora, se sustituirán los valores de todos los parámetros eléctricos:
d 2 i2 + (0.5) (10 + 15) d i2
d t2
+ (0.5) (10)(0.5) (0.5) dt
(10) (15)
+ (0.5) i = 1 d vF
(0.5) 2 0.5 d t
d 2 i2 d i2 d vF
d t2
+ 70 dt
+ 600 i2 = 2 dt
Dado que la fuente de voltaje en el intervalo t ≥ 0 será:
vF = 10 cos (30 t)
d vF
dt
= -300 sin (30 t)
y después de sustituir esta última expresión, la ecuación diferencial quedará:
d 2 i2
d t2
+ 70 ddit2 + 600 i2 = 2 [-300 sin (30 t)]
Serie 8 Problema 15.nb 3

d 2 i2 d i2
d t2
+ 70 dt
+ 600 i2 = -600 sin (30 t)
que es lo que había que demostrar.

La ecuación diferencial que determina la intensidad de corriente del inductor


superior para t > 0 es:
d 2 i2
d t2
+ 70 ddit2 + 600 i2 = -600 sin (30 t).

Obtención de la respuesta completa del voltaje del resistor de la


derecha para t > 0
Para obtener la expresión del voltaje del resistor de la derecha, R2 , será necesario determinar la
intensidad de corriente que circula por él, que es i2 . Por tanto, se resolverá la ecuación diferencial
obtenida en el inciso anterior.
Primero, se procederá a determinar las condiciones iniciales del circuito. Dado que para t < 0 el
valor de la fuente de voltaje es vF = 5 V, y se pueden considerar a los inductores como alambres, el
voltaje y la intensidad de corriente en los inductores serán:
i1 (0- ) = 0.5 A
i2 (0- ) = 0
Para t = 0+ el valor de la fuente de voltaje cambia a vF (0+ ) = 10 V, aunque las intensidades de
corriente en los inductores son los mismos pues no pueden cambiar abruptamente:
i1 (0+ ) = 0.5 A
i2 (0+ ) = 0
En el diagrama que se muestra abajo se podrá verificar el estado de las intensidades de corriente
en t = 0+ :

Por lo que la intensidad de corriente en el resistor R1 será:


iR1 (0+ ) = i1 (0+ ) - i2 (0+ )
iR1 (0+ ) = 0.5 A
y por tanto la caída de voltaje en dicho resistor será:
vR1 (0+ ) = 10 (0.5)
vR1 (0+ ) = 5 V
lo que implicará que el voltaje en el inductor L2 sea:
vL2 (0+ ) = vF (0+ ) - vR1 (0+ ) - R2 iL2 (0+ )
4 Serie 8 Problema 15.nb

vL2 (0+ ) = 10 - 5 - 0
vL2 (0+ ) = 5 V
Como:
d i2 (t)
dt
= vL2L(t)
d i2 (0+ ) vL2 (0+ )
dt
= L
d i2 (0+ ) 5
dt
= 0.5
d i2 (0+ )
dt
= 10 As

Posteriormente, se determinará la solución de la ecuación diferencial que modela al circuito,


mediante la obtención de su solución homogénea y, posteriormente la solución particular.
Para obtener la solución homogénea se requerirá la ecuación característica correspondiente:
λ2 + 70 λ + 600 = 0
cuyas raíces serán:
λ1,2 = -35 ± (35)2 - 600
λ1,2 = -35 ± 1225 - 600
λ1,2 = -35 ± 625
λ1,2 = -35 ± 25
λ1 = -10
λ2 = -60
Por consiguiente, la solución homogénea de la ecuación diferencial será:
i2,homog = C1 e-10 t + C2 e-60 t
Dado que la excitación es una sinusoidal, a solución particular será:
i2,partic = C3 sin (30 t) + C4 cos (30 t)
Luego de sustituir esta expresión en la ecuación diferencial, se obtendrá:
d2
d t2
[C3 sin (30 t) + C4 cos (30 t)] + 70 ddt [C3 sin (30 t) + C4 cos (30 t)] +
+ 600 [C3 sin (30 t) + C4 cos (30 t)] = -600 sin (30 t)
d
dt
[30 C3 cos (30 t) - 30 C4 sin (30 t)] + 70 [30 C3 cos (30 t) - 30 C4 sin (30 t)] +
+ 600 C3 sin (30 t) + 600 C4 cos (30 t) = -600 sin (30 t)
-900 C3 sin (30 t) - 900 C4 cos (30 t) + 2100 C3 cos (30 t) - 2100 C4 sin (30 t) +
+ 600 C3 sin (30 t) + 600 C4 cos (30 t) = -600 sin (30 t)
(-900 C3 - 2100 C4 + 600 C3 ) sin (30 t) + (-900 C4 + 2100 C3 + 600 C4 ) cos (30 t) = -600 sin (30 t)
(-300 C3 - 2100 C4 ) sin (30 t) + (-300 C4 + 2100 C3 ) cos (30 t) = -600 sin (30 t)
Por tanto, se podrán establecer las siguientes ecuaciones:
-300 C3 - 2100 C4 = -600
-300 C4 + 2100 C3 = 0
Se despejará C4 de la segunda expresión:
300 C4 = 2100 C3
C4 = 7 C3
y se sustituirá esta última en la primera ecuación:
-300 C3 - 2100 (7 C3 ) = -600
Serie 8 Problema 15.nb 5

-15000 C3 = -600
-600
C3 = -15 000
C3 = 0.04
y por último se sustituirá el valor de C3 en la expresión de C4 :
C4 = 0.28
Entonces, la solución particular buscada será:
i2,partic = 0.04 sin (30 t) + 0.28 cos (30 t)
Por consiguiente, la solución completa de la intensidad de corriente i2 será:
i2 (t) = C1 e-10 t + C2 e-60 t + 0.04 sin (30 t) + 0.28 cos (30 t)
Para determinar los valores de C1 y C2 será necesario sustituir las condiciones iniciales del ciruito
en la expresión anterior y en la de su derivada con respecto al tiempo. Después de sustituir
i2 (0) = 0 para t = 0 en la expresión de la solución completa, se obtendrá:
0 = C1 e-10 (0) + C2 e-60 (0) + 0.04 sin (30 (0)) + 0.28 cos (30 (0))
0 = C1 + C2 + 0.28
C1 = -C2 - 0.28
La derivada de la intensidad de corriente i2 (t) será:
d i2 (t)
dt
= -10 C1 e-10 t - 60 C2 e-60 t + 1.2 cos (30 t) - 8.4 sin (30 t)
d i2 (0- )
luego de sustituir las condiciones iniciales del circuito dt
= 0 para t = 0:
10 = -10 C1 e-10 (0) - 60 C2 e-60 (0) + 1.2 cos (30 (0)) - 8.4 sin (30 (0))
10 = -10 C1 - 60 C2 + 1.2
ahora, se sustituirá en esta última expresión la ecuación de C1 :
10 = -10 (-C2 - 0.28) - 60 C2 + 1.2
10 = 10 C2 + 2.8 - 60 C2 + 1.2
50 C2 = 4 - 10
C2 = -6
50
C2 = -0.12
por lo cual:
C1 = 0.12 - 0.28
C1 = -0.16
De donde la expresión que representa la intensidad de corriente en el inductor superior será:
i2 (t) = -0.16 e-10 t - 0.12 e-60 t + 0.04 sin (30 t) + 0.28 cos (30 t)
Finalmente, la expresión del voltaje, v, en el resistor R2 para t > 0 será:
v(t) = R2 i2 (t)
v(t) = (15) [-0.16 e-10 t - 0.12 e-60 t + 0.04 sin (30 t) + 0.28 cos (30 t)]
v(t) = -2.4 e-10 t - 1.8 e-60 t + 0.6 sin (30 t) + 4.2 cos (30 t)

La respuesta completa del voltaje en el resistor de la derecha para t > 0 es:


v(t) = -2.4 e-10 t - 1.8 e-60 t + 0.6 sin (30 t) + 4.2 cos (30 t).
6 Serie 8 Problema 15.nb

Obtención de la gráfica del voltaje en el resistor R2 para t > 0


Verificación de la solución completa del voltaje en el resistor de la derecha con el empleo de la
función DSolve de Mathematica
Antes de obtener la gráfica solicitada, se realizó la verificación de la solución completa del voltaje
en el resistor de la derecha con el empleo de la función DSolve de Mathematica, con lo que
resultó que la expresión de la intensidad de corriente en el inductor L2 así como el voltaje en el
resistor R2 fueron exactamente iguales.
A continuación se mostrará la gráfica obtenida con la expresión analítica encontrada para el
voltaje solicitado:

0.1 0.2 0.3 0.4 0.5

-2

-4

Comprabación con Proteus ISIS del funcionamiento del circuito


Para comprobar que los resultados obtenidos son los esperados, se alambró al circuito con
Proteus ISIS, cuyo diagrama se podrá observar en seguida:

Los interruptores VSWITCH abren o cierran dependiendo de si el voltaje de control aplicado es


cero o tiene un valor mayor que cero. Para este circuito, al interruptor S1 se le aplica un voltaje de
5 V para t < 0 y de cero para t > 0, y al interruptor S2 se le aplica cero para t < 0 y 5 V para t > 0, de
manera que el voltaje de la fuente corresponda al establecido en el enunciado de este problema.
Serie 8 Problema 15.nb 7

La gráfica del voltaje en el resistor R2 que se obtuvo con una gráfica analógica de Proteus ISIS
resultó ser exactamente igual a la encontrada con Mathematica, y se muestra a continuación:

Resolución del problema con Mathematica

a) Demostración de la ecuación diferencial que determina la


corriente del inductor superior para t > 0

ec1 = R1 iF[t] + vL1[t] == vF[t]


ec2 = vL2[t] + R2 i2[t] == vL1[t]
iF[t] = i1[t] + i2[t]
vL1[t] = L1 D[i1[t], t]
deriva
vL2[t] = L2 D[i2[t], t]
deriva
ec3 = ec1
ec4 = ec2
ec5 = ec3[[1]] + ec4[[1]] == ec3[[2]] + ec4[[2]]
ec6 = ec5[[1]] - L1 i1 '[t] == ec5[[2]] - L1 i1 '[t] // Expand
expande factores
ec7 = ec6[[1]][[1]] ⩵ ec6[[2]] - ec6[[1]] + ec6[[1]][[1]]
ec7[[2]] i1[t]
i1taux =
ec7[[1]]
di1taux = D[i1taux, t]
deriva
ec8 = ec4 /. i1 '[t] → di1taux
ec9 = ec8 // Expand
expande factores
L1 vF '[t] L1 vF '[t]
ec10 = - ec9[[2]] + + ec9[[1]] ==
R1 R1
8 Serie 8 Problema 15.nb

L1 L2 L2 R1 + L1 R1 + R2 L1
ec11 = i2 ''[t] + i2 '[t] + R2 i2[t] == vF '[t]
R1 R1 R1
R1 L2 R1 + L1 R1 + R2 R1 R2 R1 L1
ec12 = i2 ''[t] + i2 '[t] + i2[t] == vF '[t]
R1 L1 L2 L1 L2 R1 L1 L2
L2 R1 + L1 R1 + R2 R1 R2 1
ecDifaux = i2 ''[t] + i2 '[t] + i2[t] == vF '[t]
L1 L2 L1 L2 L2
R1 = 10
R2 = 15
1
L1 =
2
1
L2 =
2
vF[t] = 10 Cos[30 t]
coseno
dvFt = D[vF[t], t]
deriva
ecDif = ecDifaux /. vF '[t] → dvFt

b) Obtención de la respuesta completa del voltaje del resistor de la


derecha para t > 0

t0 = 0
i10m = 0.5
i20m = 0
vF0M = 10
i10M = i10m
i20M = i20m
iR10M = i10M
vR10M = R1 iR10M
vL20M = vF0M - vR10M
di2t vL2t
ec13 = ⩵
dt L2
vL20M
di20M =
L2

L2 R1 + L1 R1 + R2 R1 R2
ecCaract = λ2 + λ+ ⩵0
L1 L2 L1 L2
L2 R1 + L1 R1 + R2
bEntre2 =
2 L1 L2
bEntre2cuadrado = bEntre22
R1 R2
radicando = bEntre2cuadrado -
L1 L2
λ12 = - bEntre2 ± radicando
λ1 = - bEntre2 + radicando
λ2 = - bEntre2 - radicando
Serie 8 Problema 15.nb 9

i2homog = C1 Exp[λ1 t] + C2 Exp[λ2 t]


exponencial exponencial
i2partic = C3 Sin[30 t] + C4 Cos[30 t]
seno coseno
d2 L2 R1 + L1 R1 + R2 d R1 R2 1
ec14 = i2partic + i2partic + i2partic == dvFt
2 L1 L2 dt L1 L2 L2
dt
L2 R1 + L1 R1 + R2 R1 R2 1
ec15 = D[i2partic, t, t] + D[i2partic, t] + i2partic == dvFt
deriva L1 L2 deriva L1 L2 L2

ec16 = ec15 // Expand


expande factores
ec16[[1]][[3]] + ec16[[1]][[4]] ec16[[2]]
ec17 = == // TrigReduce
Sin[30 t] Sin[30 t] reduce funciones trigonométricas

ec16[[1]][[1]] + ec16[[1]][[2]]
ec18 = ⩵ 0 // TrigReduce
Cos[30 t] reduce funciones trigonométricas

ec19 = ec17 // Expand


expande factores
ec20 = ec18 // Expand
expande factores
ec21 = - ec20[[1]][[2]] ⩵ ec20[[1]][[1]]
ec21[[2]] C4
C4aux =
ec21[[1]]
ec22 = ec19 /. C4 → C4aux // Expand
expande factores
ec22[[2]] C3
C3Sol1 = // N
ec22[[1]] valor numérico

C4Sol1 = C4aux /. C3 → C3Sol1


i2partic1 = i2partic /. {C3 → C3Sol1, C4 → C4Sol1}

i2taux = i2homog + i2partic1


ec23 = i20M ⩵ i2taux /. t → t0
ec24 = ec23[[2]][[2]] ⩵ - ec23[[2]] + ec23[[2]][[2]]
ec24[[2]] C1
C1aux =
ec24[[1]]
di2taux = D[i2taux, t]
deriva
ec25 = di20M ⩵ di2taux /. t → t0
ec26 = ec25 /. C1 → C1aux
ec27 = ec26 // Expand
expande factores
ec28 = - ec27[[2]][[2]] == ec27[[2]][[1]] - ec27[[1]]
ec28[[2]] C2
ec29 = C2 ⩵
ec28[[1]]
C2Sol1 = ec29[[2]]
10 Serie 8 Problema 15.nb

[[ ]]
C1Sol1 = C1aux /. C2 → C2Sol1
i2t1 = i2taux /. {C1 → C1Sol1, C2 → C2Sol1}
ec30 = vt1 == R2 i2t1
ec31 = ec30 // Expand
expande factores
vt1 = ec31[[2]]

Verficación de los valores de C1 , C2 , C3 y C4 con el empleo de la función Solve de Mathematica

resp1 = Solve[{ec19, ec20, ec23, ec25}]


resuelve
C1Sol2 = C1 /. resp1[[1]]
C2Sol2 = C2 /. resp1[[1]]
C3Sol2 = C3 /. resp1[[1]]
C4Sol2 = C4 /. resp1[[1]]

c) Obtención de la gráfica del voltaje en el resistor R2 para t > 0

Verificación de la solución completa del voltaje en el resistor de la derecha con el empleo de la


función DSolve de Mathematica

resp2 = DSolve[{ecDif, i2[0] ⩵ 0, i2 '[0] ⩵ 10}, i2[t], t] // Expand


resolvedor diferencial expande factores
i2t = i2[t] /. resp2[[1]] // N
valor numérico
v = R2 i2t // Expand
expande factores
Plot[v, {t, 0, 0.5}]
representación gráfica

También podría gustarte