Research Papers

Geometric Prepatterning-Based Tuning of the Period Doubling Onset Strain During Thin-Film Wrinkling

[+] Author and Article Information
Sourabh K. Saha

Materials Engineering Division,
Lawrence Livermore National Laboratory,
PO Box 808,
Livermore, CA 94550
e-mail: saha5@llnl.gov

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received February 17, 2017; final manuscript received March 20, 2017; published online April 5, 2017. Assoc. Editor: Yihui Zhang.The United States Government retains, and by accepting the article for publication, the publisher acknowledges that the United States Government retains, a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States government purposes.

J. Appl. Mech 84(5), 051010 (Apr 05, 2017) (10 pages) Paper No: JAM-17-1089; doi: 10.1115/1.4036325 History: Received February 17, 2017; Revised March 20, 2017

Wrinkling of thin films is an easy-to-implement and low-cost technique to fabricate stretch-tunable periodic micro and nanoscale structures. However, the tunability of such structures is often limited by the emergence of an undesirable period-doubled mode at high strains. Predictively tuning the onset strain for period doubling via existing techniques requires one to have extensive knowledge about the nonlinear pattern formation behavior. Herein, a geometric prepatterning-based technique is introduced that can be implemented even with limited system knowledge to predictively delay period doubling. The technique comprises prepatterning the film/base bilayer with a sinusoidal pattern that has the same period as the natural period of the system. This technique has been verified via physical and computational experiments on the polydimethylsiloxane (PDMS)/glass bilayer system. It is observed that the onset strain can be increased from the typical value of 20% for flat films to greater than 30% with a modest prepattern aspect ratio (2·amplitude/period) of 0.15. In addition, finite element simulations reveal that (i) the onset strain increases with increasing prepattern amplitude and (ii) the delaying effect can be captured entirely by the prepattern geometry. Therefore, one can implement this technique even with limited system knowledge, such as material properties or film thickness, by simply replicating pre-existing wrinkled patterns to generate prepatterned bilayers. Thus, geometric prepatterning is a practical scheme to increase the operating range of stretch-tunable wrinkle-based devices by at least 50%.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Ahmed, S. F. , Rho, G.-H. , Lee, K.-R. , Vaziri, A. , and Moon, M.-W. , 2010, “ High Aspect Ratio Wrinkles on a Soft Polymer,” Soft Matter, 6(22), pp. 5709–5714. [CrossRef]
Bowden, N. , Brittain, S. , Evans, A. G. , Hutchinson, J. W. , and Whitesides, G. M. , 1998, “ Spontaneous Formation of Ordered Structures in Thin Films of Metals Supported on an Elastomeric Polymer,” Nature, 393(6681), pp. 146–149. [CrossRef]
Chen, X. , and Yin, J. , 2010, “ Buckling Patterns of Thin Films on Curved Compliant Substrates With Applications to Morphogenesis and Three-Dimensional Micro-Fabrication,” Soft Matter, 6(22), pp. 5667–5680. [CrossRef]
Genzer, J. , and Groenewold, J. , 2006, “ Soft Matter With Hard Skin: From Skin Wrinkles to Templating and Material Characterization,” Soft Matter, 2(4), pp. 310–323. [CrossRef]
Kim, J. B. , Kim, P. , Pegard, N. C. , Oh, S. J. , Kagan, C. R. , Fleischer, J. W. , Stone, H. A. , and Loo, Y.-L. , 2012, “ Wrinkles and Deep Folds as Photonic Structures in Photovoltaics,” Nat. Photonics, 6(5), pp. 327–332. [CrossRef]
Lin, P. C. , and Yang, S. , 2007, “ Spontaneous Formation of One-Dimensional Ripples in Transit to Highly Ordered Two-Dimensional Herringbone Structures Through Sequential and Unequal Biaxial Mechanical Stretching,” Appl. Phys. Lett., 90(24), p. 241903. [CrossRef]
Mei, Y. , Kiravittaya, S. , Harazim, S. , and Schmidt, O. G. , 2010, “ Principles and Applications of Micro and Nanoscale Wrinkles,” Mater. Sci. Eng., R, 70(3–6), pp. 209–224. [CrossRef]
Watanabe, M. , and Mizukami, K. , 2012, “ Well-Ordered Wrinkling Patterns on Chemically Oxidized Poly(Dimethylsiloxane) Surfaces,” Macromolecules, 45(17), pp. 7128–7134. [CrossRef]
Ohzono, T. , and Monobe, H. , 2012, “ Microwrinkles: Shape-Tunability and Applications,” J. Colloid Interface Sci., 368(1), pp. 1–8. [CrossRef] [PubMed]
Chiche, A. , Stafford, C. M. , and Cabral, J. T. , 2008, “ Complex Micropatterning of Periodic Structures on Elastomeric Surfaces,” Soft Matter, 4(12), pp. 2360–2364. [CrossRef]
Groenewold, J. , 2001, “ Wrinkling of Plates Coupled With Soft Elastic Media,” Physica A, 298(1–2), pp. 32–45. [CrossRef]
Jiang, H. , Khang, D.-Y. , Song, J. , Sun, Y. , Huang, Y. , and Rogers, J. A. , 2007, “ Finite Deformation Mechanics in Buckled Thin Films on Compliant Supports,” Proc. Natl. Acad. Sci., 104(40), pp. 15607–15612. [CrossRef]
Saha, S. K. , and Culpepper, M. L. , 2015, “ Design of a Compact Biaxial Tensile Stage for Fabrication and Tuning of Complex Micro- and Nano-Scale Wrinkle Patterns,” ASME J. Micro Nano-Manuf., 3(4), p. 041004. [CrossRef]
Huh, D. , Mills, K. L. , Zhu, X. , Burns, M. A. , Thouless, M. D. , and Takayama, S. , 2007, “ Tuneable Elastomeric Nanochannels for Nanofluidic Manipulation,” Nat. Mater., 6(6), pp. 424–428. [CrossRef] [PubMed]
Chung, S. , Lee, J. H. , Moon, M.-W. , Han, J. , and Kamm, R. D. , 2008, “ Non-Lithographic Wrinkle Nanochannels for Protein Preconcentration,” Adv. Mater., 20(16), pp. 3011–3016. [CrossRef]
Harrison, C. , Stafford, C. M. , Zhang, W. , and Karim, A. , 2004, “ Sinusoidal Phase Grating Created by a Tunably Buckled Surface,” Appl. Phys. Lett., 85(18), pp. 4016–4018. [CrossRef]
Yu, C. , O’Brien, K. , Zhang, Y.-H. , Yu, H. , and Jiang, H. , 2010, “ Tunable Optical Gratings Based on Buckled Nanoscale Thin Films on Transparent Elastomeric Substrates,” Appl. Phys. Lett., 96(4), p. 041111. [CrossRef]
Da, Y. , Jing, F. , Rui, M. , Xu-Lin, Z. , Yue-Feng, L. , Tong, Y. , and Hong-Bo, S. , 2015, “ Stability Improved Stretchable Metallic Gratings With Tunable Grating Period in Submicron Scale,” J. Lightwave Technol., 33(15), pp. 3327–3331. [CrossRef]
Brau, F. , Vandeparre, H. , Sabbah, A. , Poulard, C. , Boudaoud, A. , and Damman, P. , 2011, “ Multiple-Length-Scale Elastic Instability Mimics Parametric Resonance of Nonlinear Oscillators,” Nat. Phys., 7(1), pp. 56–60. [CrossRef]
Zhao, Y. , Cao, Y. , Hong, W. , Wadee, M. K. , and Feng, X.-Q. , 2015, “ Towards a Quantitative Understanding of Period-Doubling Wrinkling Patterns Occurring in Film/Substrate Bilayer Systems,” Proc. R. Soc. London A, 471(2173), p. 20140695.
Budday, S. , Kuhl, E. , and Hutchinson, J. W. , 2015, “ Period-Doubling and Period-Tripling in Growing Bilayered Systems,” Philos. Mag., 95(28–30), pp. 3208–3224. [CrossRef]
Cao, Y. , and Hutchinson, J. W. , 2012, “ Wrinkling Phenomena in Neo-Hookean Film/Substrate Bilayers,” ASME J. Appl. Mech., 79(3), p. 031019. [CrossRef]
Auguste, A. , Jin, L. , Suo, Z. , and Hayward, R. C. , 2014, “ The Role of Substrate Pre-Stretch in Post-Wrinkling Bifurcations,” Soft Matter, 10(34), pp. 6520–6529. [CrossRef] [PubMed]
Chen, Y.-C. , and Crosby, A. J. , 2014, “ High Aspect Ratio Wrinkles Via Substrate Prestretch,” Adv. Mater., 26(32), pp. 5626–5631. [CrossRef] [PubMed]
Saha, S. K. , and Culpepper, M. L. , 2016, “ Deterministic Switching of Hierarchy During Wrinkling in Quasi-Planar Bilayers,” Adv. Eng. Mater., 18(6), pp. 938–943. [CrossRef]
Bowden, N. , Huck, W. T. S. , Paul, K. E. , and Whitesides, G. M. , 1999, “ The Controlled Formation of Ordered, Sinusoidal Structures by Plasma Oxidation of an Elastomeric Polymer,” Appl. Phys. Lett., 75(17), pp. 2557–2559. [CrossRef]
Saha, S. K. , 2014, “ Predictive Design and Fabrication of Complex Micro and Nano Patterns Via Wrinkling for Scalable and Affordable Manufacturing,” Ph.D. thesis, Massachusetts Institute of Technology, Cambridge, MA.
Saha, S. K. , and Culpepper, M. L. , 2016, “ Method to Fabricate Pre-Patterned Surfaces During Manufacture of Complex Wrinkled Structures,” U.S. Patent Application No. US20160039142.
Saha, S. K. , 2017, “ Sensitivity of the Mode Locking Phenomenon to Geometric Imperfections During Wrinkling of Supported Thin Films,” Int. J. Solids Struct., 109, pp. 166–179. [CrossRef]
Hutchinson, J. W. , 2013, “ The Role of Nonlinear Substrate Elasticity in the Wrinkling of Thin Films,” Philos. Trans. R. Soc. A, 371(1993), p. 20120422.
Wang, Q. , and Zhao, X. , 2014, “ Phase Diagrams of Instabilities in Compressed Film-Substrate Systems,” ASME J. Appl. Mech., 81(5), p. 051004. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic representation of the formation and growth of wrinkles with strain during compression of (a) flat bilayers and (b) prepatterned bilayers. Onset of period doubling is delayed in prepatterned bilayers.

