应变梯度理论

1、精品文档应变梯度理论应变梯度理论是近解释材料在微米尺度下的尺寸效应现象而发展起来的一种新理论。Fleek 等6 于 1994 年在细铜丝的扭转实验中观测到微尺度下应变梯度的硬化， 其中直径 12m的无量纲扭转硬化约为直径170 m的三倍。通过对12.5 m、25 m和50 m三种厚度纯镍薄片的弯曲测试， Stolken 和 Evanslv7 于 1998 年发现镍的无量纲弯曲硬化随着薄片厚度的减小而明显增大，然而在拉伸试验中并未发现这种微尺度现象。 Chong 和 Lam8于 1999 年通过压痕实验观察到热固性环氧树脂和热塑性聚碳酸酷的无量纲硬化与应变梯度有关，材料的塑性具有微尺度效应。 M

2、cFarland 和 Colton9J 于 2005 年通过对不同厚度聚丙烯悬臂微梁的弯曲测试， 同样观测到无量纲弯曲刚度随梁厚减小而增大。 与宏观尺度相比，微尺度下结构的力学特性及行为研究主要考虑到以下两个方面(1) 尺度效应。材料不是无限可分。因此材料颗粒的固有属性将影响到微结构的力学特性。(2) 表面和界面效应。 一些在宏观尺度下常被忽略的力和现象， 在微尺度下起着重要的作用 ;而一些在宏观领域作用显著的力和现象， 在微尺度下作用微小，甚至可以忽略。例如， 微尺度下， 与特征尺寸L 的高次方成比例的惯性力、 电磁力 (L3) 等的作用相对减小， 而与尺寸的低次方成比例的粘性力、 弹性力

3、(L2) 、表面张力 (Ll) 、静电力 (L0) 等的作用相对增大。随着尺寸的减小，表面积(L2) 与体积 (L3) 之比相对增大，表面力学和物理效应将起主导作用。理论模型建立(1) 偶应力理论早在一个多世纪前， voigt12 便提出了体力偶和面力偶的概念，并建议构建考虑作用在材料微粒表面或边界上的力偶的连续模型。随后 Cosserat 兄弟 14 根据的假设建立了相关的Cosserat 理论，对应的运动方程中出现了偶应力。直到 20 世纪 60 年代左右，一些学者才开始尝试 Cosserat 理论的改进扩展工作， 他们对 Cosserat 连续体物质点的旋转施加一定约束， 并逐渐发展了一

4、种更为普遍的理论偶应力理论。 相比其它非经典连续介质理论， 偶应力理论是一种相对简单的理论。 如应变梯度理论考虑旋转梯度、 拉伸和膨胀梯度的影响， 而偶应力理论仅考虑了旋转梯度( 与偶应力共轭) 。 Ashby22 指出几何必需位错和统计储存位错是材料的塑性硬化来源，而几何必需位错产生于塑性剪切应变梯度。据此， Fleek 和Hutchinson23 及 Fleek 等 6 在偶应力理论框架上发展了一种应变梯度塑性理论(通常称为CS应变梯度塑性理论)，它是经典的J2形变或J2流动理论的推广。在理论中为了考虑旋转梯度的影响，引入了偶应力，并且服从二阶变形梯度本构率的 Clausius-Duhem

5、 热力学限制条件 24 。 这种理论不仅在模拟裂纹扩展时能消除裂纹尖端的应力奇异性25, 还能成功预测微结构力学行为中的微尺度效应。例如， Fleck 等 6 铜丝的扭转实验中证实了应变梯度硬化的存在，并应用提出的CS应变梯度塑性理论成功解释了这种微尺度现象。经典牛顿力学框架下，连续变形体的材料颗粒仅在力的作用下作平动; 在 TouPin 和 Mindiin 等学者18-21 建立的传统偶应力弹性理论中，材料颗粒不仅在力的作用下作平动，还在力偶的作用下作转动。因此，偶应力理论中的系统能量包括应力对应变和偶应力对旋转形变做的功，其中旋转形变是二阶变形梯度的反对称部分， 含有 8 个独立分量。 对

