PROYECTO INTEGRADOR DE ECUACIONES DIFERENCIALES.

Objetivo general:
Aplicar ecuaciones diferenciales en un modelo fsico (Circuito RLC).
Objetivos especficos:
Modelar un sistema fsico (Circuito RLC) mediante ecuaciones diferenciales.
Resolver el modelo por medio de ecuaciones diferenciales.
Determinar las soluciones homogneas y particular del sistema sometido a un rgimen
senoidal permanente.
Constatar las respuestas mediante un simulador.
Recursos:
Software simulador de circuitos elctricos (Multisim, Simulink, Proteus)
Software para realizar grficos de funciones (Winplot, Derive, Matlab.)
Sustento Terico

Modelado con Ecuaciones Diferenciales de Orden Superior

1. MODELOS LINEALES: PROBLEMAS CON VALORES INICIALES

Figura 1: Movimiento oscilatorio de un Resorte

Vamos a considerar varios sistemas dinmicos lineales en los que cada modelo matemtico
es una ecuacin diferencial de segundo orden con coeficientes constates junto con
condiciones iniciales especificadas en un tiempo inicial de t=0:

+
+ = (),
2

(0) = 0 , (0) = 1

g(t)= funcin de entrada, funcin de conduccin o funcin forzada.

La solucin y(t) que satisface las condiciones iniciales se llaman respuesta del sistema.
1.1. Sistemas Resorte/masa: Movimiento Libre no Amortiguado
Ley de Hooke

Figura 2: Ley de Hooke

Esta ley establece que el resorte unido a una masa ejerce una fuerza F opuesta a la
direccin de elongacin y proporcional a la cantidad de elongacin s y es expresada
en forma simple como F=ks, siendo k la constante del resorte.
El resorte se caracteriza por el numero k.

Segunda Ley de Newton

Figura 3: Resorte y Ley de Newton

Luego de unir una masa a un resorte, este se alarga una distancia s y logra una
posicin de equilibrio en la cual su peso W se equilibra mediante la fuerza
restauradora ks. El peso se define mediante W=mg, donde la masa se mide en slugs,
kilogramos o gramos y g=32 pies/s2.
La condicin de equilibrio es mg=ks . Si la masa se desplaza una cantidad x de su
posicin de equilibrio, la fuerza restauradora del resorte es k(x+s). Cabe recalcar
que esto es cuando no existe fuerzas restauradoras que actan sobre el sistema se
puede igualar la segunda ley de Newton con la fuerza neta o resultante de la fuerza
restauradora y el peso.

2
= ( + ) + = + =
2

El signo menos indica que la fuerza restauradora acta opuesta al movimiento.

Ecuacin Diferencial de un movimiento libre no amortiguado.
Dividiendo la ecuacin anterior para la masa, obtenemos:
2
+ 2 = 0
2
Esta ecuacin describe el movimiento armnico simple o movimiento libre no
amortiguado.

Ecuacin de Movimiento

La solucin de la ecuacin anterior ser:

() = 1 cos + 2
2

Periodo: T=

Frecuencia de movimiento: = = 2 ciclos por segundo

Frecuencia circular: = /

Figura 4: Movimiento Libre no Amortiguado.

1.2. Sistemas resorte/masa: Movimiento Libre Amortiguado
Ecuacin Diferencial de un Movimiento Libre Amortiguado.
Las fuerzas de amortiguamiento que actan sobre un cuerpo son consideradas
proporcionales a una potencia de la velocidad instantnea. Esta fuerza est dada por
un mltiplo constante de dx/dt. Cuando no acta otra fuerza en el sistema se
obtiene:
2

2 =

donde es una constante de amortiguamiento positiva y el signo negativo es una

consecuencia del hecho de que la fuerza de amortiguamiento acta en una direccin
opuesta al movimiento.
Dividiendo la ecuacin anterior para la masa, nos da la ecuacin diferencial del
movimiento libre amortiguado:

+ 2
+ 2 = 0
2

2= , 2 =

Donde
A continuacin, los tres casos:
CASO I: 2 2 >0

El sistema se encuentra sobreamortiguado ya que es grande respecto a la constante

