基于模擬分析探究高強(qiáng)度稀土鎂合金強(qiáng)化相形成與原子排列機(jī)制

基于模擬分析探究高強(qiáng)度稀土鎂合金強(qiáng)化相形成與原子排列機(jī)制一、引言1.1研究背景與意義在現(xiàn)代工業(yè)不斷追求高性能材料的大背景下，高強(qiáng)度稀土鎂合金憑借其獨(dú)特的優(yōu)勢，成為材料科學(xué)領(lǐng)域的研究焦點(diǎn)。鎂合金作為實(shí)際應(yīng)用中最輕的金屬結(jié)構(gòu)材料，具有密度低、比強(qiáng)度和比剛度高、阻尼性及切削加工性能好、易回收等諸多優(yōu)點(diǎn)，在航天、軍工、汽車、電子通訊等領(lǐng)域展現(xiàn)出巨大的應(yīng)用潛力，被譽(yù)為“21世紀(jì)的綠色工程材料”。然而，大部分常規(guī)鎂合金的力學(xué)性能，尤其是高溫力學(xué)性能較差，嚴(yán)重限制了其在一些對材料性能要求苛刻領(lǐng)域的廣泛應(yīng)用。為了突破鎂合金性能的瓶頸，科研人員嘗試通過添加合金化元素來改善其力學(xué)性能，其中添加稀土元素成為一種極為有效的手段。稀土元素因其獨(dú)特的核外電子排布，表現(xiàn)出特殊的物理和化學(xué)性質(zhì)，對鎂合金具有“凈化”“細(xì)化”“強(qiáng)化”“合金化”的四重作用。當(dāng)稀土元素加入鎂合金后，一方面，它能與鎂合金中的雜質(zhì)元素如鐵、鎳等形成高熔點(diǎn)的化合物，這些化合物在合金凝固過程中會(huì)浮出液面或沉淀在底部，從而凈化合金熔體，減少雜質(zhì)對合金性能的不利影響；另一方面，稀土元素可以在鎂合金凝固過程中細(xì)化晶粒，通過增加形核核心，抑制晶粒長大，使合金組織更加均勻細(xì)小，進(jìn)而顯著提高合金的強(qiáng)度和韌性。此外，稀土元素還能與鎂或其他合金化元素在合金凝固過程中形成穩(wěn)定的金屬間化合物，這些含稀土的金屬間化合物一般具有高熔點(diǎn)、高熱穩(wěn)定性等特點(diǎn)，它們呈細(xì)小化合物粒子彌散分布于晶界和晶內(nèi)，在高溫下可以釘扎晶界，抑制晶界滑移，同時(shí)阻礙位錯(cuò)運(yùn)動(dòng)，強(qiáng)化合金基體，提高合金的高溫強(qiáng)度和蠕變性能。基于稀土元素對鎂合金性能的顯著提升作用，Mg-RE系合金應(yīng)運(yùn)而生，并因其優(yōu)異的高溫拉伸性能、抗蠕變性能及良好的塑性成形能力而備受青睞，被認(rèn)為是最具應(yīng)用前景的高溫高強(qiáng)合金體系。目前，稀土鎂合金已廣泛應(yīng)用于航空航天領(lǐng)域，用于制造發(fā)動(dòng)機(jī)部件、衛(wèi)星支架、橫梁等結(jié)構(gòu)件，以滿足高溫、高強(qiáng)度的使用要求；在汽車工業(yè)中，可用于制造發(fā)動(dòng)機(jī)罩、排氣系統(tǒng)等高溫環(huán)境下的結(jié)構(gòu)件，有助于實(shí)現(xiàn)汽車的輕量化，降低能耗和排放；在軍工領(lǐng)域，因其輕質(zhì)、高強(qiáng)、耐熱的綜合性能特點(diǎn)，被廣泛應(yīng)用于導(dǎo)彈、坦克等武器裝備的制造中。盡管稀土鎂合金在實(shí)際應(yīng)用中取得了一定的成果，但目前對于其強(qiáng)化相形成及原子排列的研究仍存在許多不足。強(qiáng)化相作為決定合金力學(xué)性能的關(guān)鍵因素，其形成過程涉及復(fù)雜的物理化學(xué)變化，受到合金成分、制備工藝、熱處理?xiàng)l件等多種因素的交互影響。例如，在不同的合金成分體系中，稀土元素與其他元素之間的化學(xué)反應(yīng)和擴(kuò)散過程不同，導(dǎo)致強(qiáng)化相的種類、數(shù)量、尺寸和分布存在差異；制備工藝如鑄造、鍛造、軋制等會(huì)影響合金的凝固過程和組織結(jié)構(gòu)，進(jìn)而對強(qiáng)化相的形成產(chǎn)生影響；熱處理工藝中的溫度、時(shí)間和冷卻速度等參數(shù)的變化，也會(huì)改變強(qiáng)化相的析出行為和原子排列方式。此外，原子排列方式不僅決定了強(qiáng)化相的晶體結(jié)構(gòu)和性能，還與合金基體之間存在著復(fù)雜的界面關(guān)系，這種關(guān)系對合金的力學(xué)性能、物理性能和化學(xué)性能有著深遠(yuǎn)的影響。然而，目前對于這些復(fù)雜關(guān)系和影響機(jī)制的認(rèn)識(shí)還不夠深入，實(shí)驗(yàn)研究由于受到技術(shù)手段和觀測尺度的限制，難以全面、準(zhǔn)確地揭示強(qiáng)化相形成及原子排列的微觀過程和內(nèi)在規(guī)律。從理論研究層面來看，深入探究高強(qiáng)度稀土鎂合金中強(qiáng)化相形成及原子排列的機(jī)制，有助于豐富和完善材料科學(xué)的基礎(chǔ)理論。通過建立合理的理論模型，能夠從原子尺度上理解合金元素之間的相互作用、原子擴(kuò)散規(guī)律以及晶體結(jié)構(gòu)的演變過程，為材料的設(shè)計(jì)和優(yōu)化提供堅(jiān)實(shí)的理論依據(jù)。從實(shí)際應(yīng)用角度出發(fā)，清晰掌握強(qiáng)化相形成及原子排列與合金性能之間的關(guān)系，能夠指導(dǎo)科研人員通過調(diào)整合金成分和制備工藝，精確控制強(qiáng)化相的形成和原子排列，從而開發(fā)出性能更加優(yōu)異的稀土鎂合金材料，滿足航空航天、汽車、電子等高端領(lǐng)域?qū)Σ牧喜粩嗵岣叩男阅苄枨螅苿?dòng)相關(guān)產(chǎn)業(yè)的技術(shù)進(jìn)步和發(fā)展，具有重要的科學(xué)意義和工程應(yīng)用價(jià)值。1.2國內(nèi)外研究現(xiàn)狀在高強(qiáng)度稀土鎂合金的研究領(lǐng)域，國內(nèi)外學(xué)者圍繞強(qiáng)化相形成及原子排列展開了廣泛而深入的探索，取得了一系列具有重要價(jià)值的成果。國外方面，早期的研究重點(diǎn)主要集中在稀土元素對鎂合金微觀組織和力學(xué)性能的影響。例如，德國學(xué)者在20世紀(jì)20年代就開始了Mg-MM合金的開發(fā)工作，并將其應(yīng)用于飛機(jī)發(fā)動(dòng)機(jī)部件，開啟了稀土鎂合金的研究序幕。后續(xù)研究中，科研人員發(fā)現(xiàn)稀土元素能夠與鎂合金中的雜質(zhì)元素形成高熔點(diǎn)化合物，從而凈化合金熔體，減少雜質(zhì)對合金性能的不利影響。同時(shí)，稀土元素在鎂合金凝固過程中細(xì)化晶粒的作用也被深入研究，揭示了其通過增加形核核心、抑制晶粒長大來改善合金組織均勻性和性能的機(jī)制。在強(qiáng)化相方面，對Mg-RE系合金的時(shí)效析出強(qiáng)化相進(jìn)行了大量研究，確定了如β′′(DO19)、β′(cbco)等主要時(shí)效析出強(qiáng)化相，明確了其化學(xué)組成分別為Mg3RE和Mg15RE3。通過透射電子顯微鏡（TEM）、掃描電子顯微鏡（SEM）等微觀分析技術(shù)，對強(qiáng)化相的形貌、尺寸和分布進(jìn)行了細(xì)致觀察，發(fā)現(xiàn)這些強(qiáng)化相呈細(xì)小化合物粒子彌散分布于晶界和晶內(nèi)，在高溫下能夠釘扎晶界、阻礙位錯(cuò)運(yùn)動(dòng)，有效提高合金的強(qiáng)度和高溫性能。近年來，隨著計(jì)算機(jī)技術(shù)和計(jì)算材料學(xué)的快速發(fā)展，國外開始運(yùn)用第一性原理計(jì)算、分子動(dòng)力學(xué)模擬等方法，從原子尺度研究強(qiáng)化相的形成過程和原子排列方式，為深入理解合金化機(jī)制提供了新的視角。例如，通過第一性原理計(jì)算研究稀土元素與鎂及其他合金元素之間的相互作用能，預(yù)測強(qiáng)化相的形成可能性和穩(wěn)定性；利用分子動(dòng)力學(xué)模擬合金凝固過程中原子的擴(kuò)散和聚集行為，揭示強(qiáng)化相的形核和生長機(jī)制。國內(nèi)在高強(qiáng)度稀土鎂合金研究方面起步相對較晚，但發(fā)展迅速。在合金體系開發(fā)上，基于國外已有的研究成果，結(jié)合我國豐富的稀土資源優(yōu)勢，開發(fā)出了一系列具有自主知識(shí)產(chǎn)權(quán)的稀土鎂合金體系，如在Mg-Gd體系上形成的Mg-Gd-Y、Mg-Gd-Er等系列合金，以及在Mg-Y體系上形成的Mg-Y-Gd、Mg-Y-Nd等系列合金。在強(qiáng)化相研究中，采用多種微觀分析手段相結(jié)合的方法，對強(qiáng)化相的種類、形成條件、晶體結(jié)構(gòu)和強(qiáng)化機(jī)制進(jìn)行了系統(tǒng)研究。通過實(shí)驗(yàn)發(fā)現(xiàn)，在不同的合金成分和熱處理?xiàng)l件下，強(qiáng)化相的析出行為和原子排列方式會(huì)發(fā)生顯著變化。例如，研究發(fā)現(xiàn)通過調(diào)整合金成分和熱處理工藝，可以控制Mg-Gd-Y合金中β′′相和β′相的析出數(shù)量、尺寸和分布，從而優(yōu)化合金的力學(xué)性能。在理論研究方面，國內(nèi)學(xué)者積極開展計(jì)算材料學(xué)研究，利用模擬軟件深入探究強(qiáng)化相形成及原子排列的微觀機(jī)制，為合金的設(shè)計(jì)和優(yōu)化提供理論支持。同時(shí)，在實(shí)際應(yīng)用研究中，國內(nèi)將稀土鎂合金廣泛應(yīng)用于航空航天、汽車、電子等領(lǐng)域，通過產(chǎn)學(xué)研合作，不斷推動(dòng)稀土鎂合金的產(chǎn)業(yè)化進(jìn)程。盡管國內(nèi)外在高強(qiáng)度稀土鎂合金強(qiáng)化相形成及原子排列研究方面取得了諸多成果，但仍存在一些不足之處。從實(shí)驗(yàn)研究來看，目前對于強(qiáng)化相在復(fù)雜制備工藝和服役條件下的形成過程和演化規(guī)律研究還不夠全面。例如，在多種制備工藝協(xié)同作用下，如鑄造與鍛造、軋制相結(jié)合的復(fù)合工藝，強(qiáng)化相的形成機(jī)制尚未完全明確；在高溫、高壓、復(fù)雜應(yīng)力等多場耦合的服役環(huán)境下，強(qiáng)化相的穩(wěn)定性和原子排列的變化規(guī)律缺乏深入研究。從理論模擬角度，雖然計(jì)算材料學(xué)的發(fā)展為研究提供了有力工具，但現(xiàn)有的模擬方法和模型仍存在一定局限性。一方面，模擬過程中對一些復(fù)雜因素的考慮不夠全面，如合金元素之間的交互作用、晶體缺陷對原子擴(kuò)散和強(qiáng)化相形成的影響等；另一方面，模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的對比驗(yàn)證還不夠充分，導(dǎo)致模擬結(jié)果的可靠性和準(zhǔn)確性有待進(jìn)一步提高。此外，目前對于強(qiáng)化相形成及原子排列與合金宏觀性能之間的定量關(guān)系研究還不夠深入，難以實(shí)現(xiàn)通過精確控制微觀結(jié)構(gòu)來精準(zhǔn)調(diào)控合金性能的目標(biāo)，這在一定程度上限制了高性能稀土鎂合金的開發(fā)和應(yīng)用。1.3研究內(nèi)容與方法本研究旨在深入探究高強(qiáng)度稀土鎂合金中強(qiáng)化相形成及原子排列的微觀機(jī)制，通過理論模擬與實(shí)驗(yàn)驗(yàn)證相結(jié)合的方式，為高性能稀土鎂合金的開發(fā)和應(yīng)用提供堅(jiān)實(shí)的理論基礎(chǔ)和技術(shù)支持。具體研究內(nèi)容如下：合金成分對強(qiáng)化相形成的影響：選取典型的稀土鎂合金體系，如Mg-Gd-Y、Mg-Y-Nd等，通過改變合金中稀土元素（Gd、Y、Nd等）以及其他合金化元素（如Zn、Zr等）的含量，利用第一性原理計(jì)算和分子動(dòng)力學(xué)模擬方法，研究合金成分變化對強(qiáng)化相形成熱力學(xué)和動(dòng)力學(xué)的影響。計(jì)算不同成分下合金體系中各種可能強(qiáng)化相的形成能、生成焓等熱力學(xué)參數(shù)，分析強(qiáng)化相形成的可能性和穩(wěn)定性；模擬合金凝固過程中原子的擴(kuò)散和聚集行為，揭示合金成分對強(qiáng)化相形核和生長速率的影響規(guī)律，明確合金成分與強(qiáng)化相種類、數(shù)量之間的定量關(guān)系。制備工藝對強(qiáng)化相形成及原子排列的影響：考慮鑄造、鍛造、軋制等常見制備工藝，采用數(shù)值模擬方法研究不同制備工藝參數(shù)（如鑄造溫度、冷卻速度，鍛造的變形量、變形溫度，軋制的道次、壓下量等）對合金凝固過程、晶體生長方式以及原子擴(kuò)散路徑的影響，進(jìn)而分析這些因素如何作用于強(qiáng)化相的形成過程和原子排列方式。結(jié)合實(shí)驗(yàn)研究，利用掃描電子顯微鏡（SEM）、透射電子顯微鏡（TEM）等微觀分析手段，觀察不同制備工藝下合金中強(qiáng)化相的形貌、尺寸、分布以及原子排列特征，驗(yàn)證模擬結(jié)果的準(zhǔn)確性，建立制備工藝與強(qiáng)化相微觀結(jié)構(gòu)之間的內(nèi)在聯(lián)系。強(qiáng)化相原子排列與合金性能的關(guān)系：基于晶體結(jié)構(gòu)理論和位錯(cuò)理論，通過第一性原理計(jì)算和有限元模擬，研究強(qiáng)化相不同原子排列方式（如晶體結(jié)構(gòu)類型、晶格常數(shù)、原子占位等）對合金力學(xué)性能（如強(qiáng)度、硬度、韌性）、物理性能（如熱膨脹系數(shù)、彈性模量）和化學(xué)性能（如耐腐蝕性）的影響機(jī)制。計(jì)算強(qiáng)化相原子排列與合金基體之間的界面能、界面結(jié)合強(qiáng)度等參數(shù)，分析界面結(jié)構(gòu)對合金性能的影響；模擬在外加載荷、溫度變化、腐蝕介質(zhì)等條件下，強(qiáng)化相原子排列的穩(wěn)定性以及合金的性能演變過程，建立強(qiáng)化相原子排列與合金宏觀性能之間的定量關(guān)系模型，為通過調(diào)控強(qiáng)化相原子排列來優(yōu)化合金性能提供理論依據(jù)。在研究方法上，本研究綜合運(yùn)用多種模擬方法和實(shí)驗(yàn)技術(shù)：第一性原理計(jì)算：基于密度泛函理論，利用VASP、CASTEP等計(jì)算軟件，從電子層次出發(fā)，計(jì)算合金體系的電子結(jié)構(gòu)、晶體結(jié)構(gòu)、形成能、結(jié)合能等物理量。通過這些計(jì)算結(jié)果，分析合金元素之間的相互作用、化學(xué)鍵性質(zhì)以及強(qiáng)化相形成的熱力學(xué)驅(qū)動(dòng)力，預(yù)測強(qiáng)化相的晶體結(jié)構(gòu)和穩(wěn)定性，為分子動(dòng)力學(xué)模擬提供初始結(jié)構(gòu)和參數(shù)。分子動(dòng)力學(xué)模擬：采用LAMMPS等模擬軟件，構(gòu)建合金體系的原子模型，通過設(shè)定合適的原子間相互作用勢（如EAM勢、L-J勢等），模擬合金在不同溫度、壓力條件下的凝固過程、原子擴(kuò)散行為以及強(qiáng)化相的形核和生長過程。在模擬過程中，跟蹤原子的運(yùn)動(dòng)軌跡，統(tǒng)計(jì)原子的分布和聚集情況，獲取強(qiáng)化相形成的動(dòng)力學(xué)信息，直觀地展現(xiàn)強(qiáng)化相形成的微觀過程。實(shí)驗(yàn)研究：開展合金制備實(shí)驗(yàn)，采用熔煉、鑄造、鍛造、軋制等工藝制備不同成分和工藝條件下的稀土鎂合金樣品；利用X射線衍射（XRD）分析合金的物相組成和晶體結(jié)構(gòu)；通過掃描電子顯微鏡（SEM）和透射電子顯微鏡（TEM）觀察合金的微觀組織和強(qiáng)化相的形貌、尺寸、分布及原子排列特征；運(yùn)用能譜分析（EDS）確定合金成分和強(qiáng)化相的化學(xué)組成；進(jìn)行力學(xué)性能測試（如拉伸、硬度測試）、物理性能測試（如熱膨脹系數(shù)測量）和化學(xué)性能測試（如耐腐蝕性能測試），獲取合金的宏觀性能數(shù)據(jù)，為模擬結(jié)果的驗(yàn)證和分析提供實(shí)驗(yàn)依據(jù)。二、高強(qiáng)度稀土鎂合金概述2.1鎂合金的特性與應(yīng)用2.1.1鎂合金的基本特性鎂合金是以鎂為基礎(chǔ)加入其他元素組成的合金，其密度約為1.8g/cm3，約為鋁的2/3、鐵的1/4，是實(shí)際應(yīng)用中最輕的金屬結(jié)構(gòu)材料，這一特性使其在對重量有嚴(yán)格要求的領(lǐng)域具有顯著優(yōu)勢，如航空航天領(lǐng)域，使用鎂合金可有效減輕飛行器重量，提高飛行性能。鎂合金具有較高的比強(qiáng)度和比剛度，比強(qiáng)度是指材料的強(qiáng)度與密度之比，比剛度是指材料的剛度與密度之比。鎂合金的比強(qiáng)度明顯高于鋁合金和鋼，比剛度與鋁合金和鋼相當(dāng)，這意味著在相同強(qiáng)度或剛度要求下，鎂合金可以制成更輕薄的結(jié)構(gòu)件，同時(shí)保持良好的力學(xué)性能，滿足工程應(yīng)用中對結(jié)構(gòu)輕量化和高強(qiáng)度的需求。鎂合金還具有良好的阻尼性能，在彈性范圍內(nèi)，受到?jīng)_擊載荷時(shí)，吸收的能量比鋁合金件大，其減振性是鋁的100倍，鈦合金的300-500倍。這一特性使鎂合金在需要減震降噪的場合，如汽車發(fā)動(dòng)機(jī)支架、電子設(shè)備外殼等部件中得到廣泛應(yīng)用，能夠有效減少振動(dòng)和噪音，提高設(shè)備的穩(wěn)定性和舒適性。在加工性能方面，鎂合金熔點(diǎn)比鋁合金熔點(diǎn)低，壓鑄成型性能好，壓鑄件壁厚最小可達(dá)0.5mm，能夠適應(yīng)復(fù)雜形狀零部件的制造，且鑄造加工尺寸精度高，可進(jìn)行高精度機(jī)械加工。此外，鎂合金還具有良好的導(dǎo)電導(dǎo)熱性能，其導(dǎo)熱性略低于鋁合金及銅合金，遠(yuǎn)高于鈦合金，在電子設(shè)備散熱領(lǐng)域具有應(yīng)用潛力；同時(shí)，鎂合金具有良好的電磁屏蔽性能，能夠完全吸收頻率超過100db的電磁干擾，可用于制作3C產(chǎn)品的外殼，提供優(yōu)越的抗電磁保護(hù)作用。2.1.2鎂合金在各領(lǐng)域的應(yīng)用現(xiàn)狀在航空航天領(lǐng)域，鎂合金被廣泛應(yīng)用于制造飛機(jī)、導(dǎo)彈、飛船、衛(wèi)星上的重要構(gòu)件，如設(shè)備支架、儀器儀表殼體、操縱系統(tǒng)支座、座艙骨架、發(fā)動(dòng)機(jī)附件機(jī)匣、直升機(jī)變速箱、發(fā)動(dòng)機(jī)架、機(jī)輪輪轂等零部件，其中鎂合金鑄件超過90%。例如，在天舟八號(hào)貨運(yùn)飛船任務(wù)中，重慶大學(xué)潘復(fù)生院士團(tuán)隊(duì)陳先華教授課題組研制的大規(guī)格復(fù)雜結(jié)構(gòu)功能一體化鎂合金適配構(gòu)件獲得成功應(yīng)用，有效助力任務(wù)實(shí)施，在確保設(shè)備功能穩(wěn)定和構(gòu)件承載可靠的同時(shí)，表現(xiàn)出顯著的減重效果和性能提升。然而，傳統(tǒng)鎂合金“絕對強(qiáng)度偏低”“變形能力偏弱”和“耐蝕性能不足”等問題，限制了其在航空航天領(lǐng)域更廣泛的應(yīng)用，尤其是在一些對材料性能要求極高的關(guān)鍵部件上的應(yīng)用受到制約。在汽車工業(yè)中，鎂合金可用于制造發(fā)動(dòng)機(jī)罩、進(jìn)氣歧管、變速器殼體、車輪等零部件，有助于實(shí)現(xiàn)汽車的輕量化，降低能耗和排放。隨著汽車行業(yè)對輕量化和環(huán)保化需求的不斷增加，尤其是新能源汽車的快速發(fā)展，鎂合金在汽車領(lǐng)域的應(yīng)用前景廣闊。如一些高端汽車品牌已經(jīng)開始在部分零部件上采用鎂合金材料，以提升車輛的性能和燃油經(jīng)濟(jì)性。但鎂合金的成本相對較高，以及在高溫和腐蝕環(huán)境下的性能穩(wěn)定性問題，在一定程度上阻礙了其在汽車工業(yè)中的大規(guī)模應(yīng)用。在電子領(lǐng)域，鎂合金憑借其良好的電磁屏蔽性能、散熱性能和輕薄特性，被用于制造手機(jī)、電腦、相機(jī)等3C產(chǎn)品的外殼和內(nèi)部結(jié)構(gòu)件。例如，蘋果公司的部分筆記本電腦外殼采用了鎂合金材質(zhì)，不僅減輕了產(chǎn)品重量，還提升了產(chǎn)品的質(zhì)感和外觀。然而，鎂合金的耐蝕性較差，在潮濕環(huán)境下容易發(fā)生腐蝕，影響產(chǎn)品的使用壽命，這是其在電子領(lǐng)域應(yīng)用中需要解決的關(guān)鍵問題之一。2.2稀土元素在鎂合金中的作用2.2.1稀土元素對鎂合金組織的影響稀土元素在鎂合金中具有凈化合金熔體的重要作用。在合金熔煉過程中，稀土元素因其活潑性僅次于堿金屬和堿土金屬，能與鎂合金中常見的雜質(zhì)如氫、氧化物夾雜、硫和鐵等發(fā)生化學(xué)反應(yīng)。例如，稀土元素可以與氫反應(yīng)生成穩(wěn)定的氫化物，從而降低合金中的氫含量，減少氣孔等缺陷的產(chǎn)生；與氧化物夾雜反應(yīng)，使其轉(zhuǎn)變?yōu)楦呷埸c(diǎn)、密度大的化合物，在合金凝固過程中這些化合物會(huì)浮出液面或沉淀在底部，達(dá)到除氣精煉、凈化熔體的效果，有效提高合金的純凈度，為后續(xù)形成良好的組織結(jié)構(gòu)奠定基礎(chǔ)。在變質(zhì)和細(xì)化組織方面，稀土元素主要通過影響合金的凝固過程來實(shí)現(xiàn)。一方面，稀土元素是鎂合金的表面活性元素，在合金液凝固過程中，會(huì)在固液界面前沿富集，提高成分過冷度。根據(jù)凝固理論，成分過冷度的增加有利于在過冷區(qū)形成新的形核帶，從而產(chǎn)生大量的晶核，這些晶核在生長過程中相互競爭，抑制了晶粒的長大，最終形成細(xì)等軸晶組織。例如，在Mg-Al系合金中添加稀土元素Ce，當(dāng)Ce的加入量為0.8％時(shí)，晶粒細(xì)化效果顯著，晶粒尺寸由未細(xì)化前的約300μm下降到約20-40μm。另一方面，稀土元素的富集還能起到阻礙α-Mg晶粒長大的作用。稀土原子與鎂原子之間存在一定的相互作用，這種作用會(huì)阻礙鎂原子的擴(kuò)散和晶粒的長大，進(jìn)一步促進(jìn)了晶粒的細(xì)化。同時(shí)，稀土元素還可以細(xì)化鎂合金中的第二相，如在Mg-Al系合金中，稀土元素可以使β-Mg17Al12相的數(shù)量減少、變細(xì)，使其分布更加均勻，改善合金的微觀組織均勻性。此外，稀土元素還能改善合金的微觀組織形態(tài)和分布。在一些鎂合金中，稀土元素與其他合金化元素反應(yīng)形成新的金屬間化合物相，這些相通常具有高熔點(diǎn)、高熱穩(wěn)定性等特點(diǎn)，它們呈細(xì)小化合物粒子彌散分布于晶界和晶內(nèi)。例如，在Mg-Zn-Y系合金中，會(huì)形成Mg-Zn-Y準(zhǔn)晶相，這些準(zhǔn)晶相彌散分布在基體中，不僅可以釘扎晶界，抑制晶界滑移，還能阻礙位錯(cuò)運(yùn)動(dòng)，從而強(qiáng)化合金基體，同時(shí)也改變了合金微觀組織的形態(tài)和分布，提高了合金組織的穩(wěn)定性和均勻性。2.2.2稀土元素對鎂合金性能的提升在力學(xué)性能提升方面，稀土元素主要通過多種強(qiáng)化機(jī)制來實(shí)現(xiàn)。首先是固溶強(qiáng)化，大部分稀土元素在鎂中具有較高的固溶度。當(dāng)稀土元素固溶于鎂基體時(shí)，由于稀土元素與鎂的原子半徑和彈性模量存在差異，會(huì)使鎂基體產(chǎn)生點(diǎn)陣畸變。這種點(diǎn)陣畸變產(chǎn)生的應(yīng)力場將阻礙位錯(cuò)運(yùn)動(dòng)，從而使鎂基體得到強(qiáng)化，提高合金的強(qiáng)度。例如，在Mg-Gd合金中，Gd原子固溶到鎂基體中，引起基體的點(diǎn)陣畸變，有效提高了合金的強(qiáng)度。其次是時(shí)效沉淀強(qiáng)化，稀土元素在鎂中所具有的較高固溶度隨溫度降低而降低。當(dāng)處于高溫下的單相固溶體快速冷卻時(shí)，形成不穩(wěn)定的過飽和固溶體，經(jīng)過長時(shí)間的時(shí)效，則形成細(xì)小而彌散的析出沉淀相。這些析出相與位錯(cuò)之間發(fā)生交互作用，阻礙位錯(cuò)的滑移，從而提高合金的強(qiáng)度。如在Mg-Nd合金時(shí)效過程中，會(huì)析出Mg12Nd相，該相彌散分布在基體中，顯著提高了合金的強(qiáng)度。此外，細(xì)晶強(qiáng)化也是重要的強(qiáng)化方式之一，如前文所述，稀土元素細(xì)化晶粒，根據(jù)Hall-Petch公式，合金的強(qiáng)度隨晶粒尺寸的細(xì)化而增加，且對于密排六方結(jié)構(gòu)的鎂合金，晶粒尺寸對強(qiáng)度的影響更為顯著，通過細(xì)化晶粒，鎂合金的強(qiáng)度和韌性都能得到有效提升。在抗氧化性能方面，稀土元素能提高鎂合金的抗氧化能力。稀土元素是鎂合金熔體的表面活性元素，在熔煉過程中，稀土元素易于在合金液表面聚集，與氧反應(yīng)形成MgO、RE2O3等多元復(fù)合致密氧化物層。這層致密的氧化物膜能夠有效阻止熔體和大氣的接觸，大大提高鎂合金熔體的起燃溫度，從而增強(qiáng)合金在熔煉和使用過程中的抗氧化性能。例如，在Mg-Al合金中添加稀土元素La后，合金表面形成了更加致密的氧化膜，其抗氧化性能得到明顯改善。對于蠕變性能，稀土元素的加入可以顯著提高鎂合金的抗蠕變性能。一方面，稀土原子質(zhì)量和半徑都比鎂原子大，能夠減慢原子的擴(kuò)散速率，在高溫下，原子擴(kuò)散速率的降低有助于抑制位錯(cuò)的攀移和滑移，從而提高合金的抗蠕變性能。另一方面，稀土元素與鎂或其他合金化元素在合金凝固過程中形成的穩(wěn)定金屬間化合物，如Mg3RE、Mg12Nd等，這些化合物呈細(xì)小粒子彌散分布于晶界和晶內(nèi)，在高溫下可以釘扎晶界，抑制晶界滑移，阻礙位錯(cuò)運(yùn)動(dòng)，從而有效提高合金的高溫強(qiáng)度和蠕變抗力。如在Mg-Y合金中，添加適量的Y元素形成的Mg24Y5相，能夠有效阻礙位錯(cuò)運(yùn)動(dòng)和晶界滑移，顯著提高合金的抗蠕變性能。2.3高強(qiáng)度稀土鎂合金的強(qiáng)化機(jī)制2.3.1固溶強(qiáng)化固溶強(qiáng)化是合金強(qiáng)化的基本方式之一，其原理基于溶質(zhì)原子與溶劑原子在尺寸和性質(zhì)上的差異。當(dāng)溶質(zhì)原子固溶入溶劑晶格中時(shí)，由于溶質(zhì)原子與溶劑原子的原子半徑不同，會(huì)使晶格發(fā)生畸變。這種晶格畸變會(huì)在晶體內(nèi)部產(chǎn)生應(yīng)力場，而位錯(cuò)是晶體中的一種線缺陷，其運(yùn)動(dòng)是材料發(fā)生塑性變形的主要方式。位錯(cuò)在運(yùn)動(dòng)過程中會(huì)受到晶格畸變產(chǎn)生的應(yīng)力場的阻礙，需要消耗更多的能量才能克服這種阻礙，從而使材料的強(qiáng)度和硬度提高。例如，在金屬晶體中，當(dāng)較小的溶質(zhì)原子溶入較大的溶劑原子晶格間隙時(shí)，會(huì)產(chǎn)生間隙固溶強(qiáng)化；當(dāng)溶質(zhì)原子取代溶劑原子占據(jù)晶格節(jié)點(diǎn)位置時(shí)，會(huì)產(chǎn)生置換固溶強(qiáng)化。無論是哪種固溶方式，晶格畸變產(chǎn)生的應(yīng)力場與位錯(cuò)之間的交互作用，都會(huì)阻礙位錯(cuò)的滑移，進(jìn)而強(qiáng)化合金。在稀土鎂合金中，大部分稀土元素在鎂中具有較高的固溶度。由于稀土元素與鎂的原子半徑和彈性模量存在差異，當(dāng)稀土元素固溶于鎂基體時(shí)，會(huì)使鎂基體產(chǎn)生點(diǎn)陣畸變。以Mg-Gd合金為例，Gd原子的半徑大于Mg原子，Gd原子固溶到鎂基體中后，會(huì)導(dǎo)致鎂基體的晶格常數(shù)發(fā)生變化，產(chǎn)生晶格畸變。這種晶格畸變產(chǎn)生的應(yīng)力場會(huì)阻礙位錯(cuò)的運(yùn)動(dòng)，使得位錯(cuò)在滑移過程中需要克服更大的阻力，從而提高了合金的強(qiáng)度。研究表明，隨著Gd元素含量的增加，Mg-Gd合金的強(qiáng)度和硬度逐漸提高，這充分體現(xiàn)了固溶強(qiáng)化在稀土鎂合金中的作用效果。此外，稀土元素的固溶強(qiáng)化作用還與溶質(zhì)原子的濃度和分布有關(guān)。當(dāng)溶質(zhì)原子濃度較高且分布均勻時(shí)，產(chǎn)生的晶格畸變更加明顯，對合金的強(qiáng)化效果也更為顯著。2.3.2析出強(qiáng)化析出強(qiáng)化是指合金在時(shí)效過程中，從過飽和固溶體中析出細(xì)小彌散的第二相粒子，這些粒子阻礙位錯(cuò)運(yùn)動(dòng)和滑移，從而提高合金強(qiáng)度的強(qiáng)化機(jī)制。其具體機(jī)制主要包括以下兩種：一是位錯(cuò)繞過機(jī)制，當(dāng)位錯(cuò)運(yùn)動(dòng)遇到彌散分布的第二相粒子時(shí)，如果位錯(cuò)不能直接切過粒子，位錯(cuò)會(huì)在粒子周圍發(fā)生彎曲，隨著位錯(cuò)的不斷運(yùn)動(dòng)，位錯(cuò)環(huán)會(huì)逐漸圍繞粒子形成，位錯(cuò)則繞過粒子繼續(xù)運(yùn)動(dòng)。在這個(gè)過程中，位錯(cuò)需要克服較大的阻力，從而增加了位錯(cuò)運(yùn)動(dòng)的難度，提高了合金的強(qiáng)度。二是位錯(cuò)切過機(jī)制，當(dāng)?shù)诙嗔Ｗ优c基體之間具有良好的共格或半共格關(guān)系，且粒子尺寸較小時(shí)，位錯(cuò)可以直接切過粒子。在切過過程中，位錯(cuò)會(huì)破壞粒子與基體之間的共格關(guān)系，產(chǎn)生新的界面，同時(shí)還會(huì)引起粒子內(nèi)部的晶格畸變，這些都會(huì)消耗額外的能量，阻礙位錯(cuò)的運(yùn)動(dòng)，進(jìn)而強(qiáng)化合金。在稀土鎂合金中，析出強(qiáng)化發(fā)揮著重要作用。例如，在Mg-Nd合金時(shí)效過程中，會(huì)從過飽和固溶體中析出Mg12Nd相。這些Mg12Nd相呈細(xì)小彌散的顆粒狀分布在鎂基體中，與基體之間存在一定的共格或半共格關(guān)系。當(dāng)位錯(cuò)運(yùn)動(dòng)到Mg12Nd相粒子附近時(shí)，會(huì)受到粒子的阻礙，位錯(cuò)需要通過繞過或切過粒子的方式才能繼續(xù)運(yùn)動(dòng)，這就大大增加了位錯(cuò)運(yùn)動(dòng)的阻力，提高了合金的強(qiáng)度。研究發(fā)現(xiàn)，通過控制時(shí)效工藝參數(shù)，如時(shí)效溫度和時(shí)效時(shí)間，可以調(diào)整析出相的尺寸、數(shù)量和分布，從而優(yōu)化合金的力學(xué)性能。在較低的時(shí)效溫度下，析出相的形核速率較快，但生長速率較慢，會(huì)形成大量細(xì)小彌散的析出相，此時(shí)合金的強(qiáng)度較高；而在較高的時(shí)效溫度下，析出相的生長速率加快，尺寸逐漸增大，數(shù)量相對減少，合金的強(qiáng)度會(huì)有所下降。因此，合理控制時(shí)效工藝，能夠充分發(fā)揮析出強(qiáng)化的作用，有效提高稀土鎂合金的強(qiáng)度和硬度。2.3.3彌散強(qiáng)化彌散強(qiáng)化是通過在合金基體中引入細(xì)小彌散分布的第二相質(zhì)點(diǎn)來阻礙位錯(cuò)運(yùn)動(dòng)，從而提高合金強(qiáng)度和硬度的強(qiáng)化方式。這些彌散質(zhì)點(diǎn)通常具有高硬度、高熔點(diǎn)和熱穩(wěn)定性，它們與基體之間沒有或僅有微弱的共格關(guān)系。位錯(cuò)在遇到彌散質(zhì)點(diǎn)時(shí)，由于無法切過質(zhì)點(diǎn)，只能繞過質(zhì)點(diǎn)繼續(xù)運(yùn)動(dòng)，這就增加了位錯(cuò)運(yùn)動(dòng)的路徑長度，使位錯(cuò)運(yùn)動(dòng)變得更加困難，從而提高了合金的強(qiáng)度。與析出強(qiáng)化相比，彌散強(qiáng)化中的彌散質(zhì)點(diǎn)一般是通過外加方式引入合金中的，而析出強(qiáng)化中的析出相是在合金時(shí)效過程中從過飽和固溶體中自發(fā)析出的。此外，彌散質(zhì)點(diǎn)的尺寸和分布相對更加穩(wěn)定，在高溫下不易長大和聚集，因此彌散強(qiáng)化在提高合金高溫強(qiáng)度方面具有獨(dú)特的優(yōu)勢。在稀土鎂合金中，稀土元素與鎂或其他合金化元素在合金凝固過程中形成的穩(wěn)定金屬間化合物，如Mg3RE、Mg12Nd等，這些含稀土的金屬間化合物一般具有高熔點(diǎn)、高熱穩(wěn)定性等特點(diǎn)，它們呈細(xì)小化合物粒子彌散分布于晶界和晶內(nèi)，起到了彌散強(qiáng)化的作用。例如，在Mg-Zn-Y系合金中形成的Mg-Zn-Y準(zhǔn)晶相，這些準(zhǔn)晶相彌散分布在基體中，在高溫下可以有效釘扎晶界，抑制晶界滑移，同時(shí)阻礙位錯(cuò)運(yùn)動(dòng)。當(dāng)合金在高溫下承受外力時(shí)，位錯(cuò)運(yùn)動(dòng)到準(zhǔn)晶相粒子附近會(huì)受到阻礙，需要消耗更多的能量才能繞過粒子，從而提高了合金的高溫強(qiáng)度。研究表明，彌散強(qiáng)化可以顯著提高稀土鎂合金在高溫下的抗變形能力，使其在高溫環(huán)境下仍能保持較好的力學(xué)性能，這對于稀土鎂合金在航空航天、汽車等高溫應(yīng)用領(lǐng)域具有重要意義。2.3.4細(xì)晶強(qiáng)化細(xì)晶強(qiáng)化的原理基于晶界對滑移傳遞的阻礙作用。晶體中的位錯(cuò)在滑移過程中，當(dāng)遇到晶界時(shí)，由于晶界處原子排列不規(guī)則，位錯(cuò)難以直接穿過晶界，需要在晶界處堆積。隨著位錯(cuò)的不斷堆積，會(huì)在晶界附近產(chǎn)生應(yīng)力集中，當(dāng)應(yīng)力集中達(dá)到一定程度時(shí)，會(huì)促使相鄰晶粒內(nèi)的位錯(cuò)源開動(dòng)，從而使滑移傳遞到相鄰晶粒中。這個(gè)過程需要消耗更多的能量，增加了材料的變形難度。此外，細(xì)晶粒組織中晶界面積較大，晶界可以將材料分割成許多細(xì)小的晶粒區(qū)域，使得位錯(cuò)的滑移距離減小。根據(jù)Hall-Petch公式，合金的屈服強(qiáng)度與晶粒尺寸的平方根成反比，即晶粒尺寸越小，合金的強(qiáng)度越高。同時(shí)，細(xì)晶粒組織還能使合金在變形過程中各晶粒的變形更加均勻，減少應(yīng)力集中，從而提高合金的韌性。在稀土鎂合金中，細(xì)晶強(qiáng)化效果顯著。如前文所述，稀土元素在合金凝固過程中，通過在固液界面前沿富集，提高成分過冷度，形成新的形核帶，產(chǎn)生大量晶核，抑制晶粒長大，從而細(xì)化晶粒。以在Mg-Al系合金中添加稀土元素Ce為例，當(dāng)Ce的加入量為0.8％時(shí)，晶粒尺寸由未細(xì)化前的約300μm下降到約20-40μm。隨著晶粒的細(xì)化，合金的強(qiáng)度和韌性都得到了有效提升。研究發(fā)現(xiàn)，細(xì)晶強(qiáng)化對于密排六方結(jié)構(gòu)的鎂合金尤為重要，因?yàn)殒V合金的晶體結(jié)構(gòu)特點(diǎn)使其在變形過程中滑移系較少，而細(xì)晶粒可以增加晶界數(shù)量，為位錯(cuò)運(yùn)動(dòng)提供更多的阻礙和協(xié)調(diào)機(jī)制，從而顯著提高合金的強(qiáng)度和韌性。此外，細(xì)晶強(qiáng)化還能與其他強(qiáng)化機(jī)制（如固溶強(qiáng)化、析出強(qiáng)化等）相互配合，進(jìn)一步提升稀土鎂合金的綜合力學(xué)性能。三、模擬研究方法與模型建立3.1模擬方法選擇3.1.1分子動(dòng)力學(xué)模擬原理與優(yōu)勢分子動(dòng)力學(xué)模擬是一種基于經(jīng)典力學(xué)原理的計(jì)算機(jī)模擬方法，主要用于模擬和研究物質(zhì)在原子或分子尺度的動(dòng)態(tài)行為。其基本原理緊密依賴于牛頓第二定律和分子間相互作用勢。在模擬開始時(shí)，首先需要對系統(tǒng)中的粒子（原子或分子）進(jìn)行初始化設(shè)置，賦予它們初始位置、速度和加速度等參數(shù)。然后，依據(jù)粒子間的相互作用勢，如范德華力、庫侖力等，以及牛頓第二定律F=ma（其中F表示粒子所受的力，m為粒子質(zhì)量，a是粒子的加速度），計(jì)算每個(gè)粒子所受的力和加速度。通過這些計(jì)算得到的力和加速度，進(jìn)而更新粒子的速度和位置，其更新公式通常基于數(shù)值積分算法，如Verlet算法、leap-frog算法等。這個(gè)過程在時(shí)間上逐步推進(jìn)，不斷重復(fù)計(jì)算和更新，從而模擬出物質(zhì)在特定條件下隨時(shí)間的動(dòng)態(tài)演化過程。在研究高強(qiáng)度稀土鎂合金中，分子動(dòng)力學(xué)模擬具有顯著的優(yōu)勢。從原子尺度動(dòng)態(tài)信息獲取角度來看，它能夠直觀地展示合金在凝固過程中原子的擴(kuò)散、聚集和排列行為。通過模擬，可以清晰地觀察到稀土元素與鎂及其他合金化元素原子之間的相互作用過程，以及強(qiáng)化相在形核和生長過程中原子的動(dòng)態(tài)變化。例如，在模擬Mg-Gd-Y合金的凝固過程時(shí)，可以追蹤Gd、Y原子在鎂基體中的擴(kuò)散路徑，以及它們?nèi)绾蜗嗷ゾ奂纬珊⊥恋膹?qiáng)化相粒子，為深入理解強(qiáng)化相的形成機(jī)制提供了直接的原子尺度信息。而且，分子動(dòng)力學(xué)模擬可以在不同的溫度、壓力等條件下進(jìn)行，通過改變模擬參數(shù)，能夠研究這些外部條件對合金微觀結(jié)構(gòu)和性能的影響。比如，通過模擬不同冷卻速度下合金的凝固過程，可以分析冷卻速度對強(qiáng)化相尺寸、分布和原子排列的影響規(guī)律，這對于優(yōu)化合金制備工藝具有重要的指導(dǎo)意義。此外，相比于一些實(shí)驗(yàn)研究方法，分子動(dòng)力學(xué)模擬的計(jì)算速度相對較快，能夠在較短時(shí)間內(nèi)獲得大量的模擬數(shù)據(jù)，為理論研究提供豐富的信息資源。它可以在實(shí)驗(yàn)之前對合金的性能和微觀結(jié)構(gòu)進(jìn)行預(yù)測，幫助科研人員設(shè)計(jì)實(shí)驗(yàn)方案，減少實(shí)驗(yàn)次數(shù)和成本。同時(shí)，對于一些難以通過實(shí)驗(yàn)直接觀測的微觀過程，如原子在晶界處的擴(kuò)散、位錯(cuò)與強(qiáng)化相的交互作用等，分子動(dòng)力學(xué)模擬能夠提供有效的研究手段。3.1.2第一性原理計(jì)算方法簡介第一性原理計(jì)算是基于量子力學(xué)原理，從基本的物理和化學(xué)定律出發(fā)，對材料的原子結(jié)構(gòu)和能量進(jìn)行計(jì)算的方法，它不需要任何實(shí)驗(yàn)參數(shù)，僅依靠一些基本的物理常量，如電子質(zhì)量、光速、質(zhì)子中子質(zhì)量等，就可以計(jì)算出體系的能量和電子結(jié)構(gòu)等物理性質(zhì)。其核心思想是將多個(gè)原子構(gòu)成的體系看成是由多個(gè)電子和原子核組成的系統(tǒng)，并根據(jù)量子力學(xué)的基本原理對問題進(jìn)行最大限度的“非經(jīng)驗(yàn)性”處理。在實(shí)際計(jì)算中，第一性原理計(jì)算通常基于密度泛函理論（DFT）。DFT將電子密度作為描述多電子系統(tǒng)的關(guān)鍵參數(shù)，通過最小化材料的總能量來計(jì)算材料的電子結(jié)構(gòu)和能帶結(jié)構(gòu)。其基本方程是庫侖能和交換-相關(guān)能的和等于總能量。其中，庫侖能是電子間的相互作用能，交換-相關(guān)能是描述電子間交換作用和相關(guān)作用的能量。庫侖能可以通過庫侖定律計(jì)算得到，而交換-相關(guān)能通常需要使用一些近似方法來處理。在DFT中，電子的波函數(shù)被表示為電子密度的泛函形式，即電子的密度可以唯一地確定電子的基態(tài)波函數(shù)。由此得到的Kohn-Sham方程可以通過數(shù)值方法求解得到電子密度和能帶結(jié)構(gòu)。在研究高強(qiáng)度稀土鎂合金中，第一性原理計(jì)算在探究原子排列方面發(fā)揮著重要作用。它可以精確計(jì)算合金體系中原子之間的相互作用能，從而確定不同原子排列方式下合金的穩(wěn)定性。例如，對于含稀土強(qiáng)化相的原子排列，通過第一性原理計(jì)算可以分析不同晶體結(jié)構(gòu)、原子占位情況下強(qiáng)化相的形成能和結(jié)合能，判斷哪種原子排列方式最為穩(wěn)定，以及不同排列方式對強(qiáng)化相性能的影響。此外，第一性原理計(jì)算還能計(jì)算合金的電子結(jié)構(gòu)，如電子云分布、能帶結(jié)構(gòu)等，從電子層次揭示合金元素之間的化學(xué)鍵性質(zhì)和相互作用機(jī)制，為理解合金的物理和化學(xué)性質(zhì)提供深層次的理論依據(jù)。比如，通過分析合金的電子結(jié)構(gòu)，可以了解稀土元素與鎂及其他合金化元素之間的電子轉(zhuǎn)移和共享情況，進(jìn)而解釋合金的強(qiáng)化機(jī)制和物理性能。3.2模型構(gòu)建3.2.1確定模擬體系與原子模型本研究選擇典型的Mg-RE系合金作為模擬體系，具體以Mg-Gd-Y合金為例，其合金成分為Mg-10at%Gd-5at%Y。選擇該合金成分主要基于以下考慮：一方面，Gd和Y是稀土元素中在鎂合金研究和應(yīng)用中較為常用且具有顯著強(qiáng)化效果的元素。Gd在鎂中具有較高的固溶度，能夠通過固溶強(qiáng)化和析出強(qiáng)化有效提高鎂合金的強(qiáng)度和硬度；Y元素則可以細(xì)化鎂合金的晶粒，同時(shí)與其他元素形成穩(wěn)定的金屬間化合物，增強(qiáng)合金的高溫性能。另一方面，該成分體系在前期的實(shí)驗(yàn)和理論研究中已有一定的基礎(chǔ)，相關(guān)的實(shí)驗(yàn)數(shù)據(jù)和研究成果可為本次模擬研究提供參考和驗(yàn)證依據(jù)。在構(gòu)建原子模型時(shí)，采用周期性邊界條件構(gòu)建模擬盒子。首先，根據(jù)Mg-Gd-Y合金的晶體結(jié)構(gòu)信息，確定模擬盒子的初始原子坐標(biāo)。Mg為密排六方結(jié)構(gòu)（hcp），其晶格常數(shù)a=0.3209nm，c=0.5211nm。在模擬盒子中，按照hcp結(jié)構(gòu)的原子排列方式，放置Mg原子。然后，根據(jù)設(shè)定的合金成分，隨機(jī)替換部分Mg原子為Gd和Y原子。為了保證原子分布的隨機(jī)性和均勻性，使用隨機(jī)數(shù)生成算法來確定被替換原子的位置。例如，通過在模擬盒子內(nèi)隨機(jī)生成三維坐標(biāo)，將坐標(biāo)對應(yīng)的Mg原子替換為Gd或Y原子，直至達(dá)到設(shè)定的合金成分比例。最終構(gòu)建的模擬盒子尺寸為5nm×5nm×5nm，包含約5000個(gè)原子，該尺寸既能保證模擬體系具有足夠的代表性，又能在計(jì)算資源可承受的范圍內(nèi)進(jìn)行高效模擬。3.2.2設(shè)定模擬參數(shù)與條件在分子動(dòng)力學(xué)模擬過程中，合理設(shè)定模擬參數(shù)和條件對于獲得準(zhǔn)確可靠的模擬結(jié)果至關(guān)重要。溫度參數(shù)的設(shè)定基于合金的實(shí)際應(yīng)用場景和模擬目的。考慮到Mg-Gd-Y合金在高溫環(huán)境下的應(yīng)用需求，將模擬溫度設(shè)定為600K、700K和800K三個(gè)不同的溫度點(diǎn)。600K接近合金的實(shí)際服役溫度下限，700K為合金在一些常規(guī)高溫應(yīng)用中的典型溫度，800K則用于研究合金在接近熔點(diǎn)溫度下的微觀結(jié)構(gòu)演變。這些溫度點(diǎn)的選擇能夠全面覆蓋合金在不同工況下的溫度范圍，有助于深入研究溫度對強(qiáng)化相形成及原子排列的影響。溫度的控制采用Nose-Hoover溫控器，它通過調(diào)節(jié)系統(tǒng)的動(dòng)能來維持設(shè)定的溫度，能夠有效保證模擬體系在恒定溫度下進(jìn)行演化。壓力參數(shù)設(shè)定為1atm（101325Pa），近似于常壓環(huán)境。在實(shí)際的合金制備和應(yīng)用過程中，大部分情況處于常壓或接近常壓的條件下。因此，設(shè)定常壓條件能夠更真實(shí)地模擬合金在實(shí)際生產(chǎn)和使用過程中的狀態(tài)。通過在模擬中施加恒定的壓力，能夠保證模擬體系的體積和密度在合理范圍內(nèi)，避免因壓力變化對強(qiáng)化相形成及原子排列產(chǎn)生不必要的干擾。時(shí)間步長是分子動(dòng)力學(xué)模擬中的一個(gè)關(guān)鍵參數(shù)，它決定了模擬過程中每一步的時(shí)間間隔。經(jīng)過多次預(yù)模擬和對比分析，最終將時(shí)間步長設(shè)定為1fs（1×10?1?s）。選擇該時(shí)間步長主要基于以下考慮：一方面，時(shí)間步長不能過大，否則會(huì)導(dǎo)致模擬過程中原子的運(yùn)動(dòng)軌跡出現(xiàn)較大偏差，無法準(zhǔn)確描述原子的真實(shí)運(yùn)動(dòng)狀態(tài)，甚至可能導(dǎo)致模擬結(jié)果的不穩(wěn)定和發(fā)散；另一方面，時(shí)間步長也不能過小，過小的時(shí)間步長會(huì)顯著增加計(jì)算量和計(jì)算時(shí)間，降低模擬效率。1fs的時(shí)間步長在保證模擬精度的前提下，能夠在合理的計(jì)算時(shí)間內(nèi)完成模擬任務(wù)。在模擬過程中，總模擬步數(shù)設(shè)定為100000步，總模擬時(shí)間為100ps，這樣的模擬時(shí)長能夠充分展現(xiàn)合金在不同溫度和壓力條件下強(qiáng)化相形成及原子排列的動(dòng)態(tài)演變過程。3.3模擬流程與驗(yàn)證3.3.1詳細(xì)模擬步驟在分子動(dòng)力學(xué)模擬中，初始化過程是模擬的起點(diǎn)。首先，利用MaterialsStudio軟件搭建Mg-Gd-Y合金的初始原子模型。在構(gòu)建模型時(shí)，依據(jù)Mg-Gd-Y合金的晶體結(jié)構(gòu)信息，確定模擬盒子的尺寸和原子坐標(biāo)。Mg為密排六方結(jié)構(gòu)（hcp），其晶格常數(shù)a=0.3209nm，c=0.5211nm，按照hcp結(jié)構(gòu)的原子排列方式在模擬盒子中放置Mg原子。然后，根據(jù)設(shè)定的合金成分Mg-10at%Gd-5at%Y，使用隨機(jī)數(shù)生成算法，隨機(jī)替換部分Mg原子為Gd和Y原子，以保證原子分布的隨機(jī)性和均勻性。構(gòu)建好原子模型后，對模型進(jìn)行能量最小化處理，采用共軛梯度法等優(yōu)化算法，調(diào)整原子的位置，使體系的能量達(dá)到最低，消除初始結(jié)構(gòu)中的不合理應(yīng)力，為后續(xù)模擬提供穩(wěn)定的初始狀態(tài)。弛豫階段對于模擬體系達(dá)到平衡狀態(tài)至關(guān)重要。在能量最小化后的模型基礎(chǔ)上，將體系置于NPT系綜下進(jìn)行弛豫。NPT系綜即等溫等壓系綜，在該系綜下，體系的溫度、壓力和粒子數(shù)保持不變。溫度的控制采用Nose-Hoover溫控器，它通過調(diào)節(jié)系統(tǒng)的動(dòng)能來維持設(shè)定的溫度；壓力的控制采用Parrinello-Rahman壓控器，通過調(diào)整模擬盒子的體積來維持設(shè)定的壓力。在NPT系綜下，模擬體系經(jīng)歷一定的時(shí)間步長進(jìn)行弛豫，使體系的原子分布和結(jié)構(gòu)達(dá)到穩(wěn)定狀態(tài)。經(jīng)過多次測試和分析，確定弛豫時(shí)間為50ps，時(shí)間步長為1fs。在弛豫過程中，原子會(huì)在模擬盒子內(nèi)進(jìn)行擴(kuò)散、遷移等運(yùn)動(dòng)，體系的能量、體積和密度等參數(shù)會(huì)逐漸趨于穩(wěn)定。通過監(jiān)測體系的能量、壓力和體積等參數(shù)隨時(shí)間的變化曲線，當(dāng)這些參數(shù)在一定范圍內(nèi)波動(dòng)且波動(dòng)幅度較小時(shí)，表明體系已達(dá)到平衡狀態(tài)。數(shù)據(jù)采集與分析是模擬的關(guān)鍵環(huán)節(jié)。在體系達(dá)到平衡狀態(tài)后，切換到NVT系綜進(jìn)行數(shù)據(jù)采集。NVT系綜即正則系綜，在該系綜下，體系的粒子數(shù)、體積和溫度保持不變。在NVT系綜下，再運(yùn)行模擬一定的時(shí)間步長，如50ps，每隔一定的時(shí)間間隔（如100步）采集一次原子的位置、速度、能量等數(shù)據(jù)。利用這些采集到的數(shù)據(jù)，進(jìn)行深入的分析。例如，通過計(jì)算原子的均方位移（MSD）來研究原子的擴(kuò)散行為，均方位移的計(jì)算公式為MSD(t)=\frac{1}{N}\sum_{i=1}^{N}\left\langle\left(r_{i}(t)-r_{i}(0)\right)^{2}\right\rangle，其中N為原子總數(shù)，r_{i}(t)和r_{i}(0)分別為第i個(gè)原子在t時(shí)刻和初始時(shí)刻的位置，通過分析MSD隨時(shí)間的變化曲線，可以得到原子的擴(kuò)散系數(shù)，了解原子在不同溫度和壓力條件下的擴(kuò)散速率和擴(kuò)散機(jī)制。此外，還可以計(jì)算原子的徑向分布函數(shù)（RDF），徑向分布函數(shù)的計(jì)算公式為g(r)=\frac{V}{4\pir^{2}\DeltarN}\sum_{i=1}^{N}\sum_{j\neqi}^{N}\delta(r_{ij}-r)，其中V為模擬盒子的體積，\Deltar為徑向間距的微小增量，r_{ij}為原子i和j之間的距離，通過分析RDF曲線，可以了解原子之間的相對位置分布和原子團(tuán)簇的形成情況，揭示合金中原子的排列特征和微觀結(jié)構(gòu)信息。在第一性原理計(jì)算模擬中，同樣首先進(jìn)行初始結(jié)構(gòu)構(gòu)建。運(yùn)用VASP軟件構(gòu)建Mg-Gd-Y合金的超晶胞模型，超晶胞的大小根據(jù)研究需求和計(jì)算資源確定，例如選擇2\times2\times2的超晶胞，以充分考慮原子間的相互作用。在構(gòu)建模型時(shí)，明確各原子的坐標(biāo)和晶格參數(shù)，確保模型的準(zhǔn)確性。然后，進(jìn)行參數(shù)設(shè)置。選擇廣義梯度近似（GGA）下的PBE泛函來描述電子的交換-相關(guān)能，這種泛函在處理金屬體系時(shí)具有較好的準(zhǔn)確性。平面波截?cái)嗄茉O(shè)置為500eV，以保證計(jì)算的精度。k點(diǎn)網(wǎng)格采用Monkhorst-Pack方法生成，根據(jù)超晶胞的大小和對稱性確定合適的k點(diǎn)密度，如對于2\times2\times2的超晶胞，選擇4\times4\times4的k點(diǎn)網(wǎng)格，以充分考慮電子的動(dòng)量分布。在計(jì)算過程中，進(jìn)行自洽場（SCF）迭代計(jì)算，使體系的總能量收斂。收斂標(biāo)準(zhǔn)設(shè)置為總能量的變化小于1\times10^{-5}eV/atom，原子力小于0.01eV/?，以確保計(jì)算結(jié)果的可靠性。通過這些參數(shù)設(shè)置和計(jì)算過程，能夠準(zhǔn)確計(jì)算合金體系的電子結(jié)構(gòu)、形成能等物理量。3.3.2模擬結(jié)果的驗(yàn)證方法為了驗(yàn)證分子動(dòng)力學(xué)模擬結(jié)果的準(zhǔn)確性，將模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對比。在合金成分對強(qiáng)化相形成影響的模擬中，模擬不同合金成分下強(qiáng)化相的種類和數(shù)量，然后通過實(shí)驗(yàn)制備相同成分的Mg-Gd-Y合金樣品。利用X射線衍射（XRD）分析合金的物相組成，通過XRD圖譜中衍射峰的位置和強(qiáng)度，確定合金中存在的相，與模擬預(yù)測的強(qiáng)化相種類進(jìn)行對比。例如，若模擬預(yù)測在某一合金成分下會(huì)形成Mg5Gd相，通過XRD分析實(shí)驗(yàn)樣品，若在相應(yīng)的衍射角度出現(xiàn)Mg5Gd相的特征衍射峰，則說明模擬結(jié)果與實(shí)驗(yàn)結(jié)果相符。同時(shí)，采用掃描電子顯微鏡（SEM）和透射電子顯微鏡（TEM）觀察合金的微觀組織，測量強(qiáng)化相的尺寸和分布，與模擬得到的強(qiáng)化相尺寸和分布數(shù)據(jù)進(jìn)行對比。如在模擬中得到某一強(qiáng)化相的平均尺寸為d_1，通過TEM觀察實(shí)驗(yàn)樣品中該強(qiáng)化相的尺寸，計(jì)算其平均尺寸為d_2，若d_1和d_2在合理的誤差范圍內(nèi)，則驗(yàn)證了模擬結(jié)果在強(qiáng)化相尺寸方面的準(zhǔn)確性。在制備工藝對強(qiáng)化相形成及原子排列影響的模擬中，將模擬得到的不同制備工藝下強(qiáng)化相的形成過程和原子排列特征與實(shí)驗(yàn)結(jié)果進(jìn)行對比。例如，模擬鑄造工藝中不同冷卻速度下合金的凝固過程和強(qiáng)化相的形成，通過實(shí)驗(yàn)設(shè)置不同的冷卻速度進(jìn)行合金鑄造，然后利用SEM和TEM觀察實(shí)驗(yàn)樣品中強(qiáng)化相的形貌、尺寸和原子排列方式。若模擬預(yù)測在較快冷卻速度下強(qiáng)化相尺寸較小且分布更均勻，而在實(shí)驗(yàn)中觀察到相同冷卻速度條件下強(qiáng)化相的實(shí)際尺寸和分布情況與模擬結(jié)果一致，則表明模擬結(jié)果能夠準(zhǔn)確反映制備工藝對強(qiáng)化相形成及原子排列的影響。對于第一性原理計(jì)算模擬結(jié)果的驗(yàn)證，與其他理論計(jì)算結(jié)果進(jìn)行對比。在計(jì)算合金體系的形成能、結(jié)合能等物理量時(shí)，參考其他已有的理論計(jì)算研究成果。例如，對于Mg-Gd-Y合金體系，若其他研究采用不同的計(jì)算方法或模型得到了該合金體系中某一強(qiáng)化相的形成能為E_1，本研究通過第一性原理計(jì)算得到的形成能為E_2，對E_1和E_2進(jìn)行對比分析。若兩者數(shù)值相近，且在誤差允許的范圍內(nèi)，說明本研究的計(jì)算結(jié)果具有可靠性。同時(shí)，還可以將計(jì)算得到的電子結(jié)構(gòu)、能帶結(jié)構(gòu)等結(jié)果與相關(guān)的理論模型和已有研究進(jìn)行對比。例如，根據(jù)晶體結(jié)構(gòu)理論和電子理論，分析計(jì)算得到的能帶結(jié)構(gòu)中導(dǎo)帶和價(jià)帶的位置、寬度以及能隙大小等特征，與理論預(yù)期和其他研究結(jié)果進(jìn)行比較。若計(jì)算結(jié)果與理論模型和已有研究相符，則驗(yàn)證了第一性原理計(jì)算模擬結(jié)果在電子結(jié)構(gòu)和能帶結(jié)構(gòu)方面的準(zhǔn)確性。四、強(qiáng)化相形成的模擬結(jié)果與分析4.1強(qiáng)化相形成過程模擬4.1.1凝固過程中強(qiáng)化相的形核與長大通過分子動(dòng)力學(xué)模擬，對Mg-Gd-Y合金在凝固過程中強(qiáng)化相的形核與長大過程進(jìn)行了詳細(xì)研究。在模擬過程中，設(shè)定合金初始處于高溫液態(tài)，溫度為800K，然后以1K/ps的冷卻速度緩慢冷卻至室溫300K。模擬結(jié)果清晰地展示了強(qiáng)化相在凝固過程中的形核位置和長大過程。在凝固初期，當(dāng)溫度降低到一定程度時(shí)，在鎂基體中某些原子聚集區(qū)域開始出現(xiàn)形核現(xiàn)象。這些形核位置并非完全隨機(jī)，而是傾向于在溶質(zhì)原子（Gd和Y）相對富集的區(qū)域形成。這是因?yàn)槿苜|(zhì)原子的富集能夠降低形核的自由能壘，使得形核更容易發(fā)生。通過對模擬過程中原子分布的分析發(fā)現(xiàn)，Gd和Y原子在液態(tài)時(shí)會(huì)發(fā)生一定程度的團(tuán)簇現(xiàn)象，這些團(tuán)簇區(qū)域成為強(qiáng)化相形核的優(yōu)先位置。例如，在某一時(shí)刻的模擬快照中，可以觀察到在一個(gè)局部區(qū)域內(nèi)，Gd和Y原子形成了一個(gè)相對密集的原子團(tuán)，隨后以這個(gè)原子團(tuán)為核心，周圍的Mg原子逐漸聚集，形成了強(qiáng)化相的晶核。隨著凝固過程的繼續(xù)進(jìn)行，晶核開始長大。強(qiáng)化相的長大主要通過原子的擴(kuò)散和聚集來實(shí)現(xiàn)。在晶核周圍，液態(tài)中的原子不斷向晶核表面擴(kuò)散，并逐漸排列到晶核的晶格結(jié)構(gòu)中，使得晶核的尺寸不斷增大。從模擬結(jié)果的時(shí)間序列分析中可以看出，晶核的長大速率在初期較快，隨著時(shí)間的推移，長大速率逐漸減緩。這是因?yàn)樵诰Ш碎L大初期，晶核與周圍液態(tài)之間的濃度梯度較大，原子擴(kuò)散驅(qū)動(dòng)力較強(qiáng)，使得原子能夠快速地向晶核表面遷移。然而，隨著晶核尺寸的增大，晶核周圍的溶質(zhì)原子逐漸被消耗，濃度梯度減小，原子擴(kuò)散驅(qū)動(dòng)力減弱，同時(shí)晶核表面的原子排列逐漸趨于穩(wěn)定，也限制了原子的進(jìn)一步附著，導(dǎo)致長大速率逐漸降低。在強(qiáng)化相長大過程中，還觀察到了晶核的合并現(xiàn)象。當(dāng)兩個(gè)或多個(gè)晶核在長大過程中相互靠近時(shí)，它們會(huì)逐漸合并成一個(gè)更大的晶粒。這種晶核合并現(xiàn)象會(huì)影響強(qiáng)化相的最終尺寸和分布。通過統(tǒng)計(jì)模擬過程中不同時(shí)刻強(qiáng)化相的尺寸分布，發(fā)現(xiàn)隨著凝固時(shí)間的增加，強(qiáng)化相的平均尺寸逐漸增大，尺寸分布的范圍也逐漸變寬。在凝固后期，由于晶核合并現(xiàn)象的增多，大尺寸強(qiáng)化相的數(shù)量相對增加，而小尺寸強(qiáng)化相的數(shù)量則相應(yīng)減少。影響強(qiáng)化相形核與長大的因素眾多，其中溫度是一個(gè)關(guān)鍵因素。通過改變模擬的冷卻速度進(jìn)行對比模擬發(fā)現(xiàn)，冷卻速度越快，強(qiáng)化相的形核數(shù)量越多，但晶核的長大時(shí)間相對較短，導(dǎo)致最終形成的強(qiáng)化相尺寸較小且分布更均勻。這是因?yàn)榭焖倮鋮s使得合金液的過冷度迅速增大，增加了形核的驅(qū)動(dòng)力，從而促進(jìn)了更多晶核的形成。然而，由于冷卻速度快，原子擴(kuò)散的時(shí)間較短，晶核沒有足夠的時(shí)間長大，所以強(qiáng)化相尺寸較小。相反，冷卻速度較慢時(shí)，形核數(shù)量相對較少，但晶核有更多的時(shí)間長大，最終形成的強(qiáng)化相尺寸較大且分布不均勻。合金成分對強(qiáng)化相形核與長大也有顯著影響。在不同的Gd和Y含量下進(jìn)行模擬，結(jié)果表明，隨著Gd和Y含量的增加，強(qiáng)化相的形核數(shù)量增多，長大速度加快。這是因?yàn)檩^高的溶質(zhì)含量增加了形核的幾率，同時(shí)也為晶核的長大提供了更多的原子供應(yīng)。例如，在Mg-12at%Gd-6at%Y合金成分下，強(qiáng)化相的形核數(shù)量明顯多于Mg-10at%Gd-5at%Y合金成分，且在相同的凝固時(shí)間內(nèi)，前者形成的強(qiáng)化相尺寸更大。此外，溶質(zhì)原子之間的相互作用也會(huì)影響強(qiáng)化相的形核與長大。Gd和Y原子之間可能存在一定的相互吸引作用，使得它們更容易聚集在一起，形成穩(wěn)定的原子團(tuán)簇，從而促進(jìn)強(qiáng)化相的形核。4.1.2不同成分體系下強(qiáng)化相形成差異為了深入探究不同成分體系下強(qiáng)化相形成的差異，對比了三種不同稀土元素含量的Mg-Gd-Y合金體系，分別為Mg-8at%Gd-3at%Y、Mg-10at%Gd-5at%Y和Mg-12at%Gd-6at%Y。在凝固過程中，不同成分體系下強(qiáng)化相的形核時(shí)間存在明顯差異。模擬結(jié)果顯示，Mg-8at%Gd-3at%Y合金體系的強(qiáng)化相形核時(shí)間相對較晚，當(dāng)溫度降低到約720K時(shí)才開始出現(xiàn)明顯的形核現(xiàn)象；而Mg-12at%Gd-6at%Y合金體系在溫度降至750K左右就已經(jīng)有大量晶核形成。這是因?yàn)槿苜|(zhì)原子的含量直接影響形核的驅(qū)動(dòng)力。溶質(zhì)原子含量越高，液態(tài)合金的自由能與固態(tài)合金自由能的差值越大，形核的驅(qū)動(dòng)力也就越大，從而使得形核更容易發(fā)生，形核時(shí)間提前。在強(qiáng)化相種類方面，通過模擬過程中原子的聚集和排列情況以及對最終凝固結(jié)構(gòu)的分析發(fā)現(xiàn)，三種成分體系中均形成了Mg3(Gd,Y)和Mg24Y5等強(qiáng)化相。然而，不同成分體系下強(qiáng)化相的相對含量有所不同。在Mg-8at%Gd-3at%Y合金體系中，Mg3(Gd,Y)相的含量相對較低，約占強(qiáng)化相總體積的30%，而Mg24Y5相的含量相對較高，約占70%；隨著Gd和Y含量的增加，在Mg-12at%Gd-6at%Y合金體系中，Mg3(Gd,Y)相的含量顯著增加，達(dá)到了約50%，Mg24Y5相的含量則下降至約50%。這是由于Gd和Y含量的變化改變了合金中原子的化學(xué)勢和擴(kuò)散行為。當(dāng)Gd和Y含量較低時(shí)，Y原子更容易與Mg原子結(jié)合形成Mg24Y5相；而隨著Gd和Y含量的增加，Gd原子的作用逐漸凸顯，Gd原子與Mg原子結(jié)合形成Mg3(Gd,Y)相的幾率增大，導(dǎo)致Mg3(Gd,Y)相的含量增加。強(qiáng)化相的尺寸和分布也因成分體系的不同而有所差異。在Mg-8at%Gd-3at%Y合金體系中，強(qiáng)化相尺寸相對較小且分布較為均勻。通過對模擬結(jié)果中強(qiáng)化相尺寸的統(tǒng)計(jì)分析，其平均尺寸約為5nm，尺寸分布的標(biāo)準(zhǔn)差較小，表明強(qiáng)化相尺寸的離散度較小。這是因?yàn)檩^低的溶質(zhì)含量使得形核數(shù)量相對較少，但晶核的生長受到溶質(zhì)擴(kuò)散的限制相對較小，導(dǎo)致強(qiáng)化相在生長過程中能夠較為均勻地發(fā)展。而在Mg-12at%Gd-6at%Y合金體系中，強(qiáng)化相尺寸較大且分布不均勻。平均尺寸增大到約8nm，尺寸分布的標(biāo)準(zhǔn)差較大，存在部分尺寸較大的強(qiáng)化相團(tuán)聚現(xiàn)象。這是由于較高的溶質(zhì)含量促進(jìn)了大量晶核的形成，在晶核生長過程中，溶質(zhì)原子的擴(kuò)散距離較長，容易導(dǎo)致晶核之間的合并和團(tuán)聚，從而形成尺寸較大且分布不均勻的強(qiáng)化相。進(jìn)一步分析不同成分體系下強(qiáng)化相形成差異的原因，除了上述溶質(zhì)含量對形核驅(qū)動(dòng)力、原子化學(xué)勢和擴(kuò)散行為的影響外，溶質(zhì)原子之間的相互作用以及溶質(zhì)原子與鎂基體之間的界面能也起到了重要作用。Gd和Y原子之間的相互作用會(huì)影響它們在鎂基體中的分布和聚集方式。在高含量的Gd和Y合金體系中，Gd和Y原子之間的相互作用可能更強(qiáng)，更容易形成穩(wěn)定的原子團(tuán)簇，這些原子團(tuán)簇不僅作為形核核心促進(jìn)了更多晶核的形成，還在晶核生長過程中影響原子的擴(kuò)散和排列，導(dǎo)致強(qiáng)化相的尺寸和分布發(fā)生變化。同時(shí)，溶質(zhì)原子與鎂基體之間的界面能也會(huì)隨著成分的改變而變化。不同的界面能會(huì)影響晶核與鎂基體之間的界面穩(wěn)定性，進(jìn)而影響晶核的生長速率和形態(tài)。在成分體系變化時(shí)，界面能的改變可能導(dǎo)致晶核生長的各向異性發(fā)生變化，從而使得強(qiáng)化相的尺寸和分布呈現(xiàn)出不同的特征。4.2影響強(qiáng)化相形成的因素4.2.1溫度對強(qiáng)化相形成的影響溫度在高強(qiáng)度稀土鎂合金強(qiáng)化相形成過程中扮演著關(guān)鍵角色，其對強(qiáng)化相的形核率、生長速度和形態(tài)有著顯著且復(fù)雜的影響。通過精心設(shè)計(jì)并實(shí)施不同溫度條件下的凝固過程模擬，能夠深入剖析溫度對強(qiáng)化相形成的具體作用機(jī)制。在模擬過程中，設(shè)置了三個(gè)具有代表性的溫度點(diǎn)：600K、700K和800K，并保持其他模擬參數(shù)一致，以確保溫度成為影響強(qiáng)化相形成的單一變量。從形核率角度來看，模擬結(jié)果顯示，隨著溫度的降低，強(qiáng)化相的形核率呈現(xiàn)明顯上升趨勢。在800K時(shí)，形核率相對較低，單位時(shí)間內(nèi)形成的晶核數(shù)量較少；當(dāng)溫度降至700K時(shí)，形核率顯著提高，晶核形成速度加快；而在600K時(shí)，形核率達(dá)到最高。這是因?yàn)楦鶕?jù)經(jīng)典形核理論，形核驅(qū)動(dòng)力與過冷度密切相關(guān)，過冷度越大，形核驅(qū)動(dòng)力越大，形核率也就越高。隨著溫度降低，合金液的過冷度增大，原子的擴(kuò)散速率相對降低，原子在局部區(qū)域的聚集變得更加容易，從而增加了形核的幾率，提高了形核率。強(qiáng)化相的生長速度同樣受到溫度的顯著影響。在高溫階段，如800K時(shí)，原子具有較高的動(dòng)能，擴(kuò)散速度較快，強(qiáng)化相的生長速度相對較快。然而，隨著溫度降低至700K和600K，原子擴(kuò)散速率逐漸減慢，強(qiáng)化相的生長速度也隨之下降。這是因?yàn)閺?qiáng)化相的生長依賴于原子的擴(kuò)散和遷移，溫度降低導(dǎo)致原子的熱運(yùn)動(dòng)減弱，原子從液相中擴(kuò)散到晶核表面并排列到晶格結(jié)構(gòu)中的速度變慢，從而限制了強(qiáng)化相的生長速度。通過對模擬過程中強(qiáng)化相尺寸隨時(shí)間變化的分析，可以清晰地觀察到這種生長速度與溫度的依賴關(guān)系。在800K時(shí)，強(qiáng)化相在較短時(shí)間內(nèi)就能達(dá)到較大尺寸；而在600K時(shí)，強(qiáng)化相的生長較為緩慢，在相同時(shí)間內(nèi)尺寸相對較小。溫度還對強(qiáng)化相的形態(tài)產(chǎn)生重要影響。在高溫800K下，由于原子擴(kuò)散速度快，強(qiáng)化相在生長過程中更容易沿著各個(gè)方向均勻生長，因此傾向于形成較為規(guī)則的球形或近球形形態(tài)。隨著溫度降低到700K，原子擴(kuò)散速度有所減慢，強(qiáng)化相在不同方向上的生長速度開始出現(xiàn)差異，導(dǎo)致其形態(tài)逐漸變得不規(guī)則，可能出現(xiàn)樹枝狀或多邊形狀。當(dāng)溫度進(jìn)一步降低至600K時(shí)，原子擴(kuò)散受到更大限制，強(qiáng)化相的生長變得更加不均勻，形態(tài)更加復(fù)雜多樣，可能出現(xiàn)分枝狀或團(tuán)聚狀。這種溫度對強(qiáng)化相形態(tài)的影響是由于不同溫度下原子的擴(kuò)散各向異性以及晶核生長過程中的競爭作用導(dǎo)致的。在高溫時(shí)，原子的快速擴(kuò)散使得晶核在各個(gè)方向上的生長條件較為相似，從而形成規(guī)則形態(tài)；而在低溫下，原子擴(kuò)散的各向異性以及晶核之間對原子的競爭吸附作用，使得強(qiáng)化相在生長過程中出現(xiàn)不均勻性，導(dǎo)致形態(tài)的多樣化。4.2.2合金元素含量的作用合金元素含量在高強(qiáng)度稀土鎂合金強(qiáng)化相形成過程中發(fā)揮著至關(guān)重要的作用，其對強(qiáng)化相的數(shù)量、尺寸和分布有著顯著的影響。通過精確控制和改變稀土元素或其他合金元素的含量，利用模擬手段深入分析其對強(qiáng)化相形成的具體作用機(jī)制，對于優(yōu)化合金性能具有重要意義。在模擬研究中，以Mg-Gd-Y合金為研究對象，系統(tǒng)地改變Gd和Y元素的含量。當(dāng)Gd元素含量從8at%增加到12at%，Y元素含量從3at%增加到6at%時(shí)，模擬結(jié)果顯示，強(qiáng)化相的數(shù)量顯著增加。這是因?yàn)楹辖鹪睾康脑黾樱沟煤辖痼w系中形成強(qiáng)化相的原子數(shù)量增多，溶質(zhì)原子的富集區(qū)域也相應(yīng)增加，從而為強(qiáng)化相的形核提供了更多的核心，促進(jìn)了更多晶核的形成，導(dǎo)致強(qiáng)化相數(shù)量上升。例如，在較高的Gd和Y含量下，模擬中觀察到在鎂基體中出現(xiàn)了更多的溶質(zhì)原子團(tuán)簇，這些團(tuán)簇成為強(qiáng)化相形核的優(yōu)先位置，使得形核數(shù)量明顯增多。合金元素含量的變化對強(qiáng)化相尺寸也有明顯影響。隨著Gd和Y含量的增加，強(qiáng)化相的平均尺寸增大。在較低的合金元素含量下，由于形核數(shù)量相對較多，晶核在生長過程中競爭原子資源，導(dǎo)致每個(gè)晶核獲得的原子數(shù)量相對較少，生長受到一定限制，強(qiáng)化相尺寸較小。而當(dāng)合金元素含量升高時(shí)，雖然形核數(shù)量增加，但由于原子供應(yīng)相對充足，晶核有更多的原子可供生長，使得強(qiáng)化相能夠生長到更大的尺寸。通過對不同合金元素含量下強(qiáng)化相尺寸的統(tǒng)計(jì)分析，發(fā)現(xiàn)強(qiáng)化相的平均尺寸與合金元素含量呈現(xiàn)正相關(guān)關(guān)系。在強(qiáng)化相分布方面，合金元素含量的改變同樣產(chǎn)生顯著影響。在低含量合金元素的情況下，強(qiáng)化相在鎂基體中分布相對較為均勻。這是因?yàn)樾魏藬?shù)量相對較少，晶核在生長過程中相互之間的影響較小，能夠較為均勻地在基體中生長和分布。然而，隨著合金元素含量的增加，強(qiáng)化相的分布變得不均勻。高含量的合金元素導(dǎo)致形核數(shù)量大幅增加，在晶核生長過程中，由于原子擴(kuò)散的不均勻性以及晶核之間的相互作用，部分區(qū)域的晶核生長較快，形成較大尺寸的強(qiáng)化相團(tuán)聚體，而部分區(qū)域的強(qiáng)化相則相對較小且分布稀疏，從而導(dǎo)致強(qiáng)化相分布的不均勻性增加。例如，在高合金元素含量的模擬中，可以明顯觀察到在某些區(qū)域出現(xiàn)了強(qiáng)化相的聚集現(xiàn)象，形成了尺寸較大的強(qiáng)化相團(tuán)簇，而在其他區(qū)域強(qiáng)化相的分布則較為分散。合金元素之間的相互作用也會(huì)對強(qiáng)化相的形成產(chǎn)生影響。Gd和Y元素之間可能存在一定的相互吸引作用，在合金凝固過程中，這種相互作用使得它們更容易聚集在一起，形成穩(wěn)定的原子團(tuán)簇。這些原子團(tuán)簇不僅作為形核核心促進(jìn)了更多晶核的形成，還在晶核生長過程中影響原子的擴(kuò)散和排列。例如，當(dāng)Gd和Y元素含量發(fā)生變化時(shí)，它們之間的相互作用強(qiáng)度也可能改變，從而影響原子團(tuán)簇的形成和穩(wěn)定性，進(jìn)而對強(qiáng)化相的數(shù)量、尺寸和分布產(chǎn)生影響。此外，合金元素與鎂基體之間的界面能也會(huì)隨著合金元素含量的改變而變化。不同的界面能會(huì)影響晶核與鎂基體之間的界面穩(wěn)定性，進(jìn)而影響晶核的生長速率和形態(tài)，最終導(dǎo)致強(qiáng)化相的尺寸和分布發(fā)生變化。4.3強(qiáng)化相形成機(jī)制探討4.3.1基于模擬結(jié)果的形核機(jī)制分析根據(jù)模擬過程中對原子聚集和排列情況的詳細(xì)觀察與分析，高強(qiáng)度稀土鎂合金中強(qiáng)化相的形核機(jī)制主要涉及均勻形核和非均勻形核兩種方式，且溶質(zhì)原子的分布和相互作用在其中起到了關(guān)鍵作用。在均勻形核方面，當(dāng)合金處于液態(tài)且溫度逐漸降低時(shí)，根據(jù)經(jīng)典形核理論，體系中的原子會(huì)發(fā)生熱運(yùn)動(dòng)和擴(kuò)散。在液態(tài)合金中，由于原子的隨機(jī)熱運(yùn)動(dòng)，會(huì)在局部區(qū)域出現(xiàn)原子的短暫聚集，形成一些微小的原子團(tuán)簇。這些原子團(tuán)簇的能量狀態(tài)高于周圍液態(tài)原子，處于不穩(wěn)定狀態(tài)。隨著溫度進(jìn)一步降低，原子的擴(kuò)散速率減小，原子團(tuán)簇的穩(wěn)定性增加。當(dāng)原子團(tuán)簇的尺寸達(dá)到一定臨界值時(shí)，就形成了穩(wěn)定的晶核，即發(fā)生了均勻形核。在Mg-Gd-Y合金模擬中，通過對原子運(yùn)動(dòng)軌跡和原子團(tuán)簇形成過程的跟蹤發(fā)現(xiàn)，在液態(tài)合金中，Mg、Gd和Y原子會(huì)隨機(jī)聚集形成一些小的原子團(tuán)。在某一時(shí)刻，當(dāng)溫度降低到一定程度時(shí)，部分原子團(tuán)的尺寸逐漸增大，達(dá)到臨界尺寸后，這些原子團(tuán)就成為了均勻形核的晶核。例如，在模擬的某一溫度區(qū)間內(nèi)，通過對原子分布的統(tǒng)計(jì)分析，發(fā)現(xiàn)當(dāng)溫度降低到750K左右時(shí)，均勻形核的晶核數(shù)量開始明顯增加。這是因?yàn)樵谠摐囟认拢拥臄U(kuò)散速率和能量狀態(tài)使得原子團(tuán)簇更容易達(dá)到臨界尺寸，從而促進(jìn)了均勻形核的發(fā)生。非均勻形核在強(qiáng)化相形成過程中也起著重要作用。在實(shí)際合金體系中，存在著各種晶體缺陷，如位錯(cuò)、空位、晶界等。這些晶體缺陷處的原子排列不規(guī)則，能量較高，為形核提供了有利的位置。溶質(zhì)原子在晶體缺陷處的偏聚現(xiàn)象較為明顯。以位錯(cuò)為例，位錯(cuò)是晶體中的線缺陷，其周圍存在著應(yīng)力場。溶質(zhì)原子由于與基體原子的尺寸差異，會(huì)在應(yīng)力場的作用下向位錯(cuò)處偏聚。在Mg-Gd-Y合金中，Gd和Y原子相對Mg原子尺寸較大，更容易在位錯(cuò)處偏聚。這些偏聚的溶質(zhì)原子會(huì)降低形核的自由能壘，使得在晶體缺陷處更容易形成晶核。模擬結(jié)果顯示，在含有位錯(cuò)的區(qū)域，溶質(zhì)原子的濃度明顯高于其他區(qū)域，且在這些區(qū)域優(yōu)先形成了強(qiáng)化相的晶核。通過對晶核周圍原子排列和溶質(zhì)原子分布的分析，進(jìn)一步證實(shí)了非均勻形核的機(jī)制。在晶界處，由于晶界是不同晶粒之間的過渡區(qū)域，原子排列紊亂，能量較高，也為溶質(zhì)原子的偏聚和形核提供了條件。模擬中觀察到在晶界附近，溶質(zhì)原子會(huì)聚集形成原子團(tuán)簇，這些原子團(tuán)簇在適當(dāng)?shù)臈l件下會(huì)發(fā)展成為晶核，促進(jìn)非均勻形核的發(fā)生。溶質(zhì)原子之間的相互作用對形核過程也有顯著影響。在Mg-Gd-Y合金中，Gd和Y原子之間存在一定的相互吸引作用。這種相互作用使得它們在液態(tài)合金中更容易聚集在一起，形成相對穩(wěn)定的原子團(tuán)簇。這些原子團(tuán)簇不僅尺寸較大，而且穩(wěn)定性較高，更容易達(dá)到臨界形核尺寸，從而促進(jìn)形核的發(fā)生。通過對模擬過程中原子間相互作用能的計(jì)算和分析，發(fā)現(xiàn)Gd和Y原子之間的相互作用能為負(fù)值，表明它們之間存在吸引作用。這種吸引作用導(dǎo)致Gd和Y原子在液態(tài)合金中形成了一些較大的原子團(tuán)簇，這些原子團(tuán)簇成為了形核的核心，增加了形核的幾率。此外，溶質(zhì)原子與鎂基體原子之間的相互作用也會(huì)影響形核過程。溶質(zhì)原子與鎂基體原子之間的相互作用強(qiáng)度會(huì)影響溶質(zhì)原子在基體中的擴(kuò)散和聚集行為，進(jìn)而影響形核的驅(qū)動(dòng)力和形核速率。4.3.2生長機(jī)制的理論解釋強(qiáng)化相的生長機(jī)制涉及原子擴(kuò)散和界面能等多個(gè)理論，這些因素相互作用，共同決定了強(qiáng)化相在合金中的生長過程和最終形態(tài)。從原子擴(kuò)散理論角度來看，強(qiáng)化相的生長主要依賴于原子在液態(tài)合金中的擴(kuò)散和遷移。在合金凝固過程中，當(dāng)晶核形成后，晶核周圍的液態(tài)原子會(huì)不斷向晶核表面擴(kuò)散，并排列到晶核的晶格結(jié)構(gòu)中，從而使晶核逐漸長大。原子擴(kuò)散的驅(qū)動(dòng)力主要來自于濃度梯度和溫度梯度。在晶核周圍，由于原子不斷向晶核表面擴(kuò)散，導(dǎo)致晶核周圍液態(tài)合金中的溶質(zhì)原子濃度降低，形成了濃度梯度。這種濃度梯度促使液態(tài)合金中的原子向晶核表面擴(kuò)散，以平衡濃度差異。例如，在Mg-Gd-Y合金中，當(dāng)Mg3(Gd,Y)強(qiáng)化相晶核形成后，周圍液態(tài)合金中的Mg、Gd和Y原子會(huì)在濃度梯度的作用下向晶核表面擴(kuò)散。通過對模擬過程中原子擴(kuò)散路徑的分析，可以清晰地觀察到原子從液態(tài)合金中向晶核表面遷移的過程。同時(shí)，溫度梯度也會(huì)影響原子的擴(kuò)散。在凝固過程中，合金體系存在溫度差，高溫區(qū)域的原子具有較高的動(dòng)能，更容易向低溫區(qū)域擴(kuò)散。這種溫度梯度驅(qū)動(dòng)的原子擴(kuò)散也為強(qiáng)化相的生長提供了動(dòng)力。界面能在強(qiáng)化相生長過程中起著重要的作用。晶核與周圍液態(tài)合金之間存在著界面，界面能是指單位面積界面上的能量。界面能的大小與晶核的形狀和原子排列方式密切相關(guān)。在強(qiáng)化相生長過程中，晶核會(huì)傾向于以降低界面能的方式生長。根據(jù)晶體生長理論，當(dāng)晶核的界面能各向同性時(shí)，晶核會(huì)以球形生長，因?yàn)榍蛐尉哂凶钚〉谋砻娣e與體積比，能夠使界面能最小化。然而，在實(shí)際合金體系中，晶核的界面能往往是各向異性的，這是由于晶核不同晶面的原子排列方式和原子間相互作用不同導(dǎo)致的。在這種情況下，晶核會(huì)沿著界面能較低的晶面優(yōu)先生長，從而形成特定的晶體形態(tài)。在Mg-Gd-Y合金中，通過對強(qiáng)化相晶核不同晶面的界面能計(jì)算發(fā)現(xiàn)，某些晶面的界面能較低。模擬結(jié)果顯示，強(qiáng)化相在生長過程中，會(huì)沿著這些界面能較低的晶面優(yōu)先生長，形成了具有特定形狀的晶體，如樹枝狀或多邊形狀。這種由于界面能各向異性導(dǎo)致的晶體生長方式，使得強(qiáng)化相的形態(tài)更加復(fù)雜多樣。此外，強(qiáng)化相生長過程中還存在著溶質(zhì)原子的再分配現(xiàn)象。隨著強(qiáng)化相的生長，溶質(zhì)原子在晶核與液態(tài)合金之間不斷進(jìn)行分配。在晶核生長初期，由于晶核表面的原子對溶質(zhì)原子具有一定的吸附作用，使得晶核表面附近的溶質(zhì)原子濃度較高。隨著生長的進(jìn)行，溶質(zhì)原子在濃度梯度的作用下，會(huì)逐漸向液態(tài)合金中擴(kuò)散，導(dǎo)致晶核表面的溶質(zhì)原子濃度逐漸降低。這種溶質(zhì)原子的再分配過程會(huì)影響強(qiáng)化相的生長速率和成分均勻性。在Mg-Gd-Y合金模擬中，通過對不同生長階段強(qiáng)化相成分的分析發(fā)現(xiàn)，在生長初期，強(qiáng)化相中的Gd和Y含量相對較高，隨著生長的進(jìn)行，Gd和Y含量逐漸降低。這表明在生長過程中，溶質(zhì)原子發(fā)生了再分配，影響了強(qiáng)化相的成分和性能。同時(shí)，溶質(zhì)原子的再分配還會(huì)導(dǎo)致液態(tài)合金中成分的不均勻性，進(jìn)而影響后續(xù)晶核的形成和生長。五、強(qiáng)化相原子排列的模擬分析5.1強(qiáng)化相的晶體結(jié)構(gòu)與原子排列特征5.1.1常見強(qiáng)化相的晶體結(jié)構(gòu)類型在高強(qiáng)度稀土鎂合金中，Mg?RE和Mg??RE?是兩種典型且常見的強(qiáng)化相，它們各自具有獨(dú)特的晶體結(jié)構(gòu)類型，對合金的性能產(chǎn)生著重要影響。Mg?RE強(qiáng)化相通常呈現(xiàn)出六方晶系結(jié)構(gòu)，屬于DO??型晶體結(jié)構(gòu)。在這種結(jié)構(gòu)中，原子排列具有高度的規(guī)律性。從晶胞結(jié)構(gòu)來看，其晶胞參數(shù)具有特定的數(shù)值，a軸和b軸長度相等，且夾角為120°，c軸與a、b軸垂直。以Mg?Gd為例，其晶胞參數(shù)a=0.5203nm，c=0.8504nm。在晶胞內(nèi)，原子按照特定的方式排列，Mg原子和Gd原子通過金屬鍵相互連接，形成穩(wěn)定的晶體結(jié)構(gòu)。這種晶體結(jié)構(gòu)使得Mg?RE強(qiáng)化相具有較高的穩(wěn)定性和熱穩(wěn)定性，能夠在高溫環(huán)境下保持其結(jié)構(gòu)完整性。其獨(dú)特的原子排列方式為合金提供了良好的高溫性能，在高溫下，能夠有效阻礙位錯(cuò)運(yùn)動(dòng)，提高合金的強(qiáng)度和抗蠕變性能。Mg??RE?強(qiáng)化相則屬于正交晶系結(jié)構(gòu)，為cbco型晶體結(jié)構(gòu)。在該結(jié)構(gòu)中，晶胞參數(shù)a、b、c各不相等，且夾角均為90°。以Mg??Y?為例，其晶胞參數(shù)a=1.447nm，b=1.014nm，c=1.221nm。在晶胞內(nèi)部，Mg原子和Y原子按照特定的排列規(guī)律分布。Mg??Y?相的形成與合金成分和凝固條件密切相關(guān)。在合適的合金成分和凝固過程中，Mg和Y原子通過擴(kuò)散和聚集，逐漸形成這種具有特定原子排列的強(qiáng)化相。這種強(qiáng)化相在合金中起到了重要的強(qiáng)化作用，其原子排列方式使得位錯(cuò)在運(yùn)動(dòng)過程中難以穿過，從而提高了合金的強(qiáng)度和硬度。這些常見強(qiáng)化相的晶體結(jié)構(gòu)類型是由合金元素的原子半徑、電負(fù)性以及原子間的相互作用等多種因素共同決定的。合金元素之間的原子半徑差異會(huì)影響原子在晶體結(jié)構(gòu)中的排列方式，電負(fù)性的不同則決定了原子間化學(xué)鍵的性質(zhì)。在Mg-RE系合金中，Mg、RE等原子之間的相互作用，包括靜電相互作用、電子云的重疊等，共同促使了這些特定晶體結(jié)構(gòu)的形成。不同的晶體結(jié)構(gòu)類型導(dǎo)致強(qiáng)化相具有不同的物理和化學(xué)性質(zhì)，進(jìn)而對合金的力學(xué)性能、物理性能和化學(xué)性能產(chǎn)生顯著影響。例如，Mg?RE強(qiáng)化相的六方晶系結(jié)構(gòu)使其在某些晶向具有較高的強(qiáng)度和穩(wěn)定性，而Mg??RE?強(qiáng)化相的正交晶系結(jié)構(gòu)則賦予其在不同方向上的力學(xué)性能各向異性。這些特性對于理解合金的性能和優(yōu)化合金設(shè)計(jì)具有重要意義。5.1.2原子在晶格中的排列方式通過分子動(dòng)力學(xué)模擬，清晰地展示了原子在Mg?RE和Mg??RE?強(qiáng)化相晶格中的具體排列方式。在Mg?RE強(qiáng)化相（以Mg?Gd為例）的六方晶系DO??型結(jié)構(gòu)中，原子排列具有高度的對稱性和規(guī)律性。從晶胞的角度來看，在底面（a-b平面）上，Mg原子和Gd原子交替排列，形成六邊形的原子層。每個(gè)Mg原子周圍有6個(gè)Gd原子作為最近鄰，形成了一個(gè)穩(wěn)定的配位結(jié)構(gòu)。在c軸方向上，這些原子層按照一定的順序堆疊，相鄰原子層之間的原子位置相互錯(cuò)開，以達(dá)到最緊密堆積的狀態(tài)。通過模擬得到的原子排列快照，可以直觀地觀察到這種排列方式。在快照中，不同顏色的原子代表不同的元素，Mg原子和Gd原子在晶格中的位置一目了然。從原子間距離的分析結(jié)果來看，Mg-Gd原子間的距離約為0.30-0.32nm，這種距離既保證了原子間的相互作用強(qiáng)度，又使得晶體結(jié)構(gòu)具有一定的柔韌性。這種排列方式使得Mg?Gd強(qiáng)化相具有較高的結(jié)構(gòu)穩(wěn)定性，在高溫下能夠有效阻礙位錯(cuò)的運(yùn)動(dòng)。當(dāng)位錯(cuò)在晶體中運(yùn)動(dòng)時(shí)，遇到這種規(guī)則排列的原子結(jié)構(gòu)，需要克服較大的能量障礙，從而提高了合金的強(qiáng)度和抗蠕變性能。在Mg??RE?強(qiáng)化相（以Mg??Y?為例）的正交晶系cbco型結(jié)構(gòu)中，原子排列方式與Mg?RE強(qiáng)化相有所不同。在晶胞內(nèi)，Mg原子和Y原子分布在不同的晶格位置上。沿著a軸方向，Mg原子和Y原子交替排列形成鏈狀結(jié)構(gòu)；在b軸和c軸方向上，這些鏈狀結(jié)構(gòu)相互交織，形成復(fù)雜的三維網(wǎng)絡(luò)結(jié)構(gòu)。模擬結(jié)果顯示，在這種結(jié)構(gòu)中，Mg-Y原子間的距離在不同方向上略有差異，在a軸方向上約為0.28-0