3D Rigid body dynamics - Rod over rotating disk

[Quesadilla]

Homework Statement

A thin uniform rod is attached to an axis through its midpoint. The axis is standing on a disk rotating with constant angular speed \Omega about its symmetry axis. The rod's midpoint is located directly above the rotational axis of the disk. Let \theta denote the rod's angle with the horizontal plane. Determine how \theta varies with time.

Pictures are attached below to make the above description clearer.

Homework Equations

Tensor of inertia, I_{\alpha \beta} = \int (r^2 \delta_{\alpha \beta} - r_{\alpha} r_{\beta}) dm

Euler's equations for principal systems:

M_1 = I_{11} \dot{\omega}_1 - \omega_2 \omega_3 (I_{22} - I_{33})
M_2 = I_{22} \dot{\omega}_2 - \omega_3 \omega_1 (I_{33} - I_{11})
M_3 = I_{33} \dot{\omega}_3 - \omega_1 \omega_2 (I_{11} - I_{22})

where M_i denotes the components of the net applied moment of force and \mathbf{\omega} = \omega_1 \mathbf{e}_1 + \omega_2 \mathbf{e}_2 + \omega_3 \mathbf{e}_3 is the angular velocity vector, with \mathbf{e}_i being the principal basis vectors.

The Attempt at a Solution

I have tried several things, none of which have yielded satisfactory results. This is amongst the first 3D rigid body dynamics problems I have tried to solve.

Assume an inertial Cartesian coordinate system with the xy-plane coinciding with the horizontal plane of the disk and the z-axis pointing upward perpendicular to this plane. Now introduce an accelerated reference frame constant in the rod. Let the x'-axis be perpendicular to the rod, such that it at a time t=0 is positioned directly above the x-axis with the angle between them being \theta. In the same manner put the y'- and z'-axes, such that they make angles \theta with the y- and z-axis respectively.

In this frame we can write \Omega \mathbf{e}_z = \Omega \sin \theta \mathbf{e}_{x'} + \Omega \cos \theta \mathbf{e}_{z'}.

I am not sure how the total angular velocity vector \omega should look like now. I have tried both using the above and adding a term \dot{\theta} \mathbf{e}_{x'} to it. Neither seems to give me a result I want.

Now, I could be wrong but I think this is a principal system, so the tensor of inertia is diagonal and its diagonal elements are given by

I_{x'x'} = \frac{1}{12} ml^2
I_{y'y'} = 0 (I think, but the value shouldn't matter?)
I_{z'z'} = \frac{1}{12} ml^2

Euler's equations then gives

\frac{1}{12} ml^2 (\Omega \cos \theta \dot{\theta}+ \ddot{\theta}) = M_1 = 0.

0=0

\frac{1}{12} ml^2 (-\Omega \sin \theta \dot{\theta}) = M_3 = 0.

Here I am not happy with what I have got, so I must have gone wrong somewhere in the above. Any hints to where I am messing up are much appreciated. Thank you.

 

[Quesadilla]

I now realize one mistake I made in my original post. The paragraph about the accelerated reference frame introduced makes no sense. The origin 0' is put at the midpoint of the rod, such that the origin remains constant throughout the motion. The x- and x'-axes should coincide, and be perpendicular to rod and parallel to the xy-plane (plane of the disk). The y'- and z'- axes are then tilted an angle \theta, such that O'x'y'z' becomes a principal coordinate system for the rod.

Then I think the angular velocity vector \mathbf{\omega} should read something like
\mathbf{\omega} = \dot{\theta} \mathbf{e}_{x'} + \Omega \sin \theta \mathbf{e}_{y'} + \Omega \cos \theta \mathbf{e}_{z'}

Euler's equations then gives

\frac{1}{12}ml^2 \ddot{\theta} = 0

0= 0

-\frac{1}{12}ml^2 \Omega \dot{\theta} \sin \theta - \frac{1}{12}ml^2 \Omega \dot{\theta} \sin \theta = 0.

Maybe there should be an extra minus sign in the last equation, so that the terms will cancel out, but I can't figure out where the mistake is if there is one. Then I get \ddot{\theta} = 0 and I guess \theta(t) = \dot{\theta}_0 t + \theta_0.