基于有限元分析的壓力容器損傷陣列稀疏特征建模與定位方法

朱高亮，王志凌，鐘永騰

(1.溫州大學(xué) 機電工程學(xué)院,溫州 325035；2.金陵科技學(xué)院 機電工程學(xué)院,南京 211169)

壓力容器在工業(yè)上有著廣泛的應(yīng)用，但惡劣的工作環(huán)境使得壓力容器設(shè)備易受到破壞，進(jìn)而導(dǎo)致結(jié)構(gòu)表面或內(nèi)部受損，嚴(yán)重影響設(shè)備的運行安全，因此，利用無損檢測技術(shù)對壓力容器的工作狀態(tài)進(jìn)行監(jiān)測、檢測顯得尤為重要[1-2]。常見的無損檢測技術(shù)有射線檢測[3]、電磁檢測[4]、滲透檢測[5]、渦流檢測[6]等。與常規(guī)無損檢測技術(shù)相比，超聲導(dǎo)波具有傳播距離遠(yuǎn)、能量衰減小以及傳播面積大等優(yōu)勢，應(yīng)用較為廣泛。LUGOVTSOVA 等[7]利用超聲導(dǎo)波與傳感器的合理配置，開發(fā)了一種復(fù)合材料壓力容器健康監(jiān)測方法。劉偉成等[8]利用超聲波對海洋石油靜電脫水壓力容器進(jìn)行了檢測，結(jié)果表明超聲檢測方法可對在役壓力容器進(jìn)行快速檢測。CHEN等[9]采用電磁超聲技術(shù)對儲氫壓力容器進(jìn)行在役監(jiān)測，為氫損傷的診斷和安全性評價提供了依據(jù)。

導(dǎo)波在壓力容器中傳播具有多頻散性，借助有限元分析軟件可以對導(dǎo)波傳播行為進(jìn)行可視化操作，便于更好地了解導(dǎo)波在壓力容器中經(jīng)過損傷或缺陷處的傳播變化，提高壓力容器超聲導(dǎo)波損傷監(jiān)測的效率。PARODI等[10]通過有限元建模分析預(yù)測了導(dǎo)波多模式的傳播行為，而后簡化了模型分析來研究復(fù)合材料壓力容器中導(dǎo)波與分層的相互作用。SAUSE等[11]對壓力容器的混合多層板中的聲發(fā)射源和信號傳播進(jìn)行有限元建模，模擬了多層結(jié)構(gòu)內(nèi)的導(dǎo)波傳播。因此，有限元分析軟件可以較容易地分析導(dǎo)播的傳播，對監(jiān)測方法優(yōu)化分析，可以提高對壓力容器的監(jiān)測效率。

針對壓力容器無損檢測中較難檢測的封頭部位，筆者提出一種基于有限元分析的壓力容器損傷陣列稀疏特征建模與定位方法，通過對壓力容器的封頭結(jié)構(gòu)進(jìn)行有限元建模，采取密集型陣列傳感器布置形式，提取表面稀疏分布的損傷模擬點源波傳播的陣列響應(yīng)信號波達(dá)時間差，建立結(jié)構(gòu)表面模擬損傷點源的樣本庫信號，利用表征損傷信息的均方根值之間的相似程度來得到損傷的具體位置。最后，在封頭表面含缺陷結(jié)構(gòu)的模型上采用基于Lamb波的主動監(jiān)測技術(shù)，接收并分析含缺陷結(jié)構(gòu)的陣列響應(yīng)信號，以此對提出的方法進(jìn)行驗證。

1 陣列稀疏特征損傷定位方法

1.1 壓力容器封頭結(jié)構(gòu)有限元建模

在ABAQUS軟件中建立壓力容器封頭結(jié)構(gòu)的有限元模型，封頭材料為30GrMo，其密度為7 850 kg·m-3，彈性模量為211 GPa，泊松比為0.279。仿真結(jié)構(gòu)為類橢圓半空心球體，該球體外圓直徑為635 mm，封頭底部到外徑圓圓心頂點的高度為140 mm，外徑圓處向外設(shè)置有延伸50 mm的貼合部設(shè)計。壓力容器結(jié)構(gòu)模型如圖1所示，其中，紅色實心點表示稀疏分布在結(jié)構(gòu)表面的損傷點源，對其進(jìn)行編號，1-1表示內(nèi)圈第一個圓編號為1的損傷模擬點，順時針排布并分別編號為1-2、1-3,…,1-120，共有120個稀疏排布的損傷模擬點源，黑色實心點表示監(jiān)測損傷的主動激勵點，黃色實心圓圈為接收信號圓形陣列。激勵源位于結(jié)構(gòu)正中心位置，信號接收陣列由8個陣元組成，設(shè)置在以中心為圓心，半徑為13 mm的圓上。

圖1 壓力容器結(jié)構(gòu)模型

1.2 構(gòu)建陣列稀疏特征

Lamb波是在結(jié)構(gòu)中傳播的一種彈性波，具有多頻散和多模態(tài)的特點，選取不同板厚和激勵頻率會產(chǎn)生不同的傳播模式[12]?；贚amb波的結(jié)構(gòu)主動監(jiān)測技術(shù)就是利用驅(qū)動器在結(jié)構(gòu)表面激發(fā)主動的檢測信號，與此同時，單個傳感器或陣列傳感器在結(jié)構(gòu)同側(cè)表面接收結(jié)構(gòu)的響應(yīng)信號，對接收的Lamb波信號進(jìn)行分析，進(jìn)而達(dá)到對結(jié)構(gòu)進(jìn)行健康狀態(tài)監(jiān)測的目的。研究中利用中心陣元PZT0激發(fā)Lamb波向外傳播，Lamb波經(jīng)過表面損傷會發(fā)生散射現(xiàn)象，采用8個陣元組成的陣列傳感器接收回波信號，陣列中陣元布置及損傷回波如圖2所示。

圖2 陣列中的陣元布置及損傷回波示意

Lamb波在傳播過程中呈現(xiàn)多模態(tài)且具有頻散效應(yīng)，因此傳播信號較為復(fù)雜,而窄帶信號激勵波形不易畸變，波峰明顯，頻散效應(yīng)較小，因此選用漢寧窗調(diào)制的五波峰窄帶信號作為激勵信號，可表示為

t∈[0，n/fc]

(1)

式中：A為激勵信號幅值；fc為中心頻率，取90 kHz;S為正弦調(diào)制信號;t為傳播時間;n為正弦調(diào)制信號的波峰數(shù)。

激勵信號的時域圖和頻域圖如圖3所示。

圖3 激勵信號的時域圖與頻域圖

在有限元仿真模型中，網(wǎng)格劃分的尺寸與導(dǎo)波在結(jié)構(gòu)中傳播的最小波長有關(guān)，根據(jù)文獻(xiàn)[13]，尺寸Δd應(yīng)滿足

Δd≤λmin/10

(2)

式中：λmin為Lamb波在模型中傳播的最短波長。

時間步長的準(zhǔn)確劃分能有效減少結(jié)果的誤差，并提高計算效率，根據(jù)文獻(xiàn)[13]，有限元模型求解的時間步長Δt應(yīng)滿足不大于最高頻率結(jié)果周期的1/20，可表示為

Δt≤1/(20fmax)

(3)

式中：fmax為頻帶中頻率的最大值。

綜合考慮，筆者采用的中心頻率為90 kHz，其對應(yīng)的波長為17.47 mm，在有限元模型中網(wǎng)格劃分的尺寸為1 mm，時間步長設(shè)置為0.1 μs，滿足最小尺寸網(wǎng)格劃分和最短時間步長的要求。

壓力容器封頭某損傷點源的Lamb波傳播有限元模擬結(jié)果如圖4所示，在沖擊點施加激勵信號，分析沖擊位置的Lamb波傳播變化，收集陣列響應(yīng)信號，通過Hilbert變換[14]做包絡(luò)曲線獲取陣列直達(dá)波信號的到達(dá)時刻，將其保存。同理，將壓力容器封頭表面布置的120個稀疏點分別通過導(dǎo)波激勵，將這120個稀疏點作為樣本庫中的120個樣本，同時收集陣列傳感器處的回波信號，獲取120個稀疏點的陣列直達(dá)波信號的到達(dá)時刻，并計算這120個稀疏點的到達(dá)時間均方根值(RMS)，形成基于有限元模型的損傷稀疏特征樣本庫。均方根ΔT可表達(dá)為

圖4 壓力容器封頭某損傷點源的Lamb波傳播有限元模擬結(jié)果

(4)

