茅志穎，陳國(guó)平，張保強(qiáng)

(南京航空航天大學(xué) 航空宇航學(xué)院，南京 210016)

結(jié)構(gòu)在使用過(guò)程中，由于受外部環(huán)境等的影響容易造成損傷，會(huì)影響到結(jié)構(gòu)的力學(xué)性能，因此結(jié)構(gòu)的損傷定位成為近年來(lái)的一個(gè)熱點(diǎn)問(wèn)題。由于容易通過(guò)實(shí)驗(yàn)獲得結(jié)構(gòu)模態(tài)參數(shù)，因此，以結(jié)構(gòu)振動(dòng)性能為基礎(chǔ)的結(jié)構(gòu)損傷定位方法被廣泛使用［1－5］。任何結(jié)構(gòu)都可以被考慮成一個(gè)帶有剛度、質(zhì)量和阻尼的振動(dòng)系統(tǒng)。一旦損傷在結(jié)構(gòu)中出現(xiàn)，結(jié)構(gòu)振動(dòng)的基本參數(shù)會(huì)被改變。結(jié)構(gòu)模態(tài)參數(shù)的變化可以作為早期損傷出現(xiàn)在結(jié)構(gòu)中的信號(hào)。基于結(jié)構(gòu)振動(dòng)性能的損傷定位主要是用結(jié)構(gòu)的模態(tài)參數(shù)如固有頻率，模態(tài)阻尼和模態(tài)振型等來(lái)實(shí)現(xiàn)［6－11］。

拓?fù)鋬?yōu)化在結(jié)構(gòu)設(shè)計(jì)優(yōu)化領(lǐng)域有著廣泛的應(yīng)用，但是在結(jié)構(gòu)損傷定位上的應(yīng)用還不多見［12－13］，王云，郝際平［14］研究了金屬材料內(nèi)部微觀缺陷在幾何方程上的體現(xiàn)。本文中，利用漸進(jìn)結(jié)構(gòu)優(yōu)化法結(jié)合結(jié)構(gòu)模態(tài)參數(shù)研究了其在微小損傷結(jié)構(gòu)損傷定位應(yīng)用上的可行性。這種方法基于拓?fù)鋬?yōu)化的特質(zhì)，可以在結(jié)構(gòu)整個(gè)區(qū)域?qū)ふ覔p傷點(diǎn)，具有一定的應(yīng)用前景。

1 損傷定位的數(shù)學(xué)模型

漸進(jìn)結(jié)構(gòu)優(yōu)化法的原理是通過(guò)選取適當(dāng)?shù)臏?zhǔn)則來(lái)判定結(jié)構(gòu)中每一個(gè)單元對(duì)結(jié)構(gòu)某特性的重要程度，也可以說(shuō)是使得結(jié)構(gòu)的某一特性達(dá)到某一程度，結(jié)構(gòu)中每一個(gè)單元的貢獻(xiàn)大小。在本文中，即是指通過(guò)判斷所有單元對(duì)于結(jié)構(gòu)的幾階共振頻率和反共振頻率的貢獻(xiàn)的大小來(lái)確定結(jié)構(gòu)中最有可能損傷的區(qū)域。

損傷定位的數(shù)學(xué)模型可以表示為如下形式:

式中Q=［q1q2… qn］T為結(jié)構(gòu)單元拓?fù)鋬?yōu)化設(shè)計(jì)向量;qi∈{0 1}，分別表示單元在損傷區(qū)域和不在損傷區(qū)域兩種情況;mi為單元i的質(zhì)量。其中，ωkr和 ωka為原結(jié)構(gòu)的共振頻率和反共振頻率，和為損傷結(jié)構(gòu)的共振頻率和反共振頻率，Nr為共振頻率的約束個(gè)數(shù)，Na為反共振頻率的約束個(gè)數(shù)，ε為容許的誤差值。

結(jié)構(gòu)振動(dòng)特性由其特征方程給出，在對(duì)結(jié)構(gòu)進(jìn)行動(dòng)力學(xué)分析時(shí)，無(wú)阻尼結(jié)構(gòu)的振動(dòng)特征方程為:

式中:［K］為結(jié)構(gòu)總剛度矩陣，為N階方陣;［M］為結(jié)構(gòu)總質(zhì)量矩陣，為N階方陣;ωi為第i階固有頻率且λi=ω2i;{φi}為第i階固有振型

經(jīng)過(guò)推導(dǎo)可得結(jié)構(gòu)的i階頻率ωi的靈敏度αij為:

其中［ΔKj］和［ΔMj］為刪除第j個(gè)單元使結(jié)構(gòu)的質(zhì)量和剛度矩陣的變化。

對(duì)于一般的無(wú)阻尼結(jié)構(gòu)，其系統(tǒng)的頻響函數(shù)矩陣為:

Hts(ω)表示在s自由度激勵(lì)，在t自由度測(cè)量所獲得的頻響。其中，Hits(ω)的分子是伴隨矩陣 adj(［K］－ω2［M］)的第 ts項(xiàng) H*its(ω)，表示如下:

其中，［Kt，s］和［Mt，s］分別表示刪除第 t行和第 s列元素后的質(zhì)量矩陣和剛度矩陣。當(dāng)t=s時(shí)，稱為原點(diǎn)頻響函數(shù)，當(dāng)t≠s時(shí)，稱為跨點(diǎn)頻響函數(shù)。

結(jié)構(gòu)的反共振頻率可以通過(guò)下式的特征值問(wèn)題獲得:

其中，λa即為反共振頻率所對(duì)應(yīng)的特征值。

當(dāng)Hts(ω)為原點(diǎn)頻響函數(shù)時(shí)，共振頻率和反共振頻率交叉成對(duì)出現(xiàn)，比結(jié)構(gòu)固有頻率少一個(gè)，相當(dāng)于將原系統(tǒng)第i個(gè)自由度固定獲得的新結(jié)構(gòu)系統(tǒng)的固有頻率。而當(dāng)Hts(ω)為跨點(diǎn)頻響函數(shù)時(shí)，不存在共振頻率和反共振頻率交叉的特性。結(jié)構(gòu)各原點(diǎn)頻響函數(shù)的反共振頻率一般是不相同的，如圖1所示。共振現(xiàn)象是系統(tǒng)的全局現(xiàn)象，而反共振現(xiàn)象是系統(tǒng)的一種局部現(xiàn)象。

對(duì)于原點(diǎn)反共振頻率，結(jié)構(gòu)中刪除單元的特征值變化可以表示為:

反共振頻率特征值λai的變化值為Δjλai。

對(duì)于跨點(diǎn)反共振頻率，結(jié)構(gòu)中刪除單元的特征值的變化量可以表示為:

其中，{φail}和{φair}分別為左右特征向量。

由于，原點(diǎn)反共振頻率的特殊性，在本文中，所使用的反共振頻率都是指某點(diǎn)的原點(diǎn)反共振頻率。

圖1 不同點(diǎn)的原點(diǎn)頻響函數(shù)Fig.1 The P-FRFs of different measurement points

2 基于漸進(jìn)結(jié)構(gòu)優(yōu)化法的損傷定位策略

當(dāng)損傷很小時(shí)，搜尋的過(guò)程會(huì)比較困難。在拓?fù)鋬?yōu)化的過(guò)程中，因?yàn)樵O(shè)計(jì)范圍是結(jié)構(gòu)所有的有限元單元，所以當(dāng)單元數(shù)目較多時(shí)，尋找微小的損傷單元就會(huì)需要很多的工作量。因此在本文中，尋找微小的損傷單元時(shí)，采用一種判斷方法來(lái)逐步排除不可能的單元，同時(shí)向損傷區(qū)域逼近。這種方法可以直觀的表示為如圖2所示。

結(jié)構(gòu)的單元是否處于損傷區(qū)域的單元，是一個(gè)有可能是或者有可能非的模糊概念，因此，結(jié)合模糊數(shù)學(xué)，給每一個(gè)單元一個(gè)可能是或者可能非的程度大小一個(gè)數(shù)值衡量。在本文中采用線性隸屬函數(shù)的中間型表示結(jié)構(gòu)中每一個(gè)單元對(duì)于結(jié)構(gòu)共振頻率和反共振頻率變化的靈敏度。

其中，a，b，c，d的計(jì)算公式可以表示如下:

圖2 損傷區(qū)域逼近Fig.2 The approach strategy of damage region

當(dāng)結(jié)構(gòu)的損傷比較微小時(shí)，也即此時(shí)結(jié)構(gòu)的幾階共振頻率和反共振頻率都沒有太大的變化，此時(shí)，采用區(qū)域逼近方法的搜索過(guò)程可以表示為如圖3所示。

(1)將無(wú)損傷結(jié)構(gòu)劃分有限元單元，確定初始搜尋范圍內(nèi)的單元集合Ne;

(2)確定無(wú)損傷結(jié)構(gòu)的初始共振頻率和某點(diǎn)的原點(diǎn)反共振頻率以及損傷結(jié)構(gòu)的共振頻率和相同點(diǎn)的原點(diǎn)反共振頻率;

