陳水洲王曉書(shū)歐陽(yáng)求保張瑞

(1.上海大學(xué)計(jì)算機(jī)工程與科學(xué)學(xué)院,上海200444;2.上海交通大學(xué)材料科學(xué)與工程學(xué)院金屬基復(fù)合材料國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200240;3.上海大學(xué)材料基因組工程研究院材料信息與數(shù)據(jù)科學(xué)中心,上海200444;4.之江實(shí)驗(yàn)室,浙江杭州311100)

碳化硅顆粒(SiC particle,SiCp)增強(qiáng)鋁基復(fù)合材料具有高模量、高強(qiáng)度、低密度、低膨脹等優(yōu)異的綜合性能,在航空航天、武器裝備及交通運(yùn)輸?shù)阮I(lǐng)域有著廣泛的應(yīng)用.本工作采用原位生長(zhǎng)的方式,在SiCp表面原位自生碳納米管(carbon nanotubes,CNTs),制備出SiCp(CNT)微納復(fù)合增強(qiáng)體,并以7075Al粉末為基體,通過(guò)真空熱壓燒結(jié)工藝制備SiCp(CNT)/7075Al復(fù)合材料.通過(guò)熱處理工藝參數(shù)調(diào)節(jié)和成分設(shè)計(jì),對(duì)復(fù)合材料的力學(xué)性能進(jìn)行了調(diào)控.由于熱處理工藝參數(shù)較為復(fù)雜,包括固溶時(shí)間、固溶溫度、時(shí)效時(shí)間以及時(shí)效溫度等多個(gè)因素,因此為了實(shí)現(xiàn)復(fù)合材料力學(xué)性能的最優(yōu)化,需要進(jìn)行大量的工藝探索.機(jī)器學(xué)習(xí)在交叉學(xué)科領(lǐng)域發(fā)揮著極大的作用[1-3].通過(guò)機(jī)器學(xué)習(xí)可以探索數(shù)據(jù)的特性、預(yù)測(cè)性能,分析已有的實(shí)驗(yàn)數(shù)據(jù),篩選出對(duì)實(shí)驗(yàn)結(jié)果影響最大的特征,從而有助于大大縮短整個(gè)實(shí)驗(yàn)周期.

基于此,本工作將機(jī)器學(xué)習(xí)方法應(yīng)用到鋁基復(fù)合熱處理工藝中,對(duì)SiCp(CNT)/7075Al鋁基復(fù)合材料的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了分析,分別用嶺回歸、正則化線性回歸、彈性網(wǎng)絡(luò)回歸、支持向量機(jī)回歸(support vector regression,SVR)等多種機(jī)器學(xué)習(xí)方法,對(duì)鋁基復(fù)合材料的抗拉強(qiáng)度和延伸率兩種性能進(jìn)行了預(yù)測(cè),完成了熱處理工藝參數(shù)的優(yōu)化,并在此基礎(chǔ)上進(jìn)行了實(shí)際實(shí)驗(yàn)驗(yàn)證.

1 鋁基復(fù)合材料

顆粒增強(qiáng)鋁基復(fù)合材料具有高模量、高強(qiáng)度、低密度、低膨脹和良好的耐磨性等優(yōu)異性能,被廣泛應(yīng)用于航空航天和交通運(yùn)輸?shù)阮I(lǐng)域[4-8].目前常用的增強(qiáng)體主要包括碳化硅(SiCp)、氧化鋁(Al2O3)和碳化硼(B4C)等微米級(jí)陶瓷顆粒[9-13],其中SiCp對(duì)鋁基復(fù)合材料顯著的強(qiáng)化效果已在2xxx、5xxx和6xxx等鋁合金基體中得到較好驗(yàn)證[14-16].力學(xué)性能更優(yōu)異的基體合金往往會(huì)產(chǎn)生更強(qiáng)的復(fù)合材料,因此通常選擇高強(qiáng)度鋁合金作為結(jié)構(gòu)材料用顆粒以及增強(qiáng)鋁基復(fù)合材料的基體[17].大多數(shù)Al-Zn-Mg-Cu合金的強(qiáng)度高于2xxx、5xxx和6xxx鋁合金,如7075Al、7085Al和7055Al合金.因此,基于7系高強(qiáng)度鋁合金開(kāi)發(fā)高強(qiáng)度鋁基復(fù)合材料的研究受到了廣泛關(guān)注.

