基于完全正交化算法识别预应力桥梁现存预应力的方法技术

本发明专利技术公开了一种基于完全正交化算法的识别预应力桥梁现存预应力的方法，包括以下步骤：1)、在桥梁底面从左端依次向右布置有m个位移传感器；2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移

【技术实现步骤摘要】
基于完全正交化算法识别预应力桥梁现存预应力的方法
本专利技术属于预应力桥梁现存预应力的识别
，尤其涉及一种基于完全正交化算法的识别预应力桥梁现存预应力的方法。
技术介绍
随着我国经济的迅猛发展，公路铁路交通量均日益增加，且在今后较长的一段时间内我国交通运输需求将继续保持快速增长趋势。货物运输中大宗货物车辆荷载重型化、车速高速化的发展虽然创造了良好的经济效益，但也对我国桥梁等交通基础设施造成了潜在的危害，因此对桥梁结构的安全性、适用性与耐久性提出更高要求。如何精确评估现役桥梁的结构安全状况，避免由于桥梁结构破坏导致的人员伤害和财产损失是桥梁工程人员迫切需要解决的工程难点。预应力桥梁现存预应力的大小是评估桥梁抗裂性能和承载能力的重要指标，采用无损检测技术快速识别预应力桥梁现存预应力具有良好的工程可行性。
技术实现思路
本专利技术的目的是提供一种基于完全正交化算法由桥梁位移实时识别预应力桥梁现存预应力的方法，识别速度快且精度较高。为达到上述目的，本专利技术采用的技术方案是：一种基于完全正交化算法的识别预应力桥梁现存预应力的方法，包括以下步骤：1)、在桥梁底面从左端依次向右布置有m个位移传感器；2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1；4)、建立识别桥梁现存预应力值的系统方程T＝({S}T{S})-1{S}T{b}；5)、利用基于完全正交化算法，由系统方程初始解得到残差，由残差构造得到w维Krylov子空间，求得桥梁现存预应力值T。所述步骤2)中具体包括以下步骤：21)、建立桥梁的振动微分方程：将桥梁考虑为一单跨简支梁，假设固定荷载P(t)位于单跨简支梁上，其距离桥梁左端距离为xp，通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x,t)，则其振动微分方程：其中单跨简支梁的跨长为L，密度为ρ，截面面积为A，粘性比例阻尼为c；桥梁抗弯刚度为EI，E是材料的弹性模量，I＝bh03/12是梁截面的惯性矩，其中b是桥梁横截面宽度，h0是桥梁横截面高度；δ(x-xp)是狄拉克函数；22)、得到简支梁的模态振型函数的矩阵形式：基于模态叠加原理，假设梁的第i阶模态振型函数为Yi，梁的第i阶振幅为qi(t)，则动态位移y(x,t)的模态形式表示为：把方程(2)带入方程(1)并对方程两端的每一项均乘以Yi(x)，考虑模态正交条件并对方程中x值从0到L进行积分，得到新方程：其中，ξi和分别是第i阶模态的阻尼比和模态质量，fi(t)＝P(t)Yi(xp)是位于单跨简支梁上的模态荷载；为qi(t)的二阶导数，为qi(t)的一阶导数；简支预应力梁的模态振型函数标准化形式表示为：在时域内通过Newmark积分将方程(3)用矩阵形式表示为：其中，[I]是单位矩阵，{Q(t)}＝{q1(t),q2(t),…,qn(t)}T，{F(t)}＝{f1(t),f2(t),…,fn(t)}T，和分别为Q(t)的一阶导数和二阶导数；23)、得到桥梁竖向位移向量假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器，则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为：其中，Yi为梁的第i阶模态振型函数，qi(t)为梁的第i阶振幅，方程(6)用矩阵形式表示为：式中，是第Nm个位移传感器测得的位移向量，N是测量位移向量中包含的振型数量；则广义坐标下基于最小二乘法的振幅伪逆解表示为：式中，为矩阵的逆，[Y]T为矩阵[Y]的转置；布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为：其中，Nf为正交多项式函数的阶数，Gi为第i阶正交多项式，ai为第i阶正交多项式的系数；将方程(9)带入方程(7)并写成矩阵形式得：式中，[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵；正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到：将方程(11)带入方程(10)即可得到竖向位移向量所述的步骤3)具体为：将竖向位移向量代入方程(8)即可求得预应力桥梁的振幅{q}N×1。所述的步骤4)具体为：所述的方程(5)重写为：式中，矩阵[K′]包含待识别的桥梁预应力值T，矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi，矩阵[K]包含预应力桥梁的系统参数；在时域内求解方程(12)即可得到桥梁现存预应力值T，方程(12)在时域内可表示为：其中则待识别预应力T由最小二乘法直接解得：式中，矩阵{B}T为矩阵{B}的转置，即为包含了方程(12)所有右端项信息的向量；由方程(14)即可求得在t时刻预应力桥梁现存预应力值T(t)，当取固定荷载P(t)在t1,t2,…,tw时刻的荷载值分别为P(t1),P(t2),…,P(tw)，则方程(12)在时域内可重写为：{S}N×wTw×1＝{b}N×1(15)其中则待识别预应力T由最小二乘法直接解得：T＝({S}T{S})-1{S}T{b}(16)式中，矩阵{S}T为矩阵{S}的转置，向量{b}N×1与类似，即为包含了方程(12)所有右端项信息的向量，可由桥梁动态参数和测量位移响应推求。所述的步骤5)具体包括以下步骤：取基于桥梁动力响应识别现存预应力值系统矩阵为{S}T{S}，由系统方程初始解T0得到的残差为：r0＝{S}T{b}-{S}T{S}T0(17)取由r0构造得到的w维Krylov子空间为：记β＝‖r0‖2，则经过w步Arnoldi算法可构造一组标准正交基Vw＝[v1,v2,…vw]和一个(w+1)×w的Hessenberg矩阵Hw+1，将Hessenberg矩阵Hw的最后一行元素删除即可得到大小为w×w的Hessenberg矩阵Hw，识别现存预应力值系统矩阵和标准正交基Vw以及Hessenberg矩阵Hw有下面的关系：完全正交化算法采用正交投影理论，由标准正交基Vw和Hessenberg矩阵Hw即可在Krylov子空间解得预应力桥梁现存预应力值Tw，其表达式为：其中，β＝‖r0‖2，e1为第一个元素为1的w维单位向量，由公式(20)即可求得待识别的桥梁现存预应力值即为Tw。本专利技术具有的优点是：本专利技术提出一种基于完全正交化算法由桥梁位移实时识别预应力桥梁现存预应力的方法，通过建立预应力桥梁现存预应力与桥梁竖向动态位移响应之间的力学映射关系，实现由桥梁位移响应识别桥梁现存预应力值。为提高识别方法的迭代效率，采用完全正交化算法以节约识别时间，最终实现实时识别桥梁现存预应力值，具有迭代收敛速度快、识别时间短等优点，具有良好的工程应用价值。附图说明图1是本专利技术的方法流程图。具体实施方式如图1所示，本专利技术公开了一种基于完全正交化算法的识别预应力桥梁现存预应力的方法，包括以下步骤：(1)、在桥梁底面从左端依次向右布置有m个位移传感器；(2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移21)、建立桥梁的振动微分方程：将桥梁考虑为一单跨简支梁，假设固定荷载P(t)位于单跨简支梁上，其距离桥梁左端距离为xp，通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x,t)，则其振动微分方程：其中单跨简支梁的跨长为L，密度为ρ，截面面本文档来自技高网...

