Green-Naghdi Rate

In fields from geophysics to material engineering, scientists face a common challenge: how to accurately describe materials that undergo large movements and deformations simultaneously. Standard descriptions of change often fail in this regime, leading to models that predict physically impossible behaviors. This discrepancy arises from a fundamental requirement known as the Principle of Material Frame Indifference, or objectivity, which states that the physical laws governing a material cannot depend on the observer's motion. The simple time derivative of stress is not objective and can produce spurious, non-physical results when a material rotates.

This article tackles this foundational problem in continuum mechanics. It introduces the concept of objective stress rates, which are mathematical tools designed to correctly account for large rotations. In the first chapter, "Principles and Mechanisms," we will delve into the theory of objectivity, contrast the widely-used but flawed Jaumann rate with the more physically robust Green-Naghdi rate, and understand the deep connection between rotation, deformation, and stress. The second chapter, "Applications and Interdisciplinary Connections," will then demonstrate why this seemingly theoretical choice has critical, real-world consequences in computational simulation, materials science, and beyond. We begin our journey by exploring the core principles that demand these advanced concepts and the elegant mechanisms developed to address them.

Imagine you are trying to describe the changing shape of a football as it flies through the air, spinning rapidly. It’s a surprisingly tricky business. From your fixed position on the ground, you see a dizzying blur. The ball is deforming slightly as it travels, but most of what you see is just its rapid rotation. Now, imagine you could shrink down and ride on the surface of the ball. From this new vantage point, the world would be spinning, but the football itself would appear mostly still, with only its subtle stretching and squashing being apparent.

This simple thought experiment captures the central challenge of describing deformable materials undergoing large motions. Our physical laws, which relate how a material deforms (strain) to the internal forces it develops (stress), must be true no matter who is observing them. A physicist on the ground and a physicist on a spinning carousel must arrive at the same fundamental description of a piece of rubber being stretched. This powerful idea is known as the ​​Principle of Material Frame Indifference​​, or ​​objectivity​​.

In introductory physics, we describe change using derivatives—the rate of change of position is velocity, the rate of change of velocity is acceleration. It seems natural to describe the rate of change of stress by simply taking its time derivative, which we can write as $\dot{\boldsymbol{\sigma}}$. But here, we hit a wall.

Let's take a block of rubber that is already under some stress, say, it's being squeezed. Now, let's just rotate the block rigidly, without any additional squeezing or stretching. Physically, the stress state within the material should just rotate along with the block. An observer spinning with the block would report that the stress is constant. However, an observer in a fixed laboratory frame sees the components of the stress tensor changing because its orientation is changing. This means $\dot{\boldsymbol{\sigma}}$ is not zero!

