Cauchy's Equation of Motion

Newton's second law, $F=ma$, masterfully describes the motion of whole objects, but what governs the motion within those objects? How does one part of a bending beam or a flowing river exert force on its immediate neighbor? This question marks the critical transition from the mechanics of discrete particles to the rich and powerful world of continuum mechanics. The challenge lies in creating a law of motion that applies not just to a body, but to every single point inside it.

This article delves into the fundamental principle that answers this question: Cauchy's equation of motion. It provides a comprehensive framework for understanding how forces and motion are distributed within any continuous medium—be it solid, liquid, or gas.

We will embark on this exploration in two main parts. First, under "Principles and Mechanisms," we will deconstruct the equation itself, translating Newton's law into a local statement by introducing the crucial concepts of body forces and the stress tensor. You will learn how inertia, internal interactions, and external forces are elegantly balanced at every infinitesimal point. Following this, the "Applications and Interdisciplinary Connections" section will demonstrate the equation's immense power by surveying its diverse uses, revealing how this single law governs phenomena from the stability of bridges and the propagation of seismic waves to the dynamics of fluids and the expansion of the cosmos.

You already know Newton's second law, the famous $F=ma$. It's the bedrock of mechanics. You picture a force hitting a billiard ball, and the ball accelerates. The law works perfectly for whole objects—what we might call "chunks" of matter. But what about the matter within the chunk? When you bend a ruler, one part of the ruler pulls on the next, which pulls on the next. When a river flows, a faster-moving parcel of water drags its slower neighbor along. How does one tiny piece of a continuum—be it solid, liquid, or gas—communicate with the piece right next to it? How does the law of motion apply not just to the whole body, but to every single infinitesimal point within it?

This is the magnificent leap we are about to make. We will transform Newton's law "for a body" into a law that holds at every point in space. It is a journey from a statement about a whole volume to a precise, local equation—a theme that runs through all of modern physics.

Let’s zoom in. Way in. Imagine a tiny, imaginary cube of material buried deep inside a larger object. Perhaps it’s a cube of steel in a skyscraper's beam, or a cube of water in the ocean. What forces act on this cube? They come in two distinct flavors.

First, there are ​​body forces​​. These are mysterious, "action-at-a-distance" forces that act on the entire bulk of the cube simultaneously. Gravity is the most familiar example. Drop an apple, and every atom in the apple is pulled downward. If you were to magically double the size of our cube, the gravitational force on it would also double. These forces are proportional to the mass (or volume) of the element. We typically write the total body force on a volume as the integral of a force density, $\rho \mathbf{b}$, where $\rho$ is the mass density and $\mathbf{b}$ is the body force per unit mass. For gravity near the Earth's surface, $\mathbf{b}$ is simply the acceleration vector $\mathbf{g}$.

The second type of force is much more intimate: ​​surface forces​​. These are the contact forces exerted by the material directly touching our tiny cube. The cube to the left pushes on our left face; the cube above pulls on our top face, and so on. This is the mechanical essence of internal togetherness, what we call ​​stress​​. Unlike a body force, a surface force is proportional to the area of the face it acts on. The force per unit area is defined as stress, denoted by the symbol $\boldsymbol{\sigma}$.

Now, stress is a more subtle character than the simple pressure you learned about in high school. Pressure is a force that always pushes perpendicular to a surface. Stress can do that too—we call this ​​normal stress​​—but it can also act parallel to a surface, trying to shear it, like the force between cards in a deck when you slide the top card. This is called ​​shear stress​​. Because the force on a face depends on both its magnitude and its orientation in space, we need a mathematical object more powerful than a simple number or vector to describe it. This object is a ​​tensor​​. For now, you can think of the stress tensor $\boldsymbol{\sigma}$ as a machine that, when you tell it the orientation of a surface (via a normal vector $\mathbf{n}$), tells you the force vector $\mathbf{t}$ acting on that surface. This relationship, $\mathbf{t} = \boldsymbol{\sigma}\mathbf{n}$, is known as Cauchy's stress principle.

