樓善治 張海靜

（上海大學(xué)土木工程系，上海200072）

0 引 言

隨著計(jì)算機(jī)技術(shù)的進(jìn)步，遺傳算法發(fā)展十分迅速，有大量關(guān)于遺傳算法的專著及研究文獻(xiàn)出現(xiàn)。在框架結(jié)構(gòu)設(shè)計(jì)優(yōu)化領(lǐng)域，也有越來(lái)越多的學(xué)者運(yùn)用遺傳算法解決工程問(wèn)題。

S.Rajeev［1］把離散變量的概念用于基于標(biāo)準(zhǔn)遺傳算法的鋼筋混凝土框架優(yōu)化設(shè)計(jì)研究。Hongjin Kim［3］利用實(shí)數(shù)編碼遺傳算法對(duì)鋼筋混凝土框架優(yōu)化設(shè)計(jì)，克服了二進(jìn)制編碼個(gè)體編碼冗長(zhǎng)的缺陷。陳雨生［4］采用平面框架有限元程序和MATLAB遺傳算法工具箱結(jié)合開(kāi)發(fā)了實(shí)用的優(yōu)化系統(tǒng)程序。但在優(yōu)化程序中僅考慮了構(gòu)件的強(qiáng)度約束，結(jié)構(gòu)的整體位移約束則通過(guò)手動(dòng)調(diào)整變量值來(lái)實(shí)現(xiàn)。張琦［5］以結(jié)構(gòu)造價(jià)最小為目標(biāo)函數(shù)，通過(guò)可靠度理論，對(duì)各類約束進(jìn)行處理，優(yōu)化結(jié)果與實(shí)際工程更貼近。但其在約束處理過(guò)程中僅考慮了構(gòu)件的強(qiáng)度和剛度約束，沒(méi)有考慮結(jié)構(gòu)的整體位移約束。

遺傳算法解決工程時(shí)的一個(gè)較大缺陷是在搜索過(guò)程中需要不斷對(duì)個(gè)體進(jìn)行結(jié)構(gòu)分析，從而導(dǎo)致計(jì)算效率低下。實(shí)際工程通常結(jié)構(gòu)體型較大，優(yōu)化時(shí)設(shè)計(jì)變量個(gè)數(shù)較多，遺傳算法的種群規(guī)模也需相應(yīng)地取較大值，若對(duì)每一個(gè)個(gè)體都進(jìn)行結(jié)構(gòu)分析，會(huì)造成優(yōu)化設(shè)計(jì)的時(shí)間成本過(guò)高的問(wèn)題。本文采用構(gòu)件內(nèi)力靜定化假設(shè)，在每一輪優(yōu)化的過(guò)程中凍結(jié)構(gòu)件內(nèi)力，避免了在優(yōu)化過(guò)程中進(jìn)行大量的結(jié)構(gòu)分析，極大地提高了計(jì)算效率。

學(xué)者們?cè)趯⑦z傳算法運(yùn)用于框架結(jié)構(gòu)優(yōu)化設(shè)計(jì)的領(lǐng)域上取得了一定的成果。但在上述列舉的文獻(xiàn)中，均未在優(yōu)化程序中同時(shí)考慮構(gòu)件約束和結(jié)構(gòu)整體位移約束。本文將優(yōu)化過(guò)程分為兩部分：第一部分采用遺傳算法對(duì)構(gòu)件尺寸和混凝土強(qiáng)度等級(jí)進(jìn)行優(yōu)化設(shè)計(jì)，通過(guò)懲罰函數(shù)處理構(gòu)件約束；第二部分采用準(zhǔn)則法對(duì)不滿足結(jié)構(gòu)位移約束的最佳個(gè)體進(jìn)行位移調(diào)整。這樣在確保構(gòu)件滿足約束的同時(shí)，又能使結(jié)構(gòu)整體滿足位移約束。

現(xiàn)有的工程問(wèn)題遺傳算法優(yōu)化設(shè)計(jì)研究中，很多學(xué)者采用的還是標(biāo)準(zhǔn)的交叉、變異、選擇算子，而傳統(tǒng)的遺傳算法存在著局部尋優(yōu)能力差、迭代過(guò)程緩慢、隨機(jī)性大、易出現(xiàn)振蕩和早熟收斂等缺點(diǎn)［6］。結(jié)構(gòu)設(shè)計(jì)優(yōu)化問(wèn)題是一種多變量、多約束的復(fù)雜問(wèn)題。為了得到更好的優(yōu)化效果，本文對(duì)遺傳算法進(jìn)行改進(jìn)，通過(guò)自適應(yīng)的操作算子，在保證優(yōu)化效果的同時(shí)，可以達(dá)到加速種群進(jìn)化收斂的目的，節(jié)省了時(shí)間成本。

