孟 浩，王志剛，劉昌明

（武漢科技大學(xué)冶金裝備及其控制教育部重點實驗室，湖北 武漢 430081）

1 引言

耐火磚在生產(chǎn)過程中容易產(chǎn)生斷層、裂紋等內(nèi)部損傷[1]，損傷嚴重時會極大降低耐火磚的使用性能。若將性能低劣的耐火磚投入使用，可能會造成嚴重的安全事故。因此，對耐火磚進行損傷識別具有重大意義。目前的耐火磚損傷檢測技術(shù)應(yīng)用較多的有超聲波檢測法，聲發(fā)射檢測法，沖擊回波檢測法[2]，這些方法雖然檢測精度高，但操作復(fù)雜，成本較高，無法滿足批量檢測的要求。因此，尋求一種高效率、低成本、易操作的耐火磚損傷檢測方法迫在眉睫。

頻率不僅可以反映結(jié)構(gòu)的健康狀況，而且是結(jié)構(gòu)動力學(xué)響應(yīng)中最易測得且測量最準的特征量。利用頻率指紋進行損傷檢測具有高效性、準確性和經(jīng)濟性等優(yōu)點，因此頻率指紋被廣泛應(yīng)用于橋梁、混凝土、木材、農(nóng)作物等方面的損傷檢測中[3]，在耐火磚的損傷檢測領(lǐng)域也已有相關(guān)應(yīng)用：文獻[4]建立了一套耐火磚內(nèi)部質(zhì)量檢測系統(tǒng)，通過分析敲擊聲信號的一階固有頻率，可以自動識別受損的耐火磚；文獻[5]利用模糊聚類的方法分析耐火磚的共振頻率，識別了受損耐火磚并確定了其損傷類型。

但上述方法僅判斷了耐火磚是否受損，而沒有進一步的損傷定位或損傷定量的研究。耐火磚的模態(tài)頻率不僅與損傷程度有關(guān)，還受損傷位置的影響。文獻[6]在對結(jié)構(gòu)進行模態(tài)分析時得出結(jié)論：任意兩階頻率變化量之比僅是損傷位置的函數(shù)而與損傷程度無關(guān)。文獻[7]提出任意兩階模態(tài)頻率變化量的平方比僅是結(jié)構(gòu)損傷位置的函數(shù)。故只有確定了耐火磚的損傷位置，才能進一步確定其損傷程度。國內(nèi)外利用模態(tài)頻率進行損傷定位的研究眾多。文獻[8]提出了利用頻率變化量之比和模態(tài)應(yīng)變能變化率進行梁類結(jié)構(gòu)的損傷識別：先根據(jù)頻率變化量之比確定損傷的可疑位置；再通過損傷前后模態(tài)應(yīng)變能的變化率準確定位損傷位置。文獻[9]利用頻率變化量的平方比作為損傷的動力指紋，判定結(jié)構(gòu)的損傷位置，并采用插值法從頻率改變率中獲取剛度損失量進行損傷定量。

根據(jù)結(jié)構(gòu)的動力學(xué)響應(yīng)特性可知，高階模態(tài)頻率的能占比較低，實際測試中高階模態(tài)頻率不易測得。通過實驗可較好測得耐火磚的前三階模態(tài)頻率，根據(jù)耐火磚損傷后第三階與第一階頻率變化量的比值以及第三階與第二階頻率變化量的比值確定耐火磚的損傷位置，并根據(jù)該位置頻率變化率與損傷系數(shù)的關(guān)系曲線圖，確定耐火磚的損傷程度。

2 損傷位置的確定

2.1 理論分析

任意結(jié)構(gòu)的振動特性都可以看作是由質(zhì)量、阻尼以及剛度矩陣構(gòu)成的振動系統(tǒng)。若忽略阻尼的影響，結(jié)構(gòu)的振動方程為：

式中：K—結(jié)構(gòu)剛度矩陣；M—結(jié)構(gòu)質(zhì)量矩陣；ω—固有頻率；φ—振型。

發(fā)生損傷后，結(jié)構(gòu)的質(zhì)量和剛度會改變，從而結(jié)構(gòu)的固有頻率和振型也會發(fā)生改變，由攝動理論式（1）可知：

