袁玉全, 曾祥國, 胡燕飛

(1 四川理工學(xué)院 理學(xué)院, 四川 自貢 643000; 2 四川大學(xué) 建筑與環(huán)境學(xué)院, 四川 成都 610065)

?

含中心微裂紋α-Fe拉伸與疲勞失效機(jī)理的原子模擬

袁玉全1, 曾祥國2, 胡燕飛1

(1 四川理工學(xué)院 理學(xué)院, 四川 自貢 643000; 2 四川大學(xué) 建筑與環(huán)境學(xué)院, 四川 成都 610065)

摘要：采用分子動力學(xué)方法研究了含(010)[101]型中心裂紋的金屬α-Fe在拉伸載荷和疲勞載荷作用下裂紋擴(kuò)展的微觀機(jī)制。研究表明，含有此類裂紋的體心立方金屬鐵在單軸拉伸載荷下，其破壞機(jī)制是一個以裂尖層錯和孿晶變形共同存在、伴有裂紋鈍化，并最終以孿晶變形為主的過程。疲勞失效是以孿晶形變?yōu)橹?，全位錯、不全位錯、層錯共同作用的破壞機(jī)制，同時發(fā)現(xiàn)循環(huán)加載下 [11-1](-121)滑移系中孿晶形變比位錯更容易發(fā)生，而且孿晶的出現(xiàn)，有利于減緩疲勞裂紋的擴(kuò)展速度。含裂紋α-Fe變形與失效的過程，是一個位錯、層錯、孿晶和相變等多種微觀機(jī)制共同作用的復(fù)雜過程。

關(guān)鍵詞:分子動力學(xué); 中心裂紋; α-Fe; 拉伸; 疲勞

PACS: 78.67.Bf

材料中不可避免地存在各種微缺陷如空位、孔洞、夾雜、微裂紋等，它們會嚴(yán)重影響材料的力學(xué)性能，材料的破壞往往從這些微缺陷附近開始萌生和擴(kuò)展。因此，微缺陷附近的變形過程和失效機(jī)理吸引了工程界和學(xué)術(shù)界的廣泛關(guān)注。實驗和理論表明，材料的特性取決于材料的原子結(jié)構(gòu)與微觀結(jié)構(gòu)[1]。盡管微觀實驗儀器掃描隧道顯微鏡和原子力顯微鏡的分辨率都已經(jīng)達(dá)到原子尺寸量級，可以做到直接移動原子水平[2]，但目前要實現(xiàn)對裂紋尖端附近的變形和失效機(jī)理的動態(tài)觀察仍很困難。隨著計算機(jī)水平的發(fā)展，分子動力學(xué)模擬技術(shù)彌補了這一不足，可以展現(xiàn)材料在這一過程區(qū)的微觀動態(tài)過程。

已有大量利用分子動力學(xué)方法對材料變形和失效過程中缺陷附近微結(jié)構(gòu)演化機(jī)理的研究。Argon、Dienes等[3-4]在20世紀(jì)80年代就利用分子動力學(xué)模擬方法進(jìn)行了微裂紋的擴(kuò)展研究，但限于當(dāng)時的計算水平，其重點研究了裂紋的擴(kuò)展以及裂紋區(qū)域的脆-塑轉(zhuǎn)變過程。隨著計算水平的提高，研究體系的規(guī)模逐漸擴(kuò)大，缺陷擴(kuò)展過程中位錯、相變、孿晶等微結(jié)構(gòu)演化機(jī)制以及缺陷形式、載荷方式等對裂紋缺陷擴(kuò)展機(jī)理的影響也逐漸被研究，并且得到了一些重要結(jié)論[5-8]。近年來，分子動力學(xué)研究也擴(kuò)展到材料在疲勞載荷下失效過程的微觀機(jī)制研究領(lǐng)域[9-11]，不過這方面的工作相對較少。

