王 正，劉建雄，王 璐，謝偉云

(大連理工大學(xué)能源與動(dòng)力學(xué)院，大連116024)

改進(jìn)神經(jīng)網(wǎng)絡(luò)在疲勞短裂紋演化行為中的應(yīng)用

王 正，劉建雄，王 璐，謝偉云

(大連理工大學(xué)能源與動(dòng)力學(xué)院，大連116024)

為研究短裂紋演化行為中復(fù)雜的非線性動(dòng)力學(xué)過(guò)程，采用改進(jìn)BP神經(jīng)網(wǎng)絡(luò)算法對(duì)疲勞短裂紋的演化行為進(jìn)行表征.該方法采用遺傳算法優(yōu)化確定神經(jīng)網(wǎng)絡(luò)的權(quán)重，同時(shí)集合BP網(wǎng)絡(luò)算法的局部精確搜索和遺傳算法的宏觀搜索、全局優(yōu)化特性，可以綜合多個(gè)影響因素，反映其隱含的復(fù)雜非線性關(guān)系.通過(guò)對(duì)復(fù)雜應(yīng)力狀態(tài)下高溫低周疲勞短裂紋的試驗(yàn)研究及疲勞短裂紋密度和裂紋擴(kuò)展速率的模擬比較，表明該方法收斂速度更快、計(jì)算更精確，基于該方法建立的疲勞短裂紋演化模型合理有效.

高溫低周疲勞;短裂紋;改進(jìn)神經(jīng)網(wǎng)絡(luò);裂紋密度;裂紋擴(kuò)展速率

近年來(lái)，對(duì)疲勞短裂紋的研究已成為整個(gè)疲勞斷裂研究領(lǐng)域的重點(diǎn)之一.據(jù)統(tǒng)計(jì)，有50%～90%的工程結(jié)構(gòu)或機(jī)械零件破壞與疲勞有關(guān)，而約90%的疲勞損傷壽命都消耗在裂紋的萌生、擴(kuò)展階段［1］.雖然Paris公式［2］對(duì)滿足線彈性斷裂力學(xué)(LEFM)條件的疲勞裂紋擴(kuò)展行為的描述取得了成功，但是短裂紋的擴(kuò)展受到材料宏觀結(jié)構(gòu)和微觀組織如晶粒大小、包含物形狀及方向、晶界分布、材料成分以及局部各向異性等的強(qiáng)烈影響［3］，具有LEFM所不能描述的“異?！爆F(xiàn)象，其擴(kuò)展規(guī)律與長(zhǎng)裂紋有很大不同，因此對(duì)復(fù)雜應(yīng)力狀態(tài)下高溫低周疲勞短裂紋進(jìn)行研究就顯得十分重要.對(duì)于疲勞短裂紋擴(kuò)展規(guī)律的描述，Miller［4］、Hobson［5］、Polak［6］等都做了大量的工作，提出了各自的模型.但是對(duì)于模型中的一些參數(shù)，由于沒(méi)有十分準(zhǔn)確的物理意義，因此在實(shí)際操作中難以應(yīng)用.研究表明，短裂紋的萌生與擴(kuò)展是非線性動(dòng)力學(xué)演化的結(jié)果，如滑移帶、裂紋路徑及斷口的分形特性、損傷局部化、疲勞裂紋密度和擴(kuò)展速率的分散性等［7］.

神經(jīng)網(wǎng)絡(luò)為解決復(fù)雜的非線性動(dòng)力學(xué)系統(tǒng)問(wèn)題提供了有效的算法，由于其具有較強(qiáng)的學(xué)習(xí)、記憶、聯(lián)想、識(shí)別等功能，以及實(shí)時(shí)、并行集體運(yùn)算能力，近年來(lái)被廣泛用于研究材料科學(xué)中復(fù)雜的非線性問(wèn)題，如材料設(shè)計(jì)、材料制備工藝優(yōu)化、塑性加工、熱處理以及材料使用壽命等方面［8］.本文采用改進(jìn)神經(jīng)網(wǎng)絡(luò)對(duì)疲勞短裂紋的演化行為進(jìn)行建模，通過(guò)對(duì)疲勞短裂紋密度和擴(kuò)展速率進(jìn)行分析研究，建立疲勞短裂紋的演化預(yù)測(cè)模型.

1 改進(jìn)神經(jīng)網(wǎng)絡(luò)方法

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

