宋海洋，張文正，張豐收，呂 稀，馮志鵬

（中國核動力研究設計院 核反應堆系統(tǒng)設計技術重點實驗室，成都 610213）

新型工程產(chǎn)品和項目在研制初期一般都需要獲得相應的聲振環(huán)境預示數(shù)據(jù)，在低頻段可以采用傳統(tǒng)的有限元法或邊界元法，但在高頻段由于模態(tài)十分密集，重疊率高使得傳統(tǒng)的有限元法應用起來十分困難。20世紀60年代發(fā)展起來的統(tǒng)計能量分析方法可以克服傳統(tǒng)模態(tài)分析方法所遇到的這些困難，為解決復雜系統(tǒng)的高頻振動及噪聲問題提供了有力工具。

利用統(tǒng)計能量分析方法對工程結構的高頻振動和噪聲進行預示時，內(nèi)損耗因子和耦合損耗因子等參數(shù)的測量精度對響應預示結果的精確性有著很大的影響。目前大部分都采用實驗測量的方法來獲得結構的內(nèi)損耗因子和耦合損耗因子[1–3]，但是根據(jù)前人所做的大量實驗和所獲得的工程經(jīng)驗可知，對于大部分工程結構來說，內(nèi)損耗因子和耦合損耗因子一般都是10-2～10-4量級的小數(shù)，在實際工作中想要對這些小數(shù)進行精確運算和測量是比較困難的[4–5]，最終的測量數(shù)據(jù)通常存在誤差，而參數(shù)測量誤差的積累會導致預示出的高頻響應和響應的真實值之間存在偏差，目前尚未有相關文獻報道關于參數(shù)的測量誤差對預示結果影響的研究。

本文在考慮參數(shù)以及外載荷測量誤差的前提下，結合區(qū)間攝動方法推導出了統(tǒng)計能量分析子系統(tǒng)的模態(tài)能量區(qū)間，利用模態(tài)能量區(qū)間可進一步獲得不同頻率點處的總能量區(qū)間?？偰芰繀^(qū)間的均值代表不考慮參數(shù)和外載荷測量誤差時預示的動響應，而總能量區(qū)間上下界則反映了參數(shù)以及外載荷的測量誤差對響應預示結果的影響程度。

1 區(qū)間的定義和運算

在區(qū)間分析中，不確定性變量可以用一個封閉的和有限的區(qū)間來表示。一個區(qū)間變量xI可以由它的區(qū)間下限xl和區(qū)間上限xu來定義，具體的表達形式可以寫為[6]

本文使用上標“l(fā)”代表區(qū)間變量的下限，使用上標“u”代表區(qū)間變量的上限，同時使用上標“c”代表區(qū)間變量的均值，區(qū)間xI的均值定義為

區(qū)間xI的離差定義為

兩個區(qū)間變量xI和yI之間的運算規(guī)則如下[7]

2 統(tǒng)計能量分析功率流平衡方程

2.1 子系統(tǒng)間功率流

在統(tǒng)計能量分析理論中，子系統(tǒng)所儲存的能量通過兩種方式流出，一是子系統(tǒng)本身的能量耗散，能量被轉化成熱能等其它形式，二是能量從一個子系統(tǒng)傳遞到了另一個子系統(tǒng)。第i個統(tǒng)計能量分析子系統(tǒng)的損耗功率pid可用下式進行計算

其中f代表分析帶寬Δf內(nèi)的中心頻率，Ei代表子系統(tǒng)i在頻率點f處的總能量，ηi是子系統(tǒng)i的內(nèi)損耗因子，它表示子系統(tǒng)i自身耗散能量的能力大小。

子系統(tǒng)i到子系統(tǒng)j的功率流為

其中ηij和ηji是子系統(tǒng)間的耦合損耗因子，它們代表能量在子系統(tǒng)i和子系統(tǒng)j之間的傳遞效率。

