Stress Asymmetry in Hard Magnetic Soft Materials

This paper demonstrates that while continuum magnetomechanical formulations for hard magnetic soft materials yield different Cauchy stress symmetries depending on whether the magnetization is described in the referential or current configuration, both formulations converge to symmetric stresses with identical equilibrium behavior when the magnetization field is at its energy-minimizing state.

The Gist

Imagine you are a chef trying to bake a perfect cake. The cake is made of soft dough (the polymer) with tiny, hard chocolate chips embedded inside (the magnetic particles). This is a Hard Magnetic Soft Material.

Now, imagine you want to describe exactly how this cake behaves when you stretch it, twist it, or when a magnet pulls on the chocolate chips. To do this, scientists use complex math called "continuum mechanics."

This paper is about a very specific, tricky question in that math: Does the way we describe the "chocolate chips" change the way we calculate the forces inside the cake?

Here is the breakdown in simple terms:

1. The Two Ways to Look at the Cake

The authors point out that you can describe the position of the chocolate chips in two different ways, and both seem correct at first glance:

• The "Map" View (Referential Description): Imagine you have a blueprint of the cake before you baked it. You say, "The chocolate chip is at coordinate X on the blueprint." Even if you stretch the cake, you still refer to its original spot on the blueprint.
• The "Real-Time" View (Current Description): Imagine you are looking at the cake while it is being stretched. You say, "The chocolate chip is currently at coordinate Y in the air." As the cake stretches, this coordinate changes.

2. The Big Surprise: The Math Changes the Physics

Usually, in physics, it doesn't matter which view you use; the final answer should be the same. But this paper says: Wait a minute!

When you do the math to figure out the stress (the internal pressure or force) inside the material:

• If you use the "Map" View, the math says the forces are perfectly balanced and symmetrical (like a perfectly round wheel).
• If you use the "Real-Time" View, the math says the forces are asymmetrical (like a wobbly wheel).

The Analogy:
Think of a spinning top.

• If you describe the top's spin based on where it started on the table, the math says it's spinning perfectly straight.
• If you describe the spin based on where it is wobbling right now, the math might say it's tilting to the side.

The paper proves that simply changing your "language" (switching from the Map View to the Real-Time View) changes the calculated stress, even though the physical cake hasn't changed.

3. Why Does This Happen?

The authors explain this with a simple rule: The ingredients are coupled.

The chocolate chips (magnetization) and the dough (deformation) are stuck together. When you stretch the dough, the chips move. When you change your description of the chips, you are mathematically changing how you account for that movement.

It's like calculating the cost of a road trip.

• Method A: You calculate the cost based on the miles on your odometer before you started driving.
• Method B: You calculate the cost based on the miles on your odometer while you are driving.

If you don't adjust your formula correctly for the fact that the car is moving, Method A and Method B will give you different numbers for the "force" of the trip, even though the trip is the same.

4. The "Peace of Mind" Conclusion

So, is one method wrong and the other right? No.

The paper finds a beautiful "sweet spot":

• When the system is calm (Equilibrium): If the magnetic chips have stopped moving and settled into their most comfortable position (energy minimum), both methods give the exact same answer. The forces balance out, and the stress becomes symmetrical again.
• When the system is chaotic (Not Equilibrium): If the chips are still wiggling, spinning, or reacting to a changing magnetic field (like in a robot that is moving fast), the two methods give different answers. One says the stress is symmetrical, the other says it's not.

Why Should You Care?

This matters for engineers building soft robots or artificial muscles.

• If you are designing a robot that moves slowly and settles down, you can use either math model; they will agree.
• But if you are designing a robot that needs to snap, twist, or react instantly to a magnet, you have to be very careful which math model you pick. If you pick the wrong one, you might think the robot is stable when it's actually going to wobble or break.

The Takeaway

The paper is a warning label for scientists: "Be careful how you describe your variables." Just because two math formulas look like they are describing the same thing, they might calculate the internal forces differently. However, once everything settles down, the universe agrees, and the math works out the same.

It's a reminder that in the complex world of soft, magnetic materials, how you look at the problem changes the answer you get—unless you wait for everything to settle.

Technical Summary

1. Problem Statement

Hard Magnetic Soft Materials (HMSM) are composite materials consisting of hard magnetic particles embedded in a soft polymer matrix. They are increasingly used in robotics and soft actuators. Continuum modeling of these materials typically treats the deformation ($\mathbf{x}$) and the magnetization field ($\mathbf{m}$) as primary kinematic variables.

A critical theoretical issue has emerged regarding the symmetry of the Cauchy stress tensor in these formulations. The symmetry of the Cauchy stress is fundamentally linked to the balance of angular momentum and frame invariance (objectivity). Recent literature has debated whether the Cauchy stress in HMSM models is symmetric or asymmetric. The authors identify that this ambiguity arises from the choice of the internal variable description:

• Current Description: Magnetization defined in the deformed (spatial) configuration.
• Referential Description: Magnetization defined in the reference (material) configuration.

The paper investigates whether these energetically equivalent formulations yield identical stress tensors and, specifically, whether the symmetry of the Cauchy stress depends on this choice of variables.

2. Methodology

The authors employ a rigorous continuum mechanics framework based on variational principles and energy minimization.

