AMD系統(tǒng)主動控制控制器設(shè)計方法研究

劉 泳 伶

(四川省建筑科學(xué)研究院，四川 成都 610000)

AMD系統(tǒng)主動控制控制器設(shè)計方法研究

劉 泳 伶

(四川省建筑科學(xué)研究院，四川 成都 610000)

介紹了AMD主動控制系統(tǒng)的減震原理，基于經(jīng)典線性最優(yōu)控制算法，采用Wilson-θ法來求位移時程曲線，并利用matlab編程對一個框架剪切模型計算其在地震作用下采用主動控制前后的彈性時程反應(yīng)分析，同時進(jìn)行了對比，表現(xiàn)出了主動控制的顯著效果。

AMD，主動控制，時程分析，matlab，經(jīng)典線性最優(yōu)算法

振動主動控制系統(tǒng)主要由作動器、傳感器、控制器三部分組成，其中控制器是主動控制的核心環(huán)節(jié)，由它實現(xiàn)所需的控制律，也就是輸入與輸出的關(guān)系。因此要實現(xiàn)主動控制系統(tǒng)的最優(yōu)的控制律就需要掌握主動控制系統(tǒng)控制器的設(shè)計方法，而經(jīng)典線性最優(yōu)控制算法是目前應(yīng)用較為廣泛的一種算法。本文對經(jīng)典線性最優(yōu)控制算法進(jìn)行了介紹，利用它分析計算了采用主動質(zhì)量阻尼器AMD系統(tǒng)后的框架結(jié)構(gòu)的減震效果。

1 AMD控制系統(tǒng)的減震機理

假定將標(biāo)準(zhǔn)AMD系統(tǒng)放置在結(jié)構(gòu)頂部，圖1和圖2分別為AMD控制系統(tǒng)的計算模型與受力分析，受控結(jié)構(gòu)的運動方程如下:

(1)

(2)

(3)

當(dāng)AMD主動控制結(jié)構(gòu)體系在地震激勵下產(chǎn)生地震反應(yīng)時，AMD系統(tǒng)的驅(qū)動器驅(qū)動質(zhì)量塊，使得質(zhì)量塊產(chǎn)生運動，而從式(3)中可以看出，AMD系統(tǒng)對結(jié)構(gòu)施加的控制力u(t)與AMD系統(tǒng)質(zhì)量塊的慣性力相等，它通過彈簧、驅(qū)動器和阻尼器、作用在結(jié)構(gòu)上，通過選擇合理的控制算法確定最優(yōu)控制力，進(jìn)而達(dá)到控制和衰減結(jié)構(gòu)地震反應(yīng)的目的。

2 算例分析

如圖3所示一五層剪切型框架結(jié)構(gòu)，已知結(jié)構(gòu)的阻尼比為ξ=0.05；各樓層的質(zhì)量分別為m1=m2=m3=m5=3.5×105kg，m4=3.1×105kg；樓層的剛度分別為k1=k2=k3=k5=4.3×107N/m，k4=3.6×107N/m，本算例只在結(jié)構(gòu)底層施加一個作動器，選擇200 gal EL centro波作為激勵，利用matlab編制程序用Wilson-θ法對結(jié)構(gòu)進(jìn)行未加主動控制時的彈性時程反應(yīng)分析，然后用經(jīng)典線性最優(yōu)算法計算主動控制系統(tǒng)對結(jié)構(gòu)施加的最優(yōu)控制力，最后對比計算主動控制前后結(jié)構(gòu)在地震激勵下的位移反應(yīng)。

2.1 無主動控制時的結(jié)構(gòu)位移時程反應(yīng)

在地面運動激勵下結(jié)構(gòu)的運動方程為：

(4)

2.2 主動控制后的結(jié)構(gòu)位移時程反應(yīng)

2.2.1 運動方程

n個自由度系統(tǒng)在地面激勵下的運動方程如式(4)所示，進(jìn)行AMD主動控制后，其運動方程可表示為：

(5)

其中，BS為作動器定位矩陣；U(t)為作動器控制力向量。此算例中只在底層安裝作動器，那么BS=[1,0,0,0,0]T。

2.2.2 經(jīng)典線性最優(yōu)算法

(6)

