Spherical harmonics-based pseudo-spectral method for quantitative analysis of symmetry breaking in wrinkling of shells with soft cores

Zavodnik, Jan (author), Brojan, Miha (author)

A complete understanding of the wrinkling of compressed films on curved substrates remains illusive due to the limitations of both analytical and current numerical methods. The difficulties arise from the fact that the energetically minimal distribution of deformation localizations is primarily influenced by the inherent nonlinearities and that the deformation patterns on curved surfaces are additionally constrained by the topology. The combination of two factors – the need for dense meshes to mitigate the topological limitations of discretization in domains such as spheres where there is no spherically-symmetric discretizations, and the intensive search for minima in a highly non-convex energy landscape due to nonlinearity – makes existing numerical methods computationally impractical without oversimplifying assumptions to reduce computational costs or introducing artificial parameters to ensure numerical stability. To solve these issues, we have developed a novel (less) reduced version of shell theory for shells subjected to membrane loads, such as during wrinkling. It incorporates the linear contributions of the usually excluded tangential displacements in the membrane strain energy and thus retains the computational efficiency of reduced state-of-the-art methods while nearly achieving the accuracy of the full Kirchhoff–Love shell theory. We introduce a Galerkin-type pseudo-spectral method to further reduce computational costs, prevent non-physical deformation distribution due to mesh-induced nucleation points, and avoid singularities at the poles of the sphere. The method uses spherical harmonic functions to represent functions on the surface of a sphere and is integrated into the framework of minimizing the total potential energy subject to constraints. This robust approach effectively solves the resulting non-convex potential energy problem. Our method accurately predicts the transition between deformation modes based solely on the material and geometric parameters determined in our experiments, without the need to introduce artificial parameters for numerical stability and/or additional fitting of the experimental data.

• Language: English
• Keywords: spectral method, spherical harmonics, wrinkling, Kirchhoff–love shell, constrained optimization
• Work type: Article
• Typology: 1.01 - Original Scientific Article
• Organization: FS - Faculty of Mechanical Engineering
• Publication status: Published
• Publication version: Version of Record
• Year: 2025
• Number of pages: 19 str.
• Numbering: Vol. 433, pt. B, [art. no.] 117529
• PID: 20.500.12556/RUL-164865
• UDC: 539.3
• ISSN on article: 0045-7825
• DOI: 10.1016/j.cma.2024.117529
• COBISS.SI-ID: 214921987
• Publication date in RUL: 14.11.2024

Record is a part of a journal

• Title: Computer methods in applied mechanics and engineering
• Shortened title: Comput. methods appl. mech. eng.
• Publisher: Elsevier
• ISSN: 0045-7825
• COBISS.SI-ID: 6695685

Licences

• License: CC BY 4.0, Creative Commons Attribution 4.0 International
• Link: http://creativecommons.org/licenses/by/4.0/
• Description: This is the standard Creative Commons license that gives others maximum freedom to do what they want with the work as long as they credit the author.

Secondary language

• Language: Slovenian
• Keywords: sferične harmonike, gubanje, teorija lupin

Projects

• Funder: ARIS - Slovenian Research and Innovation Agency
• Project number: J2-9223

• Funder: ARIS - Slovenian Research and Innovation Agency
• Project number: 2-2499