基于BMO算法的结构损伤识别方法技术

本发明专利技术主要涉及鸟群遗传优化（BMO）算法在结构损伤识别这一领域的工程应用，具体是一种基于BMO算法的结构损伤识别方法，主要步骤如下：（一）通过有限单元法建立损伤结构的有限元模型，提取结构的固有频率、振型等模态参数。（二）利用损伤结构和计算结构的固有频率残差和模态确保准则构建目标函数(MAC)构建目标函数。（三）采用BMO算法优化这一目标函数，直到满足循环结束条件为止。（四）最后得到的最优解即为损伤识别结果。该方法相较于传统的灵敏度方法而言，无需借助梯度信息，利用少量的模态参数即可得到精度较高的识别结果。

【技术实现步骤摘要】

本专利技术涉及结构健康检测损伤识别
，更具体地，涉及一种利用遗传优化算法——鸟群(BMO)算法和结构的模态参数对结构的损伤进行识别的方法，是一种基于BMO算法的结构损伤识别方法。
技术介绍
近年来，我国经济迅速发展，科学技术水平也在不断提高，在许多领域都取得了很不错的成绩。与此同时，由于科技的提高和经济的发展，国家在重大土木工程诸如大跨度桥梁、大坝、大垮空间结构、超高层建筑、海洋平台等方面投入了巨额的资金。这些大型建筑大多数都要服役几十年，有的甚至达到上百年，在这么长的服役时间内，这些建筑不免要受到环境侵蚀、突变效应、载荷效应和材料老化等因素的影响，在这些因素的耦合作用下，用于构建这些大型建筑的材料和结构会逐渐的出现损伤，它们的抗力也会不断地衰减。当结构的关键部位损伤达到一定程度的时候，这些建筑就会出现无法挽回的损伤，给国家财产带来重大损失，给人身安全带来重大威胁，为了有效的预防类似情况的发生，结构损伤识别技术应运而生。结构的损伤可以发生在很多的方面，判断一个结构是否损伤也有不同的衡量参数和标准，这就导致了结构损伤识别方法的多种多样，但是大体上是两个方面：一个是基于振动的损伤识别检测方案，另一个是基于动力特性的损伤检测方案。文献“基于时域响应灵敏度分析的板结构损伤识别(振动与冲击，2015，34(4)，117～120)”提出了一种模型修正方法和灵敏度方法相结合的损伤识别新方法。该方法首先利用New-mark法获得损伤结构的时域响应，在损伤识别反问题当中，利用灵敏度分析，不断进行迭代，最终得到最后的识别结果。然而应用到时域数据时，要求测量一定时间内的响应数据，所以测量的数据点相对较多，而且这些数据很容易被噪声“污染”，进而影响方法的实际应用。
技术实现思路
为了改进现有技术的缺陷，本专利技术提出一种基于BMO算法的结构损伤识别方法，该方法采用频域数据对损伤结构进行识别并且利用BMO遗传优化算法对目标函数进行优化，得到损伤识别结果。该方法检测只需要借助前几阶模态参数就可以实现结构损伤识别，具有较高的精度。为了实现上述目的，本专利技术的技术方案为：一种基于BMO算法的结构损伤识别方法，具体步骤如下：步骤一：将结构划分为nel个单元，再利用有限单元法得到系统刚度和质量矩阵，再提取前N阶固有频率和模态。步骤二：构建损伤结构的目标函数，即待优化的目标函数。步骤三：利用BMO算法不断优化目标函数，直到满足终止条件。进一步，算法包括如下几个阶段1)初始化参数，包括算法的初始种群数量、最大迭代次数、算法中的鸟群不同繁殖方式所占的比例以及五种不同鸟类繁殖方式的算法。鸟群的繁殖方式包括以下五种，分别为：一夫一妻制的鸟类繁殖方式、一夫多妻制的鸟类繁殖方式、一妻多夫制的鸟类繁殖方式、孤雌生殖与杂交繁殖；其中一夫一妻制的鸟类后代产生方式为： x → b = x → + w × r → . × ( x → i - x → ) ]]>如果r1＞mcf，则xb(c)＝l(c)-r2×(l(c)-u(c))；其中c＝1,2,....,n其中，表示所得到的后代，w表示一个随时间变化的权重因子，表示的是一个随机向量，并且该向量中的元素取值为0到1之间的随机数，表示任意一夫一妻制的雄性鸟类，表示一夫一妻制的雌性鸟类，xb(c)是向量中对应第c个位置的数，n表示的鸟个体的总数，mcf是一个在0到1之间随机变化的突变控制因素，r1、r2是0到1之间的一个随机数，而u(c)和l(c)分别表示各单元的最大值和最小值；令一夫多妻制与一妻多夫制的鸟只产生一个后代，基因为所有的雌鸟基因混合算法，其表达式为下： x → b = x → + w × Σ j = 1 n i r → . × ( x → i j - x → ) ]]>如果r1＞mcf，则xb(c)＝l(c)-r2×(l(c)-u(c))，其中c＝1,2,....,n其中，表示所得到的后代，w表示一个随时间变化的权重因子，表示的是一个随机向量，并且该向量中的元素取值为0到1之间的随机数，为表示鸟类个体的向量，ni表示有吸引力的鸟总数，而表示的是第j个有吸引力的鸟；xb(c)是向量中对应第c个位置的数，n表示的鸟个体的总数，mcf是一个在0到1之间随机变化的突变控制因素，r1、r2是0到1之间的一个随机数，而u(c)和l(c)分别表示各单元的最大值和最小值；孤雌生殖与杂交繁殖中，各雌性都通过改变自身的基因使得后代的基因更优良，其算法给出如下：如果r1＞mcfp，则xb(i)＝x(i)+μ×(r2-r3)×x(i),否则xb(i)＝x(i)其中，x(i)为表示对应向量的初始数据，xb(i)是向量中对应第i个位置变化以后的数，n表示的鸟个体的总数，mcf是一个在0到1之间随机变化的突变控制因素，r1、r2、r3是0到1之间的一个随机数，μ表示步长，i＝1，2，…，n。2)计算种群的函数适应度值，并评价种群在损伤识别问题中，适应度的计算公式如下： f = Σ j = 1 N F w ω j 2 Δω j 2 + Σ j = 1 N M w 本文档来自技高网...

【技术保护点】
一种基于BMO算法的结构损伤识别方法，其特征在于，包括以下步骤：1)用有限元方法将结构进行简化建模，并把结构划分为nel个单元；2)提取损伤结构的NF阶频率和模态，构建目标函数f如下所示：f=Σj=1NFwωj2Δωj2+Σj=1NMwφj2(1-MACjR)]]>Δωj=|ωjC-ωjMωjM|]]>MACjR=({φjC}T{φjM})2|{φjC}|2|{φjM}|2]]>其中为结构第j阶计算得到的频率和振型，为结构第j阶测量得到的频率和振型，为权重系数；Δwj表示结构固有频率差值，代表第j阶不完整振型对应的简化模态置信准则，NF和NM分别为选用的固有频率和振型的个数；3)利用BMO算法不断优化目标函数，直到满足设定的终止条件，得到识别结果，其具体过程为：初始化参数，包括算法的初始种群数量、最大迭代次数、算法中的鸟群不同繁殖方式所占的比例以及五种不同鸟类繁殖方式的算法；在初始化后计算种群的函数适应度值，并评价种群；选出适应度为最佳，即对应折损系数最大的后代；选择一次迭代以后符合条件的所有优秀后代，用后代替代掉原来种群中不好的后代，如此循环，不断地将解进行优化；记忆目前最好的解，直到算法结束为止。...

【技术特征摘要】
1.一种基于BMO算法的结构损伤识别方法，其特征在于，包括以下步骤：1)用有限元方法将结构进行简化建模，并把结构划分为nel个单元；2)提取损伤结构的NF阶频率和模态，构建目标函数f如下所示： f = Σ j = 1 N F w ω j 2 Δω j 2 + Σ j = 1 N M w φ j 2 ( 1 ...

【专利技术属性】
技术研发人员：吕中荣，秦昌富，朱嘉健，
申请(专利权)人：中山大学，
类型：发明
国别省市：广东;44