基于有限元與深度學(xué)習(xí)的大跨度橋梁智能設(shè)計(jì)方法

■付迎春 張 龍 范一娜 趙瓊琳

（山東省交通規(guī)劃設(shè)計(jì)院集團(tuán)有限公司，濟(jì)南 250000）

斜拉橋、 懸索橋等大跨度橋梁因其力學(xué)性能好、跨越能力強(qiáng)、外形美觀而被廣泛應(yīng)用[1]。 此類橋梁大多應(yīng)用于跨海等重要程度高的區(qū)域，也為結(jié)構(gòu)的設(shè)計(jì)提出了更高的要求。 同時(shí)，大型橋梁結(jié)構(gòu)需建立大量的設(shè)計(jì)參數(shù)和力學(xué)參數(shù)，為分析計(jì)算帶來了較大的難度。 針對(duì)橋梁的設(shè)計(jì)，眾多專家學(xué)者提出了可靠的設(shè)計(jì)方法和理念。 在橋梁結(jié)構(gòu)的設(shè)計(jì)過程中，涉及大量的設(shè)計(jì)參數(shù)和力學(xué)參數(shù)，如何有效地融合這兩類關(guān)鍵參數(shù)，實(shí)現(xiàn)結(jié)構(gòu)的快速精準(zhǔn)化分析是橋梁設(shè)計(jì)亟待解決的問題[2]。 本研究提出了基于有限元與深度學(xué)習(xí)的大跨橋梁智能設(shè)計(jì)方法，將有限元分析與神經(jīng)網(wǎng)絡(luò)融合得到了設(shè)計(jì)參數(shù)評(píng)估流程，進(jìn)一步分析力學(xué)參數(shù)的可行性獲取結(jié)構(gòu)的合理設(shè)計(jì)參數(shù)，并以斜拉橋?yàn)檠芯繉?duì)象，驗(yàn)證了本研究方法的可行性。

1 橋梁設(shè)計(jì)的基本原則

隨著國(guó)家經(jīng)濟(jì)建設(shè)的不斷發(fā)展，橋梁結(jié)構(gòu)的建造技術(shù)也在不斷提升。 其中斜拉橋、懸索橋等大跨度橋梁也在發(fā)揮著重要的作用。 大跨橋梁作為連接多地區(qū)經(jīng)濟(jì)、文化交流的交通樞紐，必須在結(jié)構(gòu)設(shè)計(jì)的安全性和耐久性方面給予重視[3]。

1.1 結(jié)構(gòu)影響因素分析

橋梁的設(shè)計(jì)與施工水平直接影響著使用階段的舒適度，結(jié)構(gòu)設(shè)計(jì)參數(shù)的選擇是保證橋梁安全性和耐久性的關(guān)鍵。 影響橋梁工程結(jié)構(gòu)的主要因素為構(gòu)件的預(yù)應(yīng)力不足、截面裂縫問題和耐久性考慮不充分，3 類因素對(duì)橋梁結(jié)構(gòu)的影響見表1。

表1 3 類因素對(duì)橋梁結(jié)構(gòu)的影響

1.2 結(jié)構(gòu)可靠度設(shè)計(jì)

根據(jù)橋梁設(shè)計(jì)規(guī)范[4]，國(guó)內(nèi)常用的耐久性分析方法為基于可靠度理論的極限狀態(tài)方法。 在該方法中將影響結(jié)構(gòu)可靠性的因素作為隨機(jī)變量，統(tǒng)計(jì)結(jié)構(gòu)所承受的各類荷載。 根據(jù)結(jié)構(gòu)形式和材料特性計(jì)算結(jié)構(gòu)的抗力，最終運(yùn)用數(shù)理統(tǒng)計(jì)理論分析隨機(jī)變量的特性，獲取橋梁結(jié)構(gòu)的總體可靠度。 據(jù)此可推導(dǎo)橋梁結(jié)構(gòu)的失效方程，見式（1）。

公式中，Z（t）為結(jié)構(gòu)的可靠性，R（t）為結(jié)構(gòu)抗力隨機(jī)過程；S（t）為結(jié)構(gòu)荷載隨機(jī)過程。

2 橋梁智能設(shè)計(jì)方法

依據(jù)橋梁設(shè)計(jì)的基本原則，要將設(shè)計(jì)參數(shù)與力學(xué)參數(shù)融合以保證結(jié)構(gòu)在各類影響因素作用下的可靠性。 傳統(tǒng)的設(shè)計(jì)方法通過有限元計(jì)算不同設(shè)計(jì)參數(shù)對(duì)應(yīng)的力學(xué)參數(shù)，需要消耗大量的時(shí)間成本。 本研究將有限元計(jì)算與神經(jīng)網(wǎng)絡(luò)融合，有效降低了計(jì)算成本，同時(shí)也為類似的結(jié)構(gòu)設(shè)計(jì)提供了參考。

