3種主動控制算法的結構振動控制效果研究

樊文華，鄒萬杰

（廣西科技大學土木建筑工程學院，廣西柳州545006）

3種主動控制算法的結構振動控制效果研究

樊文華，鄒萬杰*

（廣西科技大學土木建筑工程學院，廣西柳州545006）

針對不同類型的結構，由Matlab中的相關函數命令求出相應的控制力狀態反饋增益矩陣，繼而得出結構振動的位移、加速度和控制力，比較分析了3種主動控制算法的控制效果，并給出了基于Matlab的三層框架結構的算例.分析結果表明：主動控制能夠有效地減小結構的地震響應，特別是當結構剛度較大時，3種主動控制算法控制效果均較理想；而當結構剛度較小時，各種控制算法的控制效果差異較大，LQR控制算法和模態控制算法的控制效果較好，而極點配置控制算法控制效果較差，甚至可能放大地震響應，需要優化極點配置得到適宜的狀態反饋增益矩陣以改善結構加速度響應的控制效果.

主動控制算法；結構剛度；控制效果；Matlab

0　引言

結構振動控制是振動工程領域重要的研究方向，根據控制原理不同，可分為主動振動控制、半主動振動控制和被動振動控制［1-5］.主動控制能夠根據傳感元件檢測到被控對象的振動信息，采用主動控制算法，通過作動器對被控對象施加控制力，從而減小或抑制被控對象的振動反應［1-2］.目前，主動控制相關技術方法已在航空航天領域、土木工程領域、交通運輸領域及機械工程領域的應用中初具規模［6］.與被動振動控制方法相比，主動控制方法具有更好的適應性［7-9］.目前，結構主動控制算法主要有：LQR控制算法、模態控制算法、極點配置控制算法、自適應控制算法、魯棒控制算法和智能控制算法等.LQR控制算法通過設定控制參數獲得一定的控制效果；模態控制算法利用模態坐標變化，將系統的運動方程完全解耦，將整個結構受控系統振動轉化為對各階模態控制，大大降低控制復雜度［10］；極點配置控制算法使系統的閉環極點配置到期望位置上［11-12］；自適應控制算法所采用的控制策略大致有自適應前饋控制、自校正控制和模型參考自適應控制3大類；魯棒控制算法是選擇線性反饋律使得閉環系統的穩定性或性能對于擾動具有一定的抵抗能力；模糊控制作為智能控制的一個重要分支，為解決不易或無法建模的復雜系統控制問題提供了有力手段.文中采用3種工程中最常采用的主動控制算法對不同剛度的結構振動控制效果進行了比較分析，并得出相應的結論.

1　結構振動控制系統模型

已知結構質量矩陣M、剛度矩陣K和阻尼矩陣C，其中，結構阻尼矩陣可按Rayleigh阻尼由前兩階振型阻尼比確定，即C=αM+βK.為元素均為1的列向量.采用結構層相對地面位移坐標空間建立結構運動方程，即：

2　LQR控制算法

在采用LQR控制算法時，權矩陣Q和R是2個重要的參數，但根據現有研究成果，只能依據工程經驗先確定Q和R，使其達到一定的控制效果.

設Q和R分別為：

式（3）中，α和β是待定系數.

由Matlab函數lqr可得連續狀態方程的控制力反饋增益矩陣，即：

則最優控制力：

將式（4）求解得到的狀態反饋增益矩陣G代入式（2），得受控結構狀態方程：

由Matlab函數lsim可得結構控制系統狀態反應：

式（7）中，C0為單位矩陣，D0為零矩陣，t為地震作用時間向量，y是輸出向量，它與狀態向量Z相同.

3　模態控制算法

設結構反應

由于廣義質量矩陣、剛度矩陣和阻尼矩陣均是對角矩陣，所以，式（9）可表示成n個獨立的運動方程，再將其表示為下述獨立的狀態方程：

針對每一個模態坐標方程，權矩陣

由Matlab函數lqr得到每一獨立廣義模態坐標的反饋增益矩陣：

將其組合形成廣義模態坐標狀態向量反饋增益矩陣Gc=［Gc1Gc2］.則控制力表達式：

結構控制系統反應仍利用式（7）計算.

4　極點配置控制算法

由Matlab函數eig可得到結構系統矩陣A的特征值向量，即：

期望極點可以根據減振要求的阻尼比確定.

假設G的選擇使下式成立：

式（16），式（17）中，Φj（λi）是Φ（λi）的第j列，ej是單位矩陣Ip的第j列，當（λiIn-A）-1存在時，rank［（λiIn-A）-1］=p.

如果λi（i=1，2，…，n）是不同的特征值，則狀態反饋增益矩陣可表示為：

