冉孟廷 周婭

摘 要：為解決大型橋梁工程中細(xì)部構(gòu)件力學(xué)分析時(shí)，有限元建模工作量與計(jì)算精度相矛盾的問(wèn)題，以某長(zhǎng)江大橋?yàn)槔捎么笮陀邢拊?jì)算軟件Abaqus，建立多尺度有限元模型進(jìn)行鋼混結(jié)合段受力分析，并與局部精細(xì)有限元模型和全橋精細(xì)有限元模型計(jì)算結(jié)果作對(duì)比，結(jié)果表明：多尺度有限元模型計(jì)算精度較高，能夠滿(mǎn)足工程應(yīng)用要求，并能有效降低建模工作量，提高運(yùn)算速度。

關(guān)鍵詞：軟件技術(shù)；多尺度有限元模型；建模；計(jì)算精度；橋梁工程

0 前言

工程應(yīng)用中，傳統(tǒng)力學(xué)分析以靜力學(xué)為主，并以力學(xué)方程的精確解析解為基礎(chǔ)。但隨著計(jì)算對(duì)象在結(jié)構(gòu)與材料上的復(fù)雜化，特別是在涉及到材料的粘彈塑性時(shí)，不僅計(jì)算量急劇增大，力學(xué)方程也更難以建立和求解，從而不得不簡(jiǎn)化受力條件或作出近似假設(shè)，但其求解結(jié)果往往與實(shí)際受力有所偏差。為解決復(fù)雜構(gòu)件和材料粘彈塑性的力學(xué)計(jì)算問(wèn)題，采用離散化分析從而求取近似數(shù)值解的大型有限元軟件應(yīng)運(yùn)而生，目前工程中常用的有Abaqus、Adina、Ansys和Marc等。這些大型有限元軟件針對(duì)固體、流體、氣體、磁場(chǎng)、熱力場(chǎng)以及耦合計(jì)算對(duì)象，在結(jié)構(gòu)分析、線性和非線性分析等領(lǐng)域各有所長(zhǎng)。

但對(duì)于復(fù)雜結(jié)構(gòu)的大型橋梁，采用有限元分析軟件進(jìn)行細(xì)部結(jié)構(gòu)力學(xué)分析時(shí)，也面臨建模工作量和計(jì)算精度的矛盾：若離散單元?jiǎng)澐诌^(guò)密，不但建模工作量過(guò)大，計(jì)算速度也會(huì)大大降低；若單元?jiǎng)澐诌^(guò)少，則計(jì)算精度又會(huì)大大降低，且不利于考察細(xì)部構(gòu)件的受力特征。為解決這一問(wèn)題，工程中常常采用簡(jiǎn)化邊界約束條件的局部精細(xì)有限元模型進(jìn)行分析，但其計(jì)算結(jié)果有可能與實(shí)際受力有所差異，但如果細(xì)部結(jié)構(gòu)單元進(jìn)行精細(xì)化處理，而橋梁其它部分采用稀疏網(wǎng)格進(jìn)行劃分，這樣的多尺度有限元模型不僅兼顧了細(xì)部構(gòu)件的計(jì)算精度問(wèn)題，也大大降低了建模工作量，提高計(jì)算速度。

本文以某長(zhǎng)江大橋鋼混結(jié)合段為例，采用大型有限元分析軟件Abaqus，建立局部精細(xì)有限元模型、全橋精細(xì)有限元模型和多尺度有限元模型，進(jìn)行計(jì)算精度分析，介紹多尺度有限元分析模型的應(yīng)用過(guò)程。

1工程概況

某長(zhǎng)江大橋?yàn)殡p向六車(chē)道高速公路特大橋，主跨為半漂浮體系的7跨連續(xù)雙塔混合梁斜拉橋，橋跨布置為64+2×68+608+2×68+64=1008m。邊跨采用混凝土箱梁，中跨采用鋼箱梁，鋼箱梁段索距為15.5m，混凝土梁段索距為10m和8m，塔上豎向標(biāo)準(zhǔn)索距為2.5m。索塔采用C50混凝土，塔高分別為196.7m和206.4m，索塔基礎(chǔ)采用24根直徑2.5m的鉆孔灌注樁。大橋整體布置如圖1所示，鋼混結(jié)合段細(xì)部尺寸如圖2所示。

2 有限元建模與單元?jiǎng)澐?/p>

2.1 局部精細(xì)有限元模型

根據(jù)該橋設(shè)計(jì)參數(shù)，建立局部精細(xì)有限元模型，其中鋼結(jié)構(gòu)采用S4R殼單元模擬，混凝土結(jié)構(gòu)采用實(shí)體單元C3D8R模擬，斜拉索采用桁架單元T3D2模擬。鋼結(jié)構(gòu)部分網(wǎng)格劃分較密，尺寸為200×200mm；與鋼結(jié)構(gòu)連接的混凝土結(jié)構(gòu)部分網(wǎng)格劃分尺寸保持與鋼結(jié)構(gòu)一致，為200×200×200mm，其他部分為600×600×600mm。模型共計(jì)包括55835個(gè)殼單元，14690個(gè)實(shí)體單元以及2個(gè)桁架單元。

邊界條件為：約束混凝土結(jié)構(gòu)部分截?cái)嘟缑娴娜?個(gè)方向自由度，在鋼結(jié)構(gòu)部分的截?cái)嘟缑嫣帉?duì)模型施加作用力；約束斜拉索頂部的平動(dòng)自由度，并采用降溫法對(duì)斜拉索施加張拉力。

2.2 全橋精細(xì)有限元模型

