馮彥軍 石川千 周宇 趙發(fā)剛 周徐斌

(上海衛(wèi)星工程研究所，上海 201109)

基于頻響主元的衛(wèi)星結(jié)構(gòu)損傷定位與評估方法

馮彥軍 石川千 周宇 趙發(fā)剛 周徐斌

(上海衛(wèi)星工程研究所，上海 201109)

鑒于目前國內(nèi)衛(wèi)星預(yù)、復(fù)振頻響結(jié)果的評估主要依靠設(shè)計師的經(jīng)驗分析，缺乏高效準確的量化手段，文章基于頻響函數(shù)(FRF)的主元分析(PCA)，提出了兩種針對衛(wèi)星結(jié)構(gòu)的損傷定位指標，分別是多主元綜合指標和FRF主元殘差均方根；并采用平方預(yù)測誤差(SPE)指數(shù)對衛(wèi)星結(jié)構(gòu)的整體損傷程度進行評估。以衛(wèi)星主承力部件中心承力筒為分析對象，驗證了基于衛(wèi)星結(jié)構(gòu)頻響主元的損傷定位和評估方法的有效性。文章對基于預(yù)、復(fù)振對比分析的損傷識別、量化試驗評估標準、提高衛(wèi)星研制效率有促進作用。

損傷識別；多主元綜合指標；頻響函數(shù)主元殘差均方根；平方預(yù)測誤差指數(shù)；衛(wèi)星預(yù)、復(fù)振；頻響函數(shù)

1 引言

動力學環(huán)境試驗是對衛(wèi)星等航天產(chǎn)品設(shè)計和工藝狀態(tài)適用性的考核，是研制過程中不可缺少的環(huán)節(jié)[1]。在每個滿量級的動力學環(huán)境試驗前后，應(yīng)分別進行預(yù)、復(fù)振試驗檢驗，通過對比預(yù)、復(fù)振頻響曲線獲得傳遞特性、共振頻率和各階幅值放大系數(shù)等結(jié)構(gòu)特性，識別結(jié)構(gòu)在振動試驗過程中可能發(fā)生的結(jié)構(gòu)損傷[2]。大量的振動試驗結(jié)果表明：預(yù)、復(fù)振試驗響應(yīng)曲線不一致現(xiàn)象非常普遍，典型現(xiàn)象有共振峰漂移、共振峰數(shù)量變化和共振峰幅值變化等，究其機理，主要受邊界非線性、材料非線性和結(jié)構(gòu)損傷等因素的影響[3-9]。目前對于預(yù)、復(fù)振試驗頻響結(jié)果不一致現(xiàn)象的分析，主要依靠設(shè)計師的工程經(jīng)驗，檢查手段局限于對試驗件的拆解和復(fù)查故障，準確性不高，效率低下，大大制約了衛(wèi)星的研制進度。因此，基于衛(wèi)星預(yù)、復(fù)振頻響數(shù)據(jù)的損傷識別是極具應(yīng)用價值并亟待解決的問題。