2.1 有限元分析與深度學(xué)習(xí)的融合

在橋梁結(jié)構(gòu)的設(shè)計(jì)過程中，根據(jù)有限元仿真的結(jié)果，提取設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的對(duì)應(yīng)關(guān)系，融合深度學(xué)習(xí)算法挖掘數(shù)據(jù)間的關(guān)聯(lián)規(guī)則。 針對(duì)大跨橋梁結(jié)構(gòu)的特點(diǎn)在結(jié)構(gòu)的設(shè)計(jì)階段形成設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的關(guān)聯(lián)機(jī)制，具體表述見式（2）。

公式中， f （a1，a2， …，am） 表示設(shè)計(jì)參數(shù)的集合體，a1，a2，…，am分別表示具體的設(shè)計(jì)參數(shù)，如拉索的尺寸、主梁的尺寸、拉索的間距等；g（b1，b2，…，bn）表示力學(xué)參數(shù)的集合體，b1，b2，…，bn分別表示具體的力學(xué)參數(shù)，如拉索的索力、應(yīng)力、結(jié)構(gòu)的位移等；?R表示根據(jù)結(jié)構(gòu)的技術(shù)標(biāo)準(zhǔn)及數(shù)據(jù)關(guān)聯(lián)規(guī)程則，不同設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的對(duì)應(yīng)關(guān)系。 在關(guān)聯(lián)規(guī)則的驅(qū)動(dòng)下可以實(shí)現(xiàn)對(duì)結(jié)構(gòu)的智能化找形，通過改變?cè)O(shè)計(jì)參數(shù)獲取結(jié)構(gòu)的力學(xué)參數(shù)，判斷結(jié)構(gòu)模型的可行性。

本研究將有限元分析的結(jié)果形成樣本集，輸入神經(jīng)網(wǎng)絡(luò)中進(jìn)行數(shù)據(jù)的智能化分析。 在有限元分析過程中形成設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的關(guān)聯(lián)樣本。 以該樣本為依據(jù)，將設(shè)計(jì)參數(shù)作為神經(jīng)網(wǎng)絡(luò)的輸入層，將力學(xué)參數(shù)作為神經(jīng)網(wǎng)絡(luò)的輸出層。 最終由力學(xué)參數(shù)評(píng)估設(shè)計(jì)的合理性。 有限元分析與深度學(xué)習(xí)的融合見圖1。

圖1 有限元分析與深度學(xué)習(xí)的融合

3 智能設(shè)計(jì)流程

本研究利用BP 神經(jīng)網(wǎng)絡(luò)進(jìn)行參數(shù)關(guān)聯(lián)規(guī)則的捕捉。 結(jié)構(gòu)找形過程中，采用有限元軟件進(jìn)行仿真分析。 首先選擇結(jié)構(gòu)的設(shè)計(jì)參數(shù)，調(diào)整結(jié)構(gòu)的找形工況，獲取對(duì)應(yīng)的力學(xué)參數(shù)。 結(jié)合結(jié)構(gòu)技術(shù)規(guī)程，最終通過力學(xué)參數(shù)判斷結(jié)構(gòu)找形的可行性。 在本研究中，根據(jù)工程經(jīng)驗(yàn)選取關(guān)鍵的設(shè)計(jì)參數(shù)及其具體的數(shù)值范圍。 在有限元分析過程中形成大量的樣本。樣本中包含結(jié)構(gòu)的設(shè)計(jì)參數(shù)和力學(xué)參數(shù)， 為基于BP 神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)智能找形模型的搭建提供了數(shù)據(jù)支撐。 根據(jù)大跨橋梁結(jié)構(gòu)的特點(diǎn)，本研究中所選取的結(jié)構(gòu)設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的集合體具體表述為式（3）和式（4）。

公式中，DP 表示結(jié)構(gòu)設(shè)計(jì)參數(shù)的集合，Ns表示拉索的間距，Ds表示斜拉索的直徑，Rl表示橋梁的邊跨中比，HL表示主梁的腹板長(zhǎng)度，Bl表示主梁的翼緣長(zhǎng)度，MP 表示結(jié)構(gòu)力學(xué)參數(shù)的集合，Dp表示結(jié)構(gòu)的豎向位移，σs表示拉索的應(yīng)力設(shè)計(jì)參數(shù)的選擇為仿真模型的搭建提供了基礎(chǔ)，力學(xué)參數(shù)的選擇為結(jié)構(gòu)找形的可行性提供了依據(jù)。

