Objective Stress Rates

In the study of how materials deform, a fundamental challenge arises when movements are large and complex. How do we measure the true change in stress within a material that is simultaneously being stretched and spun? This question is central to continuum mechanics and highlights a critical flaw in simplistic approaches: the standard rate of stress change is not independent of the observer's viewpoint. This failure to uphold the ​​Principle of Material Frame Indifference​​ leads to physically meaningless models, especially in scenarios involving large rotations. This article confronts this problem head-on. First, in the "Principles and Mechanisms" chapter, we will delve into the concept of objectivity, demonstrate why the simple time derivative of stress fails, and explore the elegant solution of ​​co-rotational derivatives​​ or objective stress rates. Following this theoretical foundation, the "Applications and Interdisciplinary Connections" chapter will reveal why these concepts are not just academic but are essential for the practical world of computational engineering, influencing everything from the accuracy of finite element simulations to the prediction of structural stability.

Imagine you are stirring a very thick, cold batch of taffy. It's hard work. The taffy resists being stretched and deformed—it is under ​​stress​​. But it's also being carried along by the circular motion of your spoon. If you were a tiny observer floating in the taffy, you would be spinning around. From your spinning perspective, how would you measure the "true" change in stress caused by the taffy being pulled, as opposed to the apparent change caused by your own rotation? This simple question plunges us into one of the most subtle and elegant concepts in continuum mechanics: the need for ​​objective stress rates​​.

At the heart of physics lies a beautiful principle of symmetry: the laws of nature should not depend on the observer. A material's intrinsic behavior—how it deforms in response to a force—must be the same whether it's observed from a stationary lab or from a spinning carousel. This idea is called the ​​Principle of Material Frame Indifference​​, or ​​objectivity​​. It demands that our constitutive equations, the mathematical rules that describe material behavior, must retain their form under any rigid-body motion of the observer.

Let's see what this means for stress. The stress within the taffy is described by a mathematical object called the Cauchy stress tensor, $\boldsymbol{\sigma}$. If a new observer is rotating with respect to you, described by a rotation tensor $\boldsymbol{Q}(t)$, they will measure a different stress tensor, $\boldsymbol{\sigma}^*$, whose components are simply the rotated components of your tensor: $\boldsymbol{\sigma}^* = \boldsymbol{Q}\boldsymbol{\sigma}\boldsymbol{Q}^T$. This is perfectly reasonable; you're both looking at the same physical state of stress, just from different angles.

The trouble begins when we consider how stress changes with time. We might naively think that the rate of change of stress is just its time derivative, $\dot{\boldsymbol{\sigma}} = d\boldsymbol{\sigma}/dt$. But is this rate "objective"? Let's check by taking the time derivative of the transformation rule: $\dot{\boldsymbol{\sigma}}^* = \frac{d}{dt}(\boldsymbol{Q} \boldsymbol{\sigma} \boldsymbol{Q}^T) = \dot{\boldsymbol{Q}} \boldsymbol{\sigma} \boldsymbol{Q}^T + \boldsymbol{Q} \dot{\boldsymbol{\sigma}} \boldsymbol{Q}^T + \boldsymbol{Q} \boldsymbol{\sigma} \dot{\boldsymbol{Q}}^T$ If we define the observer's spin as the skew-symmetric tensor $\boldsymbol{\Omega} = \dot{\boldsymbol{Q}}\boldsymbol{Q}^T$, this equation becomes: $\dot{\boldsymbol{\sigma}}^* = \boldsymbol{Q} \dot{\boldsymbol{\sigma}} \boldsymbol{Q}^T + \boldsymbol{\Omega}\boldsymbol{\sigma}^* - \boldsymbol{\sigma}^*\boldsymbol{\Omega}$ Look at those extra terms! The rate of stress seen by the new observer, $\dot{\boldsymbol{\sigma}}^*$, is not just the rotated version of your rate, $\boldsymbol{Q} \dot{\boldsymbol{\sigma}} \boldsymbol{Q}^T$. It has additional pieces that depend on the observer's spin $\boldsymbol{\Omega}$. This means the simple material time derivative, $\dot{\boldsymbol{\sigma}}$, is ​​not objective​​. It mixes the true, physical change in stress due to deformation with a spurious, kinematical change due to the observer's rotation. Using $\dot{\boldsymbol{\sigma}}$ in a constitutive law would be like trying to measure the growth of a plant with a ruler that shrinks and spins randomly. The numbers would be meaningless.

