基于单位脉冲相应的框架结构损伤检测.pdfVIP

基于单位脉冲相应的框架结构损伤检测 宰吕中荣，刘济科 (中山大学工学院，广东广州510275) 摘要：建立了空间框架结构的有限元振动分析方程。利用Newmark直接积分法计算了结构在单位脉冲激励下的振动响应。 将结构的局部损伤模拟为单元杨氏模量的减少．推导了加速度响应对单元杨氏模量的灵敏度。并利用此响应灵敏度进行局部 损伤识别．数值算例表明，本文方法能快速准确地识别出结构的局部损伤，并且对模拟的人工噪声不敏感。说明本文方法具 有一定的工程实用前景。 关键词：框架，脉冲激励，损伤检测，响应灵敏度 1前言 近年来，工程结构的安全性问题日益引起人们的广泛关注。基于振动响应的结构损伤识别法具有非破 坏性、方便、快速的优点，已经广泛应用到工程结构的损伤检测中【1．5】。 本文建立了多层空间框架结构的有限元分析模型。在正问题分析中，利用直接积分法计算了结构在单 位脉冲激励下强迫振动响应。在反问题中，将结构的局部损伤模拟为单元杨氏模量的减少。进一步推导了 动态响应对单元杨氏模量的灵敏度。并利用此响应灵敏度进行结构的局部损伤识别。在算例中，对结构多 个局部损伤进行了识别，同时研究了人工噪声对识别精度的影响。算例表明，本文方法对模拟的人工测量 噪声不敏感，即使在较高噪声下，也能够准确地识别出结构损伤。 2正问题 2．1有限元运动方程 对于一般的工程结构通过有限元法离散，其强迫振动方程可以写成如下形式， 【gl{a)+[C】{d)+【K】{d)={，(f)) (1) 直接积分方法来计算。 3反问题 3．1动力晌应对损伤参数的灵敏度 在反问题中，利用基于时域响应灵敏度的有限元模型修正法来进行结构的局部损伤识别。在本文中， 目，广东省科技计划项目(20078030402002)．教育部新教师基金资助项目(200805581018)． 作者简介： ·吕中荣．男，湖南人．副教授．主要从事结构损伤检测方法研究(EmajI：lvzhrhai l．sysu．edu．cn)： 刘济科，男．湖北人．教授。#婴从事结构振动理论及应用研究． ·nI-430 我们用单元杨氏模量的减少来模拟某单元上的局部损伤。 结构的动态响应可以通过数值积分方法从式(1)获得，然后对结构的损伤参数，例如第i个单元杨 氏模量E4求偏导数有 ㈨{署)+【c】{熹}+旧{砉)一警渤一口I筹渤㈣，2，．．．，n) (2) 这里，”为结构的单元总数，{熹)，{要)，{熹)分别为位移、速度和加速度对第i个单元杨氏模量的灵敏 度。由于系统动态响应已经从式(1)获得，响应对单元杨氏模量的灵敏度可以类似地从式(2)由直接积 分法得到，进一步便可形成响应灵敏度矩阵。 3．2损伤参数的识别 识别问题可以表达为：寻找杨氏模量{E)使得计算出来的响应(如加速度响应)与测量的响应最好地 匹配，也就是 【Q】{‰)={R) (3) 其中，选孥矩阵(Q】是常数矩阵，其元素为0或1，与对应的测量响应分量的自由度相吻合；％为计算 的响应，月为测量的响应。 识别方程可以表示为 {8R}=IS]{aA) (4) 的偏导数，可以从(2)式获得。 杨氏模量增量向量{万E}可以从阻尼最小平方法获得【6】 {万E)=(SrS+A，)．1Sr{aR) (5) 其中，，为单位矩阵，旯是正则化参数，本文应用L曲线方法【7】获得最优的正则化参数旯。 修jE后的杨氏模量向量为 {E)={E)+{万E) 其中，{龟}没有损伤时结构的杨氏模量向量。 4数值算例 如图l所示的五层空间框架结构，经有限元离散为17个节点，37个空间梁单元，其左端三个节点 Nl，N7，N13为固定端。设初始单元杨