基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型

王 杰 鄭楊艷 凌 祥

（南京工業(yè)大學(xué)a.機械與動力工程學(xué)院；b.能源科學(xué)與工程學(xué)院）

在航空航天、石油化工等行業(yè)中，高溫服役的金屬構(gòu)件材料性能隨著服役時間逐漸劣化，壽命縮短，提前或者突發(fā)性破壞極易造成重大經(jīng)濟(jì)損失，產(chǎn)生災(zāi)難性的后果。對高溫服役的金屬構(gòu)件進(jìn)行蠕變損傷分析、壽命評估是保證構(gòu)件安全運行的必要手段。為避免對在役構(gòu)件造成破壞，小沖孔試驗技術(shù)作為一種微型測試技術(shù)，在蠕變壽命預(yù)測方面與傳統(tǒng)蠕變試驗有相似的3個階段，對服役設(shè)備壽命預(yù)測已經(jīng)在石油、化工等領(lǐng)域被逐步應(yīng)用［1］。

20世紀(jì)90年代，PARKER J進(jìn)行了小沖孔蠕變試驗，提出小沖孔蠕變試驗技術(shù)可以用于對材料高溫蠕變特性的建立，表明了小沖孔試驗技術(shù)在高溫領(lǐng)域的可行性［2］。目前，小沖孔蠕變壽命的預(yù)測主要通過短期溫度或載荷函數(shù)進(jìn)行外推，預(yù)測材料蠕變的壽命。DOBE? F和MILICˇKA K基于Monkman-Grant模型使用最小蠕變速率和斷裂時間建立了小沖孔蠕變壽命預(yù)測公式［3］。文獻(xiàn)［4～8］的研究者基于Larson-Miller參數(shù)法，成功估算了不同材料的蠕變壽命。YANG S S等基于小沖孔能量模型推斷不同材料在恒溫下的蠕變壽命［9］。ZHENG Y Y等以20%蠕變壽命作為小沖孔蠕變壽命預(yù)測的起點，使用修正Theta法成功預(yù)測了已服役材料的壽命［10］。JEFFS S P等使用Wilshire方程根據(jù)短期蠕變數(shù)據(jù)預(yù)測長期蠕變壽命［11］，而后，HOLMSTROM B等提出一種可直接用于小沖孔蠕變測試的修正Wilshire方程［12］。近年來數(shù)據(jù)驅(qū)動的機器學(xué)習(xí)在蠕變壽命預(yù)測方面發(fā)展迅速。VENKATESH V和RACK H J采用反向神經(jīng)網(wǎng)絡(luò)并得到了很好的預(yù)測精度［13］。LIU Y等將機器學(xué)習(xí)應(yīng)用到鎳基單晶高溫合金的蠕變斷裂壽命預(yù)測［14］。

筆者建立一種基于粒子群算法（PSO）優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型，克服工程中難以獲得解析解和無法通過經(jīng)驗解描述高維變量之間非線性交互的局限性，實現(xiàn)對小沖孔蠕變壽命的預(yù)測，節(jié)約試驗成本，并與傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型進(jìn)行對比分析。

1 PSO-BP神經(jīng)網(wǎng)絡(luò)模型

1.1 BP神經(jīng)網(wǎng)絡(luò)

20世紀(jì)80年代，Rumelhart和McClelland提出了利用誤差反向傳播訓(xùn)練算法的神經(jīng)網(wǎng)絡(luò)，即BP神經(jīng)網(wǎng)絡(luò)（Back-Propagation Network），目前已成為應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一［15］。該模型使用梯度下降的搜索方法，利用反向傳播優(yōu)化神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值，最終使得輸出值和期望輸出值的均方根最小，具有很好的容錯性泛化能力［16，17］。BP神經(jīng)網(wǎng)絡(luò)模型具有3層拓?fù)浣Y(jié)構(gòu)，包括輸入層（Input Layer）、隱含層（Hidden Layer）和輸出層（Output Layer）。大量研究表明含有一個隱含層的神經(jīng)網(wǎng)絡(luò)可以實現(xiàn)對任何函數(shù)的逼近［15］，為避免模型過于復(fù)雜而引起計算量過大或預(yù)測結(jié)果不可靠的問題，構(gòu)造單隱含層的BP神經(jīng)網(wǎng)絡(luò)，其結(jié)構(gòu)如圖1所示。

圖1 單隱層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖

如圖1所示，輸入層和輸出層數(shù)據(jù)集維度分別為n、m，對應(yīng)為實際問題中自變量和因變量個數(shù)，隱含層神經(jīng)元個數(shù)即節(jié)點數(shù)為p。各層之間的連接權(quán)重分別為ωij（i=1，2，…，n；j=1，2，…，p）和νij（i=1，2，…，p；j=1，2，…，m）。

