Direct numerical simulation of buckling instability of thin films on a compliant substrate

Siavash Nikravesh,Donghyeon Ryu,Yu-Lin Shen

doi:10.1177/1687814019840470

Siavash Nikravesh, Donghyeon Ryu + Show 1 more

Open Access

https://doi.org/10.1177/1687814019840470

Copy DOI

Abstract

Surface buckling (wrinkling) driven by mechanical instability is commonly observed in thin-film structures with a compliant substrate. The resulting undulation, while sometimes undesirable, has bee...

Highlights

IntroductionFor the post-buckling step, the incompressible neo-Hookean and Arruda–Boyce constitutive models for elastomers were exploited that require the use of incompressible/hybrid elements
The simulation results are presented, with emphasis given to the wrinkling wavelength, amplitude, and critical strain
The effect of imperfection placement is discussed in section ‘‘Effect of imperfection distribution.’’ Mesh convergence analyses are presented in section ‘‘Convergence analysis and comparisons with theory,’’ along with the comparisons between simulation results and analytical solutions described in section ‘‘Brief overview of theory.’’ Section ‘‘Further discussion’’ provides additional discussion of the numerical technique

Summary

Introduction

For the post-buckling step, the incompressible neo-Hookean and Arruda–Boyce constitutive models for elastomers were exploited that require the use of incompressible/hybrid elements Applying such elements may have some benefits within the reported numerical frameworks as to simulate the higher-order deformation modes at large strains, incompressible/hybrid elements are available only in limited finite element packages such as ABAQUS.[35] These types of elements are based on relatively complex mathematical theories, and are computationally more expensive compared with regular elements, significantly lowering the efficiency for three-dimensional analysis. The present work aims at developing a simple computational modeling approach to simulate temporal evolution of buckling Having this capability enables increased accuracy and versatility in predicting wrinkle formation of a relatively stiff thin film on a thick and compliant substrate. The critical buckling stress and strain can be determined, using

À n2fÁ

Results and discussion

Conclusion