k. La solucin de la ecuacin ser:
() = (1

2 2

+ (2

2 2

Esta ecuacin describe un movimiento uniforme y no oscilatorio.

CASO II : 2 2=0
El sistema se encuentra crticamente amortiguado porque cualquier ligera
disminucin en la fuerza de amortiguamiento dara como resultado un movimiento
oscilatorio.
La solucin general de esta ecuacin seria:
() = (1 + 2 )
CASO III: 2 2<0
El sistema se encuentra subamortiguado ya que el coeficiente de amortiguamiento
es pequeo comparado con la constante del resorte.
La solucin general de esta ecuacin seria:
() = (1 2 2 + 2 2 2 )

Figura 5: Movimiento Amortiguado

1.3. Sistemas resorte/masa: Movimiento forzado

Ecuacin Diferencial de Movimiento Forzado con Amortiguamiento
Ahora tomaremos en consideracin una fuerza externa f(t) que acta sobre una masa
vibrante en un resorte. La inclusin de esta fuerza nos da la ecuacin:
2

+ 2
+ 2 = ()
2

Para resolver esta ecuacin de puede usar el mtodo de coeficientes indeterminados

o variacin de parmetros.

Figura 6: Movimiento Forzado

Trminos Transitorios y de Estado Estable
Se dice que la funcin complementaria xc(t) es un trmino transitorio o solucin
transitoria, y la funcin xp(t), la parte de la solucin que permanece despus de un
intervalo de tiempo, se llama termino de estado estable o solucin de estado estable.

Ilustracin 7: Trminos Transitorios y Estables

1.4. Circuito en serie Anlogo

CIRCUITOS LRC EN SERIE:

Figura 8: Circuito LRC

Por la segunda ley de Kirchhoff, la suma de las cadas de voltaje en el inductor,
resistor y capacitor es igual al voltaje E(t) aplicado al circuito,

1
+ + = ()

Como ecuacin diferencial de segundo orden:

2
1
+
+ = ()
2

Si E(t)=0, las vibraciones elctricas del circuito estn libres, para las ecuaciones
anteriores se dice:
El circuito es Sobreamortiguado si 2

>0

El circuito es crticamente Amortiguado si 2

El circuito es Subamortiguado si 2

=0

<0

Cuando E(t)=0 y R=0, la respuesta del circuito es Armnica simple.

2. RESONANCIA SERIE EN CIRCUITOS RLC

Definimos como resonancia al comportamiento de un circuito con elementos inductivos
y capacitivos, para el cual se verifica que la tensin aplicada en los terminales del
mismo circuito, y la corriente absorbida, estn en fase. La resonancia puede aparecer en
todo circuito que tenga elementos L y C.
Por lo tanto existir una resonancia serie y otra resonancia paralelo o en una
combinacin de ambos.
El fenmeno de resonancia se manifiesta para una o varias frecuencias, dependiendo del
circuito, pero nunca para cualquier frecuencia. Es por ello que existe una fuerte
dependencia del comportamiento respecto de la frecuencia. Deviene de ello la gran
importancia de los circuitos sintonizados, especialmente en el campo de las
comunicaciones, en lo que hace a la sintonizacin de seales de frecuencias definidas o
al "filtrado" de seales de frecuencias no deseadas.
Genricamente se dice que un circuito est en resonancia cuando la tensin aplicada y
la corriente estn en fase, el factor de potencia resulta unitario.

Figura 9: Resonancia en circuito LRC

Si trazamos el diagrama de tensiones y corrientes del circuito, se verificar que la
tensin adelantar,
Atrasar o estar en fase con la corriente. Esto resulta evidente de la expresin anterior,
en la cual, para algunas frecuencias se cumplir que:
1

>
para otras frecuencias ser:

 <

En el primer caso, se comporta el circuito en forma inductiva, en el segundo, en forma

capacitiva y, adems, para alguna frecuencia, se cumplir que:
=

Para este caso, el circuito se encontrar en resonancia, ya que la impedancia ser

resistiva pura.(tensin en fase con la corriente). Este tipo de circuito se denomina
tambin Resonante en Tensiones, dado que los mdulos de las tensiones en los
componentes reactivos, son iguales pero opuestos en fase y se cancelan.
Los diagramas fasoriales son los que se dibujan a continuacin:

Figura 10: DIagramas Fasoriales

Debe observarse que cundo, el circuito estar en resonancia, el circuito se comportar
en forma resistiva pura, mientras la impedancia ser slo la resistencia del circuito, y,
por consiguiente, la corriente ser mxima.
Frecuencia de resonancia
Se obtiene muy fcilmente, ya que la componente imaginaria de la impedancia deber
ser nula, para que el circuito se comporte como resistivo puro. Para este caso simple,
ser:
Se ve en esta ltima expresin, que la frecuencia de resonancia, ser siempre la misma
en la medida que no cambie el producto LC.

1
0

0 =
0 2 =
0 =

, 0 = 20
0 =

1
2

3. RESPUESTA FORZADA Y NATURAL DE UN SISTEMA.

La Respuesta Forzada es la solucin cuando todas las condiciones iniciales son nulas y
el sistema se encuentra sometido a la seal de entrada u(t).
Normalmente es difcil determinar. Existen distintos mtodos entre los que cabe citar el
mtodo de los coeficientes indeterminados, en el que la solucin particular depende
mucho del tipo de funcin u(t) y se encuentra tabulada.
Como se ha dicho anteriormente, la Solucin Completa para un sistema descrito por una
Ecuacin Diferencial con coeficientes constantes se obtiene sumando la Respuesta
Libre y la Respuesta Forzada.

DESARROLLO
1. Seleccione la cdula de uno de los integrantes del grupo, los tres ltimos dgitos
diferentes de cero corresponden desde el ms significativo al menos significativo a los
valores de la inductancia en mili henrios, capacitancia en mili faradios y resistencia en
ohmios.
El nmero de la cedula de identidad seleccionado fue el 010656089-9
Por lo cual los valores asumidos para el desarrollo del trabajo fueron:
Inductancia: 9mH
Capacitancia: 9mF
Resistencia: 8 ohm

2. Calcule la solucin transitoria del voltaje en el condensador del circuito serie

conformado por la inductancia y capacitancia si el condensador tiene un voltaje inicial
de 10.v.
L
C

a) Determine la frecuencia natural de oscilacin mediante la solucin transitoria.

=
=

1

1

0.009 0.009

= .
Ahora procedemos a calcular la solucin transitoria del voltaje en el condensador
Considerando que: (0) = 10
(0) = 0
Entonces:
2 1
2+ =0

0.009 + 0.009 = 0
dividimos para 0.009
+

1
=0
0.009 0.009

+ 12345.679 = 0
La ecuacin auxiliar es:
= 111.11
La solucin de la ecuacin es:
= 1(111.111) + 2cos(111.111)

Considerando la condicin:
1
= 10

Entonces (0) = 10 = 10 0.009 = 0.09

Por lo tanto
0.09 = 1(111.111(0)) + 2cos(111.111(0))
0.09 = 1
= 111.1111 cos(111.111(0)) 111.112(111.11(0))
Si () = 0 :
2 = 0
Entonces hallo la ecuacin del voltaje en el condensador
() =

0.09(111.111)
=(
)

0.009
() = (. )

b) Dibuje en un software la respuesta del voltaje en el condensador en dominio del

