Análisis de sistemas de control, PID

PID – Diseño con el controlador PI (Proporcional-Integral)

Configuración del controlador.

El PID es uno de los controladores más ampliamente utilizados en los esquemas de compensación siguientes, algunas de las cuáles se ilustran en la Figura 10-2:

  1. Compensación en serie (cascada),
  2. Compensación mediante realimentación,
  3. Compensación mediante realimentación de estado,
  4. Compensación en serie-realimentada,
  5. Compensación prealimentada,

En general, la dinámica de un proceso lineal controlado puede representarse mediante el diagrama de la Figura 10-1.

El objetivo de diseño es que las variables controladas, representadas por el vector de salida y(t), se comporten de cierta forma deseada. El problema esencialmente involucra el determinar la señal de control u(t) dentro de un intervalo prescrito para que todas las especificaciones de diseño sean satisfechas (ver Estabilidad de un sistema de control, Respuesta transitoria, Error en estado estable).

El controlador PID aplica una señal al proceso de control mostrado en la Figura 10-1, que es una combinación proporcional, integral y derivada de las señal de actuación e(t). Debido a que estos componentes de la señal se pueden realizar y visualizar fácilmente en el dominio del tiempo, los controladores PID se diseñan comúnmente empleando métodos en el dominio del tiempo.

Después de que el diseñador ha seleccionado una configuración para el controlador, debe escoger además el tipo de controlador. En general, mientras más complejo el controlador, más costoso, menos confiable y más difícil de diseñar. Por ende, en la práctica se selecciona el tipo de controlador más simple que permita cumplir con las especificaciones de diseño, lo que involucra experiencia, intuición, arte y ciencia.

Las componentes integral y derivativo de un controlador PID tienen una implicación individual en el desempeño, y sus aplicaciones requieren un entendimiento de las bases de estos elementos. Ya hemos tratado anteriormente la porción PD (PID – Diseño con el controlador PD (Proporcional-Diferencial)), ahora nos enfocamos en la porción PI.

Diseño con el controlador PI

En muchas ocasiones el controlador PD puede no ser suficiente para cumplir con las especificaciones de diseño a través del proceso de compensación. Procedemos ahora a estudiar la parte integral del controlador PID, que produce una señal proporcional a la integral con respecto al tiempo de la entrada del controlador. La Figura 10-15 muestra mediante una configuración en serie, un controlador PI que compensa un sistema de control para una planta prototipo de segundo orden.

null

La función de transferencia del controlador es:nullDonde Kp y Ki son las constantes proporcional e integral respectivamente. Mientras, la función de transferencia Gp(s) de la planta es:null

La función de transferencia directa G(s) del sistema compensado es:null

Se puede constatar que los efectos inmediatos de aplicar el controlador PI es agregar un cero simple en s=-Ki/Kp y agregar un polo simple en s=0 en la función de transferencia directa, elevando a tipo 2 el sistema, mejorando así el error en estado estable.Sin embargo, se deben seleccionar adecuadamente los valores de los parámetros Kp y Ki para evitar que el sistema pase a ser inestable. Este generalmente es el principal problema a resolver cuando se utiliza un controlador PI, seleccionar Kp y Ki para que la respuesta transitoria sea satisfactoria mientras se mejora el error en estado estable.

Cuando se diseña un controlador PI un método factible es la selección del cero s=-Ki/Kp para que esté localizado cerca del origen y lejos de los polos más significativos del sistema. Por ello los valores de Kp y Ki deben ser relativamente pequeños.

 Diseño electrónico de un PI   

Al diseñar un controlador PI con amplificadores operacionales, obtenemos el circuito de la Figura 10-16:

null
Figura 10-16

Donde:nullAl comparar esta última ecuación con la función de transferencia del PI, observamos que:null

Ejemplo

Para apreciar mejor el efecto del controlador PI, veamos el siguiente ejemplo. Supongamos que tenemos el sistema de la Figura 7-23.

null

La función de transferencia directa G(s) de este sistema viene dada por la siguiente expresión:null

Donde K es la constante del preamplificador.

Las especificaciones de diseño para este sistema son las siguientes:

nullDonde:

  • ess: Error en estado estable debido a una entrada parábola
  • Mp: Sobrepaso máximo
  • Tr: Tiempo de levantamiento
  • Ts: Tiempo de asentamiento