We are now ready to apply Newton's second law to our tiny cube. The right-hand side of $F=ma$, the $ma$ part, is the mass of the cube ($\rho$ times its volume) multiplied by its acceleration, $\mathbf{a}$. Let's call this the ​​inertial term​​.

The "F" part is the sum of all forces. We have the body force, $\rho \mathbf{b}$, acting throughout the cube's volume. And we have the surface forces from the stress acting on all six faces. Now, here is the critical insight. If the stress were perfectly uniform everywhere, the push on the left face would be exactly canceled by the push on the right face. A net force only arises if there is an imbalance—a change in stress from one side of the cube to the other. For example, if the pressure on the bottom of a submerged cube is higher than on the top, you get a net upward buoyant force.

This net force resulting from the spatial change in stress is captured by a beautiful mathematical operation called the ​​divergence of the stress tensor​​, written as $\nabla \cdot \boldsymbol{\sigma}$. It measures how much the stress is "unbalanced" at a point.

Putting it all together, Newton's law for our infinitesimal cube becomes:

$\rho \mathbf{a} = \nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b}$

This is it! This is ​​Cauchy's equation of motion​​. Take a moment to appreciate it. It is a vector equation, standing for three equations in three dimensions, and it is profoundly beautiful in its symmetry and simplicity. Let's look at what it tells us:

• ​​$\rho \mathbf{a}$​​: This is the response of the material at a point. It's inertia—the reluctance of a small parcel of mass to change its velocity.
• ​​$\nabla \cdot \boldsymbol{\sigma}$​​: This is the net force from the internal world—the combined push and pull from all of its immediate neighbors.
• ​​$\rho \mathbf{b}$​​: This is the force from the external world—the long-range pull of gravity or other fields.

The equation says that the inertia of a point-mass is overcome by the sum of the net force from its neighbors and any force from the universe at large. And notice, as a crucial sanity check, that every single term in this equation has the same physical dimensions: force per unit volume. Nature's bookkeeping is always perfect.

There is one more piece to the puzzle. We have balanced the forces to make sure our tiny cube doesn't fly off in some direction. But what about rotation? We must also balance the torques to prevent it from spinning.

Imagine the shear stresses on the four side faces of our cube in 2D. If the upward shear on the right face and the downward shear on the left face were not balanced by the shear forces on the top and bottom faces, they would create a net torque. For an infinitesimally small cube, this would lead to an infinite angular acceleration, which is physically absurd. The only way to ensure that there are no net internal torques on any infinitesimal element is for the stress tensor to be ​​symmetric​​. This means that the shear stress on the x-face in the y-direction must equal the shear stress on the y-face in the x-direction ($\sigma_{xy} = \sigma_{yx}$), and so on. This is ​​Cauchy's second law of motion​​ and it's a direct consequence of the conservation of angular momentum. This elegant symmetry is not an assumption; it is a deep truth about the nature of internal forces.

Cauchy's equation is a universal tool. Its form never changes, but its terms take on different meanings in different situations, allowing it to describe a vast range of phenomena.

Let's consider ​​statics​​, where nothing is accelerating, so $\mathbf{a} = \mathbf{0}$. The equation simplifies to an equilibrium condition: $\nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b} = \mathbf{0}$. This equation holds the secret to how a bridge supports its load or a mountain supports its own weight. Let's apply it to a fluid at rest, where the only stress is isotropic pressure, so $\boldsymbol{\sigma} = -p\mathbf{I}$ (where $\mathbf{I}$ is the identity tensor). The divergence becomes $-\nabla p$. If the only body force is gravity, $\mathbf{b} = \mathbf{g} = (0, 0, -g)$, our grand equation of motion elegantly reduces to $\nabla p = \rho \mathbf{g}$. This tells us that pressure increases with depth, which gives us the familiar high-school physics formula $\frac{dp}{dz} = -\rho g$. From a single universal principle, a familiar result appears as a simple special case!

