一种基于混合元启发式算法的钢桁架结构损伤识别方法技术

本发明专利技术公开了一种基于混合元启发式算法的钢桁架结构损伤识别方法,属于元启发式算法在钢桁架损伤识别这一工程领域的应用,主要包括四个步骤：建立钢桁架损伤结构的有限元模型,获取结构在外荷载作用下的加速度；使用混合算法计算加速度；构建钢桁架结构的目标函数；不断优化目标函数直到满足终止条件输出最优解。本发明专利技术综合两种算法的优点,兼顾全局搜索与局部搜索的平衡,使该混合算法有很好的精确性和鲁棒性；该算法采用自适应变异和交叉算子和随迭代次数而改变的动态参数,可以识别出损伤的存在,损伤的位置和损伤的程度,而且在噪声环境下,本发明专利技术方法仍可以较为精确的识别钢桁架结构的多处损伤。

A damage identification method of steel truss structure based on hybrid element heuristic algorithm

【技术实现步骤摘要】
一种基于混合元启发式算法的钢桁架结构损伤识别方法
本专利技术涉及结构损伤识别
,具体为一种基于混合元启发式算法的钢桁架结构损伤识别方法。
技术介绍
随着科学技术的进步,大量新材料、新设备、新技术的应用使工程结构朝着大型化、复杂化的方向发展,上海金融中心、北京中国尊和深圳平安金融中心等标志性建筑表明我国在土木工程领域取得了巨大进步。由于工程结构的使用寿命长达几十年,在服役期内结构易受到材料老化、超载、环境侵蚀和疲劳破坏等不利因素的耦合作用,将引起结构承载力降低和损伤累积，可能导致结构破坏引发灾难性后果。2004年法国戴高乐机场2E候机楼结构在温差的反复作用下,结构损伤不断累积,承载力逐渐降低而坍塌；2007年美国明尼阿波利斯市的一座桥梁因钢材锈蚀发生坍塌事故；2008年南方特大冰雪灾害使钢结构超过其承载能力而大量破坏。钢桁架结构是常用的工程结构形式之一,为了保证结构在服役期内正常运行,保障人们财产和人身的安全,有必要对其进行健康监测。结构损伤识别技术是健康监测的核心,是工程领域的研究热点之一。结构损伤识别方法中,基于振动的方法是利用结构的加速度、频率、振型等动态参数对损伤进行检测,具有非破坏性、方便性和有效性等优点。钢桁架结构损伤识别逆问题可以转化为约束优化问题,先根据结构的动态参数定义优化问题的目标函数,然后利用优化算法计算出目标函数的最值。近些年来,大量优化算法已经应用在钢桁架结构损伤识别问题,传统的优化方法存在明显不足,或依赖于初始值,或易陷于局部最优而非全局最优；而受自然启发的元启发式算法则表现优越,全局搜索能力更强,既无需计算梯度,又易于实现并行计算。迄今为止,研究人员提出了多种自然启发的元启发式算法解决全局优化问题,常见的有人工蜂群算法(ABC)、模拟退火算法(SA)、遗传算法(GA)、粒子群优化算法(PSO)、差分进化算法(DE)等。其中,模拟蝴蝶觅食和交配行为的蝴蝶优化算法(ButterflyOptimizationAlgorithms，简称为BOA)是一种最新被提出并受到关注的智能优化算法,主要基于蝴蝶利用嗅觉确定花蜜位置和交配对象的策略。该算法设计参数少,概念简单且易于实施,收敛速度快；但蝴蝶优化算法也存在不足,全局搜索能力较差,有时陷入局部最优,向局部最优提前收敛而无法很好地执行全局搜索等。差分进化算法利用选择、变异和交叉三步操作更新种群个体，该算法全局搜索能力强,结构简单,易与不同算法混合。单独利用蝴蝶优化算法和差分进化算法进行钢桁架的损伤识别存在无法精确识别多处损伤和对噪声敏感的问题，影响了方法的实际应用。
技术实现思路
为克服上述方法的缺陷,本专利技术提出一种混合蝴蝶优化和差分进化的元启发式算法，把差分进化算法中的变异和交叉操作加入到蝴蝶优化算法的全局搜索阶段,平衡全局搜索与局部搜索的关系,加快收敛速度,避免陷入局部最优陷阱,该混合算法可探测出损伤的存在、损伤的位置和损伤的程度，解决钢桁架损伤识别的反问题。为实现上述目的,本专利技术采用的技术方案如下：一种基于混合元启发式算法的钢桁架结构损伤识别方法,具体包括以下步骤步骤1：将钢桁架划分为M个单元，L个加速度计安装在不同节点上,利用加速度计测量结构在外荷载force作用下的响应Amea；步骤2：混合算法参数初始化,在搜索域内随机产生结构参数,假定结构损伤为刚度的变化而质量不发生变化，结构单元刚度损伤程度可用损伤因子αi,i＝1,2,...,M，αi∈[0,1]，αi＝0表示钢桁架单元没有损伤,αi＝1表示钢桁架单元完全损伤,损伤单元的整体刚度矩阵为然后利用newmark-β法计算结构的加速度Acal；步骤3：构建钢桁架结构的目标函数,即待优化的目标函数为f(x)＝||Amea-Acal||；步骤4：利用混合蝴蝶优化和差分进化的元启发式算法不断优化目标函数,直到达到最大迭代次数或满足精确要求。进一步地,步骤4中钢桁架结构损伤识别的混合蝴蝶优化和差分进化的元启发式算法包括以下几个阶段：1)混合算法参数初始化,设置参数：该混合算法中蝴蝶代表钢桁架,钢桁架个数N,钢桁架单元数即算法维数Dim,损伤因子αi(i＝1,2,...,M)的搜索空间[Lb,Ub],钢桁架的最大迭代次数Gm；在搜索域内利用rand函数随机产生N钢桁架xi(i＝1,2,...,N),Lb≤xi≤Ub；定义待优化的目标函数f(x)；2)外荷载作用在钢桁架xi,计算得到的加速度Acal带入目标函数得到适应性函数f(x),找到最优个体Bf,最优个体Bf是最接近钢桁架实际损伤的识别结果。在钢桁架损伤识别问题中,适应性函数越来越小,修正有限元模型更接近钢桁架结构,从而实现钢桁架结构的损伤识别；3)计算蝴蝶感受到的香味,蝴蝶能够产生某种与其适应性相关的香味,即当蝴蝶从一个位置移动到另一个位置时,香味的浓度会随着适应性发生变化,香味的多少可利用式(1)计算：fi＝cIa(1)其中,fi是第i只蝴蝶感受到的香味大小，即第i个钢桁架的识别函数；I是香味刺激强度,由钢桁架结构的适应性函数f(x)决定；c是感觉形态,对算法的收敛速度有很大影响,由具体的钢桁架结构决定；a是幂指数,它影响香味的吸收程度,a＝1表示香味在理想环境下传播,香味没有损失,a＝0表示香味全部被吸收,其他蝴蝶不能感受到；设置感觉形态c、幂指数a的初始值为c0、a0。4)执行混合算法的迭代操作：4.1)对每一个钢桁架判断是否r＜p,其中r为[0，1]之间的随机数,p为开关概率,使算法在普通全局搜索和密集局部搜索之间切换,开关概率p表示环境因素对蝴蝶觅食交配行为的影响程度,影响钢桁架的识别结果,取值范围为[0，1]；4.2)执行局部搜索：如果r＜p,则蝴蝶朝着最优蝴蝶位置即适应性最好的解移动,即利用式(2)执行局部搜索阶段：上式中t为钢桁架当前迭代次数,为当前钢桁架的种群,为新生成的钢桁架的种群,g*为当前所有钢桁架识别解中最好的解,r为[0，1]之间的随机数,fi是第i只蝴蝶感受到的香味浓度。4.3)执行全局搜索：如果r≥p,蝴蝶无法感觉到周围的香味,它会随机移动,即利用式(3)执行全局搜索阶段：进一步地,为提高蝴蝶算法的全局搜索能力,在全局搜索阶段分别利用式(4)和(6)执行差分进化算法中的变异和交叉两步操作：式(4)中,为钢桁架的种群中的最优个体,为当前钢桁架的种群,为变异后的钢桁架的种群,F为自适应变异算子,取值范围为[0，2],可通过(5)式计算得到,Gm为钢桁架的最大迭代次数,变异算子初始值为F0。式(6)中,第t代钢桁架种群与其变异的中间体进行交叉操作生成新的钢桁架子代个体CR为服从正态分布的交叉算子,可通过CR(t)＝normrnd(CR0,0.1)得到,其中均值为CR0,标准差为0.1。5)更新参数,设置开关概率p、感觉形态c、幂指数a为随算法的迭代而改变的动态参数,具体计算式分别为(7)、(8)、(9)：本文档来自技高网...

