next up previous contents
Next: Circular Region: Up: Corkscrew Elements Previous: Back to the Corkscrew...   Contents

Clothoid 1: $(0 \leq s_{proj} \leq a_{1}t_{1_{f}} )$

The equation for the arclength of a clothoid is used to determine $t$, the variable used for clothoid calculations.

\begin{displaymath}
t = \frac{s_{proj}}{a_{1}}
\end{displaymath}

The vertical angle through the corkscrew is found in the same manner as it is for a loop.

\begin{displaymath}
\phi = \frac{\pi t^{2}}{2}
\end{displaymath}

The projection of the corkscrew in the $x$-$z$ plane is determined by the Fresnel Integrals. The $y$ coordinate varies linearly with the arclength of the corkscrew.

\begin{eqnarray*}
x & = & a_{1} \ensuremath\int_{0}^{t} \cos(\frac{\pi}{2} u^{2...
... = & a_{1} \ensuremath\int_{0}^{t} \sin(\frac{\pi}{2} u^{2})\,du
\end{eqnarray*}



The corkscrew is then rotated downward $\lambda$ radians about the $x$ axis. See figure 2.9.
$\displaystyle x$ $\textstyle =$ $\displaystyle a_{1} \ensuremath\int_{0}^{t} \cos(\frac{\pi}{2} u^{2})\,du$ (2.84)
$\displaystyle y$ $\textstyle =$ $\displaystyle s_{proj} \tan\alpha \cos\lambda - a_{1} \sin\lambda \ensuremath\int_{0}^{t} \sin(\frac{\pi}{2} u^{2})\,du$ (2.85)
$\displaystyle z$ $\textstyle =$ $\displaystyle s_{proj} \tan\alpha \sin\lambda + a_{1} \cos\lambda \ensuremath\int_{0}^{t} \sin(\frac{\pi}{2} u^{2})\,du$ (2.86)

The calculation of the forward and radial vectors can be performed two ways. The derivatives of the final position vector (with respect to arclength) may be taken, or the derivative may be taken before the rotation is made. The forward and radial vectors would then be rotated. The latter method will be used here, as it was for a curve.

The calculation of the forward and radial vectors is nearly identical to the calculation for a loop and a curve. An abbreviated description will be given here. For a more complete description, refer to section 2.1.2.

Derivatives of the position vector with respect to arclength must be taken. The required chain rule relationships are first determined.

\begin{eqnarray*}
\ensuremath \frac{dt}{ds} & = & \frac{1}{a_{1}} \cos\alpha \\...
...ensuremath \frac{d\phi}{ds} & = & \frac{\pi t}{a_{1}} \cos\alpha
\end{eqnarray*}



The forward vector is calculated as the derivative of the position vector with respect to arclength.

\begin{eqnarray*}
f_{x} & = & \ensuremath \frac{d\left[ a_{1} \ensuremath\int_{...
...nsuremath \frac{dt}{ds} \nonumber \\
& = & \sin\phi \cos\alpha
\end{eqnarray*}



This vector is now rotated $\lambda$ radians about the $x$ axis.
$\displaystyle f_{x}$ $\textstyle =$ $\displaystyle \cos\phi \cos\alpha$ (2.87)
$\displaystyle f_{y}$ $\textstyle =$ $\displaystyle \sin\alpha \cos\lambda - \sin\phi \cos\alpha \sin\lambda$ (2.88)
$\displaystyle f_{z}$ $\textstyle =$ $\displaystyle \sin\alpha \sin\lambda + \sin\phi \cos\alpha \cos\lambda$ (2.89)

The radial vector is calculated as the derivative of the forward vector with respect to arclength.

\begin{eqnarray*}
r_{x} & = & \ensuremath \frac{d\left( \cos\phi \cos\alpha \ri...
...\phi}{ds} \\
& = & \frac{\pi t}{a_{1}} \cos\phi \cos^{2}\alpha
\end{eqnarray*}



The unit vector in the direction of the radial vector must be determined.

\begin{eqnarray*}
r_{x} & = & -\sin\phi \\
r_{y} & = & 0 \\
r_{z} & = & \cos\phi
\end{eqnarray*}



This vector is now rotated $\lambda$ radians about the $x$ axis.
$\displaystyle r_{x}$ $\textstyle =$ $\displaystyle -\sin\phi$ (2.90)
$\displaystyle r_{y}$ $\textstyle =$ $\displaystyle -\cos\phi \sin\lambda$ (2.91)
$\displaystyle r_{z}$ $\textstyle =$ $\displaystyle \cos\phi \cos\lambda$ (2.92)

The curvature is the magnitude of the (non-unit) radial vector.
\begin{displaymath}
k = \frac{\pi t}{a_{1}} \cos^{2}\alpha
\end{displaymath} (2.93)


next up previous contents
Next: Circular Region: Up: Corkscrew Elements Previous: Back to the Corkscrew...   Contents
Darla Weiss 2000-02-13