波形钢腹板PC箱梁剪力滞效应的能量变分法解

2022-03-29王社良

西安建筑科技大学学报（自然科学版） 2022年1期

高丽，周媛，王社良，刘博

(1.淮阴工学院计算机与软件工程学院，江苏淮安 223003；2.西安建筑科技大学土木工程学院，陕西西安 710055；3.长安大学理学院，陕西西安 710061；4.中铁二十局集团有限公司，陕西西安 710016)

波形钢腹板预应力混凝土(PC)箱梁桥是将薄钢板弯折成波浪形后替换普通混凝土箱梁腹板后发展成的新型组合桥梁.由于波形钢腹板高强轻质的特性，箱梁的跨度有了新的突破，动力响应弱，抗震性能大为改善.凭借波形钢腹板的褶皱效应，腹板几乎不抵抗轴向力，由抗弯效率高的混凝土土顶底板承担竖向弯矩，受力更为合理、明确[1-3].钢腹板上可支模，架设迅速，维护方便，造型优雅，在50～100 m跨径桥梁方案极具竞争力.波形钢腹板PC箱梁桥自20世纪80年代法国提出后，在欧洲、日本得到成功实践，近年来，波形钢腹板PC箱梁桥在国内也得到了迅速发展[4-5].

与普通PC箱梁一样，波形腹板PC箱梁顶板也会受剪切变形影响，横向变形不均匀，导致正应力沿宽度方向分布不均匀[6-7].已有的分析表明：波形钢腹板PC组合箱梁的剪力滞问题比普通预应力混凝土箱梁更加突出.波形腹板PC箱梁设计时不能直接套用普通PC箱梁剪力滞的研究结论[8-9].在运用能量变分法讨论PC箱梁的剪力滞效应时，许多学者拟定了二次、三次及四次曲线作为翘曲函数，其概念清晰，求解简单，取得了丰硕的成果[10-13].纵向翘曲位移函数与箱梁真实翘曲位移吻合程度决定了分析结果误差大小.

对于波形钢腹板组合箱梁的剪力滞问题，理论分析尚不全面.本文引入β阶曲线和近似余弦函数作为箱梁纵向翘曲的位移模式，推导了等截面波形钢腹板组合箱梁剪力滞系数的近似解，通过与空间有限元法计算结果比较，提出一个可靠度高，通用性强，方便实用的波形腹板PC箱梁剪力滞系数计算方法.

1 β次翘曲位移函数的变分解

基本假定：(1)由于波形钢腹板预应力混凝土箱梁在竖向荷载作用下服从“拟平截面假定”[14],可引入两个广义位移函数来描述翼板变形：

(1)

(2)

式中：u(x)为翼缘板最大转角差，其余几何参数如图1示.

图1 波形钢腹板箱梁截面参数

(2)翼缘板为应力沿厚度均匀分布的薄板，忽略板的平面外剪应变和横向弯曲变形等微小量，不考虑波形钢腹板的剪切变形.

当箱梁弯曲时，外力势能：

(3)

应变—位移几何关系为

(4)

上下翼板应变能分别为

(5)

式中：E、G分别为翼缘板的弹性模量和剪切模量.忽略翼缘板自身惯性矩，截面惯性矩为

(6)

(7)

由最小势能原理有

δ∏=0

(8)

分部积分后可得：

(9)

根据变分原理有

(10)

(11)

(12)

可得简化后的表达式:

(13)

式中：M、Q分别为计算截面的弯矩和剪力.

当等截面简支梁承受集中荷载P时，解的一般形式为

(14)

式中：u*为剪力Q(x)分布有关的特解；C1，C2为待定系数，由边界条件确定.

翼缘板纵向正应力：

(15)

固支边界：

u=0

(16)

(17)

同理，可得到其他力学图式和荷载工况下的位移和边界条件的具体表达式，最后导出应力并将剪力滞系数λ1汇总到表1.

2 近似余弦函数的变分解

余弦函数的泰勒展开公式：

(18)

由上式可以看出，余弦函数本质上是一个高阶多项式，高次幂的系数很小，一般取前三项已有足够的精度，常数项不影响应力的横向分布，可以假定翘曲位移模式为

(19)

式中：ξ,η,ζ∈[-1,1]任意参数，为保证腹板处位移的连续条件还应有：ξ+η+ζ=1

(20)

式中：

(21)

(22)

(23)

运用最小势能原理和变分原理，最后得到应力表达式及边界条件，并将剪力滞系数λ2汇总于表1：

表1 剪力滞系数实用计算公式汇总

(24)

(25)

续表1

3 算例分析

采用两种有限元分析软件ANSYS和MIDAS/FEA同时建立计算跨径为30 m的简支梁和15 m的悬臂梁的1/4空间模型，如图2所示.

图2 有限元模型

