基于混沌分形理論的金屬疲勞損傷過程的特征分析

張玉華, 李欣欣, 黃振峰, 毛漢領(lǐng)

(1.廣西大學(xué) 輕工與食品工程學(xué)院,南寧 530004； 2.廣西大學(xué) 機(jī)械工程學(xué)院,南寧 530004)

基于混沌分形理論的金屬疲勞損傷過程的特征分析

張玉華1, 李欣欣2, 黃振峰2, 毛漢領(lǐng)2

(1.廣西大學(xué) 輕工與食品工程學(xué)院,南寧 530004； 2.廣西大學(xué) 機(jī)械工程學(xué)院,南寧 530004)

為了有效地分析金屬疲勞損傷過程中裂紋擴(kuò)展與演化的特征，提出利用混沌分形理論對(duì)蘊(yùn)含疲勞裂紋擴(kuò)展規(guī)律的超聲非線性輸出信號(hào)進(jìn)行分析，用Lyapunov指數(shù)、K熵、關(guān)聯(lián)維數(shù)等特征量對(duì)金屬疲勞損傷過程的特性進(jìn)行表征和評(píng)價(jià)。采用RITEC公司的RAM-5000系統(tǒng)獲得具有不同疲勞損傷程度試件的超聲非線性輸出信號(hào)，計(jì)算疲勞過程中的Lyapunov指數(shù)、K熵、關(guān)聯(lián)維數(shù)。分析結(jié)果表明金屬疲勞損傷過程中超聲非線性輸出信號(hào)具有混沌特性；Lyapunov指數(shù)、關(guān)聯(lián)維數(shù)、K熵等特征值隨疲勞循環(huán)次數(shù)單調(diào)遞增，當(dāng)宏觀裂紋出現(xiàn)后趨于飽和，說明疲勞損傷過程中試件內(nèi)部的復(fù)雜程度和混沌特征逐漸增強(qiáng)；由疲勞過程中微裂紋的演化規(guī)律可知，混沌分形特征值可以有效地表征試件在疲勞過程中裂紋的變化規(guī)律和群體性行為。因此混沌分形理論可以有效地分析金屬疲勞損傷過程中裂紋擴(kuò)展與演化的特征，根據(jù)混沌分形特征值—疲勞循環(huán)次數(shù)的關(guān)系，為試件的疲勞壽命預(yù)測(cè)提供了一種新的分析方法。

非線性超聲；疲勞損傷；Lyapunov指數(shù)；關(guān)聯(lián)維數(shù)；K熵