將損傷識別技術(shù)應(yīng)用于衛(wèi)星預(yù)、復(fù)振不一致問題的分析，可以有效地利用試驗數(shù)據(jù)信息，并找到故障所在。基于結(jié)構(gòu)動力特性的損傷識別是近幾十年國內(nèi)研究的熱點。目前大多數(shù)損傷識別方法是通過頻響函數(shù)(Frequency Response Function,FRF)數(shù)據(jù)進行模態(tài)擬合，提取模態(tài)參數(shù)作為損傷標識量，但存在測量噪聲、模態(tài)提取誤差、數(shù)據(jù)冗余和分析難度大等問題，因此一些學者直接利用FRF中豐富的信息進行損傷識別[10-12]，或用主元分析(Principal Component Analysis,PCA)方法對FRF數(shù)據(jù)進行壓縮、提取和處理[13-16]，得到包含損傷信息的特征數(shù)據(jù)，以此構(gòu)造損傷識別指標。其中，損傷定位指標的優(yōu)劣嚴重影響損傷識別的精度：為了避免構(gòu)造不適當?shù)膿p傷定位指標，一些學者[14-16]將PCA壓縮的FRF數(shù)據(jù)直接或間接輸入到神經(jīng)網(wǎng)絡(luò)，進行損傷模式識別，但需要大量有效的實測數(shù)據(jù)來訓練神經(jīng)網(wǎng)絡(luò)，而且識別精度受測量噪聲和模型誤差的影響很大。楊彥芳等人[17]采用橢圓控制圖和T2控制圖實現(xiàn)對網(wǎng)架結(jié)構(gòu)的損傷定位，但網(wǎng)架損傷形式單一，且不能進行多損傷識別和損傷程度評估。相比于FRF，F(xiàn)RF曲率對損傷更敏感，Maia[18]、姜增國[19]和朱新圓[20]等人利用頻響函數(shù)曲率及其變形構(gòu)造損傷識別指標，對傳力形式相對單一的梁、板和桁架結(jié)構(gòu)進行了損傷識別。但現(xiàn)有的基于FRF的損傷識別理論，僅應(yīng)用于梁桿件、平板、桁架等傳力形式相對單一的特殊結(jié)構(gòu)，而衛(wèi)星結(jié)構(gòu)形式多樣，傳力路徑復(fù)雜，主承力結(jié)構(gòu)以復(fù)合材料殼為主，單純利用以上現(xiàn)有方法不能很好地識別局部損傷位置，更不能對損傷程度進行評估。

針對上述問題，本文提出了兩種線性歸一化的無量綱化指標：多主元綜合指標和FRF主元殘差均方根，分別對衛(wèi)星結(jié)構(gòu)FRF數(shù)據(jù)進行分析，實現(xiàn)結(jié)構(gòu)損傷定位，并用平方預(yù)測誤差(SquarePredictionError，SPE)指數(shù)對結(jié)構(gòu)整體損傷程度進行了評估。通過衛(wèi)星主承力筒模型在不同損傷工況下的數(shù)值仿真分析，對上述損傷定位和程度評估方法進行了驗證。

2 基于頻響主元分析的損傷識別理論

2.1 PCA方法

PCA分析步驟為：

(1)對H(ω)按列進行標準化處理，將各變量化為均值為0，方差為1的標準變量。

(2)對H(ω)的協(xié)方差矩陣S進行特征值分解：

(1)

2.2 損傷定位理論

首先對損傷識別矩陣H(ω)進行主元分析，得到包含損傷信息的主元Fn×k，然后基于其主元構(gòu)造出兩種損傷定位指標：多主元綜合指標和FRF主元殘差均方根。

2.2.1 多主元綜合指標

由多元控制理論知，當測點數(shù)n足夠大時，由任意測點j前2階主元決定的坐標(yj1,yj2)服從式(2)的橢圓控制域分布。第j點前2階主成分對應(yīng)的點距離橢圓控制邊界的距離，一定程度上代表了該點數(shù)據(jù)的異常程度，因此定義偏離橢圓邊界的相對距離為dj1，見式(3)。

(2)

(3)

(4)

(5)

則第j點的總體相對距離為dj=dj1+dj2。各測點的dj構(gòu)成向量d，經(jīng)0～1線性歸一化，即得到多主元綜合指標：

(6)

控制限(ControlLimit,CL)的確定。橢圓控制限和T2控制限對應(yīng)的相對距離dCL=0，將dCL與d一起歸一化即可得到標準化CL，見式(7)，式中L1為CL的值。d中大于L1的元素對應(yīng)的測點被認為是異常點，即位于損傷區(qū)域附近的測點。

(7)

2.2.2FRF主元殘差均方根

對損傷識別矩陣H(ω)的前k階主成分矩陣Fn×k，定義第j個測點對應(yīng)的FRF主元殘差均方根為rj(見式(8))，各測點構(gòu)成的FRF主元殘差均方根向量為R，經(jīng)0～1線性歸一化后，可得到標準化的FRF主元殘差均方根向量R，見式(9)。R中大于CL值的元素對應(yīng)的測點被認為是異常點。

(8)

(9)

