呂云卓

(大連交通大學(xué)材料科學(xué)與工程學(xué)院，遼寧 大連 116028)

?

非晶合金塑性形變基本單元模擬的研究進(jìn)展

呂云卓

(大連交通大學(xué)材料科學(xué)與工程學(xué)院，遼寧 大連 116028)

非晶合金因其獨特的長程無序、短程有序的結(jié)構(gòu)特征而呈現(xiàn)復(fù)雜多樣的宏觀塑性形變行為。目前非晶合金領(lǐng)域已經(jīng)達(dá)成一個共識：非晶合金復(fù)雜多樣的宏觀塑性形變行為具有完全相同的微觀“元過程”，是一種能夠承載剪切變形的局部原子瞬時重排過程，這一瞬時重排過程通常被稱為“流動事件”。由于“流動事件”的尺寸較小，發(fā)生速度較快，現(xiàn)有的測試手段很難同時實現(xiàn)高時間和高空間的分辨率對其進(jìn)行直接觀測，所以目前人們通常采用模擬的方法來對其進(jìn)行研究。主要介紹了目前研究非晶合金“流動事件”的模擬方法，并對這些模擬方法的優(yōu)勢和缺點進(jìn)行了評述，最后簡要展望了非晶合金塑性形變基本單元模擬研究需要重點關(guān)注的幾個問題。

非晶合金；塑性形變；基本單元；模擬方法;流動事件

1 前 言

非晶合金，又稱“金屬玻璃”，是誕生于20世紀(jì)60年代的一種新型金屬材料[1,2]。由于其獨特的長程無序、短程有序的結(jié)構(gòu)特征，使其呈現(xiàn)出一系列十分優(yōu)異的性能，例如高強度、高硬度、高耐磨性、高耐腐蝕性以及良好的軟磁性能等[3-12]。這些優(yōu)異的性能使得非晶合金在航空航天、汽車船舶、裝甲防護(hù)、精密儀器、電力、能源、電子、生物醫(yī)學(xué)等領(lǐng)域都具有廣泛的應(yīng)用前景。除此之外，特殊的原子結(jié)構(gòu)也使非晶合金呈現(xiàn)復(fù)雜多樣的宏觀塑性形變行為[13-19]。如圖1所示，Spaepen教授[20,21]發(fā)現(xiàn)，在高溫、低應(yīng)力以及低應(yīng)變率條件下，非晶合金呈現(xiàn)宏觀的均勻塑性形變；而在低溫、高應(yīng)力或者高應(yīng)變率條件下，其塑性形變僅發(fā)生在納米尺度的剪切帶形變區(qū)域內(nèi)，宏觀塑性形變呈現(xiàn)非均勻化和高度局域化。Schuh 教授[22]后來還發(fā)現(xiàn)在低溫、高應(yīng)變率條件下，高密度剪切帶也可導(dǎo)致非晶合金發(fā)生宏觀均勻塑性形變。

圖1 非晶合金形變圖[20,21]Fig.1 Deformation map for metallic glasses [20,21]

復(fù)雜多樣的宏觀塑性是決定非晶合金未來發(fā)展的關(guān)鍵性能，這是因為非晶合金的生命力取決于以其宏觀塑性為主要考慮因素的大規(guī)模工業(yè)應(yīng)用[23]。正如，人們已成功利用其在過冷液相區(qū)內(nèi)的超塑性加工出高精密、大深寬比的非晶合金零部件，并且實現(xiàn)了其在體育用品、手機配件以及醫(yī)療器械等領(lǐng)域的小范圍應(yīng)用[24,25]。然而，非晶合金的超塑性成型必須謹(jǐn)慎選擇加工溫度和保溫時間，以避免晶化所導(dǎo)致的力學(xué)性能下降。另外，極差的室溫宏觀塑性又嚴(yán)重限制了非晶合金作為結(jié)構(gòu)材料在工業(yè)領(lǐng)域的大規(guī)模應(yīng)用[26-30]。因此，為了能夠有效地指導(dǎo)人們使用和改善非晶合金復(fù)雜多樣的宏觀塑性，進(jìn)而實現(xiàn)其大規(guī)模應(yīng)用，迫切需要澄清其塑性形變機制。然而，目前人們對該問題仍然缺乏較為深入和統(tǒng)一的認(rèn)識，這主要是因為人們對其微觀原子尺度的塑性形變行為仍不清楚。

