Efficient strain-space hyperreduction in large-deformation solid mechanics

기술 요약: 대변형 고체 역학에서의 효율적인 변형률-공간 하이퍼리덕션(Strain-Space Hyperreduction)

문제 정의
변형률 공간(strain space)에서 작동하는 모델 차수 축소(Model Order Reduction, MOR) 기법은 최근 계산적 균질화, 특히 대표 체적 요소(RVE) 분야에서 우수한 성능을 입증했다. 그러나 이러한 기법의 적용은 해당 특정 맥락에 국한되어 왔다. 이를 일반적인 대변형 고체 역학 문제로 확장하는 데는 상당한 어려움이 따른다. 변위 기반 정식화와 달리, 변형률 공간 방식은 경계 조건(Dirichlet BCs)을 직접 강제할 수 없다. 왜냐하면 경계에서의 적분점 변형률을 명시적으로 알 수 없으며, 경계 값은 일반적으로 변형률이 아닌 변위로 규정되기 때문이다. 또한, 표준적인 변형률 공간 접근 방식은 주기적 경계 조건(periodic BCs)에 의존하는 경우가 많아, 임의의 경계 값을 갖는 매크로 스케일 문제에는 적용할 수 없다.

방법론
저자들은 임의의 매개변수화된 디리클레 경계 조건을 가진 문제에 변형률 공간 MOR을 적용하기 위한 일반화된 프레임워크를 제안한다. 핵심 방법론은 세 가지 구성 요소로 이루어진다:

1. 리프팅(Lifting)을 통한 경계 처리: 비제로(non-zero) 디리클레 경계 조건을 만족시키기 위해, 저자들은 각 독립 경계에 대한 단위 변위 값에 대해 오프라인으로 일관된 참조 필드(reference field)를 계산하는 리프팅 전략을 사용한다. 실제 해는 이 참조 필드(실제 경계 값에 의해 스케일링됨)와 "변동 필드(fluctuation field)"로 분해된다. 변동 필드는 구조적으로 모든 디리클레 경계에서 0이 되도록 설계되었으므로, 고유 직교 분해(POD)로부터 유도된 축소 기저(reduced basis)가 경계 조건을 사전적으로(a priori) 만족할 수 있게 한다. 이 접근 방식은 비균질 경계 조건을 축소 솔버를 위한 균질 경계 조건으로 변환한다.

2. 변형률-공간 정식화: 본 논문은 세 가지 특정 하이퍼리덕션 기법을 이 새로운 환경에 맞게 일반화한다:

   • 변형률-공간 ECM (Empirical Cubature Method): 요소 수준의 연산을 우회하여 가우스 점(Gauss points)에서 직접 축소 적분을 수행하도록 재구성되었다.
   • 변형률-공간 E3C (Empirically Corrected Cluster Cubature): 임의의 역학 문제로 일반화되었다. E3C는 변형률 공간에서 적분점과 가중치를 최적화하며, 이를 통해 원래 메쉬의 가우스 점과는 다른 지점을 선택할 수 있게 한다.
   • 변형률-공간 EMSL (Empirical Material Sampling and Linearisation): RVE를 넘어 확장되었다. EMSL은 클러스터별 기대 평균 변형률 주변에서 재료 거동을 선형화한다. 정확한 재료 법칙을 근사적으로 적분하는 다른 방법들과 달리, EMSL은 클러스터 단위의 선형 근사를 정확하게 적분하므로, 온라인 단계에서 반복적인 뉴턴-랩슨(Newton-Raphson) 솔버를 필요로 하지 않는 로드 스텝당 아핀(affine) 문제를 생성한다.

3. 비교 기준선: 이러한 변형률-공간 방법들은 구조 보존적 베이스라인 역할을 하는 표준 변위-공간 에너지 보존 샘플링 및 가중치(ECSW) 방법과 벤치마킹된다.

주요 기여

