財富管理的視角下基于NSGA-II的橋墩截面優(yōu)化設(shè)計

張士博 張寧 張世紅 張宇哲 張冠堯 練晨燁 宮中飛 趙雪鍵

[摘? 要：本文基于財富管理的視角，以客戶為中心，應(yīng)用NSGA-II建立以橋墩截面半徑和縱筋配筋率為決策變量的橋墩截面優(yōu)化設(shè)計模型，將橋墩截面優(yōu)化設(shè)計優(yōu)化問題模擬為生物進化問題，為客戶提供一套橋墩造價最小的設(shè)計方案。計算結(jié)果表明：NSGA-II方法得到的非劣解在目標(biāo)空間中分布均勻、計算速度快、收斂性和魯棒性好，該設(shè)計方法在滿足客戶對橋梁結(jié)構(gòu)性能的安全需求同時，降低了客戶投資風(fēng)險、實現(xiàn)效益最大化。

關(guān)鍵詞：橋墩;截面;NSGA-II;優(yōu)化設(shè)計

Abstract：In this paper，based on the perspective of wealth management，taking the customer as the center，the NSGA-II is applied to establish the optimal design model of pier section with the pier section radius and longitudinal reinforcement ratio as the decision variables.The optimal design problem of pier section is simulated as the biological evolution problem，and a set of design scheme with the lowest pier cost is provided for the customer.The calculation results show that the non inferior solution obtained by NSGA-II method is evenly distributed in the target space，with fast calculation speed，good convergence and robustness. The design method can meet the safety requirements of customers for the bridge structure performance，reduce the investment risk of customers and maximize the benefits.

Key words：bridge pier; section;NSGA-II;optimal design]

1引言

橋梁設(shè)計是一個復(fù)雜的過程，首先需要假設(shè)初始的截面尺寸和縱筋配筋率，然后對假設(shè)的截面配置進行分析驗算。在設(shè)計過程中需要執(zhí)行很多次迭代計算，直到假定的截面配置滿足設(shè)計規(guī)范，另外，還需要從截面配置方案中選擇滿足安全性、可靠性、設(shè)計規(guī)范及規(guī)范要求的造價、重量、美學(xué)等的最優(yōu)截面。設(shè)計師在滿足所有前面的標(biāo)準(zhǔn)的可選擇的截面配置方案中選擇最優(yōu)的截面配置是非常困難的。

本文基于財富管理的視角，以客戶為中心，首次提出將精英保留非劣排序遺傳算法（Elitist Non-dominated Sorting Genetic Algorithm，簡稱NSGA-II）應(yīng)用到橋墩截面優(yōu)化設(shè)計中來，以公路鋼筋混凝土及預(yù)應(yīng)力混凝土橋涵設(shè)計規(guī)范中圓形截面偏心受壓柱式墩計算模型為基礎(chǔ)，建立了以橋墩截面半徑和縱筋配筋率為決策變量，以抗彎強度需求能力比和橋墩的初始造價為目標(biāo)的多目標(biāo)優(yōu)化模型，并研究了該模型的具體求解方法。該方法既能夠給設(shè)計師提供更靈活的選擇空間，也可以降低了客戶投資風(fēng)險、實現(xiàn)效益最大化。

2建立圓柱橋墩截面優(yōu)化設(shè)計的數(shù)學(xué)模型

2.1設(shè)計變量

取圓柱橋墩截面的縱向鋼筋配筋率[ρ]和圓形截面半徑[r]作為設(shè)計變量。

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

目標(biāo)函數(shù)為橋墩的造價，以材料費用表示為

2.3約束函數(shù)

圓形截面鋼筋混凝土偏心受壓構(gòu)件正截面抗壓承載力約束

1）對于長細比[l0/i>17.5]的構(gòu)件，應(yīng)考慮構(gòu)件在彎矩作用平面內(nèi)的撓曲對軸向力偏心距的影響。此時，應(yīng)將軸向力對截面重心軸的偏心距[e0]乘以偏心距增大系數(shù)[η]。

2）計算系數(shù)A、B、C、D應(yīng)按下列公式計算

