jueves, 7 de junio de 2012

Problema7: Circuitos de 2º Orden

Ej7. Realizar un estudio del circuito de la figura respecto al instante t = 0:


Obtener la evolución de la corriente i(t).



La fuente de tensión del circuito se comporta cómo una función escalón, cuya definición es:


Una vez tenemos clara las condiciones del problema, vamos a abordarlo teniendo en cuenta el instante t = 0.


· t < 0.

En este instante, la fuente de tensión es nula, por lo que todos los datos que podamos obtener serán nulos, ya que en el enunciado del problema, no nos han indicado nada de que los elementos estén cargados con anterioridad (por lo que los consideramos totalmente descargados).

En modo resumen:

· ip = iL1p = 0 A
· V1 = vC1p = 0 V

Ahora, pasamos a estudiar el circuito en el instante t = 0.


· t = 0.

El circuito equivalente que tenemos es el siguiente:



Las condiciones de equivalencia:

· Condensador: Un condensador sin polaridad y de la misma capacidad pero sin ninguna carga inicial, con una fuente de tensión lineal en serie de valor y polaridad a la carga inicial del condensador.

· Inductor: Un inductor del mismo valor pero sin ninguna carga inicial, con una fuente de corriente lineal en paralelo de valor y polaridad a la carga inicial del inductor.

· NOTA: Tanto el condensador C2 cómo el inductor L2 es de la misma capacidad e inductancia respectivamente, que el condensador inicial del circuito (C1) y la inductancia inicial del circuito (L1), pero no tiene la misma polaridad y ninguna carga inicial, es por ello, que al no ser el mismo (ya que hemos empleado las condiciones de equivalencia para cuando existen elementos cargados), se ha nombrado cómo: C2 y L2.

En estas condiciones, tenemos dos comportamiento en el circuito, por un lado: El régimen transitorio (in) y por el otro: El régimen permanente (ip). La solución general será del tipo:

· i(0+) = ip + in

Así que vamos a obtener los dos por separado.


· Régimen Permanente:

Este caso, el circuito ha estado en la configuración presentada, durante un largo período de tiempo. El inductor se comportará como un corto circuito y el condensador como un circuito abierto.

El circuito equivalente es:



Por lo que, los datos que tendremos son:

· VC2p = 0 V
· ip = iL2p = V1/R1 = 12/2·103 = 6 mA


· Régimen Transitorio:

Tanto el condensador como el inductor, se irán cargando, realizamos el análisis mediante LVK:

· [Ec1] - V1 + VR1 + VL2 = 0

Donde:

· VL2 = L2Din
· VR1 = i1·R1

Si tenemos realizamos un análisis en el nodo A:

· Nodo A: i1 = in + iC2

Sustituimos:

· VC2 = VL2
· iC2 = C2DVC2 = C2DVL2= C2·L2D2 in

Ya podemos sustituir los valores obtenidos en la [Ec1]:

· L2Din + [in + C2·L2D2 in]·R1 = V1

Reordenamos y dividimos por C2·L2:

· D2 in + Din·(1/R1·C2) + (1/(L2·C2))·in = V1/R1·L2·C2

Ya tenemos nuestra ecuación diferencial de 2º Orden, cuya homogénea es:

· D2 in + Din·(1/R1·C2) + (1/(L2·C2))·in = 0

Que sabemos que debe ser del tipo:

· s2 + 2α·s + w20 = 0

Y en este caso, los parámetros α y w0 son:

· Factor de amortiguamiento ≡ α = 1/2R1·C2 = 1/(2·2·103·5·10-6) = 50 s-1
· Frecuencia natural ≡ w0 = √[ (1/L2·C2)] = √[ (1/8·10-3·5·10-6)] = 5000 rad/s

Y teniendo el valor de estos dos parámetros, podemos determinar cual es el comportamiento de nuestro circuito. En nuestro caso:

· α < w0

