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

Technische Samenvatting: Physics-Informed Single Atom Matching Pursuit voor Damage-lokalisatie in Geleidingsgolven

Probleemstelling
Structural Health Monitoring (SHM) van dunwandige componenten vertrouwt vaak op ultrasone geleidingsgolven (Ultrasonic Guided Waves, GWs) die worden gegenereerd door piëzo-elektrische transducers (PZTs). Een primaire uitdaging in dit domein is de gelijktijdige excitatie van meerdere dispersieve golfmodi (specifiek de $A_0$- en $S_0$-modi) door een enkele actuator. Terwijl deze golfpakketten voortplanten, interageren ze met grensvlakken en potentiële schade, wat resulteert in complexe signalen met overlappende en gedispergeerde bijdragen. Conventionele signaalverwerkingstechnieken kampen met aanzienlijke beperkingen:

• Dictionary-gebaseerde methoden (bijv. standaard Matching Pursuit) vereisen vooraf berekende dictionaries van alle mogelijke modale bijdragen, wat onpraktisch is voor continue structuren en de aanpasbaarheid beperkt.
• Niet-dictionary-gebaseerde methoden (bijv. EMD, VMD of data-gestuurde AutoEncoders) missen vaak een sterke verbinding met de onderliggende fysica van golfvoortplanting, wat resulteert in geëxtraheerde kenmerken die wiskundig weliswaar correct zijn, maar fysiek oninterpreteerbaar.

De kern van het probleem dat in dit werk wordt aangepakt, is de noodzaak voor een decompositiekader dat zowel computationeel efficiënt als geworteld is in de fysische principes van golfvoortplanting, waardoor directe extractie van modale golfgetalfuncties en voortplantingsafstanden mogelijk is zonder voorafgaande kennis van de specifie specificiteit van de dispersiecurves van de structuur.

Methodologie: Physics-Informed Single Atom Matching Pursuit (PISAMP)
De auteurs stellen de Physics-Informed Single Atom Matching Pursuit (PISAMP) methode voor, een signaaldecompositietechniek die gemeten signalen modelleert als vertraagde en gedispergeerde impulsresponsen van een externe excitatie. De methodologie is gebouwd op de volgende principes:

1. Fysische Formulering: De methode modelleert het signaal in het frequentiedomein. Een gemeten signaal $\hat{s}(\omega)$ wordt benaderd als een lineaire combinatie van $N$ atomen, waarbij elk atoom een golfpakket vertegenwoordigt dat een specifieke afstand $d_i$ aflegt met een specifieke golfgetalfunctie $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}$$
   Hierbij is $\hat{x}(\omega)$ de excitatie, $\alpha_i$ de amplitude (die attenuatie modelleert) en encodeert de exponentiële term zowel de voortplantingsafstand als de dispersieve effecten via de golfgetalfunctie $k_i(\omega)$.

2. Iteratieve Optimalisatie: Het probleem wordt opgelost met behulp van een gulzig (greedy), incrementeel algoritme. Voor elke term $i$ optimaliseert het algoritme gezamenlijk de amplitude $\alpha_i$, de voortplantingsafstand $d_i$ en de golfgetalfunctie $k_i(\omega)$ om de reconstructiefout in het frequentiedomein te minimaliseren.

   • Initialisatie: Het golfgetal wordt geïnitialiseerd met een niet-dispersieve schatting (bijv. de $S_0$-modus).
   • Amplitude en Afstand: Deze worden bepaald door de fout met betrekking tot de afstand te minimaliseren, waarbij het golfgetal als vaststaand wordt beschouwd uit de vorige iteratie.
   • Golfgetal Update: De golfgetalfunctie wordt bijgewerkt met behulp van een eerste-orde Taylor-linearisatie. De correctie $\Delta k_i(\omega)$ wordt benaderd met behulp van Chebyshev-polynomen van de tweede soort om de vloeiendheid en numerieke efficiëntie te waarborgen.