(3)計(jì)算無(wú)損傷結(jié)構(gòu)的單元對(duì)于某一頻率或反共振頻率的靈敏度;

(4)利用式(8)對(duì)每個(gè)單元的靈敏度求取隸屬度uiω;

圖3 微小損傷定位策略Fig.3 The strategy of tiny damage localization

(5)判斷每個(gè)單元的 uω是否小于 δ，同時(shí)更新Ne;

(6)對(duì)Ne的單元，確定可能的損傷單元;

(7)判斷是否收斂，如果收斂，則損傷搜尋停止;反之，則重復(fù)(3)～(6)步，直到搜尋到最可能的損傷區(qū)域。

3 數(shù)值算例

算例1，如圖4和圖5所示分別為一無(wú)損傷結(jié)構(gòu)和一微小的凹槽型損傷結(jié)構(gòu)，左端固支。

原結(jié)構(gòu)的共振頻率和A點(diǎn)的原點(diǎn)反共振頻率以及損傷結(jié)構(gòu)的共振頻率和A點(diǎn)的原點(diǎn)反共振頻率分別如表1所示。

表1 損傷結(jié)構(gòu)和無(wú)損傷結(jié)構(gòu)的模態(tài)參數(shù)Tab.1 Modal data of the undamaged and damaged structure

圖6 A點(diǎn)的原點(diǎn)頻響函數(shù)Fig.6 The comparison of the P-FRFs for point A

A點(diǎn)的原點(diǎn)頻響函數(shù)曲線如圖6所示。由圖6可以發(fā)現(xiàn)，結(jié)構(gòu)的共振頻率和反共振頻率都發(fā)生了較小的變化。取δ值為0.2，結(jié)構(gòu)的損傷定位使用了損傷區(qū)域逼近的方法，損傷結(jié)構(gòu)的搜尋過(guò)程如圖7所示。

圖7 損傷定位過(guò)程Fig.7 The process of the damage localization

結(jié)構(gòu)的損傷定位過(guò)程如圖7所示，在搜尋過(guò)程中，由于使用了損傷區(qū)域逼近方法，隨著搜索的進(jìn)行，其可能得損傷區(qū)域在逐步縮小范圍，需要判斷的單元也在逐漸減小，最后只剩下14個(gè)可能的損傷單元。從最后的定位結(jié)果看，采用區(qū)域逼近的方法，其損傷定位還是比較準(zhǔn)確的。

算例2，如圖8和圖9所示，分別為一無(wú)損傷結(jié)構(gòu)和一具有多處微小損傷結(jié)構(gòu)，左端固支。

原結(jié)構(gòu)的共振頻率和C、D點(diǎn)的原點(diǎn)反共振頻率以及損傷結(jié)構(gòu)的共振頻率和C、D點(diǎn)的原點(diǎn)反共振頻率分別如表2所示。

圖8 無(wú)損傷結(jié)構(gòu)Fig.8 Undamaged structure

圖9 損傷結(jié)構(gòu)Fig.9 Damaged structure

C點(diǎn)和D點(diǎn)的原點(diǎn)頻響函數(shù)曲線分別如圖10和圖11所示。由圖10和11可以發(fā)現(xiàn)，結(jié)構(gòu)的各階頻率和反共振頻率并沒有太大的變化，判斷是微小損傷，使用損傷區(qū)域逼近的方法來(lái)搜索損傷區(qū)域，取δ值為0.2，其搜尋過(guò)程如圖12所示。

表2 損傷結(jié)構(gòu)和無(wú)損傷結(jié)構(gòu)的模態(tài)參數(shù)Tab.2 Modal data of the undamaged and damaged structure

圖10 C點(diǎn)的原點(diǎn)頻響函數(shù)Fig.10 The comparison of the P-FRFs for point C

在本算例中，為了識(shí)別多處微小損傷，因此，使用了兩個(gè)點(diǎn)的反共振頻率數(shù)據(jù)，得到的定位結(jié)果很準(zhǔn)確。相比于算例1只使用了一個(gè)點(diǎn)的反共振頻率就能獲得比較準(zhǔn)確的定位結(jié)果，在本算例中，當(dāng)只使用C點(diǎn)的反共振頻率時(shí)，得到的損傷定位結(jié)果如圖13所示。

圖11 D點(diǎn)的原點(diǎn)頻響函數(shù)Fig.11 The comparison of the P－ FRFs for point D

圖12 損傷定位過(guò)程Fig.12 The process of the damage localization