1 框架結(jié)構(gòu)優(yōu)化模型

1.1 優(yōu)化變量設(shè)計(jì)

在材料單價(jià)確定的情況下，以造價(jià)為目標(biāo)函數(shù)的優(yōu)化問(wèn)題中，框架結(jié)構(gòu)的優(yōu)化變量應(yīng)為梁柱的b，h，縱向配筋值A(chǔ)s和箍筋值A(chǔ)sv。當(dāng)內(nèi)力被凍結(jié)，縱向配筋值和箍筋值都可以通過(guò)尺寸參數(shù)和內(nèi)力由規(guī)范［12-13］公式求得，也就是說(shuō)縱向配筋值A(chǔ)s和箍筋值A(chǔ)sv是關(guān)于構(gòu)件b和h的函數(shù)［7］。因此，最終的優(yōu)化設(shè)計(jì)變量?jī)H僅是梁柱的截面尺寸b和h。實(shí)際工程中為了便于施工，同層柱的混凝土等級(jí)一般相同，故本文將每層柱的混凝土等級(jí)作為變量。梁的混凝土等級(jí)采用C30不變。

1.2 目標(biāo)函數(shù)

本文以框架結(jié)構(gòu)的經(jīng)濟(jì)效益為研究目的，因此采用框架主體部分的造價(jià)為目標(biāo)函數(shù)。其中包括梁、柱混凝土，受力縱筋、箍筋以及模板。

1）單根梁造價(jià)

單根梁造價(jià)為

式中：b，h分別為截面寬度和高度；Cc為混凝土單位體積價(jià)格；Cs為鋼筋的單位重量?jī)r(jià)格；受力筋和箍筋的單價(jià)統(tǒng)一取用Cs；ρ為鋼筋的密度；Asu為上部架立筋截面面積之和；Asl為底部縱筋截面面積之和；As1，As2為梁端負(fù)彎起鋼筋截面面積；Cf為模板單位面積價(jià)格；l0為構(gòu)件長(zhǎng)度；n為整根桿件的箍筋根數(shù)；Asv為單肢箍筋截面面積；as為保護(hù)層厚度。

梁的配筋立面圖見(jiàn)圖1，配筋剖面圖見(jiàn)圖2。

圖1 梁配筋立面圖Fig.1 Elevation of beam steel bar

圖2 梁配筋剖面圖Fig.2 Profile of beam steel bar

2）單根柱造價(jià)

單根柱造價(jià)為

式中：As為縱筋面積之和；lsv為柱內(nèi)部箍筋長(zhǎng)度。

柱的立面圖見(jiàn)圖3，剖面圖見(jiàn)圖4。

圖3 柱配筋立面圖Fig.3 Elevation of column steel bar

圖4 柱配筋剖面圖Fig.4 Profile column steel bar

（3）鋼筋混凝土框架結(jié)構(gòu)的總造價(jià)為

式中，nc，nb分別為梁總數(shù)和柱總數(shù)。

1.3 約束條件

1.3.1 強(qiáng)度約束

1）梁

2）柱

（1）大偏心情況

（2）小偏壓情況（x>x）

式中：

除滿足式（7）以外，同時(shí)還需滿足式（6）的第二式。

另外，在各種情況下還需滿足斜截面承載力要求及軸壓比限制要求：

式中，αac為軸壓比限值。

框架結(jié)構(gòu)抗震等級(jí)為一級(jí)、二級(jí)、三級(jí)、四級(jí)時(shí)分別取0.65，0.75，0.85，0.90。

1.3.2 構(gòu)造約束

構(gòu)造約束參照《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》《建筑抗震設(shè)計(jì)規(guī)范》中的各項(xiàng)條文，如梁截面高寬比小于或等于4；構(gòu)件截面最大及最小配筋率、配箍率；構(gòu)件中鋼筋最大及最小間距、箍筋肢數(shù)等。

1.3.3 位移約束

每層最大彈性層間位移角

式中：Δi為第i層的層間相對(duì)位移；h i為第i層的層高。