式中：n為陣列中的陣元數(shù)目；Δti為第i個傳感器陣元接收到的回波傳播信號到達(dá)時刻。

通過已建立的損傷稀疏特征樣本庫可知，在壓力容器封頭有限元模型表面任意部位的鄰近位置都存在一個模擬的損傷點源，當(dāng)接收到實測損傷信號后，通過上述方法進(jìn)行信號處理，最終可定位到與實測損傷信號最相近的損傷點源，后續(xù)將對該稀疏特征建模定位方法進(jìn)行驗證。

1.3 損傷定位方法

對波達(dá)時間進(jìn)行均方根值處理，最后與損傷稀疏特征樣本庫進(jìn)行匹配成像，對損傷進(jìn)行定位。散射回波信號的到達(dá)時刻應(yīng)由陣列接收到回波信號的時刻與陣列接收到中心激發(fā)的直達(dá)波到達(dá)時刻相減，即

Δti=Tsi-Tdi

(5)

式中：i=1,2,…,8,為陣列中的第i個陣元；Tsi為陣列信號圖中第i個陣元回波信號到達(dá)時刻；Tdi為陣列信號圖中第i個陣元損傷激勵源的直達(dá)波信號到達(dá)時刻。

將式(5)中得到的缺陷損傷散射回波信號的到達(dá)時刻與損傷稀疏特征樣本庫中的120個稀疏特征點進(jìn)行相應(yīng)陣元接收時刻相減，得到120組包含該缺陷損傷信息的散射回波信號的差值數(shù)據(jù)，將120組數(shù)據(jù)代入式(4)中，獲得相應(yīng)的到達(dá)時間均方根值。對數(shù)據(jù)進(jìn)行均一化處理，將其歸一化至[0，1]區(qū)間中，則與樣本庫中最相近的損傷稀疏特征點位置的值歸一化至區(qū)間最小值；最后，對其進(jìn)行成像，使歸一化值數(shù)據(jù)與成像的像素點成反比關(guān)系，即可得到損傷的成像位置。

2 仿真驗證

在有限元軟件建立的壓力容器封頭模型結(jié)構(gòu)表面(130 mm,-46°)位置(即編號2-10與3-10點位之間)處模擬一個尺寸為20 mm×10 mm(長×寬)的缺陷損傷，并分析圓形陣列接收到的信號變化。壓力容器封頭表面損傷有限元模型各時刻的Lamb波傳播情況如圖5所示，根據(jù)波的傳播可以看出，Lamb波在經(jīng)過損傷處時會發(fā)生散射，損傷散射波經(jīng)過傳播被陣列接收到回波信號，收集陣列傳感器回波信號的Lamb波到達(dá)時間，對其進(jìn)行均方根值處理，然后與損傷稀疏特征樣本庫進(jìn)行匹配成像，進(jìn)而對損傷進(jìn)行定位。

圖5 壓力容器封頭表面損傷有限元模型各時刻的Lamb波傳播情況示意

壓力容器封頭損傷的成像結(jié)果如圖6所示，像素點最高的位置即為損傷發(fā)生的位置。根據(jù)成像結(jié)果可知，基于壓力容器損傷陣列稀疏特征建模的定位方法進(jìn)行成像，缺陷位于(130 mm,-45°)位置,距離誤差為20 mm，角度誤差為1°，與實際損傷位置較為符合，能夠較準(zhǔn)確地定位出缺陷的位置。

圖6 壓力容器封頭損傷的成像結(jié)果

3 結(jié)語

(1) 建立了一種基于有限元分析的壓力容器損傷陣列稀疏特征建模與定位方法，對于真實工作條件下難以通過布置傳感器進(jìn)行檢測的結(jié)構(gòu)部位，在有限元分析軟件中建立結(jié)構(gòu)模型，提取結(jié)構(gòu)模型的陣列稀疏特征，實現(xiàn)了壓力容器封頭結(jié)構(gòu)表面損傷的定位。

(2) 相比較于傳統(tǒng)的損傷成像方法，該方法簡便易實施，不涉及復(fù)雜的成像原理，提取結(jié)構(gòu)模型的陣列稀疏特征即可對其進(jìn)行監(jiān)測，可為無損檢測新技術(shù)的發(fā)展提供參考，在損傷監(jiān)測應(yīng)用中具有一定的發(fā)展前景。