斜拉桥塔-索-梁耦合面内整体动力学模型与1:1内共振影响性分析1)

陈柯帆 李 源,†, 贺拴海,† 宋一凡,† 徐珂瑶

* (长安大学公路学院,西安 710064)

† (长安大学旧桥检测与加固技术交通行业重点实验室,西安 710064)

引言

斜拉桥整体刚度低,受载荷作用产生全桥振动时,将同时存在以拉索振动模态为主导的结构局部模态和以主塔或主梁振动模态为主导的结构整体模态[1].研究表明,随着拉索参数变化,当结构局部模态频率靠近某一阶结构整体模态频率时,系统将出现明显的局部-整体模态耦合共振,此时系统能量在拉索与其他构件间周期性转换,进而产生复杂的非线性内共振[2-5].已有研究表明斜拉桥普遍存在此类特征的全桥共振[6],有的甚至已经严重影响到桥梁结构的安全性能与使用寿命[7],引起了国内外学者的广泛关注.

在部分非线性共振研究中,通过把塔、索和梁等构件等效为质量块[8],或将主梁对拉索的激励作用等效为简谐载荷[9-12],认为拉索自身模态与激励作用模态间相互独立,彼此间无影响作用.然而在斜拉桥这类塔-索-梁耦合结构的内共振中,对拉索模态的激励作用主要来源于主梁或主塔的振动模态,由于索-梁端、索-塔端存在动态协调关系,拉索的振动模态与主梁或主塔振动模态间存在相互耦合的关系,两者模态间实时变化且相互影响,其内共振机理非常复杂[13].因此,为了能更贴近工程实际结构,Fujino 等[14-15]和Gattulli 等[16-17]从试验和理论角度出发,研究并建立了考虑局部-整体模态间相互作用的索-梁耦合动力学模型,提出了振型局部化系数的计算方法,并观察到了由于整体模态和局部模态线性耦合的模态失真现象.康厚军等从连续介质力学角度出发,建立了索-梁[18-19]、索-拱[20-21]、索-塔[22]、索-曲线梁[23]等动力学模型,基于弹性体与柔性体间的非线性关系,得到了各自约化动力学控制方程,并就其非线性动力学行为进行了系列的理论与试验研究.

关于斜拉桥内共振的现有研究绝大多数以索-梁耦合结构为主[1-19,23-24],其通过假设拉索上端固结而忽略了主塔振动模态对结构整体模态的耦合作用,此时与拉索振动模态产生耦合关系的整体模态完全以主梁模态为主导.然而在目前实桥或试验观测到,全桥共振时斜拉索、主梁和主塔均发生了大幅振动,且后两者模态的参与将改变斜拉索的某些非线性动力特性[5,13,25].极少数考虑了主塔激励的研究或将其视为理想激励[26-28],或通过模态截断仅考虑了主塔的一阶模态[29],这些研究模型及结果显然不能全面揭示斜拉桥的动力学行为.因此,考虑塔-梁模态耦合作用建立的塔-索-梁耦合结构更契合工程实际结构体系,其研究结果才能更准确反映斜拉桥内共振机理.

本文基于参数质量离散方法,建立了考虑塔-梁模态耦合作用的塔-索-梁耦合动力学模型,通过有限差分法将参数体系的动力平衡微分方程转换为代数方程,采用模态拖拽法得到了结构面内运动方程和模态函数,并与数值模拟结果进行了对比验证.在此基础上,采用4～5 阶Runge-Kutta 方法编写了振动方程的数值仿真程序,重点分析和讨论了塔-梁间模态耦合作用对结构整体模态和结构1:1 内共振的影响.

1 斜拉桥面内动力学模型

忽略主塔和主梁轴向振动位移,基于参数质量体系离散方法[2,30-31],将斜拉桥主塔及主梁分别按照db,dp等间距划分为参数质量体系,以此建立固结体系斜拉桥的面内离散动力学模型如图1 所示.

图1 离散化的斜拉桥面内整体动力学模型Fig.1 The in-plane full bridge dynamic model of a cable-stayed bridge

