Research Papers

Wrinkling of a Polymeric Gel During Transient Swelling

[+] Author and Article Information
William Toh

School of Mechanical
and Aerospace Engineering,
Nanyang Technological University,
50 Nanyang Avenue,
639798, Singapore
e-mail: wtoh1@e.ntu.edu.sg

Zhiwei Ding

Institute of High Performance Computing,
1 Fusionopolis Way,
#16-16 Connexis,
138632, Singapore
e-mail: dingzw@ihpc.a-star.edu.sg

Teng Yong Ng

School of Mechanical
and Aerospace Engineering,
Nanyang Technological University,
50 Nanyang Avenue,
639798, Singapore
e-mail: mtyng@ntu.edu.sg

Zishun Liu

International Center for Applied Mechanics,
State Key Laboratory for Strength
and Vibration of Mechanical Structures,
Xi'an Jiaotong University,
No. 28, West Xianning Road,
Xi'an Shaanxi 710049, China
e-mail: zishunliu@mail.xjtu.edu.cn

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received January 22, 2015; final manuscript received April 2, 2015; published online April 30, 2015. Editor: Yonggang Huang.

J. Appl. Mech 82(6), 061004 (Jun 01, 2015) (8 pages) Paper No: JAM-15-1039; doi: 10.1115/1.4030327 History: Received January 22, 2015; Revised April 02, 2015; Online April 30, 2015

When exposed to an external solvent, a dry polymeric network imbibes the solvent and undergoes large deformation. The resulting aggregate is known as a hydrogel. This swelling process is diffusion driven and thus results in differential swelling during transient swelling. When subjected to external geometrical constraints, such as being rigidly fixed or attachment to a compliant substrate, wrinkles have been shown to appear due to mechanical instabilities. In the case of free swelling, there are no external constraints to induce the instabilities accounting for wrinkling patterns. However, during the transient swelling process, the swelling differential between the gel on the exterior and the interior causes compressive stresses and gives rise to mechanical instabilities. It is also observed that the time dependence of the swelling profile causes the wrinkles to evolve with time. In this work, we investigate this interesting phenomenon of transient wrinkle mode evolution using the finite element and state-space methods. From our simulations and prediction, we find that there is an inverse relation between critical wave number and time, which has earlier been observed in experiments.

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


