Wrinkling Simulations of Fiber-Reinforced Elastomer Sheets Under Global Tensile Loading

Numerical predictions of the wrinkling behavior of biaxially fiber-reinforced elastomer sheets are carried out under consideration of finite deformations. The Holzapfel–Gasser–Ogden material model is used to account for the anisotropic hyperelastic material behavior of the sheets, where material parameters are identified based on experimental data of tensile tests from literature. A Finite Element Method-based simulation strategy is presented to extract critical loading conditions and to access the postbuckling response using geometrical imperfections. Depending on the layup and aspect ratio of the sheets, wrinkling onset was predicted for global stretches between 10% and 25%. For sheets with fiber orientations [±45°] wrinkling is predicted at larger global stretches than for sheets with fiber orientations of [+30∘/−60∘] for the same aspect ratio. Furthermore, it is shown that short sheets have a tendency towards symmetric wrinkling patterns whereas for long sheets asymmetric wrinkles are more likely to occur. Comparison of the numerical predictions with experiments from the literature shows that the geometrical characteristics of the wrinkles, such as wavelengths and amplitudes, can be well predicted. Far into the postbuckling regime, the deviations of the predicted wrinkling amplitudes and their experimental counterparts are around 30% or less.