超弹性薄膜与可压缩基底双层结构表面失稳分析1)

黄春阳 唐山彭向和

(重庆大学航空航天学院，重庆400017)

创刊60周年专栏

超弹性薄膜与可压缩基底双层结构表面失稳分析1)

黄春阳 唐山2)彭向和3)

(重庆大学航空航天学院，重庆400017)

当上层超弹性硬质薄膜和下层可膨胀基底构成的双层结构受压时，薄膜的自由表面可通过形成褶皱降低系统能量.研究表明,上下两层的模量比不同时，上层弹性硬质薄膜将表现出不同的表面失稳模式.本文提出了一种新颖的方法可有效抑制双层软材料的表面失稳，即改变基底材料的泊松比，这种方法同时适用于不具有应变硬化的软材料.首先基于Neo-Hookean模型发展了小变形条件下双层结构表面失稳的理论模型，通过半解析的方法得到了表面失稳的临界应变；然后通过有限元计算与模拟，进一步验证了负泊松比基底可延缓表面失稳.结果表明：(1)当双层结构基底泊松比为正且趋于0.5(不可压缩)时，双层结构在较小的压缩应变下出现表面失稳；(2)当基底的泊松比为负且趋于-1时，可被压缩至46%而不出现表面失稳，即可膨胀基底能有效抑制薄膜的表面失稳.本文发展的方法及主要结果可为延展性电子器件的设计提供指导.

薄膜--基底结构,表面失稳,负泊松比,可延展性电子

引言

Biot的线性扰动分析认为，平面应变条件下，半无限超弹性材料受压时，如果施加压缩应变达到临界值εBiot=0.46，自由表面将丧失稳定性，形成波浪形褶皱[1].然而实验[24]、理论[47]以及计算[8-9]表明当施加压缩应变达到εcrease=0.35时，自由表面将形成折痕.超弹性块体在弯曲载荷下通常也可观察到折痕.折痕不同于褶皱，因为折痕临界应变小于褶皱临界应变εcrease＜εBiot.但是具有硬质薄膜和厚软基底构成的双层结构受到压缩时，则会出现波浪形褶皱以降低系统能量[1011].褶皱是自由表面整体无限小扰动，而折痕则是局部大应变下的自我接触.折痕的接触区域有尖点，不同于褶皱.这些特殊表面形貌可应用到生物仿生[12]、生物淤积[13]、表面粘附[14]和细胞工程[15]等.但是表面失稳[16,17]也会导致材料失效，比如微电子系统[1819]和复合材料的屈曲[2021].因此人们容易联想到的一个的问题是：如何抑制或者推迟表面失稳？

为回答上述问题，Cao和Hutchinson[22]通过力学理论和计算研究了超弹性硬质薄膜在超弹性软基底上的双层结构在受压下的失稳[22]，结果表明褶皱出现的临界应变εwrinkle取决于薄膜的杨氏模量Ef与基底的杨氏模量Es之比Ef/Es.例如Ef/Es=20时，εwrinkle=0.07[18]；Ef/Es=2时，εwrinkle=0.30，临界应变小于εcrease和εBiot[22].Feng和Tian[23]，Wang和Zhao[24]详细分析了受压薄膜--基底双层结构的表面失稳.考虑了薄膜与基底的杨氏模量，界面的粘结强度、薄膜厚度和基底的预拉伸等因素对表面失稳的影响.Jin和Suo[25]最近的研究表明：材料硬化可抑制表面折痕出现.随压缩应变增大，平直表面会失稳形成折痕.随着压缩应变增大和材料应变硬化，折痕将随之消失.这种抑制表面折痕的方法很大程度上依赖材料的本构行为，不具有普遍性.

本文提出了一种新颖的方法抑制双层软材料表面失稳，可适用于不具有应变硬化的软材料.理论分析表明，当基底为负泊松比材料时，表面失稳临界应变可被推迟到εBiot.为验证理论分析的结果，建立了有限元模型，研究了受压下薄膜与基体表面失稳的形变过程.有限元分析结果表明，所提出的方法能够抑制或延缓折痕、褶皱、倍褶皱或者折叠等表面失稳模式[2627].当基底材料为正泊松比时，没有同样的效应.因此，负泊松比基底是抑制双层软材料结构表面失稳的重要因素.