2.4圓形截面偏心受壓構(gòu)件承載力計算步驟和流程圖

（1）判斷墩柱長細比[l0/i]是否大于17.5，以此確定是否考慮偏心距增大系數(shù)。

（2）假定受壓區(qū)高度與圓形截面直徑的比值[ξ]。

（3）根據(jù)假定的[ξ]值帶入公式（7～10）得到計算參數(shù)A、B、C和D。

（4）將計算參數(shù)帶入公式2的右邊項，求得計算截面的設(shè)計承載力。

（5）根據(jù)求得的截面設(shè)計承載力反算偏心距[e00]，將其與設(shè)計要求的[e0]比較，若兩者相差不超過2%，則求得的截面設(shè)計承載力即為真實的承載力;否則重新假定[ξ]值，重復(fù)以上3～4，直到滿足精度要求為止。編程流程圖如下：

3精英保留非劣排序遺傳算法

NSGA-II的主要方法為：

（1）虛擬適應(yīng)度的計算

在NSGA-II算法中，個體的適應(yīng)度包括非劣解的等級和個體的虛擬適應(yīng)度.為了保持個體的多樣性、防止個體在局部堆積，NSGA-II算法首次提出了虛擬適應(yīng)度（dummy fitness）的概念，它指目標(biāo)空間上的每一點與同等級相鄰2點之間的局部擁擠距離。例如，圖1中目標(biāo)空間第[i]點的擁擠距離等于它在同一等級相鄰的點[i-1]和[i+1]在[f1]軸和[f2]軸距離的和，即由點[i-1]和[i+1]組成的矩形2個邊長之和。使用這一方法可自動調(diào)整小生境（niche），使計算結(jié)果在目標(biāo)空間比較均勻地散布，具有較好的魯棒性。

具體實現(xiàn)時，首先解碼染色體，然后按計算模型計算每個個體相應(yīng)的目標(biāo)函數(shù)，再根據(jù)目標(biāo)函數(shù)值進行非劣分層，計算每層個體的虛擬適應(yīng)度.虛擬適應(yīng)度的計算步驟如下：

5）對不同的目標(biāo)函數(shù)，重復(fù)步驟2）～4）操作。

（2）選擇運算

（3）精英策略

精英策略即保留父代中的優(yōu)良個體直接進入子代，它是遺傳算法以概率1收斂的必要條件。采用的方法是：①將父代[Pt]和子代[Qt]全部個體合成為一個統(tǒng)一的種群[Rt=Pt∪Qt]并放入進化池中，種群[Rt]的個體數(shù)成為[2N];②將種群[Rt]按非劣解等級分類并計算每一個體局部擁擠距離，依據(jù)等級的高低逐一選取個體，直至個體總數(shù)達到[N];③以此形成新一輪進化的父代種群[Pt+1]，其個體數(shù)為[N]。在此基礎(chǔ)上開始新一輪的選擇、交叉和變異，形成新的子代種群[Qt+1]。

精英保留非劣排序遺傳算法流程如圖3：

4NSGA-II用于圓柱橋墩截面優(yōu)化設(shè)計的具體實現(xiàn)

4.1實例與NSGA-II參數(shù)選擇

