胡杰

（1.中國工程物理研究院 總體工程研究所，四川 綿陽 621999；2.工程材料與結(jié)構(gòu)沖擊振動四川省重點實驗室，四川 綿陽 621999）

裝備結(jié)構(gòu)在定型并投入到實際工程應(yīng)用之前，往往都需要通過環(huán)境試驗對其環(huán)境適應(yīng)性進行考核，其目的是在試驗室中通過模擬實際環(huán)境的載荷工況或者使得產(chǎn)品產(chǎn)生與其在實際工況下等效的響應(yīng)，以考察裝備結(jié)構(gòu)是否能夠滿足各項指標(biāo)要求。其中，振動環(huán)境對大多數(shù)裝備結(jié)構(gòu)而言，都是在實際使用過程中不可避免的載荷工況，因此在環(huán)境試驗中扮演著重要的角色。

振動環(huán)境試驗中，最常采用的加載設(shè)備為振動臺，試驗中一般是在臺面或者產(chǎn)品上設(shè)置響應(yīng)控制點，使得控制點產(chǎn)生預(yù)期的響應(yīng)。目前最常見的形式為單點激勵方式，即采用單個振動臺加載。這種試驗方式應(yīng)用比較成熟[1]，諸多試驗規(guī)范和標(biāo)準也制定了相關(guān)的規(guī)定[2]。但該試驗方式只能對單個測點的響應(yīng)進行控制，產(chǎn)品上其他部位的響應(yīng)往往不受控，即只能保證單個測點的響應(yīng)達到預(yù)期目的。當(dāng)前試驗中常采用的多點平均響應(yīng)控制方式實質(zhì)上仍屬于單點控制，只不過是對某個虛擬點的控制而已，多個點的平均響應(yīng)滿足控制要求，但各個測點之間差異往往很大。

從上述分析來看，如果振動環(huán)境試驗中能夠在多個點上實現(xiàn)響應(yīng)控制，那么產(chǎn)品的響應(yīng)將能更真實地與其在實際使用環(huán)境下等效[3-5]，尤其是針對尺寸較長的結(jié)構(gòu)[6]。此時就需要借助多點激勵方式，即搭建多個加載設(shè)備（如多個振動臺）構(gòu)成的激勵系統(tǒng)。劉志華等[7]針對正弦激勵情況，提出了多點激勵正弦振動實數(shù)域控制算法的控制流程。對于振動環(huán)境中最為關(guān)注的隨機振動工況，20世紀80年代由Small-wood等[8-9]初步建立了多點隨機振動控制方法，建議根據(jù)重要性程度，對譜矩陣中的對角線元素進行排列，靠前的元素優(yōu)先收斂，但無法實現(xiàn)多點加權(quán)綜合控制。劉青林等[10]采用柔性懸掛-多點激勵試驗方法對某彈體結(jié)構(gòu)飛行過程中的面載荷進行了有效模擬。謝志勇等[11]研究了結(jié)構(gòu)在力型輸入和加速度型輸入條件下，多點激勵的振動響應(yīng)求解方法，結(jié)果表明，采用力型輸入與實驗吻合最好。韓偉等[12]應(yīng)用虛擬試驗技術(shù)在試驗前進行試驗方案設(shè)計及優(yōu)化，驗證了雙點激勵虛擬試驗的準確和可靠性。李永華等[13]則在數(shù)值模擬方法上，對多點激勵下的響應(yīng)求解開展了研究，提出了多點激勵下絕對位移直接求解的虛擬激勵法，并進一步對該算法開展了誤差的頻域分析[14]。

多點激勵的難點主要是多點激勵之間相關(guān)性（體現(xiàn)為互譜）的控制，體現(xiàn)了多個激勵載荷之間的協(xié)調(diào)配合。賀旭東等[15-16]從試驗控制的角度對互譜矩陣開展了仿真研究，基于矩陣分解提出了獲得試驗驅(qū)動信號的控制算法，但對于互譜元素非零的情況還有待進一步深入研究。葉建華[17]開發(fā)了雙振動臺隨機振動試驗控制系統(tǒng)樣機，對2個控制點的不同響應(yīng)功率譜密度分別進行了控制。祝濟之等[18]將2個響應(yīng)控制點設(shè)置在夾具與試驗件的連接處，較好地將2個振動臺的激勵傳遞到試驗件上。