式中：ΔK—剛度矩陣變化量；ΔM—質(zhì)量矩陣變化量；Δω—頻率變化量；Δφ—振型變化量。耐火磚的內(nèi)部缺陷導(dǎo)致結(jié)構(gòu)剛度明顯降低，而質(zhì)量變化并不明顯。若不考慮質(zhì)量變化，則ΔM=0；式（2）則變?yōu)椋?/p>

展開式（3），且忽略二階項，整理可得：

單元變形可以由結(jié)構(gòu)振型推導(dǎo)出來，對于某階模態(tài)i（i=1，2，3…）則有：

式中：A—結(jié)構(gòu)單元總數(shù)；m—單元號；ε—單元變形。若結(jié)構(gòu)只存在單個損傷，即損傷單元只有一個，設(shè)損傷單元為N，則：

由式（6）可以看出，在不同階數(shù)的模態(tài)頻率下，單元的變形是變化的，而單元的剛度矩陣變化量ΔK 是不變的，因此定義單元損傷系數(shù)[9]：

式中：αN—損傷系數(shù)矩陣；KN—單元損傷前的剛度—單元損傷后的剛度，則：

將式（8）帶入式（6）得：

同理，對于某階模態(tài)j（j=1，2，3…）有：

兩式相比得到損傷后任意兩階頻率變化量之比：

由式（11）可知，結(jié)構(gòu)損傷后任意兩階頻率變化量的比值僅與損傷位置有關(guān)，與損傷程度無關(guān)。

2.2 數(shù)值仿真

選取的耐火磚主要參數(shù)為：彈性模量E=21542MPa，泊松比σ=-0.53，密度 ρ=2162kg·m-3，b×h×l=115mm×60mm×230mm。耐火磚的形狀是對稱的，若采用對稱的約束形式，會推導(dǎo)出兩個損傷位置，因此將耐火磚以懸臂梁的形式進行約束。分別在耐火磚的10mm、20mm、30mm…220mm 處（0mm 處為耐火磚的自由端）制造深15mm、寬0.8mm 的貫穿性缺口模擬損傷，利用有限元數(shù)值仿真軟件對耐火磚進行模態(tài)分析，得到不同損傷位置Δω3與Δω1、Δω3與Δω2的比值。以耐火磚的損傷位置為橫坐標（橫坐標為歸一化坐標，耐火磚的自由端為橫坐標的零點），以Δω3與Δω1的比值為縱坐標繪制曲線圖，如圖1 所示。

圖1 Δω3/Δω1 隨損傷位置變化的曲線圖Fig.1 Curve Map of Δω3/Δω1 with Position Change

同樣，以耐火磚的損傷位置為橫坐標，以Δω3與Δω2的比值為縱坐標繪制曲線圖，如圖2 所示。

圖2 Δω3/Δω2 隨損傷位置變化的曲線圖Fig.2 Curve Map of Δω3/Δω2 with Position Change

通過實驗測得結(jié)構(gòu)完好時前三階頻率ω1、ω2、ω3的值，以及結(jié)構(gòu)損傷后前三階頻率的值。算出 Δω3/Δω1的值，在圖1 中作一條過點（0，Δω3/Δω1）的水平線，與曲線交于點 A1（ld1，Δω3/Δω1），該點的橫坐標 ld1即為損傷位置。或算出 Δω3/Δω2的值，在圖 2 中作一條過點（0，Δω3/Δω2）的水平線，與曲線交于點A2（ld2，Δω3/Δω2），該點的橫坐標 ld2即為損傷位置。經(jīng)仿真和實驗發(fā)現(xiàn)，當損傷位于耐火磚（0～50）mm 位置時，無法準確推導(dǎo)出損傷的位置，由于耐火磚結(jié)構(gòu)的對稱性，只需置換耐火磚的約束端與自由端，損傷位置則變換到（180～230）mm 位置之間，從而確定損傷位置。

3 損傷程度的確定

文獻[10]在利用共振頻率測試耐火磚彈性模量的研究中，推導(dǎo)了耐火磚的一階共振頻率ω、截面慣性矩J 以及彈性模量E 之間的關(guān)系：

耐火磚的截面慣性矩J=bh3/12，某處受損后，損傷處橫截面的高hd=h-h0，截面慣性矩耐火磚某處存在損傷時，此處的截面慣性矩減小，進而導(dǎo)致結(jié)構(gòu)整體的截面慣性矩減小。隨著損傷程度的加劇，結(jié)構(gòu)的整體截面慣性矩J 會隨之減小，結(jié)合式（12）可知，由于彈性模量E 為定值，頻率ω 也會隨之降低。

