Taming the T-handle continued

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.

16 thoughts on “Taming the T-handle continued

  1. 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.

    1. I noticed two things: you have opposite signs where c2 and c3 are. Perhaps you reversed the order of I1,I2,I3?
      There should be + sign after t/3 and there should be -3 under EllipticPi.

      1. The + and – reversal is Maple. It states in the FunctionAdvisor EllipticPi(-A,B,C)= -EllipticPi(A,B,C). So it swaps the signs itself. It will look further into it tonight. Would you be able to test in your Maple?

        1. The conversion is not that simple. I yet have to work it out, and then test. I think I know how to do it. And since I can’t find it online, I may devote the whole next post to this technical problem.

Leave a Reply