2
$\begingroup$

An object with known initial orientation and angular velocity $\omega_0$ is subjected to constant angular acceleration $\alpha$. Its final orientation at time $t$ is some rotation of its initial orientation. How can I find this rotation?

Angular velocity and acceleration are 3D vectors. Note, $\alpha$ is not necessarily parallel to $\omega_0$, so you can't just use the rotation $\omega_0 t+\alpha t^2$ about a common axis.

EDIT: That should be $\omega_0 t+(\alpha /2) t^2$

Thanks.

  • 0
    Thanks. You're right. But only if the axes are parallel.2011-01-16
  • 0
    $\phi = \phi_0+ \omega_{0}t + \frac{1}{2} \alpha t^2$.2011-01-16
  • 1
    Thanks Trevor, but about what axis?2011-01-16
  • 0
    Is the acceleration axis fixed to the body or fixed to the global frame? That makes a difference.2011-01-16
  • 0
    The acceleration vector is fixed in the global frame.2011-01-16
  • 0
    Hmm, I vaguely remember seeing this on Physics SE. Did this get moved perchance?2011-01-16

2 Answers 2

1

I assume you mean $\vec{\omega}_0$ and $\vec{\alpha}$ as vector quantities. Unlike rotations, angular speeds and accelerations are vectors that can be added like

$$ \mathrm{d}\vec{\omega} = \vec{\alpha}\,\mathrm{d}t $$

I suppose the end result in speed is

$$ \vec{\omega} = \vec{\omega}_0 + \vec{\alpha}\,t $$

Now for the orientations things are more complicated. They do not add up as vectors. What you have to do is decompose the orientation to three euler angles $\varphi$, $\psi$ and $\theta$ and integrate their derivatives to find the final orientation. If the 3x3 rotation matrix is $R=R_y(\varphi)\,R_x(\psi)\,R_z(\theta)$ then

$$ \vec{\omega} = \hat{j}\dot{\varphi}+R_y(\varphi)\left(\hat{i}\dot{\psi}+R_x(\psi)\,\hat{k}\dot{\theta}\right)$$

which is a result of the sequencial rotations about $y$, $x$ and $z$.

The above can be expressed as $\vec{\omega}=\mathrm{J}\,\pmatrix{\varphi,\psi,\theta}\cdot\begin{bmatrix}\dot{\varphi} \\ \dot{\psi} \\ \dot{\theta}\end{bmatrix}$ with $\mathrm{J}$ the jacobian matrix.

Next the acceleration needs to be decomposed by

$$\vec{\alpha}=\mathrm{J}\,\pmatrix{\varphi,\psi,\theta}\cdot\begin{bmatrix}\ddot{\varphi} \\ \ddot{\psi} \\ \ddot{\theta}\end{bmatrix}+\mathrm{\dot{J}}\,\pmatrix{\varphi,\psi,\theta}\cdot\begin{bmatrix}\dot{\varphi} \\ \dot{\psi} \\ \dot{\theta}\end{bmatrix}$$

and solved as

$$ \begin{bmatrix}\ddot{\varphi} \\ \ddot{\psi} \\ \ddot{\theta}\end{bmatrix} = \vec{F}\pmatrix{\varphi,\psi,\theta,\dot{\varphi},\dot{\psi},\dot{\theta}} $$

and integrated (analytically not likely, but numerically).

Good luck!

  • 0
    I think I understand, but I'm not sure about the last step. The Euler angles and their rates in F are unknown functions of t, so how do I go about integrating them? Do you mean to consider it as a system of differential equations to be solved?2011-01-16
  • 0
    @SteveB - yes the angles and rates are initial conditions for each step.2011-01-16
  • 0
    Caution when $\psi=\pi/2$ because $\vec{F}$ becomes singular.2011-01-16
0

If being considered numerically, it's going to be easier I think to just find angular velocity vectors $\omega_t=(\omega_0+\alpha t_i)$ for each timestep, construct a matrix to rotate about each vector by the appropriate amount, and multiply them together. But I was hoping for something a bit faster and more elegant. Is an algebraic answer impossible, or just difficult?