Yang, S., Khare, K., and Lin, P.-C., 2010, “Harnessing Surface Wrinkle Patterns in Soft Matter,” Adv. Funct. Mater., 20(16), pp. 2550–2564. [CrossRef]
DuPont, S. J., Jr., Cates, R. S., Stroot, P. G., and Toomey, R., 2010, “Swelling-Induced Instabilities in Microscale, Surface-Confined Poly(N-Isopropylacryamide) Hydrogels,” Soft Matter, 6(16), pp. 3876–3882. [CrossRef]
Tanaka, T., Sun, S.-T., Hirokawa, Y., Katayama, S., Kucera, J., Hirose, Y., and Amiya, T., 1987, “Mechanical Instability of Gels at the Phase Transition,” Nature, 325(6107), pp. 796–798. [CrossRef]
Guvendiren, M., Burdick, J. A., and Yang, S., 2010, “Kinetic Study of Swelling-Induced Surface Pattern Formation and Ordering in Hydrogel Films With Depth-Wise Crosslinking Gradient,” Soft Matter, 6(9), pp. 2044–2049. [CrossRef]
Peixinho, J., and Mukhopadhyay, S., 2013, “Diffusion-Mechanical Instability of a Spherical Gel,” Rencontre du Non-linéaire (Meeting of Nonlinear 2015), Paris, Mar. 18–19, pp. 116–119. http://nonlineaire.univ-lille1.fr/SNL/media/2012/CR/Peixinho_.pdf
Liu, Z. S., Swaddiwudhipong, S., Cui, F. S., Hong, W., Suo, Z., and Zhang, Y. W., 2011, “Analytical Solutions of Polymeric Gel Structures Under Buckling and Wrinkle,” Int. J. Appl. Mech., 3(2), pp. 235–257. [CrossRef]
Zhang, Y., Chen, L., Swaddiwudhipong, S., and Liu, Z. S., 2014, “Buckling Deformation of Annular Plates Describing Natural Forms,” Int. J. Struct. Stab. Dyn., 14(1), p. 1350054. [CrossRef]
Mora, T., and Boudaoud, A., 2006, “Buckling of Swelling Gels,” Eur. Phys. J. E, 20(2), pp. 119–124. [CrossRef]
Liu, Z. S., Hong, W., Suo, Z., Swaddiwudhipong, S., and Zhang, Y., 2010, “Modeling and Simulation of Buckling of Polymeric Membrane Thin Film Gel,” Comput. Mater. Sci., 49(1), pp. S60–S64. [CrossRef]
Wu, Z., Bouklas, N., and Huang, R., 2013, “Swell-Induced Surface Instability of Hydrogel Layers With Material Properties Varying in Thickness Direction,” Int. J. Solids Struct., 50(3–4), pp. 578–587. [CrossRef]
Hong, W., Zhao, X., and Suo, Z., 2009, “Formation of Creases on the Surfaces of Elastomers and Gels,” Appl. Phys. Lett., 95(11), p. 111901. [CrossRef]
Duan, Z., Zhang, J., An, Y., and Jiang, H., 2013, “Simulation of the Transient Behavior of Gels Based on an Analogy Between Diffusion and Heat Transfer,” Int. J. Appl. Mech., 80(4), p. 041017. [CrossRef]
Toh, W., Liu, Z. S., Ng, T. Y., and Hong, W., 2013, “Inhomogeneous Large Deformation Kinetics of Polymeric Gels,” Int. J. Appl. Mech., 5(1), p. 1350001. [CrossRef]
Chester, S. A., Di Leo, C. V., and Anand, L., 2014, “A Finite Element Implementation of a Coupled Diffusion-Deformation Theory for Elastomeric Gels,” Int. J. Solids Struct., 52, pp. 1–18. [CrossRef]
Zhang, J., Zhao, X., Suo, Z., and Jiang, H., 2009, “A Finite Element Method for Transient Analysis of Concurrent Large Deformation and Mass Transport in Gels,” J. Appl. Phys., 105(9), pp. 93522–93529. [CrossRef]
Wu, Z., Meng, J., Liu, Y., Li, H., and Huang, R., 2014, “A State Space Method for Surface Instability of Elastic Layers With Material Properties Varying in Thickness Direction,” ASME J. Appl. Mech., 81(8), p. 081003. [CrossRef]
Kang, M. K., and Huang, R., 2010, “A Variational Approach and Finite Element Implementation for Swelling of Polymeric Hydrogels Under Geometric Constraints,” ASME J. Appl. Mech., 77(6), p. 061004. [CrossRef]
Hong, W., Liu, Z. S., and Suo, Z., 2009, “Inhomogeneous Swelling of a Gel in Equilibrium With a Solvent and Mechanical Load,” Int. J. Solids Struct., 46(17), pp. 3282–3289. [CrossRef]
Toh, W., Ng, T. Y., Hu, J., and Liu, Z., 2014, “Mechanics of Inhomogeneous Large Deformation of Photo-Thermal Sensitive Hydrogels,” Int. J. Solids Struct., 51(25–26), pp. 4440–4451. [CrossRef]
Ding, Z., Liu, Z. S., Hu, J., Swaddiwudhipong, S., and Yang, Z., 2013, “Inhomogeneous Large Deformation Study of Temperature-Sensitive Hydrogel,” Int. J. Solids Struct., 50(16–17), pp. 2610–2619. [CrossRef]
Marcombe, R., Cai, S., Hong, W., Zhao, X., Lapusta, Y., and Suo, Z., 2010, “A Theory of Constrained Swelling of a pH-Sensitive Hydrogel,” Soft Matter, 6(4), pp. 784–793. [CrossRef]
Bouklas, N., and Huang, R., 2012, “Swelling Kinetics of Polymer Gels: Comparison of Linear and Nonlinear Theories,” Soft Matter, 8(31), pp. 8194–8203. [CrossRef]
Lee, D., Triantafyllidis, N., Barber, J. R., and Thouless, M. D., 2008, “Surface Instability of an Elastic Half Space With Material Properties Varying With Depth,” J. Mech. Phys. Solids, 56(3), pp. 858–868. [CrossRef] [PubMed]
Weiss, F., Cai, S., Hu, Y., Kyoo Kang, M., Huang, R., and Suo, Z., 2013, “Creases and Wrinkles on the Surface of a Swollen Gel,” J. Appl. Phys., 114(7), p. 073507. [CrossRef]
Cao, Y., and Hutchinson, J. W., 2011, “From Wrinkles to Creases in Elastomers: The Instability and Imperfection-Sensitivity of Wrinkling,” Proc. R. Soc. A, 468(2137), pp. 94–115. [CrossRef]
Wong, W. H., Guo, T. F., Zhang, Y. W., and Cheng, L., 2010, “Surface Instability Maps for Soft Materials,” Soft Matter, 6(22), pp. 5743–5750. [CrossRef]
Hohlfeld, E., and Mahadevan, L., 2012, “Scale and Nature of Sulcification Patterns,” Phys. Rev. Lett., 109(2), p. 025701. [CrossRef] [PubMed]
Hohlfeld, E., and Mahadevan, L., 2011, “Unfolding the Sulcus,” Phys. Rev. Lett., 106(10), p. 105702. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Schematics of a gel exposed to solvent of chemical potential μs during (a) initial state of homogeneous swelling at chemical potential of μ0 and (b) onset of wrinkling at the surface at time t

Grahic Jump Location
Fig. 2

(a) Compressive stresses experienced by a gel (Nν=0.001,χ=0.1,and μ0=-1) swelling under lateral constrains in kinetic analysis and (b) stretch ratio at various times

Grahic Jump Location
Fig. 3

(a) Discretization of a single waveform into eight elements and (b) imposition of MPC at surface nodes

Grahic Jump Location
Fig. 4

(a) Full finite element model of the gel, with thickness 1 and width of 10. (b)–(d) Different meshes with top layer mesh thickness of 0.01 H, 0.005 H, and 0.001 H, respectively. Note that meshes (b)–(d) are shortened with width to illustrate the graded mesh used in simulations.

Grahic Jump Location
Fig. 5

Time of onset of wrinkling with increasing refinement of mesh for models with 10 and 15 waves within the width of the gel

Grahic Jump Location
Fig. 6

Evolution of surface wrinkles of a hydrogel layer. Time of onset of wrinkling increases from (a) to (f).

Grahic Jump Location
Fig. 7

(a) wave number versus time and (b) wavelength versus time for gels of crosslink densities at Nν=0.001 and Nν=0.01; Flory-interaction parameter χ=0.1; and initial chemical potential of μ0=-1

Grahic Jump Location
Fig. 8

Comparison between state-space predictions for different discretization of gel layers

Grahic Jump Location
Fig. 9

Normalized wave number versus time for finite element method and state-space methods for Nν=0.01,χ=0.1,and μ0=-1



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