# Two nonholonomic particles in mutual potential

1. ### coolnessitself

32
Hi all,
1. The problem statement, all variables and given/known data
Not actually homework, but hopefully it fits well in this forum.
I have two particles that can only move in the direction of their respective heading angles. I'm trying to relate a radially-symmetric mutual potential between the particles to the acceleration of this heading angle. For example,
$$\begin{bmatrix} \dot{x}_i \\ \dot{y}_i \\ \dot{z}_i \\ \dot{\theta}_i \\ \dot{\Theta}_i \\ \dot{\phi}_i \\ \dot{\Phi}_i \end{bmatrix} = \begin{bmatrix} v\sin\theta_i\cos\phi_i \\ v\sin\theta_i\sin\phi_i \\ v\cos\theta_i \\ \Theta_i \\ f(r) \\ \Phi_i \\ g(r)\end{bmatrix}$$
for the distance r between 1 and 2. Normally we might say that $$\ddot{r} = -\nabla G(r)$$, but here in order for r to change, the heading angles must change. So my question is how to find f(r) and g(r) (the acceleration of the heading angles) so that the same behavior occurs as with $$\ddot{r} = -\nabla G(r)$$. If we require that the direction of r is constant throughout this acceleration, ie $$\dot{r} = (\ldots) \hat{r} + 0\hat{\theta} + 0\hat{\phi}$$ then it seems to me like there should be a unique solution (if one exists).

2. Relevant equations

3. The attempt at a solution
If the direction of r doesn't change,
$$(\dot{z}_2 - \dot{z}_1)r = (z_2 - z_1)\dot{r}$$ and $$(\dot{y}_2 - \dot{y}_1)(x_2-x_1) = (y_2-y_1)(\dot{x}_2 - \dot{x}_1)$$. My hope would be to then expand out $$\ddot{r}$$ until equations for $$\ddot{\theta}$$ and $$\ddot{\phi}$$ appear, then use that condition to find one solution. Something like
$$(\dot{r})^2 + r\ddot{r} = (\dot{x}_2 - \dot{x}_1)^2 + (x_2-x_1)(\ddot{x}_2 - \ddot{x}_1) + (\dot{y}_2 - \dot{y}_1)^2 + (y_2-y_1)(\ddot{y}_2 - \ddot{y}_1) + (\dot{z}_2 - \dot{z}_1)^2 + (z_2-z_1)(\ddot{z}_2 - \ddot{z}_1)$$
But after plugging in $$\dot{x}$$, etc. with phis and thetas, It doesn't seem to lead anywhere.
So I get the feeling I'm going about this incorrectly. Physics gurus, if you have any suggestions for an alternative approaches to finding this heading-angle-potential, I'd appreciate your input!

2. Science news on Phys.org