钢-混凝土组合梁静动力响应
张云龙1, 刘占莹2, 吴春利3, 王静1
1.吉林建筑大学 交通科学与工程学院,长春130118
2.长春建筑学院 交通学院,长春 130607
3.吉林大学 交通学院,长春 130022
吴春利(1978-),女,讲师,博士.研究方向:桥梁结构健康监测理论.E-mail:clwu@jlu.edu.cn

作者简介:张云龙(1975-),男,副教授,博士.研究方向:组合梁结构理论.E-mail:zyl_ql@163.com

摘要

为了准确分析钢-混凝土组合梁的静力和动力响应,基于接触理论和有限单元法的基本思想,提出了组合梁单元的合理位移函数,推导了考虑剪切滑移效应的组合梁单元刚度方程,并结合动力学经典理论,给出了组合梁无阻尼状态下考虑交界面剪切滑移效应的各阶自振频率和振型计算方法。通过具体算例分析了钢-混凝土简支组合梁在竖向荷载作用下的挠度和交界面的剪切滑移应变差沿梁长的分布形态,并与试验结果进行了对比,结果表明:有限元解与试验结果基本一致,说明该有限元解是合理、可信的。利用上述刚度矩阵,基于经典动力学理论,计算了组合梁在考虑剪切滑移效应时的频率和振型,并与完全相互作用的结果进行了对比分析,说明在进行动力响应计算时应考虑剪切滑移的影响,为实际工程的设计提供了一种可靠的理论计算方法。

关键词: 桥梁工程; 组合梁; 接触理论; 有限单元法; 动力响应
中图分类号:TU311.1 文献标志码:A 文章编号:1671-5497(2017)03-0789-07
Static and dynamic responses of steel-concrete composite beams
ZHANG Yun-long1, LIU Zhan-ying2, WU Chun-li3, WANG Jing1
1.School of Transportation Science & Engineering, Jilin Jianzhu University, Changchun 130118, China
2.School of Transportation, Changchun Architecture and Civil Engineering College, Changchun 130607, China
3.College of Transportation,Jilin University, Changchun 130022, China
Abstract

In order to accurately analyze the static and dynamic responses of steel-concrete composite beams, the reasonable element displacement function of the composite beam is derived based on the contact theory and finite element method. The stiffness equation of the composite beam element considering shear slip effect is deduced. Combined with classical dynamics theory, a method of calculating the natural frequencies and vibration mode of the composite beams without damping is established with consideration of the shear slip effect at the interface. As a case study, the deflection of a simply supported steel-composite beam under vertical load and distributed shear slip strain difference along the beam is analyzed. Results show that the finite element solution is in good agreement with the experiment results. Using stiffness matrix the natural frequencies and vibration modes of the beam considering the shear slip effect are calculated and compared with complete interaction results, which suggest that the effect of shear slip effect should be considered in dynamic response calculation. This study provides a reliable calculation method for practical engineering design.

Key words: bridge engineering; composite beam; contact theory; finite element method; dynamic response
0 引言

钢-混凝土组合梁桥是在钢结构的基础上发展而来的一种新型结构形式。将混凝土板设置在梁的受压部分, 将抗拉强度高的钢材设置在梁的受拉部分, 两种材料按这种方式结合在一起可以充分发挥各自的优势[1], 中国许多桥梁结构都采用这种形式[2]。近年来, 清华大学聂建国[3, 4, 5, 6]对钢-混凝土组合梁进行了系列研究, 推动了组合梁在实际工程中的广泛应用, 取得了良好的经济效果。由于组合梁结构的广泛应用, 对其抗剪连接件受力性能进行研究显得尤为重要。连接件在钢梁与混凝土板交界面处传递剪力, 使其成为整体, 但连接件的刚度并非无限大, 在传递剪力的过程中会产生变形, 从而导致交界面处产生滑移[7]。Ranzi等[8]分析了刚度及变形等效应对组合梁的影响, 提出了考虑滑移和不考虑滑移时的有限元分析方法。试验及研究结果表明[9, 10, 11], 组合梁的滑移效应会对结构的内力和变形产生不可忽略的影响。

在组合梁的动力响应方面, 由经典动力学理论及已有公式可以对不考虑交界面剪切滑移效应的组合梁进行模态分析[12]。目前, 我国对组合梁静力特性的研究已经较为深入, 但对其动力特性, 尤其是在考虑组合梁交界面处剪切滑移效应时的动力特性研究还处于起步阶段[13, 14, 15, 16]。东南大学王景全等[17, 18]通过室内试验分析了钢-混凝土交界面的滑移效应对组合梁动力性能的影响, 指出交界面滑移效应对组合梁自振频率的影响是不容忽视的。所以准确地对振型的形态和频率进行理论分析显得尤为重要。

综上可知, 目前钢-混凝土组合梁的设计计算方法仍停留在对经验公式的总结和特定条件下的解析解推导, 造成设计理论滞后于生产实际。由于设计方面常用的软件均采用有限元法作为计算的内核, 因此推导可靠的考虑剪切滑移效应的组合梁有限元解是在设计计算阶段亟待解决的问题。文献[7]中给出的有限元解会出现剪切锁死现象, 影响了计算公式的推广和使用。为了解决这一现象, 本文提出了组合梁有限元解合理的位移模式, 推导了一种新型有限元解, 并经过了试验验证, 为考虑剪切滑移效应的组合梁静力分析和动力响应研究奠定了理论基础。

1 试探函数的建立
1.1 基本假定

由于组合梁在受力过程中可以忽略竖向掀起效应产生的影响。二者在外力作用下能够保持共同的曲率, 由组合梁的基本性质可以做出如下基本假定:①钢梁和混凝土顶板在变形过程中能够保持平截面; ②在钢梁和混凝土顶板交界面处存在剪切滑移效应, 交界面处的剪力集度与相对水平位移差成正比; ③忽略钢梁与混凝土顶板之间的竖向掀起效应。

1.2 组合梁位移函数的建立

1.2.1 基本位移函数

图1 梁高方向的位移状态Fig.1 Displacement state of beam high direction

根据基本假定, 组合梁的基本位移函数有3个, 即混凝土顶板的轴向位移uc; 钢梁的轴向位移us; 组合梁的竖向挠度v。如图1所示, 3个基本函数可表示为:

ucx, zc=uc0x+zcφxvx=vxusx, zs=us0x+zsφx(1)

式中:uc0(x)为混凝土顶板形心处的轴向位移; us0(x)为钢梁形心处的轴向位移; v(x)为组合梁竖向位移; φ (x)为组合梁的转角, 根据弹性力学的基本原理, φ (x)可表示为:

φx=dvx/dx(2)

由基本假定可知, 混凝土顶板底面与钢梁顶面之间的轴向位移差为:

Δux=usux-ucbx=us0x-uc0x+φxzcb+zsu=us0x-uc0x+v'xzd(3)

式中:zd=zcb+zsu, zcb为混凝土顶板中性轴至组合梁交界面的距离, zsu为钢梁中性轴至组合梁交界面的距离; usu(x)为钢梁顶面的轴向位移; ucb(x)为混凝土顶板底面的轴向位移。Δ u(x)即为组合梁交界面处的剪切滑移应变差。

1.2.2 组合梁的物理方程

由基本假定可知, 混凝土顶板与钢梁交界面处的剪力集度为:

qx=kΔux=kus0x-uc0x+v'xzd(4)

式中:k为剪力连接件的剪切刚度。

图2为微元体受力状态, 根据力的平衡关系可得:

qx=-dNsdx=dNcdx(5)

图2 交界面微元体受力图Fig.2 Force diagram of interface volume element

根据材料力学的基本原理, 组合梁内力与位移之间的关系为:

Ncx=EcAcεc0x=EcAcu'c0xNsx=EsAsεs0x=EsAsu's0xM0x=EsIskx+EcIckx=B0vx(6)

式中:EcIcAc分别为混凝土顶板的弹性模量、惯性矩和横断面面积; EsIsAs分别为钢梁的弹性模量、惯性矩和横断面面积; k(x)为组合梁在x截面处的曲率, k(x)=v″(x); B0=EcIc+EsIs; Nc(x)、Ns(x)分别为混凝土顶板和钢梁的轴力, 以受压为正, 受拉为负。

1.2.3 组合梁的试探函数

设混凝土顶板的轴向位移函数为:

uc0x=c1+c2x+c3x2+c4x3(7)

根据式(5)(6)(7)可得:

us0x=-2αc3-6c4αx(8)

式中:α =EcAc/(EsAs)。

对式(8)进行积分可得:

us0x=-c3αx2-c4αx3+c5x+c6(9)

由式(4)(5)可得:

v'x=1zduc0x-1zdus0x-βzdus0x(10)

式中:β =EsAs/K

由式(7)(9)(10)可得组合梁的转角为:

v'x=1zdc2x+1+αzdc3x2+2αβzdc3+6αβzdc4x+1+αzdc4x3-1zdc5x+1zdc7(11)

对式(11)进行积分可得组合梁的竖向位移为:

vx=12zdc2x2+1+α3zdc3x3+2αβzdc3x+3αβzdc4x2+1+α4zdc4x4-12zdc5x2+1zdc7x+c8(12)

式中:c1~c8为待定系数, 设单元左侧i结点的位移分别为uciusiviφ i, 单元右侧j结点的位移分别为ucjusjvjφ j, 代入式(7)(9)(11)(12)中可解得8个待定系数, 由此可得到组合梁的位移函数为:

ue=Nδe(13)

式中:ue=[ uc0us0vφ]T为单元的位移试探函数; δ e=[uci usi vi φ i ucj usj vj φ j]T为单元的结点位移列阵; N为位移函数矩阵, 令α 1=1+α , β 1=α 1l l-2x-36α β , 则矩阵N可写成如下形式:

N=A1A2A3A4A5A6A7A8B1B2B3B4B5B6B7B8C1C2C3C4C5C6C7C8

式中:A1=1- xl+ (l-x)(l-2x)x12αβl;

A2=- (l-x)(l-2x)x12αβl;

A3=- 6zd(l-x)xα1l3;

A4=- β1(l-x)xzd12α1αβl2;

A5= xl- (l-x)(l-2x)x12αβl;

A6= (l-x)(l-2x)x12αβl;

A7= 6zd(l-x)xα1l3;

A8= β1(l-x)xzd12α1αβl2;

B1=- (l-x)(l-2x)x12βl;

B2=1- xl+ (l-x)(l-2x)x12βl;

B3= 6zd(l-x)αxα1βl3;

B4= β1(l-x)zdαx12α1αβl2;

B5= (l-x)(l-2x)αx12αβl;

B6= xl- (l-x)(l-2x)αx12αβl;

B7=- 6zd(l-x)αxα1l3;

B8=- β1(l-x)zdαx12α1αβl2;

C1= α1(l-x)2x224αβzdl;

C2=- α1(l-x)2x2+12β1-α24αβzdl;

C3=1- 3x2l2+ 2x3α1l3;

C4=x- 2x2l+ x3l2- α1(l-x)2x224αβl;

C5=- α1(l-x)2x224αβzdl;

C6= α1(l-x)2x224αβzdl;

C7= 3x2l2- 2x3α1l3;

C8=- x2l+ x3l2+ α1(l-x)2x224αβl

2 考虑剪切滑移效应的组合梁单元刚度方程

根据式(13)可得单元内任意截面的应变为:

ε1=B1δe(14)

式中:ε 1为单元的应变列阵, 可表示为:

ε1=εc0xεs0xkxT(15)

B1为应变矩阵, 可用如下矩阵表示:

B1=b11b12b13b14b15b16b17b18b21b22b23b24b25b26b27b28b31b32b33b34b35b36b37b38(16)

B1中的各元素可按式(17)求解:

bij=dnijdxb3j=d2n3jdx2(17)

式中:i=1, 2; j=1, 2, …, 8; nij为位移函数矩阵N中对应的元素。

根据式(3)(13)可得组合梁交界面处的剪切滑移量为:

ε2=Δux=B2δe(18)

式中:

B2=b41b42b43b44b45b46b47b48(19)

矩阵中各元素可采用式(20)求解:

b4i=n2i-n1i+dn3idxzd(20)

将式(15)(18)合并可写成:

ε=ε1ε2=εc0xεs0xkxΔuxT(21)

将式(16)(19)合并可写成:

B=B1B2T=b11b12b13b14b15b16b17b18b21b22b23b24b25b26b27b28b31b32b33b34b35b36b37b38b41b42b43b44b45b46b47b48(22)

基于弹性力学的基本理论, 根据式(4)(6)给出的内力与位移的关系, 可得组合梁单元的势能函数为:

π=0l12εTdx-δeTFe(23)

式中:D为结构的刚度;

D=EcAcEsAsB0k

Fe为单元结点荷载列阵, 由弹性力学基本原理可知:

Fe=0lNTqdx+jNTxPj+kdNTxdxMkl(24)

式中:q为组合梁单元受到的均布荷载; Pj为集中力; Mk为集中力矩。

根据最小势能原理, 对组合梁单元的势能函数取变分, 由式(24)可得:

δπ=0lδδeTBTDBδedx-δδeTFe=0(25)

由式(25)可知, 考虑剪切滑移效应的组合梁单元刚度方程为:

Keδe-Fe=0(26)

式中:Ke为单元刚度矩阵, 可按式(27)求解:

Ke=0lBTDBdx(27)

根据式(27)即可基于常规有限单元法基本思路求解考虑剪切滑移效应的组合梁的静力响应。

3 有限元解的理论验证

将理论分析与文献[19]给出的试验数据进行对比, 结果如图3图4所示。从图3图4可以看出:本文计算值与试验值有较好的一致性, 尤其是挠度的计算结果与试验值吻合得很好, 说明该数值解能够对组合梁在竖向荷载作用下的抗弯性能进行合理的分析和计算。

图3 试验梁的挠曲线Fig.3 Deflection curve of test beam

图4 剪切滑移应变差Fig.4 Shear slip strain difference

4 考虑剪切滑移时梁的动力响应

根据式(26)所示单元刚度矩阵和结构动力学可知, 无阻尼多自由度体系自由振动方程的推导过程如下。

Mu+Ku=0(28)

式中:MK分别为组合梁的质量和刚度矩阵。

假设多自由度体系自由运动为简谐振动, 有:

ut=φsinωt+θ(29)

式中:φ 表示体系的形状(它不随时间t而变化), 对式(29)求二次导数, 得到自由振动的加速度为:

ut-ω2φsinωt+θ=-ω2ut(30)

将式(29)(30)代入式(27)有:

-ω2sinωt+θ+sinωt+θ=0(31)

由于sin ωt+θ为任意正弦函数, 故有:

K-ω2Mφ=0(32)

式(32)有非零解的条件是:

K-ω2M=0(33)

图5 组合梁考虑剪切滑移时的竖向振型Fig.5 Vertical vibration mode of composite beam considering shear slip effect

图6 组合梁考虑剪切滑移时的轴向振型Fig.6 Axial vibration mode of composite beam considering shear slip effect

本文根据试验梁的尺寸, 在考虑剪切滑移且不考虑阻尼的情况下得到各阶竖向和轴向的自振频率和振型, 如图5图6所示(此处取前两阶的频率和振型)。

表1表2中可以看出:当考虑剪切滑移时, 组合梁的自振频率较完全相互作用有明显减小, 所以在实际进行动载测试过程中, 应考虑剪切滑移的影响。

表1 竖向自振频率和振型 Table 1 Vertical natural frequency and vibration mode
表2 轴向自振频率和振型 Table 2 Axial natural frequency and vibration mode
5 结论

(1)基于接触理论及有限元基本思想推导了考虑剪切滑移效应时的钢-混凝土组合梁的刚度方程。将计算值与试验结果进行对比分析, 验证了该理论的可靠性及合理性。

(2)得到考虑剪切滑移效应的组合梁频率和振型, 并与完全相互作用的结果进行对比分析, 说明在进行动力响应计算时应考虑剪切滑移的影响。

(3)本文给出的2结点组合梁有限元解可以完美地嵌入到目前流行的桥梁和结构设计软件中, 促进了考虑剪切滑移效应的组合梁设计理论的实际推广。

The authors have declared that no competing interests exist.

参考文献
[1] 范旭红, 石启印, 马波. 钢-混凝土组合梁的研究与展望[J]. 江苏大学学报, 2004, 25(1): 89-92.
Fan Xu-hong, Shi Qi-yin, Ma Bo. Development and perspective of steel-concrete composite beams[J]. Journal of Jiangsu University, 2004, 25(1): 89-92. [本文引用:1]
[2] 陈宝春, 牟廷敏, 陈宜言, . 我国钢-混凝土组合结构桥梁研究进展及工程应用[J]. 建筑结构学报, 2013, 34(增刊1): 1-10.
Chen Bao-chun, Mu Ting-min, Chen Yi-yan, et al. State-of-the-art of research and engineering application of steel-concrete composite beidges in China[J]. Journal of Building Strctures, 2013, 34(Sup. 1): 1-10. [本文引用:1]
[3] 聂建国, 余志武. 钢-混凝土组合梁在我国的研究及应用[J]. 土木工程学报, 1999, 32(2): 3-8.
Nie Jian-guo, Yu Zhi-wu. Research and practive of composite steel-concrete beams in China[J]. China Civil Engineering Journal, 1999, 32(2): 3-8. [本文引用:1]
[4] 聂建国, 沈聚敏. 滑移效应对钢-混凝土组合梁弯曲强度的影响及其计算[J]. 土木工程学报, 1997, 30(1): 31-36.
Nie Jian-guo, Shen Ju-min. Slip effect on strength of composite steel-concrete beams[J]. China Civil Engineering Journal, 1997, 30(1): 31-36. [本文引用:1]
[5] 聂建国, 沈聚敏, 余志武. 考虑滑移效应的钢-混凝土组合梁变形计算的折减刚度法[J]. 土木工程学报, 1995, 28(6): 11-17.
Nie Jian-guo, Shen Ju-min, Yu Zhi-wu. A reduced rigidity method for calculating deformation of composite steel-concrete beams[J]. China Civil Engineering Journal, 1995, 28(6): 11-17. [本文引用:1]
[6] 聂建国, 陶慕轩. 体外预应力钢-混凝土组合梁受力性能的研究现状与展望[J]. 工程力学, 2013, 34(增刊2): 129-141.
Nie Jian-guo, Tao Mu-xuan. Research status and perspective on externally prestressed steel-concrete composite beams[J]. Engineering Mechanics, 2013, 34(Sup. 2): 129-141. [本文引用:1]
[7] 刘寒冰, 刘天明, 张云龙. 钢-混凝土组合连续梁抗弯性能[J]. 吉林大学学报: 工学版, 2009, 39(6): 1486-1491.
Liu Han-bing, Liu Tian-ming, Zhang Yun-long. Bending resistance of steel-concrete composite continuous beam[J]. Journal of Jilin University (Engineering and Technology Edition), 2009, 39(6): 1486-1491. [本文引用:1]
[8] Ranzi G, Zona A. A steel-concrete composite beam model with partial interaction including the shear deformability of the steel component[J]. Engineering Structures, 2007, 29(11): 3026-3041. [本文引用:1]
[9] Fragiacom M, Amaido C, Macorini L. Finite element model for collapse and long-term analysis of steel-concrete composite beams[J]. Journal of Structural Engineering, 2004, 130(3): 389-497. [本文引用:1]
[10] Salari M R, Spacone E, Shing P B, et al. Nonlinear analysis of composite beams with deformable shear connectors[J]. Journal of Structural Engineering, 1998, 124(10): 1148-1158. [本文引用:1]
[11] Sebastion W M, Connel R. Nonlinear FE analysis of steel-concrete composite[J]. Journal of Structural Engineering, 2000, 126(6): 662-674. [本文引用:1]
[12] 周旺保, 蒋丽忠, 李书进. 钢-混凝土组合箱梁弯矩曲率恢复力模型研究[J]. 建筑结构学报, 2015, 36(1): 78-84.
Zhou Wang-bao, Jiang Li-zhong, Li Shu-jin. Moment-curvature restoring force model research of steel-Concrete composite box-beams[J]. Journal of Building Structures, 2015, 36(1): 78-84. [本文引用:1]
[13] 侯忠明, 夏禾, 王元清, . 钢-混凝土组合梁动力折减系数研究[J]. 振动与冲击, 2015, 34(4): 74-81.
Hou Zhong-ming, Xia He, Wang Yuan-qing, et al. Dynamic reduction coefficients for a steel-concrete composite beam[J]. Journal of Vibration and Shock, 2015, 34(4): 74-81. [本文引用:1]
[14] 戚菁菁, 蒋丽忠, 张传增, . 界面滑移、竖向掀起及剪切变形对钢-混凝土组合连续梁动力性能的影响[J]. 中南大学学报, 2010, 41(6): 2334-2343.
Qi Jing-jing, Jiang Li-zhong, Zhang Chuan-zeng, et al. Effects of interface slip, vertical uplift and shear deformation on dynamic behavior of steel-concrete composite continuous beams[J]. Journal of Central South University, 2010, 41(6): 2334-2343. [本文引用:1]
[15] Banerjee J R. Frequency equation and mode shape formulae for composite Timoshenko beams[J]. Composite Structures, 2001, 51(4): 381-388. [本文引用:1]
[16] 沈旭栋, 陈伟球, 徐荣桥. 有轴力的部分作用组合梁的动力分析[J]. 振动工程学报, 2012, 25(5): 514-520.
Shen Xu-dong, Chen Wei-qiu, Xu Rong-qiao. Dynamic analysis of partial-interactive composite beam with axial force[J]. Journal of Vibration Engineering, 2012, 25(5): 514-520. [本文引用:1]
[17] 王景全, 殷惠光, 张书兵, . 钢-混凝土组合梁考虑界面滑移的振动模态试验研究[J]. 建筑结构学报, 2016, 37(2): 142-149.
Wang Jing-quan, Yin Hui-guang, Zhang Shu-bing, et al. Experimental investigation on dynamic properties of steel-concrete composite beams considering effect of interfcial slip[J]. Journal of Building Structures, 2013, 37(2): 142-149. [本文引用:1]
[18] 张书兵, 王景全, 李明, . 考虑界面滑移效应的组合梁自振频率计算的修正折减刚度法[J]. 土木工程学报, 2015, 48(12): 41-49.
Zhang Shu-bing, Wang Jing-quan, Li Ming, et al. A modified stiffness reduction method for calculating the natural frequencies of composite beams considering the effect of interfacial slippage[J]. China Civil Engineering Journal, 2015, 48(12): 41-49. [本文引用:1]
[19] 张云龙. 体外预应力钢-混凝土组合梁结构行为的实验研究[D]. 长春: 吉林大学交通学院, 2005.
Zhang Yun-long. Test study of external prestressed composite steel-concrete beams structure[D]. Changchun: College of Transportation, Jilin University, 2005. [本文引用:1]