神經(jīng)網(wǎng)絡(luò)模型主要考慮網(wǎng)絡(luò)連接的拓?fù)浣Y(jié)構(gòu)、神經(jīng)元的特征、學(xué)習(xí)規(guī)則等.基于誤差反向傳播算法的多層前饋網(wǎng)絡(luò)(Error Back Propagation Neural Network，BPNN)是目前應(yīng)用最為廣泛最成功的網(wǎng)絡(luò)之一.在實(shí)際應(yīng)用中，80%～90%的人工神經(jīng)網(wǎng)絡(luò)模型是采用BP網(wǎng)絡(luò)［9］.

BPNN是采用Widrow-Hoff學(xué)習(xí)算法和非線性可微轉(zhuǎn)移函數(shù)的多層網(wǎng)絡(luò)，其信息處理功能是利用映射訓(xùn)練樣本(x1，y1)，(x2，y2)，…，(xk，yk)，…實(shí)現(xiàn)從n維歐氏空間的一個(gè)子集A到m維歐氏空間子集F(A)的映射，即F:AF:A?Rn→Rm.一般，BPNN是由輸入層、輸出層以及一個(gè)或多個(gè)隱層構(gòu)成，根據(jù)Kplmogorov定理，只有1個(gè)隱含層的三層BPNN可以實(shí)現(xiàn)在任意精度內(nèi)對(duì)任意非線性連續(xù)函數(shù)的逼近，如圖1所示.BPNN的具體算法過(guò)程見(jiàn)參考文獻(xiàn)［10］.

圖1 一個(gè)典型的BP網(wǎng)絡(luò)

1.2 改進(jìn)神經(jīng)網(wǎng)絡(luò)

BP算法雖然具有精確尋優(yōu)的能力，但同時(shí)存在容易陷入局部極小值、收斂速度慢和引起振蕩效應(yīng)等缺點(diǎn)，從而使訓(xùn)練精度難以提高.由于遺傳算法具有很強(qiáng)的宏觀搜索能力，且能以較大的概率找到全局最優(yōu)解，所以用來(lái)確定神經(jīng)網(wǎng)絡(luò)的權(quán)重能較好地克服BP算法的缺點(diǎn).因此，改進(jìn)神經(jīng)網(wǎng)絡(luò)算法既保存了遺傳算法“很強(qiáng)的宏觀搜索能力和良好的全局優(yōu)化能力”，又發(fā)揚(yáng)了神經(jīng)網(wǎng)絡(luò)的“較強(qiáng)的非線性逼近能力”.

遺傳算法(Genetic Algorithm，GA)的本質(zhì)是一種高效、并行、全局搜索的方法，它能在搜索過(guò)程中自動(dòng)獲取和積累有關(guān)搜索空間的知識(shí)，并自適應(yīng)的控制過(guò)程以求得最優(yōu)解.由于遺傳算法的整體搜索策略和優(yōu)化計(jì)算不依賴于梯度信息，具有很好的魯棒性，故在處理高度復(fù)雜的非線性問(wèn)題時(shí)，表現(xiàn)出了無(wú)可比擬的優(yōu)越性.

GA優(yōu)化BPNN的要素包括種群初始化、適應(yīng)度函數(shù)(fitness)、選擇操作(selection)、交叉操作(crossover)和突變操作(mutation).適應(yīng)度用來(lái)度量群體中各個(gè)體在優(yōu)化計(jì)算中可能達(dá)到、接近或有助于找到最優(yōu)解的優(yōu)良程度.適應(yīng)度函數(shù)采用:

其中，SE是神經(jīng)網(wǎng)絡(luò)的均方誤差，α是輸入層節(jié)點(diǎn)和隱含層節(jié)點(diǎn)的影響系數(shù)(0＜α≤0.01)，In-Dim、HiddenUnitNum分別是輸入層節(jié)點(diǎn)數(shù)目和隱含層節(jié)點(diǎn)數(shù)目.為避免初始權(quán)值和閾值對(duì)適應(yīng)度函數(shù)計(jì)算的影響，針對(duì)每一個(gè)體計(jì)算適應(yīng)度函數(shù)值時(shí)，均用遺傳算法對(duì)權(quán)值和閾值進(jìn)行優(yōu)化.

