混合PSO在鋼管混凝土減震框架優(yōu)化設(shè)計中的應(yīng)用研究

任鳳鳴，梁健偉

(1.廣州大學(xué)土木工程學(xué)院，廣東廣州 510006;2.廣東工業(yè)大學(xué)土木與交通工程學(xué)院，廣東廣州 510006)

鋼管混凝土結(jié)構(gòu)由于其豎向承載力高，抗震性能好，經(jīng)濟效益好等優(yōu)點［1］，在高層建筑中的應(yīng)用越來越廣泛.為了進一步提高鋼管混凝土結(jié)構(gòu)的抗震性能，已有學(xué)者對帶耗能減震裝置的鋼管混凝土結(jié)構(gòu)的抗震性能進行了研究.REN等［2］對一榀鋼管混凝土純框架和兩榀帶屈曲約束支撐的鋼管混凝土框架進行擬靜力試驗，對比得出耗能減震構(gòu)件明顯地改善了鋼管混凝土框架的抗震性能.JIA等［3］分別對相同尺寸的空心鋼管混凝土純框架和帶屈曲約束支撐的空心鋼管混凝土框架進行了擬靜力試驗，試驗結(jié)果表明帶屈曲約束支撐的空心鋼管混凝土框架具有良好的延性和耗能能力，屈曲約束支撐有效提高了結(jié)構(gòu)的整體強度、剛度和耗能能力.

對耗能減震結(jié)構(gòu)進行優(yōu)化，得到耗能減震裝置的最優(yōu)參數(shù)設(shè)置，成為了減震結(jié)構(gòu)設(shè)計亟需解決的問題之一.目前很多學(xué)者采用智能優(yōu)化算法對耗能減震結(jié)構(gòu)的優(yōu)化問題進行了研究.熊仲明等［4］采用循環(huán)布置法和遺傳算法對帶摩擦阻尼器的結(jié)構(gòu)進行優(yōu)化，通過對無阻尼器框架和兩個優(yōu)化過的減震結(jié)構(gòu)進行時程分析，結(jié)果表明遺傳算法優(yōu)化后的結(jié)構(gòu)具有更好的減震效果.曲激婷等［5］采用遺傳算法對帶黏彈性阻尼器的結(jié)構(gòu)進行阻尼器的位置優(yōu)化研究.FADI等［6］采用遺傳算法對設(shè)置了防屈曲支撐的耗能減震結(jié)構(gòu)分別進行了單目標(biāo)和多目標(biāo)優(yōu)化.

粒子群算法為智能優(yōu)化算法的一種，由于其參數(shù)較少，優(yōu)化效果好，應(yīng)用范圍越來越廣泛，但其在耗能減震結(jié)構(gòu)優(yōu)化方面的應(yīng)用甚少，本文把粒子群算法和遺傳算法進行結(jié)合，形成混合粒子群算法，并采用混合粒子群算法對帶防屈曲支撐的鋼管混凝土減震框架進行優(yōu)化分析.

1 算法原理

1.1 粒子群算法

粒子群算法由 KENNEDY等［7］于1995年共同提出，其理論主要基于大自然的鳥群和魚群的覓食行為.在鳥群尋找食物過程中，每個個體通過自身經(jīng)驗和群體經(jīng)驗更新速度和位置.

對于結(jié)構(gòu)優(yōu)化問題，設(shè)計變量的值對應(yīng)于鳥的位置，目標(biāo)函數(shù)的最優(yōu)解相當(dāng)于食物的位置.對于一個N維優(yōu)化問題，第i個設(shè)計變量表示為Xi=(xi1，xi2，…，xiN)，第i個粒子的速度表示為Vi=(vi2，vi2，…，viN)，個體的歷史最好位置表示為pbest，群體中所有個體的最優(yōu)位置表示為gbest.在每次迭代過程中，每個粒子根據(jù)速度更新公式(1)和位置更新公式(2)調(diào)整飛行的速度和個體位置.

式(1)中，w為慣性權(quán)重，其能調(diào)整上一次迭代時的速度對下一次迭代速度的影響;c1、c2為常數(shù)，稱為學(xué)習(xí)因子;rand()和Rand()是［0，1］內(nèi)的兩個不同的隨機數(shù).