6、于各向同性线弹性材料而言， 系统本构方程中除了两个经典的拉梅系数外， 还包含两个与材料微结构有关的附加常数。 在上述偶应力理论构建中， 仅用到传统的力和力矩的平衡关系， 对力偶并没有施加约束。Yang 等 28 从引入高阶平衡关系角度出发，提出一种修正偶应力理论。在添加力偶矩平衡关系后， 偶应力张量被约束成对称量， 它对与之共轭张量的曲率张量的对称部分做功， 并与应力对应变做的功一起转变为系统能量。 这种理论下的本构方程仅包含一个附加常数， 从而大大降低了非经典常数的确定难度。 Park 和 Gao29 曾使用这种新理论计算Bemoulli-Euler 微梁的弯曲，发现微梁厚度与材料内察长度相

7、当时， 呈现出明显的尺度效应， 所求得的无量纲弯曲刚度与弯曲实验测量值 28 吻合得较好。(2) 应变梯度理论应变梯度理论的基本思想是通过将高阶应变梯度和 / 或位错密度纳入支配材料行为的本构或演化方程， 来引入尺度对结构或系统的弹、 塑性变形和位错运动等力学行为的影响。 这种理论最早由 Mindlin30 提出，他将弹性体的应变能密度视为应变和它的第一、二阶导数的函数。 同时， 他也给出了一种更常用的仅包含应变和其一阶导数的简化理论， 简化后的附加变形包含了二阶变形梯度的所有18 个独立分量。比较而言，偶应力理论仅包含了二阶变形梯度中的 8 个独立分量，而应变梯度理论是一个完整的二阶梯度理论

8、。 Mindlin 为非经典连续介质力学研究提供了一种新的思路， 后人针对各种应用对其理论进行了改进和扩充。 除了弹性材料外，不少学者致力建立了塑性31-33 、弹塑性 34 、热弹性 135 等材料的应变梯度模型。 例如， 通过使用等效应变的一次和二次拉普拉斯算子表示附加的应变梯度， Aifantis等32建立了应变梯度塑性理论。Fleek等31和Gao等331则发展了另一种基于几何必需位错的应变梯度塑性理论。 Aifantis 为应变梯度理论的发展和应用做出了卓越的贡献。他和他的合作者们建立并逐步发展了模拟物体弹性、 塑性和位错动力行为的各种应变梯度理论，并就相关理论的发展、 应用及数学表

9、述给出了综述 36 。 另外， 黄克智等 37 也在他们的综述性文章中综合介绍了偶应力和应变梯度塑性理论。 除了用于描述位错组态、 材料软化和裂纹尖端附近的变形场等问题外36 ，应变梯度理论也广泛应用在微尺度效应研究中。例如，Aifantis38讨论了应变梯度弹性、 塑性理论在解释不常见微结构的标准尺寸试件或普通微结构的小尺寸试件的扭转和弯曲中的微尺度现象上的能力。 在 Mindlin30 建立的传统应变梯度弹性理论中，附加变形即引入的二阶变形梯度，它包括了 8 个独立分量的反对称部分和 10 个独立分量的对称部分在内的所有18 个独立分量。对于各向同性材料而言，二阶变形梯度对应有七个线性弹性

10、常数， 即两个拉梅系数和五个与材料微结构有关的非经典常数。 应用虚功原理得到的控制方程和边界条件也包含五个附加常数，从而能捕捉到微结构中的尺度效应。后来，Fleck 和 Hutchinson31 ， 39 重新表述了 Mindlin 的应变梯度理论， 他们将二阶变形梯度张量分解成两个独立部分，即拉伸梯度张量和旋转梯度张量。与Mindiin 的工作类似， Fleck和Hutchinson仅使用了传统平衡关系一力和力矩平衡来支配高阶应力行为。受Yang等28的工作启发，Lam等40尝试将新的高阶平衡关系应用在本构关系及控制方程的推导中。在施加附加的力偶矩平衡关系后，Lam等重新定义了高阶应变张量及

