基于可控性柱面網(wǎng)殼結(jié)構(gòu)主動控制作動器布置的優(yōu)化研究

王社良，田秀婷，代建波，趙 祥

(西安建筑科技大學(xué) 土木工程學(xué)院，西安 710055)

隨著建筑對跨度、高度和結(jié)構(gòu)形式的要求，網(wǎng)殼結(jié)構(gòu)以其突出的優(yōu)點(diǎn)已廣泛應(yīng)用于各種重要的標(biāo)志性建筑。然而，這類結(jié)構(gòu)通常柔性大、阻尼小、自振頻率較低，在動力荷載(如地震和強(qiáng)風(fēng))作用下會發(fā)生較大的振動變形、構(gòu)件破壞及一定程度的損傷積累。因此有關(guān)網(wǎng)殼結(jié)構(gòu)的抗震及減震問題更顯得突出。傳統(tǒng)的抗震設(shè)計方法依靠結(jié)構(gòu)自身的耗能能力來吸收地震能量，缺乏自我調(diào)節(jié)能力，而主動控制技術(shù)具有智能性、靈活性及減震效果明顯等優(yōu)點(diǎn)。該技術(shù)在機(jī)械、航空航天和船舶等領(lǐng)域已經(jīng)取得了廣泛應(yīng)用，而在土木工程中的研究則剛剛起步。設(shè)計適應(yīng)于土木工程結(jié)構(gòu)的主動控制作動器，如何合理地布置作動器以有效的抑制結(jié)構(gòu)振動，以及作動器參數(shù)對結(jié)構(gòu)不利因素的控制效果和影響是當(dāng)前的研究重點(diǎn)［1］。

作動器的性能在很大程度上決定了整個振動控制系統(tǒng)的性能，而作動器的性能主要取決于作動器采用何種驅(qū)動材料。超磁致伸縮材料(Giant Magnetostrictive Material，簡稱為GMM)是一種在外磁化狀態(tài)改變時，其尺寸會產(chǎn)生顯著的變化，去掉外磁場后，又恢復(fù)原來尺寸的智能材料，這類材料具有輸出力大、應(yīng)變顯著、工作頻帶寬、響應(yīng)速度快和驅(qū)動電壓低等優(yōu)點(diǎn)。利用GMM性能研制的作動器定位精度高、響應(yīng)速度快、輸出力大、設(shè)計相對簡單［2］，因此在土木工程主動控制中的應(yīng)用具有良好前景。

1 超磁致伸縮作動器結(jié)構(gòu)設(shè)計及性能

1.1 設(shè)計原理

GMM作動器就是通過外加磁場使GMM元件發(fā)生形變從而產(chǎn)生作動力與位移。勵磁繞組提供磁場，磁場的大小可以通過調(diào)節(jié)電磁鐵勵磁繞組的電流控制。外加磁場的方向應(yīng)與GMM元件產(chǎn)生變形的方向一致，當(dāng)GMM元件中的磁場達(dá)到一定值時，元件發(fā)生變形。磁場去掉以后，GMM元件基本自動恢復(fù)變形。

1.2 GMM作動器總體結(jié)構(gòu)設(shè)計

圖1是本次試驗(yàn)的GMM作動器的總體結(jié)構(gòu)示意圖與制造好的實(shí)物圖。此GMM作動器包括外套、探測線圈、偏置線圈、激勵線圈、線圈骨架、GMM元件、作動桿、連接桿、預(yù)壓碟簧和調(diào)節(jié)螺母等。通電后偏置線圈、激勵線圈提供疊加磁場，單層探測線圈用于測量磁場的大小。GMM元件在非工作狀態(tài)時就通過伸長方向上的碟簧加一足夠大的反力，可起到預(yù)緊及增加位移輸出量的作用。當(dāng)激勵線圈通以電流信號時產(chǎn)生磁場，GMM元件產(chǎn)生伸長變化且伸長效應(yīng)通過作動桿輸出給作動對象，斷掉電流信號后磁場消失，GMM元件恢復(fù)原來形狀，完成此作動器的作動效應(yīng)。

圖1 GMM作動桿總體結(jié)構(gòu)示意圖及實(shí)物圖Fig.1 GMM actuator general structure schematic and real object

1.3 GMM作動器輸出性能測試

通過電腦式伺服控制多功能力學(xué)試驗(yàn)機(jī)對制作的GMM作動桿進(jìn)行了輸出性能測試，如圖2和圖3所示。通過對不同預(yù)壓力作用下GMM作動桿的輸出特性研究發(fā)現(xiàn)，該 GMM做動桿在6MPa下能同時達(dá)到最大位移和最大力的輸出。位移輸出在0.2A ～ 1.0A 線性度較好，3.0A以后磁致伸縮量基本趨于飽和值。輸出力則在整體上呈線性增長。如圖4和圖 5所示。

