程 怡(陽(yáng)光學(xué)院 土木工程系， 福建 福州 350108)

電視塔的風(fēng)致振動(dòng)控制研究

程 怡
(陽(yáng)光學(xué)院 土木工程系， 福建 福州 350108)

風(fēng)災(zāi)是電視塔面臨的主要自然災(zāi)害之一，主動(dòng)控制方法可以有效地減小電視塔的風(fēng)致振動(dòng)。以一座空間鋼管桁架電視塔為算例，首先基于ANSYS軟件分析了該電視塔的動(dòng)力特性，然后利用諧波合成法模擬了其隨機(jī)風(fēng)場(chǎng)，同時(shí)，還提出了一種形狀記憶合金的主動(dòng)控制方法，并編制了相應(yīng)的線(xiàn)性二次型經(jīng)典最優(yōu)控制算法(LQR)，最后對(duì)比了采取控制措施前后電視塔的風(fēng)致響應(yīng)，結(jié)果表明形狀記憶合金控制方法可以有效地減少電視塔的風(fēng)致振動(dòng)。

電視塔；動(dòng)力特性；脈動(dòng)風(fēng)場(chǎng)；風(fēng)致振動(dòng)；主動(dòng)控制方法

隨著各國(guó)城市建設(shè)的發(fā)展，電視塔由原來(lái)只能作為發(fā)射電視訊號(hào)的單功能電視塔逐步向集通訊、展覽、旅游、廣告等文化活動(dòng)和商業(yè)活動(dòng)為一體的大型綜合電視塔發(fā)展，所以電視塔通常也是一個(gè)城市的地標(biāo)建筑物。電視塔一般采用自立式的塔式結(jié)構(gòu)，是高聳結(jié)構(gòu)中的一大組成部分。按結(jié)構(gòu)材料劃分，電視塔可分為三種：鋼結(jié)構(gòu)電視塔，鋼筋混凝土電視塔和鋼-混凝土組合電視塔。高度在300 m以上的電視塔主要是鋼筋混凝土結(jié)構(gòu)，而高度在300 m以下的電視塔主要采用鋼結(jié)構(gòu)。

但是由于鋼結(jié)構(gòu)截面小，高度大，剛度低，其抗風(fēng)問(wèn)題也尤為突出。由于電視塔這種高聳結(jié)構(gòu)的高度大，結(jié)構(gòu)抗彎剛度相對(duì)較柔，水平風(fēng)荷載成為塔體的控制因素，在現(xiàn)有技術(shù)條件下鋼結(jié)構(gòu)塔高度受到限制[1]。近年來(lái)，隨著高聳結(jié)構(gòu)風(fēng)致災(zāi)害的事故頻繁發(fā)生，如何減少高聳結(jié)構(gòu)在風(fēng)荷載作用下的風(fēng)致振動(dòng)成為人們關(guān)注的焦點(diǎn)，而電視塔作為一種重要的現(xiàn)代高聳結(jié)構(gòu)，因此開(kāi)展電視塔的抗風(fēng)性能研究及其風(fēng)致振動(dòng)的控制具有很重要的意義。

