On the Characterization of the Duhem Hysteresis Operator with Clockwise Input-Output Dynamics

Sintesi Tecnica: Sulla Caratterizzazione dell'Operatore di Isteresi di Duhem con Dinamica Ingresso-Uscita Oraria

Definizione del Problema
L'isteresi è un fenomeno non lineare prevalente nei sistemi fisici (ad es., attuatori piezoelettrici, materiali ferromagnetici, sistemi meccanici). Mentre i comportamenti ingresso-uscita (I/O) sono spesso caratterizzati come antiorari (CCW) o orari (CW), la letteratura esistente ha analizzato estensivamente le dinamiche CCW (associate alla passività dal derivata dell'ingresso alla derivata dell'uscita), ma manca di un quadro rigoroso per le dinamiche CW nel contesto degli operatori di isteresi di Duhem. In particolare, vi è la necessità di caratterizzare gli operatori di Duhem che esibiscono un comportamento I/O orario, dove la velocità di fornitura (supply rate) è definita dal prodotto dell'uscita e della derivata dell'ingresso ($y\dot{u}$), piuttosto che dell'ingresso e della derivata dell'uscita ($\dot{y}u$). Ciò è particolarmente rilevante per i sistemi meccanici con attrito isteretico, dove la relazione tra spostamento relativo e forza d'attrito segue spesso dinamiche CW.

Metodologia
Gli autori investigano una specifica classe di operatori di isteresi di Duhem definiti dall'equazione differenziale:
$$\dot{y}(t) = f_1(y(t), u(t))\dot{u}_+(t) + f_2(y(t), u(t))\dot{u}_-(t)$$
dove $\dot{u}_+$ e $\dot{u}_-$ rappresentano rispettivamente la parte positiva e negativa della derivata dell'ingresso.

Per caratterizzare la proprietà CW, il documento introduce tre funzioni ausiliarie derivate dalle funzioni definitorie dell'operatore $f_1$ e $f_2$:

1. Funzione di Anisteresi ($f_{an}$): Una curva dove $f_1 = f_2$, che rappresenta il luogo di punti in cui il ciclo di isteresi collassa.
2. Funzione di Traversata ($\omega_\Phi$): Descrive la traiettoria dell'operatore quando viene applicato un ingresso monotono crescente o decrescente da uno stato specifico.
3. Funzione di Intersezione ($\Lambda$): Definisce l'unico punto di intersezione tra la traiettoria di traversata e la curva di anisteresi.

Utilizzando queste funzioni, gli autori costruiscono una funzione di accumulo (storage function) $H_\curvearrowright: \mathbb{R}^2 \to \mathbb{R}^+$ definita come l'integrale della funzione di anisteresi meno l'integrale della funzione di traversata. Il nucleo della metodologia consiste nel derivare condizioni sufficienti su $f_1$ e $f_2$ (specificamente riguardo al segno di una funzione derivata $F$ rispetto alla curva di anisteresi) per garantire che la derivata temporale di questa funzione di accumulo soddisfi la disuguaglianza di dissipazione CW:
$$\frac{d}{dt}H_\curvearrowright(y(t), u(t)) \leq y(t)\dot{u}(t)$$

Contributi Chiave

1. Condizioni Sufficienti per la Dinamica CW: Il documento fornisce condizioni sufficienti esplicite sulle funzioni $f_1$ e $f_2$ di un operatore di Duhem per garantire che l'operatore esibisca dinamiche I/O orarie. Ciò è formalizzato nel Teorema 1.
2. Costruzione Esplicita della Funzione di Accumulo: A differenza dei lavori precedenti che stabilivano l'esistenza, questo articolo offre una costruzione esplicita della funzione di accumulo $H_\curvearrowright$ (Equazione 20) che soddisfa la disuguaglianza di dissipazione CW. Si dimostra che questa funzione è la "funzione di accumulo disponibile" per la velocità di fornitura $y\dot{u}$.
3. Generalizzazione alle Uscite Saturate: I risultati sono estesi nel Proposizione 2 per accomodare operatori di Duhem con uscite saturate, una caratteristica comune nei modelli fisici come il modello di Dahl.
4. Framework di Analisi della Stabilità: Il documento dimostra come la funzione di accumulo costruita possa essere utilizzata come funzione di Lyapunov per analizzare la stabilità di sistemi che incorporano l'isteresi CW.

Risultati

• Validazione Teorica: Gli autori dimostrano che se la funzione costruita $H_\curvearrowright$ soddisfa la disuguaglianza di dissipazione, l'operatore è CW. La prova si basa sulle proprietà della funzione di intersezione $\Lambda$ e sul comportamento della funzione di traversata $\omega_\Phi$.
• Applicazione al Modello di Dahl: La teoria è applicata al modello di Dahl dell'attrito. Il documento deriva esplicitamente la funzione di accumulo $H_\curvearrowright$ per il modello di Dahl (Equazione 26) e verifica che essa soddisfi la richiesta disuguaglianza di dissipazione.
• Stabilità di Sistemi del Secondo Ordine: Utilizzando la funzione di accumulo derivata, gli autori analizzano un sistema meccanico del secondo ordine con attrito isteretico descritto dal modello di Dahl. Costruendo una funzione di Lyapunov composita $V = \frac{1}{2}kx_1^2 + \frac{1}{2}mx_2^2 + H_\curvearrowright(x_3, x_1)$, dimostrano che il sistema è passivo dalla forza applicata alla velocità. Di conseguenza, una semplice legge di controllo proporzionale a feedback ($F = -dx_2$) mostra che la velocità converge asintoticamente a zero, senza richiedere la conoscenza precisa dei parametri di attrito ($\rho$ e $F_C$).

Significatività
Il documento afferma che la sua principale significatività risiede nel fornire un'estensione duale ai risultati di caratterizzazione CCW esistenti, specificamente adattata per sistemi con dinamiche I/O CW. Stabilendo un rigoroso framework di dissipatività per gli operatori di Duhem con comportamento CW, il lavoro consente l'analisi della stabilità di sistemi meccanici con attrito isteretico senza dover invertire l'operatore (il che potrebbe essere fisicamente irrealizzabile o non invertibile). Inoltre, la costruzione esplicita della funzione di accumulo offre uno strumento pratico per l'identificazione del sistema, permettendo ai ricercatori di restringere la classe di operatori di Duhem adattati ai dati di misura in base alle loro proprietà CW. I risultati colmano il divario tra la teoria astratta della dissipatività e la modellazione pratica dell'isteresi indotta dall'attrito nei sistemi di controllo.