用户登录

你在这里

2018年6月杂志俱乐部:使用中尺度离散位错动力学建模位错介导的塑性

ill.ryu's picture

生病Ryu

机械工程系

艾瑞克·琼森工程与计算机科学学院

德克萨斯大学达拉斯

1.简介

未来技术发展的紧迫挑战和机遇正在显著加速材料研究的进展。在这些努力中,纳米结构材料越来越受到关注,这不仅是因为它们在强度、韧性和热稳定性方面的改进,而且是因为它们通过减轻结构系统的重量来节省能源的潜力。

十多年来,通过在各种时间和空间尺度上对材料进行强大的模拟和建模,已经取得了显著的科学进展。利用这一关键机会需要在中尺度上支撑物质行为。材料的性质和化学反应就是在这个尺度上决定的。基于最近在基本原子长度尺度上开发和仍在研究的知识,中尺度研究将提供预测建模和仿真方面的新兴能力,有可能彻底改变新材料的发展[1,2]。

在这个杂志俱乐部,我们想讨论最近的理解和最新的建模工作,特别是使用离散位错动力学(DDD)建模。我们将跳过DDD的基本概念,DDD可以在一些书籍[3-6]和期刊上发表的优秀综述论文[7-9]中找到。然后,我想提出一些在这个领域的挑战。

2.目前对位错介导塑性的认识

塑性变形的基本机制是位错运动的剪切过程,如图1所示。通过探索位错的详细动态运动,我们可以直接从微观尺度上预测位错集体行为所产生的力学行为。以下是近年来关于位错介导的中尺度可塑性的一些研究结果。

图1。塑性变形为滑移过程[10]

2 - 1。从位错源了解变形机制

从最近的许多微观力学实验来看,“越小越强”现象在小结构中得到了充分的证明[11-13]。在已有的FCC金属尺寸依赖塑性模型中,普遍接受的解释是位错饥饿(DS)模型[14,15]和单臂源(SAS)模型[16-21]。根据DS模型,较小的样品包含较少的位错源,使得位错更容易逸出,并且需要更高的应力来产生新的位错以实现持续的塑性变形。在SAS模型中,截断源的长度在较小的样本中较短,因此需要更高的应力来操作它们。采用表面成核方案的DDD模型可以捕捉FCC微柱中震源驱动的塑性机制,如图2所示。

图2。(a) FCC微柱中直径分别为150 nm(左)和1000 nm(右)的位错微观结构。可以看出,在150nm样品中存在表面成核,在1000nm样品中存在单臂操作。(b) DD结果与SAS和SN理论模型的比较机理[16,22,23]。

在BCC金属中,位错塑性被认为主要由螺旋位错的运动控制,其迁移率低,这归因于非平面的核心结构和高的Peierls势垒[24,25]。此外,从分子动力学(MD)模拟来看,位错增殖是通过表面诱导的交叉滑移过程发生的,而孤立的螺钉位错在图像应力的辅助下可以自我增殖并产生位错源[26]。采用修正的位错迁移率,DDD模型也可以解释尺寸依赖塑性,如图3所示。

图3。直径分别为150 nm(左)和1000 nm(右)的BCC微柱中的位错微观结构。

2 - 2。加载机制取决于位错的微观结构和相应的力学行为

与体积样品相比,微米和亚微米尺度金属样品在各种加载条件下表现出较高的流动应力,包括微柱压缩/拉伸[14,15,27 -34]、金属丝扭转[35-38]、箔弯曲[39-41]以及微和纳米压痕[42-44]。由于传统的连续介质理论无法捕捉观察到的尺寸效应,因此研究不同加载条件下的位错微观结构及其相应的力学响应非常重要。如下图4所示,不同的加载条件在小尺度上产生的微观组织差异较大。基于这一观察结果,它可以限制连续介质方法,其中力学行为主要由简单的张力实验表征。

图4。FCC微柱中的位错微观结构及DDD模型[45]和MD模型[46]中的扭转。

3.挑战与展望

3 - 1。空位/杂质的影响及其与位错的相互作用

基于从间接实验和原子模型中收集的详细知识,中尺度位错动力学(DD)模型可以为研究晶体线缺陷的力学性能如何变化提供独特的机会。然而,目前尚不清楚位错是否与空位/杂质相互作用,遵循一般规则。这项研究有助于在各种类型缺陷的情况下力学性能的基本知识。

