1. 引言
波形钢腹板PC组合箱梁桥结构作为一种新型的钢–混凝土组合结构,其顶底板为钢筋混凝土结构,腹板为波形钢板,预应力采用体外或体内与体外相结合的结构形式。该结构充分利用了混凝土抗压强度高,波形钢腹板抗剪屈服强度高的优点,有效的将钢、混凝土两种材料进行了融合,提高了材料利用率,是一种较为经济、合理、高效的结构形式 [1] [2] [3] [4] 。
与传统的预应力混凝土箱梁相比,它无论是在结构形式还是受力性能方面均有诸多特点,结构方面,波形钢腹板代替混凝土腹板使自重减轻,并且由于腹板特有的折叠形状,剪切刚度要在钢板剪切模量的基础上进行折减;受力性能方面,波形钢腹板是组合箱梁的剪力的主要承担者,混凝土顶板和底板主要起到抗弯的作用,充分发挥了材料特性 [5] [6] [7] 。同时,由于混凝土顶底板又对截面抗剪有贡献,波纹腹板也对抗弯也有一定贡献,忽略它们之间的作用,虽然使得计算的截面抗力偏于安全和保守。但是在桥梁施工阶段中,对波纹钢腹板的变形的控制和计算精度要求比较高,如果计算的变形偏于保守,会使得计算结果与实际结构发生的变形产生较大的误差,给波纹钢腹板梁的安装和变形控制带来困难 [8] [9] [10] [11] 。
国内外学者对波形钢腹板组合箱梁的变形计算进行了较多研究。其中,Ezzeldin Yazeed [12] 提出基于弹性弯曲理论的推导公式,通过假定位移场,建立内力平衡方程、变形协调条件和物理方程,依据边界条件及荷载条件求出解析解;但实际中梁的边界条件、荷载状况较复杂,难以获得解析解。刘保东 [13] 等根据波形钢腹板特点,基于变分原理,结合微分方程及相关的边界条件,利用欧拉公式推导了考虑箱梁剪力滞和剪切变形的波形钢腹板变形计算公式,结果表明剪力滞对计算结果的影响较小,变形值仅增大了1.3%。李宏江、叶见曙等 [14] 基于Timoshenko梁理论,推导出了等截面波形钢腹板组合箱梁的变形计算公式,其中,混凝土顶底板仅起到抗弯作用,腹板仅起到抗剪作用,简单有效,但忽略了混凝土顶底板的抗剪能力,其计算结果偏大。在此基础上,国内外学者采用共轭梁法和等效刚度法等近似等效方法进行了改进。上述方法应用于变截面波形钢腹板箱梁时都有一定的难度。
本文基于辛普森积分公式,推导了箱梁的等效惯性矩、顶底板的等效厚度和钢腹板的等效高度的计算公式,建立了箱梁在集中荷载和自重荷载作用下的弯曲变形和剪切变形的计算公式,并且利用等截面剪应力计算公式得到了腹板承担的剪力比例。通过单箱单室变截面波形钢腹板组合箱梁算例,利用等效惯性矩法分别计算了其在集中荷载和自重荷载作用下的弯曲变形和剪切变形,并与梁单元模型和实体–壳单元模型进行了对比,结果显示集中荷载作用下等效惯性矩法计算的弯曲变形和剪切变形误差分别为1.3%和13.4%,梁单元模型计算的误差分别为4.6%和34.1%;自重荷载作用下的等效惯性矩法计算的弯曲变形和剪切变形误差分别为2.8%和15.8%,梁单元模型计算误差分别为5.5%和29.7%。最后,通过顶底板混凝土对剪切变形影响的定量分析,表明由于考虑了顶底板混凝土的抗剪能力,利用等效惯性矩方法计算的剪切变形,更接近实体–壳单元模型计算结果。本文建立的基于辛普森公式的波形钢腹板组合箱梁变形计算方法计算简便,计算精度较高,且避免了复杂积分计算和模型建立,计算效率较高。适合波施工过程中的变形分析计算。
2. 波形钢腹板组合箱梁的截面特性
2.1. 波形钢腹板轴向受力特性
波形钢腹板组合箱梁截面,是在传统的混凝土箱梁截面基础上用波形钢腹板替代混凝土腹板的组合箱梁,由于其“手风琴”构造特征如图1所示,轴力作用下,轴向变形大,故其表观弹性模量小 [15] 。波形钢腹板纵向表观弹性模量与钢板的弹性模量E0、波高h、板厚t以及波纹形状系数r有关,可按下式计算:
(1)
计算所得表观弹性模量仅为钢板模量的几百分之一,因此设计上一般忽略波形钢腹板的轴向受力问题,认为轴向力仅由混凝土顶底板承担。
2.2. 波形钢腹板组合箱梁抗弯等效截面
由于波形钢腹板纵向表观弹性模量和轴向刚度很小,可忽略腹板对截面的抗弯作用。因此,计算箱梁抗弯惯性矩时,仅考虑混凝土顶底板的作用,其等效截面如图2所示。
波形钢腹板箱梁受弯时,若不考虑波形钢腹板的作用,则混凝土顶底板的纵向正应力分布基本符合线性的特征,因此计算波形钢腹板组合箱型梁变形时采用这一“拟平截面假定” [16] [17] 。
2.3. 波形钢腹板有效剪切模量
Johnson RP [18] 等研究发现,剪力作用下波形钢板的剪切模量有所降低,提出了有效剪切模量,其计算公式为:
(2)
式中:a、b、d为波形钢腹板尺寸,其意义如图1所示;G为钢板的抗剪弹性模量,
,其中m泊松比 [19] [20] 。
3. 变截面波形钢腹板箱梁变形计算的等效惯性矩方法
与等截面梁不同,波形钢腹板变截面梁的惯性矩沿长度方向变化,等截面变形曲线公式不再适用,本文基于辛普森公式,建立适用于变截面梁的等效惯性矩方法,在此基础上分别研究集中荷载和自重荷载作用下的悬臂梁弯曲和剪切变形的等效计算方法。
3.1. 变截面波形钢腹板弯曲变形计算
3.1.1. 等效惯性矩的推导
根据材料力学的卡氏第二定理,集中荷载作用下悬臂梁自由端的端部位移的计算公式为:
(3)
式中:Ice为等效梁的等效惯性矩。
为求出等效惯性矩Ice,采用辛普森积分公式计算最大悬臂的变形值:
![](//html.hanspub.org/file/14-2750494x13_hanspub.png)
Figure 1. Cross section of corrugated steel web
图1. 波形钢腹板截面
![](//html.hanspub.org/file/14-2750494x15_hanspub.png)
Figure 2. Equivalent cross section diagram
图2. 等效截面示意图
![](//html.hanspub.org/file/14-2750494x16_hanspub.png)
Figure 3. Calculation of simple diagram by equivalent moment of inertia
图3. 等效惯性矩法计算简图
(4)
式中:
为被积分函数由(3)式得到:
(5)
将变截面悬臂波形钢腹板箱梁等分成四个节段,
表示各截面相对于中性轴的惯性矩,如图3所示。
弯矩作用下的位移可以写成:
(6)
简化上式后可得等效惯性矩Ice为:
(7)
3.1.2. 集中荷载作用下的弯曲变形公式
根据材料力学,承受弯矩、剪力共同作用的直梁,除了弯矩作用下产生弯曲变形外,剪力作用亦会产生剪切变形,即相邻截面间的竖直相互滑动 [17] 。设弯矩引起的弯曲变形为
,并取位移向下为正,x轴向右为正:
(8)
集中荷载作用下悬臂梁自由端梁截面弯矩为
,l为梁长度,P为集中荷载。对式(16)积分得:
(9)
根据边界条件
,可得
,代入式(9)即:
(10)
对式(10)积分可得:
(11)
将边界条件
,代入式(11)式得到
。故集中荷载作用下的弯曲变形公式为:
(12)
只考虑顶底板混凝土的抵抗弯矩作用,因此抗弯刚度为
,其中
为混凝土弹性模量,
为上下混凝土板对中性轴的等效惯性矩。
3.1.3. 自重荷载作用下的弯曲变形计算
变截面箱梁的自重荷载可简化为梯形形式,如图4所示。梯形的上下底长度根据梁的实际重量选定,d为上底与下底的长度的比值。
根据几何计算可得梁重ω:
(13)
式中
为定值,则悬臂梁在自重荷载作用下:
(14)
对式(14)积分得:
(15)
带入边界条件
,可以得到:
(16)
故式(16)可以写成:
(17)
对式(17)积分得:
(18)
把边界条件
,带入式(18)可以计算得到:
(19)
故集中荷载作用下的弯曲变形计算公式为:
![](//html.hanspub.org/file/14-2750494x47_hanspub.png)
Figure 4. Equivalent load diagram of self weight load
图4. 自重荷载等效示意图
(20)
3.2. 变截面波形钢腹板剪切变形计算
3.2.1. 变截面波形钢腹板截面等效
变截面波形钢腹板组合箱梁的钢腹板截面面积随着梁的长度变化,利用辛普森积分公式将其等效为等截面波形钢腹板,等效截面面积
:
(21)
式中
为钢腹板沿梁长度方向的截面高度曲线,t为钢腹板厚度。将式(21)代入辛普森积分公式,得到波形钢腹板的等效截面面积:
(22)
同理可以计算出等效截面顶板和底板厚度。
3.2.2. 变截面波形钢腹板箱梁剪应力计算公式
由于波形钢腹板箱梁特殊的受力方式,计算时一般仅考虑钢腹板的抗剪作用,偏保守与实际不符。因此,本文通过计算截面的剪应力分布,推算混凝土顶底板和波形钢腹板的剪力承担比例,从而考虑混凝土顶底板的抗剪作用。
变截面梁的等效剪应力计算公式:
(23)
τ为截面计算点的剪应力,Q为截面所受剪力,S为横截面上过所求切应力点的横线以外部分的面积对中性轴的面积矩,b为截面计算点的宽度,I为截面惯性矩。由于是组合箱梁,须考虑剪应力计算点分别位于顶板、腹板及底板三种情况。
1) 当剪应力计算点在顶板内时:
(24)
2) 当剪应力计算点在腹板内时:
(25)
3) 当剪应力计算点在底板内时:
(26)
、
分别为顶板的宽度和高度,
、
分别为底板的宽度和高度,
为截面形心,h为截面高度。波纹钢腹板的剪力承担比例:
(27)
式中:
波形钢腹板剪力承担比例,
波形钢腹板的剪应力。
剪切变形引起的剪切变形
:
(28)
式中的波形钢腹板有效剪切模量可由式(2)得到,
。由于梁截面剪应力分布并非均匀,引入应力非均匀系数K。对于波形钢腹板组合箱梁来说,剪应力沿腹板高度接近于等值分布,K取1.0。
4. 算例分析
选取悬臂长度为80 m的单箱单室变截面波形钢腹板组合箱梁为研究对象,箱梁几何尺寸如图5。箱梁结构混凝土强度为C60,弹性模量
,抗压强度标准值
,抗拉强度标准值
,泊松比0.2。波形钢腹板弯折形式为波浪形,选用Q345钢板,弹性模量
,泊松比
。自由端作用集中荷载1000 kN,自重荷载为20687.64 kN。
4.1. 有限元计算模型
采用桥梁有限元分析程序Midas Civil和Abaqus有限元分析软件分别建立梁单元和实体–壳单元模型。其中,Midas Civil梁单元模型,悬臂梁分为18个单元,如图6所示。考虑到横隔板对结构整体影响不大,忽略其影响效应。组合箱梁的模拟计算采用程序中的波形钢腹板截面。
为了验证理论计算的准确性,利用有限元软件Abaqus对该箱梁进行实体建模以及荷载作用和边界条件的模拟,三维实体模型,如图7所示。其中,采用C3D8R单元模拟混凝土结构;采用S4R壳单元模拟波形钢腹板结构;在分析计算中,假设波形钢腹板与混凝土顶底板完全连接,无相对滑移或剪切连接破坏;波形钢腹板具有足够的屈曲强度,无任何形式的屈曲破坏;不考虑钢筋以及混凝土的非线性的影响。
(a) 箱梁横截尺寸
(b) 波形钢腹板尺寸
Figure 5. Section diagram of 80 m cantilever variable cross-section corrugated steel web box girder
图5. 80 m悬臂变截面波形钢腹板箱梁截面图
4.2. 等效截面参数计算
利用建立的等效计算方法,得到悬臂波形钢腹板组合箱梁的等效顶板厚度、底板厚度、钢腹板高度和截面惯性矩,其中顶板厚度不变,按照设计块段划分选择计算截面,如表1所示。
4.3. 弯曲变形分析
利用公式(20)、(28)和表1中的惯性矩,可得该梁在集中荷载和自重荷载作用下的弯曲变形。
在集中荷载作用下两种方法计算的弯曲变形值相对于实体-壳单元模型计算的弯曲变形值误差都较
![](Images/Table_Tmp.jpg)
Table 1. Equivalent post composite box beam parameters
表1. 等效后组合箱梁参数
小,最大悬臂状态下梁单元误差为4.6%;等效惯性矩法误差为1.3%,相对于梁单元误差更小,如图8所示。在自重荷载作用下,最大悬臂状态下的误差梁单元为5.5%;等效惯性矩法误差为2.8%,如图9所示。
4.4. 剪切变形分析
为验证计算的剪切变形的准确性,与梁单元模型和实体-桥单元模型所计算的剪切变形进行比较;梁单元模型的剪切变形由考虑剪切变形的总变形值和不考虑剪切变形的弯曲变形值差值得到;实体–壳单元模型的剪切变形通过约束该梁变形时的转角得到。
计算结果如图10所示,由图可知在集中荷载作用下提出的等效方法计算的剪切变形值相对于梁单元模型计算值更加精确,随着悬臂长度的增加,等效惯性矩法计算值与实体–单元模型计算值误差逐渐减小,最大悬臂状态下等效惯性矩法的误差为13.4%,而梁单元模型计算的误差为34.1%。
自重荷载作用下的剪切变形如图11所示,由图可知利用等效惯性矩法计算的剪切变形相对于梁单元模型更加精确,最大悬臂状态下梁单元计算值误差29.7%,等效惯性矩法计算值误差为15.8%。
由于梁单元模型默认的有效抗剪面积为钢腹板截面面积,并不考虑混凝土的抗剪能力,故导致计算的剪切变形值偏大。
5. 顶底板混凝土抗剪能力对剪切变形的影响分析
该梁在自重荷载作用下各截面的波形钢腹板剪应力如图12所示,计算中假设腹板剪应力沿腹板高度均匀分布。由图可以看出,利用等截面梁计算公式计算得到的腹板剪应力相对与腹板承受全部剪力时的
![](//html.hanspub.org/file/14-2750494x84_hanspub.png)
Figure 8. Calculation of bending deformation under concentrated load by different methods
图8. 不同方法计算在集中荷载作用下的弯曲变形
![](//html.hanspub.org/file/14-2750494x85_hanspub.png)
Figure 9. Bending deformation calculated by different methods under the action of self weight load
图9. 不同方法计算在自重荷载作用下的弯曲变形
![](//html.hanspub.org/file/14-2750494x86_hanspub.png)
Figure 10. Calculation of shear deformation under concentrated load by different methods
图10. 不同方法计算在集中荷载作用下的剪切变形
![](//html.hanspub.org/file/14-2750494x87_hanspub.png)
Figure 11. Calculation of shear deformation under the action of self weight load by different methods
图11. 不同方法计算自重荷载作用下的剪切变形
![](//html.hanspub.org/file/14-2750494x88_hanspub.png)
Figure 12. Shear stress comparison of steel web
图12. 钢腹板剪应力对比
剪应力要小,这是由于考虑了混凝土顶底板的抗剪作用,因此计算结果更加精确;该梁梁端固结处截面的腹板实际承担剪力的比例86.5%,与腹板承受全部剪力时的剪切变形进行比较,如图13所示。由于腹板承受的剪力变小,故剪切变形相应减小,且与实体-壳单元模型计算值更加吻合,考虑混凝土抗剪能力
![](//html.hanspub.org/file/14-2750494x89_hanspub.png)
Figure 13. Shear deformation under different shear load ratio of the web
图13. 腹板不同剪力承载比下的剪切变形
时的剪切变形误差为15.8%,不考虑混凝土抗剪能力的剪切变形误差为33.9%。
6. 结论
1) 实体–壳单元模型可以真实反映波形钢腹板的受力情况,但建模过程复杂,计算效率不高;梁单元模型计算结果误差较大,波形钢腹板的安装对变形精度要求较高;本文基于辛普森公式计算的等效惯性矩,避免了当惯性矩沿截面变化时繁琐的积分计算,利用等效惯性矩计算得到的弯曲变形和剪切变形与实体–壳单元模型的计算结果吻合度较高。
2) 本文基于辛普森积分公式,将波形钢腹板组合箱梁的变截面形式等效为等截面,建立了计算变截面波形钢腹板组合箱梁变形的等效惯性矩法,该计算方法考虑了混凝土的抗剪能力,计算结果更贴近实际。
基金项目
国家自然科学基金面上项目(51578206, 51478159);安徽省杰出青年科学基金(1708085J06);中央高校基本科研业务费专项资金项目(PA2017GDQT0022)。