自20世紀(jì)50年代電視塔在世界各地開(kāi)始大量興建以后，國(guó)內(nèi)外學(xué)者就針對(duì)電視塔的抗風(fēng)性能及其振動(dòng)控制展開(kāi)了研究，主要集中在動(dòng)力特性、抗風(fēng)性能和控制措施三個(gè)方面：(1) 在動(dòng)力特性方面，張志強(qiáng)等[2]和劉學(xué)利等[3]均對(duì)高階振型對(duì)電視塔風(fēng)振的影響進(jìn)行了分析，結(jié)果表明高振型對(duì)鋼結(jié)構(gòu)電視塔的最大位移響應(yīng)很小，一般考慮三階或四階即可滿(mǎn)足計(jì)算精度要求，對(duì)于結(jié)構(gòu)的加速度需要考慮高振型的影響，應(yīng)該考慮八階振型才能滿(mǎn)足結(jié)構(gòu)舒適度計(jì)算的要求。但他們的研究不足之處是都將電視塔模型簡(jiǎn)化為葫蘆串式的二維模型，并未直接在原結(jié)構(gòu)進(jìn)行分析，所以得出的數(shù)據(jù)與真實(shí)數(shù)據(jù)存在一定的差別。張文福等[4]對(duì)大慶電視塔建立了有限元模型進(jìn)行模態(tài)分析，得出對(duì)于大慶電視塔高振型的影響不可忽略，鞭梢效應(yīng)顯著，天線(xiàn)的振動(dòng)在整個(gè)結(jié)構(gòu)中占主導(dǎo)地位，原因是鋼塔的質(zhì)量和剛度沿高度方向不均勻。因此，需要進(jìn)一步關(guān)注高階振型對(duì)電視塔的風(fēng)振影響。(2) 電視塔的抗風(fēng)性能研究。張志強(qiáng)[5]對(duì)合肥電視塔風(fēng)振及地震響應(yīng)的振動(dòng)控制研究分析表明合肥電視塔的地震響應(yīng)與脈動(dòng)風(fēng)荷載作用下的風(fēng)振響應(yīng)相比，電視塔的響應(yīng)以風(fēng)振響應(yīng)為主。唐振興[6]對(duì)高聳通信塔抗風(fēng)性能進(jìn)行了分析，他選取沿高度不變的Davenport譜為風(fēng)譜，用諧波疊加法來(lái)生成風(fēng)速時(shí)程，得出了在只考慮順風(fēng)向風(fēng)荷載的情況下根據(jù)規(guī)范計(jì)算較根據(jù)時(shí)程分析偏于保守。顧明等[7]對(duì)廣州新電視塔進(jìn)行了測(cè)力和測(cè)振的風(fēng)洞試驗(yàn)研究，并采用隨機(jī)振動(dòng)CQC法對(duì)該塔進(jìn)行了風(fēng)致抖振響應(yīng)及等效靜力風(fēng)荷載計(jì)算。羅書(shū)泉[8]對(duì)比分析了頻域法和時(shí)域法對(duì)高聳電視塔結(jié)構(gòu)的風(fēng)振響應(yīng)并相互驗(yàn)證了計(jì)算結(jié)果的可靠性，還分析了黏彈性阻尼器后對(duì)結(jié)構(gòu)響應(yīng)的控制效果。因此，為了研究電視塔的風(fēng)致抖振響應(yīng)，需要建立其地址區(qū)的風(fēng)場(chǎng)，然后利用時(shí)程分析方法計(jì)算結(jié)構(gòu)響應(yīng)。(3) 結(jié)構(gòu)振動(dòng)控制的形式按照外加驅(qū)動(dòng)能源的出力種類(lèi)主要可分為被動(dòng)控制、主動(dòng)控制、混合控制和半主動(dòng)控制。自20世紀(jì)50年代電視塔在世界各地開(kāi)始大量興建以后，國(guó)內(nèi)外學(xué)者就針對(duì)電視塔結(jié)構(gòu)的振動(dòng)控制展開(kāi)了研究[9-10]。國(guó)內(nèi)電視塔振動(dòng)控制的研究起步較晚，但是發(fā)展很快。王肇民[11]對(duì)高聳結(jié)構(gòu)振動(dòng)控制的理論、控制裝置設(shè)計(jì)和參數(shù)優(yōu)化進(jìn)行了系統(tǒng)的研究，并針對(duì)電視塔天線(xiàn)的控制提出了懸掛彈簧阻尼器TSD。還有其他一些學(xué)者結(jié)合電視塔振動(dòng)控制的工程實(shí)例在控制裝置設(shè)計(jì)、控制效果實(shí)測(cè)及參數(shù)優(yōu)化等方面做了許多工作[12-15]。然而，針對(duì)電視塔的主動(dòng)控制的方法和應(yīng)用相對(duì)還不成熟。

