不同晶向取向的裂纹扩展模拟
开篇：润墨网以专业的文秘视角，为您筛选了一篇不同晶向取向的裂纹扩展模拟范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！
《广西科学杂志》2016年第5期
摘要：
【目的】研究具有不同初始晶向倾角的样品在单轴拉应变作用下的纳观尺度裂纹扩展行为，了解裂纹的生长特征和扩展规律，揭示纳米级裂纹扩展机理及其对材料断裂的影响。【方法】采用晶体相场法观察不同晶向倾角下裂纹的扩展演化图及对应的应力分布图。【结果】当拉应变作用达到临界值时，无预应变的样品裂口开始起裂，并伴随着位错出现。在晶向倾角为0°、5°时，裂口在起裂时，缺口两端裂纹和裂口相连接，裂纹主要是解理脆性断裂模式扩展；在晶向倾角为10°时，裂口向左右两边各发射一个位错，位错在滑移过程中留下一系列空位，空位连通形成裂纹再与主裂口相连，裂纹主要是韧性断裂模式扩展。【结论】不同晶向取向对裂纹的扩展演化有重要影响。
关键词：
微裂纹；晶向取向；位错发射；晶体相场
0引言
【研究意义】材料的变形与破坏主要源自材料的微结构，如空位、位错、晶界与微裂纹等［1］。要想深入研究材料的变形与破坏的特性与机理，还需要开展多尺度分析，把宏观分析与微纳观分析结合起来，在更深层次上找到问题的根源和变形破坏机制［2］。加载的材料通常会表现出可逆行为，当加载除去后，变形就会消失，并且储存能在变形原子键的弹性能得以释放。然而，对于塑性变形和裂纹开裂情况，这一行为是不可逆的。材料在变形中会有能量损耗，其主要来自两个方面：一是断裂产生新表面时所损耗的，二是位错形核和发射位错时所损耗的。因此，对于存在裂纹缺陷的材料变形行为的检验，特别是当位错的形核和运动是由非弹性变形主导的情况，是非常重要的。【前人研究进展】目前已有许多断裂模拟研究，如高英俊等［3］研究韧性材料的微裂纹扩展和连通，毛鸿等［4］研究材料裂纹分叉的机理，刘晓波等［5］对铝裂纹扩展行为进行分子动力学模拟，Song等［6］用相场法研究铁电体的裂纹尖端，Abdolahi等［7］则用相场法模拟铁电体材料的断裂，张跃等［8］进行脆断微裂纹形核的原位观察。在这些模拟研究方法中，既有单尺度分子动力学模拟，也有复杂的多尺度模拟技术。分子动力学模拟有一定的局限性，即它的时间尺度主要适用于原子振动的时间尺度（10－14～10－10s），难以拓展到原子扩散的时间尺度（10－6s），而且由于分子动力学模拟强烈依赖于原子势函数的选择，这就要求所加载的应变力速率远高于实际工业生产的应变力速率（102s－1），达到107～109量级［9］。这使其对微裂纹的扩展模拟与实际存在相当距离。【本研究切入点】晶体相场方法（PFC）［10－14］是基于密度泛函理论、继承传统相场模型优势而建立的，它能够揭示晶体学结构特性以及空间尺度为原子尺度、时间尺度为扩散时间尺度下的结构演化。由于该方法能够很好地描述扩散时间尺度上的微结构演化行为［15－17］，不针对特定材料，因此更适合用于研究微纳米尺度上的裂纹扩展细节［18－19］。【拟解决的关键问题】应用PFC方法研究不同初始晶向倾角条件下，单轴拉伸应变作用的裂纹生长特征和扩展规律，揭示纳米级裂纹扩展机理及其对材料断裂的影响。
1PFC模型与方法
1．1PFC模型
基于PFC模型的液态金属研究如枝晶生长已有许多报道，但更多的是用于固态金属研究。对于固态金属材料，其原子的位置呈规则周期性排列，通过引入周期性相场变量，其局域位置的最大值对应于原子的位置；对于均匀相（液相等）中的原子分布为均匀分布，其值为常量。
1．2样品制备
本研究晶体相用二维三角点阵相表示，其中平衡相为条状相、三角相和液相（图1）。选取原子密度参量为ρ0＝0．49，温度参量γ＝－1．0（图1中B点），缺口处的参数设置为ρ1＝0．79，γ＝－1．0（图1中A点）。应用公式（3）设置单晶体结构。计算模拟区域为1024Δx×512Δy，Δx＝Δy＝π／3，其原子排列方向与y轴夹角θ为0°、5°、10°。在样品的中心位置，设置一半径r＝8的圆形缺口作为初始裂口。由于不涉及材料的物性参数，模拟所用的参数均已无量纲化处理，并将连续空间离散为正方格子，计算时采用周期性边界条件。设置的3组样品参数如表1所示。
1．3拉应变作用的施加
为观察裂纹的萌生和扩展情况，对设置的含有裂口的样品施加拉应变作用。首先，经过5×104步弛豫，得到函缺口的稳态样品（图2，其放大区域给出实验的坐标体系），然后再沿y轴方向（图3所示）施加拉应变。在变形过程中，x方向空间步长保持不变，y方向空间步长随着应变速率在每一时间步长下都有增量d＝εnΔt，其中，ε为无量纲的应变速率，n为施加拉应变的时间步长数，Δt为时间步长，Δy′＝Δy（1＋nεΔt）。
2结果与分析
2．1初始晶向倾角为0°的裂纹生长
当外加应变量达到临界应变值0．138时，裂口开始起裂（图4a），裂纹开裂方向分别向左上［1－21－］方向和右上［21－1－］方向，与y轴的夹角为60°，裂纹长度约为4～5个原子距离，两端裂尖的裂口附近都出现位错结构（图4a中的放大区域所示）。当外加应变量增加到0．162时，裂纹已经由与y轴夹角60°的向左上和右上扩展，转变为向左边水平方向和右边水平方向扩展（图4b），在裂尖的放大细节图上可以观察到裂纹尖端仍然出现位错，但左边位错方向与图4a相比发生转向，顺时针旋转60°，位错的柏氏矢量沿［2－11］方向。外加应变量继续加大，裂纹仍然沿左边水平方向和右边水平方向扩展，裂纹长度增加，左边位错矢量方向不变（图4c）。到扩展后期，外加应变量达到0．210，裂纹总体仍然沿左边水平方向和右边水平方向扩展（图4d），裂纹长度约为模拟区域的三分之二，但尖端可以看到较明显的锯齿状，右边裂尖处的位错与图4a的相比逆时针转动60°。表明裂纹尖端位错滑移方向改变，沿［1－21－］和［21－1－］方向交替变化，形成锯齿边缘结构（图4d中的放大区域所示）。整个裂纹的形状呈解理状，裂尖前端没有出现塑性变形区，也没有空洞形成和位错发射现象，表明裂纹的扩展属于脆性解理扩展。
2．2初始晶向倾角为5°的裂纹生长
开始阶段，在y向的拉应变作用下，裂口发生应力集中；当外加拉应变量达到临界应变值0．138时，裂口向两边发射滑移位错，这是典型的局部塑性变形的表现。由图5a可见，距裂口约2个原子距离的左右两边各出现一个位错，左右两边位错矢量方向相反，柏氏矢量几乎沿水平方向，且此时没有出现与裂口相连接的裂纹。当外加应变量达到0．162时（图5b），左右两边已经形成明显的裂纹并与裂口相连接，由图5b的放大细节图可知，左边的裂尖有一个位错，其柏氏矢量方向与水平方向夹角约60°。此时位错滑移受阻，裂尖应力集中，裂纹前进转向，裂纹边缘呈锯齿状，即左边裂纹扩展方式为扩展－转向－扩展。总体来说，裂纹向左下延伸与水平方向呈约30°，即［2－11］方向；右端裂纹向右边水平方向延伸。当外加应变量达到0．186时（图5c），左边裂纹向左下方扩展，扩展方式改为解理扩展，方向为［2－11］；右边裂纹扩展方式仍然为扩展－钝化－扩展，裂纹形状呈锯齿状，总体上裂纹扩展方向水平向右；右边裂尖也有一个位错，其柏氏矢量方向为向上右偏60°，方向为［101］。到扩展后期（图5d），左边裂纹位错滑行阻力较小，形成解理裂纹，解理裂纹扩展方向为［12－1］，其裂尖位错方向与图5c相同，依旧为向下水平右偏60°，裂纹右边仍然为锯齿状扩展，裂尖的刃位错的半原子面垂直向上。
2．3初始晶向倾角为10°的裂纹生长
由图6a可见，由于裂口应力集中，当外加应变量达到临界应变值0．138时，裂口向左右两侧各发射一个位错，位错滑移方向分别为方向为向上右偏沿［21－1－］方向和向下左偏沿［2－11］方向，位错在滑移过程中沿滑移线诱发生成一系列孤立的空位，空位发展成空洞，然后空洞长大连通，形成裂纹。在空位前端，左边的刃位错滑移方向沿［12－1］方向，右边的刃位错滑移方向沿［1－21－］方向（图6a，b）。当外加拉应变量继续增加时，靠近裂口的空位发展成空洞，空洞相互连接，形成一条与水平方向夹角30°的倾斜裂纹（图6b）。由图6b右下角放大图可见，在裂尖处出现应力集中，在裂尖生成新的位错。随着外加拉应变量的进一步增加，在裂口处向左下和右上延伸的过程中，裂尖应力集中，引起裂尖钝化并使裂纹扩展方向转向；远离裂口处的空位则相互连接成新的裂纹，方向与原主裂纹扩展方向相同。在拉应变作用下，远离主裂纹的空位末端处的位错沿新的方向滑移，出现裂纹分叉（图6c，d），形成二次裂纹扩展。主裂纹钝化转向的尖端也出现位错（图6d），由此可见主裂纹两端钝化转向后开始呈解理状生长，而远离裂口的空位形成的新裂纹也开始钝化转向，该新裂纹的裂尖处也出现位错，位错的方向为左侧的位错向上右偏60°，右侧的位错向下左偏60°。该二次裂纹与主裂纹的扩展方向夹角约为60°。整个区域形成几条不连通的裂纹，表明裂纹主要是韧性断裂模式扩展。
3结论
本研究采用晶体相场法，研究具有不同初始晶向倾角的样品在施加拉应变作用下的裂纹尖端扩展行为，得出如下主要结论：1）晶向角为0°、5°的样品，裂纹主要是解理脆性断裂模式扩展；晶向角为10°的样品，裂纹主要是韧性断裂模式扩展。随着晶向角的变化，裂纹从脆性向韧性裂纹扩展模式转化。2）裂纹作解理脆性扩展时，在裂尖处存在一对位错锁结构，由于位错的柏氏矢量垂直于裂纹扩展方向，这导致在裂尖处存在应力应变集中，原子键呈解理方式断开，裂纹沿光滑直线扩展；裂纹作韧性断裂扩展时，在裂尖前方不远处发生位错滑移运动，且在位错滑移轨道上产生空洞，空洞生长扩大并连通，形成裂纹，此时裂纹边缘呈锯齿形状。
参考文献：
［1］张俊善．材料强度学［M］．哈尔滨：哈尔滨工业大学出版社，2004：83－116．
［2］杨卫．宏微观断裂力学［M］．北京：国防工业出版社，1995．
［3］高英俊，罗志荣，黄礼琳，等．韧性材料的微裂纹扩展和连通的晶体相场模拟［J］．中国有色金属学报，2013，23（7）：1892－1899．
［4］毛鸿，罗志荣，黄世叶，等．材料裂纹扩展分叉机理的晶体相场法研究［J］．广西科学，2015，22（5）：499－505.
［5］刘晓波，徐庆军，刘剑．铝裂纹扩展行为的分子动力学模拟［J］．中国有色金属学报，2014，24（6）：1408－1413．
［8］张跃，王燕斌，褚武扬，等．脆断微裂纹形核的原位观察［J］．科学通报，1994，39（2）：117－120．
［10］黄礼琳，华平，王玉玲，等．凸曲率衬底外延生长界面演化的晶体相场模拟［J］．广西科学，2014，21（3）：214－246．
［13］高英俊，黄礼琳，周文权，等．高温应变下的亚晶界湮没与位错旋转机制的晶体相场模拟［J］．中国科学：技术科学，2015，45（3）：306－321．
［15］高英俊，卢成健，黄礼琳，等．晶界位错运动与位错反应过程的晶体相场模拟［J］．金属学报，2014，50（1）：110－120．
［18］高英俊，罗志荣，邓芊芊，等．韧性材料的微裂纹扩展与分叉的晶体相场模拟［J］．计算物理，2014，31（4）：471－478．
［19］郭刘洋，陈铮，龙建，等．晶体相场法研究应力状态及晶体取向对微裂纹尖端扩展行为的影响［J］．物理学报，2015，64（17）：178102．
［21］高英俊，罗志荣，黄创高，等．晶体相场方法研究二维六角相向正方相结构转变［J］．物理学报，2013，62（5）：050507．
作者:黄礼琳 叶里 胡绪志 黄创高 卢强华 高英俊 单位:广西大学物理科学与工程技术学院