11、与之共扼的高阶应力张量，并推导了相应的本构关系和应变能表述。 由于高阶平衡关系的引入， 旋转梯度的反对称部分不出现在变形能中， 与微结构有关的附加材料常数的个数由五个减少到三个。 基于所提出的新理论，Lam等40研究了微悬臂梁的弯曲问题，发现微梁的无量纲刚度与梁厚呈二次方反比关系，这与微梁的弯曲实验观测结果相吻合。(3) 微态理论微态理论是由连续介质力学大师Eringen 建立。在 1964 年， Eringen41 、 Eringen 和Suhubi42 分别提出了简单微流体和简单微弹性体理论， 他们的模型中分别考虑了微流体的局部微运动和微固体的微变形和微旋转， 并推导了对应的基本场方程、

12、边界条件和本构方程。到 1966 年， Eringen43 综合阐述了这类理论， 并将之正式命名为微态连续统力学。 这种理论把材料体看作无数变形物质点的连续集合， 每个物质点都具有有限的尺寸和内部结构。 除了经典的三个平动自由度外， 每个材料物质点还具有独立的拉伸和旋转自由度， 即允许物质 点作刚体运动和发生变形。因此，微态连续体的变形会同时产生宏观应变和微观内部应变。后来的学者发展和拓宽了该理论， 建立了弹粘塑性材料44 、 考虑热存储效应45 、 热磁祸合弹，f体46等一系列理论模型。Chen和Lee47建立了基于微态理论的所有瞬时力学变量和原子坐标及速度的联系， 并通过统计系综平均后得到

13、连续场量， 从而将微态理论和分子动力学结合起来。 微态理论已经应用于微尺度结构与材料的模拟中。 例如， Dinard48 使用基于微态理论的各向异性可压缩塑性模型， 探讨了含孔洞泡沫镍板的尺度效应， 成功预测了孔洞周围的应变集中随孔洞尺寸减小而减弱的变化趋势。微极理论微极理论实质是微态理论的一种特例。 微极理论中， 每个材料物质点除了经典的三个平动自由度外，仅添加了独立的旋转自由度，即意味着物质点可以作刚体运动，但不允许变形。Eringen49l 于 1965 年首次提出微极连续统理论， 并在同年 50 和 1967 年51 分别构建了基于微极理论的流体和弹性体模型。 相一比经典力学理论， 由

14、于附加独立自由度的引入， 微极材料会有旋转惯性矩、 体力偶和表面力偶的作用产生。 微极理论的后续研究重心在于该理论的扩展和推广。 Eringen 、 de Borst 、 Tauchert 等学者在该领域作了大量的工作，提出了用于粘弹性52 、 塑性 53 、 热弹性 54 等材料的微极理论。 微极理论同样可以描述微结构中的尺度效应。例如， McFarland 和 Colton9 使用微极弹性连续理论探讨了微结构对悬臂微梁的弯曲刚度的影响， 计算结果预测到梁的弯曲刚度随梁厚减小而增大， 这与文中实验所观测到的微尺度现象是吻合的。上述理论对每个物质点均引入了附加自由度或高阶变形，与之不同的另一种

15、理论则考虑了表面效应的影响。由于微尺度结构具有急剧增大的表面/ 体积比， 表面效应的重要性显著提高。 为了引入表面效应， Gurtin 和 Murdoch55 提出了一种表面弹性理论， 将表面看作和体不同性质、 无滑移地勃附在体上的二维弹性膜， 表面应力的出现导致了非经典的边界条件， 它和表面经典的弹性方程共同组成场方程。 对于特征尺寸大于 100nm的结构， 表面 / 体积比可以忽略， 因此表面弹性理论目前仅运用在纳米结构中。 Wang 和 Feng56-58 基于 Bemoulli-Euler 梁和 Timoshenko 梁理论，建立了一种考虑表面效应的理论模型，有力预测了纳米梁在振动和屈

