초록
이 논문에서는 반무한 고체영역의 표면에서 측정한 변위응답의 시간이력으로부터 유한요소망 연속기법을 이용해 탄성파 속도의 공간적 분포를 추정하는 역해석 문제를 소개한다. 반무한 영역에서의 역해석을 위해서는 해석 대상이 되는 유한영역의 경계에서 파동의 반사가 일어나지 않도록 하는 것이 중요하다. 이를 위해 유한영역의 경계면에 perfectly-matchedlayers(PMLs)라는 수치적 파동흡수층을 도입하였고, PML을 경계로 하는 유한영역에서 역해석 문제를 정의하였다. 이 문제를 탄성파동방정식을 구속조건으로 하는 최적화 문제로 표현하였으며, 라그랑주 승수법에 기초한 비구속 최적화 기법에 의해 탄성파속도의 최적 분포를 결정하였다. 해의 정확도와 수렴성을 높이기 위해 유한요소망 연속기법을 도입하여 점진적으로 밀도가 증가하는 유한요소망에 대해 연속적으로 역해석을 수행하였다. 1차원 예제들을 통해 유한요소망 연속기법을 이용한 역해석으로부터 탄성파속도의 분포를 정확히 추정할 수 있음을 확인하였으며, 측정 응답에 노이즈가 존재하는 경우에도 제안한 역해석 기법은 목표 탄성파속도 분포에 근사한 결과를 도출하였다.
This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.