方漂漂　鄭慧峰　喻桑桑　竺科儀　王月兵　 曹文旭

1.中國計(jì)量學(xué)院，杭州，3100182.哈帝環(huán)?？萍加邢薰?，杭州，310018

基于振動(dòng)聲調(diào)制的金屬微裂紋檢測(cè)方法

方漂漂1鄭慧峰1喻桑桑1竺科儀2王月兵1曹文旭1

1.中國計(jì)量學(xué)院，杭州，3100182.哈帝環(huán)?？萍加邢薰荆贾?，310018

摘要：為克服傳統(tǒng)線性超聲無法檢測(cè)微納級(jí)寬度裂紋的缺點(diǎn)，開展了振動(dòng)聲調(diào)制法檢測(cè)金屬微裂紋的研究。用有限元分析軟件ABAQUS開展了仿真研究，并驗(yàn)證了該方法的有效性，為實(shí)驗(yàn)參數(shù)設(shè)置提供了參考。搭建振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)，分析了不同長(zhǎng)度微裂紋的調(diào)制特性，提出了基于振動(dòng)聲調(diào)制信號(hào)的信息熵定量表征方法。結(jié)果表明，存在微裂紋的鋁板的振動(dòng)聲調(diào)制信號(hào)中的旁瓣信號(hào)明顯，且檢測(cè)信號(hào)的信息熵隨裂紋長(zhǎng)度單調(diào)增大，可作為振動(dòng)聲調(diào)制法檢測(cè)微裂紋長(zhǎng)度的特征參數(shù)。

關(guān)鍵詞：微裂紋；振動(dòng)聲調(diào)制；非線性超聲；定量表征

0引言

常規(guī)超聲檢測(cè)法由于微裂紋的反射回波強(qiáng)度低、透射率高、衰減小等原因，對(duì)微裂紋不敏感[1]。國內(nèi)外研究者發(fā)現(xiàn)，非線性超聲檢測(cè)法可以很好地檢測(cè)閉合微裂紋、復(fù)合材料分層等材料內(nèi)部的接觸性缺陷。振動(dòng)聲調(diào)制是非線性超聲檢測(cè)方法中的一種，通過低頻振動(dòng)聲波和高頻超聲波與微裂紋相互作用，可使兩列聲波產(chǎn)生非線性耦合。它對(duì)接觸缺陷有較高的靈敏度，可以用于檢測(cè)復(fù)雜結(jié)構(gòu)零件和大型構(gòu)件及其結(jié)構(gòu)遠(yuǎn)端，所以正受到越來越多的關(guān)注[2]。

目前對(duì)振動(dòng)聲調(diào)制的研究還多限于調(diào)制信號(hào)的影響因素分析、檢測(cè)結(jié)果信號(hào)處理以及損傷的定量表征，還沒有建立統(tǒng)一的定量表征參數(shù)來對(duì)微裂紋進(jìn)行損傷評(píng)價(jià)。比較常用的評(píng)價(jià)損傷的參數(shù)是旁瓣信號(hào)幅值和主頻信號(hào)幅值，Aymerich等[3]利用一階旁瓣信號(hào)幅值表征損傷程度，發(fā)現(xiàn)旁瓣信號(hào)幅值隨著損傷程度增大而增大。但鄭磊[4]通過研究發(fā)現(xiàn)，旁瓣信號(hào)幅值和主頻信號(hào)幅值并不隨著裂紋長(zhǎng)度的增加而單調(diào)增大，因此旁瓣信號(hào)幅值和主頻信號(hào)幅值不適合作為裂紋長(zhǎng)度的定量表征參數(shù)。Duffour等[5]通過實(shí)驗(yàn)發(fā)現(xiàn)，旁瓣信號(hào)幅值與主頻信號(hào)幅值之比和裂紋長(zhǎng)度也不是單調(diào)相關(guān)的。Trochidis等[6]把基于ZAM算法得到的旁瓣信號(hào)和主頻信號(hào)的幅值之比作為裂紋深度的評(píng)價(jià)參數(shù)，它比基于FFT算法的幅值比對(duì)裂紋深度變化更加敏感，但仍然不是單調(diào)相關(guān)的，在工程實(shí)際中難以作為定量表征特征。

