移动荷载作用下曲线箱梁剪力滞效应

卢海林，南子俊，张伟

武汉工程大学资源与土木工程学院，湖北 武汉 430074

移动荷载作用下曲线箱梁剪力滞效应

卢海林，南子俊，张伟

武汉工程大学资源与土木工程学院，湖北 武汉 430074

采用大型通用有限元软件ANSYS分析了移动荷载作用下曲线箱梁的剪力滞效应.分别设置了3种不同大小的移动荷载和3个不同的作用位置，研究了荷载移动到曲线箱梁L／2跨时该截面的剪力滞效应分布情况.结果表明：荷载沿箱梁中心线移动时，截面内外侧剪力滞效应分布不均匀；荷载沿外侧腹板移动时，外侧腹板处正剪力滞现象明显，内侧腹板以及翼缘板上负剪力滞现象明显；荷载沿内侧腹板移动时，内侧腹板处正剪力滞现象明显，外侧腹板以及翼缘板上负剪力滞现象明显.移动荷载大小的改变，对曲线箱梁剪力滞效应的分布情况影响较小.但随着移动荷载的加大，局部剪力滞效应会少量地增大.

曲线箱梁；剪力滞；移动荷载；有限元

0 引言

在公路立交及城市道路和桥梁工程中，曲线箱梁桥是实现各方向交通连接的必要手段.曲线箱梁桥正常作用时箱梁中存在着剪力滞效应.对曲线箱梁剪力滞效应的研究，一直以来都是桥梁专家所关心的问题.已有的研究成果多数是静载作用下的剪力滞效应，包括集中荷载和均布荷载作用下曲线箱梁剪力滞效应［1-2］.但实际工程中，桥梁结构绝大部分是受动荷载作用的.因此，开展动荷载作用下箱梁剪力滞效应的研究是十分有必要的.陈鸿鸣等［3］根据箱梁剪力滞过去的研究，总结了当前箱梁剪力研究的主要方向，包括几何和材料的非线性效应、荷载作用形式等.刘健新等［4］用变分原理推导了考虑剪力滞效应的箱梁强迫振动时的微分方程、边界条件，建立了方程解的差分格式，得出了薄壁箱梁强迫振动剪力滞效应分析的理论方法.徐勋［5］等利用有限元分析了扁平曲线箱梁考虑剪力滞效应和剪切变形影响的静动力特性，特别地对移动荷载作用下结构的弯扭耦合振动特性进行了研究.长沙理工大学的蔡汶珊［6］建立了曲线箱梁的动力分析模型，验证了一种双脊骨有限元模型的可靠性，并分析了地震作用下曲线箱梁的剪力滞效应.武汉工程大学的卢海林［7］利用大型通用有限元软件ANSYS建立了悬臂曲线箱梁模型，并分析了移动荷载速度对曲线箱梁剪力滞效应的影响.Hugo C等［8］通过大量的现场检测，研究了在行车荷载作用下某曲线箱梁桥的自振特性.本文利用有限元软件ANSYS，选取不同大小的移动荷载在不同位置移动到跨中截面时，分析跨中截面上剪力滞分布规律.

1 模型的建立

1.1 模型尺寸

模型尺寸来源于文献［1］，为一采用以有机玻璃加工而成的曲线箱梁模型.模型截面形状和尺寸如图1所示，模型空间形态示意如图2所示，图中中心线全长L＝1．047 m，曲率半径r＝2 m，圆心角角度θ＝30°.有机玻璃板厚为8 mm，上、下顶板宽分别为400 mm和192 mm，腹板高h＝100 mm.有机玻璃材料弹性模量为E＝3000MPa，泊松比μ＝0．385，密度dens＝1 180 kg／m3.

图1 箱梁横截面尺寸（mm）Fig.1 Dimensions of box girder cross－section（mm）

图2 箱梁空间示意图Fig.2 Space diagram of box girder

1.2 选择荷载

进入ANSYS瞬态分析，在节点上施加瞬间作用以模拟移动荷载，荷载步采用阶跃荷载［9］，使移动荷载中心由箱梁一端沿箱梁中心线匀速移动至另一端，完成加载过程，并进行数据采集分析.

根据文献［7］的荷载总量使用，为了使模型上剪力滞效果更加明显，分别采用总量为G1＝200 N、G2＝400 N、G3＝600 N均匀分布在4个集中力（F1、F2、F3、F4）上的移动荷载以模拟实际车辆荷载，分前后两排间距10 cm沿着顶板与腹板交接处以1 m／s的速度移动如图3所示.随后将3组移动荷载分别向内外两侧横向各偏移一次，偏移距离为10 cm（在有限元软件中大约为5个划分网格的距离），沿着外侧和内侧的腹板移动如图4所示，共得到移动荷载分别偏向内侧、沿箱梁中心线和偏向外侧移动时的九组数据，具体工况设计如表1所示.

表1 移动荷载加载工况Table1 Moving load work condition

图3 荷载中心沿箱梁中心线移动Fig.3 Load center moving along box girder midline

图4 荷载中心分别沿外侧和内侧腹板移动Fig.4 Load center moving along both inside and outside web

1.3 有限元模型

箱梁模型有限元模型如图5所示.由于梁单元不能完全满足计算精度的需求，而采用三位实体单元建模时其节点数较多较为繁杂，综合来看壳单元最为符合［10］.因此，该箱梁有限元模型采用单元类型为四节点shell63单元，模型结构一共划分为1 664个单元，1 696个节点.约束条件采用文献［10］中的简支约束条件.

图5 箱梁有限元模型Fig.5 Finite element model of box girder

2 数值结果及分析

2.1 移动荷载作用下箱梁跨中截面顶板和底板剪力滞系数

当荷载移动完毕时，进入后处理，若记整个运动时间为T、荷载步为N，将荷载加载时间调至T／2或将荷载步调至N／2时，移动荷载重心刚好经过跨中L／2截面处，此时该截面上的作用影响效果最大［11］，因此选择跨中截面为研究对象.为了研究其剪力滞效应，取出该截面上的顶板与底板节点号，读取应力数据换算为剪力滞系数后，进行整理分析，分别得出各个工况下截面剪力滞系数曲线分布比较如图6～图17所示（数据横轴正数为曲线箱梁横截面内侧），并列出顶板主要节点处的剪力滞系数值如表2所示.

图6 G＝200 N时箱梁顶板L／2截面剪力滞系数Fig.6 The shear lag coefficient of mid－section of box girder at top when G＝200 N

图7 G＝200 N时箱梁底板L／2截面剪力滞系数Fig.7 The shear lag coefficient of mid－section of box girder at bottom when G＝200 N

图8 G＝400 N时箱梁顶板L／2截面剪力滞系数Fig.8 The shear lag coefficient of mid－section of box girder at top when G＝400 N

图9 G＝400 N时箱梁底板L／2截面剪力滞系数Fig.9 The shear lag coefficient of mid－section of box girder at bottom when G＝400 N

图10 G＝600 N时箱梁顶板L／2截面剪力滞系数Fig.10 The shear lag coefficient of mid－section of box girder at top when G＝600 N

图11 G＝600 N时箱梁底板L／2截面剪力滞系数Fig.11 The shear lag coefficient of mid－section of box girder at bottom when G＝600 N

图12 不同荷载沿外侧腹板移动时顶板L／2截面剪力滞系数Fig.12 The shear lag coefficient of mid－section of box girder at top under the different load moving along outside web

图13 不同荷载沿外侧腹板移动时底板L／2截面剪力滞系数Fig.13 The shear lag coefficient of mid－section of box girder at bottom under the different load moving along outside web

图14 不同荷载沿中心线移动时顶板L／2截面剪力滞系数Fig.14 The shear lag coefficient of mid－section of box girder at top under the different load moving along center line

图15 不同荷载沿中心线移动时顶板L／2截面剪力滞系数Fig.15 The shear lag coefficient of mid－section of box girder at bottom under the different load moving along center line

图16 不同荷载沿内侧腹板移动时顶板L／2截面剪力滞系数Fig.16 The shear lag coefficient of mid－section of box girder at top under the different load moving along inside web

图17 不同荷载沿内侧腹板移动时底板L／2截面剪力滞系数Fig.17 The shear lag coefficient of mid－section of box girder at bottom under the different load moving along inside web

2.2 数据分析

从图6～图17中很容易看出，由于该模型是曲线箱梁，所以数据呈现明显的非对称性，荷载作用在哪一侧，哪一侧的剪力滞效应就会越明显，并且从图表中可以归纳得出以下结论：

a．无论荷载大小多少，荷载作用在3种不同的位置时曲线箱梁顶板剪力滞效应曲线的分布规律近似，均在腹板与顶板的交接处达到一个极大值，呈现近似M的形状.

b．从图6、图8和图10中看出，当荷载沿箱梁中心以及外侧腹板移动时，外侧顶板上出现正剪力滞现象，内侧顶板上出现负剪力滞现象，荷载沿外侧腹板移动时剪力滞效应更为明显；当荷载沿内侧腹板移动时，内侧顶板上出现正剪力滞现象，外侧顶板上出现负剪力滞现象.

c．从图12和图14中看出，当荷载沿箱梁中线以及外侧腹板移动时，外侧翼缘的端部随着荷载的加大，剪力滞系数呈减小趋势，在内侧翼缘端部均非常小；从图16中看出，荷载沿内侧腹板移动时，内外侧翼缘端部剪力滞现象同外侧.

d．从表2中的各工况下重要节点剪力滞系数值可以看出，剪力滞效应与荷载大小没有特定的单调关系，各处数值近似，表明在模型上任意位置，荷载大小对其剪力滞效应的数值影响程度不大.但当荷载加大时，荷载偏心作用情况下，剪力滞系数分化较大，在腹板与顶板交接处越来越大，而在翼缘端部越来越小.

e．荷载作用在3种不同位置对曲线箱梁底板剪力滞分布有明显影响：从图13中看出，当荷载沿外侧腹板移动时，底板上剪力滞规律由外到内越来越小；图15中，当荷载沿箱梁中线移动时，底板上剪力滞规律呈现开口向上的抛物线趋势；图17中，当荷载沿内侧腹板移动时，底板上剪力滞规律由外到内越来越大.其剪力滞数值亦同顶板，与荷载大小无明显联系，基本近似.

3 结语

本文主要采用有限元软件ANSYS进行研究分析了3种不同大小的移动荷载在结构上改变横向作用位置时，移动到曲线箱梁跨中L／2处时该截面的剪力滞效应分布规律.分析结果表明：当荷载沿箱梁中线移动时，内外侧剪力滞效应分布不均.由于偏心距较大，使得偏向外侧的顶板剪力滞效应更为明显.荷载沿两侧腹板作用时，截面上的剪力滞效应分化较为严重.荷载作用的一侧反应明显，在靠近腹板的顶板上呈现较大正剪力滞现象，而在另一侧腹板以及翼缘板上呈现负剪力滞现象.此外，不同大小的荷载作用在结构上任意位置时，随着荷载的加大，使得顶板上受力在顶板与腹板处较为集中，呈现局部应力集中现象，使得局部剪力滞效应增大，但对整体剪力滞规律的变化影响较小.

表2 各个工况下顶板主要节点剪力滞系数值Table 2 Shear lag coefficient of main roof node in each work condition

致谢

本研究得到国家自然科学基金委员会和武汉工程大学的资助，在此表示感谢！

［1］卢海林．预应力混凝土曲线箱梁剪力滞效应理论分析与试验分析［D］．天津：天津大学，2005．

LU Hai-lin．Theoretical analysis and experimental study on shear－lag effect of prestress concrete curved box beam［D］．Tianjin：Tianjin Uni versity，2005．（in Chinese）

［2］FABRIZIO Gara，GRAZIANO Leoni，LUIGINO Dezi．A beam finite element including shear lag effect for the time-dependent analysis of steel-concrete composite decks［J］．Engineering Structures，2009，31：1888－1902．

［3］陈鸿鸣，乔静宇．混凝土箱梁剪力滞研究现状与发展［J］．结构工程师，2011，27（1）：161－166．

CHEN Hong－ming，QIAO Jing－yu．State of the art about shear lag effect of concrete box girder bridges［J］．Structural Engineering，2011，27（1）：161－166．（in Chinese）

［4］刘健新，马麟，胡庆安．薄壁箱梁振动时的剪力滞效应［J］．郑州大学学报，2008，29（3）：122－128．

LIU Jian－xin，MA Lin，HU Qing－an．The research of analysis method of shear lag＇s effect on forced vibration of thin－wall box－girder［J］．Journal of Zhengzhou University，2008，29（3）：122－128．（in Chinese）

［5］徐勋，卫星，刘德军，等．扁平曲线箱梁静动力分析［J］．公路交通科技，2007，24（12）：60－65．

XU Xun，WEI Xing，LIU De－jun，et al．Analysis for static and dynamic characteristics of flat curved box beam［J］．Journal of Highway and Transportation Research and Development，2007，24（12）：60－65．（in Chinese）

［6］蔡汶珊．薄壁曲线箱梁剪力滞数值计算与实验研究［D］．长沙：长沙理工大学，2010．

CAI Wen－shan．Numerical calculate and experimental research about shear lag effect of curved thin wall box girders［D］．Changsha：Changsha University of Science＆Technology，2010．（in Chinese）

［7］卢海林，张伟，颜昌雄．移动荷载速度对悬臂曲线箱梁剪力滞效应的影响［J］．武汉工程大学学报，2014，36（3）：18－21．

LU Hai－lin，ZHANG Wei，YAN Chang－xiong．Moving load speed on shear lag effect of cantilever curved box girder［J］．Journal of Wuhan Institute of Technology，2014，36（3）：18－21．（in Chinese）

［8］HUGO C Gomeza，PAUL J Fanning，MARIA Q Fenga，et al．Testing and long-term monitoring of a curved concrete box girder bridge［J］．Engineering Structure，2011（33）：2861－2869．

［9］王新敏．ANSYS工程结构数值分析［M］．北京：人民交通出版社，2007．

WANG Xin－min．ANSYS numerical analysis of engineering structural［M］．Beijing：China Communications Press，2007．（in Chinese）

［10］李洪．曲线连续箱梁桥梁剪力滞效应分析及实用计算［D］．北京：北京交通大学，2010．

LI Hong．Analysis of the shear lag effect in curved continuous box－girder bridge and the research on practical calculating method［D］．Beijing：Beijing Jiaotong University，2010．（in Chinese）

［11］赵云安，王炎．车辆荷载作用下连续箱梁桥剪力滞效应分析［J］．铁道建筑，2013（2）：9－11．ZHAO Yun－an，WANG Yan．Analysis of the shear lag effect in curved continuous box－girder bridge under the vehicle load［J］．Railway Engineering，2013，（2）：9－11．（in Chinese）

Shear lag effect of curved box－girder under moving load

LU Hai－lin，NAN Zi－jun，ZHANG Wei
School of Resource and Civil Engineering，Wuhan Institute of Technology，Wuhan 430074，China

Large general－purpose software ANSYS was adopted to analyze the shear lag effect of curved boxgirder under the moving load．Three different sizes of moving loads and three different positions were set to investigate the distribution of shear lag effect of the cross－section when the loads moved to the mid－span of curved box－girder．The results show that the distribution of shear lag effect of both inside and outside section is uneven when the loads move along the box－girder midline；the phenomenon of positive shear lag is obvious at the outside web，the phenomenon of negative shear lag is obvious at the inside web and flange plate when the loads move along the outside web；the phenomenon of positive shear lag is obvious at the inside web，the phenomenon of negative shear lag is obvious at the web and flange plate when the loads move along the inside web．The change of the size of moving load has w influences on the distribution of shear lag effect of curved box－girder．But the partial shear lag effect increases slightly with the increase of moving load．

curved box girder；shear lag；moving load；finite element

U448.42

A

10.3969/j.issn.1674-2869.2015.05.005

1674-2869（2015）05-0023-05

本文编辑：龚晓宁

2015－03－23

国家自然科学基金（51378404）；武汉工程大学研究生教育创新基金（CX2013115）

卢海林（1965－），男，湖南安乡人，教授，博士.研究方向：桥梁结构分析、建筑结构检测与加固.