3 - 2。晶界内位错相互作用。

人们对纳米晶金属的力学性能的兴趣集中在晶粒和孪晶界在变形机制中的作用。一般来说,GBs可以阻碍位错的运动并导致硬化[47],或者当塑性为位错源限制[48]时,它们可以作为位错的现成来源。虽然已经有人尝试研究不可穿透gb在晶体材料力学行为中的作用[49-52],但仍缺乏遵循一般规律的模型[53-55],具体说明了位错与内部边界之间的各种相互作用模式,如穿透、吸收、解吸[56],如图5所示。

图5。跨孪晶边界界面介导滑移转移反应的原子模型。(a)吸收,(b)解吸,(c)直接透射[56]的原子构型。

3 - 3。位错运动离散描述的统一多尺度模型。

为了研究这些尺度下的塑性行为[57-59],一些研究小组开发了直接耦合离散位错动力学和连续介质有限元法的统一计算模型。这些模型可以同时利用DDD和FE模型;DDD框架考虑位错的动态运动,采用有限元法处理复杂边界和加载条件,并有效地考虑了位错的长程相互作用,如图6所示。

图6。DD与FE[57]的耦合原理图。

在未来的研究中,有必要建立一个通用的模型来描述位错与各种类型的缺陷(如杂质、晶粒和孪晶界以及裂纹)的相互作用,这将提高我们预测多晶固体力学性能的能力。

4.参考文献

[1]从量子到连续体:中尺度科学的机会,中尺度科学报告(2012)。

G.W. Crabtree, J.L. Sarrao,中尺度科学的机遇,MRS Bulletin 37(11)(2012) 1079-1088。

vv . Bulatov, W. Cai,位错的计算机模拟,牛津大学出版社,牛津,2013。

[4] L. Kubin,位错,中尺度模拟和塑性流动。http://public.eblib.com/choice/publicfullrecord.aspx?p=3055181, 2013)。

[5] Z. H.M.离散位错动力学导论,见:S.S. Sansour C.(编),广义连续位错理论,CISM课程与讲座,施普林格,维也纳,2012。

[6]蔡伟文,位错动力学模拟基础,在:T.G. Weinberger C.(主编),多尺度材料的纳米力学建模。《材料科学》施普林格Series in Materials Science施普林格,2016,pp. 53-87。

[7] R.S. Fertig, S.P. Baker,薄膜中位错和强度的模拟:综述,科学进展54(6)(2009)874-908。

周春华,周志强,刘志强,小尺度塑性离散位错动力学模拟,力学学报,58(5)(2010)1565-1577。

[9] D. Weygand, M. Mrovec, T. Hochrainer, P. Gumbsch, bcc金属塑性的多尺度模拟,材料研究年度评论45(1)(2015)369-390。

蔡伟文,刘伟德,刘国强。材料力学与工程学报,2016。

[11] J.R. Greer, J.T.M. De Hosson,小尺寸金属系统的塑性:内在与外在尺寸效应,科学进展56(6)(2011)654-724。

[12] O. Kraft, P.A. Gruber, R. Monig, D. Weygand,有限尺寸的塑性,年度研究报告40(2010)293-317。

[13] M.D. Uchic, P.A. Shade, D.M. Dimiduk,微米尺度的单晶压缩塑性,材料研究,39(2009)361-386。

[14] J.R. Greer, W.D. Nix,通过位错饥饿强化的纳米级金柱,Phys Rev B 73(24)(2006)。

[15] W.D. Nix, J.R. Greer, G. Feng, E.T. Lilleodden,纳米和微米长度尺度下的变形:应变梯度和位错饥饿的影响,薄固体薄膜515(6)(2007)3152-3157。

[16] T.A. Parthasarathy, S.I. Rao, D.M. Dimiduk, M.D. Uchic, D.R. Trinkle,有限样品中位错源长度随机效应对屈服强度的贡献,材料学报56(4)(2007)313-316。

饶思智,D.M. Dimiduk, T.A. Parthasarathy, M.D. Uchic, M. Tang, C. Woodward,三维离散位错模拟的非热机制,学报56(13)(2008)3245-3259。