針對(duì)以上三方面的問(wèn)題，本文以一座空間鋼管桁架電視塔為算例，首先通過(guò)模態(tài)分析得到了該結(jié)構(gòu)的自振頻率和振型等，然后建立其地址區(qū)的風(fēng)場(chǎng)，最后選取了一種形狀記憶合金的主動(dòng)控制方法，采用基于結(jié)構(gòu)振動(dòng)狀態(tài)空間模型的LQR控制算法對(duì)電視塔的風(fēng)振響應(yīng)進(jìn)行數(shù)值計(jì)算，對(duì)比得到結(jié)構(gòu)在有無(wú)控制措施下位移、速度、加速度和內(nèi)力響應(yīng)，以驗(yàn)證其減振效果。

1 電視塔的有限元分析

選取一座總高度為168 m的電視塔，重量300 t，主體為鋼管空間桁架結(jié)構(gòu)，如圖1所示。塔體分為天線(xiàn)段、塔樓、塔身3部分，標(biāo)高92 m以下為5邊形塔身，92 m～102 m為一碟形塔樓，116 m以上為4邊形天線(xiàn)段，天線(xiàn)段的邊長(zhǎng)分別為2.5 m×2.5 m(標(biāo)高116 m～135.5 m)、2.0 m×2.0 m(標(biāo)高135.5 m～157.5 m)、0.7 m×0.7 m(標(biāo)高157.5 m～168.5 m)。將該電視塔簡(jiǎn)化為15個(gè)質(zhì)點(diǎn)的模型，塔柱、橫桿和天線(xiàn)采用三維梁?jiǎn)卧猙eam 189模擬，橫桿和斜桿采用桿單元link 8來(lái)模擬，(電視塔結(jié)構(gòu)闡述如表1所述)。其中塔柱和橫桿都分別分成15層和10層來(lái)準(zhǔn)確考慮各個(gè)層的截面特性，塔底截面節(jié)點(diǎn)采用固結(jié)形式處理，電視塔整體有限元模型一共有319個(gè)節(jié)點(diǎn)，329個(gè)單元，有限元模型見(jiàn)圖2。選用Block Lanczos方法對(duì)該塔進(jìn)行了模態(tài)分析，運(yùn)用子空間迭代法求出結(jié)構(gòu)的自振頻率和模態(tài)振型，表2中列出了該塔的前6階自振頻率、周期和振動(dòng)形式。

表1 電視塔的結(jié)構(gòu)參數(shù)

由振型圖3和自振頻率可知，由于結(jié)構(gòu)的對(duì)稱(chēng)性，出現(xiàn)多個(gè)頻率相同的振型：該電視塔的第一階和第二階頻率相同，均為0.414 Hz，但是天線(xiàn)的振動(dòng)方向不同；第三階和第四階頻率相同均為0.935 Hz，但是整體結(jié)構(gòu)的振動(dòng)方向不同；第五階表現(xiàn)為整體扭轉(zhuǎn)，第六階表現(xiàn)為塔樓的振動(dòng)。這種動(dòng)力特性的特點(diǎn)表明前六階振型對(duì)電視塔的結(jié)構(gòu)響應(yīng)存在一定的影響。

圖1 電視塔實(shí)物圖

圖2 電視塔的有限元模型

2 隨機(jī)風(fēng)場(chǎng)的模擬

功率譜密度函數(shù)是平穩(wěn)隨機(jī)過(guò)程的主要數(shù)字特征，風(fēng)速脈動(dòng)分量的功率譜函數(shù)表示了紊流中各頻率的成分所作貢獻(xiàn)的大小。許多風(fēng)工程專(zhuān)家對(duì)風(fēng)功率譜進(jìn)行了研究，得到了不同形式的風(fēng)速譜表達(dá)式。其中用于結(jié)構(gòu)設(shè)計(jì)的水平脈動(dòng)風(fēng)速譜主要有Davenport譜、Kaimal譜和Harris譜等；豎向脈動(dòng)風(fēng)速譜Panofsky-McCormick譜、Lumley-Panofsky譜等[14]。