采用的單元類(lèi)型和局部精細(xì)有限元模型相同。為控制單元數(shù)量，結(jié)合段局部單元尺寸為200×200mm，其他部分為800×800mm；與鋼結(jié)構(gòu)連接的混凝土結(jié)構(gòu)部分網(wǎng)格劃分尺寸為200×200×200mm，其他部分為800×800×800mm。全橋精細(xì)有限元模型共計(jì)包括259130個(gè)殼單元，34160個(gè)實(shí)體單元及76個(gè)桁架單元。

邊界條件為：采用大橋的實(shí)際位移約束條件，約束斜拉索頂部全部平動(dòng)自由度，施加主梁在橋塔、輔助墩處的豎向及側(cè)向位移約束。斜拉索底端與主梁對(duì)應(yīng)吊點(diǎn)處通過(guò)Coupling連接，其張拉力通過(guò)降溫法施加，各索初始張拉力數(shù)值通過(guò)整橋桿系MIDAS civil模型獲得。

2.3 多尺度有限元模型

包括精細(xì)模型部分以及桿系模型部分，后者參考整橋桿系MIDAS civil模型建立，主要包括斜拉索與主梁。因?yàn)槟Ｐ头治鲋胁豢紤]橋塔、橋墩沉降產(chǎn)生的影響，所以幾何模型中不直接建立橋塔和橋墩，而將其作為邊界條件分別施加給斜拉索與主梁。

所采用單元與局部精細(xì)有限元模型和全橋精細(xì)有限元模型相同，但采用梁?jiǎn)卧狟31模擬主梁。精細(xì)模型部分網(wǎng)格尺寸與局部精細(xì)模型完全一致；桿系模型部分中，主梁的劃分尺寸為600mm，由于斜拉索部分傳力明確且不考慮拉索的幾何非線性，故使用單個(gè)桁架單元模擬一根斜拉索。多尺度有限元模型共計(jì)包括55835個(gè)殼單元，14690個(gè)實(shí)體單元，794個(gè)梁?jiǎn)卧?6個(gè)桁架單元。

多尺度有限元模型邊界條件與全橋精細(xì)有限元模型相同。使用降溫法模擬各根斜拉索的初張拉力，施加自重荷載作用。

3 計(jì)算參數(shù)與計(jì)算結(jié)果

計(jì)算時(shí)，各有限元模型中主要材料參數(shù)相同，即：

鋼材：彈性模量Es=2.06×105MPa，泊松比υ=0.3；

混凝土：彈性模量Es=3.55×104MPa，泊松比υ=0.2；

斜拉索：彈性模量Ep=1.95×105MPa。

鋼混結(jié)合段的受力計(jì)算結(jié)果如圖9 ～ 圖11所示。

鋼混結(jié)合段應(yīng)力計(jì)算結(jié)果表明，多尺度有限元模型在斜拉索吊點(diǎn)位置有明顯的應(yīng)力集中現(xiàn)象，全橋精細(xì)模型中也出現(xiàn)類(lèi)似的分布現(xiàn)象，但在局部精細(xì)模型中并未出現(xiàn)；另外，局部精細(xì)有限元模型在截?cái)嘟缑娴玫降膽?yīng)力水平相對(duì)較低，低于多尺度有限元模型計(jì)算結(jié)果，而全橋精細(xì)有限元模型的分析結(jié)果與多尺度有限元模型基本一致，從而驗(yàn)證了多尺度有限元模型分析結(jié)果的合理性。

4結(jié)語(yǔ)

本文以某長(zhǎng)江大橋鋼混結(jié)合段為例，采用大型有限元計(jì)算軟件Abaqus，通過(guò)建立局部精細(xì)有限元模型、全橋精細(xì)有限元模型和多尺度有限元模型，研究力學(xué)計(jì)算精度的差別和多尺度有限元模型的可行性，成功地解決了工程應(yīng)用技術(shù)中的棘手問(wèn)題，得出以下結(jié)論：

（1）多尺度有限元模型計(jì)算結(jié)果和精細(xì)有限元模型相近，其精度能夠滿(mǎn)足工程應(yīng)用要求；

（2）多尺度有限元模型能夠有效降低建模工作量，提高運(yùn)算速度；

（3）受邊界條件的影響，局部精細(xì)有限元模型計(jì)算精度遠(yuǎn)低于多尺度有限元模型；

（4）合理選擇現(xiàn)代信息化處理技術(shù)，能夠?yàn)楣こ碳夹g(shù)人員解決工程技術(shù)中的棘手問(wèn)題帶來(lái)事半功倍的效果。

參考文獻(xiàn)

[1]龔謙.通用有限元分析軟件ABAQUS在橋梁分析中的應(yīng)用[J].城市建設(shè)理論研究：電子版，2012（21）

[2]許麗嬌，葉平，葉惺.淺談大型有限元分析軟件Ansys[J].煤礦機(jī)械，2005（4）：49～51.

[3]趙培仲，胡芳友，黃旭仁等.有限元分析軟件在結(jié)構(gòu)損傷修理課中的應(yīng)用[J].中國(guó)信息技術(shù)教育，2015（11）：122～123.

[4]龍軍，張濤，陳居術(shù)等.復(fù)雜組合結(jié)構(gòu)溫度場(chǎng)及應(yīng)力有限元分析[J].機(jī)械設(shè)計(jì)與制造工程，2002，31（5）：55～57.

[5]王志義，錢(qián)良宏.淺談通用分析軟件對(duì)某橋的結(jié)構(gòu)分析[J].四川建材，2007，33（4）