其中，q為作動器個數(shù)；n為結(jié)構(gòu)層數(shù)。經(jīng)典線性最優(yōu)控制的二次型性能表示為：

其中，ts為地震持續(xù)時間；Q為狀態(tài)向量加權(quán)矩陣；R為控制力向量加權(quán)矩陣。 經(jīng)典線性最優(yōu)控制的目標(biāo)是在滿足運動方程(5)和狀態(tài)方程(6)的約束條件下，使目標(biāo)函數(shù)J最小，從而尋求最優(yōu)控制力U。引入Hamitonian函數(shù)，根據(jù)極值條件，可求得經(jīng)典線性最優(yōu)控制力：

U(t)=-Gz(t)

(7)

其中，G為反饋增益矩陣，G=αR-1BTP，其中P為下列Riccati方程的解：

-PA-ATP+PBR-1BTP-Q=0。

matlab內(nèi)置函數(shù)care函數(shù)提供了求解Riccati方程的方法,[P,L,G,RR]=CARE(A,B,Q,R),求出的P為Riccati方程的解，G為經(jīng)典線性最優(yōu)控制的增益矩陣。對于權(quán)矩陣Q，R選取目前還是個難題，只能通過試算的方式確定。權(quán)矩陣用于調(diào)整結(jié)構(gòu)反應(yīng)與控制力兩者之間的相對重要程度，在主動控制算法中決定了控制力和結(jié)構(gòu)反應(yīng)的大小，Q∈R2n×2n，R∈Rq×q，q為作動器的數(shù)量，而且Q和R分別是半正定矩陣和正定矩陣，最簡單的取法就是均采用對角矩陣，由于在此算例分析時只加了一個作動器，所以R是一個一階矩陣，即一個數(shù)?；赒,R矩陣上述特點，本例經(jīng)典線性最優(yōu)控制算法取Q=diag[6,6,6,6,6,6,6,6,6,6]×1013,R=2.9。在matlab中將A,B,Q,R代入care函數(shù)最終可得到增益矩陣G。

2.2.3 位移時程反應(yīng)及最優(yōu)控制力的求解

將求得的增益矩陣G代入式(4)及式(6)可得到如下的形式：

(8)

再將狀態(tài)變量z(t)代入式(8)可得到：

(9)

將G通過矩陣變換為：

(10)

(11)

3 結(jié)語

圖4～圖6分別為在結(jié)構(gòu)底層加作動器時，受控和無控結(jié)構(gòu)各層位移時程對比曲線以及底層控制力時程曲線(直線為無控時，*線為控制后)。從圖中可以明顯地看出AMD主動控制系統(tǒng)大大減小了結(jié)構(gòu)在地面激勵下的各層位移，通過有效施加控制力，起到了良好的減震效果。

[1] 周星德，姜冬菊.結(jié)構(gòu)振動主動控制[M].北京：科學(xué)出版社，2009.

[2] 周福霖.工程減震控制[M].北京：地震出版社，1997.

[3] 王小侃，馮冬青.基于matlab的LQR控制器設(shè)計方法研究[J].微計算機信息，2008,24(10)：95-97.

[4] 杜 炎.利用matlab進(jìn)行主動控制分析[J].中國科技信息，2008(2)：110-112.

[5] 徐趙東.matlab語言在建筑抗震工程中的應(yīng)用[M].北京：科學(xué)出版社，2010.

[6]AnilK.Chopra.結(jié)構(gòu)動力學(xué)理論及其在地震工程中的應(yīng)用[M].北京：高等教育出版社，2012.

[7] 徐金明.matlab實用教程[M].北京：清華大學(xué)出版社，2011.

Research on AMD system active control controller design method

Liu Yongling

(SichuanBuildingScienceResearchInstitute,Chengdu610000,China)

This paper introduced the shock absorption principle of AMD active control system, based on classical linear optimal control algorithm, using Wilson-θmethod to get the displacement time history curve, and using the matlab programming analyzed the elastic time history response to a frame shear model calculation under earthquake using active control before and after, and has carried on the contrast, showed the significant effect of active control.

AMD, active control, time history analysis, matlab, classical linear optimal algorithm

2015-01-09

劉泳伶(1987- )，女，碩士，助理工程師

1009-6825(2015)08-0065-02

TU318

A