第6卷第4期 2009年l2月 长沙理工大学学报(自然科学版) Journal of Changsha University of Science and Technology(Natural Science) VoI.6 NO.4 Dec.20e 文章编号:1672~933l(2009)04—0016—05 薄壁箱梁剪力滞效应数值计算 徐飞鸿,蔡汶珊 (长沙理工大学土木与建筑学院,湖南长沙410004) 摘要:运用有限元方法,采用板壳单元——She儿63单元,对薄壁直线箱梁和薄壁曲线箱粱剪力滞效应分 别进行了数值计算.将直线箱梁剪力滞效应的数值计算结果与变分法理论计算值及模型试验值进行了对比, 三者吻合较好。验证了本研究数值方法的正确性.在有限元理论的基础上,进一步计算了曲线箱梁在静力荷 载作用下的挠度、应力、应变及剪力滞系数值,分析了曲率半径等因素对曲线箱梁剪力滞效应的影响.汁算结 果表明,曲率半径对曲线箱梁的剪力滞效应影响较大.与直线箱梁相比,截面相同位置处的剪力滞系数随曲 率半径的减小而增大,增幅远超过5 以上.因此在曲线箱梁的设汁中应对曲率半径加以考虑. 关键词:曲线箱形粱;剪力滞;数值计算;曲率半径 中图分类号:U448.2l3;TU3l7.1 文献标识码:A Numerical calculation on shear lag effect of thin—walled box girder XU Fei—hong,CAI Wen~shan (School of Civil Engineering and Architecture,Changsha University of Science and Fechnology.Changsha 410004,China) Abstract:The numerical calculation of the shear lag effect of the thin—walled straight box girder and the thin— wailed curved box girder is carried out respectively by using the finite el、・ ement method based on the board shell unit—Shell 63 unit.The computational results agree well with the predicted value of variational method and the experimental data,which vail— dates the exactness of the numerical method.Then primary research is carried out on tom— puting the amount of deflection,stress,strain and shear lag coefficient of the thin—walled curved box girder under static load arm analyzing the influence of curvature radius on the shear lag effect of the thin—walled curved box girder.The computational results show that cur— vature radius have significant influence on the shear lag effect of the thin-walled curved box girder, the shear lag coefficient increases as the curvature radius decreasing on the same position compared with the thin-walled straight box girder,and the amplitude increases more than 5 .Therefore, the curvature radius should be considered in the design of thin—walled curved box girders. Key words:curved box girder;shear lag;numerical calculation;curvature radius 曲线箱梁在荷载作用下将发生弯扭组合变 切、形心与剪心不重合等现象.由于梁轴初曲率的 存在,这些变形具有高度的耦合性.目前,常用的 形,同时伴随有截面的畸变、剪力滞后和翘曲剪 收稿日期:2009一l2一O6 基金项目:湖南省自然科学基金资助项目(o5JJ3OlO8) 作者简介:徐飞鸿(1962一),男,湖南益阳人,长沙理工大学教授.主要从事结构工程的研究 第6卷第4期 徐飞鸿,等:薄壁箱梁剪力滞效应数值计算 l7 箱梁分析方法有空间梁格法 ]、曲线梁段 法[4 ]、空间板壳有限元法L6]、空问实体有限元法 和试验方法等.对于曲线箱梁结构的分析,解析法 的各种理论存在着各种假设或忽略一些因素,如: 用于分析曲线梁桥的单纯扭转理论[7J,它把曲线 梁桥结构当作集中在梁轴中心线处的弹性杆件来 处理,并假定梁负载后横截面不产生翘曲和畸变, 仍保持平面.该理论一般仅适用于跨长大于横截 面尺寸四倍以上的曲线梁.再如:梁格系理论 将 桥梁上部结构用一个等效梁格来代替,由于梁格 系假定横梁刚性无穷大,故该理论只适合于宽跨 比较小的曲线梁桥.因此,假设的存在不仅使理论 的应用受到局限,也影响了计算结果的正确性.试 验方法虽然可靠性强,但是成本高、周期长.随着 计算机硬件与有限元方法的发展,使得采用数值 计算方法来研究曲线箱梁剪力滞效应成为可 能 ],数值计算方法可以指导试验,并能得到比较 全面的结论. 在我国现行公路桥梁规范中,仅对直线箱梁 因剪力滞效应而引出的翼板有效宽度做了计算规 定.工程人员在进行连续曲线箱梁设计时,通常先 按直梁计算翼板的有效宽度,然后予以折减,这使 计算的精度难以得到保证.因此,对曲线箱梁的剪 力滞效应的数值计算进行深入研究很有必要-l。 . 1剪力滞效应原理 在箱形梁中,产生弯曲的剪力流通过肋板传 递给翼板,剪力在翼板上分布是不均匀的,因此剪 切变形沿翼板的分布是不均匀的,从而引起弯曲 时远离肋板的翼缘纵向位移滞后于近肋板的翼缘 纵向位移,或者反之,这就导致其弯曲正应力的横 向分布呈曲线形.这种由于翼板的剪切变形所造 成的弯曲正应力沿梁宽方向不均匀的现象称为 “剪力滞效应”[1 ,前者称为“正剪力滞”,后者称为 “负剪力滞”(如图l所示). 通常用剪力滞系数来度量剪力滞效应的变化 规律,即 、一考虑剪力滞效应所求得的正应力 / 按初等梁理论所求得的正应力 箱梁主要由上翼缘、下翼缘、腹板和横隔板四 部分构成.由于上翼缘和下翼缘的厚度远小于梁 宽和梁轴线长,腹板和横隔板的厚度也远远小于 箱梁高度和梁轴线长,所以选择的结构单元为板 壳单元——She儿63单元.它是有限元软件AN— SYS中的弹性壳单元,既具有弯曲能力,具有薄膜 效应,可以承受平面内荷载和法向荷载.该单元共 有4个节点,每个节点具有6个自由度,即分别沿 节点坐标系z,Y, 方向的平动和转动. 正斡力滞 图1 剪力滞效应 Fig.1 Shear lag effect 本研究采用文献[12]中介绍的有机玻璃简支 梁模型作为计算模型,利用Shet[63单元建立对 应的有限元模型,在有限元分析程序ANSYS中 进行求解,将数值计算结果与变分法理论计算值 及试验数据进行比较,对计算单元的选取与数值 计算方法的正确性进行验正. 2 薄壁直线箱梁的剪力滞效应分析 2.1 有限元模型的建立 有机玻璃试验梁模型跨径为80 em,梁端不 设横隔板,弹性模量E一3 000 MPa,泊松比 一 0.385,自重忽略不计.在简支直线箱梁跨中正对 腹板的上翼板处对称作用集中荷载,总量为P一 0.272 2 kN,板中面的应变取上、下测点的平均 值.梁的截面尺寸和测点布置如图2所示. 试验梁模型的建模采用直接建模法,即先建 立节点,然后依据节点建立单元.利用Shell 63单 元建立的简支直线梁有限元计算模型如图3所 示,其中,节点共有748个,单元704个. 应用Shell 63单元对箱梁剪力滞效应进行数 值计算时,单元型式的定义、单元的合理离散以及 后处理中结果数据的提取与计算非常关键. l8 ..长沙理工大学学报(自然科学版) .。r十 t t,. ,士.L .『 。 2009年12月 l 2 3 4 5 6 7 8 :- 0 9 :10● :● ● ● I1l 12 13 14 l k————————一20 8————————— 图2模型梁的横截面尺寸及测点布置图(单位:cm) Fig.2 Cross dimensions of the model beam and the distributing of metrical points(Unit:em) 图3 简支直线箱梁有限元模型 Fig.3 Finite element model of freely supported straight box girder 2.2计算结果及分析 数值计算结果与文献[13]中的试验数据及变 分法理论计算值如表l所示. 表l 简支直线箱粱跨中截面正应力计算结果 Table l The normal stress calculation result in crosssection of the span centre of freely supported straight box girder MPa 从表1可以看出,利用ANSYS有限元分析 软件所求得的数值解与变分法理论计算值及试验 数据吻合较好,其中,数值解与变分法理论计算值 的误差值在2l%以内,与试验数据的误差值不超 过15.8 ,特别在底板部分以及顶板和腹板的交 界处,本研究计算值较文献[13]更加吻合,较好地 验证了本研究数值计算方法的正确性. 为了能更清楚地反映本研究数值计算值与变 分法理论值及试验值之间的接近程度,绘出横截 面的应力分布曲线(如图4和图5所示).由图4 和图5可见,用ANSYS数值求解箱粱剪力滞效应 的计算结果精度较高,能够真实反映横截面上正 应力的分布情况. ~O.15 -0-20 塞一O.25 一0.30 皇 .35 -0-4O -0.45 O 5 J0 】5 2O 25 3O 35 4O 节点横坐标,cm 图4 直线箱梁跨中截面顶板轴向正应力分布图 Fig.4 The axial normal Stress distribution of straight box girder cross—section of the span centre at top tray _^ = j垩 节点潢坐标/cm 图5直线箱梁跨中截面底板轴向正应力分布图 Fig.5 Fhe axial normal stress distribution of straight box girder CROSS—section of the span centre at the bottom 2.3剪力滞系数的计算 根据剪力滞系数的定义,首先按初等梁理论 计算出跨中截面顶板与底板中面的正应力值: MyT O'0T一丁一 一 _x :::一0.3.O37×10 ’ 02 6 MP;…. 第6卷第4期 徐飞鸿,等:薄壁箱粱剪力滞效应数值计算 19 MyB 一T一 一 5一0.497 2 MP..O37×10 ‘ ‘ 然后将其代人剪力滞系数 的计算公式(1), 就可以得出横截面各点处的剪力滞系数,直线箱 梁剪力滞系数沿跨中横截面的分布规律如图6 所示. 强 挺 R 京 节点横坐标/cm 图6 直线箱梁跨中截面剪力滞系数沿横截面变化图 Fig.6 The variation of shear lag coefficients along cross section of straight box girder cross—section of the span centre 3薄壁曲线箱梁的剪力滞效应分析 3.1 曲率半径的影响 为了研究曲率半径对箱梁剪力滞的影响,现 选取不同曲率半径的等截面简支曲线箱梁,横截 面尺寸及测点位置与直线箱梁相同.简支曲线箱 梁的有限元模型如图7所示. 图7 简支曲线箱梁有限元模型 Fig.7 Finite element model of freely supported curved box girder 通过有限元方法数值求解,计算出了不同曲 率半径的曲线箱梁桥在受跨中对称集中荷载作用 下的剪力滞效应.半径选用2,2O,50,100,200 m 的5种情况,跨径取S一0.8 m,弹性模量E一 3 000 MPa,泊松比 一O.385.5种不同半径曲线 桥及直桥的剪力滞系数计算结果汇总如表2及图 8所示. 表2跨中横截面剪力滞系数 Table 2 The shear lag coefficient in CROSS—section of the span centre 注: 表示内侧腹板与翼板相交位置处的剪力滞系数; 表 示外侧腹板与翼板相交位置处的剪力滞系数. 1 l l l l O O O O O 轻 4 3 2 l O 9 8 7 6 5 嶷 器 U 5 l0 15 20 25 30 35 节点横坐标,cm 图8剪力滞系数随曲率半径变化曲线 Fig.8 The shear lag coefficient curve varied with the curvature radius 由计算结果可知: 1)箱梁的曲率半径越小,同一位置处的剪力 滞系数越大.如:曲率半径均为2 m的曲线梁的 内与直线梁的 在跨中截面处相差11.89%,故 半径变化对剪力滞效应的影响较大. 2)由于曲线梁内侧曲率半径小,则 外< 内, 且两者在跨中截面处差别较大,最大差距为 14.6 . 3)翼板中心处剪力滞系数不随曲率半径的变 化而变化. 4)综合考虑曲率半径对 外及 内的影响,在 曲线梁与直线梁的剪力滞系数相差小于5 时进 行控制,可以得出,只有当半径R≥300 m时,才 可以忽略半径影响直接按照直线梁进行计算. 3.2其他因素的影响 计算分析中还考虑了箱梁高宽比 ]、腹板间 距、横隔板道数及截面形式等因素对曲线箱梁剪 2O 长沙理工大学学报(自然科学版) 2009年12月 力滞效应的影响.限于篇幅,这里只给出分析结 果:高宽比变化对剪力滞效应的影响不大;较大的 腹板间距可以减小剪力滞效应的影响;横隔板的 设置在一定范围内将改变剪力滞效应沿桥纵向的 YAO Ling—seng.Curved girder[M].Beijing:China Communications Press,1989:301-319. [3] 吴西伦.弯梁桥设计[M].北京:人民交通出版社. 1990:52-75. 分布情况;腹板数量增加可使截面的受力更加均 wu Xi—lun.Curved girder bridge design[M].Bei, 匀,降低剪力滞效应的影响. 4 结论 运用有限元方法对有机玻璃简支直线箱形梁 和曲线箱形梁模型进行了数值计算与分析,可以 得出以下结论. 1)采用数值计算方法对直线箱形梁的求解结 果与变分法理论计算值及试验数据吻合较好,验 证了本研究数值计算方法的正确性,为其他箱形 梁模型剪力滞效应的数值计算提供了理论指导. 2)曲率半径对连续曲线箱梁的剪力滞效应有 较大影响.与直线箱梁相比,截面相同位置处的剪 力滞系数随着曲率半径的减小而增大,增幅远超 过5%以上. 3)曲线箱梁应力沿横向分布不均匀.引起曲 线箱梁应力沿横向分布不均匀的原因有多种,除 了曲线箱梁桥翘曲和畸变的应力以外,还有一个 非常重要的原因就是曲率沿横向变化产生纵向弯 曲正应力分布不均.由于曲率沿横向变化,使箱梁 内侧的弯曲正应力要大于外侧的纵向正应力.在 曲线箱粱桥的实际设计中,由于横坡超高的原因, 很多桥梁均将内侧腹板高度减小,而外侧腹板高 度增大,这与箱梁桥的内侧应力大于外侧应力的 特点是相反的,这需要引起工程设计界的重视. 4)当曲率半径足够大时(如:R≥300 m),曲 线箱梁可近似按直线梁来计算剪力滞系数,不会 产生较大误差. [参考文献] [1] 李惠生,张罗溪.曲梁桥结构分析[M].北京:中国铁 道出版社,1992:12卜128. LI Hui—sheng,ZHANG Luo-xi.Curved girder bridge structure analysis[M].Beijing:China Railway Press. 1992:I2I—l28. [2] 姚玲森.曲线梁[M].北京:人民交通出版社,l989: 3O卜319. jing:China Communications Press,1990:52—75. [4] 韦成龙.曾庆元.薄壁曲线箱梁考虑翘曲、畸变和剪 滞效应的空间分析[J].土木工程学报,2000,33(6): 81-87. WEI Cheng—long.ZENG Qing—yuan.A new element for thin—walled curved box girder analysis including warping distortion and shear lag effects[J].China Civil Engineering Journal,2000,33(6):81—87. [5] 吴幼明.岳珠峰,吕震宙.薄壁曲线箱梁剪力滞计算 的有限段方法[J].物理学报,2009(6):4 002—4 009. Wu You—ming,Yue Zhu-feng,Lv Zhen-zhou.Finite segment method for calculating the shear lag in thin- walled curved box girders[J].Acta Physica Sinica, 2009(6):4 002—4 009. [6] 张峰,叶见曙 徐向锋.基于连续体退化组合壳单元 的连续箱梁非线性分析[J].公路交通科技.2007,24 (4):89—94. ZHANG Feng,YE Jian—shu.XU Xiang—fcng.Nonlin— ear analysis of prestressed concrete continuous box girders based on CB degenerated combined shell[J]. Journal of Highway and Transportation Research and Development,2007,24(4):89—94. [7] McManus P E.Nasir C A.Culver C G.Horizontally curved girders state—of-the—art[J].JS D ASCE。1969. 95(5):823-829. [8] 高岛春生.曲线梁桥[M].北京:中国建筑工业出版 社.1994:28—31. GAODAO Chun-sheng.Curved girder bridge[M]. Beijing:China Building Industry Press,l994:28—31. [9] 欧阳永金.大宽跨比连续钢箱梁桥的剪力滞效应研 究[J].世界桥梁,2009(1):29—32. OUYANG Yong-jin.Investigation of shear lag effect of continuous steel box girder bridge with great width—to-span ratio[J].World Bridges,2009(1):29— 32. [io] 尤文刚.杨海科,周家有.桥梁工程中的剪力滞效应 [J].北方交通.2008(6):100—101. YOU Wen—gang.YANG Hal—ke,ZH()U Jia—you. Shear lag effect in bridge engineering[J].Northern Communications,2008(6):i00-i01. (下转至第26页) 26 长沙理工大学学报(自然科学版) 2009年l2月 [5] 夏禾.徐幼麟.阎全胜.大跨度悬索桥在风与列车荷 journal of the China Railway Society。2004,26(3): 载同时作用下的动力响应分析[J].铁道学报,2002, 7l一75. 24(4):83-91. [93 C S Cai,S R Chen.Framework of vehicle-bridge-wind XIA He,XU You 1in.YAN Quan-sheng.Dynamic dynamic analysis[J].Journal of Wind Engineering response of longspan suspension bridge to high wind and Industrial Aerodynamics,2004(92):579-607. and running train[J].Journal of the China Railway [1O] Chen S R.Dynamic performance of bridges and ve— Society,2002,24(4):83—91. hicles under strong windED].Baton Rouge:Depart— [63 Xu Y L.Xia H。Yan Q S.Dynamic response of SHS— ment of Civil and Environmental Engineering of pension bridge to high winds and running train[J]. Louisiana State University,2004. Journal of Bridge Engineering.2003(8):46—55. [儿] 韩万水.风一汽车一桥梁系统空问耦合振动研究 [73 李永乐.风一车~桥系统非线性空间耦合振动研究 [D].上海:同济大学.2006. [D].成都:西南交通大学,2003. HAN Wan-shui.Three—dimensional coupling vibra— LI Yong—le.Nonlinear three—dimensional coupling vi— lion of wind—vehicle bridge system[D].Shanghai: bration of wind vehicle-bridge system[D].Chengdu: Tongji University.2006. Southwest Jiaotong University,2003. [12] 韩万水,陈艾荣.随机车流下的风一汽车一桥梁系 [8] 李永乐,廖海黎,强士中.车桥系统气动特性的节段 统空间耦合振动研究[J].土木工程学报.2008.41 模型风洞实验研究[J].铁道学报,2004.26(3):7I一 (9):97—102. 75. HAN Wan-shui。CHEN Ai—rong.Three—dimensional LI Yong—le,LIAO Hai—li,QIANG Shi—zhong.Study coupling vibration of wind—・vehicle—・bridge systems on aerodynamic characteristics of the vehicle—bridge under random traffic flow[J].China Civil Engineer— system by the section model wind tunnel test[J]. ing Journal,2008.4l(9):97—102. (上接第2O页) 应分析中的应用[j].甘肃科技纵横.2006(35):144一 In]廖志宏.ANSYS在箱梁剪力滞效应分析中的应用 l46. [J].河南建材.2009(3):68—69. SHEN Guo—shun.LIU Shi—zhong.The application of LIAO Zhi—hong.The application of ANSYS tO anal— ANSYS tO analysis of shear lag effect of thin-wall ysis of shear lag effect of box-girder[J].Henan box—girder[J].Scientific and Technical Information Building Materials,2009(3):68-69. of Gansu,2006(35):144—146. [12]张士铎,邓小华,王文州.箱形薄壁梁剪力滞效应 [143李长凤.杜文学.周莉.曲线箱梁剪力滞效应及其影 [M].北京:人民交通出版社.1998:l2l-122. 响因素rJ].黑龙江科技学院学报.2008(3):98—102. Z H ANG Shi-duo.DENG Xiae-hua.WANG Wen-zhou. LI Chang—feng.DU Wen-xue.ZHOU Li.Analysis of Shear lag effect of box girder with thin-wail[M].Bei— shear lag and effect factors in curved box girder[J]. jing:China Communications Press.1998:12l—l22. Journal of Heilongjiang Institute of Science and [13] 申国顺,刘世忠.ANSYS软件在薄壁箱梁剪力滞效 Technology,2008(3):98—102.