CL的確定，可以通過兩種方法確定：一種是通過統(tǒng)計的方法，統(tǒng)計學中常用3σ準則來排除數(shù)據(jù)中偏差較大的數(shù)據(jù)，本文中用L2=E(R)+2D(R)確定CL值，E(R)為R的均值，D(R)為R的標準差；另一種方法，考慮到具體實際試驗中由于噪聲、激勵能量不穩(wěn)定及其它干擾因素，即使對于同樣工況的同一試件，其實測頻響數(shù)據(jù)的主元殘差均方根一般不為零。因此，以完好結(jié)構(gòu)相同條件下的連續(xù)兩次實測頻響數(shù)據(jù)的主元殘差均方根向量中最大元素為CL。若某工況下某測點對應(yīng)的FRF主元殘差均方根大于CL，則被認為是異常點。

2.3 損傷程度評估理論

(10)

式中：I為單位矩陣；SPE指數(shù)定義為Hd(ω)在RS投影的方差，表征了有損結(jié)構(gòu)實測數(shù)據(jù)對無損結(jié)構(gòu)主元模型的偏離程度，見式(11)，E為SPE指數(shù)的值。因此，當結(jié)構(gòu)無損傷(或損傷程度較小)時，Hd(ω)的主要信息落在主元子空間(PrincipleComponentSubspace,PCS)，SPE指數(shù)很小，主要為噪聲；當結(jié)構(gòu)有損傷(或損傷程度較大)時，Hd(ω)在RS的投影就會增加，SPE指數(shù)也增加，主要為結(jié)構(gòu)損傷引起的異常信息。因而可通過SPE指數(shù)來評估結(jié)構(gòu)的整體損傷程度。

(11)

3 衛(wèi)星承力筒的頻響主元分析算例

3.1 算例模型

為了驗證上述提出的損傷識別方法，采用衛(wèi)星典型主承力部件承力筒進行數(shù)值驗證。如圖1所示，承力筒筒體柱段外徑Ф1114mm，內(nèi)徑Ф1090mm，高3260mm，主體采用碳纖維/環(huán)氧-鋁蜂窩夾層結(jié)構(gòu)復(fù)合材料，橫向加強框和縱向桁條均為碳纖維材料。采用Pro/E軟件進行三維實體建模，Patran軟件進行有限元建模。筒體部分采用殼單元，利用三明治夾芯等效理論模擬等效蜂窩夾層復(fù)合材料；面板安裝法蘭和縱向桁條采用梁單元。筒體部分單元總數(shù)為10 000，周向和縱向各劃分100個網(wǎng)格節(jié)點，編號從星箭連接環(huán)端至上端加強框依次增大。

圖1 衛(wèi)星承力筒模型Fig.1 Satellite bearing cylinder model

通常認為結(jié)構(gòu)損傷并不引起質(zhì)量的損失，而只引起剛度損失。因此，本文采用剛度折減法模擬損傷，即對損傷部位的材料彈性模量折減來模擬局部損傷。設(shè)計9個工況，工況1為完好結(jié)構(gòu)，用于和有損結(jié)構(gòu)形成對比；工況2～5用來驗證損傷定位方法；工況3及工況6～9用于驗證損傷程度評估方法，見表1。每種工況下?lián)p傷區(qū)域均為沿軸向5行、周向25列的125個剛度折減單元的連續(xù)區(qū)域，如圖2。承力筒從下至上，每周5個測點，共11周，最下面一周測點編號為1～5，編號依次增大，最上面一周測點編號為51～55，在損傷區(qū)域附近布置4個測點，一個測點位于損傷區(qū)域中心，其它三個測點位于損傷區(qū)域邊緣，圖3是承力筒的平面展開圖。各工況對應(yīng)的損傷區(qū)域附近測點編號見表2。

表 1 工況設(shè)計Table 1 Design of working condition

注：h表示損傷區(qū)域距離筒體下端的距離占筒體高的百分比，損傷程度s為剛度折減率

圖2 損傷區(qū)域Fig.2 Damage region

圖3 測點和損傷區(qū)域位置(沿承力筒母線展開)Fig.3 Measuring point and damage region location (unfold drawing along cylinder generatrix line)

