基于频比偏移系数的改进线性模型振动法索力测量方法
【技术领域】
[0001] 本发明涉及一种振动法索力测量方法，属于结构工程技术领域。
【背景技术】
[0002] 振动法索力测量是索支承桥梁、体外预应力桥梁检测拉索受力状态的主要手段。 其机理即拉索索力与其固有振动频率之间存在稳定的函数关系。振动法索力测量普遍被认 为精度不足的方法。因为，在拉索长度相对较短的情况下，即l〇m以内时，边界条件及自身 刚度的影响较为显著。其索力-频率关系较弦模型变得复杂。
[0003] 为了解决这一问题，目前已有许多学者提出了考虑边界条件、有效长度、截面抗弯 刚度的方法。这些方法的核心都在各阶频率与索力之间建立考虑上述因素的最小二乘方 法。从本质上讲，索力识别同其他因素的共同识别是一个基于样本观测的逆问题。当同时 识别的参数较多和样本数量有限时，逆问题面临多解和解不稳定的问题，且效率较低。
[0004] 索力识别逆问题的实质是由样本值找到其他频率值处对应的索力。现有的基于线 性模型的振动法拉索索力测量方法较目前其他方法更为简洁且精度较高。但是将多阶频率 分别考虑，对于每一阶频率分别计算线性回归系数。即基于线性模型的振动法拉索索力测 量方法将各阶频率的样本值视为相互独立，没有充分考虑各阶频率之间的联系。
【发明内容】
[0005] 本发明的目的是提供一种基于频比偏移系数的改进线性模型振动法索力测量方 法，以充分考虑各阶频率之间的联系构建线性回归模型，实现精度较高的索力估计值测量。
[0006] 本发明为解决上述技术问题采取的技术方案是：
[0007] -种基于频比偏移系数的改进线性模型振动法索力测量方法，所述方法的实现过 程为：
[0008] 在桥梁调索施工中记录待测拉索的2级或以上的张拉力，并测试所述待测拉索的 相应受力水平下振动频率；通过频比偏移系数确定的频率分布特征，导出名义一阶频率，识 别索力关于名义一阶频率的改进线性模型；在运营阶段测试待测拉索的振动频率，利用所 述改进线性模型识别待测拉索的索力；
[0009] 具体实现过程为：
[0010] 步骤一、根据设计资料，确定待测拉索的设计索力的范围为[Tmin，T_];
[0011] 步骤二、在调索施工过程中，在范围内施加不同的张拉力下，分别标定 待测拉索的索力Ti，通过拉索振动自功率谱频谱图识别出拉索m阶自振频率数据[f n-fim]，标定次数为2次或以上；下角标i用于区分不同的张拉力；
[0012] 步骤三、计算标定索力水平下的各高阶频率对应的频比偏移系数
[0014] fk表示第k阶振动频率；带下角标i表示测量数据，没带下角标i表示通式；
[0015] 步骤四、标定数据的频比偏移系数均值式关于k的线性回归
[0017] 得到频比偏移系数回归值a、b均为线性回归系数；
[0018] 步骤五、计算名义一阶频率J；
[0020] 步骤六、建立索力T关于名义一阶频率平方的线性模型
[0022] A、B为线性模型的系数；
[0023] 计算索力凡与/；；的线性相关系数P ;
[0025] 式中，Cov( ?)为统计学中的协方差函数，将向量1\与做为随机变量，计算其协 方差；D(_)为统计学中的方差函数，计算1\与的方差；
[0026] 步骤七、计算系数A，B ;
[0027] 若线性相关系数P >0. 95则说明凡与力存在显著线性关系，线性模型 7'= B成立，通过最小二乘法乘拟合确定系数A, B，
[0030] 式中，E(〇为统计学中的期望函数；
[0031] 步骤八、在桥梁运营阶段测试拉索振动频率计算测定频率的名义一阶频率 j及其#，由标定索力确定的频比偏移系数回归值I,计算名义一阶频率j ,将振动频率 代入线性回归模型r=j/f+5,得到索力估计值：
[0032] 当标定次数为2次时（即i = 1，2)，可以采用以下简化公式直接计算A，B，
[0035] 本发明的有益效果是：
[0036] 本发明方法在2级拉力水平下张拉拉索，在每级索力水平下识别出1至5阶频率。 本发明引入频比偏移系数，根据频比偏移系数对索力水平不敏感的性质，建立频比偏移系 数关于频率阶数的线性回归关系.用回归后的频比偏移系数，将高阶频率转化为名义一阶 频率.建立索力与名义一阶频率的线性回归模型，实现参数识别。由识别系数后的线性模 型和各阶固有振动频率即可实现拉索索力识别。本发明提出了频比偏移系数，反映了拉索 物理模型的固有性质，综合各阶频率信息可以得到精度较高的索力测量值。基于本方法的 索力识别，在不同的索力水平下预测精度不变，降低了索力测试中的不确定性。
[0037] 本发明所述的基于频比偏移系数的改进线性模型振动法索力测量方法，引入了 "频比偏移系数"以描述拉索振动模型的频率分布特征，引入了"名义一阶频率"以集成多阶 频率信息。
[0038] 理想的弦模型认为高阶频率与一阶基频存在倍频关系：
[0040] fk表示k阶频率，k = 1，2…，f i表示1阶频率。而实测数据表明，由于拉索抗弯 刚度以及边界条件的影响，振动频率不满足倍频关系：
[0042] 实际拉索与理想的弦模型的差异可以理解为，真实高阶频率fk是在基频倍数kf : 基础上的偏移。定义频比偏移系数为
[0044] 通过实验数据统计发现，Xk在不同的索力水平下基本稳定，并且关于阶数k呈线 性关系。附图4可以看出，A k在不同的索力水平下数值基本一致，且关于阶数k线性变化。 即对索力水平不敏感，相比自振频率fk随着索力增大而增大，所以，可以认为X k直接 反映了拉索物理模型的固有特征。
[0045] 据此，考虑利用将多次测量的频率数据统一到一个参数一一名义一阶频率：
[0047] 建立索力T与名义一阶频率平方的线性关系。
[0049] 式中A，B为线性回归系数。根据实验数据统计，以上线性关系仅与拉索物理模型 相关，对索力水平不敏感。本发明提出了结合施工过程的调索数据，利用2级张拉力及其对 应频率数据对确定改进模型系数A，B，在桥梁服役阶段将自振频率对应的名义一阶频率代 入改进模型即可得到精度较高的索力估计值。
【附图说明】
[0050] 图1为基于频比偏移系数的改进线性模型振动法索力测量方法操作流程图；
[0051] 图2为实施例中拉索振动自功率谱频谱图（张拉索力水平T = 180. 7kN);
[0052] 图3改进线性模型与索力识别图；
[0053] 图4频比偏移系数Xk在不同索力水平下线性统计特征图。
【具体实施方式】
[0054] 下面结合附图，以一根型号S4的拉索的张拉测试试验为例，对本方法进行进一步 详细说明。
[0055] 本例拉索为S4型，截面包含4束4X7根钢丝，间隙及外部由聚乙烯（PE)包裹。索 长 1 = 11. 57m，线密度 m = 5. 59kg/m，截面积 A = 556mm2,极限索力 Tlim= 1041. 60kN，弹性 模量E = 2. 0X 105MPa。截面抗弯惯性矩I是梁的力学参数，索的截面不是连续的整体，不 存在理论意义上抗弯刚度的性质。
[0056] 根据以上试验资料，对本方法的实施进行印证。
[0057] 1.根据设计资料，确定拉索工作状态的索力范围[Tmin，Tmax];