12 Vector Valued Functions

12.5 The Arc Length Parameter and Curvature

In normal conversation we describe position in terms of both time and distance. For instance, imagine driving to visit a friend. If she calls and asks where you are, you might answer “I am 20 minutes from your house,” or you might say “I am 10 miles from your house.” Both answers provide your friend with a general idea of where you are.

Currently, our vector-valued functions have defined points with a parameter t, which we often take to represent time. Consider Figure 12.5.1(a), where r(t)=t2-t,t2+t is graphed and the points corresponding to t=0, 1 and 2 are shown. Note how the arc length between t=0 and t=1 is smaller than the arc length between t=1 and t=2; if the parameter t is time and r is position, we can say that the particle traveled faster on [1,2] than on [0,1]. margin: -2246246t=0t=1t=2r(t)xy (a)-2246246s=0s=1s=2s=3s=4s=5s=6r(s)xy (b) Figure 12.5.1: Introducing the arc length parameter. Λ

Now consider Figure 12.5.1(b), where the same graph is parameterized by a different variable s. Points corresponding to s=0 through s=6 are plotted. The arc length of the graph between each adjacent pair of points is 1. We can view this parameter s as distance; that is, the arc length of the graph from s=0 to s=3 is 3, the arc length from s=2 to s=6 is 4, etc. If one wants to find the point 2.5 units from an initial location (i.e., s=0), one would compute r(2.5). This parameter s is very useful, and is called the arc length parameter.

How do we find the arc length parameter?

Start with any parameterization of r. We can compute the arc length of the graph of r on the interval [0,t] with

arc length =0tr(u)du.

We can turn this into a function: as t varies, we find the arc length s from 0 to t. This function is

s(t)=0tr(u)du. (12.5.1)

This establishes a relationship between s and t. Knowing this relationship explicitly, we can rewrite r(t) as a function of s: r(s). We demonstrate this in an example.

Example 12.5.1 Finding the arc length parameter

Let r(t)=3t-1,4t+2. Parameterize r with the arc length parameter s.

SolutionUsing Equation (12.5.1), we write

s(t)=0tr(u)du.

We can integrate this, explicitly finding a relationship between s and t:

s(t) =0tr(u)du
=0t32+42du
=0t5du
=5t.

Since s=5t, we can write t=s/5 and replace t in r(t) with s/5:

r(s)=3(s/5)-1,4(s/5)+2=35s-1,45s+2.

Clearly, as shown in Figure 12.5.2, the graph of r is a line, where t=0 corresponds to the point (-1,2). What point on the line is 2 units away from this initial point? We find it with r(2)=1/5,18/5.

margin: 12-1-2246t=0t=1s=1s=2s=3s=4s=5s=0xy Figure 12.5.2: Graphing r in Example 12.5.1 with parameters t and s. Λ

Is the point (1/5,18/5) really 2 units away from (-1,2)? We use the Distance Formula to check:

d=(15-(-1))2+(185-2)2=3625+6425=4=2.

Yes, r(2) is indeed 2 units away, in the direction of travel, from the initial point.

Things worked out very nicely in Example 12.5.1; we were able to establish directly that s=5t. Usually, the arc length parameter is much more difficult to describe in terms of t, a result of integrating a square-root. There are a number of things that we can learn about the arc length parameter from Equation (12.5.1), though, that are incredibly useful.

First, take the derivative of s with respect to t. The Fundamental Theorem of Calculus (see Theorem 5.4.1) states that

dsdt=s(t)=r(t). (12.5.2)

Letting t represent time and r(t) represent position, we see that the rate of change of s with respect to t is speed; that is, the rate of change of “distance traveled” is speed, which should match our intuition.

The Chain Rule states that

drdt =drdsdsdt
r(t) =r(s)r(t).

Solving for r(s), we have

r(s)=r(t)r(t)=T(t), (12.5.3)

where T(t) is the unit tangent vector. Equation (12.5.3) is often misinterpreted, as one is tempted to think it states r(t)=T(t), but there is a big difference between r(s) and r(t). The key to take from it is that r(s) is a unit vector. In fact, the following theorem states that this characterizes the arc length parameter.