統(tǒng)計能量分析中的耦合損耗因子和模態(tài)密度存在一個互易關系

其中di表示子系統(tǒng)i的模態(tài)密度。

2.2 功率流平衡方程

假設一個統(tǒng)計能量分析系統(tǒng)被劃分為k個相互耦合的子系統(tǒng)，如圖1所示

圖1 k個相互耦合的統(tǒng)計能量分析子系統(tǒng)

由式(5)和式(6)可得子系統(tǒng)i的功率流平衡方程

其中pi代表子系統(tǒng)i的外部激勵源對子系統(tǒng)i的輸入功率。為子系統(tǒng)i的能量Ei的變化率。本文研究的是穩(wěn)態(tài)振動情況所以即子系統(tǒng)能量不隨時間而變化。

利用式(7)給出的互易原理，將所有子系統(tǒng)的功率流平衡方程寫成矩陣形式如下

其中N是一個包含模態(tài)密度di(i=1,2,…,k) 、內(nèi)損耗因子ηi(i=1,2,…,k)和耦合損耗因子ηij(i,j=1,2,…,k)的k×k階系統(tǒng)矩陣，f是分析帶寬Δf內(nèi)的中心頻率，為計算方便令Pi=pi/πp f，ei(i=1,2,…,k)稱為模態(tài)能量指的是在分析帶寬內(nèi)每階模態(tài)所具有的能量。統(tǒng)計能量分析功率平衡方程是在分析帶寬內(nèi)各階模態(tài)所具有的能量相等的假設下建立的，子系統(tǒng)i在各頻率點處總能量Ei和模態(tài)能量關系為

3 基于區(qū)間攝動法的穩(wěn)態(tài)響應區(qū)間計算

3.1 子系統(tǒng)模態(tài)能量攝動分析

利用式(9)，子系統(tǒng)的模態(tài)能量向量可以通過系統(tǒng)矩陣N的逆陣和功率向量P的乘積來獲得

根據(jù)攝動理論，當系統(tǒng)矩陣N存在一個攝動量ΔN，功率向量P存在一個攝動量ΔP時，攝動法的平衡方程可以表示為

e+Δe能被進一步寫為

當矩陣范數(shù) ‖N-1+ΔN‖小于1時，可以利用紐曼級數(shù)將(N+ΔN)-1展開如下

將式(14)代入到式(13)中得到

忽略式(15)的高階微量，可以得到模態(tài)能量向量的攝動量為

3.2 帶有測量誤差參數(shù)的區(qū)間表示

用實驗測量的方式來獲取結構的內(nèi)損耗因子、耦合損耗因子和輸入功率時，測量結果通常存在一定量的誤差，下面采用式(1)的區(qū)間形式來表示帶有測量誤差的參數(shù)和輸入功率

內(nèi)損耗因子、耦合損耗因子和輸入功率的測量值是式(17)的區(qū)間均值即，假設測量誤差不超過±s%，則式(17)可以進一步寫成

相比于內(nèi)損耗因子和耦合損耗因子的測量精確性，模態(tài)密度的測量結果相對精確很多。再者任何區(qū)間分析方法都存在區(qū)間擴張問題，即過多的區(qū)間相乘或相除運算會導致最終的計算結果不精確。綜合以上兩方面原因，將子系統(tǒng)的模態(tài)密度處理成常數(shù)而不是區(qū)間變量。

3.3 子系統(tǒng)模態(tài)能量區(qū)間

本節(jié)將帶有測量誤差的內(nèi)損耗因子、耦合損耗因子和輸入功率(即式(17))代入到統(tǒng)計能量分析的功率流平衡方程，最終可以得到子系統(tǒng)模態(tài)能量區(qū)間，區(qū)間均值是不考慮參數(shù)和輸入功率誤差時預示的模態(tài)能量，而模態(tài)能量區(qū)間的上下限反映了參數(shù)和輸入功率的測量誤差對預示結果的影響程度和大小。將式(17)代入到式(9)中得區(qū)間形式的功率平衡方程為

