pith. sign in

arxiv: 2606.28036 · v1 · pith:6BNYV3KLnew · submitted 2026-06-26 · 🧮 math.NA · cs.NA

A robust mixed finite element formulation for third medium contact

M. Vorwerk , J. Schr\"oder , P. Wriggers This is my paper

Pith reviewed 2026-06-29 03:24 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords third medium contactmixed finite elementsauxiliary fieldgradient regularizationlarge deformationself-contactfinite element stabilization
0
0 comments X

The pith

Continuous low-order auxiliary fields stabilize third medium contact even with first-order displacement elements.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a mixed finite element formulation for third medium contact that replaces explicit contact constraints with a compliant fictitious medium. It introduces an independent auxiliary field resembling the deformation gradient, coupled to the physical deformation gradient by a penalty term while a gradient term on the auxiliary field supplies regularization. This mechanism avoids direct evaluation of higher-order displacement derivatives. Numerical benchmarks demonstrate that continuous low-order auxiliary approximations remain stable for large deformations, progressive self-contact, and severe compression even when the displacement field uses first-order elements.

Core claim

Treating a deformation-gradient-like field as an independent unknown in the third medium, coupled by penalty to the physical deformation gradient and regularized by a gradient contribution acting on the auxiliary field, yields stable third medium contact simulations that do not require higher displacement derivatives and remain effective with low-order displacement approximations.

What carries the argument

Auxiliary deformation-gradient-like field treated as independent unknown with penalty coupling to the physical deformation gradient and gradient regularization on the auxiliary field.

If this is right

  • Continuous low-order auxiliary fields suffice for stabilization when the displacement field is approximated by first-order finite elements.
  • Element-wise discontinuous auxiliary fields remain local to each element and can be eliminated by static condensation without increasing the global system size.
  • The formulation remains effective under progressive self-contact and severe third-medium compression.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The local condensation option for discontinuous fields may preserve computational efficiency when scaling to large three-dimensional problems.
  • Avoiding explicit higher derivatives could simplify integration of third medium contact into existing low-order finite element libraries.

Load-bearing premise

The penalty coupling of the auxiliary field to the physical deformation gradient together with the gradient regularization on the auxiliary field is sufficient to eliminate instabilities in the third medium across the tested benchmarks without introducing new numerical artifacts.

What would settle it

A large-deformation self-contact benchmark in which instabilities or new artifacts appear in the third medium despite the introduction of the continuous low-order auxiliary field and its penalty-gradient regularization.

Figures

Figures reproduced from arXiv: 2606.28036 by J. Schr\"oder, M. Vorwerk, P. Wriggers.

