[논문 리뷰] A Mixed Finite Element Method for Multi-Cavity Computation in Incompressible Nonlinear Elasticity

Weijie Huang, Zhiping Li|arXiv (Cornell University)|2018. 01. 08.

Advanced Numerical Methods in Computational Mathematics인용 수 3

한 줄 요약

이 논문은 비압축성 비선형 탄성에서 다중 캐비티 성장을 계산하기 위해 $Q_2$-$P_1$ 및 $P_2^+$-$P_1$ 등가요소를 조합한 혼합 유한요소법을 제안한다. 이 방법은 이론적으로 잠김이 없는 것으로 증명되고 수렴성이 입증되었으며, 수치 실험을 통해 하중 감소 시 비대칭에서 대칭 구성을 향한 전이를 보이는 새로운 분기 현상이 발견되었다. 이 현상은 임계 변형률 임계값 $ olambda_c ∈ (1.175, 1.180)$ 에서 발생한다.

ABSTRACT

A mixed finite element method combining an iso-parametric $Q_2$-$P_1$ element and an iso-parametric $P_2^+$-$P_1$ element is developed for the computation of multiple cavities in incompressible nonlinear elasticity. The method is analytically proved to be locking-free and convergent, and it is also shown to be numerically accurate and efficient by numerical experiments. Furthermore, the newly developed accurate method enables us to find an interesting new bifurcation phenomenon in multi-cavity growth.

연구 동기 및 목표

• 비압축성 비선형 탄성 재료에서 큰 변형과 비압축성 제약 조건으로 인해 발생하는 도전 과제들에 대비해 다중 캐비티 발생 현상을 시뮬레이션할 수 있는 견고한 수치적 방법을 개발하는 것.
• 특히 캐비티 형성 부근에서 흔히 발생하는 잠김 효과를 극복하기 위해 표준 유한요소법을 비압축성 탄성에 적용할 때 발생하는 문제를 해결하는 것.
• 다중 캐비티 문제에 대해 수치적 안정성과 수렴성을 보장하는 혼합 유한요소 형식을 수립하는 것.
• 다양한 하중 조건 하에서 다중 캐비티 해의 분기 행동을 수치적으로 조사하여, 특히 축대칭과 비축대칭 캐비티 구성을 전환하는 현상을 밝혀내는 것.
• 제안된 방법이 복잡한 변형 패tern, 특히 에너지 기반의 캐비티 형태 선택을 정확히 포착할 수 있음을 이론적 및 수치적 증거로 제시하는 것.

제안 방법

• 이 방법은 두 가지 다른 등가요소를 사용하는 혼합 유한요소 형식을 적용한다: 변위와 압력에 대해 $Q_2$-$P_1$를, 변위장의 향상된 근사에 대해 $P_2^+$-$P_1$를 사용한다.
• 혼합 형식은 비압축성 제약 조건 $\det \nabla \mathbf{u} = 1$ 을 약한 형태로 강제하기 위해 압력 유사 라그랑주 승수 $p$ 를 도입하여 강한 강제가 필요 없도록 한다.
• 비선형 오일러-라그랑즈 방정식은 체계의 안정성과 수렴성을 확보하기 위해 사다리꼴 뉴턴 방법을 사용하여 변분 문제의 이차형식에서 유도된다.
• 이전의 이중 매개변수 혼합 요소에 대한 기법을 확장하여 이론적으로 잠김이 없는 것으로 증명되었으며, 적절한 정규성 조건 하에서 최적 수렴성을 보장한다.
• 저장 에너지 밀도 $W(F) = \mu|F|^s + d(\det F)$ 를 사용하며, 비압축성 조건 하에서 상수이지만 수렴성 증명에서 핵심적인 역할을 한다.
• 변형률과 행렬식의 $L^2$ 및 $L^1$ 오차 노름을 통해 수치적 수렴성을 검증하였으며, 최적 수렴 속도를 확인하였다.

실험 결과

연구 질문

• RQ1 다중 캐비티 해가 비대칭에서 대칭 구성을 전환하는 임계 하중 매개변수 $\lambda_c$ 는 무엇인가?
• RQ2 축대칭과 비축대칭 캐비티 해 간의 에너지 차이 $E_s - E_r$ 는 변형이 증가함에 따라 어떻게 변화하는가?
• RQ3 비압축성 비선형 탄성 하에서 다중 캐비티 시스템에서 한 캐비티가 다른 캐비티보다 더 크게 성장하는 데 영향을 주는 요소는 무엇인가?
• RQ4 고도로 비대칭적인 변형을 겪는 다중 캐비티 문제에 대해 혼합 유한요소법이 잠김이 없고 수렴성 있는가?
• RQ5 초기 조건은 다중 캐비티 성장에서 최종 캐비티 형태에 어떤 영향을 미치며, 에너지 최소화는 안정된 구성을 선택하는 데 어떤 역할을 하는가?

주요 결과

• 제안된 혼합 유한요소 방법은 이론적으로 잠김이 없고 수렴성이 입증되었으며, $L^2$ 및 $L^1$ 오차에 대해 수치적으로 최적 수렴 속도를 확보하였다.
• 수치 실험 결과, $\lambda > 1.18$ 일 경우 축대칭 해와 오른쪽 우세 해 간의 에너지 차이 $E_s - E_r$ 가 급격히 증가하여 비대칭 구성을 선호하는 경향을 보였다.
• 비축대칭 해에서 두 캐비티의 부피 비율 $v_2/v_1$ 는 $\lambda > 1.18$ 일 경우 급격히 증가하여 캐비티 성장의 강한 비대칭성을 나타냈다.
• 임계 전환점 $\lambda_c$ 는 $(1.175, 1.180)$ 의 구간에 위치하며, 하중 감소에 따라 비대칭 해를 선호하는 시스템에서 대칭적이고 잘 분리된 구성을 향해 전환된다.
• 이 방법은 이론적 예측과 일치하는 에너지 기반의 안정된 구성을 자동으로 선택하는 것을 성공적으로 포착하였다. Henao & Serfaty 의 다중 캐비티 거동 이론적 예측과 일치한다.
• 새로운 분기 현상이 발견되었다: 하중 매개변수 $\lambda$ 가 $\lambda_c$ 를 감소하면서 지나갈 때, 다수의 안정된 형태(비대칭 또는 융합된 캐비티)에서 대칭적이고 잘 분리된 구조로의 전환이 발생한다.