基于改進蝙蝠算法的鋼筋混凝土箱梁橋截面優(yōu)化設(shè)計研究

基于改進蝙蝠算法的鋼筋混凝土箱梁橋截面優(yōu)化設(shè)計研究袁冰玲

橋梁截面尺寸參數(shù)轉(zhuǎn)變會影響構(gòu)造的整體重量及受力性能，在保證構(gòu)造受力合理的前提下優(yōu)化截面參數(shù)，可節(jié)省本錢。文章基于改良的蝙蝠算法，結(jié)合有限元模型，對某已建鋼筋混凝土箱梁橋截面進展優(yōu)化設(shè)計，比照分析優(yōu)化前后參數(shù)、撓度等數(shù)據(jù)變化狀況。討論結(jié)果說明：改良蝙蝠算法改良了原始算法多樣性缺乏、易于陷入局部最有解的缺陷;截面設(shè)計優(yōu)化后，箱梁總體構(gòu)造重量削減了9.2%，構(gòu)造剛度有些許降低，撓度結(jié)果影響較小;優(yōu)化后的橋梁整體構(gòu)造安全牢靠，使用階段承載力滿意要求。

鋼筋混凝土箱梁橋;截面優(yōu)化;改良蝙蝠算法;有限元模型;優(yōu)化設(shè)計

U442.5A351224

0引言

20世紀(jì)我國橋梁建立興起，受限于計算工具以及施工工藝，橋梁設(shè)計安全指標(biāo)較大，這就導(dǎo)致了橋梁構(gòu)造笨重以及材料的鋪張，目前有諸多學(xué)者采納各種方式對橋梁構(gòu)造進展優(yōu)化設(shè)計。劉明慧[1]采納掌握變量法針對鋼板-混凝土組合梁橋進展優(yōu)化設(shè)計，對于高跨比、寬厚比以及橫向聯(lián)系的設(shè)計值給出了合理化建議;魯業(yè)紅[2]運用多目標(biāo)模糊算法，結(jié)合T梁橋工程實例進展構(gòu)造優(yōu)化，節(jié)省了造價本錢，跨中截面彎矩降低，剛度提升13.3%;孫潔等[3]運用ANSYS有限元模型，進展多變量優(yōu)化分析，得出波浪鋼腹板尺寸以及箱梁截面尺寸數(shù)值;燕松波等[4]針對PC斜拉橋主梁箱形截面缺乏之處，通過優(yōu)化截面高度、頂?shù)装搴穸纫约皩挾鹊葏?shù)提高橋梁受力性能及局部穩(wěn)定性，并給出了相關(guān)截面參數(shù)值取值的合理化建議。

本文引入改良的蝙蝠優(yōu)化算法，結(jié)合有限元分析模型，依托某已建鋼筋混凝土箱梁橋，以混凝土總重為目標(biāo)函數(shù)，對其截面參數(shù)進展優(yōu)化分析，得出各變量優(yōu)化數(shù)值并進展優(yōu)化前后的結(jié)果比照，同時給出合理化建議，為同類橋梁設(shè)計供應(yīng)肯定參考。

1蝙蝠算法

自然界中蝙蝠在黑夜中能夠精確進展捕食或避開障礙物的阻擋，主要是通過自身主動發(fā)出和承受聲波來進展推斷，2023年Yang[5]依據(jù)蝙蝠這種捕食或避開障礙物的這種行為提出了蝙蝠算法。蝙蝠算法中各參數(shù)分別模擬的是蝙蝠聲音中的響度、放射率、頻率以及本身的位置和速度等。

算法中，t時刻下，蝙蝠個體i進展全局搜尋時的公式如下所示：

fi=fmin+fmax-fminβ

vti=vt-1i+xt-1i-x*（1）

xti=xt-1i+vti

式中：fmax、fmin——頻率的上、下限;

β——隨機變量，其值介于0～1之間;

x*——計算當(dāng)前狀態(tài)下個體處于的最優(yōu)位置。

蝙蝠算法在進展全局搜尋的過程中，為了保證整體的尋優(yōu)效果，會以肯定的概率進展局部搜尋，其更新公式為：

xnew=x*+ε·At（2）

式中：ε——勻稱分布于[-1，1]的隨機變量;

At——時刻t下，蝙蝠種群的平均響度。

在算法進展過程中，響度以及放射率都要隨著算法的運行不斷更新，其公式為：

Ait=αAt-1i

rti=r0i1-exp[-γ（t-1）]00（3）

從式（1）～（3）可以看出，隨著算法的運行，響度漸漸趨向于0，rti趨向于r0i。其中，rti主要作為局部搜尋的推斷條件。

2蝙蝠算法的改良

2.1種群多樣性改良

針對種群多樣性改良方面，此處引入反向?qū)W習(xí)[6]的概念，其原理是比照原始解和反向解的優(yōu)劣，擇優(yōu)者作為下一代個體，其初始化如下。