S.I. Rao, D.M. Dimiduk, M. Tang, T.A. Parthasarathy, M.D. Uchic, C. Woodward,估算微米尺寸单晶中单端位错源的强度,Philos Mag 87(30)(2007) 4777-4794。

[19]吴k.s., A.H.W. Ngan, micropillars中Schmid定律的分解,scriptmater 59(7)(2008) 796-799。

J. a . El-Awady, M. Wen, N.M. Ghoniem,最薄弱环节机制在控制微柱塑性中的作用,J机械物理固体57(1)(2009)32-50。

[21] F. Mompiou, M. Legros, A. Sedlmayr, D.S. Gianola, D. Caillard, O. Kraft,亚微米Al纤维的源增强,材料学报60(3)(2012)977-983。

刘毅,蔡玮,高洪杰,面心立方微柱塑性变形的随机行为,物理学报,95(2015):176-183。

[23] W.D. Nix, s - w。李,表面位错形核控制的微柱塑性,Philos Mag 91(7-9)(2011) 1084-1096。

[24]蔡伟,V.V. Bulatov,位错动力学模拟中的移动规律,材料科学与工程,A 387-389(2004)。

[25] Po G.,崔宇,D. Rivera, D. Cereceda, T.D. Swinburne, J. Marian, N. Ghoniem, bcc金属的唯象位错迁移律,学报119(2016)123-135。

[26] C.R. Weinberger, W. Cai,表面控制的金属微柱位错倍增,中国科学院学报105(38)(2008)14304-14307。

[27] J.R. Greer, W.C. Oliver, W.D. Nix,在没有应变梯度的微米尺度下的金力学性能的尺寸依赖,材料学报53(6)(2005)1821-1830。

[28] M.D. Uchic, D.M. Dimiduk, J.N. flordo, W.D. Nix,样品尺寸对强度和晶体塑性的影响,科学305(5686)(2004)986-989。

[29] D. Kiener, W. Grosinger, G. Dehm, R. Pippan,对尺寸依赖性晶体塑性的进一步理解:微型单晶铜样品的原位拉伸实验,材料学报56(3)(2008)580-592。

[30]吴桂生,颜永华,铝微柱塑性的随机性质,物理学报56(8)(2008)1712-1720。

王志强,王志强,王志强,纯镍微晶的单滑移行为,材料工程学报,53(5)(2005),466 - 466。

D.M. Dimiduk, C. Woodward, R. LeSar, M.D. Uchic,晶体塑性无尺度间歇流动,科学312(5777)(2006)1188-1190。

[33] C.A. Volkert, E.T. Lilleodden,亚微米Au柱变形的尺寸效应,Philos Mag 86(33-35)(2006) 5567-5579。

[34] D. Kiener, C. Motz, T. Schoberl, M. Jenko, G. Dehm,微米尺度下铜力学性能的测定,机械工程学报8(11)(2006)1119-1125。

[35] N.A. Fleck, G.M. Muller, M.F. Ashby, J.W. Hutchinson,应变梯度塑性理论与实验,金属学报42(2)(1994)475-487。

[36] D.J.邓斯坦,B. Ehrler, R. Bossis, S. Joly, K.M.Y. P'ng, A.J. Bushby,弹性极限和应变硬化在扭转,物理学报103(15)(2009)。

刘丹波,何彦敏,张波,张志平,胡培平,丁慧敏,细金属丝扭转过程中尺寸效应的进一步研究,力学学报,41(2013):30-52。

刘丹波,何彦敏,张波,张志平,胡培平,丁慧敏,微米尺度金属丝循环扭转的异常塑性,物理学报,32(12)(2013)。

[J.S. Stolken, A.G. Evans,一种测量塑性长度尺度的微弯曲试验方法,材料学报46(14)(1998)559 -5115.]

[40]王志强,王志强,纳米薄膜的应变梯度效应,物理学报,51(11)(2003)3053-3061。

[41] C. Motz, D. Weygand, J. Senger, P. Gumbsch,三维离散位错动力学模拟与实验的比较,物理学报56(9)(2008)1942-1955。

[42]马Q., D.R. Clarke,银单晶硬度的尺寸相关,J Mater Res 10(4)(1995) 853-863。