盡管如此，目前非晶合金領(lǐng)域已經(jīng)達(dá)成一個共識：非晶合金復(fù)雜多樣的宏觀塑性形變行為具有完全相同的微觀“元過程”[31-33]。那就是，其塑性形變的基本單元是一種能夠承載剪切變形的局部原子瞬時重排過程，這一瞬時重排過程通常被稱為“流動事件”[34,35]?；诖?，從單個“流動事件”演化到宏觀塑性形變的可能物理圖像是：在外部熱力作用下，若“流動事件”在非晶合金內(nèi)部的任何位置都等概率發(fā)生，那么該合金的每個部分都參與塑性形變，則該合金在宏觀上表現(xiàn)為均勻塑性形變；若“流動事件”在非晶合金內(nèi)部的一些局域條帶區(qū)域內(nèi)更容易發(fā)生，那么這些條帶區(qū)域?qū)a(chǎn)生較為顯著的塑性形變，進(jìn)而形成剪切帶，則該合金在宏觀上表現(xiàn)為非均勻塑性形變。由此可見，揭示“流動事件”的物理本質(zhì)以及建立其與宏觀塑性形變之間的聯(lián)系，是澄清非晶合金塑性形變機制的關(guān)鍵。

2 非晶塑性形變基本單元的理論模型

為了描述這一“流動事件”的流變過程，人們已經(jīng)提出了許多理論模型。如圖2a所示，最早在該問題上取得突破的是Spaepen教授[20,36,37]提出的自由體積模型(Free Volume Model)，該理論模型本質(zhì)上是將“流動事件”看作是材料內(nèi)部較大自由體積區(qū)域單個原子類似擴散的跳躍。之后，Argon教授[38]、Falk教授[39]和Langer教授[40]提出了剪切轉(zhuǎn)變區(qū)模型(Shear Transformation Zone (STZ) Model)，該理論模型認(rèn)為“流動事件”不是由單個原子的躍遷完成，而是由多個原子參與的局部結(jié)構(gòu)重排過程, 如圖2b。除此之外，還有汪衛(wèi)華院士[41,42]提出的流變單元模型(Flow Unit Model)，這里提到的流變單元是非晶合金基體中原子排列相對疏松的納米尺度區(qū)域，它們具有較低的彈性模量和強度，原子流動性較高，在外界溫度和應(yīng)力的作用下會發(fā)生塑性“流動事件”, 如圖2c。另外還有戴蘭宏研究員與蔣敏強副研究員提出的拉伸轉(zhuǎn)變區(qū)模型(Tension Transformation Zone (TTZ) Model)[43,44]，該模型把“流動事件”描述為主要在拉應(yīng)力作用下局部原子協(xié)同運動導(dǎo)致的原子團簇體積膨脹, 如圖2d。

圖2 描述“流動事件”的理論模型：(a) 自由體積模型[20,36,37]，(b) 剪切轉(zhuǎn)變區(qū)模型[38-40]，(c) 流變單元模型[41,42]，(d) 拉伸轉(zhuǎn)變區(qū)模型[43,44]Fig.2 Theoretical models used to describe the “flow event”: (a) free volume model[20,36,37], (b) shear transformation zone (STZ) model[38-40], (c) flow unit model[41,42] and (d) tension transformation zone (TTZ) model [43,44]

