复合材料力学小论文

2011年秋季学期复合材料力学

课程小论文

《微纳米尺度复合材料界面强度的研究》

班级：

姓名：

评分：

微纳米尺度复合材料界面强度的研究

课题背景：

复合材料是一大类新型材料，其强度高、刚度大、质量轻，并具有抗疲劳、减振、耐高温、可设计等一系列优点，近40年来，在航空航天、能源、交通、建筑、机械、信患，生物、医学和体育等工程和部门日益得到广泛的应用。

随着微纳米结构加工技术的发展，人类已经能够在微电子器件、微纳米机电系统中实现多种功能。这些期间与系统通常包含沉积于衬底上的部件，比如薄膜、细线或管状结构、岛状结构或点状结构等。它们的三维尺度非常小，其特征尺寸往往处于微米甚至纳米量级。这些器件的构成材料多样，比如陶瓷、金属间化合物、金属、聚合物等。在这些微纳米尺寸的结构中存在的双相材料界面，由于界面两侧材料物性不同导致的变形不匹配，往往会造成界面处应力集中。因此，界面是微纳米结构极易发生失效的位置之一。界面脱黏、剥离、分层破坏是这类材料的制作及其制成品的使用中的主要失效行为。这些界面失效经常会导致器件报废，直接影响其成品率和使用寿命。因此，为了保证所研制的或正在使用中的微纳米器件或系统的可靠性与耐久性，非常有必要对这些材料及其结构中的界面的强度特性进行研究与分析评价。

论文内容：

本文主要有两大内容：其一，从实验的角度展现在微纳米尺度下复合材料界面端裂纹启裂的现象，并对实验结果作简要分析；其二，

基于分子动力学理论，利用分子动力学模拟软件LAMMPS和Abaqus 模拟微纳米尺度下复合材料界面端裂纹启裂的强度准则。

1 实验

对于微纳米材料与结构而言，比较突出的问题是界面裂纹启裂问题。因为在微纳米系统中，结构或器件的一点失效往往就意味着整个系统的失效。特别是在表面与界面的交汇处，即界面端，由于变形失配引起的应力奇异性使得界面分层裂纹常于此萌生。需要指出的是，随着部件尺寸逐渐缩小，奇异应力主导区的尺寸也在相应减小。对于微纳米部件，根据线弹性应力分析，应力奇异区域一般仅有数十到数百纳米大小。此时，经典断裂力学的概念能否使用值得商榷。

在物理意义上，界面指的是结合材料的结合部，一般不是一个理想的面，存在所谓的界面层、界面相或界面中间层。界面材料的特性不仅复杂，而且具有一定的随机性。但是，界面层厚度非常小，因此，在研究界面的宏观力学性能时，一般将界面结合部理想化为一个面，即，界面。从力学分析与评价的角度来看，这种处理可以避免被界面相得复杂性所困扰，以获取界面、结合材料整体的强度特性参数等。因此，这里将界面定义为材料内的物性间断或不连续处，界面本身不具有厚度，在界面两侧，材料的物性截然不同。

实验装置如图所示，图中的长度量纲为mm，由于采用了夹层结构和较为刚性的悬臂梁，实验中可以有效地抑制薄膜塑性变形与断裂的产生。在远离薄膜材料界面端的悬臂梁左端面处，机械载荷可以被方便地施加并传递到界面端处。

图1 图2

结论：

实验结果观察到载荷-位移曲线在分层破坏之前该曲线基本上是一条直线，在分层破坏发生的瞬间，载荷突然下降。而且分层裂纹一旦启裂，试样旋即破断为两部分，裂纹几乎没有稳态扩展的阶段，表现出明显的脆性破坏特征。用显微镜观察，破断面非常均匀、平整，裂纹启裂于Cu/TiN界面端，并一直沿该界面向前扩展，直至整个试样断为两部分。裂纹尖端存在较弱的奇异性，经计算，界面端应力奇异指数为0.02，如果考虑残余应力的影响之后，界面端附近的应力奇异指数约等于0.18，相应的界面分层启裂应力强度临界值为

0.18

2.8MPa m

，此即Cu/TiN薄膜界面端的分层启裂特性。

2 分子动力学模拟

采用分子动力学模拟方法，可以对结合材料的界面破坏过程进行模拟，从而获得结合材料的界面应力和界面破坏之间的关系。界面破坏可以分为奇异应力场作用下的破坏，和界面应力集中引起的破坏两种。从微观上来说，界面裂纹的产生是一个离散的原子结合对的分离过程，因此，采用分子动力学方法进行数值模拟，应能获得一些关于破坏准则的本质性的认识，有助于我们建立宏观的评价准则。最早的关于裂纹扩展的原子模拟出现于二十世纪七十年代早期，之后被进一

步地应用于脆性材料的断裂、位错的形成和运动、复合材料的断裂特性等方面。

界面破坏包含两个过程，即界面裂纹启裂以及界面裂纹的扩展，本文只模拟界面裂纹的启裂。目前已有一些关于宏观的界面裂纹的破坏准则。本文采用经典的分子动力学和Morse对势，模拟在不同远场载荷作用下，结合材料界面附近原子的运动，从而获得结合材料界面应力和界面破坏之间的关系，进而提出一个结合材料界面破坏的准则。对应于界面破坏的两个过程，本文建立分析模型:界面附近存在初始裂纹(界面启裂)。虽然在分子动力学模拟中我们对界面作了高度的简化，但模拟结果表明，分子动力学模拟仍可给宏观的界面破坏准则提供诸如破坏准则形式等方面的参考。

2.1 分析模型

考虑结合材料受远场应力的作用，为计入不同应力状态的影响，考虑远场斜面力（角度可变）的作用，如图3(a)所示我们只需对相应界面附近的一个微小区域进行分子动力学模拟计算，即图3(a)中的MD(Molecular Dynamics)区，但必须考虑与其他宏观区域的相互作用。由于本文的重点不在分子动力学分析的精度，而在破坏条件的形式方面，故直接采用由有限元分析的结果，作为边界条件施加到MD 区边缘的分子上的方法。分析中材料考虑为铜(Cu)和钨(W)组成的结合材料，为了简单，不考虑界面结合部的缺陷、位错等的影响，并假定钨的晶格常数与铜一致。平面应变有限元分析模型整体尺寸为

nm

18⨯。

2.7

100⨯，分子动力学模拟区域尺寸为nm

200

nm

nm

2.2 分子动力学模拟

两种母材的初始构型均为面心立方晶体(fcc)，在厚度方向上采用周期性边界条件，总共模拟4000个原子。

本文用经典的Morse 对势来描述原子间的相互作用：

02()()()[2]ij ij o a r r a r r ij r D e e ----Φ=- (1)

其中Φ是原子相互作用势,ij r 是原子i 和j 之间的距离,势能参数D 、α、o r 分别表示了原子相互作用的强度、范围和势能最小值的位置。

图3

但是，界面两侧不同种原子之间的相互作用，是与各自的母材内部的势函数不同的，为了方便，相应的势能参数采用以下假定进行计算：

⎪⎪⎪⎩

⎪⎪⎪⎨⎧+=+==1221210122112212112/2ln )()(21)(101αλλαααr D D D ……………………(2) 其中i i i r αλ/2ln 0-=，2,1=i ；下标1、2分别表示母材铜和钨。由于上式