Now for ​​dynamics​​, where acceleration is key. Think of a skyscraper swaying in an earthquake. The term $\rho \mathbf{a}$ is now all-important. The French mathematician d'Alembert had a wonderfully clever idea. Why not just move the inertia term to the other side of the equation? $\nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b} - \rho \mathbf{a} = \mathbf{0}$ Look at this! It has the exact form of a statics problem. We can think of dynamics as a "quasi-static" problem where we have simply added a new, fictitious body force, $-\rho \mathbf{a}$, which we call the ​​inertial force​​. It’s the "force" you feel pushing you back in your seat when a car accelerates. This brilliant trick allows us to use all the mathematical tools of statics to solve complex dynamics problems.

Cauchy's equation provides the universal law of the interior. But to solve a real-world problem, we need to specify what's happening at the edges. These are the ​​boundary conditions​​. For any piece of the surface of our body, we have two fundamental choices:

1. ​​Prescribe the Position:​​ We can clamp a part of the boundary, forcing its displacement to be a known value, $\mathbf{u} = \bar{\mathbf{u}}$. This is a displacement or Dirichlet boundary condition.
2. ​​Prescribe the Force:​​ We can apply a known force or traction to a part of the boundary, $\boldsymbol{\sigma}\mathbf{n} = \bar{\mathbf{t}}$. This is a traction or Neumann boundary condition.

Here is a supremely important point: you cannot, in general, do both on the same piece of boundary. Why? Because the physics has its own logic. If you specify the displacement of the boundary, Cauchy's equation will calculate for you the unique stress field throughout the body that is consistent with that state. This, in turn, determines the reaction force at the boundary. You are no longer free to demand that the force be something else. It would be like telling your thermostat to maintain the room at 70 degrees, and also demanding that the furnace run for exactly 5 minutes every hour. The system can’t obey both commands if they are contradictory. You set one, and the system determines the other.

Finally, for a dynamic problem, we also need to know the state of the system at the starting gun. This means we must specify the initial displacement and the initial velocity of every point in the body. With the governing equation, the boundary conditions, and the initial state all laid out, the fate of the continuum is sealed. Its entire subsequent motion, the dance of stresses and strains within it, is uniquely determined. This is the astonishing predictive power of physics, all stemming from a thoughtful application of Newton's laws to a speck of continuous matter.

You might be thinking, after our journey through the elegant architecture of Cauchy’s equation of motion, "This is all very beautiful, but what is it for?" It’s a fair question. The true power and beauty of a physical law lie not just in its logical coherence, but in its ability to reach out and describe the world around us. And in this, Cauchy’s equation is a supreme artist.

This single statement of momentum balance, a sort of $F=ma$ for continuous stuff, is a master key that unlocks secrets across an astonishing spectrum of science and engineering. It governs the silent, immense forces that hold a bridge against the wind, the vibrations that travel through the Earth after an earthquake, the strange flow of toothpaste from a tube, and even the motion of galaxies in our expanding universe. The equation itself remains the same; what changes is the physical context and the material’s unique “personality,” which we describe with a constitutive law. Let’s take a tour of this vast landscape and see the master key at work.

Much of modern civilization is built—quite literally—on our ability to predict how materials respond to forces. We want our dams to hold, our engines to run, and our airplanes to fly. At the heart of this predictive power lies Cauchy’s equation in its simplest form: the state of equilibrium.

When things are not accelerating, the right-hand side of the equation is zero, leaving us with a statement of perfect balance: $\nabla \cdot \boldsymbol{\sigma} + \rho\mathbf{b} = \mathbf{0}$. This says that for any tiny cube of material, the net "tugging and pushing" from its neighbors ($\nabla \cdot \boldsymbol{\sigma}$) must exactly cancel out any body forces like gravity (the $\rho\mathbf{b}$ term). If they didn't, the cube would start to move.