【技术保护点】
1.一种基于混合元启发式算法的钢桁架结构损伤识别方法，其特征在于：具体包括以下步骤：/n步骤1：将钢桁架划分为M个单元，L个加速度计安装在不同节点上,利用加速度计测量结构在外荷载force作用下的响应A

【技术特征摘要】
1.一种基于混合元启发式算法的钢桁架结构损伤识别方法，其特征在于：具体包括以下步骤：
步骤1：将钢桁架划分为M个单元，L个加速度计安装在不同节点上,利用加速度计测量结构在外荷载force作用下的响应Amea；
步骤2：混合算法参数初始化,在搜索域内随机产生结构参数,假定结构损伤为刚度的变化而质量不发生变化，结构单元刚度损伤程度可用损伤因子αi,i＝1,2,...,M，αi∈[0,1]，αi＝0表示钢桁架单元没有损伤,αi＝1表示钢桁架单元完全损伤,损伤单元的整体刚度矩阵为然后利用newmark-β法计算结构的加速度Acal；
步骤3：构建钢桁架结构的目标函数,即待优化的目标函数为f(x)＝||Amea-Acal||；
步骤4：利用混合蝴蝶优化和差分进化的元启发式算法不断优化目标函数,直到达到最大迭代次数或满足精确要求。