工況損傷區(qū)域附近測點編號21,2,6,73、6、7、8、916,17,21,22431,32,36,37546,47,51,52

3.2 損傷定位驗證

仿真驗證的總體流程為：采用MSC/Nastran軟件對承力筒進行工況1～9下20～2000 Hz的頻響分析，加速度激勵施加在星箭連接環(huán)基礎(chǔ)上，利用PCL提取55個測點的FRF數(shù)據(jù)。利用Matlab軟件編制PCA程序，提取FRF主元和對應(yīng)的貢獻率，構(gòu)造出兩種損傷定位指標，驗證算法可行性。

3.2.1 構(gòu)造損傷識別矩陣

以工況1(完好結(jié)構(gòu))和不同損傷工況2～5(損傷結(jié)構(gòu))下各測點的FRF曲線幅值差為變量，以每一測點得到FRF數(shù)據(jù)為樣本，組成55行(測點數(shù))900列(譜線數(shù))組成損傷識別矩陣H(ω)。

3.2.2 求各階主元及相應(yīng)的貢獻率

利用工況2～5下對應(yīng)的結(jié)構(gòu)損傷識別矩陣，求其協(xié)方差陣的900個特征值和特征向量，提取前10階主元，對應(yīng)的貢獻率如圖4和表3。

圖4 前10階主元貢獻率Fig.4 First ten principle components contribution rate

工況2345前10階累積貢獻率0.90400.91060.89660.9292

從表3和圖4可以看出，各工況下的前10階主元累積貢獻率大于89%，明顯大于其它各階成分，說明前10階主元包含了損傷識別矩陣900個變量的絕大部分信息。因而僅利用前10階主元所包含的數(shù)據(jù)信息進行分析損傷識別是可行的。

3.2.3 兩種損傷定位指標識別結(jié)果及分析

利用多主元綜合指標對工況2～5四種損傷狀態(tài)進行識別，置信度為95%。識別結(jié)果見圖5。

圖5 多主元綜合指標仿真結(jié)果Fig.5 Simulation of multi-principal component comprehensive index

歸一化的FRF主元殘差均方根法，表征了各測點間的相對損傷程度。損傷程度從0～1，0代表無損傷，1代表損傷程度最嚴重。根據(jù)2σ準則，各工況相應(yīng)的CL如圖6所示。

圖6 FRF主元殘差均方根指標損傷識別結(jié)果Fig.6 Simulation of FRF principle component residual mean square root

該兩種損傷定位指標識別得到的異常點，即落在CL以外的測點編號，匯總對比后，如表4所示。

表4 兩種損傷定位指標識別結(jié)果對比Table 4 Contrast of identification result between the two damage location indexes

對以上分析得到以下幾點結(jié)果：

(1)多主元綜合指標識別結(jié)果存在錯判(測點11、41)，而FRF主元殘差均方根無錯判。4個工況下，實際損傷區(qū)域測點為16個，前者識別到13個異常點，有2個錯判，準確率達84.6%；而后者識別到11個異常點，無錯判，準確率達100%。兩種指標基本能準確識別到損傷位置。測點11、41之所以被錯判，可能是該測點距離損傷區(qū)域比較近(測點11位于工況2損傷區(qū)域軸向正上方，測點41位于工況5損傷區(qū)域軸向正下方)，軸向傳力特性比較明顯，且遠離節(jié)點，響應(yīng)較大，因此該測點受損傷區(qū)域的影響較大，表現(xiàn)為異常點。

(2)多主元綜合指標沒有識別出測點1、16、31、46，存在漏判。原因可能是測點1、16、31、46位于損傷區(qū)域中心，其它測點位于損傷區(qū)域和完好區(qū)域的過渡區(qū)，結(jié)構(gòu)特性發(fā)生突變，更容易被檢測出來。而利用FRF主元殘差均方根指標，工況3、4下能夠識別到損傷區(qū)域中心的測點。因此，綜合參考兩種損傷識別指標，在識別到的異常點附近進行檢查校核，能夠有效防止漏判情況。

