Taming the T-handle continued

Posted on 2017/02/13 by arkajad

This is a simple continuation from the last post “Taming the T-handle“. We ended up with the equation

(1)   \begin{equation*}\frac{d\psi(t)}{dt}=c_1+\frac{c_2}{1+c_3\,\sn^2(Bt,m)},\end{equation*}

where

(2)   \begin{equation*}c_1=\frac{L}{I_3},\quad c_2=L\left(\frac{1}{I_1}-\frac{1}{I_3}\right),\quad c_3=\frac{I_3(I_2-I_1)}{I_1(I_3-I_2)}.\end{equation*}

To find \psi(t), let us assume that \psi(0)=0, we need to integrate

    \[\psi(t)=c_1 t+c_2\int_0^t \frac{1}{1+c_3\,\sn^2(Bs,m)}\,ds.\]

Setting Bt=u we transform the above into

    \[\psi(t)=c_1 t+\frac{c_2}{B}\int_0^{Bt} \frac{1}{1+c_3\,\sn^2(u,m)}\,du.\]

Searching the net we can find the integral of this type, for instance on the Wolfram’s page about elliptic integrals of the third order we can find:

    \[\Pi(n;z|m)=\int_0^{F(z|m)}\,\frac{1}{1-n\,\sn(u|m)^2}\, du.\]

Then we set F(z|m)=Bt, therefore z= \am(Bt,m) – see Jacobi amplitude- realism or cubism,
so that

    \[\psi(t)=c_1 t+\frac{c_2}{B}\,\Pi(-c_3;\am(Bt,m)|m).\]

In principle we could now substitute the values of c_1,c_2,c_3 and consider our task essentially done. But that would be not very prudent! The point is that we have the parameter m and, depending on the case, this parameter can have value: <1, 1, and >1. The case m=1 is very special and should be treated separately. No elliptic functions are needed in this case. Moreover not every software used for calculations and simulations will know what to do with m>1, and even if it pretends to know, it may happen that in this domain the software is not sufficiently tested and debugged. It is for this reason that we now consider the three cases separately.

The case of d<1/I_2, i.e. m<1.

This is the most straightforward case. From Taming the T-handle, Eq. (10), we have

(3)   \begin{equation*}B=L\sqrt{\frac{(1-dI_1)(I_3-I_2)}{I_1I_2I_3}},\end{equation*}

(4)   \begin{equation*}m=\frac{(dI_3-1)(I_2-I_1)}{(1-dI_1)(I_3-I_2)}.\end{equation*}

Thus

    \[\frac{c_2}{B}=(I_3-I_1)\sqrt{\frac{I_2}{I_1I_3(I_3-I_2)(1-dI_1)}},\]

and we are done.

The case of d=1/I_2, i.e. m=1.

This is a tricky case. For m=1 we have (see Jacobi Elliptic Functions, Eqs. (33)-(35), and this comment)

(5)   \begin{eqnarray*} \cn(u,1)&=&\frac{1}{\cosh u},\\ \sn(u,1)&=&\tanh u,\\ \dn(u,1)&=&\frac{1}{\cosh u}. \end{eqnarray*}

To find \psi(t) we need to integrate

    \[\psi(t)=c_1 t+c_2\int_0^t \frac{1}{1+c_3\,\tanh^2(Bs,m)}\,ds.\]

We can go to Wolfram Alpha, and ask it to integrate it. It will, but we will not like the answer. It looks very-very ugly!

In principle we could try to simplify it, but it is better to make one step back in order to make two steps forward! Let us go back to Eq. (1) and write it as just one fraction:

(6)   \begin{equation*}\frac{d\psi(t)}{dt}=c_1+\frac{c_2}{1+c_3\,\tanh^2 Bt} =\frac{c_1+c_2+c_1c_3\tanh^2 Bt}{1+c_3\tanh^2 Bt}.\end{equation*}

Now, we write \tanh^2 u=\sinh^2 u /\cosh^2 u=(\cosh^2u -1)/\cosh^2 u, and get

(7)   \begin{equation*}\frac{d\psi(t)}{dt}=\frac{-c_1c_3+(c_1+c_2+c_1c_3)\cosh^2 Bt}{-c_3+(1+c_3)\cosh^2 Bt}.\end{equation*}

So, we need to integrate the function

(8)   \begin{equation*}\frac{d\psi}{dt}=\frac{\alpha+\beta \cosh^2Bt}{\gamma+\delta\cosh^2 Bt}.\end{equation*}

We can go again to Wolfram Alpha and what we get this time is so much nicer!

All we have to substitute the constants. It is a mechanical task, so we can use REDUCE, so we use REDUCE or any other software capable of symbolic computation. I used Mathematica. Here is the result, followed by a complete summary of the case:

(9)   \begin{equation*}\psi(t)=\frac{Lt}{I_2}+L\arctan\left( \sqrt{\frac{(I_2-I_1)I_3}{I_1(I_3-I_2)}}\tanh\sqrt{\frac{(I_2-I_1)(I_3-I_2)}{I_1I_2^2I_3}}t \right)\end{equation*}

