Sign in to follow this  
telcom_un

curvature of an arc

Recommended Posts

Hi everybody,

I got a problem and I don't know how to figure it out.

In short I have an arc r(t)=[x(t), y(t), z(t)]
i'm gonna calculate the radius of this arc. As we know the curvature is the inverse of radius. If we consider the following formula:
K(t)=||r'(t) x r''(t)|/||r'(t)||^3
when ' means the first derivative to time and '' means the second derivative to the time and the the X is cross product.

So by having this equation and considering one simple example: A circle at the origin of a 2D space:

x(t)= R (1-t^2)/(1+t^2)
y(t)=R 2t/(1+t^2)
z(t)=0
when tE[0,1] and R is the radius. (EQ1)

By the formula of curvature K(t) is expected to be 1/R but when I compute it strangely for t=0 I get 4*R and for the t=1 I get R and for the time between two margins the K(t) is a non-linear distribution from 4R to R.

I have no idea why it happens. It's very simple example and I am wondering whether I do a silly mistake in my calculations or it's just a wrong assumption in my calculation;.

I created one Maple file to do so. And the result is flawless in sense of numerical and parametric calculation I guess.

In addition, by selecting another form of a circle equation by sinusoidal ones. I get better consistency but it's not even the answer. So if I select

x(t)=R.cos(t)
y(t)=R.sin(t)
z(t)=0

t [0,pi/2]
(EQ2)

The comuted curvature for the whole range of t is R!!! it's wrong because we expect 1/R not R.

So I looked for any possible example. Finally I found one in a literature. They change the parameters and the new circle equation is:

x(t)=R.cos(t/R)
y(t)=R.sin(t/R)
z(t)=0
t [0,pi/2]
(EQ3)

And the result is correct we get the curvature by this equation 1/R!!

So I'm really confused why three different answers I get by the same formula when the resulting geometrical shape is the same. Please help me out thank you in advance.

all the best

Share this post


Link to post
Share on other sites
With the second equation, I get r'(t)=[-Rsin(t), Rcos(t), 0] and r''(t)=[-Rcost(t),-Rsin(t),0], taking the cross product gives [0,0,R^2.sin^2(t)+R^2.cos^2(t)] which is equal to [0,0,R^2] which has a magnitude of R^2. r'(t) has a magnitude of R, so when I divide ||r'(t) X r''(t)|| by the cube of that magnitude I get 1/R. How are you computing the cross product?

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

Sign in to follow this