本文约定下标“p,c,b”分别代表塔、索和梁的相关参数,其中r和j分别表示对塔段和梁段的计数(r∈[1,R],j∈[1,J]),i表示拉索的计数(i∈[1,I]),定义Ci#和Bj#分别表示小里程至大里程方向的第i根拉索和第j个梁段,Pr#表示竖直由上至下方向的第r个塔段,而Pri#和Bji#对应拉索锚固处的塔段及梁段;θci表示拉索与主梁大里程方向夹角.需要说明的是,为方便计算,将P0#的质量归于P1#,其余边界质点无位移,故忽略其质量对结构整体运动的影响.对Ci#及其连接的Pri#和Bji#开展受力分析如图2 所示.

图2 Ci#振动子系统受力分析示意图Fig.2 The forced diagram of the Ci# dynamic subsystem

图中zpr,xci和xbj分别表示Pr#,Ci#和Bj#在各自轴向的坐标;mpr和mbj分别表示Pr#和Bj#的质量,另有mci表示Ci#的单位长度质量;wci,vci和uci分别简写自wci(xci),vci(xci,t) 和uci(xci,t),其中wci为拉索的静平衡线形,vci为拉索的横向振动位移,uci为拉索的轴向振动位移;s0ci为Ci#静平衡长度,sci为Ci#振动弧长;vpr和vbj分别简写自vpr(zpr,t),vbj(xbj,t),分别表示Pr#的横向振动位移及Bj#的竖向振动位移.βb(j-1,j)表示Bj-1#和Bj#面内运动的转角,Fb(j-1,j)与Nb(j-1,j)分别表示为Bj-1#与Bj#间的剪力与轴力,同理如βp(r-1,r),Fp(r-1,r),Np(r-1,r)之于Pr#.试验表明[32],张紧弦的低阶模态是自由振动时的主导模态.若仅考虑拉索的低阶模态,则其轴向惯性效应可以忽略(准静态近似),此时拉索的轴向动力学将通过静态缩聚的形式被其横向动力学所控制[33].在此基础上,假设拉索质量沿弦向均匀分布,建立拉索面内静平衡方程及静平衡附近的横向非线性振动方程分别如式(1)和式(2)所示

式中,g表示重力加速度,取9.806 m/s2;Tci和τci分别表示Ci#在切向的初始索力与索力动增量,其与弦向的初始索力(Hci)与索力动增量间(hci)存在近似关系[7,19]

拉索横向运动的单位动应变表达式为

对式(4)沿xci方向积分,代入式(3)并简化后,可以得到拉索弦向的索力动增量表达式为

式中,lci表示Ci#上下端锚固点的弦向距离,Lci表示Ci#静平衡长度;Uci表示Ci#轴向振动伸长量.fci为Ci#的垂度,表示静平衡状态下的索中点与弦向中点的横向距离[19,24,31].Ci#轴向振动伸长量由其上下边界处连接的Pri#和Bji#运动分量构成

为了使对称拉索的两端边界条件符号一致,认为塔段的右侧位移、梁段的向上位移、拉索沿轴线顺时针横向振动位移为正.基于此,拉索振动时边界条件为

根据牛顿定律和图2 的受力分析,并考虑到狄拉克函数性质(δ),分别建立任意Pr#横向及Bj#竖向动力平衡方程

式中,参数上标“·”表示对时间t求偏导;ηpk和ηbk分别表征固结体系斜拉桥中主塔及主梁运动通过索-梁固结点相互影响作用,下标“k”表示与该点相关参数,下同.对于任一塔段而言,Pr#左右侧轴力表达式为

对于任一梁段而言,Bj#左右侧轴力表达式为

任意相邻的塔段或梁段间存在平衡关系

式中,参数上标“'”表示对轴向坐标求偏导;Mpr和Mbj分别表示Pr#和Bj#处的弯矩;Epr和Ebj分别表示Pr#和Bj#的弹性模量;Ipr表示Pr#的面内横向弯曲惯性矩;Ibj表示Bj#的面内竖向弯曲惯性矩.假设质量体系分布较密,相邻梁段间的相对位移较小,采用差分法将对位移的偏微分多项式转换为差分代数方程[2,31]

