I’d like to know specifically about an elegant way of deriving a second derivative of an orientation quaternion from a torque and a moment of inertia matrix, if one is available.
The straight forward, but somewhat round-about way that I can think of is to solve for the derivative of the orientation of a rigid body in terms of Euler angles via classical methods. Then I would convert those angular displacement vectors into unit quaternions and continually compose the unit quaternions.
I have been trying to find a formulation for mapping net torques as a function of time to orientation quaternions as a function of time, but have been unsuccessful.