基于Timoshenko 梁理论的钢-混组合梁动力刚度矩阵法

孙琪凯，张 楠，刘 潇

(北京交通大学土木建筑工程学院，北京 100044)

钢-混组合梁因其自重较轻、承载力高、刚度大等优点而得到了越来越广泛的应用。常见的钢-混组合梁结构形式是由柔性剪力键连接钢梁和混凝土子梁而成，这种结构既能发挥钢材受拉又能发挥混凝土受压的材料特点[1-3]。

钢-混组合梁的动力学研究已比较常见[4-7]。早期，基于Euler-Bernoulli 梁理论，Girhammar 等[8]推导了考虑界面剪切滑移效应的钢-混组合梁的运动平衡微分方程，指出了界面剪切滑移效应对组合梁自振频率的折减效应是不容忽视的。Huang 和Su[9]提出了影响组合梁频率折减量的两个无量纲参数，组合连接系数和上下层截面抗弯刚度比。Wu 等[10]推导了梁端轴力作用下组合梁自振频率计算公式，以简支梁为参考，给出了悬臂、简支-固支和固支-固支等三种典型边界条件下自振频率的近似表达式。侯忠明等[11-12]提出了“动力折减系数”的概念，包括“刚度折减系数”和“频率折减系数”两部分。从理论和实验的角度讨论了“动力折减系数”随剪力键刚度的变化规律。孙琪凯等[13]从能量法的角度给出了有限元计算公式，可用于分析轴向变抗弯刚度的钢-混组合梁的动力特性。对于细长梁的低阶频率的分析，Euler-Bernoulli梁理论计算精度一般可以满足工程要求。但是，对于短粗梁或者进行高阶频率分析时，如果继续忽略钢-混组合梁剪切变形和转动惯量的影响，则计算误差会增大，不再满足工程要求。因此，部分学者[14-17]将Timoshenko 梁理论用于组合梁动力分析。Xu 等[14]推导了Timoshenko组合梁运动微分方程，说明了组合梁动力分析时，剪切变形和转动惯量是不可忽略的。Lin等[15-16]给出了刚度矩阵和质量矩阵，从有限元法的角度求解了组合梁的动力性能。以上研究中虽然考虑了剪切变形和转动惯量的影响，但是均假设混凝土子梁与钢梁的转角相等，这与组合梁的实际运动状态是不符的。因为混凝土子梁与钢梁的剪切角与各自的剪切模量有关，而两者的剪切模量具有很大的差异。Nguyen 等[17]给出了考虑混凝土子梁和钢梁剪切角不同时的组合梁运动微分方程，得到了频率解析解。说明了当假设子梁转角相同时，组合梁的剪切刚度被高估，造成计算所得频率仍然高于实际频率。但是该计算模型中忽略了考虑转动惯量的影响，且无法分析沿梁轴向变抗弯刚度的钢-混组合梁的动力性能[13]。

动力刚度矩阵法在无滑移的组合梁结构分析中得到了广泛的应用[18-20]。该方法推导过程中没有引入力或位移形函数，因此是精确解而非近似解。由于可以采用与静力刚度法类似的方式组装动力刚度矩阵，因此，可被用于分析沿梁轴向变抗弯刚度的组合梁的动力特性。目前在考虑剪切滑移的组合梁动力分析方面，Li 等[21]、Bao 等[22]和Sun 等[23]学者给出了基于Euler-Bernoulli 梁理论的钢-混组合梁的动力刚度矩阵。如前所述，其无法满足高阶振型和短粗梁的频率分析要求。Wang 等[24]基于Timoshenko 梁理论及子梁等转角假设给出了6 个自由度的钢-混组合梁的动力刚度矩阵。子梁等转角假设过高的估计了结构的剪切刚度，造成频率计算结果仍然偏大。

本文基于Timoshenko 梁理论提出了用于分析钢-混组合梁动力特性的动力刚度矩阵法。该方法中，假设混凝土子梁和钢梁具有独立的剪切角，其大小与各自的剪切模量相关。该假设更加符合钢-混组合梁的实际运动状态，从而具有更高的计算精度。而且由于矩阵推导过程中没有使用力场或位移场近似形函数，因此计算结果是准确的。最后，采用本文提出的方法分析了不同跨高比、剪力键刚度时，剪切变形和转动惯量对钢-混组合梁自振频率的影响。

1 运动微分方程

1.1 基本假设

图1 钢-混组合梁构造图Fig. 1 Structural drawing of steel-concrete composite beam

本文中基于Timoshenko 梁理论分析钢-混组合梁的动力特性时，基本假设如下：

1) 本文只研究钢-混组合梁平面内的弯曲振动，且基于小变形假设。

2) 混凝土子梁和钢梁的剪切角是独立的。

3) 混凝土子梁和钢梁之间可以沿水平向相对滑动，而不能竖向掀起脱离。

4) 钢混结合面上的剪力全部由剪力键承担，且剪切滑移量与剪力键承受的剪力成正比关系，剪力键可以等效为连续分布的弹簧，剪力键的等效剪切刚度为K。

1.2 建立运动方程

根据以上假设，混凝土子梁与钢梁之间的剪切滑移量ucs(图2)可以表示为：

图2 剪切滑移量示意图Fig. 2 Schematic diagram of shear slip

式中：uc、us分别为混凝土子梁和钢梁中性轴处的轴向位移；θc、θs分别为混凝土子梁和钢梁的截面转角。

由Hamilton 原理，钢-混组合梁的运动学问题可以表示为如下：

式中，T、U和Ucs分别为结构的动能、应变能和界面剪切滑移势能。

考虑转动惯量的影响，因此，动能T为：

式中，w为结构的竖向位移。应变能U为：

界面剪切滑移势能Ucs为：

把式(3)～式(5)代入式(2)，并进行变分，可得钢-混组合梁的运动方程式：

式中，W、Ucs、Θs和Θc为相应的振型函数；ω为结构的自振频率；φ为初始相位角。

把式(8)代入式(6)，可得：

联立式(9)中的四个式子，可以得到组合梁运动微分方程的最终形式：

式中：

2 动力刚度矩阵法

式(10)为8 阶常微分方程，其特征方程为：

式(11)是一个一元四次方程，由费拉里法可以获得方程的四个根，如下所示：

式中，Ai、Bi、Ci和Di为振型待定系数。把式(14)代入式(9)可得待定系数之间的关系如下：

式中：

式中：

由式(14)可得，位移向量ue与待定系数向量a之间的关系，如下：

把式(14)～式(16)代入式(7)，可得力向量pe与待定系数向量a的关系式，如下：

式中：

由式(17)和式(18)可得：

在中国古代两件事最为重要，一是祭祀，一是战争。祭祀排在战争的前面，因为只有祭祀才能获得神灵的佑助，而能不能获得神灵的佑助是部落能不能得到发达的根本。从文献资料，中国远在尧舜禹三帝之时，就已经有以乐舞祭神的形式了。而实物的佐证就是良渚出土了玉钺。良渚时代当是大早于禹的尧舜时代，有学者甚至认为，良渚部落就是大舜的部落。从良渚反山发掘的十几座古墓来看，每墓出土玉钺只有一件，可见它极为珍贵。

式中，Ke=Me即为钢-混组合梁的8 自由度动力刚度矩阵。

按照与静力刚度法相同的过程，组装动力刚度矩阵，构建整体刚度矩阵Kg，则有：

式中，pg和ug为整体力向量和整体位移向量。

工程中常见的边界条件有：自由(F)、简支(S)和固支(C)三种，对应的边界条件为：

图3 给出了计算组合梁自振特性时的流程图。

图3 动力刚度矩阵法计算流程图Fig. 3 Calculation flow chart of dynamic stiffness matrix method

3 算例分析

3.1 理论验证

采用文献[13]中的钢-混组合实验梁2 对本文中的理论进行验证，结构尺寸见图4 所示。实验梁2 的横截面是由1700 mm×300 mm 的矩形混凝土子梁和550 mm×450 mm×28 mm 的工字钢梁组成，两者之间采用直径为22 mm 的剪力钉连接。梁长为8.5 m，跨径为8 m。沿梁长方向按照剪力键的疏密程度，划分为5 个区段[13]：(2.0+3×1.5+2.0) m；相应的剪力键刚度K为：(2204+1291+461+1291+2204) MPa，梁材料见表1 所示。

表1 梁材料参数Table 1 Material parameters of composite beam

图4 实验梁构造图 /mmFig. 4 Structural drawing of experimental beam

具体试验过程和ANSYS 有限元软件建模过程均已在文献[13]中详述，此处不在赘述。

图5 给出了ANSYS 计算结果的前5 阶振型云图如下所述。

图5 前5 阶振型Fig. 5 First five modes

由表2 可得：

表2 实验梁2 自振频率分析结果对比表Table 2 Comparison of eigenfrequencies obtained by different methods for experimental beam 2

1) 试验测试、ANSYS 计算和本文理论计算的前3 阶频率基本一致。说明文中理论可用于分析钢-混组合梁的自振特性。

2) Euler-Bernoulli 组合梁模型和子梁同剪切角假设的Timoshenko 组合梁模型的频率计算结果明显大于试验测试、ANSYS 计算和本文理论，而且阶数越高，误差越明显。第5 阶频率Euler-Bernoulli组合梁模型误差达48.2%，子梁同剪切角假设的Timoshenko 组合梁模型达22.6%。说明钢-混组合梁自振特性分析时，不可忽略剪切变形的影响；且不可简单的假设混凝土子梁与钢梁剪切角相等。

3.2 影响因素分析

本节在讨论不同剪力键刚度、跨高比时，剪切变形和转动惯量对钢-混组合梁自振频率的影响。以图3 中的钢-混组合梁为参考，改变部分结构或材料参数，讨论相关影响因素，最终得到跨高比为何时，可以忽略剪切变形和转动惯量的影响。工程应用中，前3 阶竖向频率比较重要，因此本节主要针对前3 阶竖向频率进行讨论。

3.2.1 剪力键刚度分析

图6 可以明显得出：

图6 前3 阶自振频率随剪力键刚度的变化情况Fig. 6 Effect of shear connector stiffness on first three eigenfrequencies

1) 相对于本文模型计算结果，频率阶数越高，α 越小(剪力键刚度越大)，Euler-Bernoulli 组合梁模型结果误差越大。

2) 当α＜10-3(剪力键刚度很大，剪切滑移量很小)或者α＞1(剪力键刚度很小，上下层间几乎无剪切约束)时，两种计算模型的计算结果误差基本不变。

综上所述，在以下讨论跨高比的影响时，选取组合连接系数α=10-3，并选取第3 阶为控制频率进行分析。

3.2.2 跨高比的影响

图7 给出了钢-混组合梁前3 阶频率随跨高比的变化规律。图7 中，ωE-B为采用Euler-Bernoulli组合梁模型(文献[23])计算所得频率；ωT为采用本文计算模型(Timoshenko 组合梁模型)所得频率。

图7 表明：

图7 前3 阶自振频率随跨高比的变化情况Fig. 7 Effect of depth-to-span on first three eigenfrequencies

1) 随着跨高比的增大，Euler-Bernoulli 梁理论的计算误差逐渐减小。且跨高比越小，这种变化越明显。

2) 跨高比对高阶频率的影响明显大于低阶频率。对于1 阶频率，跨高比大于10 后，两种理论计算结果误差即小于5%；对于2 阶频率，跨高比应大于18；对于3 阶频率，应大于25。

4 结论

本文基于Timoshenko 梁理论，提出了一种分析钢-混组合梁动力性能的新方法—动力刚度矩阵法。通过实验梁模型对方法的适用性进行了验证，并讨论了高跨比、剪力键刚度对频率计算结果的影响，主要结论如下：

(1) 本文提出的动力刚度矩阵法适用于分析钢-混组合梁的动力特性。相比于Euler-Bernoulli 组合梁计算模型和子梁同剪切角假设的Timoshenko 组合梁计算模型，本文的方法具有更高的计算精度。

(2) 钢-混组合梁动力特性分析时，尤其是进行高阶频率分析时，不可忽略剪切变形和转动惯量的影响；而且不能简单的假设混凝土子梁与钢梁具有相同的剪切角。

(3) 相对于Timoshenko 组合梁理论，分析的频率越高或剪力键刚度越大，Euler-Bernoulli 组合梁模型计算结果误差越大。

(4) 跨高比越大，Euler-Bernoulli 梁理论计算结果误差越小。若控制误差在5%以内，则对于第1 阶频率，跨高比需大于10；第2 阶，需大于18；第3 阶，需大于25。