🔢 mathematics

A robust mixed finite element formulation for third medium contact

This paper introduces a robust mixed finite element formulation for third medium contact that utilizes an auxiliary-field stabilization with penalty coupling to regularize the problem, demonstrating effectiveness across various interpolation schemes and benchmark scenarios involving large deformations and self-contact.

Original authors: M. Vorwerk, J. Schröder, P. Wriggers

Published 2026-06-29

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Vorwerk, J. Schröder, P. Wriggers

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to simulate how two soft, squishy objects (like rubber blocks) bump into each other on a computer. In the world of engineering simulations, this is notoriously difficult.

The Old Problem: The "Hard Wall" Approach
Traditionally, computer programs treat contact like two people trying to walk through a solid wall. The software has to constantly ask, "Are they touching? Are they overlapping? Stop! Go back!" It has to search for the exact moment of impact and enforce strict rules. If the objects slide, twist, or fold in complex ways (like a rubber band snapping back onto itself), the computer gets confused, the math breaks, and the simulation crashes. It's like trying to herd cats while they are running through a maze.

The New Idea: The "Magic Sponge"
This paper introduces a clever trick called Third Medium Contact. Instead of treating the space between the objects as empty air or a hard wall, the authors fill that gap with a "fictitious medium"—think of it as a super-soft, invisible magic sponge.

Before contact: The sponge is so squishy that the objects can move through it without feeling any resistance. It's like moving through warm air.
During contact: As the objects push together, the sponge gets squeezed. The more you squeeze it, the harder it gets to compress. Eventually, it becomes as hard as a rock, pushing the objects apart.
The Result: The objects never actually "touch" in the computer's eyes; they just compress the sponge between them. This avoids all the tricky "search and stop" rules, making the simulation much smoother and less likely to crash.

The New Problem: The "Sponge Collapse"
There was a catch with this sponge idea. If the objects squished the sponge too hard, the computer's math grid (the mesh) would get distorted, twisted, and eventually collapse, causing the simulation to fail again. It was like trying to squeeze a sponge until it turns into a flat, useless pancake.

The Solution: The "Ghost Guide"
The authors propose a new stabilization method to keep the sponge from collapsing. They introduce a second, invisible field they call a "Ghost Guide" (technically an auxiliary field).

Here is the analogy:
Imagine the sponge is made of tiny, independent tiles. If you squeeze the sponge, some tiles might get squished while their neighbors stay flat. This creates a jagged, unstable mess.

The Ghost Guide is like a smooth, flexible sheet laid over the tiles.
The tiles are loosely tied to this sheet. If a tile tries to squish too much compared to its neighbor, the sheet pulls it back into line.
Crucially, this sheet doesn't need to be a complex, high-tech material. It can be simple and low-resolution (using "low-order" math), which makes the computer calculation fast and cheap.

Why This Matters
The paper tests this idea with several scenarios:

Two Blocks: Pushing two blocks together. The new method allows them to get very close without the math breaking, even if the "sponge" is made of very simple, blocky math tiles.
The Box: A frame that bends and touches itself inside a box. The old methods often failed here because the "sponge" got too twisted. The new "Ghost Guide" keeps the sponge stable, allowing the frame to fold and touch itself naturally.
The C-Shape: A curved beam that bends until its ends touch. The new method handles the twisting and turning smoothly, matching the results of much more complex, older methods but without the headaches.

The Bottom Line
The authors found that if you make the "Ghost Guide" continuous (connected smoothly from one piece to the next), it acts like a strong stabilizer, preventing the simulation from crashing even when the objects are squished into extreme shapes. If you make the guide "discontinuous" (broken into separate pieces per tile), it's easier to compute but doesn't stop the collapse as well.

In short, this paper gives engineers a robust, crash-proof way to simulate objects rubbing, sliding, and folding into each other by replacing hard contact rules with a smart, self-stabilizing "magic sponge" and a simple "Ghost Guide" to keep the math from falling apart.

Technical Summary: A Robust Mixed Finite Element Formulation for Third Medium Contact

Problem Statement
Classical finite element contact formulations rely on explicit discretization of contact interfaces, enforcing constraints via penalty, Lagrange multipliers, or barrier methods. While mature, these approaches face significant algorithmic burdens in scenarios involving large sliding, self-contact, severe mesh distortion, or topology optimization where contact boundaries are unknown a priori. The "third medium contact" (TMC) method addresses this by replacing explicit constraints with a highly compliant fictitious medium filling the space between potentially contacting bodies. This medium remains soft before contact but stiffens drastically under compression, transmitting contact forces through the continuum.

However, a central challenge in TMC is stabilization. The medium must be soft enough to avoid spurious forces pre-contact yet stable enough to withstand extreme compression and shear without mesh distortion. Existing stabilization strategies often require higher-order finite elements (to compute second derivatives of displacement), introduce global auxiliary fields that increase system size, or rely on specialized discretizations. There is a need for a robust stabilization mechanism compatible with simple, computationally efficient first-order finite elements that avoids direct evaluation of second displacement derivatives.

