【大学课件】金属材料形变与强化机制的分子动力学研究.ppt
1,金属材料形变与强化机制的分子动力学研究,http:/,2,裂尖形变机制的分子动力学模拟,镍基单晶高温合金沉淀强化机制研究,位错形核和运动(韧脆转变)形变孪晶 相变和新晶粒形核 原子-连续模型的尺寸行为,NiAl,Fe,界面微结构(MD)蠕变和筏化机制(MD+FE+)形状和尺寸以及体积分数的影响;纳米尺寸效应行为(MD+FE+CM+MMM),http:/,3,金属材料强度和裂纹扩展行为与裂尖变形机制及微结构演化过程密切相关,位错形核和运动 形变孪晶 相变和新晶粒形核,铁中裂尖形变机制的分子动力学模拟,http:/,4,Abraham et al.,PNAS,MD,http:/,5,0.1KIC/300步 0.1KIC/800步 0.1KIC/4000步,KI=1.6KIC,Brittle fracture in bcc-iron,(a)010,KI=1.11KIC(b)110,KI=1.27KIC(c)110,KI=1.49KIC(d)010,KI=1.9KIC,Anisotropy for cleavage fracture,Stacking fault(SF),Twinning,http:/,6,*Guo Y-F,Wang C-Y,Zhao D-L.Mater Sci Eng A,2003,349(1-2):29*Y-F Guo,C-Y Wang,Y-S Wang.Phil.Mag.Lett,2004,84(12):763断裂理论基础,范天佑著,科学出版社,2003,p360-366,Brittle to ductile transition,KI=2.2KIC,5K,100K,200K,300K,http:/,7,铁中裂尖形变机制的分子动力学模拟,裂尖相变,相变和新晶粒形核及生长机制,对裂纹扩展影响,A:100;B:110;C:111,http:/,8,F-S N-body potential,Get initial cracks according to continuum mechanics,Relax the crack by MD,Observe the structure evolution of the crack tip,different loads,different types of crack,bcc-Iron,fixed-displacement boundary condition,4000 steps,510-15s,37nm-37nm,Plane strain condition,Length of crack:3.7nm,procedure,http:/,9,Anisotropic linear elastic continuum theory,Displacement field u,Initial crack,http:/,10,G.C.Sih,H.Liebowitz(ed.),in Fracture:An Advanced Treatise(vol.2),Academic,New York,1968,p67.,Crack tip field for plane strain condition,http:/,11,Structure evolution of the crack tip for different types of cracks,A:100 crack B:110 crack C:111 crack,KI=1.9 KIC,the stress-induced phase transformation or recrystallisation can only occurs at KI 1.8 KIC,http:/,12,B:110 crack,C:111 crack,partial dislocations emission on successive planes,bcc-hcp,http:/,13,B:110 crack,C:111 crack,90 steps,110 steps,140 steps,40 steps,80 steps,140 steps,http:/,14,bcc,hcp,Phase transformation,a:2.8665 to 2.733,b:4.0525 to 4.7336,http:/,15,A:100 crack,new grain nucleated with a different orientation from the original bcc crystal(90),boundary between the twin strip and the new grain is a symmetrical tilt boundary with the rotation angle of 70.5,http:/,16,New grain nucleation,new grain grows with the growth of two twin regions,and the lengths of the new grain in and directions are always equal to the thickness of the twin region,and,20 steps,60 steps,140 steps,http:/,17,Distribution of shear stress,A:100 crackB:110 crack,A,B,C,C:111 crack,symmetric,unsymmetric,http:/,18,铁中裂尖形变机制的分子动力学模拟,裂尖孪晶形成和新晶粒形核是低温裂尖变形的两种主要方式；新晶粒形成包括相变和非相变两种形式。伴随裂尖结构演化过程，裂尖应力集中得到松弛。裂尖相变过程中新相为hcp结构。由于hcp结构能量高于bcc结构能量(0.11eV/atom)，相变过程只能在外载很高时才能发生。,界面结构 相变机制 裂纹沿裂尖新晶粒晶界扩展,Y-F Guo,Y-S Wang,D-L Zhao.Acta Materialia(in press),http:/,19,Phase transformation in NiAl,B2,110 crack plane is the only cleavage plane in B2 NiAl(MD,experiments),http:/,20,110 crack,100 crack,B2-L10,http:/,21,110 crack,100 crack,http:/,22,原子-连续模型研究断裂问题时的尺寸相关行为,裂尖应力奇异性,http:/,23,原子弛豫对位移场的影响,010100裂纹,垂直于裂纹表面方向：UY，UX随距离r变化曲线,垂直于裂纹表面方向：dy，dx 随距离r变化曲线,http:/,24,原子区尺寸对裂尖场弛豫结果的影响,如果所取原子区尺寸大于裂尖非线性区半径，体系大小对于弛豫结果没有明显影响,http:/,25,裂纹类型对裂尖场弛豫结果的影响,010101裂纹,http:/,26,原子区尺寸对裂尖区能量分布及裂纹扩展长度的影响,KI=1.6KIC；LY为原子区在Y方向长度,http:/,27,原子区尺寸对裂尖结构演化的影响,LY=29,KI=1.2KIC,KI=1.6KIC,KI=1.9KIC,LY=57,KI=1.9KIC,KI=1.2KIC,KI=1.6KIC,KI=1.1KIC,LY=140,KI=1.6KIC,KI=1.9KIC,http:/,28,原子区尺寸对裂尖结构演化的影响,KI=1.9KIC(a)LY=140；(b)LY=280；(c)LY=360,Y-F Guo,Y-C Gao.Computational Materials Science,2006,36:432-439.,http:/,29,结论,原子弛豫对原子位置影响较大的区域在距离裂尖约100 的范围之内,原子弛豫后的裂尖离散区半径约为100 如果原子区直径大于300，可以运用固定位移边界条件有效模拟bcc-Fe中裂纹的低温脆性解理扩展行为 连续模型在小于300（30nm）的范围内不再适用,原子尺度应力,模型本征,http:/,30,移动边界技术,虚线为移动前的边界，实线为移动后的边界；a为扩展前裂纹长度，a为裂尖移动后的长度,固定位移边界只能保证边界区位移连续，不能保证边界应力连续,用原子-连续模型模拟裂纹稳态扩展：运用应力边界条件运用固定位移边界条件时采用移动边界技术,http:/,31,裂纹稳态扩展,（a）加载500步后的裂尖组态（b）弛豫4000步后的裂尖组态,初始载荷1.2KIC,每弛豫500步根据裂尖移动的距离更新边界,载荷达到1.4KIC后可以观察到裂纹稳态扩展,KI=1.4KIC,http:/,32,镍基单晶高温合金沉淀强化机制研究,强化相,软(基体相)+硬(沉淀相/)无晶界(相界;共格性好)纳米尺寸效应行为,结构特点,http:/,33,强化机制,固溶强化,沉淀强化,/数量(70或更高)、尺寸(0.3-0.5m)、形状(立方)，及本身的固溶强化程度等,筏化现象,/,在高温施加应力条件下，特别是1100K以上，沉淀相会发生定向粗化形成筏状，这种筏化行为直接影响了镍基高温合金的蠕变疲劳寿命，并且断裂面通常沿筏化方向发生(外应力的幅值和方向、晶格错配度、弹性常数差),http:/,34,研究兴趣,界面微结构(MD);蠕变和筏化机制(MD+FE+)形状和尺寸以及体积分数的影响;纳米尺寸效应行为(MD+FE+CM+MMM),Size?Temperature?Stress?,http:/,35,W-P Wu，Y-F Guo,Y-S Wang.Transactions of Nonferrous Metals Society of China(accepted),the critical external tensile stress of about 180Mpa is required for bowing a creep dislocation in the horizontal matrix channel at 950,tensile stress of 300Mpa,http:/,36,THANKS FOR YOUR ATTENTION!,http:/,