(3)FRF殘差均方根方法在工況2漏判情況較其它工況嚴重，漏判率達50%。原因可能是，承力筒約束方式為固支，1、2測點位于底端，響應(yīng)較小，信噪比太低，不容易被識別出來。故在實際測試中，對于遠離約束端和節(jié)點位置的響應(yīng)較大的測點，更容易被識別出來。

(4)損傷定位可以采用“先軸向、再周向”原則。從圖5和圖6中可以看出，從軸向和周向兩個維度看，越靠近損傷區(qū)域，定位指標數(shù)值越大。如圖6(b)中，測點6、11、16、21、26、31位于損傷區(qū)域軸向正上(下)方，比同一高度其它測點指標更大；而16、17、18、19、20位于損傷區(qū)域周向左右，比同一周向其它測點指標更大。位于承力筒上的損傷區(qū)域附近或與損傷區(qū)域同周向角度、同軸向高度的測點，兩種損傷定位指標較高，更有利于被識別出來。

3.3 損傷程度評估

工況3和6～9代表結(jié)構(gòu)損傷區(qū)域不同的材料剛度折減率，用于驗證損傷程度評估算法，求得對應(yīng)的SPE指數(shù)，如表5和圖7。

表5 不同損傷程度下SPE指數(shù)Table 5 SPE indexes under different damage degrees

圖7 不同損傷程度下SPE變化趨勢Fig.7 SPE tendency at various damage degree

從圖7中和表5中可以看出，工況1下(無損結(jié)構(gòu))頻響數(shù)據(jù)的對應(yīng)的SPE指數(shù)為3 219.1，隨著剛度折減率逐漸增加，結(jié)構(gòu)損傷程度逐漸加重，從10%增加到80%，對應(yīng)的SPE指數(shù)也逐漸增加，從3 340.4逐漸到10 474.0。因此，可以將工況1對應(yīng)的SPE指數(shù)設(shè)為閾值，大于該值說明損傷發(fā)生，偏離該值越大，說明結(jié)構(gòu)整體損傷程度越嚴重。由此可以看出，SPE指數(shù)既可以作為結(jié)構(gòu)整體損傷程度的評估，又可以作為結(jié)構(gòu)損傷開始的預(yù)示。

4 結(jié)論

本文以航天器振動試驗中預(yù)、復(fù)振頻響不一致問題為背景，將主元分析和損傷識別理論應(yīng)用于預(yù)、復(fù)振頻響數(shù)據(jù)的分析，提出了兩種基于頻響主元的損傷定位指標，分別是多主元綜合指標和FRF主元殘差均方根指標，并用SPE指數(shù)對結(jié)構(gòu)整體損傷程度進行了評估。通過衛(wèi)星典型主承力部件承力筒進行了數(shù)值仿真，驗證了上述損傷識別方法的有效性，得到了以下結(jié)論：

(1) 測點利用頻響數(shù)據(jù)構(gòu)造損傷識別矩陣，其前10階主元累積貢獻率達到了90%左右，包含了FRF數(shù)據(jù)矩陣的大部分信息，因此利用前10階主元進行損傷識別是可行的。

(2) 兩種損傷識別指標均能基本準確識別出損傷區(qū)域的大致位置。多主元綜合指標可以識別到損傷區(qū)域邊緣及其附近遠離節(jié)點、響應(yīng)較大的測點，但無法識別到損傷區(qū)域內(nèi)部的測點；FRF主元殘差均方根指標可以準確識別到損傷區(qū)域邊緣，并能夠識別到部分損傷區(qū)域內(nèi)部的測點。實際識別時，應(yīng)該綜合考慮兩種指標，減小錯判和漏判的可能，提高損傷定位效率。

(3)位于承力筒上的損傷區(qū)域附近或與損傷區(qū)域同周向角度、同軸向高度的測點，兩種損傷定位指標較高，更利于被識別出來。因此實際識別時，可以采取“先軸向、后周向”的損傷定位原則，先定位出損傷區(qū)域所在的軸向高度，再定位出損傷區(qū)域的周向角度。