從上述分析來看，在目前多點激勵的研究中，一般考慮的是激勵數(shù)量與響應(yīng)控制點數(shù)量相同的情況，控制的響應(yīng)點越多，等效程度越高，試驗效果越好，但要求的加載設(shè)備也越多。實際情況下，搭建由多個加載設(shè)備構(gòu)成的協(xié)調(diào)激勵系統(tǒng)往往較為困難。如四點響應(yīng)控制中，搭建4個振動臺的激勵系統(tǒng)并非易事，且協(xié)調(diào)加載要求也很高，難度也更大。此時，若能夠通過少量的加載設(shè)備實現(xiàn)多點響應(yīng)控制將具有重要的實際應(yīng)用意義。其相應(yīng)的控制載荷是否存在，本文針對該問題，并結(jié)合實際試驗情況，對控制載荷的存在性從理論上進行分析。

1 目標(biāo)響應(yīng)與控制載荷間傳遞關(guān)系分析

經(jīng)過有限元離散化，頻域內(nèi)結(jié)構(gòu)的動力學(xué)方程為：

式（1）中：M、C、K分別為結(jié)構(gòu)的質(zhì)量、阻尼和剛度矩陣；Xa(ω)、Xv(ω)、Xd(ω)分別為結(jié)構(gòu)的加速度、速度和位移響應(yīng)；F(ω)為結(jié)構(gòu)受到的外載荷。

振動環(huán)境試驗中，常關(guān)注正弦振動和隨機振動兩種工況，且考慮到實際試驗中常關(guān)注結(jié)構(gòu)的加速度響應(yīng)，正弦振動下，結(jié)構(gòu)加速度的諧響應(yīng)為：

式（2）中：H(ω)為載荷與加速度響應(yīng)之間的傳遞函數(shù)矩陣。

記c為所關(guān)注的多點響應(yīng)集合，其維數(shù)為m；f為載荷施加的節(jié)點集合，其維數(shù)為n。則關(guān)注點c的加速度諧響應(yīng)為：

其中，Hc(ω)為H(ω)中與c所對應(yīng)的行和f對應(yīng)的列所構(gòu)成的矩陣，其維數(shù)為m×n。

對于隨機振動而言，關(guān)注點響應(yīng)的功率譜密度矩陣Pc(ω)為：

式中：dω為頻率分辨率[19]，為圓頻率。

由于諧響應(yīng)分析是隨機振動的基礎(chǔ)，可由諧響應(yīng)推導(dǎo)出隨機振動的響應(yīng)，因此在討論多點隨機振動下的控制載荷存在性時，可通過多點諧波激勵分析開展研究，具有等價性。

2 控制載荷存在性分析

將式（4）寫成完備矩陣形式，如式（6）所示，其中元素的虛部反映了響應(yīng)的相位信息。

記傳遞函數(shù)矩陣Hc的增廣矩陣Hz為：

另記Hc和Hz的秩分別為Rc和Rz，以下基于矩陣論中增廣矩陣對解的判定研究，討論多點響應(yīng)等效時控制載荷的存在性問題。

2.1 同時控制響應(yīng)幅值和相位

1）Rc=Rz=n，即傳遞函數(shù)矩陣的秩與其增廣矩陣的秩相等，且為滿秩。對實際工程結(jié)構(gòu)而言，一般意味著響應(yīng)控制點數(shù)量與激勵數(shù)量相等，此時式（6）存在唯一解[20]，控制載荷可由式（8）計算獲得。

將式（6）的虛部進一步改寫為相位角的形式，則有：

式中：α和β分別為激勵與響應(yīng)的相位角。

2）Rc=Rz＜n，對于工程結(jié)構(gòu)的傳遞函數(shù)矩陣而言，一般意味著響應(yīng)控制點數(shù)量小于控制載荷數(shù)量，此時式（6）存在無窮解[20]。

將式（6）改寫為：

式（10）中：Ha的維數(shù)為m×m，Hb的維數(shù)為m×(n?m)。將式（10）展開為：

當(dāng)Ff,b=0時，式（11）就退化到m=n時的情況。需要指出的是，Ha可由H矩陣中任意m列元素構(gòu)成。這種組合方式很多，相比m=n時的情況，提供了更豐富的加載組合方式。因此，控制載荷也不唯一，有利于多點激勵時響應(yīng)控制的實現(xiàn)。

3）Rc＜Rz，對于工程結(jié)構(gòu)的傳遞函數(shù)矩陣而言，一般意味著響應(yīng)控制點數(shù)量大于控制載荷數(shù)量，即期望通過少量的加載設(shè)備實現(xiàn)對多個點的響應(yīng)的控制，此時式（6）的解往往不存在[20]。

2.2 僅控制響應(yīng)幅值