金屬試件在服役過程中，疲勞斷裂是主要的失效形式。在外界載荷作用下，試件內(nèi)部微觀結(jié)構(gòu)表現(xiàn)為位錯(cuò)群的產(chǎn)生、駐留滑移帶和微裂紋的形成和擴(kuò)展，宏觀裂紋的產(chǎn)生，直至試件斷裂失效[1]。(在早期力學(xué)性能退化階段，由于位錯(cuò)和滑移帶等微觀缺陷，單一頻率的超聲波在金屬內(nèi)傳播時(shí)波形發(fā)生畸變，產(chǎn)生高次諧波，利用二次諧波幅值對(duì)材料的早期疲勞損傷進(jìn)行檢測(cè)和評(píng)價(jià)。大量學(xué)者開展了許多相關(guān)研究，Shah等[2]通過實(shí)驗(yàn)探究線性和非線性參數(shù)對(duì)材料疲勞的敏感程度，隨著疲勞次數(shù)的增加，相對(duì)非線性系數(shù)明顯增加，但衰減系數(shù)、聲速等線性參數(shù)變化不大。吳斌、稅國(guó)雙等[3-4]探究了聲學(xué)非線性參數(shù)和材料疲勞性能的關(guān)系。Matlack等[5]利用瑞利波探究汽車零件17-4PH不銹鋼的熱輻射損傷，隨著試件老化時(shí)間的增加，試件的聲學(xué)非線性參數(shù)逐漸減小，該結(jié)果與試件的微觀結(jié)構(gòu)演化規(guī)律相一致。陳振華等[6]利用高次諧波檢測(cè)薄板點(diǎn)焊焊核的微小缺陷，表明非線性超聲檢測(cè)具有很高的檢測(cè)能力。但是這些研究通過對(duì)疲勞過程中的超聲非線性輸出信號(hào)進(jìn)行頻譜分析，由二次諧波和基波幅值探究試件的非線性特性。由頻譜分析的原理可知，該方法適合處理平穩(wěn)的線性信號(hào)[7]。而在疲勞損傷過程中，微裂紋個(gè)數(shù)、尺寸、擴(kuò)展路徑等一直在發(fā)生變化[8]，其微觀結(jié)構(gòu)高度不均勻且十分復(fù)雜，超聲非線性輸出信號(hào)是反映金屬疲勞損傷、并經(jīng)過復(fù)雜傳播的非平穩(wěn)時(shí)間序列，如何有效地提取超聲非線性輸出信號(hào)，合理的表征疲勞損傷過程中裂紋的群體行為和演化規(guī)律是非線性超聲檢測(cè)技術(shù)的關(guān)鍵問題。

混沌和分形理論的發(fā)展為深入分析非線性時(shí)間序列的內(nèi)在機(jī)理提供理論基礎(chǔ)，廣泛應(yīng)用于機(jī)械設(shè)備故障診斷和特征提取的研究中。孫自強(qiáng)[9]利用混沌分形理論對(duì)風(fēng)力機(jī)械傳動(dòng)系統(tǒng)的振動(dòng)數(shù)據(jù)進(jìn)行分析，對(duì)相關(guān)的故障診斷和狀態(tài)預(yù)測(cè)進(jìn)行研究，并把分維數(shù)可以作為風(fēng)電運(yùn)行狀態(tài)以及不同故障識(shí)別的一種有效手段；Takuma等[10]計(jì)算剪切加工過程中聲發(fā)射信號(hào)的關(guān)聯(lián)維數(shù)和最大Lyapunov 指數(shù)，并對(duì)刀具的磨損狀態(tài)和剩余壽命進(jìn)行評(píng)估。然而這些研究都是對(duì)系統(tǒng)的振動(dòng)信號(hào)或與振動(dòng)有關(guān)的信號(hào)進(jìn)行特性分析，振動(dòng)信號(hào)包含了較多的環(huán)境噪聲，且數(shù)據(jù)量比較大，對(duì)系統(tǒng)的特征分析和提取帶來了一定的困難。而超聲非線性輸出信號(hào)是反映系統(tǒng)特性的另一種時(shí)間序列，受環(huán)境影響較??；超聲波在試件中傳播與缺陷發(fā)生非線性相互作用，輸出信號(hào)能夠更直觀地反映系統(tǒng)的內(nèi)在特性，因此，在對(duì)試件疲勞損傷特征進(jìn)行分析時(shí)，超聲非線輸出信號(hào)可以有效地替代振動(dòng)信號(hào)。Wang等[11-12]探究了復(fù)雜應(yīng)力狀態(tài)下高溫低周疲勞表面短裂紋行為，并用計(jì)盒維數(shù)對(duì)裂紋的群體行為進(jìn)行表征。由于不同的混沌特征值可以從不同的角度反映非線性系統(tǒng)的特征，為了充分、全面地分析金屬疲勞損傷過程中裂紋擴(kuò)展與演化特征，對(duì)超聲非線性輸出信號(hào)的多組特征值進(jìn)行分析，深入探究金屬疲勞過程的特性。

本文利用RITEC公司的RAM-5000系統(tǒng)獲得具有不同疲勞損傷程度試件的超聲非線性輸出信號(hào)，在相空間重構(gòu)的基礎(chǔ)上，利用小數(shù)量法、最大似然法和GP算法估算超聲非線性輸出信號(hào)的Lyapunov指數(shù)、K熵、關(guān)聯(lián)維數(shù)，探究金屬疲勞損傷過程中混沌分形特征值的演化規(guī)律，通過對(duì)混沌分形特征值的分析深入地理解金屬材料疲勞過程的裂紋擴(kuò)展和演化規(guī)律，通過建立混沌分形特征值-疲勞損傷的關(guān)系，為金屬構(gòu)件疲勞壽命的預(yù)測(cè)提供了一種新的分析方法。

1 混沌分形特征值估算方法

混沌吸引子是混沌系統(tǒng)最顯著的特征，它表示系統(tǒng)最終會(huì)做規(guī)律的運(yùn)動(dòng)，產(chǎn)生規(guī)則、有形的軌跡。這種軌跡在經(jīng)過類似拉伸、折疊后轉(zhuǎn)化成與時(shí)間相關(guān)的序列，表現(xiàn)出混亂、復(fù)雜的特性。相空間重構(gòu)技術(shù)利用合適的嵌入維數(shù)m和延遲時(shí)間τ，重構(gòu)m維的相空間Y(ti)，恢復(fù)吸引子的特性。通過對(duì)吸引子特征量——Lyapunov指數(shù)、K熵、關(guān)聯(lián)維數(shù)等的分析，深入了解非線性系統(tǒng)的內(nèi)在特征[10]。

1.1Lyapunov指數(shù)

混沌運(yùn)動(dòng)最基本的特征是對(duì)初值條件極為敏感，兩個(gè)初始臨近點(diǎn)產(chǎn)生的軌道，隨著時(shí)間的推移按指數(shù)方式分離，Lyapunov指數(shù)就是研究這一特征的量。根據(jù)超聲非線性輸出信號(hào)的特性，計(jì)算Lyapunov指數(shù)采用小數(shù)據(jù)量法，它是利用軌道上每對(duì)初始臨近點(diǎn)的平均發(fā)散速率估計(jì)的。若時(shí)間序列{x1,x2,…,xN}，當(dāng)潛入維數(shù)為m，時(shí)間延遲為τ時(shí)，重構(gòu)相空間為

Yi=(xi,xi+τ,…xi+(m-1)τ)òRm，i=1,2,…,M

(1)

在相空間中，尋找給定軌道上點(diǎn)Yj的最近鄰點(diǎn)Yj′，并限制短暫分離，即

j-j′>p

(2)

P為時(shí)間序列的平均周期。每對(duì)初始臨近點(diǎn)，經(jīng)過i個(gè)離散步后的距離dj(i)，

dj(i)=Yj+i-Yj′+i,i=1,2,…min(M-j,M-j′)

(3)

對(duì)于每個(gè)離散時(shí)間步i，求出所有j的lndj(i)的平均值y(i)，

(4)

式中:q為非零dj(i)的個(gè)數(shù)，并用最小二乘法對(duì)i-y(i)作回歸直線，該直線的斜率即為最大Lyapunov指數(shù)。

1.2K熵

Kolmogorov熵(K熵)是吸引子的一個(gè)重要特征值，它對(duì)系統(tǒng)的混亂程度進(jìn)行整體的度量[13]。在混沌系統(tǒng)中，由于非線性系統(tǒng)的初值敏感性使軌道的指數(shù)發(fā)散，因此K熵為一正值，其值越大，信息的損失速率越大，系統(tǒng)越復(fù)雜。本文采用最大似然法估算時(shí)間序列的K熵。

1.3關(guān)聯(lián)維數(shù)

分形維數(shù)是描述吸引子的另一個(gè)重要參數(shù)，可以對(duì)系統(tǒng)的非線性行為進(jìn)行定量描述，關(guān)聯(lián)維數(shù)就是其中的一種，它表示吸引子相空間結(jié)構(gòu)的復(fù)雜程度，因此，它可以很好地反映時(shí)間序列中蘊(yùn)含的系統(tǒng)信息。

重構(gòu)m維相空間中有N個(gè)點(diǎn)，如果兩點(diǎn)yj和yj間的歐式距離dij小于給定值臨界值r，稱有關(guān)聯(lián)的分量，有關(guān)聯(lián)的矢量對(duì)數(shù)所占的比例稱為關(guān)聯(lián)積分Cn(r)：

(5)

其中，θ(x)為Heaviside單位函數(shù)，其定義為：

(6)

由于關(guān)聯(lián)積分與給定臨界值存在如下關(guān)系：

(7)

則D就是關(guān)聯(lián)維數(shù)，其計(jì)算式如下

(8)

用最小二乘法擬合雙對(duì)數(shù)曲線，該線的最佳擬合直線段區(qū)為無標(biāo)度區(qū)，無標(biāo)度區(qū)的斜率就是關(guān)聯(lián)維數(shù)。

2 試驗(yàn)系統(tǒng)與方法

2.1非線性超聲測(cè)試系統(tǒng)

圖1是非線性超聲測(cè)試系統(tǒng)原理圖，利用RAM-5000主機(jī)產(chǎn)生15個(gè)周期的正弦脈沖信號(hào)，經(jīng)衰減器和兩級(jí)低通濾波器后，驅(qū)動(dòng)中心頻率為5 MHz的激勵(lì)換能器，超聲波在試件中傳播與缺陷發(fā)生非線性作用，并由中心頻率為10 MHz的換能器接收，直接傳輸至RAM-5000主機(jī)的receiver2通道。帶通濾波器和前置放大器專門處理接收信號(hào)中微弱的二次諧波信號(hào)，并將處理后的信號(hào)傳輸至receiver1通道。采用黃油作耦合劑，同時(shí)也要保證激勵(lì)和接收換能器在同一軸線上，避免能量損失。利用示波器采集接收換能器的輸出時(shí)域信號(hào)(超聲非線性輸出信號(hào))，并對(duì)信號(hào)進(jìn)行256次平均，提高信號(hào)的信噪比。

圖1 非線性超聲測(cè)試系統(tǒng)原理圖

2.2試件制備

試件規(guī)格如圖2所示，材料選用45#鋼，其長(zhǎng)、寬、厚度分別為130 mm、30 mm、15 mm，用線切割技術(shù)在試件中部加工5 mm的機(jī)械缺口，線切割鉬絲直徑為0.18 mm，試件表面進(jìn)行拋光處理。試件共6根，分別編號(hào)1,2,3,4,5和6。在三點(diǎn)彎曲疲勞試驗(yàn)中，兩支點(diǎn)的跨度為120 mm，支點(diǎn)1距試件左邊緣的距離為5 mm，支點(diǎn)和加載位置如圖2所示。

三點(diǎn)彎曲疲勞試驗(yàn)在長(zhǎng)春新科PX20高頻疲勞試驗(yàn)機(jī)上完成，其參數(shù)設(shè)置如表1所示。

首先對(duì)試件5和6進(jìn)行疲勞試驗(yàn)，測(cè)得疲勞循環(huán)壽命分別為1.78×105和1.85×105次，并根據(jù)經(jīng)驗(yàn)確定該批試件的疲勞循環(huán)壽命為1.80×105次。對(duì)試件1、2、3、4進(jìn)行中斷疲勞試驗(yàn)，當(dāng)加載至相應(yīng)的疲勞次數(shù)時(shí)(疲勞循環(huán)壽命的百分比為10%、20%、30%、40%、50%、60%、70%、80%、90%)，暫停試驗(yàn)機(jī)，對(duì)試件進(jìn)行超聲激勵(lì)試驗(yàn)，獲得超聲非線性輸出信號(hào)，然后再對(duì)試件繼續(xù)進(jìn)行疲勞試驗(yàn)，直至試件失效。每次超聲試驗(yàn)均重復(fù)進(jìn)行6次，以減小隨機(jī)誤差。

圖2 試件/mm

設(shè)備參數(shù)設(shè)定值加載類型正弦交變應(yīng)力加載頻率自適應(yīng)最大加載力/kN-1．3最小加載力/kN-13

3 試驗(yàn)結(jié)果與分析

3.1非線性超聲測(cè)試系統(tǒng)性能分析

在非線性超聲測(cè)試系統(tǒng)中，試驗(yàn)儀器、耦合劑、換能器等都會(huì)產(chǎn)生高次諧波，對(duì)測(cè)量結(jié)果造成影響。

當(dāng)激勵(lì)信號(hào)頻率和超聲波傳播距離一定時(shí)，基波幅值的平方和二次諧波幅值成正比，據(jù)此可以判定試驗(yàn)系統(tǒng)的可靠性[14]。在其他條件相同的情況下，對(duì)試件3和4進(jìn)行超聲激勵(lì)試驗(yàn)，激勵(lì)電壓分別為50 V、60 V、70 V、80 V、90 V、100 V，分析基波幅值平方與二次諧波幅值的關(guān)系，實(shí)驗(yàn)結(jié)果如圖3所示。

圖3 基波幅值和二次諧波幅值平方的關(guān)系

3.2超聲非線性輸出信號(hào)的混沌分形特征值計(jì)算

圖4是示波器采集的超聲非線性輸出信號(hào)，該信號(hào)包含了兩端的暫態(tài)過程和中間的穩(wěn)態(tài)過程，為了抑制暫態(tài)過程對(duì)試驗(yàn)結(jié)果的影響，截取信號(hào)中間2 048個(gè)采樣點(diǎn)作為混沌分形分析的時(shí)間序列。

計(jì)算吸引子的特征值時(shí)，首先對(duì)時(shí)間序列進(jìn)行相空間重構(gòu)，將其拓展到高維空間，恢復(fù)吸引子的動(dòng)力特性，本文采用C-C法對(duì)超聲非線性輸出信號(hào)進(jìn)行相空間重構(gòu)。以試件4疲勞循環(huán)次數(shù)為30%的6組超聲非線性輸出信號(hào)為例，在最佳延遲時(shí)間和嵌入維數(shù)的基礎(chǔ)上，利用小數(shù)量法、最大似然法、GP算法分別計(jì)算吸引子的Lyapunov指數(shù)、K熵和關(guān)聯(lián)維數(shù)，將六組特征值的平均值作為該疲勞循環(huán)次數(shù)下的混沌分形特征值。圖5是Lyapunov指數(shù)和關(guān)聯(lián)維數(shù)的計(jì)算結(jié)果。

圖4 超聲非線性輸出信號(hào)

(a) Lyapunov指數(shù)

(b) 關(guān)聯(lián)維數(shù)

表2是試件4在不同疲勞損傷程度下的混沌分形特征值。采用上述方法依次計(jì)算試件1,2,3的混沌分形特征值。

從表2可以看出，在不同的疲勞損傷程度下，Lyapunov指數(shù)均大于0，說明疲勞損傷過程中超聲非線性輸出信號(hào)具有混沌特性。隨著疲勞循環(huán)次數(shù)的增加，Lyapunov指數(shù)、K熵和關(guān)聯(lián)維數(shù)均不斷增加，說明試件內(nèi)部的混沌特性不斷增強(qiáng)。

表2 試件4的混沌分形特征值

3.3混沌分形特征值的演化規(guī)律

為進(jìn)一步探究金屬疲勞損傷過程中混沌分形特征值的演化規(guī)律，對(duì)各特征值進(jìn)行歸一化處理。將未經(jīng)過疲勞加載試件的特征值記為λ0，經(jīng)過疲勞加載試件的特征值記為λ′，利用λ′/λ0進(jìn)行歸一化處理，可得到試件歸一化Lyapunov指數(shù)、歸一化關(guān)聯(lián)維數(shù)和歸一化K熵與疲勞循環(huán)壽命之間的關(guān)系，如圖6所示。

隨著疲勞循環(huán)次數(shù)的增加，混沌分形特征值逐漸變大，說明試件內(nèi)部的復(fù)雜性和混亂程度不斷增加。通過對(duì)比分析可知，四個(gè)試件的Lyapunov指數(shù)、關(guān)聯(lián)維數(shù)和K熵隨疲勞循環(huán)次數(shù)的變化趨勢(shì)大致相同，在疲勞循環(huán)的早期階段，混沌分形特征值單調(diào)遞增，在疲勞循環(huán)壽命為70%之后，各混沌分形特征值趨于穩(wěn)定。在三組變化曲線中Lyapunov指數(shù)和關(guān)聯(lián)維數(shù)的增加趨勢(shì)較為明顯，且在疲勞循環(huán)壽命為10%時(shí)，關(guān)聯(lián)維數(shù)的增幅約為17%，說明關(guān)聯(lián)維數(shù)對(duì)試件疲勞損傷的初期階段較為敏感。

圖6 歸一化特征值與疲勞循環(huán)壽命的關(guān)系

3.4疲勞損傷過程中微裂紋的演化規(guī)律

由于試件具有5 mm深的機(jī)械缺口，所以缺口處局部應(yīng)力較大。在循環(huán)載荷作用下，缺口的尖端形成應(yīng)力集中區(qū)域，會(huì)有大量的微裂紋萌生并伴隨有擴(kuò)展的現(xiàn)象，試件內(nèi)部結(jié)構(gòu)會(huì)隨著微裂紋的出現(xiàn)變復(fù)雜，因此試件的非線性特性增強(qiáng)，混沌分形特征值增大。由于關(guān)聯(lián)維數(shù)是系統(tǒng)混亂程度的度量，所以在大量的微裂紋萌生時(shí)，關(guān)聯(lián)維數(shù)有一個(gè)較大的增幅。隨著循環(huán)載荷的進(jìn)行，微裂紋繼續(xù)萌生和擴(kuò)展，并有一定數(shù)量的微裂紋合并，少量的裂紋擴(kuò)展較快，發(fā)展為主導(dǎo)裂紋，混沌分形特征值繼續(xù)增加。此階段中試件的損傷狀態(tài)已經(jīng)較為明顯，多條主導(dǎo)裂紋持續(xù)擴(kuò)展，且合體現(xiàn)象加速，形成致使試件失效的主裂紋。而在主裂紋的擴(kuò)展過程中，試件內(nèi)部不再有新的裂紋形成。

試件在疲勞循環(huán)壽命為70%左右時(shí)已出現(xiàn)明顯的宏觀裂紋，該宏觀裂紋由主裂紋發(fā)展而來，當(dāng)疲勞循環(huán)壽命在70%之后，宏觀裂紋尺寸進(jìn)一步變大，裂紋的數(shù)量和擴(kuò)展路徑的復(fù)雜度并未發(fā)生大的變化，所以混沌分形特征值在在疲勞循環(huán)壽命為70%時(shí)達(dá)到最大，而在疲勞循環(huán)的末期趨于穩(wěn)定狀態(tài)。因此，混沌分形特征值的變化可以有效地表征試件在疲勞損傷過程中裂紋的變化規(guī)律，揭示金屬疲勞損傷過程的非線性特性。

金屬試件在外界載荷作用下，混沌分形特征值逐漸增加，代表了裂紋形態(tài)的復(fù)雜程度及試件內(nèi)部應(yīng)力復(fù)雜程度的增加，即表面無序性的增長(zhǎng)。由于Lyapunov指數(shù)可以表征系統(tǒng)的混沌程度，關(guān)聯(lián)維數(shù)和K熵是系統(tǒng)混亂程度和復(fù)雜度的度量，而材料內(nèi)部結(jié)構(gòu)的復(fù)雜度反映了加載條件、所處環(huán)境和微觀結(jié)構(gòu)等多種因素，因此混沌分形特征值的演化規(guī)律可以表征材料的疲勞損傷過程，并可對(duì)試件的疲勞壽命進(jìn)行預(yù)測(cè)。

4 結(jié) 論

(1)在不同的疲勞損傷程度下，超聲非線性輸出信號(hào)的Lyapunov指數(shù)均大于0，說明蘊(yùn)含疲勞裂紋擴(kuò)展規(guī)律的超聲非線性輸出信號(hào)具有混沌特性。

(2)混沌分形特征值隨疲勞循環(huán)次數(shù)單調(diào)遞增，當(dāng)宏觀裂紋出現(xiàn)之后，三個(gè)特征值趨于飽和，說明在疲勞損傷過程中，材料內(nèi)部的混亂程度不斷增加。且關(guān)聯(lián)維數(shù)對(duì)疲勞損傷的初期階段較為敏感。

(3)通過分析微裂紋的演化特征可知，混沌分形特征值隨疲勞循環(huán)次數(shù)的變化曲線可以有效地表征疲勞損傷過程中微裂紋的演化規(guī)律和群體性行為。因此混沌分形理論可以有效地提取金屬疲勞損傷過程的特征，通過對(duì)疲勞過程中混沌分形特征值的分析，利用特征值的演化規(guī)律可以對(duì)試件的疲勞壽命進(jìn)行預(yù)測(cè)。

[1] 萬楚豪，剛鐵，劉斌. 非線性超聲脈沖反轉(zhuǎn)法在鋁合金焊縫疲勞壽命預(yù)測(cè)中的應(yīng)用[J]. 焊接學(xué)報(bào), 2015(2): 27-30.

WAN Chuhao, GANG Tie, LIU Bin. Nonlinear ultrasonic evaluation of fatigue life of aluminum alloy welded joint based on pulse-inversion technique[J]. Transactions of The China Welding Institution, 2015(2): 27-30.

[2] SHAH A A, RIBAKOV Y, ZHANG C. Efficiency and sensitivity of linear and non-linear ultrasonics to identifying micro and macro-scale defects in concrete[J]. Materials & Design, 2013, 50: 905-916.

[3] WU B, YAN B, HE C. Nonlinear ultrasonic characterizing online fatigue damage and in situ microscopic observation[J]. Transactions of Nonferrous Metals Society of China, 2011, 21(12): 2597-2604.

[4] SHUI G, WANG Y, GONG F. Evaluation of plastic damage for metallic materials under tensile load using nonlinear longitudinal waves[J]. NDT & E International, 2013, 55: 1-8.

[5] MATLACK K H, BRADLEY H A, THIELE S, et al. Nonlinear ultrasonic characterization of precipitation in 17-4PH stainless steel[J]. NDT & E International, 2015, 71: 8-15.

[6] 陳振華，史耀武，趙海燕，等. 微小缺陷的非線性超聲檢測(cè)及其成像技術(shù)[J]. 聲學(xué)學(xué)報(bào), 2010, 35(1): 9-13.

CHEN Zhenhua, SHI Yaowu, ZHAO Haiyan,et al. Nonlinear ultrasonic testing and imaging for tiny flaw[J]. Acta Acustica, 2010, 35(1): 9-13.

[7] 董雷. 基于混沌理論的非線性聲學(xué)特性研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2011.

[8] 王璐，王正，宋希庚，等. 基于分形理論的復(fù)雜應(yīng)力狀態(tài)下高溫低周疲勞表面短裂紋行為研究[J]. 機(jī)械工程學(xué)報(bào), 2011, 47(14): 49-53.

WANG Lu, WANG Zheng, SONG Xigeng, et al. Study on behavior of surface short cracks for low cycle at high temperature and complex stress state based on fractal theory[J]. Journal of Mechanical Engineering, 2011, 47(14): 49-53.

[9] 孫自強(qiáng). 基于混沌分形理論的大型風(fēng)電機(jī)械故障診斷研究[D]. 沈陽(yáng):沈陽(yáng)工業(yè)大學(xué), 2013: 114.

[10] TAKUMA M, SHINKE N, NISHIURA T, et al. Acoustic emission evaluation systems of tool life for shearing of piano and stainless steel wires[J]. Journal of Acoustic Emission, 2006, 24.

[11] WANG L, WANG Z, XIE W, et al. Fractal study on collective evolution of short fatigue cracks under complex stress conditions[J]. International Journal of Fatigue, 2012, 45: 1-7.

[12] 王魁. 基于分形理論的復(fù)雜應(yīng)力狀態(tài)下高溫低周疲勞短裂紋群體行為研究[D]. 大連: 大連理工大學(xué), 2011.

[13] 劉飛. 基于關(guān)聯(lián)維數(shù)和Kolmogorov熵的轉(zhuǎn)子振動(dòng)故障模式判別[D]. 沈陽(yáng): 沈陽(yáng)航空工業(yè)學(xué)院, 2009.

[14] 鄺文川. 基于非線性超聲縱波的高溫蠕變損傷檢測(cè)與評(píng)價(jià)研究[D]. 上海:華東理工大學(xué), 2011.

Featureanalysisofmetalfatiguedamageprocessesbasedonthechaosandfractaltheory

ZHANGYuhua1,LIXinxin2,HUANGZhenfeng2,MAOHanling2

(1. Institute of Light Industrial and Food Engineering, Guangxi University, Nanning 530004, China； 2. College of Mechanical Engineering, Guangxi University, Nanning 530004, China)

To effectively extract the characteristics of crack extension and evolution in the process of fatigue damages, the chaos and fractal theory were introduced to analyze the output signals of the receiver in a nonlinear ultrasonic system. The ultrasonic nonlinear output signals were gathered by the RAM-5000 system and analysed by virtue of the chaos and fractal theory. Then the Lyapunov exponent, correlation dimension and K entropy were calculated to evaluate the nonlinear characteristics of fatigue damages. The study shows that the ultrasonic nonlinear output signals of metal fatigue damages are of chaotic characters. The Lyapunov exponent, correlation dimension and K entropy monotonously increase with the development of fatigue damages, and reach a stable state when the macro cracks appear. The chaotic characteristic values can effectively characterize the extension and evolution of cracks and fully reveal the nonlinearity of metal fatigue processes. The results provide a new analysis method for the fatigue life analysis based on the relationship between the chaotic characteristic values and fatigue life.

nonlinear ultrasonic; fatigue damage; Lyapunov exponent; correlation dimension; K entropy

TG115.28

A

10.13465/j.cnki.jvs.2017.21.012

國(guó)家自然科學(xué)基金(51365006;51445013)

2016-05-17 修改稿收到日期：2016-09-04

張玉華 女,博士生,1987年4月生

毛漢領(lǐng) 男,博士,教授,1963年11月生,E-mail:maoh179@gxu.edu.cn