对于塔-梁固结斜拉桥,索-梁固结点处离散段同时为桥塔和主梁的一部分,塔与梁振动模态间存在相互影响,且通过该点进行传递.截取索-梁节点受力示意图如图3 所示.

图3 塔-梁固结点处的受力分析示意图Fig.3 The forced diagram of the pylon-beam connection joint

如图3 所示,根据离散模型的假设,该点处横向及纵向振动位移为零,而其水平、竖向的剪力与轴力表达式间分别存在相互作用

图1 中,主塔上端自由、下端固结,主梁左右侧均为简支端,其边界条件分别为

式中,e+,e-分别表示大小里程边界.考虑边界条件式(27)～式(30)后,整理式(15)～式(26)并分别代入式(13) 和式(14),可以得到面内的Pr#横向及Bj#竖向运动方程.

2 结构整体模态定义与模态分析

2.1 运动方程与模态函数

为便于找到系统的内共振形式,采用模态拖拽法定义第n阶斜拉桥塔、索、梁自由振动位移表达式[31,34]

式中,加粗符号表示矩阵或向量(下同),基本形式列于附录A.式(31)是一个含有多维向量(塔/梁)和连续函数(索) 的混合表达式,其中,元素为φpr(n),φbj(n)的向量分别表示塔、梁在第n阶面内振动模态下的振动模态振型向量;qnpr简写自qpr(t)(n),表示与时间相关的主塔第n阶振动模态形状变化因子,同理如qnci与qnbj;元素为fpi(xci),fbi(xci)表示与边界条件式(10)和式(12)关联的模态拖拽函数向量,定义其基本形式为[34]

对于图1 所示的主塔与主梁离散参数质量体系,Pr#和Bj#的形状变化因子已包括了各自各阶振动模态下的振型函数.因此,Pr#和Bj#的振动位移表达式可以简写为

式中,Apr(n)(zpr)表示Pr#在结构第n阶振动模态下的振幅,同理如Aci(n)(xci)和Abj(n)(xbj).为简化表达,选取三角函数作为式(31)中拉索振动方程的振型基函数[19,26-27]

基于此,可以得到Ci#横向振动方程为

将上式连同式(34)和式(35)分别代入整合后的式(2)、式(13)和式(14),使用Galerkin 方法进行模态截断可得图1 所示固结斜拉桥的面内整体运动方程.考虑到张紧弦的低阶模态是其自由振动时的主导模态,为简化计算,在此仅考虑了拉索一阶振动模态,并将qnci简写为qci,得到结构方程如下所示

式中,对角矩阵Ωpr,Ωnci和Ωbj是主塔、拉索与主梁振动模态的特征对角矩阵,其主元参数量纲与频率相同.其中,Ωpr和Ωbj表达式为

式中,矩阵S,D,C,G分别表征了剪力效应、轴力效应、拉索的弹性支承效应与重力效应对Pr#或Bj#振动模态的影响,各参数表达式详见附录B.从式(38)～式(40)可以看出,Pr#,Ci#和Bj#的运动方程存在耦合项,表明离散的塔段或梁段运动并不彼此独立,其通过剪力与轴力效应相互影响,共同参与结构的整体运动.此外,线性化后的式(38)～式(40)实质是关于结构固有振动模态频率ωN的K阶齐次超越方程

式中,下标“N” 的参数表示结构面内固有振动模态参数;EN表示结构整体运动的特征矩阵,QN表示结构形状变化因子向量,其表达式如下所示

式(44) 中,除塔、索和梁特征矩阵外,存在Πbj,3等耦合项系数,表征了塔、索、梁等构件振动模态在结构振动模态下的耦合作用.式(43)存在非零解的前提须EN行列式为零.因此,对EN进行特征值求解可得结构的振动模态参数

式中,ωN(K)表示结构前K阶固有振动模态频率.其中第n阶频率ωN(n)对应的特征向量AN(n)表示该阶振动模态下的塔、索和梁构件振动模态的振幅向量.将结构第n阶固有振动模态频率ωN(n)和对应的特征向量AN(n)代入式(34)、式(35)和式(37)后,可以得到塔、索、梁的振动位移表达式.由于本文仅考虑了拉索的一阶频率,则式(46)的K个解中包含了C1#～C4#的1 阶局部模态频率,其余为结构的整体振动模态频率,此时

若将主梁和拉索无限细分,且使用拉索的前N阶模态进行截断,则式(46)有N个解,此时K=N,按照ωN(k)数值大小升序排列即为结构的第1 阶至第N阶固有振动模态频率.

2.2 基于FEM 的斜拉桥整体模态验证

参考中国西北地区一座斜拉桥设计参数,建立一座单塔四索斜拉桥的动力学简化模型,构件参数如表1 和表2 所示.

表1 拉索参数Table 1 Parameters of cables

表2 主塔及主梁参数Table 2 Parameters of the pylon and beam

为辨别结构模态,引入模态局部化程度系数Λc(n),Λp(n)和Λb(n)[1,24-25]以定义构件的振动模态在第n阶整体模态下的模态参与度

基于上式获取的结构模态参数仅包含结构整体模态的特征值与特征向量.将本文方法(EMP)与有限元方法(FEM)得到的G1st～G5th整体模态振型及频率汇总如表3.

表3 两种方法得到的结构体系整体模态振型Table 3 Global modal shapes obtained via these two methods

表3 中,两种方法得到的结构整体模态振型几乎一致,表4 中模态频率误差随离散间距d取值减小而减小,这与式(21)～式(24)中有限差分法应用假设相关.采用本文方法求解复杂塔-索-梁结构体系模态参数,在误差可接受范围内确定离散间距d及结构特征矩阵(式(44))元素的取值,借助MATLAB,Excel 等工具即可高效率获取结构体系固有振型及频率,无需进行大量细化的有限元建模分析.本文后续研究为避免较大计算量,离散间距取值为4 m.

表4 两种方法得到的结构体系整体模态频率Table 4 Global modal frequencies obtained via these two methods

2.3 主塔振动模态对整体模态的影响

表3 显示,在本文塔-索-梁耦合模型中,结构的整体模态不再完全以主塔或主梁的振动模态为主,而是存在由两者共同参与的混合整体模态,如G1st,G4th等.为开展进一步研究,建立了3 个动力学模型,如表5 所示.

表5 3 种索-支撑简化动力学模型Table 5 Three dynamic models reduced from the cable-support structures

表5 中DM2#被广泛采用于索-梁模态耦合的相关非线性振动研究[2,15-18,26-27,31].采用本文方法对表5 中的3 个模型进行特征值求解,剥离其整体模态频率汇总如图4 所示.

图4 3 个结构的整体模态频率Fig.4 The global modal frequencies of these three structures

图4 中Bnth表示以梁第n阶局部模态为主的整体模态;Pnth表示以塔第n阶局部模态为主的整体模态.上图表明,DM1#的整体模态由塔和梁的振动模态组成.在数值大于10 Hz 的高阶模态中,DM1#整体模态频率数值上分别与DM2#中梁的振动模态频率和DM3#中塔的振动模态频率相等,表明塔-梁振动模态间的耦合作用无法影响结构的高阶整体模态频率;而在低阶整体模态中,三者模型的模态频率数值上存在一定差异.为进一步研究,取图4 中阴影部分的G1st～G9th阶模态参数如图5 所示.

图5 G1st～G9th 结构的整体模态频率变化曲线Fig.5 The trends of the structural global modal frequencies of G1st～G9th

图5 中在G1st,G3rd,G4th,G5th,G7th和G8th阶等低阶反对称整体模态中,DM1#与其他两个模型得到的模态频率存在较大误差,且局部化程度系数显示主塔振动模态参与了以主梁振动模态为主导的整体模态运动,使得该阶整体模态由完全以主梁模态为主导的完全整体模态转为塔-梁振动模态耦合的混合整体模态.而对于低阶的对称整体模态(如G2nd,G6th)、高阶整体模态(G9th及以上),塔-梁模态耦合作用对结构整体模态没有明显影响,此时塔或梁的振动模态完全主导了该阶整体模态.这表明,在考虑了主塔振动模态以后,斜拉桥等塔-索-梁耦合结构的整体模态可以进一步被细分为完全整体模态和混合整体模态,其细分流程如图6所示.

图6 斜拉桥面内固有振动模态的细分流程Fig.6 The reduced process of the in-plane natural vibration modes of cable-stayed bridges

上述研究结果表明,主塔振动模态参与度决定了结构整体模态的动力特性,在斜拉桥这类塔-索-梁耦合模型中,考虑主塔振动模态对研究因局部-整体模态耦合激励产生的内共振现象具有重要意义.

3 斜拉桥1:1 内共振分析

3.1 模态翻转与1:1 内共振动力特性

为排除起振阈值影响以更清晰观测内共振响应变化规律[6],本文以“1:1”内共振为主要研究对象.在考虑拉索前2 阶振动模态的基础上,从纯数值角度变化斜拉索索力与垂度参数,汇总结构第5 阶固有振动频率(N5th)至结构第9 阶固有振动频率(N9th)变化规律如图7 所示.

图7 结构N5th～N9th 整体模态频率随C1#的垂度与索力变化曲线Fig.7 The trends of the global modal frequencies from N5th to N9th with the sag and force of C1#

图7 中,随着C1#垂度及索力参数变化,当结构局部模态频率临近某一阶结构整体模态频率时,两者模态振型发生了快速且连续的交换(Veering 现象[1]).此时局部-整体模态间发生耦合,结构能量在拉索与系统内其他构件间进行周期性地转换,容易激励产生1:1 内共振.研究表明,内共振完全由系统非线性特性决定[35].为进一步研究这种现象并验证本文结构的非线性动力特性,引入系数σci表示Ci#拉索局部频率(ωci)增量,Θci(n)表示Ci#拉索局部模态频率与结构第n阶整体模态频率(ωG(n))的靠近程度[24,31]

采用4～5 阶Runge-Kutta 积分方法对式(38)～式(40)进行数值仿真.为清晰展示系统构件在达到稳态后一个振动周期内的动力特性,提高仿真效率,本文的数值仿真过程设定固定步长为0.001 s,迭代时间为100 s.在塔顶(P1#)、梁小里程端点(B1#)附近同时设置了0.02 m 的初始位移,改变σci至分别与第G6th与G8th模态频率满足“1:1”比例关系.由图5 可知,G6th为以主梁模态为主导的整体模态,而G8th为以主塔模态为主导的混合模态.对获取的拉索振动响应曲线进行快速傅里叶分析(FFT),汇总拉索中点振动响应及对应频谱图如图8 所示.

图8 局部-整体模态耦合激励产生“1:1”内共振现象Fig.8 The phenomenon of the 1:1 internal resonance excited by the local-global modal interactions

图8(a) 与图8(b) 显示,C1#分别与G6th和G8th整体模态“1:1”耦合并产生了“拍”特性明显的内共振.这表明,拉索局部模态与HGM 和EGM (定义见图6)两者“1:1”耦合都可以激励拉索产生“1:1”内共振.而这样“拍”特性明显的内共振现象同样验证了本文推导方程及数值仿真方法可以有效进行非线性动力响应分析.

此外,对比图8(a)与图8(b)后可得,在同等的初始条件下,由两种不同主导模态激励产生的内共振响应的“拍”频及幅值不同,表明拉索的内共振动力特性随主导模态变化而变化.实际上在塔-索-梁耦合模型中,除局部-整体模态存在交换现象外,还存在主导模态间的交换现象.以表5 中DM1#和DM2#为例,变化主塔单位长度质量而使主梁构造参数保持不变,计算DM#中G7th模态参数如图9 所示.

图9 结构G7th 模态参数随主塔单位长度质量变化规律Fig.9 The trend of the G7th-order modal property varied with the mass per unit length of the pylon

上图显示,随着塔单位长度质量变化至6.72×103kg/m 附近,DM1#模型的G7th主导模态发生了快速且连续的交换,由以主梁模态为主导的HGM 变化为以主塔模态为主导的HGM,各构件的模态局部化程度系数也随之发生变化.此时,DM1#模型中主导模态的变化将进一步影响拉索内共振响应的动力特性,其变化规律与图8 中一致,在此不过多赘述.

3.2 HGM 耦合内共振的参与度分析

不同于以往EGM 耦合内共振的相关研究,在塔-索-梁耦合系统中,当某阶HGM 与某根拉索的局部模态“1:1”耦合时,主塔和主梁的振动模态将同时参与“1:1”内共振.为研究以主梁模态为主导的HGM内共振动力特性,设定了1#研究工况.

1#研究工况(RC1#): 仅考虑拉索1 阶模态,针对DM1 与DM2 模型,在主梁小里程桩号附近增加0.02 m 初始位移,改变σc1至分别与G7th,G9th模态频率满足“1:1”关系.

汇总RC1#工况条件下的C1#振动响应曲线如图10 所示.

图10 RC1#工况下拉索时程曲线Fig.10 The time histories of C1# under RC1#

图10(a)和图10(b)中,DM1#模型的G7th是以主梁振动模态为主的混合整体模态,主塔振动模态参与了整体模态运动,同时也影响了拉索的内共振动力特性.图10(c)和图10(d)中,G9th是完全以主梁振动模态为主的整体模态,此时DM1#与DM2#系统拉索的内共振时程曲线动力特性相同.为进一步明晰混合模态内共振时的能量传递路径,采用无相位滤波法分离RC1#中C1#及其相连的塔段(P1#)和梁段(B2#)在主共振频率附近的振动响应信号,汇总如表6所示.

表6 DM1#中发生G7th 和G9th 内共振时,C1#振动子系统构件在各自主共振频率附近的分离响应信号Table 6 The separated signals around the main resonance frequency in the dynamic sub-system of C1# when the internal resonance of G7th-order or G9th-order occurring in DM1#

表6 中,前25 s 的分离信号存在滤波延迟,以后75 s 的滤波信号为研究对象.当C1#局部模态与G7th混合模态“1:1”耦合时,P1#,C1#和B2#在主共振频率附近的响应信号此起彼伏、互相耦合,表明共振时能量在三者间来回转换.当与G9th混合模态“1:1”耦合时,C1#与B2#的振动响应曲线互相耦合,而P1#与C1#的振动特性相同且其响应幅值与B2#存在数量级差异,表明P1#的振动响应仅来源于塔-索端点耦合而未参与内共振,此时系统能量在拉索与主梁间来回转换.采用上述同样的滤波方法,获取在各自主共振频率下运行时刻t1的结构振型如图11 所示.

图11 DM1#和DM2#模型在RC1#工况下主共振频率对应的结构振动形态Fig.11 The structural dynamic configuration corresponding to the main resonance frequency in DM1# and DM2# under RC1#

图11(a) 显示,由主塔第3 阶面内模态和主梁第5 阶面内竖向模态共同构成的混合整体模态激发了塔-索-梁耦合内共振,此时主塔或主梁振动模态参与度系数Λp和Λb与原G7th的局部化程度系数误差较小.而在图11(b)～图11(d)中,由主梁的第5 阶或第6 阶面内竖向模态为主的整体模态激发了索-梁耦合内共振,拉索与主梁间存在能量交换作用,而主塔未参与整体运动.为研究主塔振动模态对内共振的影响,设定了2#研究工况.

2#研究工况(RC2#): 针对DM1 与DM3 模型,在塔顶增加0.02 m 初始位移,改变σc1至分别与G8th,G11th模态频率满足“1:1”关系.

为简化表达,汇总RC2#工况条件下的C1#振动响应曲线,及其在各自主共振频率下运行时刻t1的结构振动形态分别如图12 所示.

图12 RC2#工况下C1#时程曲线与对应的结构振动形态Fig.12 The time histories of C1# and the corresponding the structural dynamic configuration under RC2#

图12 结果再次验证了在HGM 中,非主导的主梁振动模态将参与由主塔振动模态主导的内共振能量转换,并将改变拉索内共振响应曲线的“拍”频及振幅,此时各构件的整体模态局部化程度系数是重要的量化因素.

3.3 HGM 耦合内共振的影响性分析

为研究非主导的主梁振动模态在混合整体模态内共振下的影响效应,设定了3#研究工况.

3#研究工况(RC3#): 针对DM1 与DM2 模型,在主梁小里程桩号附近增加0.02 m 初始位移,分别改变σc1,σc2,σc3和σc4至与以主梁振动模态为主导的整体模态(G1th,G2th,G4th,G6th,G7th,G9th,G10th)频率满足1:1 比值关系.

拉索的振动响应曲线(Aci)能一定程度上反映内共振时拉索动能变化规律,定义分别表示在EGM 和HGM 模型中拉索振动响应的最大值.同时为尽量减少多自由度系统内复杂组合内共振引起的响应误差,采用无相位滤波法分离各子工况下主共振频率附近的拉索响应,汇总其响应振幅最大值曲线如图13 所示.

图13 DM1#和DM2#模型在RC3#工况下拉索主共振频率附近的振幅最大值Fig.13 The maximum response of cables in DM1# and DM2# under RC3#

图13 中,对于不同阶次的整体模态内共振,DM1#模型中的拉索共振响应最大值与DM2#的计算结果存在较大差异,表明非主导的主塔振动模态参与并影响了HGM 模态内共振,其影响效应随阶次变化而不同.其中,对称布置的C1#与C4#(或C2#和C3#)基础构件参数、锚固位置对于整体模态振型有效质量相同,因此对称索的各阶振幅最大值曲线趋势、误差数值基本一致.对于同一侧的C1#与C2#(或C3#与C4#)各阶振幅最大值曲线变化趋势相似度较高但差值不同,表明各阶非主导主塔模态对于任意拉索内共振的影响效应类型一致,但影响效应大小随拉索变化而变化.对于两种模型计算C1#(图13-A点) 或C4#(图13-B点) 结果差值较大的G6th内共振,结合对应的响应与频谱图可得,DM1#中较长的C1#或C4#在整体模态激励作用下发生了组合内共振,故而其主共振频率附近的模态响应弱于DM2#计算结果.为研究非主导的主塔振动模态在混合整体模态内共振下的影响效应,设定了4#研究工况.

4#研究工况(RC4#): 针对DM1 与DM4 模型,在塔顶增加0.02 m 初始位移,分别改变σc1,σc2,σc3和σc4至与以主塔振动模态为主导的整体模态(G3rd,G5th,G8th,G11th)频率满足1:1 比值关系.

同样采用无相位滤波法分离各子工况下主共振频率附近的拉索响应,汇总其响应振幅最大值曲线如图14 所示.

图14 结果与图13 结果规律基本一致,即非主导模态的影响效应随阶次变化而变化,且基本以DM3#模型的计算结果较大,表明在相同初始条件下,非主导构件模态的参与会减弱拉索的内共振响应,且对于各阶的影响效应也随拉索锚固位置与基础参数的不同而不同.为了能更进一步明晰HGM内共振模式的影响效应,定义ηci(n)表示Ci#在第n阶EGM(AciE)和HGM(AciH)两者模型中振动响应最大值的误差率,如下

ηci(n)＞ 0 表明非主导模态的影响效应为抑制作用,反之为激励作用.为了对比EGM 和HGM 两种模式内共振下的影响效应,汇总C1#和C2#在G1st～G11th中两种模式的响应误差率如图15 所示.

图15 同侧拉索在两类GM 下振动响应最大值的误差率Fig.15 The error coefficient of same-side cables’ maximum response in these two GMs

图15 中,除15(a)中的G6th外,EGM 模态下的影响系数数值上接近零,表明EGM 模态下的内共振响应不受非主导构件的振动模态影响,此时拉索模态与主导的主梁或主塔模态“1:1”耦合并激励产生剧烈内共振.而在HGM 模态下,两索的影响系数变化趋势基本相同,且除G1st外基本表现为抑制作用,此时非主导模态参与了内共振并分走了一部分系统能量,所以拉索的内共振响应幅值随之减小.

4 结论

本文基于离散的参数质量体系,考虑了塔-索、索-梁间的动态关系与斜拉索初始垂度、大位移振动引起的几何非线性影响,建立了新的塔-索-梁耦合面内整体动力学模型.通过有限差分法代数转换了参数体系动力平衡方程的偏微分多项式,采用模态拖拽法获得了结构运动方程与模态函数,精细化分析了塔-梁模态耦合作用对整体模态的影响,首次提出了完全整体模态和混合整体模态两种模式,考虑拉索一阶模态并通过数值仿真开展了混合整体模态下的1:1 内共振影响性分析,得到结论如下.

(1)塔-梁模态耦合作用对本文结构的低阶对称及高阶整体模态无影响作用,此时结构的整体模态完全由主塔或主梁振动模态主导,为完全整体模态;对结构的低阶反对称整体模态具有显著影响,此时主塔以及主梁振动模态共同参与了该阶整体模态,为混合整体模态.

(2)拉索局部模态与两种整体模态“1:1”耦合都将激励产生内共振,而主塔主导或主梁主导的两种拉索内共振响应动力特性不同.在同一阶混合整体模态下,当主导模态的局部化程度系数靠近非主导模态局部化程度系数时,该阶整体模态的主导模态将出现快速且连续的交换,而拉索振动响应的动力特性也将随之改变.

(3)构件的模态局部化程度系数是两种模式内共振参与度的重要量化因素.在完全整体模态内共振中,非主导模态局部化程度系数接近0,因此系统能量转换仅发生在拉索与主导模态间,非主导模态不参与内共振;而在混合整体模态内共振中,系统能量转换发生在塔、索和梁局部模态间,主导与非主导模态的结构振型参与度与各自局部化程度系数一致.

(4)针对本文固结体系四索结构的内共振分析发现,由于各阶整体模态能量总数恒定,混合整体模态内共振下非主导的结构模态参与了系统能量转换并将改变拉索振动响应的动力特性,其影响效应随整体模态阶次与拉索构件参数变化而变化.

(5)本文建立的塔-索-梁耦合动力学模型为深入开展斜拉桥大系统非线性动力特性研究提供了一种可靠途径,下一步将结合构件几何参数、各阶整体模态频率和振型以更准确研究和定义模态能量,进而深入开展斜拉桥内共振的参数敏感性分析.

附录A

式(31)中,vpr,vci,vbj,φpr(n),φci(n),φbj(n),qpr,qpr,qpr,fpi(xci),fpi(xci),fpi(xci)等皆为形式相同的列向量,为避免赘述,在此仅展示vpr

依据式(1)和拉索边界条件,可得拉索的静平衡线形为二次抛物线线形,其表达式为[19,24,26-27]

附录B

式(38)～式(40)中,q相关参数表示为相同形式的R维、I维、J维的列向量,为避免赘述,仅取如下所示

除 Πpr,1,Πci,1和 Πbj,1外,其余 Π 相关参数矩阵表示为相同形式的R阶、I阶、J阶的对角矩阵.为避免赘述,仅取Πpr,1,Πci,1和 Πbj,1如下所示

Ωci为拉索的局部模态对角矩阵,其主元形式为

相关参数矩阵的主元具体形式如下所示

为简化表达,定义Ppm,Ppm,n,γpi,γbi,μpi和μbi表达式分别为

Spr与Sbj,Dpr与Dbj,Cpr与Cbj形式相同,仅角标(p,r或b,j)不同,限于篇幅,在此仅展示Spr,Dpr,Cpr为

塔的重力影响矩阵Gpr为