Conservation of Momentum and Equations of Motion

In the study of solid mechanics, the primary objective is to understand how objects deform and move under the influence of various forces. To predict these behaviors accurately, we rely on a set of fundamental physical principles. Among these, the Law of Conservation of Momentum stands as a cornerstone. It serves as the theoretical bridge that allows us to transition from Newton’s classical laws—originally formulated for discrete particles—to the complex, continuous mathematical models required to describe the motion of solids.

The Concept of Momentum in a Continuum

In classical particle mechanics, momentum is a straightforward product of mass and velocity ($\mathbf{P} = m\mathbf{v}$). However, when we shift our focus to continuum mechanics, we no longer deal with isolated points. Instead, we consider a continuous medium occupying a specific volume $\Omega$.

To describe the momentum of such a region, we must account for the distribution of mass and velocity throughout the entire volume. The total momentum $\mathbf{P}$ of a material volume $\Omega$ is expressed as the volume integral:

$$\mathbf{P} = \int_{\Omega} \rho \mathbf{v} , d\Omega$$

where $\rho$ represents the material density and $\mathbf{v}$ is the velocity vector field. According to the principle of conservation of momentum, the time rate of change of this total momentum must be equal to the sum of all external forces acting on the system:

$$\frac{d\mathbf{P}}{dt} = \mathbf{F}_{ext}$$

Deriving the Local Equations of Motion

While the integral form above provides a global view of the system, engineering analysis requires local (differential) equations that describe the state of the material at any specific point. To derive these, we must categorize the external forces acting on the continuum into two distinct types:

  1. Surface Forces (Traction): These forces act on the boundary $\Gamma$ of the volume. Examples include pressure, friction, or contact loads. The force per unit area on the boundary is referred to as the traction vector $\mathbf{t}$.
  2. Body Forces: These forces act throughout the entire volume $\Omega$, such as gravity, magnetism, or centrifugal forces. The force per unit volume is denoted by $\mathbf{b}$.

1. The Integral Balance of Forces

The total external force $\mathbf{F}_{ext}$ is the sum of the surface and body force contributions:

$$\mathbf{F}{ext} = \int{\Gamma} \mathbf{t} , d\Gamma + \int_{\Omega} \rho \mathbf{b} , d\Omega$$

2. Applying Cauchy’s Stress Principle

To relate the surface traction $\mathbf{t}$ to the internal state of the material, we invoke Cauchy’s Stress Principle. This principle states that the traction vector at any point on a surface is a function of the internal stress tensor $\boldsymbol{\sigma}$ and the unit outward normal vector $\mathbf{n}$ to that surface:

$$\mathbf{t} = \boldsymbol{\sigma} \cdot \mathbf{n}$$

By substituting this into the surface force integral and applying Gauss’s Divergence Theorem, we can convert the surface integral into a volume integral:

$$\int_{\Gamma} \sigma_{ij} n_j , d\Gamma = \int_{\Omega} \frac{\partial \sigma_{ij}}{\partial x_j} , d\Omega$$

3. The Cauchy Equation of Motion

By combining these components into the momentum balance equation, we obtain:

$$\int_{\Omega} \rho \frac{d\mathbf{v}}{dt} , d\Omega = \int_{\Omega} \left( \nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b} \right) , d\Omega$$

Since this equality must hold for any arbitrary volume $\Omega$ within the continuum, the integrands themselves must be equal. This leads us to the fundamental Cauchy Equation of Motion:

$$\rho \left( \frac{\partial v_i}{\partial t} + v_j \frac{\partial v_i}{\partial x_j} \right) = \frac{\partial \sigma_{ij}}{\partial x_j} + \rho b_i$$

In most structural engineering applications involving small deformations and relatively low velocities, the convective term ($v_j \frac{\partial v_i}{\partial x_j}$) is negligible. The equation then simplifies to a more manageable form relating acceleration to stress and body forces:

$$\rho \frac{\partial^2 u_i}{\partial t^2} = \sigma_{ij,j} + \rho b_i$$

where $u_i$ is the displacement vector and $\sigma_{ij,j}$ represents the divergence of the stress tensor.

Physical Interpretation of the Terms

The equation of motion is essentially a dynamic balance of three distinct physical phenomena:

  • Inertial Term ($\rho \ddot{u}_i$): This represents the resistance of the material to acceleration. It is the "mass-effect" that governs how much force is required to change the state of motion.
  • Internal Stress Gradient ($\sigma_{ij,j}$): This term represents the net force resulting from the spatial variation of internal stresses. If the stress is uniform, the gradient is zero; if there is a gradient, it creates a driving force that moves the material.
  • Body Force Term ($\rho b_i$): This accounts for the direct influence of external fields (like gravity) acting on the mass of the material element.

Analysis of Special Cases

Static Equilibrium

When a body is at rest or moving at a constant velocity, the acceleration is zero ($\ddot{u}_i = 0$). The equation of motion reduces to the Equations of Equilibrium:

$$\frac{\partial \sigma_{ij}}{\partial x_j} + \rho b_i = 0$$

This is the fundamental governing equation for all statics and structural stability analyses.

Wave Propagation in Elastic Media

In dynamic scenarios involving elastic materials, the stress is related to strain via a constitutive law (e.g., $\sigma = \mathbf{C} : \varepsilon$). When substituted into the equation of motion, the system describes how disturbances travel through the medium. For a simple one-dimensional elastic rod, this results in the wave equation:

$$\rho \frac{\partial^2 u}{\partial t^2} = E \frac{\partial^2 u}{\partial x^2}$$

This equation allows us to calculate the speed of sound (stress waves) within a solid, defined by $c = \sqrt{E/\rho}$.

Practical Application: Impact on a 1D Rod

Consider a practical example: an elastic rod of length $L$ and cross-section $A$ is fixed at one end and subjected to a sudden impact force $F(t)$ at the other. To model this:

  1. Assumptions: We neglect gravity ($\rho b_i = 0$) and consider only axial motion.
  2. Governing Equation: The motion is governed by $\rho \frac{\partial^2 u}{\partial t^2} = \frac{\partial \sigma_{xx}}{\partial x}$.
  3. Constitutive Relation: Using Hooke's Law, $\sigma_{xx} = E \frac{\partial u}{\partial x}$.
  4. Resulting Model: Combining these yields the wave equation $\rho \frac{\partial^2 u}{\partial t^2} = E \frac{\partial^2 u}{\partial x^2}$.

By solving this partial differential equation, engineers can predict the time-varying stress distribution and the displacement of the rod, which is critical for designing components that must withstand sudden mechanical shocks.

Conclusion

The transition from the global Law of Conservation of Momentum to the local Cauchy Equations of Motion is one of the most vital logical leaps in mechanics. It transforms a general physical principle into a powerful mathematical tool capable of describing the complex interplay between inertia, internal stress, and external loads. Whether analyzing the static load on a bridge or the high-speed impact of a projectile, the underlying logic remains the same: the conservation of momentum dictates the motion of the world.