2 基于遺傳算法的作動器位置優(yōu)化

在主動控制中作動器的結(jié)構(gòu)形式以及配置方案不僅對控制效果有著重要的影響，而且直接決定了結(jié)構(gòu)控制系統(tǒng)的穩(wěn)定性、可控性、以及控制的有效性。作動器位置合理布置可以在較小耗能情況下達(dá)到較好的控制效果;反之，不合理的布置會產(chǎn)生控制系統(tǒng)的硬件成本增加、復(fù)雜性增加、控制效果降低等弊端，甚至?xí)?dǎo)致系統(tǒng)不穩(wěn)定［3］。作動器的優(yōu)化位置本質(zhì)上是選擇合適的優(yōu)化算法使得某一指標(biāo)函數(shù)取得全局最優(yōu)［4］。

遺傳算法是模擬自然界中生物進(jìn)化過程發(fā)展起來的隨機(jī)全局搜索和優(yōu)化方法，是一種高效、并行、全局搜索的方法。它使用適者生存的原則，在潛在的解決方案種群中逐次產(chǎn)生一個近似最優(yōu)的方案［5］。遺傳算法主要包括優(yōu)化參數(shù)編碼、初始種群的確定和遺傳操作三個步驟。

2.1 柱面網(wǎng)殼結(jié)構(gòu)模型

桿件加入作動器后有限元動力模型為:

其中:Me，Ke，De分別為單元質(zhì)量矩陣、剛度矩陣、耦合剛度矩陣。

為了簡化推導(dǎo)柱面網(wǎng)殼結(jié)構(gòu)的動力學(xué)方程，做以下假設(shè):(1)結(jié)構(gòu)振動完全處于彈性范圍內(nèi)。(2)結(jié)構(gòu)的振動模態(tài)不發(fā)生不耦合。受控結(jié)構(gòu)的動力學(xué)方程可表示為［6］:

不考慮外擾力的作用，對于實(shí)際的物理模型，低階模態(tài)相比高階模態(tài)顯得更加重要，為降低模型維數(shù)我們做模態(tài)截斷處理，即保留低階模態(tài)，忽略高階模態(tài)。只考慮控制系統(tǒng)前nc階主模態(tài)。根據(jù)模態(tài)疊加原理:

其中:［φ］c={φ1，…，φnc}為前nc階控制模態(tài)的特征向量矩陣，{p}為模態(tài)坐標(biāo)向量。

其中:pi為第i個模態(tài)坐標(biāo)，ωi為第i階固有頻率，ξi為第i個模態(tài)的阻尼比，{fi}為模態(tài)控制力向量。

2.2 優(yōu)化準(zhǔn)則

如果作動器的布置位置可以保證產(chǎn)生最大的模態(tài)控制力，則可保證作動效率也就是可控性較高，這個位置就是最優(yōu)的位置。

令B=［φc］T［D］，由式(5)可以得到:

其中:A=BTB為正定對稱矩陣。

可以證明如下不等式成立［7］

從式(6)可以看出，如果A的所有特征值均集中分布且很大，則可以保證系統(tǒng)作動效率較高。這樣就可以定義一個作動器位置優(yōu)化準(zhǔn)則［8］。

優(yōu)化準(zhǔn)則乘積前一項(xiàng)考慮了A所有的特征值的算術(shù)平均值，后一項(xiàng)考慮了A所有特征值的幾何平均值，兩項(xiàng)乘積最大可以保證A所有特征值分布集中且最大。在提高作動效率的同時還可以保證各個作動器的控制作用更加均衡。

圖6 遺傳算法優(yōu)化流程圖Fig.6 Flow chart of genetic algorithm optimization

2.3 優(yōu)化步驟

柱面網(wǎng)殼結(jié)構(gòu)作動器位置遺傳算法優(yōu)化流程見圖6。

(1)編碼和初始種群:采用二進(jìn)制碼對個體編碼。用網(wǎng)架結(jié)構(gòu)的桿件總數(shù)r來代表作動器可能布置的位置，若基因值為1，表示作動器放置在相應(yīng)的位置，若基因值為0，表示相應(yīng)位置沒有布置作動器。通過隨機(jī)方法生成初始種群。

(2)計算適應(yīng)度:選擇個體目標(biāo)函數(shù)的值作為其適應(yīng)度，適應(yīng)度值的大小決定了個體生存的機(jī)會。

(3)選擇操作:是對群體中的個體進(jìn)行優(yōu)勝劣汰操作的過程。根據(jù)各個個體適應(yīng)度值，按照一定規(guī)則使適應(yīng)度較高的個體復(fù)制到下一代。這使群體中個體的適應(yīng)度值不斷接近最優(yōu)解。

