Stability Analysis and Controller Design for a Linear System with Duhem Hysteresis Nonlinearity

Resumen Técnico: Análisis de Estabilidad y Diseño de Control para un Sistema Lineal con No linealidad de Histéresis de Duhem

Planteamiento del Problema
Este artículo aborda el análisis de estabilidad y el diseño de control para interconexiones de retroalimentación entre un sistema lineal y un operador de histéresis de Duhem. La histéresis es una no linealidad prevalente en diversos sistemas, incluyendo actuadores piezoeléctricos, materiales ferromagnéticos y fricción mecánica. Si bien existen varios modelos (p. ej., Preisach, histéresis de zona muerta [backlash], elasto-plástico), este trabajo se centra en el operador de Duhem debido a su capacidad para capturar dinámicas complejas de entrada-salida (I/O). Específicamente, los autores investigan sistemas donde la planta lineal y el operador de histéresis exhiben dinámicas de I/O de sentido antihorario (CCW) o de sentido horario (CW). El desafío central es establecer condiciones de estabilidad y diseñar controladores para estos sistemas interconectados sin requerir un conocimiento paramétrico preciso del operador de histéresis, basándose en su lugar en las propiedades estructurales de I/O.

Metodología
Los autores emplean un enfoque basado en la disipatividad, utilizando funciones de almacenamiento (funciones de Lyapunov) derivadas de las características específicas de I/O de los subsistemas.

1. Caracterización de las Dinámicas:

   • Dinámicas CCW (Antihorarias): Definidas por la desigualdad $\liminf_{T\to\infty} \int_0^T \langle \dot{y}(t), u(t) \rangle dt > -\infty$. Para los operadores de Duhem, esto se caracteriza por una función de almacenamiento $H^\circleleft$ que satisface $\frac{d}{dt}H^\circleleft \leq \dot{y}u$.
   • Dinámicas CW (Horarias): Definidas por la desigualdad dual $\liminf_{T\to\infty} \int_0^T y(t)^T \dot{u}(t) dt > -\infty$. Para los operadores de Duhem, esto se caracteriza por una función de almacenamiento $H^\circleright$ que satisface $\frac{d}{dt}H^\circleright \leq y\dot{u}$.
   • El artículo revisa las condiciones sobre las funciones de Duhem $f_1$ y $f_2$ (específicamente involucrando una función de anhistéresis $f_{an}$ y funciones de trayectoria) que garantizan estas propiedades CCW o CW.

2. Análisis de Estabilidad de las Interconexiones:
   El artículo analiza cuatro casos distintos de interconexión entre un sistema lineal $P$ y un operador de Duhem $\Phi$ (resumidos en la Tabla I del texto original):

   • Caso (a): Tanto $P$ como $\Phi$ son CCW (Retroalimentación positiva). La estabilidad se demuestra utilizando una función de Lyapunov compuesta que involucra la forma cuadrática del sistema lineal y la función de almacenamiento de Duhem $H^\circleleft$. Las condiciones dependen del límite de sector de la función de anhistéresis $f_{an}$.
   • Caso (b): $P$ es CW y $\Phi$ es CCW (Retroalimentación negativa). La estabilidad se establece mediante Desigualdades Matriciales Lineales (LMI) que son independientes de los parámetros específicos de $\Phi$, basándose únicamente en su propiedad CCW.
   • Caso (c): $P$ es CCW y $\Phi$ es CW (Retroalimentación negativa). La estabilidad se demuestra mediante la propiedad S-CCW del sistema lineal y la función de almacenamiento de Duhem $H^\circleright$.
   • Caso (d): Tanto $P$ como $\Phi$ son CW (Retroalimentación positiva). La estabilidad se analiza utilizando representaciones extendidas de espacio de estados y LMI, requiriendo acotaciones sobre las funciones $f_1$ y $f_2$.

3. Diseño de Controladores:
   Basándose en los teoremas de estabilidad, los autores proponen una metodología de diseño de control. El controlador $C$ se diseña para asegurar que el sistema lineal en cascada (planta + controlador) satisfaga las condiciones de LMI o algebraicas específicas requeridas para el caso de interconexión correspondiente. Esto permite la estabilización sin necesidad de invertir el operador de histéresis.

Contribuciones Clave y Resultados

• Condiciones de Estabilidad Suficientes: El artículo proporciona condiciones de estabilidad suficientes y rigurosas para sistemas de bucle cerrado que involucran histéresis de Duhem. Estas condiciones se categorizan según la naturaleza CCW/CW de la planta y del operador de histéresis.
• Construcción de Lyapunov: Se proporcionan construcciones explícitas de funciones de Lyapunov para cada caso de interconexión, combinando la energía del sistema lineal con las funciones de almacenamiento ($H^\circleleft$ o $H^\circleright$) del operador de Duhem.
• Síntesis de Controladores: Se presenta un algoritmo sistemático para diseñar controladores para plantas lineales con actuadores o sensores histeréticos. El método se basa en la resolución de LMI o desigualdades algebraicas, evitando la necesidad de un modelado preciso de la histéresis o de la compensación inversa.
• Validación Numérica: Los resultados teóricos se validan mediante ejemplos numéricos que involucran un sistema masa-amortiguador-resorte con actuadores histeréticos tanto CCW (tipo piezoeléctrico) como CW (tipo fricción). Las simulaciones demuestran que los controladores diseñados estabilizan con éxito los sistemas, llevando las velocidades a cero y las posiciones a valores constantes.

Significancia y Afirmaciones
Los autores afirman que este trabajo ofrece un paradigma diferente para controlar sistemas con histéresis. A diferencia de los enfoques tradicionales que a menudo requieren definir un operador de histéresis inverso para cancelar la no linealidad, este método aprovecha las propiedades disipativas inherentes (CCW o CW) de la histéresis. En consecuencia, el diseño de control no requiere información precisa sobre los parámetros del operador de histéresis, solo sus dinámicas cualitativas de I/O. El artículo sostiene que este enfoque de disipatividad puede extenderse a plantas no lineales que exhiban dinámicas CCW, ampliando la aplicabilidad de los resultados más allá de los sistemas lineales. El trabajo está motivado por y se construye sobre resultados previos relacionados con la interconexión de sistemas de imagen negativa (negative imaginary) y sistemas CCW, extendiendo estos conceptos para incluir dinámicas CW y la estructura específica del operador de Duhem.