16、曲行为中的微尺度效应。实验研究方法微尺度材料的力学性能测试主要包括弹性模量、 泊松比、 残余应力屈服强度、 疲劳强度和断裂强度等参数， 常用的测量方法有拉伸测试法、弯曲测试法、 纳米压入法、 鼓膜法和共振频率测试法等59拉伸试验法(1) 拉伸测试法拉伸测试是测量弹性模量、 泊松比、 屈服强度和断裂强度等参数的最直接的方法。 一般的拉伸测试装备主要包括加载装置、 力传感器、位移传感器、 机械框架和夹具五部分。 载荷和位移是拉伸法测量的主要内容， 所需材料参数通过绘出的包括塑性变形在内的拉伸应力一应变曲线得到。早在1955年，Eisner60就对直径1 m进行了拉伸试验。最近，Chasiotis

17、和Knauss61设计了一种新的拉伸实验装备测得表面微加工多晶硅的弹性模量和拉伸强度， 他们使用原子力显微镜获得变形试样的表面形貌后， 再通过数字图像相关法（DIC） 来确定应变。 在测量 FCc薄膜力学特性的纯拉伸实验中，Espinosa 等 62 发现金、铜和铝试样的屈服应力比试样厚度减小得快，显示出明显的尺度效应。(2) 弯曲测试法弯曲测试法是微尺度材料测试领域比较常用的方法之一。 与拉伸法相比， 弯曲法可以很容易地使用光学显微镜测量，如原子力显微镜（AFM）、扫描隧道显微镜（STM）和力调制显微镜（FMM）等Pearson等631于1957年通过弯曲实验研究了硅丝（直径为20 m）的屈

18、服应力和断裂强度。在新的测量仪器出现后， Espinosa 等 64 使用原子力显微镜、纳米硬度计和 Mirau 型干涉显微镜对金薄膜进行了三点弯曲实验，测得了杨氏模量、屈服应力和残余应力。前文中已提到，已有学者在金属和聚合物材料的弯曲实验中观测到尺度效应 ( 如无量纲弯曲刚度随特征尺寸的减小而增大)7-9 。(3) 纳米压入法纳米压入法具有极高的位移分辨率和加载精度， 通过分析所记录的加载与卸载过程中的载荷一位移曲线来确定材料参数， 主要是材料硬度。 这种方法分为纳米压痕和划痕技术， 其中纳米硬度计是常见的实验装备。 Oliver 和 Pharr6s 使用 Berkovich 压头进行压痕实

19、验并测得熔融石英、钠钙玻璃和单晶铝等六种材料的弹性模量和硬度。 Jardret 等66 项 使用 Berkovich 压头对一些金属和高分子材料展开了划痕实验，指出其中的技术要点。采用扫描隧道显微镜(STM)， Stemashenk 。等67 对单晶钨薄膜所做的压痕实验发现，当压痕深度由10 m 到 1 m 时，材料硬度值急剧增大，表现出强烈的尺度效应。(4) 鼓膜法鼓膜法通常被称之为两轴拉伸试验， 是最早用于研究薄膜力学性能的技术之一。 薄膜凸起高度是主要测量物理量之一， 最早是采用指针仪等机械法进行测量， 随后发展为超声厚度测量，目前则使用精度更高的激光干涉测量68 。 由于金属具有良好的

20、延展性， 早期的鼓膜实验主要集中在对金属薄膜的研究上。 Beam69 最早于 1959 年使用鼓膜技术对金、银薄膜的力学性能进行测量。后来的学者们也逐渐对非金属材料产生兴趣，以矩形薄膜为对象,Tabataa等70测得了 LPCV必晶硅和PECVtM化硅组成的复合薄膜的内应力和弹性模量。(5) 共振频率法微结构的谐振频率与材料的弹性模量、 剪切模量和泊松比等材料特性参数密切相关， 通过检测谐振频率来计算材料参数是一种较早采用的动态测试方法。 该方法采用非接触光学测量技术， 使用激光多普勒干涉仪、 迈克尔逊干涉仪等仪器测量振动位移， 再进一步确定微结构的振动频率。早在1979 年， Peterse

