I have done this calculation several times over the years, so I might as well write it down in detail, in case it may be of use to someone else.
We are interested in the curvature C=1/R of a planar curve y=f(x) at a given point A, where R is the curvature radius at that particular point, defined with respect to the curvature center O (intersection of the normals raised to the curve in A and its infinitesimal neighbor B.)
The angle subtending AB is: dα=ds/R⇒C=dαds
The length of the curve element AB is: ds=√dx2+dy2⇒dsdx=√1+f′(x)2
The derivative of f is directly related to the angle α: f′(x)=dydx=tanα⇒α=arctandydx=arctan[f′(x)]⇒dαdx=11+f′(x)2f″(x)
Putting together the three relations above yields:
C=dαds=f″(x)[1+f′(x)2]3/2
No comments:
Post a Comment