tiempo.

Grafica de solucin de ecuacin

c) En la grfica realizada en 2.b determine el tiempo entre dos cruces por cero con
la funcin descendente prximos. A este tiempo se lo llama periodo (T). Vea
figura 2.

Figura 2: Periodo de una onda.

El tiempo entre los dos cortes ascendentes y descendentes es :

= .
d) Calcule la frecuencia (f=1/T) segn la grfica.
=

1
0.056

= .

e) Compare los dos valores obtenidos en los puntos 2.a y 2.d.

Frecuencia natural = 111.111Hz
Frecuencia segn la grfica = 17.85

f) Arme en un simulador electrnico el circuito.

9mH

9mF

g) Grafique mediante un osciloscopio en el simulador la forma de onda del voltaje

sobre el condensador.
XSC1
Ext T rig
+
_
B

9mH

9mF

Grafica de simulacin

h) Compare las grficas obtenidas en los puntos 2.b y 2.g.

Grafica obtenida en el simulador

Grafica obtenida en derive

La graficas obtenida mediante la solucin dela ecuacin diferencial del voltaje sobre el
condensador es la misma que la obtenida en el simulador.

3. Incluya al circuito anterior la resistencia formando un circuito RLC serie en el que

el voltaje inicial en el condensador es 10.V