1.2 遺傳算法

遺傳算法是由 HOLLAND［8］于1975年提出，其基本的理論是基于進化論中的“適者生存”和遺傳學(xué)的機理.

對于N維優(yōu)化問題，第i個設(shè)計變量Xi=(x1，x2，…xn)表示為一個染色體，xn為染色體中的基因.在優(yōu)化過程中，通過對染色體進行選擇、交叉和變異3種操作，不斷迭代，得到最優(yōu)解.其中選擇是指選擇適應(yīng)度值好的個體復(fù)制到下一代;交叉是把兩個配對的染色體的部分基因進行交換的過程;變異是指根據(jù)變異概率隨機對染色體中的某個基因進行變異，得到一個新個體.

1.3 混合粒子群算法

遺傳算法具有全局搜索能力強、魯棒性好［9］等優(yōu)點，粒子群算法的局部搜索能力可以通過調(diào)整慣性權(quán)重的值得到增強.把遺傳算法和粒子群算法結(jié)合，組成混合粒子群算法，在優(yōu)化前期采用遺傳算法，優(yōu)化后期采用粒子群算法，使混合粒子群算法具有前期全局搜索能力強，后期局部搜索能力好的特點.采用混合粒子群算法進行優(yōu)化，在優(yōu)化前期由于算法的全局搜索能力強，防止算法陷入局部最優(yōu)，在優(yōu)化后期，算法較好的局部搜索能力有利于得到最優(yōu)解.

混合粒子群算法的優(yōu)化流程如下:

(1)確定算法參數(shù)，包括種群數(shù)目、交叉概率、變異概率和代溝等;

(2)初始化種群;

(3)采用遺傳算法進行優(yōu)化，計算每個個體的適應(yīng)度值，對每個染色體執(zhí)行選擇、交叉和變異等遺傳操作，得到新的種群;

(4)判斷是否達(dá)到轉(zhuǎn)換算法的條件，滿足則跳到步驟(5)，否則返回步驟(3);

(5)把遺傳算法的個體傳遞給粒子群算法，采用粒子群算法進行優(yōu)化;

(6)計算每個個體的適應(yīng)度值，根據(jù)公式(1)和(2)更新粒子的速度和位置;

(7)判斷是否滿足結(jié)束條件，滿足則輸出最優(yōu)解，否則返回步驟(6).

2 數(shù)學(xué)模型

本文采用混合粒子群算法對帶防屈曲支撐的鋼管混凝土減震框架進行優(yōu)化，其中設(shè)計變量為各層設(shè)置的防屈曲支撐的核心單元的橫截面積，約束條件為結(jié)構(gòu)的層間位移角滿足規(guī)范規(guī)定的層間位移角限值和構(gòu)件內(nèi)力小于極限承載力，見公式(3)～(5).目標(biāo)函數(shù)為結(jié)構(gòu)所有防屈曲支撐的核心單元的橫截面積的總和，通過優(yōu)化得到最小值，見公式(6).

式中，Xi為第i層的防屈曲耗能支撐的核心單元的橫截面積，Amax為防屈曲耗能支撐的核心單元橫截面積的上限，θi為第i層的層間位移角，θlim為層間位移角限值，F(xiàn)為結(jié)構(gòu)構(gòu)件內(nèi)力，F(xiàn)lim為結(jié)構(gòu)構(gòu)件極限承載力，f(X)為所有防屈曲耗能支撐的核心單元橫截面積之和，n為結(jié)構(gòu)的總層數(shù).

3 算例

3.1 結(jié)構(gòu)概況

某10層鋼管混凝土框架，平面尺寸24 m×24 m，柱距8 m，層高4 m.結(jié)構(gòu)抗震設(shè)防烈度為8度，設(shè)計基本地震加速度為0.30 g，結(jié)構(gòu)處于Ⅱ類場地，設(shè)計地震分組為第二組.結(jié)構(gòu)平面圖見圖1，結(jié)構(gòu)構(gòu)件材料和尺寸見表1.

表1 結(jié)構(gòu)參數(shù)Table 1 Parameters of the structure

圖1 鋼管混凝土框架示意圖Fig.1 Figure of the CFST frame