21、n 等71 通过测量悬臂微梁的横向共振频率得到氧化硅薄膜的弹性模量。在压电陶瓷或电磁驱动器激振下，Mazza等72使用激光多普勒干涉仪对馍和铁镍合金微梁进行了振动实验，测得材料的弹性模量、弹性极限和失效应力。 Chen 等 73使用电场驱动的谐振实验装置， 在 Zno 纳米丝的弹性模量测量中观察到尺度效应， 发现当试样直径小于120nm 时，弹性模量随直径的减小而显著增大。(6) Fleck N A, Muller G M ， Ashby M F, etal. Strain gradient Plasticity: theory andexperiment. Acta Metall. Mate

22、r. ， 1994 ， 42:475 一 4877Stolken J S， Evans A G.A microbend test method for measuring the Plasticity length scale. Aeta Metall. Mater. ， 1998 ， 46:5109 一 51158Chong A C M， LamD C C. Strain gradient plasticity effete in indentation hardness of polymers. Journal of Materials Research， 1999， 14(10):410

23、341109Mcfarland A W, Colton J S. Role of material microstructure in Plate stiffnesswith relevanee to microcantilever sensors. J. Micromeeh. Microeng. ， 2005 ， 15:106010610 魏悦广 . 机械微型化所面临的科学难题一尺度效应 . 世界科技研究与发展，2000 ， 22(2):57 一 6111 林谢昭 . 微机电系统的尺度效应及其影响 . 机电产品开发与创新， 2005 ， 18(5):28 一 3012Voigt W.Theo

24、retisehe studien uber die elastieitatsvethaltnissederkrystalle.Abh.5欢迎下载 。精品文档Ges.WissGottingen , 1887 , 34:3 511893, 33:257 25813 Duhem R lliEoric mathematique de la lumiere, Revue des Questions Scientifiques, 1893,33:257-25814 Cosserat Et Cosserat F. fheoric des corps deformablcs. Paris: A* Ilemif

25、inn et Fils, 190915 Eringen A Cf Kafadar C B.微极场论，戴天民译.南京：江苏科学技术出版社, 1982|16 Edelen D G B.非局部场论.戴天民译 南京：江苏科学技术出版社,198117 EringenAC.非局部微极场论.戴天民译.南京：江苏科学技术出版社, 198218 Toupin R A. Elastic materials with couple stresses. Arch+ Ration. Meeh. Anal., 1962, 11:385-41419J Mindlin R D, Tiersten H E Effects of

26、 couple-stresses in linear elasticity Arch. Ration. Meeh. Anal., 1962,】：415y4820J Mindlin R D, Influence of couple-strcsscs on stress concenirations. Exp. Meeh” 1963, 3:1-7f21 Koiter W T. Couple-stresses in the theory of elasticity: I and II, Proc. K. Ned. Akad Wtl. B, 1964.67:17M422 Ashby M F. The

27、defbrinalioii uf plastically iioj * homo Ltieou s alloys, Phil.-Mag., 1970,21:39X2423 Fleck N A, Hutchinson J W. A phenomenological thcor)r for strain gradient cflfccts in plasticity, X Medu Phys. Sol ids, 1993,41: 1825-1857f24| Acharya A. Shawki T G I hcrmodynamic restrict Sens on constitutive equa

28、tions for &eco nd -de form atio n-grad i ent inelastic behavior. J. Meeh. Phys. Solids, 1995, 43: 1751-177225 Xia Z C, Hu(chinsnn J W Crock tip fields in strain gradient plasticity. Journal of the Mechanics and Physics of Sol ids, 1996,10): 1621164826J Nowacki W. Couple slresses in the theory of the

29、rmite last icily. In： Farkus H, Sedov L I, Eds, Proc.lUTAM Symposia, Vienna: Springer 1966t 259278f27j Ristinmaa M, Vecchi M. Use of couple-stress theory in clasto-plasticiij; Comput.Methods AppL Mh. E口gre，* 】9,6,136; 205-22428 Yang F, Chong ACM, Lam 1) C C, cl aL Ccupk: struts bayjd strain gradient