盡管這些理論模型加深了對非晶合金塑性形變基本“流動事件”的理解，但是迄今為止，仍然有許多與這一“流動事件”相關(guān)的基本科學(xué)問題難以解釋或存在爭議。例如，非晶合金塑性形變的基本單元到底是什么? 塑性形變基本單元的形狀、平均尺寸、激活能等基本特征是什么? “流動事件”之間在時間和空間上具有怎樣的相互作用規(guī)律? 溫度、應(yīng)力以及應(yīng)變率等外加條件如何影響“流動事件”的激活與分布? “流動事件”如何跨尺度演化成剪切帶? 等等。若想徹底解決這些有關(guān)“流動事件”的問題，最直接和最可靠的途徑就是利用原位觀測手段，獲得在外力作用下非晶合金原子結(jié)構(gòu)演化的清晰圖像和具體細(xì)節(jié)，從而直接捕捉“流動事件”的發(fā)生。然而，由于非晶合金的原子尺寸較小，運動速度較快，現(xiàn)有的測試手段很難同時實現(xiàn)高時間和高空間的分辨率。比如，先進(jìn)的電子顯微鏡可以在高空間分辨率的條件下觀測到非晶合金的微觀原子結(jié)構(gòu)信息, 但卻不具備足夠高的時間分辨率來捕獲其微觀原子結(jié)構(gòu)隨時間的變化，從而無法直接觀察到動態(tài)“流動事件”，以及其從微觀原子尺度到宏觀塑性形變尺度的演化過程。因此，盡管原位觀測非晶合金的原子結(jié)構(gòu)演化過程是揭示“流動事件”本質(zhì)的有效途徑，但是，目前在實驗方法上仍然面臨巨大的挑戰(zhàn)。

3 非晶合金塑性形變的模擬研究方法

采用模擬的方法可突破上述實驗方法的限制，能夠幫助人們深刻認(rèn)識非晶合金的塑性形變基本“流動事件”。目前，用來模擬研究非晶合金的微觀塑性形變形變的方法可以分為以下兩大類(如圖3)[48,49]：

圖3 模擬研究非晶合金微觀塑性形變行為的方法[48,49]Fig.3 Simulation methods for micro plastic deformation behaviors of metallic glasses[48,49]

3.1 計算機模擬

原子尺度的計算方法是最早被用來模擬研究非晶合金微觀塑性形變行為的計算機模擬方法。原子尺度計算機模擬方法不但可以構(gòu)造出非晶合金的結(jié)構(gòu)，并且能夠精確追蹤每一個原子的運動，從而可以給出有關(guān)塑性形變基本單元的微觀細(xì)節(jié)[45-47]。原子尺度計算機模擬方法又可以分為分子靜力學(xué)方法和分子動力學(xué)方法。其中，分子靜力學(xué)方法不考慮溫度對研究體系的影響，將熱能從體系移除，從而熱起伏對體系沒有影響，那么體系狀態(tài)的改變僅是外部條件所致。利用分子靜力學(xué)方法研究非晶合金微觀塑性形變行為時通常對體系緩慢施加應(yīng)變，非晶體系的應(yīng)力應(yīng)變曲線最初表現(xiàn)出彈性特征，也就是應(yīng)力和應(yīng)變呈線性關(guān)系。然而，經(jīng)過較大的應(yīng)變之后，非晶體系的應(yīng)力會發(fā)生突降，這是非晶體系內(nèi)部結(jié)構(gòu)發(fā)生突然改變所致。當(dāng)發(fā)現(xiàn)應(yīng)力突降現(xiàn)象時，非晶體系內(nèi)部原子的位移分布呈現(xiàn)出四極關(guān)聯(lián)特征[48,49]，如圖4所示。

不同于分子靜力學(xué)方法，分子動力學(xué)方法考慮了熱能對體系的影響，這使得可以利用該方法研究在一定溫度下非晶體系的結(jié)構(gòu)隨時間的演化。例如，Deng等[50]利用分子動力學(xué)方法研究了非晶體系的微觀塑性形變行為。研究發(fā)現(xiàn)，類似于STZ轉(zhuǎn)變的“流動事件”通常會降低體系的自由能，并且是體系因剪切而引起膨脹的主要載體。Falk和Langer[51,52]同樣在其分子動力學(xué)模擬非晶合金的微觀塑性形變行為中發(fā)現(xiàn)STZ轉(zhuǎn)變的“流動事件”。他們還發(fā)現(xiàn)自由體積的產(chǎn)生和湮滅與“流動事件”的發(fā)生密切相關(guān)。

圖4 非晶體系內(nèi)部原子位移分布的四極關(guān)聯(lián)特征[48,49]Fig.4 Quadrupolar feature of atomic motion distribution for amorphous system[48,49]