3. Schaalcalibratie en Damage-lokalisatie: Een kritieke uitdaging die is geïdentificeerd, is de schaal-onbepaaldheid van de exponentiële term, waarbij het product $k_i(\omega)d_i$ constant blijft, zelfs als $k_i$ en $d_i$ invers met elkaar worden geschaald. Om dit op te lossen:

   • Gebruikt de methode het directe pad van de actuator naar de sensor in een onbeschadigde structuur als referentie om een globale schalfactor te definiëren, waardoor numerieke afstanden worden omgezet in fysieke eenheden.
   • Er is een gemodificeerd foutcriterium geïntroduceerd dat specifiek gericht is op het identificeren van "first arrival" golfpakketten, die vaak een verwaarloosbare amplitude hebben ten opzichte van het totale signaal, maar cruciaal zijn voor lokalisatie.
   • Zodra de fysieke afstanden zijn vastgesteld, wordt de locatie van de schade bepaald met behulp van Elliptische Lokalisatie. De totale padlengte van de actuator naar de schade en vervolgens naar de sensor wordt gemodelleerd als de som van twee Euclidische afstanden. Een non-lineair kleinste kwadratenprobleem (opgelost via Gauss-Newton of Levenberg-Marquardt algoritmen) wordt gebruikt om de schade-coördinaten te vinden die het beste passen bij de geschatte padlengtes van meerdere sensoren.

Belangrijkste Resultaten
Het artikel valideert de PISAMP-methode door middel van numerieke onderzoeken:

• Signaalbenadering: In een academische casus betreffende een isotrope plaat bereikte PISAMP een reconstructiefout van 4% met 6 submodi, waarbij het succesvol onderscheid maakte tussen niet-dispersieve ($S_0$) en dispersieve ($A_0$) modi.
• Complexe Structuren: Toegepast op experimentele gegevens van een Airbus A380 fan cowl structuur (een composietplaat met verstevigingen), behaalde de methode een reconstructiefout van 20% met 40 submodi, waarbij reflecties en materiaaldemping effectief werden afgehandeld.
• Damage-lokalisatie: In een 2D vierkante plaat simulatie met 30 verschillende schadeconfiguratiesen bereikte de voorgestelde aanpak (combinatie van PISAMP afstandschatting met Elliptische Lokalisatie) een gemiddelde relatieve fout van 2,35% voor de x-coördinaat en 2,26% voor de y-coördinaat. De methode identificeerde de locatie van de schade succesvol met behulp van enkel de verstrooide golfcomponenten.

Betekenis en Claims
De auteurs stellen dat de belangrijkste nieuwheid van dit onderzoek ligt in de integratie van de fysica van golfvoortplanting direct in een compact signaaldecompositiekader. In tegen tegenstelling tot puur data-gestuurde of numeriek intensieve technieken, implementeert PISAMP sterke fysische restricties om een laag-dimensionale, computationeel efficiënte representatie te produceren.

De betekenis van het werk is tweeledig:

1. Interpreteerbaarheid: Het maakt directe identificatie mogelijk van fysiek betekenisvolle kenmerken, specifiek modale golfgetalfuncties en voortplantingsafstanden, zonder dat een vooraf berekende dictionary van alle mogheden nodig is.
2. Lokalisatievermogen: Het biedt een robuuste, interpreteerbare methodologie voor de lokalisatie van schade in dunne plaatstructuren door de gedecomposeerde signaalkenmerken direct te koppelen aan geometrische padlengtes.

Het artikel merkt bescheiden op dat de huidige implementatie beperkt is tot vlakke oppervlakken waar Euclidische afstand van toepassing is. Het identificeert de uitbreiding naar gebogen geometrieën (die geodetische afstandberekeningen vereisen) en validatie met real-world experimentele gegevens als noodzakelijke stappen voor toekomstig onderzoek.