然而,SiCp在7xxx鋁合金中的強(qiáng)化效果遠(yuǎn)不如2xxx鋁合金[18-21],甚至SiCp在SiC/Al-Zn-Mg-(Cu)復(fù)合材料中表現(xiàn)出弱化作用,即復(fù)合材料的強(qiáng)度低于鋁合金[22-24].Lee等[18]發(fā)現(xiàn),增強(qiáng)體與基體間的不匹配會(huì)導(dǎo)致增強(qiáng)體周?chē)a(chǎn)生嚴(yán)重的應(yīng)力集中,從而使該區(qū)域的基體過(guò)早失效,最終導(dǎo)致復(fù)合材料屈服強(qiáng)度的降低.Min等[25]認(rèn)為,弱的基體和界面脫黏可能會(huì)導(dǎo)致SiCp的強(qiáng)化效果無(wú)法充分體現(xiàn).因此,制備出致密的復(fù)合材料并緩解SiCp周?chē)膽?yīng)力集中,對(duì)于充分實(shí)現(xiàn)SiCp的強(qiáng)化效果十分重要.另外,Al-Zn-Mg-Cu合金中的主要析出相MgZn2也將對(duì)復(fù)合材料的力學(xué)性能產(chǎn)生顯著影響.Ma等[26]在SiC/7075Al復(fù)合材料中發(fā)現(xiàn),SiCp可以加速時(shí)效過(guò)程.這說(shuō)明SiCp可以為沉淀相提供更多的優(yōu)先成核位點(diǎn),降低固溶體原子擴(kuò)散的活化能.SiCp引入的位錯(cuò)和晶格畸變也會(huì)加速部分析出相的異常生長(zhǎng).此外,Mg元素的添加也可以提升復(fù)合材料的力學(xué)性能[17].

綜上所述,為了制備出力學(xué)性能優(yōu)異的SiCp/7075復(fù)合材料,本工作通過(guò)在SiCp原位自生CNTs的方式制備出SiCp(CNT)微納復(fù)合增強(qiáng)體,旨在緩解增強(qiáng)體顆粒周?chē)膽?yīng)力集中.之后,采用真空熱壓燒結(jié)工藝將其與7075Al基體進(jìn)行復(fù)合,制備出致密的微納復(fù)合構(gòu)型復(fù)合材料.最后,結(jié)合熱處理工藝的優(yōu)化以及Mg元素的添加,制備出強(qiáng)塑性匹配的SiCp(CNT)復(fù)合增強(qiáng)7075鋁基復(fù)合材料SiCp(CNT)/7075Al.同時(shí),為了實(shí)現(xiàn)復(fù)合材料力學(xué)性能的最優(yōu)化,本工作通過(guò)機(jī)器學(xué)習(xí)進(jìn)行復(fù)合材料成分(如Mg元素含量)設(shè)計(jì)、熱處理工藝參數(shù)的探索和設(shè)計(jì),并加以實(shí)驗(yàn)驗(yàn)證,大幅度縮短了實(shí)驗(yàn)周期.本方法可以為開(kāi)發(fā)新的鋁基復(fù)合材料以及鋁基復(fù)合材料工藝設(shè)計(jì)提供技術(shù)支撐.

2 機(jī)器學(xué)習(xí)對(duì)鋁基復(fù)合材料的性能預(yù)測(cè)和逆向設(shè)計(jì)

2.1 數(shù)據(jù)描述和預(yù)處理

