UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. 6. ESTABILIDAD DE SISTEMAS LINEALES Existen varias definiciones de estabilidad de acuerdo al problema que se quiera estudio de invariantes analizar. sistemas en el Para el lineales tiempo generalmente se suele escuchar los términos de Neutralmente estable, Asintóticamente Asintóticamente estable estable (local), (globa) e inestable. En general el estudio de la estabilidad de un sistema tiene que ver con que el comportamiento del sistema (respuesta) siempre permanezca dentro de límites aceptables. Es decir que “exista control”. En general el análisis de la estabilidad depende del sistema que se esté tratando de analizar. Es así que ciertos resultados son únicamente válidos para sistemas lineales, mientras que otros criterios más sofisticados son aplicables a cualquier tipo de sistema. El estudio de la estabilidad de sistemas lineales e invariantes en el tiempo se la puede realizar tanto en el dominio del tiempo como en el dominio de la frecuencia. Las siguientes secciones estudian las distintas definiciones y herramientas necesarias para el estudio de la estabilidad de estos sistemas. 66 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. 6.1 Herramientas para el Análisis de Estabilidad Para estabilidad es suficiente con determinar la localización de los polos o de las raíces del sistema en lazo cerrado. En la actualidad la determinación de los valores de las raíces es bastante sencilla debido a la gran cantidad de métodos numéricos disponibles en computadoras y calculadores. Sin embargo, tiempo atrás esto era más complicado por lo que se recurría al uso de resultados algebraicos y métodos semigráficos que permiten determinar la localización de los polos del sistema sin necesidad de encontrar sus valores. En general se destacan cuatro métodos de análisis, 1. Criterio de Routh-Hurwitz: Es un método algebráico que proporciona información sobre la localización de las raíces de un polinomio. 2. Lugar Geométrico de las Raíces (LGR): Es un método semigráfico el cual permite determinar como afecta el cambio de una constante en la localización de las raíces de un polinomio. 3. Criterio de Nyquist: Es un método semigráfico que provee información sobre el número de polos y ceros de la función de transferencia en lazo cerrado que están en el semiplano derecho 67 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. del plano s mediante la observación del comportamiento de la traza de Nyquist. 4. Criterio de Bode: El criterio de Bode es en general un criterio integral que define la estabilidad de un sistema a través del valor de esta. Se suele utilizar la traza de Bode para cuantificar la estabilidad de un sistema (Estabilidad Relativa) de manera similar al uso de la traza de Nyquist. El uso de todos estos métodos y criterios no se detallaran todos en esta sección, solamente recordaremos y analizaremos el Criterio de RouthHurwitz. Criterio de Routh-Hurwitz Edward John Routh (1831-1907) ganó en 1877 el premio Adams por su método de analizar la estabilidad de sistemas dinámicos de cualquier orden, siguiendo estudios de James Clerk Maxwell. El matemático alemán Adolf Hurwitz (18591919) propuso su método para asegurar la estabilidad de turbinas de vapor en 1895, siguiendo estudios de Aurel Boreslav Stodola. 68 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. En 1911 se demostró que ambos criterios de estabilidad eran equivalentes. El método de Routh-Hurwitz está formado por dos componentes: El criterio de John Hurwitz y la tabulación de Routh. El criterio de Hurwitz establece que dado un polinomio de coeficientes reales como el que se presenta a continuación, P(s) s n an1s n1 ... a1s ao para que P(s) no tenga raíces con partes reales positivas es necesario que se cumplan las siguientes dos condiciones: Todos los coeficientes tengan el mismo signo Ninguno de los coeficientes sea igual a cero Si las dos condiciones anteriores se cumplen no necesariamente las raíces estarán en el semiplano izquierdo. Es decir si una de las dos condiciones anteriores no se cumple, entonces necesariamente existirán raíces en el semiplano derecho, pero si se cumplen puede como no puede que existan raíces en el semiplano derecho. La condición de suficiencia se la obtiene utilizando la tabulación de Routh, que establece que: a6 s 6 a5 s 5 a4 s 4 a3 s 3 a2 s 2 a1s1 a0 0 69 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Para que no existan raíces en el semiplano derecho todos los coeficientes de la primera columna de la tabulación de Routh deben ser del mismo signo. El número de cambios de signo en los elementos de la primera columna es igual al número de raíces con partes reales positivas. La tabulación de Routh se la hace de la siguiente manera para, por ejemplo, un polinomio de sexto orden como el siguiente: Ejemplo: Realizar el análisis de la localización de las raíces de los siguientes polinomios: 1. P(s) s 3 4s 2 s 6 Como existe un coeficiente con signo distinto a los demás se sabe que existirá al menos una raíz en el semiplano derecho, por lo que el sistema es inestable. Se procede entonces a la tabulació de Routh, 70 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Como se ve existen 2 cambios de signo, por lo que hay dos raíces en el semiplano derecho. 2. P(s) 2s 4 s 3 3s 2 5s 10 Todos los coeficientes tienen el mismo signo, por lo que es necesario realizar la tabulación de Routh. Se observan dos cambios de signo en la primera columna, por lo que hay dos raíces en el semiplano derecho. 3. P(s) s 4 s 3 s 2 s K Todos los coeficientes son mayores que cero y debemos asegura que también lo sea. Si entonces se cumplen las condiciones necesarias, pero es necesario realizar la tabulación para ver que rangos debe tener 71 . UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. En la tercera fila se ha reemplazado el valor de 0 por un valor muy pequeño lógicamente mayor que cero, esto se hace para evitar divisiones para cero. Como se ve en el resultado para cualquier valor de K, el sistema siempre será inestable dado que la expresión K2 siempre será negativa. 4. P(s) s 5 2s 4 2s 3 4s 2 11s 10 Todos los coeficientes tienen el mismo signo, por lo que es necesario realizar la tabulación de Routh. Existen 2 cambios de signo, por lo que el sistema es inestable. EJERCICIOS 1. Analice la estabilidad en lazo abierto y lazo cerrado para una planta dada por T ( s) 72 1 s 3s 1 2 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. 2. Sea un proceso físico cuya ecuación característica es: s 4 K s3 s 2 s 1 0 Determinar por el Criterio de Routh el rango de valores de K para que el sistema sea estable. 3. Sea un proceso físico cuya ecuación característica es: s 3 3408.3 s 2 1204000 s 1.5 107 K 0 Determinar por el Criterio de Routh el rango de valores de K para que el sistema sea estable. 4. Sea un proceso físico cuya ecuación característica es: s 3 3 K s 2 (k 2) s 4 0 Determinar por el Criterio de Routh el rango de valores de K para que el sistema sea estable. 5. Sea un proceso físico cuya ecuación característica es: s5 4 s 4 8 s3 8 s 2 7 s 4 0 Determinar por el Criterio de Routh si el sistema es estable. 73 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. 7. DISEÑO DE UN SISTEMAS DE CONTROL USANDO REALIMENTACIÓN DE VARIABLES DE ESTADO 7.1. DISEÑO DE SERVOSISTEMAS Un servosistema es un tipo de sistema de control en el que, en estado de funcionamiento normal, la variable de entrada está variando con el tiempo. A dicha variable se le suele denominar señal de referencia. Dentro de los servosistemas encontramos un caso particular, este es el servomecanismo: servosistema en el que la variable controlada es una posición mecánica. Ahora vamos a analizar el método de asignación de polos para el diseño de servosistemas de tipo 1. Estudiaremos servosistemas de tipo 1 para plantas con un integrador y sin integradores. Diremos que una planta con un integrador es aquella en la que su ecuación característica posee un factor . Por ejemplo o también se limitará a sistemas que tengan una señal de control . El método escalar (tipo escalón) y una salida y también escalar. Supóngase que la planta con integrador se define mediante: x Ax Bu y Cx Du 74 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. donde x es el vector de estado para la planta (vector de dimensión x ) u =señal de control (escalar) y =señal de salida (escalar) A=matriz de coeficientes constantes B=matriz de coeficientes constantes C=matriz de coeficientes constantes Mediante una elección adecuada es posible seleccionar la salida igual a una de las variables de estado. La siguiente figura muestra una configuración general de este servosistema, donde y x1 y se supone que la entrada de referencia r es una función escalón. En este sistema se utiliza el siguiente esquema de control mediante realimentación del estado. En donde u Kx k1r y K [k1 , k2 , k3 ,..., kn ] 75 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Si suponemos que la entrada de referencia se aplica en . Asi para , la dinámica del sistema se puede describir mediante: x Ax Bu ( A BK ) x Bk1r Se diseñará el servosistema de modo que los polos en lazo cerrado se localicen en las posiciones asintóticamente estable, deseadas. y() El sistema diseñado tenderá al valor constante será r y u () tenderá a cero. En el estado estacionario tendremos: x() ( Ax Bu) x() Bk1r () Se tiene que r () r (t ) r para . Restando ambas ecuaciones se tiene que: x(t ) x() ( A BK )[ x(t ) x()] definimos x(t ) x() e(t ) , y la ecuación se convierte en: e ( A BK )e . Esta ecuación describe la dinámica de error. Ahora el diseño del servosistema de tipo 1 se convierte aquí en el diseño de un sistema regulador asintóticamente estable tal que tiende a cero, para cualquier condición inicial Luego si el sistema definido al principio del todo es completamente controlable, entonces sólo tenemos que especificar los valores propios deseados para la matriz . La matriz hallar mediante Matlab como le explicare más adelante. 76 se puede UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Ejemplo 1: Diseñar un servosistema de tipo 1 cuando la función de transferencia de la planta tiene un integrador. La función de transferencia de la planta se obtiene mediante: T ( s) Observe un factor Y ( s) 1 1 3 U ( s) s( s 1)(s 2) s 3s 2 2s en el denominador, esto es un integrador. Consideremos los polos deseados en lazo cerrado en la posición s 2 2 j 3 , s 2 2 j 3 . Se supone que la entrada de referencia es la función escalón. Se pide obtener también la respuesta a un escalón unitario del sistema diseñado. Analizaremos su estabilidad usando el criterio de Ruth y por simulación. Como hay un coeficiente nulo, según el criterio de Ruth decimos que el sistema es inestable. Analizando el sistema mediante Simulink de Matlab tenemos: Planta con un integrador 77 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Respuesta de la planta ante una entrada tipo escalón unitario Realizaremos el diseño utilizando Matlab. Determinaremos primero las matrices A, B, C y D de las variables de estado y analizaremos su cotrolabilidad. >>num=[1] >>den=[1 3 2 0] >>[A, B, C, D]=tf2ss(num,den) >> M=ctrb(A,B) >>rank(M) Con esto notamos que el rango de la matriz de controlabilidad es 3 por lo tanto el sistema es controlable. Encontrando las ecuaciones de estado, construimos el diagrama de bloques y observamos también que la respuesta es completamente inestable. 78 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Planta en variables de estado Respuesta de la planta ante una entrada tipo escalón unitario Encontremos ahora la matriz K de controlabilidad: >>J=[ 2 2 j 3 22j 3 10 ] >> K=acker(A,B,J) En el vector J se ha añadido un polo en -10 para completar el orden, esto polo no afectará puesto que se ha seleccionado por lo menos unas 5 veces el dominante. Esta es la matriz K de realimentación del estado. 79 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. El diagrama de bloques con la realimentación y la respuesta ante una entrada tipo escalón es: Servosistema tipo 1. Planta con un integrador Respuesta ante entrada tipo escalón unitario Ejemplo 2: Diseñemos un servosistema de tipo 1 cuando la planta no tiene integrador. 80 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Supongamos que se tiene el sistema dado por: Deseamos unos polos en s 1 j 3 y en s 1 j 3 . El estudiante puede verificar que se trata de un sistema sin integradores encontrando su función de transferencia. Se trata primero de generamos una matriz ampliada de controlabilidad de la forma Y encontramos su rango. Si el rango es n+1, diremos que es completamente controlable. >>P=[A B; -C 0] >>rank(P) Como puede notar es completamente controlable. Ahora encontremos la matriz K. 81 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. Debemos ahora ampliar las matrices A y B. Esto se requiere para que queden de orden y respectivamente y puedan ser operadas: y >>A1=[A zeros(4,1); -C 0] >>B1=[B; 0] Asignemos los polos, por ejemplo: >>J=[-1+j*sqrt(3) -1-j*sqrt(3) -5 -5 -5] Note usted se añaden otros 3 polos alejados por lo menos 5 veces de la parte real de los polos deseados importantes para que se pueda operar y no tengan efecto. Escribimos ahora en Matlab: >> K1=acker(A1, B1, J) En este caso la matriz K da un valor más. Aquí debemos tomar Nótese que se debe tomar con signo contrario. 82 UNIVERSIDAD SANTIAGO DE CALI FACULTAD INGENIERIAS PROGRAMA DE INGENIERIA ELECTRONICA CONTROLES II Profesor: MSc. JAVIER A. MURILLO M. El sistema en variables de estado y con los valores de realimentación se muestra en la figura siguiente. La realimentación debe ser negativa pero como cada valor de es negativo, se le ha cambiado el signo al sumador. El valor se puede también introducir al bloque con su signo: Servosistema tipo 1. Planta sin integrador La respuesta ante un escalón unitario es de la forma: Respuesta ante entrada tipo escalón unitario 83
© Copyright 2024