刘 超， 黄钰豪， 高 展

(同济大学 土木工程学院, 上海 200092)

波形钢腹板桥梁是一种新型的组合结构桥梁，波形钢质量轻、抗剪性能优良，利用波形钢腹板替代混凝土腹板[1]，可以降低工程造价，同时波形钢腹板的褶皱效应也提高了预应力的效率[2].目前中国很多大桥和特大桥都采用波形钢腹板形式，如南昌朝阳大桥、山西运宝大桥等.然而，由于桥梁结构形式较为复杂，腹板和顶、底板的受力特性尚未被研究透彻.国内外学者对波形钢腹板截面抗剪做了许多研究.在剪力分配方面，Abbas等[3]和Cafolla[4]通过试验发现，混凝土顶、底板主要承担截面弯矩，波形钢腹板主要承担截面剪力.李立峰等[5]认为，波形钢腹板承担剪力在80%以上.侯立超等[6]进行了单箱三室波形钢腹板预应力混凝土箱梁静力试验研究，得出了单箱三室波形钢腹板剪力分配特点.周绪红等[7]和聂建国等[8]通过试验对比，分别给出了考虑局部屈曲强度的半经验半理论式.Driver等[9]通过波形钢腹板I形梁的局部屈曲和整体屈曲试验数据分析发现，基于板壳屈曲的理论高估了波形钢腹板梁桥的抗剪承载力，并提出了一种低边界的抗剪承载力计算式.在不考虑腹板屈曲的情况下，苏俭等[10]推导了变截面波形钢腹板竖向剪应力和水平剪力的实用计算方法，并以南京滁河大桥为例，进行了有限元模型计算分析，验证了计算方法的正确性.现有的计算方法主要有以下三种：① 将截面总剪力除以波形钢腹板面积得到剪应力[11-12]，中国和日本的规范均采用这种方法进行计算，但该方法忽略了顶、底板的承剪比例，认为剪力仅由波形钢腹板承担，对于钢腹板来说偏保守，对于混凝土顶、底板来说偏不安全；② 利用传统的计算截面静矩的方法来计算剪应力[13]，该方法仅仅考虑由剪力产生的剪应力，而忽略了弯矩和轴力产生的剪应力，对于等截面桥梁来说是合适的，但对于变截面桥梁来说准确性不足；③ 在方法②的基础上，同时考虑弯矩和轴力产生的剪应力[10,14]，该方法计算结果较为准确，但目前尚未推广至单箱多室和变宽桥梁的计算.随着波形钢腹板桥梁的跨径越做越大，单箱多室与变宽桥梁越来越常见.因此，针对单箱多室波形钢腹板变宽变高截面剪应力计算理论的空白，以及目前规范中计算方法的不足，提出了一种变宽变高波形钢腹板截面剪应力计算理论，认为剪力、弯矩和轴力均会产生剪应力，并且后两者仅在变截面时产生剪应力.结合有限元分析，给出了顶、底板和腹板的剪力分配比例，为波形钢腹板桥梁抗剪设计提供借鉴.

1 剪应力计算方法

1.1 计算基本假定

根据波形钢腹板截面的受力特点，给出如下计算基本假定：① 在计算的整个过程中，变形前后截面始终保持平面或近似平面，平截面假定始终成立；② 钢和混凝土为各向同性材料，均处于线弹性受力阶段，在计算过程中结构始终处于弹性阶段，尚未进入塑形，所有变形均为小变形；③ 波形钢腹板不承受轴力，轴力全由混凝土顶、底板承担，而且由于波形钢腹板厚度较小，不考虑弯矩对其造成的轴向应力；④ 不考虑波形钢腹板和顶、底板的轴向滑移； ⑤ 计算时不考虑畸变和扭转影响.

1.2 波形钢腹板变截面计算理论

图1a为一般波形钢腹板箱梁截面图.箱梁腹板数为n，上顶板宽度为b1，厚度为t1，下顶板宽度为b2，厚度为t2，箱梁截面高度为h，形心距上顶板距离为yc，任意高度处截面距上顶板距离为y0，坐标轴原点取在上顶板左上角，向右为z轴正方向，向下为y轴正方向.图1b为距上顶板距离y0处截取的箱梁截面微元图，沿桥轴向长度为dx，左截面受到的轴向合力为F，右截面受到的轴向合力为F+dF，y0处断面存在的剪应力为τ，由于微元长度很小，因此可以认为在dx长度内，剪应力大小不变，均为τ.化简轴向平衡方程后有

(1)

式中:b(y)为y高度处截面总宽度.

根据剪应力互等定理，式(1)所求剪应力τ与截面剪应力相等，所以由式(1)所求得的τ即可代表截面剪应力.在图1b中，轴向合力F可以看成是轴向力FN和由弯矩产生的弯曲轴向力FM的线性叠加，即：

F=FN+FM

(2)

(3)

(4)

式中：N为截面上轴力；M为截面上弯矩；σyN为轴力在截面上产生的正应力；σyM为弯矩在截面上产生的正应力；Ay0为坐标0～y0截面总面积；A为截面总面积；I为截面惯性矩.

a 波形钢腹板箱梁截面

b 截面微元区段

对于式(4)，可作如下变化：

(5)

式中：y*为坐标0～y0截面形心坐标；Sy0为坐标0～y0截面的静矩.

将式(3)～(5)代入式(2)，可得

(6)

截取全截面轴向的一段小微元，如图2所示.图2中，q为外荷载，α为轴线与水平线的夹角，Q为剪力.对右截面轴线处取距，由力矩平衡可知

(7)

(8)

将式(6)代入式(1)，结合式(8)，有

图2 轴向微元区段

(9)

(10)

(11)

式(11)即为材料力学中截面剪应力的传统计算式.在变截面波形钢腹板箱梁截面中，剪应力由三个部分组成，分别是由弯矩造成的剪应力τM、由轴力造成的剪应力τN和由剪力造成的剪应力τQ，前两者仅存在于变截面剪应力计算中.τ可表示为

τ=τM+τQ+τN

(12)

A=b1t1+b2t2

(13)

(14)

(15)

将式(13)对x求偏导，可得

(16)

将式(14)对x求偏导，可得

(17)

(18)

(1)y0处于上翼缘板内(0

(19)

(2)y0处于腹板内(t1≤y0≤h-t2)

(20)

(3)y0处于下翼缘板内(h-t2

Ay0=b1t1+b2(y0+t2-h)，

(21)

将式(13)～(21)代入式(10)即可求得截面上任意位置的剪应力大小.

2 有限元模型

为了验证上述计算模型的正确性，本文利用Abaqus软件建立了悬臂梁模型并进行有限元仿真.梁长20 m，截面为变宽变高波形钢腹板截面，有限元模型如图3所示.自由端截面上顶板长6.00 m，宽0.25 m，下顶板长4.00 m，宽0.25 m，为了与工程相仿，腹板厚度取0.02 m，截面总高度取2.75 m；固定端截面上顶板长8.0 m，宽0.4 m，下顶板长6.0 m，宽0.5 m，腹板厚度0.02 m，截面总高度4.0 m，中间截面各参数由自由端向固定端线性变化.截面示意图如图4a所示，波形钢腹板示意图如图4b所示.在悬臂梁自由端施加竖向集中力1 000 kN，同时在自由端截面顶、底板施加轴向均布荷载，合力为1 000 kN.在悬臂梁取自由端靠里1 m位置为截面A(为排除集中力附近应力集中干扰)，中跨为截面B，固定端为截面C，作为计算截面.

图3 有限元模型

a 箱梁截面

b 波形钢腹板

模型中采用的材料为钢和混凝土，使用线性应力-应变关系.由前述假定认为波形钢腹板与混凝土顶、底板间不存在相对滑移，故在有限元模型里，将接触面间所有自由度耦合，并遵从共同变形原则.同时，定义混凝土和钢为各向同性材料.模型网格划分采用扫略划分，为了保证精度，波形钢腹板和混凝土顶、底板均采用空间六面体C3D8网格，共计345 523个节点和210 279个单元.各参数设置如表1所示.

表1 悬臂梁参数设置

3 结果分析

3.1 有限元正确性验证

提取有限元模型(FEM)中梁的内力和挠度，与解析解进行对比.通过解析的方法进行模型挠度的计算.根据欧拉梁理论，弯矩产生的挠度计算式为

(22)

考虑变截面情况，分母中惯性矩为x的函数，沿梁轴向进行积分，分别计算出距固定端1 m，2 m，3 m，…，20 m间距处梁挠度.

当杆件在两相邻的横截面处有一对垂直于杆轴，但方向相反的横向力作用时，杆件发生的变形为该两截面沿横向力方向发生的相对错动，此变形称为剪切变形.波形钢腹板在受力过程中发生明显的剪切变形，从而引起波形钢腹板梁变形明显增大.进一步增加剪切变形产生的挠度计算，根据铁木辛柯梁理论，剪切应变能

(23)

式中：G为剪切模量.

由于目前尚未有完善的变截面变宽波形钢腹板剪应力计算理论，τ按照本文所提出理论计算，对剪切应变能求偏微分可得剪力产生挠度，如下所示：

(24)

将有限元挠度计算结果和解析解计算结果列于图5中.悬臂端有限元挠度计算结果为6.60 mm，解析解挠度计算结果为6.27 mm，解析解与有限元解误差为5%，证明了有限元模型的正确性.

图5 各种挠度计算方法对比

3.2 整体应力分析

悬臂梁剪应力云图如图6所示.从图6可以看出：① 悬臂梁在自由端竖向集中力和轴向均布力作用下，剪应力从自由端向固定端递减，一个原因是截面由自由端向固定端逐渐变高、变宽，从而使截面的惯性矩I增大，导致剪应力减小，另一个原因是弯矩和轴力产生的剪应力对总剪应力有影响；② 波形钢腹板褶皱处剪应力较小，平直处剪应力较大，轴向应力分布不仅呈现出从自由端向固定端逐渐减小的趋势，在局部还呈现出大—小—大的分布；③ 顶、底板的剪应力远远小于腹板的剪应力，主要原因是顶、底板的b过大，是腹板的上百倍，从而导致很小的剪应力；④ 波形钢腹板剪力分配不均匀，中腹板的剪应力小于边腹板,不考虑连接处局部效应，中腹板最大剪应力为6.14 MPa，边腹板最大剪应力为7.70 MPa，相差约20%.

图6 悬臂梁剪应力云图

3.3 截面剪应力分析

对于该悬臂梁，给出本文提出的变截面计算理论(VCM)(见式(9))、传统计算理论(TCM)(见式(11))和有限元(FEM)方法的计算结果用以精度对比.三个截面的计算结果如图7所示，其中腹板截面应力取边腹板的应力.

a 截面A

b 截面B

c 截面C

图7 截面剪应力分布

Fig.7 Sectional shear stress distribution

从图7可以看出：① 从截面A到截面B再到截面C，腹板和顶、底板的剪应力逐渐衰减；② 顶、底板剪应力相比腹板来说很小；③ 本文提出的变截面计算理论和有限元模型计算结果误差很小，最大误差不足5%,而传统计算理论由于忽略弯矩和轴力产生的剪应力，误差较大，误差从自由端到固定端递增，最大可达200%.由于波形钢腹板的特殊性，从数值模拟结果来看，传统计算理论由于误差较大，已不再适用于计算，本文提出的变截面计算理论误差很小，对变宽变高波形钢腹板截面适用性很强.

3.4 顶、底板剪力分配

目前，大部分设计均不考虑顶、底板的抗剪能力，认为所有剪力由波形钢腹板承受，但实际上由于混凝土顶、底板具有部分承剪能力，这样的设计对波形钢腹板来说过于保守，对于混凝土顶、底板来说偏不安全.对悬臂梁截面A、B和C的腹板应力进行积分，可得腹板承受的剪力大小，再除以总剪力Q即可得腹板承剪比例，如下所示：

(20)

式中：ρw为腹板承剪比例；τw为腹板剪应力.顶、底板承剪比例为(1-ρw).截面A、B和C各部分承剪比例如图8所示.

离自由端越远，悬臂梁受到的弯矩越大，弯矩对顶、底板和腹板剪应力产生的影响越大.对于有限元计算结果，波形钢腹板承剪比例在自由端达到了近87.1%，后逐渐减小，在固定端达到谷值42.3%，而顶、底板在自由端承剪比例为12.9%，后逐渐增大，在固定端达到峰值57.7%.这说明传统假定仅仅在无弯矩或弯矩很小区域适用，在弯矩较大区域，由于顶、底板承剪比例很大，因此忽视顶、底板承剪是不合理的.有限元模型计算结果和本文提出的变截面计算理论的计算结果很相近，说明了本文计算理论的有效性.

图8 剪力分配比例

3.5 剪力、轴力、弯矩剪应力分布

根据前述计算式，变截面波形钢腹板梁某一截面上钢腹板的剪应力可以看成由以下三个部分组成：弯矩剪应力τM，轴力剪应力τN，剪力剪应力τQ.利用式(20)，计算出悬臂梁截面A、B和C的钢腹板剪应力分配比例，如图9所示.

图9 剪应力分配比例

从图9可以看出：① 由弯矩和轴力产生的剪应力方向和剪力产生的剪应力方向相反，而截面总剪应力为三部分的叠加，所以总剪应力数值小于剪力产生的剪应力；② 从自由端到固定端，剪力剪应力逐渐减小，弯矩剪应力逐渐增大；③ 从自由端到固定端，弯矩和轴力产生的剪应力与剪力剪应力的比值从0.096 0增大到0.615 6，说明传统计算方法只计算剪力剪应力是不正确的.

4 结论

(1)提出了一种针对变宽变高波形钢腹板桥梁的剪应力计算方法，认为剪力、弯矩和轴力均会产生剪应力，并且后两者仅仅在变截面时产生剪应力.本文方法计算结果与有限元模型计算结果吻合程度高，可为波形钢腹板桥梁的抗剪设计提供依据.

(2)针对变截面波形钢腹板桥梁，传统剪应力计算理论适用性较差.在腹板断面上，弯矩和轴力产生的剪应力方向与剪力剪应力方向相反，并且在弯矩较大时，弯矩和轴力产生的剪应力与剪力剪应力的比值可达0.6以上，从而导致实际剪应力数值小于剪力剪应力，这也是传统剪应力计算理论误差较大的原因.

(3)波形钢腹板截面顶板剪应力很小，底板剪应力相对顶板而言较大，主要来源于底板受到的压应力的竖向分力.在无弯矩或弯矩很小区域认为仅由腹板承剪是可以接受的，但当弯矩较大时，顶、底板承剪比例能达到50%以上，因此应考虑顶、底板的承剪能力.传统计算理论、计算模型和规范认为剪力全部由钢腹板承担而不考虑混凝土顶、底板作用，容易造成波形钢腹板浪费，顶、底混凝土板不足的情况.

(4)对单箱双室波形钢腹板桥梁进行有限元分析，发现腹板间存在剪应力分配，中腹板剪应力较小，边腹板剪应力较大.本文提出的计算理论对边腹板吻合程度高，在一定程度上提高了结构安全性.