盡管原子尺度的計算機模擬方法可以給出非晶合金塑性形變微觀“流動事件”的具體細(xì)節(jié)，然而，由于受計算能力的限制，原子尺度的計算方法僅能模擬一元或二元體系的非晶合金，其構(gòu)建的原子結(jié)構(gòu)比真實多元非晶合金的原子結(jié)構(gòu)簡單的多。除此之外，原子尺度的計算方法無法跨尺度研究非晶合金微觀與宏觀塑性形變行為之間的聯(lián)系。例如，非晶合金的原子運動和“流動事件”發(fā)生的時間尺度大約是10-12～10-11s，而剪切帶形成的時間尺度大約是10-5～10-3s；再如，非晶合金的原子和“流動事件”的尺寸大約是10-10～10-9m，遠(yuǎn)小于剪切帶的厚度10-8～10-7m[53]。所以，利用原子尺度的計算模擬方法難以得到從微觀“流動事件”演化成宏觀塑性流動的整個過程。

近年來，Schuh教授[53,54]研究組針對上述問題，采用介觀尺度的連續(xù)介質(zhì)有限元模擬方法，成功地模擬了非晶合金的宏觀塑性形變過程。如圖5所示，他們利用有限元模擬方法模擬了非晶合金剪切帶形核和擴展的整個過程。圖中每一個紅色的圓圈代表一個“流動事件”的激活，這里“流動事件”被設(shè)定為是STZ。可以看到，隨著累積應(yīng)變的不斷增大，“流動事件”的發(fā)生越來越多，并集中發(fā)生在一個條帶區(qū)域，最終形成剪切帶。

盡管介觀尺度的連續(xù)介質(zhì)有限元模擬方法可以成功模擬非晶合金的宏觀塑性形變過程，但是，該模擬方法的理論框架是建立在把剪切轉(zhuǎn)變區(qū)(STZ)認(rèn)定為“流動事件”的基礎(chǔ)之上，而忽略了“流動事件”的具體細(xì)節(jié)。所以，利用介觀尺度的連續(xù)介質(zhì)有限元模擬方法無法得到“流動事件”的微觀信息。

綜上可見，原子尺度計算機模擬方法和觀尺度連續(xù)介質(zhì)有限元模擬方法在模擬研究非晶合金塑性形變基本單元上都有各自的優(yōu)勢，然而，巨大的時間和空間跨度使得利用計算機模擬方法難以獲得從微觀“流動事件”演化到宏觀塑性形變的全尺度物理圖像。

圖5 有限元模擬剪切帶形成過程[53,54]Fig.5 Finite element modeling the formation process of shear band[53,54]

3.2 可視化實物仿真

可視化實物仿真方法是指利用尺寸較大的實物粒子模擬真實的原子，構(gòu)建出與真實合金結(jié)構(gòu)相同的仿真結(jié)構(gòu)，然后利用該仿真結(jié)構(gòu)的形變過程可視化研究實際合金的形變行為。最早利用該方法研究非晶合金塑性形變的是Argon教授[55]，如圖6所示，他以氣泡作為模擬金屬原子的微觀粒子，利用兩種尺寸不同的氣泡搭建出類似非晶合金的二維無序結(jié)構(gòu)，然后通過對該無序排列的氣泡筏施加剪切力來原位觀察單個氣泡的運動過程，并依據(jù)此結(jié)果提出了STZ模型[38]。

以膠體體系為模型，模擬金屬材料的微觀結(jié)構(gòu)、動力學(xué)行為以及力學(xué)行為等，是研究金屬材料基本科學(xué)問題的新穎方法。所謂膠體是一種分散質(zhì)粒子的直徑介于粗分散體系和溶液之間的一類均勻分散體系，它是微米級或納米級的膠體微球粒子分散在液體中。在外場作用下，膠體粒子可以組裝成具有特定有序結(jié)構(gòu)的晶體或者具有無序結(jié)構(gòu)的非晶體。近些年，人們以膠體為模型，已經(jīng)成功地解決和驗證了晶體的形核率、位錯動力學(xué)、點缺陷擴散方式，以及非晶體的力學(xué)行為等一系列以往難以解決的基本科學(xué)問題[56-59]。采用同種電荷不同尺寸的膠體粒子構(gòu)建大尺寸多元膠體玻璃，是可視化模擬研究非晶合金塑性“流動事件”相關(guān)問題的理想模型。這是因為：

