下の記事の続きです
前の記事の最後に内円の中にさらに円がある内々トロコイドの動画を紹介しました。
本記事ではその内々トロコイドの方程式を求め方を説明します。
方程式の求め方
基準となる位置を定めるため図1を初期位置とします。外円(Of)に対する内円を円1(O1)、円1の中の円を円2(O2)と呼びます。

円1が外円に沿って転がり、円2が外円に沿って転がると図2のようになります。

初期位置で外円と円1が接していた点の内、外円の円周上の点を Pf 、円1の円周上の点を P1 、 ∠O2OfPf を θ1 、 ∠O2O1P1 を θ2 、外円、円1、円2の半径をそれぞれ Rf,R1,R2 とします。また、前の記事と同様に横軸が実部、縦軸が虚軸とした複素平面で考えます。
このとき、 O1 の位置は
O1=(Rf−R1)exp(θ1i)(1)
O2の位置はθ2=0の場合、前記事の点Dと等価なため、前記事の式(5)より
O2=O1+(R1−R2)exp((1−R1Rf)θ1i)(2)
θ2=0 の場合、これに θ2 分が加わるので、
O2=O1+(R1−R2)exp(((1−R1Rf)θ1+θ2)i)(3)
点 D の位置は、円2が転がらないで移動したと仮定した場合、
D=O2+Rdexp(((1−R1Rf)θ1+θ2)i)(4)
となります。
この式に円2の転がりを考慮します。
円2が図2まで転がったとき、円1と円2が接した長さは θ′=∠AO2D とすると、
R1θ2=R2θ′(5)
これより θ′ は
θ′=R2R1θ2(6)
となり、 θ′ は θ2 と反対方向に回転するため、点Dの位置は式(4)、(6)より
D==O2+Rdexp(((1−R1Rf)θ1+θ2−θ′)i)O2+Rdexp(((1−R1Rf)θ1+(1−R2R1)θ2)i)(7)
となります。
点Dが初期位置に戻る条件
曲線の始点と終点を繋げたい場合は点 D が初期位置に戻るように θ1,θ2 を定める必要があります。本節では点Dが初期位置に戻るための θ1,θ2 の条件を考察します。
円1が初期位置に戻る条件は
θ1=2πN1(N1∈Z)(8)
このとき、円2が元の位置に戻る条件は式(4)より
(1−R1Rf)θ1+θ2(N2∈Z)=2πN2(9)
式(8)を代入して
(1−R1Rf)θ1+θ2=2πN2(10)
θ2 の式に直すと、
θ2=(R1Rf−1)2πN1+2πN2(11)
このとき点Dが初期位置に戻る条件は式(7)より
(1−R1Rf)θ1+(1−R2R1)θ2=2πNd(12)
式(8)、(11)を代入して
2πNd(Nd∈Z)===(1−R1Rf)2πN1+(1−R2R1)(−(1−R1Rf)2πN1+2πN2)(1−1+R2R1)(1−R1Rf)2πN1+(1−R2R1)2πN2R2R1−Rf2πN1+(1−R2R1)2πN2(13)
よって、点Dが初期位置に戻る条件は
((R1−Rf)N1−R1N2)modR2=0N1,N2∈Z(14)
であり、このときの θ1,θ2 は
θ1θ2==2πN1(R1Rf−1)2πN1+2πN2(15)
となります。
例として Rf=11,R1=5,R2=2,N1=3 の場合
((5−11)3−5∗N2) =(−18−5N2)mod2=0(16)
より、 N2 が偶数の場合に条件を満たすことが分かります。
ここで N2=14 とすると図3のようになり冒頭紹介した動画のサムネの曲線になります。

おわりに
内々トロコイドの方程式を見ると規則的な式になっているので、N階の内々トロコイドも計算できそうです。
ただ、初期位置に戻るように描こうとするとNが増えるほどに条件を求めるのが難しくなりそうなので、Nを大きくする場合は始点と終点を繋げずに描くのが良さそうです。
終わり。