金屬單質(zhì)鐵常溫下以α-Fe形式存在，為體心立方bcc結(jié)構(gòu)，因金屬鐵在結(jié)構(gòu)中有著廣泛的用途，學(xué)者們對其微觀尺度上的失效和破壞機(jī)理做了大量的研究。但目前很多關(guān)于用分子動力學(xué)研究α-Fe的工作都僅僅考慮材料在某一特定晶向下的失效機(jī)理[12-13]，盡管曹莉霞等[14]用F-S多體勢[15]模擬研究了沿幾種典型晶體取向的α-Fe裂紋在外載作用下的形變及相關(guān)缺陷的形成和發(fā)展過程，但是作者沒有考慮α-Fe在疲勞載荷作用下材料變形與失效的微觀機(jī)理。而現(xiàn)實工況中，結(jié)構(gòu)受疲勞載荷作用而失效的情況相當(dāng)普遍，因而這一研究有著重要意義，Uhnáková等[16-18]研究了室溫下帶邊裂紋單晶的α-Fe在Ⅰ、Ⅱ、Ⅲ型疲勞加載下，裂紋擴(kuò)展和位錯發(fā)射機(jī)理，但又沒有考慮不同晶向的影響，對裂紋擴(kuò)展過程中一些重要的現(xiàn)象也未詳盡描述。

為了能更全面地理解材料的微觀斷裂行為,我們系統(tǒng)地建立了幾種典型的α-Fe裂紋分子動力學(xué)模型，采用大規(guī)模分子動力學(xué)模擬程序LAMMPS[19]，全面系統(tǒng)地研究了預(yù)制初裂紋單晶α-Fe拉伸與疲勞破壞失效的微觀機(jī)理，揭示了材料不同晶向?qū)ξ⒘鸭y擴(kuò)展方式、位錯形核與發(fā)射、孿晶、相變等微觀過程的影響機(jī)制。本文主要報道(010)[101]型裂紋的破壞機(jī)制。

1模型和方法

標(biāo)準(zhǔn)的實驗和理論工作都指出，體心立方鐵的解理面是(100) 和(110)面，主要的滑移系是{110}〈111〉和{112}〈111〉[20-21]，為此，建立了4種不同取向的裂紋模型。對于一個特定的裂紋取向,用裂紋所在的平面和裂紋前沿方向這兩個量來進(jìn)行標(biāo)定，即(裂紋所在平面)[裂紋前沿方向]的形式，本系列研究所取的4種裂紋構(gòu)型分別為(010)[101]、(010)[001]、(0-11)[100]、(01-1)[011]，本文將針對(010)[101]型中心微裂紋的破環(huán)機(jī)制進(jìn)行詳細(xì)研究。

圖1為模擬中所用的幾何構(gòu)型,a0為中心裂紋寬度，a0/w=0.05。W、H、t分別表示x、y、z三方向的模型尺寸。位移載荷沿y方向施加，在整個模擬過程中，將頂部和底部各1.25個晶格常數(shù)高度的原子在y方向剛化。x和z方向取為周期性邊界條件，y方向取為非周期性邊界條件。

圖1　模型幾何構(gòu)型

按圖1建立α-Fe的4種分子動力學(xué)裂紋模型，其晶體取向與坐標(biāo)軸對應(yīng)關(guān)系為x:[10-1]；y:[010]；z:[101]。模型尺寸大小為W=32.43 nm，H=22.93 nm，t=2.03 nm，模擬體系原子總數(shù)達(dá)到128 800個。

原子間相互作用采用Mendelev[22]擬合的嵌入原子勢，為了驗證該勢函數(shù)的可靠性，文章計算了α-Fe的一些重要物理量，并與實驗和其他工作所獲得的參數(shù)對比，如表1和表2所示。

表1　晶格常數(shù)和彈性常數(shù)計算與實驗值比較

通過Mendelev[22]勢函數(shù)下的材料晶格常數(shù)和彈性常數(shù)與實驗值的比較，表1表明，除C11以外，其他都符合的很好。表2還表明，與文獻(xiàn)[20，24]一致，(110)和(100)面的表面能更小，(110)表面最小，最有可能成為α-Fe的最優(yōu)解理面，符合實驗規(guī)律[20-21]。該勢函數(shù)在各項物理性能上有極高的可靠性，可以用于本文對α-Fe裂紋的拉伸和疲勞的變形與失效機(jī)制研究工作。