Theorem 12.5.1 Arc Length Parameter

Let r(s) be a vector-valued function. The parameter s is the arc length parameter if, and only if, r(s)=1.

Curvature

margin: 12312ABxy (a)12312ABxy (b) Figure 12.5.3: Establishing the concept of curvature. Λ

Consider points A and B on the curve graphed in Figure 12.5.3(a). One can readily argue that the curve curves more sharply at A than at B. It is useful to use a number to describe how sharply the curve bends; that number is the curvature of the curve.

We derive this number in the following way. Consider Figure 12.5.3(b), where unit tangent vectors are graphed around points A and B. Notice how the direction of the unit tangent vector changes quite a bit near A, whereas it does not change as much around B. This leads to an important concept: measuring the rate of change of the unit tangent vector with respect to arc length gives us a measurement of curvature.

Definition 12.5.1 Curvature

Let r(s) be a vector-valued function where s is the arc length parameter. The curvature κ of the graph of r(s) is

κ=dTds=T(s).

If r(s) is parameterized by the arc length parameter, then

T(s)=r(s)r(s)andN(s)=T(s)T(s).

Having defined T(s)=κ, we can rewrite the second equation as

T(s)=κN(s). (12.5.4)

We already knew that T(s) is in the same direction as N(s); that is, we can think of T(s) as being “pulled” in the direction of N(s). How “hard” is it being pulled? By a factor of κ. When the curvature is large, T(s) is being “pulled hard” and the direction of T(s) changes rapidly. When κ is small, T(s) is not being pulled hard and hence its direction is not changing rapidly.

We use Definition 12.5.1 to find the curvature of the line in Example 12.5.1.

Example 12.5.2 Finding the curvature of a line

Use Definition 12.5.1 to find the curvature of r(t)=3t-1,4t+2.

SolutionIn Example 12.5.1, we found that the arc length parameter was defined by s=5t, so r(s)=3s/5-1,4s/5+2 parameterized r with the arc length parameter. To find κ, we need to find T(s).

T(s) =r(s)(recall this is a unit vector)
=3/5,4/5.
Therefore
T(s) =0,0
and
κ =T(s)=0.

It probably comes as no surprise that the curvature of a line is 0. (How “curvy” is a line? It is not curvy at all.)

While the definition of curvature is a useful mathematical concept, it is nearly impossible to use most of the time; writing r in terms of the arc length parameter is generally very hard. Fortunately, there are other methods of calculating this value that are much easier. There is a trade-off: the definition is “easy” to understand though hard to compute, whereas these other formulas are easy to compute though it may be hard to understand why they work.

Theorem 12.5.2 Formulas for Curvature

Let C be a smooth curve on an open interval I in the plane or in space.

  1. (a)

    If C is defined by y=f(x), then

    κ=|f′′(x)|(1+(f(x))2)3/2.
  2. (b)

    If C is defined as a vector-valued function in the plane, r(t)=x(t),y(t), then

    κ=|xy′′-x′′y|((x)2+(y)2)3/2.
  3. (c)

    If C is defined in space by a vector-valued function r(t), then

    κ=T(t)r(t)=r(t)×r′′(t)r(t)3=a(t)N(t)v(t)2.
  • Proof


    We’ll prove statement 3; statements 1 and 2 are applications that we leave to the exercises. By the chain rule and then Equation (12.5.2),

    dT(t)dt=dT(s)dsdsdt=dT(s)dsr(t),

    so that

    κ=dT(s)ds=dT(t)dtr(t)=dT(t)dtr(t)=T(t)r(t).

    Now, r(t)=r(t)T(t) so that

    r′′(t) =(r(t))T(t)+r(t)T(t)  and
    r(t)×r′′(t) =r(t)T(t)×r(t)T(t)=r(t)2T(t)×T(t)

    Because T=1, Theorem 12.2.5 tells us that T and T are orthogonal. This means that

    r(t)×r′′(t) =r(t)2T(t)T(t)=r(t)2T(t),and
    κ =T(t)r(t)=r(t)×r′′(t)r(t)3.

    Theorem 12.4.2 tells us that r′′(t)=aTT(t)+(r′′(t)N(t))N(t). Since r(t)=v(t) and T(t) are parallel, their cross product is zero and

    κ=v(t)×[aTT(t)+(r′′(t)N(t))N(t)]v(t)3=v(t)×(a(t)N(t))N(t)v(t)3.

    Since v(t) and N(t) are orthogonal, the norm of their cross product is the product of their norms, and

    κ=v(t)(a(t)N(t))N(t)v(t)3=a(t)N(t)v(t)2.