因此，本文根據(jù)振動(dòng)聲調(diào)制的理論，結(jié)合有限元仿真分析結(jié)果，搭建了振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)。以鋁板為檢測(cè)對(duì)象，研究不同激勵(lì)信號(hào)與不同長(zhǎng)度裂紋之間的調(diào)制結(jié)果，實(shí)現(xiàn)了振動(dòng)聲調(diào)制檢測(cè)技術(shù)對(duì)裂紋的有效檢測(cè)，并研究了信息熵與裂紋長(zhǎng)度之間的關(guān)系，為振動(dòng)聲調(diào)制的定量表征提供參考。

1振動(dòng)聲調(diào)制的檢測(cè)原理

振動(dòng)聲調(diào)制的基本原理如圖1所示[4]，在檢測(cè)過程中，向被測(cè)材料同時(shí)輸入高頻(ω1)和低頻(ω2)聲波，如果材料中沒有缺陷，則振動(dòng)聲調(diào)制信號(hào)為兩列聲波的線性疊加；如果存在缺陷，材料內(nèi)部空間在兩列聲波作用下會(huì)發(fā)生不均勻現(xiàn)象，缺陷處的質(zhì)點(diǎn)振動(dòng)速度改變，從而引起波的傳播速度發(fā)生畸變，產(chǎn)生以頻率ω1±ω2為主的非線性調(diào)制信號(hào)。

圖1　振動(dòng)聲調(diào)制原理

為了更加具體地描述振動(dòng)聲調(diào)制，下面將通過公式推導(dǎo)來解釋振動(dòng)聲調(diào)制的物理機(jī)理。

在一維條件下的各向同性體中，波動(dòng)方程(保留二階非線性項(xiàng))為[2]

(1)

式中，c為波速；β為非線性系數(shù)；u為質(zhì)點(diǎn)振動(dòng)位移；x為傳播距離；t為時(shí)間。

利用微擾近似理論，得到方程的近似解：

u(x,t)=u0(x,t)+βu1(x,t)

(2)

式中，u0為線性位移；u1為非線性位移。

暫不考慮衰減，則振動(dòng)聲調(diào)制信號(hào)的線性位移為

u0(x,t)=A1cosω1τ+A2cosω2τ

(3)

τ=t-x/c

式中，A1、A2分別為低頻聲波幅值和高頻聲波幅值。

假設(shè)非線性位移為

u1(x,t)=xh(τ)

(4)

式中，h(τ)為待定函數(shù)。

先將式(3)、式(4)代入式(2)，再代入式(1)，最后得到：

u(x,t)=u0(x,t)+βu1(x,t)=A1cosω1t+

βA1A2k1k2xsinω1tsinω2t/2

(5)

k1=ω1/c k2=ω2/c

從式(5)可以看出，兩種信號(hào)與微缺陷相互作用，產(chǎn)生了二次諧波信號(hào)和旁瓣信號(hào)，并且旁瓣信號(hào)的幅值正比于非線性系數(shù)和基波信號(hào)幅值。

2有限元分析

為了給實(shí)驗(yàn)研究提供指導(dǎo)依據(jù)，先用ABAQUS軟件對(duì)鋁板進(jìn)行振動(dòng)聲調(diào)制仿真，分析振動(dòng)聲調(diào)制超聲非線性特性，驗(yàn)證振動(dòng)聲調(diào)制對(duì)微小裂紋檢測(cè)的有效性。

2.1模態(tài)分析

根據(jù)Straka等[7]的研究，當(dāng)振動(dòng)信號(hào)頻率選用鋁板的模態(tài)頻率時(shí)，振動(dòng)聲調(diào)制現(xiàn)象最明顯，旁瓣信號(hào)幅值較大。所以用ABAQUS分別對(duì)尺寸為100mm×40mm，裂紋長(zhǎng)度為0～20mm的鋁板作模態(tài)分析。設(shè)置鋁板密度為2700kg/m3，彈性模量為72GPa，泊松比為 0.35。為簡(jiǎn)化計(jì)算，采用二維平面模型，裂紋設(shè)在鋁板模型中心，并使用線性攝動(dòng)分析步，特征值求解器選擇Lanczos，特征值數(shù)目設(shè)為30，模態(tài)分析結(jié)果如表1所示。

