基于静力位移数据桥梁结构损伤识别方法.docVIP

基于静力位移数据的桥梁结构损伤识别方法 摘 要：当前桥梁营运状态健康监测越来越受到重视，近十多年来国内外学者一直在寻找能适用于复杂结构的整体探伤评估方法，也使得结构的损伤识别在大跨度桥梁健康监测领域中成为一个研究热点[1]。本文在总结与研究当前国内外桥梁结构损伤识别理论的基础上，提出了一种基于静态位移测量数据结合有限元计算的结构损伤识别方法，与动力损伤指纹法相比，静态方法具有更高的识别精度与可操作性。数值模拟分析证明了该方法的有效性。 关键词：健康监测；损伤识别；影响矩阵；有限元； 一、前言 对于土木工程结构，尤其是对于以混凝土为主要材料的大跨度桥梁，其模态频率的变化对于结构局部损伤的敏感性极低，这一点已经在理论计算与工程实践中被证明[2]，而且基于模态振型的动力损伤指标[3]也往往因为模态振型测量的精度无法保证而失去作用，因此可以认为动力损伤指标(包括频率、振型等)并不适合于大跨混凝土桥梁结构的局部损伤探测。相反，静力位移测量的精度相对较高，为此，很多学者致力于通过静力测量数据进行结构探伤的研究[4、5]。本文在国内外现有的研究成果基础上，针对现有理论的不足，探索并提出了一种基于静力位移测量数据进行结构损伤识别的简单实用方法。 二、基于静力位移数据的损伤识别方法的研究 对基于静力位移测量数据的结构损伤识别方法的研究兴起于上个世纪90年代初，也是得益于现代控制理论中系统识别理论的快速发展。事实上，静力损伤识别方法也属于系统识别的范畴，它的目的就是根据系统的输入与响应来辨识系统内部的未知参数。 静力法损伤识别的基础就是有限元法的静力平衡方程： (1) 其中，[K]为结构的刚度矩阵，{u}为位移向量，{f}为外部荷载向量。当结构发生局部损伤时，在有限元模型中表现为刚度矩阵[K]发生了变化，假设外部荷载没有发生变化，从而引起位移{u}发生相应的变化，这个过程可以表示为一个静力摄动方程： (2) 其中，ΔX表示变量X对应的摄动量。 将(2)式展开后并结合(1)式便可以得到： (3) 对于梁单元结构，[ΔK]是由杨氏弹性模量E或者截面的惯性矩I发生损失后的单元刚度矩阵变化量的集成，因此[ΔK]可以看作是损伤单元刚度矩阵变化量的线性叠加，文献[4]将(3)式写成： (4) (5) 其中，[Kj]是由第j个单元的单元刚度矩阵扩展成的整体刚度矩阵；{ud}是发生损伤后的位移；αj表示第j个单元的损伤程度，并且有0≤αj1，j=1,2,…,p。并将(4)式表示成： (6) 其中， (7) (8) 至此，(6)式很容易求解，得到每个单元损伤程度的估计。 然而在实际测量中，我们很难在真实结构上测量与有限元模型的对应的全部自由度，而只能测量其中一部分，或者说我们在真实结构上测得的静力数据是不完备的，这对于复杂结构尤其明显。文献[5]中采用了自由度静态缩聚的办法使得不完备的实测数据能够应用于损伤识别。将基本平衡方程(1)式写成分块矩阵的形式： (9) 其中，下标a对应已测的自由度，b对应未测自由度。同时进一步得到： (10) 因此平衡方程可以写成： (11) 然而，我们可以发现由于[kbb]-1这一项的存在，导致[Ka]成为刚度系数的非线性函数。因此当采用自由度静态凝聚后，(6)式与(7)式中的线性叠加无法实现。文献[5]对缩聚后的刚度矩阵[Ka]变换形式后作了一阶Taylor级数展开： (12) 式中，p便是刚度参数，这一项就是前面(2)式中的[ΔK]，显然当采用静态缩聚后，[ΔK]也是关于刚度参数的非线性函数。 既然无法实现(6)式中的线性叠加，也就难以对(8)进行求解。而文献[5]中所提到的方法需要对整体刚度矩阵进行求导而得到关于刚度参数的敏感度矩阵，这种方法较为困难。本文结合结构影响矩阵概念[6]提出一种简单实用的求解方法。 对于上述问题可采用刚度参数－位移影响矩阵来求解。经过静态缩聚后，将基本平衡方程写成： (13) 显然，{ua}是关于{p}的非线性函数，于是原问题变成了一