[논문 리뷰] Lie group variational integrators for rigid body problems using quaternions

Xuefeng Shen, Melvin Leok|arXiv (Cornell University)|2017. 05. 11.

Numerical methods for differential equations참고 문헌 13인용 수 3

한 줄 요약

이 논문은 S³의 내재된 리 군 구조를 활용하여 제약 조건 없이 단위 허수를 사용하는 강체 역학을 위한 리 군 변분 적분법을 개발한다. 이 방법은 장기간에 걸쳐 심플렉틱성, 에너지, 단위 허수 노름, 각운동량을 유지하며 항공우주, 로봇공학, 컴퓨터 그래픽스 분야에서 계산 효율성과 기하학적 정확성을 제공한다.

ABSTRACT

Rigid body dynamics on the rotation group have typically been represented in terms of rotation matrices, unit quaternions, or local coordinates, such as Euler angles. Due to the coordinate singularities associated with local coordinate charts, it is common in engineering applications to adopt the unit quaternion representation, and the numerical simulations typically impose the unit length condition using constraints or by normalization at each step. From the perspective of geometric structure-preserving, such approaches are undesirable as they are either computationally less efficient, or interfere with the preservation of other geometric properties of the dynamics. In this paper, we adopt the approach used in constructing Lie group variational integrators for rigid body dynamics on the rotation group to the representation in terms of unit quaternions. In particular, the rigid body dynamics is lifted to unit quaternions, and the Lie group structure of unit quaternions is used to represent tangent vector intrinsically, thereby avoiding the use of a Lagrange multiplier. A Lie group variational integrator in the unit quaternion representation is derived, and numerical results are presented.

연구 동기 및 목표

• 제약 조건이나 정규화에 의존하는 전통적인 단위 허수 방법의 한계를 해결하여 기하학적 구조 유지 성능을 향상시키기 위해.
• 강체 역학을 위한 전역적으로 유효하고 구조를 유지하는 수치적 적분법을 단위 허수 표현을 사용하여 개발하기 위해.
• 국소 좌표계나 라그랑주 승수를 사용하지 않고도 심플렉틱성, 에너지, 단위 노름, 각운동량과 같은 핵심 기하학적 불변량을 유지하기 위해.
• 항공우주 및 컴퓨터 그래픽스와 같이 허수를 표준으로 사용하는 분야에서 강체 운동의 효율적이고 안정적인 장기 시뮬레이션을 가능하게 하기 위해.
• 기하학적 적분 이론과 공학적 실무 간 격차를 해소하기 위해 널리 사용되는 단위 허수 형식에 리 군 변분 적분법을 적응시키기 위해.

제안 방법

• 강체 역학을 SO(3)에서 단위 허수 군 S³로 상향 변환하여, 허수 곱셈에 의해 정의되는 리 군으로 간주한다.
• 접선 벡터를 S³의 리 대수(실수 공간 R³로 식별됨)를 통해 내재적으로 표현함으로써, 명시적 제약 조건이나 정규화가 필요 없도록 한다.
• 지수 사상과 군 구조를 사용하여 Q × Q 위에 이산 라그랑지안을 구성함으로써, 이산 방정식이 연속 시스템의 기하학적 성질을 그대로 이어받도록 보장한다.
• 이산 오일러-라그랑ż 방정식과 이에 상응하는 심플렉틱 사상은 이산 레지오트르 변환을 통해 유도되며, 운동량과 심플렉틱 구조를 유지한다.
• 각 시간 단계에서 비선형 방정식을 해결하기 위해 브로이든 유형의 방법을 구현함으로써 정확한 야코비안 계산을 피한다.
• 제약 조건이나 투영을 통해가 아니라 군 구조 자체를 통해 단위 허수 조건이 내재적으로 유지되도록 보장한다.

실험 결과

연구 질문

• RQ1 제약 조건이나 정규화에 의존하지 않고 단위 허수를 사용하여 리 군 변분 적분법을 강체 역학에 적용할 수 있는가?
• RQ2 이러한 방법이 장기간의 통합 시간 동안 심플렉틱성, 에너지, 허수의 단위 노름을 유지하는가?
• RQ3 에너지 및 노름 유지 성능 측면에서 기존의 허수 기반 방법과 비교해 볼 때 이 적분법의 성능은 어떠한가?
• RQ4 군 작용에 대한 불변성 덕분에 각운동량 및 기타 보존량을 유지할 수 있는가?
• RQ5 기하학적 정확성을 유지하면서도 매트릭스 기반 리 군 적분법보다 계산 비용이 낮은가?

주요 결과

• 제안된 적분법은 1000 단위 시간 동안 오차 수준 10⁻¹⁴ 이내로 단위 허수 노름을 유지하여 내재된 제약 조건 만족을 입증한다.
• 에너지 오차는 1000 단위 시간 동안 10⁻⁷ 수준에서 유한하고 안정적으로 유지되어 심플렉틱성에 기인한 장기적인 에너지 안정성을 확인한다.
• 총 각운동량 오차는 약 10⁻¹¹ 수준에서 유한하게 유지되어 군 불변성에 기인한 이산 노이터 정리에 의해 보존됨을 확인한다.
• 계산 효율성이 확보된다: 허수는 4개의 스칼라로 표현되며(행렬의 9개보다 적음), 허수 곱셈은 16회의 곱셈으로 수행되며(행렬 곱셈의 27회보다 적음).
• 라그랑주 승수를 사용하지 않으며 정규화 단계도 필요로 하지 않아 계산 오버헤드가 감소하고 기하학적 구조가 유지된다.
• 수치 결과는 중심 중력장 하에서 강체를 시뮬레이션할 때 에너지나 노름의 드리프트 없이 뛰어난 장기 안정성과 정확성을 보여준다.