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

기술 요약: 유도파 손상 국지화를 위한 물리 기반 단일 원자 매칭 퍼서스(PISAMP)

문제 정의
박판 구조물의 구조 건전성 모니터링(SHM)은 압전 변환기(PZT)에 의해 생성된 초음파 유도파(GW)에 의존하는 경우가 많습니다. 이 분야의 주요 과제는 단일 액추에이터에 의해 생성되는 다중 분산 파동 모드(특히 $A_0$ 및 $S_0$ 모드)를 동시에 여기시키는 것입니다. 생성된 파동 패킷은 경계면과 잠재적 손상과 상호작용하며, 이 과정에서 중첩되고 분산된 기여도를 가진 복잡한 신호를 형성합니다. 기존의 신호 처리 기술은 다음과 같은 상당한 한계에 직면해 있습니다:

• 사전 구축된 딕셔너리 기반 방법(예: 표준 매칭 퍼서스)은 가능한 모든 모드 기여도를 포함하는 딕셔너리를 미리 계산해야 하므로, 연속적인 구조물에 적용하기가 비실용적이며 적응성을 제한합니다.
• 비딕셔너리 기반 방법(예: EMD, VMD 또는 데이터 기반 오토인코더)은 파동 전파의 근저에 있는 물리 법칙과의 연결성이 부족하여, 추출된 특징이 수학적으로는 타당할지라도 물리적으로는 해석 불가능한 결과를 초래할 수 있습니다.

본 연구에서 다루는 핵심 난제는, 구조물의 특정 분산 곡선에 대한 사전 지식 없이도 파동수의 기능(wavenumber function)과 전파 거리를 직접 추출할 수 있도록, 물리적 전파 원리에 기반하면서도 계산 효율적인 신호 분해 프레임워크를 구축하는 것입니다.

방법론: 물리 기반 단일 원자 매칭 퍼서스 (PISAMP)
저자들은 측정된 신호를 외부 여기(excitation)의 지연 및 분산된 임펄스 응답으로 모델링하는 신호 분해 기법인 물리 기반 단일 원자 매칭 퍼서스(PISAMP) 방법을 제안합니다. 이 방법론은 다음 원칙을 바탕으로 구축되었습니다:

1. 물리적 정식화: 이 방법은 주파수 영역에서 신호를 모델링합니다. 측정된 신호 $\hat{s}(\omega)$는 $N$개의 원자(atom)의 선형 결합으로 근사화됩니다. 여기서 각 원자는 특정 전파 거리 $d_i$와 특정 파수 함수 $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}$$
   여기서 $\hat{x}(\omega)$는 여기 신호이고, $\alpha_i$는 진폭(감쇠 모델링)이며, 지수 항은 파수 함수 $k_i(\omega)$를 통해 전파 거리와 분산 효과를 모두 인코딩합니다.

2. 반복적 최적화: 이 문제는 탐욕적(greedy)이고 점진적인 접근 방식을 사용하여 해결됩니다. 각 항 $i$에 대해, 알고리즘은 재구성 오차를 최소화하기 위해 진폭 $\alpha_i$, 전파 거리 $d_i$, 그리고 파수 함수 $k_i(\omega)$를 공동으로 최적화합니다.

   • 초기화: 파수는 비분산 추정치(예: $S_0$ 모드)를 사용하여 초기화됩니다.
   • 진폭 및 거리: 파수를 이전 반복 단계의 값으로 고정한 상태에서, 거리에 따른 오차를 최소화함으로써 진폭과 거리를 결정합니다.
   • 파수 업데이트: 파수 함수는 1차 테일러 선형화를 사용하여 업데이트됩니다. 보정값 $\Delta k_i(\omega)$는 매끄러움과 수치적 효율성을 보장하기 위해 제2종 체비쇼프 다항식(Chebyshev polynomials of the second kind)을 사용하여 근사화됩니다.

3. 스케일 교정 및 손상 국지화: 지수 항의 스케일 부정성(scale indeterminacy), 즉 $k_i(\omega)$와 $d_i$가 서로 반비례하게 스케일링되어도 곱 $k_i(\omega)d_i$가 일정하게 유지되는 문제를 해결해야 합니다. 이를 위해:

   • 본 방법은 손상이 없는 구조에서의 액추에이터-센서 간 직접 경로를 기준으로 사용하여, 수치적 거리를 물리적 단위로 변환하는 글로벌 스케일링 인자를 정의합니다.
   • 전체 신호에 비해 진폭은 미미하지만 국지화에 매우 중요한 "최초 도착(first arrival)" 파동 패킷을 구별하기 위해 수정된 오차 기준을 도입합니다.
   • 물리적 거리가 확립되면, **타원형 국지화(Elliptical Localization)**를 사용하여 손상 위치를 결정합니다. 액추에이터에서 손상을 거쳐 센서로 이어지는 총 경로 길이는 두 유클리드 거리의 합으로 모델링됩니다. 여러 센서로부터 추정된 경로 길이에 가장 잘 부합하는 손상 좌표를 찾기 위해 비선형 최소제곱 문제(Gauss-Newton 또는 Levenberg-Marquardt 알고리즘을 통해 해결)가 사용됩니다.

주요 결과
본 논문은 수치적 조사를 통해 PISAMP 방법의 유효성을 검증하였습니다:

• 신호 근사: 등방성 판재를 대상으로 한 학술적 사례에서, PISAMP는 6개의 서브모드를 사용하여 4%의 재구성 오차를 달 dengan 했으며, 비분산($S_0$) 모드와 분산($A_0$) 모드를 성공적으로 구분했습니다.
• 복잡한 구조물: 에어버스 A380 팬 카울 구조(보강재가 있는 복합재 판재)의 실험 데이터에 적용했을 때, 이 방법은 40개의 서브모드를 사용하여 20%의 재구성 오차를 달성했으며, 반사 및 재료 감쇠를 효과적으로 처리했습니다.
• 손상 국지화: 30가지의 서로 다른 손상 구성이 포함된 2D 정사각형 판 시뮬레이션에서, 제안된 방식(PIMAP 거리 추정과 타원형 국지화의 결-합)은 x-좌표에서 2.35%, y-좌표에서 **2.26%**의 평균 상대 오차를 달성했습니다. 이 방법은 산란파 성분만을 사용하여 손상 위치를 성공적으로 식별했습니다.

의의 및 주장
저자들은 본 연구의 주요 독창성이 물리적 파동 전파를 압축된 신호 분해 프레임워크에 직접 통합했다는 점에 있다고 주장합니다. 순수 데이터 기반 방식이나 수치 집약적인 기법과 달리, PISAMP는 강력한 물리적 제약 조건을 내장하여 저차원의 효율적인 표현을 생성합니다.

본 연구의 의의는 두 가지 측면에서 나타납니다:

1. 해석 가능성: 모든 가능한 모드를 위한 사전 계산된 딕셔너리 없이도, 물리적으로 의미 있는 특징인 모드 파수 함수와 전파 거리를 직접 식별할 수 있게 합니다.
2. 국지화 능력: 분해된 신호 특징을 기하학적 경로 길이와 직접 연결함으로써, 박판 구조물에서의 손상 국지화를 위한 견고하고 해석 가능한 방법론을 제공합니다.

논문은 현재 구현된 방식이 유클리드 거리가 적용되는 평면 표면에 국한되어 있음을 겸허히 명시하고 있습니다. 향-후 연구로서 측지선 거리(geodesic distance) 계산이 필요한 곡면 기하학으로의 확장과 실제 실험 데이터에 대한 검증이 필요함을 밝히고 있습니다.