表1　模態(tài)分析結(jié)果(模態(tài)頻率)　Hz

從表1可知：在相同模態(tài)階次下，不同長(zhǎng)度裂紋鋁板的模態(tài)頻率變化很小。所以下面對(duì)有無裂紋鋁板的振動(dòng)聲調(diào)制仿真中，都選擇相同的低頻振動(dòng)頻率。有限元仿真時(shí)，為減小計(jì)算量，采用了二維平面模型，所以鋁板模型的質(zhì)量會(huì)比實(shí)際小，導(dǎo)致其固有頻率比實(shí)際值偏大，但在同一階次下模態(tài)分析結(jié)果的變化趨勢(shì)并沒有變化，所以仍然可以為實(shí)驗(yàn)操作提供參考。

2.2有限元建模及仿真

根據(jù)振動(dòng)聲調(diào)制原理，在鋁板模型左端邊界取2個(gè)點(diǎn)(高頻信號(hào)和低頻信號(hào)的輸入點(diǎn))，并施加正弦應(yīng)力波，使材料內(nèi)部質(zhì)點(diǎn)產(chǎn)生振動(dòng)，來模擬超聲波的激勵(lì)效果。由于裂紋在高低頻應(yīng)力波的共同作用下發(fā)生反復(fù)張開閉合，閉合的裂紋面將相互作用，所以需要將裂紋面定義為自接觸，接觸作用屬性為法向硬接觸。對(duì)于分析步設(shè)置，由于振動(dòng)聲調(diào)制是復(fù)雜非線性問題，所以選用顯式動(dòng)態(tài)分析。根據(jù)多次仿真經(jīng)驗(yàn)以及ABAQUS作業(yè)管理器監(jiān)控器輸出的信息，把分析步總時(shí)間設(shè)為5 ms，固定增量步時(shí)間設(shè)置為5 ns。低頻應(yīng)力波頻率設(shè)為4123 Hz，幅值為3 μm，高頻應(yīng)力波頻率設(shè)為50 kHz，幅值為1.5 μm。仿真結(jié)果如圖2、圖3所示。

(a)時(shí)域圖　　　　　　(b)頻域圖圖2　無裂紋鋁板仿真信號(hào)時(shí)域圖與頻譜圖

(a)時(shí)域圖　　　　　　(b)頻域圖圖3　有裂紋鋁板仿真信號(hào)時(shí)域圖與頻譜圖

從圖2、圖3可以看出：有缺陷鋁板的時(shí)域信號(hào)相比于無缺陷鋁板的時(shí)域信號(hào)，出現(xiàn)了明顯的調(diào)制。從頻域信號(hào)觀察，當(dāng)鋁板中無缺陷時(shí)，高頻信號(hào)周圍未出現(xiàn)明顯的旁瓣信號(hào)，但當(dāng)材料中存在長(zhǎng)為5 mm，寬為10 nm的微小裂紋時(shí)，在低頻信號(hào)的倍頻處出現(xiàn)了明顯的二次諧波A，在高頻信號(hào)周圍出現(xiàn)了頻率為(50 000±4123)Hz的旁瓣信號(hào)B、C。這與上文的振動(dòng)聲調(diào)制物理機(jī)理研究結(jié)果一致，表明振動(dòng)聲調(diào)制檢測(cè)法可以對(duì)鋁板中的微小缺陷進(jìn)行有效檢測(cè)。

2.3檢測(cè)結(jié)果影響因素分析

對(duì)于振動(dòng)聲調(diào)制檢測(cè)信號(hào)，調(diào)制現(xiàn)象越明顯，旁瓣信號(hào)的幅值就越大，測(cè)量結(jié)果準(zhǔn)確度越高。在仿真研究過程中發(fā)現(xiàn)，除了低頻信號(hào)的頻率會(huì)影響旁瓣信號(hào)的幅值外，基波信號(hào)幅值也會(huì)對(duì)旁瓣信號(hào)幅值有影響。

