Wrinkling of Randomly Heterogeneous Film-Substrate Systems

Abstract

Wrinkling instabilities in stiff films on compliant substrates are strongly affected by spatial fluctuations in film stiffness. We develop a homogenized instability theory for one-dimensional film--substrate systems with random bending stiffness. The heterogeneous stability equation is reformulated as a Lippmann--Schwinger equation for the curvature field, and a strong-contrast expansion is derived using a local--nonlocal kernel decomposition and a cavity-field formulation. Truncation at third order yields an effective polarizability and a Dyson-type dispersion relation for predicting the critical load and wavenumber. The stiffness is modeled as an exponentially mapped Gaussian random field, allowing the required two- and three-point connected statistics to be obtained analytically. The theory is validated against generalized eigenvalue calculations and Fourier spectral simulations. Increasing stiffness contrast lowers the critical load and shifts the instability toward higher wavenumbers, producing shorter wrinkles. At weak contrast, the threshold follows the universal scaling $N_c^{(0)}-N_c\sim\varepsilon^2$, whereas at moderate and strong contrast the third-order approximation is more accurate than the second-order theory. Wavelength selection is controlled by the ratio of the dominant material wavelength $λ^\ast$ to the harmonic-mean reference wavelength $λ_H$. For $λ^\ast/λ_H<1$, the harmonic-mean model accurately predicts the wrinkle wavelength. The framework provides a mechanics-based tool for reliability assessment and design of statistically heterogeneous film--substrate systems.