• 변형률-공간 MOR의 일반화: 본 연구는 변형률 공간 MOR 기법(ECM, E3C, EMSL)을 계산적 균질화에서 일반적인 대변형 고체 역학 및 임의의 디리클레 경계 조건으로 성공적으로 확장하였다.
• 새로운 경계 처리: 리프팅 기반의 변동 필드 정식화를 도입하여, 경계 변형률 데이터 없이도 변형률 공간에서 임의의 디리클레 경계 조건을 부과할 수 있게 하였다.
• 포괄적인 벤치마킹: 저자들은 역 매개변수 추정, 가변 재료 매개변수, 그리고 복잡한 하중 경로(인장, 전단 및 혼합 하중)를 포함하는 두 가지 초탄성 예제를 통해 변형률-공간 방식과 변위-공간 ECSW 방식 간의 엄격한 비교를 제공한다.

결과
방법론은 타원형 구멍이 있는 두 가지 초탄성 판 문제에 대해 테스트되었다:

1. 예제 1 (중등도 변형): 최대 20% 변형률을 갖는 인장 하중 하의 판.
2. 예제 2 (대변형 및 복잡한 하중): 동일한 기하학적 구조에 인장(최대 50%), 전단(최대 25%) 및 혼합 하중을 가한 경우로, 더 미세한 메쉬와 높은 계산 비용을 수반한다.

성능 결과:

• 정확도 대비 실행 시간: 두 예제 모두에서 변형률-공간 방법들은 변위-공간 ECSW에 비해 정확도와 실행 시간 사이의 트레이드오프 측면에서 일관되게 우수한 성능을 보였다.
• 가속도(Speedups):
  • E3C: 충분한 적분점이 사용되었을 때 높은 정확도(예: 오차 <1%)를 유지하면서, 전체 모델(full-order model) 대비 약 10,000배(예제 1) 및 30,000배(예제 2)의 속도 향상을 달야냈다. E3C는 주어진 계산 예산 내에서 많은 적분점을 사용할 경우 변형률-공간 방법들 중 일반적으로 가장 높은 정확도를 제공했다.
  • EMSL: 탁월한 실행 시간 효율성을 입증하여, 최대 100,000배(예제 1) 및 300,000배(예제 2)의 가속도를 달성했다. EMSL은 온라인 및 오프라인 예산이 매우 제한적일 때 뛰어난 성능을 보였으며, 종종 E3C보다 속도 면에서 우수했으나, 고비선형 영역에서는 많은 적분점을 가진 E3C에 비해 약간 낮은 정확도를 보이기도 했다.
  • ECM: ECSW보다는 성능이 좋았으나, 테스트된 특정 지표에서는 E3C 및 EMSL에 비해 뒤처지는 경향을 보였다.
• 오프라인 비용: E3C는 반복적인 최적화로 인해 더 비싼 오프라인 훈련 단계를 갖지만, 저자들은 이것이 관리 가능한 수준이며 훈련 데이터를 생성하는 데 걸리는 시간보다 현저히 적다고 언급했다. EMSL은 가장 낮은 오프라인 비용을 갖는다.

의의 및 주장
본 논문은 제안된 방법론이 계산적 균질화를 넘어 변형률-공간 MOR 연구를 정당화한다고 주장한다. 저자들은 이러한 기법들이 다음과 같은 이점을 제공한다고 강조한다:

• 우수한 효율성: 높은 정확도를 유지하면서도 온라인 실행 시간을 극적으로 단축한다(수십 배에서 수만 배).
• 비침투성(Non-Intrusiveness): 이 방법들은 변위, 변형률, 응력 데이터(블랙박스 솔버에서 쉽게 추출 가능)와 재료 루틴만을 필요로 한다. 요소 수준의 연산을 우회하며, 이산화 상세 내용(예: 요소 유형)이나 잔차 벡터(residual vectors)에 대한 지식을 요구하지 않는다.
• 적용 가능성: 임의의 디리클able 경계 조건을 처리하는 데 성공함으로써, 이러한 효율적인 기법들을 매크로 스케일 구조 문제에 적용하는 데 있어 주요 장벽을 제거하였다.

저자들은 고비선형 문제(예: 탄소성)와 EMSL에서의 클러스터별 변형률 추정 등의 과제가 여전히 남아 있지만, 입증된 성능은 이러한 방법들이 현재 계산 비용이 너무 높아 실행이 어려운 제어, 최적화 및 재료 모델 피팅 분야에서 새로운 공학적 응용을 가능하게 할 것임을 시사한다고 결론지었다.