在研究材料的表面失稳时，先前的理论[1,4,7,10]和计算[89,22,24]通常假设材料为不可压缩的高分子弹性体(泊松比为0.5).但软材料或多或少具有一定的压缩性.可压缩性与温度T及材料的微观结构有关，例如高分子聚合物在温度T高于自身玻璃化温度Tg时，往往不可压缩；但是当T＜Tg时，往往具有可压缩性.Greaves等[28]对此做了详细的综述，分子动力学模拟也证明高分子材料的泊松比与温度有关[29].对于各项同性材料，泊松比ν满足-1＜ν＜0.5.一些研究者通过设计具有凹角或者弧形的网状结构，得到具有负泊松比性质的材料[2832].

1 表面失稳理论模型

如图1表示，考虑薄膜--基体双层结构.上层为硬质薄膜，下层为基底.为考虑泊松比效应，采用可压缩Neo-Hookean模型[29,33]，其自由能函数为

假设上层薄膜初始厚度为h，下层基底初始厚度为H.可压缩基底可通过微观结构设计得到[32]，比如孔洞.即使存在微观结构，依然假设其为均质超弹性材料，以方便分析.薄膜和基底材料均假设为可压缩的Neo-Hookean材料.一般情况下，考虑无穷厚基底，假设h/H=100.在下面的符号中，下角标“s”和“f”分别表示基底和薄膜.

图1 薄膜--基底双层结构示意图Fig.1 Schematic ofa thin film-substrat bilayer

在大变形和平面应变条件下，一种半解析方法可用于研究多层结构的表面失稳问题.它主要包括两个分析步骤：第一步，统一施加沿x1方向的预拉伸λ;第二步，在当前构型下，施加任意波长的扰动.

在Lagrangian构型中，增量问题通过下式描述

式中Div表示在初始构型下的散度算子.通常在即时构型中构造失稳问题，引入变换[35]

平衡方程又可以写为

式中div表示即时构型下的散度算子.

取以下形式的扰动

式中K为波数.

假设薄膜与基底不存在界面分离和相对滑移，即在薄膜和基底界面处[36]

其中[·]表示跳跃算子，定义为

上标“+”和“-”表示薄膜与基体界面的上方和下方.

在平面应变条件下，连续性条件可表示为

表面失稳问题的边界条件为：

上表面

下表面

求解平衡方程可得由待定未知系数表示的v1和v2的解.按照图1中讨论的双层结构，连续条件和边界条件共有8个方程和8个待定未知系数.存在非平凡解的条件为8个方程的系数矩阵行列式等于零.由此可得表面失稳的临界应变.在小变形条件下，可得下式表示的表面失稳的临界应变[37]

其中,A0为上层薄膜褶皱的振幅，n为褶皱的数量，H为基底厚度.

本文考虑了两种不同的薄膜和基底模量比：Ef/Es=1和Ef/Es=20.

假设上层薄膜不可压，即νf=0.5，基底泊松比νs在-1到0.5之间变化.图2显示了预测的临界失稳应变和基底泊松比的关系.其中定义临界应变εcrit=1-λcrit，结果与Biot[1]，Cao和Hutchnison[22]的结果一致.从图2中可见，当Ef/Es=1，对于不同的νs，褶皱出现的临界应变均为0.46，与Biot基于均匀不可压超弹性材料预测结果相同[1].当Ef/Es=20时，褶皱出现的临界应变随基底泊松比νs减小而单调增加.当基底泊松比为正时，泊松比对于褶皱出现的影响可忽略；当基底的泊松比从0到-1变化时，褶皱出现临界应变剧烈增加(从0.08增加到0.46).当νs接近于-1时，临界应变趋近εBiot=0.46.尽管只展示了Ef/Es=1和Ef/Es=20两种情形，但详尽的计算表明，使用其他模量比的结果也与此规律相符.综上所述，受压情况下，基底为负泊松比材料可有效抑制和推迟表面失稳发生，临界应变的上限εBiot=0.46.

