Representación en Variables de Estado de un Sistema Masa-Resorte-Amortiguador.

 

Las variables de estado son la herramienta más poderosa de la Ingeniería de Control Moderna, ya que no está limitada a sistemas lineales como sí o está el método hasta ahora visto, La Transformada de Laplace.

Las variables de estado en el caso del sistema masa-resorte-amortiguador de la Figura 8, nos permitirá reescribir un sistema de segundo orden en un sistema de primer orden. El siguiente material fue obtenido del video: State-Space Representation

Figura 8

Seleccionando nuevamente el desplazamiento como la coordenada generalizada, la ecuación de movimiento del sistema es la siguiente:

El objetivo es expresar esta ecuación en una forma equivalente que tiene la siguiente forma:

Aquí el vector es un Vector de Estado, y X1, X2, son variables de estado que sustituyen a la original variable generalizada X y, más importante, a sus derivadas. El describir el sistema en forma de matrix, ofrecerá la enorme ventaja de utilizar el poder de las computadoras para procesar información y ejecutar análisis de datos presentados en forma matricial (Matrix Algebra).

Las ecuaciones encerradas en círculos amarillos muestran como la primera forma de escribir es la forma compacta de escribir las ecuaciones para y.

El primer paso es definir las variables de estado:

Este procedimiento nos permite obtener de inmediato la primera ecuación de estado :

…..por tanto

El segundo paso consiste en forzar al coeficiente que acompaña al orden más alto, el coeficiente líder, a ser igual a la unidad. Para ello, en nuestro caso, se divide la ecuación de movimiento original entre m (y en general, entre el valor que ocupe ese lugar):

En el tercer paso se despeja la derivada de mayor orden:

El cuarto paso consiste en sustituir las derivadas de la variable original por sus ya asignadas variables de estado:

Y así hemos encontrado la segunda ecuación de estado:

….

Y así hemos completado el objetivo. La ecuación de movimiento original puede ser expresado como variables de estado en la siguiente forma:

Written by: Larry Francis Obando – Technical Specialist

Escuela de Ingeniería Eléctrica de la Universidad Central de Venezuela, Caracas.

Escuela de Ingeniería Electrónica de la Universidad Simón Bolívar, Valle de Sartenejas.

Escuela de Turismo de la Universidad Simón Bolívar, Núcleo Litoral.

Contact: Ecuador (Quito, Guayaquil, Cuenca)

WhatsApp: 00593984950376

email: dademuchconnection@gmail.com

Copywriting, Content Marketing, Tesis, Monografías, Paper Académicos, White Papers (Español – Inglés)

Dinámica de un Sistema Masa-Resorte-Amortiguador

Descargar en PDF: Dinámica de un Sistema Masa-Resorte-Amortiguador

Elementos básicos de un sistema mecánico.

Los elementos básicos de todo sistema mecánico son la masa, el resorte y el amortiguador. El estudio del movimiento en sistemas mecánicos se corresponde con el análisis de sistemas dinámicos. En robótica, por ejemplo, la palabra Forward Dynamic se refiere a lo que le sucede a los actuadores cuando le aplicamos a los mismos ciertas fuerzas y torques.

La masa, el resorte, el amortiguador, son actuadores elementales de un sistema mecánico.

En consecuencia, para controlar el robot es necesario conocer muy bien la naturaleza del movimiento de un sistema masa-resorte-amortiguador.

Además, este sistema elemental se presenta en numerosos campos de aplicación, de allí la importancia de su análisis. De nuevo, en robótica, cuando se habla de Inverse Dynamic, se habla sobre el cómo hacer que el robot se mueva de una manera deseada, cuáles fuerzas y torques debemos aplicar sobre los actuadores para que nuestro robot se mueva de una manera particular.

Antes de realizar el Análisis Dinámico de nuestro sistema masa-resorte-amortiguador, debemos obtener su modelo matemático. Éste es el primer paso a ejecutar por toda persona que pretenda conocer a profundidad la dinámica de un sistema, especialmente el comportamiento de sus componentes mecánicos.

Iniciaremos nuestro estudio con el modelo de un sistema masa-resorte.

Esto es conveniente por el motivo siguiente. Todos los sistemas mecánicos presentan una naturaleza en su movimiento que le impulsa a oscilar, como cuando un objeto pende de un hilo en el techo y con la mano lo empujamos. O un zapato sobre una plataforma con resortes. Es bueno saber qué función matemática es la que mejor describe ese movimiento.

