基于市場機制的控制策略在磁流變阻尼器上的應用

半主動裝置的應用在上部結構之間配置隔震層，可以有效減少地震能量向上部結構的轉移。作為一種行之有效的結構被動控制方法,它改變了傳統的依賴增加結構強度、剛度和延性的抗震設計。在隔震結構中以橡膠隔震支座為代表的隔震技術發(fā)展最為成熟且應用最廣泛。盡管應用普通橡膠隔震支座可以延長結構自振周期,避開地震作用優(yōu)勢周期的分布范圍,減小上部結構遭受的地震作用,但將導致隔震層產生過大的位移。尤其在強震作用下,變形超過規(guī)定范圍影響正常使用甚至造成碰撞而破壞。為防止隔震層發(fā)生過度變形,在隔震結構中設置阻尼器是必不可少的,如:滯回衰減型阻尼器、黏性衰減型阻尼器等。然而,用被動式阻尼器來減小隔震層的位移是以增加上部結構的變形和加速度為代價的,而且被動式阻尼器特性固定,無法保證在任意隨機的環(huán)境激勵下都能獲得良好的控制效果。因此,學者們提出在隔震層附加半主動裝置,以結構的半主動隔震來實現有效限制隔震層位移同時確保上部結構振動反應在規(guī)定的范圍。半主動裝置中尤以磁流變阻尼器(MagnetorheologicalDamper,MRD)為代表。它是一種性能優(yōu)良的可控智能式變阻尼控制裝置,通過調節(jié)磁場強度可以迅速地改變磁流變液的力學阻尼特性,具有能耗小、阻尼力可調范圍寬、響應速度快等優(yōu)點。將MRD與普通橡膠隔震支座結合組成的智能隔震系統,不僅能夠減小上部結構的地震反應,而且能夠有效地抑制隔震層的大變形,是一種性能優(yōu)異的混合控制系統。同時學者們應用不同的控制算法,如瞬時最優(yōu)控制算法、神經網絡控制算法、模糊控制算法等等,對此控制系統進行了數值分析?；谑袌鰴C制的控制(Market-BasedControl,MBC)策略應用自由經濟中的市場競爭機制模擬復雜的控制系統,它充分利用市場的價格機制來分配控制系統中有限的控制能量,可以看成是一種在虛擬的市場環(huán)境中完成控制能量最優(yōu)分配的策略。與之相應的算法即為MBC算法。本文將MBC理論應用于MRD-隔震結構的半主動控制中,提出了相應的控制策略,實現對MRD電壓的實時控制,從而保證控制力的實時最優(yōu)調節(jié)。通過一個五層隔震結構的數值分析驗證這種方案的有效性,顯示此策略能很好地應用于MRD-隔震結構的半主動控制中。1mrd-flun結構模型1.1t、kxn具有n個自由度MRD-隔震結構(如圖1所示,n=6),在地震作用和MRD阻尼力作用下的運動方程可表示為:Μ¨x(t)+C˙x(t)+Κx(t)=-Μ?！g(t)-BsfΜRD(t)(1)Mx¨(t)+Cx˙(t)+Kx(t)=?MΓx¨g(t)?BsfMRD(t)(1)式中,M、C和K分別為系統的質量、阻尼和剛度矩陣;x、˙xx˙和¨xx¨分別為位移、速度和加速度向量;¨xx¨g為地面加速度,Γ為n維單位列向量;fMRD為MRD提供的阻尼力,Bs為相應的位置矩陣。將運動方程(1)寫成狀態(tài)方程形式:˙Ζ(t)=AΖ(t)+BfΜRD(t)+D¨xg(t)(2)Z˙(t)=AZ(t)+BfMRD(t)+Dx¨g(t)(2)其中,Ζ(t)={x(t)¨x(t)}Z(t)={x(t)x¨(t)};A=[0Ι-Μ-1Κ-Μ-1C]A=[0?M?1KI?M?1C];B=[0-Μ-1Bs]B=[0?M?1Bs];D=[0-Γ]D=[0?Γ]。1.2黏滯力及bouc-wn滯變模型MRD所用的磁流變流體在外磁場作用下能在幾毫秒時間內完成從流動性能良好的線性黏滯流體狀態(tài)到具有一定剪切屈服強度的半固體狀態(tài)的轉變。目前,國內外學者針對MRD提出了不同的力學模型。本文采用修正的Bouc-Wen模型,即去掉了Bouc-Wen模型中的線性彈簧單元,認為MRD產生的阻尼力是黏滯力和Bouc-Wen滯變阻尼力之和,其阻尼力表示為如下形式:fΜRD=c0¨xb+αzΜRD(3)[ΗS1*2/3]˙zΜRD=-γ|˙x|zΜRD|˙z|n-1-β˙xb|˙zΜRD|n+A˙xb(4)fMRD=c0x¨b+αzMRD(3)[HS1*2/3]z˙MRD=?γ|x˙|zMRD|z˙|n?1?βx˙b|z˙MRD|n+Ax˙b(4)式中:¨xx¨b為隔震層的速度,zMRD為MRD的滯回變量;γ、β、n和A為控制卸載時力與速度曲線的線性性狀和層服前至層服后漸變段光滑度的參數;而參數c0與α為與電壓u線性相關的量,可表示為:c0=c0a+c0bu(5)α=αa+αbu(6)c0=c0a+c0bu(5)α=αa+αbu(6)其中u為產生相應阻尼力時MRD上的電壓值。由于受到回路電阻及電感等因素的影響,電壓u會滯后于根據控制算法得出的線圈輸入電壓值υ,這種滯后關系可用一階微分方程表示:˙u=-η(u-υ)(7)u˙=?η(u?υ)(7)2md結構模型對MRD進行半主動控制時,首先應由主動控制算法,如線性二次型經典最優(yōu)控制算法(LQR)等,計算出結構所需要的最優(yōu)控制力,然后求得需要加在MRD上的電壓大小,進而使其產生趨于主動最優(yōu)控制力的阻尼力。本文采用MBC策略結合修正的Clipped-optimal控制算法完成MRD線圈輸入電壓值υ的實時控制。2.1基于權變的md-隔震結構的能量分配算法經濟學(economics)研究的是在資源稀缺的給定條件下,消費者、生產廠家和政府部門為實現各自目標而進行的選擇。在不需政府干預的自由市場經濟中,作為市場主體的銷售商(seller,賣方)和消費者(consumer,買方)在價格準則作用下分配稀缺資源來實現各自利益的最大化。消費者與銷售商之間的買賣行為都是在一定的價格條件下進行的,而市場的總供給和總需求相互影響決定著這個價格。結構控制系統與自由市場經濟體系之間具有一定的相似性,兩者都要實現稀缺資源的合理有效分配。在自由市場經濟中,銷售商擁有的商品是待分配的稀缺資源,而在結構控制系統中,如MRD-隔震系統中,系統所需的控制能量(以電壓形式促成MRD出力)是待分配的稀缺資源。Lynch等將MBC理論應用到結構振動控制領域,其研究表明采用MBC算法對結構進行振動控制可以獲得與基于LQR算法相當的控制效果?；袅稚葢肕BC算法對半主動的TLCD進行控制也表明采用此方法不僅減振效果好,而且節(jié)約能量。用市場機制來模擬結構控制系統時,虛擬市場中的供需函數并沒有嚴格的形式。在選取供給函數時,主要考慮的是市場價格,價格越高,供應量越大。而在選取需求函數時,主要考慮的是結構的反應(如位移與速度)和控制能量的價格,即:結構反應越大,為減小結構振動所需購買的能量越大;價格上漲,結構趨向于購買較少的能量。在每一個控制時間段內市場中的總供需達到均衡,形成控制能量的均衡價格(如圖2所示),系統在此價格下實現利潤與效用的最大化。為提高計算效率,在MBC策略中,供給函數一般取為線性函數形式而需求函數模型可取不同的形式,如文獻和需求函數采用了線性形式而文獻和采用了冪函數形式。本文需求函數采用另一種形式即指數函數形式,其優(yōu)點在于市場中能量價格為零時,結構系統只能得到有限的能量(即可能獲得的最大能量),防止虛擬市場中價格較小時控制能量的“過度”購買,且該形式也符合經濟學中需求函數的一般形式。由于在MRD-隔震結構中,隔震層的剛度一般都較小。當結構遭受地震作用時,由隔震層來承受地震動引起的位移運動,而上部結構可看作近似平動。附加MRD的作用在于消耗地震能量,對上部結構的整體平動影響較小。應用MBC策略時,鑒于MRD-隔震結構這種“平動”的特點,采用隔震層的反應來確定該控制系統所需的控制能量,即有:供給函數:QS=ζ·p(8)需求函數:QD=|μ?xb+λ?˙xb|?e-cp(9)式中,ζ為反映能量源供給的常數;xb和˙xb分別表示隔震層相對地面的位移與速度,μ≥0,λ>0為相應的權系數;c為需求調節(jié)系數,可取為1。當市場達到均衡時:Qs=QD(10)可求得均衡價格p。則可求得系統所受的最優(yōu)控制力,即正比于以均衡價格購得的能量:fc=Κ?(μ?xb+λ?˙xb)?e-cp(11)其中,K>0為控制力增益系數,可按MRD的出力界限及所要達到的控制效果進行選取。2.2智能加壓ld控制對MRD控制系統,其半主動控制律即為加壓方式。加壓方式問題是關系到MRD能否充分體現其連續(xù)可調等優(yōu)點的關鍵,直接影響到系統的控制效果。到目前為止,國內外學者對MRD提出了許多加壓方式,如恒定加壓式、Heaviside函數式、動力逆模型加壓、智能加壓式等等。本文采用修正的Clipped-optimal控制算法,使施加于MRD上的電壓可在0與最大電壓Vmax之間變化。該方法形式簡單而效果顯著。υ=VcΗ({fc-fΜRD}?fΜRD)(12)Vc={φfcVmax,for|fc|≤fΜRD,max,for|fc|>fΜRD,max(13)其中:H(·)為Heaviside階躍函數;fMRD,max為MRD所能產生的最大阻尼力;φ=Vmax/fMRD,max,為電壓與力的關系系數。3計算值的例和分析3.1n4d種加速度以圖1所示的剪切型MRD-隔震結構為例,模型基本參數參見文獻。隔震層處附加的MRD可實現最大阻力為50kN,具體參數參見文獻。施加電壓變化范圍為0～5V,MRD阻尼出力可調倍數約為24。時程分析時選取三條最不利地震動:F4(ElCentro,N69W,1979)、F5(Taft,N21E,1952)和N2(Gengma,S00E,1988),加速度峰值統一調整為400Gal,相當于烈度為8度地區(qū)大震情況。進行三種工況分析,即:非隔震(No-Iso.)、MRD-隔震LQR半主動(LQR-Semi)和MRD-隔震MBC半主動(MBC-Semi)。3.2lqr控制策略表1列出了各種工況下結構反應的峰值。從結果可以看出,對于本算例,在保證隔震層相對地面位移反應峰值相當的情況下,應用本文提出的MBC策略對上部結構層間位移的控制效果要略好于基于經典最優(yōu)LQR控制策略;同時,從結構的加速度峰值反應可以看到,MBC策略要明顯好于LQR控制策略。圖3給出了隔震層相對地面位移與頂層加速度反應的時程曲線,也可知在控制時域內兩種策略在控制隔震層相對地面位移相當的情況下,MBC策略在結構加速度方面表現出的控制效果要好于LQR控制策略。圖4給出了各工況下的MRD出力時程曲線(其為連接各仿真計算時刻點上MRD出力值所得到,該連續(xù)曲線上的穿零點僅意味著其兩側計算時刻點上MRD出力方向發(fā)生變化,并不代表該時刻MRD出力為零,實際MRD出力也不可能達到零值),從時程曲線上看在控制時域內MBC策略控制下MRD總出力要比LQR策略下的小。這意味著要達到相當的控制目標,MBC策略所需的控制能量輸入要小于LQR控制策略的能量輸入。由于MRD的變阻尼特性,它提供的控制力將以阻尼力的形式體現,且該阻尼力與MRD的相對速度有近似線性的關系。圖5為各地震動輸入下,MRD在LQR算法與MBC算法下的出力與相對速度之間的關系圖?？芍狹BC算法比LQR算法能更好地反映出MRD出力的這種阻尼特性。從算法本身來看,LQR算法是基于系統全狀態(tài)反饋的最優(yōu)控制算法,算法的有效性取決于全部狀態(tài)能否實時準確反饋。當結構自由度較多時,在線實時計算往往較長,容易出現時滯問題。而MBC算法所需應用的信息僅僅是隔震層的信息,相當于對隔震層反應進行靜態(tài)離散的Pareto最優(yōu)控制,進而實現全結構系統的最優(yōu)控制。其所需的信息量少,實時計算時間短,可以很好地消除時滯影響。在對MRD-隔震結構多維控制情況下,可將本文提出的MBC算法稍加改進,用兩個相互獨立的供需方程來分別計算X和Y方向的最優(yōu)控制力需求,便能實現對偏心MRD-隔震結構的有效控制。