基于单元模态应变能的梁结构分级损伤识别方法研究.pdf

欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟收稿日期：基金项目：河北建筑工程学院校级研究生创新基金（）作者简介：吴文斌（），男，硕士研究生，从事防灾减灾领域的研究。通讯作者：杨海军（），男，教授，博士，主要从事结构优化设计、防灾减灾领域的研究。犱 狅 犻：基于单元模态应变能的梁结构分级损伤识别方法研究吴文斌杨海军，马磊蒋亚贤伍莎（河北建筑工程学院，河北 张家口 ；张家口市工程力学分析重点实验室，河北 张家口 ）摘要：针对目前利用单元模态应变能对结构进行损伤识别时精度过低、计算量大的问题，提出一种基于单元模态应变能变化率分级损伤识别方法。首先通过计算单元的各阶平均模态应变能变化率，初步判断损伤位置或损伤疑似区间；然后在损伤疑似区域内施加质量块，通过计算单元的各阶平均模态应变能变化率，进行损伤二次定位；最后利用遗传算法对损伤位置的损伤程度进行量化分析，确定损伤位置的损伤程度。关键词：分级损伤识别；模态应变能；附加质量；遗传算法；有限元分析中图分类号：犜 犝文献标识码：犃 引言结构的损伤识别主要包括两个方面，一是结构损伤位置的判断，二是结构损伤程度的确定对于损伤位置的判断，刘文光等推导了弹性薄板的单元模态应变能，并构造损伤识别指标对弹性薄板进行了损伤定位分析，得出单元模态应变能可以定位结构的损伤位置，但这种损伤识别方法确定的损伤是一个区域，真实的损伤位置在这一区域内，还需要通过优化算法对这一区域内的所有单元进行分析，判断出精确损伤位置与对应的损伤程度对于损伤程度的确定，通常是通过建立损伤结构的优化数学模型，结合损伤检测技术来进行的，遗传算法因其良好的寻优能力被广泛应用这一领域缪炳荣提出一种利用单元模态应变能与优化算法结合的损伤识别方法，该方法通过单元模态应变能确定结构的损伤区域，再结合遗传算法组合优化的思路对损伤区域的精确损伤位置与损伤程度进行判断虽然能够最终确定损伤位置与损伤程度，但是初步定位的损伤区域过大，导致遗传算法的初始变量过多，计算量较大因此，基于目前损伤识别中存在的定位不精准、计算量大的问题，本文基于单元模态应变能和遗传算法，通过附加质量块的方法进行精确损伤定位以较小的计算量确定损伤程度，实现结构的分级损伤识别 理论公式及原理分析由文献 知，结构损伤后，由于损伤位置与损伤程度未知，结构自身力学属性的改变量无法确定，单元刚度矩阵难以确定，因此可用损伤前的单元刚度矩阵近似代替损伤后的单元刚度矩阵，研究表明，这样处理后的单元模态应变能对损伤更为敏感因此，损伤前后的单元模态应变能及单元模态应变能变化率可表示为：犕 犛 犈犻 犼犻 犼犜犓犼 犻 犼（）犕 犛 犈犻 犼犻 犼犱犜犓犼 犻 犼犱（）第 卷 第期 年月河 北 建 筑 工 程 学 院 学 报犑 犗 犝 犚 犖 犃 犔犗 犉犎 犈 犅 犈 犐犐 犖 犛 犜 犐 犜 犝 犜 犈犗 犉犃 犚 犆 犎 犐 犜 犈 犆 犜 犝 犚 犈犃 犖 犇犆 犐 犞 犐 犔犈 犖 犌 犐 犖 犈 犈 犚 犐 犖 犌犞 狅 犾 犖 狅 犑 狌 狀 犕 犛 犈 犆 犚犻 犼狘犕 犛 犈犻 犼犱犕 犛 犈犻 犼狘犕 犛 犈犻 犼（）式中：犕犛 犈犻 犼、犕犛 犈犻 犼犱分别表示损伤前后的第犼个单元的犻阶单元模态应变能，犕犛 犈 犆 犚犻 犼表示第犼个单元的犻阶单元模态应变能变化率，犻 犼 表示第犼个单元的犻阶位移列向量，犓犼 表示第犼个单元的单元刚度矩阵考虑到高阶模态时，由于噪声的影响，结构的振型频率将与实际情况偏差较大，因此，用单元的平均单元模态应变变化率作为损伤指标，这样将降低噪声的影响，更趋于合理化，公式如下：犕 犛 犈 犆 犚犼犿犿犻狘犕 犛 犈犻 犼犱犕 犛 犈犻 犼狘犕 犛 犈犻 犼（）结合式（）、（）、（）可知，结构的单元模态应变能变化率受单元的节点位移列向量犻 犼 的影响，孙虎等的研究表明，在结构的损伤区域施加质量块，结构的模态振型将发生改变，从而引起相应区域单元模态应变能变化率的改变，因此可以通过施加质量块的方式来达到精确定位损伤的目的在精确定位损伤位置的基础上，以损伤单元的单元模态应变能变化率为目标函数，利用遗传算法对损伤程度进行量化分析，判断损伤位置的损伤程度 模型算例 建立有限元模型及损伤工况以两端均为固定端的梁结构为研究对象，进行结构的损伤定位与损伤程度的判断在 中用 梁单元进行有限元建模，弹性模量为 ，泊松比为 ，密度为 建立长为 ，截面为 的矩形截面梁，划分 个单元，有限元模型如图，梁结构单元及节点编号如图图梁结构有限元模型图梁单元及节点编号分布相比于无损伤工况，改变预设损伤位置的弹性模量，达到模拟损伤的效果设定四个损伤工况进行分析研究，见表：表预设损伤工况表设定损伤单元预设损伤程度工况，工况，工况，工况，确定损伤位置 初步损伤定位分别对损伤工况与未损伤工况进行五阶模态分析，根据模态参数，结合公式（）计算不同损伤工况下各个单元的平均模态应变能变化率，分布规律如图第期吴文斌 等基于单元模态应变能的梁结构分级损伤识别方法研究（）工况（）工况（）工况（）工况图各单元平均模态应变能变化率各工况下的损伤疑似单元见表：表各工况损伤单元及疑似损伤区域预设损伤单元确定的损伤单元损伤疑似区域工况，工况，工况，工况，由表可知，通过单元的各阶平均单元模态应变能变化率可以得出相应的损伤位置（工况、中均有确定的损伤位置）或者损伤疑似区域（个工况都有相应的损伤疑似区域），在损伤疑似区域内各单元的各阶平均单元模态应变能变化率不明显，即可能存在损伤单元，也可能不存在，因此要对该区域进行损伤二次定位 二次损伤定位分别对各工况损伤疑似区域的中间节点施加与整个损伤疑似区域质量相同的质量块，各工况质量块施加情况如下：（）对于工况，损伤疑似区域包含、号节点，在中间号节点施加与损伤疑似区域质量相等的质量块；（）对于工况，损伤疑似区域包含、号个节点，在中间、节点分别施加与损伤疑似区域质量相等的质量块；（）对于工况，损伤疑似区域包含、号个节点，在号节点施加与损伤疑似区域质量相等的质量块；（）对于工况，损伤疑似区域包含、号个节点，在中间、节点分别施加与损伤疑似区域质量的质量块质量块数值模拟施加方法：通过定义 质量单元，赋予相关节点质量属性，达到施加相应位河 北 建 筑 工 程 学 院 学 报第 卷置质量块的目的为控制单一变量，在分析种损伤工况时，还要同时对无损伤结构的对应位置施加节点质量，然后再次进行模态分析，只需提取损伤疑似区域及邻近单元（左右各一个单元，以便于损伤疑似区域的数据进行对比）的平均单元模态应变能变化率进行分析，所得结果如图（）工况（）工况（）工况（）工况图各单元平均单元模态应变能变化率由图知，在工况、各自的损伤疑似区域内，所确定的损伤位置分别为：工况，、号单元；工况，、号单元；工况，号单元；工况，号单元 最终损伤定位结果通过两次分析，最终确定损伤位置，见表：表两次损伤定位结果工况预设损伤单元一次识别二次识别（对损伤疑似区域加质量块）最终识别损伤单元，（），号单元 ，（）号单元为损伤疑似区域 ，（）号单元为疑似损伤区域号单元为损伤单元，（），号单元为损伤单元 号单元为损伤疑似区域，由表知，对于一次定位出现的损伤疑似区域施加质量块进行二次损伤定位，两次总的损伤定位结果与预设损伤位置相同，定位准确 确定损伤程度在确定损伤位置的情况下，对损伤程度进行判断，用遗传算法对损伤程度进行判断，设定初始种群为 ，交叉概率 ，变异概率 ，利用 与 联合仿真优化，以损伤单元在一次定位时的单元模态应变能为损伤量化值各工况损伤单元在一次损伤定位时的平均单元模态应变能见表：第期吴文斌 等基于单元模态应变能的梁结构分级损伤识别方法研究表损伤单元在一次定位时的平均单元模态应变能损伤单元一次定位时的平均单元模态应变能变化率工况，工况，工况，工况，设遗传算法中变量对应的损伤单元的单元模态应变能变化率为犻（犻，），则各工况的适应度函数表达式为式（）（）：犳 犻 狋狘（）（）狘（）犳 犻 狋狘（）（）狘（）犳 犻 狋狘（）（）狘（）犳 犻 狋狘（）（）（）狘（）式中：犳 犻 狋犳 犻 狋分别为工况的适应度函数值本算例迭代终止条件设为适应度函数值小于等于 ，各工况迭代进化曲线如图所示图迭代进化曲线个工况下利用遗传算法进行损伤程度的确定，得到各工况的迭代进化曲线终止代数、适应度最小值、对应的设计变量及损伤程度见表：表各工况损伤程度的量化分析表工况损伤单元终止代数适应度最小值对应的设计变量损伤程度（预设损伤程度），代 （预设），代 （预设），代 （预设），代 （预设）由表知，在确定损伤位置的情况下，利用遗传算法对损伤位置的损伤程度进行分析，在一定的遗传迭代后，能够判断出损伤位置损伤程度，并且与预设损伤情况相比较，该方法所识别的损伤程度误差均小于 ，接近预设损伤情况 结论通过以梁结构为研究对象，提出一种以单元模态应变能变化率为基础，通过施加质量块方法进行损伤精确定位，再结合遗传算法进行损伤程度判断的分级损伤识别方法，得出以下结论：（）在损伤位置判断方面，利用单元模态应变能变化率进行损伤区域的初步确定，再通过施加质量块的方法进行损伤位置的二次确定，最终得到结构的精确损伤位置河 北 建 筑 工 程 学 院 学 报第 卷（）在损伤程度确定方面，在损伤位置确定的基础上，利用遗传算法仅对损伤程度进行分析计算，解决了利用遗传算法进行损伤程度分析时由于虚假损伤造成计算量过大的问题参考文献 刘文光，颜龙，郭隆清基于模态应变能法的弹性薄板损伤识别噪声与振动控制，（）：缪炳荣，张盈，黄仲，张哲，杨树旺利用模态应变能变化率的结构损伤识别优化方法 振动工程学报：龚篧基于曲率模态的梁桥结构损伤识别方法研究重庆交通大学，张斌，苏云磊，范康基于模态应变能变化率结构损伤定位法的改进水利与建筑工程学报，（）：邓旭华，唐雪松基于模态应变能和 神经网络的混凝土框架结构损伤识别交通科学与工程，（）：孟海平，蔺新艳，杨健辉结构损伤定位中模态应变能法的改进广西大学学报（自然科学版），（）：姜绍飞，胡春明基于模态指标与数据融合的钢管混凝土拱桥损伤识别振动与冲击，（）：孙虎，狄生奎，杜祝遥，王立宪，项长生基于附加质量的单点模态振型参数损伤识别研究中外公路，（）：犛 狋 狌 犱 狔狅 狀犵 狉 犪 犱 犲 犱犱 犪 犿 犪 犵 犲 犾 狅 犮 犪 狋 犻 狅 狀犫 犪 狊 犲 犱狅 狀犲 犾 犲 犿 犲 狀 狋犿 狅 犱 犪 犾 狊 狋 狉 犪 犻 狀犲 狀 犲 狉 犵 狔犠犝 犠 犲 狀 犫 犻 狀，犢 犃 犖 犌犎 犪 犻 犼 狌 狀，犕犃犔 犲 犻，犑 犐 犃 犖 犌犢 犪 狓 犻 犪 狀，犠犝犛 犺 犪（，；，）犃 犫 狊 狋 狉 犪 犮 狋：，；，；，犓 犲 狔狑 狅 狉 犱 狊：；第期吴文斌 等基于单元模态应变能的梁结构分级损伤识别方法研究