(1)同種電荷膠體粒子是模擬金屬原子的理想模型粒子。膠體粒子尺寸在微米級或納米級，盡管膠體粒子的尺寸較大，但是膠體粒子具有與金屬原子完全相同的熱布朗運動[60]。另外，同種電荷膠體粒子之間的相互作用與金屬原子之間的相互作用相同。例如，同種電荷膠體粒子之間在相距較遠(yuǎn)時存在著吸引勢，而隨著膠體粒子間距的減小，粒子之間排斥勢逐漸趨于主導(dǎo)[61]，這與原子之間的相互作用十分相似。此外，膠體固體的摩爾彈性常數(shù)[62]和摩爾融化潛熱[63]在數(shù)值上與原子固體一致，這說明膠體體系中粒子間的作用能與原子體系原子間的作用能處于同一數(shù)量級[64]。

(2)利用不同尺寸的膠體粒子搭建出的多元膠體玻璃與實際多元非晶合金的微觀結(jié)構(gòu)相同。例如，Miracle教授[65,66]曾成功利用膠體粒子堆積出與二元Cu-Zr非晶合金和三元Ca-Mg-Zn 非晶合金具有相同微觀結(jié)構(gòu)的膠體玻璃。利用該膠體玻璃得到的原子團簇形狀、配位數(shù)、堆垛密度以及徑向分布函數(shù)等原子結(jié)構(gòu)信息，與相應(yīng)非晶合金的實驗和計算機模擬結(jié)果一致。

(3)大尺寸膠體玻璃可以建立不同尺度間的貫通，能夠原位觀察從微觀“流動事件”演化到宏觀塑性形變的整個過程。相比于金屬原子，膠體粒子尺寸較大(納米級或微米級)，運動速度較為緩慢，利用光學(xué)顯微鏡可以直接觀察膠體體系的微觀結(jié)構(gòu)變化和追蹤膠體粒子的運動軌跡[67,68]。因此，在大尺寸膠體玻璃中既可以觀察到單個微觀粒子的運動行為，進(jìn)而捕捉“流動事件”，也可以觀察到剪切帶的萌生過程或宏觀均勻塑性流動。例如，最近Spaepen教授和Schall教授[69,70]利用硬球膠體Si粒子構(gòu)建出一元硬球膠體玻璃，并通過該玻璃的剪切實驗成功地觀察到類似于剪切轉(zhuǎn)變區(qū)的局域粒子重排過程，如圖7所示。

上述可視化實物仿真方法可以給出“流動事件”存在的直接實驗證據(jù)，但是，所構(gòu)建的仿真結(jié)構(gòu)與實際合金的原子結(jié)構(gòu)仍有一定差距。例如，氣泡不具有金屬原子的熱布朗運動；硬球膠體粒子之間不存在金屬原子之間的相互吸引作用；仿真結(jié)構(gòu)都是粒子尺寸相同的一元體系，與真實復(fù)雜的非晶合金體系相差較大。盡管如此，若能選用合適的膠體粒子構(gòu)建出與非晶合金原子結(jié)構(gòu)相同的仿真結(jié)構(gòu)，那么，可視化實物仿真方法仍然是研究非晶合金塑性“流動事件”的最直接和最可靠的手段。

圖7 膠體玻璃中觀察到的“剪切轉(zhuǎn)變”事件[69,70]Fig.7 “Shear transformation”event observed in hard-sphere colloidal glass[69,70]

4 結(jié) 語