鋁基復(fù)合材料的制備過(guò)程如下:首先,采用化學(xué)氣相沉積法在SiCp表面進(jìn)行CNTs的原位生長(zhǎng),制備SiCp(CNT)微納復(fù)合增強(qiáng)體;然后,以SiCp(CNT),Mg粉以及Al粉末作為原材料,采用混粉、真空熱壓及熱擠壓的粉末冶金工藝,制備出致密的SiCp(CNT)/7075Al復(fù)合材料;之后,通過(guò)調(diào)控固溶溫度、固溶時(shí)間、時(shí)效溫度以及時(shí)效時(shí)間4個(gè)熱處理工藝參數(shù)以及Mg粉的添加量,進(jìn)而調(diào)控SiCp(CNT)/7075Al復(fù)合材料的力學(xué)性能;最后,制備出強(qiáng)塑化的鋁基復(fù)合材料.

硬度測(cè)試采用HVS-30P型顯微維氏硬度計(jì),測(cè)試時(shí)的使用載荷為4 903 N,加載時(shí)間為10 s.拉伸性能的測(cè)試樣品取自擠壓后的棒材,然后根據(jù)GB/T 228.1—2010對(duì)棒材拉伸試樣的規(guī)定,通過(guò)機(jī)械加工的方法制得.拉伸測(cè)試在Z100 TEW-100 KN型萬(wàn)能材料拉伸試驗(yàn)機(jī)上進(jìn)行,拉伸速率為0.5 mm/min,標(biāo)距為25 mm.

本實(shí)驗(yàn)共獲得57條樣本數(shù)據(jù).根據(jù)復(fù)合材料的制備過(guò)程提取每條數(shù)據(jù)的特征,包括化學(xué)成分、熱處理工藝和力學(xué)性能.化學(xué)成分包括材料的基體部分(Al91Zn5.2Mg2.3Cu1.5、Mg)和增強(qiáng)體部分(SiCp、CNTs).熱處理工藝主要包括對(duì)樣品進(jìn)行固溶過(guò)程和時(shí)效過(guò)程的工藝參數(shù),即固溶溫度、固溶時(shí)間、時(shí)效溫度以及時(shí)效時(shí)間.力學(xué)性能主要包括抗拉強(qiáng)度和延伸率.57條樣本數(shù)據(jù)經(jīng)過(guò)數(shù)據(jù)預(yù)處理,篩除無(wú)效數(shù)據(jù)后剩余22條.對(duì)每個(gè)特征進(jìn)行統(tǒng)計(jì),結(jié)果如表1所示.

表1 鋁基復(fù)合材料22條樣本數(shù)據(jù)中單個(gè)特征的統(tǒng)計(jì)量Table 1 Numerical statistics of individual features of 22 aluminum matrix composites samples

對(duì)待預(yù)測(cè)的抗拉強(qiáng)度和延伸率的數(shù)據(jù)分布進(jìn)行分析,結(jié)果如圖1所示,其中點(diǎn)的分布越靠近紅色直線,表示越接近正態(tài)分布.顯然,2個(gè)性能指標(biāo)均接近正態(tài)分布,不需要對(duì)其再做進(jìn)一步的處理.

圖1 抗拉強(qiáng)度及延伸率的頻數(shù)分布圖和Q-Q圖Fig.1 Frequency distribution plots and Q-Q plots of tensile strength and elongation

2.2 特征分析

2.2.1 基于Shapley值的特征重要性分析

Shapley值方法用于分析每個(gè)特征對(duì)預(yù)測(cè)目標(biāo)的貢獻(xiàn)程度[27].本工作分析了不同特征對(duì)抗拉強(qiáng)度和延伸率這2個(gè)性能指標(biāo)的影響程度,結(jié)果如圖2所示.圖2(a)和圖2(c)分別表示單個(gè)特征對(duì)抗拉強(qiáng)度和延伸率的Shapley值,即單個(gè)特征在預(yù)測(cè)抗拉強(qiáng)度和延伸率上的貢獻(xiàn),其中的每個(gè)點(diǎn)代表著一個(gè)樣本,點(diǎn)的顏色越紅,表明特征值越大,貢獻(xiàn)值也越大.圖2(b)和圖2(d)分別表示單個(gè)特征對(duì)抗拉強(qiáng)度和延伸率Shapley值計(jì)算后的絕對(duì)值.