[43] W.D. Nix,薄膜的机械性能,金属Trans A 20(11)(1989) 2217-2245。

[44] N.A. Stelmashenko, M.G. Walls, L.M. Brown, Y.V. Milman, W和Mo取向单晶的微压痕研究,金属学报41(10)(1993)2855-2865。

刘毅,蔡玮,高洪杰,面心微柱塑性变形的各向异性,机械工程学报,68(1)(2016):366 - 366。

蔡伟杰,蔡伟杰,金属纳米线在扭转作用下的取向相关塑性:扭转边界的形成和Eshelby扭转,纳米工程10(1)(2010)139-142。

杨波,杨波,H. Vehoff,纳米硬度对压痕尺寸和晶粒尺寸的影响——位错和晶界相互作用的局部检验,材料学报55(3)(2007)849-856。

[48]李志伟,刘志伟,金薄膜纳米压痕行为的微观结构尺度效应,金属学报54(6)(2006)1583-1593。

范洪东,李志强,李志强,双孪生形变对多晶镁硬化响应的影响,材料工程学报,32(2015):366 - 366。

范洪东,范洪东,王晓明,王晓明,超细晶粒镁的晶粒尺寸效应,材料工程,37(2015):357 - 357。

范洪东,李志辉,黄明生,双晶/单晶微柱压缩过程中间歇塑性响应的进一步研究,材料工程,66(10)(2012)813-816。

[52] d . s Balint, V.S. Deshpande, A. Needleman, E. Van der Giessen,多晶体流动强度对晶粒尺寸依赖的离散位错塑性分析,国际力学学报24(12)(2008)2149-2172。

[53] J. Kacher, B.P. Eftink, B. Cui, I.M. Robertson,晶界位错相互作用,固体与材料科学18(4)(2014)227-243。

[54] T.R. Bieler, P. Eisenlohr, C. Zhang, H.J. Phukan, M.A.滑移转移的压边、晶界和界面,固体与材料科学,18(4)(2014)212-226。

[55] D.E. Spearot, M.D. Sangid,从原子模拟对晶界滑移传输的见解,固体与材料科学18(4)(2014)188-195。

[56]朱涛,李俊杰,A. Samanta, H.G. Kim, S. Suresh,纳米结构金属的界面塑性控制应变率敏感性和延展性,物理学报,104(9)(2007)3031-3036。

[57] A. Vattré, B. Devincre, F. Feyel, R. Gatti, S. Groh, O. Jamond, A. Roos,用三维位错动力学和有限元方法模拟晶体塑性:离散-连续模型,J机械物理固体63(2014)491-505。

刘志林,刘小明,庄宗忠,尤晓春,亚微米-纳米尺度下晶体塑性的多尺度计算模型,中国有色金属学报,25(8)(2009)1436-1455。

[59] A. Vattré, B. Devincre, A. Roos,镍基单晶高温合金塑性变形的方向依赖性:离散-连续模型模拟,材料学报58(6)(2010)1938-1951。

附件 大小
形象的图标图1 116.97 KB
形象的图标图2 - 1 103.14 KB
形象的图标图2 - 2 132.2 KB
形象的图标图3 189.36 KB
形象的图标图4 270.96 KB
形象的图标图5 276.26 KB
形象的图标图6 182.69 KB

评论

ill.ryu's picture

我期待听到你们的问题和评论。

Shailendra的照片

亲爱的生病,

感谢您对中尺度塑性建模的深刻概述。正如你所指出的,其中一个突出的挑战是位错-界面相互作用的建模。为这种交互开发一套规则的障碍是什么?随机性在其中起作用吗?在并行建模方法中,使用原子离散化来对界面区域建模,而使用DD/均质连续体来对大块区域建模,这方面是否有努力?

谢谢,

Shailendra

ill.ryu's picture

回答你的第一个问题为这种交互开发一套规则,它主要是由于表征晶界和位错的参数较多。具体来说,晶界的几何形状可以用5个自由度(dof)来描述:3是两个晶粒之间的相对取向差,2是晶界平面法线的方向。同样,位错也需要5个自由度:3个Burgers向量,2个方向。此外,还涉及到材料的固有性质,如核心结构和层错能等。我们还需要考虑入射角的影响。