Figure 1
Figure 1. Figure 1: One-dimensional illustration of the auxiliary-field regularization for two linear finite ele￾ments. a) Prescribed displacement field, b) gradient of u, c) auxiliary field with a continuous interpola￾tion, and d) auxiliary field with a discontinuous in￾terpolation. This example illustrates why a low-order discretization can still generate a non￾vanishing gradient of the auxiliary field Θ, and why this mecha… view at source ↗
Figure 2
Figure 2. Figure 2: Reference configuration of the two-block contact benchmark. The upper and lower solid blocks are separated by a third medium layer, which transfers the contact forces during compression. Dirichlet boundary conditions are prescribed on the top and bottom boundaries. The lower boundary is fixed with u = 0, while the upper boundary is displaced in case 1 until u = [0, −60]T and in case 2 until u = [10, −60]T … view at source ↗
Figure 3
Figure 3. Figure 3: Displacement-gap curves for the two-block benchmark under vertical compression. The gap measure is defined by u o 2 − u u 2 , where u o 2 and u u 2 denote the vertical displacements of the marked reference nodes. Results are shown for different values of the third-medium scaling parameter γ given in a) as the full loading path and b) as a magnified perspective of the contact regime [PITH_FULL_IMAGE:figure… view at source ↗
Figure 4
Figure 4. Figure 4: Final deformed configurations of the two-block benchmark under vertical com￾pression for different values of the third-medium scaling parameter γ. The gray overlay indicates the reference configuration. The trend observed in [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Displacement-gap curves for the two-block benchmark under combined horizontal and vertical loading. The gap measure is defined by u o 2 − u u 2 , where u o 2 and u u 2 denote the vertical displacements of the marked reference nodes in [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Final deformed configurations of the two-block benchmark under combined hor￾izontal and vertical loading for different values of the third-medium scaling parameter γ. The gray overlay indicates the reference configuration. The final deformed configurations for case 2 are shown in [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Reference and deformed configurations of the self-contact-within-a-box bench￾mark for two penalty parameters and the Tu 1T Θ 1 approximation. The solid frame is shown in orange, while the enclosed third medium is shown in yellow. Graphics a) and c) use γ = 10−5 , αr = 102 and pΘ = 10−2 ; graphics b) and d) use γ = 10−5 , αr = 102 and pΘ = 10−1 . Panels a) and b) show the reference configurations, whereas c… view at source ↗
Figure 8
Figure 8. Figure 8: Reference configuration of the C-shaped box benchmark for finite deformation contact. Orange indicates the physical solid domain, and yellow indicates the third medium. Blue dots mark the clamped boundary, while the violet dot marks the point where the prescribed displacement is applied. Since the penalty coupling is mesh-dependent, the parameter pΘ is adjusted for each dis￾cretization rather than kept fix… view at source ↗
Figure 9
Figure 9. Figure 9: Deformation history of the C-shaped benchmark for a calibrated sequence of in￾creasingly refined meshes using the continuous Tu 1T Θ 1 discretization. The penalty parameter pΘ is adjusted for each mesh. Snapshots correspond to the load parameters λload = 0.4, λload = 0.7 and λload = 1.0 [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Deformation history of the C-shaped benchmark using the continuous Tu 1T Θ 1 dis￾cretization. Snapshots are shown for increasing load parameters λload. All deformations are plotted without displacement scaling. The gray overlay indicates the reference configuration. A more detailed view of the deformation process is given in [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Reaction-force-displacement response of the C-shaped benchmark. The third medium contact formulation is compared with a classical node-to-segment contact formula￾tion. a) Full loading path. b) Magnified view of the contact regime [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of third-medium deformations in the C-shaped benchmark at λload = 0.2. The left column shows the proposed formulation with continuous auxiliary￾field interpolations, the middle column shows the corresponding element-wise discontinuous variants, and the right column shows the deformation-gradient averaging approach of Fal￾tus et al. [18]. The upper row shows triangular discretizations using a) T… view at source ↗
Figure 13
Figure 13. Figure 13: shows the deformed configuration from two perspectives. During loading, the up￾per part of the frame bends into the box, the third medium is compressed between the approaching solid surfaces, and the deformation remains regular. Overall, the example confirms that the auxiliary-field stabilized third medium formulation carries over directly to three-dimensional finite deformation contact problems. a) b) [… view at source ↗
read the original abstract

Third medium contact provides a smooth continuum alternative to classical contact algorithms by replacing explicit contact constraints with a highly compliant fictitious medium. In this work, an auxiliary-field stabilization is introduced in which a deformation-gradient-like field is treated as an independent unknown in the third medium and coupled to the physical deformation gradient by a penalty term. A gradient contribution acting on the auxiliary field provides the regularization mechanism without requiring a direct evaluation of higher displacement derivatives. Linear and quadratic interpolation spaces are investigated, including continuous and element-wise discontinuous auxiliary-field approximations. The numerical results show that continuous low-order auxiliary fields provide an effective gradient-type stabilization of the third medium, even when the displacement field is approximated by first-order finite elements. For element-wise discontinuous auxiliary fields, the additional unknowns remain local to each element and can be eliminated locally by static condensation, so that the global system does not necessarily contain additional auxiliary degrees of freedom. Benchmark problems involving large deformation, progressive self-contact and severe third-medium compression are used to assess the formulation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper introduces a mixed finite-element formulation for third-medium contact in which a deformation-gradient-like auxiliary field is introduced as an independent unknown within the fictitious medium and coupled to the physical deformation gradient by a penalty term; a gradient regularization is applied directly to the auxiliary field. Continuous and element-wise discontinuous approximations are considered for both linear and quadratic spaces. Numerical experiments on large-deformation, progressive self-contact and severe-compression benchmarks indicate that continuous low-order auxiliary fields stabilize the third medium effectively even when the displacement field is discretized with first-order elements; discontinuous auxiliary fields can be eliminated by static condensation.

Significance. If the reported stabilization holds under the stated penalty and regularization choices, the formulation supplies a practical continuum-based route to contact that avoids explicit constraint enforcement and higher-order displacement derivatives. The observation that low-order continuous auxiliary fields suffice with first-order displacements is computationally attractive and could simplify implementation in existing low-order codes. The local condensation property for discontinuous fields is a further implementation advantage.

minor comments (3)
  1. The abstract states that continuous low-order auxiliary fields 'provide an effective gradient-type stabilization' but does not indicate the range of penalty parameters over which this holds or whether any post-hoc tuning was required; a brief statement in §4 or §5 on parameter sensitivity would strengthen the claim.
  2. The description of the auxiliary-field interpolation spaces ('linear and quadratic') is not accompanied by the corresponding polynomial degrees for the displacement field in the same sentence; clarifying the exact combination tested in each benchmark (e.g., P1 displacement + P1 continuous auxiliary) would improve readability.
  3. No mention is made of the precise form of the gradient regularization term (e.g., whether it is the full gradient of the auxiliary field or a projected version); adding the weak-form expression in §3 would eliminate ambiguity.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of our work on auxiliary-field stabilization in mixed finite elements for third medium contact and for recommending minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The provided abstract and context present the auxiliary-field stabilization as a direct extension of standard mixed finite-element methods, with the central claim (effectiveness of continuous low-order auxiliary fields for gradient-type regularization) resting on the formulation's penalty coupling and gradient term rather than any self-definitional reduction, fitted prediction, or load-bearing self-citation. No equations, ansatzes, or uniqueness theorems are quoted that collapse to the paper's own inputs by construction. The derivation chain is therefore self-contained against external benchmarks of mixed FE theory.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The formulation rests on standard continuum mechanics and finite-element discretization assumptions plus the new auxiliary field and its penalty coupling; no independent evidence for the auxiliary field's effectiveness outside the paper is supplied in the abstract.

free parameters (1)
  • penalty parameter for auxiliary-field coupling
    The strength of the penalty term that couples the auxiliary field to the physical deformation gradient is a tunable parameter whose value must be chosen to achieve stabilization.
axioms (1)
  • standard math Standard assumptions of finite element discretization, continuum mechanics, and penalty regularization hold for the mixed formulation.
    The approach builds directly on established mixed finite-element theory without stating new foundational lemmas.
invented entities (1)
  • auxiliary deformation-gradient-like field no independent evidence
    purpose: Serves as an independent unknown providing gradient regularization without explicit higher-order derivatives of displacement.
    New field introduced in the stabilization scheme; no external falsifiable prediction is given in the abstract.

pith-pipeline@v0.9.1-grok · 5704 in / 1372 out tokens · 48431 ms · 2026-06-29T03:24:40.639226+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 26 canonical work pages

  1. [1]

    Wriggers.Computational Contact Mechanics

    P. Wriggers.Computational Contact Mechanics. Springer, Berlin, Heidelberg, 2 edition, 2006. doi: 10.1007/978-3-540-32609-0

    work page doi:10.1007/978-3-540-32609-0 2006

  2. [2]

    Laursen.Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfa- cial Phenomena in Nonlinear Finite Element Analysis

    T.A. Laursen.Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfa- cial Phenomena in Nonlinear Finite Element Analysis. Springer, Berlin, Heidelberg, 2002. doi: 10.1007/978-3-662-04864-1

    work page doi:10.1007/978-3-662-04864-1 2002

  3. [3]

    Yastrebov.Numerical Methods in Contact Mechanics

    V.A. Yastrebov.Numerical Methods in Contact Mechanics. Wiley, London, 2013. ISBN 9781848215191. doi: 10.1002/9781118647974

    work page doi:10.1002/9781118647974 2013

  4. [4]

    Rus and M.T

    D. Rus and M.T. Tolley. Design, fabrication and control of soft robots.Nature, 521(7553):467–475,

  5. [5]

    doi: 10.1038/nature14543. M. Vorwerk, J. Schr¨ oder, and P. Wriggers23

    work page doi:10.1038/nature14543

  6. [6]

    Frederiksen, O

    A.H. Frederiksen, O. Sigmund, and K. Poulios. Topology optimization of self-contacting structures. Computational Mechanics, 73:967–981, 2024. doi: 10.1007/s00466-023-02396-7

    work page doi:10.1007/s00466-023-02396-7 2024

  7. [7]

    doi:10.1016/j.cma.2003.10.010 Sharpa Robotics

    M.A. Puso and T.A. Laursen. A mortar segment-to-segment contact method for large deformation solid mechanics.Computer Methods in Applied Mechanics and Engineering, 193(6–8):601–629, 2004. doi: 10.1016/j.cma.2003.10.010

    work page doi:10.1016/j.cma.2003.10.010 2004

  8. [8]

    Sauer and L

    R.A. Sauer and L. De Lorenzis. A computational contact formulation based on surface po- tentials.Computer Methods in Applied Mechanics and Engineering, 253:369–395, 2013. doi: 10.1016/j.cma.2012.09.002

    work page doi:10.1016/j.cma.2012.09.002 2013

  9. [9]

    Bluhm, O

    G.L. Bluhm, O. Sigmund, and K. Poulios. Internal contact modeling for finite strain topology optimization.Computational Mechanics, 67:1099–1114, 2021. doi: 10.1007/s00466-021-01974-x

    work page doi:10.1007/s00466-021-01974-x 2021

  10. [10]

    Wriggers, J

    P. Wriggers, J. Korelc, and P. Junker. A third medium approach for contact using first and second order finite elements.Computer Methods in Applied Mechanics and Engineering, 436:117740, 2025. doi: 10.1016/j.cma.2025.117740

    work page doi:10.1016/j.cma.2025.117740 2025

  11. [11]

    Wriggers, J

    P. Wriggers, J. Schr¨ oder, and A. Schwarz. A finite element method for contact using a third medium. Computational Mechanics, 52:837–847, 2013. doi: 10.1007/s00466-013-0848-5

    work page doi:10.1007/s00466-013-0848-5 2013

  12. [12]

    T. Bog, N. Zander, S. Kollmannsberger, and E. Rank. Normal contact with high order finite elements and a fictitious contact material.Computers and Mathematics with Applications, 70:1370–1390, 2015. doi: 10.1016/j.camwa.2015.04.020

    work page doi:10.1016/j.camwa.2015.04.020 2015

  13. [13]

    Kruse, N

    R. Kruse, N. Nguyen-Thanh, P. Wriggers, and L. De Lorenzis. Isogeometric frictionless contact analysis with the third medium method.Computational Mechanics, 62:1009–1021, 2018. doi: 10.1007/s00466-018-1547-z

    work page doi:10.1007/s00466-018-1547-z 2018

  14. [14]

    Huang, N

    J. Huang, N. Nguyen-Thanh, and K. Zhou. An isogeometric-meshfree coupling approach for contact problems by using the third medium method.International Journal of Mechanical Sciences, 148: 327–336, 2018. doi: 10.1016/j.ijmecsci.2018.08.031

    work page doi:10.1016/j.ijmecsci.2018.08.031 2018

  15. [15]

    G. L. Bluhm, O. Sigmund, and K. Poulios. Inverse design of mechanical springs with tailored nonlinear elastic response utilizing internal contact.International Journal of Non-Linear Mechanics, 157:104552, 2023. doi: 10.1016/j.ijnonlinmec.2023.104552

    work page doi:10.1016/j.ijnonlinmec.2023.104552 2023

  16. [16]

    Dalklint, J

    A. Dalklint, J. Alexandersen, A.H. Frederiksen, K. Poulios, and O. Sigmund. Topology optimization of contact-aided thermo-mechanical regulators.International Journal for Numerical Methods in Engineering, 126(2):e7661, 2025. doi: 10.1002/nme.7661

    work page doi:10.1002/nme.7661 2025

  17. [17]

    Frederiksen, O

    A.H. Frederiksen, O. Rokoˇ s, K. Poulios, O. Sigmund, and M.G.D. Geers. Adding friction to third medium contact: A crystal plasticity inspired approach.Computer Methods in Applied Mechanics and Engineering, 432:117412, 2024. doi: 10.1016/j.cma.2024.117412

    work page doi:10.1016/j.cma.2024.117412 2024

  18. [18]

    Frederiksen, A

    A.H. Frederiksen, A. Dalklint, O. Sigmund, and K. Poulios. Improved third medium formulation for 3D topology optimization with contact.Computer Methods in Applied Mechanics and Engineering, 436:117595, 2025. doi: 10.1016/j.cma.2024.117595

    work page doi:10.1016/j.cma.2024.117595 2025

  19. [19]

    Faltus, M

    O. Faltus, M. Amato, and M. Hor´ ak. Deformation gradient averaging regularization for third medium contact.Computer Methods in Applied Mechanics and Engineering, 458:119072, 2026. doi: 10.1016/j.cma.2026.119072

    work page doi:10.1016/j.cma.2026.119072 2026

  20. [20]

    Dahlberg, F

    V. Dahlberg, F. Sj¨ ovall, A. Dalklint, and M. Wallin. A rotation-based approach to third medium contact regularization.Computer Methods in Applied Mechanics and Engineering, 453:118801, 2026. doi: 10.1016/j.cma.2026.118801

    work page doi:10.1016/j.cma.2026.118801 2026

  21. [21]

    In: Bue, A.D., Canton, C., Pont-Tuset, J., Tommasi, T

    P. Wriggers, J. Korelc, and BB. Xu. Third medium contact: An overview. In Y. Bazilevs and K. Takizawa, editors,Advances in Computational Fluid-Structure Interaction and Flow Simulation. Modeling and Simulation in Science, Engineering and Technology, Springer, 2026. doi: 10.1007/978- 3-032-25894-6 15. A robust mixed finite element formulation for third med...

    work page doi:10.1007/978- 2026

  22. [22]

    Wriggers

    P. Wriggers. A third medium approach for thermo-mechanical contact based on low order ansatz spaces.Finite Elements in Analysis and Design, 255:104522, 2026. doi: 10.1016/j.finel.2026.104522

    work page doi:10.1016/j.finel.2026.104522 2026

  23. [23]

    Xu and P

    B.-B. Xu and P. Wriggers. Stabilization-free virtual element method for 2D third medium contact.Computer Methods in Applied Mechanics and Engineering, 450:118611, 2026. doi: 10.1016/j.cma.2025.118611

    work page doi:10.1016/j.cma.2025.118611 2026

  24. [24]

    B.-B. Xu, T. Xue, and P. Wriggers. Three-dimensional third medium contact model for hyperelastic contact and pneumatically actuated systems.Journal of the Mechanics and Physics of Solids, 213: 106617, 2026. doi: 10.1016/j.jmps.2026.106617

    work page doi:10.1016/j.jmps.2026.106617 2026

  25. [25]

    von Zabiensky, D.R

    M. von Zabiensky, D.R. Jantos, and P. Junker. A fast and robust third medium contact approach using the neighbored element method.Finite Elements in Analysis and Design, 255:104489, 2026. doi: 10.1016/j.finel.2025.104489

    work page doi:10.1016/j.finel.2025.104489 2026

  26. [26]

    Faltus, M

    O. Faltus, M. Hor´ ak, M. Doˇ sk´ aˇ r, and O. Rokoˇ s. Third medium finite element contact formulation for pneumatically actuated systems.Computer Methods in Applied Mechanics and Engineering, 431: 117262, 2024. doi: 10.1016/j.cma.2024.117262

    work page doi:10.1016/j.cma.2024.117262 2024

  27. [27]

    Wriggers, J

    P. Wriggers, J. Korelc, and P. Junker. First order finite element formulations for third medium contact.Computational Mechanics, 76:829–845, 2025. doi: 10.1007/s00466-025-02628-y. M. Vorwerk, J. Schr¨ oder, and P. Wriggers25 A Appendix T able 4:Residual gap, total Newton iteration count and achieved load parameter for the self-contact-within-a-box benchma...

    work page doi:10.1007/s00466-025-02628-y 2025