【技术保护点】
1.一种基于完全正交化算法识别预应力桥梁现存预应力的方法，其特征在于：包括以下步骤：1)、在桥梁底面从左端依次向右布置有m个位移传感器；2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移

【技术特征摘要】
1.一种基于完全正交化算法识别预应力桥梁现存预应力的方法，其特征在于：包括以下步骤：1)、在桥梁底面从左端依次向右布置有m个位移传感器；2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1；4)、建立识别桥梁现存预应力值的系统方程T＝({S}T{S})-1{S}T{b}；5)、利用基于完全正交化算法，由系统方程初始解得到残差，由残差构造得到w维Krylov子空间，求得桥梁现存预应力值T。2.如权利要求1所述的基于完全正交化算法识别预应力桥梁现存预应力的方法，其特征在于：所述步骤2)中具体包括以下步骤：21)、建立桥梁的振动微分方程：将桥梁考虑为一单跨简支梁，假设固定荷载P(t)位于单跨简支梁上，其距离桥梁左端距离为xp，通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x，t)，则其振动微分方程：其中单跨简支梁的跨长为L，密度为ρ，截面面积为A，粘性比例阻尼为c；桥梁抗弯刚度为EI，E是材料的弹性模量，I＝bh03/12是梁截面的惯性矩，其中b是桥梁横截面宽度，h0是桥梁横截面高度；δ(x-xp)是狄拉克函数；22)、得到简支梁的模态振型函数的矩阵形式：基于模态叠加原理，假设梁的第i阶模态振型函数为Yi，梁的第i阶振幅为qi(t)，则动态位移y(x，t)的模态形式表示为：把方程(2)带入方程(1)并对方程两端的每一项均乘以Yi(x)，考虑模态正交条件并对方程中x值从0到L进行积分，得到新方程：其中，ξi和分别是第i阶模态的阻尼比和模态质量，fi(t)＝P(t)Yi(xp)是位于单跨简支梁上的模态荷载；为qi(t)的二阶导数，为qi(t)的一阶导数；简支预应力梁的模态振型函数标准化形式表示为：在时域内通过Newmark积分将方程(3)用矩阵形式表示为：其中，[I]是单位矩阵，{Q(t)}＝{q1(t)，q2(t)，...，qn(t)}T，[C]＝diag(2ξiωi)，{F(t)}＝{f1(t)，f2(t)，...，fn(t)}T，和分别为Q(t)的一阶导数和二阶导数；23)、得到桥梁竖向位移向量假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器，则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm，t)表示为：其中，Yi为梁的第i阶模态振型函数，qi(t)为梁的第i阶振幅，方程(6)用矩阵形式表示为：式中，是第Nm个位移传感器测得的位移向量，N是测量位移向量中包含的振型数量；则广义坐标下基于最小二乘法的振幅伪逆解表示为：式中，为矩阵的逆，[Y]T为矩阵[Y]的转置；布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj，t)由广义正交多项式拟合为...

【专利技术属性】
技术研发人员：姚建斌，姚文志，杜培荣，刘明堂，陈俊峰，郑辉，
申请(专利权)人：华北水利水电大学，
类型：发明
国别省市：河南,41