選擇操作有輪盤賭法、競(jìng)標(biāo)賽法等多種方法，為避免局部最優(yōu)解主導(dǎo)整個(gè)迭代過(guò)程，以及全局最優(yōu)值在形成之前就被淘汰，算法中優(yōu)先選擇適應(yīng)度值高于10%范圍內(nèi)的一代直接復(fù)制進(jìn)入下一代，其余的染色體通過(guò)輪盤賭法選擇.

對(duì)于BP神經(jīng)網(wǎng)絡(luò)初始權(quán)值和閾值的優(yōu)化，交叉操作采用算術(shù)交叉算子，其方法如下:

式中:p'1、p'2為交叉操作后的新個(gè)體;p1、p2為一組配對(duì)的兩個(gè)個(gè)體;r1、r2是［0，1］間的隨機(jī)數(shù).

為了避免算法迭代后期出現(xiàn)種群過(guò)早收斂，同時(shí)保持種群的多樣性，需要對(duì)新生代進(jìn)行突變操作.設(shè)p是父代，子染色體p'=p+ε，其中ε是［0，1］間的隨機(jī)數(shù).

傳統(tǒng)的交叉、突變操作基于一定的隨機(jī)性，但事實(shí)上自然界遺傳演化是動(dòng)態(tài)的，因此引入一種改進(jìn)的表達(dá)式:

式中:f是個(gè)體適應(yīng)度;fave是種群平均適應(yīng)度; fmax、fmin分別是種群適應(yīng)度的最大、最小值.

利用改進(jìn)神經(jīng)網(wǎng)絡(luò)建立疲勞短裂紋演化行為模型的步驟如下:

1)導(dǎo)入數(shù)據(jù)并歸一化，得到P_train、P_test、T _train.歸一化的具體作用是歸納統(tǒng)一樣本的統(tǒng)計(jì)分布性，歸一化后的數(shù)據(jù)分布在［-1，1］之間，目的在于加快網(wǎng)絡(luò)的收斂性;

2)創(chuàng)建、訓(xùn)練及仿真單BP神經(jīng)網(wǎng)絡(luò);

3)初始化參數(shù)，如popsize，Generationmax;

4)初始化種群，計(jì)算個(gè)體適應(yīng)度f(wàn);

5)GA優(yōu)化BPNN:根據(jù)適應(yīng)度值選擇和復(fù)制種群;計(jì)算交叉概率pcrossover，進(jìn)行交叉操作;計(jì)算突變概率pmutation，進(jìn)行突變操作;

6)計(jì)算種群適應(yīng)度值fitness，找到最優(yōu)個(gè)體Chrombest;

7)generation=generation+1，if generation≤Generationmax進(jìn)入下一代，返回執(zhí)行步驟(5); else停止迭代;

8)提取新訓(xùn)練集/測(cè)試集數(shù)據(jù)p_train、p_test，創(chuàng)建優(yōu)化BP神經(jīng)網(wǎng)絡(luò)，訓(xùn)練、仿真測(cè)試，得到預(yù)測(cè)輸出數(shù)據(jù)，對(duì)得到的數(shù)據(jù)進(jìn)行反歸一化，并與實(shí)驗(yàn)數(shù)據(jù)相比較.

2 疲勞短裂紋演化行為模型

2.1 疲勞短裂紋演化行為描述

材料疲勞短裂紋的演化受諸多因素的影響，如材料特性、載荷特性以及環(huán)境因素等.對(duì)于疲勞裂紋密度，可以用下式描述:

式中:N為載荷循環(huán)次數(shù);Δσ為外加應(yīng)力幅;R為應(yīng)力比;Kc為材料斷裂韌度;d為材料的微結(jié)構(gòu)尺度;C為環(huán)境影響因素.

對(duì)于疲勞短裂紋擴(kuò)展速率，可用下式表示:

式中，a為載荷循環(huán)N次時(shí)的裂紋長(zhǎng)度.

董聰?shù)龋?］和劉長(zhǎng)虹等［11］研究表明，疲勞裂紋密度和疲勞短裂紋擴(kuò)展速率與先前的演化歷史相關(guān)，是非線性動(dòng)力學(xué)系統(tǒng)的描述.其中f(a，△σ，R，Kc，d，C…)在一定條件下表現(xiàn)出非線性自相關(guān)混沌的特性.

