Physics-Informed Single Atom Matching Pursuit: Guided-Waves Wavenumbers and Propagation Distance Estimation for Damage Localization in Structural Health Monitoring

Sintesi Tecnica: Ricerca di Corrispondenza a Singolo Atomo Informata dalla Fisica per la Localizzazione di Danni tramite Onde Guidate

Enunciato del Problema
Il monitoraggio dello stato strutturale (SHM) di componenti a parete sottile si basa spesso su onde guidate ultrasoniche (GW) generate da trasduttori piezoelettrici (PZT). Una sfida primaria in questo dominio è l'eccitazione simultanea di molteplici modi di onde dispersive (specificamente i modi $A_0$ e $S_0$) mediante un singolo attuatore. Man mano che questi pacchetti d'onda si propagano, interagiscono con i bordi e potenziali danni, creando segnali complessi con contributi sovrapposti e dispersi. Le tecniche convenzionali di elaborazione del segnale affrontano limitazioni significative:

• Metodi basati su dizionario (ad esempio, il Matching Pursuit standard) richiedono dizionari pre-calcolati di tutti i possibili contributi modali, il che è impraticabile per strutture continue e limita l'adattabilità.
• Metodi non basati su dizionario (ad esempio, EMD, VMD o AutoEncoder basati sui dati) spesso mancano di un forte collegamento con la fisica sottostante della propagazione delle onde, risultando in caratteristiche estratte che sono matematicamente solide ma fisicamente non interpretabili.

La difficoltà centrale affrontata in questo lavoro è la necessità di un framework di decomposizione che sia al contempo computazionalmente efficiente e fondato sui principi fisici della propagazione delle onde, consentendo l'estrazione diretta delle funzioni del numero d'onda modale e delle distanze di propagazione senza la conoscenza pregressa delle specifiche curve di dispersione della struttura.

Metodologia: Physics-Informed Single Atom Matching Pursuit (PISAMP)
Gli autori propongono il metodo Physics-Informed Single Atom Matching Pursuit (PISAMP), una tecnica di decomposizione del segnale che modella i segnali misurati come risposte all'impulso ritardate e disperse di un'eccitazione esterna. La metodologia si basa sui seguenti principi:

1. Formulazione Fisica: Il metodo modella il segnale nel dominio della frequenza. Un segnale misurato $\hat{s}(\omega)$ è approssimato come una combinazione lineare di $N$ atomi, dove ogni atomo rappresenta un pacchetto d'onda che viaggia una specifica distanza $d_i$ con una specifica funzione del numero d'onda $k_i(\omega)$:
   $$\hat{s}(\omega) \approx \sum_{i=1}^{N} \alpha_i(d_i) \hat{x}(\omega) e^{-j k_i(\omega) d_i}$$
   Qui, $\hat{x}(\omega)$ è l'eccitazione, $\alpha_i$ è l'ampiezza (che modella l'attenuazione), e il termine esponenziale codifica sia la distanza di propagazione che gli effetti dispersivi tramite la funzione del numero d'onda $k_i(\omega)$.

2. Ottimizzazione Iterativa: Il problema viene risolto utilizzando un approccio greedy e incrementale. Per ogni termine $i$, l'algoritmo ottimizza congiuntamente l'ampiezza $\alpha_i$, la distanza di propagazione $d_i$ e la funzione del numero d'onda $k_i(\omega)$ per minimizzare l'errore di ricostruzione nel dominio della frequenza.

   • Inizializzazione: Il numero d'onda viene inizializzato utilizzando una stima non dispersiva (ad esempio, il modo $S_0$).
   • Ampiezza e Distanza: Questi vengono determinati minimizzando l'errore rispetto alla distanza, trattando il numero d'onda come fisso rispetto all'iterazione precedente.
   • Aggiornamento del Numero d'Onda: La funzione del numero d'onda viene aggiornata utilizzando una linearizzazione di Taylor del primo ordine. La correzione $\Delta k_i(\omega)$ è approssimata utilizzando polinomi di Chebyshev di secondo tipo per garantire fluidità e un'efficienza numerica.

3. Calibrazione della Scala e Localizzazione del Danno: Una sfida critica identificata è l'indeterminazione di scala del termine esponenziale, dove il prodotto $k_i(\omega)d_i$ rimane costante anche se $k_i$ e $d_i$ vengono scalati inversamente. Per risolvere questo:

   • Il metodo utilizza il percorso diretto attuatore-sensore in una struttura non danneggiata come riferimento per definire un fattore di scala globale, convertendo le distanze numeriche in unità fisiche.
   • Viene introdotto un criterio di errore modificato per identificare specificamente i pacchetti d'onda di "primo arrivo", che spesso hanno un'ampiezza trascurabile rispetto al segnale totale ma sono cruciali per la localizzazione.
   • Una volta stabilite le distanze fisiche, la posizione del danno viene determinata tramite Localizzazione Ellittica. La lunghezza totale del percorso dall'attuatore al danno e poi al sensore è modellata come la somma di due distanze euclidee. Un problema di minimizzazione dei quadrati non lineari (risolto tramite algoritmi Gauss-Newton o Levenberg-Marquardt) viene utilizzato per trovare le coordinate del danno che meglio si adattano alle lunghezze di percorso stimate da più sensori.

Risultati Chiave
Il documento valida il metodo PISAMP attraverso investigazioni numeriche:

• Approssimazione del Segnale: In un caso accademico riguardante una piastra isotropa, PISAMP ha raggiunto un errore di ricostruzione del 4% utilizzando 6 submodi, distinguendo con successo tra i modi non dispersivi ($S_0$) e quelli dispersivi ($A_0$).
• Strutture Complesse: Applicato a dati sperimentali da una struttura di carter di ventola di un Airbus A380 (una piastra composita con irrigidimenti), il metodo ha raggiunto un errore di ricostruzione del 20% utilizzando 40 submodi, gestendo efficacemente riflessioni e smorzamento del materiale.
• Localizzazione del Danno: In una simulazione di una piastra quadrata 2D con 30 diverse configurazioni di danno, l'approccio proposto (combinando la stima della distanza PISAMP con la Localizzazione Ellittica) ha ottenuto un errore relativo medio del 2,35% per la coordinata x e del 2,26% per la coordinata y. Il metodo ha identificato con successo le posizioni del danno utilizzando solo le componenti delle onde diffuse.

Significatività e Rivendicazioni
Gli autori affermano che la principale novità di questo studio risiede nell'integrare la fisica della propagazione delle onde direttamente in un framework di decomposizione del segnale compatto. A differenza delle tecniche puramente basate sui dati o numericamente intensive, PISAMP incorpora forti vincoli fisici per produrre una rappresentazione a bassa dimensionalità ed efficiente dal punto di vista computazionale.

La significatività del lavoro è duplice:

1. Interpretabilità: Consente l'identificazione diretta di caratteristiche fisicamente significative, specificamente le funzioni del numero d'onda modale e le distanze di propagazione, senza richiedere un dizionario pre-calcolato di tutti i modi possibili.
2. Capacità di Localizzazione: Fornisce una metodologia robusta e interpretabile per la localizzazione dei danni in strutture a piastra sottile, collegando direttamente le caratteristiche del segnale decomposto alle lunghezze dei percorsi geometrici.

Il documento nota modestamente che l'implementazione attuale è limitata a superfici planari dove si applica la distanza euclidea. Identifica l'estensione a geometrie curve (che richiedono il calcolo della distanza geodetica) e la validazione con dati sperimentali reali come passi necessari per la ricerca futura.