[논문 리뷰] Stability conditions for the explicit integration of projection based nonlinear reduced-order and hyper reduced structural mechanics finite element models

C. Bach, Linze Song외 2명|arXiv (Cornell University)|2018. 06. 29.

Model Reduction and Neural Networks참고 문헌 41인용 수 5

한 줄 요약

이 논문은 구조역학에서 갈레르킨 투영 기반 비선형 축소차수모델(ROM)이 명시적 중심차분 시간적분을 사용할 경우, 전체차수모델에 비해 더 큰 또는 동일한 안정적인 시간스텝 크기를 가지는 엄밀한 이론적 분석을 제공한다. 핵심 기여는 기저 절단에 의한 고유값 분리가 수치적 안정성을 향상시켜 더 큰 시간스텝과 더 빠른 시뮬레이션을 가능하게 한다는 수학적 증명이다. 이는 안정성을 훼손하지 않고도 성능 향상을 이룬다.

ABSTRACT

Projection-based nonlinear model order reduction methods can be used to reduce simulation times for the solution of many PDE-constrained problems. It has been observed in literature that such nonlinear reduced-order models (ROMs) based on Galerkin projection sometimes exhibit much larger stable time step sizes than their unreduced counterparts. This work provides a detailed theoretical analysis of this phenomenon for structural mechanics. We first show that many desirable system matrix properties are preserved by the Galerkin projection. Next, we prove that the eigenvalues of the linearized Galerkin reduced-order system separate the eigenvalues of the linearized original system. Assuming non-negative Rayleigh damping and a time integration using the popular central difference method, we further prove that the theoretical linear stability time step of the ROM is in fact always larger than or equal to the critical time step of its corresponding full-order model. We also give mathematical expressions for computing the stable time step size. Finally, we show that under certain conditions this increase in the stability time step even extends to some hyper-reduction methods. The findings can be used to compute numerical stability time step sizes for the integration of nonlinear ROMs in structural mechanics, and to speed up simulations by permitting the use of larger time steps.

연구 동기 및 목표

• 명시적 적분을 사용할 경우 축소차수모델이 전체차수모델보다 더 큰 안정적인 시간스텝을 가지는 것으로 관측된 현상을 설명하기 위해.
• 명시적 시간적분 하에서 갈레르킨 투영 기반 비선형 ROM의 수치적 안정성을 엄밀히 분석하기 위해.
• ROM이 전체차수모델의 안정성과 동일하거나 이를 초월하는 수학적 조건을 설정하기 위해.
• 이 안정성 향상 효과가 DEIM 및 ECSW와 같은 초차수화 기법으로도 확장되는지 조사하기 위해.
• 축소차수 기저(ROB)가 시스템 고유값과 안정성 시간스텝에 어떤 영향을 미치는지 정량화하기 위해.

제안 방법

• 갈레르킨 투영 기반 ROM의 선형화된 연립방정식을 유도하고, 이를 통해 유도된 행렬의 고유값 구조를 분석한다.
• 선형화된 ROM 시스템의 고유값이 전체차수모델의 고유값과 교차함으로써 안정성 향상이 보장됨을 증명한다.
• 중앙차분법을 시간적분에 적용하고 고유값의 범위에 기반한 이론적 안정성 조건을 유도한다.
• 비음성 레일라이즈 감쇠 및 M-직교 또는 반직교 축소기저가 안정성에 미치는 영향을 분석한다.
• DEIM 및 ECSW와 같은 초차수화 기법이 대칭성, 정정규성, 고유값 교차성의 유지 여부를 평가한다.
• 행렬 이론과 스펙트럼 분석을 활용하여 ROM에서 안정적인 시간스텝 크기를 위한 닫힌 형태의 수식을 유도한다.

실험 결과

연구 질문

• RQ1 왜 갈레르킨 투영 기반 ROM이 명시적 적분을 사용할 경우, 전체차수모델보다 더 큰 안정적인 시간스텝을 허용하는가?
• RQ2 선형 안정성 시간스텝이 ROM에서 전체차수모델보다 크거나 같아지는 수학적 조건은 무엇인가?
• RQ3 DEIM 또는 ECSW를 사용한 초차수화 모델에서도 안정성 향상 효과가 어느 정도 유지되는가?
• RQ4 축소차수 기저(ROB)의 선택이 고유값 분포와 ROM의 안정성에 어떤 영향을 미치는가?
• RQ5 시스템 행렬의 어떤 구조적 성질(예: 대칭성, 정정규성)이 ROM에서 안정성 유지 또는 향상에 필수적인가?

주요 결과

• 선형화된 갈레르킨 ROM 시스템의 고유값은 전체차수모델의 고유값과 교차하며, 이로 인해 ROM의 최대 고유값이 전체모델의 최대 고유값보다 크지 않음을 보장한다.
• 비음성 레일라이즈 감쇠 및 반직교 또는 M-직교 축소기저 조건 하에서, 이론적 선형 안정성 시간스텝은 항상 전체차수모델의 것보다 크거나 같게 된다.
• 안정성 향상의 근본 원인은 기저의 적절한 직교 절단에 의한 고유값 분리이다.
• DEIM를 통한 초차수화는 고유값 교차성과 행렬 대칭성을 유지하지 못해 안정성 저하가 발생할 수 있다.
• ECSW는 행렬 대칭성과 정반정규성을 유지하지만 고유값 교차성은 반드시 보장되지 않아 중간 수준의 안정성 행동를 보인다.
• 이 발견들은 온라인 ROM 시뮬레이션에서 더 큰 시간스텝 사용을 가능하게 하여 안정성을 훼손하지 않고도 계산 성능을 직접적으로 향상시킬 수 있다.