기하학적으로 정확한 원형 단면 로드의 교차면 변형을 포함한 리 그룹 변분 적분기

초록

본 논문은 원형 단면을 가진 코스테라트 로드에 평면 교차면 변형을 도입하고, 이를 연속적인 라그랑지안 형태로 유도한다. 이후 리 그룹 변분 적분기(LGVI) 기법을 적용해 회전 및 단면 팽창·수축을 동시에 보존하는 이산 모델을 구축한다. 제안된 모델은 부피 보존, 회전 구성 보존, 에너지 오차의 유한한 경계 유지 등 기하학적 구조를 유지하면서도 높은 정확도를 보인다. 다양한 초기 조건에 대한 수치 실험을 통해 모델의 물리 재현성과 장기 안정성을 검증하였다.

상세 분석

이 연구는 기존 코스테라트 로드 모델이 단면을 강체로 가정하는 한계를 극복하고, 원형 단면이 변형 후에도 원형을 유지한다는 가정 하에 평면 변형을 허용한다. 핵심은 국부적인 팽창 계수 ε(S,t) 를 도입해 현재 아크 길이 s 와 초기 아크 길이 S 사이의 비율을 정의하고, 질량 보존을 이용해 단면 면적 A(S) = Ā(S) ε(S,t) 으로 표현함으로써 부피 보존을 자연스럽게 구현한다. 이 팽창 계수는 연속 방정식에 새로운 변형 자유도를 제공하면서도, 전단·굽힘·비틀림과 같은 기존 변형 모드와 결합되어 전체 변형 텐서를 재정의한다.

라그랑지안은 전역적인 운동 에너지와 변형 에너지로 구성된다. 전이 운동 에너지는 질량이 ε 에 의존하지 않으므로 기존 식과 동일하게 유지된다. 반면 회전 운동 에너지는 관성 모멘트 J(S,t) 가 단면 면적의 제곱에 비례함에 따라 ε 의 제곱에 비례하게 변한다. 이는 원형 단면의 두 번째 모멘트 I₁=I₂= A² 4π, I₃= A² 2π 에 ε 을 삽입함으로써 얻어진다. 따라서 회전 관성은 축방향 연신에 따라 동적으로 변하며, 이는 변분 적분 단계에서 정확히 반영된다.

이산화 과정에서는 구성 변수 (xₖ,Rₖ) 와 속도 변수 (γₖ,ωₖ) 를 사용해 시간 간격 h 에 대한 디스크리트 라그랑지안을 구성한다. 리 그룹 구조 SO(3) 위에서의 변분 원리를 적용하면, Cayley 변환을 이용한 Rₖ₊₁ = Rₖ cay(h Ωₖ) 형식의 업데이트 식이 도출된다. 이 식은 회전 행렬이 수치적으로 SO(3) 에 남아 있음을 보장한다. 또한, 변분 원리 자체가 심플렉틱성을 내포하므로, 전체 시스템은 에너지와 선형·각운동량을 장기적으로 보존한다.

수치 실험에서는 (1) 자유 진동, (2) 외부 하중에 의한 비틀림·연신, (3) 복합적인 초기 변형을 가진 경우를 다루었다. 모든 실험에서 부피 보존( A·ds = Ā·dS )이 수치적으로 10⁻⁸ 이하의 오차로 유지되었으며, 에너지 진동은 시간에 따라 유한한 경계 내에서 안정적으로 유지되었다. 특히, 팽창 계수를 포함하지 않은 기존 LGVI와 비교했을 때, 단면 변형이 큰 경우(연신률 ≈ 30 %)에 현저히 더 정확한 변형 및 응력 분포를 제공한다는 점이 강조된다.

한계점으로는 단면이 원형을 유지한다는 가정이 실제 비원형 혹은 비대칭 변형을 다루기 어렵게 만든다. 또한, 선형 응력‑변형 관계와 포아송 비를 명시적으로 고려하지 않아, 큰 압축·인장 상황에서 정확도가 떨어질 가능성이 있다. 향후 연구에서는 비원형 단면, 비선형 재료 모델, 그리고 3차원 변형을 포함한 일반화된 팽창 텐서를 도입함으로써 모델의 적용 범위를 확대할 수 있을 것이다.