考慮主梁損傷的簡支箱梁剪力滯效應研究

時元緒，鄔曉光，鄧淇元，李院軍

考慮主梁損傷的簡支箱梁剪力滯效應研究

時元緒，鄔曉光，鄧淇元，李院軍

(長安大學 公路學院，陜西 西安 710064)

針對主梁在發生損傷情況下截面剛度將發生變化，為計算在役簡支箱梁發生主梁損傷情況下的考慮剪力滯效應的彎曲正應力，基于階梯形剛度模型撓度計算方法，并引入一個主梁抗彎剛度折減系數ζ，將考慮主梁損傷后的計算撓度和主梁抗彎剛度折減系數代入剪力滯變分方程中，得到一種主梁發生損傷后考慮剪力滯效應的簡支箱梁的彎曲正應力公式，并分析主梁不同剛度對上翼緣彎曲正應力的影響。為驗證此方法的正確性，將用本文計算方法得到的翼緣板彎曲正應力與實體有限元結果、實橋數據進行對比。研究結果表明：此方法的計算結果與實體有限元結果吻合度較高，與實測值也基本吻合。得出隨著主梁剛度不斷折減，上翼緣應力值增幅不斷加大，當主梁剛度折減到60%時，上翼緣板應力值增幅均大于5%。

橋梁工程；簡支箱梁橋；剪力滯效應；有限元；變分法；階梯形剛度模型

箱形截面翼緣板與腹板在對稱豎向荷載作用下，基本不會產生相對變形，但腹板和翼緣板之間產生的縱向剪力流將引起剪切應力與相應的翼緣板的翹曲。由于寬翼緣板寬度方向上剪切變形的非均勻分布，引起受荷載時遠離腹板部分的縱向位移滯后于靠近腹板處的縱向位移，這一現象被稱為“剪力滯后現象”[1]。近幾年來，對剪力滯問題的研究主要集中在基于最小勢能原理的變分法[2?4]；基于剪切變形規律的有限梁段法[5?7]。在翹曲位移函數的研究中，較多研究不同翹曲位移函數對剪力滯效應的影響[8]。但是，對于在役橋梁出現裂縫等損傷情況下的剪力滯效應分析較少，其中，曹國輝等[9]基于換算截面法、運用變分原理分析開裂簡支箱梁剪力滯效應；何志剛等[10]分析了配筋以后混凝土開裂狀態下的剪力滯效應。本文基于階梯形剛度模型撓度計算方法，考慮主梁損傷對在役簡支箱梁剪力滯效應的影響，運用變分原理，計算在役簡支箱梁發生主梁損傷情況下的剪力滯解析解。提出一種適合于在役簡支箱梁簡便、實用的剪力滯計算方法，并分析了主梁發生不同程度損傷對翼緣板彎曲應力的影響。

1 基于階梯型剛度模型計算主梁損傷撓度

基于文獻[11]，雖然裂縫的隨機性較大，但從統計學觀點來看，裂縫的某些特征又具有一定規律性。

式中：cri為各條裂縫的高度；cri為各條裂縫的間距；cr為裂縫的開裂范圍；cri為各條裂縫開展的程度；為裂縫數量。

由式(2)，式(3)和彎矩與曲率的關系，最終求得簡支梁開裂截面的損傷剛度：

式中：為在役橋梁某種工況下外力對開裂截面換算截面形心軸產生的彎矩。

根據梁段開裂裂縫特性，其開裂梁段彎曲裂縫與剛度區間劃分如圖1和圖2所示。計算橋梁損傷情況下的考慮剪力滯的彎曲正應力公式時，可根據梁段不同程度損傷進行損傷剛度區間計算。

圖1 梁開裂后彎曲裂縫示意圖

圖2 梁開裂后區間劃分示意圖

如圖2所示，從左邊開始計算，由材料力學相關知識可知：

撓度與彎矩的關系：

考慮主梁損傷剛度折減以后，通過分段積分即可求得開裂后各段撓度值：

2 主梁抗彎剛度折減系數ζi

基于鄔曉光等[12]的研究，由于開裂帶來的剛度折減效應，在荷載作用下主梁撓度會逐漸增加，且撓度值可由2部分組成，分別為梁體無損傷情況下產生的撓度值和梁體發生損傷時產生的附加撓度值。其中以剛度分配系數來表征梁體剛度折減的程度，計算如式(9)和式(10)所示。

式中：為梁體出現損傷時，某一區段實際產生的撓度值，其結果可由式(5)，(6)，(7)和(8)計算得到。

因此，可用撓度值的增加來反映梁體剛度的折減情況：

式中：EI為梁體階梯剛度中的剛度；為梁體未發生損傷時的剛度。

由式(11)可知，

由(10)式，便可最終計算出主梁抗彎剛度折減系數：

在式(13)中，一方面：和可以從上述理論可以計算得到；另一方面，可以通過在役損傷橋梁的荷載試驗中的撓度測得，可以通過剛通車橋梁的荷載試驗中的撓度測得。通過荷載試驗中的撓度值進行計算，更加方便準確。

3 損傷情況下的箱梁剪力滯變分方程

在箱梁彎曲的剪力滯分析中，將考慮主梁損傷后的計算撓度w和主梁剛度折減系數引入到箱梁剪力滯變分方程中，推導考慮主梁損傷情況下的箱梁剪力滯彎曲正應力方程。其箱梁截面尺寸示意如圖3所示。

圖3 箱梁截面尺寸

式中：()為考慮主梁開裂的撓度；考慮(,)為梁的縱向位移；()為翼緣板剪切轉角的最大差值；為箱室凈寬的一半；h為截面形心到頂、底板的距離。

梁受到彎曲時的外力勢能：

式中：()為考慮主梁開裂的撓度。

腹板及上下翼板的應變能為：

式中：V為腹板應變能；V為上翼板應變能；V為下翼板應變能；I為腹板抗彎慣性矩；t為上翼緣板厚度；t為下翼緣板厚度；為主梁抗彎剛度折減系數；為考慮主梁開裂的撓度

式(19)中的第3式為邊界條件，并令

由式(19)和(20)可得：

式(21)的通解為

由式(19)第一式得到

由式(23)可知，橋梁損傷以后的撓度公式中考慮了主梁抗彎剛度折減系數。

在考慮剪力滯之后的翼板彎曲正應力為：

由式(24)可知，主梁發生損傷情況下的考慮剪力滯影響的翼緣板彎曲正應力比主梁正常情況下的彎曲正應力公式中多一個主梁抗彎剛度折減系數，且翼緣板的縱向位移函數中的中也多一個主梁抗彎剛度折減系數。

4 算例分析

4.1 實體有限元分析

基于羅旗幟[13]中的簡支箱梁有機玻璃模型，在箱梁跨中附近出現損傷，本文采用3組不同的剛度，表征橋梁發生3種不同程度的損傷，其剛度值為80%，70%和60%。剛度變化示意圖如圖4。

圖4 簡支箱梁剛度變化示意圖

采用實體有限元法(FEA)及本文提出的計算方法分別主梁翼緣板的剪力滯效應進行求解，并分別與實體有限元法計算結果進行對比，以驗證本文引入主梁抗彎剛度折減系數的考慮剪力滯的彎曲正應力方程的正確性。有機玻璃模型簡支箱梁材料彈性模量取=3 000 MPa，泊松比為0.385，橋梁跨徑=0.8 m，上翼緣共有8個計算點，每個計算點間距為20 mm，其結構、荷載布置形式及計算點如圖5所示。

采用Midas FEA實體有限元軟件分別建立3組不同剛度的模型，并建立1組完全沒有損傷的模型作為對照，模型中翼緣板劃分成邊長3 mm的四面體實體單元，腹板劃分成邊長2.5 mm的四面體實體單元，模型共劃分177 149個單元，有限元模型如圖6所示。

單位：mm

圖6 實體有限元模型

由4組簡支箱梁有限元模型得到4組撓度值，取=420 mm處截面的撓度，并將撓度值代入式(13)中得到主梁抗彎剛度折減系數，將撓度值與主梁抗彎剛度折減系數匯總如表1所示。

表1 主梁抗彎剛度折減系數

由表1可得出：本文用抗彎剛度折減系數公式計算得到的剛度值與模擬剛度值誤差較小，結果較為吻合，證明本文抗彎剛度折減系數公式的正確性與可靠性。

由本文提出的計算公式對上述4組簡支箱梁上翼緣應力進行計算，取=420 mm處截面的上翼緣。其計算過程如下：

當

當

根據邊界條件，

可得：

代入式(26)和(27)可得到

將本文計算方法算得的=420 mm截面頂板上8個計算點的正應力絕對值與實體有限元計算結果絕對值進行對比，如表2所示。

由表2可知，在3組不同折減剛度情況下，實體有限元的結果與本文計算值誤差較小，應力值較為吻合，說明本文方法計算結果的正確性。另外，隨著剛度不斷折減，上翼緣應力值也不斷增加，將不同折減剛度相對于未折減剛度的應力增加幅值匯總如表3所示。

表2 簡支箱梁(x=420 mm)上翼緣應力值(MPa)

表3 不同折減剛度相對于未折減剛度應力增幅

由表3可知，主梁發生損傷后，上翼緣應力值遠離腹板增幅較大，靠近腹板增幅較小。當主梁剛開始損傷時，即剛度折減到80%，總體應力增幅均較小，可以忽略不計；而隨著損傷進一步加大，其增幅越來越顯著；當剛度折減到60%時，上翼緣腹板附近應力增幅均大于5%，對于結構產生不利影響，該應力增幅不可忽略。因此，當主梁剛度折減到80%時，此時主梁損傷程度不大，需要及時修復主梁損傷部位，以避免應力值繼續增大，導致橋梁 損壞。

4.2 實橋數據分析

本文以2018年9月15日在貴州省威寧縣迤那鎮實測的一座裝配式簡支箱梁橋的荷載試驗數據來驗證本文計算方法的正確性。試驗橋基本概況：全橋由3跨組成，兩邊跨為空心板，中跨為箱梁，其跨徑為30 m，共6片主梁，每片單梁高為1.6 m，寬為3.3 m，橋寬20 m，混凝土等級為C50，在橋梁外觀檢測中發現箱梁跨中附近出現損傷，具體損傷情況見表4。中跨中梁應變測點布置如圖7所示。

表4 貴州省威寧縣林場中橋主梁病害調查統計表

注：部件編號中第1個數字表示第幾跨，第2個數字為第幾片梁；為裂縫長度，為裂縫寬度；′為裂縫平均長度，′為裂縫平均寬度。

單位：mm

本文荷載數據采用荷載試驗的中載工況，該工況采用4輛總重為400 kN的三軸載重汽車，所以載工況汽車總荷載為1 600 kN，作用形式為集中荷載，具體汽車荷載布置位置如圖8所示。

根據剛成橋時的撓度與橋梁發生損傷的撓度計算簡支箱梁中梁的抗彎剛度折減系數，即

并根據中梁的撓度與各主梁的總撓度之比，求出作用在中梁上的集中荷載占中載工況中總荷載的比例，根據比例可將中載工況中的總荷載向中梁進行分配，以確定中梁所受的集中荷載，即

式中：為總荷載；y為第片梁的撓度；3為中梁所受總荷載。

單位：mm

圖8 中載工況加載布置圖

Fig. 8 Loading layout for medium load conditions

已知抗彎剛度折減系數與所受荷載，通過本文計算式(32)便可以計算中梁跨中截面的上翼緣正應力值。

將本文計算公式算得的結果與實測應力進行對比，結果匯總如表4所示。

表4 本文計算應力與實測應力對比(MPa)

由表4可知，中梁在跨中1 m范圍內出現裂縫，主梁發生損傷情況下，本文方法計算值與實測數據誤差小于10%，數據基本吻合，證明本文方法具有一定的可靠性。

5 結論

1) 隨著主梁損傷程度越大，即主梁剛度越小，上翼緣板應力越來越大，應力增幅越來越顯著，當主梁剛度折減到80%時，需要及時采取措施修復主梁損傷部位。

2) 引入了主梁抗彎剛度折減系數的概念，該系數為小于1的數，可由主梁各區段實際所測撓度值來反映，實際撓度值越大，主梁抗彎剛度折減系數越小，反之則越大。

3) 將提出的主梁抗彎剛度折減系數代入到剪力滯變分方程中，最終得到主梁發生損傷情況下的考慮剪力滯效應影響后的彎曲正應力方程，為計算橋梁發生損傷后翼緣板彎曲正應力提供了一種 思路。

4) 代入提出的主梁抗彎剛度折減系數所計算出的應力值與有限元結果、實測數據基本吻合，說明本文提出的計算方法對于計算主梁發生損傷情況的剪力滯效應有一定的可靠性。

[1] LIN Zhibin, ZHAO Jian. Revisit of AASHTO effective flange-width provisions for box girders[J]. Journal of Bridge Engineering, 2011, 16(6): 881?889.

[2] QIN X X, LIU H B, WANG S J, et al. Symplectic analysis of the shear lag phenomenon in a T-Beam[J]. Journal of Engineering Mechanics, 2015, 141(5): 1279? 1284.

[3] ZHOU Shijun. Shear lag analysis in prestressed concrete box girders[J]. Journal of Structural Engineering, 2011, 16(4): 500?512.

[4] ZHANG Yuanhai. Improved finite-segment method for analyzing shear lag effect in thin-walled box girders[J]. Journal of Structural Engineering, 2012, 138(10): 1279? 1284.

[5] 周朋, 藺鵬臻. 變截面連續箱梁剪力滯分析的有限梁段法[J]. 鐵道科學與工程學報, 2018, 15(4): 948?955. ZHOU Peng, LIN Pengzhen. Finite segment method of shear lag analysis of continuous box girders with variable height[J]. Journal of Railway Science and Engineering, 2018, 15(4): 948?955.

[6] 周朋, 藺鵬臻, 何志剛. 基于剪切變形規律的剪力滯分析的有限梁段法[J]. 鐵道科學與工程學報, 2018, 15(1): 103?109. ZHOU Peng, LIN Pengzhen, HE Zhigang. Finite segment method of shear lag analysis of box girders based on flange-slab shear deformation law[J]. Journal of Railway Science and Engineering, 2018, 15(1): 103?109.

[7] 張元海, 李琳, 林麗霞, 等. 以附加撓度作為廣義位移時薄壁箱梁剪力滯效應的梁段有限元分析[J]. 土木工程學報, 2013, 46(10): 100?107. ZHANG Yuanhai, LI lin, LIN lixia, et al. Beam-segment finite element analysis on shear lag effect of thin-walled box girder adopting additional deflection as generalized displacement[J]. China Civil Engineering Journal, 2013, 46(10): 100?107

[8] 肖軍, 李小珍, 劉德軍, 等. 一種剪滯翹曲位移函數的解析構造法[J]. 哈爾濱工業大學學報, 2017, 49(3): 162?167. XIAO Jun, LI Xiaozhen, LIU Dejun, et al. An analytical construction method of wraping displacement function of shear lag[J]. Journal of Harbin Institute of Technology, 2017,49(3): 162?167.

[9] 曹國輝, 方志. 變分原理分析開裂簡支箱梁剪力滯效應[J]. 計算力學學報, 2007, 24(6): 853?858. CAO Guohui, FANG Zhi. Shear lag effect of cracked simply supported box girderanalyzed by variational principle[J]. Chinese Journal of Computational Mechanics, 2007, 24(6): 853?858.

[10] 何志剛, 劉應龍, 周朋, 等. 混凝土箱梁開裂的剪力滯效應分析[J]. 鐵道建筑, 2017, 57(9): 25?30. HE Zhigang, LIU Yinglong, ZHOU Peng, et al. Analysis on shear lag effect of concrete box girder after cracking[J]. Railway Engineering, 2017, 57(9): 25?30.

[11] 董水英. 預應力混凝土簡支梁開裂后撓度模型試驗研究[D]. 西安: 長安大學, 2010. DONG Shuiying. Model experimental study of deflection of the cracked prestressed concrete beam[D]. Xi’an: Chang’an University, 2010.

[12] 鄔曉光, 錢若霖, 蘇佩. 考慮剛度損傷的裝配式T梁橋橫向分布計算[J]. 公路交通科技, 2018, 35(6): 61?66. WU Xiaoguang, QIAN Ruolin, SU Pei. Calculation of lateral distribution of prefabricated T beam bridge considering stiffness damage[J]. Highway Traffic Technology, 2018, 35(6): 61?66.

[13] 羅旗幟. 基于能量原理的薄壁箱梁剪力滯理論與試驗研究[D]. 長沙: 湖南大學, 2005. LOU Qizhi. Theory and model test studies of shear lag in thin walled box girders based on energy principle[D]. Changsha: Hunan University, 2005.

Research on shear lag effect of simply supported box girder considering main girder damage

SHI Yuanxu, WU Xiaoguang, DENG Qiyuan, LI Yuanjun

(School of Highway, Chang’an University, Xi’an 710064, China)

The stiffness of section will change under the condition of the main girder damage. In order to calculate the bending normal stress considering shear lag effect in the case of damaged main girder of simply supported box girder in service, this paper presented a deflection calculation method based on ladder stiffness model and introduced a reduction coefficientζof bending stiffness of main girder considering the damage of main girder. The calculated deflection and reduction coefficient of bending stiffness of the main girder considering damage of the main girder were substituted into the shear lag variational equation, and a formula of bending normal stress of simply supported box girder considering shear lag effect after damage of the main girder was obtained, and the influence of different bending stiffness of the main girder on the bending normal stress of the upper flange was analyzed. In order to verify the correctness of this method, bending stress of flange plate calculated by method of this paper was compared with the results of solid finite element method and real bridge data. The results show that the calculated results of this paper’s method are in good agreement with those of the solid finite element method and with the measured values. In addition, this paper draws the conclusion that with the reduction of the stiffness of the main girder, the increase of the stress value of the upper flange is increasing. When the stiffness of the main girder is reduced by 60%, the increase of the stress value of the upper flange plate is more than 5%.

bridge engineering; simply supported box girder bridge; shear lag effect; finite element; variational method; stepped stiffness model

U448.21

A

1672 ? 7029(2019)08? 1998 ? 08

10.19713/j.cnki.43?1423/u.2019.08.017

2018?11?13

山西省交通運輸廳科技資助項目(2017-1-37)

鄔曉光(1961?)，男，湖北英山人，教授，博士，從事橋梁設計理論與橋梁檢測加固；E?mail：wxgwst.cn@126.com

(編輯 蔣學東)