橋梁結(jié)構(gòu)在大學(xué)生結(jié)構(gòu)設(shè)計競賽中的優(yōu)化分析

黃文虎

(安徽理工大學(xué)土木建筑學(xué)院，安徽 淮南 232001)

1 概述

全國大學(xué)生結(jié)構(gòu)設(shè)計競賽是由浙江大學(xué)在2005年首先發(fā)起，緊接著全國各大高校紛紛參與，距今已舉辦過十四屆。該比賽不僅提高了全國各地大學(xué)生的創(chuàng)新意識，而且比賽成果還會在各個領(lǐng)域得到充分的應(yīng)用，體現(xiàn)了良好的工程價值[1]。該比賽舉辦以來，賽題遍布工程的各個領(lǐng)域，參賽高校越來越多，難度越來越大，給賽題舉辦方、參賽高校帶來了巨大的挑戰(zhàn)，充分的展現(xiàn)了比賽的價值和魅力。

2 主要研究內(nèi)容及意義

2.1 研究內(nèi)容

本文以2021年第十四屆全國大學(xué)生結(jié)構(gòu)設(shè)計競賽賽題為例，闡述有限元軟件ANSYS在大學(xué)生結(jié)構(gòu)設(shè)計競賽中的應(yīng)用。利用軟件對結(jié)構(gòu)模型進(jìn)行計算、分析、優(yōu)化并與實際模型進(jìn)行對比分析，得出結(jié)論。

2.2 研究意義

結(jié)構(gòu)設(shè)計競賽從剛開始幾十所高校參賽，到目前的一百多所高校；從最初高層建筑結(jié)構(gòu)設(shè)計，到雙跨兩車道結(jié)構(gòu)設(shè)計、發(fā)電塔結(jié)構(gòu)設(shè)計、體育場懸挑屋蓋結(jié)構(gòu)設(shè)計……到現(xiàn)在變參數(shù)橋梁結(jié)構(gòu)設(shè)計；從最初用巴西白卡紙、蠟線、透明紙、白膠、鉛發(fā)絲線、桐木條等制作，到現(xiàn)在用竹皮、竹條、502膠水等制作；從最初的靜荷載到現(xiàn)在的動荷載；從最初的定參數(shù)設(shè)計到現(xiàn)在變參數(shù)設(shè)計；從最初的試驗到現(xiàn)在的試驗加數(shù)值模擬……無一不體現(xiàn)了結(jié)構(gòu)設(shè)計大賽的不斷發(fā)展，參賽難度也在逐一攀升，對學(xué)生以及老師都提出了巨大的挑戰(zhàn)[2]。難度的上升也導(dǎo)致競賽不能夠僅僅從試驗和經(jīng)驗就能夠很好的掌握比賽，而是對數(shù)值模擬仿真提出了要求。以往的結(jié)構(gòu)設(shè)計競賽一般采用的有限元軟件是邁達(dá)斯公司開發(fā)的Midas Gen和Midas Civil兩款軟件[3]。前者主要優(yōu)勢在于建筑結(jié)構(gòu)分析，后者更強調(diào)橋梁結(jié)構(gòu)的分析計算。容易上手，分析簡單，可以自動生成原始計算書，出圖美觀等都是它們的優(yōu)點，但它們也具有一定的缺陷，無論是Midas Gen還是Midas Civil，它們的建模過程相對復(fù)雜，甚至還得借助其他建模軟件，如CAD。然而隨著結(jié)構(gòu)設(shè)計競賽的不斷發(fā)展或許在軟件上也會有更替。ANSYS軟件在科研中應(yīng)用較多，本文就用此來對結(jié)構(gòu)進(jìn)行優(yōu)化計算。這符合競賽的發(fā)展趨勢，必然會給比賽帶來不一樣的新體驗。

3 賽題分析

上海交通大學(xué)作為本屆比賽的東道主在結(jié)構(gòu)設(shè)計競賽中一直位于第一梯隊，同時也是本次賽題的出題方，給這次比賽也帶來了不一樣的氣氛。下面就對本次賽題進(jìn)行分析，因賽題可查且十分復(fù)雜，因此只簡單描述與本文相關(guān)重點內(nèi)容。具體完整賽題需要讀者自行查找，不多贅述。部分內(nèi)容如下：

如圖1所示，賽題要求設(shè)計并制作一座橋梁結(jié)構(gòu)模型，結(jié)構(gòu)自行確定。要求該模型在能夠承受分散作用的豎向集中荷載以及橋面移動荷載的同時模型的變形滿足一定要求[4]。

其中,A1,A2,B1,B2,C1,C2,D1,D2八個豎向荷載施加點(荷載質(zhì)量為待定參數(shù))；橋下凈空頂標(biāo)高Hmin以及3軸支座標(biāo)高V2均為待定參數(shù)。

因賽題復(fù)雜，待定加載參數(shù)多樣，下面取其中一組分析。

取3軸支座頂面標(biāo)高V2=-160 mm，橋下凈空頂標(biāo)高最小值Hmin=-150 mm，各加載點要加載的砝碼重量分別為：GA1=10 kg，GA2=8 kg，GB1=7 kg，GB2=4 kg，GC1=11 kg，GC2=12 kg，GD1=6 kg，GD2=9 kg。本文只討論一級加載情況，也即八個加載點分別施加上述對應(yīng)荷載，加載完成后模型不破壞(即桿件受到的應(yīng)力滿足竹材的順紋抗壓強度30 MPa、順紋抗拉強度60 MPa的規(guī)定)以及位移測點豎向位移不超過10 mm(下文將其簡化為模型任意一點豎向最大位移不超過10 mm)。

4 模型結(jié)構(gòu)選型

在進(jìn)行結(jié)構(gòu)設(shè)計時，結(jié)構(gòu)受傳力路徑對結(jié)構(gòu)的變形、強度等力學(xué)性能影響很大。所以良好的傳力路徑，是進(jìn)行設(shè)計的前提[5]。

根據(jù)結(jié)構(gòu)設(shè)計大賽竹材抗拉壓能力強，抗彎能力相對較弱的特點[6]，綜合分析，確定橋梁模型結(jié)構(gòu)如圖2所示，并對各構(gòu)件進(jìn)行命名。

5 構(gòu)件截面的選擇

在確定好結(jié)構(gòu)整體受力體系之后，接下來最重要的就是確定組成整體結(jié)構(gòu)的單個桿件的力學(xué)特性了。

結(jié)構(gòu)構(gòu)件的截面形狀對其力學(xué)性能有很大的影響。以結(jié)構(gòu)設(shè)計大賽為例，如圖3所示，圓形截面制作相對復(fù)雜，但承受軸力的能力最大，抗彎性能則相對較差。方形截面抗彎性能相對圓形較好，而抗壓抗拉則相對較差，制作簡單，一般均采用方形截面。方形截面有很多種制作方法，圖3中第2～第5種均為方形截面。其中第3種截面制作相對于第2種而言，容易產(chǎn)生扭轉(zhuǎn)現(xiàn)象，而后兩種雖然性能大大提高，但材料用量較多，用的也相對較少[7]。結(jié)構(gòu)設(shè)計大賽追求的是在能承受荷載的前提下，材料用量最小，所以大多采用第2種矩形截面。本文采取第2種矩形截面進(jìn)行有限元模擬。

6 有限元分析及優(yōu)化

6.1 有限元法

有限元法又稱有限單元法。該方法在各個工程領(lǐng)域的應(yīng)用十分廣泛，是一種數(shù)值模擬的計算方法。

1952年一家波音公司項目小組在進(jìn)行機翼變形的相關(guān)研究中，意外發(fā)現(xiàn)用很多塊三角形構(gòu)件拼裝機翼，此時算出的機翼變形數(shù)值與實際情形大致吻合，他們果斷判斷這個發(fā)現(xiàn)不尋常，于是很快將這種計算方法寫成論文。論文中將此法命名“直接剛度法”，這就是有限元法的鼻祖。

1965年“有限元”這個詞第一次出現(xiàn)在了大家的視野。有限元法的本質(zhì)就是將實際結(jié)構(gòu)離散為有限個單元，這與直接剛度法大致相似。在滿足工程精度的前提下進(jìn)行分析[8]。

6.2 ANSYS的出現(xiàn)和發(fā)展

市面上的有限元軟件浩如煙海，Midas，Abaqus，SAP-5，ADINA等等。其中ANSYS(Analysis System)是美國匹茲堡大學(xué)John Swanson教授在1970年開發(fā)。經(jīng)過大大小小幾十次的更新?lián)Q代，現(xiàn)在的ANSYS軟件集眾多分析模塊于一身。在機械、土木、航空航天、能源、交通等各個領(lǐng)域都有應(yīng)用，逐步成為世界上最有影響力的有限元軟件之一[9]。

6.3 優(yōu)化分析模塊

6.3.1 優(yōu)化設(shè)計

所謂“優(yōu)化設(shè)計”就是在很多套設(shè)計方案中選出最優(yōu)設(shè)計的一種方法。而最優(yōu)設(shè)計方案指的是在滿足所有的設(shè)計需求前提下，花費人力、物力、資金等最小。對于結(jié)構(gòu)設(shè)計大賽而言，最優(yōu)設(shè)計就是在滿足結(jié)構(gòu)強度、變形、穩(wěn)定性等基本要求下所需的材料質(zhì)量最小，在材料密度相同時就是材料體積最小。

6.3.2 優(yōu)化變量的確定

在ANSYS優(yōu)化分析模塊會涉及幾種優(yōu)化變量，包括設(shè)計變量、狀態(tài)變量以及目標(biāo)函數(shù)等[10]。

設(shè)計變量，即自變量，也就是我們想優(yōu)化的量。對于本文來講就是桿件的截面尺寸，在結(jié)構(gòu)體系和構(gòu)件的截面形狀確定之后，就是使構(gòu)件的截面尺寸盡可能達(dá)到最優(yōu)。

狀態(tài)變量，隨設(shè)計變量的改變而改變，即因變量，也就是設(shè)計需要滿足的條件。對于本文來講就是，第一級加載豎向位移UY≤10 mm，加載的任意階段桿件壓應(yīng)力不超過竹材順紋抗壓強度30 MPa，拉應(yīng)力不超過竹材順紋抗拉強度60 MPa。

目標(biāo)函數(shù)，是設(shè)計變量(自變量)的函數(shù)，是想盡量減少的量。對于本文來講就是結(jié)構(gòu)材料質(zhì)量盡量最小。因結(jié)構(gòu)使用同一種材料，密度相同，也就是材料體積盡量小。

各桿件初始設(shè)計值如表1所示。

表1 桿件截面參數(shù)

6.3.3 ANSYS模型建立

1)單元類型選擇。采用Structural Mass/Beam 2 node 188單元進(jìn)行分析，該單元可以較好的模擬實際情形[11]。

2)定義桿件初始設(shè)計變量。A=10 mm，B=15 mm，C=15 mm，D=10 mm，E=2 mm。

3)定義材料屬性。材料彈性模量6 GPa，泊松比0.32，密度為0.8 g/cm3。

4)ANSYS模型建立以及網(wǎng)格劃分見圖4。

5)施加邊界條件。位移約束：1柱、2柱端部共6個點施加固定端約束。力約束：如前所述施加8個加載點的豎向集中荷載(取重力加速度g=10 N/kg)。圖5為模型邊界條件施加圖。

6)初始設(shè)計變量下的計算結(jié)果。進(jìn)行求解，部分結(jié)果如表2所示。

表2 初始設(shè)計變量下的結(jié)果

由表2可知，最大豎向節(jié)點位移DY=16.023 mm，大于賽題要求10 mm。即初始設(shè)計變量不滿足設(shè)計要求。2柱的壓應(yīng)力最大，為10.015 MPa，小于賽題規(guī)定順紋抗壓強度30 MPa；拉條拉應(yīng)力最大為31.622 MPa小于賽題規(guī)定60 MPa，說明結(jié)構(gòu)強度滿足設(shè)計要求。材料用量為237 757 mm3。乘以材料密度0.8 g/mm3，理論上結(jié)構(gòu)重量190.205 6 g，實際大約重量乘以系數(shù)1.2(制作膠水等會加大模型重量，每個人系數(shù)不同，僅供參考)，為228.246 72 g。

圖6為每根桿件的體積(質(zhì)量)分布圖。

由圖6可知，模型的質(zhì)量主要集中在2柱，其次是梁，再者是1柱，最后是拉條。

7)優(yōu)化計算。定義設(shè)計變量控制范圍為：A：5 mm～15 mm，B:10 mm～20 mm，C:10 mm～20 mm，D:5 mm～15 mm，E:1 mm～3 mm。

定義狀態(tài)變量：DY=ABS(UY)<10 mm，即控制最大豎向位移不超過10 mm；桿件最大壓應(yīng)力小于30 MPa，最大拉應(yīng)力小于60 MPa，即滿足材料順紋抗壓強度不超過30 MPa，順紋抗拉強度不超過60 MPa。

優(yōu)化方法采用子問題近似法，迭代次數(shù)為50。

進(jìn)行迭代運算，得到最優(yōu)結(jié)果如下：A=14.090 mm，B=16.145 mm，C=19.977 mm，D=14.980 mm，E=2.987 9 mm。

最優(yōu)設(shè)計變量下的結(jié)果如表3所示。

表3 最優(yōu)設(shè)計變量下的結(jié)果

在最優(yōu)設(shè)計變量下，結(jié)構(gòu)最大變形為DY=10.102 mm(計算收斂于10 mm)，小于初始設(shè)計16.023 mm，滿足賽題規(guī)定10 mm,證明ANSYS優(yōu)化計算效果明顯。模型桿件最大拉應(yīng)力為30.867 MPa，小于初始設(shè)計31.622 MPa，小于材料順紋抗拉強度60 MPa。2柱壓應(yīng)力最大，為6.785 MPa，小于初始設(shè)計10.015 MPa，小于材料順紋抗壓強度30 MPa。模型體積一共為284 290.293 mm3，乘以材料密度0.8 g/mm3，理論上結(jié)構(gòu)重量227.992 g。實際大約重量乘以系數(shù)1.2(制作膠水等會加大模型重量，每個人系數(shù)不同，僅供參考)，為273.590 g。

最優(yōu)設(shè)計變量下材料體積分布如圖7所示，與初始設(shè)計相比，梁的質(zhì)量分布明顯提高，這也是豎向位移明顯降低的原因。與實際情況相符，也體現(xiàn)了ANSYS軟件在結(jié)構(gòu)設(shè)計大賽中應(yīng)用的可行性。

圖8為優(yōu)化前、優(yōu)化后和賽題要求數(shù)據(jù)對比圖。

從圖8中可以看出無論是材料的壓應(yīng)力還是拉應(yīng)力，其數(shù)值優(yōu)化后均小于優(yōu)化前，并小于賽題規(guī)定數(shù)值。而結(jié)構(gòu)最大變形更是從優(yōu)化前的不符合賽題要求變?yōu)闈M足賽題要求。

圖9為生成的設(shè)計變量對迭代次數(shù)的敏感性變化圖。

從圖9中可以看出參數(shù)E對迭代次數(shù)最不敏感，參數(shù)B，C對迭代次數(shù)最敏感。即同等條件下梁高度(參數(shù)C)和2柱寬度(參數(shù)B)比改變拉條尺寸(參數(shù)E)對控制模型變形，減少桿件應(yīng)力效果更好。提高了模型設(shè)計效率，與實際情況吻合。

7 模型制作與加載

實際模型(如圖10所示)制作需要把制作誤差、材料強度損失等其他不確定因素考慮在內(nèi)，因此制作時給予模型一定安全度。加載情況(如圖11所示)與模擬情況類似。

8 結(jié)論

通過ANSYS的優(yōu)化計算與實際模型加載情況進(jìn)行比較得到以下結(jié)論：

1)ANSYS的優(yōu)化計算可以幫助參賽選手更快的設(shè)計出耗材少又滿足賽題要求的結(jié)構(gòu)，具有一定的參考價值，給競賽的發(fā)展帶來了新思路。

2)ANSYS的應(yīng)用減少了模型制作加載次數(shù)，提高了效率。

3)該軟件的計算結(jié)果可以指導(dǎo)實際模型的制作，方便參賽選手理解加載過程以及結(jié)構(gòu)整體受力情況，以便進(jìn)一步進(jìn)行優(yōu)化設(shè)計。