根據(jù)BP 神經(jīng)網(wǎng)絡(luò)的架構(gòu)和結(jié)構(gòu)設(shè)計(jì)的特征，將各類設(shè)計(jì)參數(shù)歸為神經(jīng)網(wǎng)絡(luò)的輸入層，各類力學(xué)參數(shù)歸為輸出層。 通過有限元中設(shè)置設(shè)計(jì)參數(shù)與找形工況生成對(duì)應(yīng)的力學(xué)參數(shù)。 將有限元分析得出的力學(xué)參數(shù)作為判斷神經(jīng)網(wǎng)絡(luò)準(zhǔn)確性的重要依據(jù)。 根據(jù)設(shè)計(jì)參數(shù)，輸入模型中獲取預(yù)測(cè)的力學(xué)參數(shù)。 通過對(duì)比兩種途徑獲取的力學(xué)參數(shù)，不斷調(diào)整預(yù)測(cè)模型的權(quán)重，最終形成高效可靠的找形預(yù)測(cè)模型。

4 案例應(yīng)用

4.1 基本情況

在橋梁智能設(shè)計(jì)方法的驅(qū)動(dòng)下，本研究以斜拉橋?yàn)榉治鰧?duì)象進(jìn)行了案例驗(yàn)證。 本研究所選取的斜拉鋼索橋跨度為200 m，見圖2。

圖2 斜拉鋼索橋的結(jié)構(gòu)形式

根據(jù)設(shè)計(jì)規(guī)范和本研究所提出的智能設(shè)計(jì)方法，在設(shè)計(jì)過程中，主要考慮的設(shè)計(jì)參數(shù)及其取值范圍見表2，本研究重點(diǎn)考察結(jié)構(gòu)的豎向位移是否超限，根據(jù)橋梁設(shè)計(jì)規(guī)范，最大位移的不得超過跨度的1/400。 根據(jù)式（1）計(jì)算最不利荷載工況下橋體的最大豎向位移。

表2 設(shè)計(jì)參數(shù)及其取值范圍

4.2 結(jié)構(gòu)受力性能預(yù)測(cè)

在有限元中進(jìn)行設(shè)計(jì)參數(shù)的輸入，并在荷載工況下獲取結(jié)構(gòu)的豎向位移。 有限元中設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的計(jì)算見圖3。根據(jù)位移云圖可知，結(jié)構(gòu)的最大豎向位移發(fā)生于邊索與橋墩之間，應(yīng)加強(qiáng)該部位的剛度。

圖3 有限元中設(shè)計(jì)參數(shù)與力學(xué)參數(shù)的計(jì)算

根據(jù)設(shè)計(jì)參數(shù)的取值，在荷載作用下形成了大量的樣本。 將設(shè)計(jì)參數(shù)與力學(xué)參數(shù)形成的關(guān)聯(lián)樣本作為數(shù)據(jù)集訓(xùn)練神經(jīng)網(wǎng)絡(luò)。 在訓(xùn)練找形預(yù)測(cè)模型過程中，當(dāng)mse≤6.5×10-4時(shí)，模型停止迭代。 本研究分別探析了5 類設(shè)計(jì)參數(shù)對(duì)豎向位移的影響。 在訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)中獲取了不同設(shè)計(jì)參數(shù)對(duì)應(yīng)的豎向位移值，同時(shí)在有限元模型中獲取相同設(shè)計(jì)參數(shù)對(duì)應(yīng)的力學(xué)參數(shù)。 部分樣本中豎向位移的預(yù)測(cè)值與實(shí)際值的對(duì)比見圖4。

圖4 部分樣本中豎向位移的預(yù)測(cè)值與實(shí)際值的對(duì)比

在分析過程中發(fā)現(xiàn)，本研究提出的設(shè)計(jì)方法具有較高的擬合度。 由于分析過程中將物理問題轉(zhuǎn)化為數(shù)學(xué)問題，減少了邊界條件設(shè)置、工況分析等步驟，顯著降低了計(jì)算時(shí)間。 通過數(shù)據(jù)分析對(duì)比，本研究提出的智能設(shè)計(jì)方法在計(jì)算效率方面得到了顯著提升。 處理相同數(shù)量的樣本，對(duì)比有限元分析的時(shí)間（TE）和神經(jīng)網(wǎng)絡(luò)分析的時(shí)間（TM）獲取計(jì)算時(shí)間的節(jié)省量（θ），見式（5）。