表2　α-Fe不同晶面的表面能

*本工作計算值。

模擬溫度限制在100 K，時間積分步長取為0.001 ps，采用Velocity-Verlet算法求解運動方程。模擬開始前，對系統(tǒng)進(jìn)行能量最小化并在NVE系綜下弛豫，接著在NPT系綜下，沿y方向上下對稱地施加拉伸載荷以及拉伸和壓縮載荷，以模擬產(chǎn)生拉伸荷載和循環(huán)荷載，整個模擬計算采用大規(guī)模分子動力學(xué)模擬程序LAMMPS[19]。為了盡可能接近實際加載情況，本文采用沿拉伸方向線性分配拉伸速度的方法，給模型施加需要的速度，如下圖2a所示，模擬過程每隔200步統(tǒng)計記錄相關(guān)數(shù)據(jù)。

在分子動力學(xué)計算中，基于計算效率的考慮，不宜采用較低的應(yīng)變率，參考其他分子動力學(xué)模擬工作[25]，在模擬過程中給模型施加3×109s-1的應(yīng)變率。加載過程如下：對模型首先進(jìn)行單向拉伸，同時觀察構(gòu)型和應(yīng)力應(yīng)變的變化情況，以得到裂紋擴(kuò)展的臨界應(yīng)變εcri，從而確定最大循環(huán)加載應(yīng)變εmax，循環(huán)加載過程中，取最小應(yīng)變與最大應(yīng)變比例關(guān)系為εmin/εmax=0.75來確定εmin，然后，按圖2b所示加載循環(huán)載荷，以模擬材料疲勞失效過程。

圖2拉伸和循環(huán)加載示意圖

Fig.2Schematic diagram of uniaxial tensile and fatigue loading研究發(fā)現(xiàn)各類裂紋構(gòu)型擴(kuò)展的臨界應(yīng)變和最大應(yīng)變值與裂紋取向密切相關(guān)，差別較大。對于本文報道的(010)[101]裂紋，由于裂紋呈現(xiàn)嚴(yán)重鈍化，所以無法確定其臨界應(yīng)變，多次試驗后取εmax=0.062 98。

2模擬結(jié)果和分析

詳細(xì)分析了預(yù)制(010)[101]中心裂紋單晶α-Fe在拉伸和循環(huán)加載兩種典型外載荷下的變形與失效機(jī)理。

圖3為(010)[101]裂紋構(gòu)型拉伸時的應(yīng)力-應(yīng)變曲線和疲勞載荷下的裂紋長度與循環(huán)周次的關(guān)系，在后面所有的構(gòu)型圖中，我們的觀察平面都取為x-y平面。

2.1單向拉伸下材料變形與失效機(jī)理分析

從圖3a可以看出，拉伸過程中，曲線的變化形式豐富，因此材料破壞過程中必然包含重要的物理機(jī)制，而這些信息，也唯有通過原子尺度上的微觀分析才能得到。

圖3　載荷方向的拉伸應(yīng)力-時間步關(guān)系(a)和裂紋長度與循環(huán)周次關(guān)系(b)

圖4　(010)[101]裂紋拉伸過程典型時間步

如圖4c當(dāng)加載到22 500步時，滑移面上的原子堆垛順序又發(fā)生變化，剪切帶中的原子結(jié)構(gòu)再次發(fā)生轉(zhuǎn)變，逐漸轉(zhuǎn)變回bcc結(jié)構(gòu)，并出現(xiàn)孿晶現(xiàn)象，盡管此時裂紋仍在擴(kuò)展，但應(yīng)力下降速度明顯減慢，接著圖3a應(yīng)力應(yīng)變曲線中出現(xiàn)平臺，此平臺對應(yīng)于加載過程中孿晶的出現(xiàn)，繼續(xù)加載，伴隨著孿晶的“擴(kuò)展”，直到25 500步，剪切帶中的孿晶基本轉(zhuǎn)換完畢，裂紋鈍化加劇，如圖4d。此時圖3a中應(yīng)力曲線才繼續(xù)下降。