從圖2(a)中可以看出:與抗拉強(qiáng)度呈明顯正相關(guān)的特征有Mg元素含量、時(shí)效時(shí)間;與抗拉強(qiáng)度呈明顯負(fù)相關(guān)的特征有CNTs含量、固溶時(shí)間、時(shí)效時(shí)間.從圖2(b)中可以看出:對(duì)抗拉強(qiáng)度預(yù)測(cè)貢獻(xiàn)最大的前4個(gè)特征依次為Mg元素含量、基底材料Al91Zn5.2Mg2.3Cu1.5含量、CNTs含量和時(shí)效時(shí)間.從圖2(c)中可以看出:與延伸率呈正相關(guān)的特征有基底材料Al91Zn5.2Mg2.3Cu1.5含量;與延伸率呈負(fù)相關(guān)的特征有固溶時(shí)間和CNTs含量.從圖2(d)中可以看出:對(duì)延伸率預(yù)測(cè)貢獻(xiàn)最大的前4個(gè)特征依次為基底材料Al91Zn5.2Mg2.3Cu1.5含量、時(shí)效時(shí)間、固溶時(shí)間、CNTs含量.

圖2 不同特征對(duì)抗拉強(qiáng)度及延伸率的Shapley值分析Fig.2 Shapley values analysis of different features with tensile strength and elongation

2.2.2 基于樹(shù)形模型的特征重要性分析

除了使用Shapley值方法,本工作還基于隨機(jī)森林算法對(duì)特征重要性進(jìn)行了分析,結(jié)果如圖3所示.

圖3 基于隨機(jī)森林模型的特征重要性分析Fig.3 Feature importance analysis with Random Forest model

從圖3(a)可以看出,針對(duì)抗拉強(qiáng)度,特征重要性的前4位依次是Mg元素含量、基底材料Al91Zn5.2Mg2.3Cu1.5含量、CNTs含量、時(shí)效時(shí)間.從圖3(b)可以看出,針對(duì)延伸率,特征重要性的前4位依次是基底材料Al91Zn5.2Mg2.3Cu1.5含量、SiCp含量、時(shí)效時(shí)間、固溶時(shí)間.此結(jié)果和2.2.1節(jié)Shapley值分析結(jié)果吻合.

2.3 集成框架

本工作使用多個(gè)機(jī)器學(xué)習(xí)算法構(gòu)建集成框架,如圖4所示.整個(gè)框架的機(jī)器學(xué)習(xí)算法包括:基于線性理論的Ridge回歸、Lasso回歸和ElasticNet回歸;基于支持向量機(jī)原理的SVR;基于最近鄰理論的K-Neighbors回歸;基于樹(shù)形模型的隨機(jī)森林回歸;基于Boosting方法的Xgboost回歸和LightGBM回歸.

圖4 集成多個(gè)機(jī)器學(xué)習(xí)算法的模型調(diào)參和模型選擇框架Fig.4 Framework of model parameters tuning and selection based on the multi-machine learning algorithms

整個(gè)集成框架的實(shí)驗(yàn)流程如下:首先,輸入經(jīng)過(guò)預(yù)處理的樣本數(shù)據(jù);然后,通過(guò)格搜索(GridSearch)的方式進(jìn)行調(diào)參,并通過(guò)字典的方法將完整的最優(yōu)模型(best_estimator_)、調(diào)參的結(jié)果(best_params )以及測(cè)試集上的評(píng)分(test_score)保存下來(lái);最后,將上述結(jié)果輸出為CSV文件,用于后續(xù)的實(shí)際實(shí)驗(yàn)驗(yàn)證和材料逆向設(shè)計(jì).每個(gè)模型的核心超參數(shù)以及簡(jiǎn)要說(shuō)明如表2所示.

表2 機(jī)器學(xué)習(xí)模型核心超參數(shù)及簡(jiǎn)要說(shuō)明Table 2 Machine learning model core hyperparameters and brief descriptions