算例為某市立交橋，該橋為連續(xù)曲梁橋，鋼筋混凝土獨柱橋，計算模型取圓柱形橋墩高度[H=20m]，上部結(jié)構(gòu)集中質(zhì)量[W=1450t]，偏心距[e0=0.32m]墩柱混凝土采用[C30]，軸心抗壓強度設(shè)計值[fcd]取13.8MPa;主筋采用HRB335級，抗壓強度設(shè)計值[f'sd]取280MPa，沿截面周邊均勻布置?；炷痢⒖v筋的單價分別取為700元/m3、3400元/t。鋼的密度78.5kN/m3。

在本算例中，以抗彎強度需求能力比和橋墩的初始造價分別作為多目標(biāo)優(yōu)化問題的目標(biāo)函數(shù)，建立多目標(biāo)優(yōu)化模型如下：

考慮到圓柱橋墩截面優(yōu)化設(shè)計的具體要求，目前情況下，以圓柱橋墩截面優(yōu)化設(shè)計必須滿足正截面抗壓承載力為約束要求，為此可以采用罰函數(shù)法，構(gòu)造圓柱橋墩截面優(yōu)化設(shè)計目標(biāo)函數(shù)，即：

以橋墩截面半徑[r]和縱筋配筋率[ρ]為優(yōu)化參數(shù)，則待優(yōu)化的決策變量有[x1，x2]2個，即[X=x1，x2T]，采用實數(shù)編碼形式。優(yōu)化參數(shù)的邊界約束條件如下

其他參數(shù)選擇：種群規(guī)[pop=200]，最大進化代數(shù)[gen=300]，隨機聯(lián)賽規(guī)模[N=2]，交叉概率[Pc=0.9]，變異概率[Pm=0.1]。

4.2結(jié)果與分析

（1）NSGA-II算法經(jīng)過300次遺傳迭代后的結(jié)果如圖4，從圖4中的非劣解中設(shè)計師能夠在廣泛的范圍內(nèi)更靈活地選擇截面設(shè)計方案。

（2）由圖4可知，抗彎強度需求能力比越小，橋墩的初始造價越大，即抗彎強度越大，橋墩的初始造價越大。

（3）本文所用的方法可以滿足在不同的抗彎強度需求和初始造價提供更加靈活的選擇空間。

5結(jié)論與展望

通過將NSGA-II算法對規(guī)則鋼筋混凝土橋梁圓形單柱橋墩進行截面優(yōu)化設(shè)計，可以得到以下結(jié)論：

（1）經(jīng)過NSGA-II算法優(yōu)化之后可以獲得大量在設(shè)計空間中廣泛分布的Pareto最優(yōu)設(shè)計變量和最優(yōu)目標(biāo)函數(shù)序列，從而使得設(shè)計師能夠在廣泛的范圍內(nèi)更靈活地進行截面設(shè)計，而不僅僅局限于傳統(tǒng)意義上造價最小的設(shè)計方案。

（2）本文所用的算法只是對圓柱形橋墩截面的優(yōu)化，對于具體的橋墩抗震截面設(shè)計優(yōu)化還未涉及，還有待于作進一步研究。

（3）本文基于財富管理的視角下，以客戶為中心，為客戶提供一系列橋墩造價最小的設(shè)計方案，降低了客戶投資風(fēng)險、實現(xiàn)效益最大化。

參考文獻

[1]中交公路規(guī)劃設(shè)計院.JYDD62-2004公路鋼筋混凝土及預(yù)應(yīng)力混凝土橋涵設(shè)計規(guī)范[S].北京：人民交通出版社，2004.

[2]彭波.圓柱形橋墩截面設(shè)計[J].中國水運，2010，10（02）：151-152.

[3]石志源.橋梁結(jié)構(gòu)電算[M].北京：人民交通出版社，2010.

[4]雷英杰，張善文，李繼武，等.MATLAB遺傳算法工具箱及應(yīng)用[M].西安：電子科技大學(xué)出版社.2006.

[5]張士博，郄志紅，吳鑫淼，等.基于精英保留非劣排序遺傳算法的玉米非充分灌溉制度優(yōu)化[C].Chinese Control and Decision Conference，（CCDC08）2008，Vol.6.

[6]柳春光，張士博，柳英洲.基于全壽命抗震性能的近海橋梁結(jié)構(gòu)多目標(biāo)優(yōu)化設(shè)計方法[J].大連理工大學(xué)學(xué)報，2015，55（01）：39-46.

[7]柳春光，張士博，柳英洲.基于NSGA-Ⅱ的規(guī)則橋梁圓形單柱RC橋墩全壽命造價優(yōu)化設(shè)計[J].水利與建筑工程學(xué)報，2015，13（02）：1-6.

[8]中華人民共和國交通運輸部.JTG/TB02-01-2008公路橋梁抗震設(shè)計細則[S].北京：人民交通出版社，2008.