這里NI稱作區(qū)間系統(tǒng)矩陣，它是一個包含模態(tài)密 度di(i=1,2,…,k)，內(nèi)損耗因子區(qū)間和耦合損耗因子區(qū)間的矩陣，是子系統(tǒng)i的模態(tài)能量區(qū)間。

用區(qū)間均值和不確定區(qū)間來表示式(19)中的區(qū)間系統(tǒng)矩陣和區(qū)間輸入功率向量

將式(20)和式(21)代入到式(19)中得到

ΔNI可以看作是圍繞Nc的一個攝動，而ΔPI可以看作是圍繞Pc的一個攝動，這樣就可以將式(26)看作是一個攝動形式的方程即

將式(30)代入到式(26)中得到

根據(jù)式(31)可以獲得子系統(tǒng)模態(tài)能量的區(qū)間向量為

基于式(19)的區(qū)間功率流平衡方程，式(32)中的Nc和ΔNI可以表示為

假設Nc的逆矩陣為

將式(34)和式(35)代入到式(32)中得

將式(36)和式(37)相加，最終可以得到任意一個子系統(tǒng)t(t=1,2,…,k)模態(tài)能量區(qū)間的上界和下界的計算公式

4 實驗驗證

圖2是兩塊垂直相連的復合板結構，本節(jié)引用文獻[8]針對圖2典型板結構的振動實驗結果，并將該振動實驗結果與利用區(qū)間方法預示得到的響應結果進行對比。

圖2 實驗現(xiàn)場圖[8]

右側板1的長度為400 mm，寬度為500 mm，左側板2的長度為600 mm，寬度為500 mm，復合板整體厚度為1.2 mm。兩塊復合板被劃分為兩個統(tǒng)計能量分析子系統(tǒng)，板1為子系統(tǒng)1，板2為子系統(tǒng)2。兩個板分別布置了10和15個加速度傳感器，每個板布置5個激振點。依次對各個激振點激勵并測量激振力和激振點的加速度響應，對于穩(wěn)態(tài)實驗數(shù)據(jù)，由阻抗頭測量的加速度信號的互譜密度函數(shù)計算穩(wěn)態(tài)輸入功率，由各板加速度信號的自譜密度函數(shù)的平均計算子系統(tǒng)的總能量，然后將子系統(tǒng)總能量和輸入功率的比值作為歸一化的總能量。

表1 兩塊板的內(nèi)損耗因子和耦合損耗因子

假設表1中數(shù)據(jù)的測量誤差不超過±8%，那么根據(jù)式(17)和式(18)的定義，可以將所有帶有測量誤差的參數(shù)都寫成區(qū)間形式，其中區(qū)間的均值就是實際的參數(shù)測量值。

表2 用區(qū)間形式表示帶有測量誤差(不超過±8%)的參數(shù)

當統(tǒng)計能量分析的參數(shù)區(qū)間確定以后，需要驗證一下是否滿足式(30)的收斂條件，即矩陣范數(shù)小于1。

當分析帶寬的中心頻率為640 Hz時，可得

經(jīng)計算可得兩板的總能量區(qū)間如圖3和圖4所示。

圖3和圖4給出了由三種方式得到的板振動總能量：

(1)在考慮所有內(nèi)損耗因子測量誤差、耦合損耗因子測量誤差以及外載荷測量誤差的情況下，由本文提出的區(qū)間方法從理論上預示出的板總能量區(qū)間值；

(2)在不考慮任何內(nèi)損耗因子測量誤差、耦合損耗因子測量誤差以及外載荷測量誤差的情況下，由非區(qū)間的傳統(tǒng)方法從理論上預示出的板總能量確定值(即區(qū)間均值)；

圖3 板1在各頻率點上的總能量區(qū)間

