균열 유발 비선형 탄성의 FDEM 시뮬레이션

초록

본 연구는 유한-이산 요소법(FDEM)을 2차원으로 구현하여, 수평 및 경사 균열을 가진 고체 내에서 비선형 탄성 파동 전파를 직접 시뮬레이션한다. 균열면 접촉이 발생할 정도의 충분한 진폭이 가해질 때만 비선형 현상이 나타나며, 정상·전단 접촉이 비선형 탄성 생성에 미치는 영향을 정량적으로 분석한다. 결과는 FDEM가 비선형 초음파 NDT에서 균열 메커니즘을 해석하는 강력한 도구임을 보여준다.

상세 분석

이 논문은 비선형 초음파 비파괴검사(NDT)에서 핵심적인 메커니즘인 균열에 의한 비선형 탄성 현상을 수치적으로 재현하기 위해, 유한 요소법(FEM)과 이산 요소법(DEM)의 장점을 결합한 FDEM을 2차원 형태로 구현하였다. 기존 FEM 기반 모델은 균열면의 접촉·마찰을 정확히 표현하기 어려웠지만, FDEM에서는 균열 양측의 유한 요소를 별도의 이산 요소로 취급해 정상(N) 및 전단(T) 상호작용을 명시적으로 계산한다. 이를 통해 균열면이 닫히거나 미끄러질 때 발생하는 비선형 응답을 직접 추적할 수 있다.

모델링 단계에서 저자들은 먼저 균열이 없는 완전 연속체와 수평·경사 균열을 각각 포함한 시뮬레이션 시나리오를 설계하였다. 파동은 고정된 경계에서 사인 파형으로 입력되며, 진폭을 단계적으로 증가시켜 접촉 임계점을 탐색한다. 결과는 진폭이 임계값 이하일 때는 선형 전파만 관찰되지만, 임계값을 초과하면 균열면이 접촉하면서 고조파와 하모닉 성분이 급격히 증폭되는 비선형 현상이 나타난다. 특히, 정상 접촉(N‑contact)은 주로 압축 파동에 의해 유도되며, 전단 접촉(T‑contact)은 전단 변형이 큰 경우에만 활성화된다. 두 접촉 메커니즘이 동시에 작용할 때는 고조파 스펙트럼이 복합적으로 변형되어, 비선형 신호의 민감도가 크게 향상된다.

또한, 경사 균열의 경우 접촉 면적과 접촉 방향이 변하기 때문에, 동일한 진폭에서도 비선형 응답이 수평 균열보다 더 빨리 발생한다는 점을 확인하였다. 이는 실제 구조물에서 균열의 기하학적 특성이 비선형 초음파 신호에 미치는 영향을 정량화하는 데 중요한 통찰을 제공한다.

수치적 정확성을 검증하기 위해 저자들은 시간 스텝과 격자 크기에 대한 수렴 시험을 수행했으며, 결과는 격자 해상도가 충분히 높을 경우 고조파 진폭이 안정적으로 수렴함을 보여준다. 또한, DEM 부분에서 사용된 접촉 법칙(노멀 스프링·댐퍼와 쿠션 마찰 모델)이 비선형 응답에 미치는 영향을 파라미터 스터디를 통해 분석하였다. 마찰 계수가 증가하면 전단 접촉에 의한 고조파가 억제되는 반면, 노멀 스프링 강도가 커질수록 접촉 임계 진폭이 감소한다는 결론을 도출하였다.

이러한 결과는 FDEM이 균열면 접촉·마찰을 물리적으로 정확히 모델링함으로써, 비선형 초음파 NDT에서 관측되는 고조파·하모닉 신호의 원인을 해석하고, 균열의 형태·크기·방향을 역추정하는 데 활용될 수 있음을 시사한다. 특히, 실험적 NDT 데이터와의 정량적 매칭을 위해서는 모델 파라미터(재료 탄성계수, 접촉 스프링·댐퍼, 마찰 계수 등)를 실험적으로 보정하는 절차가 필요하다는 점도 강조한다.