已有研究表明，影響疲勞裂紋演化行為的因素有:應(yīng)力幅、應(yīng)力比、材料的斷裂韌性、材料微結(jié)構(gòu)尺度和環(huán)境影響因素等.疲勞短裂紋演化行為與各影響因素的關(guān)系，不是簡(jiǎn)單的線性迭加關(guān)系，而是一種十分復(fù)雜的非線性映射關(guān)系.但是通常的疲勞短裂紋演化行為研究，一般是改變單一變量，而保持其他影響因素不變，從中找出疲勞短裂紋演化行為與該變量的關(guān)系，然后采用類似的方法再研究疲勞短裂紋演化行為與其他因素的關(guān)系.這種方法雖然能反映出疲勞短裂紋演化行為隨某一因素變化的趨勢(shì)，但不能反映出實(shí)際的非線性映射關(guān)系.

BPNN就是利用已有的映射關(guān)系作為學(xué)習(xí)樣本，根據(jù)學(xué)習(xí)法則，尋求最佳權(quán)集以實(shí)現(xiàn)正確輸出.這個(gè)權(quán)集就代表了輸出的近似解析關(guān)系與真正的映射關(guān)系之間的誤差.

2.2 高溫低周疲勞短裂紋演化試驗(yàn)

應(yīng)用改進(jìn)神經(jīng)網(wǎng)絡(luò)，首先需根據(jù)問(wèn)題確定網(wǎng)絡(luò)輸入層、隱含層及輸出層，然后用疲勞試驗(yàn)數(shù)據(jù)加以訓(xùn)練、預(yù)測(cè).試驗(yàn)所用材料為20鋼(質(zhì)量分?jǐn)?shù)0.20%C)，材料力學(xué)性能參數(shù)如表1所示.試樣采用帶圓環(huán)缺口的光滑圓棒，其尺寸如圖2所示.由于不同尺寸缺口的存在，試樣在簡(jiǎn)單軸向加載情況下在缺口處產(chǎn)生不同應(yīng)力分量比的三維應(yīng)力場(chǎng).缺口表面經(jīng)金相砂紙和拋光機(jī)重點(diǎn)拋光，以避免裂紋在表面缺陷處優(yōu)先產(chǎn)生.

表1 20鋼的力學(xué)性能參數(shù)

圖2 試樣尺寸

試驗(yàn)采用MTS Landmark 100KN材料疲勞試驗(yàn)機(jī)，應(yīng)變控制，循環(huán)特性為R=-1的三角波，加載頻率為0.5 Hz，試驗(yàn)溫度400和500℃.采用中斷試驗(yàn)，對(duì)試件缺口中心表面處裂紋進(jìn)行圖像采集.

試驗(yàn)過(guò)程中發(fā)現(xiàn)，疲勞短裂紋一般是在試樣缺口中心附近或試樣表面缺陷處首先出現(xiàn)，在擴(kuò)展過(guò)程中多呈Z形分布.為研究疲勞裂紋擴(kuò)展的方向路徑，以加載方向?yàn)榛鶞?zhǔn)，忽略裂紋本身彎曲而以基本走向?yàn)闇?zhǔn)對(duì)裂紋角度進(jìn)行測(cè)量.圖3顯示了不同缺口尺寸和不同應(yīng)變幅條件下裂紋平均角度的變化情況，其中Nf為循環(huán)應(yīng)力范圍下降到初始應(yīng)力范圍的75%的循環(huán)次數(shù).

圖3 裂紋平均角度圖

由圖3可以看出，疲勞短裂紋角度數(shù)據(jù)比較分散，裂紋平均角度大部分處于75°以上，壽命分?jǐn)?shù)對(duì)裂紋平均角度影響不大.試樣在不同壽命分?jǐn)?shù)下，裂紋平均角度在裂紋擴(kuò)展過(guò)程中是波動(dòng)變化的.對(duì)于不同缺口形狀和不同載荷條件，短裂紋平均角度均大于70°，這說(shuō)明了短裂紋分布具有較強(qiáng)的角度依存性.王璐等［12］研究發(fā)現(xiàn)，高溫低周條件下大多裂紋沿與疲勞表面的應(yīng)力軸接近垂直的方向優(yōu)先萌生與擴(kuò)展.試驗(yàn)統(tǒng)計(jì)結(jié)果與其基本一致，這也進(jìn)一步說(shuō)明了短裂紋的萌生與擴(kuò)展是有方向選擇性的.