圖4 板2在各頻率點上的總能量區(qū)間

(3)由實驗測量得到的板總能量。

圖3和圖4中的橫坐標選取的是統(tǒng)計能量分析理論1/3倍頻程的幾個中心頻率點，在統(tǒng)計能量分析理論中用中心頻率點處的能量來代表分析帶寬內(nèi)能量的一個平均值，比如中心頻率640 Hz的下限和上限頻率分別是596 Hz和716 Hz，那么在統(tǒng)計能量分析中以640 Hz點處的能量代表從596 Hz到716 Hz帶寬內(nèi)的一個平均能量。由圖3和圖4可見板結構的總能量區(qū)間包含著總能量的真實值(即實驗測量值)。

為了更具體地衡量參數(shù)和外載荷的測量誤差對穩(wěn)態(tài)響應預示結果的影響程度，定義

基于圖3和圖4給出的考慮參數(shù)和載荷測量誤差時預示得到的板總能量區(qū)間上下限以及不考慮任何參數(shù)和載荷測量誤差時預示得到的板總能量的確定值(即區(qū)間均值)，現(xiàn)給出由參數(shù)和外載荷的測量誤差所導致的穩(wěn)態(tài)響應的誤差值如表3。

表3 由參數(shù)和載荷的±8%測量誤差所導致的穩(wěn)態(tài)響應誤差

5 結語

對于復雜的工程結構在應用穩(wěn)態(tài)統(tǒng)計能量分析方法預示高頻段動響應時，往往并沒有考慮到參數(shù)和外載荷的測量誤差，導致最終的動響應預示結果不夠精確，給后續(xù)的結構設計帶來不準確的信息，所設計出來的結構也會存在一定的安全隱患，在實際惡劣的聲振環(huán)境下破壞或者失效，因此在進行動響應預示時有必要考慮到參數(shù)以及外載荷的測量誤差對預示結果的影響。本文考慮到穩(wěn)態(tài)統(tǒng)計能量分析模型中損耗因子和外載荷的不確定性，將帶有測量誤差的內(nèi)損耗因子和耦合損耗因子等參數(shù)以區(qū)間變量形式表示。通過求解統(tǒng)計能量分析區(qū)間功率流平衡方程獲得了每個子系統(tǒng)的總能量區(qū)間，通過這個總能量區(qū)間可以得到參數(shù)和外載荷的測量誤差對預示響應的影響程度和大小，為后續(xù)的結構安全和可靠性設計提供參考。

[1]PANUSZKA R,WICIAK J,IWANIEC M.Experimental assessment of coupling loss factors of thin rectangular plates[J].Archives ofAcoustics,2005,30(4):533-551.

[2]EWING M S,DANDE H.Damping loss factor estimation for coupled plates using experimental transient statistical energy analysis[C].52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference,Denver,United states,2011:1-7.

[3]PEREIRA R,ARENAS J P,ZUMELZU E.Comparison of four test methods to measure damping properties of materials by using piezoelectric transducers[J].Materials and Design,2011,32(4):2423-2428.

[4]姚德源，王其政.統(tǒng)計能量分析原理及其應用[M].北京：北京理工大學出版社，1995.

[5]LYON R H,DEJONG R G.Theory and application of statistical energy analysis(Second Edition)[M].The MIT Press,Cambridge.1995.

[6]郭書祥,呂震宙.區(qū)間運算和靜力區(qū)間有限元[J].應用數(shù)學和力學,2001,22(12):1249-1254.

[7]DEGRAUWE D,LOMBAERT G,DE ROECK G.Improving interval analysis in finite element calculations by meansofaffine arithmetic[J].Computersand Structures,2010,88(3-4):247-254.

[8]張紅亮.統(tǒng)計能量分析參數(shù)辨識及模型修正方法研究[D].哈爾濱：哈爾濱工業(yè)大學，2013.