圖13所示的是只使用了共振頻率和C點(diǎn)一個(gè)點(diǎn)的反共振頻率的定位結(jié)果，從結(jié)果看，可以發(fā)現(xiàn)，此時(shí)得到的定位結(jié)果并不是很準(zhǔn)確，由此可以說(shuō)明，為了得到盡可能準(zhǔn)確的結(jié)果，足夠的模態(tài)數(shù)據(jù)，是非常必要的。

圖13 沒有使用D點(diǎn)反共振頻率的定位結(jié)果Fig.13 The result without the anti-resonances of the point D

4 結(jié)論

本文中以結(jié)構(gòu)模態(tài)參數(shù)主要是結(jié)構(gòu)共振頻率和結(jié)構(gòu)某點(diǎn)的反共振頻率為基礎(chǔ)，結(jié)合漸進(jìn)結(jié)構(gòu)優(yōu)化法，討論了其在微小損傷結(jié)構(gòu)的損傷定位中應(yīng)用的可行性。文章提出了針對(duì)微小損傷結(jié)構(gòu)的損傷區(qū)域逼近方法，并提出了合適的結(jié)合損傷區(qū)域逼近的定位策略。針對(duì)本文所提的微小損傷結(jié)構(gòu)，給出了的算例分別對(duì)單一的微小損傷和多處微小損傷的結(jié)構(gòu)進(jìn)行討論。算例表明，漸進(jìn)結(jié)構(gòu)優(yōu)化法在已經(jīng)測(cè)得結(jié)構(gòu)損傷前后共振頻率和反共振頻率的基礎(chǔ)上，在結(jié)構(gòu)整個(gè)區(qū)域?qū)ふ覔p傷點(diǎn)，即使對(duì)于微小的損傷結(jié)構(gòu)也可以有很好的搜尋效果，可以在工程上推廣應(yīng)用。

［1］ Doebling S W，F(xiàn)arrar C R，Prime M B，et al.Damage identification and health monitoring of structural and mechanical systems from changes in their vibration characteristics:a literature review［J］.Technical Report，Research Report LA－13070－MS，Los Alamos National Laboratory，1996.

［2］ Zou Y，Tong L，Steven G P.Vibration-based modeldependentdamagedelamination identification and health monitoring for composite structures-a review［J］.Journal of Sound and Vibration，2000，230(2):357－378.

［3］ Doebling S W，F(xiàn)arrar C R，Prime M B.A summary review of vibration-based damage identification methods［R］.Technical Report，Los Alamos National Laboratory，1998.

［4］ Sohn H，F(xiàn)arrar C，Hemez F，et al.A review of structural health monitoring literature［R］.Technical Report，Report LA-13976-MS，Los Alamos National Laboratory，2001.

［5］ Carden E，F(xiàn)anning P.Vibration based condition monitoring:a review［J］.Structural Health Monitoring，2004，3(4):355－377.

［6］ Shahdin A，Morlier J，Gourinat Y.Correlating low energy impactdamagewith changesin modalparameters:a preliminary study on composite beams［J］.Structural Health Monitoring，2009，8(6):523－536.

［7］ Shahdin A，Morlier J，Gourinat Y.Significance of low energy impact damage on modal parameters of composite beams by design of experiments［J］.Seventh International Conference on Modern Practice in Stress and Vibration Analysis，MPSVA 2009，2009.

［8］ Tracy J，Pardoen G.Effect of delamination on the natural frequencies of composite laminates［J］.Journal of Composite Materials，1989，23(12):1200－1215.

［9］ Adams R D，Cawley R.The localization of defects in structures from measurements of natural frequencies［J］.Journal of Strain Analysis，1979，14:49－57.

［10］ Richardson M H，Mannan M A.Correlating minute structural faults with changes in modal parameters［J］.Proceedings of SPIE，International Society of Optical Engineering，1993，1923(2):893－898.

［11］ Zhang Z，Hartwig G.Relation of damping and fatigue damage of unidirectional fibre composites［J］.International Journal of Fatigue，2004，24:713－738.

［12］ Lee J S，Kim J E，Kim Y Y.Damage detection by the topology design formulation using modal parameters［J］.International Journal for Numerical Methods in Engineering，2007，69:1480－1498.

［13］ Niemann H，Morlier J，Shahdin A，et al.Damage localization using experimental modal parameters and topology optimization［J］.Mechanical Systems and Signal Processing，2010，24:636－652.

［14］ 王 云，郝際平.含微缺陷金屬材料損傷理論的幾何拓?fù)洌跩］.建筑鋼結(jié)構(gòu)進(jìn)展，2009，11(1):49－53.