為了能夠有效地指導(dǎo)人們使用和改善非晶合金復(fù)雜多樣的宏觀塑性，進(jìn)而實現(xiàn)其大規(guī)模應(yīng)用，迫切需要澄清其微觀塑性形變機制。然而由于現(xiàn)有實驗方法的限制，對該問題的認(rèn)識仍然處于盲人摸象的階段。目前主要通過模擬的方法來研究非晶合金的微觀塑性形變行為。原子尺度的和介觀尺度的計算機模擬方法在研究該問題上都有各自的優(yōu)勢，然而，巨大的時間和空間跨度使得利用計算機模擬方法難以獲得從微觀“流動事件”演化到宏觀塑性形變的全尺度物理圖像。可視化實物仿真方法可以給出“流動事件”的跨尺度演化過程，但是，所構(gòu)建的仿真結(jié)構(gòu)與實際非晶合金的原子結(jié)構(gòu)仍有一定差距。今后的模擬研究需要重點關(guān)注以下幾個問題：

對于原子尺度的計算機模擬，①得到的塑性形變基本“流動事件”到底是什么？②非晶合金的原子團簇、結(jié)構(gòu)有序度、自由體積等微觀結(jié)構(gòu)與“流動事件”的發(fā)生具有怎樣的內(nèi)在聯(lián)系？

對于介觀尺度的連續(xù)介質(zhì)有限元計算機模擬，③“流動事件”之間在時間和空間上是否具有關(guān)聯(lián)性？④若“流動事件”之間的時空關(guān)聯(lián)存在，它們在時間和空間上的自組織演化規(guī)律是什么？

對于可視化實物仿真方法，⑤如何構(gòu)建出與實際多元非晶合金結(jié)構(gòu)完全相同的大尺寸仿真結(jié)構(gòu)？⑥在此基礎(chǔ)之上，重點研究從單個“流動事件”演化成剪切帶的整個過程。

References

[1] Greer A L.Science[J], 1995, 267(5206): 1947-1953.

[2] Wang W H, Dong C, Shek C H.MaterialsScienceandEngineering:R:Reports[J], 2004, 44(2): 45-89.

[3] Wang Shenghai(王勝海)，Yang Chuncheng(楊春成), Bian Xiufang(邊秀房).MaterialsReview(材料導(dǎo)報) [J], 2012, 26(1): 88-93.

[4] Lin B, Bian X, Wang P,etal.MaterialsScienceandEngineering:B[J], 2012, 177(1): 92-95.

[5] Liu Fengjuan(劉鳳娟), Pang Shujie(逄淑杰), Li Ran(李 然),etal.RareMetalMaterialsandEngineering(稀有金屬材料與工程)[J], 2008, 37(A04): 733-736.

[6] Inoue A.Engineering[J], 2015, 1(2): 185-191.

[7] Song Wenli(宋溫麗), Song Kaikai(宋凱凱), Liu Zengqian(劉增乾),etal.MaterialsChina(中國材料進(jìn)展)[J], 2014，5:300-311.

[8] Axinte E M, Chirileanu M P I.RecentPatentsonMaterialsScience[J], 2012, 5(3): 213-221.

[9] Chen Weirong(陳偉榮), Wang Yingmin(王英敏), Qian Jianbing(羌建兵),etal.Tribology(摩擦學(xué)學(xué)報)[J], 2003, 23(1): 14-17.

[10]Chen M.NPGAsiaMaterials[J], 2011, 3(9): 82-90.

[11]Huang Caiyun(黃彩云), Sheng Qi(諶 祺), Liu Lin(柳 林).ActaMetallurgicaSinica(金屬學(xué)報)[J], 2010, 46(6): 681-686.

[12]Ding Hongyu(丁紅瑜), Yao Kefu(姚可夫).RareMetalMaterialsandEngineering(稀有金屬材料與工程)[J], 2014, 43(7): 1787-1792.

[13]Dai L H.ShearBandinginBulkMetallicGlasses.AdiabaticShearlocalization:FrontiersandAdvances[M].Waltham: Elsevier,2012: 311-361.

[14]Sun B A, Wang W H.ProgMaterSci[J], 2015, 74:211.

[15]Sun B A, Wang W H.ProgressinMaterialsScience[J], 2015, 74: 211-307.

[16]Zhang Zhefeng(張哲峰), Wu Fufa(伍復(fù)發(fā)), Fan Jitang(范吉堂),etal.ScienceChina(中國科學(xué))[J], 2008, 38(4): 349-372.

