基于共轭残差平方算法识别预应力桥梁现存预应力的方法技术

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

【技术实现步骤摘要】
基于共轭残差平方算法识别预应力桥梁现存预应力的方法
本专利技术属于预应力桥梁现存预应力的识别
，尤其涉及一种基于共轭残差平方算法的识别预应力桥梁现存预应力的方法。
技术介绍
从受力性能的角度来看，预应力混凝土结构不仅具有自重轻、耐久性强、抗疲劳性强的特点，而且能够有效增强结构或构件的抗剪能力，近年来，预应力混凝土结构在水利和土木工程等领域得到了广泛应用。据不完全统计，在现代桥梁结构中，70％以上的桥梁结构都采用了预应力技术。但在超载型大宗货物车辆反复作用下，预应力桥梁混凝土结构不可避免的存在一定程度的损伤和破坏。预应力桥梁的现存预应力值可较直观的反应桥梁目前的受力状态和结构性能，但我国目前仍没有成熟的无损检测桥梁现存预应力方法。因此，提出可应用于工程实践的预应力桥梁现存预应力值无损检测方法是我国桥梁工程人员迫切需要解决的问题。
技术实现思路
本专利技术的目的是提供一种基于共轭残差平方算法由桥梁位移实时识别预应力桥梁现存预应力的方法，识别速度快且精度较高。为达到上述目的，本专利技术采用的技术方案是：一种基于共轭残差平方算法的识别预应力桥梁现存预应力的方法，包括以下步骤：1)、在桥梁底面从左端依次向右布置有m个位移传感器；2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1；4)、建立识别桥梁现存预应力值的系统方程T＝({S}T{S})-1{S}T{b}；5)、利用基于共轭残差平方算法，由系统方程初始解得到残差，可得到共轭残差平方算法的迭代解，进而求得桥梁现存预应力值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，[C]＝{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)通过定义共轭残差平方算法的两个常量参数α,β和四个辅助向量d,e,h,l以及一个搜索方向向量H可得到共轭残差平方算法的迭代解为：Tk+1＝Tk+αk(ek+hk)(18)其中，αk为常量参数α的第k步迭代，可表示为：式中，系统方程残差r的第k步迭代解为：rk＝rk-1-αk-1(dk-1+lk-1)(20)式中，搜索方向向量H的第k步迭代解为：Hk＝dk+βk-1(lk-1+βk-1Hk-1)(21)式中，四个辅助向量d,e,h,l的第k步迭代解为：dk＝{S}T{S}rk+βk-1lk-1ek＝rk+βk-1hk-1hk＝ek-αkHklk＝dk-αk{S}T{S}Hk(22)其中，βk为常量参数β的第k步迭代，可表示为：将公式(19)-(23)代入公式(18)，当Tk+1-Tk≤ε满足识别精度误差要求ε时，即可由共轭残差平方算法求得待识别的桥梁现存预应力值即为Tk+1，最终实现由桥梁动态位移响应实时识别现存预应力值。本专利技术具有的优点是：本专利技术提出一种基于共轭残差平方算法由桥梁位移实时识别预应力桥梁现存预应力的方法，通过建立预应力桥梁现存预应力与桥梁竖向动态位移响应之间的力学映射关系，实现由桥梁位移响应识别桥梁现存预应力值。为提高预应力值的识别效率，对预应力识别系统矩阵采用共轭残差平方算法以节约识别时间，最终实现实时识别桥梁现存预应力值，具有识别速度快且精度较高等优点，具有良好的工程应用价值。附图说明图1是本专利技术的方法流程图。具体实施方式如图1所示，本专利技术公开了一种基于共轭残差平方算法的识别预应力桥梁现存预应力的方法，包括以下步骤：(1)、在桥梁底面从左端依次向右布置有m个位移传感器；(2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移21)、建立桥梁本文档来自技高网...

【技术保护点】
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)、利用基于共轭残差平方算法，由系统方程初始解得到残差，可得到共轭残差平方算法的迭代解，进而求得桥梁现存预应力值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)由广义正交多项式拟合为：其中，Nf为正交多项式函数的阶数，Gi为第i阶正交多项式，ai为第i阶正交多项式的系数；将方...

【专利技术属性】
技术研发人员：陈震，陈贡联，陈爱玖，王慧，潘丽云，何伟，刘世明，裴松伟，吕亚军，姚贤华，
申请(专利权)人：华北水利水电大学，
类型：发明
国别省市：河南,41