鋼桁腹式混凝土組合箱梁的扭轉(zhuǎn)效應(yīng)分析

楊霞林， 于小芹， 張?jiān)?/p>

(1.蘭州交通大學(xué) 土木工程學(xué)院，蘭州 730070； 2.山東省建筑工程質(zhì)量檢驗(yàn)檢測中心有限公司，濟(jì)南 250031)

1 引 言

鋼桁腹式混凝土組合箱梁是一種較新型的鋼-混組合箱梁結(jié)構(gòu)，具有現(xiàn)場施工便捷、自重輕、避免腹板開裂、通透性好以及橋型優(yōu)美等優(yōu)點(diǎn)。我國已有的該種橋型最為典型的是南京江山橋，已建成的該種橋型還有國道107寶安段匝道以及水碾堡特大橋。目前，國內(nèi)外學(xué)者對該類結(jié)構(gòu)的受力特點(diǎn)進(jìn)行了較為深入的研究與探索，但關(guān)于扭轉(zhuǎn)性能方面的研究涉及很少[1-5]，王彤等[6]利用換算薄壁箱梁法思想及箱梁結(jié)構(gòu)理論，提出了桁腹式組合桁架彎曲變形、扭轉(zhuǎn)和畸變的計(jì)算方法。國內(nèi)外對其他類型的鋼-混組合箱梁的扭轉(zhuǎn)性能研究比較成熟[7-9]，故本文以薄壁箱梁扭轉(zhuǎn)理論為基礎(chǔ)，結(jié)合現(xiàn)有鋼-混組合箱梁的扭轉(zhuǎn)研究方法對鋼桁腹式混凝土組合箱梁進(jìn)行扭轉(zhuǎn)效應(yīng)的初步探究。

首先，按照剪切變形相等的原則將不連續(xù)鋼桁腹桿轉(zhuǎn)化為連續(xù)的薄壁混凝土腹板，并推導(dǎo)其各項(xiàng)扭轉(zhuǎn)幾何特性值。然后，基于薄壁箱梁扭轉(zhuǎn)理論，推導(dǎo)出組合箱梁閉口斷面的混凝土頂?shù)装搴蛽Q算鋼腹板的扭轉(zhuǎn)翹曲應(yīng)力表達(dá)式，進(jìn)而推導(dǎo)出組合箱梁約束扭轉(zhuǎn)控制微分方程；利用初參數(shù)法求解微分方程，并分析出翹曲雙力矩以及扭轉(zhuǎn)翹曲正應(yīng)力隨梁跨的變化規(guī)律。最后，通過有限元模擬分析，將有限元值和理論值進(jìn)行比較，以驗(yàn)證理論計(jì)算方法的可靠性與準(zhǔn)確性。

為便于開展相應(yīng)的分析與研究，本文引入以下基本假定。 (1) 斜腹桿為二力桿，僅承受軸力，不考慮失穩(wěn)； (2) 上部荷載作用下的彎矩由頂?shù)装宄袚?dān)，腹桿僅承受剪力且受力均勻； (3) 組合箱梁在偏載作用下發(fā)生扭轉(zhuǎn)時(shí)，忽略鋼桁腹桿自身的扭轉(zhuǎn)； (4) 符合小變形假定，且各個(gè)桿件的受力均在彈性范圍內(nèi)； (5) 橫斷面的縱向應(yīng)力分布符合擬平截面假定。

2 鋼桁腹桿縱向表觀彈性模量

由于鋼桁腹式混凝土組合箱梁的腹桿不連續(xù)，且腹桿與混凝土頂?shù)装宀牧喜灰恢拢虼诵枰惹蠼怃撹旄箺U的縱向表觀彈性模量。

相鄰兩腹桿相接節(jié)點(diǎn)通過剪力鍵鉗固于混凝土承托中，因此混凝土翼板所受的力能夠平順地傳遞到鋼桁腹桿，可將混凝土頂?shù)装逯g的鋼桁腹桿看作連續(xù)折線。任取一組折線單元如圖1所示，每組折線單元以兩個(gè)彎矩零點(diǎn)為界限[3]。

利用圖乘法計(jì)算單元在軸向力作用下的縱向變形Δ1。圖2為乘法受力圖。

圖1 腹桿折線單元

圖2 腹桿折線單元受力

折線單元縱向變形Δ1為