图2 临界应力εcrit随基底泊松比νs的变化Fig.2 Criticalstrainεcritatonsetof surfacew rinklesvs.Poisson’s ratio of substrateνs

2 表面失稳的有限元分析

本节将应用有限元模型进一步验证负泊松比基底可延缓表面失稳的结论.本文所有有限元计算均采用商用软件ABAQUS完成[33].基于之前的工作建立有限元模型[29,38]，与Cao和Hutchinson[22]的工作类似.有限元模型采用CPE4H单元，单元总数为20000.分析主要分两步完成：(1)线性扰动分析；(2)后屈曲分析.第一步采用线性扰动分析，得出基础失稳模态；第二步引入微小缺陷(约为0.05h)乘以基础失稳模态，对初始有限元网格进行扰动，进行后屈曲分析.

采用与前节的理论分析相同的边界条件，其他参数Ef/Es=20，νf=0.5.图3显示了不同压缩应变下双层结构表面形貌.首先讨论基底泊松比νs=0.4的结果.周期褶皱出现临界应变为εcrit=0.085，如图3(a)所示，与理论分析相一致.表面正弦褶皱的幅值逐步增大.当压缩应变接近0.18时，出现二次分岔.随着压缩应变进一步增加到0.24，倍褶皱会出现.与Cao和Hutchinson所观察到的现象相符.若基底的泊松比νs=-0.9，在压缩应变为0.28时，薄膜表面并没有产生褶皱或者折痕，如图3(b)所示.当压缩应变到0.37时，表面出现很小幅值的褶皱.压缩应变0.37已经大于均匀材料中出现折痕的临界应变值εcrease=0.35[4,7].随着压缩应变增加到0.40，表面折痕会在薄膜表面形成，但其深度很小，很难用肉眼观察到.对比图2所示的理论分析结果，可见有限元模拟结果与理论分析结果几乎一致.

图3 不同压缩应变下薄膜--基底结构表面形貌Fig.3 Surfacemorphology of thin film-substrat system atdi ff erent compressive strain

3 结论

目前，高分子聚合物基底被广泛用于提高金属和半导体薄膜的延展性[3943].比如铜薄膜可在Kapton基底上达到超过50%的应变[39]；银薄膜可延展至46%的应变[40].这些研究可帮助设计可延展性电子产品[42].本文提出了一种新的方法在压缩条件下抑制或者推迟双层基底与薄膜的表面失稳.负泊松比基底与超弹性薄膜构成的双层结构，上层薄膜可压缩至46%而未出现表面失稳.考虑到压缩和拉伸同等重要性，该方法为延展性电子产品设计提供了一种新思路.

1 BiotMA.Surface instability of rubber in compression.Flow,Turbulence and Combustion,1963,12(2)：168-182

2 GentAN,Cho IS.Surface instabilities in compressed orbent rubber blocks.RubberChemistry and Technology,1999,72(2)：253-262

3 Ghatak A,Das AL.Kink instability of a highly deformable elasticcylinder.PhysicalReview Letters,2007,99(7)：076101

4 Hong W,Zhao X,Suo Z.Formation of creases on the surfaces of elastomersand gels.Applied Physics Letters,2009,95(11)：111901

5 Needleman A,Tvergaard V,Van derGiessen E.Indentation of elastically soft and plastically compressible solids.Acta Mechanica Sinica,2015,31(4)：473-480

6 Huang SQ,Feng XQ.Spinodalsurface instability of softelastic thin films Acta Mechanica Sinica,2008,24(3)：289-296

7 Hohlfeld E,Mahadevan L.Unfolding the sulcus.Physical Review Letters,2011,106(10)：105702

8 Cao Y,Hutchinson JW.From w rinkles to creases in elastomers：the instability and imperfection-sensitivity of w rinkling//Proc.R.Soc.A.The RoyalSociety,2012,468(2137)：94-115

9 Wong WH,Guo TF,Zhang YW,etal.Surface instabilitymaps for softmaterials.SoftMatter,2010,6(22)：5743-5750

10 Mei H,Huang R,Chung JY,et al.Bucklingmodes of elastic thin film on elastic substrates.Applied Physics Letters,2007,90(15)：151902

