下の記事の続きです
前の記事の最後に内円の中にさらに円がある内々トロコイドの動画を紹介しました。
本記事ではその内々トロコイドの方程式を求め方を説明します。
方程式の求め方
基準となる位置を定めるため図1を初期位置とします。外円($O_f$)に対する内円を円1($O_1$)、円1の中の円を円2($O_2$)と呼びます。
図1円1が外円に沿って転がり、円2が外円に沿って転がると図2のようになります。
図2初期位置で外円と円1が接していた点の内、外円の円周上の点を$P_f$、円1の円周上の点を$P_1$、$\angle O_2 O_fP_f$を$\theta_1$、$\angle O_2 O_1 P_1$を$\theta_2$、外円、円1、円2の半径をそれぞれ$R_f, R_1, R_2$とします。また、前の記事と同様に横軸が実部、縦軸が虚軸とした複素平面で考えます。
このとき、$O_1$の位置は
\begin{eqnarray}
O_1 = (R_f – R_1)\exp(\theta_1 i)
\end{eqnarray}
(1)
$O_2$の位置は$\theta_2=0$の場合、前記事の点$D$と等価なため、前記事の式(5)より
\begin{eqnarray}
O_2 &=& O_1 + (R_1 – R_2) \exp((1 – \frac{R_f}{R_1}) \theta_1 i) \\
\end{eqnarray}
(2)
$\theta_2 \ne 0$の場合、これに$\theta_2$分が加わるので、
\begin{eqnarray}
O_2 &=& O_1 + (R_1 – R_2) \exp(((1 – \frac{R_f}{R_1}) \theta_1 + \theta_2) i) \\
\end{eqnarray}
(3)
点$D$の位置は、円2が転がらないで移動したと仮定した場合、
\begin{eqnarray}
D &=& O_2 + R_d \exp(((1 – \frac{R_f}{R_1}) \theta_1 + \theta_2) i) \\
\end{eqnarray}
(4)
となります。
この式に円2の転がりを考慮します。
円2が図2まで転がったとき、円1と円2が接した長さは$\theta’=\angle AO_2D$とすると、
\begin{eqnarray}
R_1 \theta_2 = R_2 \theta’
\end{eqnarray}
(5)
これより$\theta’$は
\begin{eqnarray}
\theta’ = \frac{R_1}{R_2} \theta_2
\end{eqnarray}
(6)
となり、$\theta’$は$\theta_2$と反対方向に回転するため、点Dの位置は式(4)、(6)より
\begin{eqnarray}
D &=& O_2 + R_d \exp(((1 – \frac{R_f}{R_1}) \theta_1 + \theta_2 – \theta’)i) \\
&=& O_2 + R_d \exp(((1 – \frac{R_f}{R_1}) \theta_1 + (1 – \frac{R_1}{R_2})\theta_2)i) \\
\end{eqnarray}
(7)
となります。
点$D$が初期位置に戻る条件
曲線の始点と終点を繋げたい場合は点$D$が初期位置に戻るように$\theta_1, \theta_2$を定める必要があります。本節では点$D$が初期位置に戻るための$\theta_1,\theta_2$の条件を考察します。
円1が初期位置に戻る条件は
\begin{eqnarray}
\theta_1 = 2 \pi N_1 \\
(N_1 \in\mathbb{Z})
\end{eqnarray}
(8)
このとき、円2が元の位置に戻る条件は式(4)より
\begin{eqnarray}
(1-\frac{R_f}{R_1}) \theta_1 + \theta_2 &=& 2 \pi N_2 \\
(N_2 \in\mathbb{Z})
\end{eqnarray}
(9)
式(8)を代入して
\begin{eqnarray}
(1-\frac{R_f}{R_1}) \theta_1 + \theta_2 &=& 2 \pi N_2 \\
\end{eqnarray}
(10)
$\theta_2$の式に直すと、
\begin{eqnarray}
\theta_2 = (\frac{R_f}{R_1}-1) 2 \pi N_1 + 2 \pi N_2\\
\end{eqnarray}
(11)
このとき点$D$が初期位置に戻る条件は式(7)より
\begin{eqnarray}
(1 – \frac{R_f}{R_1}) \theta_1 + (1 – \frac{R_1}{R_2})\theta_2 = 2\pi N_d
\end{eqnarray}
(12)
式(8)、(11)を代入して
\begin{eqnarray}
2\pi N_d &=& (1 – \frac{R_f}{R_1}) 2\pi N_1 + (1 – \frac{R_1}{R_2})(-(1-\frac{R_f}{R_1}) 2 \pi N_1 +2 \pi N_2)\\
&=& (1 – 1 + \frac{R_1}{R_2})(1 – \frac{R_f}{R_1}) 2\pi N_1 + (1 – \frac{R_1}{R_2})2 \pi N_2 \\
&=& \frac{R_1 – R_f}{R_2} 2\pi N_1 + (1 – \frac{R_1}{R_2})2 \pi N_2 \\
(N_d \in\mathbb{Z})
\end{eqnarray}
(13)
よって、点$D$が初期位置に戻る条件は
\begin{eqnarray}
((R_1-R_f)N_1 -R_1 N_2)\mod R_2 = 0 \\
N_1, N_2 \in \mathbb{Z}\\
\end{eqnarray}
(14)
であり、このときの$\theta_1,\theta_2$は
\begin{eqnarray}
\theta_1 &=& 2 \pi N_1 \\
\theta_2 &=& (\frac{R_f}{R_1}-1) 2 \pi N_1 + 2 \pi N_2\\
\end{eqnarray}
(15)
となります。
例として $R_f=11, R_1=5, R_2=2, N_1=3$の場合
\begin{eqnarray}
((5-11)3 -5*N_2)\ &=& (-18 – 5 N_2) \mod 2 = 0\\
\end{eqnarray}
(16)
より、$N_2$が偶数の場合に条件を満たすことが分かります。
ここで$N_2=14$とすると図3のようになり冒頭紹介した動画のサムネの曲線になります。
図3おわりに
内々トロコイドの方程式を見ると規則的な式になっているので、N階の内々トロコイドも計算できそうです。
ただ、初期位置に戻るように描こうとするとNが増えるほどに条件を求めるのが難しくなりそうなので、Nを大きくする場合は始点と終点を繋げずに描くのが良さそうです。
終わり。
