주파수 영역 BEM 효율 향상으로 빠른 과도 탄성동역학 분석

초록

본 논문은 주파수 영역 빠른 경계요소법(Fast BEM)에 지수 윈도우와 주파수 도메인 윈도잉, 그리고 해석 초기값 외삽 기법을 결합해 행렬 조건을 개선하고, Gibbs 진동을 억제하며, 연속 선형 시스템 해결 속도를 높이는 세 가지 전략을 제시한다. 0.7백만 자유도 규모의 사례 실험을 통해 첫 번째와 세 번째 전략이 계산 시간을 크게 단축하고, 두 번째 전략이 시간 영역 응답의 정확성을 확보함을 입증한다.

상세 분석

본 연구는 기존 주파수 영역 BEM이 가진 두드러진 장점인 고속 푸리에 변환 기반 연산 효율성을 유지하면서, 수치적 안정성과 정확성을 동시에 향상시키는 방법론을 제시한다. 첫 번째 전략은 지수 윈도우 함수의 감쇠 파라미터(η)를 크게 설정함으로써 복소수 주파수 평면에서 시스템 행렬의 스펙트럼을 압축하고, 결과적으로 조건수가 크게 개선된다. 이는 반복 해석 시 GMRES·CG와 같은 Krylov 서브스페이스 방법의 수렴 속도를 가속화한다. 다만, 큰 η값은 시간 영역 복원 과정에서 급격한 진동(Gibbs 현상)을 야기할 위험이 있다. 이를 보완하기 위해 두 번째 전략으로 주파수 도메인 윈도잉(예: Hanning, Blackman)을 적용한다. 윈도잉은 고주파 성분을 부드럽게 감쇠시켜 역변환 시 진동을 완화하고, 실제 물리적 응답과의 일치도를 높인다. 세 번째 전략은 연속적인 주파수 샘플에 대해 이전 해의 정보를 활용해 초기 해를 외삽하는 방법이다. 구체적으로, 선형 또는 다항식 회귀를 이용해 현재 주파수에 대한 초기값을 예측함으로써, 각 선형 시스템을 풀 때 필요한 반복 횟수를 현저히 감소시킨다. 이러한 세 가지 전략을 결합하면, 대규모(수백만 자유도) 문제에서도 메모리 사용량은 기존 FFT 기반 BEM과 동일하게 유지하면서도 전체 계산 시간이 30~50% 정도 단축되는 효과를 얻는다. 또한, 시간 영역 결과는 오차가 1% 이하로 수렴하여 실무 적용에 충분한 정확성을 제공한다.