@jbriggs444 That's the sort of thing I had in mind as a way of formalising a description of the rotational motion of a rigid body in terms of all the points - effectively regarding it as an uncountable collection of point particles, all maintaining fixed distances from one another.
jbriggs444 said:
If I am not clumsily mistaken, it should be immediately clear that P(t,i,j) = (rt+icos(t)+jsin(t),r+jsin(t)-icos(t)) is such a function.
I don't think you're mistaken, but it's not
immediately clear to me. However it looks promising, so I'll try to perform the service of working through the calcs to show this is the case.
But, being inherently lazy, I'm going to change to local polar coordinates because I think that will be easier. Let the time-##t## location of a point on the wheel have polar coordinates ##r(t),\theta(t)## in a coordinate system whose origin is the axle and whose x-axis is parallel to the floor. The frame for those polar coordinates moves as the axle does.
We also use a global, stationary frame of reference, with Cartesian coordinates, to describe the absolute location of a point on the wheel. Let the origin be the initial axle position and the x-axis be parallel to the floor.
Each point particle ##P## in the wheel has a locus that is a function ##f_P## from ##\mathbb R## to ##\mathbb R^2## giving the points global coordinates at each time. Let's define functions ##r_P,\theta_P,x_p,y_P## that give, respectively, the local polar radial coordinate, local polar angular coordinate, global Cartesian ##x## coordinate and global Cartesian ##y## coordinate of the particle in terms of time.
To keep things simple, consider a wheel that is translating horizontally at a constant speed ##v##, and at the same time rotating around its axle with angular velocity ##\omega##. It's straightforward to extend this to non-constant motion, but we'd need to introduce integrals, which would be a bit messy.
Applying the definitions of rotation and translation, we see that the locus functions are:
\begin{align*}
r_P(t)&=r_P(0)\\
\theta_P(t)&=\theta_P(0) + \omega t\\
x_P(t) &= vt + r_P(t) \cos \theta_P(t)= vt + r_P(0) \cos (\theta_P(0)+\omega t)\\
y_P(t) &= r_P(t) \sin \theta_P(t)= r_P(0) \sin (\theta_P(0)+\omega t)
\end{align*}
Then the time-##t## distance between two particles ##P## and ##Q## will be
$$\sqrt{(x_P(t)-x_Q(t))^2 + (y_P(t)-y_Q(t))^2}$$
This will remain constant iff its square remains constant, which means its square has a zero time derivative. The time derivative of the square is:
\begin{align*}
\frac d{dt} [ (x_P(t)&-x_Q(t))^2 + (y_P(t)-y_Q(t))^2 ]
=
\frac d{dt} [ (vt + r_P(0) \cos (\theta_P(0)+\omega t)) - (vt + r_Q(0) \cos (\theta_Q(0)+\omega t)))^2 \\
&\quad\quad\quad\quad+ ((r_P(0) \sin (\theta_P(0)+\omega t)) - (r_Q(0) \sin (\theta_Q(0)+\omega t)))^2 ]\\
&=
\frac d{dt} [ (r_P(0) \cos (\theta_P(0)+\omega t) - r_Q(0) \cos (\theta_Q(0)+\omega t)))^2 \\
&\quad\quad\quad\quad+ ((r_P(0) \sin (\theta_P(0)+\omega t)) - (r_Q(0) \sin (\theta_Q(0)+\omega t)))^2 ] \\
&=
\frac d{dt} [ (r_P(0)^2 +r_Q(0)^2
-2 r_P(0) r_Q(0) \cos (\theta_P(0)+\omega t) \cos (\theta_Q(0)+\omega t)\\
&\quad\quad\quad\quad-2 r_P(0) r_Q(0) \sin (\theta_P(0)+\omega t) \sin (\theta_Q(0)+\omega t)\\
&=
-2 r_P(0) r_Q(0)
\frac d{dt} (
(\cos \theta_P(0) \cos \omega t - \sin \theta_P(0) \sin \omega t)
(\cos \theta_Q(0) \cos \omega t - \sin \theta_Q(0) \sin \omega t)\\
&\quad\quad\quad\quad+
(\sin \theta_P(0) \cos \omega t - \cos \theta_P(0) \sin \omega t)
(\sin \theta_Q(0) \cos \omega t - \cos \theta_Q(0) \sin \omega t)
)\\
&=
-2 r_P(0) r_Q(0)
\frac d{dt} (
\sin \theta_P(0) \sin \theta_Q(0) + \cos \theta_P(0) \cos \theta_Q(0))\\
&=0
\end{align*}
So the shape does indeed remain rigid.
To forbid slipping, we set ##v=\omega R## where ##R## is the wheel radius. That makes the instantaneous speed of the lowest point of the wheel in the global Cartesian frame always zero.
On reflection, what this proves is not something about physics, but rather the mathematical fact that rotation and translation are isometries - ie transformations that, when applied to a set of points, preserve the shape and size of that set.
To incorporate physics into the analysis, we need to introduce one or more forces and - using the rigidity as a set of constraint forces - apply d'Alembert's principle (I think). That's beyond the scope of this post, which is already too long.