因此,影响交互机制的因素很多,因此绘制一个简单的规则来确定各种类型的交互是非常具有挑战性的。

从分析和计算两方面对相互作用机理进行了许多研究,但仍缺乏一个普遍的规律。我希望特性转化也会发挥作用,但我不太了解它的细节。关于随机性的其他意见/建议非常感谢。

其次,我同意模型界面的并行建模方法是处理复杂相互作用的好方法,因为连续介质方法将有助于考虑大变形(包括旋转),也会产生准确的远场应力计算。然而,它仍然需要一个“规则”,我相信,这个规则只能来自具有适当势能的原子建模。

李晓燕的照片

亲爱的生病,

非常感谢你发布如此令人印象深刻的话题。如图4所示,DD模拟所揭示的变形机制与MD模拟所揭示的变形机制非常相似。但DD模拟与MD模拟在长度尺度上存在明显差异。有什么方法可以将MD和DD连接起来吗?

最好的

小燕

ill.ryu's picture

谢谢你的意见。DD与MD的位错结构非常相似,并且基于max得到简单的形核,这是非常有趣的。解析了不同取向微柱的剪切应力。然而,DD的位错网络与MD的结果并不完全相同。例如,在<100>定向的岩柱中,MD显示为矩形网络,而DD显示为十字形位错。我认为这两种模型都有局限性,无法得到“真正的”结构。也就是说,MD中很高的应变率会使位错活动局限在单个滑移面上,DD中位错的微观结构是由形核规则形成的。我认为需要进行原位TEM实验来观察真正的位错结构,但在这种规模下进行扭转实验也具有挑战性。我期待看到各向异性位错微观结构的实验证明。

对于MD和DD不同的时间尺度,通过改进的加载条件(即所谓的截止塑性应变率法),DD结果可以得到与稳态实验非常相似的应力-应变关系,而MD对高应变率有内在的限制。为了刻画个体位错迁移和相互作用机制,DD模型需要有原子输入。因此,MD和DD之间的伸缩桥接将采用分层耦合的方式处理,而不是由于时间尺度的不同而采用并行的方式。

哈罗德·s·帕克的照片

亲爱的生病了:

谢谢你提出这个有趣的日记俱乐部。既然你提到了MD的众所周知的时间尺度问题,我想提到一些最近在这一领域所做的关于金属塑性的工作。具体来说,我们最近展示了,使用势能景观探测技术,我们可以捕捉双晶金属纳米线中实验观察到的速率依赖的延性到脆性转变:

https://pubs.acs.org/doi/abs/10.1021/acs.nanolett.7b04972

同时还预测了单晶金属纳米线在MD无法达到的时间尺度上出现的新的超塑性蠕变行为:

https://pubs.acs.org/doi/abs/10.1021/acsnano.8b02199

也有其他小组在使用类似的时间尺度桥接技术对不同材料系统的预测建模方面做了出色的相关工作。例如Sid Yip和非晶态固体的合作者(http://www.pnas.org/content/114/52/13631),以及Pradeep Sharma研究各种材料系统(电池、金属等)。https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.1.055401).

ill.ryu's picture

亲爱的哈罗德

感谢您的评论和分享非常相关的论文。MD中的应变率是一个长期存在的挑战,许多努力已经解决了这一局限性。看到您最近的工作和其他人的工作非常有趣,这些工作可以在MD建模中实现物理上合理的时间持续时间。

另一种有趣的解决MD中时间步长的方法是由d.e. Shaw Research (DESRES)设计和建造的用于生物分子模拟的特殊用途超级计算机。无需开发新的时间步长算法,他们就可以进行长时间跨度的MD模拟(https://www.nature.com/articles/nphys3553),它使我们能够研究分子的详细运动,而不丢失原子的“抖动和摆动”运动,理查德·费曼在他开创性的1963年物理学讲座中说。

李晓燕的照片

亲爱的生病,

非常感谢你的详细解释,这为多尺度建模提供了重要的见解。你提出的分层耦合可能有效地连接DD和MD,这将促进从原子水平到位错水平对塑性的理解。

订阅“2018年6月期刊俱乐部:使用中尺度离散位错动力学建模位错介导的塑性”的评论

更多的评论

辛迪加

订阅辛迪加