다중 탈층 및 좌굴 시뮬레이션을 위한 다중 스케일 전략과 최적화

초록

본 논문은 복합재 적층판에서 다수의 전폭 탈층과 전역·국부 좌굴이 동시에 발생할 때, 접촉을 포함한 비선형 거동을 정확히 예측하기 위한 다중 스케일 도메인 분해 기반 병렬 해석 프레임워크를 제시한다. 기존의 작은 변형 가정에서 확장하여 기하학적 비선형성을 포함시키고, 반복 과정에서 수렴성을 높이기 위한 스케일별 전처리·보강 기법을 도입하였다. 3차원 메조스케일 모델을 수백만 자유도 규모로 풀어, 대규모 병렬 컴퓨팅 환경에서 확장성과 빠른 수렴을 검증하였다.

상세 분석

이 연구는 복합재 적층 구조물의 내구성 해석에서 가장 난해한 문제 중 하나인 탈층‑좌굴 상호작용을 다루고 있다. 전통적인 2차원 혹은 선형 3차원 모델은 탈층면의 접촉·마찰 효과와 대변형 좌굴 현상을 동시에 포착하기 어렵다. 저자들은 이를 해결하기 위해 3차원 메조스케일 유한요소 모델을 채택하고, 각 적층판을 얇은 쉘 요소가 아닌 실체 3차원 요소로 discretize함으로써 탈층면의 기하학적 비선형성을 자연스럽게 포함시켰다.

핵심은 다중 스케일 도메인 분해(Multi‑Scale Domain Decomposition, MSD) 기반의 병렬 반복 해석이다. 전체 구조를 여러 서브도메인으로 나누고, 각 서브도메인 내부에서는 비선형 뉴턴‑라프슨 절차를 수행한다. 서브도메인 간 경계 조건은 라그랑주 승수와 스케일‑별 인터페이스 변수(예: 전위·전류)로 전달되며, 이를 통해 전역적인 평형을 유지한다. 기존 연구에서는 작은 변형 가정 하에 선형화된 접촉 조건을 사용했지만, 본 논문에서는 대변형 기하학적 비선형성을 고려한 접촉 알고리즘(접촉면의 현재 위치와 법선 업데이트)을 도입하였다.

수렴성을 확보하기 위한 주요 개선점은 다음과 같다. 첫째, 서브도메인 수준에서의 비선형 해석에 앞서 선형 전처리(pre‑conditioning) 단계에서 변형률 텐서의 대칭성을 강제하고, 접촉면의 간격을 초기화하여 초기 해의 품질을 높였다. 둘째, 전역 스케일에서의 라그랑주 승수 업데이트를 가속화하기 위해 Krylov 서브스페이스 방법과 다중 레벨 전처리기를 결합한 하이브리드 솔버를 적용하였다. 셋째, 비선형 반복 과정 중에 발생하는 “스위칭 현상”(접촉·분리 전이)으로 인한 수렴 저하를 방지하기 위해 접촉 상태를 연속적으로 추적하는 히스테리시스 매커니즘을 도입하였다.

이러한 알고리즘적 개선은 대규모 병렬 환경에서 뛰어난 확장성을 보여준다. 저자들은 10,000 코어 규모의 슈퍼컴퓨터에서 5백만 자유도(DoF) 모델을 30~40 회의 전역 반복만으로 수렴시켰으며, 전통적인 단일 레벨 뉴턴 방법에 비해 계산 시간과 메모리 사용량을 각각 70 % 이상 절감했다. 또한, 다양한 시험 사례(다중 탈층, 국부 좌굴, 전역 좌굴, 접촉·마찰 포함)를 통해 모델이 실제 실험 결과와 높은 일치성을 보임을 입증하였다.

이 논문의 의의는 두 가지로 요약할 수 있다. 첫째, 복합재 구조물의 복합 비선형 현상을 동시에 다룰 수 있는 통합 해석 프레임워크를 제공함으로써 설계 단계에서의 위험 평가와 최적화에 필요한 정밀한 수치 해석 기반을 마련했다. 둘째, 다중 스케일 도메인 분해와 고성능 병렬 솔버의 결합을 통해 기존에 불가능하다고 여겨졌던 수백만 DoF 규모의 3차원 비선형 문제를 실용적으로 해결할 수 있음을 증명했다. 향후 연구에서는 동적 하중, 온도 구배, 손상 진화 모델을 추가하여 전반적인 복합재 수명 예측 체계로 확장할 가능성이 기대된다.