• General Framework: They define two energy functionals, $E[\mathbf{x}, \mathbf{n}]$ (based on current internal variable $\mathbf{n}$) and $E_0[\mathbf{x}, \mathbf{n}_0]$ (based on referential internal variable $\mathbf{n}_0$). These are related by a change of variables $\mathbf{n} = \hat{\mathbf{n}}(\mathbf{F}, \mathbf{n}_0)$, where $\mathbf{F}$ is the deformation gradient.
• Frame Invariance Analysis: They apply rigid body rotations ($\mathbf{Q} \in SO(3)$) to the current configuration to derive conditions for frame invariance. By differentiating the energy invariance conditions with respect to the rotation parameter, they derive the Noether-type conditions for angular momentum balance.
• Variational Derivatives: The stress tensors are derived as partial derivatives of the energy density with respect to the deformation gradient ($\mathbf{F}$), keeping the internal variable fixed.
  • For the referential formulation: $\mathbf{P}_0 = \frac{\partial W_0}{\partial \mathbf{F}}|_{\mathbf{n}_0}$.
  • For the current formulation: $\mathbf{P} = \frac{\partial W}{\partial \mathbf{F}}|_{\mathbf{n}}$.
• Case Studies: The theory is validated and illustrated using two specific material systems:
  1. Liquid Crystal Elastomers (LCEs): A simpler system with a director vector as the internal variable, including gradient terms.
  2. Hard Magnetic Soft Materials (HMSM): The primary focus, incorporating electromagnetic (EM) energy terms and Maxwell stresses.

3. Key Contributions

The paper makes several fundamental theoretical contributions:

• Demonstration of Stress Asymmetry: The authors prove that energetically equivalent formulations (related by a change of variables) generally yield different Cauchy stress tensors. Specifically, the choice of the internal variable (current vs. referential) dictates the symmetry of the stress.
• Symmetry Conditions:
  • Formulations based on the referential description ($\mathbf{n}_0$) naturally produce a symmetric Cauchy stress ($\text{skw}(\boldsymbol{\sigma}_0) = 0$).
  • Formulations based on the current description ($\mathbf{n}$) generally produce an asymmetric Cauchy stress ($\text{skw}(\boldsymbol{\sigma}) \neq 0$) unless specific equilibrium conditions are met.
• Equivalence at Equilibrium: A crucial finding is that while the stress tensors themselves differ, their divergences (which govern the balance of linear momentum) are identical if and only if the internal variable is at its energy-minimizing equilibrium configuration (i.e., $\frac{\partial W}{\partial \mathbf{n}} = 0$).
• Clarification of "Lagrangian" vs. "Eulerian": The authors clarify that they are not switching between Lagrangian and Eulerian descriptions of the motion (all quantities remain Lagrangian). Instead, the distinction lies in how the internal variable transforms under rigid rotation: the referential description is invariant, while the current description rotates with the body.

4. Key Results

A. General Mathematical Derivation

The relationship between the two stress tensors is derived as:
$$\mathbf{P}_0 - \mathbf{P} = \left( \frac{\partial W}{\partial \mathbf{n}} \right) \otimes \mathbf{n}_0$$
Consequently, the skew-symmetric part of the Cauchy stress in the current formulation is:
$$\text{skw}(\boldsymbol{\sigma}) = -J^{-1} \text{skw}\left( \left( \frac{\partial W}{\partial \mathbf{n}} \right) \otimes \mathbf{n} \right)$$
This term vanishes only when $\frac{\partial W}{\partial \mathbf{n}} = 0$ (equilibrium).

B. Liquid Crystal Elastomers (LCEs)

In the LCE example, the director $\mathbf{n}$ is the internal variable.

• The referential formulation yields a symmetric stress.
• The current formulation yields an asymmetric stress proportional to the torque on the director.
• When the director is in equilibrium (no net torque), the stresses coincide, and the Cauchy stress becomes symmetric.

C. Hard Magnetic Soft Materials (HMSM)

For HMSM, the magnetization $\mathbf{m}$ is the internal variable.

• Referential Formulation: The Maxwell stress derived from the referential magnetization $\mathbf{m}_0$ is symmetric.
• Current Formulation: The Maxwell stress derived from the current magnetization $\mathbf{m}$ is generally asymmetric.
• Physical Implication: The asymmetry in the current formulation represents the internal magnetic torque acting on the material. If the magnetization is not in equilibrium (e.g., during dynamic evolution via Landau-Lifshitz-Gilbert dynamics), this torque manifests as an asymmetric stress tensor.

5. Significance and Implications

• Resolution of Controversy: The paper resolves recent debates in the literature regarding stress symmetry in HMSM. It clarifies that both symmetric and asymmetric stress formulations are mathematically valid but correspond to different choices of internal variables.
• Dynamic vs. Static Behavior:
  • Static/Equilibrium: If the internal variable (magnetization) is relaxed to equilibrium, the choice of formulation does not affect the force balance ($\text{div } \boldsymbol{\sigma}$), and the stress is symmetric.
  • Dynamic/Non-Equilibrium: If the internal variable is evolving (e.g., magnetization lagging behind the field), the formulations diverge. The current formulation captures the internal torque explicitly as stress asymmetry, while the referential formulation absorbs this into the definition of the stress.
• Modeling Guidance: For numerical simulations and theoretical analysis, the choice of formulation depends on the physics being modeled. If the focus is on the evolution of the internal variable (dynamics), the current formulation may be more intuitive but requires handling asymmetric stresses. If the focus is on static equilibrium, the referential formulation offers the mathematical convenience of symmetric stress tensors.
• Frame Invariance: The work reinforces that frame invariance does not mandate a unique stress tensor; rather, it imposes constraints on the energy density that lead to different stress expressions depending on the kinematic variables chosen.

In summary, Güner et al. demonstrate that stress asymmetry is not an error but a feature of formulations using current internal variables, representing the internal torque. This asymmetry vanishes only when the internal variable reaches equilibrium, at which point all formulations converge to the same force balance and symmetric stress state.