Consider the following textbook problem:
An electric car of mass $m = 3500 kg$ has stored $E = 100 kWh$ in its battery. The rolling resistance coefficient is $\mu_R = 0.8$. How far can the car travel at most?
The numerical solution is easy enough to obtain, you simply equate the work done by the rolling friction with the initial energy:
$$ F_R \cdot \Delta x = E \Rightarrow \Delta x = \frac{E}{F_R} = \frac{E}{\mu_R m g} = 13.11\ km.$$
However, upon scrutiny this derivation reveals itself to be flawed: the term $F_R \cdot \Delta x$ merely describes the pseudo-work done by the rolling friction, which a priori is not related in any fixed way to the energy of the car; of course, the pseudo-work done by the resulting force describes the change in center-of-mass kinetic energy per unit of time, however, the resulting force would not only include the rolling friction but also the static friction which is needed to even accelerate the car.
Calculation from first principles
If you start from first principles, you would have the following. The car is a system of particles of mass $m_i$, $i = 1,\ldots,N$, which interact according to some potential $U = U(\overrightarrow{x_1},\ldots,\overrightarrow{x_N})$, and which moreover is acted upon some external forces $\overrightarrow{F_i}$, so that the equations of motion would be
$$ m_i \ddot{\overrightarrow{x_i}} = -\frac{\partial U}{\partial \overrightarrow{x_i}} + \overrightarrow{F_i}.$$
Letting the energy of the system be $$ E = \sum_i \frac{1}{2}m_i \dot{\overrightarrow{x_i}}^2 + U,$$
it then follows from the equations of motion that the change of this energy per unit of time is $$ \frac{dE}{dt} = \sum_i \overrightarrow{F_i} \cdot \dot{\overrightarrow{x_i}},$$ i.e. the work done by the external forces alone.
Moreover, in this special situation, we know that there are two types of external forces acting, the static friction accelerating the car and the rolling friction "slowing down" the car; thus we can write
$$ \overrightarrow{F_i} = \overrightarrow{F_i^s} + \overrightarrow{F_i^r}.$$ These forces are unknown, however, we do know that static friction does no work since it acts upon mass points which are at rest at all times. Therefore we find that
$$\frac{dE}{dt} = \sum_i \overrightarrow{F_i^r}\cdot \dot{\overrightarrow{x_i}}.$$
Up until this point, everything is unproblematic.
The problem
However if we want to continue the calculation, we are troubled by the fact that the forces $\overrightarrow{F_i^r}$ are not known; what is known is only the resulting force
$$ \sum_i \overrightarrow{F_i^r} = - \mu_R m g \frac{\overrightarrow{v_S}}{|\overrightarrow{v_S}|},$$ where $v_s = \dot{\overrightarrow{x_S}}$ is the velocity of the center of mass.
And there is an even a more serious problem. If you think about it, it's not reasonable to assume that the rolling friction should do any work at all; in particular, there's no reason to assume that the work done equals the pseudo-work done by the resulting force, even approximately.
Frictional forces do no work
To see this, consider a block, starting with some initial velocity, sliding on a table. Of course, what happens is that the block slows down and stops eventually due to the frictional force. However that does not allow to conclude that the table does negative energy on the block, as the situation is more or less symmetric; if you don't agree, consider two rotating, equal discs pressed against each other by thermally isolated rollers instead.
What should really happen of course is that the frictional force converts ordered kinetic energy into unordered kinetic energy, i.e. heat, inside the same system, having the effect of both slowing down and heating up the block.
So a more realistic assumption is actually to assume that frictional forces do no work at all.
Dissipational energy to the rescue?
So to summarize we are troubled by two things, first that the work done by the rolling friction is unknown since the forces constituting the resulting frictional force are unknown, and second that that work should actually be zero and in particular is unrelated to the distance travelled by the center of mass.
Ignoring the first problem for the moment, we can remedy the second by including another term in the expression for the energy which makes up for any work done by the frictional forces: $$ E' \equiv \sum_i \frac{1}{2}m_i \dot{\overrightarrow{x_i}}^2 + U - \int dt\ \sum_i \overrightarrow{F_i^r} \cdot \dot{\overrightarrow{x_i}}; $$ this of course implies $\frac{dE'}{dt} = 0$.
This new energy term is problematic since it not only depends on positions or velocities at a given time $t$ but the whole history of the system up until $t$; moreover, it's also not clear what it should represent because the heat would seem to be already included in the kinetic energy term.
However, if we press on, we can make the ad-hoc assumption that the frictional forces are given as $$ \overrightarrow{F_i^r} = -\frac{\partial P}{\partial \overrightarrow{v_i}},$$ in terms of a dissipation function $P = P(\overrightarrow{v_1},\ldots,\overrightarrow{v_N}) = P(\dot{\overrightarrow{x_1}}, \ldots, \dot{\overrightarrow{x_N}}).$ Assuming that $P$ is homogeneous of degree $k$, we can compute the work as $$ \sum_i \overrightarrow{F_i^r} \cdot \overrightarrow{v_i} = \sum_i -\frac{\partial P}{\partial \overrightarrow{v_i}} \cdot \overrightarrow{v_i} = -k P.$$
We would then get that
$$ \frac{dE}{dt} = -k P,$$
which for
$$ P = P(|\overrightarrow{v_S}|) = \mu_R m g |\overrightarrow{v_S}|$$
would imply ($k = 1$) first of all that $\frac{dE}{dt} \leq 0$ and more precisely that $$ dE = -\mu_R m g |dx_S|,$$ which then leads to the same result that we get using the naive calculation if we assume some boundedness of the energy $E$ that excludes the dissipational term, i.e. if we assume some bound on $U$.
My Question
It would seem that the introduction of a "dissipational energy" term resolves the problem posed by the fact that--as introspection suggests--frictional forces shouldn't do work, i.e. should not transfer energy into or out of the system, and that the ansatz for the frictional force resolves the problem that the individual forces making up the rolling friction are not known.
However this "solution" seems rather ad hoc and there are at least three problems, some of which have been already indicated above:
- The dissipational energy term seems redundant if heat is interpreted as kinetic energy.
- Dissipational energy is a function of the history of the system.
- There is no reason to assume that the frictional forces are of the form assumed above.
So my question then would be: What's a proper and realistic treatment of this problem in Newtonian mechanics, and how would it differ from the approach above?