분할 표면 기반 이방성 성장 쉘 요소

초록

본 논문은 Loop 서브디비전 표면을 이용한 얇은 쉘 유한요소 모델을 제시한다. Kirchhoff‑Love 이론에 이방성 평면 성장 텐서를 도입해 큰 변형과 성장 현상을 동시에 해석할 수 있다. 여러 정형 로드 테스트와 성장 시뮬레이션을 통해 계산 효율성과 적용 범위를 검증한다.

상세 분석

이 연구는 얇은 쉘 구조 해석에 있어 두 가지 핵심 기술을 결합한다. 첫째, Loop 서브디비전 표면을 기반으로 한 쉘 요소는 C² 연속성을 자연스럽게 제공한다는 점에서 전통적인 4‑노드 혹은 6‑노드 쉘 요소보다 기하학적 정확도가 뛰어나다. 서브디비전은 정점 주변의 인접 삼각형을 재귀적으로 세분화하면서 새로운 정점을 삽입하고, 기존 정점의 위치를 가중 평균한다. 이 과정에서 얻어지는 매끄러운 곡면은 Kirchhoff‑Love 이론이 요구하는 중간면의 연속성을 완벽히 만족한다.

둘째, 성장 이론을 쉘 모델에 통합하기 위해서는 변형 텐서를 두 부분으로 분해하는 multiplicative decomposition, F = Fᵉ · Fᵍ 를 채택한다. 여기서 Fᵍ는 순수한 이방성 평면 성장 텐서이며, Fᵉ는 탄성 변형을 담당한다. 논문은 이 분해를 2‑차원 표면 좌표계에 직접 적용함으로써, 성장에 따른 잔여 응력과 곡률 변화를 정확히 계산한다. 특히, 성장 텐서가 비등방성일 경우(예: 한 방향으로만 팽창) 쉘의 곡률이 급격히 변하면서 경계 불안정이 발생하는 메커니즘을 수식적으로 도출한다.

수치 구현 측면에서는, 서브디비전 기반 쉘 요소의 자유도는 정점당 3개의 변위와 3개의 회전(또는 곡률)으로 구성된다. 이때, 성장 텐서는 각 정점에 할당된 이중 텐서 필드로 저장되며, 시간 적분은 변형률-속도 연관식에 기반한 뉴턴‑라프슨 반복으로 수행한다. 요소 강성 행렬은 기존 Kirchhoff‑Love 쉘 행렬에 성장 텐서의 미분항을 추가함으로써 얻어지며, 이는 행렬의 대칭성을 유지한다.

검증 사례로는 (1) 원형 판의 균일 성장, (2) 구형 쉘의 축방향·방사방향 이방성 성장, (3) 경계에서 발생하는 주름 형성 등이 제시된다. 모든 사례에서 서브디비전 쉘 요소는 기존 4‑노드 쉘 요소 대비 2~3배 빠른 수렴 속도와 높은 정확도를 보였다. 특히, 구형 쉘 내부에 고정된 구멍에 성장하는 시트와, 반대로 구멍이 수축하면서 고정된 시트를 압축하는 과정을 역학적으로 동등함을 수치적으로 확인함으로써, 성장‑압축 문제 사이의 대칭성을 입증하였다.

이러한 결과는 생물학적 조직 성장, 부풀림형 소프트 로봇, 그리고 잔류 응력에 민감한 얇은 구조물 설계 등 다양한 분야에 바로 적용 가능함을 시사한다.