30、 theory for elasticity, Int. J. Solids Struct, 2002,39:273U2745129Park S K (San X L. Bemoullj-Euler beam model based on a modified couple stressiheon. J. Micromech, Mkroeng 2006,16; 2355235930 Mindlin R D. Second gradient of strain and surface tension in linear elasticity. Ini. J, Solids Struct, 196

31、5, k 417-4S831 Fleck N A, Hutchinson J W. Strain gradient plasticity. Adv. App. Meeh, 1997, 33: 295T6132 Aifantis E C On the mkrojrTuctural origin of certain inelastic models J, Mater. Engng/rcchnd 1984,106: 326-33()33 Gao H, Huang 丫 Nix W D, et aL Mcchanism-based stniin gradient plasikky-I. J. Vkxh

32、. Phys. Solids, 1 妙9, 47:】239-126334 Zerviis A, Papanaslasioti P? VardnulaLi I. A finite element displacement 命rmulation for |jraditni eJasiuplasdcity. International Journal fbr Numerical Methods in Engineering, 2001,50(6): 1369138835 Ahnudia G Fiioozbakhsha K. First strain gradient theory of thcrmo

33、elasiicity. In也ma1沁nal Journal of Solids md Siructures.339-34536 Aifajilis E C Update on a class of gradient thcoj ies. Median ks of Materials, 2003, 35： 2528037黄克智，邱信明，妾汉卿.应变梯度理论的新进展(一L偶应力理论和SG理 论.机械强度，1999, 21(2): 81473g： Aifantis E C Slrain gradient interpretation of size effects Inlemational Jou

34、rnal of Fracture, 1W9, 95(13): 299-314391 Fleck N A： Hutchinson J W. A reformulation of strain gradient plasticit)； Mcch. rhys. Solid% 加)1,49: 2245-227140 Lani D C C, Yang E Chung A C M, Cl id. Exprimerits 即id theory in strain gradient eJastic晚上 Meeh. Phys. Solid. 2003. 51; 1477-1 SOS|41 Eringen A C

35、. Simple niicrofluids. Int. L Engrg. ScL, 1964,2(2): 205-21742 Eriiigen A C Suhubi E R. Nonlinear theory of simple microelastic solidsl & H.Int. J Engng. Sd.T 1964,2(2): 11(203,389-404431Ering如 A C. Mechanics of micromorphic maierials. In: Gortlcr H, Ed. Proc 11出ini. Qjng匚 AppL Meeh Berlin: Springer

36、, 1966. 131-13844 Sansour C. A unified concepl of clastic-viscoplastic Cosserat and micromorphic continua J, Phys J V Pmc., 199H.今 34l348451 Eniicn A C Theory cf inicroniorpliic matciiah with inemoiy. Ink L Engng Sci., 1972, 10(7):623 64146 Eringen A C Cnntinunm theory of micromorphic electromagneti

37、c thermoeliKtic solids. Int J. En部g. Sci.t 2003,41(7): 653-665147 Chen Y ?, Lee J D. Connecting molecular dynamics to micromorphic theory (I), InstanUneous and averaged mechanical variables & (II). Balance laws. Physics A, 2003,322: 35979243 Dilhrd I： Forest S, Jenny P. Micrunwrpliic comitiuutD iMuk

38、lling of the dcibrmation and fracture behaviour of nickel fbams European Jouma of Mechanics A/Sdids, 2006,不：526-54949 Eringcn A C. Theory of micropolar continuuni. Dev. Meeh, 1965, 3(1: 23-4050 Efingcn A Ci Theory of nikropolai fluids. J. Math. Meeh, 1996f 16(1)二 IT851 Eritigen A C Theory of micropolar elasticity. In; T teboiiz H, Ed. Fracture. An advanced treatise vol. 11. New York: Academic Preg 1Q6N, 621-72952 Eringen A C. Linear theory of micropolar viscoelasticity Int. J. Engng. Sei” 1967, 5(2)： 191-20455 de Borst R. A generalisation of Jz-flow ikeojy lor polar continua Com