Análisis de sistemas dinámicos.
Para comenzar el análisis de los sistemas dinámicos, veamos un sistema conocido y sencillo, el sistema masa-resorte.
Figura 5: Sistema masa-resorte.
¿De qué elementos se compone este sistema?
Masa y resorte.
¿Qué fuerzas actúan sobre el sistema?
Fuerza gravitacional y del resorte.
Aplicando la segunda ley de Newton, se obtiene la ec. de mov.
\[
m \ddot{y}=\sum Fuerzas
\]
es decir,
\begin{equation}
m \ddot{y}+ky=mg
\label{ecMR}
\end{equation}
De la Figura 5 se puede ver que,
\begin{equation}
y=x+\delta
\label{ecdelta}
\end{equation}
además note que al colocarse la masa el resorte sufre un desplazamiento de su posición en $y=x$ hasta su nueva posición $x+\delta$, por lo tanto,
\[
k \delta=mg
\]
sustituyendo la ecuación (\ref{ecdelta}) en (\ref{ecMR}) se obtiene,
\begin{equation}
m \ddot{y}+k(x+\delta)=mg
\end{equation}
La ecuación (\ref{ecdelta}) se deriva con respecto al tiempo y al ser $\delta$ constante se obtiene,
\[
\dot{y}=\dot{x}+0
\]
\[
\ddot{y}=\ddot{x}
\]
con lo cual se nota que $\ddot{y}=\ddot{x}$ y por lo tanto la ecuación (\ref{ecMR}) puede escribirse de la forma,
\begin{equation}
m \ddot{x}+kx=0
\label{ecMRx}
\end{equation}
Supóngase que la masa mostrada en la Figura 5 se jala hacia abajo y después se suelta, la distancia que se movió hacia abajo se conoce como condición inicial (C.I.) del experimento, también puede darse una velocidad inicial, y en ambos casos, la masa oscilará y se moverá periódicamente, a esto se le conoce como vibración libre.
Figura 6: Sistema masa-resorte con condiciones iniciales.
Ya se tiene la ecuación de movimiento de un sistema masa resorte, dada en \ref{ecMRx}, ahora se debe resolver esta ecuación para poder conocer y analizar el comportamiento de este sistema.
Para este efecto, hay varios métodos, se estudiará la solución por forma exponencial, por lo tanto, suponga que,
\begin{equation}
x(t)=Ce^{\lambda t}
\label{exp1}
\end{equation}
donde $C$ y $\lambda$ son constantes, si se deriva la ecuación (\ref{exp1}) con respecto al tiempo, se obtienen,
\begin{eqnarray}
\dot{x}(t)&=& C\lambda e^{\lambda t}\\
\ddot{x}(t)&=& C \lambda^{2} e^{\lambda t}
\label{exp2}
\end{eqnarray}
sustituyendo (\ref{exp1}) y (\ref{exp2}) en la ecuación (\ref{ecMRx}) se obtiene,
\[
m C \lambda^{2} e^{\lambda t}+kCe^{\lambda t}=0
\]
que al ser divida por $ Ce^{\lambda t}$, se obtiene,
\begin{equation}
m \lambda^{2}+k =0.
\label{eccar}
\end{equation}
A la ecuación (\ref{eccar}) se le conoce como la ecuación característica del sistema masa-resorte. Y se pueden obtener las raíces del sistema al despejar a $\lambda$,
\[
\lambda^{2}=-\frac{k}{m}
\]
las raíces del sistema son,
\[
\lambda_{1}=i \sqrt{\frac{k}{m}} \ \ \ \ \lambda_{2}=-i \sqrt{\frac{k}{m}}
\]
sustituyendo estas dos raíces en la ecuación (\ref{exp1})
\begin{equation}
x(t)=C_{1}e^{i\left(\sqrt{\frac{k}{m}}\right)t}+C_{2}e^{-i\left(\sqrt{\frac{k}{m}}\right)t}
\label{solexp}
\end{equation}
se nota que esta es la solución general del sistema masa-resorte y tiene dos constantes arbitrarías $C_{1}$ y $C_{2}$. Por lo cual ahora se debe conocer el valor de estas constantes para obtener la solución de este sistema en función de los valores conocidos de los elementos que componen al sistema masa-resorte.
Para esto, recuerde que una función exponencial, también puede expresarse en términos de funciones seno y coseno,
\begin{eqnarray*}
e^{i \omega t} &=& \cos \omega t + i \sin \omega t \\
e^{-i \omega t} &=& \cos \omega t – i \sin \omega t
\end{eqnarray*}
al sustituir esto en la ecuación (\ref{solexp}), se tiene,
\begin{eqnarray*}
x(t)&=& C_{1} \left( \cos \sqrt{\frac{k}{m}}t + i \sin \sqrt{\frac{k}{m}}t\right) + C_{2} \left( \cos \sqrt{\frac{k}{m}}t – i \sin \sqrt{\frac{k}{m}}t\right)\\
&=& i (C_{1}- C_{2}) \sin \sqrt{\frac{k}{m}}t + (C_{1} + C_{2}) \cos \sqrt{\frac{k}{m}}t \\
&=& A \sin \sqrt{\frac{k}{m}}t + B \cos \sqrt{\frac{k}{m}}t
\end{eqnarray*}
donde,
\[
A=i(C_{1}-C_{2}), \ \ \ \ B=(C_{1}+C_{2})
\]
Por lo tanto la solución general mostrada en (\ref{solexp}) también puede expresarse en la forma,
\begin{eqnarray}
x(t) &=& A \sin \left(\sqrt{\frac{k}{m}}t\right) + B \cos \left(\sqrt{\frac{k}{m}}t\right)
\label{solexp2}
\end{eqnarray}
donde es importante recordar que las constantes A y B depende de las condiciones iniciales, es decir de la posición y velocidad que tienen en el tiempo $t=0$, por lo tanto la posición inicial se denota como $x(0)$ y la velocidad inicial como $\dot{x}(0)$.
Conocer el valor de estas constantes es sencillo, solo se necesita saber lo que sucede en el sistema en el tiempo $t=0$, por lo tanto sustituimos este valor en la ecuación (\ref{solexp2}),
\begin{eqnarray*}
x(0) &=& A \sin \left(\sqrt{\frac{k}{m}}(0)\right) + B \cos \left(\sqrt{\frac{k}{m}}(0)\right)
\end{eqnarray*}
y se obtiene que,
\begin{equation}
B=x(0)
\label{CB}
\end{equation}
Por lo tanto el valor de $B$ es igual a la posición inicial $x(0)$ del sistema.\\
El valor de $A$ puede obtenerse ahora al utilizar la velocidad inicial del sistema $\dot{x}(0)$, por lo cual, es necesario derivar la ecuación (\ref{solexp2}) con respecto al tiempo,
\begin{eqnarray}
\dot{x}(t)&=& A \sqrt{\frac{k}{m}} \cos \left(\sqrt{\frac{k}{m}}t\right) – B \sqrt{\frac{k}{m}}\sin \left(\sqrt{\frac{k}{m}}t\right)
\label{dersol2}
\end{eqnarray}
en donde nuevamente para conocer la velocidad que el sistema tenía en el tiempo $t=0$, es necesario sustituir este tiempo en la ecuación (\ref{dersol2}),
\begin{eqnarray*}
\dot{x}(0)&=& A \sqrt{\frac{k}{m}} \cos \left(\sqrt{\frac{k}{m}}(0)\right) – B \sqrt{\frac{k}{m}}\sin \left(\sqrt{\frac{k}{m}}(0)\right)
\end{eqnarray*}
y se obtiene,
\begin{eqnarray*}
\dot{x}(0)&=& A \sqrt{\frac{k}{m}}
\end{eqnarray*}
lo cual nos indica que,
\begin{eqnarray}
A &=& \sqrt{\frac{m}{k}}\dot{x}(0)
\label{CA}
\end{eqnarray}
la constante $A$ depende de la masa, el resorte y la velocidad inicial que tenga el sistema.
Sustituyendo el valor de (\ref{CB}) y (\ref{CA}) en la ecuación (\ref{solexp2}), tenemos la solución general del sistema en términos que depende solamente de elementos conocidos como su masa, resorte y condiciones iniciales,
\begin{eqnarray}
x(t)&=& \dot{x}(0) \sqrt{\frac{m}{k}}\sin \sqrt{\frac{k}{m}}t + x(0) \cos \sqrt{\frac{k}{m}}t
\label{solci}
\end{eqnarray}
Al movimiento periódico que describe esta ecuación se le denomina movimiento armónico simple.
De la ecuación (\ref{solci}) se puede notar que la frecuencia con la que se moverá el sistema está dada por,
\[
\sqrt{\frac{k}{m}}
\]
esta frecuencia se conoce como la frecuencia natural del sistema y se representa por,
\begin{equation}
\omega_{n} =\sqrt{\frac{k}{m}}
\label{fnatural}
\end{equation}
Note que la ecuación (\ref{ecMRx}) puede dividirse entre la masa $m$ y quedar de la forma,
\begin{equation}
\ddot{x}+\frac{k}{m}x=0
\end{equation}
por lo tanto puede reescribirse como,
\begin{equation}
\ddot{x}+\omega_{n}^{2}x=0
\end{equation}
Sistema masa-resorte-amortiguador.
Hasta ahora se analizó el sistema masa-resorte, en el cual se puede ver un movimiento oscilatorio, debido a que no hay ningún elemento que absorba la energía del sistema.
Figura 7: Sistema masa-resorte-amortiguador.
El elemento que absorbe la energía del sistema se conoce como amortiguador y se introduce ahora en el sistema mostrado por la figura 7, se analiza el efecto que tiene en este nuevo sistema, al cual se le conoce como sistema masa-resorte-amortiguador (MRA).
La ecuación de movimiento para un sistema MRA, está dada por,
\[
\sum Fuerzas=m\ddot{x}+kx+b\dot{x},
\]
esta ecuación describe el movimiento pero también es el modelo matemático del sistema,
\begin{equation}
m\ddot{x}+b\dot{x}+kx=f.
\label{MRAE1}
\end{equation}
Nuevamente para encontrar la solución del sistema, suponga que,
\begin{equation}
x(t)=C_{1}e^{s_{1}t}+C_{2}e^{s_{2}t}+C_{3},
\label{MRAS1}
\end{equation}
recuerde que para encontrar la solución de este sistema, primero se deben obtener la primera y segunda derivada de (\ref{MRAS1}), después se sustituirán en (\ref{MRAE1}) y se agruparán para encontrar la solución.
La primer derivada de (\ref{MRAS1}) es,
\begin{equation}
\dot{x}(t)=\lambda_{1}C_{1}e^{\lambda_{1}t}+\lambda_{2}C_{2}e^{\lambda_{2}t}
\label{der1}
\end{equation}
mientras que la segunda derivada de (\ref{MRAS1}) es,
\begin{equation}
\ddot{x}(t)=\lambda_{1}^{2}C_{1}e^{\lambda_{1}t}+\lambda_{2}^{2}C_{2}e^{\lambda_{2}t}
\label{der2}
\end{equation}
Sustituyendo (\ref{der1}) y (\ref{der2}) en (\ref{MRAS1}) y reagrupando términos, se tiene,
\[
(m\lambda_{1}^{2}+c\lambda_{1}+k)C_{1}e^{\lambda_{1}t}+(m\lambda_{2}^{2}+c\lambda_{2}+k)C_{2}e^{\lambda_{2}t}+kC_{3}=f
\]
Donde se puede observar que, la solución del sistema puede obtenerse de la siguiente forma,
\[m s_{1}^{2}+bs_{1}+k=0
\]
\[m s_{2}^{2}+bs_{2}+k=0
\]
\[C_{3}=\frac{f}{k}
\]
si $f=0$, entonces se puede encontrar la ecuación característica de este sistema sin ser afectado por una fuerza externa,
\[
ms^{2}+bs+k=0
\]
y las raíces son,
\[
s=\frac{-b \pm \sqrt{b^{2}-4mk}}{2m}=-A \pm Bi
\]
si las raíces son complejas, entonces,
\[
A=-\frac{b}{2m} \ \ \ B=\sqrt{\frac{k}{m}-\left(\frac{b}{2m}\right)^{2}}
\]
Entonces la solución del sistema es,
\[
x(t)=e^{-At}\left[x(0)\cos(Bt)+\frac{\dot{x}(0)+Ax(0)}{B}\sin(Bt)\right]
\]
Donde la oscilación de la respuesta libre del sistema está dada por B y se le conoce como frecuencia natural amortiguada y se denota con el símbolo $\omega_{d}$.
\begin{equation}
\omega_{d}=\sqrt{\frac{k}{m}-\left(\frac{b}{2m}\right)^{2}}
\label{fnamor}
\end{equation}
Note que en esta ecuación se tiene la frecuencia natural vista en \ref{fnatural} y además un término que incluye el coeficiente de amortiguamiento, y este valor puede hacer que la frecuencia natural amortiguada sea real, sea cero o sea imaginario ( un valor negativo dentro de la raíz cuadrada), es decir,
- si b (coeficiente de amortiguamiento) es cero, entonces,
\[
\omega_{n}=\omega_{d}
\] - si b es grande, entonces, $\omega_{d}=0$ o imaginario, por lo tanto las raíces serán reales y no habrá oscilación.
El valor de b para el cual $\omega_{d}=0$ es,
\[ b=2\sqrt{mk}\] - Este valor se conoce como amortiguamiento crítico, ya que es el límite donde la respuesta libre del sistema puede o no oscilar.
Si $b<2\sqrt{mk}$ el sistema oscila.
Si $b>2\sqrt{mk}$ la respuesta del sistema es exponencial.
Razón de amortiguamiento.
Por lo tanto puede encontrarse una razón de amortiguamiento, entre el coeficiente de amortiguamiento y el valor que hace cero la frecuencia natural amortiguada, esa razón es,
\[ \zeta=\frac{b}{2\sqrt{mk}} \]
y dependiendo el valor del mismo, se sabrá si el amortiguamiento es grande, pequeño o crítico.
° Críticamente amortiguado:
Raíces repetidas. $b=2\sqrt{mk}.$ Es decir $\zeta=1$.
°Sobreamortiguado:
Cuando $b>2\sqrt{mk}$, Raíces distintas y reales. Es decir $\zeta>1$. Comportamiento exponencial.
°Subamortiguado:
Cuando $b<2\sqrt{mk}$, Raíces complejas. Es decir $\zeta<1$. Se presentan oscilaciones.