a) Encuentre la solucin de la ecuacin diferencial que define el voltaje sobre la

resistencia.

2
1
+
+ = ()
2

1
0.009 2 + 8
+
=0

0.009
0.009

+
8
+ 111.11 = 0
2

Se usa el mtodo de coeficientes constantes homognea para poder encontrar la ecuacin

auxiliar requerida:
0.009 2 + 8 + 111.11 = 0
Se obtiene los valores de r con la formula
2 4
1.2 =
2
1 = 14.44
2 = 874.44
Resultado de la ecuacin es:

() = 1 14.44 + 2 874.44

Obtenemos el valor de q con esta ecuacin:

1
10 =

0.009
= 0.09
= 0
=

Con esto seran los valores iniciales

(0) = 0.09

 (0) = 0
Derivada de la ecuacin obtenida
() = 14.44 1(5.35 107 ) 874.44 2(1.72 10380 )

Reemplazamos los valores de frontera

0.09 = 1 14.440 + 2 874.440
0.09 = 1 + 2
0.09 2 = 1
Reemplazamos los valores de frontera en la otra ecuacin
0 = 14.44 1(5.35 107 )0 874.44 2(1.72 10380 )0
0 = 14.44 1 874.44 2
Reemplazamos la ecuacin de c1
0 = 14.44 (0.09 2) 874.44 2
0 = 1.2996 860 2
1.2996
2 =
860
2 = 0.0015
1 = 0.09 + 0.0015
1 = 0.0915
Solucin de la ecuacin diferencial
() = 0.091 14.44 0.0015 874.44
Voltaje en funcin de la resistencia

=
= [1.32144 14.44 + 1.3116 874.44 ] 8
= . . + . .

b) Dibuje en un software la respuesta obtenida en el literal 3.a

Grafica de solucin de ecuacin

c) Arme el circuito en un simulador y grafique con el osciloscopio el voltaje en la
resistencia.
R
8

L3
9mH
C3
9mF
IC=10V

Grafica de simulacin

d) Compare las grficas de los puntos 3.b y 3.c

Las dos graficas son idnticas y muestran el comportamiento de la tensin de la
resistencia en funcin de tiempo, debido a que inicialmente el condensador est
cargado entonces este hace de fuente que alimenta al circuito pero como no es una
fuente constante a medida que transcurre el tiempo el condensador se descarga
entonces todo el sistema quedara sin alimentacin por ende la tensin en la
resistencia durara hasta que el condensador se descargue.

e) Qu tipo de amortiguamiento tiene el sistema?

Es un sistema sobreamortiguado

4. Incluya en serie a todos los elementos una fuente de voltaje senoidal de 10.V.

a) Configure la fuente a la frecuencia natural del sistema obtenida en el literal

2.a.
R1
8
V2
10 Vpk
111.11 Hz
0

L1
9mH

C1
9mF

b) Calcule y dibuje tanto en el simulador con el osciloscopio como en un software

para grficas el voltaje en la resistencia 4.a.

Tenemos la ecuacin que define los voltajes en los respectivos elementos:

2
1
2 +
+ = ()

Sustituimos los valores con los datos dados

2

1
0.009 2 + 8
+
= 10()

0.009
Ahora dividimos para 0.009
2
8
1
10
+
+

=
()
2 0.009 0.009 0,009
0,009
2

+ 888.88
+ 12345.67 = 1111.11()
2

La ecuacin auxiliar para resolver resta ecuacin diferencial seria:

2 + 888.88 + 12345.67 = 0
Las soluciones de esta ecuacin son:
1 = 14.11
2 = 874.8
Entonces la solucin complementaria estar dada por la siguiente expresin:
= 1 1 + 2 2
= . .
Ahora se debe hallar la solucin particular la cual es:
= () + cos()
= () ()
= 2 () 2 ()
Ahora sustituimos en la ecuacin diferencial inicial
2 ( ) + 2 ( )
+ 888.88( ( ) ( )) + 12345.67
( ( ) + ( )) = 1111.11 ( )
Sustituimos = 2
En este caso tomamos la frecuencia natura = 111.11
Entonces = 2 111.11 = 698.124

