A robust mixed finite element formulation for third medium contact

기술 요약: 제3매질 접촉을 위한 강건한 혼합 유한 요소 정식화

문제 정의
고전적인 유한 요소 접촉 정식화는 접촉 인터페이스의 명시적 이산화를 통해 페널티, 라그랑주 승수 또는 장벽 방법을 사용하여 제약 조건을 강제하는 데 의존한다. 이러한 방식은 성숙하였으나, 큰 슬라이딩, 자기 접촉(self-contact), 심한 격자 왜곡, 또는 접촉 경계가 사전에 알려지지 않은 위상 최적화 시나리오에서는 상당한 알고리즘적 부담을 초래한다. "제3매질 접촉(Third Medium Contact, TMC)" 방법은 명시적인 제약 조건 대신, 잠재적으로 접촉할 수 있는 물체들 사이의 공간을 채우는 매우 유연한 가상 매질을 도입함으로써 이 문제를 해결한다. 이 매질은 접촉 전에는 부드러운 상태를 유지하다가 압축 시 급격히 단단해지며, 연속체를 통해 접촉력을 전달한다.

그러나 TMC의 핵심 과제는 안정화이다. 매질은 접촉 전에는 스퓨리어스 힘(spurious forces)을 피할 만큼 충분히 부드러워야 하며, 동시에 극심한 압축과 전단에도 격자 왜곡 없이 견딜 수 있을 만큼 충분히 안정적이어야 한다. 기존의 안정화 전략은 종종 2차 변위 미분을 계산하기 위해 고차 유한 요소를 요구하거나, 시스템 크기를 증가시키는 전역 보조 장(auxiliary fields)을 도입하거나, 특수한 이산화를 사용해야 한다. 따라서 2차 변위 미분의 직접적인 평가를 피하면서도 단순하고 계산 효율적인 1차 유한 요소와 호환되는 강건한 안정화 메커니즘이 필요하다.

방법론
저자들은 제3매질을 위한 **보조 장 안정화(auxiliary-field stabilization)**를 도입하는 혼합 유한 요소 정식화를 제안한다. 이 방법의 핵심은 다음과 같다:

1. 보조 장 도입: 제3매질 내에 변형 구배(deformation-gradient)와 유사한 장 $\Theta$를 도입하며, 이를 물리적 변위장 $u$와 분리된 독립적인 미지수로 취급한다.
2. 약한 결합(Weak Coupling): $\Theta$는 페널티 항($W^p_{tm} \propto |\Theta - F|$)을 통해 물리적 변형 구배 $F$와 약하게 결합된다. 이는 $\Theta$가 물리적 변형을 추적하도록 보장한다.
3. 구배 정규화(Gradient Regularization): 물리적 변형 구례의 구배($\nabla F$)를 직접 페널티를 주는 대신(이는 변위의 2차 미분을 요구함), 보조 장의 구배($\nabla \Theta$)를 페펜티를 준다. 정규화 에너지는 $W^r_{tm} \propto |\nabla \Theta|^2$로 정의된다.
4. 혼합 이산화: 정식화는 $u$와 $\Theta$의 독립적인 보간을 허용한다. 저자들은 $\Theta$에 대한 연속(continuous) 및 요소별 불연속(element-wise discontinuous) 근사를 조사하였다.
   • 연속 $\Theta$: 노드 값이 요소 간에 공유되어 요소 간 결합을 생성한다.
   • 불연속 $\Theta$: 자유도가 각 요소에 국부적으로 존재한다. 이들은 **정적 축소(static condensation)**를 통해 제거될 수 있으며, 이를 통해 전역 시스템이 변위 자유도만을 유지하도록 한다.

주요 기여

• 저차 호환성: 이 정식화는 변위의 2차 미분을 계산하지 않고도 구배 유형의 정규화를 달라는 특징을 가지며, 이를 통해 표준 1차(선형) 및 2차 유한 요소와 호환된다.
• 안정화 메커니즘: 논문은 저차 변위 요소(여기서 $F$는 요소별로 일정한 값임)의 경우, $\Theta$의 연속 보간이 필수적임을 밝힌다. 이는 $F$의 요소 간 점프를 $\Theta$의 비영(non-vanishing) 구배로 전달하여 정규화를 활성화한다. 반면, 불연속 $\Theta$ 장은 이러한 요소 간 구배를 생성할 수 없어 저차 요소에 대한 안정화 효과를 약화시킨다.
• 정적 축소: 불연속 보조 장의 경우, 추가되는 미지수들을 국부적으로 축소할 수 있어 전역 시스템 크기를 증가시키지 않으면서도 국부적 안정화를 유지할 수 있음을 입증하였다.