3 結(jié)果分析

3.1 疲勞短裂紋密度

疲勞短裂紋在萌生和擴(kuò)展過(guò)程中，同時(shí)還存在裂紋的合體，在整個(gè)裂紋演化過(guò)程中，裂紋密度的變化體現(xiàn)了短裂紋的萌生和合體狀況.從試驗(yàn)統(tǒng)計(jì)結(jié)果發(fā)現(xiàn)，短裂紋密度并不是線性變化的，而是呈現(xiàn)先快速增長(zhǎng)，然后逐漸減小再趨于平穩(wěn)的趨勢(shì).這主要是因?yàn)樵趬勖捌诩s30%，即微觀結(jié)構(gòu)短裂紋(MSC)階段，受材料微觀組織影響，材料損傷主要以裂紋萌生為主.而后期由于主裂紋的形成，損傷主要以裂紋擴(kuò)展為主，伴隨發(fā)生裂紋合體等現(xiàn)象，這也預(yù)示著裂紋群體性逐漸減弱，逐漸演化為單條疲勞裂紋行為.當(dāng)循環(huán)次數(shù)達(dá)到一定程度后，出現(xiàn)了明顯的裂紋密度飽和現(xiàn)象.

在疲勞裂紋演化行為模型中，對(duì)裂紋密度的神經(jīng)網(wǎng)絡(luò)訓(xùn)練，選取外加應(yīng)力幅、缺口半徑、應(yīng)變幅、溫度、壽命分?jǐn)?shù)作為輸入節(jié)點(diǎn)變量，以裂紋密度n(N)作為輸出層節(jié)點(diǎn)變量，即網(wǎng)絡(luò)輸入層的神經(jīng)元數(shù)(InDim)取5，輸出層神經(jīng)元(OutDim)取1.由于隱含層的神經(jīng)元個(gè)數(shù)(HiddenUnitNum)選取沒(méi)有明確的公式，為防止過(guò)擬合［13］，在遵循保證神經(jīng)網(wǎng)絡(luò)收斂的情況下個(gè)數(shù)盡量少的原則，通過(guò)逐漸增大數(shù)值選取最佳值.根據(jù)多次實(shí)驗(yàn)確定隱含層神經(jīng)元HiddenUnitNum=17較為合適.

模型訓(xùn)練時(shí)，學(xué)習(xí)率η=0.05，動(dòng)量因子α= 0.95，目標(biāo)誤差為ε=1×10-4.遺傳算法的參數(shù)選取:種群規(guī)模50，遺傳代數(shù)100.訓(xùn)練中采用8個(gè)樣本對(duì)，共112組數(shù)據(jù)，隨機(jī)選取98組數(shù)據(jù)為學(xué)習(xí)樣本訓(xùn)練神經(jīng)網(wǎng)絡(luò)，其余14組數(shù)據(jù)作為檢測(cè)樣本驗(yàn)證模擬網(wǎng)絡(luò)性能.

為比較模型的性能，也采用單BP網(wǎng)絡(luò)對(duì)模型進(jìn)行了訓(xùn)練.單BP神經(jīng)網(wǎng)絡(luò)模型采用與改進(jìn)神經(jīng)網(wǎng)絡(luò)相同的參數(shù)，以及相同的訓(xùn)練、測(cè)試數(shù)據(jù).兩種神經(jīng)網(wǎng)絡(luò)模型的裂紋密度對(duì)比見(jiàn)圖4.

圖4 神經(jīng)網(wǎng)絡(luò)模型的裂紋密度對(duì)比圖

在單BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練中，出現(xiàn)了收斂時(shí)間過(guò)長(zhǎng)、收斂速度過(guò)慢的情況，而在引入改進(jìn)算法以后，非常明顯地改進(jìn)了收斂速度過(guò)慢的問(wèn)題.由圖4可以看出，單BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)相對(duì)誤差較大，而在采用改進(jìn)BP神經(jīng)網(wǎng)絡(luò)后，預(yù)測(cè)值與試驗(yàn)值吻合較好，基本能反映出疲勞短裂紋密度與各主要影響因素間的非線性映射關(guān)系，且具有較高的精度.

3.2 疲勞短裂紋擴(kuò)展速率