模擬電視塔的風(fēng)速時(shí)程時(shí)，由于電視塔高度上各點(diǎn)的風(fēng)速不同，互譜密度函數(shù)需要考慮相位角，風(fēng)譜密度矩陣(復(fù)數(shù)矩陣)有可能不正定。一般情況下脈動(dòng)風(fēng)的互譜密度函數(shù)為：

(1)

由于電視塔節(jié)點(diǎn)較多，需簡(jiǎn)化電視塔的模擬區(qū)域，本文沿高度變化將電視塔從下到上分為15段，各段中心點(diǎn)的實(shí)際高度見(jiàn)表3。用每段中心點(diǎn)的風(fēng)速時(shí)程代表該段內(nèi)各節(jié)點(diǎn)的風(fēng)速時(shí)程。截止圓頻率：12.56，頻率等分?jǐn)?shù)：1024，采樣時(shí)距：0.25 s，相干系數(shù)：1，采樣時(shí)間：600 s，粗糙高度：0.01 m。對(duì)于實(shí)際邊界層風(fēng)場(chǎng)中高度Zj處的風(fēng)速UZj：

(2)

其中，Hg和α分別為梯度風(fēng)高度和平均風(fēng)速剖面的冪函數(shù)指數(shù)，對(duì)于該電視塔所處的B類(lèi)地貌，HG=500 m和α=0.16。Ug為實(shí)際邊界層風(fēng)場(chǎng)中梯度風(fēng)速。根據(jù)電視塔所在地氣象部門(mén)提供的資料，當(dāng)?shù)?0年重現(xiàn)期的基本風(fēng)速為29.6 m/s。根據(jù)式(2)得到不同高度處的平均風(fēng)速，如表3所示。

表3 電視塔架模擬點(diǎn)高度及各高度處平均風(fēng)速

本文利用MATLAB軟件，采用諧波合成法編寫(xiě)了模擬電視塔架結(jié)構(gòu)多點(diǎn)的脈動(dòng)風(fēng)速程序，相關(guān)參數(shù)如下：(1) 假定該電視塔架所在地10 m高處的設(shè)計(jì)風(fēng)速為30 m/s；(2) 時(shí)間步長(zhǎng)0.25 s，采樣時(shí)間600 s；(3) 選取能反映高度變化的Davenport譜為目標(biāo)功率譜。經(jīng)過(guò)模擬計(jì)算分別得到了各段中心點(diǎn)處的脈動(dòng)風(fēng)速時(shí)程曲線(xiàn)，如圖4所示。

圖3 前六階的結(jié)構(gòu)振型圖

為了驗(yàn)證生成的風(fēng)速時(shí)程符合要求，圖5是風(fēng)速時(shí)程的驗(yàn)證結(jié)果，由生成的脈動(dòng)風(fēng)功率譜和目標(biāo)功率譜的符合程度以及生成的脈動(dòng)風(fēng)互相關(guān)函數(shù)和目標(biāo)函數(shù)的符合程度可以看出此方法生成的風(fēng)速時(shí)程滿(mǎn)足要求。

由生成的風(fēng)速時(shí)程可以根據(jù)下面的公式得出風(fēng)荷載時(shí)程(F)：

(3)

(4)

P=P0+p

(5)

F=PAus

(6)

式中：P0、p、P分別代表平均風(fēng)壓、脈動(dòng)風(fēng)壓、總風(fēng)壓；ρ代表空氣密度；us代表體形系數(shù)。利用上面的計(jì)算方法，可得風(fēng)荷載時(shí)程曲線(xiàn)，也即得到該電視塔的風(fēng)場(chǎng)，如圖6所示。

