마찰 접촉 구조물의 형태 및 위상 최적화를 위한 위상장 방법

초록

본 논문은 탄성체와 강체 사이의 마찰 접촉을 기술하는 비선형·비볼록 헤미변분 부등식 모델을 기반으로, 에너지형 목적함수에 대한 형태 미분을 유도하고 정규화 기법을 통해 일반 목적함수의 감도 분석을 수행한다. 이를 토대로 경계 변분법을 이용한 형태 최적화와, Allen‑Cahn 기반 1차·2차 정규화 위상장, 그리고 위상 도함수를 결합한 세 가지 위상 최적화 알고리즘을 제안한다. 수치 실험을 통해 제안 방법들의 정확성과 효율성을 검증한다.

상세 분석

이 연구는 마찰 접촉 문제를 기술하는 헤미변분 부등식(hemivariational inequality, HVI)의 수학적 특성을 깊이 파고든다. HVI는 접촉면에서 비스무스하고 비볼록한 마찰 퍼텐셜 jτ 을 서브다이버전스로 표현함으로써 전통적인 변분 부등식보다 실제 마찰 현상을 더 정밀히 모델링한다. 그러나 비볼록성 때문에 존재·유일성 증명과 감도 분석이 어려워, 저자들은 Clarke의 일반화된 방향미분과 서브다이버전트를 활용해 물질 미분(material derivative)과 Eulerian 파생을 엄밀히 정의한다. 특히 에너지형 목적함수에 대해서는 정규화 없이도 직접적인 Eulerian 미분을 도출하고, 일반 목적함수에 대해서는 비스무스 항을 매끄러운 근사함수로 정규화한 뒤, 정규화 파라미터가 0으로 갈 때의 극한을 통해 원문 문제와의 일치성을 asymptotic analysis 로 증명한다.

감도 분석 결과는 경계 변분법(boundary variational method)으로 형태 최적화 문제를 풀 수 있게 해준다. 여기서는 변형된 영역의 경계 변분을 직접 계산하고, Newton 방법을 이용해 비선형 방정식을 효율적으로 해결한다.

위상 최적화 측면에서는 세 가지 위상장 알고리즘을 제시한다. 첫 번째는 Allen‑Cahn 방정식에 기반한 1차 흐름으로, 목적함수의 에너지 형태를 L²‑gradient flow 로 전개한다. 두 번째는 목적함수 자체에 2차 정규화(term ∥∇φ∥²)를 추가해 보다 부드러운 인터페이스를 유지하면서도 수치적 안정성을 확보한다. 세 번째는 기존 위상장에 위상 도함수(topological derivative)를 결합해, 초기 설계에서 작은 구멍을 삽입하거나 제거하는 형태의 급격한 위상 변화를 효율적으로 탐색한다. 각 알고리즘은 HVI의 정규화된 형태에 적용되며, Newton‑Krylov 스키마와 적응형 메쉬 재구성을 통해 고차원 3D 문제까지 확장 가능하도록 설계되었다.

수치 실험에서는 2D 및 3D 접촉 구조물에 대해 부피 제약 하의 최소 에너지 설계, 마찰 계수 변동에 따른 최적 형상 변화, 그리고 위상 도함수를 이용한 구멍 생성/소멸 사례를 제시한다. 결과는 제안된 방법들이 기존 변분 기반 최적화에 비해 더 정확한 마찰 반응을 재현하고, 위상 변화가 자유로운 설계 공간에서도 수렴성이 뛰어남을 보여준다. 전반적으로 이 논문은 HVI라는 난해한 수학적 모델을 실제 엔지니어링 설계에 적용할 수 있는 체계적인 이론·알고리즘 프레임워크를 제공한다는 점에서 큰 의의를 가진다.