由式(5)可見，旁瓣信號(hào)幅值正比于非線性系數(shù)和基波信號(hào)幅值。高頻信號(hào)幅值分別為1 μm、1.5 μm、2 μm，低頻信號(hào)幅值從2 μm遞增到5 μm，高頻信號(hào)幅值、低頻信號(hào)幅值與旁瓣信號(hào)幅值之間的關(guān)系如圖4所示。其中低頻信號(hào)頻率為4123 Hz，高頻信號(hào)頻率為50 kHz。

圖4　低頻信號(hào)幅值與旁瓣信號(hào)幅值關(guān)系(高頻信號(hào)幅值為1 μm、1.5 μm、2 μm)

從圖4可以看出：高頻信號(hào)幅值越大，旁瓣信號(hào)幅值越大；旁瓣信號(hào)幅值與低頻幅值成正比，但在低頻信號(hào)幅值超過4.5 μm后增長(zhǎng)變緩。所以可以選取較大的激勵(lì)信號(hào)幅值,以獲得更好的振動(dòng)聲調(diào)制檢測(cè)效果。但同時(shí)要避免低頻信號(hào)幅值過大導(dǎo)致鋁板結(jié)構(gòu)變形嚴(yán)重，使得振動(dòng)聲調(diào)制信號(hào)信噪比降低。

高頻信號(hào)頻率和接收信號(hào)采樣位置同樣會(huì)對(duì)檢測(cè)結(jié)果有影響：高頻信號(hào)的頻率處于鋁板的模態(tài)頻率時(shí)，振動(dòng)聲調(diào)制效果最好，但其對(duì)結(jié)果的影響沒有低頻信號(hào)的影響明顯；采樣點(diǎn)和裂紋處于水平位置并且距離裂紋較近時(shí)，檢測(cè)到的非線性特性最好。

3實(shí)驗(yàn)研究

振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)如圖5所示，主要包括： 400 mm×200 mm×2 mm的薄鋁板，φ25 mm×2 mm、中心頻率為40 kHz的厚度伸縮型壓電陶瓷片，信號(hào)發(fā)生器DG1022U、示波器DS2102、功率放大器JYH-200M及工控機(jī)等。利用信號(hào)發(fā)生器產(chǎn)生一路高頻信號(hào)和一路低頻信號(hào)，通過功率放大器放大后，將兩路信號(hào)分別輸入到2個(gè)壓電陶瓷片產(chǎn)生超聲激勵(lì)和振動(dòng)聲激勵(lì)，再由另1個(gè)壓電陶瓷片接收振動(dòng)聲調(diào)制信號(hào)，最后由示波器采集傳送給計(jì)算機(jī)進(jìn)行數(shù)據(jù)處理并顯示。

圖5　振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)示意圖

本實(shí)驗(yàn)中，需要疲勞實(shí)驗(yàn)機(jī)對(duì)鋁板進(jìn)行周期性加載，形成微裂紋，但因?qū)嶒?yàn)室設(shè)備有限，所以采用過盈配合的方式模擬微裂紋[8]。先在鋁板中利用線切割加工一個(gè)公差等級(jí)為H6的圓形槽，然后在槽中塞入同樣材料、公差等級(jí)為n5的塞片，把塞片與槽的接觸面作為裂紋，如圖 6所示，塞片的直徑作為裂紋的長(zhǎng)度。

圖6　有裂紋鋁板示意圖

實(shí)驗(yàn)前，先要確定鋁板的模態(tài)頻率，用上述振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)對(duì)鋁板模態(tài)進(jìn)行測(cè)試，發(fā)現(xiàn)在1000～5000 Hz頻率范圍內(nèi)存在兩處振動(dòng)幅值較大的頻率(1200 Hz、3500 Hz)。通過仿真得知，在同一模態(tài)階次下，裂紋對(duì)鋁板的模態(tài)頻率影響不大，所以可將無裂紋鋁板的模態(tài)頻率作為低頻振動(dòng)聲信號(hào)的頻率。