[17]Wu Y, Liu W H, Wang X L,etal.AppliedPhysicsLetters[J], 2014, 104(5): 051910.

[18]Qiao J, Jia H, Liaw P K.MaterialsScienceEngineeringRReports[J], 2016, 100:1-69.

[19]Huang Y, Khong J C, Connolley T,etal.InternationalJournalofPlasticity[J], 2014, 60: 87-100.

[20]Spaepen F.ActaMetallurgica[J], 1977, 25(4): 407-415.

[21]Spaepen F.ScriptaMaterialia[J], 2006, 54(3): 363-367.

[22]Schuh C A, Lund A C, Nieh T G.ActaMaterialia[J], 2004, 52(20): 5879-5891.

[23]Axinte E.Materials&Design[J], 2012, 35(35): 518-556.

[24]Inoue A, Takeuchi A.ActaMaterialia[J], 2011, 59(6): 2243-2267.

[25]Shen Jun(沈 軍), Wang Gang(王 剛), Sun Jianfei(孫劍飛),etal.ActaMetallurgicaSinica(金屬學(xué)報)[J], 2004, 40(5): 518-522.

[26]Zhang Z F, He G, Eckert J,etal.PhysicalReviewLetters[J], 2003, 91(4): 045505.

[27]Perepezko J H, Imhoff S D, Chen M W,etal.ProceedingsoftheNationalAcademyofSciences[J], 2014, 111(11): 3938-3942.

[28]Liu Y, Wang Y M, Liu L.ActaMaterialia[J], 2015, 92: 209-219.

[29]Hui Xier(惠希爾)，Chen Guoliang(陳國良).BulkAmorphousAlloy(塊體非晶合金)[M]. Beijing: Chemical Industry Press, 2007.

[30]Kumar G, Neibecker P, Liu Y H,etal.NatureCommunications[J], 2013, 4(2): 1536.

[31]Schuh C A, Hufnagel T C, Ramamurty U.ActaMaterialia[J], 2007, 55(12): 4067-4109.

[32]Lu Z, Jiao W, Wang W H,etal.PhysicalReviewLetters[J], 2014, 113(4): 045501.

[33]Gu L, Xu L, Zhang Q,etal.ScientificReports[J], 2015, 5: 9122.

[34]Wang Z, Sun B A, Bai H Y,etal.NatureCommunications[J], 2014, 5: 5823.

[35]Jiang Minqiang(蔣敏強).MaterialsChina(中國材料進(jìn)展)[J], 2014, 33(5): 257-264.

[36]Spaepen F, Turnbull D.ScriptaMetallurgica[J], 1974, 8(5):563-568.

[37]Spaepen F.DefectsinAmorphousMetals[M]. Amsterdam: North-Holland, 1981: 146-162.

[38]Argon A S.ActaMetallurgica[J], 1979, 27(1): 47-58.

[39]Argon A S, Shi L T.ActaMetallurgica[J], 1983, 31(4): 499-507.

[40]Langer J S.PhysicalReviewE[J], 2004, 70(4): 041502.

[41]Wang Weihua(汪衛(wèi)華).ScienceChina(中國科學(xué))[J], 2014, 44(4): 396-405.

[42]Wang Z, Wen P, Huo L S,etal.AppliedPhysicsLetters[J], 2012, 101(12): 121906.

[43]Jiang M Q, Ling Z, Meng J X,etal.PhilosophicalMagazine[J], 2008, 88(3): 407-426.

[44]Jiang F, Jiang M Q, Wang H F,etal.ActaMaterialia[J], 2011, 59(5): 2057-2068.

[45]Falk M L, Maloney C E.TheEuropeanPhysicalJournalB[J], 2010, 75(4): 405-413.

[46]Mayr S G.PhysicalReviewLetters[J], 2006, 97(19): 195501.

[47]Peng H L, Li M Z, Wang W H.PhysicalReviewLetters[J], 2011, 106(13): 135503.

[48]Tanguy A, Leonforte F, Barrat J L.TheEuropeanPhysicalJournalE[J], 2006, 20(3): 355-364.

