框架结构损伤定位比例柔度矩阵分解法试验研究.docVIP

框架结构损伤定位的比例柔度矩阵分解法试验研究 安永辉，欧进萍 (大连理工大学土木工程学院，辽宁 大连116023) 摘要：本文从结构，。在实验室一个6层集中质量框架模型，基于振动台试验结果均表明：能准确地对结构损伤进行定位工程应用奠定, TU317.1 引　　言 结构损伤识别的核心问题之一是寻找与结构动力特性密切相关且对损伤敏感的特征参数[1]。结构损伤部位的刚度、强度等力学性能降低，从而导致动力特性发生改变。不同的损伤引起的结构模态参数的改变也不相同，因此可以基于结构模态参数的改变来识别结构损伤[2-3]。柔度矩阵由模态频率和振型构建而成，也集成了这两个模态参数的优点，在基于振动的损伤识别方法中有着较多应用[4-6]。 损伤识别关键要把握两点：第一，第二，数学工具来处理选定的，建立与损伤的联系。这两点直接关系到损伤的精度。①从结构，基于直接识别的模态参数构建结构模态柔度矩阵；③然后对损伤前后的柔度矩阵差进行LU分解；④最后引入以往在模态振型上使用的曲率法。结构损伤前后的加速度不确定因素模拟研究了在基底白噪声激励下的单多损伤工况在University of Illinois at Urbana-Champaign设计建造了一个6层剪切型集中质量框架模型，基于。LU分解的损伤定位方法 提出的方法基于比例柔度矩阵的LU分解和曲率法，下面给予详细介绍。 1.1 方法的理论基础 1.1.1 比例柔度矩阵的构建 如方程1[14]所示，柔度矩阵F可用动力测试模态参数（模态频率和振型）来表示， (1) 其中是第i阶质量归一化振型，是第i阶模态频率，n 是自由度数。方程1中可以看出，柔度矩阵可以仅仅用前几阶低阶振型近似得到较好的估计[14]。记和为从测试数据中基于模态参数识别技术ERA[15]（针对自由振动响应）或NExT[16]-ERA技术（针对环境激励响应）直接识别得到的模态频率和振型。要想得到质量归一化振型，需要至少一个测点上同时测试输入和输出；但是对于环境激励等随机激励来说，测试输入是很难甚至不可能的。本文通过构建比例柔度矩阵来代替柔度矩阵。段忠东等[11-12]提出了构建结构比例柔度矩阵的方法，指出比例柔度矩阵与真正的柔度矩阵仅相差一个比例因子。任意情况下构建比例柔度矩阵的方法详见文献[11-12]。当质量矩阵中各测点位置对应元素相等时（如文中结构或常见的质量均匀分布的结构），本文提出了一种简便的构建比例柔度矩阵的方法。 结构质量矩阵为对角矩阵，当质量矩阵中各测点位置对应元素相等时，假设其质量矩阵为 (2) 其中为常数，I为维数与M相同的单位阵。假设直接识别得到的振型与质量归一化振型之间的关系为 (3) 其中是常数。结构的质量归一化振型可通过下式求得 (i = 1, 2, …, n) (4) 将式(2)、(3)带入式(4)可得 (i = 1, 2, …, n) (5) 所以可求得 (6) 上式中只有常数未知，可假设=1，记为假设=1情况下的虚拟的，则二者之间的关系为 (7) 将式(3)带入式(1)可得结构的真实的柔度矩阵F为 (8) 假设质量矩阵M = I，即=1时结构的比例柔度矩阵为 (9) 比较式(8)和(9)，可得两个柔度矩阵之间的关系为 (10) 可知两者之间只相差一个比例因子，并且与文献[11-12]指出的该比例常数为结构的第一阶模态质量这一结论保持一致。 本文基于比例柔度矩阵Fp来进行损伤定位。损伤前后柔度矩阵的改变记作, (11) 其中和分别是无损和损伤结构的比例柔度矩阵。结构中有损伤时，其柔度矩阵改变。下一步是寻找一个指标可以将柔度矩阵差中对应损伤的那部分信息凸显出来。 1.1.2 比例柔度矩阵差的LU 分解 基于LU分解，一个方阵A总是可以分解为一个下三角矩阵L和一个上三角矩阵U的乘积。以一个维数为m×m的矩阵A（A=）为例， (12) 得到U 矩阵后，对于m×m维的U矩阵，定义向量UU为矩阵U的每列元素的平方和，即 (13) 1.1.3 基于向量UU曲率的损伤指标 曲率模态法由Pandey等[1]于1991年提出，之后研究者基于该法做了研究。曲率为: (14) 其中v为模态位移，h为数值模型中的单元长度。 目前，曲率方法已经成为损伤识别领域的一个基本数学工具。类似地，本文定义向量UU的“曲率”如下： (15) 其中C是节点处的曲率，是DI如下，