基于线性模型的振动法拉索索力测量方法
【技术领域】
[0001] 本发明涉及一种振动法拉索索力测量方法，属于结构工程技术领域。
【背景技术】
[0002] 斜拉桥、悬索桥及吊杆拱桥均通过拉索将桥面系及车辆荷载传至主要受力构件， 统称索支承桥梁。拉索索力振动法测量是索支承桥梁检测及施工监控的重要技术手段。
[0003] 振动法的索力测试在运用时需要解决两个关键问题：1)拉索自振频率的识别，2) 索力-自振频率对应关系的确定。在现有的测量仪器及分析手段条件下，测定频率的精度 可达到0. 005Hz。因此，频率法索力测试的精确性主要取决于索力与自振频率关系的准确 性，即模型的精度。
[0004] 索力与自振频率的对应关系受诸多因素的影响，如拉索的长度、垂度、抗弯刚度以 及边界条件等。索力与自振频率的关系仅在两种简化模型下存在显式对应关系：
[0007] 上式中，T表示索力（N)，m表示拉索线密度（kg/m)，1表示拉索长度（m)，fk表示 k阶自振频率（Hz)，EI表示索截面抗弯刚度（Ν·πι2)，π表示圆周率常数。
[0008] 实际的工作中，索力与频率的关系难以满足以上两种模型假设。其中，长索受垂度 影响，产生非线性效应；短索受边界条件、截面抗弯刚度的影响，基于（1)、（2)式计算结果 比实际索力值偏差很大。因此，传统振动法索力测量的主要技术难点在于参数的识别，线密 度m，计算长度1。，索截面抗弯刚度ΕΙ，边界刚度等。其中计算长度1。表示拉索实际长度除 去边界约束及过渡段影响后的力学计算长度，一般与实际长度有一定差别。
[0009] 现有的测量方法对拉索的物理参数及边界条件等因素分情况进行考虑，得到多种 特定情况下的索力与自振频率的对应关系，导出的改进方法对（2)式的系数进行修正，总 体基本思路为：
[0010] a)、建立反映拉索相对长度的无量纲参数，划分长索，中索，短索；
[0011] b)、识别若干物理参数，如m，1。，EI，边界条件；
[0012] c)、基于特定低阶振动频率（1-3阶）推导计算公式。
[0013] 现有索力测量方法的主要问题有：
[0014] a)、必须对长索短索区别处理，根据无量纲参数分段给出计算式，分段计算；
[0015] b)、依赖于特定低阶频率，如前3阶频率或1阶基频。实际索力测试中，受场地条 件限制，对于较长拉索，加速度传感器难以布置到拉索中点，而是靠近桥面，这样测得的动 力响应以高阶频率为主，低阶频率精度不足将影响索力计算精度。同时没能有效利用高阶 频率包含的信息。
[0016] c)、实测信息只有频率一种物理量，要据此识别多个参数及边界刚度，会造成参数 多解或识别精度下降；
[0017] d)、每一种方法都有特定的物理模型假设，因而适用范围受限，同一种方法难以兼 顾不同型号拉索及不同的边界条件。
[0018] 在工程实践中，往往一座桥梁的拉索系统包含多种长度，直径规格，边界条件，空 间垂度，各不相同。这就导致了基于特定物理模型假设的索力测试方法，难以对多种拉索索 力统一识别。此外，不同的测量方法的参数定义不同，混合使用操作困难。
【发明内容】
[0019] 本发明的目的是提供一种基于线性模型的振动法拉索索力测量方法，以解决现有 技术中基于特定物理模型假设的索力测试方法难以对多种拉索索力统一识别，不同的测量 方法的参数定义不同，混合使用操作困难，影响测量精度的问题。
[0020] 本发明为解决上述技术问题采取的技术方案是：
[0021] -种基于线性模型的振动法拉索索力测量方法，在桥梁调索施工中记录待测拉索 的2级或以上的张拉力，并测试所述待测拉索的相应受力水平下振动频率，从而识别线性 模型的系数以确定线性模型；在运营阶段测试待测拉索的振动频率，利用所确定的线性模 型识别待测拉索的索力；
[0022] 所述方法的实现过程为：
[0023] 步骤一、根据设计资料，确定待测拉索的设计索力的范围为[T_，T_];
[0024] 步骤二、在调索施工过程中，在[Τ_，Τ_]范围内施加不同的张拉力下，分别标定 待测拉索的索力T1与及其对应的m阶振动频率f u…fini数据，
[0025] 标定次数为2次或以上；下角标i用于区分不同的张拉力；
[0026] 步骤三、根据待测拉索的索力T1与及其对应的m阶振动频率f u…fini数据拟合线 性模型的系数Λ-,或，k e [l，m]，线性模型的表达式为
[0028] 4表示通过4估计的索力值，
[0029] fk表示第k阶振动频率；带下角标i表示测量数据，没带下角标i表示通式；
[0030] 步骤四、计算索力1\与
^的线性相关系数Pk;
[0032] 式中，Cov(·)为统计学中的协方差函数，将向量1\与
做为随机变量，计算其协 方差；D( ·)为统计学中的方差函数，计算1\与
的方差；
[0033] 步骤五、拟合系数為，4 :
[0034] 若线性相关系数P，0. 95则说明1\与
存在显著线性关系，式（1)表达的线性 模型成立，通过最小二乘法乘拟合确定系数,
[0037] 式中，Ε( ·)为统计学中的期望函数；
[0038] 步骤六、在桥梁运营阶段测试待测拉索的振动频率匕…圪；
[0039] 步骤七、将步骤四获得的振动频率匕…圪、已经确的Ji,氧系数代入线性模型中，得 到m个索力估计值：
[0041 ] 步骤八、取均值可以得到的识别索力数值
[0042] 在步骤五中，
[0043] 当标定次数为2次时（即i = 1，2)，可以采用以下简化公式直接计算反，
[0046] 本发明的有益效果是：
[0047] 本发明方法在2级拉力水平下分别张拉拉索，在每级应力水平下使用加速度传感 器测量环境激励或人工激励作用下加速度信号，将其转换为自功率频谱图，在图中识别出1 至5阶频率。将2级拉力水平下各阶频率与拉力代入线性模型，识别每阶频率对应的线性 模型系数，实现模型参数识别。由识别后的线性模型和各阶固有振动频率即可实现拉索索 力识别。本发明基于线性模型方法可以得到精度较高的索力测量值，回避了传统振动法索 力测量物理参数的识别问题，降低了索力测试中的模型不确定性。
[0048] 本发明克服了现有技术的不足，提出了利用施工索力张拉的时机，基于2级索力 与频率实测数据，建立了一种基于线性模型的新的索力测量方法。
[0049] 实际索支承桥梁在施工过程中，往往需要对张拉索力进行调整，对拉索多级别分 批张拉。施工调索过程的千斤顶张拉力为真实索力，测试相应索力水平下的拉索振动频率， 这些数据即为未来的索力识别依据。
[0050] 业已证明，对于特定拉索，其第k阶频率的平方与拉力呈线性关系：
[0052] 式中4,4为线性回归系数。实测表明，当索力发生变化时4,為比较稳定，这就使 得基于线性模型测量索力成为可能。这一方法的最大优势在于回避了对具体物理参数的识 另|J。在新建桥梁或桥梁加固等具备张拉索力的条件下，利用2级张拉力及其对应频率数据 对回归系数為求解，从