2차 변형도 그라디언트 빔을 위한 새로운 약형 적분 요소

초록

본 논문은 8차 미분 방정식으로 기술되는 2차 변형도 그라디언트 Euler‑Bernoulli 빔 이론에 대해, 라그랑주와 헤르미트 보간을 이용한 두 종류의 약형 적분(Weak‑Form Quadrature) 요소를 제안한다. 내부에서는 변위만 자유도로 두고, 경계에서는 변위·기울기·곡률·삼차 미분을 추가 자유도로 배치한다. 가우스‑로바토‑레전드(GLL) 점을 요소 노드와 수치 적분에 동시에 활용하고, 가중계수 행렬을 수정하는 새로운 절차를 제시한다. 제안된 요소는 굽힘, 자유진동, 안정성 해석에서 높은 정확도와 효율성을 보인다.

상세 분석

이 연구는 기존 1차 변형도 그라디언트 빔 이론(6차 방정식)에서 한 단계 더 나아가, 2차 변형도 그라디언트 이론이 요구하는 8차 미분 방정식을 직접 다루는 최초의 약형 적분 요소를 개발했다는 점에서 학술적 의의가 크다. 핵심 아이디어는 ‘약형(weak‑form)’ 접근을 차용해 변위와 그 고차 도함수들을 변분 원리에서 자연스럽게 포함시키는 것이다. 이를 위해 저자들은 두 가지 보간 전략을 선택했다. 첫 번째는 전통적인 라그랑주 다항식으로 내부 노드에서 변위만을 근사하고, 경계에서는 추가 자유도를 별도로 연결한다. 두 번째는 헤르미트 보간을 이용해 변위·기울기·곡률·삼차 미분을 한 번에 표현함으로써 경계 자유도의 구현을 간소화한다.

가중계수 행렬의 수정 절차는 특히 주목할 만하다. 일반적인 차분법에서는 고차 도함수에 대한 정확한 가중치를 얻기 어려우나, 저자들은 기존 1차·2차·3차·4차 가중계수를 재정의하고, 경계 자유도에 대응하는 행과 열을 0 또는 기존 계수의 복사 형태로 채워 넣음으로써 8차 미분 연산을 정확히 구현한다. 이렇게 구성된 행렬을 GLL 점의 가중치와 결합해 요소 강성, 기하 강성, 질량 행렬을 구하고, 외력 벡터까지 일관되게 계산한다.

수치 실험에서는 균일 빔의 정적 굽힘, 자유진동 고유값, 압축 하중에 의한 좌굴 해석을 수행했으며, 전통적인 유한요소(FEM)와 기존 차분 기반 DQ 방법과 비교해 오차가 현저히 감소함을 확인했다. 특히 경계 조건이 복합적인 경우(예: 클램프‑프리, 클램프‑클램프)에도 정확한 결과를 얻어, 비고전적 자유도가 많은 고차 연속체 모델에 대한 적용 가능성을 입증했다.

한계점으로는 현재 1차원 빔 전용 구현에 머물러 있어, 2차원 판이나 쉘으로 확장하려면 보간 함수와 가중계수 재구성이 필요하다. 또한 GLL 점을 고정적으로 사용함에 따라 비균일 메쉬나 국부적 정밀도 향상이 제한될 수 있다. 향후 연구에서는 적응형 메쉬와 고차 보간을 결합하거나, 비선형 및 동적 하중 상황에 대한 확장도 기대된다.