圖4 各層點(diǎn)的風(fēng)速時(shí)程曲線(xiàn)

圖5 結(jié)構(gòu)135.5 m處的風(fēng)速和風(fēng)荷載時(shí)程曲線(xiàn)以及兩種功率譜的對(duì)比

圖6 風(fēng)荷載時(shí)程曲線(xiàn)

風(fēng)荷載可以寫(xiě)成：

(7)

(8)

(9)

(10)

高聳結(jié)構(gòu)考慮風(fēng)與結(jié)構(gòu)耦合的運(yùn)動(dòng)方程是一個(gè)復(fù)雜的非線(xiàn)性動(dòng)力方程組，本文采用通用有限元程序ANSYS和數(shù)學(xué)軟件MATLAB，針對(duì)結(jié)構(gòu)的幾何非線(xiàn)性，選用Newton-Raphson迭代法和Wilson-theta法等直接積分法進(jìn)行數(shù)值求解。

3 形狀記憶合金的主動(dòng)控制方法和算法

選用形狀記憶合金(SMA)作為控制器的智能材料，控制器由彈簧與SMA絲相組合，如圖7所示。因?yàn)镾MA具有相容性好、變形量大、加熱后驅(qū)動(dòng)力強(qiáng)、響應(yīng)慢，兼具感應(yīng)和驅(qū)動(dòng)功能，比較適合制作成傳感器/驅(qū)動(dòng)器應(yīng)用到土木工程中的結(jié)構(gòu)振動(dòng)控制中[15-16]。

圖7 形狀記憶合金控制方法

對(duì)于安裝了形狀記憶合金控制措施的電視塔整體結(jié)構(gòu)，在風(fēng)荷載和控制力作用下，結(jié)構(gòu)的閉環(huán)控制運(yùn)動(dòng)方程可表示為：

(11)

式中：F(t)為外部荷載作用力向量；U(t)為結(jié)構(gòu)上的控制力向量；Ds為風(fēng)荷載位置矩陣； Bs為控制器的控制力作用位置矩陣。

為了尋求最優(yōu)控制，使系統(tǒng)由任意給定的初始狀態(tài)定義Z(t0)=Z0到終端狀態(tài)Zf時(shí)，可設(shè)系統(tǒng)的二次型性能指標(biāo)為：

(12)

式中：性能指標(biāo)中的第一項(xiàng)是末值項(xiàng)，又稱(chēng)為終端代價(jià)，實(shí)際上是根據(jù)對(duì)控制的要求，在指定的時(shí)間點(diǎn)tf到來(lái)時(shí)，系統(tǒng)的終端狀態(tài)Zf接近預(yù)期狀態(tài)的程度。第二項(xiàng)中ZTQZ表示對(duì)結(jié)構(gòu)振動(dòng)響應(yīng)的控制，UTRU表示對(duì)消耗控制能量的懲罰，Q和R是兩個(gè)決定控制力和結(jié)構(gòu)反應(yīng)大小的相應(yīng)維數(shù)的權(quán)矩陣。由式(12)表示的目標(biāo)函數(shù)可以看出，權(quán)矩陣Q越大，結(jié)構(gòu)的反應(yīng)越小，控制力越大；R越小，控制力越大，結(jié)構(gòu)的反應(yīng)越小。

系統(tǒng)狀態(tài)最優(yōu)控制問(wèn)題即為在無(wú)限時(shí)間[t0，∞]內(nèi)，尋找最優(yōu)控制U(t)，使系統(tǒng)從初始狀態(tài)Z0轉(zhuǎn)移到零狀態(tài)附近，并使式(12)定義的性能泛函取極小值。因此，系統(tǒng)狀態(tài)最優(yōu)控制的問(wèn)題在數(shù)學(xué)上描述為：

(13)

(14)