La ecuacin diferencial inicial seria:

698.1242 (698.124) + 698.1242 (698.124)
+ 888.88(698.124(698.124) 698.124 (698.124))
+ 12345.67((698.124) + (698.124))
= 1111.11(698.124)
Agrupando y trminos podemos encontrar los coeficientes A y B

= 0.0008746749872
= 0.001086157559]
Ahora sustituimos en la ecuacin particular:
= 0.0008746749871 (698.1247194 ) 0.001086157559
(698.1247194 )
La solucin de la ecuacin diferencial es:
() = +
() = 1 1 + 2 2 + 0.0008746749871 (698.1247194 )
0.001086157559 (698.1247194 )
Con la condiciones de iniciales q(0)=0
0 = 1 10 + 2 2 0 + 0.0008746749871 (698.1247194 0)
0.001086157559 (698.1247194 0)
Tenemos:
1 + 2 = 0.001086157559
Ahora derivamos la q(t) con respecto a t
() = 14,111 14.11 + 874.8 2 (1.200016537 10380 ) 0.6106322299
(698.1247194 ) 0.7582734411 (698.1247194
Sustituyendo la condiciones (0) = 0
Se tiene:
14.11 1 + 874.8 2 = 0.6106322299
Resolviendo las ecuaciones se obtiene C1 yC2
1 + 2 = 0.001086157559

14.11 1 + 874.8 2 = 0.6106322299

1 = 0.0003944955821
2 = 0.0006916619767]
Ahora sustituimos en la solucin de la ecuacin:
= 0.0003944955821 ^( 14.11 ) + 0.0006916619766
(1.200016537 10380 ) 0.001086157559 (698.1247194
) + 0.0008746749871 (698.1247194 )
La derivada

= 0.005566332663 14.11 + 0.6050658971 (1.200016537 10380 )

0.6106322299 (698.1247194 ) 0.7582734411

(698.1247194 )
Ahora pa encontrar el voltaje en la resistencia segn la ley de ohm
=
= . . + . (. ) .
(. ) . (. )

Grafica de solucin de ecuacin

Grafica de simulacin
c) Repita 4.a y 4.b configurando la frecuencia de la fuente a la mitad de la frecuencia
natural.
Para este punto repetimos los clculos realizados en el punto a y b entonces:

2 ( ) + 2 ( )
+ 888.88( ( ) ( )) + 12345.67
( ( ) + ( )) = 1111.11 ( )
Sustituimos = 2
En este caso tomamos la mitad de la frecuencia natura = 55.5
Entonces = 2 55.5 = 349.06
2

349.06 (349.06) + 349.06 (349.06)

+ 888.88(349.06(349.06) 698.124 (349.06))
+ 12345.67((349.06) + (349.06)) = 1111.11(349.06)
Agrupando y trminos podemos encontrar los coeficientes A y B

= 0.001304773440
= 0.003017118410]
= 0.001304773440 (349.0309438) 0.003017118410
(349.0309438)

() = +
() = 1 1 + 2 2 + 0.001304773440 (349.0309438)
0.003017118410 (349.0309438)
Con la condiciones de iniciales q(0)=0
0 = 1 10 + 2 2 0 + 0.001304773440 (349.0309438 0)
0.003017118410 (349.0309438 0)
Tenemos:
1 + 2 = 0.003017118410
Ahora derivamos la q(t) con respecto a t
() = 14,11 1 14.11 + 874.8 2 (1.200016537 10380 ) 0.4554063052
(349.0309437 ) 1.053067686 (349.0309437 )
Sustituyendo la condiciones (0) = 0
14.11 1 + 874.8 2 = 0.4554063051
Resolviendo las ecuaciones se obtiene C1 yC2
1 + 2 = 0.003017118410
14.11 1 + 874.8 2 = 0.4554063051
1 = 0.002537462826
2 = 0.0004796555836]
Ahora sustituimos en la solucin de la ecuacin:
= 0.002537462826 14.11 + 0.0004796555835 (1.200016537 10380 )
0.003017118410 (349.0309437 ) + 0.001304773440
(349.0309437 )

La derivada

= 0.03580360047 14.11 + 0.4196027044 (1.200016537 10380 )

0.4554063050 (349.0309436 ) 1.053067685