(4)交叉操作:采用多斷點(diǎn)交叉法產(chǎn)生新的后代，每個基因?qū)?yīng)一個優(yōu)化位置信息。設(shè)交叉概率為pc，初始種群大小為pop，則隨機(jī)選擇出n=round(pop×pc/2)對染色體作為雙親。

(5)變異操作:交叉可以擴(kuò)大搜索空間，有效防止未成熟的收斂而陷入局部最優(yōu)解。設(shè)變異概率為pm，對每一個基因，隨機(jī)產(chǎn)生一個［0，1］之間均勻分布的數(shù)a，若a≥pm，則該基因進(jìn)行變異操作。

對于一般的遺傳算法而言，運(yùn)行參數(shù)的取值范圍為:初始種群的大小pop取10～100，交叉概率pc取0.4～0.99，變異概率pm取0.000 1 ～0.1。

3 算例分析

3.1 柱面網(wǎng)殼模型

一正放四角錐雙層柱面網(wǎng)殼，跨度B=38 m，長度L=40 m，矢跨比F/B=2/7，周邊固定。桿件鋼材為Q235，彈性模量為 2.1 ×106N/m2，泊松比為 0.3，密度為7 800 kg/m3。結(jié)構(gòu)共有57個節(jié)點(diǎn)，158根桿件，共12個支座，如圖7所示。

3.2 優(yōu)化模型

利用MATLAB軟件建立優(yōu)化模型，進(jìn)行遺傳算法進(jìn)行優(yōu)化時，只能求得適應(yīng)度函數(shù)的最小值，須對上面求得的性能指標(biāo)進(jìn)行調(diào)整，從而得到適應(yīng)度函數(shù)［9］:

圖7 柱面網(wǎng)殼有限元模型Fig.7 Cylindrical shells structure finite element model

其中:J0為調(diào)整前的適應(yīng)度值，p為當(dāng)前布置的作動器數(shù)目，m為預(yù)期布置的作動器數(shù)目數(shù)學(xué)模型:

3.3 作動器優(yōu)化

本算例分別選取 10、20、30、40、50、70、90、110、130、158個作動器情況，結(jié)合遺傳算法，分別進(jìn)行位置優(yōu)化。遺傳算法中的參數(shù)設(shè)置選擇如下:初始種群的大小pop取50個，最大代數(shù)取200代，交叉概率pc取0.8，變異概率pm取 0.01。

根據(jù)柱面網(wǎng)殼的特點(diǎn)和ANSYS的分析，在地震作用下結(jié)構(gòu)的最大位移發(fā)生在跨中附近。本算例在結(jié)構(gòu)中輸入的是一條400 gal的x方向EL-CENTRO波，根據(jù)柱面網(wǎng)殼的特點(diǎn)，在受到x方向地震波的時候，結(jié)構(gòu)中的節(jié)點(diǎn)x方向位移最大，z方向次之，故選取13節(jié)點(diǎn)x方向位移在不同作動器數(shù)目布置下的控制效果進(jìn)行對比，從而對作動器的數(shù)目進(jìn)行初步確定。

為了對比結(jié)構(gòu)在布置作動器和未布置作動器兩種情況下結(jié)構(gòu)在外部激勵作用下的反應(yīng)控制效果，引入一個控制效果系數(shù)β:

其中:X未為未布置作動器時的結(jié)構(gòu)反應(yīng)，X有為布置作動器時的結(jié)構(gòu)反應(yīng)。

圖8為結(jié)構(gòu)13節(jié)點(diǎn)x方向位移在不同數(shù)目作動器下控制效果的對比。從圖8可以看出:通過布置不同數(shù)目的作動器，對結(jié)構(gòu)中13節(jié)點(diǎn)均起到控制效果。隨著作動器數(shù)目增加，對13節(jié)點(diǎn)x方向的控制效果越好。作動器數(shù)目較少時，增加作動器數(shù)量對控制效果提高比較大，當(dāng)增加到一定數(shù)目，增加作動器數(shù)量對控制效果的提高較為緩慢。綜合考慮控制效果、硬件成本、系統(tǒng)復(fù)雜性以及穩(wěn)定性，結(jié)構(gòu)上布置40個作動器較為適宜。

選取布置40個作動器情況，結(jié)合遺傳算法進(jìn)行位置優(yōu)化設(shè)計，作動器位置優(yōu)化結(jié)果見表1位置優(yōu)化結(jié)果表明:由于柱面網(wǎng)殼結(jié)構(gòu)是完全對稱的，作動器的位置多對稱集中在結(jié)構(gòu)的兩端，也就是結(jié)構(gòu)各階振型變形最大的地方。從具體布置位置來看，作動器位置主要集中在斜腹桿件，其次是上、下弦桿件，并主要布置在上弦縱向桿件和下弦橫向桿件。