式（18）中，Γ=［Φj1（λ1）Φj2（λ2）…Φjn（λn）］∈Rn×n，e=［ej1ej2…ej3］∈Rp×n.

將式（18）得出的狀態反饋增益矩陣G代入式（2），由Matlab函數eig得到閉環控制系統的極點為

與設定的極點相同.

結構控制系統反應仍利用式（7）計算.

5　算例

某3層框架結構，設結構層質量為mi=4×105kg，i=1，2，3［5］，假定結構前二階振型阻尼比ζ1=ζ2=0.05.結構的激勵為EI Centro（NS，1940）地震波，地震輸入峰值為200 gal.為比較各種不同控制算法在不同類型結構（剛性和柔性結構）中的控制效果，結構的層間抗側剛度ki分別取不同的值，采用上述3種主動控制算法計算結構控制系統反應及控制力.為衡量結構各種控制算法的減震效率，定義別為無控及施加控制力時的地震反應，δi為第i層減震效率，δi越大，減震效率越高.

對于普通結構，可取ki=p×107N／m，p=2，…，10［5］.結構振動響應見圖1～圖2（只給出結構抗側剛度ki=7×107N／m時，第3層的位移、加速度，其他層類似），減震效率見圖3～圖4.

通過計算分析，從圖1～圖2中看出，無控時結構的第3層位移和加速度比用3種主動控制算法控制時的第3層位移和加速度都要大（其他層類似）.從圖3～圖4中可看出，主動控制能夠顯著地減小結構的地震響應.

圖1　3種主動控制算法和無控時第3層層間位移Fig.1 Three active control algorithms and the third layer displacement without control（ki=7×107N·m-1）

圖2　3種主動控制算法和無控時第3層加速度Fig.2 Three active control algorithms and the third layer acceleration without control（ki=7×107N·m-1）

圖3　3種主動控制算法對結構總加速度的減震效率Fig.3 Three active control algorithm s on the total acceleration of the structure of the dam ping efficiency

圖4　3種主動控制算對結構總位移的減震效率Fig.4 Three active control algorithm s on the total displacement of the structure of the dam ping efficiency

對于剛性結構，可取ki=p×108N／m，p=2，…，10［5］.結構振動響應見圖5～圖6（只給出結構抗側剛度ki=7×108N／m時，第3層的位移、加速度，其他層類似），減震效率見圖7～圖8.從圖中可看出，當結構剛度較大時，3種主動控制算法控制效果均較理想，甚至基本一致.

對于柔性結構，可取ki=p×106N／m，p=2，…，10［5］.結構振動響應見圖9～圖10（只給出結構抗側剛度ki=7×106N／m時，第3層的位移、加速度，其他層類似），其中極點配置控制算法中極點配置方法按本文式（19）確定，減震效率見圖11～圖14.

從圖9～圖14中得出，當結構較柔、剛度較小時，3種控制算法對結構位移響應均具有良好的控制效果，對于加速度響應，LQR控制算法和模態控制算法仍具備一定的控制效果，而極點配置控制算法使結構第3層加速度的減震效率為負值，說明當結構剛度較小時，極點配置控制算法非但不能減小結構地震響應，甚至會放大結構的加速度響應，且總位移減震效率也明顯低于其他2種主動控制算法；因此，對柔性結構來說，采用式（19）極點配置方法的極點配置控制算法并不適宜，而對于結構振動控制來說，一般對加速度的控制比對位移的控制更加迫切，此時，極點配置控制算法中采用文獻［13］提出的極點優化配置方法，得到適宜的狀態反饋增益矩陣，可以顯著改善結構加速度響應控制效果，如圖15～圖16所示.

圖5　有控和無控時剛性結構第3層層間位移Fig.5 The third layer displacement of rigid structure with control and without control（ki=7×108N·m-1）

圖6　有控和無控時剛性結構第3層加速度Fig.6 The third layer acceleration of rigid structure with control and without control（ki=7×108N·m-1）

圖7　有控時剛性結構總加速度的減震效率Fig.7 Dam ping efficiency of the total acceleration of the rigid structure with control

圖8　有控時剛性結構總位移的減震效率Fig.8 Dam ping efficiency of the total displacement of the rigid structure with control

圖9　有控和無控時柔性結構第3層層間位移Fig.9 The third layer displacement of flexible structure with control and without control（ki=7×106N·m-1）

圖10　有控和無控時柔性結構第3層加速度Fig.10 The third layer acceleration of rigid structure with control and without control（ki=7×106N·m-1）

圖11　有控時柔性結構第3層層間位移的減震效率Fig.11 Damping efficiency of the third layer displacement of flexible structure with control

圖12　有控時柔性結構第3層加速度的減震效率Fig.12 Damping efficiency of the third layer acceleration of flexible structure with control