So, how do we fix this? The answer is a stroke of genius. If the problem is rotation, let's subtract it out! We need to invent a new kind of derivative, an ​​objective stress rate​​, that measures the change in stress from the perspective of an observer who is "riding along" with the material, co-rotating with it at every instant. This is a ​​corotational rate​​.

The most intuitive of these is the ​​Jaumann stress rate​​. The idea is to use the material's own local rate of rotation, called the ​​spin tensor​​ $\boldsymbol{W}$, to define the corotating frame. The spin tensor is simply the skew-symmetric (or rotational) part of the velocity gradient, $\boldsymbol{W} = \frac{1}{2}(\boldsymbol{L} - \boldsymbol{L}^T)$. The Jaumann rate is then defined as: $\overset{\nabla}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W}$ The terms we've added, known as the Lie bracket or commutator, are the "correction" that precisely accounts for the rate of change of the stress tensor's components due to the material's rotation. It's a beautiful piece of mathematical engineering. To see it work, let's revisit the test case of a pre-stressed body undergoing a pure rigid rotation. In this case, there is no deformation, only spin, so the velocity gradient is just $\boldsymbol{L} = \boldsymbol{W}$. The stress tensor physically rotates with the body, which means its time derivative is exactly $\dot{\boldsymbol{\sigma}} = \boldsymbol{W}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{W}$. Plugging this into the Jaumann rate definition gives: $\overset{\nabla}{\boldsymbol{\sigma}} = (\boldsymbol{W}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{W}) - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W} = \boldsymbol{0}$ It works perfectly! The objective rate is zero, telling us correctly that no new stress is being generated by the deformation—because there isn't any deformation. The stress is just being passively rotated. This satisfies our physical intuition and the principle of objectivity.

If only it were so simple. It turns out that the notion of "the material's rotation" is not unique. Is it the instantaneous spin of the velocity field, $\boldsymbol{W}$? Or is it the rotation of the underlying material structure itself, described by the rotation tensor $\boldsymbol{R}$ from the polar decomposition of the deformation gradient ($\boldsymbol{F} = \boldsymbol{RU}$)? Different choices for the corotating frame lead to a whole "zoo" of different objective rates.

Besides the ​​Jaumann rate​​, other famous members of this zoo include:

• The ​​Green-Naghdi rate​​, which uses the spin of the material rotation tensor $\boldsymbol{R}$. This rate is considered by many to be more physically tied to the material's actual orientation.
• The ​​Truesdell rate​​, which isn't purely corotational but is derived by relating the rate of change back to the reference configuration. It has a fundamentally different structure.

For very small deformations, these different rates give nearly identical results. Differences only become significant when we have large rotations combined with deformation. And the differences can be dramatic! A famous example is the prediction of stress in a simple shearing motion. If we build a simple rate-based (​​hypoelastic​​) material model using the Jaumann rate, it predicts that the shear stress will oscillate in an unphysical way as the shear increases. A model using the Green-Naghdi rate, however, predicts a much more reasonable, monotonic increase in stress. This shows that the choice of objective rate is not just an academic debate; it has real, measurable consequences for predicting how a material will behave.

This brings us to a much deeper question. The whole idea of elasticity is rooted in the storage and release of energy, like in a spring. When you stretch a truly elastic material and then let it go back to its original shape, the net work done should be zero. The energy you put in is stored and then given back. Materials with this property are called ​​hyperelastic​​, and their behavior can be derived from a single scalar function called the ​​stored energy potential​​, $\psi$. In this elegant picture, the stress is simply a derivative of the potential with respect to a measure of strain (e.g., $\boldsymbol{S} = 2 \frac{\partial \psi}{\partial \boldsymbol{C}}$, where $\boldsymbol{C}$ is the right Cauchy-Green strain tensor).

Crucially, in hyperelasticity, you don't need objective stress rates at all! Objectivity is automatically satisfied by defining the potential $\psi$ as a function of an objective strain tensor like $\boldsymbol{C}$, which is itself invariant to the observer's rotation. The stress at any moment is determined solely by the deformation at that moment, not by the path taken to get there.