用振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)分別對(duì)無裂紋和裂紋長(zhǎng)度為10 mm的鋁板進(jìn)行檢測(cè)。由信號(hào)發(fā)生器、功率放大器和2個(gè)壓電陶瓷片產(chǎn)生振動(dòng)聲激勵(lì)和超聲激勵(lì)，其頻率分別為1200 Hz和40 kHz，電壓為100 V和50 V。再由壓電陶瓷片接收振動(dòng)聲調(diào)制信號(hào)，并用示波器采集傳送給計(jì)算機(jī)，得到圖7、圖8所示的振動(dòng)聲調(diào)制信號(hào)。

(a)時(shí)域圖　　 　　　　(b)頻域圖圖7　無裂紋鋁板振動(dòng)聲調(diào)制信號(hào)

(a)時(shí)域圖　　 　　　　(b)頻域圖圖8　有裂紋鋁板振動(dòng)聲調(diào)制信號(hào)

從圖7、圖8可以看出：無裂紋鋁板的振動(dòng)聲調(diào)制信號(hào)沒有明顯的旁瓣信號(hào)，有微裂紋鋁板的檢測(cè)結(jié)果中出現(xiàn)了非常明顯的一階旁瓣信號(hào)(頻率為(40±1.2)kHz)和二階旁瓣信號(hào)(頻率為(40±2.4)kHz)，二階旁瓣信號(hào)幅值小于一階旁瓣信號(hào)幅值。所以，可以根據(jù)有無旁瓣信號(hào)產(chǎn)生來定性判斷是否存在裂紋，但還不能對(duì)裂紋的尺寸進(jìn)行定量評(píng)價(jià)。

4微裂紋的定量表征方法

裂紋尺寸的定量表征一直是研究難題，由于振動(dòng)聲調(diào)制物理機(jī)理復(fù)雜，所以一直沒有確定統(tǒng)一的表征參數(shù)[9]。信息熵可以用來描述信號(hào)的不規(guī)則度和復(fù)雜性，系統(tǒng)的混沌程度和復(fù)雜程度越大，熵值越大。一般故障信號(hào)和非故障信號(hào)相比，故障信號(hào)的復(fù)雜度較大，所以可以利用信息熵進(jìn)行診斷[10]。如果概率分布p(xi)(i=1，2，…，n)，記為p1、p2、…、pn，則熵函數(shù)為

(6)

其中，底數(shù)B取不同值，信息熵有不同的量綱。B=2時(shí)，量綱為bit(比特)；B=e時(shí)，量綱為nat(納特)，本文取B=e。

通過經(jīng)驗(yàn)?zāi)Ｊ椒纸?Empiricalmodedecomposition，EMD)方法把有裂紋鋁板振動(dòng)聲調(diào)制時(shí)域信號(hào)分解為8個(gè)固有模態(tài)函數(shù)(intrinsicmodefunction，IMF)，如圖9所示。然后把所得的IMF1分量作Hilbert變換得到其瞬時(shí)頻率和瞬時(shí)幅值，如圖10所示。從圖10可以看出：瞬時(shí)幅值隨著時(shí)間周期性變化，瞬時(shí)頻率以頻率40kHz為中心上下跳動(dòng)，說明振動(dòng)聲調(diào)制中存在頻率和幅值兩種調(diào)制。在有裂紋的鋁板中，裂紋與兩列聲波發(fā)生頻率和幅值調(diào)制，會(huì)增加振動(dòng)聲調(diào)制信號(hào)的復(fù)雜程度，所以熵值會(huì)比無裂紋的鋁板的調(diào)制信號(hào)大。因此，在振動(dòng)聲調(diào)制中可利用信息熵進(jìn)行損傷判斷。

圖9　振動(dòng)聲調(diào)制信號(hào)EMD分解結(jié)果