Grahic Jump Location
Fig. 2

Demonstration of delayed onset of period doubling via prepatterning. (a)–(c) Fabricated wrinkles on top of flat and prepatterned bilayers. Finite element modeling of wrinkles on (d) flat and (e) prepatterned bilayers. The prepattern amplitude, strain, and natural period were 168 nm, 0.249, and 2.0 μm for both experiments and simulations. Scale bars are 20 μm long for (a), 5 μm for (b) and (c), and 2.5 μm for (d) and (e).

Grahic Jump Location
Fig. 3

Wrinkle patterns on three different bilayer samples, each with a pair of flat and prepatterned regions. (a)–(c) Period-doubled wrinkles on flat regions and (d)–(f) single-period wrinkles on the corresponding prepatterned regions. λn, εp, and ε for (a) and (d) were 2.3 μm, 0.04, and 0.255; for (b) and (e) were 2.2 μm, 0.07, and 0.249; and for (c) and (f) were 2.3 μm, 0.06, and 0.245. Scale bars are 10 μm long.

Grahic Jump Location
Fig. 4

(a) Finite element modeling of wrinkle formation during compression of a prepatterned bilayer. Color information (green-to-red) represents increasing first principal strain in the system. (b) Growth in amplitude of wrinkles with increasing strain. The prepattern amplitude, natural period, thin-film thickness, and Young’s moduli ratio were 168 nm, 2.0 μm, 42 nm, and 1690.4. (For interpretation of the references to color in this figure, the reader is referred to the web version of this article)