(1)

式中E0為鋼桁腹桿初始彈性模量(即鋼材彈性模量),I0為鋼桁腹桿截面慣性矩，A0為鋼桁腹桿截面面積，P為鋼桁腹桿所受軸力，h為鋼桁腹桿豎向高度，θ為相鄰兩鋼桁腹桿之間的夾角。

取一塊等效平鋼板代替鋼桁腹桿[6]，假設(shè)此平鋼板厚度為ts，則

ts=E0/Gs·(A0hm/L3)

(2)

式中Gs為鋼材剪切模量，m為折線單元縱向長度的一半，L為鋼桁腹桿長度。

在作用力P作用下，等效平鋼板的縱向位移為

(3)

式中Ex為鋼桁腹桿縱向表觀彈性模量，As為等效鋼腹板的截面積。

由Δ1=Δ2，聯(lián)立式(1～3)，得鋼桁腹桿縱向表觀彈性模量Ex為

(4)

3 扭轉(zhuǎn)荷載分解

圖3所示為鋼桁腹式混凝土組合箱梁的橫斷面，斷面上作用有豎向偏載時(shí)會(huì)發(fā)生扭轉(zhuǎn)變形，一般而言，將該偏心荷載分解為一對對稱荷載與一對反對稱荷載[10]。組合箱梁在反對稱荷載作用下，截面會(huì)產(chǎn)生剛性扭轉(zhuǎn)與畸變變形。為阻止組合箱梁截面發(fā)生畸變變形，在截面施加一虛擬的對角支撐，使其只發(fā)生剛性扭轉(zhuǎn)變形[9]，如圖3所示。

圖3 反對稱荷載的分解

根據(jù)靜力平衡條件，得扭轉(zhuǎn)荷載為[11]

(5)

4 鋼桁腹式混凝土組合箱梁的扭轉(zhuǎn)分析

4.1 扭轉(zhuǎn)幾何特性

分析鋼桁腹式組合箱梁的約束扭轉(zhuǎn)翹曲應(yīng)力，需先計(jì)算鋼桁腹式組合箱梁的截面扭轉(zhuǎn)幾何特性(包括扭心、輔助扇性坐標(biāo)、主扇性坐標(biāo)、主扇性慣性矩、主扇性靜矩和廣義扇性靜矩)。

圖4的點(diǎn)1為混凝土頂板與鋼桁腹桿交接處參考計(jì)算點(diǎn)，點(diǎn)2為頂部翼緣板端部參考計(jì)算點(diǎn)，點(diǎn)3為底板與鋼桁腹桿交接處參考計(jì)算點(diǎn)。輔助極點(diǎn)選取頂板中心點(diǎn)B處，計(jì)算控制點(diǎn)1，2和3輔助扇性坐標(biāo)ωB、主扇性坐標(biāo)ωA、主扇性靜矩Sω以及扭心位置R和主扇性慣性矩Iω。各參數(shù)(寬度、厚度和傾角)均在圖4中標(biāo)明。

圖4 鋼桁腹式混凝土組合箱梁截面

點(diǎn)1，2和3的輔助扇性坐標(biāo)表達(dá)式為

(6)

扭心R到頂板中心的距離d為

d=-Iω B x/Iy=h/[a3ta+c3tb+(b3-c3)tc+

(7)

式中Iω B x為輔助扇性坐標(biāo)圖與x坐標(biāo)圖圖乘所得，Iy為組合箱梁截面對x軸的慣性矩，

點(diǎn)1，2和3的主扇性坐標(biāo)表達(dá)式為

(8)

主扇性慣性矩可通過將主扇性坐標(biāo)圖自乘求得，

(9)

式中nc為混凝土彈性模量與鋼桁腹桿縱向表觀彈性模量之比，即nc=Ec/Ex。

主扇性靜矩計(jì)算式為

(10)

4.2 自由扭轉(zhuǎn)

組合箱梁橫斷面扭轉(zhuǎn)變形的截面位移可以表示為[10]

(11)

根據(jù)式(11)，可得組合箱梁橫截面的正應(yīng)力和剪應(yīng)力為

(12a)

(12b)

式中v為泊松比。

純扭轉(zhuǎn)翹曲位移為

(13)

組合箱梁的自由扭轉(zhuǎn)微分方程為

(14)

式中Mz為集中扭矩，mz為分布扭矩，Id為扭轉(zhuǎn)慣性矩。

(15)

式中mG=Gc/Gs，式(15)表示若將換算鋼腹板進(jìn)一步等效換算為混凝土腹板，可將其按混凝土薄壁箱梁計(jì)算Id。Id通常由閉口箱室與開口懸臂板迭加求得[12]，而在計(jì)算組合箱梁的Id時(shí)，開口部分的扭轉(zhuǎn)剛度很小，為計(jì)算簡便，一般忽略不計(jì)。

4.3 約束扭轉(zhuǎn)

根據(jù)薄壁箱梁設(shè)計(jì)理論，需要引入位移量β′(組合箱梁翹曲位移自由度)對翹曲位移計(jì)算進(jìn)行修正，使得組合箱梁發(fā)生約束扭轉(zhuǎn)后仍能保持變形后的連續(xù)。

翹曲位移表達(dá)式為

(16)

翹曲正應(yīng)力的表達(dá)式可以表示為

(17)

對于鋼桁腹式混凝土組合箱梁，頂?shù)装搴蛽Q算鋼腹板由于其材料不同，各自翹曲正應(yīng)力計(jì)算式也不相同。

混凝土翼板約束扭轉(zhuǎn)正應(yīng)力計(jì)算式為

(18)

式中Ec為混凝土的彈性模量。

換算鋼腹板的扭轉(zhuǎn)翹曲正應(yīng)力計(jì)算式為

(19)

式中Ex為鋼桁腹桿的縱向表觀彈性模量。

翹曲正應(yīng)力σω對應(yīng)的翹曲雙力矩Bω是一對大小相等方向相反的力矩，計(jì)算式為

(20)

式中Iω為主扇性慣性矩，

(21)

式中A1，A2和A3為混凝土頂板和底板截面面積以及換算腹板截面面積。

由式(18～20)，翹曲正應(yīng)力表示為

(22)

約束扭矩Mz可由翹曲雙力矩Bω表示為

(23)

5 約束扭轉(zhuǎn)微分方程及其參數(shù)解

5.1 約束扭轉(zhuǎn)微分方程

雙力矩Bω、扭矩Mz以及翹曲正應(yīng)力σω的求解需要先求出撓曲變形系數(shù)β′，然后通過β和θz的關(guān)系求解θz。

截面總扭矩為

(24)

根據(jù)式(24)可得β和θz的關(guān)系為

(25)

(26)

Mz為集中扭矩，mz為分布扭矩，兩者之間的關(guān)系為

dMz/dz=-mz

(27)

對式(26)求導(dǎo)并聯(lián)立式(27)，經(jīng)整理得β的微分方程為

(28)

同理可推導(dǎo)出θz的微分方程為

(29)

5.2 初參數(shù)法求解微分方程

規(guī)定z=0處，四個(gè)初參數(shù)中θ0為該處角位移，β′0為翹曲位移，Bω 0為扭轉(zhuǎn)翹曲雙力矩，M0為扭矩，微分方程的齊次方程的通解為

θz=C1+C2z+C3sinhkz+C4coshkz

(30)

(31)

Bω=-GId(C3sinhkz+C4coshkz)

(32)

Mz=C2GId

(33)

初參數(shù)解為

(34)

(35)

(36)

Mz=M0

(37)

在式(34～37)中，組合箱梁的邊界條件決定4個(gè)初參數(shù)的取值，初參數(shù)解適用于無外界荷載作用下的組合箱梁，一般而言，對于作用有外荷載的組合箱梁，應(yīng)在初參數(shù)解的基礎(chǔ)上補(bǔ)充相應(yīng)的影響項(xiàng)。

如圖5所示，當(dāng)組合箱梁跨內(nèi)有扭矩時(shí)，式(34～37)分別為

(38)

(39)

(40)

(41)

式中符號=a表示只有z>a時(shí)的數(shù)值才計(jì)入算式，同樣=b表示只有z>b時(shí)的數(shù)值才計(jì)入算式，當(dāng)z>c時(shí)，對應(yīng)的積分上限也會(huì)發(fā)生改變。

5.3 邊界條件

邊界條件不同，角位移θz以及撓曲變形系數(shù)β′的取值也不相同[13]。

(1) 固定端θz=0，β′=0

(2) 簡支端θz=0，β″=0

5.4 作用集中扭矩的簡支組合箱梁扭轉(zhuǎn)分析

在簡支組合箱梁中，距梁端a處作用一集中扭矩，如圖6所示，根據(jù)邊界約束條件θz=0，β′=0，z=0時(shí)，θz=0，Bω 0=0，可得微分方程解為

當(dāng)z

(42)

當(dāng)z≥a時(shí)，

Mz=-(a/l)T

(43)

圖6 簡支組合箱梁作用集中扭矩

6 算例分析

6.1 扭轉(zhuǎn)幾何特性計(jì)算

以某鋼桁腹式混凝土組合箱梁簡支梁為例[14]，其上部結(jié)構(gòu)采用單跨(35 m)等截面鋼桁腹式預(yù)應(yīng)力混凝土組合箱梁結(jié)構(gòu)，主梁為單箱單室截面，頂板寬8.5 m，底板寬4.8 m，梁高2.3 m?；炷涟宀捎肅50混凝土，彈性模量為3.45×104MPa，鋼桁腹桿采用Q345C級鋼管，規(guī)格為Φ351×16，鋼桁腹桿水平傾角為67° 左右，節(jié)間距為1.95 m。組合箱梁橫斷面和縱斷面分別如圖7和圖8所示。為簡化計(jì)算，分別將頂板和底板的板厚均分，將鋼桁腹桿轉(zhuǎn)換成薄壁混凝土腹板，橫斷面簡化為如 圖9 所示截面。按上述理論推導(dǎo)式(6～10)和式(15)求得扭轉(zhuǎn)幾何特性，結(jié)果見表1及圖10～圖13。

圖7 組合箱梁橫斷面(單位：m)

圖8 組合箱梁縱斷面(單位：m)

圖9 組合箱梁計(jì)算模型及橫向加載位置(單位：mm)

表1 扭轉(zhuǎn)幾何特性

圖10 輔助扇性坐標(biāo)ωB(單位：m2)

圖11 主扇性坐標(biāo)ω(單位：m2)

圖12 主扇性靜矩單位：m4)

圖13 廣義扇性靜矩單位：m4)

6.2 結(jié)果分析

依據(jù)《公路橋規(guī)》[15]中車道荷載的情況對該橋施加P=330 kN的偏載，偏心距e=2075 mm，不計(jì)橋梁自重。橫向加載如圖14所示。

以此算例為研究對象，探究翹曲雙力矩隨著集中扭矩分別作用在簡支鋼桁腹式混凝土組合箱梁上l/4和3l/4跨以及跨中截面處的變化，圖15所示為翹曲雙力矩沿梁縱向的變化曲線?？梢?，翹曲雙力矩在集中扭矩作用處達(dá)到最大值，并且衰減速度相當(dāng)快，距離集中扭矩作用處越遠(yuǎn)，翹曲雙力矩就越小，梁端處為0。

以跨中作用集中扭矩為例，繪出點(diǎn)1(左側(cè)腹桿與混凝土頂板交點(diǎn))、點(diǎn)2(混凝土頂板左端點(diǎn))與點(diǎn)3(左側(cè)腹桿與混凝土底板交點(diǎn))翹曲正應(yīng)力沿梁縱向的變化曲線，如圖16所示，圖中應(yīng)力以拉為正，以壓為負(fù)。集中扭矩作用處的翹曲雙力矩很大，使得該處組合箱梁截面產(chǎn)生了很大的翹曲應(yīng)力，而距離集中扭矩作用處越遠(yuǎn)，組合箱梁的翹曲雙力矩越小，翹曲正應(yīng)力也越小，梁端處為0。

圖14 橫向加載位置

圖15 翹曲雙力矩沿梁縱向的變化曲線

圖16 扭轉(zhuǎn)翹曲正應(yīng)力沿梁縱向的變化曲線

6.3 對比分析

為驗(yàn)證以上理論計(jì)算方法的可靠性與準(zhǔn)確性，運(yùn)用ANSYS進(jìn)行有限元數(shù)值模擬計(jì)算，并進(jìn)一步與理論計(jì)算結(jié)果進(jìn)行對比分析。有限元模型橫向和縱向尺寸構(gòu)造如圖7和圖8所示，選用SOLID185單元模擬頂板和底板，BEAM188單元模擬鋼桁腹桿[16]，不考慮腹桿與混凝土的相對滑移，鋼桁腹桿與混凝土板的剛性連接通過建立約束方程進(jìn)行自由度耦合。簡支鋼桁腹式混凝土組合箱梁的邊界一端為活動(dòng)鉸支座，約束梁的豎向和橫向位移，另一端為固定鉸支座，約束梁的豎向和縱向和橫向位移。模型的相關(guān)數(shù)據(jù)列入表2。腹桿傾斜的角度為67° 左右，兩相鄰腹桿的最大距離為1.95 m。對該模型跨中截面處施加如圖17所示的一對反對稱荷載，使其滿足僅表現(xiàn)出扭轉(zhuǎn)與畸變的疊加效應(yīng)，而不包含彎曲效應(yīng)。

根據(jù)上述所推導(dǎo)的理論計(jì)算方法求取算例中鋼桁腹式混凝土組合箱梁的跨中截面處和0.1l跨截面處1點(diǎn)(腹桿與頂板連接點(diǎn))、2點(diǎn)(頂板左端點(diǎn))以及3點(diǎn)(腹桿與底板連接點(diǎn))的扭轉(zhuǎn)翹曲正應(yīng)力,疊加畸變翹曲正應(yīng)力(計(jì)算過程本文略)后與有限元計(jì)算結(jié)果進(jìn)行比較，列入表3。

圖17 反對稱荷載施加

表2 有限元模型的相關(guān)參數(shù)

表3 翹曲正應(yīng)力有限元數(shù)值與理論數(shù)值的對比(單位：kPa)Tab.3 Comparison between finite element value and theoretical value(unit：kPa)

由表3可知，有限元數(shù)值與本文方法計(jì)算的理論數(shù)值相差并不大，差值百分比在10%以內(nèi)，可見本文計(jì)算方法合理可行。

7 結(jié) 語

(1) 通過選取合適的輔助極點(diǎn)，推導(dǎo)出組合箱梁橫截面輔助扇性坐標(biāo)ωB、扭心R、主扇性坐標(biāo)ωA、主扇性慣性矩Iω及主扇性靜矩Sω等扭轉(zhuǎn)幾何特性表達(dá)式，并通過算例求解出扭轉(zhuǎn)幾何特性值，繪出輔助扇性坐標(biāo)圖、主扇性坐標(biāo)圖、主扇性靜矩以及廣義扇形靜矩圖。

(2) 對鋼桁腹式混凝土組合箱梁進(jìn)行自由扭轉(zhuǎn)與約束扭轉(zhuǎn)分析，推導(dǎo)出組合箱梁約束扭轉(zhuǎn)微分方程，并通過初參數(shù)法求解出轉(zhuǎn)角、翹曲位移、扭轉(zhuǎn)翹曲雙力矩以及扭矩的一般公式，以跨中作用集中扭矩的簡支組合箱梁為例，賦予邊界條件，求解出四個(gè)參數(shù)的表達(dá)式。

(3) 根據(jù)算例，分析翹曲雙力矩以及扭轉(zhuǎn)翹曲正應(yīng)力隨梁跨變化，翹曲雙力矩在集中扭矩作用處達(dá)到最大值，并且衰減速度很快，使得該處箱梁截面的翹曲應(yīng)力達(dá)到最大值，箱梁的翹曲雙力矩在遠(yuǎn)離集中扭矩作用處幾乎為0，翹曲正應(yīng)力也幾乎為0。

(4) 建立有限元模型施加荷載并求解出翹曲正應(yīng)力，利用本文所述方法求取跨中截面和0.1l跨截面點(diǎn)1，2和3的扭轉(zhuǎn)翹曲正應(yīng)力，與畸變翹曲正應(yīng)力疊加后與有限元值作比較，結(jié)果表明有限元數(shù)值與本文方法計(jì)算的理論數(shù)值差值百分比在10%以內(nèi)，可見本文推導(dǎo)的理論計(jì)算方法合理可行。

(5) 本文沒有考慮預(yù)應(yīng)力對扭轉(zhuǎn)效應(yīng)的影響，實(shí)際工程中鋼桁腹式混凝土組合箱梁常在上下翼板布置預(yù)應(yīng)力筋，故需進(jìn)一步研究預(yù)應(yīng)力產(chǎn)生的扭轉(zhuǎn)效應(yīng)。