A differentiable variational model for structural self-contact and fracture
File(s)
Published version
Author(s)
Ciceri, Mirko
Aveline, Charlie
Thillaithevan, Dilaksan
Hewson, Robert
Santer, Matthew
Type
Journal Article
Abstract
Numerical modelling of structural self-contact and crack propagation presents significant challenges due to the inherently discontinuous and non-differentiable
nature of the underlying physical phenomena. Traditional contact models demand explicit definition and tracking of contact points, while fracture models often rely on predefined crack initiation sites, sharp interfaces, and re-meshing. This study introduces a novel framework that overcomes these limitations within a unified and numerically stable variational formulation. The contact phenomenon is described through the hyperelastic third medium contact model and fracture is represented by a phase field. Structures are embedded in a third
medium, that stiffens under compression, enabling the transfer of forces between structural members. Crack propagation occurs in regions in which it is energetically favourable for the system to evolve toward a fully damaged state, specifically where the critical energy release rate is exceeded. Careful treatment is required when coupling the two phenomena, particularly concerning the void material behaviour. This work presents an efficient and differentiable numerical model that captures both nonlinear phenomena within a unified framework. This framework will allow designers and engineers to efficiently analyse complex nonlinear structural behaviours, previously requiring separate models that involved pre-defined crack initiation sites and contact points. Lastly, the differentiable nature of the model facilitates straightforward future integration into
topology optimisation pipelines, providing designers the ability to intentionally design for and leverage self-contact interactions and material failure as functional,
performance-enhancing features.
nature of the underlying physical phenomena. Traditional contact models demand explicit definition and tracking of contact points, while fracture models often rely on predefined crack initiation sites, sharp interfaces, and re-meshing. This study introduces a novel framework that overcomes these limitations within a unified and numerically stable variational formulation. The contact phenomenon is described through the hyperelastic third medium contact model and fracture is represented by a phase field. Structures are embedded in a third
medium, that stiffens under compression, enabling the transfer of forces between structural members. Crack propagation occurs in regions in which it is energetically favourable for the system to evolve toward a fully damaged state, specifically where the critical energy release rate is exceeded. Careful treatment is required when coupling the two phenomena, particularly concerning the void material behaviour. This work presents an efficient and differentiable numerical model that captures both nonlinear phenomena within a unified framework. This framework will allow designers and engineers to efficiently analyse complex nonlinear structural behaviours, previously requiring separate models that involved pre-defined crack initiation sites and contact points. Lastly, the differentiable nature of the model facilitates straightforward future integration into
topology optimisation pipelines, providing designers the ability to intentionally design for and leverage self-contact interactions and material failure as functional,
performance-enhancing features.
Date Issued
2026-02-19
Date Acceptance
2026-01-28
Citation
Engineering with Computers, 2026, 42
ISSN
0177-0667
Publisher
Springer
Journal / Book Title
Engineering with Computers
Volume
42
Copyright Statement
© The Author(s) 2026. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
License URL
Identifier
10.1007/s00366-026-02285-6
Subjects
Phase field
Fracture
Third medium contact
Publication Status
Published
Article Number
ARTN 52