(4) SPE指數(shù)可以表征結(jié)構(gòu)的整體損傷程度。隨著損傷程度的加深，SPE指數(shù)逐漸增大。以無損結(jié)構(gòu)得到的SPE指數(shù)為CL，既可以預(yù)示結(jié)構(gòu)發(fā)生損傷，也可以評估結(jié)構(gòu)的整體損傷程度。

本文提出的方法對衛(wèi)星部件(承力筒)的損傷定位和評估有效，是否適用于更復(fù)雜的衛(wèi)星系統(tǒng)還有待進一步驗證。該方法對基于預(yù)、復(fù)振對比分析的損傷識別、量化試驗評估標準及提高衛(wèi)星研制效率有促進作用。

References)

[1] 柯受全,金恂叔.衛(wèi)星環(huán)境工程和模擬試驗[M].中國宇航出版社,1996

Ke Shouquan,Jin Xunshu. Satellite environmental engineering and simulation test[M]. Beijing: China Astronautics Press,1996 (in Chinese)

[2]向樹紅. 航天器力學環(huán)境試驗技術(shù)[M]. 北京: 中國科學技術(shù)出版社,2008: 18-58

Xiang Shuhong. Spacecraft mechanics environment test technology[M]. Beijing: China Science and Technology Press,2008:18-58 (in Chinese)

[3]Nobukatsu O kuizumi,M C Natori. Nonlinear vibrations of a satellite truss structure with gaps[C]//45thAIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics & Materials Confer. Washington D.C.: AIAA,2004

[4]衛(wèi)洪濤,孔憲仁,王本利,等. 非線性連接結(jié)構(gòu)對一個典型衛(wèi)星頻率漂移的影響[J]. 航天器環(huán)境工程,2012,29(3): 297-302

Wei Hongtao,Kong Xianren,Wang Benli,et al. Research about the influence of nonlinear joints on a typical satellite frequency shift[J]. Spacecraft Environment Engineering,2012,29(3): 297-302 (in Chinese)

[5]陳昌亞,宋漢文,王德禹,等.衛(wèi)星振動頻率漂移現(xiàn)象與非線性參數(shù)識別[J]. 上海交通大學學,2005,39(7): 1197-1200

Chen Changya,Kong Xianren,Wang Deyu,et al. The natural frequency shift of satellite vibration test and parameter identification of nonlinear in satellite structure[J]. Journal of Shanghai Jiaotong University,2005,39(7): 1197-1200 (in Chinese)

[6]吉桂秀. 含多分層復(fù)合材料層合板和加筋板的自由振動特性研究[D]. 大連: 大連理工大學,2007

Ji Guixiu. Analysis of free vibration characteristics for multi-delaminated composite plates and stiffened plates[D].Dalian: Dalian University of Technology,2007 (in Chinese)

[7]Anderson J R,Van Noord J L,Soulas G C. Environmental testing of the NEXT PM1 ion engine[R]. Washington D.C.: AIAA,2007

[8]Dervan J,Robertson B,Staab L,et al. Development testing and subsequent failure investigation of a spring strut mechanism[C]// Proceedings of the 42nd Aerospace Mechanisms Symposium. Greenbelt: NASA Goddard Space Flight Center,2014: 1-14

[9]王東升,周桐,張志旭. 淺析振動試驗中結(jié)構(gòu)非線性引發(fā)的現(xiàn)象[J]. 航天器環(huán)境工程,2007,24(5): 314-317

Wang Dongsheng,Zhou Tong,Zhang Zhixu. A simple phenomenon analysis for structure nonlinearity in vibration tests[J]. Spacecraft Environment Engineering,2007,24(5): 314-317 (in Chinese)

[10] Carden E P. Vibration based condition monitoring: a review[J]. Structural Health Monitoring,2004,3(4):355-377

[11]Lee U,Shin J. A frequency response function-based structural damage identification method[J]. Computers & Structures,2002,80(2):117-13

[12]Maia N M M,Silva J M M,Almas E A M,et al. Damage detection in structures: from mode shape to frequency response functionmethods[J]. Mechanical Systems & Signal Processing,2003,17(3):489-498

