Difference between revisions of "Motion in Space"

From Department of Mathematics at UTSA
Jump to navigation Jump to search
(Created page with "* [https://openstax.org/books/calculus-volume-3/pages/3-4-motion-in-space Motion in Space], openstax")
 
 
(6 intermediate revisions by 2 users not shown)
Line 1: Line 1:
* [https://openstax.org/books/calculus-volume-3/pages/3-4-motion-in-space Motion in Space], openstax
+
== Curvature ==
 +
 
 +
=== Terminology ===
 +
Before we start discussing curvature, there are some important vectors and concepts we need to be at least aware of.
 +
 
 +
==== The unit tangent vector ====
 +
In the differentiation section of this chapter, we discussed the derivatives of a vector function. We know that <math>\mathbf{v}(t)=\mathbf{r}'(t)</math> at <math>t=t_0</math> is tangent to the curve <math>\mathbf{r}(t)</math> at <math>t=t_0</math>. <math>\mathbf{r}'(t)</math> is called the tangent vector. The unit tangent vector, however, eliminates the aspect of magnitude because it is defined as:<blockquote><math>\mathbf{T}(t)=\frac{\mathbf{r}'(t)}{|\mathbf{r}'(t)|}</math></blockquote>As we can see, the magnitude of the unit tangent vector is always <math>1</math>. We can imagine that <math>\mathbf{r}(t)</math> as the displacement of a particle with respect to time. So, the unit tangent vector can be perceived as the direction of the velocity of the particle with respect to time. It can also be perceived as the direction of the tangential acceleration of the particle with respect to time. We will discuss motion in space in the next section, but this is a useful method to intuitively understand some vectors.
 +
 
 +
==== The unit normal vector ====
 +
The unit normal vector is defined as<blockquote><math>\mathbf{N}(t)=\frac{\mathbf{T}'(t)}{|\mathbf{T}'(t)|}</math></blockquote>The unit normal is orthogonal to the unit tangent because since <math>|\mathbf{T}(t)|=1</math>, we can get that:<blockquote><math>\frac{d}{dt}|\mathbf{T}(t)|^2=0=\frac{d}{dt}[\mathbf{T}(t)\cdot\mathbf{T}(t)]=2\mathbf{T}'(t)\cdot\mathbf{T}(t)</math>
 +
 
 +
 
 +
<math>\Leftrightarrow \mathbf{T}'(t)\cdot\mathbf{T}(t)=0</math></blockquote>This means that <math>\mathbf{T}'(t)</math> is orthogonal to <math>\mathbf{T}(t)</math>. Therefore, <math>\mathbf{N}(t)</math> is orthogonal to <math>\mathbf{T}(t)</math>. We can imagine that the unit normal vector is the direction of the normal acceleration of the particle with respect to time.
 +
 
 +
==== The binormal vector ====
 +
The binormal vector is defined as<blockquote><math>\mathbf{B}(t)=\mathbf{T}(t)\times\mathbf{N}(t)</math></blockquote>The binormal vector is perpendicular to both the unit tangent and the unit normal because of the properties of the cross product. The magnitude of the binormal is always 1 because<blockquote><math>|\mathbf{B}(t)|=|\mathbf{T}(t)||\mathbf{N}(t)|\sin\theta=1</math></blockquote>
 +
 
 +
==== The normal plane, osculating plane, and the osculating circle ====
 +
 
 +
* The normal plane is the plane determined by the normal and binormal vectors <math>\mathbf{N}\text{ and }\mathbf{B}</math>. The normal plane consists of all lines that are orthogonal to the tangent vector <math>\mathbf{T}</math>.
 +
* The osculating plane is the plane determined by the unit tangent and unit normal <math>\mathbf{T}\text{ and }\mathbf{N}</math>. It is the plane that comes closest to containing the part of the curve near a point where <math>t=t_0</math>.
 +
* The osculating circle is the circle that lies in the osculating plane towards the direction of <math>\mathbf{N}</math> with a radius <math>r=\frac{1}{\kappa}</math> (the inverse of the curvature, which we will immediately discuss after this). It best describes how the curve behaves near the point where <math>t=t_0</math> because it shares the same tangent, normal, and curvature at that point.
 +
These concepts are very important in the branch of differential geometry and in its applications to the motion of spacecraft.
 +
 
 +
=== Curvature ===
 +
The curvature of a curve at a given point is a measure of how quickly the curve changes direction at that point. We define it to be the magnitude of the rate of change of the unit tangent with respect to arc length. We use arc length so that the curvature will be independent of parametrization.<blockquote>Suppose that a space curve has vector function <math>\mathbf{r}</math>, unit tangent vector <math>\mathbf{T}</math>, and arc length <math>s</math>. The curvature of this curve is: <math>\kappa=\biggl|\frac{d\mathbf{T}(s)}{ds}\biggr|</math>.</blockquote>There are two other ways to express the curvature. We can express curvature in terms of <math>t</math> instead of <math>s</math> by utilizing the chain rule (recall that <math>\frac{dL}{dt}=|\mathbf{r}'(t)|</math>):<blockquote><math>\kappa=\biggl|\frac{d\mathbf{T}}{ds}\biggr|=\Biggl|\frac{\frac{d\mathbf{T}}{ds}\frac{ds}{dt}}{\frac{ds}{dt}}\Biggr|=\Biggl|\frac{\frac{d\mathbf{T}(t)}{dt}}{\frac{ds}{dt}}\Biggr|=\frac{|\mathbf{T}'(t)|}{|\mathbf{r}'(t)|}</math></blockquote>The third way is more complicated to deduce, but it is often more convenient to apply because it only requires <math>\mathbf{r}(t)</math> and its derivatives.<blockquote><math>\kappa(t)=\frac{|\mathbf{r}'(t)\times\mathbf{r}''(t)|}{|\mathbf{r}'(t)|^3}</math></blockquote>And now for the proof for this theorem:<blockquote>According to the definition of the unit tangent vector, we know that <math>\mathbf{r}'=|\mathbf{r}'|\mathbf{T}</math>. So the second derivative of <math>\mathbf{r}</math> should be:<blockquote><math>\begin{align}
 +
\mathbf{r}'' & = (|\mathbf{r}'|\ \mathbf{T})' \\
 +
& = |\mathbf{r}''|\ \mathbf{T}+|\mathbf{r}'|\ \mathbf{T}' \quad\text{the product rule}\\
 +
\end{align}</math></blockquote>Now we calculate <math>\mathbf{r}'\times\mathbf{r}''</math>.<blockquote><math>\begin{align}
 +
\mathbf{r}'\times\mathbf{r}'' & = \bigl(|\mathbf{r}'|\mathbf{T}\bigr)\times\bigl(|\mathbf{r}''|\ \mathbf{T}+|\mathbf{r}'|\ \mathbf{T}' \bigr)\quad\text{substitution} \\
 +
& = \bigl(|\mathbf{r}'|\mathbf{T}\times|\mathbf{r}''|\mathbf{T}\bigr)+\bigl(|\mathbf{r}'|\mathbf{T}\times|\mathbf{r}'|\mathbf{T}'\bigr)\quad\text{distribution} \\
 +
& = |\mathbf{r}'||\mathbf{r}''|(\mathbf{T}\times\mathbf{T})+|\mathbf{r}'|^2(\mathbf{T}\times\mathbf{T}')\quad\text{rearrangement} \\
 +
& = |\mathbf{r}'|^2(\mathbf{T}\times\mathbf{T}') \quad\text{realizing that }\mathbf{T}\times\mathbf{T}=0 \\
 +
\end{align}</math></blockquote>And then we calculate <math>|\mathbf{r}'\times\mathbf{r}''|</math>.<blockquote><math>\begin{align}
 +
|\mathbf{r}'\times\mathbf{r}''| & = |\mathbf{r}'|^2|\mathbf{T}\times\mathbf{T}'|\quad\text{substitution} \\
 +
& = |\mathbf{r}'|^2|\mathbf{T}||\mathbf{T}'|\sin\theta \quad\text{the magnitude for the cross product} \\
 +
& = |\mathbf{r}'|^2|\mathbf{T}'| \quad\text{when you realize }|\mathbf{T}|=1\text{ and }\mathbf{T}\perp\mathbf{T}' \\
 +
\end{align}</math></blockquote>We rearrange the equation into:<blockquote><math>|\mathbf{T}'|=\frac{|\mathbf{r}'\times\mathbf{r}''|}{|\mathbf{r}'|^2}</math></blockquote>Since <math>\kappa=\frac{|\mathbf{T}'|}{|\mathbf{r}'|}</math>, we can substitute <math>|\mathbf{T}'|</math> with <math>\frac{|\mathbf{r}'\times\mathbf{r}''|}{|\mathbf{r}'|^2}</math> and get: <math>\kappa=\frac{|\mathbf{r}'\times\mathbf{r}''|}{|\mathbf{r}'|^3}</math></blockquote>Here is a little summary on ways to calculate the curvature.<center>
 +
{| class="wikitable"
 +
!Definition
 +
!With parametrization with respect to <math>\mathbf{t}</math> 
 +
!In terms of <math>\mathbf{r}(t)</math> and its derivatives
 +
|-
 +
|align=center|<math>\kappa=\biggl|\frac{d\mathbf{T}}{ds}\biggr|</math>
 +
|align=center|<math>\kappa=\frac{|\mathbf{T}'|}{|\mathbf{r}'|}</math>
 +
|align=center|<math>\kappa=\frac{|\mathbf{r}'\times\mathbf{r}''|}{|\mathbf{r}'|^3}</math>
 +
|}</center>
 +
 
 +
== Motion in space ==
 +
 
 +
=== Velocity and acceleration ===
 +
Remember in 2-dimensional calculus, we mentioned that a particle with displacement function <math>f(x)</math> has velocity <math>f'(x)</math> and acceleration <math>f''(x)</math>. In vector functions, the definition is basically the same. Suppose a particle moves through space so that its position vector at time <math>t</math> is <math>\mathbf{r}(t)</math>, its velocity function and acceleration function are:<blockquote><math>\mathbf{v}(t)=\lim_{h\rightarrow 0}\frac{\mathbf{r}(t+h)-\mathbf{r}(t)}{h}=\mathbf{r}'(t)\quad</math> and <math>\quad\mathbf{a}(t)=\lim_{h\rightarrow 0}\frac{\mathbf{v}(t+h)-\mathbf{v}(t)}{h}=\mathbf{v}'(t)=\mathbf{r}''(t)</math></blockquote>Simply put: <math>\mathbf{r}''(t)=\mathbf{v}'(t)=\mathbf{a}(t)</math>.
 +
 
 +
 
 +
The speed of the particle ignores the direction. It is the magnitude of the velocity vector: <math>|\mathbf{v}(t)|</math>. The distance traveled by the particle from <math>t=a \text{ to } b</math> is <math>\int_a^b|\mathbf{v}(t)|dt</math>, which is also the formula for the arc length.
 +
 
 +
 
 +
With the help of the Fundamental Theorem of Calculus, we can deduce the velocity function and position function when we know the particle's acceleration.<blockquote><math>\mathbf{v}(t)=\mathbf{v}(t_0)+\int_{t_0}^t\mathbf{a}(u)du\quad</math> and <math>\quad\mathbf{r}(t)=\mathbf{r}(t_0)+\int_{t_0}^t\mathbf{v}(u)du</math></blockquote>
 +
 
 +
=== Tangential and normal acceleration ===
 +
We can split the acceleration vector into two components: the tangential acceleration <math>a_T</math> and the normal acceleration <math>a_N</math>. The tangential acceleration faces the same direction as the unit tangent vector (<math>\mathbf{T}</math>), and the normal acceleration faces the same direction as the unit normal vector (<math>\mathbf{N}</math>). Since both <math>\mathbf{T}</math> and <math>\mathbf{N}</math> are unit vectors, the acceleration vector can be written as the sum of two vectors:<blockquote><math>\mathbf{a}=a_T\mathbf{T}+a_N\mathbf{N}</math></blockquote>Our goal is to figure out how to describe the two components.<blockquote>Recall that <math>\mathbf{T}=\frac{\mathbf{v}}{|\mathbf{v}|}</math>, thus<blockquote><math>\mathbf{v}=|\mathbf{v}|\mathbf{T}</math></blockquote>Now we differentiate both sides of the equation,<blockquote><math>\begin{align}
 +
\mathbf{v}' & = \bigl(|\mathbf{v}|\mathbf{T}\bigr)' \\
 +
& = |\mathbf{v}|'\mathbf{T}+|\mathbf{v}|\mathbf{T}'\\
 +
\end{align}</math></blockquote>So we get <math>\mathbf{a}=|\mathbf{v}|'\mathbf{T}+|\mathbf{v}|\mathbf{T}'</math>. <blockquote>Recall that <math>\kappa=\frac{\big|\mathbf{T}'\big|}{\big|\mathbf{r}'\big|}=\frac{\big|\mathbf{T}'\big|}{\big|\mathbf{v}\big|}</math>, thus <math>\big|\mathbf{T}'\big|=\kappa|\mathbf{v}|</math>.</blockquote><blockquote>Recall that <math>\mathbf{N}=\frac{\mathbf{T}'}{\big|\mathbf{T}'\big|}</math>, thus <math>\mathbf{T}'=\big|\mathbf{T}'\big|\mathbf{N}=\kappa|\mathbf{v}|\mathbf{N}</math>.</blockquote>We substitute <math>\mathbf{T}'</math> for <math>\kappa|\mathbf{v}|\mathbf{N}</math> to yield: <math>\mathbf{a}=|\mathbf{v}|'\mathbf{T}+\kappa|\mathbf{v}|^2\mathbf{N}</math></blockquote>That leaves us with:<blockquote><math>a_T=|\mathbf{v}|'\quad</math> and <math>\quad a_N=\kappa|\mathbf{v}|^2</math></blockquote>
 +
 
 +
<math>\blacksquare</math>
 +
 
 +
Of course, it will be more convenient if those components can be written in terms if <math>\mathbf{r}(t)</math> and its derivatives. Suppose that <math>\theta</math> is the angle between <math>\mathbf{T}</math> and <math>\mathbf{a}</math>, then we can write <math>a_T,a_N</math> like this:<blockquote><math>\begin{align}
 +
a_T & = |\mathbf{a}(t)|\cos\theta \\
 +
& = \frac{|\mathbf{v}(t)||\mathbf{a}(t)|\cos\theta}{|\mathbf{v}(t)|} \quad\text{algebraic manipulation}\\
 +
& = \frac{\mathbf{v}(t)\cdot\mathbf{a}(t)}{|\mathbf{v}(t)|}\quad\text{dot product}\\
 +
& = \frac{\mathbf{r}'(t)\cdot\mathbf{r}''(t)}{|\mathbf{r}'(t)|}\quad\text{in terms of }\mathbf{r}(t)\text{ and its derivatives}\\
 +
\end{align}</math><math>\begin{align}
 +
a_N & = |\mathbf{a}(t)|\sin\theta \\
 +
& = \frac{|\mathbf{v}(t)||\mathbf{a}(t)|\sin\theta}{|\mathbf{v}(t)|} \quad\text{algebraic manipulation}\\
 +
& = \frac{|\mathbf{v}(t)\times\mathbf{a}(t)|}{|\mathbf{v}(t)|}\quad\text{cross product}\\
 +
& = \frac{|\mathbf{r}'(t)\times\mathbf{r}''(t)|}{|\mathbf{r}'(t)|}\quad\text{in terms of }\mathbf{r}(t)\text{ and its derivatives}\\
 +
\end{align}</math></blockquote>
 +
 
 +
<math>\blacksquare</math>
 +
 
 +
To sum it up:<blockquote><math>a_T=|\mathbf{v}|'=\frac{\mathbf{r}'(t)\cdot\mathbf{r}''(t)}{|\mathbf{r}'(t)|}\quad</math> and <math>\quad a_N=\kappa|\mathbf{v}|^2= \frac{|\mathbf{r}'(t)\times\mathbf{r}''(t)|}{|\mathbf{r}'(t)|}</math></blockquote>
 +
 
 +
 
 +
==Resources==
 +
* [https://en.wikibooks.org/wiki/Calculus/Vector_Functions Vector Functions], Wikibooks: Calculus
 +
* [https://openstax.org/books/calculus-volume-3/pages/3-4-motion-in-space Motion in Space], OpenStax
 +
 
 +
==Licensing==
 +
Content obtained and/or adapted from:
 +
* [https://en.wikibooks.org/wiki/Calculus/Vector_Functions Vector Functions, Wikibooks: Calculus] under a CC BY-SA license

Latest revision as of 10:17, 2 November 2021

Curvature

Terminology

Before we start discussing curvature, there are some important vectors and concepts we need to be at least aware of.

The unit tangent vector

In the differentiation section of this chapter, we discussed the derivatives of a vector function. We know that at is tangent to the curve at . is called the tangent vector. The unit tangent vector, however, eliminates the aspect of magnitude because it is defined as:

As we can see, the magnitude of the unit tangent vector is always . We can imagine that as the displacement of a particle with respect to time. So, the unit tangent vector can be perceived as the direction of the velocity of the particle with respect to time. It can also be perceived as the direction of the tangential acceleration of the particle with respect to time. We will discuss motion in space in the next section, but this is a useful method to intuitively understand some vectors.

The unit normal vector

The unit normal vector is defined as

The unit normal is orthogonal to the unit tangent because since , we can get that:


This means that is orthogonal to . Therefore, is orthogonal to . We can imagine that the unit normal vector is the direction of the normal acceleration of the particle with respect to time.

The binormal vector

The binormal vector is defined as

The binormal vector is perpendicular to both the unit tangent and the unit normal because of the properties of the cross product. The magnitude of the binormal is always 1 because

The normal plane, osculating plane, and the osculating circle

  • The normal plane is the plane determined by the normal and binormal vectors . The normal plane consists of all lines that are orthogonal to the tangent vector .
  • The osculating plane is the plane determined by the unit tangent and unit normal . It is the plane that comes closest to containing the part of the curve near a point where .
  • The osculating circle is the circle that lies in the osculating plane towards the direction of with a radius (the inverse of the curvature, which we will immediately discuss after this). It best describes how the curve behaves near the point where because it shares the same tangent, normal, and curvature at that point.

These concepts are very important in the branch of differential geometry and in its applications to the motion of spacecraft.

Curvature

The curvature of a curve at a given point is a measure of how quickly the curve changes direction at that point. We define it to be the magnitude of the rate of change of the unit tangent with respect to arc length. We use arc length so that the curvature will be independent of parametrization.

Suppose that a space curve has vector function , unit tangent vector , and arc length . The curvature of this curve is: .

There are two other ways to express the curvature. We can express curvature in terms of instead of by utilizing the chain rule (recall that ):

The third way is more complicated to deduce, but it is often more convenient to apply because it only requires and its derivatives.

And now for the proof for this theorem:

According to the definition of the unit tangent vector, we know that . So the second derivative of should be:

Now we calculate .

And then we calculate .

We rearrange the equation into:

Since , we can substitute with and get:

Here is a little summary on ways to calculate the curvature.

Definition With parametrization with respect to In terms of and its derivatives

Motion in space

Velocity and acceleration

Remember in 2-dimensional calculus, we mentioned that a particle with displacement function has velocity and acceleration . In vector functions, the definition is basically the same. Suppose a particle moves through space so that its position vector at time is , its velocity function and acceleration function are:

and

Simply put: .


The speed of the particle ignores the direction. It is the magnitude of the velocity vector: . The distance traveled by the particle from is , which is also the formula for the arc length.


With the help of the Fundamental Theorem of Calculus, we can deduce the velocity function and position function when we know the particle's acceleration.

and

Tangential and normal acceleration

We can split the acceleration vector into two components: the tangential acceleration and the normal acceleration . The tangential acceleration faces the same direction as the unit tangent vector (), and the normal acceleration faces the same direction as the unit normal vector (). Since both and are unit vectors, the acceleration vector can be written as the sum of two vectors:

Our goal is to figure out how to describe the two components.

Recall that , thus

Now we differentiate both sides of the equation,

So we get .

Recall that , thus .

Recall that , thus .

We substitute for to yield:

That leaves us with:

and

Of course, it will be more convenient if those components can be written in terms if and its derivatives. Suppose that is the angle between and , then we can write like this:

To sum it up:

and


Resources

Licensing

Content obtained and/or adapted from: