鋼-混結合梁有限元分析技術研究及應用

盧 飛 ,關 薇

(中國電建集團西北勘測設計研究院有限公司，西安 710065)

文章編號：1006—2610(2015)04—0068—04

鋼-混結合梁有限元分析技術研究及應用

盧 飛 ,關 薇

(中國電建集團西北勘測設計研究院有限公司，西安 710065)

漢中市西二環三塔斜拉-懸索組合橋梁，結構空間構造錯綜復雜，設計難度大。文章利用ANSYS邊界耦合方式建立一種符合實際情況的結合梁有限元分析模型，研究大跨度橋梁工程中結合梁段在柔索組合體系下鋼箱梁的受力特性；利用ANSYS的參數化設計語言編制了分析文件和優化控制文件，分析解決可能存在局部應力過大和局部穩定等關鍵性問題，經計算獲得結合梁的最優截面形式，從而為復雜的空間結構設計提供有益的參考。

鋼-混結合梁;ANSYS;參數化;約束方程

0 前 言

漢中市西二環大橋為三塔斜拉-自錨式懸索組合體系橋，大橋主橋上部結構主要由2個對稱的副塔和1個主塔組成的邊跨自錨段、懸索區段和斜拉索區段3部分組成，如圖1所示主橋跨徑組合為：邊跨自錨段(25 m)+懸索段(90 m)+斜拉索區段(2×162.5 m)+懸索段(90 m)+邊跨自錨段(25 m)，橋梁主橋全長555 m。

結合梁是懸索段承受主纜吊桿拉力主要結構，全橋設置左右對稱2段結合梁，每段長63 m，梁高2.5 m，梁頂寬40 m，梁底寬29 m，懸臂端4.45 m，見圖2。

圖1 大橋主橋立面圖 單位：cm

結合梁采用單箱4室結構，結合梁內鋼錨箱與吊桿相對應，每個節段左、右半幅邊箱室中各設置1榀縱橋向鋼錨梁，鋼錨箱固定于鋼錨梁之中，且鋼錨箱角度與吊桿空間角度相同；鋼錨箱采用井格狀布置，豎向40 mm厚鋼板與其垂直50 mm厚底板組成，鋼錨梁單側為2塊豎向32 mm厚槽形鋼板，采用豎向聯系鋼板連接，左右兩單側再與鋼錨箱垂直連接，頂、底采用鋼板連接；鋼錨梁兩端垂直焊接在結合梁鋼箱橫隔板上,鋼錨梁、鋼錨箱連接詳見圖3、4。

鋼箱結合梁由于承受主纜吊桿拉力，同時傳遞與邊跨自錨段的剪力，結構受力復雜。常規的結構計算軟件無法準確模擬其受力狀態。

運用ANSYS既可合理模擬斜拉-懸索組合橋結合梁空間受力狀態，又可將其成果直接應用其它復雜斜接-懸索橋梁工程，具有很高的實用價值。

圖2 結合梁段立面布置圖 單位：mm

圖3 鋼錨梁、鋼錨箱連接透視圖

圖4 鋼錨梁、鋼錨箱橫橋向剖面構造圖 單位：mm

1 鋼箱結合梁參數化建模枝術

西二環結合梁段內布置主纜吊桿鋼錨梁、鋼錨箱，結合梁段在荷載作用下，吊桿力通過鋼錨箱傳遞到鋼錨梁，再由鋼錨梁傳遞到結合梁橫隔板上，受力狀態復雜；結合梁自身又布置錯綜復雜的腹板、橫隔板、小縱梁、加勁板等，空間結構復雜。設計階段為解決結合梁段復雜受力狀態、空間結構，需要不斷地調整模型的幾何結構、載荷、邊界條件等，從而對結合梁段進行結構優化，實現結合梁段整體受力及空間結構的最佳。這個調試過程可能很不經濟。利用ANSYS的APDL語言采用參數化建模，不僅提高了模型建模的效率，而且便于模型結構修改、設計優化。

1.1 APDL的參數定義和流程控制

APDL的一個重要特征就是可實現參數化建模，參數是指變量和數組，變量參數有2種類型：數值型和字符型；數組參數有3種：數值型、字符型和表[1]。

利用APDL建模時變量的定義和賦值有6種途徑：① 利用*SET命令；② 利用賦值符號“=”；③ 利用菜單路徑Utility Menu>Parametres>Scalar Parametres；④ 在啟動時利用驅動命令；⑤ 利用*GET命令及其等效函數提取ANSYS數據庫數據；⑥ 利用*ASK命令。

APDL語言中的流程控制語句主要有：① *GO無條件分支；② *IF-*ELSEIF-*ELSE-*ENDIF條件分支；③ *DO-*ENDDO循環；④ *DOWHILE循環；⑤ *REPEAT重復一個命令，增加一個或多個命令參數。其中*DO-*ENDDO語句最為常用，它類似于C語言的for循環語句，在建模過程中遇到的多次重復操作可采用此語句，如重復移動工作平面、重復建立線面等。

本文采用PSPad作為APDL開發工具，它能提供語法高亮顯示、自動糾錯、拼寫檢查等功能，十分方便。

1.2 訪問數據庫、編號控制技術

在ANSYS中，每一個元素在創建時都會被自動分配一個編號，每一種類型的元素如關鍵點、線、面、節點和單元等都是單獨編號的，不同類型的編號互不相關。因此，當需要對命令流改寫時，點、線、面、體等元素的編號就會發生變化，特別當模型較為復雜時，建模過程非常容易出錯[2]。

利用APDL提供訪問ANSYS數據庫中各種數據的命令，如環境數據、目錄路徑、當前工作名、模型數據、結果數據等。既可控制元素的編號，又能避免元素編號變化時建模出錯的問題。

ANSYS中，*GET命令的功能十分強大，幾乎可以提取數據庫中的任何數據，包括任何對象(點、線、面、節點、單元等)的相關數據信息以及處理器的設置或狀態數據信息等。

*GET使用格式：*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM。其中，Par是賦值的參數名；Entity是被提取對象關鍵字；可以是NODE、ELEM、KP、LINE、AREA和VOLU等；ENTNUM是實體編號(為0時指全體實體)；Item1是某實體對應的項目名，在ANSYS Commands Reference中有完整說明。

實例：

(1) *GET,NNUM,KP,,COUNT !當前激活的關鍵點數目賦給NNUM;

(2) *GET,NMAX,NODE,,NUM,MAX !當前激活的節點較大的編號賦給NMAX;

(3) *GET,vtot,SSUM,ITEM,VOLUME !提取總體積并存貯到變量vtot中。

多使用*GET等訪問數據庫的命令，可使APDL的編寫更加簡潔。

在采用APDL進行參數化編程時，通過NUMCMP,KP 命令可壓縮所有關鍵點編號，并結合*GET命令獲取關鍵點數目，使得創建線段的命令代碼中不需要關鍵點編號。在建立其他實體模型中均可利用類似方法實現。

1.3 建立參數化模型

建模前首先對結構進行合理的簡化，去掉模型結構上的圓角、倒角及對結構受力不影響的部分；其次根據APDL參數化語言，并通過定義參數化設計變量，建立模型的幾何約束和尺寸間約束關系，實現僅由幾個參數驅動整個模型的建立。建成后模型如圖5。

圖5 結合梁整體有限元模型圖

2 鋼錨梁與鋼混結合梁耦合連接技術

鋼錨梁是主纜吊桿的主要受力結構，承受荷載很大。它與結合梁的連接模擬方式直接影響模型的計算結果。ANSYS通常的做法：將結合梁橫隔板與鋼錨梁耦合部位進行切割，使鋼錨梁與結合梁有共同的邊界，從而實現節點耦合。由于結合梁內結構復雜，采用切割的方法易導致結合梁出現壞單元，計算結果不收斂。無法準確模擬受力狀態。

為準確模擬結合梁與鋼錨梁的耦合，利用約束方程，把結合梁節點的自由度與鋼錨梁一個或者多個節點或單元的自由度通過單元形函數聯系起來，從而建立結合梁與鋼錨梁邊界單元與節點的耦合關系，將兩者受力狀態聯系起來[3-4]。ANSYS提供通過建立臨近區，即Adjacent Regions建立約束方程。通過選取鋼錨梁邊界處較密的節點與結合梁隔板處較稀的單元，把鋼錨梁的節點與結合梁橫隔板的單元自由度耦合起來， 耦合單元容差取0.35。最終建立1 937個節點耦合約束方程，如圖6、7。

圖6 鋼錨梁約束方程模型圖

圖7 鋼錨梁與結合梁整體約束方程模型圖

3 計算模型

3.1 網格劃分

結合梁及鋼錨梁、鋼錨箱均為鋼板焊接的組合鋼結構。用SHELL63單元劃分網格，SHELL63為4節點、6自由度三維彈性殼單元，既具有彎曲能力又具有膜力，可承受平面內荷載和法向荷載[5]。用映射網格劃分控制單元尺寸，單元尺寸設為50。經過劃分網格，共42 386個單元，43 006個節點。

單元的材料屬性設置選用線彈性，各向同性材料。鋼結構彈性模量E=2.1×105MPa，剪切模量G=8.1×104MPa，泊松比μ=0.3，質量密度ρ=7 850 kg/m3。其中結合梁鋼結構選用Q345qD，鋼錨梁、鋼錨箱選用Q420qD。Q345qD抗拉、抗壓強度基本容許應力均為σ0=196 MPa，容許剪應力τ=120 MPa。Q420qD抗壓、抗拉強度基本容許應力均為σ0=216 MPa，容許剪應力τ=140 MPa。

3.2 模型荷載及邊界條件

圖8 完成邊界條件加載后圖

結合梁段兩端部截面底緣采用全自由度約束，順橋向對稱中線采用順橋向約束，順橋向非對稱面采用豎向約束。橫橋向對稱中線采用橫橋向約束。因此，可準確模擬結合梁段在荷載作用下的受力及變形狀態吊桿力通過錨箱傳遞到錨梁，再由錨梁傳遞到橫隔板上。由整體模型知，吊桿力N=1 730 000 N，考慮實際情況將吊桿力轉化為均布力施加在錨墊板底面，其值q=17 610 385 N/m2，完成邊界條件加載后如圖8。

4 計算結果與分析

有限元計算前首先對求解控制器進行設置，其中模型求解控制選項選擇靜力分析；采用自動時間步技術；求解器選用PCG雅可比共軛梯度法，收斂準則采用位移收斂。計算結果如下：

(1) 計算結果表明，鋼箱梁拉應力均值約20 MPa左右，主拉應力峰值約185 MPa，小于抗拉強度設計值；壓應力均值約30 MPa左右，主壓應力峰值189 MPa，小于抗壓強度設計值，均滿足要求。

(2) 鋼錨箱最大主拉應力為62.8 MPa，最大主壓應力為89.6 MPa。

(3) 鋼錨梁主拉應力均值約70 MPa左右，主拉應力峰值257 MPa；主壓應力均值為50 MPa左右，主壓應力峰值約223 MPa。鋼錨梁最大主拉應力及主壓應力均出現在鋼錨梁與橫隔板相交的角點位置，范圍極小，屬應力集中現象，焊接構造能有效消除或減小應力集中現象。

(4) 橫隔板主拉應力均值約50 MPa左右，主拉應力峰值445 MPa；主壓應力均值約47 MPa左右，主壓應力峰值425 MPa。同鋼錨梁處應力集中現象相同，橫隔板最大主拉應力及主壓應力均出現在鋼錨梁與橫隔板相交的角點位置，范圍極小，屬應力集中現象，采取相應的焊接構造能有效消除或減小應力集中現象。

5 結 語

(1) 對于比較復雜的模型，通過APDL語言建模可提高效率，避免了模型復雜時，由于刪除或布爾操作導致大量元素編號變化時建模容易出錯的問題。

(2) 通過約束方程可準確地實現對復雜結構的有限元節點耦合，從而合理地模擬復雜結構構件之間的接點連接過程。

(3) 通過對鋼-混結合梁鋼結構有限元分析，提出結構應力，為結構設計提供參考依據；同時也為類似結構分析提供有益的參考。

[1] 張勝民.基于有限元軟件的結構分析[M].北京:清華大學出版社,2003.

[2] 周寧.ANSYS-APD高級工程應用實例分析與二次開發[M].北京：中國水利水電出版社,2007.[3] 吳慶鳴,陳永強.基于約束函數法的熱力耦合有限元分析[J].計算力學學報,2008，(2)：183-187.[4] 張效松,葉天麒.非連續邊界元-有限元耦合方法分析[J].固體力學學報,1998，(4)：336-340.

[5] Release 10.0 Documentation for ANSYS[CP/DK].ANSYS Inc, 2005.

Study and Application of Finite Element Analysis Technology on Steel-Concrete Composite Girder

LU Fei, GUAN Wei

(Northwest Engineering Co., Ltd., Xi'an 710065,China)

The suspended cable-stayed composite girder with three towers on the West 2nd Ring Road in Hanzhong City features complicated structure and difficult design. In this paper, the ANSYS boundary coupling mode is applied to build one finite element analysis model of the composite girder which satisfies the actual conditions to study the action characteristics of the steel box girder of the composite girder section under the flexible cable composite system in the large-span girder bridge projects. The parameterization design language of ANSYS is applied to compile the analysis documents and optimize the control documents as well. The key issues such as the potential local extreme stress and the local stability are analyzed. The optimum section shape of the composite girder is obtained through calculation. Accordingly, the beneficial reference is provided to the design of the complicated structures.Key words:steel-concrete composite girder; ANSYS; parameterization; constraining formula

2014-05-14

盧飛(1986- )，男，河南省開封市人，工程師，從事輔企設計工作.

TU393.3

A

10.3969/j.issn.1006-2610.2015.04.018