1.2 PSO算法

粒子群優(yōu)化算法（Particle Swarm Opitimization）是由Eberhart和Kennedy受到鳥類群體覓食過程的啟發(fā)，通過群體中個體之間的競爭和協(xié)同尋找最優(yōu)值所創(chuàng)立的一種群體自適應(yīng)搜索算法［18］。

種群由粒子構(gòu)成，位置和速度是每個粒子所包含的基本信息。設(shè)種群規(guī)模為n，目標(biāo)搜索的空間維數(shù)為D，則第i個粒子在D維空間中位置表示為Xi=（xi1，xi2，…，xiD），速度表示為Vi=（vi1，vi2，…，viD），其中i=1，2，…，n。粒子會在迭代中修正自身的位置和速度，直至找到全局范圍內(nèi)的最優(yōu)位置，其中單個粒子最優(yōu)位置為Pibest，粒子群最優(yōu)位置為Pgbest。粒子位置和速度的更新依據(jù)為：

1.3 PSO優(yōu)化BP算法

BP神經(jīng)網(wǎng)絡(luò)在預(yù)測時容易產(chǎn)生收斂速度慢、易陷入局部極值及學(xué)習(xí)過程中易發(fā)生振蕩等問題，采用PSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)，根據(jù)粒子位置和速度修正BP神經(jīng)網(wǎng)絡(luò)中的閾值和權(quán)值，以期提高BP神經(jīng)網(wǎng)絡(luò)的運行速度和泛化能力。筆者將粒子群優(yōu)化算法（PSO）與BP神經(jīng)網(wǎng)絡(luò)結(jié)合，其流程如圖2所示。

圖2 PSO-BP流程圖

2 基于PSO-BP的小沖孔蠕變壽命預(yù)測模型

2.1 模型建立

將PSO優(yōu)化BP神經(jīng)網(wǎng)絡(luò)應(yīng)用于小沖孔蠕變壽命預(yù)測中，建立如圖3所示的小沖孔蠕變壽命的預(yù)測模型，模型共分為5個模塊：試驗數(shù)據(jù)采集模塊、特征數(shù)據(jù)提取模塊、數(shù)據(jù)處理模塊、基于PSO-BP的小沖孔蠕變壽命預(yù)測模塊和數(shù)據(jù)接收模塊。

圖3 小沖孔蠕變壽命的預(yù)測模型結(jié)構(gòu)

首先用試驗數(shù)據(jù)采集模塊對小沖孔試驗進(jìn)行實時數(shù)據(jù)采集，其次用特征數(shù)據(jù)提取模塊提取小沖孔試驗中的溫度、載荷和對應(yīng)條件下的破斷時間，然后數(shù)據(jù)處理模塊將所提取的特征數(shù)據(jù)采用MATLAB提供的mapminmax函數(shù)進(jìn)行歸一化處理，接著基于PSO-BP的小沖孔蠕變壽命預(yù)測模塊將歸一化處理后的部分?jǐn)?shù)據(jù)進(jìn)行網(wǎng)絡(luò)訓(xùn)練并將剩余數(shù)據(jù)用于驗證測試，最后由數(shù)據(jù)接收模塊接收壽命預(yù)測模塊所輸出的數(shù)據(jù)。

2.2 試驗數(shù)據(jù)

針對小沖孔蠕變破斷時間的預(yù)測問題，選取JEFFS S P和LANCASTER R J對第2代鎳基單晶高溫合金CMSX-4在不同溫度和載荷下所做的小沖孔蠕變試驗數(shù)據(jù)進(jìn)行分析［19］。CMSX-4擁有較高的高溫強度、優(yōu)異的蠕變與疲勞抗力、良好的抗氧化性、抗熱腐蝕性、組織穩(wěn)定性和使用可靠性，已經(jīng)成為先進(jìn)航空發(fā)動機高壓渦輪葉片和地面燃?xì)廨啓C葉片的關(guān)鍵材料［20］。

CMSX-4的合金成分及其含量見表1。小沖孔蠕變試驗分別在950、1 050、1 150 ℃3種溫度、50～400 N載荷范圍下進(jìn)行，共獲得17組蠕變試驗數(shù)據(jù)（圖4）。

表1 CMSX-4的合金成分及其含量wt%

圖4 不同載荷、溫度下小沖孔蠕變破斷時間

2.3 參數(shù)設(shè)置

為了保證壽命預(yù)測模型的可靠性，隨機選取上述試驗中14個不同溫度、載荷條件下的小沖孔破斷時間作為模型訓(xùn)練樣本，剩余3個作為模型測試樣本，需要注意的是，3個測試樣本應(yīng)盡可能選擇不同溫度，并進(jìn)行多次重復(fù)試驗。

