[논문 리뷰] A finite element implementation of surface elasticity at finite strains using the deal.II library

Andrew McBride, Ali Javili외 2명|arXiv (Cornell University)|2015. 06. 02.

Elasticity and Material Modeling참고 문헌 1인용 수 4

한 줄 요약

이 논문은 표면 탄성의 완전 비선형 유한요소 구현을 제안하며, 개방소스 deal.II 라이브러리를 사용하여 유한 변형에서 수행한다. 총 라그랑주 형식을 사용하고 체적 및 표면 영역의 정확한 탄성 기여도를 포함한다. 이 구현은 유한 변형 영역에서 표면 효과를 정확하게 모델링할 수 있게 하며, 벤치마크 검증과 함께 나노스케일 응용을 위한 전체 소스 코드를 제공한다.

ABSTRACT

The potentially significant role of the surface of an elastic body in the overall response of the continuum can be described using the mature theory of surface elasticity. The objective of this contribution is to detail the finite element approximation of the underlying governing equations (both in the volume and on its surface) and their solution using the open-source finite element library deal.II. The fully-nonlinear (geometric and material) setting is considered. The nonlinear problem is solved using a Newton--Raphson procedure wherein the tangent contributions from the volume and surface are computed exactly. The finite element formulation is implemented within the total Lagrangian framework and a Bubnov-Galerkin spatial discretization of the volume and the surface employed. The surface is assumed material. A map between the degrees of freedom on the surface and on the boundary of the volume is used to allocate the contribution from the surface to the global system matrix and residual vector. The deal.II library greatly facilitates the computation of the various surface operators, allowing the numerical implementation to closely match the theory developed in a companion paper. Key features of the theory and the numerical implementation are elucidated using a series of benchmark example problems. The full, documented source code is provided.

연구 동기 및 목표

• 기존 수치적 처리 방식이 무한소 변형에 국한되어 있는 점을 고려해, 유한 변형 영역에서 표면 탄성을 위한 유한요소 구현을 개발하는 것.
• 총 라그랑주 프레임워크와 바뷰노프-갈레르킨 공간 이산화를 사용하여 기하학적 및 재료적 일致성을 갖춘 완전 비선형, 기하학적으로 일致한 형식을 제공하는 것.
• 정확한 탄성 기여도 행렬 계산을 통한 표면 및 체적 기여도 통합을 통해 나노구조물에서의 표면 효과를 정확하게 모델링할 수 있도록 하는 것.
• 공개 소스로 제공되며, 병렬 처리를 위한 TBB 지원을 통해 복잡한 기하구조에서의 성능을 확보한 완전 문서화된 구현을 제공하는 것.
• deal.II의 코드미니멀리티-일(_codimension-one_) 기능을 활용하여 이론적 표면 탄성과 실용적 수치 시뮬레이션을 연결하는 것.

제안 방법

• 유한요소 형식은 총 라그랑주 프레임워크에 기반하며, 체적 및 표면 영역의 공간 이산화에 바뷰노프-갈레르킨 방법을 사용한다.
• 표면 기여는 자유도 매핑을 통해 체적 영역의 경계에 매핑되어 표면 및 체적 자유도 간의 일관된 결합을 가능하게 한다.
• 비선형 시스템은 체적 및 표면 약한 형식에서 정확하게 계산된 탄성 기여도를 포함한 뉴턴-랩소프슨 방법으로 해결된다.
• deal.II 라이브러리의 코드미니멀리티-일 루틴을 사용하여 3차원 공간에 임bed된 2차원 다각형에서의 표면 기하학 및 연산자를 처리한다.
• 표면는 네오-후크형 자유 에너지를 갖는 재료 막으로 모델링되며, 유체적 거동를 위해 표면 장력을 포함하도록 확장된다.
• 공유 메모리 환경에서 TBB 라이브러리를 사용하여 병렬 처리를 구현하였으며, 분산 메모리 확장 지원은 현재 검토 중이다.

실험 결과

연구 질문

• RQ1 유한 변형 영역에서 일관된 유한요소 형식을 사용하여 표면 탄성을 정확하게 구현할 수 있는 방법은 무엇인가?
• RQ2 정확한 탄성 기여도 행렬 계산이 비선형 표면탄성 문제에서 수렴성과 안정성 확보에 미치는 역할은 무엇인가?
• RQ3 deal.II 라이브러리의 코드미니멀리티-일 기능을 활용하여 3차원 연속체에서 복잡한 표면 기하학을 어떻게 모델링할 수 있는가?
• RQ4 메esh 세분화에 따라 표면 대 체적 비율이 증가함에 따라 발생하는 수치적 과제는 무엇이며, 이를 어떻게 완화할 수 있는가?
• RQ5 복잡한 표면 구조를 가진 실질적인 나노구조물의 효율적 시뮬레이션을 위해 이 구현은 어느 정도 병렬화될 수 있는가?

주요 결과

• 정확한 탄성 기여도를 통해 고정밀도로 유한 변형 표면 효과를 성공적으로 캡처하였으며, 표면 장력 및 곡률 의존성 변형을 포함한 벤치마크 문제를 통해 검증되었다.
• 뉴턴-랩소프슨 스킴에서 정확한 탄성 기여도를 사용함으로써 이차 수렴성과 비선형 반복의 안정성이 보장되었다.
• deal.II의 코드미니멀리티-일 기능은 복잡하고 곡률이 있는 표면에서도 정확하고 효율적인 표면 연산자 계산을 가능하게 하였다.
• 완전하고 문서화된 소스 코드가 공개되어 있어 재현성 확보 및 새로운 재료 모델 또는 응용 분야로의 확장이 용이하다.
• TBB를 통한 병렬 처리로 대규모 문제에서 성능 향상이 뚜렷했으며, 분산 환경에서의 추가 스케일링 가능성도 보였다.
• 복잡한 표면 기하학을 가진 실질적인 나노구조물의 시뮬레이션도 가능하여, 학술적 벤치마크를 넘는 실용적 유용성을 입증하였다.