對(duì)于裂紋擴(kuò)展速率的神經(jīng)網(wǎng)絡(luò)訓(xùn)練，選取裂紋長(zhǎng)度、外加應(yīng)力幅、缺口半徑、應(yīng)變幅、溫度、壽命分?jǐn)?shù)作為輸入節(jié)點(diǎn)變量，以裂紋擴(kuò)展速率da/ dN作為輸出層節(jié)點(diǎn)變量.輸入層神經(jīng)元InDim= 6，輸出層神經(jīng)元OutDim=1，根據(jù)多次實(shí)驗(yàn)確定隱含層神經(jīng)元HiddenUnitNum=12較為合適.

BP神經(jīng)網(wǎng)絡(luò)采用與短裂紋密度模型相同的三層結(jié)構(gòu)，其余參數(shù)設(shè)置與裂紋密度模型相同.采用168組數(shù)據(jù)，選取其中的137組作為學(xué)習(xí)樣本訓(xùn)練神經(jīng)網(wǎng)絡(luò)，其余14組數(shù)據(jù)作為檢測(cè)樣本，驗(yàn)證網(wǎng)絡(luò)性能.訓(xùn)練結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比見(jiàn)圖5.

圖5 神經(jīng)網(wǎng)絡(luò)模型的裂紋擴(kuò)展速率對(duì)比圖

由圖5可知，使用改進(jìn)BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)數(shù)據(jù)與試驗(yàn)數(shù)據(jù)吻合較好.疲勞短裂紋擴(kuò)展速率的演化過(guò)程呈現(xiàn)出由慢變快又減速的變化趨勢(shì)，表明了短裂紋擴(kuò)展速率具有自身特有的規(guī)律.出現(xiàn)如圖5所示的變化趨勢(shì)，是因?yàn)榱鸭y萌生后，在3～5個(gè)晶粒尺度［14］時(shí)期內(nèi)，其生長(zhǎng)行為受材料微觀組織的強(qiáng)烈影響，其中以晶界對(duì)疲勞短裂紋生長(zhǎng)的阻礙作用表現(xiàn)最為顯著.當(dāng)裂紋擴(kuò)展至晶界時(shí)，裂紋尖端應(yīng)力集中在一定程度上得到釋放，使擴(kuò)展受阻而減慢，因此其擴(kuò)展速率平緩;隨著循環(huán)的進(jìn)行，當(dāng)裂紋擴(kuò)展到脫離微觀組織影響時(shí)，裂紋擴(kuò)展速率又加快.

通過(guò)上述模型實(shí)驗(yàn)分析，該神經(jīng)網(wǎng)絡(luò)具有較強(qiáng)的通用性，綜合考慮了影響疲勞短裂紋演化行為的諸多因素，只要用于學(xué)習(xí)的樣本含有全面的信息，網(wǎng)絡(luò)初始結(jié)構(gòu)、特征參數(shù)選取合理恰當(dāng)，就可獲得良好的預(yù)測(cè)結(jié)果，從而減少試驗(yàn)數(shù)量，節(jié)約試驗(yàn)成本.

4 結(jié)論

1)通過(guò)高溫低周疲勞短裂紋試驗(yàn)，結(jié)果表明疲勞短裂紋角度分散，具有方向選擇性，其與加載方向角度大于75°.

2)采用改進(jìn)BP神經(jīng)網(wǎng)絡(luò)，對(duì)疲勞短裂紋演化行為建立了預(yù)測(cè)模型，模擬了疲勞短裂紋密度和擴(kuò)展速率，訓(xùn)練結(jié)果表明，采用此模型可以在綜合各因素試驗(yàn)條件下獲得比較精確的結(jié)果，該方法是處理疲勞短裂紋數(shù)據(jù)和描述疲勞裂紋演化行為較好的工具，有利于推廣.

［1］ SUBRA S.Fatigue of Materials［M］(2nd ed).London:Cambridge University Press，1998.

［2］ PARIS P C，ERDOGAN F.A critical analysis of crack propagation laws［J］.Trans ASME Ser D，1963，85: 528-534.

［3］ HOBSON P D.The growth of short fatigue cracks in a medium carbon steel［D］.Sheffield:University of Sheffield，1985:468-492.

［4］ MILLER K J.Microstructural short cracks in fatigue［C］//H C VAN ELST，A.BAKKER.Proceedings of fracture control of engineering structures，Amsterdam: Engineering Materials Advisory Services Ltd，1986: 2149-2167.

