Eikonal Blog


Vector bases for 3D kinematics

Filed under: mathematics — Tags: — sandokan65 @ 16:02

Relevant sets of basis vectors

There are four relevant sets of orthonormal basis vectors (=triads) here:

  • {\cal B}_1 :\equiv \{\underline{e}_x,\underline{e}_y,\underline{e}_x\} – orts corresponding to the Cartesian coordinates \{x,y,z\},
  • {\cal B}_3 :\equiv \{\underline{n},\underline{p},\underline{q}\}, defined by:
    • \underline{n} :\equiv \frac1{r}\underline{r} = (\cos(\theta)\cos(\varphi),\cos(\theta)\sin(\varphi),\sin(\theta)) ,
    • \underline{p} :\equiv \frac{\partial}{\partial \theta} \underline{n} = (-\sin(\theta)\cos(\varphi),-\sin(\theta)\sin(\phi),\cos(\theta)),
    • \underline{q} :\equiv \frac1{\cos(\theta)}\frac{\partial}{\partial \varphi} \underline{n} = (-\sin(\varphi),\cos(\phi),0),

    Note that here: \underline{n}\times\underline{p} = - \underline{q}, ie. that this basis is left-handed.

  • {\cal B}_3 :\equiv \{\underline{n},\underline{m},\underline{k}\} – radius-based basis, defined by:
    • \underline{n} :\equiv \frac1{r}\underline{r},
    • \underline{m} :\equiv \hbox{Normalized}(\underline{n}'),
    • \underline{k} :\equiv \underline{n}\times\underline{m}.
  • {\cal B}_4 :\equiv \{\underline{T},\underline{O},\underline{B}\} – tangent-based basis (“frame”), defined by:
    • \underline{T} :\equiv \hbox{Normalized}(\underline{r}') = \frac1{v}\underline{v} = \frac{r'}{v} \underline{n} + \frac{r \Phi'}{v} \underline{m},
    • \underline{N} :\equiv \hbox{Normalized}(\underline{T}')
    • \underline{B} :\equiv \underline{T}\times\underline{N}.

    where \underline{v} :\equiv \underline{r}' = r' \underline{n} + r \Phi' \underline{m} and v(t) = \sqrt{r'^2+r^2\Phi'^2}.
    Note that \underline{N} \ne \underline{n}.

Triad kinematic formulas

{\cal B}_2:

  • \underline{n}' = \theta' \underline{p} + \cos(\theta)\varphi' \underline{q},
  • \underline{p}' = -\theta' \underline{n} - \sin(\theta)\varphi' \underline{q},
  • \underline{q}' = -\cos(\theta) \theta' \underline{n} + \sin(\theta)\varphi' \underline{p}.

{\cal B}_3:

  • \underline{n}' = \Phi' \underline{m},
  • \underline{m}' = -\Phi' \underline{n} + \Theta' \underline{k},
  • \underline{k}' = -\Theta' \underline{m}.

Note that \Phi and \Theta are not the same as \varphi and \theta. They are defined by these equations.

A vector \underline{A} = A_n\underline{n}+A_m\underline{m}+A_k\underline{k} evolves according to:
\underline{A}' = [A'_n - \Phi' A_m]\underline{n} + [A'_m + \Phi' A_n - \Theta' A_k]\underline{m} + [A'_k + \Theta' A_m]\underline{k}.

{\cal B}_4 is guided by the Frenet–Serret formulas:

  • \frac{d}{ds}\underline{T} = \kappa \underline{N},
  • \frac{d}{ds}\underline{N} = -\kappa  \underline{T} + \tau \underline{B},
  • \frac{d}{ds}\underline{B} = -\tau \underline{N}.

where s(t) :\equiv \int_0^t dt' v(t').

Relations between triads

{\cal B}_3 in terms of {\cal B}_2:

  • \underline{n} =  \underline{n},
  • \underline{p} = \frac{\theta'}{\Phi'} \underline{m} + \cos(\theta)\frac{\varphi'}{\Phi'} \underline{k},
  • \underline{q} = \cos(\theta) \frac{\varphi'}{\Phi'} \underline{m} - \frac{\theta'}{\Phi'} \underline{k}.

{\cal B}_2 in terms of {\cal B}_3:

  • \underline{n} =  \underline{n},
  • \underline{m} = \frac{\theta'}{\Phi'} \underline{p} + \cos(\theta)\frac{\varphi'}{\Phi'} \underline{q},
  • \underline{k} = \cos(\theta) \frac{\varphi'}{\Phi'} \underline{p} - \frac{\theta'}{\Phi'} \underline{q}.

Note: Here

  • \Phi'^2 :\equiv \theta'^2 + \cos(\theta)^2 \varphi'^2,
  • \Theta' :\equiv \sin(\theta) \varphi'  \left(1+\frac{\theta'^2}{\Phi'^2}\right) + \cos(\theta) \frac{\theta''\varphi'-\theta'\varphi''}{\Phi'^2}.

Note: if we replace t with \Phi, then

  • \left(\frac{d\theta}{d\Phi}\right)^2 + \cos(\theta)^2 \left(\frac{d\varphi}{d\Phi}\right)^2=1,
  • \frac{d\Theta}{d\Phi} = \sin(\theta) \left(\frac{d\varphi}{d\Phi}\right)  \left[1+\left(\frac{d\theta}{d\Phi}\right)^2\right] + \cos(\theta) \left[\frac{d^2\theta}{d\Phi^2}\frac{d\varphi}{d\Phi}-\frac{d\theta}{d\Phi}\frac{d^2\varphi}{d\Phi'^2}\right].


Leave a Comment »

No comments yet.

RSS feed for comments on this post. TrackBack URI

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

Create a free website or blog at WordPress.com.

%d bloggers like this: