周藝，邱勇，何文云，郭金楠，阮艷平

（1云南農(nóng)業(yè)大學(xué) 水利學(xué)院，云南昆明 650201；2云南省玉溪市水利電力勘測(cè)設(shè)計(jì)院，云南玉溪 653100）

0 引言

傳統(tǒng)工程結(jié)構(gòu)的設(shè)計(jì)主要依靠工程人員的經(jīng)驗(yàn)和力學(xué)概念，并未形成系統(tǒng)的方法[1]。蔡新等[2]提及早期優(yōu)化設(shè)計(jì)多采用基于直覺(jué)的準(zhǔn)則法，缺少嚴(yán)格的理論依據(jù)，求得的待定解不一定是真實(shí)解。隨著計(jì)算機(jī)的發(fā)展，基于數(shù)學(xué)和力學(xué)方法的結(jié)構(gòu)優(yōu)化設(shè)計(jì)，已經(jīng)由尺寸優(yōu)化發(fā)展到形狀優(yōu)化、拓?fù)鋬?yōu)化等階段[2]。

拓?fù)鋬?yōu)化能在工程結(jié)構(gòu)設(shè)計(jì)初始階段為設(shè)計(jì)者提供一個(gè)新穎高效的概念性設(shè)計(jì)，使結(jié)構(gòu)在布局上實(shí)現(xiàn)最優(yōu)方案[3]。陳艾榮等[4]通過(guò)實(shí)例，介紹了應(yīng)用拓?fù)鋬?yōu)化技術(shù)進(jìn)行橋梁結(jié)構(gòu)找型的方法，并展示了結(jié)構(gòu)衍化過(guò)程及優(yōu)化結(jié)果；張宏等[3]基于密度懲罰法，提出了一種連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化問(wèn)題的優(yōu)化準(zhǔn)則算法，在橋梁概念設(shè)計(jì)中得到了理想的優(yōu)化結(jié)果。

就大學(xué)生結(jié)構(gòu)設(shè)計(jì)競(jìng)賽而言，夏雨等[5]對(duì)紙質(zhì)結(jié)構(gòu)模型進(jìn)行了極限承載力的實(shí)驗(yàn)及結(jié)構(gòu)優(yōu)化；陳慶軍等[6]提出了如何在一個(gè)限定空間內(nèi)，在滿足強(qiáng)度及剛度的前提下，找出一個(gè)最優(yōu)或者較優(yōu)的空間桿系結(jié)構(gòu)；王文利等[7]在結(jié)構(gòu)設(shè)計(jì)初始階段經(jīng)過(guò)理論分析和不斷的模型試驗(yàn)，嘗試了近50個(gè)模型之后得到滿足賽題要求的結(jié)構(gòu)。

總的來(lái)說(shuō)，拓?fù)鋬?yōu)化能夠在概念設(shè)計(jì)階段通過(guò)對(duì)結(jié)構(gòu)進(jìn)行整體找型，避免橋梁設(shè)計(jì)和制作階段的盲目試載，節(jié)約大量時(shí)間并控制成本。

1 賽題簡(jiǎn)介

賽題所要求的橋梁結(jié)構(gòu)設(shè)計(jì)，應(yīng)能夠承受分散作用的豎向集中荷載以及橋面移動(dòng)荷載，亦即在限定的空間尺寸和加載條件下實(shí)現(xiàn)變參數(shù)（支座頂面標(biāo)高、橋下凈空頂標(biāo)高和豎向加載荷載）橋梁結(jié)構(gòu)預(yù)定功能。

1.1 結(jié)構(gòu)模型設(shè)計(jì)要求

模型設(shè)計(jì)時(shí)，需要考慮的空間尺寸限制條件（圖1）如下：（1）平面上橋梁結(jié)構(gòu)設(shè)計(jì)須規(guī)避給定的4個(gè)陰影區(qū)域；（2）為保證通航要求，橋下凈空頂部標(biāo)高H取值范圍為-150～-50mm，以50mm階梯取值；（3）②軸支座頂面標(biāo)高V1=-160mm；（4）③軸支座頂面標(biāo)高V2取值范圍為-160～140mm，以75mm階梯取值。

圖1 加載裝置圖

模型加載過(guò)程中，B軸和B1B2交點(diǎn)的位移測(cè)點(diǎn)最大變形不得超過(guò)10mm。結(jié)構(gòu)所使用的竹材力學(xué)指標(biāo)見(jiàn)表1，模型設(shè)計(jì)或制作可供選擇的竹材幾何截面尺寸見(jiàn)表2。

表1 竹材力學(xué)指標(biāo)

表2 竹材規(guī)格

1.2 加載方式

8個(gè)豎向加載點(diǎn)處的荷載取值范圍為40～130N（取10N的倍數(shù)），各加載點(diǎn)荷載值不能重復(fù)。

（1）一級(jí)加載，8個(gè)加載點(diǎn)分別懸掛各自的待定荷載。

（2）二級(jí)加載，包括兩個(gè)步驟：步驟一，保持C1、C2加載點(diǎn)靜載不變，左側(cè)加載點(diǎn)或右側(cè)加載點(diǎn)中任選一個(gè)加載點(diǎn)將該加載點(diǎn)的所有荷載值轉(zhuǎn)移至另一側(cè)的任一加載點(diǎn)上；步驟二，將第一步移入荷載點(diǎn)上的所有荷載值，全部轉(zhuǎn)移至該加載點(diǎn)的B軸對(duì)稱點(diǎn)，或者移至第一步移出荷載點(diǎn)的B軸對(duì)稱點(diǎn)上。

（3）三級(jí)加載，保持上一級(jí)靜載作用，施加移動(dòng)荷載50N。

2 參數(shù)取值與方案設(shè)計(jì)

待定參數(shù)的加入，不僅賦予橋梁結(jié)構(gòu)設(shè)計(jì)競(jìng)賽更多的靈活性，還提供了方案選擇的諸多可能。

2.1 支座頂面標(biāo)高取值

支座頂面標(biāo)高V2變參數(shù)設(shè)計(jì)取值范圍為-160～140mm，按照75mm階梯進(jìn)行 取 值，有-160mm、-85mm、-10mm、65mm、140mm等5種情況。前3種支座標(biāo)高均低于橋面，與橋身和支座相連的桿件主要受壓；后2種支座高度均高于橋面，與橋身和支座相連的桿件主要受拉。為充分考慮與橋身和支座相連桿件的受壓和受拉情況，分別選取支座標(biāo)高低于橋面和高于橋面的三組取值-160mm、-85mm和140mm進(jìn)行對(duì)比分析。

2.2 橋下凈空頂部標(biāo)高取值

橋下凈空頂部標(biāo)高H取值范圍為-150～-50mm，以50mm階梯取值，有-150mm、-100mm和-50mm等3種情況。

2.3 加載取值與方案設(shè)計(jì)

考慮變參數(shù)橋梁結(jié)構(gòu)受不同荷載組合的影響，橋梁結(jié)構(gòu)在滿足第一級(jí)加載的情況下，需要針對(duì)荷載施加方案中第二級(jí)加載方案進(jìn)行探討。經(jīng)過(guò)分析，選取以下兩組工況（表3）：工況一和工況二的一級(jí)加載均為最大荷載組合，其中工況一的荷載基本對(duì)稱布置，最大荷載130N布置在自由端，承擔(dān)一定的扭矩，工況二的荷載偏心布置，最大荷載130N布置在自由端，加大偏心影響；兩種工況下二級(jí)加載的移動(dòng)荷載點(diǎn)的選擇均考慮最不利荷載（扭矩）組合影響；在此基礎(chǔ)上施加移動(dòng)荷載50N完成三級(jí)加載。

表3 加載點(diǎn)取值方案(A1～D2）（單位：N）

3 結(jié)構(gòu)優(yōu)化

變參數(shù)橋梁結(jié)構(gòu)優(yōu)化一般可分為3個(gè)層次：拓?fù)鋬?yōu)化、形狀優(yōu)化、尺寸優(yōu)化[8]，分別對(duì)應(yīng)實(shí)際設(shè)計(jì)中的三個(gè)不同階段：概念設(shè)計(jì)、基本設(shè)計(jì)和詳細(xì)設(shè)計(jì)。拓?fù)鋬?yōu)化是在給定的設(shè)計(jì)空間中，尋求滿足一定約束條件的最優(yōu)構(gòu)件或材料分布形式，優(yōu)化的三要素包括設(shè)計(jì)變量、目標(biāo)函數(shù)和約束條件，根據(jù)不同的問(wèn)題建立相應(yīng)的數(shù)學(xué)模型，尋找求解設(shè)計(jì)變量的合適方法[9]。

通過(guò)拓?fù)鋬?yōu)化可以對(duì)變參數(shù)橋梁結(jié)構(gòu)進(jìn)行整體找型，使得設(shè)計(jì)部分能夠不受已有橋梁結(jié)構(gòu)體系的思維束縛，從而得出不拘一格的結(jié)構(gòu)形式[4]；形狀優(yōu)化是用來(lái)發(fā)掘工程系統(tǒng)構(gòu)件的合理內(nèi)外邊界形狀的；尺寸優(yōu)化則是以結(jié)構(gòu)桿件外形或者孔洞形狀為優(yōu)化對(duì)象，力求得到滿足工程安全要求的最輕結(jié)構(gòu)。

3.1 優(yōu)化模型建立

考慮到橋面板上部空間范圍內(nèi)需要滿足移動(dòng)荷載要求，為了賦予結(jié)構(gòu)更多的可選設(shè)計(jì)空間，模型上部結(jié)構(gòu)采用左右對(duì)稱的板殼結(jié)構(gòu)，下部結(jié)構(gòu)采用實(shí)體結(jié)構(gòu)，如圖2所示。

圖2 結(jié)構(gòu)初始模型圖

3.2 結(jié)構(gòu)拓?fù)鋬?yōu)化

通過(guò)拓?fù)鋬?yōu)化得到初始結(jié)構(gòu)最不利荷載的位置為D1、D2點(diǎn)處（圖3），經(jīng)過(guò)找型后的模型以板面結(jié)構(gòu)和桿件結(jié)構(gòu)為主：橋面板采用桿件結(jié)構(gòu)，上部結(jié)構(gòu)采用板殼，下部結(jié)構(gòu)采用板殼和桿件結(jié)構(gòu)相結(jié)合的形式。有限元模型的單元類(lèi)型采用梁?jiǎn)卧?，網(wǎng)格數(shù)量為677129，殼單元和實(shí)體單元均采用剛性連接。

圖3 結(jié)構(gòu)找型后模型圖

3.2.1 支座頂面標(biāo)高下的拓?fù)鋬?yōu)化

考慮結(jié)構(gòu)保留率50%時(shí)，對(duì)不同支座頂面標(biāo)高=-160mm、-85mm和140mm進(jìn)行拓?fù)鋬?yōu)化，得到的刪除率依次為55.95%、55.54%、52.67%（圖4）?？紤]到支座頂面標(biāo)高對(duì)結(jié)構(gòu)重量的影響以及冗余部分刪除情況，支座頂面標(biāo)高V2取值-85mm時(shí)，優(yōu)化后得到的結(jié)構(gòu)質(zhì)量相對(duì)更小。

圖4 支座頂面標(biāo)高優(yōu)化后結(jié)構(gòu)圖

3.2.2 橋下凈空拓?fù)鋬?yōu)化

當(dāng)支座頂面標(biāo)高取值V2=-85mm時(shí)，橋下凈空取值只能考慮H=-50mm（H=-100mm和H=-150mm無(wú)法滿足通航要求）。由于右側(cè)支座兩側(cè)荷載原因，橋下凈空拓?fù)鋬?yōu)化結(jié)果見(jiàn)圖5。

圖5 橋下凈空優(yōu)化后結(jié)構(gòu)圖

3.2.3 加載取值拓?fù)鋬?yōu)化

綜合考慮對(duì)稱荷載和不對(duì)稱荷載作用情況，保留率為50%時(shí)，上部結(jié)構(gòu)能夠刪除冗余部分，甚至形成桿件，但部分區(qū)域仍然存在板殼部分或者桿件截面面積還具有優(yōu)化的空間。經(jīng)過(guò)多次拓?fù)鋬?yōu)化后，得到滿足要求的不同方案結(jié)構(gòu)保留率僅為26.50%、27.90%（圖6），此時(shí)對(duì)稱荷載作用下的結(jié)構(gòu)質(zhì)量為2.205kg，不對(duì)稱荷載作用下的結(jié)構(gòu)質(zhì)量為2.286kg。相比初始結(jié)構(gòu)，優(yōu)化后的質(zhì)量下降84.86%、83.25%。

圖6 加載取值優(yōu)化后結(jié)構(gòu)圖

綜合考慮，選取工況二進(jìn)行結(jié)構(gòu)設(shè)計(jì)能夠同時(shí)滿足工況一加載要求，故后續(xù)優(yōu)化設(shè)計(jì)僅對(duì)工況二進(jìn)行討論。

3.3 結(jié)構(gòu)形狀優(yōu)化

由于實(shí)際工程的復(fù)雜性，拓?fù)鋬?yōu)化的結(jié)果與實(shí)際需求相比，還遠(yuǎn)遠(yuǎn)不能達(dá)到工程應(yīng)用要求[3]。形狀優(yōu)化通過(guò)調(diào)整結(jié)構(gòu)內(nèi)邊界形狀來(lái)改善結(jié)構(gòu)性能，其主要用來(lái)發(fā)掘構(gòu)件的合理內(nèi)外邊界形狀[2]?；诮Y(jié)構(gòu)所用材料要求，形狀優(yōu)化后的結(jié)構(gòu)圖均為桿件，如圖7所示。

從圖7可以看出，形狀優(yōu)化后的結(jié)構(gòu)桿件集中在右側(cè)支座附近，和結(jié)構(gòu)加載的主要受力情況吻合。

圖7 形狀優(yōu)化后結(jié)構(gòu)圖

3.4 結(jié)構(gòu)尺寸優(yōu)化

尺寸優(yōu)化以構(gòu)件外形為優(yōu)化對(duì)象，對(duì)構(gòu)件的截面尺寸進(jìn)行調(diào)整，以適應(yīng)不同加載需求?？紤]不對(duì)稱荷載作用，按照表2所給竹材規(guī)格對(duì)圖7的結(jié)構(gòu)進(jìn)行尺寸優(yōu)化，得到圖8所示的結(jié)構(gòu)，優(yōu)化流程具體如圖9所示。

圖8 尺寸優(yōu)化后結(jié)構(gòu)圖

圖9 優(yōu)化流程圖

結(jié)構(gòu)剛度的優(yōu)化結(jié)果見(jiàn)表4。

表4 結(jié)構(gòu)剛度優(yōu)化結(jié)果

表4反映了結(jié)構(gòu)優(yōu)化過(guò)程中剛度的變化：結(jié)構(gòu)從初始狀態(tài)的實(shí)體結(jié)構(gòu)到支座標(biāo)高下的拓?fù)鋬?yōu)化過(guò)程中，模型剛度在增大（變形減?。坏珡牟贿B續(xù)的拓?fù)浣Y(jié)構(gòu)到竹材桿件結(jié)構(gòu)，由于結(jié)構(gòu)質(zhì)量繼續(xù)減少，桿件結(jié)構(gòu)的剛度也出現(xiàn)減小（變形增大），故需要進(jìn)一步驗(yàn)證最終優(yōu)化后的結(jié)構(gòu)是否能夠滿足承載要求。

3.5 結(jié)構(gòu)內(nèi)力及應(yīng)力應(yīng)變計(jì)算

為了判斷結(jié)構(gòu)強(qiáng)度、剛度以及穩(wěn)定性是否滿足要求，需要對(duì)拓?fù)鋬?yōu)化后的結(jié)構(gòu)（圖8）進(jìn)行內(nèi)力分析和應(yīng)力應(yīng)變計(jì)算，結(jié)果如圖10所示。

圖10 三級(jí)荷載作用下模型內(nèi)力及應(yīng)力應(yīng)變圖

其中：A桿為軸向拉伸Nmax=654.9N，B桿為軸向壓縮Nmax=656.8N；C桿為最大剪力Qmax=4.7N，D桿為最小剪力Qmin=4.4N；C桿為最大扭矩MTmax=181.3N，D桿為最小扭矩MTmin=262.5N；C桿為最大正彎矩Mmax=504.5N·mm，D桿為最大負(fù)彎矩Mmax=851.1N·mm；A桿為最大拉應(yīng)力σmax=43.4MPa，E桿為最大壓應(yīng)力σmax=29.6MPa；A2點(diǎn)處最大位移值εmax=-11.270mm。

根據(jù)以上計(jì)算結(jié)果，列出結(jié)構(gòu)在不同工況作用下的桿件最大應(yīng)力及應(yīng)變值，如表5所示。

表5 結(jié)構(gòu)應(yīng)力應(yīng)變計(jì)算結(jié)果

由表5可知，結(jié)構(gòu)在不同工況情況下，強(qiáng)度均能滿足要求。位移最大值出現(xiàn)在A2加載點(diǎn)，其位移εmax達(dá)到了-11.618mm和-11.270mm，位移測(cè)點(diǎn)的變形εmax僅為-2.642mm和-2.823mm。

對(duì)于受壓桿件結(jié)構(gòu)，即便強(qiáng)度滿足要求，但如果變形過(guò)大，也會(huì)導(dǎo)致其出現(xiàn)失穩(wěn)破壞。因此還需要結(jié)合桿件截面尺寸和內(nèi)力組合情況進(jìn)行穩(wěn)定驗(yàn)算，表6給出了橋梁結(jié)構(gòu)在工況二情況下主要桿件的允許承載力Pcr。

表6 主要桿件穩(wěn)定性驗(yàn)算結(jié)果

結(jié)合表5和表6可以得出：經(jīng)過(guò)拓?fù)鋬?yōu)化后，受力桿件和整體結(jié)構(gòu)均能滿足強(qiáng)度、剛度以及穩(wěn)定性要求。

4 結(jié)語(yǔ)

對(duì)于變參數(shù)橋梁結(jié)構(gòu)，通過(guò)拓?fù)鋬?yōu)化可以方便地對(duì)結(jié)構(gòu)進(jìn)行整體找型：（1）在滿足預(yù)定功能時(shí)左右對(duì)稱最不利荷載情況下的結(jié)構(gòu)保留率為26.50%時(shí)，優(yōu)化后的冗余部分刪除率為76.57%；（2）左右不對(duì)稱最不利荷載情況下的保留率為27.90%時(shí)，優(yōu)化后冗余部分刪除率為69.57%。亦即刪除冗余部分過(guò)程中，不同荷載分布的作用影響不大。

基于拓?fù)鋬?yōu)化結(jié)果和賽題所給材料進(jìn)行形狀優(yōu)化，得到斜拉式的橋梁上部結(jié)構(gòu)，再結(jié)合竹材規(guī)格進(jìn)行尺寸優(yōu)化得到最優(yōu)截面：（1）左右對(duì)稱最不利加載情況下，結(jié)構(gòu)模型質(zhì)量為290.2g，桿件最大拉應(yīng)力σmax=34.5MPa、最大壓應(yīng)力σmax=-27.9MPa，最大位移ε=-12.002mm（測(cè)試點(diǎn)位移滿足要求）；（2）左右不對(duì)稱最不利加載情況下，結(jié)構(gòu)模型質(zhì)量為294.2g，桿件應(yīng)力和位移值同樣能夠滿足要求。表明滿足不對(duì)稱加載的變參數(shù)橋梁結(jié)構(gòu)也能夠滿足對(duì)稱加載要求。

綜上所述，拓?fù)鋬?yōu)化能夠在橋梁概念設(shè)計(jì)階段通過(guò)整體找型功能，結(jié)合形狀和尺寸修改，得到滿足強(qiáng)度、剛度以及穩(wěn)定性要求的結(jié)構(gòu)，從而避免盲目試載，節(jié)約時(shí)間和控制成本。