[13]Zang C,Imregun M. Structural damage detection using artificial neural networks and measured FRF data reduced via principal component projection[J]. Journal of Sound and Vibration,2001,242(5): 813-827

[14]Ni Y Q,Zhou X T,Ko J M. Experimental investigation of seismic damage identification using PCA-compressed frequency response functions and neural networks[J]. Journal of Sound & Vibration,2006,290(1-2):242-263

[15]Dackermann U,Li J,Samali B,et al. Structural damage identification utilising PCA-compressed frequency response functions and neural network ensembles[C]//Proceedings of 20th Australasian Conference on the Mechanics of Structures and Materials. Toowoomba: Faculty of Engineering and Surveying,2008:803-809

[16]Nguyen V V,Dackermann U,Li J,et al. Damage identification of a concrete arch beam based on frequency response functions and artificial neural networks[J]. Electronic Journal of Structural Engineering,2015,14(1):75-84

[17]楊彥芳,宋玉普,紀衛(wèi)紅. 基于實測頻響函數(shù)主成分的在役網(wǎng)架損傷識別方法[J]. 振動與沖擊,2007,26(9): 128-132

Yang Yanfang,Song Yupu,Ji Weihong. A damage identification method for applied truss structure based on measured frequency response functions principle component[J]. Journal of Vibration and Shock,2007,26(9): 128-132 (in Chinese)

[18]Maia N M M,Silva J M M,Sampaio R P C. Localization of damage using curvature of the frequency res-ponse functions[J]. Proceedings of SPIE,1997,3089: 942-946

[19]姜增國,張楨. 基于應(yīng)變頻響函曲率的結(jié)構(gòu)損傷識別[J]. 建筑科學與工程學報,2009,26(4): 40-43

Jiang Zengguo,Zhang Zhen. Structural damage identification based on SFRF curvature[J]. Journal of Architecture and Civil Engineering,2009,26(4): 40-43 (in Chinese)

[20]朱新圓,阿肯江·托呼提. 鋼桁架結(jié)構(gòu)損傷識別的頻響函數(shù)曲率法[J]. 實驗室研究與探索,2014,33(6): 40-44

Zhu Xinyuan,Akenjiang·Tuohuti. A damage identification study of steel truss structure based on frequency response function[J]. Research and Exploration in Laboratory,2014,33(6): 40-44 (in Chinese)

[21]Johnson R A,Wichern D W.實用多元統(tǒng)計分析[M]. 北京：清華大學出版社，2001

Johnson R A,Wichern D W. Practical multivariate statistical analysis[M]. Beijing: Tsinghua University Press,2001 (in Chinese)

(編輯：張小琳)

Methods of Damage Location and Evaluation for Satellite Structure Based on Frequency Response Principle Component

FENG Yanjun SHI Chuanqian ZHOU Yu ZHAO Fagang ZHOU Xubin

(Shanghai Institute of Satellite Engineering,Shanghai 201109,China)

Considering that current evaluation methods for satellites’ pre and post vibration frequency response mainly depend on designers’ experience,lacking efficient quantization methods,this paper first puts forward two damage location indexes for satellite structure based on principle component analysis (PCA) of frequency response function (FRF). The two indexes are mulit-principal component comprehensive index and FRF principle component residual mean square root. Then,square prediction error (SPE) index is adopted to evaluate the whole damage degree for satellite structure. Finally,by taking satellite bearing cylinder as analysis objective, the efficiency of the damage location and evaluation methods based on PCA is validated. This paper brings positive effect on damage identification based on pre and post vibration analysis,quantization of experiment evaluation standard and improvement of satellite manufacture.

damage identification; multi-principle component comprehensive index; FRF principle component residual mean square root; SPE index; satellite pre and post vibration; frequency res-ponse function

2016-06-27；

2016-07-12

國家重大航天工程

馮彥軍，男，碩士研究生，從事衛(wèi)星結(jié)構(gòu)設(shè)計和振動測試工作。Email:feng_yan_jun@126.com。

V416.2；V414.6

A

10.3969/j.issn.1673-8748.2016.04.018