하이브리드 NURBS 강화 유한요소법 육면체 메쉬

초록

본 논문은 NURBS 기반 경계 표현을 3차원 육면체 메쉬에 결합한 하이브리드 유한요소법(NEFEM‑HEX)을 제안한다. 경계층에는 NURBS‑강화 요소를, 내부 영역에는 1차 라그랑주 요소를 사용하고, 특수 적분 규칙과 안정적인 보간 연산자를 도입한다. h‑refinement와 knot insertion을 동시에 적용할 수 있으며, Poisson 문제에 대한 a priori 오류 추정과 수치 실험을 통해 이론적 수렴성을 확인한다.

상세 분석

이 연구는 기존의 등변형 유한요소법(isoparametric FEM)과 이소지오메트릭 분석(IGA)의 장점을 융합한 새로운 하이브리드 프레임워크를 제시한다. 먼저, 열린 볼록 영역 Ω의 경계 ∂Ω를 단일 NURBS 패치 P로 모델링하고, 이를 기반으로 경계층 Ω_B와 내부 영역 Ω_int를 구분한다. 경계층 요소는 NURBS 기반의 유리형 변환 ˜F_Q를 통해 정확한 곡면을 재현하고, 내부 요소는 전통적인 선형 등변형 변환을 사용한다. 이때, 경계층 요소는 최대 다섯 개의 곡면을 가질 수 있으며, 하나는 NURBS 면, 나머지는 내부 곡면으로 구성된다.

핵심 기여는 두 가지이다. 첫째, NURBS‑강화 요소에 대한 특수 적분 규칙을 설계하여 곡면상의 정확한 수치 적분을 가능하게 한다. 기존 NEFEM에서는 삼각형·테트라hedral 메쉬에 적용했으나, 본 논문은 육면체 메쉬에 맞춰 3차원 부피 적분을 효율적으로 수행한다. 둘째, 보간 연산자를 정의하고 그 안정성 및 근사성을 정리함으로써, h‑refinement(전통 FEM의 격자 세분화)와 knot insertion(IGA의 매듭 삽입)을 동시에 적용할 수 있는 이론적 기반을 제공한다. 이는 경계층에서 NURBS 매개변수와 내부 격자 간의 일관성을 유지하면서도, 전체 해석의 자유도를 효율적으로 제어한다는 점에서 의미가 크다.

오차 분석에서는 H^1 노름에 대한 a priori 추정식을 도출하였다. 경계층에서 NURBS 기반 보간이 제공하는 고차 정확도와 내부 선형 요소의 전형적인 1차 수렴률이 결합되어, 전체 해석이 기대되는 최적 수렴률을 달성함을 증명한다. 특히, 경계가 정확히 표현되므로 기하학적 오차가 해석 오차에 미치는 영향을 제거할 수 있다.

수치 실험으로는 3차원 Poisson 문제를 선택하고, 다양한 격자 정밀도에서 오류 수렴 곡선을 제시하였다. 결과는 이론적 추정과 일치하며, 특히 경계 근처에서 전통 등변형 FEM보다 현저히 낮은 오류를 보였다. 또한, knot insertion을 통한 국소 정밀도 향상이 h‑refinement과 동일한 수렴 효과를 제공함을 확인했다.

이와 같이 NEFEM‑HEX는 CAD와 CAE 사이의 인터페이스 문제를 완화하면서, 복잡한 곡면을 정확히 다루는 동시에 전통 FEM의 구현 용이성을 유지한다는 장점을 가진다. 다만, 텐서곱 형태의 NURBS 기반 함수가 로컬 메쉬 정밀도에 제약을 주는 점, 그리고 고차 매시와 다중 매듭 관리가 추가적인 구현 복잡성을 야기한다는 한계도 논의된다. 향후 연구에서는 계층적 B‑스플라인이나 T‑스플라인을 도입해 로컬 정밀도 제어를 개선하고, 경계층과 내부 영역 사이의 비선형 결합을 확장하는 방향이 제시된다.