2
$\begingroup$

Well I stumbled upon the next differential equation system: $$ \begin{cases} \dot{\omega}_1 + a \omega_2 \omega_3 = t_1 \\ \dot{\omega}_2 + b \omega_1 \omega_3 = t_2 \\ \dot{\omega}_3 + c \omega_1 \omega_2 = t_3\end{cases}$$ with $a,b,c \in \mathbb{R}$. If necessary you can assume $t_2 = t_3 = 0$. I don't know if I ever solved this or not, but I do not remember how to do it.

  • 0
    These are the euler equations for a rotating top under the action of a torque. Instead of $\omega$'s as variables it's used the Euler angles and usually for the restricted case of a symmetric top under constant force. The lagrangian approax is too preferred (as it offers almost from the start a first integral for the system, the energy of the top). I recognized vaguely them when answer the question believing they could be integrated with not much difficulty. It's not so:2017-06-24
  • 0
    [Here](http://hepweb.ucsd.edu/ph110b/110b_notes/node36.html) you have the relation between the Euler angles and the $\omega$'s and [here](http://www.srl.caltech.edu/phys106/p106a00/topic8.pdf) (page 7) the Euler equations.2017-06-24

1 Answers 1

1

It's a no linear system, needing its own trick to solve, if any. I'll try first the case with $t_2=t_3=0$. It's possible to separate variables.

$$\begin{cases} \dot{\omega}_1 + a \omega_2 \omega_3 = t_1 \\ \dot{\omega}_2 + b \omega_1 \omega_3 =0 \\ \dot{\omega}_3 + c \omega_1 \omega_2 =0 \end{cases}$$

$$\frac{\mathrm d\omega_3/\mathrm dt}{\mathrm d\omega_2/\mathrm dt}= \frac{c\omega_2}{b\omega_3}\implies\frac{\mathrm d\omega_3}{\mathrm d\omega_2}=\frac{c\omega_2}{b\omega_3}\implies\omega_3\mathrm d\omega_3=(c/b)\omega_2\mathrm d\omega_2$$

$$\omega_3^2=(c/b)\omega_2^2+c_1$$

$$\frac{\mathrm d\omega_2/\mathrm dt}{\mathrm d\omega_1/\mathrm dt}= \frac{-b\omega_1\omega_3}{t_1-a\omega_2\omega_3}\implies\frac{\mathrm d\omega_2}{\mathrm d\omega_1}=\frac{-b\omega_1\sqrt{(c/b)\omega_2^2+c_1}}{t_1-a\omega_2\sqrt{(c/b)\omega_2^2+c_1}}\implies$$

$$\implies\frac{t_1-a\omega_2\sqrt{(c/b)\omega_2^2+c_1}}{\sqrt{(c/b)\omega_2^2+c_1}}\mathrm d\omega_2=-b\omega_1\mathrm d\omega_1$$

$$(a/b)\omega_2^2-\frac{2t_1\ln\left|\sqrt{\frac{c}{bc_1}}\omega_2+\sqrt{1+\frac{c}{bc_1}\omega_2^2}\right|}{\sqrt{c/b}}+c_2=\omega_1^2$$

And now, things get increasingly "messy" when you try to recover the t dependence:

$$\dot\omega_2=-b\left((a/b)\omega_2^2-\frac{2t_1\ln\left|\sqrt{\frac{c}{bc_1}}\omega_2+\sqrt{1+\frac{c}{bc_1}\omega_2^2}\right|}{\sqrt{c/b}}+c_2\right)^{1/2}\left((c/b)\omega_2^2+c_1\right)^{1/2}$$

At least, it's separable too. I am not able to discuss the assimptotic behavior neither the orbits, if any, the critical points...

Anyway, I hope this helps.

  • 0
    thank you very much for your attention2017-02-24