2 自適應(yīng)遺傳算法及優(yōu)化流程

2.1 自適應(yīng)遺傳算法

遺傳算法是借鑒了進(jìn)化生物學(xué)中的一些現(xiàn)象而發(fā)展起來(lái)的一種啟發(fā)式算法，這些現(xiàn)象包括遺傳、突變、自然選擇以及雜交等。傳統(tǒng)遺傳算法是僅使用選擇算子、交叉算子和變異算子的最基本的遺傳算法，其控制參數(shù)是事先確定且在遺傳進(jìn)化過(guò)程中保持不變的。如何最優(yōu)給定有關(guān)控制參數(shù)，以便使遺傳算法的性能最優(yōu)，本身就是一個(gè)極其復(fù)雜的優(yōu)化問(wèn)題，要從理論上嚴(yán)格解決這一問(wèn)題還很困難。自適應(yīng)遺傳算法通過(guò)對(duì)遺傳參數(shù)的自適應(yīng)調(diào)整，大大提高了遺傳算法的收斂精度，加快了收斂速度［7］。

2.1.1 適應(yīng)度函數(shù)

適應(yīng)度是遺傳算法中評(píng)判個(gè)體優(yōu)劣的標(biāo)準(zhǔn)。適應(yīng)度函數(shù)值越大的個(gè)體越接近最優(yōu)解，能遺傳到下一代的概率也越大。反之則概率越小。選取合適的適應(yīng)度函數(shù)可以保證優(yōu)秀個(gè)體能夠有機(jī)會(huì)使優(yōu)良特性得以遺傳。本文采用的適應(yīng)度函數(shù)為

式中，G(x)是自適應(yīng)懲罰函數(shù)，在2.1.3節(jié)會(huì)詳細(xì)介紹。

2.1.2 選擇、交叉、變異算子

1）自適應(yīng)交叉算子

交叉概率和變異概率的取值對(duì)遺傳算法的性能有很大的影響。交叉概率Pc值過(guò)大，新個(gè)體將快速出現(xiàn)，種群多樣性更強(qiáng)，但容易破壞優(yōu)秀個(gè)體的基因，造成優(yōu)秀個(gè)體流失。Pc值過(guò)小，容易陷入局部最優(yōu)解。傳統(tǒng)遺傳算法中Pc的取值往往需要多次試算來(lái)決定。為了克服上述困難，本文采用動(dòng)態(tài)的自適應(yīng)交叉算子，根據(jù)待交叉的兩個(gè)個(gè)體的適應(yīng)值來(lái)計(jì)算交叉概率，提高搜索速度的同時(shí)也保護(hù)優(yōu)秀個(gè)體不受破壞［8］。自適應(yīng)交叉概率的計(jì)算公式為

式中：favg為當(dāng)前種群的平均適應(yīng)值；fmax為當(dāng)前種群中的最佳適應(yīng)值；f′為待交叉的兩個(gè)個(gè)體中適應(yīng)度更高者的適應(yīng)值。Pc1，Pc2為相關(guān)常數(shù)，可取Pc1=0.9，Pc2=0.6。

當(dāng)待交叉的兩個(gè)個(gè)體中適應(yīng)值較大者的適應(yīng)值不小于平均適應(yīng)值時(shí)，調(diào)整減小交叉概率Pc；當(dāng)待交叉的兩個(gè)個(gè)體的適應(yīng)值均小于平均適應(yīng)值時(shí)，交叉概率Pc等于原先設(shè)定的初始交叉概率Pc1。

2）自適應(yīng)變異算子

變異概率Pm過(guò)大會(huì)造成盲目搜索，Pm過(guò)小則難以產(chǎn)生新個(gè)體，種群多樣性差。本文同樣采用動(dòng)態(tài)的自適應(yīng)交叉算子，自適應(yīng)變異概率的計(jì)算公式為

式中：f為待變異個(gè)體的適應(yīng)值；Pm1，Pm2為相關(guān)常數(shù)，可取Pm1=0.1，Pm2=0.01。當(dāng)待變異個(gè)體的適應(yīng)值較大者的適應(yīng)值不小于平均適應(yīng)值時(shí)，按式（14）調(diào)整減小變異概率Pm；當(dāng)待變異個(gè)體的適應(yīng)值均小于平均適應(yīng)值時(shí)，Pm等于原先設(shè)定的初始變異概率Pm1。

3）選擇算子

