Non-smooth optimization meets automated material model discovery
Flaschel M, Hastie T, Kuhl E (2026)
Publication Type: Journal article
Publication year: 2026
Journal
Book Volume: 563
Article Number: 115043
DOI: 10.1016/j.jcp.2026.115043
Abstract
Automated material model discovery has gained significant traction in recent years, as it disrupts the tedious and time-consuming cycle of iteratively calibrating and modifying manually designed models. Non-smooth L 1-norm regularization is the backbone of automated model discovery; however, the current literature on automated material model discovery offers limited insights into the robust and efficient minimization of non-smooth objective functions. In this work, we examine the minimization of functions of the form f(w)+α∥w∥1, where w are the material model parameters, f is a metric that quantifies the mismatch between the material model and the observed data, and α ≥ 0 is a regularization parameter that determines the sparsity of the solution. We investigate both the straightforward case where f is quadratic and the more complex scenario where it is non-quadratic or even non-convex. Importantly, in contrast to previous works on automated material model discovery, we do not only focus on methods that solve the sparse regression problem for a given value of the regularization parameter α , but propose methods to efficiently compute the entire regularization path, facilitating the selection of a suitable α . Specifically, we present four algorithms and discuss their roles for automated material model discovery in mechanics: First, we recapitulate a well-known coordinate descent algorithm that solves the minimization problem assuming that f is quadratic for a given value of α , also known as the LASSO. Second, we discuss the algorithm LARS, which automatically determines the critical values of α , at which material parameters in w are set to zero. Third, we propose to use the proximal gradient method ISTA for automated material model discovery if f is not quadratic, and fourth, we suggest a novel pathwise extension of ISTA for computing the regularization path. Many of these algorithms have not yet been applied to problems in computational material mechanics. We demonstrate the applicability of all algorithms for the automated discovery of incompressible hyperelastic material models from uniaxial tension and simple shear data.
Authors with CRIS profile
Involved external institutions
How to cite
APA:
Flaschel, M., Hastie, T., & Kuhl, E. (2026). Non-smooth optimization meets automated material model discovery. Journal of Computational Physics, 563. https://doi.org/10.1016/j.jcp.2026.115043
MLA:
Flaschel, Moritz, Trevor Hastie, and Ellen Kuhl. "Non-smooth optimization meets automated material model discovery." Journal of Computational Physics 563 (2026).
BibTeX: Download
United States (USA) (US)