Consider the design of a thick-walled pressure vessel or a pipe, essential for everything from power plants to plumbing. By applying Cauchy's equilibrium equation in the cylindrical symmetry of the pipe, we can precisely calculate the distribution of radial and hoop (circumferential) stresses. This allows an engineer to know the point of maximum stress and ensure the material is strong enough to withstand the internal pressure without failing. It’s a textbook case of turning a fundamental physical law into a tool for safety and design.

The principle extends from simple pipes to the complex structures that define our skylines. Take a cantilever beam, like a balcony or a bridge segment, supporting its own weight. If you apply the local equilibrium equation to every infinitesimal piece of the beam and then integrate—or "sum up"—these effects over the beam's cross-section, you don't get a mess. Instead, a beautiful simplification occurs: the detailed continuum equations miraculously transform into the classic engineering formulas for shear force and bending moment that structural engineers use every day. This isn't a coincidence; it’s a profound demonstration of how macroscopic engineering theories are built upon the bedrock of continuum mechanics.

Of course, not everything is static. What about objects in motion, like a spinning flywheel or a turbine disk in a jet engine? Here, the inertial term $\rho\mathbf{a}$ comes alive. For a disk rotating at a constant angular velocity $\omega$, every piece of the material is constantly accelerating towards the center. This centripetal acceleration, $\mathbf{a} = -\omega^2 r \mathbf{e}_r$, acts like an effective body force. Cauchy's equation handles this seamlessly. By solving the equation of motion with this inertial term, engineers can calculate the stress distribution caused by rotation, preventing the disk from flying apart at high speeds.

When forces are not in balance, things move, and often, they oscillate. This is where Cauchy's equation sings. It is the parent equation for nearly all waves that travel through solid matter. The term $\rho \ddot{\mathbf{u}}$—density times acceleration—is the engine of dynamics.

