Control predictivo por modo deslizante para robots manipuladores
←
→
Transcripción del contenido de la página
Si su navegador no muestra la página correctamente, lea el contenido de la página a continuación
REVISTA INGENIERÍA UC. Vol. 11, No 2, 39-47, 2004 Control predictivo por modo deslizante para robots manipuladores Winston García-Gabín, Darine Zambrano Postgrado en Automatización e Instrumentación, Escuela de Ingeniería Eléctrica, Facultad de Ingeniería. Universidad de los Andes, Mérida-Venezuela Email: winstong@ula.ve, darine@ula.ve Resumen En este artículo se presenta el desarrollo de un nuevo controlador predictivo basado en modelo. El controlador predictivo por modo deslizante propuesto combina la técnica de control predictivo con el control por modo deslizante. Éste esquema presenta las ventajas del control predictivo en cuanto a desempeño y las de robus- tez debidas al control por modo deslizante. La robustez del controlador propuesto es evaluada por medio del gráfi- co de robustez y es comparada con la de un controlador predictivo clásico. El desempeño es mostrado por la simulación de un modelo no lineal de un brazo manipulador con fricción, y ha sido comparado con el control pre- dictivo y el control por modo deslizante. Palabras clave: Control predictivo, control por modos deslizantes, robot manipulador, fricción. Sliding mode predictive control for manipulator robots Abstract This article shows the synthesis of a new model based predictive controller. The proposed sliding mode predictive controller combines predictive control with sliding mode control. The sliding mode predictive control shows has the good performance of model predictive control model and the robustness of sliding mode control. The robustness is evaluated using a robustness graphic, and compared a classic predictive control. The performance is shown by simulation of a nonlinear manipulator robot model with friction, and compared with the model predictive control and the sliding mode control. Keywords: Predictive control, sliding mode control, manipulator robot, friction 1. INTRODUCCIÓN El Control Predictivo (Model Predictive Con- las ventajas que presenta el control predictivo basado trol, MPC) se desarrolló a finales de los setenta y ha en modelos se puede destacar [1]: ajuste intuitivo y tenido un desarrollo considerable desde entonces. El simple, puede ser utilizado para el control de sistemas término Control Predictivo no designa a una estrategia con dinámicas complejas, el caso multivariable puede de control particular sino a un conjunto de métodos de tratarse fácilmente, compensación interna del retardo control que hacen uso explícito de un modelo del pro- del sistema, permite la compensación de las perturba- ceso para obtener la señal de control minimizando una ciones medibles, la incorporación de restricciones es función objetivo. Estos métodos de control conllevan conceptualmente simple durante el proceso de dise- a controladores que tienen básicamente la misma es- ño del controlador, permite la utilización de las refe- tructura y los mismos elementos [1]. El MPC es una rencias futuras cuando son conocidas, que es de espe- técnica bastante popular tanto en el ámbito académico cial utilidad en el área de robótica. Es una metodolo- como industrial, donde ha sido implementado en gía abierta basada en ciertos principios que permite numerosas aplicaciones mostrando buen desempeño; futuras extensiones. Los distintos algoritmos de con- se pueden encontrar aplicaciones en la industria de trol predictivo difieren en el tipo de modelo utilizado cemento, desecadoras, brazos robóticos [2, 3]. Entre para representar al proceso y a las perturbaciones, y la Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004 39
Control predictivo por modo deslizante para robots manipuladores función objetivo considerada. Una vez que el sistema es atrapado en la su- perficie, donde la estructura y los parámetros han sido La idea básica de un MPC es calcular una establecidos por el diseñador, la dinámica a lazo ce- secuencia de futuras señales de control de tal manera rrado viene determinada completamente por las ecua- que minimice una función de costo sobre un horizon- ciones que la definen, y es independiente de las per- te de predicción. El índice a minimizar es el valor turbaciones en los parámetros del sistema, de esta esperado de una función del error cuadrático entre la manera se logra una excelente robustez. salida futura de la planta y la secuencia de valores de la referencia futura, más una función cuadrática que En la Figura 1, se presenta la interpretación pondera el esfuerzo de control. Con el objetivo de gráfica del control por modo deslizante, partiendo de implementar el MPC, se utiliza un modelo del sistema un punto inicial la ley de control lo conduce a la su- para predecir las salidas futuras de la planta, sobre la perficie de deslizamiento (modo de alcanzabilidad), base de los valores pasados de las entradas y salidas una vez allí se desliza atrapado por la superficie hasta del sistema. Es por ello que el modelo del sistema el valor deseado (modo deslizamiento). constituye un elemento fundamental en el desempeño del controlador, y es conveniente seleccionar un mo- delo lo más representativo posible del sistema real con el objeto de lograr una estimación exacta de las salidas futuras del sistema. La formulación convencional del MPC em- plea un modelo lineal para realizar la estimación de las salidas futuras de la planta, es por ello que las no linealidades constituyen un serio obstáculo para el buen rendimiento del controlador, puesto que el mo- delo lineal deja de producir estimaciones exactas de la dinámica futura del sistema al alejarse del punto de operación donde se obtuvo el modelo lineal. El Control por Modo Deslizante (Sliding Mo- Figura 1. Superficie de deslizamiento. de Control (SMC)) es un caso especial de Control por Estructura Variable (Variable Structure Control En el control de brazos manipuladores dos de Control (VSC)), desarrollado en la Unión Soviética las dificultades más relevantes son las no linealidades en la década de los 60, [4]. Las técnicas de SMC han y la fricción [5]. La fricción aparece por la interfase sido empleadas en diversos sistemas, por las ventajas física entre dos superficies en contacto. En el control desde el punto de vista de robustez que presenta en su de brazos manipuladores los errores en estado esta- aplicación para el control de sistemas no lineales, va- cionario y errores en el seguimiento de trayectoria riantes en el tiempo y sistemas con incertidumbre. son causados principalmente por la fricción entre los eslabones del brazo [6]. La fricción es comúnmente La metodología de diseño de un controlador modelada como una combinación de fricción Cou- por modos deslizantes [4], implica dividir el pro- lomb donde el par producido es de valor constante blema en dos subproblemas de menor dimensión: Pri- cambiando el signo en función del signo de la veloci- mero establecer la dinámica de la superficie de desli- dad y fricción viscosa donde el par producido es pro- zamiento, que debe de ser estable y puede ser lineal, porcional a la velocidad del eslabón en movimiento diseñando entonces, la superficie que es la encargada [5]. El modelado de la fricción ha sido estudiado por de regir el comportamiento del sistema cuando se en- varios autores entre otros [7, 8]. La fricción es alta- cuentra en el modo de deslizamiento. Luego se diseña mente no lineal y produce ciclos límites y mal desem- una ley de control que garantice que las trayectorias peño en los controladores. Al aplicar controladores cerca de la superficie de deslizamiento conduzcan PD a sistemas con fricción se presentan errores de hacia la superficie. estado estacionario, este se puede mejorar al aumentar la ganancia pero esto tiende a producir 40 Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004
García-Gabín y Zambrano inestabilidad en el sistema a lazo cerrado. Con la ac- forma más apropiada por medio del modelo ción integral, el PID puede mejorar el error en estado CARIMA. estacionario pero produce un ciclo límite [9]. Para el desarrollo del controlador es considera- En vista de las ventajas destacadas en ambos da la siguiente superficie deslizante métodos, resulta atractivo realizar una combinación entre la técnica de control predictivo y el control por S t + j = P ( q − 1 ) ( y t + j − w t + j ) + Q ( q − 1 ) ∆ u t + j −1 (2) modo deslizante, de manera de desarrollar un contro- lador con las ventajas en cuanto a desempeño dadas El algoritmo de control predictivo consiste en por el controlador predictivo y la robustez del control obtener una secuencia que minimiza la función de por modo deslizante, para ser aplicado a brazos mani- costo (3) [10]. Esta minimización produce ∆ut, ∆ut+1, puladores con fricción. Esta combinación ha presenta- ..., ∆ut+Nu, pero solamente ∆ut es aplicado como cam- do buenos resultados para sistemas de fase no mínima bio en la señal de control. En el instante de tiempo [10]. Ha sido aplicado a modelos de procesos quími- t+1 es calculada una nueva secuencia, esta metodolo- cos con características no lineales y de fase no gía es llamada horizonte deslizante. mínima [11]. ∑ [S ] + ∑ λ [∆u ] N2 Nu J (u, t ) = 2 2 El artículo se organiza de la siguiente manera: t+ j t + j −1 (3) j = N1 j =1 en la sección 2 se presenta el desarrollo del controla- dor predictivo por modo deslizante (Sliding Mode Predictive Control , SMPC). En la sección 3 se pre- Donde, senta la prueba de robustez para el controlador pro- N1 : Inicio del horizonte de predicción puesto. En la sección 4 se muestran las pruebas de N2 : Final del horizonte de predicción desempeño de los controladores por medio de la simu- Nu : Horizonte de control lación de un modelo no lineal de un brazo manipula- λ : Factor de ponderación del esfuerzo de control dor con fricción Coulomb y viscosa. Finalmente son presentadas las conclusiones. Con el objetivo de resolver el problema de la minimización de (3), yt+j se calcula como la j-ésima salida predicha para j= N1, ..., N2 basado en la infor- 2. CONTROL PREDICTIVO POR MODO mación conocida en el instante t y los valores futuros DESLIZANTE de los incrementos de control, la siguiente ecuación diofántica es considerada A la mayoría de los sistemas es posible aproxi- marlos, luego de una linealización sobre un punto de C ( q −1 ) = E j ( q −1 ) A ( q −1 ) ∆ + q −1 F ( q −1 ) (4) operación, por medio de: ( ) ( ) A q −1 y t = B q −1 u t −1 + ( ) C q −1 ξt (1) Combinando el modelo de la planta (1), y la ∆ ecuación diofántica (4), es obtenida la siguiente ecua- Donde: ción de predicción, yt : Variable controlada Fj EjB u t : Variable manipulada yˆ t + j = yt + u t + j −1 (5) ξ t : Ruido blanco de media cero C C A(q-1): Polinomio moníco En esta expresión, yˆ t + j es una función de los B(q-1): Polinomio que contiene los ceros del modelo valores conocidos en el instante de tiempo t y de las C(q-1): Polinomio moníco futuras señales de entrada de control que aún no han D: 1-q-1 sido calculadas. Empleando una segunda ecuación diofántica para distinguir el pasado y el futuro de las Este modelo es conocido como modelo CARI- señales de control se tiene, MA (Controller Auto-Regressive Integrated Moving- Average); varios autores [1] sostienen que diversas E j ( q − 1 ) B ( q − 1 ) = G j ( q − 1 ) C ( q − 1 ) + q − j Γ j ( q − 1 ) (6) aplicaciones industriales pueden ser representadas de Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004 41
Control predictivo por modo deslizante para robots manipuladores La siguiente expresión de la predicción de la ⎡ g0 0 L 0 ⎤ salida es obtenida, ⎢ g ⎢ 1 g0 L 0 ⎥⎥ yˆ t + j = G j ∆u t + j −1 + Γ j u tf−1 + F j y tf (7) ⎢ M M O M ⎥ (17) G=⎢ ⎥ Donde ∆u tf y ytf son los valores filtrados ⎢ g N u −1 g Nu −2 L g0 ⎥ de y, ∆u t ⎢ M M M ⎥ ⎢ ⎥ ( ) u tf = C −1 q −1 ∆u t (8) ⎢⎣ g N 2 −1 g N2 −2 L g N 2 − N u ⎥⎦ y tf = C −1 (q )y −1 t (9) Sustituyendo (16) en (2) se tiene Finalmente, (7) puede ser escrita como, N2 ∑ ⎡⎣ P ( q ) ( Gu + f − w ) + Q ( q 2 J (u, t ) = )u ⎤⎦ yˆ t + j = G j (q )∆u −1 −1 −1 (10) t + j −1 + yˆ t + j t j = N1 Donde yˆ t + j t es la respuesta libre predicha de yˆ t + j Nu 2 (18) + ∑ λ ⎡⎣ ∆ut + j −1 ⎤⎦ asumiendo que los futuros incrementos de control lue- j =1 go del tiempo t-1 serán cero, ( ) yˆ t + j t = Γ j q −1 u tf−1 + F j q −1 ytf ( ) (11) Reagrupando los términos de la respuesta libre N2 Sustituyendo Ej de (4) en (6), produce J (u, t ) = ∑ [ PGu + P ( f + P j = N1 * y + Q *u ) − B ( A∆ ) = G j + q − j Γ j C −1 + q − j BF j ( A∆C ) (12) −1 −1 Nu Pw + Q ( q −1 )u ]2 + ∑ λ ⎡⎣ ∆ ut + j −1 ⎤⎦ 2 (19) j =1 Definiendo el vector f, como las predicciones de la respuesta libre, Donde: ⎡ p0 0 L 0 ⎤ f = [yˆ t +1 t , yˆ t + 2 t ,..., yˆ t + N 2 t ] T (13) ⎢ p P=⎢ 1 p0 L 0 ⎥ ⎥ (20) ⎢M M O M ⎥ el vector de los futuros incrementos de control, ⎢ ⎥ ⎢⎣ 0 p Nc L p0 ⎥⎦ [ u~ = ∆u t , ∆u t +1 , L, ∆u t + N u −1 ] T (14) ⎡ p1 L p Nc −1 p Nc 0⎤ 0 P = ⎢⎢ p2 L 0 ⎥⎥ y el vector de salidas predichas del sistema, * (21) L p Nc L [ yˆ = yˆ t +1 , yˆ t + 2 , L, yˆ t + N 2 ] T (15) ⎢⎣ 0 0 0 0 0 0 ⎥⎦ A partir de la predicción (11) la relación ⎡ q0 0 L 0 ⎤ entrada-salida de la planta puede ser escrita como la ⎢q q0 L 0 ⎥ ecuación vectorial, Q=⎢ 1 ⎥ (22) ⎢M M O M ⎥ yˆ = Gu~ + f (16) ⎢ ⎥ ⎢⎣ 0 q Nc L q0 ⎥⎦ Donde la matriz G esta compuesta por los ⎡ q1 L q Nc −1 q Nc 00⎤ parámetros de la respuesta escalón gi del modelo de la Q = ⎢⎢ q2 * L q Nc L L 0 ⎥⎥ (23) planta. ⎢⎣ 0 0 0 0 0 0 ⎥⎦ 42 Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004
García-Gabín y Zambrano ambos parámetros se encuentra la frontera de la re- ⎡λ 0 L 0 ⎤ gión de estabilidad. ⎢0 λ L 0 ⎥ λ=⎢ ⎥ (24) En la gráfica mostrada en la Figura 3 se en- ⎢M M O M ⎥ ⎢ ⎥ cuentra dibujada una ventana delimitada por los pun- ⎣ 0 0 L λ⎦ tos (2, 0), que se corresponde con un aumento del do- ble de la ganancia, (0.5, 0) que se corresponde con Cuando no hay restricciones en la señal de una disminución a la mitad de la ganancia, y (0, 0.5), control, la minimización de la ecuación cuadrática (2) (0, 2) para variaciones de la relación de la frecuencia puede ser resuelta como un problema de álgebra li- del polo. Todas las combinaciones de relación de ga- neal: nancia y relación de constante de tiempo con un factor de dos son acotadas dentro de la ventana de robustez. ∆u = K ( w − f ) (25) Si el límite de estabilidad del controlador se encuentra −1 fuera de la ventana, el controlador se considera robus- K = ⎡ ( PG + Q ) ( PG + Q ) + λ ⎤ ( PG + Q ) T T P (26) to. ⎣ ⎦ 3. PRUEBA DE ROBUSTEZ En está sección se compara la robustez del SMPC contra un MPC clásico. La robustez es evalua- da utilizando el gráfico de robustez [12]. En el gráfico de robustez el límite de estabilidad para el sistema en lazo cerrado es dibujado en coordenadas logarítmicas cuando la constante de tiempo del proceso y la ganan- cia cambian. Los ejes de la gráfica son la relación de ganancia (ganancia del proceso dividido entre la ga- nancia con la cual fue sintonizado el controlador ini- cialmente) con respecto a la relación de la frecuencia del polo del sistema (frecuencia del polo del proceso Figura 2. Ajuste inicial de los controladores. dividido entre frecuencia del polo del modelo con el cual fue sintonizado el controlador inicialmente). La gráfica de robustez se construyó ajustando inicialmente los controladores para el modelo de la La construcción del gráfico de robustez de ecuación (27) de manera tal que tuviesen el mismo esta forma, constituye una modificación a la propuesta sobrepaso como se observa en la Figura 2. Se observa de Shinskey [12] en la cual se representa la relación como con el mismo sobrepasó el tiempo de asenta- de ganancia versus relación retardo. miento del SMPC es menor que el del MPC. 1 El gráfico tiene una región de estabilidad y Gs = (27) otra de inestabilidad. Para encontrar el límite de la s ( s + 1) región de estabilidad, el controlador es sintonizado inicialmente en base a un grupo de parámetros del Utilizando el gráfico de la Figura 3 se obtie- modelo del proceso representados en el punto central nen los siguientes valores del Índice de Robustez (IR) (1,1) de las coordenadas. La ganancia del sistema es (MPC = 0,9 y SMPC = 2,7). Los resultados de esta aumentada hasta que el sistema en lazo cerrado pre- sección muestran como el SMPC propuesto presenta sente una oscilación sostenida al margen de la inesta- una mayor robustez que el MPC convencional, tanto bilidad, ese valor de la relación entre la ganancia para cambios en la ganancia como para cambios en la respecto a la ganancia inicial origina un punto en la frecuencia del polo del modelo. gráfica. De igual forma la constante de tiempo del modelo del sistema es cambiada respecto al valor ini- cial. Realizando en forma simultanea variaciones de Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004 43
Control predictivo por modo deslizante para robots manipuladores La Figura 5, presenta la señal de control para los cambios de posición, de la Figura 4. Posición (Grados) 100 Referencia 80 Posición 60 40 20 Figura 3. Gráfico de robustez. 0 0 1 2 3 4 5 6 7 4. EJEMPLO DE APLICACIÓN Tiempo (s) Con el objetivo de estudiar el desempeño del Figura 4. Trayectoria de posición del brazo manipulador. controlador, se utilizó el modelo no lineal de brazo manipulador con fricción Coulomb y fricción viscosa, Señal de Control 6 la ecuación que rige su dinámica viene determina por: 4 τ = m l 2θ&& + vθ + c sgn (θ& ) + m l g Cos (θ ) 2 (28) 0 -2 Donde, -4 τ : Par aplicado al sistema -6 θ : Posición del brazo m: Masa del sistema (3 kg) -8 0 1 2 3 4 5 6 7 l: Longitud del brazo (0,3 m) Tiempo (s) v: Coeficiente de fricción (0,6 kg m2/s) Figura 5. Señal de control para la trayectoria de posición. c: Coeficiente de fricción viscosa (0,3 kg m2/s) g: Aceleración de gravedad (9,81 m/s2) 4.1 Pruebas de desempeño Sustituyendo los valores se tiene: En esta sección se muestra un estudio del τ = 0,27θ&& + 0,60θ& + 0,30 sgn (θ& ) + 8.83Cos(θ ) comportamiento del sistema a lazo cerrado, para el (29) que se lleva el modelo dinámico no lineal a puntos de operación donde se observa un comportamiento bas- donde se observa el carácter no lineal del sistema. tante oscilatorio en los controladores, siendo el objeti- vo principal contrastar el comportamiento de los con- En la Figura 4, se observa el comportamiento troladores ante cambios en el sistema una vez que los del brazo manipulador aplicando cambios de posición controladores han sido ajustados. Se presentan com- entre 0° y 90°, para diferentes velocidades de cambio paraciones entre el controlador predictivo por modo de posición. deslizante propuesto contra un MPC, en una formula- ción convencional [2, 3], y contra un SMC que puede Se puede observar como el controlador mantie- ser aplicado a sistemas integrantes[13, 14]. ne sin oscilaciones la posición del brazo manipulador, en la referencia para diferentes posiciones y velocida- En la Figura 6 se muestra el desempeño de des de cambio de posición, en todo el rango de opera- SMPC vs. MPC ajustados con los mismos parámetros ción del sistema. de sintonización, para un cambio de referencia en la posición del brazo manipulador, se observa como el 44 Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004
García-Gabín y Zambrano SMPC presenta un menor tiempo de respuesta y un cambio la masa del sistema de 3 kg a 9 kg, dejando mayor amortiguamiento del efecto producido por la los valores de sintonización del controlador con el fricción en el brazo manipulador. valor de ajuste inicial. Posición (Grados) 160 Posición (Grados) 70 GPC MPC SMPC 160 140 SMPCGPC+SM SMC 60 140 120 SMC 50 MPC 120 Po 100 100 sic ió 80 40 n (G 80ra 60 30 60 do s) 40 20 40 20 10 20 0 0 0 -20 0 0.5 1 1.5 2 2.5 3 3.5 4 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 0 0,5 1,0 1,5 2,0 2,5 Tiempo (s) 3,0 3,5 4,0 Tiempo (s) Tiempo (s) Figura 6. Cambio de referencia SMPC vs. MPC Figura 8. Cambio de referencia SMPC vs. MPC, SMC (igual ajuste). Observando en la Figura 8, la posición del bra- Ajustando con los mismos parámetros el SMPC zo manipulador se encuentra que el MPC es inestable y SMC se muestra el desempeño para un cambio de para el cambio de masa realizado (3 kg a 9 kg ). El referencia en la posición del brazo manipulador. En la SMPC y SMC logran estabilizar la posición en su Figura 7 se observa cómo el control predictivo por valor de referencia, pero el SMPC presenta un com- modo deslizante presenta un menor tiempo de res- portamiento menos oscilatorio y con un menor tiempo puesta, menor sobreoscilación, y similar amortigua- de estabilización. miento del efecto producido por la fricción en el brazo Señal de Control manipulador en estado estacionario. 800 Posición (Grados) 600 SMPC MPC 80 400 70 SMC SMPC 200 60 0 50 -200 40 -400 30 -600 20 -800 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 10 Tiempo (s) 0 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 Figura 9. Señal de control SMPC vs. MPC Tiempo (s) La Figura 9 muestra la señal de control SMPC Figura 7. Cambio de referencia SMPC vs. SMC vs. MPC, se ve como el MPC presenta un comporta- miento inestable y el SMPC estabiliza el sistema. Para realizar una mejor comparación entre el desempeño de los controladores, se realizaron diferen- En la Figura 10 se observa cómo la señal de tes ajustes en los mismo, de manera de lograr en con- control producida por el SMC, presenta oscilación de diciones nominales la misma dinámica en la posición alta frecuencia y una mayor amplitud en comparación del brazo manipulador para un cambio de referencia. con la señal de control producida por el SMPC. Luego con el objetivo de estudiar la robustez del con- trolador ante cambios en el brazo manipulador se Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004 45
Control predictivo por modo deslizante para robots manipuladores Señal de Control ñal de control respecto al SMPC. 500 800 Señal de Control 300 600 SMPC MPC 100 400 0 200 -100 0 -300 -200 -400 -500 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 -600 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 Tiempo (s) Tiempo (s) Figura 10. Señal de control SMC Figura 12. Señal de control SMPC vs. MPC Como segunda prueba para estudiar la robustez Señal de Control del controlador, se cambiaron los coeficientes de fric- 600 ción del brazo manipulador aumentándolo en un 40 %. La Figura 11 muestra la posición del brazo ma- 400 nipulador para estas nuevas condiciones, con los con- 200 troladores ajustados con los valores nominales. 0 Posición (Grados) -200 160 GPC 140 SMC GPC-SM MPC -400 160 SMC 140 120 -600 120 Po 100 0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 sic ión 80 Tiempo (s) 100 (G ra 80do 60 s) Figura 13. Señal de control SMC 60 40 40 20 20 5. CONCLUSIONES 0 SMPC 0 -20 0 0 0,5 1,0 1,5 0.5 1 1.5 2 2,0 2,5 Tiempo (s) 2.5 3 3.5 3,0 3,5 4 4,0 Se presentó el desarrollo de un controlador ba- sado en la combinación de las técnicas de control pre- Tiempo (s) dictivo y control por modo deslizante. El controlador propuesto presenta la robustez característica del con- Figura 11. Cambio de referencia SMPC vs. MPC, SMC trol por modo deslizante y el buen desempeño de los controladores predictivos. Se observa como el MPC presenta un compor- tamiento inestable a diferencia del SMPC y el SMC La robustez del controlador ha sido valorada que logran estabilizar el sistema, sin embargo el SMC por medio de la técnica del gráfico de robustez, al que produce una sobreoscilación un 60 % mayor que el se le ha realizado una modificación a la versión origi- SMPC. nal de manera de estudiar las variaciones de ganancia y cambios en, la frecuencia del polo del modelo. Se ha La Figura 12 presenta la señal de controla para observado como según esta evaluación el controlador MPC y SMPC, se observa como el MPC presenta un propuesto es considerado robusto, y se ha mostrado la comportamiento inestable. poca robustez del MPC convencional. El desempeño de SMPC se probó por medio de un modelo no lineal La Figura 13 muestra la señal de control del de un brazo manipulador con fricción Coulomb y vis- SMC donde se puede observar el efecto de oscilacio- cosa, ha presentado un buen desempeño con el contro- nes de alta frecuencia y una mayor amplitud de la se- lador ajustado en condiciones nominales. 46 Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004
García-Gabín y Zambrano El comportamiento de los controladores se ha [8] Canudas de Wit, C., Astrom, K. J., & Lischinsky, comparado para cambios de referencia, cambios en la P. (1995). “A new model for control of systems masa del brazo manipulador y en la fricción del with friction”. IEEE Transactions on Automatic sistema. Se observó que el SMPC presenta un mejor Control, 40(3), pp. 419-425. comportamiento en cuanto a tiempo de estabilización [9] Radclittle, C., & Southward, S. (1990). “A y sobreoscilación. Se encontró que el SMPC presenta property of stick-slip friction models which una mejor robustez que la propuesta convencional de promotes limit cycle generation”. Proceedings of MPC. Las señales de control del SMPC son más the American Control Conference, San Diego, suaves, eliminando los cambios bruscos en la señal de California, USA, pp. 1198-1203. control del SMC. [10] García-Gabín, W., E.F. Camacho (2003) “Sliding mode model based predictive control for non REFERENCIAS minimum phase systems”. European Control Conference. Cambridge, UK. [1] Camacho, E.F., Bordons, C. (1999). “Model [11] García-Gabín, W. (2003) “Control predictivo por Predictive Control”. Springer-Verlag. London. modo deslizante aplicado a un reactor no lineal [2] Clarke, D., C. Mohtadi and P. Tuffs (1987a). con respuesta inversa”. IV Congreso de “Generalized predictive Control-Part I. Basic Automatización y Control. Mérida, Venezuela. Algorithm”. Automatica. Vol. 23, pp.137-148. [12] Shinskey, F. G. (1990) “How good are our [3] Clarke, D., C. Mohtadi and P. Tuffs (1987b). controller in absolute performance and “Generalized predictive Control-Part II. robustness?”. Measurement and Control, Vol. Extensions and Interpretations”. Automatica. Vol. 23, pp. 114-121. 23, pp. 149-160. [13] Camacho, O. E., García-Gabín, W., and Rojas, R. [4] Utkin, V. Sliding (1992). “Sliding Modes in D., (1998), “Sliding Mode Control: A Robust Control Optimization”. Springer-Verlag. Moscow. Approach to Integrating Systems with [5] Craig J. (1989). “Introduction to robotics”. Deadtime”, Proceedings of IEEE-ICCDCS, Addison-Wesley. New York. Margarita, Venezuela, pp. 401–406. [6] Jee-Hwan Ryu, Jinil Song, Dong-Soo Kwon. [14] Camacho O.; Rojas R. “A General Sliding Mode (2001). “A nonlinear friction compensation Controller For Nonlinear Chemical Processes”. method using adaptive control and its practical ASME J. of Dynamic Systems, Measurements application to an in-parallel actuated 6-DOF and Control. (2000), 122, pp. 650-655. manipulator”. Control Engineering Practice (9), pp. 159-167. [7] Armstrong-Helouvry, B., Dupont, P. E., & Canudas de Wit, C. (1994). “A survey of analysis tools and compensation methods for control of machines with friction”. Automatica, 30(7), pp. 1198-1203. . Rev. INGENIERÍA UC. Vol. 11, No 2, Agosto 2004 47
También puede leer