模擬結(jié)果與Ján Bosansky等[26]在實驗中觀察的現(xiàn)象一致,即α-Fe在變形與破壞過程中，裂尖層錯和孿晶的形成是共同存在、相互伴隨的過程，而且最終以孿晶變形為主，如圖4e所示。

2.2循環(huán)載荷下的疲勞失效分析

圖5　(010)[101]裂紋在各循環(huán)結(jié)束應(yīng)變

由于在第一個循環(huán)中，就已經(jīng)有變形孿晶出現(xiàn)，所以在以后的循環(huán)過程中，由于孿晶形變和位錯的發(fā)射而使裂紋尖端的應(yīng)力得到一定釋放，裂紋緩慢平穩(wěn)地發(fā)生塑性擴(kuò)展擴(kuò)展，見圖3b。顯然,孿晶的出現(xiàn)有利于減緩疲勞裂紋的開裂擴(kuò)展速度。

3 結(jié)論

本文詳細(xì)分析單晶α-Fe典型的(010)[101]裂紋構(gòu)型在I型單向拉伸和循環(huán)載荷作用下的變形與失效微觀機(jī)理，得到以下主要結(jié)論：

(1)具有(010)[101]裂紋的單晶α-Fe的拉伸破壞機(jī)理是一個以裂尖層錯和孿晶形成共同存在、相互伴隨的過程，而且最終以孿晶變形為主。

(3)具有(010)[101]裂紋的單晶α-Fe單向拉伸載荷作用時形成孿晶前還存在從bcc到fcc再到bcc的相變過程，而循環(huán)加載時α-Fe中直接形成孿晶結(jié)構(gòu)。

因此α-Fe變形與失效的過程，是一個與裂紋方向、加載方式相關(guān)的位錯、層錯、孿晶和相變等多種微觀機(jī)制共同作用的復(fù)雜過程。

參考文獻(xiàn):

[1] SHACKELFORD J F. Introduction to materials science for engineers[M]. New Jersey: Prentice Hall, 2004.

[2] EIGLER D M, SCHWEIZER E K. Positioning single atoms with a scanning tunnelling microscope[J].Nature,1990, 344(6266):524-526.

[3] ARGON A S, YIP S. Molecular dynamics simulation of crack tip processes in alpha-iron and copper[J]. Journal of Applied Physics, 1983, 54(9): 4864-4878.

[4] DIENES G J, PASKIN A. Molecular dynamic simulations of crack propagation[J]. Journal of Physics and Chemistry of Solids, 1987, 48(11): 1015-1033.

[5] ZHANG Y W, WANG T C, TANG Q H. Simulation of nucleation and emission of dislocations by molecular dynamics method[J]. Journal of Applied Physics, 1995, 77(6): 2393-2399.

[6] KUCHEROV L, TADMOR E B. Twin nucleation mechanisms at a crack tip in an hcp material: Molecular simulation[J]. Acta Materialia, 2007, 55(6): 2065-2074.

[7] QI H G, GUO Y F, TANG X Z, et al. Atomistic simulation of the structural evolution in magnesium single crystal under c-axis tension[J]. Acta Metallurgica Sinica(English Letters), 2011,24(6):487-494 .

[8] YUAN Y Q, CHEN H Y, ZENG X G, et al. MD simulation on evolution of micro structure and failure mechanism around interactional voids in pure Al[J]. Applied Mechanics and Materials, 2014, 444/445:183-190.

[9] CHANG W J, FANG T H. Influence of temperature on tensile and fatigue behavior of nano scale copper using molecular dynamics simulation[J]. Journal of Physics and Chemistry of Solids, 2003, 64(8): 1279-1283.