(10)   \begin{eqnarray*}A_1&=&L\sqrt{\frac{I_1(I_3-I_2)}{I_2(I_3-I_1)}},\\ A_2&=&L,\\ A_3&=&L\sqrt{\frac{I_3(I_2-I_1)}{I_2(I_3-I_1)}},\\ B&=&\frac{L}{I_2}\sqrt{\frac{(I_2-I_1)(I_3-I_2)}{I_1 I_3}},\end{eqnarray*}

(11)   \begin{eqnarray*} L_1(t)&=&\frac{A_1}{\cosh Bt},\\\ L_2(t)&=&A_2 \tanh Bt,\\\ L_3(t)&=&\frac{A_3}{\cosh Bt},\end{eqnarray*}

(12)   \begin{equation*} L_p(t)=\sqrt{L_1(t)^2+L_2(t)^2},\end{equation*}

(13)   \begin{equation*}Q_0(t)=\begin{bmatrix}\frac{L_1(t)L_3(t)}{LL_p(t)}&\frac{L_2(t)L_3(t)}{LL_p(t)}&-\frac{L_p(t)}{L}\\-\frac{L_2(t)}{L_p(t)}&\frac{L_1(t)}{L_p(t)}&0\\\frac{L_1(t)}{L}&\frac{L_2(t)}{L}&\frac{L_3(t)}{L}\end{bmatrix},\end{equation*}

(14)   \begin{equation*}Q_1(t)=\begin{bmatrix}\cos\psi(t)&-\sin\psi(t)&0\\\sin\psi(t)&\cos\psi(t)&0\\0&0&1\end{bmatrix},\end{equation*}

(15)   \begin{equation*}Q(t)=Q_1(t)Q_0(t),\end{equation*}

Changing signs of any two of three components (L_1,L_2,L_3) we can obtain altogether four different evolutions of a painted T_handle. But that, and also the case of m>1 must wait for the next post.

But here are four pictures of the painted T-handle, rotated according to the above algorithm, at t=10, where the three last cases correspond to flipped sign of (L_1,L_2,L_3) as follows: (1,2),(2,3),(1,3)

Update

You can download my experimental Mathematica cdf file that can be run using free CDF Player. It shows the code and the animation for the four cases mentioned above.

Print 🖨 PDF 📄 eBook 📱
This entry was posted in Attitude matrix, Dzhanibekov effect, Elliptic functions, Mathematica. Bookmark the permalink.

16 Responses to Taming the T-handle continued

  1. Bjab says:
    2017/02/14 at 09:56

        \[\Pi(n;z|m)=\int_0^{F(z|m)}\,\frac{1}{1-n\,\sn(u|m)^2}\. du.\]

    ->

        \[\Pi(n;z|m)=\int_0^{F(z|m)}\,\frac{1}{1-n\,\sn(u|m)^2} du.\]

    (The dot over d is unneeded I suppose. (In many places.))

  2. Bjab says:
    2017/02/14 at 10:10

    From we have ->
    ?

  3. Bjab says:
    2017/02/14 at 10:28

    so e can ->
    so we can

  4. Bjab says:
    2017/02/14 at 10:33

    sign of (L_1,L_2,L3) ->
    sign of (L_1,L_2,L_3)

  5. Bjab says:
    2017/02/14 at 10:34

    In (13): L_p ->
    L_p(t)

  6. Bjab says:
    2017/02/14 at 10:53

    So now there is a hypothesis that:
    velocity of the end of the leg of T-handle is always perpendicular to the plane of T-handle.

  7. Ronan says:
    2017/02/19 at 03:16

    I am struggling a bit with
    {\it c1}\,t+{\frac {{\it c2}\,{\it EllipticPi} \left( -{\it c3},{\it JacobiAM} \left( t,m \right) ,m \right) }{B}} in Maple.
    When I enter the values putting d=49/100
    I get
    t/3-{\frac {20\,\sqrt {51}\sqrt {2}\sqrt {3}}{153}{\it EllipticPi} \left( 3,{\it JacobiAM} \left( t,{\frac {\sqrt {47}\sqrt {51}}{51}} \right) ,{\frac {\sqrt {47}\sqrt {51}}{51}} \right) }
    My m is the square root of your m because of the way Maple handles Jacobi functions.

    Setting t=0 have tested up to 50 the function evaluates to -.82693+3.6654*I

    I have had some success by getting an amimation to run for d=1/2. I will post it when I am more sure about my code/understanding.

  8. Ronan says:
    2017/02/20 at 01:31

    I found this don’t know if it helps. Scroll down to section 2.6. https://www.researchgate.net/publication/262689147_Elliptic_functions_and_elliptic_integrals_for_celestial_mechanics_and_dynamical_astronomy
    mentions Mathematica and its variation on the notation.
    Aside, how do I shorten these links to one word?

Leave a Reply