收稿日期：2022-07-29""" 修回日期：2022-11-14

基金項目：國家自然科學(xué)基金青年項目（No.12102290；12201450）；山西省科技重大專項計劃“揭榜掛帥”資助項目（No.202101120401007）；山西省關(guān)鍵核心技術(shù)和共性技術(shù)研發(fā)攻關(guān)專項資助項目（No.2020XXX017）；山西省應(yīng)用基礎(chǔ)研究計劃面上青年項目（No.201901D211067）

通信作者：劉鵬。 E-mail：liupeng01@tyut.edu.cn

引用格式：

范朝珠，李志強(qiáng)，劉鵬，等. 改進(jìn)的基于貝葉斯框架和Lamb波的復(fù)合材料損傷定位方法［J］.應(yīng)用力學(xué)學(xué)報，2024，41（3）：499-508.

FAN Chaozhu，LI Zhiqiang，LIU Peng，et al. An improved method of damage localization of composites based on Bayesian frames and Lamb waves［J］.Chinese journal of applied mechanics，2024，41（3）：499-508.

文章編號：1000-4939（2024）03-0499-10

摘" 要：研究了一種改進(jìn)的在貝葉斯框架下將橢圓軌跡法與概率成像法進(jìn)行概率集成的重構(gòu)算法，通過研究Lamb波在復(fù)合材料層合板上的傳播特性，對重構(gòu)算法進(jìn)行改進(jìn)，解決了原有算法對于傳感器網(wǎng)絡(luò)邊緣的分層損傷定位不準(zhǔn)確的問題。本研究還以曲面復(fù)合材料層合板為例，進(jìn)行了幾種不同分層損傷工況下的有限元分析，提出了權(quán)重系數(shù)β，為該方法應(yīng)用在曲面復(fù)合材料層合板上提供了可靠的解決方法。結(jié)果表明，改進(jìn)的重構(gòu)算法可以實現(xiàn)對傳感器網(wǎng)絡(luò)邊緣分層損傷和曲面復(fù)合材料板分層損傷的精準(zhǔn)定位，定位的絕對誤差小于1cm。

關(guān)鍵詞：復(fù)合材料；分層損傷；貝葉斯；損傷成像；信號處理

中圖分類號：TP212" 文獻(xiàn)標(biāo)志碼：A

DOI：10.11776/j.issn.1000-4939.2024.03.002

An improved method of damage localization of composites based on Bayesian frames and Lamb waves

FAN Chaozhu1，LI Zhiqiang2，3，LIU Peng2，3，WANG Han2

（1.School of Mechanical and Vehicle Engineering，Taiyuan University of Technology，030024 Taiyuan，China;2.School of Aeronautics and Astronautics，Taiyuan University of Technology，030024 Taiyuan，China;3.National

Experimental Teaching Demonstration Center for Mechanics，Taiyuan University of Technology，030024 Taiyuan，China）

Abstract：In this paper，an improved reconstruction algorithm that integrates the elliptical trajectory method and the probabilistic imaging method under the Bayesian framework is studied，and the reconstruction algorithm is improved by studying the propagation characteristics of lamb waves on the composite laminate，which solves the problem that the original algorithm is inaccurate in the positioning of hierarchical damage at the edge of the sensor network.In this paper，the finite element analysis under several different layering damage conditions is carried out by taking the curved composite laminate as an example，and the weight coefficient is proposed to provide a reliable solution for the application of the method on the curved composite laminate.The results show that the improved reconstruction algorithm can achieve accurate positioning of layered damage at the edge of the sensor network and layered damage of curved composite plates，and the absolute error of positioning is less than 1cm.

Key words：composite; delaminated damage; Bayesian; damage imaging; signal processing

航空航天領(lǐng)域廣泛應(yīng)用到了復(fù)合材料，而復(fù)合材料結(jié)構(gòu)在生產(chǎn)與使用過程中可能會發(fā)生損傷，例如，纖維斷裂、基體開裂、層合板分層等［1］。為能及時發(fā)現(xiàn)這些可能存在的損傷，并判斷其位置，確定其程度，結(jié)構(gòu)健康監(jiān)測技術(shù)（structure health monitoring，SHM）作為一種在線、實時、快速的檢測方法成為當(dāng)前研究熱點(diǎn)。目前，基于Lamb波的檢測方法大多是根據(jù)結(jié)構(gòu)響應(yīng)信號的某些特征參數(shù)（速度、相位、幅值等）來評價結(jié)構(gòu)的健康狀況［2］REF_Ref102038116＼r＼h＼*MERGEFORMAT。經(jīng)過多年發(fā)展，該方法已經(jīng)可以判斷結(jié)構(gòu)是否發(fā)生損傷以及損傷發(fā)生的大致位置。但是想要達(dá)到定量的識別，精確地判斷出損傷發(fā)生的位置，還需要進(jìn)一步的研究。目前常用的SHM方法大致分為以下幾類。

1）基于飛行時間的橢圓軌跡法

基于飛行時間（time of flight，ToF）的方法是較先研究出來的一種方法，它的基本原理是通過獲得損傷散射波的波達(dá)時間、波速、傳播距離來計算出損傷發(fā)生的位置。劉洋等［3］REF_Ref102038137＼r＼h＼*MERGEFORMAT提出了一種橢圓軌跡算法對損傷進(jìn)行檢測。后來LIU等［4］REF_Ref102038157＼r＼h＼*MERGEFORMAT介紹了各向異性復(fù)合材料板的色散特性，提出了基于各向異性速度的橢圓軌跡法。HAMEED等［5］REF_Ref102038181＼r＼h＼*MERGEFORMAT在橢圓軌跡算法的基礎(chǔ)上，提出了一種基于歸一化Lamb波信號連續(xù)小波變換（continuous wavelet transform，CWT）的損傷估計方法，利用基于聯(lián)合和妥協(xié)融合方案的損傷指數(shù)構(gòu)建綜合損傷圖像。XU等［6］REF_Ref102038198＼r＼h＼*MERGEFORMAT提出了一種基于橢圓環(huán)的非線性Lamb波不完全群速度匹配的三維結(jié)構(gòu)微裂紋定位與成像方法。對于橢圓軌跡法來說最為嚴(yán)重的缺陷是存在識別盲區(qū)，因此，需要更多的PZT傳感器來產(chǎn)生更多橢圓軌跡來精確定位損傷。

2）基于成像算法的損傷識別

基于成像算法的識別方法可以把損傷通過成像的方式呈現(xiàn)出來，為損傷的定量識別和位置確定提供了便利。目前成像的方法包括：時間反轉(zhuǎn)法、多重信號分類法、RAPID（reconstruction algorithm for probabilistic inspection of damage）

層析成像法。HAY等［7］REF_Ref102038220＼r＼h＼*MERGEFORMAT在研究醫(yī)學(xué)CT成像的基礎(chǔ)上，提出了一種基于概率損傷的成像算法。HUANG等［8］REF_Ref102038238＼r＼h＼*MERGEFORMAT提出了一種利用多徑反射Lamb波進(jìn)行結(jié)構(gòu)損傷檢測的方法。QIU等［9-10］REF_Ref102038274＼r＼h＼*MERGEFORMAT提出了一種 Lamb 波層析成像的方法。XU等［11］REF_Ref102038300＼r＼h＼*MERGEFORMAT在此基礎(chǔ)上采用最小方差無失真響應(yīng)算法對稀疏系數(shù)進(jìn)行處理，從而生成圖像的Lamb波成像法。馮勇明等［12］REF_Ref102038317＼r＼h＼*MERGEFORMAT利用時間反轉(zhuǎn)法對復(fù)合材料損傷進(jìn)行檢測，并利用 RAPID 算法進(jìn)行損傷成像。對于RAPID法來說，其精度和分辨率高度依賴于傳感器陣列的密度。

3）基于貝葉斯融合法的損傷識別

橢圓軌跡法和RAPID方法的靈敏度和分辨率取決于傳感器陣列的密度和布局。因此，迫切需要一種既能提高損傷成像的靈敏度和分辨率，又能減少上述局限性的方法。SU等［13-15］REF_Ref102038339＼r＼h＼*MERGEFORMAT首次引入了概率密度函數(shù)來對傳統(tǒng)的基于ToF的損傷定位方法進(jìn)行改進(jìn)。劉龍等［16］提出了基于支持向量機(jī)的分步識別方法。YAN等［17］提出了將貝葉斯系統(tǒng)識別方法應(yīng)用于利用離散Lamb波的ToF測量的損傷定位問題。ZHANG等［18］采用多稀疏貝葉斯學(xué)習(xí)（M-SBL）策略進(jìn)行損傷成像，同時，將包括損傷反射波在內(nèi)的多個殘差信號分解為基于位置分量的稀疏矩陣。HUO等［19］提出了一種基于貝葉斯框架的橢圓軌跡法和快速缺陷概率檢測重建算法的概率集成方法，該方法考慮了多種損傷敏感特征的合理整合，提高了給定傳感器陣列的可靠性和魯棒性，但是在對傳感器網(wǎng)絡(luò)邊緣的損傷定位精度還有待提高。

本研究對貝葉斯框架得到的重構(gòu)算法進(jìn)行改進(jìn)后進(jìn)行損傷定位，解決了該方法在傳感器網(wǎng)絡(luò)邊緣定位不準(zhǔn)確的問題，并將其應(yīng)用到某機(jī)型后機(jī)身的復(fù)合材料蒙皮進(jìn)行分層損傷檢測。

1" 基于貝葉斯框架的損傷定位方法

在貝葉斯框架下，將橢圓軌跡法和RAPID法融合得到一種新的分層損傷識別方法。這種方法的基本思路是將RAPID法的概率密度函數(shù)作為先驗分布，再利用橢圓軌跡法的原理構(gòu)造似然函數(shù)，最后將兩者在貝葉斯框架下進(jìn)行概率融合得到其的后驗分布，便可以通過一組實測數(shù)據(jù)實現(xiàn)對損傷位置進(jìn)行參數(shù)估計，提高了損傷檢測的準(zhǔn)確性和有效性。

1.1" 橢圓軌跡法

在板間形成的Lamb波經(jīng)過損傷時會產(chǎn)生反射波，利用橢圓軌跡法進(jìn)行損傷檢測的關(guān)鍵在于對反射波飛行時間進(jìn)行測量，如圖1所示。

根據(jù)測得的反射波飛行時間，便可以估算出Dad+Dds的距離，那么損傷可能的位置就在以P1和P2為焦點(diǎn)的橢圓上，表示為

Dad+Dds=T×vg=" xd-xa2+yd-ya2+

xs-xa2+ys-ya2

（1）

式中：（xd，yd）是損傷坐標(biāo)；信號發(fā)射器P1和接收器P2的坐標(biāo)為（xa，ya）和xs，ys；vg是Lamb波的群速度；T是P2接收到損傷反射信號與直達(dá)信號的時間差。

將被檢測對象表面的每一個坐標(biāo)點(diǎn)視作一個像素點(diǎn)，可以得到每一個像素點(diǎn)（x，y）的損傷指數(shù)p（x，y），即

p（x，y）=∑Npi=1exp-T′i（x，y）-Tiτ0

（2）

式中：Np為傳感器路徑的總數(shù)；Ti是從實測的散射信號中得到的T；T′i（x，y）是結(jié)合損傷位置參數(shù)（xd，yd）理論計算得到的T；τ0是損傷因子，取0.01ms。損傷指數(shù)p（x，y）最大值對應(yīng)的坐標(biāo)點(diǎn)即為預(yù)測的最有可能發(fā)生損傷的位置。

1.2" 損傷概率重建法

RAPID方法的基本原理（圖2）是發(fā)射器-接收器的路徑上的損傷會對信號產(chǎn)生變化，而其變化幅度會隨著損傷與直接路徑之間距離的增加而減少。將發(fā)射器信號與接收器信號的相關(guān)系數(shù)作為損傷對該路徑信號的影響因子I，定義為

I=∑Kk=1Xk-μxYk-μy" ∑Kk=1Xk-μx2" ∑Kk=1Yk-μy2

（3）

式中：Xk為未損傷板接收器的信號；Yk為損傷板接收器的信號；μx、μy為其各自的均值；K為數(shù)據(jù)記錄的點(diǎn)數(shù)。

假設(shè)在被檢測對象的任意位置（x，y）處可以用一個含位置信息的損傷指數(shù)Pi（x，y）來表示損傷概率的大?。?0］REF_Ref102038479＼r＼h＼*MERGEFORMAT，基于線性減少的橢圓分布假設(shè)［21］REF_Ref102038497＼r＼h＼*MERGEFORMAT，在被檢測對象的任意位置（x，y）處的損傷概率可以通過以下公式估計，即

Px，y=∑Npi=1Pix，y

=∑Npi=1Ai-1β-1R+ββ-1

（4）

其中

R=RCx，y，xai，yai，xsi，ysiRC（x，y，xai，yai，xsi，ysi）lt;ββ，RC（x，y，xai，yai，xsi，ysi）gt;β，

RCx，y，xai，yai，xsi，ysi

=da+dsdas

=" x-xai2+y-yai2" xai-xsi2+yai-ysi2+

x-xsi2+y-ysi2" xai-xsi2+yai-ysi2

式中：RC為成像點(diǎn)（x，y）到發(fā)射器位置（xa，ya）的距離da和到接收器中心位置（xs，ys）的距離ds之和與發(fā)射器位置到接收器位置距離das的比值；Pi（x，y）為第i條傳感器路徑上存在損傷的概率估計；Ai=1-I，為第i條傳感器路徑信號差異系數(shù)，即損傷指標(biāo)；β為大于1的尺寸參數(shù)，控制著傳感器路徑影響區(qū)域的大小，本研究取1.07。

在基于Lamb波的RAPID方法中，由多個傳感器組成的傳感器網(wǎng)絡(luò)會形成多條發(fā)射器-接收器路徑。通過對每條發(fā)射器-接收器路徑確定的橢圓區(qū)域的概率進(jìn)行疊加，最終獲得待檢測的整個區(qū)域內(nèi)全部像素點(diǎn)出現(xiàn)損傷的概率。

1.3" 貝葉斯概率融合

首先在橢圓軌跡法的基礎(chǔ)上，利用實測的Ti與根據(jù)橢圓軌跡法求得的理論T′ixd，yd的差值構(gòu)造似然函數(shù)。

hθ|D=1" 2πσ2exp-Ti-T′ixd，yd22σ2

（5）

式中：θ為未知參數(shù)向量（xd，yd，σ）；D為實測T（T1，T2，…，TNp），Np是發(fā)射器-接收器路徑總數(shù)。

聯(lián)合先驗分布由RAPID法中的xd，yd的概率密度函數(shù)與上式中σ的概率密度函數(shù)得到

pxd，yd=1Np∑Npi=11-Iiβ-Ri（x，y）β-1（6）

根據(jù)貝葉斯理論便可以得到參數(shù)向量θ的后驗概率密度函數(shù)為

pθ|D∝pD|θ·pθ

=1σNp+1∑Npi=11-Iiβ-Ri（x，y）β-1·" exp-∑Npi=1Ti-T′ixd，yd22σ2（7）

從上式可以看到，在貝葉斯框架下得到的后驗分布是非常復(fù)雜的，這種情況下通常采用馬爾科夫鏈蒙特卡羅（MCMC）來對未知參數(shù)θ的分布進(jìn)行估計，即采用MCMC方法從后驗概率密度函數(shù)中提取樣本，利用樣本對θ的分布進(jìn)行預(yù)測［22］REF_Ref102038528＼r＼h＼*MERGEFORMAT。

1.4" 馬爾科夫鏈蒙特卡羅

面對僅知其表達(dá)式的復(fù)雜分布，它的均值、最大值、最小值等無法用常規(guī)方法求出，甚至對其進(jìn)行采樣都格外棘手。馬爾科夫鏈蒙特卡羅（MCMC）采樣方法允許通過構(gòu)造一個馬爾科夫鏈，并借助已知的標(biāo)準(zhǔn)分布，使得以某種接受概率采樣的馬爾科夫鏈的平穩(wěn)分布為目標(biāo)復(fù)雜分布。

假設(shè)狀態(tài)序列：…xt－1，xt，xt+1…，有

P（xt+1xi）=P（xt+1xt）（8）

則這組序列叫馬爾科夫鏈。當(dāng)有一個初始概率分布π使得：πP=π 時，那么π是該馬爾科夫鏈的平穩(wěn)分布，P是其狀態(tài)轉(zhuǎn)移矩陣。

狀態(tài)轉(zhuǎn)移矩陣P必須滿足馬爾科夫鏈的細(xì)致平衡方程，即

π（i）P（i，j）=π（j）P（j，i）（9）

任取Q使得P（i，j）=Q（i，j）α（i，j），即通過對一個隨機(jī)給定的矩陣Q進(jìn)行采樣，該樣本被留下來的概率為α。

根據(jù)馬爾科夫鏈的細(xì)致平衡方程有

α（i，j）=minπ（j）Q（j，i）π（i）Q（i，j），1（10）

由此，便可以根據(jù)目標(biāo)分布π，初始分布Q 來對目標(biāo)分布的未知參數(shù)進(jìn)行采樣。

2" 方法的改進(jìn)

目前，貝葉斯融合法對金屬板損傷定位的結(jié)果較好，但是，在復(fù)合材料層合板上的定位結(jié)果不太理想，特別是在傳感器網(wǎng)絡(luò)的邊緣區(qū)域會出現(xiàn)比較明顯的定位偏差［19］REF_Ref102038453＼r＼h＼*MERGEFORMAT。本研究基于目前貝葉斯融合法對于傳感器網(wǎng)絡(luò)邊緣分層損傷定位不準(zhǔn)確的不足，通過對復(fù)合材料層合板不同方向Lamb波的傳播速度進(jìn)行研究，對貝葉斯融合法進(jìn)行了改進(jìn)。

2.1" 復(fù)合材料對Lamb波的影響

本研究使用ABAQUS來對Lamb在復(fù)合材料層合板上的傳播進(jìn)行研究。該層合板是由4層復(fù)合材料鋪設(shè)而成，每一層的厚度為0.32mm，總厚度為1.28mm，鋪設(shè)角度為［45°/-45°/0°/90°］s對稱鋪層［23］。其材料性能如表1所示。

層合板尺寸為450mm×450mm×1.28mm，傳感器的布置如圖3所示。其中T為發(fā)射器，周圍R1～R5為接收器。

激勵信號是由漢寧窗調(diào)制的中心頻率為160kHz的窄帶五峰波正弦信號（圖4）。激勵信號通過在傳感器位置施加載荷來進(jìn)行模擬，仿真總時長為500μs，采樣頻率為100MHz。

以R1為例，對響應(yīng)信號（圖5）進(jìn)行小波變換后便可以得到Lamb波在該方向上的傳播時間（圖6），進(jìn)而可以求出Lamb的傳播速度（圖7）。

根據(jù)各方向Lamb的傳播速度，對成像算法進(jìn)行改進(jìn)，將不同方向Lamb波的速度與信號的傳播路徑對應(yīng)起來，以此來提高定位精度。

2.2" 改進(jìn)算法的驗證

使用有限元模擬方法，對相同尺寸和材料屬性的復(fù)合材料平板分層損傷進(jìn)行識別定位。

如圖8所示，板上放置4個傳感器網(wǎng)絡(luò)P1～P4，傳感器選擇壓電換能器，坐標(biāo)原點(diǎn)位于平板的左下角，各個傳感器的位置坐標(biāo)為P1（100，350）、P2（350，350）、P3（100，100）、P1（350，100），有限元網(wǎng)格為1mm×1mm×0.32mm的八節(jié)點(diǎn)六面體單元。由4個傳感器組成的監(jiān)測區(qū)域如圖9所示，共有6條不同的發(fā)射器-接受器路徑，通過對比不同位置損傷的定位結(jié)果來證明改進(jìn)方法的優(yōu)勢。

本研究設(shè)置了2個損傷工況，損傷大小都是半徑為10mm的圓形區(qū)域。各工況損傷位置的中心坐標(biāo)為：工況一（165，225）mm、工況二（290，300）mm。

對獲得到的響應(yīng)信號處理流程如圖10所示：無損傷平板上得到健康信號，損傷平板上得到損傷信號（圖10a）；獲得健康信號與損傷信號的相關(guān)系數(shù)DIi；將損傷信號與健康信號做差得到散射信號（圖10b）；將散射信號做連續(xù)小波變換（CWT）［24］REF_Ref102038589＼r＼h＼*MERGEFORMAT，T即為激勵信號與散射信號CWT后的第一個峰值時間（圖10c）。

得到損傷工況一每條路徑的相關(guān)系數(shù)I與T列于表2。

利用表2中的損傷特征以及貝葉斯后驗分布得到的損傷位置的概率密度函數(shù)，采用MCMC抽樣程序從后驗分布中生成樣本。在本研究中，所有從貝葉斯后驗中提取的變量的初始值都是使用極大似然估計獲得的。本例中共抽取了60000個變量樣本，將前10000個樣本設(shè)置為老化數(shù)據(jù)，共得到50000個樣本用來估計損傷發(fā)生的位置。圖11（a）（b）為估計的位置參數(shù)（x，y）的直方圖，根據(jù)直方圖可以構(gòu)造一個二維正態(tài)分布（圖11c）來表示估計的損傷位置，損傷識別結(jié)果如圖12（c）所示。

改進(jìn)前和改進(jìn)后兩種損傷工況下的定位結(jié)果如圖12所示。圖中黑色圓圈代表實際損傷發(fā)生的位置，紅色區(qū)域表示預(yù)測的損傷位置。兩種方法預(yù)測的損傷位置的中心坐標(biāo)見題注。圖12（a）和（d）為平板損傷模型，圖12 （b）和（e）為方法改進(jìn)前的損傷定位結(jié)果，圖12 （c）和（f）為方法改進(jìn)后的損傷定位結(jié)果，可以看出與之前的方法相比，改進(jìn)后的方法在損傷位置的定位上產(chǎn)生了更為準(zhǔn)確的結(jié)果。

3" 機(jī)身復(fù)合材料蒙皮分層損傷檢測

在實際應(yīng)用中，相對于平板而言，曲面層合板的應(yīng)用更為廣泛。特別是在航空航天領(lǐng)域復(fù)合材料扮演了重要的角色，許多民用和軍用飛機(jī)都大量使用了曲面復(fù)合材料層合板作為蒙皮，以發(fā)揮其高比強(qiáng)度、低重量、耐腐蝕和優(yōu)異的抗疲勞性能。

然而，該損傷定位方法與波的傳播路徑息息相關(guān)，將該方法用在曲面復(fù)合材料板上會造成嚴(yán)重的定位誤差。因此，本研究提出了權(quán)重系數(shù)β解決曲面層合板的成像問題。

3.1" 有限元模型的建立

某型號無人機(jī)全機(jī)采用復(fù)合材料蒙皮，本研究以該型號無人機(jī)的后機(jī)身頂部的機(jī)身蒙皮為例進(jìn)行分層損傷定位，檢測區(qū)域的尺寸為450mm×450mm×1.28mm，該蒙皮模型由4層復(fù)合材料鋪設(shè)而成，每一層的厚度為 0.32mm ，總厚度為1.28mm，鋪設(shè)角度為［45°/-45°/0°/90°］s。

該機(jī)型后機(jī)身復(fù)合材料蒙皮模型的建立如圖13［25］所示REF_Ref102038551＼r＼h＼*MERGEFORMAT。根據(jù)成像區(qū)域截取對應(yīng)的幾何模型，損傷區(qū)域以及傳感器如圖14所示。模型中上部紅色區(qū)域為預(yù)制分層損傷區(qū)域。本次研究設(shè)置了4種不同的損傷位置，分層位置都在第2層和第3層之間，損傷大小都是半徑為10mm的圓形區(qū)域。各損傷工況的損傷位置中心坐標(biāo)為損傷一（225，350）mm、損傷二（225，100）mm、損傷三（305，305）mm、損傷四（150，225）mm。

3.2" 損傷位置監(jiān)測

獲得完好模型和損傷模型響應(yīng)信號后，通過對信號進(jìn)行處理得到6條路徑各自的相關(guān)系數(shù)I與T。因為被檢測模型是曲面，而成像算法是在平面上進(jìn)行理論計算得到的，因此，需要對曲面上的數(shù)據(jù)進(jìn)行處理，即將實測T乘以一個權(quán)重系數(shù)β，將其投影到平面上。β的定義為

β=L′L（11）

式中：L′是發(fā)射器與接收器在成像平面距離；L是發(fā)射器與接收器在實際曲面的距離。各路徑的β如表3所示。處理后的T和損傷因子I如表4所示。

損傷識別結(jié)果如圖15所示。采用該方法識別的損傷位置中心坐標(biāo)為（233，345）mm，實際損傷的中心位置坐標(biāo)為（225，350）mm，其識別結(jié)果具有很高的準(zhǔn)確性。

對其余3個損傷工況進(jìn)行同樣的有限元分析來驗證該方法的有效性。各損傷工況的損傷特征如表5所示。

各損傷工況的損傷位置識別結(jié)果如圖16所示。將設(shè)定損傷的位置坐標(biāo)與求得的位置坐標(biāo)進(jìn)行比較，結(jié)果誤差為

E=" x2-x12+y2-y12（12）

式中：（x1，y1）為實際的損傷位置坐標(biāo);（x2，y2）為預(yù)測的損傷位置坐標(biāo)。誤差分析如表6所示。

由分析結(jié)果可以得出，該方法對于復(fù)合材料蒙皮分層損傷的檢測結(jié)果的定位誤差在可接受范圍之內(nèi)，由此表明該方法有一定的準(zhǔn)確性與可靠性。

4" 結(jié)" 論

本研究基于貝葉斯框架和Lamb波的復(fù)合材料結(jié)構(gòu)損傷定位方法，對貝葉斯框架下得到的重構(gòu)算法進(jìn)行改進(jìn)，并將其應(yīng)用到曲面復(fù)合材料層合板分層損傷進(jìn)行識別。該方法是在貝葉斯框架下將橢圓軌跡法和RAPID法進(jìn)行融合，形成新的有關(guān)損傷位置變量的概率密度函數(shù)，通過研究復(fù)合材料層合板的各向異性對Lamb波的影響對定位算法進(jìn)行改進(jìn)，最終利用改進(jìn)算法對位置變量進(jìn)行抽樣，從而識別損傷發(fā)生的位置。本研究還通過有限元仿真對4種不同損傷位置下的曲板模型進(jìn)行了損傷定位，對識別結(jié)果進(jìn)行了誤差分析，最終得到的誤差均在可接受范圍之內(nèi)，驗證了該模型在不同損傷位置的有效性，結(jié)果表明該方法能夠獲得準(zhǔn)確、可靠的檢測效果。通過研究解決了以下問題。

1）Lamb波在復(fù)合材料層合板上不同方向上的傳播速度的差異會影響傳感器網(wǎng)絡(luò)邊緣區(qū)域的定位精度，改進(jìn)后的方法獲得了更為可靠地定位結(jié)果。

2）對于實際應(yīng)用更為廣泛的曲面板來說，權(quán)重系數(shù)β的提出將曲面板上獲得的成像數(shù)據(jù)投影在了成像平面上，使數(shù)據(jù)和成像方法得以匹配，為該方法的實際應(yīng)用提供了可靠的解決方案。

參考文獻(xiàn)：

［1］" 郭雙喜，李雪芹.基于非線性彈簧界面的復(fù)合材料分層損傷預(yù)測［J］.應(yīng)用力學(xué)學(xué)報，2019，36（4）：958-965.

GUO Shuangxi，LI Xueqin.Delamination prediction in laminate composites based on nonlinear spring interface elements［J］.Chinese journal of applied mechanics，2019，36（4）：958-965（in Chinese）.

［2］" 王建強(qiáng)，余龍，周蘇楓，等.基于認(rèn)知技術(shù)的Lamb波結(jié)構(gòu)損傷檢測［J］.應(yīng)用力學(xué)學(xué)報，2017，34（6）：1102-1107.

WANG Jianqiang，YU Long，ZHOU Sufeng，et al.Lamb wave structure damage detection based on cognitive technology［J］.Chinese journal of applied mechanics，2017，34（6）：1102-1107（in Chinese）.

［3］" 劉洋，王乘，張鋒.基于壓電晶體傳感器的橢圓損傷定位算法［J］.傳感技術(shù)學(xué)報，2005，18（3）：596-600.

LIU Yang，WANG Cheng，ZHANG Feng.Piezoelectric wafer-based damage localization using ellipses method［J］.Chinese journal of sensors and actuators，2005，18（3）：596-600（in Chinese）.

［4］" LIU Y，HONG X B，ZHANG B.A novel velocity anisotropy probability imaging method using ultrasonic guided waves for composite plates［J］.Measurement，2020，166：108087.

［5］" HAMEED M S，LI Z，ZHENG K H.Damage detection method based on continuous wavelet transformation of Lamb wave signals［J］.Applied sciences，2020，10（23）：8610.

［6］" XU J C，ZHU W J，XIANG Y X，et al.Localization and imaging of micro-cracks using nonlinear Lamb waves with imperfect group-velocity matching［J］.Applied sciences，2021，11（17）：8069.

［7］" HAY T R，ROYER R L，GAO H D，et al.A comparison of embedded sensor Lamb wave ultrasonic tomography approaches for material loss detection［J］.Smart materials and structures，2006，15（4）：946-951.

［8］" HUANG L P，ZENG L，LIN J，et al.Baseline-free damage detection in composite plates using edge-reflected Lamb waves［J］.Composite structures，2020，247：112423.

［9］" QIU L，YUAN S F，ZHANG X Y，et al.A time reversal focusing based impact imaging method and its evaluation on complex composite structures［J］.Smart materials and structures，2011，20（10）：105014.

［10］HETTLER J，TABATABAEIPOUR M，DELRUE S，et al.Guided wave tomography based inspection of CFRP plates using a probabilistic reconstruction algorithm［C］//The 6th International Conference on Emerging Technologies in Nondestructive Testing.Brussels：［s.n.］，2015.

［11］XU C B，YANG Z B，ZUO H，et al.Minimum variance Lamb wave imaging based on weighted sparse decomposition coefficients in quasi-isotropic composite laminates［J］.Composite structures，2021，275：114432.

［12］馮勇明，楊建元，陳衛(wèi)東，等.基于概率成像的復(fù)合材料結(jié)構(gòu)健康監(jiān)測系統(tǒng)［J］.電子測量與儀器學(xué)報，2013，27（3）：253-257.

FENG Yongming，YANG Jianyuan，CHEN Weidong，et al.Structural health monitoring system for composite matrial based on probabilistic damage imaging［J］.Journal of electronic measurement and instrumentation，2013，27（3）：253-257（in Chinese）.

［13］SU Z Q，YE L，LU Y.Guided Lamb waves for identification of damage in composite structures：a review［J］.Journal of sound and vibration，2006，295（3/4/5）：753-780.

［14］ZHOU C，SU Z Q，CHENG L.Probability-based diagnostic imaging using hybrid features extracted from ultrasonic Lamb wave signals［J］.Smart materials and structures，2011，20（12）：125005.

［15］YU L，SU Z Q.Application of kernel density estimation in Lamb wave-based damage detection［J］.Mathematical problems in engineering，2012，2012：406521.

［16］劉龍，黃海，孟光.基于支持向量機(jī)的結(jié)構(gòu)損傷分步識別研究［J］.應(yīng)用力學(xué)學(xué)報，2007，24（2）：313-317.

LIU Long，HUANG Hai，MENG Guang.Multi-step structure damage identification approach based on support vector machine［J］.Chinese journal of applied mechanics，2007，24（2）：313-317（in Chinese）.

［17］YAN G.A Bayesian approach for damage localization in plate-like structures using Lamb waves［J］.Smart materials and structures，2013，22（3）：035012.

［18］ZHANG H，HUA J D，GAO F，et al.Efficient Lamb-wave based damage imaging using multiple sparse Bayesian learning in composite laminates［J］.NDT amp; E international，2020，116：102277.

［19］HUO H D，HE J J，GUAN X F.A Bayesian fusion method for composite damage identification using Lamb wave［J］.Structural health monitoring，2021，20（5）：2337-2359.

［20］ZHOU C，SU Z Q，CHENG L.Quantitative evaluation of orientation-specific damage using elastic waves and probability-based diagnostic imaging［J］.Mechanical systems and signal processing，2011，25（6）：2135-2156.

［21］NG C T，VEIDT M.A Lamb-wave-based technique for damage detection in composite laminates［J］.Smart materials and structures，2009，18（7）：074006.

［22］JASRA A.Exploring Monte Carlo methods［J］.Materials today，2012，15（1）：69.

［23］鐵瑛，徐平，李成.不同鋪層角含孔復(fù)合材料板的損傷研究［J］.應(yīng)用力學(xué)學(xué)報，2014，31（1）：156-161.

TIE Ying，XU Ping，LI Cheng.Damage study of laminate containing a hole with different ply orientation angles［J］.Chinese journal of applied mechanics，2014，31（1）：156-161（in Chinese）.

［24］陳換過，姜節(jié)勝，閆云聚，等.小波分析在復(fù)合材料層合結(jié)構(gòu)損傷檢測中的應(yīng)用［J］.應(yīng)用力學(xué)學(xué)報，2004，21（4）：140-142.

CHEN Huanguo，JIANG Jiesheng，YAN Yunju，et al.Wavelet analysis in damage detection for laminated composite structure［J］.Chinese journal of applied mechanics，2004，21（4）：140-142（in Chinese）.

［25］袁浩凡.基于振動的FRP弧形層合板中分層損傷的無損檢測研究［D］.廣州：廣州大學(xué)，2016.

（編輯" 李坤璐）