本文根據(jù)選用適應(yīng)度比例法進(jìn)行選擇。適應(yīng)度比例法，又稱輪盤賭選擇法，是將個(gè)體的相對(duì)適應(yīng)度作為該個(gè)體被選擇的概率，某個(gè)體的相對(duì)適應(yīng)度越高，它被選擇的概率就越大。

2.1.3 約束處理

本文采用懲罰函數(shù)處理約束問(wèn)題。懲罰函數(shù)是將約束條件的違反度作為懲罰項(xiàng)加到目標(biāo)函數(shù)中，從而構(gòu)造出帶參數(shù)的增廣目標(biāo)函數(shù)。個(gè)體違反約束的程度越大，其對(duì)應(yīng)的懲罰項(xiàng)的值和增廣目標(biāo)函數(shù)值越大，個(gè)體適應(yīng)度值越小，從而達(dá)到使違反約束程度大的個(gè)體遺傳到下一代的概率減小的“懲罰”目的。這種約束處理的思想是把一系列的約束優(yōu)化問(wèn)題轉(zhuǎn)化成無(wú)約束的優(yōu)化問(wèn)題求解。

在約束條件下求函數(shù)最小值的數(shù)學(xué)表達(dá)式為

式中：f（x）為目標(biāo)函數(shù)；Pm1和Pm1分別為不等式約束條件和等式約束條件；x?Ω?S,S?Rn是搜索空間，Ω是可行域，即滿足所有約束條件的空間。

懲罰函數(shù)方法把問(wèn)題轉(zhuǎn)化為min G(x)=f(x)+cφ(x)，其中，G（x）稱為懲罰函數(shù)，c為懲罰因子，φ(x)為懲罰項(xiàng)。懲罰項(xiàng)個(gè)體到可行域的距離，個(gè)體x到第j個(gè)約束條件的距離可以表示為

傳統(tǒng)的懲罰函數(shù)中，懲罰系數(shù)c的選取比較困難，c過(guò)大會(huì)降低對(duì)不可行域的搜索，因而失去由不可行解提供的一些有價(jià)值的的信息；c過(guò)小會(huì)導(dǎo)致約束條件被忽略。為了解決上述問(wèn)題，蔡海鸞［10］提出了一種改進(jìn)的自適應(yīng)懲罰函數(shù)，利用當(dāng)前群體總的可行解的比例來(lái)權(quán)衡目標(biāo)函數(shù)和約束違反程度，實(shí)現(xiàn)對(duì)懲罰系數(shù)的動(dòng)態(tài)調(diào)整。假設(shè)懲罰系數(shù)函數(shù)為

式中：ρ為可行解在種群中所占的比例；α是一個(gè)調(diào)整的參數(shù)，α>0，可?。?，10］之間的一個(gè)數(shù)。

綜上可得自適應(yīng)懲罰函數(shù)為

本文用自適應(yīng)懲罰函數(shù)處理構(gòu)件約束。層間位移角約束采用下文介紹的位移調(diào)整策略來(lái)實(shí)現(xiàn)。

2.2 準(zhǔn)則法位移調(diào)整策略

在完成一次遺傳算法后需要對(duì)獲得的最優(yōu)個(gè)體進(jìn)行位移驗(yàn)算，如果不滿足位移規(guī)范的要求，則從該個(gè)體出發(fā)，在其附近搜索滿足位移約束的新個(gè)體，并令其替換為本輪整體優(yōu)化的最優(yōu)個(gè)體［11］。

本文結(jié)合傳統(tǒng)設(shè)計(jì)時(shí)調(diào)整不滿足層間位移角的思路，提出了一種基于準(zhǔn)則法的簡(jiǎn)化調(diào)整策略，使得在遺傳算法最優(yōu)個(gè)體附近快速找到一個(gè)滿足位移約束的個(gè)體。

傳統(tǒng)結(jié)構(gòu)設(shè)計(jì)實(shí)踐中，當(dāng)框架結(jié)構(gòu)不滿足位移約束時(shí)，有兩種常見(jiàn)的調(diào)整思路：

（1）放大本層不滿足位移要求方向的梁尺寸，一般只放大梁高尺寸；

（2）放大本層與不滿足位移要求方向的同向柱邊尺寸。

