0% encontró este documento útil (0 votos)
5 vistas6 páginas

Resolución de Ecuación RLC en Circuito

Este documento presenta la resolución de la ecuación diferencial de un circuito RLC en serie. Se calcula la corriente i(t) para t>0 y t<0. Para t<0, la corriente es i(t) = e-9t(cos(4.36t) - 4.82sen(4.36t)). Para t>0, cuando el interruptor se abre, la corriente disminuye exponencialmente. El autor nota algunas diferencias entre las gráficas de Matlab y LTSpice, posiblemente debido a errores en la simul

Cargado por

Luis Santacruz
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
5 vistas6 páginas

Resolución de Ecuación RLC en Circuito

Este documento presenta la resolución de la ecuación diferencial de un circuito RLC en serie. Se calcula la corriente i(t) para t>0 y t<0. Para t<0, la corriente es i(t) = e-9t(cos(4.36t) - 4.82sen(4.36t)). Para t>0, cuando el interruptor se abre, la corriente disminuye exponencialmente. El autor nota algunas diferencias entre las gráficas de Matlab y LTSpice, posiblemente debido a errores en la simul

Cargado por

Luis Santacruz
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

IMT

Tarea

Resolución de la Ecuación diferencial de un circuito RLC en serie


o en paralelo.

Luis Omar Herrera Santacruz

A01700973
Dr. Alejandro Rafael Vega Salinas

Santiago de Querétaro, Qro.

19/04/2022
CIRCUITO POR RESOLVER:
Para el circuito se busca calcular i(t) para t>0

La fuente tiene un valor de 10 V y el interruptor se abre al momento de que t=0, primero


sacamos la ecuación diferencial y tomaremos el momento de i
𝑉𝑓 = 𝑉𝑟 + 𝑉𝑐 + 𝑉𝑙
Al poner los valores nos queda la siguiente ecuación:
𝑡
𝑑𝑖 𝐿 1
𝑉𝑓 = 𝐿 ∗ + 𝑉𝑐(0+ ) + ∗ ∫ 𝑖 𝑑𝑡 + 𝑅𝑖
𝑑𝑡 𝐶 𝑡𝑜
Al derivar nos queda la siguiente ecuación
𝑑𝑖 𝐿2 𝑖𝑐 𝑑𝑖
0=𝐿∗ 2
+0+ +𝑅∗
𝑑𝑡 𝐶 𝑑𝑡
Quitamos la L dividiendo todo entre L y nos queda la siguiente ecuación
𝑑𝑖 𝐿2 𝑖𝑐 𝑅 𝑑𝑖
0= + + ∗
𝑑𝑡 2 𝐿𝐶 𝐿 𝑑𝑡
Y usando los valores de w y α para dejar nuestra ecuación de forma directa

𝑑𝑖 𝐿2 𝑑𝑖
0= 2
+ 2α ∗ + 𝑤2 ∗ 𝑖
𝑑𝑡 𝑑𝑡

Sacando los valores de w y α para saber de que forma será nuestra respuesta
1
𝑤= = 10
√𝐿𝐶
𝑅
𝛼 = 2𝐿 = 9
Al obtener la ecuación obtenemos:
𝑑𝑖 𝐿2 𝑑𝑖
+ 18 ∗ + 100 𝑖 = 0
𝑑𝑡 2 𝑑𝑡
Como sabemos que α>w, podemos decir que es una respuesta sub amortiguada

Para t<0
El circuito queda de la siguiente manera

Sacamos primero los valores iniciales


10𝑉
𝑖𝐿 (0− ) = = 1𝐴
10Ω
𝑉𝑐 (0− ) = 1𝐴 ∗ 6Ω = 6V
Al saber que se va a obtener una función sub amortiguada sabemos que
𝑖(𝑡) = 𝑒 −𝛼𝑡 (𝐴1 cos( 𝑤𝑑 ∗ 𝑡) + 𝐴2 sin(𝑤𝑑 ∗ 𝑡))
Sacando Wd tenemos que es igual a