수치 결과
정식화는 다음과 같은 벤치마크 문제들을 통해 평가되었다:

• 이중 블록 압축: 제3매질 스케일링 파라미터 $\gamma$가 잔여 간극(접촉 정확도)과 비선형 수렴 강건성 사이의 트레이드오프를 조절함을 입증하였다. 작은 $\gamma$는 더 타이트한 접촉을 제공하지만 더 많은 뉴턴 반복(Newton iterations)을 요구한다. 이 방법은 압축과 전단이 결합된 상황에서도 안정성을 유지했다.
• 박스 내 자기 접촉:
  • 보간 비교: 연속 보조 장($T_u^1 T_\Theta^1$, $T_u^2 T_\Theta^2$)은 견고한 접촉 폐쇄를 제공하였다. 혼합 쌍인 $T_u^2 T_\Theta^1$(2차 변위, 선형 보조 장)은 $T_u^2 T_\Theta^2$와 유사한 결과를 제공하면서도 더 적은 보조 자유도를 사용하였다.
  • 불연속 vs 연속: 불연속 보조 장은 특히 선형 요소에서 구배 항이 소멸하거나 요소 간 결합에 실패하여 정규화 효과가 감소함을 보여주었다.
  • 변형 구배 평균화와의 비교: Faltus 등[18]의 방법과 비교했을 때, 제안된 혼합 정식화는 훨씬 적은 하중 단계와 뉴턴 반복을 필요로 하여 우수한 비선ular 강건성을 입증하였다. 혼합 정식화가 더 큰 대수 시스템을 가짐에도 불구하고, 개선된 수렴성이 비용을 상쇄하였다. 다만, 매우 미세한 격수($>10^6$ 개의 미지수)에서는 보조 자유도의 비용이 지배적이 되었다.
• C-자형 박스: 정식화는 큰 회전, 점진적인 자기 접촉, 심한 제3매질 압축을 성공적으로 포착하였다. 힘-변위 응답은 고전적인 노드-세그먼트 접촉 정식화와 밀접하게 일치하였으며, 접촉 시 부드러운 전이를 보였다.
• 3D 확장: 3차원 자기 접촉 예제는 이 방법이 평면 변형 설정 너머에도 적용 가능함을 확인시켜 주었다.

의의 및 주장
본 논문은 제안된 보조 장 안정화가 고전적인 접촉 알고리즘 및 기존의 TMC 정규화에 대한 강건한 저차 호환 대안을 제공한다고 주장한다. 변위의 2차 미분을 직접 평가하는 것을 피함으로써, 복잡한 접촉 문제를 포함한 대변형 상황에서 단순하고 계산 비용이 저렴한 유한 요소(선형 및 2차)를 사용할 수 있게 한다.

저자들은 이 방법의 효과가 보조 장의 연속성과 본질적으로 연결되어 있음을 강조한다. 즉, 저차 변위 요소를 위한 구배 기반 안정화 메커니즘을 생성하기 위해서는 연속 보간이 필수적이다. 이 정식화는 접촉 경계가 알려지지 않았거나 진화하는 경우, 예를 들어 접촉 보조 메커니즘, 소프트 로보틱스, 또는 위상 최적화와 같은 분야에 특히 적합하며, 여기서 연속체 기반 접근 방식은 명시적인 접촉 탐색 알고리즘이나 활성 집합(active set) 전략의 필요성을 제거한다. 결론적으로, 이 방법은 더 큰 국부 대수 시스템을 대가로 전역적 강건성을 향상시키지만, 고차 요소나 특수한 격자 기술의 복잡성 없이 가혹한 접촉 시나리오를 시뮬레이션할 수 있는 실행 가능한 경로를 제시한다.