D'Alembert's Principle and Generalized Forces
Statement
For a system of N particles subject to ideal (workless) constraints, the total virtual work of the applied forces together with the inertial "reversed effective forces" vanishes for every virtual displacement consistent with the instantaneous constraints: Σi (Fi(a) − ṗi) · δri = 0. Projected onto generalized coordinates qj, this is equivalent to d/dt(∂T/∂q̇j) − ∂T/∂qj = Qj, where Qj = Σi Fi(a) · ∂ri/∂qj is the generalized force.
Why it matters
D'Alembert's principle is the bridge from vectorial Newtonian mechanics to the analytical mechanics of Lagrange. Its decisive move is to eliminate the unknown constraint forces from the equations of motion entirely: because ideal constraints do no virtual work, they simply drop out of the balance, leaving only quantities we can actually compute — applied forces and kinetic energy.
It also reframes dynamics as a statement of virtual work, i.e. a variational-flavoured condition, rather than a set of vector equations. This is the conceptual seed of the whole Lagrangian and Hamiltonian program, and it lets us choose coordinates adapted to the constraints (angles, arc-lengths, mode amplitudes) instead of Cartesian components tied together by awkward reaction forces.
Assumptions
Derivation
Result
Reading. The applied forces and the reversed inertial forces are in balance against every allowed virtual displacement; projected onto each independent coordinate this becomes an equation of motion in which the generalized inertia (the T-terms) equals the generalized force. Constraint reactions never appear. When the applied forces are derivable from a potential, Qj=−∂V/∂qj and the pair collapses to the Euler–Lagrange equations with L=T−V.
Units check. Every term of Σ(F−ṗ)·δr is force×length = N·m = J (virtual work). In the coordinate form, ∂T/∂q̇j is a generalized momentum with units J·(unit of q̇j)−1=J·s·(unit of qj)−1; its time derivative and ∂T/∂qj both carry units J·(unit of qj)−1, matching Qj, so that Qjδqj is an energy. For a Cartesian q this is a force (N); for an angular q it is a torque (N·m).
Limiting cases
- Statics (all q̇j=q̈j=0): T≡0, so the equations reduce to Qj=0 — the principle of virtual work for equilibrium. D'Alembert is virtual work with inertia added.
- No constraints, Cartesian coordinates (qj→x,y,z): ∂ri/∂qj are unit basis vectors, Qj=Fj, T=½m|v|2, and each equation is exactly mr̈=F. Newton is recovered as a special case.
- Conservative applied forces: Qj=−∂V/∂qj with V=V(q,t); the equations become d/dt(∂L/∂q̇j)−∂L/∂qj=0, L=T−V.
- Single free particle, one coordinate: reduces to the familiar 1-D Newton equation, confirming the projection introduces nothing new when there is nothing to constrain.
Breaks when
- Constraint forces do work (non-ideal constraints). Sliding (kinetic) friction, rolling with slip, or magnetic-brake reactions give Σifi·δri≠0. Step 4 fails; the reactions no longer cancel and must be modelled explicitly as applied forces or via a Rayleigh dissipation function.
- Non-holonomic, non-integrable constraints. When constraints restrict velocities but not configurations (e.g. a vertical coin rolling without slipping), the δqj in step 13 are not independent, so the brackets cannot be set to zero individually; one must retain Lagrange multipliers λa for each constraint one-form.
- Rheonomic constraints treated with real displacements. If one carelessly uses dri (with δt≠0) instead of the instantaneous virtual displacement, moving constraints inject spurious work Σfi·(∂ri/∂t)dt and the cancellation is lost.
- Relativistic or field-momentum systems. With p=γmv or momentum stored in fields, the identity ṗ·δr no longer reduces to derivatives of the Newtonian T=½Σmṙ2; a covariant action principle is required instead.
Failure modes
- Confusing δr with dr. Students use the actual trajectory displacement (which includes the constraint's own motion in time) rather than the frozen-time virtual displacement, so constraint work fails to cancel.
- Including constraint forces in Qj. Adding tension, normal reactions, or rod forces into the generalized force double-counts what D'Alembert was designed to remove; only genuinely applied forces (gravity, springs, drives) enter Qj.
- Forgetting the −∂T/∂qj term. In curvilinear coordinates T depends on q (e.g. T=½m(ṙ2+r2θ̇2)); dropping ∂T/∂q loses the centrifugal/Coriolis pseudo-forces.
- Setting each bracket to zero for non-holonomic systems. The independence of δqj is only guaranteed for holonomic constraints; applying step 13 to a rolling constraint gives wrong equations.
- Miscounting degrees of freedom. Using more coordinates than independent DOF (without multipliers) makes the δqj dependent and the resulting equations inconsistent.
- Differentiating T before substituting the constraint. One must express T in the independent q,q̇ first; substituting constraints into T prematurely can silently drop terms.
Discussion
D'Alembert's genius was to convert a dynamics problem into a statics problem. By reversing the inertial force −ṗi and treating it on the same footing as an applied force, every instant of motion is recast as an equilibrium under the combined system of applied, constraint, and inertial forces. The principle of virtual work — already trusted for statics — then applies, and the ideal-constraint hypothesis does the real work: it annihilates precisely the forces we can never easily compute. What survives is a scalar (energy) statement, not a set of coupled vector equations, and scalars transform trivially under changes of coordinate. This is why we may adopt angles, arc-lengths, or normal-mode amplitudes freely.
The connection to the rest of analytical mechanics is direct. For conservative forces the principle becomes the Euler–Lagrange equations and, by Legendre transform, Hamilton's equations. It is closely related to, but logically distinct from, Hamilton's principle δ∫L dt=0: Hamilton's principle is an integral (global, over a path) variational statement, whereas D'Alembert's is a differential (local, at each instant) one. One can derive Hamilton's principle from D'Alembert's by integrating over time and discarding endpoint terms, which is why some texts call D'Alembert's principle the "central equation" of mechanics.
A subtlety worth stressing is the geometric status of the virtual displacement. The set {δri} lives in the tangent space to the constraint surface at fixed t; the vectors ∂ri/∂qj are a basis for that tangent space. D'Alembert's principle asserts that the vector (F(a)−ṗ) is orthogonal to this tangent space, i.e. it lies purely in the space spanned by the constraint reactions (the normal bundle). For non-holonomic systems the admissible virtual displacements form a distribution that is not the tangent space of any surface, and orthogonality to it no longer permits component-wise vanishing — hence the multipliers. This is the precise geometric reason holonomic and non-holonomic constraints behave so differently, a distinction obscured if one works only with components.
Common misconceptions. D'Alembert's principle is not merely "F=ma with the terms moved to one side" — the physics is the vanishing of constraint virtual work, which is an independent postulate about the nature of ideal constraints. Nor is the "inertial force" −ṗ a real force in any frame; it is a bookkeeping device (an effective force) that becomes a genuine fictitious force only after one commits to a non-inertial frame. Finally, the principle is not restricted to statics-like slow motion: it holds instant by instant for arbitrary accelerated motion.
Worked examples
Reading. The tension never entered the calculation — D'Alembert removed it automatically. Units check. (kg)(m/s2)/(kg)=m/s2. ✓
Reading. The rod tension (a large radial force) is invisible here because it does no virtual work; only the tangential gravity torque drives θ. Units check. (m/s2)/m × dimensionless = s−2=rad/s2. ✓
Problems
- A bead of mass m slides on a frictionless horizontal wire; a horizontal applied force F(t) acts along the wire. Using x as generalized coordinate, write Qx and the equation of motion, and identify why the normal reaction is absent.
Solution
The wire exerts a normal (vertical) reaction f⊥δr (which is horizontal), so f·δr=0 — ideal constraint, dropped by D'Alembert. T=½m ẋ2, so d/dt(∂T/∂ẋ)=m ẍ, ∂T/∂x=0. The only applied force along the wire is F(t): Qx=F(t). Equation of motion: m ẍ = F(t). The normal reaction never appears because it is orthogonal to every admissible virtual displacement. - For the plane pendulum of Worked Example 2, take the small-angle limit and find the period for l=0.50 m. State the approximation used and give a numeric answer.
Solution
For small θ, sinθ≈θ, so θ̈=−(g/l)θ, simple harmonic with ω=√(g/l). Period T=2π√(l/g)=2π√(0.50/9.81)=2π√(0.05097)=2π(0.2258)=1.42 s. Approximation: sinθ≈θ (error <1% for θ<14°); the amplitude-dependence of the true period is thereby neglected. - A particle moves in a plane described by polar coordinates (r,θ) with no constraints, under a central applied force F=Frêr. Using T=½m(ṙ2+r2θ̇2), derive the two equations of motion and show angular momentum is conserved.
Solution
Generalized forces: with r=rêr, ∂r/∂r=êr and ∂r/∂θ=rêθ, so Qr=Fr, Qθ=0 (central force has no θ-component).
r-equation: ∂T/∂ṙ=m ṙ, d/dt=m r̈; ∂T/∂r=m rθ̇2. Thus m r̈ − m rθ̇2 = Fr (the −m rθ̇2 is the centrifugal term).
θ-equation: ∂T/∂θ̇=m r2θ̇, d/dt(m r2θ̇)=0; ∂T/∂θ=0. So d/dt(m r2θ̇)=0 ⇒ L=m r2θ̇ is conserved. The vanishing of Qθ is exactly the statement that a central force exerts no torque about the origin. - A block of mass m rests on a frictionless wedge of mass M that is free to slide on a frictionless floor; the incline angle is α. Using generalized coordinates X (wedge position) and s (block's displacement along the incline relative to the wedge), set up QX, Qs and the kinetic energy, then write the two equations of motion. (Do not solve fully; reduce to two coupled equations.)
Solution
Let the incline descend to the right at angle α. Block position: xb=X+s cosα, yb=−s sinα. Velocities: ẋb=Ẋ+ṡcosα, ẏb=−ṡsinα.
T=½M Ẋ2+½m[(Ẋ+ṡcosα)2+(ṡsinα)2]=½M Ẋ2+½m[Ẋ2+2Ẋṡcosα+ṡ2].
Applied force is gravity on the block, (0,−mg). QX=(0,−mg)·∂rb/∂X=(0,−mg)·(1,0)=0. Qs=(0,−mg)·∂rb/∂s=(0,−mg)·(cosα,−sinα)=mg sinα.
X-equation: d/dt[(M+m)Ẋ+m ṡcosα]=0 ⇒ (M+m)Ẍ+m s̈cosα=0 (horizontal momentum conserved, since QX=0).
s-equation: ∂T/∂ṡ=m(Ẋcosα+ṡ), ∂T/∂s=0, so m(Ẍcosα+s̈)=mg sinα. These are the two coupled equations; note the constraint (normal) force between block and wedge has dropped out. - State precisely, with one worked reason each, two physical situations in which the step "each bracket vanishes because the δqj are independent" (Derivation step 13) is invalid, and say what replaces it.
Solution
(i) Non-holonomic constraint — e.g. a vertical disc rolling without slipping on a plane. The rolling condition relates velocities (ẋ=Rφ̇cosθ, etc.) but is not integrable to a relation among coordinates, so the δqj satisfy differential relations Σjaajδqj=0 and are not independent. Replacement: introduce Lagrange multipliers λa, giving d/dt(∂T/∂q̇j)−∂T/∂qj=Qj+Σaλaaaj. (ii) Redundant coordinates — if one uses n coordinates for a system with fewer independent DOF (a holonomic constraint g(q,t)=0 still imposed), then Σj(∂g/∂qj)δqj=0 makes the δqj dependent. Replacement: again a multiplier λ times ∂g/∂qj, or eliminate one coordinate using g=0 to restore independence. In both cases the naive "set each bracket to zero" gives wrong equations of motion.