本工作將上述8個(gè)機(jī)器學(xué)習(xí)算法集成到材料基因工程數(shù)據(jù)平臺(tái),在平臺(tái)前端頁(yè)面上傳待建模數(shù)據(jù)后,輸入不同機(jī)器學(xué)習(xí)算法所需參數(shù),點(diǎn)擊提交任務(wù),平臺(tái)即可返回結(jié)果.具體而言,前端可分為參數(shù)定義部分和結(jié)果展示部分.參數(shù)定義部分要求用戶輸入數(shù)據(jù)的特征部分和預(yù)測(cè)目標(biāo),以及對(duì)應(yīng)模型參數(shù)的開(kāi)始值、結(jié)束值和搜索步長(zhǎng).點(diǎn)擊“提交”按鈕,后臺(tái)即開(kāi)始計(jì)算并返回結(jié)果,同時(shí)顯示不同模型的最佳參數(shù)和評(píng)價(jià)得分.平臺(tái)首頁(yè)如圖5所示,超參優(yōu)化和最優(yōu)模型選擇如圖6所示.

圖5 材料基因工程數(shù)據(jù)平臺(tái)首頁(yè)Fig.5 Portal of Materials Genome Engineering Data Platform

圖6 超參自動(dòng)優(yōu)化和最優(yōu)模型選擇軟件運(yùn)行界面Fig.6 Software interface of hyperparameters automatic optimization and optimal model selection

在訓(xùn)練過(guò)程中,使用K折交叉驗(yàn)證(本實(shí)驗(yàn)取K=10)進(jìn)行模型評(píng)分,即將原始數(shù)據(jù)集隨機(jī)劃分成K份,以其中的K-1份作為訓(xùn)練集,剩下的1份作為測(cè)試集.最后以模型在測(cè)試集上的表現(xiàn)作為評(píng)價(jià)標(biāo)準(zhǔn).本工作使用的評(píng)價(jià)函數(shù)為平均絕對(duì)百分比誤差(mean absolute percentage error,MAPE),

式中:n表示數(shù)據(jù)樣本的數(shù)量;i表示模型在第i個(gè)樣本上的預(yù)測(cè)結(jié)果.

采用1-MAPE計(jì)算評(píng)分,各模型在測(cè)試集上對(duì)抗拉強(qiáng)度和延伸率的預(yù)測(cè)評(píng)分如表3所示.可以看出:在對(duì)延伸率的預(yù)測(cè)上,SVR模型獲得最高評(píng)分;在對(duì)抗拉強(qiáng)度的預(yù)測(cè)上,7個(gè)模型的評(píng)分都超過(guò)了0.986且相差很小.考慮到SVR模型在小數(shù)據(jù)處理上的優(yōu)勢(shì),因此兩個(gè)任務(wù)都選擇SVR回歸模型進(jìn)行逆向設(shè)計(jì),并通過(guò)實(shí)際實(shí)驗(yàn)來(lái)進(jìn)行驗(yàn)證.

表3 不同機(jī)器學(xué)習(xí)模型在兩個(gè)性能預(yù)測(cè)上的評(píng)分(1-MAPE)Table 3 Scores of different machine learning models on two performance predictions(1-MAPE)

2.4 模型驗(yàn)證和逆向設(shè)計(jì)

本工作采用逆向設(shè)計(jì)方法來(lái)驗(yàn)證本模型的有效性.逆向設(shè)計(jì)的目的是對(duì)化學(xué)成分和熱處理工藝進(jìn)行優(yōu)化,制備出高性能的鋁基復(fù)合材料.首先,根據(jù)2.2節(jié)的特征重要性分析,預(yù)選出基底材料Al91Zn5.2Mg2.3Cu1.5含量、Mg元素含量、SiCp含量、CNTs含量、固溶時(shí)間和時(shí)效時(shí)間6個(gè)特征.然后,根據(jù)7075鋁基復(fù)合材料的設(shè)計(jì)需求,主要調(diào)節(jié)基底材料Al91Zn5.2Mg2.3Cu1.5和Mg元素的含量,并選擇固溶時(shí)間和時(shí)效時(shí)間兩個(gè)工藝特征做細(xì)粒度建模,列舉所有特征的優(yōu)化區(qū)間或取值,得到如表4所示的特征取值范圍和步長(zhǎng).之后,對(duì)上述結(jié)果進(jìn)行排列組合,產(chǎn)生眾多虛擬樣本點(diǎn).最后,利用2.3節(jié)選出的最優(yōu)模型預(yù)測(cè)虛擬樣本的抗拉強(qiáng)度和延伸率.結(jié)合對(duì)鋁基復(fù)合材料抗拉強(qiáng)度和延伸率的綜合指標(biāo)要求,從眾多虛擬樣本中篩選出3條性能表現(xiàn)較好的虛擬樣本,其對(duì)應(yīng)的成分、工藝和預(yù)測(cè)性能如表5所示.

