가우시안 커널을 이용한 비대칭 호환 결합 기반 페리다이내믹스 모델

초록

본 논문은 기존 결합 기반 페리다이내믹스(BB‑PD)의 비대칭 호환성 문제와 경계 셀 부피 계산의 복잡성을 해결하기 위해, 무한히 정의된 가우시안 커널을 도입한 새로운 BB‑PD 모델(BBPD‑GK)을 제안한다. 가우시안 커널을 직접 메쉬프리 이산화하면 별도의 호환성 보정이나 부피 보정 없이도 연속 탄성 해와 일치함을 증명하고, 에너지 등가성을 통해 물성계수를 도출한다. 또한 손상 모델과 선형화 버전을 제시하고, 2차원 매끄러운 문제, 균열이 있는 판, Kalthoff‑Winkler 고속 충격 실험 등에서 수치 검증을 수행하였다.

상세 분석

본 연구는 결합 기반 페리다이내믹스(BB‑PD)의 근본적인 한계인 “비대칭 호환성(asymptotic compatibility)” 부재와 경계 영역에서의 부피 보정(volume correction) 문제를 근본적으로 재구성한다는 점에서 혁신적이다. 기존 BB‑PD는 유한 반경 δ 내에서만 상호작용을 허용하는 유계(kernel) 함수를 사용했으며, 이로 인해 격자 간격을 무한히 미세화하더라도 비국소 효과가 완전히 사라지지 않아 Navier 방정식과의 일치성이 보장되지 않았다. 또한, 경계 근처 입자들의 상호작용 영역이 부분적으로만 겹치기 때문에 정확한 부피 계산이 어려워, QWJ, PA‑AC, VCA 등 복잡한 보정 알고리즘이 필요했다.

논문은 이러한 문제를 “무한히 정의된 가우시안 커널” G(ξ)= (2πσ²)⁻¹ exp(−|ξ|²/σ²) 로 교체함으로써 해결한다. 가우시안은 이론적으로 무한히 퍼져 있지만, σ가 작아질수록 급격히 감소하므로 실용적으로는 적절한 절단 반경을 두어 유한 영역만을 계산한다. 핵심은 가우시안 커널 자체가 이미 거리 가중치를 포함하고 있어, 입자 간 부피를 별도 보정할 필요가 없다는 점이다.

에너지 등가성 분석을 통해 2차원 경우 β = 3E/(4σ⁴), 3차원 경우 β = 4E/(5σ⁴) 라는 스프링 상수를 도출하였다. 여기서 E는 탄성계수이며, σ는 가우시안 분산이다. 이러한 파라미터는 기존 BB‑PD의 μ(ξ)와 달리 연속 탄성 에너지와 정확히 일치하도록 설계되었다.

손상 모델은 기존 BB‑PD와 동일하게 임계 스트레치 s_c 를 정의하고, 손상 변수 φ(x,t)=1−∫μ dV/∫dV 로 표현한다. 가우시안 커널을 이용한 파괴 에너지 G₀는 해석적으로 적분되어 s_c = s · 8√(2π) G₀/(45 E h σ) 로 얻어진다. 이는 기존에 복잡한 수치 적분이 필요했던 절차를 간단한 폐쇄식으로 대체한다.

선형화 과정에서는 g(η,ξ)≈C(ξ)η 로 전개하고, C(ξ)= (3E/(8πσ⁶)) exp(−|ξ|²/σ²) (ξ⊗ξ) 를 얻는다. 다중 차원 가우시안의 4차 모멘트를 이용해 비국소 연산자를 전통적인 라플라시안 형태인 ∇·(D:ε) 로 변환함으로써, 선형 BBPD‑GK가 정확히 클래식 탄성 방정식과 동등함을 증명한다.

수치 실험에서는 메쉬프리 방식으로 직접 이산화했으며, 매끄러운 변형 문제에서 수렴률이 2차 이상을 보였고, 균열이 있는 판에서 손상 전파가 실험적 파열 패턴과 일치했다. 특히 Kalthoff‑Winkler 고속 충격 실험 재현에서 파동 전파와 균열 각도가 기존 PD 보정 기법보다 더 정확하게 예측되었다.

결과적으로, 가우시안 커널 기반 BBPD‑GK는 (1) 비대칭 호환성을 자연스럽게 만족, (2) 부피 보정 없이도 정확한 힘 전달을 구현, (3) 선형 및 비선형 손상 모델을 일관되게 제공, (4) 계산 비용을 절감하면서도 높은 정확도를 유지한다는 장점을 가진다. 이는 차후 다중 물리·다중 스케일 시뮬레이션에서 페리다이내믹스를 보다 실용적인 도구로 전환시키는 중요한 전환점이 될 것으로 기대된다.