Peeter
- 303
- 3
teodorakis said:if i could undersand the peeter's derivation
Okay, how about this for the angular velocity part (first derivatives). Does this make sense (if it doesn't then it doesn't make much sense to do the next step which is taking second derivatives for acceleration).
Starting point is taking the derivative of:
<br /> \mathbf{r} = r \mathbf{\hat{r}}<br />
<br /> \mathbf{v} = \mathbf{r}' = r' \mathbf{\hat{r}} + r \mathbf{\hat{r}}'<br />
It can be shown (see for example, Salus and Hille, "Calculus") that the unit vector derivative can be expressed using the cross product:
<br /> \mathbf{\hat{r}}' = \frac{1}{r} \left(\mathbf{\hat{r}} \times \frac{d\mathbf{r}}{dt}\right) \times \mathbf{\hat{r}}.<br />
Now, one can express r' in terms of \mathbf{r} as well as follows:
<br /> \left(\mathbf{r} \cdot \mathbf{r}\right)' = 2 \mathbf{v} \cdot \mathbf{r} = 2 r r'.<br />
Thus the derivative of the vector magnitude is part of a projective term:
<br /> r' = \mathbf{\hat{r}} \cdot \mathbf{v}.<br />
Putting this together one has velocity in terms of projective and rejective
components along a radial direction:
<br /> \mathbf{v} = \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right) \mathbf{\hat{r}} + \left(\mathbf{\hat{r}} \times \mathbf{v}\right) \times \mathbf{\hat{r}}.<br />
Now \boldsymbol{\omega} = \frac{\mathbf{r} \times \mathbf{v}}{r^2} term is what we call the angular velocity. The magnitude of this
is the rate of change of the angle between the radial arm and the direction of rotation. The direction of this
cross product is normal to the plane of rotation and encodes both the rotational plane and the direction of the
rotation. Putting these together one has the total velocity expressed radially:
<br /> \mathbf{v} = \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right) \mathbf{\hat{r}} + \boldsymbol{\omega} \times \mathbf{r}.<br />