生成根本初始種群x0，種群數(shù)量為No，此時i個體的j維重量表達(dá)式為：

x0ij=xjmin+xjmax-xjminrj（4）

式中：xjmax——j維變量的上限;

xjmin——j維變量的下限。

解出反向解x～0ij：

x～0ij=xjmin+xjmax-x0ij（5）

比照分析x0i與x～0i優(yōu)劣，選取最優(yōu)解作為計算的初始種群。

2.2算法過程中參數(shù)優(yōu)化

在蝙蝠算法的過程中，其速度公式?jīng)]有肯定優(yōu)化措施，這會使得搜尋時間增加，同時無法保證尋優(yōu)質(zhì)量，不易收斂。為了避開上述狀況的發(fā)生，此處選擇引入權(quán)重系數(shù)ω與調(diào)整因子c[7]。更新后的蝙蝠飛行速度公式為：

vt+1i=ωvti+x*-xtifi+xi*-xtic

ω=ωmax-ωmax-ωtNmax-NtN+ωe

ωti1.2×ωti…iffxti≤fxt-1i

0.2×ωti…else（6）

式中：權(quán)重ω∈（0，1）;

ωi——權(quán)重初值;

ωmax——權(quán)重最大值;

ωe——權(quán)重變化值。

調(diào)整因子c表達(dá)式為：

c=cmin+cmax-cminω-ωmaxωmax-ωmin（7）

式中，調(diào)整因子c屬于區(qū)間（1，2）。調(diào)整因子可以讓距離食物比擬遠(yuǎn)的蝙蝠權(quán)重降低，使得這些遠(yuǎn)離目標(biāo)的蝙蝠不會使得種群有趨向局部最優(yōu)的傾向。其速率更新為：

rti=1-1-r0iγt-1（8）

2.3算法流程

（1）初始化種群。生成臨時種群，依據(jù)反向?qū)W習(xí)原理得出反向種群，評價兩個種群中的優(yōu)者形成算法初始種群。

（2）依據(jù)根本算法更新蝙蝠速度，依據(jù)目標(biāo)函數(shù)計算蝙蝠個體適應(yīng)度，得出當(dāng)前最優(yōu)蝙蝠。若滿意停頓運算條件則退出，不滿意則根據(jù)式（6）和式（7）進展權(quán)重與因子計算，轉(zhuǎn)變飛行速度。

（3）挑出蝙蝠種群中一個個體，依據(jù)式（3）和式（8）算出其速率和響度，若其速率大于當(dāng)前最優(yōu)蝙蝠，響度小于最優(yōu)個體，則次挑出個體為當(dāng)前最優(yōu)。

（4）除去第三步中的個體外再選擇一個進展比照，若響度小于上述的當(dāng)前最優(yōu)個體，則進展替代，更新最優(yōu)解。

（5）若當(dāng)前最優(yōu)解滿意要求，則終止，不滿意，則返回其次步進展循環(huán)計算，直至滿意要求。

3有限元模型的建立

本文選取湖南某已建大路路段鋼筋混凝土箱梁橋，橋梁跨徑（17+23+17）m，截面為單箱雙室截面，截面高1.3m，橋面寬13.5m，全橋采納C50混凝土材料，鋼筋標(biāo)號采納HRB335，橋梁設(shè)計荷載為汽車-超20。其截面如下頁圖1所示。

依據(jù)圖紙以及優(yōu)化變量可以建立相應(yīng)的MidasCivil有限元模型，變量轉(zhuǎn)變只需修改模型中的截面尺寸即可。全橋共離散成35個節(jié)點、34個單元，約束為一個支座約束全部平動方向，另外三個支座約束DY、DZ方向，二期荷載通過單元荷載施加，有限元模型如圖2所示。

4算法優(yōu)化及結(jié)果分析

4.1蝙蝠算法優(yōu)化模型構(gòu)建

本文擬將截面頂、底板高程，內(nèi)、外腹板厚度及高程作為優(yōu)化變量，分別表示為X1、X2、X3、X4、X5[8]，此處考慮實際橋梁的行車寬度，沒有將截面寬度作為優(yōu)化變量，因此進展截面承載力驗算時截面鋼筋排布間距根據(jù)原設(shè)計擬定，依據(jù)圖紙截面幾何特性，以橋梁構(gòu)造總重作為優(yōu)化目標(biāo)函數(shù)，此優(yōu)化目標(biāo)實行無量綱化，其表達(dá)式為：

G=1.5+2.6X4+2.6X3+8.5X1+8.5X2+[（1.3-X1-X2-X5）/2+0.15]+[0.1+（1.3-X1-X2-X5）]×0.2（9）

約束條件主要針對優(yōu)化變量的取值范圍，以及有限元模型中的驗算結(jié)果，其各數(shù)值取值范圍為：

xi∈x0i-5，x0i+5

σtp≤0.4ftk

Vd≤1/γ00.51×10-3fcu，kbh

fp≤L/600（10）