Let’s start with the simplest case: a one-dimensional elastic bar. If you strike one end, a compression wave travels down its length. By applying Cauchy’s law to an infinitesimal slice of the bar, along with a simple constitutive law stating that stress is proportional to strain (Hooke's Law), we arrive directly at the classic one-dimensional wave equation: $\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}$. The wave speed $c$ is determined entirely by the material's stiffness and density, $c = \sqrt{E/\rho}$. This is a beautiful result: the complex dynamics of countless atoms are captured by a single emergent property—the speed of sound.

This idea scales up magnificently to our entire planet. The Earth is, to a good approximation, a giant elastic sphere. When an earthquake occurs, it doesn't just shake the ground locally; it sends waves reverberating through the planet's interior. Cauchy's equation for elastodynamics, which combines momentum balance with the laws of elasticity, is the governing principle. For an isotropic material (one whose properties are the same in all directions), the theory predicts two fundamental types of waves that can travel through the bulk of the material.

• ​​Primary (P) waves:​​ These are longitudinal waves, where the particles of the medium oscillate back and forth in the same direction as the wave is traveling, just like sound waves. They are the fastest seismic waves, with a speed $c_P = \sqrt{(\lambda+2\mu)/\rho}$.
• ​​Secondary (S) waves:​​ These are transverse waves, where particles oscillate perpendicular to the direction of wave travel, like flicking a rope. They are slower, with a speed $c_S = \sqrt{\mu/\rho}$.

The existence of these two distinct wave speeds is not an assumption; it is a direct mathematical consequence of Cauchy's equation of motion. And it has a profound implication: the shear modulus $\mu$ appears in the S-wave speed, but not so simply in the P-wave speed. Since liquids have no shear stiffness ($\mu=0$), they cannot support S-waves. By observing that S-waves from an earthquake on one side of the Earth do not reach the other side, seismologists deduced that the Earth must have a liquid outer core—a stunning discovery about a place we can never visit.

This leads to an even more powerful idea: the inverse problem. Not only can we predict wave speeds from material properties, but we can also infer material properties by measuring wave speeds. By timing the arrival of P- and S-waves at seismograph stations around the globe, scientists can map out the density and elastic moduli ($\lambda$ and $\mu$) of the Earth’s mantle and core, painting a detailed picture of our planet's inner structure.

The world, however, is rarely so simple. Many materials, from wood to modern composites to the crystals in rock, are anisotropic—their properties depend on direction. Cauchy's equation, coupled with a more complex anisotropic constitutive law, rises to the challenge. When we seek plane wave solutions in such a medium, the equation transforms into an algebraic problem called the Christoffel eigenproblem. Its solution reveals that for a given direction of travel, there can be up to three different wave speeds with distinct polarizations, a richer and more complex behavior that is essential for materials science, crystallography, and more accurate seismic modeling.

The true universality of Cauchy’s equation becomes apparent when we step outside the familiar world of elastic solids. The equation describes any continuum, as long as we can define its stress.

What about fluids? The famous Navier-Stokes equations, which govern everything from water flowing in a pipe to the air over an airplane wing, are nothing more than a specific instance of Cauchy’s equation of motion. The only difference is the constitutive law: for a simple fluid, stress is related not to the amount of strain, but to the rate of strain. If we swap out the constitutive law for that of a more exotic substance, like a Bingham plastic which behaves like a solid until a certain yield stress is exceeded, Cauchy's equation still holds. It can then describe the flow of mud, wet cement, or even toothpaste squeezing from a tube, predicting where the material will flow and where it will remain as a rigid "plug".

The interconnectedness of physics provides an even wider stage. Physical processes are rarely isolated. What happens when a body heats up? It expands, creating internal strains and stresses. Conversely, deforming a body can generate heat. This is the domain of thermoelasticity, where mechanics and thermodynamics are inextricably linked. The full theory requires two coupled equations: Cauchy’s equation for momentum balance, where the stress now includes a term due to temperature changes, and the heat equation for energy balance. Solving this coupled system is crucial for designing structures that undergo significant temperature variations, such as in jet engines or spacecraft re-entering the atmosphere.

The mathematical language of Cauchy's equation, in its most elegant form, is that of tensors, which allows it to be expressed independently of any chosen coordinate system. This gives it tremendous power. Imagine a fluid flowing not on a flat plane, but constrained to a curved surface, like a soap film on a wireframe. The fundamental principle of momentum balance still applies, but to write it down, we must use the tools of differential geometry, like covariant derivatives and Christoffel symbols, which account for the curvature of the space. The resulting equation looks more complex, but its physical heart is unchanged. It’s a beautiful glimpse into how physical laws retain their integrity even when the stage they play on is bent and warped—a prelude to the concepts of Einstein’s General Relativity.

Perhaps the most breathtaking application of all takes us to the largest possible scale: the cosmos itself. On an astronomical scale, galaxies and gas clouds can be treated as a fluid, governed by Cauchy's equation. But the universe is expanding. The distance between any two points is stretching over time. To make sense of motion in this context, cosmologists use "comoving coordinates" that expand along with the universe. If we start with Cauchy’s equation in our familiar physical coordinates and painstakingly transform it into this expanding frame, a remarkable new term appears. The equation for the "peculiar" velocity of a galaxy (its motion relative to the general expansion) gains a term that acts like a drag force, $-2H(t)\mathbf{u}$, where $H(t)$ is the Hubble parameter describing the rate of cosmic expansion. This "Hubble friction" isn't a real dissipative force; it is a purely inertial effect—the consequence of viewing a local law of motion from the perspective of a globally expanding spacetime. The momentum of a particle gets diluted as space itself stretches.

From the mundane to the magnificent, from designing a pipe to describing the cosmos, Cauchy's equation of motion stands as a testament to the unifying power of physics. It reminds us that the same fundamental principles choreograph the dance of atoms, the shaking of planets, and the stately procession of the galaxies.