Lo primero que vamos a hacer es analizar el error en estado estable ess del sistema antes de compensar, y ver que tanto cumple o no con el primer requerimiento de diseño (Para un repaso de este tema ver: Error en estado estable de un sistema de control).

Para una entrada parabólica debemos trabajar con la constante de aceleración Ka:nullEsto significa un error infinito para una entrada parabólica:Para mejorar el error en estado estable incorporamos un controlador PI en la trayectoria directa del sistema, el cual ahora tendrá la siguiente función de transferencia:nullAl aumentar la tipología del sistema de tipo 1 a tipo 2 de inmediato mejoramos el error en estado estable. Ahora, el ess debido a la entrada parábola será una constante:null

Es decir:

Este ejercicio ya lo trabajamos con el controlador PD. Para ese caso seleccionamos un valor de K=181.17 (Ejemplo 1 – Diseño de un controlador PD (Proporcional-Diferencial) – Matlab)

Nos permitimos la libertad de considerar el mismo valor para K en este caso con el fin de mantener la respuesta transitoria bajo condiciones aceptables y conocidas al aplicar el PD del ejemplo 1. De ser necesario ajustaremos el valor de K más adelante. Podemos apreciar que para lograr un ess según se especifica en el diseño, mientras más grande sea K, más pequeño tendrá que ser Ki, algo que puede ser conveniente. Con los valores de y ess podemos calcular un primer valor aproximado de Ki para cumplir con los requerimientos:

Lo siguiente que haremos será analizar la estabilidad del sistema debido a que de ello depende enormemente la selección de los parámetros Kp y Ki. Aplicaremos el criterio de Routh-Hurwitz (para un repaso ver Estabilidad de un sistema de control) para calcular los valores límites de los mencionados parámetros de manera tal que el sistema permanezca estable. Para ello requerimos de la ecuación característica que surge de la función de transferencia del sistema en lazo cerrado Gce(s):

null

La ecuación característica del sistema es:

null

Con esta ecuación aplicamos el criterio de Routh-Hurwitz. De esta manera descubrimos que el sistema es estable para el siguiente rango de valores:

null

Este resultado indica que el Ki/Kp del controlador no puede estar muy cerca del cero, por lo que no resulta conveniente un valor tan bajo tal como Ki=0.002215. Otro criterio para seleccionar Ki/Kp es que resulta conveniente seleccionar el cero añadido por el controlador en s=-Ki/Kp para que esté localizado cerca del origen y lejos de los polos más significativos del sistema. Mediante el lugar geométrico de las raíces en Matlab podemos ver cuáles son los polos más significativos de la ecuación característica, suponiendo una relación aceptable entre los parámetros Ki y Kp desde el punto de vista del estudio de estabilidad anterior pero relativamente cerca del origen, digamos Ki/Kp=10, manteniendo Ki constante y variando Kp (Para un repaso del tema ver: El lugar geométrico de las raíces con Matlab) Primero, arreglamos la ecuación característica en su forma 1+G(s)H(s):

Notar que de esta manera el controlador PI está añadiendo un zero en s=-10.

Aplicamos el siguiente comando en Matlab para obtener el lugar geométrico de las raíces de este sistema compensado:

>> s=tf(‘s’)

>> sys=(815265*(s+10))/(s^3+361.2*s^2)

>> rlocus(sys)

Y obtenemos:

Como se puede observar, el polo más significativo de la ecuación característica está ubicado en s=-361. Por tanto, el criterio que debemos utilizar para seleccionar s=-Ki/Kp es:nullCon este resultado, la función de transferencia directa G(s):nullQuedaría simplificada como:null

El término Ki/Kp sería despreciable comparado con la magnitud de s cuando s asume valores a lo largo del lugar geométrico de las raíces que se corresponde con un factor de amortiguamiento relativo conveniente de 0.7<ζ< 0.1. Luego, un zero en cero anula un polo en cero. El sobrepaso máximo debe ser igual o menor al 5%. Esto significa que se desea un factor de amortiguamiento relativo ζ aproximado al siguiente valor:

null

null

Con la ayuda del lugar geométrico podemos ubicar los polos que se corresponden con este valor de ζ:

De acuerdo con la gráfica, el valor de Kp requerido para obtener este factor de amortiguamiento es:

Y por tanto:

También observamos en la gráfica anterior que si Kp=0.0838, entonces las raíces de la ecuación característica (los polos del sistema) están en s1=-175+j184 y en s2=-175-j184 Si miramos alrededor de esas raíces podemos notar que el zero añadido por el controlador en s=-10 está muy cerca del origen comparado con los polos de s1 y s2, prácticamente cancelando un polo en el origen, ratificando así la aproximación que hicimos anteriormente para la función de transferencia directa de este sistema luego de la compensación:

Así:

Podemos observar la respuesta del sistema al escalón unitario de acuerdo con estos resultados parciales y la comparación de los sistemas compensado y sin compensar, mediante la siguiente simulación:

>> Ga=(815265)/(s*(s+361.2))

>> Gd=(68319)/(s*(s+361.2))

>> Gce1=feedback(Ga,1)

>> Gce2=feedback(Gd,1)

>> step(Gce1,Gce2)

La gráfica anterior, con el sistema después de la compensación en línea roja, muestra que el PI mejora el error en estado estable y reduce el sobrepaso, pero a expensas de aumentar significativamente el tiempo de levantamiento. La gráfica también nos muestra que el sobrepaso máximo es de 5%, por tanto se cumple el requerimiento. Es necesario notar que se puede seleccionar otra relación para Ki y Kp que cumpla con el requerimiento y aún mejoren el sobrepaso, por ejemplo Ki/Kp=5, Ki/Kp=2. Sólo hay que prestar atención al tema de la estabilidad del sistema. Con los valores calculados, volvemos a calcular el sobrepaso, y vemos como quedan el tiempo de levantamiento tr y el tiempo de asentamiento ts. El siguiente comando nos facilita el valor de ζ y ωn, el factor de amortiguamiento relativo y la frecuencia natural respectivamente (Para un repaso ver:Respuesta Transitoria de un Sistema de Control).

>> damp(Gce2)

Pole: -1.81e+02 + 1.89e+02i /   -1.81e+02 + 1.89e+02i

Damping: 6.91e-01

Frequency: 2.61e+02

  •  Máximo Sobrepaso (MP):

  • Tiempo de levantamiento (Tr):

>> Gd=(68319)/(s*(s+361.2))

>> sys=feedback(Gd,1)

>> step(sys)

  • Tiempo de asentamiento (Ts):null

Cosa que también podemos constatar en la gráfica de respuesta al escalón unitario generada anteriormente:

Podemos concluir que el sistema compensado cumple con los requerimientos del diseño, aunque se sobrepasa un poco en el tiempo de asentamiento. Este último se puede reducir levemente con una relación Ki/Kp=2pero a expensas de acercarse demasiado a la zona inestable del sistema.

ANTERIOR: PID – Diseño con el controlador PD (Proporcional-Diferencial)

Escrito por: Larry Francis Obando – Technical Specialist – Educational Content Writer.

Mentoring Académico / Empresarial / Emprendedores

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

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.

Contacto: España +34 633129287

Caracas, Quito, Guayaquil, Cuenca

WhatsApp: +34 633129287 +593998524011

Atención Inmediata !!

Twitter: @dademuch

FACEBOOK: DademuchConnection

email: dademuchconnection@gmail.com

Atención:

Si lo que Usted necesita es resolver con urgencia un problema de “Sistema Masa-Resorte-Amortiguador” (encontrar la salida X(t), gráficas en Matlab del sistema de 2do Orden y parámetros relevantes, etc.), o un problema de “Sistema de Control Electromecánico” que involucra motores, engranajes, amplificadores diferenciales, etc…para entregar a su profesor en dos o tres días, o con mayor urgencia…o simplemente necesita un asesor para resolver el problema y estudiar para el próximo examen…envíeme el problema…Yo le resolveré problemas de Sistemas de Control, le entrego la respuesta en digital y le brindo una video-conferencia para explicarle la solución…incluye además simulación en Matlab.

Relacionado:

Respuesta Transitoria de un Sistema de Control

Estabilidad de un sistema de control

Simulación de Respuesta Transitoria con Matlab – Introducción

 

 

3 comentarios en “PID – Diseño con el controlador PI (Proporcional-Integral)”

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión /  Cambiar )

Google photo

Estás comentando usando tu cuenta de Google. Cerrar sesión /  Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión /  Cambiar )

Conectando a %s