［5］ MILLER K J，DE LOS RIOS E R.The Behavior of Short Fatigue Cracks［M］.London:Mechanical Engineering Publications，1986:441-459.

［6］ OBRTLIK K，POLAK J，HAJEK M，et al.Short fatigue crack behavior in 316L stainless steel［J］.Int J Fatigue，1997，19:471-475.

［7］ 董聰，何慶芝.微裂紋演化過(guò)程中分岔與混沌現(xiàn)象的描述及若干問(wèn)題探討［J］.力學(xué)進(jìn)展，1994，24 (1):106-116.DONG Cong，HE Qing-zhi.Description and discussion on bifurcation and chaos in evolution of microcrack［J］.Advances in Mechanics，1994，24(1):106-116.

［8］ 徐強(qiáng)，張幸紅，韓杰才，等.人工神經(jīng)網(wǎng)絡(luò)在材料科學(xué)中的應(yīng)用與展望［J］.材料科學(xué)與工藝，2005，13 (4):352-356.XU Qiang，ZHANG Xing-hong，HAN Jie-cai，et al.Application and prospect of artificial neural networks in materials science［J］.Materials Science＆Technology，2005，13(4):352-356.

［9］ 陳新，黃洪鐘，孫道恒.應(yīng)用神經(jīng)網(wǎng)絡(luò)的疲勞裂紋演化規(guī)律的描述［J］.鐵道學(xué)報(bào)，2000，22(3):33-37.CHEN Xin，HUANG Hong-zhong，SUN Dao-heng.Application of neural networks to the evolution of fatigue cracks［J］.Journal of The China Railway Society，2000，22(3):33-37.

［10］YU Shi-wei，GUO Xiu-fu，ZHU Ke-jun，et al.A neurofuzzy GA-BP method of seismic reservoir fuzzy rules extraction［J］.Expert Systems with Applications，2010，37:2037-2042.

［11］ 劉長(zhǎng)虹，陳虬，丁宇，等.疲勞短裂紋擴(kuò)展中的混沌現(xiàn)象［J］.大型鑄鍛件，2001(3):10-14.LIU Chang-hong，CHEN Qiu，DING Yu，et al.The chaotic behaviors in the development of fatigue short crack［J］.Heavy Castings and Forgings，2001(3):10-14.

［12］ 王璐，王正，于淼.高溫低周疲勞表面短裂紋合體與干涉行為的實(shí)驗(yàn)研究及數(shù)值模擬［J］.機(jī)械強(qiáng)度，2008，138(04):642-646.WANG Lu，WANG Zheng，YU Miao.Experimental study and numerical simulation on coalescence and interference of short cracks for low cycle fatigue at high temperature［J］.Journal of Mechanical Strength，2008，138(04):642-646.

［13］ KAASTRA I，BOYD M.Designing a neural network for forecasting financial and economic time series［J］.Neurocomputing，1996:10(3):215-236.

［14］ 呂文閣.疲勞短裂紋形成和擴(kuò)展的劃分定義探討［J］.機(jī)械工程材料，1999，23(6):18-19.Lü Wen-ge.Study on the definition of division of short fatigue crack initiation and propagation［J］.Materials for Mechanical Engineering，1999，23(6): 18-19.

Application of improved back-propagation neural network to short fatigue crack evolution

WANG Zheng，LIU Jian-xiong，WANG Lu，XIE Wei-yun
(School of Energy and Power Engineering，Dalian University of Technology，Dalian 116024，China)

To research the complicated nonlinear dynamics process of the short crack evolution behavior，a way that improves back-propagation neural network aiming at evolution of short fatigue crack is shown in this paper.This method optimizes the weight of the BP network，and aggregates the characteristics of the local precise search of the BP network and the global optimization of the improved genetic algorithm，which integrates more factors and reflects complicated relation.Comparing the results of the experiment of short fatigue crack for low cycle under complex stress at high temperature with the simulation results of improved back-propagation neural network，it is proved that the method is feasible，accurate and converged quickly.

high temperature low cycle fatigue;short crack;improved back-propagation neural network;crack density;crack propagation rate

TB303 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1005-0299(2012)06-0045-05

2011-09-13

國(guó)家自然科學(xué)基金資助項(xiàng)目(50771024).

王 正(1960-)，男，教授.

王 正，E-mail:zhengw@dlut.edu.cn.

(編輯 程利冬)