[49]Maloney C E, Lematre A.PhysicalReviewE[J], 2006, 74(1): 016118.

[50]Deng D, Argon A S, Yip S.PhilosophicalTransactionsoftheRoyalSocietyofLondonA:Mathematical,PhysicalandEngineeringSciences[J], 1989, 329(1608): 613-640.

[51]Falk M L, Langer J S.PhysicalReviewE[J], 1998, 57(6): 7192.

[52]Conner R D, Li Y, Nix W D,etal.ActaMaterialia[J], 2004, 52(8): 2429-2434.

[53]Homer E R, Schuh C A.ActaMaterialia[J], 2009, 57(9): 2823-2833.

[54]Li L, Homer E R, Schuh C A.ActaMaterialia[J], 2013, 61(9): 3347-3359.

[55]Argon A S, Kuo H Y.MaterialsScienceandEngineering[J], 1979, 39(1): 101-109.

[56]Lu Y, Lu X, Qin Z,etal.JournalofNon-CrystallineSolids[J], 2015, 420: 34-37.

[57]Auer S, Frenkel D.Nature[J], 2001, 409(6823): 1020-1023.

[58]Lu Y, Lu X, Qin Z,etal.PhysicalReviewE[J], 2016, 94(1): 012606.

[59]Nagamanasa K H, Gokhale S, Sood A K,etal.NaturePhysics[J], 2015, 11(5): 403-408.

[60]Lu Y, Li M, Rahman A,etal.ScriptaMaterialia[J], 2014, 90(1): 21-24.

[61]Grier D G.Nature[J], 1998, 393(6686): 621-623.

[62]Crandall R S, Williams R.Science[J], 1977, 198(4314): 293-295.

[63]Williams R, Crandall R S, Wojtowicz P J.PhysicalReviewLetters[J], 1976, 37(6): 348.

[64]Murray C A, Grier D G.AmericanScientist[J], 1995, 83(3): 238-245.

[65]Kramb R C, Ward L T, Jensen K E,etal.ActaMaterialia[J], 2013, 61(6): 2025-2032.

[66]Kramb R C, Ward L T, Jensen K E,etal.ActaMaterialia[J], 2013, 61(18): 6911-6917.

[67]Lu Y, Lu X, Qin Z,etal.ScriptaMaterialia[J], 2016, 115: 104-107.

[68]Zheng Z, Wang F, Han Y.PhysicalReviewLetters[J], 2011, 107(6): 065702.

[69]Schall P, Weitz D A, Spaepen F.Science[J], 2007, 318(5858): 1895-1899.

[70]Jensen K E, Weitz D A, Spaepen F.PhysicalReviewE[J], 2014, 90(4): 042305.

(編輯 吳 琛)

Review of Simulation Research on the Fundamental Unit for the Plastic Deformation of Metallic Glass

LV Yunzhuo

(School of Material Science and Engineering,Dalian Jiaotong University, Dalian 116028, China)

Due to their unique structural characteristics such as long-range disorder and short-range order, metallic glasses exhibit complex and diverse macroscopic plastic deformation behaviors. The field of amorphous alloy has currently reached a consensus: complex and diverse macroscopic plastic deformation behaviors of metallic glasses have the same microscopic "unit process", which is a local rearrangement of atoms undergoing inelastic shear distortion to accommodate plastic deformation, commonly referred to as "flow event". Since the size of the "flow event" is too small and its occurrence is too fast, the existing testing methods, which are difficult to achieve high time and spatial resolutions at the same time, are hard to directly observe the "flow event". Thus, researchers often adopt simulation methods to carry on the research of the "flow event". In this paper, we briefly review these simulation methods as well as their advantages and disadvantages. Finally, we identify several important points that deserve further investigation in the simulation research of the fundamental unit for the plastic deformation of metallic glass.

metallic glass; plastic deformation; fundamental unit; simulation methods; flow event

2016-08-23

國家自然科學(xué)基金資助項目(51401041, 51671042)

呂云卓，男，1985年生，博士，副教授，Email： luyz@djtu.edu.cn

10.7502/j.issn.1674-3962.2017.04.06

TG139.8

A

1674-3962(2017)04-0289-06