對鋼管混凝土框架進行時程分析，向結(jié)構(gòu)輸入峰值加速度為110 gal的EL-centro波，得到結(jié)構(gòu)各層的最大層間位移角，見表3(表3中第二列數(shù)據(jù)).從表3可知，優(yōu)化前鋼管混凝土框架的第二、三層的層間位移角超過《鋼管混凝土結(jié)構(gòu)技術(shù)規(guī)程》(CECS28:2012)對鋼管混凝土框架最大層間位移角1/300的要求［10］，第四到第六層的層間位移角較大，為了提高結(jié)構(gòu)的抗震性能，減少結(jié)構(gòu)的層間位移角，在圖1(a)中軸線1和軸線4兩榀框架的中間跨對稱布置防屈曲支撐，并采用混合粒子群算法對單向地震波作用下的鋼管混凝土減震框架進行優(yōu)化.

3.2 算法參數(shù)

混合粒子群算法的參數(shù)分為兩部分，其中遺傳算法部分的種群數(shù)為40，交叉概率為0.7，變異概率為0.01，代溝為0.95;粒子群算法部分的種群目為10，慣性權(quán)重由0.9遞減至0.4，學(xué)習(xí)因子c1=c2=0.5.總的迭代次數(shù)為250次，其中前40次迭代采用遺傳算法，迭代40次以后采用粒子群算法進行優(yōu)化.

3.3 優(yōu)化結(jié)果

采用混合粒子群算法對帶防屈曲支撐的鋼管混凝土減震框架進行優(yōu)化的收斂過程見圖2.優(yōu)化后結(jié)構(gòu)各層的防屈曲支撐的核心單元的總橫截面積見表2.

圖2 鋼管混凝土減震框架優(yōu)化收斂曲線Fig.2 Convergence curves of the optimization process of the CFST frame with buckling-restrained braces

從圖2可見，混合粒子群算法的前期收斂速度較快，在優(yōu)化的后期，混合粒子群算法主要進行局部搜索，迭代到220次時，算法得到最終解.從表3可知，優(yōu)化后的結(jié)構(gòu)的底部樓層的防屈曲支撐核心單元橫截面積較大，頂部樓層的橫截面積較小.

為了對比優(yōu)化效果，基于總的核心單元橫截面積相等的前提下，分別對優(yōu)化前后的帶防屈曲支撐的鋼管混凝土減震框架進行時程分析.其中優(yōu)化前的鋼管混凝土減震框架的各層的防屈曲支撐的核心單元的橫截面積為0.003 36 m2.

優(yōu)化前后的帶防屈曲支撐的鋼管混凝土減震框架的各層最大層間位移角見表3(表3中第3和第4列數(shù)據(jù)).

表2 優(yōu)化結(jié)果Table 2 Optimization results

表3 結(jié)構(gòu)層間位移角Table 3 Inter-story drift angle of the structure

從表3可見，優(yōu)化后的鋼管混凝土減震框架底部樓層的層間位移角比優(yōu)化前小，頂部樓層的層間位移角比優(yōu)化前大.這是由于優(yōu)化后的結(jié)構(gòu)底部樓層的防屈曲支撐的核心單元的橫截面積較大，頂部樓層的橫截面積較小.從結(jié)構(gòu)的整體來看，無控框架的最大層間位移角在底部樓層.因此，進行優(yōu)化計算后，減少了結(jié)構(gòu)的最大層間位移角.

優(yōu)化前后鋼管混凝土減震框架的柱彎矩值見表4.表中數(shù)值為圖1中軸線1和軸線B交點處的柱彎矩值.

從表4可知，由于優(yōu)化后的鋼管混凝土減震框架底部樓層的防屈曲支撐的核心單元的橫截面積較大，頂部樓層的橫截面積較小，因此，優(yōu)化后的帶防屈曲支撐的鋼管混凝土減震框架底部樓層的柱彎矩值比優(yōu)化前小，頂部樓層的柱彎矩值比優(yōu)化前大.從結(jié)構(gòu)整體來看，無控框架底部樓層的柱彎矩值比頂部樓層大，因此，優(yōu)化后的結(jié)構(gòu)的防屈曲支撐的核心單元橫截面積分配更合理.