基于以上兩種思路，本文設(shè)計(jì)了如圖5所示的調(diào)整位移不滿足約束的搜索策略：首先，把遺傳算法計(jì)算出的最優(yōu)個(gè)體分別按照上述（1），（2）思路進(jìn)行模型調(diào)整；然后，比較兩種思路調(diào)整后的個(gè)體適應(yīng)值，取個(gè)體適應(yīng)值最大的替換遺傳算法的優(yōu)化結(jié)果。

圖5 位移不滿足約束時(shí)的調(diào)整策略Fig.5 Adjustment strategy when the displacement does not meet the constraint

2.3 整體優(yōu)化流程

本文采用分部?jī)?yōu)化思想來(lái)處理框架結(jié)構(gòu)的尺寸和內(nèi)力的非線性問(wèn)題。整體優(yōu)化分為兩步執(zhí)行，第一步凍結(jié)內(nèi)力做結(jié)構(gòu)尺寸的優(yōu)化，第二步在第一步得到的新最優(yōu)個(gè)體的基礎(chǔ)上，更新模型內(nèi)力，再次執(zhí)行第一步的操作，反復(fù)這個(gè)過(guò)程直到造價(jià)滿足收斂條件。

第一步優(yōu)化又分兩個(gè)級(jí)別。第一級(jí)為不考慮結(jié)構(gòu)剛度對(duì)內(nèi)力的影響，在梁柱構(gòu)件約束和層間約束下執(zhí)行遺傳算法；第二級(jí)為在遺傳算法優(yōu)化得到的最佳個(gè)體后，將變量回代入有限元分析軟件進(jìn)行受力分析并輸出層間位移，分析結(jié)果若不滿足位移約束，則執(zhí)行2.2節(jié)的位移調(diào)整策略不斷更新ETABS模型，直到滿足位移約束。整個(gè)優(yōu)化過(guò)程的執(zhí)行流程圖見(jiàn)圖6。

圖6 鋼筋混凝土框架結(jié)構(gòu)整體優(yōu)化流程Fig.6 Overall optimization process of reinforced concrete frame structure

3 框架結(jié)構(gòu)優(yōu)化實(shí)例

某8層鋼筋混凝土框架結(jié)構(gòu)，左右對(duì)稱，結(jié)構(gòu)層平面圖見(jiàn)圖7。初始模型各變量參數(shù)輸入值見(jiàn)表2。結(jié)構(gòu)設(shè)防烈度為7度（0.1g），地震分組為第一組。樓面恒荷載輸入為2.5 kN/m2，活荷載輸入為2.0 kN/m2，基本風(fēng)壓為0.55 kN/m2，地面粗糙程度為B類，場(chǎng)地類別為Ⅱ類。

圖7 標(biāo)準(zhǔn)層平面圖（單位：mm）Fig.7 Standard floor plan（Unit：mm）

圖8 為該框架模型優(yōu)化過(guò)程中的造價(jià)曲線。整體優(yōu)化在第5次達(dá)到收斂精度要求，總耗時(shí)為18分33秒，框架優(yōu)化后的造價(jià)收斂于428 833元，與原始模型相比減少了近25%的造價(jià)。表1為優(yōu)化前后結(jié)構(gòu)層間位移角，表2為優(yōu)化前后各變量的值的對(duì)比?？梢?jiàn)結(jié)構(gòu)層間位移角均滿足約束要求。

圖8 結(jié)構(gòu)整體優(yōu)化造價(jià)下降曲線Fig.8 Cost decline curve of overall structural optimization

表1 框架結(jié)構(gòu)優(yōu)化后的層間位移角Table 1 Inter story displacement angle of frame structure after optimization

表2 框架結(jié)構(gòu)優(yōu)化前后數(shù)據(jù)對(duì)比Table 2 Data comparison of frame structure around optimization

4 結(jié) 論

本文采用自適應(yīng)遺傳算法，在內(nèi)力靜定化假設(shè)的前提下對(duì)鋼筋混凝土框架結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)。內(nèi)力靜定化假設(shè)避免了在優(yōu)化過(guò)程中進(jìn)行大量的結(jié)構(gòu)分析，極大地提高了計(jì)算效率。自適應(yīng)遺傳算法通過(guò)對(duì)遺傳參數(shù)的自適應(yīng)調(diào)整，大幅度提高了遺傳算法的收斂精度，加快了收斂速度。通過(guò)算例實(shí)踐，驗(yàn)證了本文設(shè)計(jì)的框架優(yōu)化程序的可行性和實(shí)用性。在大大提高工作效率的同時(shí)，也能保證優(yōu)化效果良好。