式(14)描述的控制問(wèn)題是泛函條件極值問(wèn)題，可采用變分法來(lái)求解。構(gòu)造哈密爾頓函數(shù)H為：

(15)

式中：λT為拉格朗日乘子矢量。

使得哈密爾頓函數(shù)H取得極小值的必要條件為：

協(xié)態(tài)方程：

(16)

(17)

控制方程：

(18)

將式(15)代入式(17)，即可得到狀態(tài)方程。

將式(15)代入式(18)，便得到：

(19)

由于R為正定矩陣，所以：

U(t)=-R-1BTλ(t)

(20)

由上式便可以確定最優(yōu)控制，為了尋求最優(yōu)反饋還需要把U(t)與狀態(tài)向量Z(t)聯(lián)系起來(lái)。由式(15)和式(16)，得到：

(21)

式(20)確定的U(t)是λ(t)的線(xiàn)性函數(shù)。為了使U(t)能由狀態(tài)反饋實(shí)現(xiàn)，應(yīng)該建立λ(t)與Z(t)之間的線(xiàn)性變換關(guān)系。設(shè)：

λ(t)=P(t)Z(t)

(22)

將式(22)代入式(20)，得：

U(t)=-R-1BTP(t)Z(t)

(23)

下面介紹求解P(t)的方法：

將式(23)代入式(21)，得式(21)等號(hào)的左端為：

(24)

而式(21)等號(hào)的右端為：

-QZ(t)-ATλ(t)=-[Q+ATP(t)]Z(t)

(25)

將式(24)、式(25)代入式(21)，整理得：

(26)

按照控制的意義，當(dāng)t→∞時(shí)系統(tǒng)狀態(tài)充分接近零狀態(tài)，但還不是零狀態(tài)。這樣，可以得到：

P(t)t→∞=0

(27)

上式即為式(26)求解P(t)時(shí)應(yīng)滿(mǎn)足的邊界條件。

式(26)稱(chēng)為Riccati(里卡蒂)方程，P(t)稱(chēng)為Riccati矩陣函數(shù)，Riccati方程是一個(gè)非線(xiàn)性的矩陣微分方程。在求解式(25)Riccati方程時(shí)考慮到狀態(tài)反饋Z(t)的任意性，式(26)可以簡(jiǎn)化為：

(28)

由上式求出P(t)后，則可得到最優(yōu)控制力U(t)。

令：

G=R-1BTP(t)

(29)

則最后得到的最優(yōu)控制的表達(dá)式為：

U(t)=-R-1BTP(t)Z(t)=-GZ(t)

(30)

式中：G為最優(yōu)狀態(tài)反饋增益矩陣。P(t)為Riccati方程的解。

由式(30)可以看出，最優(yōu)控制U(t)取決于狀態(tài)反饋Z(t)和反饋增益G?；诖?，可以獲得結(jié)構(gòu)響應(yīng)的最優(yōu)的控制值。

4 控制效果對(duì)比分析

本文在該電視塔的每層質(zhì)點(diǎn)上安裝一個(gè)SMA控制器，但是根據(jù)表2中動(dòng)力特性分析，天線(xiàn)段先于整體結(jié)構(gòu)出現(xiàn)X、Z方向的振動(dòng)，故需要在塔頂天線(xiàn)段的X、Z方向分別安裝一個(gè)控制器。本文將ANSYS建立的模型調(diào)入MATLAB中，利用MABLAB建立了電視塔的風(fēng)荷載和控制器的控制矩陣，如式(11)中右邊荷載項(xiàng)，編寫(xiě)了SMA作動(dòng)器的LQR控制算法，最后在MATLAB進(jìn)行數(shù)值求解并輸出計(jì)算結(jié)果。

