Find wrapping angle of helix on a torus

774 Views Asked by At

I need some help in calculating the wrapping angle of a spiral helix wrapped on a torus with constant angle against all the meridians of the torus. The wrapping angle (or the angle measured around and/or against the torus circular cross-section See here a (13,6)tori knot) always remains constant as the helix curve spirals around the torus. In layman's terms means it spirals with the same angle always. if I were to know the arc length of one such turn, my solution will be trivial: $$Arclen = {2\pi R_2\over cos(wrapping\_angle)} .(1)$$ I know this formula is counterintuitive to be true and valid regardless of R1 (torus major radius), but it is. I have an algorithm to test it out and it is true beyond any shadow of a doubt. It took me a while to figure it out, but indeed it does seem to suggest that the helix is on a cylinder, but it's on a curved cylinder (a torus). It doesn't make sense, but the formula is correct. I have a way to test it out. I'm trying to see that the tangent of any point on the curve is constant (derivative=0?? i don't know) so it projects(unwinds) as a straight line in a 2D plane. All i know is that equation 1 is correct. From here i can find my wrapping angle, where R2 is torus minor radius. This formula is valid only for 1 turn spiral helix and ONLY if I have constant wrapping angle. The thing is I don't know the arc length of one turn, but I do know R1, R2 and the azimuth angle of 1 turn (or 1 turn step-angle or the angle measured around torus central axis see it here). I want to find the arc length of 1 turn of the helix as a function of step-angle and using such formula in combination with the one I already have, I can find my wrapping angle. I don't know how to integrate but I see something in my mind along these lines: $$Arclen={ \int_0^{stp} \pi {{ \int_0^{2\pi}(R_1-R_2cos(\theta) )d\theta }\over 180} }\beta d\beta =something .(2)$$ I don't know what this is. I believe it's some sort of a double integration in polar coordinates(?) of infinitesimal circle-arcs which i know nothing about. stp is my step_angle, $\theta$ is the infinitesimal angle made by R2 in respect to equatorial plane of torus, and $\beta$ is supposed to be the step_angle chopped down into infinitesimal bits. I know equation (2) is wrong, but let's say it's not, then i write: $$equation(1)=equation(2) => wrapping\_angle=arccos({{2\pi R_2} \over something})$$ Since equation(2) is wrong, I also see in my mind something like this: $$Arclen=\int_0^{stp} (\sqrt{(\pi [\int_0^{2\pi}(R_1-R_2cos(\theta)) d\theta]*{\beta \over 180})^2+(\int_0^{2\pi}{{\pi R_2\alpha} \over180}d\alpha)^2})d\beta . (3)$$ I don't know how to integrate. Please don't ask me where i got these equations. Maybe they can provide a place to start for you to help me. I just need a little bit of help to correct this mess. Can you help?

2

There are 2 best solutions below

9
On

The curves you have described do not have constant torsion (in the usual technical sense of differential geometry). The condition that the curve make a constant angle with the meridian circles on the torus leads to a very complicated parametric equation. For any value of the constant $c$, here is such a curve on the torus obtained by rotating a circle of radius $1$ around a circle of radius $2$:

$$\alpha(t) = \big((2+\cos t)\cos\theta(t),(2+\cos t)\sin\theta(t),\sin t\big),\quad \text{with }\theta(t)= c\arctan\Big(\dfrac{\tan(t/2)}{\sqrt3}\Big).$$

You could replace $2+\cos t$ with $2+\cos at$ for some constant $a$. I haven't had the fortitude to work everything out in that generality.

Now, for the constant value $c$ appearing in the formula above, we find that the speed of $\alpha(t)$ is $\sqrt{1+c^2}$, and the arclength of one trip around the torus will be $$L=\int_0^T \|\alpha'(t)\|dt = T\sqrt{1+c^2},$$ with $T$ given by $\int_0^T \theta'(t)\,dt = 2\pi$. So this gives us $$2\pi = \theta(T) = c\arctan\Big(\dfrac{\tan(T/2)}{\sqrt3}\Big), \quad\text{so ... drumroll ...} \\ T = 2\arctan\Big(\sqrt3\tan\big(\frac{2\pi}c\big)\Big).$$

Finally, the arclength is $$L=(\sqrt{1+c^2})T = 2\sqrt{1+c^2}\arctan\Big(\sqrt3\tan\big(\frac{2\pi}c\big)\Big).$$ I should add that your wrapping angle is given by $$\cos(\text{wrapping angle}) = \frac1{\sqrt{1+c^2}}.$$

0
On

Ted Shifrin answer to my original question is incomplete and I can't use it the way it is. The almost-complete version is like this:

$$\alpha(t)=\begin{pmatrix} x=\Big(R_1+R_2\cos(t)\Big)\cos\biggr(c\arctan\Big(\dfrac{(R_1-R_2)\tan(t/2)}{\sqrt{R_1^2-R_2^2}}\Big)\biggr) \\y=\Big(R_1+R_2\cos(t)\Big)\cos\biggr(c\arctan\Big(\dfrac{(R_1-R_2)\tan(t/2)}{\sqrt{R_1^2-R_2^2}}\Big)\biggr) \\z=R_2\sin(t) \end{pmatrix}$$

I have injuries for scratching my head trying to more or less do semi-educated guess what extra terms goes where by trail and error (https://en.wikipedia.org/wiki/Trial_and_error) That's because I'm an engineer not a mathematician which should make sense to as why I did it like this.

Using "consecrated" math terminology from here: https://en.wikipedia.org/wiki/Toroidal_and_poloidal (I hope it's ok to use the word "consecrated" in the meaning to "to devote to a purpose with deep dedication" (according to Merriam-Webster dictionary) english is not my native language so I'm asking for leniency)

Let there be the term : "toroidal direction step angle" or the angle done in toroidal direction by only 1 poloidal turn of my curve (a STEP in the right way - literally).

Therefore I can't show mathematically, but have proven undoubtedly in numerous tests and experiments that: $$c=\frac{\text{step angle°}}{180} $$

$$\text{1 poloidal turn wire length}=\frac{2\pi R_2}{\cos(\text{wrapping angle})} $$

$$\text{step angle°}=\frac{360R_2\tan(\text{wrapping angle})}{\sqrt{R_1^2-R_2^2}}$$

$$\text{wrapping angle°}=\arctan\biggr(\frac{\text{step angle}\sqrt{R_1^2-R_2^2}}{360R_2}\biggr)$$ Again, emphasis on how these equations were determined: semi-educated guesses, trial and error a dose of madness and insanity and lots of dedication and perseverance. These formulas never failed me so far for any value of $R_1$ , $R_2$, c, step angle or wrapping angle.

Of course for the answer to be complete, I would have to replace $R_1+R_2\cos(t)$ with $R_1+R_2\cos(at)$ for some constant "a" representing how many turns in the poloidal direction my curve does in the same step angle, the final goal being to generate a special family of closed curves tori knots (with curve tangent(aka velocity vector) always constant against all meridians and parallels of the torus)

I was asked to provide more context to my question: https://www.researchgate.net/publication/311742883_Visualization_in_3D_of_Dynamics_of_Toroidal_Helical_Coils_in_quest_of_optimum_designs_for_a_Concordian_Mandala

I would like to add some sentences (which most likely will be deleted by someone):

  1. Nothing less than a miracle or a fluke for a professor like Ted Shifrin answering my question.

  2. My question was written with sloppy terminology.

  3. It would have been amazing to have seen the math he did behind the scenes for the parametrization of such curves.

  4. Can I hope for yet another miracle?

Why is speaking like this forbidden on this forum? Reading between the lines 1st sentence shows gratitude towards altruism, respect and awe. 2nd- introspection, constructive self-criticism and honesty. 3rd mild regret, call for compassion and desire. and 4th expectation and hope Why speaking (in and with) these terms illegal on this forum? Aren't these aspects extremely important in any form of communication? As a newbie on this forum I'm very confused... Yes... stick to the math and disregard and delete everything else...monstrously strict and rigid.