有限差分法在斜拉桥拉索内力测试中应用.docVIP

有限差分法在斜拉桥拉索内力测试中的应用 陈水生1 温世游1 谢彩文2 （1 华东交通大学土木建筑学院 南昌 330013） （2 江西省公路开发总公司 南昌 ） 摘 要：使用有限差分法将拉索自由振动偏微分方程转化为常微分方程，用状态空间法求解广 义特征值，从而计算系统的模态频率。将计算得到的模态频率与测试得到的模态频率比较，进 而修正拉索张拉力，直到计算频率与测试频率近似相等，最后修正的张拉力则为拉索内力。通 过对实际工程的测试结果分析表明，本文方法与传统计算方法的不同，说明了本文算法的通用、 准确、可靠性。 关键词：桥梁工程；有限差分法；斜拉索；测试技术；模态频率；中图分类号：U443.38;O329 0 前言 斜拉索是斜拉桥的主要承重构件之一，也是为斜拉桥主梁提供弹性支撑的构件，拉索内力的大小直接影响拉索自身及主梁的安全。所以，如何准确的测试拉索内力，是保障大桥安全运营极为重要的研究课题。对于正在使用的斜拉索，使用振动测试是首选方法，通过振动响应信号，由频谱分析得到拉索的振动基频，然后通过拉索振动频率与内力的关系求得内力。影响拉索内力的准确性有两个方面：其一是准确的测定拉索的振动频率，这涉及测试设备的质量和频谱分析的准确性；其二是振动频率与拉索内力关系方程的准确性。目前广泛使用的是由标准弦方程得到的内力与频率方程[1]，这对拉索垂度很小的斜拉索是适用的，而对拉索垂度较大的斜拉索将有较大的误差。大量研究表明，拉索振动频率与拉索垂度、抗弯刚度有关。M.Pacheco[2]及陈[3]的研究均表明，随拉索垂度的增大拉索振动频率增大。所以，为了准确计算拉索振动频率，考虑垂度及抗弯刚度的影响是必要的。 本文综合考虑索的垂度及抗弯刚度，使用有限差分方法将拉索自由振动方程离散，使用复特征值求解方法，计算拉索振动频率。由环境激励测试拉索振动响应，通过频谱分析得到索的振动基频。对钱塘江三桥进行了现场测试，使用本文方法计算了拉索内力，并与传统方法进行了比较。 1 斜拉索自由振动方程及特征值求解 不计索的轴向振动，由牛顿定律并考虑索的垂度及抗弯刚度，可以得出斜拉索面内自由振动方程[4]： （1） 其中，为拉索在Y方向的振动位移，为索的自重在Y方向产生的垂度曲线，为初始切向拉力，为索的振动拉伸而产生的附加切向动拉力， 为索的弧长坐标，为索的单位长度质量，为索的倾斜角度，为索的抗弯刚度。 认为拉索抗弯刚度沿索长为常数，由几何变换关系[5]，方程(1)可表示为： (2) 其中：分别为X方向静、动张拉力，在整个索长L范围内为常数。令的解表示为如下分离变量形式： (3) 其中可表示为： (4) 由线性阻尼系统理论可知，可表示为： (5) 为无阻尼索的面内自振频率，为对应的模态阻尼比。方程(2)中的线性表达式为[5]： (6) (7) 其中为跨中垂度在Y方向的值。 由(3)～(7)各式，方程(2)可记为： 图 4-1 斜拉索差分离散模型 用中心差分法求解，将索均分为段，每段长度。将方程(8)中的微分项化为各结点位移形式。对于非边界结点，有差分表达式： (9a) 式中 对边界结点（铰节点）： (10a) (10b) 将的积分转化为离散形式： (11) 由方程(8)可以得出如下特征方程： (12) (13) 为对角矩阵，对角元素为。，、、分别是抗弯刚度、静态张拉力、动态张拉力提供的等效刚度矩阵。 为求解(12)式的特征值问题，引进状态向量： (14) 则(12)式可以转化为以下广义特征值问题： (15) 其中： (16) 则由(15)式可以求得对共轭特征值，再由(5)式可以求出各阶模态无阻尼自振频率。 2 斜拉桥拉索自振频率测试及内力分析 由方程(8)和(13)式可知，特征分析得到的自振频率与拉索内力是对应的，即某一内力值，可以计算出拉索对应的各阶自振频率。使用频谱分析及参数识别技术，由环境振动测试得到的振动加速度信息，识别出拉索振动频率。 根据识别出的拉索振动频率，由标准弦振动频率方程得出一个力： 　　　　　 (17) 式中：为拉索的轴向内力；为拉索的阶振动频率；为拉索