[논문 리뷰] A consistent mixed-dimensional coupling approach for 1D Cosserat beams and 2D surfaces in 3D space

Ivo Steinbrecher, Nora Hagmeyer외 2명|arXiv (Cornell University)|2022. 10. 28.

Fluid Dynamics and Vibration Analysis인용 수 4

한 줄 요약

이 논문은 3D 고체의 2D 표면과 1D 코스세라 빔을 일관된 모르타르 기반 유한요소 방법으로 결합하는 방법을 제안하며, 일관된 표면 법선장과 회전 텐서 표현을 통해 위치 및 회전 연속성을 강제한다. 이 방법은 각운동량 보존을 보장하고 복잡한 빔-표면 상호작용을 정확하고 견고하게 시뮬레이션할 수 있으며, 스텐트-유체-구조 상호작용 예제를 통해 법선 응력 최대 1.82 kN/m² 및 비틀림 응력 최대 1.29 kN/m² 예측을 통해 성능을 입증하였다.

ABSTRACT

The present article proposes a novel computational method for coupling arbitrarily curved 1D fibers with a 2D surface as defined, e.g., by the 2D surfaces of a 3D solid body or by 2D shell formulations. The fibers are modeled as 1D Cosserat continua (beams) with six local degrees of freedom, three positional and three rotational ones. A kinematically consistent 1D-2D coupling scheme for this problem type is proposed considering the positional and rotational degrees of freedom along the beams. The positional degrees of freedom are coupled by enforcing a constant normal distance between a point on the beam centerline and a corresponding point on the surface. This strategy requires a consistent description of the surface normal vector field to guarantee fundamental mechanical properties such as conservation of angular momentum. Coupling of the rotational degrees of freedom of the beams and a suitable rotation tensor representing the local orientation within a solid volume has been considered in a previous contribution. In the present work, this coupling approach will be extended by constructing rotation tensors that are representative of local surface orientations. Several numerical examples demonstrate the consistency, robustness and accuracy of the proposed method. To showcase its applicability to multi-physics systems of practical relevance, the fluid-structure interaction example of a vascular stent is presented.

연구 동기 및 목표

• 3D 공간 내 1D 코스세라 빔과 2D 고체 표면 간 일관된 혼합 차원 결합 접근법을 개발하는 것.
• 빔-표면 인터페이스에서 위치 및 회전 자유도 모두에 대해 운동학적 일관성을 확보하는 것.
• 결합 형식에서 각운동량 보존과 같은 기본 기계적 성질을 유지하는 것.
• 유연한 동맥 벽과 상호작용하는 혈관 스텐트와 같은 복잡한 빔 시스템의 정확한 시뮬레이션을 가능하게 하는 것.
• 기존의 빔-고체 부피 결합 방법을 진정한 1D-2D 결합으로 확장하여 위치 및 회전 제약 조건을 모두 포함하는 것.

제안 방법

• 빔 중심선과 고체 표면 간의 결합 제약 조건을 약한 형태로 강제하기 위해 모르타르 유형의 유한요소 방법을 사용한다.
• 위치 결합은 빔 중심선 점과 고체 표면 사이의 일정한 법선 거리를 유지함으로써 강제되며, 이는 일관된 표면 법선 벡터장이 필요하다.
• 운동 자유도는 국소 표면 방향을 나타내는 표면 삼중체 장을 통해 구성된 회전 텐서를 사용하여 결합된다.
• 제약 조건을 강제하기 위해 라그랑주 승수를 사용하고, 이산 시스템의 안정성을 확보하기 위해 가중치가 부여된 페널티 정규화를 적용한다.
• 빔 모델은 각 점에 대해 세 가지 이동도와 세 가지 회전 자유도를 가진 기하학적으로 정확한 코스세라 빔 이론에 기반한다.
• 기존의 1D-3D 빔-고체 부피 방법에서 유도된 결합 형식을 진정한 1D-2D 빔-표면 결합으로 확장하여, 내장된 빔을 가진 박판 구조의 효율적이고 정확한 모델링을 가능하게 한다.

실험 결과

연구 질문

• RQ1 빔 중심선과 고체 표면 사이의 일정한 법선 거리 강제가 1D-2D 빔-표면 결합에서 각운동량 보존을 보장하는가?
• RQ2 표면 법선 장을 일관되게 포함할 경우, 표준 메쉬 바인딩 기법에 비해 위치 결합의 정확도와 안정성은 어떻게 향상되는가?
• RQ3 1D-2D 문제에서 회전 결합을 위한 국소 표면 방향을 표현하기 위해 일관된 표면 삼중체 장을 구성할 수 있는가?
• RQ4 위치 및 회전 결합의 병합이 혈관 스텐트와 같은 복잡한 빔 시스템의 기계적 거동에 미치는 영향은 무엇인가?
• RQ5 제안된 방법이 환자 특화 혈관 스텐트 모델에서 전역적인 유체-구조 상호작용 효과를 정확히 포착할 수 있는가?

주요 결과

• 오직 일관된 위치 결합 변형만이 —현재 표면 법선 장의 평가가 필요함— 각운동량 보존을 보장하고 물리적으로 올바른 결과를 도출한다.
• 특히 큰 변형 하에서, 일관된 방법은 표준 메쉬 바인딩 기법에 비해 기계적 일관성과 견고성 면에서 뛰어나다.
• 제안된 방법은 스텐트 빔과 동맥 벽 사이의 상호작용 선력(선형 응력)을 성공적으로 포착하였으며, 최대 법선 응력은 1.8179·10³ N/m², 최대 비틀림 응력은 1.28899·10³ N/m²로 추정되었다.
• 수치 예제들은 복잡한 빔 기하구조, 특히 스텐트 유사 구조에 대해 이 방법의 견고성과 정확성을 확인하였다.
• 다물리적 유체-구조 상호작용 예제는 이 방법이 탄성 변화와 혈류역학에 미치는 영향을 모델링할 수 있음을 보여주었다.
• 결합 방법은 인터페이스 힘의 신뢰성 있는 정량화를 가능하게 하여, 스텐트 성능과 재협착 위험 예측의 장기적 분석을 지원한다.