Grahic Jump Location
Fig. 5

Effect of strain on the ratio of the strain energy in the base for two competing wrinkle modes. The smaller period (λ) was 2.25 μm, and the prepattern amplitudes were 0.084 times the periods.

Grahic Jump Location
Fig. 6

Effect of prepattern size on the period doubling onset strain quantified using finite element simulations. Onset strain was evaluated for different natural periods by separately varying film thickness (h) and moduli ratio (η). The natural periods were 1 μm, 2.25 μm, and 4.5 μm. The parameter ho was 47.5 nm and ηo was 1690.4.

Grahic Jump Location
Fig. 7

Increased in amplitude and aspect ratio due to increase in prepattern size. (a) Amplitude versus strain. (b) Aspect ratio versus strain. The natural period was 2.25 μm, the film thickness was 47.5 nm, and the moduli ratio was 1690.4 for the two bilayers.

Grahic Jump Location
Fig. 8

Effect of prepattern size on the maximum aspect ratio of wrinkles. Maximum aspect ratio is the aspect ratio at the period doubling onset strain. Aspect ratio was evaluated for different natural periods by separately varying film thickness (h) and moduli ratio (η). The natural periods were 1 μm, 2.25 μm, and 4.5 μm. The parameter ho was 47.5 nm and ηo was 1690.4.

Grahic Jump Location
Fig. 9

Finite element model of (a) prepatterned bilayer and (b) flat bilayer for the case when the prepattern and natural periods are not identical. The natural period of the bilayers was 2.25 μm, and the prepattern period and amplitude were 3 μm and 36 nm. The onset strain for the flat bilayer is 0.19, whereas the onset strain for the prepatterned bilayer is 0.225. The strain in the two bilayers is 0.205. Color information (green-to-red) represents increasing first principal strain in the system. Scale bars are 2 μm long. (For interpretation of the references to color in this figure, the reader is referred to the web version of this article)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In