11 Sun JY,Xia S,Moon MW,et al.Folding w rinkles of a thin sti ff layeron a softsubstrate//Proc.R.Soc.A.The RoyalSociety,2012,468(2140)：932-953

12 Kim J,Yoon J,Hayward RC.Dynam ic display of biomolecular patterns throughan elastic creasing instability ofstimuli-responsivehydrogels.Nature Materials,2010,9(2)：159-164

13 Shivapooja P,Wang Q,Orihuela B,etal.Bioinspired surfacesw ith dynamic topography foractive controlofbiofouling.Advanced Materials,2013,25(10)：1430-1434.

14 Chan EP,Karp JM,Langer RS.A“self-pinning”adhesive based on responsive surface w rinkles.Journal of Polymer Science Part B:PolymerPhysics,2011,49(1)：40-44

15 Saha K,Kim J,Irw in E,et al.Surface creasing instability of soft polyacrylam ide cell culture substrates.Biophysical Journal,2010,99(12)：L94-L96

16 Xin F,Lu TJ.Acoustomechanical constitutive theory for softmaterials.Acta Mechanica Sinica,2016,32(5)：828-840

17 HuiCY.Crack buckling in soft gels under compression.Acta Mechanica Sinica,2012,28(4)：1098-1105

18 Krylov S,Ilic BR,Schreiber D,etal.The pull-in behavior of electrostatically actuated bistablem icrostructures.Journal ofM icromechanicsand Microengineering,2008,18(5)：055026

19 Li Y,Wang XS,Meng XK.Buckling behavior of metal fil/substrate structure under pure bending.Applied Physics Letters,2008,92(13)：131902

20 Bolotin VV.Delaminationsin compositestructures：Itsorigin,buckling,grow th and stability.Composites Part B:Engineering,1996,27(2)：129-145

21 Hu Y,Hiltner A,Baer E.Buckling in elastomer/plastic/elastomer 3-layer films PolymerComposites,2004,25(6)：653-661

22 Cao Y,Hutchinson JW.W rinkling phenomena in neo-Hookean fil/substrate bilayers.JournalofApplied Mechanics,2012,79(3)：031019

23 Feng PL,Lu TJ.A three-layer structuremodel for the e ff ect of a softmiddle layer on Lovewaves propagating in layered piezoelectric systems.Acta Mechanica Sinica,2012,28(4)：1087-1097

24 Wang Q,Zhao X.Phase diagrams of instabilities in compressed film-substrat systems.JournalofApplied Mechanics,2014,81(5)：051004

25 Jin L,Suo Z.Smoothening creases on surfaces of strain-sti ff ening materials.Journal of the Mechanics and Physics of Solids,2015,74：68-79

26 Li Z,Zhou Z,Li Y,et al.E ff ectof cyclic loading on surface instability of silicone rubber under compression.Polymers,2017,9(4)：148

27 Tang S,Gao B,Zhou Z,et al.Dimension-controlled formation of crease patternson softsolids.SoftMatter,2017,13(3)：619-626

28 GreavesGN,Greer AL,LakesRS,etal.Poisson’s ratio andmodern materials.NatureMaterials,2011,10(11)：823-837

29 Tang S,LiY,LiuWK,etal.Surface ripplesofpolymeric nanofiber under tension：the crucial role of Poisson’s ratio.Macromolecules,2014,47(18)：6503-6514

30 M ilton GW.Compositematerialsw ith Poisson’s ratios close to—1.JournaloftheMechanicsand Physics ofSolids,1992,40(5)：1105-1137

31 Grima JN,A lderson A,Evans KE.Auxeticbehaviour from rotating rigid units.Physica Status Solidi(b),2005,242(3)：561-575

32 BabaeeS,Shim J,Weaver JC,etal.3D Softmetamaterialsw ithnegative Poisson’s ratio.Advanced Materials,2013,25(36)：5044-5049

33 Hibbitt,Karlsson,Sorensen.ABAQUS：Theory Manual.Hibbitt,Karlsson&Sorensen,1997

34 Belytschko T,LiuWK,Moran B,etal.Nonlinear finit elements for continuaand structures.New York：JohnWiley&Sons,2013

