基于残余力向量的结构损伤识别的两步法.docxVIP

基于残余力向量的结构损伤识别的两步法 模型修正法是结构损伤识别方法的一种非常重要的方法。它通过修正损伤后的结构模型以使其满足新的测量结果。这类方法又可细分为很多种方法, 如残余力向量法, 最小范数修正法, 最小秩修正法, 灵敏度法, 特征结构分配法等。Doebling 推广了最小秩修正法, 直接采用单元体刚度参数进行最小秩修正。Lim等人在最小范数修正法基础上, 提出了最佳逼近特征向量法。其基本思想是:先假设仅某一处单元体发生损伤, 对每个单元体建立一个最佳逼近特征向量, 考察该向量与测量振型之间的接近程度来确定可能损伤位置, 在此基础上再计算损伤程度。 本文在以上两种方法的基础上, 提出了一种新的损伤识别两步法。首先通过计算各单元体的损伤定位标准 (DLAC, damage location assurance criterion) 值来判定可能损伤单元体, 再利用最佳逼近的思想来精确定位并计算损伤程度。以一个8自由度的质量—弹簧系统为例验证了本文方法的可行性。算例表明该方法仅需一阶模态参数便可有效进行损伤识别。 1 损伤前后系统刚度矩阵的变化 假设损伤前后质量矩阵不变, 则结构系统损伤后的特征方程为 其中,M表示质量矩阵,Ku,Kd分别表示损伤前后系统的刚度矩阵, ΔK表示损伤前后刚度矩阵的变化, 它们均为n阶方阵,λdj, ?dj分别表示损伤后第j个特征值和振型。 方程 (2) 代入 (1) 中, 整理可得 记向量bj= (Ku-λdjM) ?dj, 方程 (3) 可改写为 方程 (4) 即为残余力方程, 其中,bj称为第j个残余力向量。 2 受伤识别的两个阶段法 2.1 刚度联系向量与残余力向量 损伤前后刚度变化矩阵ΔK可表示为 其中,αi为第i个单元体的刚度损伤参数,αi=0表示单元体未损伤,αi=1表示单元体完全损伤。 A称为刚度联系矩阵 (n×N维,n为自由度数目,N为单元体数目) , 且A=[a1a2…aN], 其中ai为第i个单元体的刚度联系向量。A在损伤前后不变。 方程 (5) 可进一步展开成 方程 (6) 代入 (4) 可得 方程 (7) 中aTi?dj为一个数值, 令 方程 (8) 代入 (7) , 有 由方程 (9) 可知, 残余力向量是各单元体刚度联系向量的线性组合。假设仅第i个单元体损伤, 即仅αi≠0, 其它N-1个刚度损伤参数均为零, 方程 (9) 变为 定义单元体的刚度联系向量ai与残余力向量bj间的夹角余弦值的平方为损伤定位标准DLAC, 即 若DLACij=1, 则损伤仅由第i个单元体引起, 为单个损伤情况。对于多个损伤, 不存在DLACij=1的情况, 此时, 可依次计算各单元体的DLAC值, 该值较大者即为可能发生损伤的单元体。 2.2 未损伤单元体的计算 不失一般性, 设可能损伤单元体为 (以各单元体的刚度联系向量表示)a1,a2, …,am, 此时方程 (9) 变为 记L1=[a1,a2, …,am], 方程 (12) 变为 由方程 (13) , 对bj建立其最佳逼近向量 其中, 上标“+”表示广义逆。定义向量bj,ba1j1ja的偏差为 其中, ‖·‖2表示2_范数。当d1j为0时, 表明可能损伤单元体没有漏选, 实际损伤的单元体必包含在a1,a2, \:.am中。 利用如下的筛选方法来去掉未损伤单元体, 保留实际损伤单元体: ① 假设a1未发生损伤, 先去掉a1, 可能损伤单元体组合形成矩阵L2=[a2,a3, …,am], 则方程 (13) 变为 相应最佳逼近向量为 计算偏差为 若d2j仍很小 (近似为0) , 表明a1未发生损伤, 下轮可能损伤单元体组合中不再保留它;若d2j较大 (相对于0而言) , 则表明a1为实际发生损伤单元, 下轮计算中仍应保留它。 这里不妨假设a1未发生损伤, 将其从组合中剔除。 ② 再假设a2未发生损伤, 去掉a2, 可能损伤单元体组合形成矩阵 同理可计算d3j, 根据它的值大小来决定a2是否损伤及下轮计算中是否保留。 依此类推, 最终可将未发生损伤单元体逐次剔除, 剩下的就是实际损伤单元体。 ③ 将最终所得损伤单元体组合成矩阵Lr, 各对应βi组合成向量{β}, 由方程 (12) , 有 于是可解得 方程 (20) 联立方程 (8) 即可求得各损伤程度值αi。 3 计算 以图1所示的8自由度质量—弹簧系统为例, 来验证上述方法的可行性。 系统的刚度联系矩阵为 其中 系统损伤前刚度矩阵和质量矩阵分别为 3.1 非邻单元体dlac1的值 设第3个弹簧刚度损失20%, 即损伤后k3=0.8, 其它弹簧不变。仅计算第一阶模态, 可得 各单元体的DLACi1值列于表1。 由表1可见, DLAC31=1, 而其相邻单元DLAC21=DLAC41=0.2