圖13　有控時柔性結構總加速度的減震效率Fig.13 Dam ping efficiency of the total acceleration of the flexible structure with control

圖15　有控時柔性結構總加速度的減震效率Fig.15 Dam ping efficiency of the total acceleration of the flexible structure with control

圖16　有控時柔性結構第3層加速度的減震效率Fig.16 Damping efficiency of the third layer acceleration of flexible structure with control

6　結論

本文采用3種主動控制算法，計算了結構振動控制系統反應，比較分析了3種主動控制算法的振動控制效果，得到以下結論：

1）主動控制能夠有效地減小結構的地震響應，當結構剛度較大時，3種主動控制算法計算的控制效果均較理想，甚至基本一致；

2）當結構剛度較小、為柔性結構時，3種主動控制算法對結構的控制效果并不相同.整體而言，采用LQR控制算法和模態控制算法得到的結構振動響應較小，減震效率較高，控制效果較好，而極點配置控制算法減震效率明顯低于其他2種主動控制算法，甚至會放大結構加速度響應，這時需要優化極點配置，以得到適宜的狀態反饋增益矩陣，當采用文獻［13］所述的極點配置方法時，可以使結構加速度響應控制效果得到顯著改善.

［1］PREUMONTA.Vibration Control of Active Structures：A n Introduction［M］.Berlin：Springer Science＆Business Media，2012.

［2］HANSENC，SNYDERS，QIUX，et al.Active Control of Noise and Vibration［M］.Florida：CRC Press，2012.

［3］邱賢波.阻尼器減震結構隨機地震響應及等效阻尼比分析［D］.柳州：廣西工學院，2012.

［4］葛新廣，李創第，陸運軍.帶懸吊質量擺減振系統的建筑結構隨機風振響應的解析解法［J］.廣西工學院學報，2009，20（4）：6－10.

［5］歐進萍.結構振動控制——主動、半主動和智能控制［M］.北京：科學出版社，2003.

［6］CASCIATIF，RODELLAR J，YILDIRIM U.Active and Semi－Active Control of Structures-Theory and Applications：A Review of Recent Advances［J］.Journal of Intelligent Material Systems and Structures，2012，23（11）：1181－1195.

［7］張春良.振動主動控制及應用［M］.哈爾濱：哈爾濱工業大學出版社，2011.

［8］KORKMAZS.A Review of Active Structural Control：Challenges for Engineering Informatics［J］.Computers and Structures，2011，89（23）：2113－2132.

［9］EL－KHOURY O，ADELI H.Recent Advances on Vibration Control of Structures under Dynamic Loading［J］.Archives of Computational Methods in Engineering，2013，20（4）：353－360.

［10］侯軍芳，李榮利，白鴻柏.主動隔振平臺模態區間參數模型的獨立模態控制方法［J］.振動工程學報，2011，24（6）：670－675.

［11］MOTTERSHEAD JE，T EHRANI M G，J AMES S，et al.Active Vibration Suppression by Pole－Z ero Placement Using Measured Receptances［J］.Journal of Sound and Vibration，2008，311（3）：1391－1408.

［12］ZHANG J，OUYANG H，YANG J.Partial Eigenstructure Assignment for Undamped Vibration Systems Using Acceleration and Displacement Feedback［J］.Journal of Sound and Vibration，2014，333（1）：1－12.

［13］杜永峰，李春鋒，李慧.結構主動控制最優極點配置算法研究［J］.振動與沖擊，2012，31（19）：95－98，103.

（學科編輯：黎婭）

Research on the structural vibration control effects of three active control algorithms

FAN Wen-hua,ZOU Wan-jie*
(School of Civil Engineering and Architecture,Guangxi University of Science and Technology, Liuzhou 545006,China)

According to the different types of structure,the corresponding control force state feedback gain matrix is obtained by the correlation function commands in Matlab,then the displacement,acceleration and control force of the structure vibration are obtained,and the control effects of the three active control algorithms are compared.And given a calculation example of three layer frame structure based on Matlab,the results show that the active control can effectively reduce the seismic response of the structure.Especially when the structure stiffness is large,the control effects of the three active control algorithms are ideal.When the structure stiffness is small,the control effects of various control algorithms are relatively large.The control effects of LQR control algorithm and Modal control algorithm are better,and the pole placement algorithm is less effective,and it may even amplify seismic response.In this case,we need to optimize the pole placement to obtain an adequate state feedback gain matrix to improve the structural acceleration control effects.

active control algorithm;structural stiffness;control effect;Matlab

TU31

A

2095-7335（2016）03-0069-07

10.16375/j.cnki.cn45-1395/t.2016.03.013

2015-12-31

廣西自然科學基金（2011GXNSFA018032）資助.

＊鄒萬杰，博士，副教授，研究方向：結構振動控制，E-mail：705450818＠qq.com.