Pero resulta que las oscilaciones de nuestro ejemplos no son infinitas. Existe una fuerza de roce que amortigua el movimiento. En el caso del objeto que cuelga de un hilo es el aire, un fluido. Por lo que luego de estudiar el caso de un sistema ideal masa-resorte, sin amortiguación, pasaremos a considerar dicha fuerza de roce y añadir a la función ya encontrada un nuevo factor que describa el decaimiento del movimiento.

 

Sistema Masa-Resorte.

 

Fuente: Física. Robert Resnick

La dinámica de un sistema se representa en primer lugar mediante un modelo matemático compuesto por ecuaciones diferenciales. En el caso de el sistema masa-resorte, dicha ecuación es la siguiente (en construcción):

Esta ecuación se conoce como Ecuación de Movimiento de un Oscilador Armónico Simple. Veamos de donde se deriva.

Si nuestra intención es obtener una fórmula que describa la fuerza que ejerce un resorte en contra del desplazamiento que lo estira o lo encoge, la mejor manera es visualizando la energía potencial que se inyecta al resorte cuando tratamos de estirarlo o encogerlo. La siguiente gráfica describe cómo se comporta esta energía en función del desplazamiento horizontal:

A medida que la masa m de la figura anterior, sujeta al extremo del resorte como se muestra en la figura 5, se aleja del punto X=0 de relajación del resorte en sentido positivo o negativo, la energía potencial se acumula y aumenta en forma parabólica, llegando a un valor superior de energía E del sistema en la máxima elongación o compresión del resorte. La ecuación matemática que en la práctica describe mejor esta forma de curva, incorporando una constante K para la propiedad física del material que aumenta o disminuye la inclinación de dicha curva, es la siguiente:

La fuerza se relaciona con la energía potencial de la siguiente manera:

Por lo tanto:

Tiene sentido ver que F(x) es inversamente proporcional al desplazamiento de la masa m. Porque está claro que si estiramos el resorte, o lo encogemos, esta fuerza se opone a dicha acción, intentando devolver al resorte a su posición relajada o natural. Por ello se le llama fuerza de restitución. La ecuación anterior es conocida en la academia como La Ley de Hooke, o ley de la fuerza para resortes. La siguiente es una gráfica representativa de dicha fuerza, en relación con la energía como se ha venido mencionando, sin intervención de fuerzas de roce (amortiguación), por lo que se le conoce como Oscilador Armónico Simple. Es importante recalcar la relación proporcional entre desplazamiento y fuerza, pero con pendiente negativa, y que, en la práctica, es más compleja, no lineal.

Fuente: Física. Robert Resnick

Para una análisis animado del resorte, corto, sencillo pero contundente, recomiendo observar los videos: Potential Energy of a Spring, Restoring Force of a Spring

AMPLITUDE AND PHASE: SECOND ORDER II (Mathlets)

Sistema MRA

Amplitude-and-Phase-2nd-Order-II

He realizado un resumen de los textos originales consultados para analizar las ecuaciones de los elementos que se consideran en este documento: masa, resorte, amortiguador.

Regresando a la Figura 5:

Acudimos a la Segunda Ley de Newton:

Esta ecuación nos dice que la sumatoria vectorial de todas las fuerzas que actúan sobre el cuerpo de masa m, es igual al producto del valor de dicha masa por su aceleración adquirida debido a dichas fuerzas. Recordando que la aceleración es la segunda derivada del desplazamiento y aplicando la Segunda Ley de Newton a nuestro sistema resorte-masa obtenemos la siguiente ecuación:

Arreglando un poco las cosas, obtenemos la ecuación que queríamos obtener desde un principio:

Esta ecuación representa La Dinámica de un Sistema Masa-Resorte ideal.

A parte de la figura 5, otra forma común de representar este sistema es mediante la configuración siguiente:

Fuente: Dinámica de Sistemas. Katsuhiro Ogata

En este caso debemos considerar la influencia del peso en la sumatoria de fuerzas que actúan sobre el cuerpo de masa m. El peso P está determinado por la ecuación P=m.g, donde g es el valor de la aceleración del cuerpo en caída libre.

Si se jala la masa hacia abajo y luego se suelta, actúa la fuerza de restitución del resorte, provocando una aceleración ÿ en el cuerpo de masa m. Obtenemos la siguiente relación aplicando Newton:

Si implícitamente consideramos la deflexión estática, es decir, si realizamos las medidas a partir del nivel de equilibrio de la masa colgando del resorte sin moverse, entonces podemos obviar y descartar la influencia del peso P en la ecuación. Si hacemos y=x, obtenemos de nuevo la ecuación:

Sistema Masa-Resorte-Amortiguador

 

Si no existiera ninguna fuerza de roce, el oscilador armónico simple oscilaría infinitamente. En la realidad, la amplitud de la oscilación disminuye gradualmente, un proceso conocido como amortiguación, descrito gráficamente a continuación:

Fuente: Física. Robert Resnick

El desplazamiento de un movimiento oscilatorio se grafica contra el tiempo, y su amplitud se representa mediante una función sinusoidal amortiguada por un factor exponencial decreciente que en la gráfica se manifiesta como una envolvente. La fuerza de fricción que actúa en el Movimiento Armónico Amortiguado es proporcional a la velocidad en la mayoría de los casos de interés científico. Dicha fuerza tiene la forma F = bV, donde b es una constante positiva que depende de las características del fluido que ocasiona la fricción, entre otras cosas. Esta fricción, también conocida como Fricción Viscosa, se representa mediante un diagrama que consiste en un pistón y un cilindro lleno de aceite:

La manera más popular de representar un sistema masa-resorte-amortiguador es mediante una conexión en serie como la siguiente:

Figura 6

Fuente: Física. Robert Resnick

 

Así como la siguiente:

Fuente: Dinámica de Sistemas. Katsuhiro Ogata

En ambos casos se obtiene el mismo resultado al aplicar nuestro método de análisis. Considerando la figura 6, podemos observar que es la misma configuración mostrada en a figura 5, pero agregando el efecto del amortiguador. Aplicando la segunda Ley de Newton a este nuevo sistema, obtenemos la siguiente relación:

Esta ecuación representa La Dinámica de un Sistema Masa-Resorte-Amortiguador.

Transformada de Laplace de un Sistema Masa-Resorte-Amortiguador

Una solución para la ecuación 37 se presenta a continuación:

Fuente: Física. Robert Resnick

La ecuación 38 muestra claramente lo que se había observado con anterioridad. Una ejemplo puede simularse en Matlab mediante el siguiente procedimiento:

Tcontinuo

La forma de la curva del desplazamiento en un sistema masa-resorte-amortiguador está representada por una sinusoide amortiguada por un factor exponencial decreciente. Esta solución puede derivarse mediante el método tradicional para resolver ecuaciones diferenciales. Sin embargo, dicho método es poco práctico cuando nos encontramos con sistemas más complicados como el siguiente:

Figura 7

Fuente: Control System Engineering. Norman Nise.

Surge entonces la propuesta de un método más práctico para hallar la dinámica de los sistemas y facilitar el posterior análisis de su comportamiento mediante simulación computarizada.

En la ecuación 37 no es fácil despejar x(t), que en ese caso es la función de interés. Tampoco puede representarse una ecuación diferencial en forma de Diagrama de Bloques que es el lenguaje más utilizado por los ingenieros para modelar sistemas, haciendo de lo complejo un objeto visual más fácil de entender y analizar. Esto conduce al primer objetivo para un método más práctico. El primer paso es separar claramente la función de salida c(t), la función de entrada r(t) y la función del sistema, alcanzando una representación como la siguiente:

La Transformada de Laplace permite alcanzar este objetivo de una manera rápida y rigurosa. Consiste en cambiar la función del dominio del tiempo al dominio de la frecuencia mediante la siguiente ecuación:

Fuente: Control System Engineering. Norman Nise.

La ventaja del método radica en que transforma derivadas en sumas y restas, luego, mediante asociaciones, podemos despejar la función de interés aplicando las simples reglas del algebra. Además, no es necesario aplicar la ecuación 2.1 a la función f(t) cuando se dispone de tablas que de antemano ya nos indican la transformada de funciones que se utilizan con gran frecuencia. En el caso de nuestros elementos básicos para un sistema mecánico, es decir: masa, resorte y amortiguador, contamos con la siguiente tabla:

Es decir, aplicamos un diagrama de fuerzas para cada unidad de masa del sistema, sustituimos la expresión de cada fuerza en tiempo por su equivalente en frecuencia (que en la tabla se denomina Impedancia, haciendo analogía entre sistemas mecánicos y sistemas eléctricos) y aplicamos la propiedad de superposición (cada movimiento se estudia por separado y luego se suma el resultado).

Considerando la figura 7, el resultado sería el siguiente:

La Transformada de Laplace llama a la función del sistema Función de Transferencia, cuya definición depende de cual es la función de entrada y cual la salida. Por ejemplo, para la figura 7 nos interesa conocer la Función de Transferencia G(s)=X2(s)/F(s).

Arreglando en forma matricial las ecuaciones del movimiento obtenemos lo siguiente:

Las ecuaciones 2.118a y 2.118b muestran un patrón que siempre se cumple y se puede aplicar para cualquier sistema masa-resorte-amortiguador:

La consecuencia inmediata del método anterior es que facilita enormemente obtener las ecuaciones del movimiento para un sistema masa-resorte-amortiguador, al contrario de lo que sucede con las ecuaciones diferenciales. Además, podemos llegar rápidamente a la solución exigida. En el caso de nuestro ejemplo:

donde

que son resultados que se obtienen aplicando las reglas del Algebra Lineal, lo que concede un gran poder computacional al método de Transformada de Laplace.

Caso Rotacional

Hasta ahora se ha considerado solamente el caso traslacional. En el caso de que el desplazamiento sea rotacional, la siguiente tabla resume la aplicación de la transformada de Laplace en ese caso:

Para ilustrar su uso consideramos el siguiente ejemplo:

Las siguientes figuras ilustran la manera cómo realizar el diagrama de fuerzas para este caso:

De esta manera, el resultado se obtiene a continuación:

Siendo:

Observamos que de nuevo se cumple que:

Representación en Variables de Estado de un Sistema Masa-Resorte-Amortiguador.

Las variables de estado son la herramienta más poderosa de la Ingeniería de Control Moderna, ya que no está limitada a sistemas lineales como sí o está el método hasta ahora visto, La Transformada de Laplace.

Las variables de estado en el caso del sistema masa-resorte-amortiguador de la Figura 8, nos permitirá reescribir un sistema de segundo orden en un sistema de primer orden. El siguiente material fue obtenido del video: State-Space Representation

                                                                      Figura 8

 

Seleccionando nuevamente el desplazamiento como la coordenada generalizada, la ecuación de movimiento del sistema es la siguiente:

El objetivo es expresar esta ecuación en una forma equivalente que tiene la siguiente forma:

Aquí el vector es un Vector de Estado, y X1, X2, son variables de estado que sustituyen a la original variable generalizada X y, más importante, a sus derivadas. El describir el sistema en forma de matrix, ofrecerá la enorme ventaja de utilizar el poder de las computadoras para procesar información y ejecutar análisis de datos presentados en forma matricial (Matrix Algebra).

Las ecuaciones encerradas en círculos amarillos muestran como la primera forma de escribir es la forma compacta de escribir las ecuaciones para y.

El primer paso es definir las variables de estado:

Este procedimiento nos permite obtener de inmediato la primera ecuación de estado :

…..por tanto

El segundo paso consiste en forzar al coeficiente que acompaña al orden más alto, el coeficiente líder, a ser igual a la unidad. Para ello, en nuestro caso, se divide la ecuación de movimiento original entre m (y en general, entre el valor que ocupe ese lugar):

En el tercer paso se despeja la derivada de mayor orden:

El cuarto paso consiste en sustituir las derivadas de la variable original por sus ya asignadas variables de estado:

Y así hemos encontrado la segunda ecuación de estado:

….

Y así hemos completado el objetivo. La ecuación de movimiento original puede ser expresado como variables de estado en la siguiente forma:

Dinámica de un Sistema Masa-Resorte-Amortiguador

Los libros de donde extraje imágenes, ecuaciones e información, que recomiendo para iniciar este estudio, son los siguientes:

  1. Robert Resnick, tomo1
  2. Dinamica_de_Sistemas, Katsuhiko Ogata
  3. Control Systems Engineering, Norman Nise
  4. Sistemas de Control Automatico, Benjamin Kuo
  5. Ingenieria de Control Moderna, 3° ED. – Katsuhiko Ogata

Escrito por Larry Francis Obando – TSU

Escuela de Ingeniería Eléctrica de la Universidad Central de Venezuela, CCs.

Escuela de Ingeniería Electrónica de la Universidad Simón Bolívar, Valle de Sartenejas.

Escuela de Turismo de la Universidad Simón Bolívar, Núcleo Litoral.

Contact:

WhatsApp: 00593984950376

email: dademuchconnection@gmail.com

Copywriting, Content Marketing, Tesis, Monografías, Paper Académicos, White Papers (Español – Inglés)