本模型構(gòu)造單隱含層網(wǎng)絡(luò)，隱含層的神經(jīng)元數(shù)為通過尋優(yōu)獲得的最小均方誤差所對應(yīng)的最佳節(jié)點數(shù)。訓(xùn)練迭代的最大次數(shù)為1 000，學(xué)習(xí)速率為0.01，期望誤差為0.001。PSO算法的粒子群規(guī)模n=10，慣性權(quán)重ω=0.8，進(jìn)化次數(shù)為30，學(xué)習(xí)因子c1=c2=1.99445。

3 結(jié)果及對比分析

采用MATLAB編程模擬神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過程，得到如圖5所示的算法學(xué)習(xí)性能。

圖5 算法學(xué)習(xí)性能

圖5a所示的BP神經(jīng)網(wǎng)絡(luò)的算法性能中，前7步3條誤差線幾乎重合，7步之后逐漸分開，迭代次數(shù)為105時得到的訓(xùn)練誤差最小，為3.6746×10-5；圖5b所示的PSO-BP神經(jīng)網(wǎng)絡(luò)的算法性能中，前5步3條誤差線幾乎重合，5步之后逐漸分開，迭代次數(shù)為11時得到的訓(xùn)練誤差最小，為1.2779×10-5；兩種預(yù)測模型對比可知，基于PSO優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)收斂速度塊，均方誤差小且整體性能好。

為了檢驗基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型的準(zhǔn)確性和穩(wěn)定性，采用多次隨機抽取測試樣本的方法對模型的預(yù)測結(jié)果進(jìn)行誤差分析，預(yù)測結(jié)果見表2，預(yù)測結(jié)果對比如圖6所示。

表2 預(yù)測結(jié)果與誤差分析

圖6 破斷時間預(yù)測結(jié)果對比

由圖6a可知，PSO-BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型的預(yù)測值相較于BP神經(jīng)網(wǎng)絡(luò)更接近試驗值，預(yù)測結(jié)果更加準(zhǔn)確。結(jié)合表2、圖6b可知，PSO-BP預(yù)測模型最大相對誤差為0.329 3，平均相對誤差為0.102 2；傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)的最大相對誤差為0.399 9，平均相對誤差為0.182 5。考慮到蠕變試驗因為實際試驗條件的不確定性，有學(xué)者指出蠕變試驗的破斷時間偏差在20%以內(nèi)均可接受［21］。在第6次測試中，BP神經(jīng)網(wǎng)絡(luò)預(yù)測值為247.27 h，相對誤差為0.186 7；PSO-BP神經(jīng)網(wǎng)絡(luò)預(yù)測值為276.97 h，相對誤差為0.329 3，與實際值208.36 h誤差更為顯著，可能是蠕變試驗的誤差，也可能是模型自身仍存在一些問題。但整體情況而言，基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型預(yù)測結(jié)果符合實際應(yīng)用情況，且相較于傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)泛化能力更勝一籌。

4 結(jié)論

4.1 采用粒子群優(yōu)化算法（PSO）優(yōu)化了BP神經(jīng)網(wǎng)絡(luò)的權(quán)值、閾值，根據(jù)不同溫度和載荷下小沖孔蠕變破斷時間，建立基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型，解決了BP神經(jīng)網(wǎng)絡(luò)收斂速度慢、易陷入局部極值的問題，提高了小沖孔蠕變壽命的預(yù)測精度。

4.2 用PSO-BP神經(jīng)網(wǎng)絡(luò)小沖孔蠕變壽命預(yù)測模型多次隨機選取測試樣本進(jìn)行預(yù)測，破斷時間預(yù)測結(jié)果與試驗值的平均相對誤差為0.102 2，BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型的平均相對誤差為0.182 5，基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型預(yù)測精度較BP神經(jīng)網(wǎng)絡(luò)更高，穩(wěn)定性更好，且預(yù)測相對誤差除第6次試驗，其余均在蠕變試驗破斷時間偏差范圍內(nèi)，達(dá)到很好的預(yù)期，符合實際要求。

4.3 筆者建立的基于PSO-BP神經(jīng)網(wǎng)絡(luò)的小沖孔蠕變壽命預(yù)測模型預(yù)測精度整體較好，但其慣性權(quán)重、學(xué)習(xí)因子等參數(shù)還需要依據(jù)經(jīng)驗進(jìn)行選取，缺乏相應(yīng)理論的指導(dǎo)，在第6次測試中誤差顯著，如何選取預(yù)測模型參數(shù)提高預(yù)測的準(zhǔn)確性和穩(wěn)定性仍需要進(jìn)一步研究。