2026/07/31 by Xinyu Xing, Liyu Zhong, Feng Deng +1
Physics and Astronomy · #cond-mat.mtrl-sci
arxiv created 2026/07/31 · arxiv updated 2026/08/04
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 Nc(0)-Nc∼ε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 λ^∗ to the harmonic-mean reference wavelength λH. For λ^∗/λ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.