根據(jù)智能設(shè)計(jì)方法獲取了保證結(jié)構(gòu)位移不超限的設(shè)計(jì)參數(shù)集合。 拉索的間距為12 m，斜拉索的直徑為150 mm，橋梁的邊跨中比為0.8，主梁的腹板長(zhǎng)度為2 000 mm，主梁的翼緣長(zhǎng)度為2 000 mm。

4.3 參數(shù)化分析

本研究在神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)模型中進(jìn)行了結(jié)構(gòu)設(shè)計(jì)參數(shù)的分析，分別獲取了多種設(shè)計(jì)參數(shù)取值下的力學(xué)響應(yīng)。 根據(jù)表3 中5 種設(shè)計(jì)參數(shù)的取值，分析了各類參數(shù)對(duì)結(jié)構(gòu)豎向位移的影響。 通過參數(shù)化分析可以有效獲取提高結(jié)構(gòu)安全性能的措施。

4.3.1 拉索的間距

根據(jù)規(guī)范要求，采用密距布置方案，拉索間距的取值范圍是8～24 m。 在本研究中，控制其他設(shè)計(jì)參數(shù)不變，設(shè)置不同的拉索間距，豎向位移的變化見圖5。由圖5 可知，適當(dāng)降低拉索的間距可以保證結(jié)構(gòu)豎向位移不超過限制。 然而隨著間距的降低，對(duì)豎向位移的影響明顯降低。 因此，斜拉索的間距應(yīng)保持在一定的范圍內(nèi)，過于密集對(duì)豎向位移的影響不明顯且提高了建造成本。

圖5 拉索間距對(duì)豎向位移的影響

4.3.2 拉索的直徑

本研究拉索間距的取值范圍是50～300 mm。 在本研究中，控制其他設(shè)計(jì)參數(shù)不變，設(shè)置不同的拉索直徑，豎向位移的變化見圖6。 由圖6 可知，適當(dāng)增加拉索的直徑可以保證結(jié)構(gòu)豎向位移不超過限制。 然而隨著拉索直徑的增加，對(duì)豎向位移的影響明顯降低。 因此，斜拉索的直徑應(yīng)保持在一定的范圍內(nèi)，尺寸過大對(duì)豎向位移的影響不明顯且提高了建造成本。

圖6 拉索直徑對(duì)豎向位移的影響

4.3.3 邊跨中比

本研究邊跨中比的取值范圍是0.5～1.0。在本研究中，控制其他設(shè)計(jì)參數(shù)不變，設(shè)置不同的邊跨中比，豎向位移的變化見圖7。 由圖7 可知，適當(dāng)減小邊跨中比可以保證結(jié)構(gòu)豎向位移不超過限制。 然而隨著邊跨中比的減小， 對(duì)豎向位移的影響明顯降低。 因此，邊跨中比應(yīng)保持在一定的范圍內(nèi)，比例過小對(duì)豎向位移的影響不明顯且出現(xiàn)較長(zhǎng)拉索提高了建造成本。

圖7 邊跨中比對(duì)豎向位移的影響

4.3.4 主梁截面

本研究主梁截面由腹板長(zhǎng)度和翼緣寬度表示。在本研究中，控制其他設(shè)計(jì)參數(shù)不變，設(shè)置不同的主梁截面，豎向位移的變化見圖8。 由圖8 可知，主梁的截面提供了橋體的剛度，對(duì)豎向位移的影響最大。 因此，當(dāng)橋梁變形超限時(shí)，增大主梁的截面尺寸是最有效的維護(hù)措施。

圖8 主梁截面對(duì)豎向位移的影響

5 結(jié)論

大型橋梁的設(shè)計(jì)過程中，由于涉及眾多構(gòu)件和參數(shù)，存在著設(shè)計(jì)的效率低和精度不足的問題。 本研究提出了基于有限元與深度學(xué)習(xí)的大跨橋梁智能設(shè)計(jì)方法。 在研究過程中，明確了橋梁設(shè)計(jì)的基本原則，分析了結(jié)構(gòu)性能的影響因素，并提出了結(jié)構(gòu)耐久性設(shè)計(jì)的關(guān)鍵指標(biāo)。 通過輸入設(shè)計(jì)參數(shù)可以精準(zhǔn)快速地獲取結(jié)構(gòu)的力學(xué)響應(yīng)，為結(jié)構(gòu)的設(shè)計(jì)提供了新思路。 本研究以斜拉橋的設(shè)計(jì)為例進(jìn)行了研究方法的應(yīng)用，所提出的智能設(shè)計(jì)方法相較于傳統(tǒng)有限元方法效率得到顯著提升，且具有較高的計(jì)算精度。