라미네이트 복합판의 기하 비선형 등가기하학적 해석: 고차 전단 변형 이론 기반 IGA

초록

본 논문은 고차 전단 변형 이론(HSDT)을 적용한 라미네이트 복합판의 기하 비선형 거동을 등가기하학적 해석(IGA)으로 효율적으로 분석한다. B‑spline·NURBS 기반의 연속성 높은 형상함수를 이용해 전단 변형을 정확히 표현하고, 전체 라그랑주 접근법과 von‑Kármán 변형 가정을 결합해 큰 변형을 모델링한다. 다양한 재료·배향 조합에 대한 수치 검증을 통해 제안 방법의 정확성과 효율성을 입증한다.

상세 분석

이 연구는 복합재 라미네이트 판의 비선형 거동을 고차 전단 변형 이론(HSDT)과 등가기하학적 해석(IGA)을 결합함으로써 기존 유한요소법(FEM)의 한계를 극복한다는 점에서 의미가 크다. HSDT는 전통적인 1차 전단 변형 이론(CLT)에서 요구되는 전단 보정계수를 필요로 하지 않으며, 전단 변형 에너지를 정확히 계산한다. 이는 특히 두께‑길이비가 큰 복합판에서 전단 효과가 지배적일 때 중요한데, 전단 보정계수의 추정 오류가 전체 해석 정확도를 크게 저하시킬 수 있기 때문이다.

IGA는 B‑spline·NURBS 기반의 고차 연속 형상함수를 사용한다. HSDT는 변위장에 대해 2차 연속성(C¹)이 요구되는데, 전통 FEM에서는 이를 만족시키기 위해 추가 자유도나 특수 요소(예: Hermite 요소)를 도입해야 한다. 반면 IGA는 기본적으로 높은 차수와 연속성을 제공하므로, 별도의 보강 없이도 HSDT 모델을 직접 구현할 수 있다. 이는 구현 복잡성을 크게 낮추고, 메쉬 정밀도에 대한 의존성을 감소시킨다.

비선형성은 전체 라그랑주 접근법을 채택하고, 변형률 텐서는 von‑Kármán 가정을 적용해 2차 비선형 항을 포함한다. 이는 중·대변형 영역에서의 응력‑변형 관계를 정확히 포착한다. 논문은 뉴턴‑라프슨 반복법을 이용해 비선형 방정식을 해결하고, 수렴성을 확보하기 위해 적절한 초기값과 스텝 크기 제어 전략을 제시한다.

수치 검증에서는 등방성 금속판, 직교 이방성(orthotropic) 판, 교차 적층(cross‑ply) 및 각도 적층(angle‑ply) 복합판을 대상으로 선형·비선형 정적 굽힘, 자유진동, 그리고 압축‑굽힘(플라터) 문제를 풀었다. 결과는 전통 FEM(예: 8‑node Lagrange 요소)과 비교했을 때, 동일 메쉬에서 더 높은 정확도를 보이며, 메쉬를 절반으로 줄여도 오차가 크게 증가하지 않는다. 특히 전단 변형이 큰 두께‑비율 높은 판에서 HSDT‑IGA 조합이 전단 보정계수를 사용한 CLT‑FEM보다 월등히 우수함을 확인했다.

또한, 층간 전이 효과와 각도 적층에 따른 비선형 경계조건 변화도 자연스럽게 반영된다. 이는 설계 단계에서 복합재 구조물의 안전성 평가와 최적화에 직접 활용될 수 있다. 마지막으로 계산 비용 측면에서, 고차 연속성 덕분에 자유도는 다소 증가하지만, 수렴 속도가 빨라 전체 시뮬레이션 시간이 기존 FEM 대비 30~40% 단축되는 것으로 보고된다.

이와 같이, 논문은 HSDT와 IGA의 시너지 효과를 통해 복합재 라미네이트 판의 기하 비선형 해석을 보다 정확하고 효율적으로 수행할 수 있음을 입증한다. 향후 동적·충격 해석, 비선형 재료 모델(예: 플라스틱·손상)과의 결합, 그리고 3‑D 복합 구조에 대한 확장 가능성도 제시한다.