[논문 리뷰] Dual Quaternion Variational Integrator for Rigid Body Dynamic Simulation

Jiafeng Xu, Karl Henning Halse|arXiv (Cornell University)|2016. 11. 02.

Numerical methods for differential equations참고 문헌 19인용 수 9

한 줄 요약

이 논문은 이중 허수 벡터를 사용하여 모든 6자유도에서 강체 운동을 통합적으로 모의하는 심플렉틱 일차 리 군 적분기인 이중 허수 벡터 변분 적분기(DQVI)를 소개한다. 기하적 구조, 에너지, 운동량을 유지함으로써 DQVI는 기준점과 6×6 관성 행렬이 임의로 정의된 경우에도 큰 시간 간격으로 고정밀 실시간 시뮬레이션을 가능하게 한다.

ABSTRACT

We introduce a symplectic dual quaternion variational integrator(DQVI) for simulating single rigid body motion in all six degrees of freedom. Dual quaternion is used to represent rigid body kinematics and one-step Lie group variational integrator is used to conserve the geometric structure, energy and momentum of the system during the simulation. The combination of these two becomes the first Lie group variational integrator for rigid body simulation without decoupling translations and rotations. Newton-Raphson method is used to solve the recursive dynamic equation. This method is suitable for real-time rigid body simulations with high precision under large time step. DQVI respects the symplectic structure of the system with excellent long-term conservation of geometry structure, momentum and energy. It also allows the reference point and 6-by-6 inertia matrix to be arbitrarily defined, which is very convenient for a variety of engineering problems.

연구 동기 및 목표

• 모의 과정에서 강체 시스템의 심플렉틱 구조, 에너지, 운동량을 유지하는 기하학적 적분기를 개발하는 것.
• 특수 에울러리안 군 SE(3)를 사용하여 강체 역학에서 이동과 회전을 분리하지 않고 통합적으로 다루는 것.
• 높은 정밀도와 수치적 안정성을 유지하면서도 큰 시간 간격으로 실시간 시뮬레이션을 가능하게 하는 것.
• 더 넓은 공학적 적용 가능성을 위해 기준점과 6×6 관성 행렬을 임의로 정의할 수 있도록 하는 것.
• 다단계 적분 및 잠재력장과 비보존력이 있는 시스템으로의 확장을 위한 기반을 마련하는 것.

제안 방법

• SE(3)에서 강체 자세를 표현하기 위해 이중 허수 벡수를 사용하여 기하학적 특이점 없이도 회전과 이동을 단일하고 압축된 표현으로 통합한다.
• 재투영 또는 제약 조건 없이 리 군 위에서 시스템이 자연스럽게 진화하도록 보장하는 일차 리 군 변분 적분기를 적용한다.
• 리 대수 se(3)의 지수 사상 표현을 위해 6개의 독립 변수(Φ와 Ψ)를 사용하여 구성이 항상 SE(3) 위에 유지되도록 한다.
• 변분 원리를 통해 이산 동역학 방정식을 유도하여 연속 시스템의 기하학적 구조를 유지한다.
• 초기 추정치로 증분 자세 변화를 활용한 뉴턴-라프슨 방법을 사용하여 유도된 비선형 방정식을 해결한다.
• 잔차 방정식의 애널리틱한 야코비안 표현식을 사용하여 수렴 속도를 빠르고 안정적으로 보장한다.

실험 결과

연구 질문

• RQ1 이동과 회전을 분리하지 않고 전체 6자유도 강체 운동을 하나의 리 군으로 다룰 수 있는 변분 적분기를 구성할 수 있는가?
• RQ2 리 군 변분 적분기에서 이중 허수 벡수를 사용할 경우 기존 방법에 비해 에너지, 운동량, 기하학적 구조의 보존 성능이 향상되는가?
• RQ3 기준점과 6×6 관성 행렬을 임의로 정의한 경우에도 큰 시간 간격에서 높은 정밀도와 안정성을 유지할 수 있는가?
• RQ4 제안된 초기 추정 전략을 사용한 뉴턴-라프슨 솔버의 수렴 속도는 얼마나 빠른가?
• RQ5 미래 연구에서 비보존력과 잠재력장을 포함한 시스템으로의 확장을 이 방법이 어떻게 가능하게 하는가?

주요 결과

• DQVI 방법은 시스템의 심플렉틱 구조를 유지하여 장기적으로 에너지와 운동량을 매우 잘 보존한다.
• 기준점과 6×6 관성 행렬을 임의로 정의할 수 있어 해양 차량의 추가 질량과 저항성 등 공학 문제에 대한 적용 가능성이 높아진다.
• 적절한 초기화 조건에서 뉴턴-라프슨 솔버는 최소 3회 반복 내에 기계 정밀도에 도달하여 효율적인 실시간 시뮬레이션을 가능하게 한다.
• 증분 자세 변화를 기반으로 한 초기 추정은 매개변수화의 유효성을 보장한다(|Φ| < 1), 이는 시간 간격을 조절하여 180도 이내의 회전을 유지한다.
• 질량 중심에 기준점을 두고 순수 이동만 발생하는 경우, 뉴턴-라프슨 솔버는 1단계 전진 오일러 방법으로 축소되며, 이는 극한에서의 일致성을 확인한다.
• DQVI는 SE(3)를 전체로 간주하여 회전 및 이동 역학을 분리하지 않고 처리하는 최초의 리 군 변분 적분기이다.