1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

A Singularity also with Euler-Rodrigues parametrisation?

  1. Mar 2, 2017 #1
    Trying to develop a multibody dynamics software of my own (just to understand the nitty-gritty details of such stuff), I chose the Lagrangian equations approach, with the Euler-Rodrigues parametrisation (quaternion) for 3D rotation as it is supposed to remove the gimbal locking singularity of conventional Euler angles. When angular position is to be part of the generalized coordinates, I take b, c and d and leave a as a dependant variable computed from the normalisation condition: ##a^2= b^2+c^2+d^2 ## . I take the positive value of a, by convention. Sounds good but... The 3x3 transformation from angular velocity vector ## \begin{pmatrix} u \\ v \\ w \end{pmatrix} ## to vector ## \begin{pmatrix} \dot b \\ \dot c \\ \dot d \end{pmatrix} ## is singular for any rotation such that a=0 (determinant is a/8), which a π rd rotation about any axis. This singularity is a potential problem for the calculation of the inverse transform ## \frac {\partial \Omega} {\partial \dot q}## and also ## \frac {\partial^2 \Omega} {\partial \dot q \partial q}## that are necessary for computing the Lagrange equations.
    Looks like another sort of gimbal-locking case! Did I miss something :confused: ? Or how to circumvent this problem?
    Thanks for any clarification!
  2. jcsd
  3. Mar 3, 2017 #2
    Evidently you're normalization introduces singularities in the Euler-Rodriguez paramemterization for rotations greater than##\pi## .I suggest you Google: O. Buchan, L. Trained "The Vectorial Paramemterization of Rotation". Try implementing the scheme found in Appendix A.
  4. Mar 4, 2017 #3
    Excellent! I should have spotted this reference before posting since I know some other publications of Olivier Bauchau (relevant to non linear elasticity modelling of slender beams).
    Thanks a lot for support!
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted