[Search for users] [Overall Top Noters] [List of all Conferences] [Download this site]

Conference rusure::math

Title:Mathematics at DEC
Moderator:RUSURE::EDP
Created:Mon Feb 03 1986
Last Modified:Fri Jun 06 1997
Last Successful Update:Fri Jun 06 1997
Number of topics:2083
Total number of notes:14613

1228.0. "Looking for formula for a "cubic spiral"." by QB::QBARRY (Jeff Barry, Dtn 297-6756. Let it go.) Wed Apr 25 1990 17:38

Hi,

I'm looking for the formula for something I've seen refered to as a "cubic
sprial".  I don't know if it goes by any other names.  The context in which
I've seen this mentioned is in civil engineering: A cubic sprial is the
transition curve used to minimize "lurch" in a vehicle moving from a
straight path to a constant radius curve.

If anyone could give me the formula or a reference to it, I'd appreciate it.

Jeff
T.RTitleUserPersonal
Name
DateLines
1228.2ALLVAX::JROTHIt's a bush recording...Thu Apr 26 1990 14:4336
   The problem is that a straight line tangent to a circular arc
   has only first order geometric continuity.  That is, the curvature
   undergoes a discontinuous jump, and this is felt as a sudden
   force.

   Any curve which joins the circular arc and the straight line which
   has 3-rd order geometric continuity could smooth out the jump
   in acceleration.  A cubic curve in the shape of a spiral is the
   minimum order polynomial curve which would accomplish this.

   You could also think of it in terms of Fourier analysis - the idea
   would be to have no harmonics below a certain order in the path
   by "low pass filtering" the path.

   Here's a possible way.  Suppose you approach a circle of radius
   1 from below at a speed of w units per second, meeting it at (1,0)
   at t = 0, and thereafter follow the circle with velocity w radians per
   second.

   then
	x(t) = 1,		y(t) = wt		t < 0
	x(t) = cos(wt),		y(t) = sin(wt)		t >= 0

	x'(t) = 0,		y'(t) = w		t < 0
	x'(t) = -w*sin(wt),	y'(t) = w*cos(wt)	t >= 0

	x''(t) = 0,		y''(t) = 0		t < 0
	x''(t) = -w^2*cos(wt),	y''(t) = -w^2*sin(wt)	t >= 0

    You can see a jump in the second x derivative at t = 0; simply
    specify an acceleration that smoothly joins the discontinuity
    over some time interval [-tau,0] so that the path still meets the circle.
    It could even be a simple linear ramp in x acceleration...  upon
    integration this would give a cubic curve.

    - Jim