沿x 軸方向?qū)⒛突鸫u平均劃分成23 個單元，如圖3 所示。其中耐火磚的右端為自由端。

圖3 單元劃分圖Fig.3 Unit Division

在第一個單元的中間位置制造寬為0.8 mm、深為h0的貫穿性缺口模擬損傷，通過有限元數(shù)值仿真軟件對耐火磚的損傷進行分析，計算出未受損時耐火磚的一階模態(tài)頻率ω1以及不同損傷深度時耐火磚的一階模態(tài)頻率引用損傷系數(shù)α=h0/h 描述損傷程度，以損傷系數(shù)α 為橫坐標，頻率變化率為縱坐標，根據(jù)分析結(jié)果繪制出單元1 的頻率變化率與損傷系數(shù)的關(guān)系曲線圖。利用相同的方法分別繪制出第2-23 個單元頻率變化率與損傷系數(shù)的關(guān)系曲線圖，曲線圖的變化趨勢基本相同，其中單元15（（150～160）mm 位置）的頻率變化率與損傷系數(shù)的關(guān)系曲線圖，如圖4 所示。

圖4 單元15 的頻率變化率-損傷系數(shù)變化曲圖Fig.4 Curve of Frequency Variation with Damage at Fifteenth Units

4 實驗驗證

在耐火磚的155 mm 處（位于第15 個單元），制造一個寬0.8 mm、深20 mm 的貫穿性缺口，使用敲擊錘敲擊耐火磚，接收并存儲敲擊產(chǎn)生的聲信號。利用MATLAB 軟件對聲信號進行FFT 變換，提取耐火磚的前三階模態(tài)頻率。分析結(jié)果，如表1 所示。

表1 分析結(jié)果Tab.1 Analysis Results

在圖1 中作一條過點（0，2.33）的水平線，與曲線交于點A1，如圖5 所示。在圖2 中作一條過點（0，3.75）的水平線，與曲線交于點A2，如圖6 所示。由圖5 可知A1的橫坐標在0.67 附近，對應(yīng)的耐火磚損傷位置ld1為154.1 mm。由圖6 可知A2的橫坐標在0.66 附近，對應(yīng)的耐火磚損傷位置ld2為151.8 mm。根據(jù)ld1和ld2可以確定損傷位于第15 個單元。在圖 4 中作一條過點（0，14.5）的水平線，交曲線于點D，如圖7 所示。點D 對應(yīng)的橫坐標誤差為故此方法可以較為準確的推斷耐火磚的損傷位置和損傷程度。

圖5 Δω3/Δω1 隨損傷位置變化的曲線圖Fig.5 Curve Map of Δω3/Δω1 with Position Change

圖6 Δω3/Δω2 隨損傷位置變化的曲線圖Fig.6 Curve Map of Δω3/Δω2 with Position Change

圖7 單元15 的頻率變化率-損傷系數(shù)變化曲線圖Fig.7 Curve of Frequency Variation with Damage at Fifteenth Units

5 結(jié)論

首先從理論上推導(dǎo)了損傷導(dǎo)致的任意兩階頻率變化量之比僅是損傷位置的函數(shù)而與損傷程度無關(guān)，依此繪制出耐火磚頻率變化量的比值與損傷位置的關(guān)系曲線圖，根據(jù)頻率變化量的比值確定了耐火磚的損傷位置；再繪制出該位置第一階頻率變化率與損傷系數(shù)的關(guān)系曲線圖，通過頻率變化率確定了耐火磚的損傷程度。并和實驗結(jié)果做了對比，結(jié)果表明：（1）通過耐火磚損傷后的第一階與第二階頻率變化量的比值可以確定損傷位置。（2）確定損傷位置后，根據(jù)該位置第一階頻率變化率與損傷系數(shù)的關(guān)系曲線圖可以確定損傷程度；（3）對位于耐火磚（0～50）mm 位置的損傷，直接采用該方法不能獲取損傷位置，須置換耐火磚的約束端與自由端，此時損傷位于耐火磚（180～230）mm 位置，即可確定損傷位置。