(349.0309436 )

Ahora para encontrar el voltaje la resistencia segn la ley de ohm

VR = I R
= . . + . (. ) .
(. ) . (. )

Grafica de solucin de ecuacin

Grafica de simulacin

d) Repita 4.a y 4.b configurando la frecuencia de la fuente al doble de la frecuencia

natural.

2 ( ) + 2 ( )
+ 888.88( ( ) ( )) + 12345.67
( ( ) + ( )) = 1111.11 ( )
Sustituimos = 2
En este caso tomamos el doble de la frecuencia natural = 222.22
Entonces = 2 22.22 = 1396.12
2

1396.12 (1396.12) + 1396.12 (1396.12)

+ 888.88(31396.12(1396.12) 1396.12 (1396.12))
+ 12345.67((1396.12) + (1396.12))
= 1111.11(1396.12)
Agrupando y trminos podemos encontrar los coeficientes A y B
= 0.0004044821660
= 0.0002558809283]

= 0.0004044821660 (1396.249438 ) 0.0002558809283

(1396.249438 )
() = +
() = 1 1 + 2 2 + 0.0004044821660 (1396.249438 )
0.0002558809283 (1396.249438 )
Con la condiciones de iniciales q(0)=0
0 = 1 1 0 + 2 2 0 + 0.0004044821660 (1396.249438 0)
0.0002558809283 (1396.249438 0)

Tenemos:
1 + 2 = 0.0002558809283
Ahora derivamos la q(t) con respecto a t
() = 14,11 1 14.11 + 874.8 2 (1.200016537 10380 ) ) 0.5647579969
(1396.249438 ) 0.3572736023 (1396.249438 )

Sustituyendo la condiciones (0) = 0

14.11 1 + 874.8 2 = 0.5647579969
Resolviendo las ecuaciones se obtiene C1 yC2
1 + 2 = 0.0002558809283
14.11 1 + 874.8 2 = 0.5647579969
1 = 0.0003960930892
2 = 0.0006519740173
Ahora sustituimos en la solucin de la ecuacin:
= 0.0003960930891 14.11 + 0.0006519740173 (1.200016537 10380 )
0.0002558809283 (1396.249438 ) + 0.0004044821660
(1396.249438 )

La derivada

= 0.005588873487 14.11 0.5703468703 (1.200016537 10380 )

+ 0.5647579969(349.0309436 )0.3572736023
(349.0309436 )
Ahora para encontrar el voltaje la resistencia segn la ley de ohm
VR = I R
= . . . (. ) + .
(. ) + . (. )

Grafica de solucin de ecuacin

Grafica de simulacin
e) Cul de las tres frecuencias presenta un voltaje en la resistencia ms alto?
En la frecuencia menor es decir la mitad de la frecuencia natural
f) A qu se debe este fenmeno?
Esto se debe a que mientras menor se la frecuencia el periodo ser menor entonces
le fuente deber suministrar la misma cantidad de voltaje pero en un tiempo menor
entonces para cumplir el ciclo se deber extender la onda en amplitud es decir el
voltaje ser mayor.

CONCLUSIONES
En un circuito RLC se pueden dar 3 casos de amortiguamiento que son:
subamortiguado, sobreamortiguado y crticamente amortiguado, dependiendo de la
relacin entre los componentes del circuito.
El comportamiento de un circuito RLC subamortiguado se asemeja al comportamiento
mecnico de un resorte, al liberarlo despus de aplicarle una fuerza externa.
Como pudimos ver las ecuaciones diferenciales de segundo orden nos ayudan a resolver
con
ayuda
de
los
mtodos
aprendidos
un
circuito
RLC.
Podemos ver que la solucin de las ecuaciones, es decir la grfica nos sale idntica a la
grfica obtenida en el osciloscopio. Gracias a estas ecuaciones podemos hacer un
estudio mas a fondo de los problemas a estudiar..

BIBLIOGRAFA:
[1] ZILL DENNIS, Matemticas Avanzadas para Ingeniera, Edit. Mcgraw-Hill, Mxico.
4ta edicin. 2012.
[2] NAGLE, SAFF,SNIDER, Ecuaciones Diferenciales y problemas von valores en la
frontera, / Edit. PEARSON, 2005 ,4 Edicin