𝑤𝑑 = √𝑤 2 − 𝛼 2 = 4.36 𝑟𝑎𝑑/𝑠
Sacando la ecuación nos quedaría de la siguiente manera:
𝑖(𝑡) = 𝑒 −9𝑡 (𝐴1 cos(4.36 ∗ 𝑡) + 𝐴2 sin (4.3 ∗ 𝑡)
Sabiendo que i(t) = i (0):

1𝐴 = 𝑒 −9𝑡 (𝐴1 cos(4.36 ∗ 𝑡) + 𝐴2 sin (4.3 ∗ 𝑡)


Despejando con t = 0 nos da que el valor de A1 es:
1A = 1(A1 (1) + 0) = 1 = A1
Entonces tenemos que A1 = 1
Para sacar el valor de A2 hacemos la integral, pero para ello sacamos el valor de la derivada de
la corriente en iL el cual sería:
𝑑𝑖𝐿(0+) 𝐴
= −30
𝑑𝑡 𝑠
Por lo que al sacar la derivada de toda la ecuación nos daría lo siguiente
𝐴
−30 = 9 𝑒 −9𝑡 (−𝐴1 4.36 sin(4.36 ∗ 𝑡) + 𝐴2 ∗ 4.36 cos(4.36 ∗ 𝑡))
𝑠
𝐴
−30 = −9 ∗ 𝐴1 + 4.36 ∗ 𝐴2
𝑠
Al despejar A2
−30 + 9 ∗ 𝐴1 −30 + 9
𝐴2 = = = −4.82
4.36 4.36

Por lo tanto, la ecuación nos daría:


Para t>0

𝒊(𝒕) = 𝒆−𝟗𝒕 ∗ (𝒄𝒐𝒔(𝟒. 𝟑𝟔 ∗ 𝒕) − (𝟒. 𝟖𝟐 ∗ 𝒔𝒆𝒏(𝟒. 𝟑𝟔 ∗ 𝒕)))

Grafica en Matlab
Grafica de LTSPICE
Lp
Como conclusión tengo una duda con lo diferente de las gráficas, esto debido a que en LTSPICE
no se si coloque de manera correcta la ecuación de la fuente de voltaje, además de que no se
puede iniciar una simulación de LTSPICE antes de cero, además de que el switch no se si tenga
mucha influencia en el proceso, sin embargo, si podemos ver como al principio la corriente que
va por el circuito es alta y como va disminuyendo debido a que se desconecta la fuente de
voltaje.

Common questions

Con tecnología de IA

The switch's behavior impacts analysis by introducing a marked transition: at t=0, the source of excitation is disconnected, transforming the circuit to a natural response. This affects current and voltage evolution since only stored energy in the inductor and capacitor affects subsequent dynamics, needing careful evaluation of initial conditions and transient analysis using the derived differential equation. The abrupt disconnection initiates diminishing oscillations, governed by derived parameters like 𝛼 and 𝑤, reflecting the circuit's inherent characteristics without external excitation .

The differential equation for a series RLC circuit when the switch opens at t=0 is derived by first expressing the voltage sum as 𝑉_𝑓= 𝑉_𝑟+ 𝑉_𝑐+ 𝑉_𝑙, which gives 𝑉_𝑓= 𝐿∗𝑑𝑖/𝑑𝑡 + 1/𝐶∫𝑖 𝑑𝑡 + 𝑅𝑖. Differentiating results in 0 = 𝐿∗𝑑^2𝑖/𝑑𝑡^2 + 𝑅∗𝑑𝑖/𝑑𝑡 + 𝑖/𝐶. Simplifying by dividing by L gives 0 = 𝑑^2𝑖/𝑑𝑡^2 + (𝑅/𝐿)∗𝑑𝑖/𝑑𝑡 + (1/𝐿𝐶)𝑖. With 𝛼=𝑅/2𝐿 and 𝑤=1/√𝐿𝐶, it transforms to 0 = 𝑑^2𝑖/𝑑𝑡^2 + 2α∗𝑑𝑖/𝑑𝑡 + 𝑤^2∗𝑖. Given 𝛼 > 𝑤, indicating a sub-amortized response, the differential equation becomes 𝑑^2𝑖/𝑑𝑡^2 + 18∗𝑑𝑖/𝑑𝑡 + 100𝑖=0 .

A sub-amortized response is indicated when the damping factor 𝛼 is greater than the natural frequency 𝑤 of the circuit, i.e., 𝛼 > 𝑤. In the context of the RLC circuit presented, where 𝛼=9 and 𝑤=10, it is represented in the solution as an exponentially decaying oscillation, given by 𝑖(𝑡) = 𝑒^−9𝑡(𝐴₁ cos(4.36∗𝑡) + 𝐴₂ sin(4.36∗𝑡)). The angular frequency of oscillation 𝑤𝑑 is √(𝑤²−𝛼²) = 4.36 .

Exponential decay is integral because it characterizes energy dissipation over time as oscillations in the circuit reduce. Mathematically, with a sub-amortized RLC response, it is expressed by the e^−9t term in i(t) = 𝑒^−9t(𝐴₁ cos(4.36t) + 𝐴₂ sin(4.36t)), where 9 exemplifies the decay rate determined by the damping factor 𝛼. This decay reflects the energy lost through resistive components, revealing how quickly the influence of initial conditions wanes, leading to a cessation of oscillation in steady state .

Matlab is utilized by graphing the calculated time-domain response function of the RLC circuit, allowing for visual analysis of the circuit's transient behavior. One challenge mentioned is the discrepancy between Matlab and LTSPICE simulations, possibly due to setting the equation of the voltage source incorrectly or issues initiating simulations in LTSPICE before t=0. Differences in results can lead to confusion if simulations do not properly reflect the theoretical model due to initial conditions or software constraints .

An engineer might be uncertain about the initial conditions in an LTSPICE simulation due to limitations within the software, such as not being able to correctly set the initial voltage source or configure the circuit to reflect pre-switch settings accurately. Moreover, initiating a simulation at t=0 might not capture pre-switch dynamics, introducing discrepancies between expected transient behavior and simulated results. This could lead to misalignment in observed current or voltage waveforms compared to theoretical expectations .

The value of 𝑤d (damped natural frequency) is crucial because it represents the frequency of oscillation in the underdamped response of the RLC circuit. It is calculated as 𝑤d = √(𝑤² - 𝛼²) and dictates how quickly the oscillations decrease in amplitude. Knowing 𝑤d is important for accurately predicting waveform behavior over time, such as the frequency of oscillatory components in i(t). In this scenario, the calculated 𝑤d is 4.36 rad/s, reflecting the damping effect on oscillation frequency .

Calculating initial values like 𝑖𝐿(0−) and 𝑉𝑐(0−) is essential for accurately describing the circuit's initial state before the switch changes position at t=0. They provide necessary conditions for solving the differential equation, ensuring continuity and correct initial behavior of currents and voltages. In the given scenario, 𝑖𝐿(0−) was determined to be 1A, and 𝑉𝑐(0−) was 6V, which establish immediate values around which the circuit response evolves once the system dynamics change .

To solve for constants A1 and A2 in i(t) = 𝑒^−9𝑡(𝐴₁ cos(4.36∗𝑡) + 𝐴₂ sin(4.36∗𝑡)), we use boundary conditions. Setting t=0 gives i(0)=1A, leading to A1=1. To find A2, the derivative of the current immediately after the switch, 𝑑𝑖𝐿(0+)/𝑑𝑡 = −30 A/s, is equated to the derived expression: 𝑑𝑖(𝑡)/𝑑𝑡 at t=0, resulting in an equation in terms of A₁ and A₂. Solving -30 = -9A₁ + 4.36A₂ with A₁ = 1, we find A₂ = -4.82 .

Discrepancies between Matlab and LTSPICE simulations can challenge the validation process of a circuit model by indicating possible errors in initial condition settings or the implementation of mathematical models in the tools. Such differences may arise from inaccuracies in setting voltage source equations, limitations in handling time-zero conditions, or unforeseen interaction issues in software, thus prompting further investigation into system setup and assumptions. Analyzing these disparities aids in refining the theoretical model or adjusting simulation configurations for more accurate real-world predictions .

También podría gustarte