pith. sign in

arxiv: 1907.00447 · v1 · pith:WKDWT45Lnew · submitted 2019-06-30 · 🧮 math.AP · cond-mat.mtrl-sci

Energy minimising configurations of pre-strained multilayers

Miguel de Benito Delgado , Bernd Schmidt This is my paper

Pith reviewed 2026-05-25 12:20 UTC · model grok-4.3

classification 🧮 math.AP cond-mat.mtrl-sci
keywords pre-strainmultilayersvon Karman theoryenergy minimizationphase transitionspontaneous curvatureplate theoriesasymptotics
0
0 comments X

The pith

A family of interpolated von Kármán functionals for pre-strained plates converges to cylindrical minimisers for small θ and spherical ones for large θ.

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

The authors construct a continuous family of effective von Kármán energies for thin pre-strained multilayers that interpolates between the linearised Kirchhoff and linearised von Kármán regimes. They prove rigorous convergence of the minimisers to the known explicit forms in the two asymptotic limits of the interpolation parameter. Numerical computations for intermediate values reveal a sharp transition between cylindrical and spherical configurations at a critical finite value of the parameter. This provides a mathematical description of how the strength of pre-strain selects the optimal shape in the critical nonlinear regime.

Core claim

The effective von Kármán functionals I^θ_vK admit explicit minimisers that are cylindrical when θ is small and spherical when θ is large, with a phase transition occurring at some critical θ in between, as confirmed by analysis in the limits and by numerics in the interior.

What carries the argument

The one-parameter family of von Kármán functionals (I^θ_vK) that interpolates between the linearised plate theories with spontaneous curvature induced by the pre-strain.

If this is right

  • The minimisers converge to the explicit cylindrical solutions of the linearised Kirchhoff theory as θ approaches 0.
  • The minimisers converge to the explicit spherical solutions of the linearised von Kármán theory as θ approaches infinity.
  • Numerical minimisation indicates the transition occurs abruptly at a specific critical value of θ.
  • The hierarchy of plate theories depends on the strength of the pre-strain in the multilayers.

Where Pith is reading between the lines

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

  • The critical θ may mark a point where the competing effects of bending and stretching balance in a new way for pre-strained structures.
  • Similar interpolation techniques could be used to study transitions in other asymptotic regimes of thin structure theories.
  • Experimental realisation of pre-strained multilayers with tunable parameters could test the predicted shape transition.

Load-bearing premise

The pre-strain must be such that the three-dimensional energy reduces to one of the effective two-dimensional plate models with spontaneous curvature in the hierarchy.

What would settle it

Failure to observe convergence to the predicted explicit minimisers when numerically minimising the interpolated functional for very small or very large θ, or absence of a transition in the shape of minimisers as θ varies.

Figures

Figures reproduced from arXiv: 1907.00447 by Bernd Schmidt, Miguel de Benito Delgado.