表4 柱彎矩值Table 4 Moment of the columns

優(yōu)化前后帶防屈曲支撐的鋼管混凝土減震框架的層間剪力對比見圖3.

圖3 結(jié)構(gòu)層間剪力對比Fig.3 Comparison of story shear of the structure

從圖3可見，優(yōu)化后的鋼管混凝土減震框架的各層層間剪力均有不同程度的降低.

總體而言，優(yōu)化后的帶防屈曲支撐的鋼管混凝土減震框架的最大層間位移角、最大柱彎矩值和層間剪力均得到降低，分布更為合理，優(yōu)化后結(jié)構(gòu)的抗震性能進一步得到提高.

4 結(jié)論

本文采用混合粒子群算法對帶防屈曲支撐的鋼管混凝土減震框架優(yōu)化進行了研究，通過對算例進行優(yōu)化，得出以下結(jié)論:

(1)混合粒子群算法對地震波作用下的鋼管混凝土減震框架進行參數(shù)優(yōu)化是可行的.

(2)把遺傳算法和粒子群算法進行結(jié)合，該新型的混合粒子群算法具有前期全局搜索能力強，后期局部搜索能力好的特點.

(3)采用混合粒子群算法對地震波作用下的鋼管混凝土減震框架進行優(yōu)化，能得到結(jié)構(gòu)各層最優(yōu)的防屈曲支撐的核心單元的橫截面積，改善了結(jié)構(gòu)的抗震性能，為鋼管混凝土減震框架優(yōu)化設(shè)計提供一個新方法.

［1］ 鐘善桐.鋼管混凝土統(tǒng)一理論——研究與應(yīng)用［M］.北京:清華大學(xué)出版社，2006.ZHONG S T.Unified theory of CFST［M］.Beijing:Tsinghua University Press，2006.

［2］ REN F M，ZHOU Y，ZHANG J B，et al.Experimental study on seismic performance of CFST frame structures with energy dissipation devices［J］.J Constr Steel Res，2013，90:120-132.

［3］ JIA M M，LU D G，GUO L H，et al.Experimental research and cyclic behavior of buckling-restrained braced composite frame［J］.J Constr Steel Res，2014，95:90-105.

［4］ 熊仲明，王佳，郭瑜.基于耗能減震結(jié)構(gòu)耗能裝置的優(yōu)化設(shè)置的研究［J］.西安建筑科技大學(xué)學(xué)報:自然科學(xué)版，2011，43(2):166-171.XIONG Z M，WANG J，GUO Y.Research on optimization setting based on passive energy dissipation structural device［J］.J Xi'an Univ Arch Tech:Nat Sci Edi，2011，43(2):166-171.

［5］ 曲激婷，李宏男.粘彈性阻尼器在結(jié)構(gòu)減震控制中的位置優(yōu)化研究［J］.振動與沖擊，2008，27(6):87-91.QU J T，LI H N.Optimal placement of viscoelastic dampers for passive response control［J］.J Vibrat Shock，2008，27(6):87-91.

［6］ FARHAT F，NAKAMURA S，TAKAHASHI K.Application of genetic algorithm to optimization of buckling restrained braces for seismic upgrading of existing structures［J］.Comput Struct，2009，87(1/2):110-119.

［7］ KENNEDY J，EBERHART R.Particle Swarm Optimization［C］∥Proceedings of the IEEE International Conference on Neural Networks，Piscataway，NJ，USA:1995:1942-1948.

［8］ HOLLAND J H.Adaptation in natural and artificial systems［M］.Ann Arbor:University of Michigan Press，1975.

［9］ 邊霞，米良.遺傳算法理論及其應(yīng)用研究進展［J］.計算機應(yīng)用研究，2010，27(7):2425-2429.BIAN X，MI L.Development on genetic algorithm theory and its applications［J］.Appl Res Comput，2010，27(7):2425-2429.

［10］CECS28:2012，鋼管混凝土結(jié)構(gòu)技術(shù)規(guī)程［S］.北京:中國計劃出版社.CECS28:2012，Technical specification for concrete-filled steel tubular structures［S］.Beijing:China Planning Press.