擬定參數(shù)設(shè)置如下：種群大小為80，最大迭代次數(shù)400，最小頻率和最大頻率分別為0和1，最大響度1，響度衰減因子0.8，最小脈沖速率0.4，脈沖增加因子0.8，ωi=1，ωe=1，ωmax=1，學(xué)習(xí)因子cmin=1.2。優(yōu)化流程采納2.3節(jié)中算法流程。

4.2優(yōu)化結(jié)果分析

依據(jù)改良蝙蝠算法編制相應(yīng)程序，結(jié)合MidasCivil有限元分析軟件，對于鋼筋混凝土箱梁截面參數(shù)進展優(yōu)化設(shè)計，結(jié)果如表1所示。

依據(jù)表1數(shù)據(jù)可以看出，優(yōu)化后總體尺寸數(shù)據(jù)小于原設(shè)計尺寸，優(yōu)化后總體重量相較于原設(shè)計削減了約9.2%。優(yōu)化后，一次成橋計算狀態(tài)下，中跨最大位移為8.21mm，邊跨最大位移為6.23mm，優(yōu)化前后位移比照如圖3所示，可以看出優(yōu)化后中跨、邊跨位移都有所增長，說明構(gòu)造整體重量雖有減輕，剛度卻有肯定減弱，但是總體來看位移仍是可以承受的程度。

在優(yōu)化后需要針對橋梁整體構(gòu)造的承載力量進展分析，模型采納雙車道布置，荷載等級為汽車-超20，依據(jù)橋涵通用設(shè)計標(biāo)準(zhǔn)進展荷載組合，得出其承載力結(jié)果如圖4所示。

依據(jù)圖4承載力驗算結(jié)果可知，抗彎、抗裂驗算均有充裕值，說明橋梁截面仍有優(yōu)化余地，此處主要受限于行車需要，截面寬度并沒有作為優(yōu)化變量。由圖4（b）可以看出斜截面抗剪承載力充裕缺乏，這是由于構(gòu)造抗剪不僅與抗剪鋼筋布置有關(guān)，還與截面腹板厚度有關(guān)，所以在進展優(yōu)化時，受限于抗剪承載力需要，腹板厚度不行削減過多。

5結(jié)語

本文依據(jù)蝙蝠算法根本原理，引入了反向解以及權(quán)重參數(shù)對算法進展改良，并依托某鋼筋混凝土橋梁工程背景，對于其截面參數(shù)進展優(yōu)化設(shè)計，得出如下結(jié)論：

（1）改良蝙蝠算法彌補了原始算法種群多樣性缺乏、易于陷入局部最優(yōu)解的缺陷，使得算法收斂速度更快，結(jié)果更為精確。

（2）運用改良蝙蝠算法后，箱梁尺寸得以優(yōu)化，總體構(gòu)造重量削減了9.2%，位移邊、中跨總體增加較小，構(gòu)造剛度有所降低。

（3）優(yōu)化后的橋梁整體構(gòu)造安全牢靠，使用階段承載力滿意要求，同時在同類橋梁優(yōu)化過程中需留意腹板厚度的選取，滿意構(gòu)造抗剪要求。

[1]劉明慧.鋼板-混凝土組合梁橋截面優(yōu)化討論[D].西安：西安科技大學(xué)，2023.

[2]魯業(yè)紅.橋梁設(shè)計中多目標(biāo)模糊優(yōu)化求解分析[J].大路工程，2023，43（6）：117-120.

[3]孫潔，劉磊，彭益.波浪鋼腹板簡支結(jié)合箱梁的截面優(yōu)化[J].鐵道標(biāo)準(zhǔn)設(shè)計，2023（4）：64-65，72.

[4]燕松波，李偉平，蔡鎖德.PC斜拉橋主梁箱形截面優(yōu)化設(shè)計討論[J].大路與汽運，2023（3）：154-157，186.

[5]YangXS.ANewMetaheuristicBat-InspiredAlgorithm[J].StudiesinComputationalIntellcgence，2023，284：65-74.

[6]龔雪嬌，郝東光，朱瑞金.基于改良蝙蝠算法的水火電力系統(tǒng)短期優(yōu)化調(diào)度[J].水力發(fā)電，2023，46（8）：84-87，91.

[7]閻震.基于改良蝙蝠算法的裝配式混凝土構(gòu)造優(yōu)化討論[D].邯鄲：河北工程大學(xué)，2023.

[8]邵世斌.鋼筋混凝土箱型橋梁斷面優(yōu)化設(shè)計要點分析[J].工程建立與設(shè)計，2023（7）：155-157.

猜你喜愛響度箱梁種群調(diào)頻播送響度掌握的方法及技巧科技傳播(2023年5期)2023-04-07由種群增長率反向分析種群數(shù)量的變化中學(xué)生物學(xué)(2023年8期)2023-03-01種群數(shù)量變化中的“率”和“速率”中學(xué)生物學(xué)(202