Al ser la frecuencia natural mayor que el factor de amortiguamiento, el comportamiento de nuestro circuito es: Respuesta Subamortiguada, cuya solución a la respuesta homogénea es:

· in(t) = e-α·t·[k1·cos(wd·t) + k2·sin(wd·t)]

Dónde:

· Frecuencia natural amortiguada ≡ wd = (w20 - α2)1/2 = (50002 - 502)1/2 ≈ 4999.749994 rad/s

Por lo tanto, la solución completa es, la compuesta por la respuesta permanente más la respuesta transitoria:

· i(t) = ip + in = 0.006 + e-50·t·[k1·cos(4999.749994·t) + k2·sin(4999.749994·t)]

Ya tenemos la estructura de cómo será la evolución de la corriente al abrirse el interruptor de nuestro circuito. Nos falta obtener los valores de los parámetros.

Para ello, vamos a aplicar las condiciones iniciales de los componentes que forman nuestro circuito.


· Inductor:

La corriente no puede cambiar bruscamente en un instante concreto. El inductor de nuestro circuito, no presenta ninguna carga de intensidad inicial:

· iL2(0-) = 0 A

Tenemos todos los datos para poder trabajar:

· i(t) = 0.006 + e-50·t·[k1·cos(4999.749994·t) + k2·sin(4999.749994·t)]
· iL2(0-) = iL2(0+) = 0 A

· NOTA: El inductor L2 está totalmente descargado.

Por lo tanto:

· 0.006 + e-50·t·[k1·cos(4999.749994·t) + k2·sin(4999.749994·t)] = 0

En t = 0:

· 0.006 + e-50·0·[k1·cos(4999.749994·0) + k2·sin(4999.749994·0)] = 0

Operamos:

· k1 = - 0.006

Ya tenemos el valor del primer parámetro, ahora vamos a por el del segundo parámetro, para ello, analizamos las condiciones iniciales del siguiente componente.


· Condensador:

La tensión no puede cambiar bruscamente en un instante concreto. El condensador de nuestro circuito, no presenta ninguna carga de tensión inicial:

· vC2(0-) = 0 V.

Como suele ser habitual, existen varias maneras de abordar esta parte, para este caso, la más fácil es mediante la propia definición de la tensión del condensador:

· [Ec2] VC2 = VL2 = L2Di(t)

Donde:

· i(t) = 0.006 + e-50·t·[- 0.006·cos(4999.749994·t) + k2·sin(4999.749994·t)]

Obtenemos la primera derivada de la corriente del circuito:

· Di(t) = - 50e-50·t·[- 0.006·cos(4999.749994·t) + k2·sin(4999.749994·t)] + e-50·t·[29.9985·sin(4999.749994·t) + 4999.749994·k2·cos(4999.749994·t)]

En el instante t = 0:

· Dv(0) = - 50·(- 0.006) + 4999.749994·k2 = 0.3 + 4999.749994·k2

Tenemos en cuenta las condiciones del condensador:

· vC2(t-) = vC2(t+) = 0 V

· NOTA: El inductor C2 está totalmente descargado.

Sustituimos en [Ec2]:

· 0 = L2·[0.3 + 4999.749994·k2]

Sustituimos el valor del inductor y despejamos el parámetro que nos interesa obtener:

· k2 = - 0.3/4999.749994 ≈ - 0.00006

Ya hemos obtenido ambos parámetros, por lo que, una expresión para la corriente del inductor (y la que nos piden obtener en el problema) es:

· i(t) = 0.006 + e-50·t·[- 0.006·cos(4999.749994·t) - 0.00006·sin(4999.749994·t)], t ≥ 0


Por lo tanto, en modo resumen, una expresión para la evolución de la corriente, dadas las condiciones que nos expone el enunciado del problema, es:



Cuya representación gráfica es:



Os dejo una simulación realizada en LTSpice:

Problema 7: Transitorios de 2º Orden
Problema 7
Problema 7
Problema 7Enlace

0 comentarios: