韓家山,曹翁愷,顧海龍,陳新培
(洛陽雙瑞特種裝備有限公司,河南 洛陽 471000)
作為連接橋梁上下部結(jié)構(gòu)的“關(guān)節(jié)”,橋梁支座可以將橋梁上部結(jié)構(gòu)中反力和變形可靠地傳遞給橋梁下部結(jié)構(gòu),是橋梁結(jié)構(gòu)中的一個(gè)重要組成部分[1]。隨著我國鐵路、公路交通運(yùn)輸?shù)陌l(fā)展,大跨度橋梁越來越多,橋梁支座的承載力也在不斷刷新紀(jì)錄。目前,支座最大承載力已達(dá)2×106kN,支座質(zhì)量已高達(dá)百噸。橋梁支座過重不僅會(huì)增加建設(shè)成本,而且會(huì)增加現(xiàn)場施工的難度。因此,在保證使用性能的前提下,對(duì)橋梁支座進(jìn)行輕量化設(shè)計(jì),對(duì)于減少建設(shè)成本和降低施工難度具有重要的實(shí)際意義。目前,針對(duì)橋梁支座輕量化的研究,大多是從支座結(jié)構(gòu)設(shè)計(jì)參數(shù)的角度出發(fā),予以形狀優(yōu)化和尺寸優(yōu)化,而有關(guān)支座最優(yōu)拓?fù)湫问降奶接懸娭谖墨I(xiàn)的研究還很少[2-4]。
拓?fù)鋬?yōu)化能在給定的空間結(jié)構(gòu)中得到優(yōu)化的形狀及材料分布形式,通過在結(jié)構(gòu)區(qū)域內(nèi)形成一些孔洞而去除結(jié)構(gòu)中多余的材料,以達(dá)到結(jié)構(gòu)最優(yōu)。與形狀優(yōu)化和尺寸優(yōu)化相比,拓?fù)鋬?yōu)化能夠在概念設(shè)計(jì)階段就提供一種合理而又富有啟發(fā)性的結(jié)構(gòu)形式[5-7],在一些零部件的輕量化設(shè)計(jì)中已得到廣泛應(yīng)用。本文以結(jié)構(gòu)拓?fù)鋬?yōu)化方法為理論基礎(chǔ),基于HyperWorks中Optistruct優(yōu)化軟件,以一座大橋豎向承載力是 17 500 kN 的多向球型支座下座板為研究對(duì)象,以剛度最大化為設(shè)計(jì)目標(biāo),對(duì)支座下座板進(jìn)行拓?fù)鋬?yōu)化設(shè)計(jì),旨在為橋梁支座輕量化設(shè)計(jì)提供一種新的設(shè)計(jì)思路。
根據(jù)支座初始設(shè)計(jì)尺寸建立支座下座板三維模型,并在HyperMesh中對(duì)模型進(jìn)行前處理。采用四面體進(jìn)行網(wǎng)格劃分,建立的模型共包含 120 863 個(gè)單元和 27 539 個(gè)節(jié)點(diǎn)。材料彈性模量為 210 GPa,密度為7.8×103kg/m3,泊松比為0.3。根據(jù)支座工作原理可知,支座下座板主要承受來自中座板的豎向荷載,荷載通過中座板球面?zhèn)鬟f給支座下座板,然后再傳遞給墩臺(tái)。根據(jù)支座實(shí)際工況,對(duì)下座板模型施加邊界條件,即對(duì)下座板球面施加F=17 500 kN的集中力。為保證集中力均勻分布,將球面上節(jié)點(diǎn)進(jìn)行多點(diǎn)耦合,下座板底面進(jìn)行全約束。荷載及約束模型見圖1。
圖1 荷載及約束模型
對(duì)初始設(shè)計(jì)模型進(jìn)行分析,計(jì)算其變形及應(yīng)力,下座板變形及應(yīng)力云圖如圖2所示??芍?,支座下座板最大變形為0.05 mm,最大變形發(fā)生在下座板的球面區(qū)域,由上到下變形逐漸變小。最大應(yīng)力為38.42 MPa,最大應(yīng)力出現(xiàn)在球面中心,由中心向周邊逐漸減小。另外,在外圓柱面根部也有較大應(yīng)力區(qū)域。
圖2 下座板變形及應(yīng)力云圖
拓?fù)鋬?yōu)化作為一種優(yōu)化設(shè)計(jì)方法,同樣具有優(yōu)化設(shè)計(jì)3要素:目標(biāo)函數(shù)、設(shè)計(jì)變量及約束條件。本文以支座下座板體積為約束條件,以支座下座板剛度最大化(最小柔順度)為設(shè)計(jì)目標(biāo),對(duì)支座下座板進(jìn)行輕量化設(shè)計(jì)。在體積約束下的結(jié)構(gòu)最大化剛度問題[8]可以描述為
Findx={x1,x2,…,xi},i=1,2,…,N
(1)
式中:x為設(shè)計(jì)變量;xi為單元設(shè)計(jì)變量,包括1和一個(gè)很小的值xmin;C(x)為單元結(jié)構(gòu)柔順度;U為位移矩陣;K為總剛度矩陣;ui為單元位移矩陣;ki為單元?jiǎng)偠染仃?;N為結(jié)構(gòu)單元數(shù)量;p為材料懲罰因子;V*為初始結(jié)構(gòu)總體積;Vi為單元體積。
Optistruct軟件中采用的是SIMP材料插值模型。在SIMP方法中對(duì)于設(shè)計(jì)域內(nèi)的單元,其彈性模量可引入如下關(guān)系式[8]
(2)
式中,E0為材料的真實(shí)楊氏模量。
對(duì)于不需要去除的單元,由于設(shè)計(jì)變量xi=1,因此其楊氏模量就是材料的真實(shí)楊氏模量。對(duì)于需要去除的單元,設(shè)計(jì)變量xi=xmin,再加上材料懲罰因子指數(shù),此時(shí)該單元的楊氏模量與真實(shí)楊氏模量相比差了幾個(gè)數(shù)量級(jí)(取決于xmin和p值),是一個(gè)很小的數(shù),這樣的單元一般被稱為軟單元或空單元。最后通過推導(dǎo)單元靈敏度公式得到單元的靈敏度值[8]為
(3)
根據(jù)單元靈敏度值對(duì)單元進(jìn)行增刪操作,經(jīng)過一定次數(shù)的迭代和循環(huán)便可得到最優(yōu)的拓?fù)湫问健?/p>
考慮支座最初拓?fù)湫问轿粗?,首先要根?jù)支座實(shí)際工況確定支座邊界條件。支座下座板通過地腳螺栓與墩臺(tái)相連來傳遞梁部的水平荷載。因此,支座下座板地腳螺栓孔周圍區(qū)域必須進(jìn)行約束,而下座板底面中間區(qū)域?yàn)椴淮_定設(shè)計(jì)區(qū),可不進(jìn)行約束。因此可以將底面進(jìn)行全約束或部分約束,在2種不同的邊界條件下分別進(jìn)行優(yōu)化設(shè)計(jì)。底面模型見圖3。考慮到球型支座需具備傳力、滑移等正常功能,將支座下座板分為設(shè)計(jì)區(qū)域與非設(shè)計(jì)區(qū)域,分別對(duì)應(yīng)圖中綠色部分和黃色部分。非設(shè)計(jì)區(qū)域在優(yōu)化過程中,只參與結(jié)構(gòu)的有限元分析,但不進(jìn)行單元的更新操作。
圖3 底面模型
在Optistruct軟件中,分別定義支座下座板的目標(biāo)函數(shù)、約束條件、設(shè)計(jì)變量等參數(shù),然后進(jìn)行拓?fù)鋬?yōu)化計(jì)算。支座下座板在進(jìn)行拓?fù)鋬?yōu)化后保留了密度值在0.3以上的單元后的最終結(jié)構(gòu),底面非全約束模型密度圖視角1和視角2分別見圖4(a)和圖4(b)。
圖4 底面非全約束模型密度圖
理論上,模型中應(yīng)該只存在2種單元,即密度值為0的單元和密度值為1的單元。密度值為0的單元代表空單元,對(duì)結(jié)構(gòu)貢獻(xiàn)度最小,可以刪除。密度值為1的單元代表實(shí)單元,對(duì)結(jié)構(gòu)貢獻(xiàn)度最大,應(yīng)該保留。實(shí)際上,即使設(shè)定了材料懲罰因子對(duì)這些單元進(jìn)行懲罰,結(jié)構(gòu)中依然會(huì)存在較多中間密度單元,仍不能使它們的密度值向0或向1逼近。根據(jù)經(jīng)驗(yàn),一般將密度值在0.3以上的單元進(jìn)行保留,低于0.3的單元進(jìn)行刪除。底面全約束模型密度圖見圖5。
圖5 底面全約束模型密度圖
對(duì)比2種模型密度圖可知,底面非全約束模型在結(jié)構(gòu)優(yōu)化過程中,支座下座板中間部分材料對(duì)結(jié)構(gòu)貢獻(xiàn)度較低,通過逐步刪除該部分材料形成了一種底板中空的結(jié)構(gòu)形式,達(dá)到了結(jié)構(gòu)的剛度最優(yōu)。底面全約束模型通過逐漸刪除外圓材料,以外圓中縮形式達(dá)到了支座下座板的剛度最優(yōu)。
經(jīng)過拓?fù)鋬?yōu)化形成的密度圖形式相對(duì)復(fù)雜,不能直接用于實(shí)際生產(chǎn),還必須進(jìn)行二次設(shè)計(jì)??紤]中空設(shè)計(jì)會(huì)造成比較大的墩臺(tái)應(yīng)力,因此只對(duì)底面全約束模型的拓?fù)浣Y(jié)果進(jìn)行二次設(shè)計(jì)。
根據(jù)拓?fù)鋬?yōu)化結(jié)果密度圖并考慮可鑄造性能及可加工性,對(duì)支座下座板進(jìn)行二次設(shè)計(jì)。設(shè)計(jì)中,主要對(duì)下座板外圓直徑進(jìn)行了一定程度的減小,二次設(shè)計(jì)后支座下座板三維模型見圖6。
圖6 二次設(shè)計(jì)后支座下座板三維模型
對(duì)優(yōu)化后模型再次進(jìn)行有限元分析,以驗(yàn)證模型是否滿足強(qiáng)度和剛度要求。優(yōu)化后模型的變形及應(yīng)力云圖見圖7。可知,拓?fù)鋬?yōu)化后支座下座板最大變形(0.15 mm)發(fā)生在下座板上部圓周區(qū)域;最大應(yīng)力(82.22 MPa)發(fā)生在下座板外圓的中間部位。
圖7 優(yōu)化后模型的變形及應(yīng)力云圖
支座下座板優(yōu)化前后參數(shù)對(duì)比如表1所示??芍ё伦鍍?yōu)化后質(zhì)量減小了21.2%,最大位移由0.05 mm增加到0.15 mm,最大應(yīng)力由38.42 MPa增加到82.22 MPa。根據(jù)GB 50017—2003《鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范》[9]可知,下座板最大變形不超過0.5 mm,容許應(yīng)力為150 MPa,優(yōu)化后結(jié)構(gòu)強(qiáng)度和剛度滿足要求。通過拓?fù)鋬?yōu)化,在減小結(jié)構(gòu)剛度與強(qiáng)度的情況下達(dá)到了輕量化的目的。
表1 支座下座板優(yōu)化前后參數(shù)對(duì)比
本文以橋梁支座為研究對(duì)象,利用結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)思想,在2種不同的邊界條件下對(duì)下座板進(jìn)行了輕量化設(shè)計(jì),得到以下結(jié)論:
1)對(duì)于底面非全約束模型,結(jié)構(gòu)在優(yōu)化過程中,通過分析單元靈敏度判定支座下座板中間部分材料對(duì)結(jié)構(gòu)的貢獻(xiàn)度較低,通過逐步刪除該部分材料,得到了一種底板中空的結(jié)構(gòu)形式。該結(jié)構(gòu)雖然不能滿足墩臺(tái)應(yīng)力要求,但是相對(duì)于初始設(shè)計(jì)模型,其材料消耗大大降低。同時(shí),該結(jié)構(gòu)形式給設(shè)計(jì)者預(yù)留了較大的設(shè)計(jì)空間,為支座實(shí)現(xiàn)更多功能(如測力、隔振)提供了更多可能。
2)對(duì)于底面全約束模型,結(jié)構(gòu)在優(yōu)化過程中,通過逐漸刪除外圓部分材料,支座下座板形成了一種外圓中縮的結(jié)構(gòu)形式。
3)通過對(duì)外圓中縮結(jié)構(gòu)二次構(gòu)型及再分析,底面全約束模型拓?fù)鋬?yōu)化得到的支座下座板結(jié)構(gòu)重量減小了21.2%,驗(yàn)證了拓?fù)鋬?yōu)化理論在橋梁支座輕量化中應(yīng)用的可行性,為其他橋梁支座結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了思路和依據(jù)。
[1]莊軍生.橋梁支座[M].北京:中國鐵道出版社,2015.
[2]楊春平,胡宇新,郭紅峰,等.GPZ(Ⅱ)盆式橡膠支座結(jié)構(gòu)優(yōu)化與試驗(yàn)研究[J].鐵道建筑,2014,54(3):22-23.
[3]張常勇,鐘鐵毅,顧正偉,等.鐵路連續(xù)梁橋鉛芯橡膠支座優(yōu)化設(shè)計(jì)研究[J].鐵道學(xué)報(bào),2012,34(5):97-100.
[4]劉岳兵.鐵路橋梁支座數(shù)值仿真研究[D].成都:西南交通大學(xué),2010.
[5]左孔天.連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化理論與應(yīng)用研究[D].武漢:華中科技大學(xué),2004.
[6]郭中澤,張衛(wèi)紅,陳裕澤.結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)綜述[J].機(jī)械設(shè)計(jì),2007,24(8):1-6.
[7]彭細(xì)榮,隋允康.對(duì)連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化合理模型的再探討[J].固體力學(xué)學(xué)報(bào),2016,37(8):181-184.
[8]BENDSOE M P,SIGMUND O.Topology Optimization:Theory,Method and Applications[M].New York:Springer Berlin Heidelberg,2003.
[9]中華人民共和國建設(shè)部,中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 50017—2003鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范[S].北京:中國計(jì)劃出版社,2003.