結(jié)構(gòu)模型是一個(gè)具有6n個(gè)自由度的電視塔結(jié)構(gòu)(n為結(jié)構(gòu)的節(jié)點(diǎn)數(shù))，在電視塔中安裝有16個(gè)SMA控制器。采用上述的基于結(jié)構(gòu)振動(dòng)狀態(tài)空間模型的LQR控制算法，對(duì)有無(wú)采用SMA控制器的主動(dòng)控制方法的電視塔風(fēng)振響應(yīng)進(jìn)行了數(shù)值模擬計(jì)算。為了分析主動(dòng)控制措施對(duì)高聳結(jié)構(gòu)電視塔的風(fēng)致響應(yīng)影響，本文主要對(duì)比分析了是否考慮SMA控制器作用下該塔重要部位的位移響應(yīng)、加速度響應(yīng)和內(nèi)力響應(yīng)，響應(yīng)最大值的結(jié)果見(jiàn)表4和表5。

表4 該塔重要位置處的風(fēng)致響應(yīng)最大值對(duì)比

表4和表5表明，SMA控制器可以有效地減少電視塔的風(fēng)致響應(yīng)：對(duì)于位移響應(yīng)， 塔底處的控制效果最大，達(dá)到了21.7%，塔頂?shù)奈灰埔矎?.945 m減少到0.813 m；57 m處速度響應(yīng)的控制效果最大，達(dá)到了17.65%，塔底處加速度響應(yīng)的效果為11.3%；對(duì)于塔底的剪力和彎矩響應(yīng)，主動(dòng)控制效果分別為5.0%和48.4%，特別是彎矩的最大值從8690 kN/m減小為4 486.7 kN/m。

表5 該塔塔底內(nèi)力響應(yīng)對(duì)比

5 結(jié) 論

本文首先建立了一座高聳電視塔的有限元模型和風(fēng)場(chǎng)，然后基于一種形狀記憶合金主動(dòng)控制方法編寫(xiě)了LQR控制算法程序，對(duì)比計(jì)算了該電視塔結(jié)構(gòu)在有無(wú)該主動(dòng)控制措施作用下結(jié)構(gòu)的風(fēng)振響應(yīng)。主要的結(jié)論如下：

(1) 利用ANSYS建立了電視塔的三維有限元模型，通過(guò)比較得出的前六階頻率和振型，并說(shuō)明高階振型不可忽略。

(2) 討論了脈動(dòng)風(fēng)速時(shí)程的模擬方法，針對(duì)電視塔的結(jié)構(gòu)幾何特性，并借助于MATLAB軟件，應(yīng)用諧波合成法編寫(xiě)風(fēng)速程序模擬了各層質(zhì)點(diǎn)的脈動(dòng)風(fēng)速時(shí)程。根據(jù)模擬風(fēng)荷載的方法，為對(duì)電視塔的風(fēng)振響應(yīng)及主動(dòng)控制研究提供了風(fēng)荷載。

(3) 介紹了形狀記憶合金主動(dòng)控制方法和編制了相應(yīng)的線(xiàn)性二次型(LQR)經(jīng)典最優(yōu)控制算法，然后對(duì)有無(wú)主動(dòng)控制作用下該電視塔風(fēng)振響應(yīng)進(jìn)行數(shù)值計(jì)算。

(4) 對(duì)比分析了電視塔在有無(wú)控制下的結(jié)構(gòu)位移、速度、加速度和內(nèi)力響應(yīng)。結(jié)果顯示，形狀記憶合金控制方法可以有效地減小結(jié)構(gòu)的風(fēng)致響應(yīng)，對(duì)于高聳結(jié)構(gòu)的風(fēng)致振動(dòng)控制是一種較好的主動(dòng)控制方法。此外，對(duì)于該控制系統(tǒng)的穩(wěn)定性、可靠性尚需進(jìn)一步深入研究。

[1] 張相庭.結(jié)構(gòu)風(fēng)工程:理論規(guī)范實(shí)踐[M].北京：中國(guó)建筑工業(yè)出版社,2006.

[2] 張志強(qiáng),李?lèi)?ài)群,何建平,等.高振型對(duì)合肥電視塔風(fēng)振響應(yīng)的影響分析[J].特種結(jié)構(gòu),2004,21(3):42-49.