If we were to write a naive constitutive law like "the rate of stress is proportional to the rate of deformation" (e.g., $\dot{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$), this law would fail spectacularly. For our rigid rotation, the rate of deformation $\boldsymbol{D}$ is zero, so the law would predict $\dot{\boldsymbol{\sigma}} = \boldsymbol{0}$. But as we saw, $\dot{\boldsymbol{\sigma}}$ is not zero. Our law is broken. It predicts that the stress doesn't rotate with the material, leading to the creation of ​​spurious stresses​​ that don't physically exist. This is not a minor error; it’s a fundamental violation of objectivity. The ordinary time derivative is tainted by the observer's frame of reference; it's simply the wrong tool for the job.

So, how do we fix this? The solution is as elegant as our football analogy suggests: we stop trying to describe the deformation from a fixed frame and instead do our calculations in a frame that rotates with the material. This moving, spinning frame of reference is called a ​​corotational frame​​.

By jumping onto this "mathematical carousel," we conceptually untangle the rigid spinning of the material from its pure deformation. In this special frame, the only changes we see are the real stretches and shears. Here, a simple time derivative makes sense again. The overall procedure becomes:

1. Take the current stress tensor.
2. Rotate it into the corotational frame.
3. Calculate how it changes in this frame due to pure deformation.
4. Rotate the result back to our fixed laboratory frame.

When this procedure is expressed mathematically, it gives rise to what is known as an ​​objective stress rate​​, often denoted by a symbol like $\overset{\diamond}{\boldsymbol{\sigma}}$. All such objective rates take a general form:

Here, $\dot{\boldsymbol{\sigma}}$ is the problematic simple derivative, and the extra terms, $\boldsymbol{\sigma}\boldsymbol{\Xi} - \boldsymbol{\Xi}\boldsymbol{\sigma}$, are the correction. They precisely counteract the non-objective parts of $\dot{\boldsymbol{\sigma}}$. The tensor $\boldsymbol{\Xi}$ is a ​​spin tensor​​; it's a skew-symmetric tensor that describes the instantaneous rate of rotation of the corotational frame. The whole expression $\overset{\diamond}{\boldsymbol{\sigma}}$ can be beautifully interpreted as the time derivative of the stress as seen by an observer in the corotating frame, but expressed in the coordinates of the laboratory frame.

This brings us to a crucial question: how, exactly, do we define the "rotation of the material"? For a rigid body like a carousel, it's obvious. But for a squishy, deforming continuum, different parts might be spinning differently. This ambiguity leads to different definitions of the corotational frame, and thus, different objective stress rates. Let’s meet the two most famous ones.

​​1. The Jaumann Rate:​​ This is perhaps the most straightforward approach. It defines the spin of the corotational frame using the local ​​vorticity tensor​​, $\boldsymbol{W}$. You can picture this as the spin of an infinitesimally small paddle wheel dropped into the "flow" of the deforming material. It measures the average instantaneous rotation of the material at a point. The resulting objective rate, the ​​Jaumann rate​​, is built using this spin $\boldsymbol{\Xi}=\boldsymbol{W}$.

​​2. The Green-Naghdi Rate:​​ This rate is based on a deeper, more powerful mathematical idea: the ​​polar decomposition​​. This theorem states that any deformation, described by the deformation gradient tensor $\boldsymbol{F}$, can be uniquely decomposed into a pure stretch ($\boldsymbol{U}$) followed by a pure rigid rotation ($\boldsymbol{R}$). Think of it like a recipe for getting from the original shape to the final shape: first stretch the material along certain axes, then rotate the whole stretched shape into its final orientation. The equation is simply $\boldsymbol{F} = \boldsymbol{R}\boldsymbol{U}$.

The ​​Green-Naghdi rate​​ defines its corotational frame using this very rotation, $\boldsymbol{R}$. It tracks the orientation of the material's underlying structure as it deforms. The spin it uses, $\boldsymbol{\Omega} = \dot{\boldsymbol{R}}\boldsymbol{R}^T$, is the rate of change of this material rotation.

Now for the critical insight: these two spins, the vorticity $\boldsymbol{W}$ and the material spin $\boldsymbol{\Omega}$, are ​​not the same​​ in general! This fundamental difference has profound consequences.

So, which definition is "better"? Physics often judges theories by testing them in simple, well-understood scenarios. A classic test for material models is ​​simple shear​​, which is what happens when you spread soft butter on toast or slide a deck of cards.

When we use a simple elastic model with the ​​Jaumann rate​​ to simulate a large amount of simple shear, something unphysical happens. Instead of the shear stress building up smoothly, it begins to oscillate, wiggling up and down. In more complex materials, it can even incorrectly predict the generation of stresses in directions where there should be none. It’s as if the model is getting "dizzy" by tracking the local fluid-like vorticity $\boldsymbol{W}$ instead of the true orientation of the material.

The ​​Green-Naghdi rate​​, on the other hand, excels here. Because it is tied to the more fundamental material rotation $\boldsymbol{R}$, it typically produces a smooth, monotonic stress response in simple shear, free of the spurious oscillations. It provides a more physically faithful account of what the material is actually doing.

Interestingly, the choice of rate dramatically affects the rate at which the principal stress axes rotate. The difference between the two rates is a term that dictates how the direction of the internal forces evolves, and the Jaumann rate can make them rotate in a non-physical way.

Why does the Jaumann rate behave so poorly? The issue lies in a deep concept called ​​integrability​​. A truly elastic material, like an ideal spring, has a perfect "memory." The work you do to deform it is stored as potential energy, and you get all of it back when you release it. The final energy depends only on the final state of deformation, not the specific path taken to get there. This is called a ​​path-independent​​ or ​​hyperelastic​​ response.

A rate-based, or ​​hypoelastic​​, law is not guaranteed to have this property. It turns out that a simple hypoelastic law using the Jaumann rate is ​​not integrable​​—it doesn't correspond to any underlying potential energy function. The spurious oscillations in simple shear are a direct symptom of this pathology; the model is creating and destroying energy in unphysical cycles.

The Green-Naghdi rate is a significant improvement, but it too is not perfectly integrable for all types of materials. The search for a "perfectly" integrable rate leads to the ​​logarithmic rate​​, which is connected to a specific measure of strain called Hencky strain. This rate yields a fully path-independent elastic model.

So, where does that leave the Green-Naghdi rate? It represents a crucial step up the ladder of physical fidelity from the simpler Jaumann rate. It provides a robust, physically intuitive, and computationally effective way to ensure objectivity while avoiding the most glaring pathologies. For many practical engineering problems, especially in metals where elastic strains remain small, the differences between these rates are higher-order effects that can be negligible. But in the world of large deformations—in soft robotics, biomechanics, and materials science—understanding this elegant piece of mechanics is not just an academic exercise, but a necessity for building models that truly reflect the beautiful and complex reality of the physical world.

Now that we have grappled with the principles of objectivity and the spinning, stretching nature of deformation, we might be tempted to ask, "So what?" Does the choice between one objective rate and another—say, the Jaumann versus the Green-Naghdi—really matter outside the ivory towers of theoretical mechanics? The answer, it turns out, is a resounding yes. This choice is not merely a mathematical subtlety; it is a profound physical statement with far-reaching consequences that ripple through computational engineering, materials science, and even geophysics. To appreciate this, we must leave the clean world of pure theory and venture into the messy, demanding realm of real-world problems. This is where the true character and utility of a concept like the Green-Naghdi rate are revealed.

In physics, we often have what you might call a "hydrogen atom"—a problem so simple in its setup that its solution, or lack thereof, reveals deep truths about our theories. In solid mechanics, one such problem is simple shear. Imagine shearing a thick book by pushing the top cover sideways. It's a motion we see everywhere, from the grinding of tectonic plates to the cutting of metal.

What happens if we subject a simple, elastic block to this motion and model it with a rate equation? You would expect the force required to increase smoothly as the shear increases. However, if we use the Jaumann rate to define our stress evolution, something utterly bizarre can happen. For large amounts of shear, this model predicts that the shear stress will oscillate! It is as if you are pushing the book with a steady hand, yet it begins to alternately resist you less, then more, then less again, all on its own. This is not some esoteric numerical error; it is a direct prediction of the model. This prediction of "spurious shear oscillations" is, of course, physically nonsensical for a simple elastic material.

Why does this happen? The Jaumann rate is based on the continuum spin, $\boldsymbol{W}$, which for simple shear, is constant. It essentially treats the material as a fluid-like continuum that is spinning. But a solid has an internal structure, a lattice, which is stretching along certain axes. The Green-Naghdi rate, by being tied to the rotation of these principal stretch axes, $\boldsymbol{\Omega} = \dot{\boldsymbol{R}}\boldsymbol{R}^T$, paints a much more faithful picture. In a hypoelastic model based on the Green-Naghdi rate, the predicted stress response in simple shear is smooth and monotonic, just as our intuition demands. The simple shear test, therefore, acts as a crucial proving ground, showing that the Green-Naghdi rate’s perspective—focusing on the material's own rotation—is a more robust foundation for building models of reality.

The most profound impact of these ideas is felt in computational mechanics, the field dedicated to simulating the behavior of materials and structures using computers. Here, the Green-Naghdi rate is not just a 'better' choice; it is often an essential one for obtaining physically meaningful results.

The Menace of Cyclic Loading: Ratcheting and Phantom Energy

Imagine designing a component for an aircraft engine or a building in an earthquake zone. It will be subjected to thousands of cycles of loading and unloading. We need to be able to predict its behavior with absolute confidence. Let's consider a symmetric, back-and-forth simple shear cycle. We apply a certain amount of shear, then reverse it to return to the starting point.

Here again, the Jaumann rate leads us astray in the most dramatic fashion. In a simulation of plastic deformation under this symmetric cyclic loading, a model using the Jaumann rate can predict ratcheting: the material begins to accumulate a net shear stress, drifting further and further with each cycle, even though the applied deformation is perfectly balanced. This is a catastrophic failure of the model, as it would incorrectly predict premature failure where none should exist.

Even more disturbingly, in a purely elastic cycle, the Jaumann model can predict a non-zero net work done. This means the model seems to create (or destroy) energy from nothing over a closed loop, forming an apparent hysteresis loop where there should be none. This is a direct violation of the laws of thermodynamics.

A formulation based on the Green-Naghdi rate, however, behaves beautifully. Because it is built upon the clean separation of rigid rotation from pure stretch (the polar decomposition), it correctly predicts a stable, non-ratcheting stress response in a symmetric cyclic test. For an elastic material, it correctly predicts zero net work over a closed cycle. This robustness is precisely why corotational formulations, of which the Green-Naghdi rate is a prime example, are a cornerstone of modern computational plasticity software for analyzing fatigue and cyclic loading.

The Heart of the Simulation: Plastic Flow and Computational Speed

When a material deforms plastically, like a piece of metal being stamped into a car door, the direction in which the material flows is governed by the state of stress. As we've seen, different objective rates lead to different predicted stress evolutions under large rotations. Consequently, the choice of objective rate directly influences the predicted shape and final state of the formed part. For industries where precision is paramount, this is a critical detail.

Furthermore, the choice of rate has a very practical impact on the efficiency of the simulation itself. In the Finite Element Method (FEM), the governing equations are solved incrementally. At each step, the computer has to solve a massive system of linear equations, a process that is vastly more efficient if the system's "tangent matrix" is symmetric. Formulations based on the Green-Naghdi rate (or its close cousin, the logarithmic rate) can be constructed to produce a symmetric tangent matrix for elastic-plastic materials [@problem_id:2676218, @problem_id:2568886]. Using a non-integrable rate like the Jaumann rate generally produces a non-symmetric matrix, significantly increasing the computational cost and complexity. In essence, choosing a physically well-founded rate like Green-Naghdi not only gives you a better answer, it helps the computer find that answer much, much faster.

The significance of objective rates extends beyond engineering simulations into the fundamental principles that govern all matter.

From Solids to Fluids: The Response to Pressure

Consider a material under immense pressure, but with no shear stress—a state we call hydrostatic. This could be water at the bottom of the ocean or rock deep within the Earth's crust. What happens if this material begins to flow and deform? Will shear stresses immediately appear, or will the stress remain purely pressure-like?

The answer depends on your choice of objective stress rate! A constitutive model based on the Jaumann or Green-Naghdi rate predicts that if you start with a hydrostatic stress, any arbitrary deformation will initially produce a purely hydrostatic rate of stress. The deviatoric (shear) stresses only build up later. In contrast, another perfectly valid objective rate, the Truesdell rate, predicts that shear deformation will instantaneously generate shear stresses from a hydrostatic state. This reveals a deep connection: the choice of objective rate is tied to fundamental assumptions about how a material's shear and volumetric responses are coupled. This has profound implications in geophysics and high-pressure fluid dynamics, where the behavior of materials under extreme compression is a central question.

Stability and the Arrow of Time

Finally, we arrive at the deepest connection of all: thermodynamics. Any valid material model must obey the laws of physics, most notably the second law of thermodynamics, which states that entropy, or disorder, in an isolated system can only increase. For a material undergoing plastic deformation, this manifests as plastic dissipation—the irreversible work that is converted into heat. This dissipation must always be non-negative. A material cannot spontaneously organize itself and give back energy during plastic flow.

This principle is embodied in Drucker's Stability Postulate. For this physical law to be meaningful, it must be stated in a way that is independent of the observer. This is achieved by framing the plastic dissipation rate, $\mathcal{D}$, as the product of work-conjugate quantities: the Cauchy stress $\boldsymbol{\sigma}$ and the plastic rate of deformation $\boldsymbol{d}^{p}$, such that $\mathcal{D} = \boldsymbol{\sigma} : \boldsymbol{d}^{p} \ge 0$.

The choice of objective rate (Jaumann, Green-Naghdi, etc.) resides in the elastic part of the model. It determines the stress path the material takes, but it does not change the fundamental law of dissipation itself. However, a "bad" choice of rate (like Jaumann) can lead to such unphysical stress paths that it becomes difficult to construct models that behave sensibly in all conditions. The robust, physically intuitive nature of the Green-Naghdi rate, rooted in the material's own deformation, provides a much sounder foundation upon which to build models that are guaranteed to respect the fundamental laws of thermodynamics.

In the end, the journey through the applications of the Green-Naghdi rate teaches us a beautiful lesson. The quest for an objective stress rate is a quest for the right perspective. It's about finding the reference frame—not one fixed in space, nor one spinning with the continuum, but one that rides along with the material itself, rotating and stretching in harmony with its internal structure. By adopting this perspective, the physics simplifies, paradoxes vanish, and we are left with a powerful and elegant tool for understanding and predicting the complex dance of matter in motion.