(a)瞬時(shí)幅值　　 　　　　(b)瞬時(shí)頻率圖10　有裂紋鋁板振動(dòng)聲調(diào)制信號(hào)瞬時(shí)幅值、瞬時(shí)頻率

為了研究信息熵能否對(duì)裂紋長(zhǎng)度進(jìn)行定量表征，在保持其他參數(shù)不變的情況下，使振動(dòng)聲激勵(lì)電壓分別為50 V、70 V、100 V，用振動(dòng)聲調(diào)制系統(tǒng)對(duì)裂紋長(zhǎng)度為0～15 mm的鋁板進(jìn)行檢測(cè)。然后對(duì)所得的振動(dòng)聲調(diào)制檢測(cè)信號(hào)頻域進(jìn)行采樣，并計(jì)算采樣數(shù)組的概率分布，再代入式(6)求取熵值，得到的結(jié)果如圖11所示。從圖11可知：在不同電壓的低頻振動(dòng)聲激勵(lì)下，調(diào)制信號(hào)的信息熵都隨著裂紋長(zhǎng)度的增加而線性單調(diào)增大；低頻動(dòng)信號(hào)幅值越大，振動(dòng)聲調(diào)制現(xiàn)象越明顯，所以在同樣裂紋長(zhǎng)度下，隨著低頻振動(dòng)信號(hào)幅值的增大，信息熵也增大。

圖11　信號(hào)熵值與裂紋長(zhǎng)度的關(guān)系

為了驗(yàn)證這種關(guān)系的正確性，在低頻信號(hào)的頻率1200 Hz、電壓100 V，高頻信號(hào)的頻率40 kHz、電壓50 V情況下，對(duì)裂紋長(zhǎng)度為2 mm、5 mm、8 mm、11 mm、14 mm的鋁板分別進(jìn)行振動(dòng)聲調(diào)制檢測(cè)，經(jīng)計(jì)算得到信息熵，如表2所示。根據(jù)圖11的擬合直線，查詢到對(duì)應(yīng)的裂紋長(zhǎng)度，經(jīng)計(jì)算可得所測(cè)數(shù)據(jù)的相對(duì)測(cè)量誤差均小于3%。因此，可利用信息熵對(duì)裂紋長(zhǎng)度進(jìn)行定量評(píng)估。

表2　信息熵對(duì)裂紋長(zhǎng)度的定量表征結(jié)果

5結(jié)論

(1)本文通過公式推導(dǎo)分析了振動(dòng)聲調(diào)制法檢測(cè)微裂紋的機(jī)理，發(fā)現(xiàn)振動(dòng)聲調(diào)制信號(hào)中會(huì)出現(xiàn)正比于材料非線性的諧波信號(hào)，并通過有限元分析軟件ABAQUS驗(yàn)證了該方法的有效性。

(2)搭建了一套振動(dòng)聲調(diào)制檢測(cè)系統(tǒng)，對(duì)鋁板中的微小裂紋進(jìn)行了檢測(cè)，分析其調(diào)制特性。研究發(fā)現(xiàn)，裂紋的存在會(huì)增加信號(hào)的復(fù)雜程度。因此，可將調(diào)制信號(hào)的信息熵作為損傷評(píng)價(jià)指標(biāo)。

參考文獻(xiàn)：

[1]吳斌，顏丙生，李佳銳，等.鎂合金疲勞早期非線性超聲在線檢測(cè)實(shí)驗(yàn)研究[J].聲學(xué)學(xué)報(bào)，2011，36(5)：527-533.

Wu Bin，Yan Binsheng，Li Jiarui，et al.Experimental Study of Early Fatigue Damage Online Measurement in Magnesium Using Nonlinear Ultrasonic[J].Acta Acustica，2011，36(5)：527-533.

[2]胡海峰.板狀金屬結(jié)構(gòu)健康監(jiān)測(cè)的非線性超聲理論與關(guān)鍵技術(shù)研究[D].長(zhǎng)沙：國防科學(xué)技術(shù)大學(xué)，2011.

[3]Aymerich F，Staszewski W J.Impact Damage Detection in Composite Laminates Using Nonlinear Acoustics[J].Composites Part A：Applied Science and Manufacturing，2010，41(9)：1084-1092.

