Abstract

Unit quaternion representation is widely used in flight simulation to overcome the limitations of the standard numerical ordinary-differential-equations (ODEs) based on three-parameters rotation variables (such as Euler angels), as they may impose kinematic singularities during aircraft's attitude reconstruction. However, these benefits do not come without a price, since the classical way of integrating rotational quaternions includes solving of differential-algebraic equations (DAEs) that requires post-integration numerical stabilization of the additional algebraic constraint enforcing the quaternion unitary norm. This can pose a problem in the case of longer flight simulations since improper numerical treatment of the quaternion-normalization constraint may induce numerical drift into the simulation results. As a remedy, the proposed novel algorithm circumvents DAE problem of quaternion integration by shifting update-integration-process from configuration manifold to the local tangential level of the incremental rotations (reducing thus integration to standard three ODEs problem at tangential Lie algebra level). This can be done due to the isomorphism of the Lie algebras of the rotational SO(3) group and the configuration manifold unit quaternion Sp(1) group. Besides avoiding DAE formulation by reducing integration process to standard three ODEs problem, the proposed algorithm also exhibits numerical advantages as it is discussed in the presented example.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.