[10] POTIRNICHE G P, HORSTEMEYER M F, JELINEK B, et al. Fatigue damage in nickel and copper single crystals at nano scale[J]. International Journal of Fatigue, 2005, 27(10): 1179-1185.

[11] TANG T, KIM S, JORDON J B, et al. Atomistic simulations of fatigue crack growth and the associated fatigue crack tip stress evolution in magnesium single crystals[J]. Computational Materials Science, 2011, 50(10): 2977-2986.

[12] ANNA M， GLENN E B. Ductile-brittle behavior of (001)[110]nano-cracks in bcc iron[J]. Materials Science and Engineering A , 2004, 387:414-418.

[13] 吳映飛,王崇愚, 郭雅芳.體心立方鐵中裂紋擴(kuò)展的結(jié)構(gòu)演化研究[J]. 自然科學(xué)進(jìn)展，2005, 15(2):206-211.

[14] 曹莉霞,王崇愚. α-Fe 裂紋的分子動力學(xué)研究[J]. 物理學(xué)報，2007, 56(1): 413-422.

[15] FINNIS M W, SINCLAIR J E. A simple N-body potential for transition metals[J]. Philosphical Magazine A, 1984, 50(1): 45-55.

[19] PLIMPTON S J. Fast parallel algorithms for short-range molecular dynamics[J]. Journal of Computational Physics, 1995, 117(1): 1-19.

[20] 馮端.金屬物理學(xué)[M]. 北京：科學(xué)出版社，1999.

[21] HU S Y, LUDWIG M, KIZLER P,et al. Atomistic simulations of deformation and fracture of α-Fe[J]. Modelling and Simulation in Materials Science and Engineering, 1998, 6(5): 567 - 586

[22] MENDELEV M, HAN S, SROLOVITZ D, et al. Development of new interatomic potentials appropriate for crystalline and liquid iron[J]. Philosophical Magazine, 2003, 83(35): 3977-3994.

[23] LEESE J, LORD A E J. Elastic stiffness coefficients of single-crystal iron from room temperature to 500 ℃[J]. Journal of Applied Physics, 1968,39(8):3986-3988.

[25] POTIRNICHEA G P, HORSTEMEYERA M F, JELINEKA B, et al. Fatigue damage in nickel and copper single crystals at nano scale[J]. International Journal of Fatigue, 2005, 27(10):1179-1185.

〔責(zé)任編輯 李博〕

Atomistic simulation of failure mechanism for α-Fe with central crack under uniaxial tensile and fatigue loading

YUAN Yuquan1， ZENG Xiangguo2, HU Yanfei1

(1 School of Science, Sichuan University of Science & Engineering, Zigong 643000, Sichuan, China;2 College of Architecture and Environment, Sichuan University, Chengdu 610065, Sichuan, China)

Abstract:Molecular dynamics simulation is performed for the failure mechanism for α-Fe with (010)[101] centered pre-crack defect under uniaxial tensile and fatigue loading. The results show that the failure mechanism under uniaxial tensile loading is a process with crack blunting, stacking faults and twins formation. While under fatigue loading, the failure mechanism is twins deformation based, coupled with dislocation, partial dislocation, stacking faults.Furthermore the twin deformation occurs easier than the dislocation emission in [11-1](-121) slip system. The study indicate that the deformation and failure process of α-Fe with crack defect is a complicated process, which is the result that many mechanisms act together, such as dislocation, stacking faults, twin and phase transition.

Keywords:molecular dynamics simulation; central crack; α-Fe; tension; fatigue

中圖分類號:O483

文獻(xiàn)標(biāo)志碼:A

基金項目：國防重點基金(B1520132013-1); 四川理工學(xué)院科研項目(2015RC41;2015RC44); 四川省教育廳科研項目(15ZB0207)

收稿日期：2015-08-01

doi:10.15983/j.cnki.jsnu.2016.02.221

文章編號：1672-4291(2016)02-0021-06

第一作者： 袁玉全，男，副教授，博士，研究方向為計算固體力學(xué)和原子分子物理。E-mail：yuquan_yuan@suse.edu.cn