表1 40個作動器位置優(yōu)化結(jié)果Tab.1 location optimization results of 40 actuators

根據(jù)布置40個作動器的位置優(yōu)化結(jié)果，應(yīng)用LQR控制算法［10］，選取模型結(jié)構(gòu)13節(jié)點(diǎn)進(jìn)行控制與未控制的位移時程曲線對比分析，如圖9、10所示。

由圖9、10可知，加入作動器進(jìn)行主動控制后，柱面網(wǎng)殼結(jié)構(gòu)13節(jié)點(diǎn)x方向和z方向的位移有較大的控制，說明采用GMM作動器進(jìn)行柱面網(wǎng)殼結(jié)構(gòu)振動主動控制是可行的而且效果非常明顯。

4 結(jié)論

本文針對網(wǎng)殼結(jié)構(gòu)的振動控制，采用GMM材料為核心元件制作了主動控制作動器并進(jìn)行了輸出性能測試，最后通過算例分析了作動器在網(wǎng)殼結(jié)構(gòu)中的數(shù)量及位置優(yōu)化問題以及地震反應(yīng)控制效果，結(jié)論如下:

(1)根據(jù)GMM變形機(jī)理和磁控特性設(shè)計制作了GMM作動器。經(jīng)輸出性能測試可知此作動器非常適宜用于網(wǎng)殼結(jié)構(gòu)振動主動控制。

(2)分析了作動器數(shù)目與控制效果系數(shù)的關(guān)系，隨著作動器數(shù)目的增加，控制效果系數(shù)增加的速度減小，最后趨于平緩。利用控制效果系數(shù)對作動器數(shù)目進(jìn)行初步優(yōu)化，可以讓作動器對結(jié)構(gòu)的控制效果、控制系統(tǒng)效率以及硬件成本達(dá)到相對完美的平衡。

(3)遺傳算法是一種全局搜索方法，運(yùn)用遺傳算法成功的優(yōu)化了作動器在空間結(jié)構(gòu)中的位置，效果明顯，收斂快，大大節(jié)省了計算時間。位置優(yōu)化結(jié)果表明作動器多集中在柱面網(wǎng)殼的兩端，多布置在斜腹桿上。

(4)基于LQR主動控制算法，將GMM作動器集成到原桿件中進(jìn)行了結(jié)構(gòu)振動主動控制的數(shù)值模擬，并對控制效果進(jìn)行了分析，結(jié)果表明GMM作動桿可有效地減小結(jié)構(gòu)的動力反應(yīng)，為結(jié)構(gòu)的主動振動控制提供了一種新方法。

［1］Cho S W，Kim B W，Jung H J.Implementation of modal control for seismically excited structures using magnetorheological dampers［J］.Journal of Engineering Mechanics，2005，131(2):177-184.

［2］陳定方，舒 亮，盧全國.超磁致伸縮致動器建模與控制仿真［J］.機(jī)械工程學(xué)報，2007，43(8):12-16.

［3］徐 斌，歐進(jìn)萍，姜節(jié)勝.壓電智能桁架作動器配置的優(yōu)化設(shè)計［J］.機(jī)械科學(xué)與技術(shù)，2007，26(11):1511-1516.

［4］潘 繼，陳龍祥，蔡國平.柔性板壓電作動器的優(yōu)化位置與主動控制實(shí)驗(yàn)研究［J］.振動與沖擊，2010，29(2):117-120.

［5］雷英杰，張善文.MATLAB遺傳算法工具箱及應(yīng)用［M］.西安電子科技出版社，2005:2-3.

［6］歐進(jìn)萍.結(jié)構(gòu)振動控制——主動、半主動和智能控制［M］.北京:科學(xué)出版社，2003:39-44.

［7］孫東昌，張洪華，等.壓電智能梁振動控制中致動片優(yōu)化布置與實(shí)驗(yàn)［J］.中國空間科學(xué)技術(shù)，1992，12(6):46-52.

［8］曾 光，李東旭.空間智能桁架作動器/傳感器位置優(yōu)化中的遺傳算法應(yīng)用［J］.宇航學(xué)報，2007，28(2):461-464.

［9］張東昱.雙層柱面網(wǎng)殼主動控制及其優(yōu)化研究［D］.北京:北京工業(yè)大學(xué)，2001:27.

［10］容 里.壓電作動器設(shè)計及雙層球面網(wǎng)殼結(jié)構(gòu)振動主動控制研究［D］.杭州:浙江大學(xué)，2008:57.