So, can our rate-based ​​hypoelastic​​ models (like $\overset{\nabla}{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$) be integrated to find an underlying energy potential? Is the work done path-independent? For most choices of objective rates, including the popular Jaumann rate, the stunning answer is ​​no​​. You can deform the material in a closed loop (e.g., a cycle of simple shear) and find that the model predicts a net creation or destruction of energy. This is because the "stress power one-form" is not mathematically "exact"—it cannot be written as the total derivative of a potential function.

There is, however, a beautiful exception that proves the rule. There exists a special objective rate, the ​​logarithmic rate​​, which is integrable when used with a constant elasticity tensor. It corresponds to a hyperelastic model where the energy is a quadratic function of the logarithmic (or Hencky) strain. This reveals a profound and subtle unity between kinematics (the choice of rate) and thermodynamics (the existence of a potential).

From the simple problem of describing a stirred pot of taffy, our journey has led us through the fundamental principle of objectivity, the clever invention of corotational rates, and the practical consequences of choosing among them. Finally, it has shown us the deep connection between these rate-based models and the more fundamental, energy-based framework of hyperelasticity. It's a testament to how in physics, a careful and honest analysis of a seemingly simple observation can unfold into a rich and beautiful theoretical structure, one that is not only intellectually satisfying but also essential for the practical engineering of our world—from designing car tires to simulating the behavior of biological tissues.

Now that we have grappled with the principle of material frame indifference and the elegant mathematical tools known as objective stress rates, you might be wondering, "Where does the rubber meet the road?" It's a fair question. This is where physics ceases to be an abstract exercise and becomes a powerful lens through which we can understand, predict, and engineer the world. The need for an objective rate isn't some esoteric footnote for theorists; it is the very key that unlocks our ability to accurately simulate the complex dance of materials under large deformations. Let's embark on a journey to see where this key fits.

In an introductory mechanics course, you likely never encountered a "Jaumann rate" or a "Green-Naghdi rate." And for good reason! When things bend and twist by only a tiny amount, the world is a much simpler place. Imagine a small block of steel being gently loaded. The rotations are minuscule. In this regime, the spurious stress changes induced by rotation are incredibly small compared to the stress changes caused by actual stretching and straining. A careful scaling argument shows that the rotational effects are of a higher order of smallness, and we can, with a clear conscience, ignore them. The simple material time derivative, $\dot{\boldsymbol{\sigma}}$, works just fine.

But the real world is not always so gentle. Think of a car crash, the forging of a turbine blade, or the slow, immense folding of tectonic plates. Here, rotations can be enormous, even if the material itself is only stretching slightly. In these scenarios, the "negligible" rotational terms in the stress rate can become as large, or even larger, than the terms due to physical straining. To ignore them would be to live in a fantasy world where simply spinning an object could magically generate real stresses. Our physical models would fail catastrophically. Thus, the objective stress rate is born out of necessity—it is the physicist's tool for "subtracting out the merry-go-round" to see only the true, physical evolution of stress.

The most profound impact of objective stress rates is in the realm of computational mechanics, the art and science of simulating physical phenomena on a computer. Modern engineering and science rely on powerful tools like the Finite Element Method (FEM) and the Material Point Method (MPM) to model everything from the safety of a bridge to the flow of a glacier. These methods work by breaking a problem down into small, discrete steps in time. At each step, a fundamental question must be answered: given a small change in deformation, how does the stress change?

This is precisely where our objective rates come into play. A naive, rate-based model of a material's elastic response, known as a ​​hypoelastic model​​, would state that the objective stress rate is proportional to the rate of deformation: $\widehat{\boldsymbol{\sigma}} = \mathbb{C} : \boldsymbol{D}$. But which objective rate, $\widehat{\boldsymbol{\sigma}}$, should we choose?

It turns out that the choice matters immensely. The conceptually simplest and most traditional choice is the ​​Jaumann rate​​, which uses the continuum spin tensor $\boldsymbol{W}$ to define its corotating frame. For a long time, it was the workhorse of computational plasticity. Yet, it harbors a subtle flaw. When subjected to a simple, continuous shearing motion—like sliding a deck of cards—a hypoelastic model using the Jaumann rate can predict unphysical oscillating normal stresses. This is a jarring symptom that the model, despite being "objective," is not thermodynamically consistent. It doesn't correspond to a true stored energy potential, and it can predict that energy is created or destroyed in purely elastic cycles.

This deficiency led to a fascinating intellectual journey in mechanics. One alternative, the ​​Green–Naghdi rate​​, uses a spin derived from the polar decomposition of the deformation gradient ($\boldsymbol{F}=\boldsymbol{R}\boldsymbol{U}$). This aligns the rate more closely with the material's actual rotational history and cures some of the Jaumann rate's pathologies.

The ultimate resolution, however, came from a different perspective entirely. Instead of patching a rate-based law, why not build a model that is correct from the ground up? This leads to ​​hyperelastic-plastic models​​. In this framework, the elastic part of the response is derived from a true stored energy function, $\psi$, which depends only on objective measures of elastic strain. Objectivity is baked into the very foundation of the model. Stress is no longer the integral of a rate; it is a direct function of the current deformed state. Consequently, explicit objective stress rates are no longer needed for the elastic law. This approach is more complex to implement but is physically and thermodynamically superior. The ​​logarithmic rate​​ can be seen as the rate-based formulation that is energetically consistent with this hyperelastic viewpoint, forming a "best of both worlds" approach that combines the structure of a rate-form update with the thermodynamic rigor of a potential-based model.

The application of these ideas is not confined to the traditional grid-based FEM. In the ​​Material Point Method (MPM)​​, a powerful technique for simulating problems with massive deformations like landslides, explosions, and fluid-solid interaction, the same principles hold. Stresses are carried by moving particles, and their state must be updated at each time step using information from a background grid. To do this correctly in the face of the large rotations and flows that MPM is designed for, the particle's stress must be updated using an objective rate computed from the local velocity gradient. Again, the principle of objectivity proves to be a universal requirement for any valid simulation of continuum mechanics.

The principle of material frame indifference is not just a rule for simple elastic solids. Its reach extends across the vast landscape of material behaviors.

Consider a metal being bent back and forth. It gets harder to deform. This phenomenon, called ​​hardening​​, is an internal memory of the material's plastic deformation history. In many advanced models, this memory is stored in internal variables, such as a "backstress" tensor $\boldsymbol{\alpha}$ for kinematic hardening. But what is this backstress? It is a tensorial quantity representing a shift of the center of the yield surface in stress space. As the material element rotates, the backstress must rotate with it. Therefore, just like the Cauchy stress $\boldsymbol{\sigma}$, the rate of change of the backstress $\boldsymbol{\alpha}$ must also be objective to avoid creating a fictitious internal state change from pure rotation. The principle applies not just to what we see on the outside (stress) but also to the hidden internal state of the material.

Let's venture even further, into the world of ​​viscoelasticity​​—the realm of polymers, biological tissues, and even the Earth's mantle. These materials exhibit a fascinating combination of solid-like elasticity and fluid-like viscosity. A classic way to model them is the generalized Maxwell model, which pictures the material as a parallel set of springs and dashpots, each with a different relaxation time. When we formulate such a model for large deformations, the same problem reappears. Each elastic element (spring) in the model stores stress, and its state must be described in a frame-indifferent way. Modern theories do this using a multiplicative decomposition of the deformation, where the evolution of stress in each "Maxwell branch" is governed by an objective rate equation. This allows us to build sophisticated models of rheology that are valid for the intense processing conditions of polymers or the slow, powerful creep of geological formations, and correctly reduce to the familiar Prony-series representation in the small-strain limit. From ductile metals to gooey polymers, the same fundamental principles provide the framework for a correct description.

Perhaps the most dramatic application of these ideas lies in the prediction of structural stability. When does a column buckle? When does an airplane wing begin to flutter? These are questions of life and death, and their answers are hidden in the mathematics of the system's response to a small perturbation.

In a computational simulation, this response is governed by the ​​algorithmic tangent stiffness matrix​​, $\mathbf{K}_T$. This matrix is the discrete version of the material's tangent modulus; it dictates how the internal forces change in response to a small change in displacement. The stability of the structure is tied to the eigenvalues of this matrix. A static instability, like buckling, occurs when an eigenvalue becomes zero.

Here is the beautiful connection: the choice of objective stress rate directly influences the structure of this crucial matrix. A thermodynamically consistent formulation, such as a hyperelastic one (or a hypoelastic one using a work-conjugate pair like the logarithmic rate), leads to a symmetric tangent stiffness matrix $\mathbf{K}_T$. This is computationally convenient, but more importantly, it implies that instabilities will be static (buckling).

However, if one uses a less ideal formulation, like a hypoelastic law with the Jaumann rate, the resulting consistent tangent matrix $\mathbf{K}_T$ is generally non-symmetric. This non-symmetry is not just a numerical inconvenience; it is a signature of a non-conservative system. A system with a non-symmetric tangent can exhibit dynamic instabilities, or ​​flutter​​, where oscillations grow exponentially in time. A standard stability analysis designed for symmetric systems would completely miss such a failure mode.

Think about that! The subtle choice of how to define a "rate of change" for stress in the presence of rotation has profound consequences, determining not just the accuracy of a stress calculation, but the very nature of the predicted structural collapse—the difference between a graceful buckling and a catastrophic, violent flutter. It is a stunning example of how a seemingly small detail in our fundamental physical description can have macroscopic, system-level repercussions. The journey from the abstract principle of frame indifference to the concrete prediction of a wing's flutter is a testament to the power, unity, and inherent beauty of mechanics.