[发明专利]基于截断奇异值分解法识别预应力桥梁现存预应力的方法

权利要求书

1.一种基于截断奇异值分解法识别预应力桥梁现存预应力的方法，其特征在于：包括以下步骤：

1)、在桥梁底面从左端依次向右布置有m个位移传感器；

2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移

3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}_N×1；

4)、建立识别桥梁现存预应力值的系统方程T＝({S}^T{S})^-1{S}^T{b}；

5)、利用基于截断奇异值分解法，先进行分解，然后引入滤波因子和修改截断参数得到改进的截断奇异值分解法，求得桥梁现存预应力值T。

2.如权利要求1所述的基于截断奇异值分解法识别预应力桥梁现存预应力的方法，其特征在于：所述步骤2)中具体包括以下步骤：

21)、建立桥梁的振动微分方程：

将桥梁考虑为一单跨简支梁，假设固定荷载P(t)位于单跨简支梁上，其距离桥梁左端距离为x_p，通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x，t)，则其振动微分方程：

其中单跨简支梁的跨长为L，密度为ρ，截面面积为A，粘性比例阻尼为c；桥梁抗弯刚度为EI，E是材料的弹性模量，I＝bh₀³/12是梁截面的惯性矩，其中b是桥梁横截面宽度，h₀是桥梁横截面高度；δ(x-x_p)是狄拉克函数；

22)、得到简支梁的模态振型函数的矩阵形式：

基于模态叠加原理，假设梁的第i阶模态振型函数为Y_i，梁的第i阶振幅为q_i(t)，则动态位移y(x，t)的模态形式表示为：

把方程(2)带入方程(1)并对方程两端的每一项均乘以Y_i(x)，考虑模态正交条件并对方程中x值从0到L进行积分，得到新方程：

其中，ξ_i和分别是第i阶模态的阻尼比和模态质量，f_i(t)＝P(t)Y_i(x_p)是位于单跨简支梁上的模态荷载；为q_i(t)的二阶导数，为q_i(t)的一阶导数；

简支预应力梁的模态振型函数标准化形式表示为：

在时域内通过Newmark积分将方程(3)用矩阵形式表示为：

其中，[I]是单位矩阵，{Q(t)}＝{q₁(t)，q₂(t)，...，q_n(t)}^T，[C]＝diag(2ξ_iω_i)，{F(t)}＝{f₁(t)，f₂(t)，...，f_n(t)}^T，和分别为Q(t)的一阶导数和二阶导数；

23)、得到桥梁竖向位移向量

假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器，则在模态坐标下距简支梁左端支座x_m处的第m个位移传感器实测竖向位移y(x_m，t)表示为：

其中，Y_i为梁的第i阶模态振型函数，q_i(t)为梁的第i阶振幅，方程(6)用矩阵形式表示为：

式中，是第N_m个位移传感器测得的位移向量，N是测量位移向量中包含的振型数量；

则广义坐标下基于最小二乘法的振幅伪逆解表示为：

式中，为矩阵的逆，[Y]^T为矩阵[Y]的转置；

布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(x_j，t)由广义正交多项式拟合为：

其中，N_f为正交多项式函数的阶数，G_i为第i阶正交多项式，a_i为第i阶正交多项式的系数；

将方程(9)带入方程(7)并写成矩阵形式得：

式中，[A]，[G]分别为正交多项式的系数矩阵和正交多项式矩阵；正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到：

将方程(11)带入方程(10)即可得到竖向位移向量