[논문 리뷰] An Eulerian approach to the simulation of deformable solids: Application to finite-strain elasticity

Ken Kamrin, Jean‐Christophe Nave|ArXiv.org|2009. 01. 24.

Elasticity and Material Modeling참고 문헌 12인용 수 11

한 줄 요약

이 논문은 초탄성 고체에서 대변형을 시뮬레이션하기 위한 새로운 오일러 기반 수치 방법을 제안한다. 이 방법은 물질 입자의 초기 위치에서 현재 위치로의 이동을 추적하는 기준 지도(field)를 사용하며, 라그랑주 메시 없이도 정확하고 2차 정확도를 갖는 정적 및 동적 유한변형 탄성 시뮬레이션을 가능하게 한다. 분석적 해와의 검증 결과 수렴 속도는 1.97~1.99를 기록하였다.

ABSTRACT

We develop a computational method based on an Eulerian field called the "reference map", which relates the current location of a material point to its initial. The reference map can be discretized to permit finite-difference simulation of large solid-like deformations, in any dimension, of the type that might otherwise require the finite element method. The generality of the method stems from its ability to easily compute kinematic quantities essential to solid mechanics, that are elusive outside of Lagrangian frame. This introductory work focuses on large-strain, hyperelastic materials. After a brief review of hyperelasticity, a discretization of the method is presented and some non-trivial elastic deformations are simulated, both static and dynamic. The method's accuracy is directly verified against known analytical solutions.

연구 동기 및 목표

• 라그랑주 유한요소에 의존하지 않고 고체 유사 물질의 대변형 변형을 시뮬레이션하기 위한 완전한 오일러 프레임워크를 개발하는 것.
• 고체의 물성 법칙을 오일러 기반 속도 형태로 변환하는 데서 발생하는 수학적 복잡성을 해결하기 위해 운동학적 정밀도를 유지하는 기준 지도를 도입하는 것.
• 오일러 설정에서 변형률 계열과 고체역학의 핵심 양들을 정확하게 계산할 수 있도록 하여 이전 방법에서 흔히 발생하는 근사치를 피하는 것.
• 비틀림 없는 정적 및 동적 탄성 변형 문제, 특히 최대 36% 변형률을 가진 진행하는 압축파를 포함한 비정질 문제에 대해 방법을 검증하는 것.
• 탄성-플라스틱성 및 유체-고체 상호작용과 같은 더 복잡한 재료 모델로의 확장을 위한 기반을 마련하는 것.

제안 방법

• 이 방법은 현재 공간 위치 $\mathbf{x}$ 에서 초기 물질 좌표 $\mathbf{X}$ 로 되돌리는 기준 지도 필드 $\boldsymbol{\xi}(\mathbf{x},t)$ 를 사용하여 물질의 역사적 경로를 오일러 방식으로 추적한다.
• 변형률 계열 $\mathbf{F}$ 는 기준 지도의 기울기에서 유도되며, $\mathbf{F} = \nabla_{\mathbf{x}} \boldsymbol{\xi}^{-1}$ 로 표현되어 고체역학에서 필수적인 운동학적 양에 직접 접근할 수 있다.
• 기준 지도의 진화는 속도장에서 유도된 운반 방정식 $\partial_t \boldsymbol{\xi} + \mathbf{v} \cdot \nabla \boldsymbol{\xi} = 0$ 으로 제어되며, 이는 유체 유사한 이류와 일관성을 유지한다.
• 공간 및 시간 이산화에 두 번째 차수 정확도를 갖는 유한차분 스킴을 사용하고, 안정성과 정확도를 유지하기 위해 룬게-쿠타 시간 적분기를 적용한다.
• 기본적으로 초탄성성과 같은 물성 법칙은 기준 지도에서 유도된 변형률 계열을 사용하여 적용되며, 오일러 기준에서 응력 계산이 가능해진다.
• 경계 조건은 기준 지도를 통해 적용되며, 압력이 없는 유체 상과의 상호작용을 통해 자유 표면 조건을 구현하는 것이 제안된다.

실험 결과

연구 질문

• RQ1 라그랑주 메시 없이도 근사치 없이 순수 오일러 프레임에서 대변형 탄성 변형을 정확하게 시뮬레이션할 수 있는가?
• RQ2 고체 유사 물질에 대해 오일러 설정에서 변형률 계열과 같은 운동학적 양을 신뢰성 있게 계산할 수 있는 방법은 무엇인가?
• RQ3 제안된 방법이 정적 및 동적 변형에 대해 분석적 해와 비교했을 때 수렴 속도와 정확도는 어느 정도인가?
• RQ4 기준 지도 설정은 최대 36% 변형률을 가진 진행하는 압축파와 같은 복잡한 대변형 문제를 처리할 수 있는가?
• RQ5 이 방법은 탄성-플라스틱성 및 유체-고체 상호작용과 같은 더 복잡한 물성 거동으로 확장될 수 있는가?

주요 결과

• 이 방법은 이격된 두 최소 그리드 간에 $\xi$ 에 대해 1.97, 속도에 대해 1.99의 수렴 속도를 보이며, 변위 및 속도 필드 모두에서 전역적으로 2차 정확도를 달성한다.
• 진행하는 압축파의 수치적 해는 조밀한 격자조차도 매우 정확하게 분석적 해와 일치하며, 이는 방법의 강건성과 정확도를 확인한다.
• 파동 펄스 내 최대 압축 변형률은 약 36%에 이르며, 이는 방법이 대규모 유한변형 변형을 처리할 수 있음을 보여준다.
• 기준 지도를 통해 오일러 설정에서 변형률 계열과 기타 운동학적 변수를 정확히 계산할 수 있으며, 물성 법칙 설정 시 일반적으로 발생하는 근사치가 필요 없어진다.
• 이 방법은 정적 및 동적 탄성 변형을 매우 높은 정확도로 시뮬레이션하며, 알려진 분석적 해와의 검증을 통과하였다.
• 초기 결과에 따르면, 이 방법은 향후 탄성-플라스틱성 및 유체-고체 상호작용 시뮬레이션으로의 확장에 매우 유망한 것으로 보인다.