총라그랑지 기반 대변형 다물체 동역학을 위한 유한요소 해석 프레임워크

초록

본 논문은 총라그랑지(Total Lagrangian) 관점에서 변형체의 대변형·대회전·대변위를 동시에 다룰 수 있는 유한요소(FEA) 프레임워크를 제안한다. 절대노드좌표형식(ANCF) 빔·쉘·테트라헤드론 요소들을 통합하고, St‑Venant‑Kirchhoff 및 Mooney‑Rivlin 하이퍼탄성 모델, Kelvin‑Voigt 점성 감쇠를 포함한다. 또한 다물체 시스템에 대한 위치·속도 제약을 증강라그랑지(augmented Lagrangian) 방식으로 처리한다.

상세 분석

이 연구는 총라그랑지(TL) 접근법을 기반으로, 변형도와 응력 모두를 기준구조(reference configuration)에서 정의한다는 점에서 기존의 업데이트 라그랑지(UL) 방식과 차별화된다. TL‑FEA는 변형도 그라디언트 F 를 핵심 변수로 삼아, 모든 요소(ANCF 빔, ANCF 쉘, 10‑노드 테트라헤드론)에서 동일한 수식적 틀을 적용한다. 특히, 전통적인 ANCF 해석에서 위치 r 를 S(X) e(t) 형태로 표현하던 것을, r = N(t) s(X) 형태로 재배열함으로써 F = N(t) H(X) ( H = ∂s/∂X )라는 명시적 구조를 얻는다. 이는 요소마다 다른 형태의 변형그라디언트를 동일한 행렬 연산으로 처리하게 해, 구현 복잡성을 크게 낮춘다.

요소별로는 2‑노드 3243 빔(12 DOF/노드), 4‑노드 3443 쉘(12 DOF/노드), 그리고 10‑노드 2차 테트라(T10) 고차원 요소를 제시한다. 각 요소는 전역 노드 행렬 N(t) 와 스칼라 형태함수 s(u) 를 통해 위치와 그라디언트를 동시에 보간한다. 특히, T10 요소는 베시안 좌표와 2차 라그랑지 형태함수를 사용해 C⁰ 연속성을 보장하고, 5‑점 가우시안 적분으로 정확한 질량·내부력 행렬을 얻는다.

재료 모델은 선형 St‑Venant‑Kirchhoff과 비선형 Mooney‑Rivlin을 포함한다. 두 모델 모두 내부 힘 f_int 과 일관된 탄성 탄성 행렬 K_tan 을 변형그라디언트 F 에 대한 미분으로 도출한다. 또한, Kelvin‑Voigt 점성 모델을 도입해 C_visc = η · \dot{F} 형식의 점성 행렬을 추가함으로써 고변형 동역학에서 수치적 안정성을 확보한다.

다물체 시스템의 제약은 위치 c(r,A,t)=0 형식으로 정의하고, 증강라그랑지(augmented Lagrangian) 방법을 적용한다. 라그랑지 승수 λ와 페널티 파라미터 ρ를 이용해 제약을 1차·2차 항으로 결합하고, 가상작업 원리를 통해 전체 방정식 M\ddot{q}+C\dot{q}+Kq+G=F_ext+J^Tλ 을 얻는다. 이때 J 은 제약의 야코비안이며, 반복적인 뉴턴‑라프슨 혹은 준-뉴턴 방식으로 λ와 q 를 동시에 업데이트한다.

전체 프레임워크는 대변형·대회전·다물체 상호작용을 하나의 통합된 수치 해법으로 구현함으로써, 기존에 각각 별도 구현되던 빔·쉘·솔리드 요소와 접촉·마찰 모델을 일관되게 결합한다는 점이 큰 장점이다. 또한, 변형그라디언트 중심의 설계는 새로운 요소 추가와 복합 재료 모델 확장을 용이하게 만든다.