以上討論都是從數(shù)學(xué)角度出發(fā)，是嚴格意義上的響應(yīng)等效分析，要求響應(yīng)功率譜密度矩陣Pc(ω)中的所有元素都得到控制，每個元素實部和虛部都等效。從動力學(xué)意義上來說，則表現(xiàn)為多個測點的響應(yīng)幅值以及多個測點之間的相位關(guān)系都是固定的，對控制精度要求較高，實現(xiàn)難度較大。目前在工程應(yīng)用中的實際情況是只控制測點響應(yīng)的幅值。以隨機振動試驗為例，即只控制式（5）中功率譜密度矩陣Pc(ω)的對角線元素，并不控制多個控制點之間響應(yīng)的相位，即Pc(ω)中的非對角線元素（互譜），多個控制點之間的相位可以是隨機的，因此等效程度有所放寬，此時控制載荷的存在性結(jié)論將發(fā)生變化：

1）Rc=Rz=n時，即在只控制響應(yīng)自譜的情況下，加載設(shè)備數(shù)量與響應(yīng)控制點數(shù)量相同時，由于式（9）中的β具有任意性，因此控制載荷Ff并不唯一，這種特性便于多點激勵時響應(yīng)控制的實現(xiàn)。

2）Rc=R＜n時，即在只控制響應(yīng)自譜的情況下，加載設(shè)備數(shù)量大于響應(yīng)控制點數(shù)量時，與上述分析類似，控制載荷仍然不唯一。

3）Rc＜Rz時，即在只控制響應(yīng)自譜的情況下，加載設(shè)備數(shù)量小于響應(yīng)控制點數(shù)量的情況。實際試驗中，往往也只能基于有限的加載設(shè)備搭建多維振動試驗系統(tǒng)，那么在不控制多點響應(yīng)之間相關(guān)性這種相對弱化的條件下，控制載荷與等效響應(yīng)幅值之間的傳遞關(guān)系如何？若能夠通過少量的加載設(shè)備實現(xiàn)對多個點響應(yīng)的控制，將具有重要的工程應(yīng)用意義。

以m=3、n=2，即兩點激勵、三點響應(yīng)控制為例，對控制載荷從數(shù)值求解的角度進行分析。為便于分析，取頻率分辨率dω=2π，即1 Hz，得到3個目標(biāo)響應(yīng)點的響應(yīng)功率譜密度Pcm(m=1、2、3)的表達式為[21]：

其中，θ=α1?α2，為兩點激勵之間的相位差。式（12）表明，兩點激勵時，影響3個控制點響應(yīng)的載荷參數(shù)有3個，即2個振動臺的臺面響應(yīng)功率譜和，以及它們之間的相位差θ。式（12）也表明了響應(yīng)功率譜與控制載荷之間的非線性關(guān)系，難以參考非齊次線性方程組解的存在性理論，說明控制載荷不存在。

文[21]中采用優(yōu)化算法對梁結(jié)構(gòu)在兩點激勵、三點響應(yīng)控制的控制載荷進行了數(shù)值模擬，結(jié)果表明，控制載荷在某些頻段內(nèi)并不存在，而在某些頻段范圍內(nèi)控制載荷仍然可以獲得，能夠滿足控制要求，凸顯了問題的復(fù)雜性。

上述分析是從數(shù)學(xué)意義上討論控制載荷的存在性，但對于實際工程而言，當(dāng)加載位置的選擇點不合理。如處于某階模態(tài)的振型節(jié)點附近，響應(yīng)控制點位于振型值較大的位置時，在該模態(tài)頻率附近要實現(xiàn)多點響應(yīng)控制將十分困難，或者說對推力的要求往往是非常大的，導(dǎo)致無法實現(xiàn)。以下通過某典型結(jié)構(gòu)的多種工況的數(shù)值仿真對該情況進行分析說明。

3 數(shù)值模擬

3.1 仿真模型

如圖1所示筒狀結(jié)構(gòu)，材料為鋼，長度為2000 mm，小端內(nèi)徑為300 mm，大端內(nèi)徑為536.8 mm，筒壁厚8 mm，質(zhì)量約280 kg。試驗中采用大小兩端柔索懸掛模擬自由邊界條件，響應(yīng)控制點為A1—A3，可根據(jù)不同的控制需求選擇控制點組合，多點激勵可供選擇的加載點為B1—B5。載荷類型為力，坐標(biāo)方向如圖1所示，軸向為x方向，橫向為y方向，響應(yīng)控制和加載的方向均為y方向。

圖1 筒結(jié)構(gòu)多點激勵示意Fig.1 Multi-point excitation sketch of tube shape structure

研究中設(shè)置的多種仿真工況見表1。其中，諧波載荷和隨機振動載荷下的目標(biāo)響應(yīng)控制曲線如圖2所示。諧波載荷激勵頻率范圍為5~1000 Hz；隨機振動下，激勵頻率范圍為10~1000 Hz，頻率分辨率為0.5 Hz，且3個控制點的響應(yīng)功率譜密度曲線設(shè)置一致。