表4 鋁基復(fù)合材料逆向設(shè)計(jì)樣本取值范圍Table 4 Range of values for the inverse design samples of aluminum matrix composites

表5 鋁基復(fù)合材料逆向設(shè)計(jì)得到的最優(yōu)虛擬樣本Table 5 Optimal virtual samples obtained by inverse design of aluminum matrix composites

根據(jù)上述模型預(yù)測(cè)結(jié)果選出具有較好抗拉強(qiáng)度和延伸率的虛擬樣本.在考慮抗拉強(qiáng)度最優(yōu)的前提下選擇虛擬樣本2,并實(shí)際制作鋁基復(fù)合材料,其中基底材料Al91Zn5.2Mg2.3Cu1.5占比84%、Mg元素占比1%、SiCp占比14.5%、CNTs占比0.5%.依次通過(guò)時(shí)長(zhǎng)40 min、溫度為470°C的固溶操作和時(shí)長(zhǎng)15 h、溫度為120°C的時(shí)效操作,獲得該鋁基復(fù)合材料的性能測(cè)試結(jié)果,如圖7所示.可以看出:該復(fù)合材料的抗拉強(qiáng)度和延伸率在真實(shí)實(shí)驗(yàn)條件下分別達(dá)到了647.0 MPa和3.31%,而模型預(yù)測(cè)結(jié)果為617.48 Mpa和2.98%;2個(gè)物理性能的MAPE依次為4.56%和9.97%,均滿足了材料制備的誤差范圍要求.這說(shuō)明利用機(jī)器學(xué)習(xí)預(yù)測(cè)建模并逆向設(shè)計(jì)出性能最優(yōu)的材料,是切實(shí)可行的.

圖7 SiCp(0.5CNT)/7075Al-1.0Mg復(fù)合材料15 h時(shí)效拉伸曲線Fig.7 Tensile curve of SiCp(0.5CNT)/7075Al-1.0Mg composite at 15 h aging

3 結(jié)束語(yǔ)

本工作采用真空熱壓燒結(jié)工藝制備了致密的SiCp(CNT)/7075Al復(fù)合材料,并通過(guò)機(jī)器學(xué)習(xí)分析方法高效完成了鋁基復(fù)合材料熱處理工藝的優(yōu)化,實(shí)現(xiàn)了SiCp(CNT)/7075Al復(fù)合材料的強(qiáng)塑性匹配制備.在熱處理優(yōu)化過(guò)程中,結(jié)合8種常見(jiàn)的機(jī)器學(xué)習(xí)算法,對(duì)SiCp(CNT)/7075Al鋁基復(fù)合材料的抗拉強(qiáng)度和延伸率進(jìn)行預(yù)測(cè),用Shapley模型和樹(shù)形模型兩種手段分析了材料特征的重要程度,逆向設(shè)計(jì)出了鋁基復(fù)合材料,并進(jìn)行了實(shí)驗(yàn)驗(yàn)證.結(jié)果表明:在470°C固溶40 min,120°C時(shí)效15 h的熱處理工藝下,SiCp(0.5CNT)/7075Al-1.0Mg復(fù)合材料抗拉強(qiáng)度和延伸率的預(yù)測(cè)值為617.48 MPa和2.98%,實(shí)驗(yàn)值為647.0 MPa和3.31%;2個(gè)物理性能的MAPE誤差依次為4.56%和9.97%,可以滿足實(shí)際應(yīng)用需求.綜上可知,本工作所采用的數(shù)據(jù)驅(qū)動(dòng)方法,對(duì)設(shè)計(jì)高性能的鋁基復(fù)合材料有一定指導(dǎo)意義.