2.根据权利要求1所述的一种基于混合元启发式算法的钢桁架结构损伤识别方法，其特征在于：步骤4中钢桁架结构损伤识别的混合蝴蝶优化和差分进化的元启发式算法包括以下几个阶段：
1)混合算法参数初始化,设置参数：该混合算法中蝴蝶代表钢桁架,钢桁架个数N,钢桁架单元数即算法维数Dim,损伤因子αi(i＝1,2,...,M)的搜索空间[Lb,Ub],钢桁架的最大迭代次数Gm；在搜索域内利用rand函数随机产生N钢桁架xi(i＝1,2,...,N),Lb≤xi≤Ub；定义待优化的目标函数f(x)；
2)外荷载作用在钢桁架xi,计算得到的加速度Acal带入目标函数得到适应性函数f(x),找到最优个体Bf,最优个体Bf是最接近钢桁架实际损伤的识别结果；在钢桁架损伤识别问题中,适应性函数越来越小,修正有限元模型更接近钢桁架结构,从而实现钢桁架结构的损伤识别；
3)计算蝴蝶感受到的香味,蝴蝶能够产生某种与其适应性相关的香味,即当蝴蝶从一个位置移动到另一个位置时,香味的浓度会随着适应性发生变化,香味的多少可利用式(1)计算：
fi＝cIa(1)
其中,fi是第i只蝴蝶感受到的香味大小，即第i个钢桁架的识别函数；I是香味刺激强度,由钢桁架结构的适应性函数f(x)决定；c是感觉形态,对算法的收敛速度有很大影响,由具体的钢桁架结构决定；a是幂指数,它影响香味的吸收程度,a＝1表示香味在理想环境下传播,香味没有损失,a＝0表示香味全部被吸收,其他蝴蝶不能感受到；设置感觉形态c、幂指数a的初始值为c0、a0；

【专利技术属性】
技术研发人员：周宏元，张广才，王小娟，倪萍禾，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11