Figure 1
Figure 1. Figure 1: The partition of ω 0 into bodies and arms. ∇y is constant in the bodies (colored) and along each of the straight lines making up the arms (white). The existence of covered domains for isometric immersions y ∈ W1,2 is shown in [28, Corollary 1.2]. Proposition 1 Let v ∈ W 2,2 sh (ω) and x0 ∈ ω. There exists a neighbourhood U of x0 such that, if ∇2v 6= 0 a.e. in U, then for a suitable ε > 0 there exist maps γ… view at source ↗
Figure 2
Figure 2. Figure 2: Final configurations after gradient descent starting with a flat disk [PITH_FULL_IMAGE:figures/full_fig_p028_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Mean principal strains (left) and symmetry (right) of the minimiser [PITH_FULL_IMAGE:figures/full_fig_p029_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Initial (left) and final (right) states starting with skewed paraboloid. [PITH_FULL_IMAGE:figures/full_fig_p029_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Mean principal strains (left) and symmetry (right) of the minimiser [PITH_FULL_IMAGE:figures/full_fig_p029_5.png] view at source ↗
read the original abstract

We investigate energetically optimal configurations of thin structures with a pre-strain. Depending on the strength of the pre-strain we consider a whole hierarchy of effective plate theories with a spontaneous curvature term, ranging from linearised Kirchhoff to von K\'arm\'an to linearised von K\'arm\'an theories. While explicit formulae are available in the linearised regimes, the von K\'arm\'an theory turns out to be critical and a phase transition from cylindrical (as in linearised Kirchhoff) to spherical (as in von linearised K\'arm\'an) configurations is observed there. We analyse this behavior with the help of a whole family $(\mathcal{I}^{\theta}_{\rm vK})_{\theta \in (0,\infty)}$ of effective von K\'arm\'an functionals which interpolates between the two linearised regimes. We rigorously show convergence to the respective explicit minimisers in the asymptotic regimes $\theta \to 0$ and $\theta \to \infty$. Numerical experiments are performed for general $\theta \in (0,\infty)$ which indicate a stark transition at a critical value of $\theta$.

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 manuscript derives a hierarchy of effective 2D plate models with spontaneous curvature for pre-strained thin multilayers, obtained from 3D nonlinear elasticity depending on pre-strain strength. Explicit minimizers are identified in the linearised Kirchhoff and linearised von Kármán regimes. A one-parameter family of von Kármán functionals (I^θ_vK) is introduced that interpolates between these regimes. Rigorous convergence of minimizers to the explicit ones is established as θ → 0 and θ → ∞. Numerical experiments for intermediate θ indicate a phase transition from cylindrical to spherical configurations at a critical value of θ.

Significance. If the convergence theorems hold, the work supplies a rigorous variational justification for the effective models across regimes and clarifies the critical role of the von Kármán scaling via the interpolating family. The explicit minimizers and the observed transition are of direct interest to the mechanics of morphing and pre-strained structures. The combination of Gamma-convergence-style analysis with targeted numerics strengthens the contribution to the literature on dimension reduction with pre-strain.

minor comments (3)
  1. [§4] §4 (or the section containing the definition of I^θ_vK): the precise dependence of the spontaneous curvature term on the pre-strain parameters should be stated explicitly, including any scaling assumptions that distinguish the linearised Kirchhoff from the linearised von Kármán limit.
  2. [Numerical experiments] Numerical section: the discretization scheme (finite elements, mesh refinement strategy) and the precise functional form used for the numerical minimization of I^θ_vK are described only briefly; adding these details would improve reproducibility.
  3. [§2] The statement that the pre-strain 'permits reduction to the stated hierarchy' (abstract and §2) would benefit from a short remark on the admissible class of pre-strain tensors that allow the spontaneous curvature to appear at the indicated orders.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the positive assessment, including the recommendation for minor revision. No major comments are listed in the report, so there are no specific points requiring point-by-point rebuttal or changes to the analysis. We are happy to address any minor editorial matters in the revised version.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper performs a rigorous variational analysis of pre-strained multilayers, deriving a hierarchy of effective plate models (linearised Kirchhoff, von Kármán, linearised von Kármán) from 3D elasticity and proving Gamma-convergence of an interpolating family of functionals to explicit minimisers as θ → 0 and θ → ∞. No load-bearing step reduces by definition or construction to its own inputs, no parameters are fitted and then relabelled as predictions, and no uniqueness claims rest on self-citations. The derivation is self-contained mathematical analysis whose central results are externally falsifiable via the stated convergence statements and numerical observations.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on abstract; main assumptions concern the form of pre-strain allowing dimension reduction to plate models with spontaneous curvature and the validity of the chosen scaling regimes.

axioms (1)
  • domain assumption Pre-strain of appropriate strength permits reduction to effective plate theories with spontaneous curvature term
    Invoked to justify the hierarchy ranging from linearised Kirchhoff to von Kármán.

pith-pipeline@v0.9.0 · 5724 in / 1137 out tokens · 39000 ms · 2026-05-25T12:20:11.276154+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    M. S. Alnaes, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. Richardson, J. Ring, M. E. Rognes, and G. N. Wells. The FEniCS Project Version 1.5. Archive of Numerical Software, 3(100), 2015. 27

    work page 2015

  2. [2]

    Babuˇ ska and J

    I. Babuˇ ska and J. Pitk¨ aranta. The plate paradox for hard and soft simple support. SIAM Journal on Mathematical Analysis , 21(3):551–576, 1990. 22

    work page 1990

  3. [3]

    Babuˇ ska and M

    I. Babuˇ ska and M. Suri. On Locking and Robustness in the Finite Element Method. SIAM Journal on Numerical Analysis , 29(5):1261–1293, 1992. 22

    work page 1992

  4. [4]

    S. Bartels. Approximation of Large Bending Isometries with Discrete Kirch- hoff Triangles. SIAM Journal on Numerical Analysis , 51(1):516–525, 2013. 19

    work page 2013

  5. [5]

    S. Bartels. Numerical Methods for Nonlinear Partial Differential Equations, volume 47 of Springer Series in Computational Mathematics . Springer International Publishing, 2015. 20, 21, 26

    work page 2015

  6. [6]

    S. Bartels. Numerical solution of a F¨ oppl–von K´ arm´ an model.SIAM J. Numer. Anal., 55(3):1505–1524, 2017. 19

    work page 2017

  7. [7]

    Bartels, A

    S. Bartels, A. Bonito, and R. H. Nochetto. Bilayer plates: Model reduc- tion, Γ-convergent finite element approximation, and discrete gradient flow. Communications on Pure and Applied Mathematics , 70(3):547–589, 2017. 2

    work page 2017

  8. [8]

    S. C. Brenner, M. Neilan, A. Reiser, and L.-Y. Sung. A c0 interior penalty method for a von K´ arm´ an plate.Numerische Mathematik, 135(3):803–832,

    work page

  9. [9]

    S. C. Brenner and L. R. Scott. The Mathematical Theory of Finite Element Methods. Number 15 in Texts in Applied Mathematics. Springer New York, New York, NY, 3rd edition, 2008. 21, 23

    work page 2008

  10. [10]

    de Benito Delgado

    M. de Benito Delgado. Implementation of a nonlinear Kirchhoff plate model. https://bitbucket.org/mdbenito/nonlinear-kirchhoff, 2017. 27

    work page 2017

  11. [11]

    de Benito Delgado and B

    M. de Benito Delgado and B. Schmidt. A hierarchy of multilayered plate models. In preparation, 2019. 2, 3, 4, 5, 6, 8, 9, 12

    work page 2019

  12. [12]

    A. I. Egunov, J. G. Korvink, and V. A. Luchnikov. Polydimethylsilox- ane bilayer films with an embedded spontaneous curvature. Soft Matter , 12(1):45–52, 2016. 2 30

    work page 2016

  13. [13]

    Finot and S

    M. Finot and S. Suresh. Small and large deformation of thick and thin- film multi-layers: Effects of layer geometry, plasticity and compositional gradients. Journal of the Mechanics and Physics of Solids , 44(5):683 – 721,

    work page

  14. [14]

    Mechanics and Physics of Layered and Graded Materials. 2

    work page

  15. [15]

    L. B. Freund. Substrate curvature due to thin film mismatch strain in the nonlinear deformation range. J. Mech. Phys. Solids , 48(6-7):1159–1174,

    work page

  16. [16]

    The J. R. Willis 60th anniversary volume. 2

    work page

  17. [17]

    Friesecke, R

    G. Friesecke, R. D. James, and S. M¨ uller. A hierarchy of plate models derived from nonlinear elasticity by Γ-convergence. Archive for Rational Mechanics and Analysis, 180(2):183–236, 2006. 3, 10

    work page 2006

  18. [18]

    K. Greff, A. Klein, M. Chovanec, F. Hutter, and J. Schmidhuber. The Sacred Infrastructure for Computational Research. Proceedings of the 16th Python in Science Conference , pages 49–56, 2017. 27

    work page 2017

  19. [19]

    Grossmann, H.-G

    C. Grossmann, H.-G. Roos, and M. Stynes. Numerical Treatment of Partial Differential Equations . Universitext. Springer Berlin Heidelberg, Berlin, Heidelberg, 2007. 21, 22

    work page 2007

  20. [20]

    Grundmann

    M. Grundmann. Nanoscroll formation from strained layer heterostructures. Applied Physics Letters, 83:2444–2446, 2003. 2

    work page 2003

  21. [21]

    P. Hornung. Approximation of flat W 2,2 isometric immersions by smooth ones. Archive for Rational Mechanics and Analysis , 199(3):1015–1067,

    work page

  22. [22]

    Howell, G

    P. Howell, G. Kozyreff, and J. Ockendon. Applied Solid Mechanics. Num- ber 43 in Cambridge Texts in Applied Mathematics. Cambridge University Press, 2008. 19

    work page 2008

  23. [23]

    C. S. Kim and S. J. Lombardo. Curvature and bifurcation of mgo-al 2o3 bi- layer ceramic structures. Journal of Ceramic Processing Research, 9(2):93– 96, 2008. 2

    work page 2008

  24. [24]

    Kluyver, B

    T. Kluyver, B. Ragan-Kelley, F. P´ erez, B. Granger, M. Bussonnier, J. Fred- eric, K. Kelley, J. Hamrick, J. Grout, S. Corlay, P. Ivanov, D. Avila, S. Ab- dalla, and C. Willing. Jupyter Notebooks – a publishing format for repro- ducible computational workflows. In F. Loizides and B. Schmidt, editors, Positioning and Power in Academic Publishing: Players, A...

    work page 2016

  25. [25]

    Mallik and N

    G. Mallik and N. Nataraj. Conforming finite element methods for the von K´ arm´ an equations.Advances in Computational Mathematics , 42(5):1031– 1054, 2016. 19

    work page 2016

  26. [26]

    Mallik and N

    G. Mallik and N. Nataraj. A nonconforming finite element approxima- tion for the von K´ arm´ an equations.ESAIM: Mathematical Modelling and Numerical Analysis, 50(2):433–454, 2016. 19

    work page 2016

  27. [27]

    C. B. Masters and N. Salamon. Geometrically nonlinear stress-deflection relations for thin film/substrate systems. International Journal of Engi- neering Science, 31(6):915 – 925, 1993. 2 31

    work page 1993

  28. [28]

    C. Ortner. Γ-Limits of Galerkin Discretizations with Quadrature. Tech- nical Report 04/26, Oxford University Computing Laboratory, Numerical Analysis Group, 2004. 19

    work page 2004

  29. [29]

    Paetzelt, V

    H. Paetzelt, V. Gottschalch, J. Bauer, H. Herrnberger, and G. Wagner. Fabrication of iii–v nano- and microtubes using movpe grown materials. physica status solidi (a) , 203(5):817–824, 2006. 2

    work page 2006

  30. [30]

    M. R. Pakzad. On the Sobolev space of isometric immersions. Journal of Differential Geometry, 66(1):47–69, 2004. 9, 10

    work page 2004

  31. [31]

    G. Prathap. Finite element analysis as computation. 2001. 22

    work page 2001

  32. [32]

    Quaglino

    A. Quaglino. Membrane Locking in Discrete Shell Theories . Doctoral dis- sertation, Georg-August-Universit¨ at G¨ ottingen, G¨ ottingen, 2012. 22

    work page 2012

  33. [33]

    Salamon and C

    N. Salamon and C. B. Masters. Bifurcation in isotropic thinfilm/substrate plates. International Journal of Solids and Structures , 32(3):473 – 481,

    work page

  34. [34]

    Special topics in the theory of elastic: A volume in honour of Professor John Dundurs. 2

    work page

  35. [35]

    B. Schmidt. Minimal energy configurations of strained multi-layers. Calcu- lus of Variations and Partial Differential Equations , 30(4):477–497, 2007. 2

    work page 2007

  36. [36]

    B. Schmidt. Plate theory for stressed heterogeneous multilayers of finite bending energy. Journal de Math´ ematiques Pures et Appliqu´ ees, 88(1):107– 122, 2007. 2, 10

    work page 2007

  37. [37]

    O. G. Schmidt and K. Eberl. Thin solid films roll up into nanotubes. Nature, 410(9):168, 2001. 2

    work page 2001

  38. [38]

    V. R. Subramanian. Omniboard: a web-based dashboard for Sacred, 2018. 27 32

    work page 2018