We practice using these formulas.

Example 12.5.3 Finding the curvature of a circle

Find the curvature of a circle with radius r, defined by c(t)=rcost,rsint.

SolutionBefore we start, we should expect the curvature of a circle to be constant, and not dependent on t. (Why?)

We compute κ using the second part of Theorem 12.5.2.

κ =|(-rsint)(-rsint)-(-rcost)(rcost)|((-rsint)2+(rcost)2)3/2
=r2(sin2t+cos2t)(r2(sin2t+cos2t))3/2
=r2r3=1r.

We have found that a circle with radius r has curvature κ=1/r.

Example 12.5.3 gives a great result. Before this example, if we were told “The curve has a curvature of 5 at point A,” we would have no idea what this really meant. Is 5 “big” — does is correspond to a really sharp turn, or a not-so-sharp turn? Now we can think of 5 in terms of a circle with radius 1/5. Knowing the units (inches vs. miles, for instance) allows us to determine how sharply the curve is curving.

Let a point P on a smooth curve C be given, and let κ be the curvature of the curve at P. A circle that:

  • passes through P,

  • lies on the concave side of C,

  • has a common tangent line as C at P and margin: 12312ABxy Figure 12.5.4: Illustrating the osculating circles for the curve seen in Figure 12.5.3. Λ

  • has radius r=1/κ (hence has curvature κ)

is the osculating circle, or circle of curvature, to C at P, and r is the radius of curvature. Figure 12.5.4 shows the graph of the curve seen earlier in Figure 12.5.3 and its osculating circles at A and B. A sharp turn corresponds to a circle with a small radius; a gradual turn corresponds to a circle with a large radius. Being able to think of curvature in terms of the radius of a circle is very useful. (The word “osculating” comes from a Latin word related to kissing; an osculating circle “kisses” the graph at a particular point. Many useful ideas in mathematics have come from studying the osculating circles to a curve.)

margin: -10-5510xy Figure 12.5.5: Examining the curvature of y=x2. Λ
Example 12.5.4 Finding curvature

Find the curvature of the parabola defined by y=x2 at the vertex and at x=1.

SolutionWe use the first formula found in Theorem 12.5.2.

κ(x) =|2|(1+(2x)2)3/2
=2(1+4x2)3/2.

At the vertex (x=0), the curvature is κ=2. At x=1, the curvature is κ=2/(5)3/2. So at x=0, the curvature of y=x2 is that of a circle of radius 1/2; at x=1, the curvature is that of a circle with radius (5)3/2/25.59. This is illustrated in Figure 12.5.5. At x=3, the curvature is 0.009; the graph is nearly straight as the curvature is very close to 0.

Example 12.5.5 Finding curvature

Find where the curvature of r(t)=t,t2,2t3 is maximized.

SolutionWe use the third formula in Theorem 12.5.2 as r(t) is defined in space. We leave it to the reader to verify that

r(t) =1,2t,6t2,
r′′(t) =0,2,12t,and
r(t)×r′′(t) =12t2,-12t,2.

Thus

κ(t) =r(t)×r′′(t)r(t)3
=12t2,-12t,21,2t,6t23
=144t4+144t2+4(1+4t2+36t4)3
margin: -2-11212xy (a)
(fullscreen)
(b)
Figure 12.5.6: Understanding the curvature of a curve in space. The top is κ as a function of t; the bottom is r(t) as a function of t. Λ

While this is not a particularly “nice” formula, it does explicitly tell us what the curvature is at a given t value. To maximize κ(t), we should solve κ(t)=0 for t. This is doable, but very time consuming. Instead, consider the graph of κ(t) as given in Figure 12.5.6(a). We see that κ is maximized at two t values; using a numerical solver, we find these values are t±0.189. In part (b) of the figure we graph r(t) and indicate the points where curvature is maximized.

Curvature and Motion

Let r(t) be a position function of an object, with velocity v(t)=r(t) and acceleration a(t)=r′′(t). In Section 12.4 we established that acceleration is in the plane formed by T(t) and N(t), and that we can find scalars aT and aN such that

a(t)=aTT(t)+aNN(t).

Theorem 12.4.2 gives formulas for aT and aN:

aT=ddt(v(t))andaN=v(t)×a(t)v(t).

We understood that the amount of acceleration in the direction of T relates only to how the speed of the object is changing, and that the amount of acceleration in the direction of N relates to how the direction of travel of the object is changing. (That is, if the object travels at constant speed, aT=0; if the object travels in a constant direction, aN=0.)

In Equation (12.5.2) at the beginning of this section, we found s(t)=v(t). We can combine this fact with the above formula for aT to write

aT=ddt(v(t))=ddt(s(t))=s′′(t).

Since s(t) is speed, s′′(t) is the rate at which speed is changing with respect to time. We see once more that the component of acceleration in the direction of travel relates only to speed, not to a change in direction.

Now compare the formula for aN above to the formula for curvature in Theorem 12.5.2:

aN=v(t)×a(t)v(t)andκ=r(t)×r′′(t)r(t)3=v(t)×a(t)v(t)3.

Thus

aN=κv(t)2=κ(s(t))2 (12.5.5)

This last equation shows that the component of acceleration that changes the object’s direction is dependent on two things: the curvature of the path and the speed of the object.

Imagine driving a car in a clockwise circle. You will naturally feel a force pushing you towards the door (more accurately, the door is pushing you as the car is turning and you want to travel in a straight line). If you keep the radius of the circle constant but speed up (i.e., increasing s(t)), the door pushes harder against you (aN has increased). If you keep your speed constant but tighten the turn (i.e., increase κ), once again the door will push harder against you.

Putting our new formulas for aT and aN together, we have

a(t)=s′′(t)T(t)+κv(t)2N(t).

This is not a particularly practical way of finding aT and aN, but it reveals some great concepts about how acceleration interacts with speed and the shape of a curve.

Example 12.5.6 Curvature and road design

The minimum radius of the curve in a highway cloverleaf is determined by the operating speed, as given in the table in Figure 12.5.7. For each curve and speed, compute aN.

Solutionmargin: Operating Speed (mph) Minimum Radius (ft) 35 310 40 430 45 540 Figure 12.5.7: Operating speed and minimum radius in highway cloverleaf design. Λ Using Equation (12.5.5), we can compute the acceleration normal to the curve in each case. We also need to convert each speed from “miles per hour” to “feet per second”.

35mph, 310ft aN=1310ft(35mi1hr×1hr3600s×5280ft1mi)2=8.50ft/s2
40mph, 430ft aN=1430ft(40mi1hr×1hr3600s×5280ft1mi)2=8.00ft/s2
45mph, 540ft aN=1540ft(45mi1hr×1hr3600s×5280ft1mi)2=8.07ft/s2

Note that each acceleration is similar; this is by design. Considering the classic “Force = mass × acceleration” formula, this acceleration must be kept small in order for the tires of a vehicle to keep a “grip” on the road. If one travels on a turn of radius 310ft at a rate of 50mph, the acceleration is double, at 17.35ft/s2. If the acceleration is too high, the frictional force created by the tires may not be enough to keep the car from sliding. Civil engineers routinely compute a “safe” design speed, then subtract 5-10mph to create the posted speed limit for additional safety.

We end this chapter with a reflection on what we’ve covered. We started with vector-valued functions, which may have seemed at the time to be just another way of writing parametric equations. However, we have seen that the vector perspective has given us great insight into the behavior of functions and the study of motion. Vector-valued position functions convey displacement, distance traveled, speed, velocity, acceleration and curvature information, each of which has great importance in science and engineering.

Exercises 12.5

 

Terms and Concepts

  1. 1.

    It is common to describe position in terms of both              and/or             .

  2. 2.

    A measure of the “curviness” of a curve is             .

  3. 3.

    Give two shapes with constant curvature.

  4. 4.

    Describe in your own words what an “osculating circle” is.

  5. 5.

    Complete the identity: T(s)=¯N(s).

  6. 6.

    Given a position function r(t), how are aT and aN affected by the curvature?

Problems

In Exercises 7–10., a position function r(t) is given, where t=0 corresponds to the initial position. Find the arc length parameter s, and rewrite r(t) in terms of s; that is, find r(s).

  1. 7.

    r(t)=2t,t,-2t

  2. 8.

    r(t)=7cost,7sint

  3. 9.

    r(t)=3cost,3sint,2t

  4. 10.

    r(t)=5cost,13sint,12cost

In Exercises 11–22., a curve C is described along with 2 points on C.

  1. (a)

    Using a sketch, determine at which of these points the curvature is greater.

  2. (b)

    Find the curvature κ of C, and evaluate κ at each of the 2 given points.

  1. 11.

    C is defined by y=x3-x; points given at x=0 and x=1/2.

  2. 12.

    C is defined by y=1x2+1; points given at x=0 and x=2.

  3. 13.

    C is defined by y=cosx; points given at x=0 and x=π/2.

  4. 14.

    C is defined by y=1-x2 on (-1,1); points given at x=0 and x=1/2.

  5. 15.

    C is defined by r(t)=cost,sin(2t); points given at t=0 and t=π/4.

  6. 16.

    C is defined by r(t)=cos2t,sintcost; points given at t=0 and t=π/3.

  7. 17.

    C is defined by r(t)=t2-1,t3-t; points given at t=0 and t=5.

  8. 18.

    C is defined by r(t)=tant,sect; points given at t=0 and t=π/6.

  9. 19.

    C is defined by r(t)=4t+2,3t-1,2t+5; points given at t=0 and t=1.

  10. 20.

    C is defined by r(t)=t3-t,t3-4,t2-1; points given at t=0 and t=1.

  11. 21.

    C is defined by r(t)=3cost,3sint,2t; points given at t=0 and t=π/2.

  12. 22.

    C is defined by r(t)=5cost,13sint,12cost; points given at t=0 and t=π/2.

In Exercises 23–26., find the value of x or t where curvature is maximized.

  1. 23.

    y=16x3

  2. 24.

    y=sinx

  3. 25.

    r(t)=t2+2t,3t-t2

  4. 26.

    r(t)=t,4/t,3/t

In Exercises 27–30., find the radius of curvature at the indicated value.

  1. 27.

    y=tanx, at x=π/4

  2. 28.

    y=x2+x-3, at x=π/4

  3. 29.

    r(t)=cost,sin(3t), at t=0

  4. 30.

    r(t)=5cos(3t),t, at t=0

In Exercises 31–34., find the equation of the osculating circle to the curve at the indicated t-value.

  1. 31.

    r(t)=t,t2, at t=0

  2. 32.

    r(t)=3cost,sint, at t=0

  3. 33.

    r(t)=3cost,sint, at t=π/2

  4. 34.

    r(t)=t2-t,t2+t, at t=0

  1. 35.
    For Theorem 12.5.2, use part 3 to prove part 2: If r(t)=x(t),y(t) is a vector-valued function in the plane, then κ=|xy′′-x′′y|((x)2+(y)2)3/2.
  2. 36.
    For Theorem 12.5.2, use part 2 or 3 to prove part 1: If y=f(x), then κ=|f′′(x)|(1+(f(x))2)3/2.
Omni CMS