# Scientific Papers of Josiah Willard Gibbs, Volume 2/Chapter I

I.
ON THE FUNDAMENTAL FORMULÆ OF DYNAMICS.
[American Journal of Mathematics, vol. ii. pp. 49-64, 1879.]
Formation of a new Indeterminate Formula of Motion by the Substitution of the Variations of the Components of Acceleration for the Variations of the Coordinates in the usual Formula.
The laws of motion are frequently expressed by an equation of the form
| $\sum [(X-m{\ddot {x}})\delta x+(Y-m{\ddot {y}})\delta y+(Z-m{\ddot {z}})\delta z]=0,$ | (1) |
| $x,y,z$ | its rectangular coordinates, |
| ${\ddot {x}},{\ddot {y}},{\ddot {z}}$ | the second differential coefficients of the coordinates with respect to the time, |
| $X,Y,Z$ | the components of the forces acting on the particle, |
| $\delta x,\delta y,\delta z$ | any arbitrary variations of the coordinates which are simultaneously possible, and |
| $\sum$ | a summation with respect to all the particles of the system. |
It is evident that we may substitute for $\delta x,\delta y,\delta z$ any other expressions which are capable of the same and only of the same sets of simultaneous values.
Now if the nature of the system is such that certain functions $A,B$, etc. of the coordinates must be constant, or given functions of the time, we have
| $\sum \left({\frac {dA}{dx}}\delta x+{\frac {dA}{dy}}\delta y+{\frac {dA}{dz}}\delta z\right)=0,$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (2) |
| $\sum \left({\frac {dB}{dx}}\delta x+{\frac {dB}{dy}}\delta y+{\frac {dB}{dz}}\delta z\right)=0,$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (2) |
| etc. | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (2) |
These are the equations of condition, to which the variations in the general equation of motion (1) are subject. But if $A$ is constant or a determined function of the time, the same must be true of ${\dot {A}}$ and ${\ddot {A}}$. Now
* ${\dot {A}}=\sum \left({\frac {dA}{dx}}{\dot {x}}+{\frac {dA}{dy}}{\dot {y}}+{\frac {dA}{dz}}{\dot {z}}\right)$
* and ${\ddot {A}}=\sum \left({\frac {dA}{dx}}{\ddot {x}}+{\frac {dA}{dy}}{\ddot {y}}+{\frac {dA}{dz}}{\ddot {z}}\right)+H,$
$H$
$A$
| $\delta {\ddot {A}}=\sum \left({\frac {dA}{dx}}\delta {\ddot {x}}+{\frac {dA}{dy}}\delta {\ddot {y}}+{\frac {dA}{dz}}\delta {\ddot {z}}\right)=0,$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (3) |
| and, in like manner, | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (3) |
| $\delta {\ddot {B}}=\sum \left({\frac {dB}{dx}}\delta {\ddot {x}}+{\frac {dB}{dy}}\delta {\ddot {y}}+{\frac {dB}{dz}}\delta {\ddot {z}}\right)=0.$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (3) |
| etc. | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (3) |
Comparing these equations with (2), we see that when the accelerations of the particles are regarded as subject to the variation denoted by $\delta$, but not their positions or velocities, the possible values of $\delta {\ddot {x}},\delta {\ddot {y}},\delta {\ddot {z}}$ are subject to precisely the same restrictions as the values of math>\delta x, \delta y, \delta z</math>, when the positions of the particles are regarded as variable. We may, therefore, write for the general equation of motion
| $\sum [(X-m{\ddot {x}})\delta {\ddot {x}}+(Y-m{\ddot {y}})\delta {\ddot {y}}+(Z-m{\ddot {z}})\delta {\ddot {z}}]=0,$ | (4) |
$\delta$
| $\delta x=0,$ | $\delta y=0,$ | $\delta z=0,$ | $\scriptstyle {\left.{\begin{matrix}\ \\\ \end{matrix}}\right\}\,$ (5) |
| $\delta {\ddot {x}}=0,$ | $\delta {\ddot {y}}=0,$ | $\delta {\ddot {z}}=0.$ | $\scriptstyle {\left.{\begin{matrix}\ \\\ \end{matrix}}\right\}\,$ (5) |
We have so far supposed that the conditions which restrict the possible motions of the systems may be expressed by equations between the coordinates alone or the coordinates and the time. To extend the formula of motion to cases in which the conditions are expressed by the characters $\leqq$ or $\geqq$, we may write
| $\sum [(X-m{\ddot {x}})\delta {\ddot {x}}+(Y-m{\ddot {y}})\delta {\ddot {y}}+(Z-m{\ddot {z}})\delta {\ddot {z}}]\leqq 0.$ | (6) |
$\delta {\ddot {x}},\delta {\ddot {y}},\delta {\ddot {z}},$
$\delta x,\delta y,\delta z,$
| $\sum [(X-m{\ddot {x}})\delta x+(Y-m{\ddot {y}})\delta y+(Z-m{\ddot {z}})\delta z]\leqq 0,$ | (7) |
This may be illustrated by a simple example.
Let it be required to find the acceleration of a material point, which, at a given instant, is moving with given velocity on the ​frictionless surface of a body (which it cannot penetrate, but which it may leave), and is acted on by given forces. For simplicity, we may suppose that the normal to the surface, drawn outward from the moving point at the moment considered, is parallel to the axis of $X$ and in the positive direction. The only restriction on the values of $\delta x,\delta y,\delta z$ is that
* $\delta x\geqq 0.$
* ${\ddot {x}}\geqq {\frac {X}{m}},{\ddot {y}}\geqq {\frac {Y}{m}},{\ddot {z}}\geqq {\frac {Z}{m}}\cdot$
${\ddot {x}}$
${\ddot {x}}$
$N$
${\ddot {x}}$
* ${\ddot {x}}\geqq N.$
Let us now apply formula (6) to the same problem. Since ${\ddot {x}}$ cannot be less than $N$,
* if ${\ddot {x}}=N,\,\,\delta {\ddot {x}}\geqq 0.$
$\delta {\ddot {x}}$
${\ddot {x}}>N$
$\delta {\ddot {x}}$
| if ${\ddot {x}}=N,$ | ${\ddot {x}}\geqq {\frac {X}{m}};$ |
| but if ${\ddot {x}}>N,$ | ${\ddot {x}}={\frac {X}{m}}:$ |
—that is (since ${\ddot {x}}$ cannot be less than $N$), that ${\ddot {x}}$ shall be equal to the greater of the quantities $N$ and ${\frac {X}{m}}$, or to both, if they are equal,—and that
* ${\ddot {y}}={\frac {Y}{m}},\,\,{\ddot {z}}={\frac {Z}{m}}\cdot$
${\ddot {x}},{\ddot {y}},{\ddot {z}}$
The following considerations will show that what is true in this case is also true in general, when the conditions to which the system ​is subject are such that certain functions of the coordinates cannot exceed certain limits, either constant or variable with the time. If certain values of $\delta {\ddot {x}},\delta {\ddot {y}},\delta {\ddot {z}}$ (with unvaried values of $x,y,z$, and ${\dot {x}},{\dot {y}},{\dot {z}}$) are simultaneously possible at a given instant, equal or proportional values with the same signs must be possible for $\delta x,\delta y,\delta z$ immediately after the instant considered, and must satisfy formula (1), and therefore (6), in connection with the values of ${\ddot {x}},{\ddot {y}},{\ddot {z}},X,Y,Z$ immediately after that instant. The values of ${\ddot {x}},{\ddot {y}},{\ddot {z}}$, thus determined, are of course the very quantities which we wish to obtain, since the acceleration of a point at a given instant does not denote anything different from its acceleration immediately after that instant.
For an example of a somewhat different class of cases, we may suppose that in a system, otherwise free, ${\dot {x}}$ cannot have a negative value. Such a condition does not seem to affect the possible values of $\delta x$, as naturally interpreted in a dynamical problem. Yet, if we should regard the value of $\delta x$ in (7) as arbitrary, we should obtain
* ${\ddot {x}}={\frac {X}{m}},$
$\delta x$
$\delta x\geqq 0$
* ${\ddot {x}}\geqq {\frac {X}{m}}.$
$\delta x$
${\ddot {x}}$
${\dot {x}}=0$
$\delta x\geqq 0$
${\ddot {x}}$
The application of formula (6), in problems of this kind, presents no difficulty. From the condition
* ${\dot {x}}\geqq 0,$
* we obtain, first, if${\dot {x}}=0,$ ${\ddot {x}}\geqq 0,$
* then, if${\dot {x}}=0$and${\ddot {x}}=0,$ $\delta {\ddot {x}}\geqq 0,$
$\delta {\ddot {x}}$
* ${\dot {x}}=0,$${\ddot {x}}=0$and${\ddot {x}}\geqq {\frac {X}{m}};$
* or${\ddot {x}}={\frac {X}{m}}\cdot$
${\dot {x}}=0,{\ddot {x}}$
${\frac {X}{m}}$
${\ddot {x}}={\frac {X}{m}}$
In cases of this kind also, in which the function which cannot exceed a certain value involves the velocities (with or without the coordinates), one may easily convince himself that formula (6) is always valid, and always sufficient to determine the accelerations with the aid of the conditions afforded by the constraints of the system.
But instead of examining such cases in detail, we shall proceed to consider the subject from a more general point of view.
Comparison of the New Formula with the Statical Principle of Virtual Velocities.—Case of Discontinuous Changes of Velocity.
Formula (1) has so far served as a point of departure. The general validity of this, the received form of the indeterminate equation of motion, being assumed, it has been shown that formula (6) will be valid and sufficient, even in cases in which both (1) and (7) fail. We now proceed to show that the statical principle of virtual velocities, when its real signification is carefully considered, leads directly to formula (6), or to an analogous formula for the determination of the discontinuous changes of velocity, when such occur. This will be the case even if we start with the usual analytical expression of the principle
| $\textstyle \sum \displaystyle (X\delta x+Y\delta y+Z\delta z)\geqq 0,$ | (8) |
$\delta x,\delta y,\delta z$
$x,y,z$
$\delta x,\delta y,\delta z$
Now when we suppose a point to have a certain position, specified by $x,y,z$, at a certain time, its position at that time is no longer a subject of hypothesis or of question. It is its future positions which ​form the subject of inquiry. Its position in the immediate future is naturally specified by
* $x+{\dot {x}}dt+{\tfrac {1}{2}}{\ddot {x}}dt^{2}+etc.,$$y+{\dot {y}}dt+{\tfrac {1}{2}}{\ddot {y}}dt^{2}+etc.,$$z+{\dot {z}}dt+{\tfrac {1}{2}}{\ddot {z}}dt^{2}+etc.,$
$\delta x,\delta y,\delta z$
$x,y,z$
${\dot {x}},{\dot {y}},{\dot {z}}$
Let us first consider the most important case, in which there is no discontinuous change of velocity. In this case, ${\dot {x}},{\dot {y}},{\dot {z}}$ are not to be regarded as variable (by $\delta$), and the variations of the above expressions are represented by
* ${\tfrac {1}{2}}\delta {\ddot {x}}dt^{2},$${\tfrac {1}{2}}\delta {\ddot {y}}dt^{2},$${\tfrac {1}{2}}\delta {\ddot {z}}dt^{2},$
$\delta x,\delta y,\delta z$
${\tfrac {1}{2}}dt^{2}$
But if the circumstances are such that there is (or may be) a discontinuity in the values of ${\dot {x}},{\dot {y}},{\dot {z}}$ at the instant considered, it is necessary to distinguish the values of these xpressions before and after the abrupt change. For this purpose, we may apply ${\dot {x}},{\dot {y}},{\dot {z}}$ to the original values, and denote the changed values by ${\dot {x}}+\Delta {\dot {x}},{\dot {y}}+\Delta {\dot {y}},{\dot {z}}+\Delta {\dot {z}}$. The value of $x$ at a time very shortly subsequent to the instant considered, will be expressed by ${\dot {x}}+({\dot {x}}+\Delta x)dt+etc.,$, in which we may regard $\Delta {\dot {x}}$ as subject to the variation denoted by $\delta$. The variation of the expression is therefore $\delta \Delta {\dot {x}}dt$. Instead of $-m{\ddot {x}}$, which expresses the reaction against acceleration, we need in the present case $-m\Delta {\dot {x}}$ to express the reaction against the abrupt change of velocity. A reaction against such a change of velocity is, of course, to be regarded as infinite in intensity in comparison with reactions due to acceleration, and ordinary forces (such as cause acceleration) may be neglected in comparison. If, however, we conceive of the system as acted on by impulsive forces (i.e., such as have no finite duration, but are capable of producing finite changes of velocity, and are measured numerically by the discontinuities of velocity which they produce in the unit of mass), these forces should be combined with the reactions due to the discontinuities of velocity in the general formula which determines these discontinuities. If the impulsive forces are specified by $X,Y,Z$, the formula will be
| $[(X-m\Delta {\dot {x}})\delta \Delta {\dot {x}}+(Y-m\Delta {\dot {y}})\delta \Delta {\dot {y}}+(Z-m\Delta {\dot {z}})\delta \Delta {\dot {z}}]\geqq 0.$ | (9) |
​The reader will remark the strict analogy between this formula and (6), which would perhaps be more clearly exhibited if we should write ${\frac {d{\dot {x}}}{dt}},{\frac {d{\dot {y}}}{dt}},{\frac {d{\dot {z}}}{dt}}$ for ${\ddot {x}},{\ddot {y}},{\ddot {z}}$ that formula.
But these formulæ may be established in a much more direct manner. For the formula (8), although for many purposes the most convenient expression of the principle of virtual velocities, is by no means the most convenient for our present purpose. As the usual name of the principle implies, it holds true of velocities as well as of displacements, and is perhaps more simple and more evident when thus applied.
If we wish to apply the principle, thus understood, to a moving system so as to determine whether certain changes of velocity specified by $\Delta {\dot {x}},\Delta {\dot {y}},\Delta {\dot {z}}$ are those which the system will really receive at a given instant, the velocities to be multiplied into the forces and reactions in the most simple application of the principle are manifestly such as may be imagined to be compounded with the assumed velocities, and are therefore properly specified by $\delta \Delta {\dot {x}},\delta \Delta {\dot {y}},\delta \Delta {\dot {z}}.$ The formula (9) may therefore be regarded as the most direct application of the principle of virtual velocities to discontinuous changes of velocity in a moving system.
In the case of a system in which there are no discontinuous changes of velocity, but which is subject to forces tending to produce accelerations, when we wish to determine whether certain accelerations, specified by ${\ddot {x}},{\ddot {y}},{\ddot {z}},$ are such as the system will really receive, it is evidently necessary to consider whether any possible variation of these accelerations is favored more than it is opposed by the forces ​and reactions of the system. The formula (6) expresses a criterion of this kind in the most simple and direct manner. If we regard a force as a tendency to increase a quantity expressed by ${\ddot {x}}$, the product of the force by $\delta {\ddot {x}}$ is the natural measure of the extent to which this tendency is satisfied by an arbitrary variation of the accelerations. The principle expressed by the formula may not be very accurately designated by the words virtual velocities, but it certainly does not differ from the principle of virtual velocities (in the stricter sense of the term), more than this differs from that of virtual displacements,—a difference so slight that the distinction of the names is rarely insisted upon, and that it is often very difficult to tell which form of the principle is especially intended, even when the principle is enunciated or discussed somewhat at length.
But, although the formulæ (6) and (9) differ so little from the ordinary formulæ, they not only have a marked advantage in respect of precision and accuracy, but also may be more satisfactory to the mind, in that the changes considered (to which $\delta$ relates), are not so violently opposed to all the possibilities of the case as are those which are represented by the variations of the coordinates. Moreover, as we shall see, they naturally lead to various important laws of motion.
Transformation of the New Formula.
Let us now consider some of the transformations of which our general formula (6) is capable. If we separate the terms containing ​the masses of the particles from those which contain the forces, we have
| $\textstyle \sum \displaystyle (X\delta {\ddot {x}}+Y\delta {\ddot {y}}+Z\delta {\ddot {z}})-\textstyle \sum \displaystyle [{\tfrac {1}{2}}m\delta ({\ddot {x}}^{2}+{\ddot {y}}^{2})+{\ddot {z}}^{2}]\leqq 0,$ | (10) |
$u$
| $\textstyle \sum \displaystyle (X\delta {\ddot {x}}+Y\delta {\ddot {y}}+Z\delta {\ddot {z}})-\delta \textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\leqq 0.$ | (11) |
$X\delta x$
$P\delta p$
$p$
$P\delta {\ddot {p}}$
$\delta {\ddot {x}},\delta {\ddot {y}},\delta {\ddot {z}}$
$\delta x,\delta y,\delta z$
$x,y,z$
Hence, if the moments of all the forces of the system are represented by the sum
* ${\mathfrak {S}}(Pdp),$
| ${\mathfrak {S}}(P\delta {\ddot {p}})-\delta \textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\leqq 0$ | (12) |
$V$
* $\delta {\ddot {V}}-\delta \textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\leqq 0,$
| or$\delta [{\ddot {V}}-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})]\leqq 0.$ | (13) |
$X,Y,Z$
$P$
$\delta$
| $\delta \textstyle \sum \displaystyle (X{\ddot {x}}+Y{\ddot {y}}+Z{\ddot {z}}-{\tfrac {1}{2}}mu^{2})\leqq 0,$ | (14) |
| or$\delta [{\mathfrak {S}}(P{\ddot {p}})-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\leqq 0.$ | (15) |
$\delta {\dot {X}}=0,\delta {\dot {Y}}=0,\delta {\dot {Z}}=0$
general formula may be written
| $\delta \left[{\frac {d}{dt}}\textstyle \sum \displaystyle (X{\ddot {x}}+Y{\ddot {y}}+Z{\ddot {z}})-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\right]\quad \leqq 0,$ | (16) |
| or$\delta \left[{\frac {d}{dt}}{\mathfrak {S}}(P{\dot {p}})-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})\right]\quad \leqq 0.$ | (17) |
$\delta$
${\dot {p}},{\dot {x}},{\dot {y}},{\dot {z}}$
$\delta$
${\ddot {p}}',{\ddot {x}}',{\ddot {y}}',{\ddot {z}}'$
* ${\begin{aligned}{\mathfrak {S}}[P({\ddot {p}}+{\ddot {p}}')]&-\textstyle \sum \displaystyle [{\tfrac {1}{2}}m{({\ddot {x}}+{\ddot {x}}')^{2}+({\ddot {y}}+{\ddot {y}}')^{2}+({\ddot {z}}+{\ddot {z}}')^{2}}]\\&-{\mathfrak {S}}(P{\ddot {p}})+\textstyle \sum \displaystyle [{\tfrac {1}{2}}m({\ddot {x}}^{2}+{\ddot {y}}^{2}+{\ddot {z}}^{2})]\\&={\mathfrak {S}}(P{\ddot {p}}')-\textstyle \sum \displaystyle [m({\ddot {x}}{\ddot {x}}'+{\ddot {y}}{\ddot {y}}'+{\ddot {z}}{\ddot {z}}')]-\textstyle \sum \displaystyle [{\tfrac {1}{2}}m({\ddot {x}}'^{2}+{\ddot {y}}'^{2}+{\ddot {z}}'^{2})].\end{aligned}}$
${\ddot {p}}',{\ddot {x}}',{\ddot {y}}',{\ddot {z}}'$
$\delta {\ddot {p}}',\delta {\ddot {x}}',\delta {\ddot {y}}',\delta {\ddot {z}}'$
* ${\mathfrak {S}}(P{\ddot {p}}')-\textstyle \sum \displaystyle [m({\ddot {x}}{\ddot {x}}'+{\ddot {y}}{\ddot {y}}'+{\ddot {z}}{\ddot {z}}')]\leqq 0.$
To show the same with regard to (16) and (17), we have only to observe that the quantities affected by $\delta$ in these formulae differ from those affected by the same symbol in (14) and (15) only by the terms
* $\textstyle \sum \displaystyle ({\dot {X}}{\dot {x}}+{\dot {Y}}{\dot {y}}+{\dot {Z}}{\dot {z}})$and${\mathfrak {S}}({\dot {P}}{\dot {p}}),$
When the forces are determined by the configuration (with or without the time), the principle may be enunciated as follows: The accelerations in the system are always such that the acceleration of the rate of work done by the forces diminished by one-half the sum of the products of the masses of the particles by the squares of their accelerations has the greatest possible value.
The formula (17), although in appearance less simple than (15), not only is more easily enunciated in words, but has the advantage that the quantity ${\frac {d}{dt}}{\mathfrak {S}}(P{\dot {p}})$ is entirely determined by the system with its forces and motions, which is not the case with ${\mathfrak {S}}(P{\ddot {p}})$. The value of the latter expression depends upon the manner in which we choose to represent the forces. For example, if a material point is revolving in a circle under the influence of a central force, we may write either $X{\ddot {x}}+Y{\ddot {y}}+Z{\ddot {z}}$ or $R{\ddot {r}}$ for $P{\ddot {p}}$, $R$ and $r$ denoting respectively the force and radius vector. Now $X{\ddot {x}}+Y{\ddot {y}}+Z{\ddot {z}}$ is manifestly unequal to $R{\ddot {r}}$. But $X{\dot {x}}+Y{\dot {y}}+Z{\dot {z}}$ is equal to $R{\ddot {r}}$, and ${\frac {d}{dt}}(X{\dot {x}}+Y{\dot {y}}+Z{\dot {z}})$ is equal to ${\frac {d}{dt}}(R{\dot {r}})$.
It may not be without interest to see what shape our general formulæ will take in one of the most important cases of forces dependent upon the velocities. If a body which can be treated as a point is moving in a medium which presents a resistance expressed by ​any function of the velocity, the terms due to that resistance in the general formula of motion may be expressed in the form
* $\delta \left[\phi (v){\frac {\dot {x}}{v}}{\ddot {x}}+\phi (v){\frac {\dot {y}}{v}}{\ddot {y}}+\phi (v){\frac {\dot {z}}{v}}{\ddot {z}}\right],$
$v$
$\phi (v)$
* ${\frac {{\dot {x}}{\ddot {x}}}{v}}+{\frac {{\dot {y}}{\ddot {y}}}{v}}+{\frac {{\dot {z}}{\ddot {z}}}{v}}={\frac {dv}{dt}}={\dot {v}}.$
* $\delta [\phi (v){\dot {v}}],$
* or,$\delta {\frac {d}{dt}}f(v),$
$f$
$\phi$
Discontinuous Changes of Velocity.—Formula (9), which relates to discontinuous changes of velocity, is capable of similar transformations.
* If we set$w^{2}=\Delta {\dot {x}}^{2}+\Delta {\dot {y}}^{2}+\Delta {\dot {z}}^{2},$
| $\delta \textstyle \sum \displaystyle (X\Delta {\dot {x}}+Y\Delta {\dot {y}}+Z\Delta {\dot {z}}-{\tfrac {1}{2}}mw^{2})\leqq 0,$ | (10) |
$X,Y,Z$
${\mathfrak {S}}({\text{P}}dp)$
${\text{P}}$
| $\delta [{\mathfrak {S}}({\text{P}}\Delta {\dot {p}})-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mw^{2})]\leqq 0.$ | (19) |
$\delta$
Deduction of other Properties of Motion.
The principles which have been established furnish a convenient point of departure for the demonstration of various properties of motion relating to maxima and minima. We may obtain several such properties by considering how the accelerations of a system, at a given instant, will be modified by changes of the forces or of the constraints to which the system is subject. Let us suppose that the forces $X,Y,Z$ of a system receive the increments $X',Y',Z'$, in consequence of which, and of certain additional constraints, which do not produce any discontinuity in the velocities, the components of acceleration ${\dot {x}},{\dot {y}},{\dot {z}}$ receive the increments ${\ddot {x}},{\ddot {y}},{\ddot {z}}$. The expression
| ${\begin{aligned}\textstyle \sum \displaystyle [(X+X')({\ddot {x}}+{\ddot {x}}')&+(Y+Y')({\ddot {y}}+{\ddot {y}}')+(Z+Z')({\ddot {z}}+{\ddot {z}}')\\&-{\tfrac {1}{2}}m{({\ddot {x}}+{\ddot {x}}')^{2}+({\ddot {y}}+{\ddot {y}}')^{2}+({\ddot {z}}+{\ddot {z}}')^{2}}]\end{aligned}}$ | (20) |
${\ddot {x}},{\ddot {y}},{\ddot {z}}$
| $\textstyle \sum \displaystyle [(X+X'){\ddot {x}}+(Y+Y'){\ddot {y}}+(Z+Z'){\ddot {z}}-{\tfrac {1}{2}}m({\ddot {x}}^{2}+{\ddot {y}}^{2}+{\ddot {z}}^{2})],$ | (21) |
| $\textstyle \sum \displaystyle [X{\ddot {x}}'+Y{\ddot {y}}'+Z{\ddot {z}}'-m({\ddot {x}}{\ddot {x}}'+{\ddot {y}}{\ddot {y}}'+{\ddot {z}}{\ddot {z}}')],$ | (22) |
| and$\textstyle \sum \displaystyle [X'{\ddot {x}}'+Y'{\ddot {y}}'+Z'{\ddot {z}}'-{\tfrac {1}{2}}m({\ddot {x}}^{'2}+{\ddot {y}}^{'2}+{\ddot {z}}^{'2})].$ | (23) |
${\ddot {x}}',{\ddot {y}}',{\ddot {z}}'$
* $\textstyle \sum \displaystyle [X\delta {\ddot {x}}+Y\delta {\ddot {y}}+Z\delta {\ddot {z}}-m({\ddot {x}}\delta {\ddot {x}}+{\ddot {y}}\delta {\ddot {y}}+{\ddot {z}}\delta {\ddot {z}})]=0$
$\delta {\ddot {x}},\delta {\ddot {y}},\delta {\ddot {z}}$
${\ddot {x}}',{\ddot {y}}',{\ddot {z}}'$
${\ddot {x}}+{\ddot {x}}',{\ddot {y}}+{\ddot {y}}',{\ddot {z}}+{\ddot {z}}'$
$-{\ddot {x}}',-{\ddot {y}}',-{\ddot {z}}'$
* $u^{'2}={\ddot {x}}^{'2}+{\ddot {y}}^{'2}+{\ddot {z}}^{'2},$
| $\textstyle \sum \displaystyle (X'{\ddot {x}}'+Y'{\ddot {y}}'+Z'{\ddot {z}}'-{\tfrac {1}{2}}mu^{'2}).$ | (24) |
${\mathfrak {S}}(Qdq)$
$q$
* $\textstyle \sum \displaystyle (X'{\dot {x}}+Y'{\dot {y}}+Z'{\dot {z}})={\mathfrak {S}}(Q{\dot {q}}).$
${\frac {d^{2}q}{dt^{2}}}$
${\ddot {q}}$
${\ddot {q}}+{\ddot {q}}'$
* $\textstyle \sum \displaystyle [{\dot {X}}'{\dot {x}}+{\dot {Y}}'{\dot {z}}+{\dot {Z}}'{\dot {z}}+X'({\ddot {x}}+{\ddot {x}}')+Y'({\ddot {y}}+{\ddot {y}}')+Z'({\ddot {z}}+{\ddot {z}}')]$ $={\mathfrak {S}}[{\dot {Q}}{\dot {q}}+Q({\ddot {q}}+{\ddot {q}}')];$
$\textstyle \sum \displaystyle (X'{\ddot {x}}'+Y'{\ddot {y}}'+Z'{\ddot {z}}')$
${\mathfrak {S}}(Q{\ddot {q}}')$
| ${\mathfrak {S}}(Q{\ddot {q}}')-\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2})$ | (25) |
It will be observed that the condition which determines these relative accelerations is of precisely the same form as that which determines absolute accelerations.
An important case is that in which new constraints are added but no new forcea The relative accelerations are determined in this case by the condition that $\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2})$ is a minimum. In any case of motion, in which finite forces do not act at points, lines or surfaces, we may first calculate the accelerations which would be produced if there were no constraints, and then determine the relative accelerations due to the constraints by the condition that $\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2})$ is a minimum. This is Gauss's principle of least constraint.
Again, in any case of motion, we may suppose $u$ to denote the acceleration which would be produced by the constraints alone, and $u'$ the relative acceleration produced by the forces; we then have
* $\textstyle \sum \displaystyle [m({\ddot {x}}{\ddot {x}}'+{\ddot {y}}{\ddot {y}}'+{\ddot {z}}{\ddot {z}}')]=0,$
$u''$
* $\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2})+\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2})=\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{''2}).$
* ${\mathfrak {S}}(Q\delta {\ddot {q}}')-\textstyle \sum \displaystyle [m({\ddot {x}}'\,\delta {\ddot {x}}'+{\ddot {y}}'\,\delta {\ddot {y}}'+{\ddot {z}}'\,\delta {\ddot {z}}')]=0,$
$\delta {\ddot {q}}',\delta {\ddot {x}}',\delta {\ddot {y}}',\delta {\ddot {z}}'$
${\ddot {q}}',{\ddot {x}}',{\ddot {y}}',{\ddot {z}}',$
* ${\mathfrak {S}}(Q{\ddot {q}}')=2\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2}).$
Particular Equations of Motion.
From the general formula (12), we may easily obtain particular equations which will express the laws of motion in a very general form.
Let $d\omega _{1},d\omega _{2},$ etc. be infinitesimals (not necessarily complete differentials) the values of which are independent, and by means ​of which we can perfectly define any infinitesimal change in the configuration of the system; and let
* ${\dot {\omega }}_{1}={\frac {d\omega _{1}}{dt}},\,\,{\dot {\omega }}_{2}={\frac {d\omega _{2}}{dt}}\,\,{\text{etc.,}}$
$d\omega _{1},d\omega _{2}$
$dt$
* ${\ddot {\omega }}_{1}={\frac {d{\dot {\omega }}_{1}}{dt}},\,\,{\ddot {\omega }}_{2}={\frac {d{\dot {\omega }}_{2}}{dt}}\,\,{\text{etc.}}$
$U=\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{2}).$
$U$
${\dot {\omega }}_{1},{\dot {\omega }}_{2}$
${\ddot {\omega }}_{1},{\ddot {\omega }}_{2}$
$\delta$
* $\delta U={\frac {dU}{d{\ddot {\omega }}_{1}}}\delta {\ddot {\omega }}_{1}+{\frac {dU}{d{\ddot {\omega }}_{2}}}\delta {\ddot {\omega }}_{2}+{\text{etc.}}$
$p$
* ${\ddot {p}}={\frac {dp}{d{\omega }_{1}}}{\dot {\omega }}_{1}+{\frac {dp}{d{\dot {\omega }}_{2}}}{\dot {\omega }}_{2}+{\text{etc.,}}$
${\frac {dp}{d\omega _{1}}}$
$dp$
$d\omega _{1}$
$d\omega _{2}$
${\frac {dp}{d\omega _{2}}}$
$P$
* ${\mathfrak {S}}(P{\dot {p}})=\Omega _{1}{\dot {\omega }}_{1}+\Omega _{2}{\dot {\omega }}_{2}+{\text{etc.,}}$
$\Omega _{1}={\mathfrak {S}}\left(P{\frac {dp}{d\omega _{1}}}\right),$
$\Omega _{2}={\mathfrak {S}}\left(P{\frac {dp}{d\omega _{2}}}\right),$
${\text{etc.}}$
$t$
$\delta$
* ${\mathfrak {S}}(P\delta {\ddot {p}})=\Omega _{1}\delta {\ddot {\omega }}_{1}+\Omega _{2}\delta {\ddot {\omega }}_{2}+{\text{etc.}}$
| $\left(\Omega _{1}-{\frac {dU}{d{\ddot {\omega }}_{1}}}\right)\delta {\ddot {\omega }}_{1}+\left(\Omega _{2}-{\frac {dU}{d{\ddot {\omega }}_{2}}}\right)\delta {\ddot {\omega }}_{2}+{\text{etc.}}\geqq 0.$ | (26) |
$V$
| $\left({\frac {dV}{d\omega _{1}}}-{\frac {dU}{d{\ddot {\omega }}_{1}}}\right)\delta {\ddot {\omega }}_{1}+\left({\frac {dV}{d\omega _{2}}}-{\frac {dU}{d{\ddot {\omega }}_{2}}}\right)\delta {\ddot {\omega }}_{2}+{\text{etc.,}}$ | (27) |
${\frac {dV}{d\omega _{1}}}$
$dV$
$d\omega _{1}$
$d\omega _{2}$
​If the variations $\delta \omega _{1},\delta \omega _{2}$, etc. are capable both of positive and of negative values, we must have
| ${\frac {dU}{d{\ddot {\omega }}_{1}}}=\Omega _{1},$${\frac {dU}{d{\ddot {\omega }}_{2}}}=\Omega _{2},$etc., | (29) |
| or${\frac {dU}{d{\ddot {\omega }}_{1}}}={\frac {dV}{\omega _{1}}},$${\frac {dU}{d{\ddot {\omega }}_{2}}}={\frac {dV}{\omega _{2}}},$etc. | (29) |
$d\omega _{1},d\omega _{2}$
$d\omega _{1},d\omega _{2},d\omega _{3}$
$\Omega _{1},\Omega _{2},\Omega _{3}$
$U$
* ${\begin{aligned}2U&=(a+b+c)({\dot {\omega }}_{1}^{2}+{\dot {\omega }}_{2}^{2}+{\dot {\omega }}_{3}^{2})^{2}-({\dot {\omega }}_{1}^{2}+{\dot {\omega }}_{2}^{2}+{\dot {\omega }}_{3}^{2})(a{\dot {\omega }}_{1}^{2}+b{\dot {\omega }}_{2}^{2}+c{\dot {\omega }}_{3}^{2})\\&+2(b-c){\dot {\omega }}_{2}^{2}{\dot {\omega }}_{3}^{2}{\dot {\omega }}_{1}^{2}+2(c-a){\dot {\omega }}_{3}^{2}{\dot {\omega }}_{1}^{2}{\dot {\omega }}_{2}^{2}+2(a-b){\dot {\omega }}_{1}^{2}{\dot {\omega }}_{2}^{2}{\dot {\omega }}_{3}^{2}\\&+(b+c){\ddot {\omega }}_{1}^{2}+(c+a){\ddot {\omega }}_{2}^{2}+(a+b){\ddot {\omega }}_{3}^{2},\end{aligned}}$
$a,b$
$c$
$a+b,b+c,c+a$
| ${\frac {dU}{d{\ddot {\omega }}_{1}}}=(b-c){\dot {\omega }}_{2}{\dot {\omega }}_{3}+(b+c){\ddot {\omega }}_{1},$ | ${\frac {dU}{d{\ddot {\omega }}_{2}}}=(c-a){\dot {\omega }}_{3}{\dot {\omega }}_{1}+(c+a){\ddot {\omega }}_{2},$ |
| ${\frac {dU}{d{\ddot {\omega }}_{3}}}=(a-b){\dot {\omega }}_{1}{\dot {\omega }}_{2}+(a+b){\ddot {\omega }}_{3};$ ||
and the equations of motion are
* ${\ddot {\omega }}_{1}={\frac {(c-b){\dot {\omega }}_{2}{\dot {\omega }}_{3}+\Omega _{1}}{c+b}},$ ${\ddot {\omega }}_{2}={\frac {(a-c){\dot {\omega }}_{3}{\dot {\omega }}_{1}+\Omega _{2}}{a+c}},$ ${\ddot {\omega }}_{3}={\frac {(b-a){\dot {\omega }}_{1}{\dot {\omega }}_{2}+\Omega _{3}}{b+a}}\cdot$
* ↑ The failure of the formula (7) in this ease is rather apparent than real; for, although the formula apparently allows to ${\ddot {x}}$, at the instant considered, a value exceeding both $N$ and ${\frac {X}{m}}$, it does not allow this for any interval, however short. For if ${\ddot {x}}<N$ the point will immediately leave the surface, and then the formula requires that${\ddot {x}}={\frac {X}{m}}\cdot$
* ↑ Even in Statics, the principle of virtual velocities, as distinguished from that of virtual displacements, has a certain advantage in respect of its evidence. The demonstration of the principle in the first section of the Mécanique Analytique, if velocities had been considered instead of displacements, would not have been exposed to an objection, which has been expressed by M. Bertrand in the following words: "On a objecté, avec raison, à cette assertion de Lagrange l'example d'un point pesant en équilibre au sommet le plus élevé d'une courbe; il est évident qu'un déplacement infiniment petit le ferait descendre, et, pourtant, ce déplacement ne se produit pas." (Mécanique Analytique, troiséme édition, tome 1, page 22, note de M. Bertrand.) The value of $z$ (the height of the point above a horizontal plane) can certainly be diminished by a displacement of the point, but the value of i is not affected by any velocity given to the point. The real difficulty in the consideration of displacements is that they are only possible at a time subsequent to that in which the system has the configuration to which the question of equilibrium relates. We may make the interval of time infinitely short, but it will always be difficult, in the establishing of fundamental principles, to treat a conception of this kind (relating to what is possible after an infinitesimal interval of time) with the same rigor as the idea of velocities or accelerations, which, in the cases to which (9) and (6) respectively relate, we may regard as communicated immediately to the system.
* ↑ It may have seemed to some readers of the Mécanique Analytique—a work of which the unity of method is one of the most striking oharaoteristios, and that to whioh its universally recognized artistic merit is in great measure due—that the treatment of dynamical problems in that work is not entirely analogous to the treatment of statical problems. The statical question, whether a system will remain in equilibrium in a given configuration, is determined by Lagrange by considering all possible motions of the system and inquiring whether there is any reason why the system should take any one of them. A similar method in dynamics would be based upon a comparison of a proposed motion with all other motions of which the system is capable without violating its kinematical conditions. Instead of this, Lagrange virtually reduces the dynamical problem to a statical one, and considers, not the possible variations of the proposed motion, but the motions which would be possible if the system were at rest. This reduction of a given problem to a simpler one, which has already been solved, is a method which has its advantages, but it is not the characteristic method of the Mécanique Anaiytique. That which most distinguishes the plan of this treatise from the usual type is the direct application of the general principle to each particular case. The point is perhaps of small moment, and may be differently regarded by others, but it is mentioned here because it was a feeling of this kind (whether justified or not) and the desire to express the formula of motion by means of a maximum or minimum condition, in which the conditions under which the maximum or minimum subsists should be such as the problem naturally affords (Gauss's principle of least constraint being at the time unknown to the present writer, and the conditions under which the minimum subsists in the principle of least action being such that that is hardly satisfactory as a fundamental principle), which led to the formula proposed in this paper.
* ↑ As an illustration of the significance of this limitation, we may consider the condition afforded by the impenetrability of two bodies in contact. Let us suppose that if subject only to the original forces and constraints they would continue in contact, but that, under the influence of the additional forces and constraints, the contact will cease. The impenetrability of the bodies then ceases to be operative as a constraint. Such cases form an exception to the principle which is to be established. But there are no exceptions when all the original constraints are expressed by equations.
* ↑ This principle may be derived very directly from the general formula (6), or vice versa, for $\textstyle \sum \displaystyle ({\tfrac {1}{2}}mu^{'2})$ may be put in the form $\textstyle \sum \displaystyle \left[{\tfrac {1}{2}}m\left\{\left({\ddot {x}}-{\frac {X}{m}}\right)^{2}+\left({\ddot {y}}-{\frac {Y}{m}}\right)^{2}+\left({\ddot {z}}-{\frac {Z}{m}}\right)^{2}\right\}\quad \right],$ the variation of which, with the sign changed, is identical with the first member of (6).