模型梁为单箱单室波形钢腹板箱梁，梁高3 m，顶板宽5.5 m，悬臂长1.25 m，底板宽3.5 m.顶底板板厚0.4 m.端部设厚度为0.9 m的端隔板.波形钢腹板厚为10 mm，波长0.9 m，直幅宽0.25 m，斜幅宽0.2 m，波高0.12 m.混凝土强度C55，弹模为3.5×104MPa，泊松比为0.2，采用带中间节点的高次实体单元模拟.钢腹板材质为Q345c优质低碳钢，弹模2.06×102GPa，泊松比0.3，选用shell63板单元建模，并与混凝土顶底板及端隔板共节点.在对称面施加对称约束，简支梁的支座位置采用面上铰约束，悬臂梁的固定端约束全部自由度.暂不计入预应力单元的影响.由于混凝土顶底板在钢混共节点区域刚度不均匀，受力复杂，恰好又是最大剪力滞系数的敏感位置，网格务需精细、中正，竖向至少4层单元.模型采用柔度较好的20个节点单元solid95，共计41 752个单元和198 917个节点.

表2给出了集中荷载20 kN下简支梁和悬臂梁典型截面剪力滞系数.

表2 集中荷载20 kN下典型截面剪力滞系数

从表2可以看出，当引入的曲线从2次增加至20次，简支梁的最大剪力滞系数理论值由1.20增加到1.29，有限元值为1.17.悬臂梁的最大剪力滞系数从1.20上升至1.30，有限元值为1.23.翘曲位移函数的阶次愈大，最大剪力滞系数的计算结果也愈大.可见，翘曲位移函数的阶次太高会导致结果偏离真实值，采用三次抛物线是个比较保守，富余量又不大的选择.

表3给出了10 kN/m均布荷载作用下简支梁和悬臂梁典型截面剪力滞系数.

表3 均布荷载10 kN/m下典型截面剪力滞系数

简支梁在均布荷载下的最大剪力滞系数为1.02，且对曲线阶次并不敏感，工程中可以忽略跨中剪力滞引起的应力集中效应.悬臂梁在均布荷载下，剪力滞系数为1.23左右，比集中荷载下更为显著，因此，悬臂施工的薄壁箱梁必须重视剪力滞因素引起的不均匀现象.

由图3可见，变分法求解波形钢腹板预应力混凝土箱梁的最大剪力滞系数可以很好地逼近有限元值，但该法将翼缘板做了拟平面应力假定，自由端存在一定的误差，并且从有限元分析可以看出翼板悬挑段与箱室内段的剪力滞系数并非完全对称，可以考虑拟定多条翘曲函数来分别计算.

图3 简支梁集中荷载下跨中截面剪滞系数分布

图4 悬臂梁根部剪力滞系数横向分布

悬臂梁在均布荷载下的最大剪力滞系数值1.34，二次曲线的变分值为1.38.然而，在集中荷载下最大剪力滞系数为1.23，高于二次曲线的求解值1.20，可以考虑采用4次曲线来解答，对于不同的荷载工况，也可以考虑采用不同阶次的曲线.

4 超静定桥梁的剪力滞解答

由最小势能原理推出超静定波形钢腹板组合箱梁桥的剪力滞系数存在两个难点：一是联立求解困难，公式繁杂，应用不便；二是超静定箱梁翼缘板有效宽度的变化影响弯矩M(x)沿跨长方向的分布.可以结合表1的实用公式应用肢解法和叠加法[15]，使分析得以简化.

连续梁桥，斜拉桥，连续刚构等超静定桥梁在各种荷载工况下都存在多个弯矩为0的反弯点，可以人为地在该点将桥梁拆解成若干根简支梁或悬臂梁，最后运用表1的实用公式间接计算超静定结构的剪力滞系数，这是肢解法.

在变分法求解过程中，仅假定了纵向翘曲位移的横向分布，纵向和竖向位移依然满足连续性要求和边界条件，考虑到约束的本质是对桥梁施加相应的支座反力，可以用等代力替代赘余约束，将超静定桥梁简化为静定形式，根据线弹性阶段的应力叠加原理，将叠加原理推广到箱梁剪力滞的求解中，这是叠加法.计算截面I的剪力滞系数表达式为：

式中：M为I截面真实弯矩；Mi为基本体系在单个工况下I截面处的弯矩；λi为基本体系单个工况下I截面的剪力滞系数；λ为I截面真实剪力滞系数；

以第3节中的模型梁截面为基本截面，分别建立3 m×30 m的三跨连续梁和连续刚构实体有限元模型，得到工况1：中跨跨中和边跨跨中同时施加集中荷载10 kN时典型截面的剪力滞效应，见表4；工况2：全桥均布荷载10 kN/m下，典型截面的剪力滞效应，见表5.

表4 工况1典型截面剪力滞系数

表5 均布荷载下典型截面剪力滞系数

从表4和表5可以看出，肢解法、叠加法和有限元结果基本吻合，虽然有些误差，但相对误差不大于3.4%，足以满足工程精度的要求.

5 剪力滞效应的参数分析

改变有限元模型的桥宽，建立宽跨比为0.1～1的10组波形钢腹板预应力悬臂箱梁模型，分析在悬臂端集中荷载和均布荷载工况下的剪力滞效应.

从图5可以看出，宽跨比是影响箱梁剪力滞效应的重要几何参数之一.悬臂梁在相同荷载工况下，顶底板的剪力滞系数基本一致.均布荷载工况的剪力滞效应对宽跨比的改变更为敏感.当宽跨比由0.1递增到1时，均布荷载下的最大剪力滞系数相应由1.19上升至2.67；集中荷载下的最大剪力滞从1.07浮升至2.00.因此，本文按照最小二乘法给出悬臂梁根部剪力滞系数y和宽跨比x的关系函数.