정확한 선탐색만으로 단계적 수렴을 확보하는 위상장 연성 파괴 계산

초록

변분 위상장 모델은 비볼록 에너지 함수를 이용한 지역 제한 최소화 문제를 제시한다. 이산화된 형태에서는 두 미지 변수에 대한 에너지의 별도 볼록성을 활용한 교대 최소화가 일반적으로 사용된다. 각 서브문제가 성공적으로 해결될 경우 이 방법은 이론적으로 수렴이 보장된다. 그러나 에너지 함수의 강한 비선형성은 하나 또는 양쪽 서브문제에서 반복 수렴 실패를 초래할 수 있다. 본 논문에서는 이분법을 기반으로 한 정확한 선탐색 알고리즘을 제안한다. 특정 조건 하에서 이 알고리즘은 각 서브문제에 대한 뉴턴 방법의 전역 수렴을 보장하고, 따라서 교대 최소화 스킴을 통해 에너지의 임계점을 성공적으로 찾을 수 있다. 다양한 변형 에너지 분해와 두·세 차원에서 비가역성 제약을 적용하는 두 가지 전략을 이용한 여러 벤치마크 테스트를 통해 접근법의 견고성을 입증하고, 일반적으로 사용되는 다른 선탐색 알고리즘과의 효율성을 비교 평가한다.

상세 요약

위상장 모델은 파손 메커니즘을 연속적인 스칼라 필드인 손상 변수로 전환함으로써 파괴 현상을 수치적으로 다루는 강력한 프레임워크이다. 이 모델의 핵심은 전체 시스템의 자유 에너지를 손상 변수와 변위 변수에 대한 함수로 정의하고, 이를 최소화함으로써 물리적으로 의미 있는 균형 상태를 찾는 것이다. 에너지 함수는 일반적으로 변형 에너지와 손상 연관 항으로 구성되며, 손상 변수에 대한 비선형성 및 비볼록성이 존재한다. 이러한 특성 때문에 직접적인 전역 최소화는 계산적으로 비현실적이며, 대신 교대 최소화(Alternate Minimization, AM) 방식이 널리 사용된다. AM은 두 변수 집합을 번갈아 고정하고 각각에 대해 볼록 최적화 문제를 푸는 절차로, 각 단계에서 수렴이 보장되는 장점이 있다. 하지만 실제 구현에서는 두 서브문제 모두가 강한 비선형성을 포함하고 있어, 뉴턴‑라프슨 방법 등 2차 최적화 기법을 적용할 때 라인 서치(line search)의 선택이 수렴에 결정적인 영향을 미친다. 부적절한 라인 서치는 단계적 감소가 충분히 이루어지지 않아 뉴턴 단계가 발산하거나, 최소한 허용 오차 내에 도달하지 못하게 만든다.

본 논문이 제시하는 ‘정확한 선탐색(exact line search)’ 알고리즘은 이분법(bisection) 기반으로, 주어진 뉴턴 방향에 대해 에너지 감소를 보장하는 최적의 스텝 길이를 수치적으로 찾는다. 이 방법은 전통적인 백트래킹(line search with Armijo condition)이나 Wolfe 조건 기반 기법과 달리, 스텝 길이 후보를 반복적으로 반으로 나누어 에너지 감소가 확실히 일어나는 구간을 좁혀간다. 중요한 점은, 특정 수학적 가정(예: 에너지 함수의 연속 미분 가능성, 라인 검색 구간 내의 단조 감소성)이 충족될 경우, 이 알고리즘이 뉴턴 방법의 전역 수렴을 보장한다는 이론적 증명이다. 따라서 교대 최소화 전체 흐름에서도 각 서브문제가 반드시 수렴하게 되며, 전체 알고리즘이 안정적으로 임계점(critical point)에 도달한다.

실험 부분에서는 두 가지 주요 변형 에너지 분해—볼츠만-피에르스(Volumetric–Deviatoric) 분해와 분해된 전형적인 에너지 형태—를 적용하고, 비가역성(irreversibility) 제약을 강제하는 두 전략(예: 히스테리시스 연산자와 라그랑주 승수 방식)을 비교하였다. 2차원 및 3차원 베엔치마크(예: 단일 균열 전파, 복합 균열 네트워크)에서 제안된 정확한 선탐색은 기존의 백트래킹이나 Wolfe 기반 라인 서치에 비해 반복 횟수와 총 계산 시간을 현저히 감소시켰으며, 특히 강한 비선형 구간에서의 발산 현상을 완전히 억제하였다. 이는 실제 엔지니어링 시뮬레이션에서 파손 전파를 정확히 예측하고, 계산 비용을 절감하는 데 큰 의미가 있다.

요약하면, 이 논문은 위상장 파괴 해석에서 흔히 발생하는 수렴 문제를 ‘정확한 선탐색’이라는 비교적 간단한 알고리즘으로 해결함으로써, 기존 방법론 대비 이론적 보증과 실용적 효율성을 동시에 제공한다는 점에서 학술적·산업적 가치를 동시에 지닌다.