Motion in Space

This worksheet is a continuation of ParaCurves , and focuses on the motion of a particle along a parametric curve

x ( t ) = x ( t ) i + y ( t ) j + z ( t ) k = ( x ( t ), y ( t ), z ( t )).

in space. Consider, for instance, the first curve in ParaCurves , the circular helix

x ( t ) = ( ).

Its veclocity and acceleration vectors are

v ( t ) = ( ) and a ( t ) = ( ).

Since for any value of
t the speed of a particle moving along the curve is

||
v ( t ) || = = ,

the formula for the unit tangent vector at any time t is

T ( t ) = ( ).

To find the unit normal vector N ( t ), first differentiate T ( t ), and then divide by its length, which here is simply the constant . The calculation therefore yields

[ T ( t )] = ( ), so N ( t ) = i j .

The following Maple routine plots the helix, along with its unit tangent and normal vectors at the point corresponding to
t = /2. It draws the unit tangent vector by plotting the segment of the tangent line

x = x ( t ) = x ( /2) + t T ( /2) = (0, 1, /2) + t ( )

at the point
P (0, 1, /2) in the direction of the unit tangent vector T = ( ) for t ranging from 0 to 1. ( P corresponds to t = 0, and t = 1 gives the end of T ( /2) drawn from initial point P .) It uses the same approach for N ( /2): starting from P it plots the normal line

x = x ( t ) = x ( /2) + t N ( /2) = (0, 1, /2) + t ( ).

for
t between 0 and 1. (Again, the unit normal vector ends at the point corresponding to t = 1.)

> with (plots):
curve := spacecurve( [cos(t), sin(t), t], t = 0..2*Pi, color = black, axes=boxed, numpoints = 150, labels = ["x", "y", "z"] ):
utan := spacecurve( [-t/sqrt(2), 1, Pi/2 + t/sqrt(2)], t = 0..1, color = blue, axes = boxed, numpoints = 150 ):
unorm := spacecurve( [0, 1 - t, Pi/2], t = 0..1, color = green, axes = boxed, numpoints = 150 ):
xaxis := spacecurve( [t, 0, 0, t = -1..1, color = magenta] ):
yaxis := spacecurve( [0, t, 0, t = -1..1, color = magenta] ):
zaxis := spacecurve( [0, 0, t, t = -1..6, color = magenta] ):
labx := textplot3d([1.2, 0, -.2, `x`], color = magenta):
laby := textplot3d([0,1.2, -.2, `y`], color = magenta):
labz := textplot3d([0, 0, 6.5, `z`], color = magenta):
display(curve, utan, unorm, xaxis, yaxis, zaxis, labx, laby, labz);

The Moving Trihedral. The binormal vector has formula B ( t ) = T ( t ) ´ N ( t ). For each value of t , the three vectors T ( t ), N ( t ) and B ( t ) define a local right-handed coordinate system with base point ( x ( t ), y ( t ), z ( t )). This leads to the term moving trihedral for the array of these three vectors emanating from the point x (t). For the circular helix x ( t ) = ( ) , the binormal vector is

B ( t ) = ( )

=
(1, 0, 1) for t = /2.

The next routine uses the technique of the preceding one to generate a plot that includes
B ( /2). It thus gives shows the moving trihedral at the point x ( /2) = (0, 1, /2).

> with (plots):
curve := spacecurve( [cos(t), sin(t), t], t = 0..Pi, color = black, axes=boxed, numpoints = 150, labels = ["x", "y", "z"] ):
utan := spacecurve( [-t/sqrt(2), 1, Pi/2 + t/sqrt(2)], t = 0..1, color = blue, axes = boxed, numpoints = 150 ):
unorm := spacecurve( [0, 1 - t, Pi/2], t = 0..1, color = green, axe BookFinder.com: Search for New, Out of Print and Used Books

 Search | About | Interact | Help
 The Open Marketplace for Books Online Books in other languages...EnglishFrenchGermanItalian

Tangental and Normal Components of Acceleration. To decompose a (1) into its tangential and normal components at the point P (1, 1, 1) corresponding to t = 1, just use the foregoing calculaations:

(1) = [ (a á v)/(v á v)] v = [22/14] (1, 2, 3) = (11/7, 22/7, 33/7).

Then (with pointed brackets enclosing the coordinates of vectors to accommodate Maple's syntactical requirements in mathematical expressions),

(1) = (1) = = .

From this, it follows that

N (1) = [ 1/|| (1) || ] (1) = .

The next routine plots the accleration and its tangential and normal components at t = 1.

> with (plots):
curve := spacecurve( [t, t^2, t^3], t = 0..2, color = black, axes=boxed, numpoints = 150, labels = ["x", "y", "z"] ):
accel := spacecurve( [1, 1 + 2*t, 1 + 6*t], t = 0..1, color = green, axes = boxed, numpoints = 150 ):
asubt := spacecurve( [1 + 11*t/7, 1 + 22*t/7, 1 + 33*t/7], t = 0..1, color = blue, axes = boxed):
asubn := spacecurve( [1 - 11*t/7, 1 - 8*t/7, 1 + 9*t/7], t = 0..1, color = blue, axes = boxed ):
parall1 := spacecurve( [18/7 - 11*t/7, 29/7- 8*t/7, 40/7 + 9*t/7], t = 0..1, color = cyan, axes = boxed ):
parall2 := spacecurve( [-4/7 + 11*t/7, -1/7 + 22*t/7, 16/7 + 33*t/7], t = 0..1, color = cyan, axes = boxed ):
xaxis := spacecurve( [t, 0, 0, t = -1..1, color = magenta] ):
yaxis := spacecurve( [0, t, 0, t = -1..1, color = magenta] ):
zaxis := spacecurve( [0, 0, t, t = -1..6, color = magenta] ):
labx := textplot3d([1.2, 0, -.2, `x`], color = magenta):
laby := textplot3d([0,1.2, -.2, `y`], color = magenta):
labz := textplot3d([0, 0, 6.5, `z`], color = magenta):
display(curve, accel, asubt, asubn, parall1, parall2, xaxis, yaxis, zaxis, labx, laby, labz);

>

Calculation of Curvature from the Normal Component of Acceleration. The decomposition of a into its tangential and normal components provides a quick and often convenient way to calculate the curvature K . This avoids ugly differentiating or computation of a cross product. Above the normal component of acceleration is

K (1) N (1) = 14 K (1) = ,

It follows immediately then that

14 K (1) = , and so K (1) = = .

As you can check, that agrees with the result of using the formula from Theorem 3.10 in Section 2.3.