錢 淼

(中鐵大橋勘測(cè)設(shè)計(jì)院集團(tuán)有限公司 武漢 430050)

波形鋼腹板組合箱梁橋充分利用混凝土承壓、鋼腹板承剪的特性，自重輕、結(jié)構(gòu)受力明確，是一種經(jīng)濟(jì)、合理的橋梁結(jié)構(gòu)[1-3]。目前，國內(nèi)外學(xué)者認(rèn)為波形鋼腹板組合箱梁橋的豎向彎曲剪應(yīng)力主要由鋼腹板承擔(dān),且在設(shè)計(jì)波形鋼腹板箱梁橋時(shí)，也只考慮鋼腹板承受豎向剪力，但實(shí)際上混凝土結(jié)構(gòu)也承受一部分豎向剪力[4-6]，精細(xì)化設(shè)計(jì)應(yīng)該考慮混凝土結(jié)構(gòu)承受部分豎向剪力。目前單箱多室箱梁腹板的豎向剪力分配計(jì)算難度大，所以有必要研究波形鋼腹板箱梁豎向剪應(yīng)力計(jì)算理論和開發(fā)相應(yīng)的軟件。

1 有限元計(jì)算算法

將文獻(xiàn)[7-9]的計(jì)算方法進(jìn)行改進(jìn)，使其適用于多材料結(jié)構(gòu)截面的橫向剪應(yīng)力計(jì)算，其受力模型，見圖1。

圖1 梁橫向受力圖

懸臂梁，自由端分別受X和Y方向的剪力Vx和Vy作用，采用翹曲函數(shù)表示截面的位移。

(1)

(2)

φ(x,y)-a1y-a2x

(3)

式中：u、v和w分別為X、Y、Z方向的位移分量；φ(x，y)為翹曲函數(shù)；L為梁?jiǎn)卧L(zhǎng)度。

(4)

式中：Vx和Vy分別為X和Y方向的剪力；a1、a2、a3、a4為常數(shù)，與梁的邊界條件相關(guān)。

(5)

(6)

(7)

(8)

橫向力作用下，桿的總勢(shì)能為

Ep=U-W

(9)

式中：U為應(yīng)變能，

(10)

其中：υz為鋼材泊松比。

(11)

(12)

W為橫向力作用下的外力做功。

(13)

文獻(xiàn)[7]、[10]將截面網(wǎng)格劃分為等參單元，令等參單元形函數(shù)為H(ξ，η)，ξ和η為自然坐標(biāo)，單元內(nèi)翹曲函數(shù)可以由形函數(shù)和節(jié)點(diǎn)翹曲函數(shù)獲得。

(14)

對(duì)翹曲函數(shù)求導(dǎo)

(15)

利用最小勢(shì)能原理δEp，由式(9)、式(10)和式(13)得截面在受剪時(shí)的平衡方程。

(16)

將式(15)代入(16)，得到梁在受到橫向力作用時(shí)的平衡方程表達(dá)式

Keψi=fs

(17)

式中：ψi為單元節(jié)點(diǎn)的翹曲函數(shù)值向量，是待求未知量；Ke為單元?jiǎng)偠染仃嚕籪s為在橫向力作用下的廣義力向量。根據(jù)下式計(jì)算

(18)

yCy)]detJdξdη

(19)

鋼腹板豎向剪應(yīng)力計(jì)算公式為

(20)

式中：υz為鋼材泊松比；Ez鋼材彈性模量；EBi為混凝土與鋼材彈性模量比；GBi為混凝土與鋼材剪切模量比。Ixx和Iyy為截面在過形心的x和y方向的換算截面慣性矩；Ixy為換算截面慣性積；A為換算截面面積；As為實(shí)際截面面積；Az為鋼材實(shí)際面積；Ai為混凝土實(shí)際面積。換算截面是將混凝土換算成鋼材。

2 計(jì)算實(shí)例

實(shí)例的波形鋼腹板箱形截面見圖2，為單箱三室箱形組合截面，頂?shù)装鍨殇摻罨炷?，鋼腹板?6 mm。受豎向剪力1 000 kN作用，求鋼腹板剪應(yīng)力。

圖2 波形鋼腹板箱梁截面(單位：m)

2.1 本文算法計(jì)算結(jié)果

根據(jù)本文算法編寫的軟件計(jì)算了頂板懸臂板寬度分別為4.24，1.74和0.24 m的3種截面豎向剪應(yīng)力分布，邊腹板豎向剪應(yīng)力分別為4.258，4.264，4.250 MPa，中腹板豎向剪應(yīng)力分別為4.242，4.256，4.266 MPa。根據(jù)計(jì)算結(jié)果頂板懸臂板的寬度越窄邊腹板分配的剪力越小，3種截面的邊中腹板剪力分配比分別為1.004，1.002，0.996。

2.2 ANSYS計(jì)算結(jié)果

根據(jù)有限元軟件ANSYS計(jì)算圖2截面，采用實(shí)體單元SOLID45建立100 m長(zhǎng)的懸臂梁，在自由端加1 000 kN豎向剪力，取懸臂梁中部查看豎向剪應(yīng)力，中腹板最大豎向剪應(yīng)力4.227 MPa，邊腹板最大豎向剪應(yīng)力為4.311 MPa。

2.3 簡(jiǎn)化算法計(jì)算結(jié)果

簡(jiǎn)化算法基本假定：豎向剪力全部由鋼腹板承受，且剪應(yīng)力沿高度均勻分布。根據(jù)簡(jiǎn)化算法，鋼腹板剪應(yīng)力為5.351 MPa。

3 結(jié)論

本文算法只需要截面幾何信息、材料信息和豎向剪力就可以獲得剪應(yīng)力分布結(jié)果，其值與有限元軟件ANSYS實(shí)體模型計(jì)算的剪應(yīng)力結(jié)果一致，且比簡(jiǎn)化算法計(jì)算結(jié)果小20%。

根據(jù)實(shí)例分析發(fā)現(xiàn)，雖然腹板的翼緣越寬其分配的剪力越大，但變化很小。按簡(jiǎn)化算法設(shè)計(jì)鋼腹板，明顯過于保守,會(huì)造成材料的浪費(fèi)。本文算法也適用于其它多材料的組合梁橫向剪應(yīng)力計(jì)算。