桥梁频率不仅是反映桥梁整体结构完整性的重要指标,也是桥梁动力特性的重要参数,对桥梁承载力的评估具有重要参考价值[1-2]。现行测量频率方法应用最多的是桥梁动载试验,通过在桥梁上安装多组加速度传感器,利用有线或无线采集系统记录桥梁的动力响应,该方法需要中断交通且繁琐费力[3]。
根据车桥耦合理论和间接测量法原理进行频率识别是一种有效的识别手段,车辆驶过桥梁时,车辆与桥梁相互作用和影响的问题为车桥耦合问题。其研究自1825年人类有了第一条铁路就已开始[4-5],在20世纪60年代以前,主要集中在简单移动荷载作用下铁路桥梁的车桥耦合动力分析研究上,20世纪60年代以后,随着高等级公路桥梁的大量修建,公路桥梁的车桥耦合问题越来越吸引桥梁界科研人员的关注,并且随着计算机的出现和有限元软件的应用,使得建立复杂的车桥耦合振动模型成为可能。Yang等[6]首次提出基于动力测试间接测量桥梁频率的方法,用单自由度车模型模拟检测车辆,检测车辆通过桥梁时,其垂向的惯性力和弹簧力经车轮以接触力形式传递到桥梁造成桥梁结构的振动,由于耦合作用,桥梁振动将再次影响到检测车的垂向振动,并为其主要振动源。因此,只要分析检测车垂直方向的响应,即可得出桥梁特征参数信息,从而规避了车桥耦合时变系统有载频率变化的限制[7]。其主要贡献是利用单轴四分之一车辆通过某一简支梁桥推导出了解析解,该公式结果与Yang等[8]提出的利用车辆与桥梁耦合单元计算出的数值解吻合,理论验证了车辆竖向加速度解析解中包含了桥梁频率的基本信息。Lin等[9]通过现场试验验证了该方法的有效性,说明了间接测量法识别桥梁频率的可行性。在此基础上,陈上有等[10]用车轮加弹簧-阻尼器-簧上质量模拟移动小车,根据车辆的加速度响应的频谱分析,成功地识别出简支梁桥的一阶频率;王希等[11]进行了在桥梁平整状态下由车辆振动信息识别简支梁桥固有频率的研究,可以准确识别桥梁的第一阶基频;Malekjafarian等[12]利用列车不同车厢的振动信息识别出了简支梁桥的基频; Kong等[13]利用单自由度小车识别出某简支梁桥的前三阶频率;Yang等[14]在桥梁间接测量法的回顾文章中阐明了间接测量法应用的注意事项以及单自由度车辆模型的局限性。
前人的研究多集中于利用简支梁桥和单自由度车辆模型进行数值模拟分析,且主要是桥梁基频的识别。笔者选用更加贴近实际的双轴半车车辆平面模型,并选用实际工程中的三跨斜交梁桥建模,运用间接测量理论识别出了桥梁的前3阶频率。在20世纪60年代以前,车桥耦合多是现场试验研究,总结出经验公式或理论来指导桥梁设计;60年代以后,车桥动力相互作用的研究方法多采用试验与理论相结合的方法,即用试验的方法校核理论的正确性,用校核过的、正确的理论方法进行车桥耦合振动的仿真分析,研究结构各参数对振动的影响。笔者主要运用后者方法,基于车桥耦合理论与间接测量法原理,通过有限元数值模拟方法得到理论解与试验数值,其基本流程:1)利用Ansys有限元软件建立桥梁模型,通过模态分析模块得到桥梁前3阶频率的理论解;2)建立车桥耦合模型,利用Ansys有限元软件的瞬态动力分析模块,模拟车辆驶过桥梁的过程,通过采用不同的参数得到车体的加速度时程后利用峰值拾取法等手段得到识别频率;3)将这两步得到的频率进行误差对比和验证分析。桥梁间接测量法不同于传统的桥梁频率测量方法,将加速度传感器放置于测试车辆上,不用在桥上架设各种仪器,且桥上任意移动车辆对频率识别无不利影响[15],因此具有较大的应用前景。
依据间接测量法理论[6],单轴四分之一测量车辆模型匀速驶过等截面简支梁桥,量测车质量为mv、刚度为kv、车辆移动速度为v。桥梁总长度为L,单位长度质量为m,截面抗弯刚度为EI,车桥振动简图如图 1所示,桥梁因车辆作用导致振动产生坐标x处挠度为u(x)。
采用有限元方法将桥梁离散化,桥梁振动方程可以表示为
车辆的运动方程为
经推导得到车辆竖向加速度解析解为
式中:
式中:${{\ddot q}_{\rm{v}}}(t)$为在时间t时刻车辆的竖向加速度;Δst, n为在车辆作用下桥梁第n阶模态产生的静位移;Sn为无量纲速度参量;ωb, n为桥梁的第n阶自振频率;ωv为车辆的竖向振动频率。
车辆竖向加速度时程响应主要由4种类型的频率余弦波所叠加组成[10],分别为车辆的固有频率ωv、驱车频率2nπv/L、桥梁的左移频率ωb, n-nπv/L和右移频率ωb, n+nπv/L。然后通过快速傅里叶变换等频谱分析技术,桥梁的n阶频率信息可以从车辆的振动加速度响应中提取出来。
运用ANSYS中的瞬态动力学分析方法进行车桥耦合振动数值模拟计算分析[16],利用耦合命令约束车轮与相关桥面点的竖向位移,编制APDL车桥耦合振动命令流计算车辆和桥梁的动态响应,然后利用车辆响应识别出桥梁前3阶自振频率。
工程实例为跨径组合为28 m+45 m+28 m的混凝土连续梁桥,斜交角为30°,桥面宽度为8 m,上部结构为单箱单室主梁。全桥主梁采用C50预应力混凝土,弹性模量为EC=3.45×104 MPa,为模拟桥梁的倾斜支承边界条件,桥梁中横梁和端横梁均模拟为刚性横梁。桥梁横截面如图 2所示。
桥梁的相关结构参数如表 1所示。
应用beam4单元建立桥梁的有限元计算模型,桥梁有限元模型如图 3所示。
应用有限元程序计算得到桥梁前3阶理论频率如表 2所示。
车辆模型采用双轴半车模型[17],车辆模型示意如图 4所示。
车辆模型具体参数如表 3所示。
双轴半车模型相比于单轴1/4车辆模型,增加了车体质量点的转动刚度,因此,多了一个转动频率,利用ANSYS软件建立车辆模型进行频率分析得到车辆模型的频率为:车辆竖向振动频率ωv为0.75 Hz、转动频率ωθ为1.49 Hz。
车辆频率信息也可由式(4)、式(5)计算。
式中:Iv为车体点头刚度。ki为车辆前轴垂向刚度,kj为车辆后轴垂向刚度,mv为车辆重量。
基于车桥耦合理论建立车桥耦合有限元模型进行数值模拟分析,由于高阶桥梁频率对车辆动力响应贡献较小,间接测量法的识别效果较差[18],现对桥梁只进行前3阶桥梁频率识别并进行相应的识别精度误差分析。
车辆匀速驶过该3跨连续斜交梁桥,选取6个不同车速,分别是10、20、30、40、50、60 km/h。在车辆质量为4 t的条件下计算得到车辆加速度时程数据,然后进行快速傅里叶变换得到频谱图,识别结果的频谱如图 5所示,识别得到的桥梁频率数值如表 4所示,识别频率与理论计算频率的相对误差如表 5所示。
从图 5可以看出,频谱图中共有5种频率信息,前3种频率是已知的车辆转动频率1.49 Hz、驱车频率和车辆竖向振动频率0.75 Hz,且随着车速的增大,驱车频率越来越大。后3种频率信息分别对应前3阶桥频,可以看到桥频都为双峰(桥梁左移频率ωb, n-nπv/L、右移频率ωb, n+nπv/L),并且随着车速的增大,双峰间距越来越远,峰值愈加不明显,可以认为双峰的均值为对应的桥梁频率。
通过表 5误差分析可得,速度对桥梁基频识别结果基本无影响,对高阶识别较不利。速度在20 km/h及以下时,识别精度较高,高车速时无法有效识别桥梁的频率信息。
为了分析不同的车桥质量比mv/(mL)对桥梁频率识别精度的影响,在10 km/h的车速下,选取不同的车桥质量比进行频率识别,结果如图 6所示。
从图 6中可以看出,每条车辆加速度频谱曲线都有6个峰值,它们分别是车辆转动频率1.67 Hz,10 km/h时驱车频率2πv/L为0.08 Hz,车辆竖向振动频率ωv因车辆质量不相等而数值不同,车辆竖向频率随着车重增加而逐渐减少,这3个频率信息可以通过事先测试车辆而已知,并且数值相对于桥梁频率一般较小,所以可以较明显地区分出后面3个峰值为桥梁的前3阶频率信息。观察上述加速度时程数据频谱图,可以轻松地排除已知的驱车频率、车辆转动频率和车辆竖向频率,最后识别出桥梁前3阶频率信息。
从表 6和表 7可以看出,随着车桥质量比的增大,其识别出的频率值和识别频率误差基本无变化,仍能较好识别出桥梁的前3阶频率。通过图 6看出,随着车桥质量比的增大,识别出的桥梁频率峰值逐渐不明显,变得难以辨识。因此,为了获得较好的识别效果,建议选用小的车桥质量比。
改变桥梁的斜交角度,分析间接测量法提取桥梁频率的识别效果,选取8个不同的桥梁斜交角度,角度分别是0°、10°、20°、30°、40°、50°、60°、70°。在车重为4 t、车速为10 km/h的条件下进行车桥耦合振动计算得到车辆加速度时程数据,然后进行快速傅里叶变换,得到车辆加速度频谱图,识别结果如图 7所示。
从表 8中可以看出,斜交桥不同的斜交角度会对桥梁的刚度产生影响,随着斜交角度的增大,基频增高。
从图 7和表 9中可以看出,间接测量法对所有不同斜交角度桥梁的识别效果均较好,与理论频率的最大误差不超过4%,均满足工程精度要求。相对而言,斜交角度在70°以下时识别精度较好,对于斜交角度在50°时,2阶桥频和3阶桥频由于太过接近而混为一簇大峰,使得2阶和3阶桥频峰值识别困难,仔细辨别后仍能从频谱图中识别出桥梁频率相关峰值信息;斜交角度在70°以上时,由于斜边支撑线与桥梁轴线接近平行,识别效果变差。但斜交角度在70°以上的桥梁在实际工程中几乎不存在,因此,可以认为斜交角度对桥梁频率的间接识别精度几乎没有影响。
桥面不平整度是在实际应用中无法避免的一个重要影响因素。Chang等[19]提出了使用两个单自由度车模型匀速通过粗糙桥面,然后对得到的两个加速度时程曲线分别进行快速傅里叶变换,将得到的两个频谱图相减,以达到消除桥面不平整度对间接测量法的影响。Elhattab等[20]提出可以使用两个单自由度车模型匀速通过粗糙桥面,得到前后两车在过桥时段的加速度时程,然后选取它们在桥面同一位置时的加速度数值相减,使用加速度差值曲线进行桥梁模态参数识别,进而消除桥面不平整度对识别结果的影响,笔者采用后者方法进行频率识别。桥面不平整度参照《车辆振动输入与桥面平度表示方法》(GB/T 7031—1986)模拟,采用简谐波叠加法(三角级数法),生成101 m长桥面C级不平整度,如图 8所示。
选取桥面C级不平整度进行分析,用前后两辆相同1/4单轴车辆模型匀速以10 km/h驶过桥梁,两辆车间距1 m,车辆质量为4 t,刚度kv=170 kN/m,该单轴车辆的自振频率为1.038 Hz,在桥梁同一位置对前后车体加速度进行加速度差值,对差值数据进行处理,识别桥梁前3阶频率,与单车识别出的频谱图进行对比,结果如图 9所示。
从图 9中可以看出,单辆车无法识别粗糙桥梁的频率,两车差值方法可以准确识别出粗糙桥梁的前3阶频率,1阶频率识别误差为0.41%;2阶频率识别误差为1.30%;3阶频率识别误差为2.12%。可以看出C级桥面粗糙条件下差值法能够较好地解决粗糙度对桥梁频率识别的干扰。
车速低于20 km/h时能够较好地识别出斜交梁桥的前3阶频率,车速较高时无法识别桥梁的频率信息;相对较小的车桥质量比对桥梁频率识别有利;斜交梁桥不同的斜交角度基本不影响桥梁频率识别的精度,桥面粗糙条件下采用差值法仍能较好地识别桥频。数值模拟表明,间接测量法对于不规则斜交梁桥频率仍有较好的识别效果,可进一步推进移动传感间接测量方法在桥梁检测工程中的应用。
本方法的优点是省时省力、不中断交通;缺点是用弹簧来模拟非线性的车辆减振系统可能仍不够精确;现场试验时的有色噪声也是重要的干扰因素,因此,在工程中的应用还需要进一步的实桥试验研究。