[4]鄭磊.管道閉合裂紋振動(dòng)聲調(diào)制檢測(cè)方法研究[D].北京：北京工業(yè)大學(xué)，2012.

[5]Duffour P，Morbidini M，Cawley P.Comparison between a Type of Vibro-acoustic modelat Ion and Damping Measurement as ND Techniques[J].NDT & E International，2006，39(2)：123-131.

[6]Trochidis A，Hadjileontiadis L，Zacharias K. Analysis of Vibroacoustic Modulations for Crack Detection：a Time-frequency Approach Based on Zhao-Atlas-Marks Distribution[J].Shock and Vibration，2014，2014(10)：1-8.

[7]Straka L，Yagodzinskyy Y，Landa M. Detection of Structural Damage of Aluminum Alloy 6082 Using Elastic Wave Modulation Spectroscopy[J].NDT & E International，2008，41(7)：554-563.

[8]Tirosh J，Catz E.Mixed-mode Fracture Angle and Fracture Locus of Materials Subjected to Compressive Loading[J].Engineering Fracture Mechanics，1981，14(1)：27-38.

[9]Gurbatov S N，Rudenko O V，Saichev A L，et al.Waves and Structures in Nonlinear Nondispersive Media[M].Beijing：Higher Education Press，2011.

[10]劉國華，吳志根.引入信息熵理論的砼結(jié)構(gòu)損傷動(dòng)力識(shí)別新思路[J].振動(dòng)與沖擊，2011，30(6)：162-171.

Liu Guohua，Wu Zhigen.New Thought on Dynamic Identification Technology for Damage Detection of RC Structures by Introducing Information Entropy Theory[J].Journal of Vibration and Shock，2011，30(6)：162-171.

(編輯張洋)

Metal Micro Crack Detection Based on Vibro-acoustic Modulation Method

Fang Piaopiao1Zheng Huifeng1Yu Sangsang1Zhu Keyi2Wang Yuebing1Cao Wenxu1

1.China Jiliang University, Hangzhou, 310018 2.Hadi Environmental Protection Technology Company, Hangzhou, 310018

Abstract:In order to overcome the shortages that traditional linear ultrasonic testing could not detect the micro-cracks, VAM was studied. The VAM was simulated with finite element software ABAQUS to demonstrate the validity of the proposed method. The parameter setting of experimental testing system was also depended on the simulation results. VAM experimental testing system was established to investigate the nonlinear characteristic for micro cracks with variable lengths, and a quantitative characterization method was proposed based on information entropy of modulation signals. The results show that side-band signals are obvious when there is a micro crack in aluminum plate. The information entropy of modulation signals increases monotonically with increasing micro crack lengths and can be used as the characteristic parameter of crack length using VAM method.

Key words:micro crack；vibro-acoustic modulation (VAM)； nonlinear ultrasonic； quantitative characterization

收稿日期：2015-06-17

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(11474259)；浙江省自然科學(xué)基金資助項(xiàng)目(LY14E050013，LY15E050012)；浙江省公益技術(shù)應(yīng)用研究項(xiàng)目(2014C31109)；浙江省“儀器科學(xué)與技術(shù)”重中之重學(xué)科開放基金資助項(xiàng)目；工程車輛安全性設(shè)計(jì)與可靠性技術(shù)湖南省重點(diǎn)實(shí)驗(yàn)室資金資助項(xiàng)目(KF1601)

中圖分類號(hào)：TB523

DOI：10.3969/j.issn.1004-132X.2016.11.015

作者簡(jiǎn)介：方漂漂，女，1991年生。中國計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院碩士研究生。主要研究方向?yàn)槌暉o損檢測(cè)。鄭慧峰，男，1981年生。中國計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院副教授。喻桑桑，女，1989年生。中國計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院碩士研究生。竺科儀，男，1981年生。哈帝環(huán)?？萍加邢薰竟こ處?。王月兵，男，1963年生。中國計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院教授。曹文旭，女，1991年生。中國計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院碩士研究生。