필터링 재료 특성으로 FFT 기반 수치 균질화 개선

초록

본 논문은 이산 푸리에 변환(FFT) 기반 수치 균질화 방법에 재료 특성 필터링을 도입함으로써 격자 정밀도와 계산 시간을 크게 절감하고, 국부장에 나타나는 인위적 진동을 감소시키는 새로운 접근법을 제시한다. Voigt·Reuss 필터와 비교해 제안된 “2‑layer” 필터는 강성 포함체와 연성 매트릭스, 혹은 그 반대 상황 모두에서 일관된 성능 향상을 보인다.

상세 분석

FFT 기반 수치 균질화는 1990년대 초반부터 복합재료의 탄성 거동을 효율적으로 해석하는 도구로 자리 잡았다. 기존 연구에서는 연산 속도 향상을 위해 전처리 기법, 고차원 스펙트럼 보정, 그리고 변형된 이산 그린 연산자 등을 도입했지만, 재료 파라미터 자체를 조정하는 접근은 상대적으로 적었다. 본 논문은 이러한 공백을 메우기 위해 “재료 특성 필터링”이라는 개념을 도입한다. 구체적으로는 각 격자 셀에 할당된 탄성 텐서를 주변 이웃 셀들의 평균값으로 대체하는 전통적인 Voigt(등변형)와 Reuss(등응력) 필터를 검토하고, 두 개의 연속된 층을 가정한 “2‑layer” 모델을 새롭게 제안한다.

2‑layer 필터는 포함체와 매트릭스 사이에 존재하는 계면을 두 개의 등가 층으로 근사함으로써, 실제 미세구조에서 발생하는 응력·변형률 집중 현상을 보다 정확히 포착한다. 이때 각 층의 유효 탄성 상수는 체적 평균법과 전단 연속성을 동시에 만족하도록 계산된다. 실험적으로는 2‑D 및 3‑D 이중상(두 물질) 구조를 대상으로 격자 해상도를 32×32에서 256×256까지 변화시키며, 전통적인 FFT 솔버와 필터링된 버전을 비교하였다. 결과는 두드러졌다. 2‑layer 필터를 적용하면 동일한 정확도를 유지하면서 최소 30 % 이상의 격자 수를 감소시킬 수 있었으며, 이는 곧 계산 시간과 메모리 요구량의 비례적 감소로 이어졌다. 또한, 국부 응력·변형률 필드에서 관찰되는 고주파 진동(스파키) 역시 현저히 억제되었다.

흥미로운 점은 Voigt와 Reuss 필터가 미세구조에 따라 성능이 크게 달라지는 반면, 2‑layer 필터는 강성 포함체가 연성 매트릭스에 삽입된 경우와 그 반대 경우 모두에서 일관된 개선을 보였다는 것이다. 이는 2‑layer 모델이 계면 효과를 보다 물리적으로 반영하기 때문으로 해석된다. 논문은 이러한 결과가 선형 탄성뿐 아니라 열전도, 전기 전도 등 다른 물리적 현상에도 동일하게 적용될 가능성을 제시한다.

요약하면, 재료 특성 필터링은 기존 FFT 기반 균질화 기법에 간단히 추가할 수 있는 전처리 단계이며, 특히 복합재료와 같이 상이한 강성을 가진 구성요소가 혼재하는 경우에 큰 효율성을 제공한다. 향후 연구에서는 다중상(multi‑phase) 시스템, 비선형 재료, 그리고 비정형 격자에 대한 확장 가능성을 탐색할 필요가 있다.