圖2 各控制點的目標(biāo)響應(yīng)曲線Fig.2 Target control point response: a) harmonic excitation;b) random vibration

表1 控制載荷數(shù)值模擬工況Tab.1 Numerical simulation working condition of control load

3.2 仿真結(jié)果

工況①、②下，控制載荷的幅值和相位角的頻譜曲線如圖3所示。為便于分析，考慮2個控制點響應(yīng)之間相位差為0，此時，控制載荷是唯一的，不同的加載點選擇，所需要的控制力譜也不同。另外也可以看出，工況②下所需的控制載荷幅值要比工況①下大得多，這說明在試驗中，合理地選擇加載點能夠充分利用設(shè)備推力加載能力，以達到試驗?zāi)康摹?/p>

圖3 工況①、②的控制載荷Fig.3 Control load of working condition ①,②: a) amplitude;b) phase angle

工況③下，僅控制A1、A3兩點諧響應(yīng)幅值譜時的控制載荷曲線如圖4所示，圖中結(jié)果只是其中的一種解的情況。

圖4 工況③的控制載荷Fig.4 Control load of working condition ③: a) amplitude; b)phase angle

諧波激勵下，對A1、A2、A3的幅值和相位均控制時的控制載荷如圖5所示，3個響應(yīng)點之間的相位差設(shè)置為0。圖6則給出了僅控制幅值譜時的控制載荷，與工況③類似，也只是多解中的一種情況。

圖5 工況④的控制載荷Fig.5 Control load of working condition ④: a) amplitude; b)phase angle

圖6 工況⑤的控制載荷Fig.6 Control load of working condition ⑤: a) amplitude; b)phase angle

隨機載荷下，兩點激勵、三點響應(yīng)控制的控制載荷如圖7所示。只控制三點響應(yīng)的幅值，數(shù)值模擬中采用遺傳優(yōu)化算法，對兩點激勵的幅值以及它們之間的相位差進行優(yōu)化，同時，考慮了實際隨機振動試驗中控制容差為±3 dB。可以看出，在諸多頻段范圍內(nèi)，仍然能找到滿足控制容差的控制載荷。

圖7 工況⑥的控制載荷Fig.7 Control load of working condition ⑥: a) amplitude; b)phase angle

圖8 中，123 Hz附近的控制響應(yīng)的超差比較嚴重，主要表現(xiàn)為A1和A3點響應(yīng)較大，而A2點的響應(yīng)較小。通過模態(tài)分析，可以觀察到在該頻率點處存在第六階模態(tài)，頻率為121.6 Hz，其振型如圖9所示。A1、A2、A3各點的振型值分別為4.82×10?4、3.38×10?5、4.99×10?4，可見也正好表現(xiàn)為A1、A3點響應(yīng)大，A2的響應(yīng)小的特點，造成在模態(tài)頻率附近難以控制。

圖9 第六階模態(tài)振型Fig.9 The sixth mode shape

從圖8中也可以看出，在575 Hz左右，控制精度較高，而在123 Hz處，控制超差嚴重。將式（12）進一步改寫為：

圖8 工況⑥控制載荷產(chǎn)生的響應(yīng)與目標(biāo)響應(yīng)對比Fig.8 Comparison of control load caused response with target response in working condition⑥

其中：

進一步得到：

式（14）中，m=1、2、3，根據(jù)式（15），給出了在123 Hz和400 Hz處式（12）數(shù)值解的分布曲面，如圖10所示。

圖10 不同頻率點上控制載荷解的分布Fig.10 Solution distribution of control load at different frequencies

由圖10可見，在123 Hz處，3個解的分布曲面無交點，即無精確控制解，優(yōu)化過程中超差也較為嚴重。在575 Hz處，3個解的分布曲面有交點，有精確控制解，優(yōu)化仿真過程中，偏差很小。

4 結(jié)語

從結(jié)構(gòu)動力學(xué)的載荷與響應(yīng)之間的傳遞關(guān)系出發(fā)，基于諧響應(yīng)分析，將傳遞函數(shù)矩陣及其增廣矩陣的秩的大小關(guān)系對比作為控制載荷存在性的判據(jù)。從廣義上給出了多點激勵下，同時控制響應(yīng)幅值和相位時控制載荷的存在性，并結(jié)合典型結(jié)構(gòu)進行了控制載荷的數(shù)值模擬。結(jié)合實際振動試驗中只控制響應(yīng)幅值的情況，給出了僅控制幅值時控制載荷存在性結(jié)論，尤其是針對利用少量加載設(shè)備實現(xiàn)多點響應(yīng)幅值控制的情況，其控制載荷在某些頻段上仍然存在。