[3] 劉學(xué)利，金建新，王肇民.鋼結(jié)構(gòu)電視塔風(fēng)振反應(yīng)的高振型影響[J].四川建筑科學(xué)研究，1998(1):25-27.

[4] 張文福，孫曉剛.大慶電視塔有限元分析[J].油氣田地面工程，2006，25(12):32.

[5] 張志強(qiáng)．合肥電視塔風(fēng)振及地震響應(yīng)的振動(dòng)控制研究[D].南京:東南大學(xué)，2003.

[6] 唐振興.高聳通信塔抗風(fēng)性能分析[D].長(zhǎng)沙:湖南大學(xué)，2011.

[7] 周日亙毅，顧 明，朱樂(lè)東.廣州新電視塔模型測(cè)力風(fēng)洞試驗(yàn)及風(fēng)致響應(yīng)研究Ⅱ:風(fēng)致響應(yīng)分析[J].土木工程學(xué)報(bào)，2009，42(7)：14-19.

[8] 羅書(shū)泉.高聳電視塔的風(fēng)振響應(yīng)及控制分析[D].西安：西安建筑科技大學(xué)，2010.

[9] Harris R I. The nature of wind, the modern design of wind-sensitive structures[C]//Construction Industry Research and Information Association， London: 1971.

[10] Roorda J. Tendon control in tall structures[J]. Journal of Structural Division, 1975,101(3):505-521.

[11] 王肇民.鋼結(jié)構(gòu)電視塔振動(dòng)控制及其工程實(shí)踐[J].特種結(jié)構(gòu)，1998，15(3):1-5.

[12] 陸 飛,程文嚷,李?lèi)?ài)群.南京電視塔風(fēng)振主動(dòng)控制的實(shí)施方案研究[J].東南大學(xué)學(xué)報(bào)，2002，32(5):799-803.

[13] 何敏娟，馬人樂(lè).336m黑龍江廣播電視塔懸掛水箱振動(dòng)控制及其實(shí)測(cè)[J].建筑結(jié)構(gòu)學(xué)報(bào)，2001,22(1):39-41.

[14] 汪大洋，孫作玉，趙桂峰.高層結(jié)構(gòu)風(fēng)振失舒率的模糊概率評(píng)價(jià)[J].水利與建筑工程學(xué)報(bào)，2013,11(4):28-33.

[15] 馬 駿.工業(yè)廠房風(fēng)場(chǎng)數(shù)值模擬和抗風(fēng)優(yōu)化[J].水利與建筑工程學(xué)報(bào),2015,13(1):187-192.

Active control of Wind-induced Vibration of a Television Tower

CHENG Yi

(CollegeofCivilEngineering,YangoCollege,Fuzhou,Fujian350108,China)

Wind is one of natural disasters for television towers, and active control method could effectively mitigate the wind-induced vibration of the television tower. A television tower was taken as an example in this study, the dynamic characteristic of the television tower was firstly analyzed based on the ANSYS software. Subsequently, the stochastic wind field was simulated by the harmonic synthesis. Also, a shape memory alloy (SMA) control method was selected and the numerical simulation of a classical linear optimum LQR control algorithm was established. Finally, the control effects of wind-induced responses of the tower structure with and without the SMA control method were compared. The results shows that the control method can effectively reduce the wind-induced vibration of tower structures.

television tower; dynamic characteristic; fluctuating wind field; wind-induced vibration; active control method

10.3969/j.issn.1672-1144.2016.06.044

2016-07-29

福建省中青年教師教育科研項(xiàng)目(JAT160627)

程 怡(1986—)，女，福建福州人， 講師，主要從事建筑結(jié)構(gòu)風(fēng)工程方面的研究工作。E-mail: 84342523@qq.com

TU311.3

A

1672—1144(2016)06—0227—08