Methodology
The authors propose a mixed finite element formulation that introduces an auxiliary-field stabilization for the third medium. The core of the method involves:

Auxiliary Field Introduction: A deformation-gradient-like field, denoted as $\Theta$ , is introduced in the third medium and treated as an independent unknown, decoupled from the physical displacement field $u$ .
Weak Coupling: The auxiliary field $\Theta$ is weakly coupled to the physical deformation gradient $F$ via a penalty term ( $W^p_{tm} \propto |\Theta - F|^2$ ). This ensures $\Theta$ tracks the physical deformation.
Gradient Regularization: Instead of penalizing the gradient of the physical deformation gradient ( $\nabla F$ ), which would require second derivatives of displacement, the formulation penalizes the gradient of the auxiliary field ( $\nabla \Theta$ ). The regularization energy is defined as $W^r_{tm} \propto |\nabla \Theta|^2$ .
Mixed Discretization: The formulation allows for independent interpolation of $u$ $u$ and $\Theta$ $Θ$ . The authors investigate both continuous and element-wise discontinuous approximations for $\Theta$ $Θ$ .
- Continuous $\Theta$ : Nodal values are shared between elements, creating inter-element coupling.
- Discontinuous $\Theta$ : Degrees of freedom are local to each element. These can be eliminated via static condensation, ensuring the global system retains only displacement degrees of freedom.

Key Contributions

Low-Order Compatibility: The formulation achieves gradient-type regularization without computing second derivatives of the displacement field, making it compatible with standard first-order (linear) and second-order finite elements.
Mechanism of Stabilization: The paper elucidates that for low-order displacement elements (where $F$ is element-wise constant), a continuous interpolation of $\Theta$ is essential. It transfers element-wise jumps in $F$ into non-vanishing gradients of $\Theta$ , thereby activating the regularization. In contrast, discontinuous $\Theta$ fields cannot generate these inter-element gradients, weakening the stabilization effect for low-order elements.
Static Condensation: For discontinuous auxiliary fields, the method demonstrates that the additional unknowns can be condensed out locally, preventing an increase in the global system size while retaining local stabilization.

Numerical Results
The formulation was assessed through several benchmark problems:

Two-Block Compression: Demonstrated that the third-medium scaling parameter $\gamma$ controls the trade-off between residual gap (contact accuracy) and nonlinear convergence robustness. Smaller $\gamma$ yields tighter contact but requires more Newton iterations. The method remained stable under combined compression and shear.
Self-Contact within a Box:
- Interpolation Comparison: Continuous auxiliary fields ( $T_u^1 T_\Theta^1$ , $T_u^2 T_\Theta^2$ ) provided robust contact closure. The mixed pair $T_u^2 T_\Theta^1$ (quadratic displacement, linear auxiliary) offered results comparable to $T_u^2 T_\Theta^2$ with fewer auxiliary degrees of freedom.
- Discontinuous vs. Continuous: Discontinuous auxiliary fields showed reduced regularization effectiveness, particularly for linear elements where the gradient term vanished or failed to couple elements.
- Comparison with Deformation-Gradient Averaging: Compared to the method by Faltus et al. [18], the proposed mixed formulation required significantly fewer load steps and Newton iterations, demonstrating superior nonlinear robustness. While the mixed formulation had a larger algebraic system, the improved convergence compensated for the cost in moderate mesh sizes, though the cost of auxiliary DOFs became dominant at extremely fine meshes ( $>10^6$ unknowns).
C-Shaped Box: The formulation successfully captured large rotations, progressive self-contact, and severe third-medium compression. The force-displacement response closely matched classical node-to-segment contact formulations, with smooth transitions into contact.
3D Extension: A three-dimensional self-contact example confirmed the method's applicability beyond plane strain settings.

Significance and Claims
The paper claims that the proposed auxiliary-field stabilization provides a robust, low-order-compatible alternative to classical contact algorithms and existing TMC regularizations. By avoiding direct evaluation of second displacement derivatives, it enables the use of simple, computationally cheap finite elements (linear and quadratic) for complex contact problems involving large deformations and self-contact.

The authors emphasize that the method's effectiveness is intrinsically linked to the continuity of the auxiliary field; continuous interpolations are necessary to generate the gradient-based stabilization mechanism for low-order displacement elements. The formulation is particularly suited for applications where contact boundaries are unknown or evolve, such as in contact-aided mechanisms, soft robotics, and topology optimization, where the continuum-based approach eliminates the need for explicit contact search algorithms and active set strategies. The work concludes that while the method trades a larger local algebraic system for improved global robustness, it offers a viable path for simulating severe contact scenarios without the complexity of higher-order elements or specialized mesh technologies.

Technical Summary: A Robust Mixed Finite Element Formulation for Third Medium Contact

More like this