35 Bertoldi K,Gei M.Instabilities in multilayered soft dielectrics.Journal of the Mechanics and Physics of Solids,2011,59(1)：18-42

36 Li Y,Wang XS,Fan Q.E ff ects of elastic anisotropy on the surface stability of thin fil/substrate system.International JournalofEngineering Science,2008,46(12)：1325-1333

37 Tang S,Li Y,Yang Y,etal.The e ff ectofmechanical-riven volumetric change on instability patterns of bilayered soft solids.Soft Matter,2015,11：7911-7919

38 Tang S,Greene MS,Peng XH,etal.Chain confinemen drives the mechanical properties of nanoporous polymers.EPL(Europhysics Letters),2014,106(3)：36002

39 Lu N,Wang X,Suo Z,et al.Metal film on polymer substrates stretched beyond 50%.Applied Physics Letters,2007,91(22)：221909

40 Xiao J,Carlson A,Liu ZJ,et al.Stretchable and compressible thin film of sti ff materials on compliantwavy substrates.Applied Physics Letters,2008,93(1)：013109

41 Kim DH,Ahn JH,ChoiWM,etal.Stretchableand foldable silicon integrated circuits.Science,2008,320(5875)：507-511

42 Rogers JA,Someya T,Huang Y.Materials and mechanics for stretchableelectronics.Science,2010,327(5973)：1603-1607

43 Xu F,LuW,Zhu Y.Controlled 3D buckling of silicon nanow ires for stretchableelectronics.Acs Nano,2010,5(1)：672-678

STUDY OFSURFACE INSTABILITY ABOUT HYPERELASTIC FILMSON AUXETIC SUBSTRATESUNDER COMPRESSION1)

Huang Chunyang Tang Shan2)Peng Xianghe3)
(College ofAerospace Engineering，Chongqing University，Chongqing 400017，China)

When abilayerstructure consisting ofa thin sti ff fil and a thick compliantsubstrate subjected to compressive deformation,its free surfacewould bew rinkled tom inim ize the energy of the system,and di ff erentw rinkle patternsmay appear fordi ff erentratiosof themodulusof the fil to thatof thesubstrate.In thisarticle,wedeveloped anovelapproach to suppress the surface instability of such bilayermaterialsunder severe compression by adjusting the Poisson’s ratio of the substrates.Thisapproach isalso applicable to thebilayer consisting ofa softsubstrateand a fil w ith elasticmodulus sim ilar to thatof the substrate.We developed an analyticalapproach for surface instability of the bilayer based on Neo-Hookeanmodel in the case of small deformation,and obtained the critical strain of the bilayerw ith a semi-analytical method.Then,we used finit elementapproach(FEA)to illustrate that the instability of the thin fil can be delayed if the substrate hasa negative Poisson’s ratio.We showed that：(1)when the Poisson’s ratio of the substrate is positive and close to 0.5(nearly incompressible),the surface instabilitymay occur to the bilayer system ata very small compressive strain;(2)if the Poisson’s ratio of thesubstrate isnegativeand close to-1,the fil can be compressed up to 46%w ithout occurence surface instability.The approach developed and the results obtained in this article imply a greatpotentialofauxeticmaterialsused to enhance the compressibility of thin films which can provideguidance for the design of lam inate ductileelectronic devices.

thin film-substrat system,surface instability,negative Poisson’s ratio,stretchableelectronics

O34

A

10.6052/0459-1879-17-161

2017-05-08收稿，2017-05-31录用,2017-05-31网络版发表.

1)国家自然科学基金资助项目(11472065).

2)唐山，教授，主要研究方向：跨尺度数值模拟与力学表征.E-mail：shan tang 0917@163.com

3)彭向和，教授，主要研究方向：跨尺度数值模拟与力学表征.E-mail：xhpeng@cqu.edu.cn

黄春阳,唐山,彭向和.超弹性薄膜与可压缩基底双层结构表面失稳分析.力学学报,2017,49(4)：758-762

Huang Chunyang,Tang Shan,Peng Xianghe.Study of surface instability abouthyperelastic film on auxetic substratesunder compression.Chinese JournalofTheoreticaland Applied Mechanics,2017,49(4)：758-762