钟昌廷, 李 刚
(大连理工大学 工程力学系 工业装备结构分析国家重点实验室,大连 116024)
结构可靠性分析是衡量结构不确定性和保证结构安全的重要手段。现有的可靠性分析方法主要包括抽样方法[1]、一阶可靠性方法[2]、代理模型[3]和矩方法[4]等。一阶可靠性方法简单高效,依照是否需要梯度信息,可分为梯度类方法和非梯度类方法。梯度类优化方法中,HL-RF迭代方法最具代表性[5,6],但在求解强非线性问题时会导致收敛困难。为了克服这一缺陷,研究者们提出了诸多改进策略。Yang等[7]提出了基于混沌理论的可靠度指标求解方法,并分析了HL-RF方法在寻优过程中出现震荡和分叉现象的机理。Gong等[8]提出了改进的有限步长方法。李刚等[9]提出了一种基于多目标优化策略的逆可靠度分析方法。夏雨等[10]提出了基于曲面搜索的方法求解结构可靠指标。上述基于梯度的方法取得了较大进展,但高维非线性隐式问题仍然是一大挑战,仍然可能面临不收敛的问题。
在非梯度类方法中,基于群智能算法的寻优策略逐渐受到研究者的关注。群智能优化算法是一类具有随机搜索能力的优化方法,具有无须梯度信息、多候选解和并行计算等特性,在工程结构优化中得到了广泛的应用[11,12]。目前,研究者已开发了多种群智能优化算法,用于改善一阶可靠性方法的寻优能力,如粒子群算法[13]、差分进化算法[14]、蛙跳算法[15]和碰撞体优化算法[16]等。钟昌廷等[17]提出了基于樽海鞘群优化算法和哈里斯鹰优化算法[18]的改进一次二阶矩法,对高维可靠度问题取得了一定的效果,但对高维工程算例的研究较少。因此需要考虑结构可靠度分析的特点,对现有群智能优化算法的搜索机制进行改进,提高其在高维非线性结构可靠度问题中的收敛能力。
海鸥优化算法是一种新型群智能优化算法[19],具有较好的求解优化问题的能力。Li等[20]将海鸥优化算法用于求解冷热电联供系统优化问题,Dhiman提出了[21]多目标海鸥优化算法求解多目标优化问题。然而,海鸥优化算法在结构可靠性分析中存在早熟和局部收敛等现象。本文引入生物地理优化算法中的迁移和变异策略[22],提出生物地理-海鸥群混合优化算法迭代求解结构可靠度指标,可有效求解高维非线性可靠度问题,并通过框架和网架工程结构验证了所提方法的有效性。本文提出基于生物地理-海鸥群优化的改进一阶可靠度方法,研究生物地理-海鸥群混合优化策略处理高维结构可靠度问题和搜索全局可靠度指标的性能,并采用两个高维工程算例验证了所提方法的有效性。
结构可靠性分析的目标是评估结构在不确定性下的失效概率,其表达式为
(1)
式中Pf为失效概率,fX(x1,…,xn)为基本随机变量x的概率密度函数,g(x) 为极限状态函数,小于0时的区域定义为失效域。由于式(1)为多维积分,难以直接求解,通常采用近似方法进行失效概率评估。一阶可靠度方法的主要任务为搜索最可能失效点,确定最优可靠度指标,在寻优过程中需要求解一个约束优化问题[5],
s.t.g(u)=0
(2)
式中u为标准正态空间内最可能失效点的位置向量,β为可靠度指标,定义为标准正态空间内极限状态面上离原点最近的距离。当随机变量不服从正态分布时,可采用当量正态化方法处理非正态随机变量[6],得到随机变量的均值和标准差。可以看出,式(2)是一个约束优化问题,其可靠度指标的求解质量受优化算法性能的影响很大。随着问题复杂程度的增加,对于优化算法的性能要求也越来越高,需要进一步开发具有更高性能的优化算法来求解高维非线性可靠度问题。
Dhiman等[19]提出一种新型的群智能优化算法,即海鸥优化算法SOA(Seagull optimization algorithm),其灵感来源于海鸥群体的迁徙和觅食行为。海鸥的迁徙行为主要包括躲避碰撞、朝最优位置方向和接近最优位置三种。在迁徙过程中,海鸥成群结伴地移动,为了防止海鸥之间碰撞误伤,海鸥相互保持一定距离,直至到达最佳位置。在攻击猎物时,空中飞行呈螺旋式运动,且盘旋半径逐渐减小直至捕食成功。根据上述行为,建立了海鸥优化算法的粒子位置更新机制。
海鸥在迁徙过程中,位置更新表达式为
Ds(t)=|AXi-2rA2(Xbest-Xi)|
(3)
式中Ds(t)为海鸥的更新位置,Xbest和Xi分别为海鸥群最优位置和海鸥当前位置,r为(0,1)的随机数,A为与迭代步t和最大迭代步Tmax相关的调节参数,计算式为A=2(1-t/Tmax)。
海鸥在攻击过程中,位置更新表达式为
Xi=Ds(t)exp(-2bL)cos(2πL)+Xbest
(4)
式中b为螺旋形状相关常数,L通过调节参数A计算得到,计算式为L=r(A-1)-1。
海鸥优化算法概念简单,易于执行,但求解结构可靠度问题时存在早熟和陷入局部最优的现象。其原因在于,海鸥优化算法的探索能力较差,一旦种群中的最优解陷入局部最优,其他个体解也跟随种群内最优解方向进行移动,其变异性较差。
为了改善海鸥优化算法的缺陷,本文引入生物地理优化算法BBO(Biogeography-based optimization)[22]的栖息地物种迁移和变异策略来改进海鸥优化算法的寻优能力。生物地理学优化算法中,栖息地视为候选解,而栖息地适应度指数视为适应度值。其位置更新策略如下。
在迁移过程中,迁入率和迁出率为描述栖息地变化过程中的两个参数,其计算公式为
λi=λmax(1-si/smax),μi=μmax(si/smax)
(5)
式中λi和μi分别为第i个栖息地的迁入率和迁出率,λmax和μmax分别为迁入率和迁出率的最大值。si和smax分别为栖息地的物种数量和最大值。通过计算迁移过程中的参数来确定解集更新分量,融合解集信息。
在变异过程中,其栖息地的突变概率与物种概率为反向关系,物种概率Pk计算公式为
(6)
式中栖息地Hk突变概率mk为
mk=mmax(1-Pk/Pmax)
(7)
式中mmax为栖息地的最大突变概率,Pmax为物种概率的最大值。
生物地理学优化算法的探索能力较强,可以弥补海鸥群优化算法探索能力的不足。本文结合两者算法的特点,提出了生物地理-海鸥群优化算法BBSOA(Biogeography-based seagull optimization algorithm),引入栖息地物种迁移和变异策略来增加种群的随机搜索能力,同时利用海鸥优化算法中较强的开发能力,对种群中每个个体生成随机数,并根据概率值判断更新策略,从而进行位置和适应度值更新。在混合群智能优化算法中,生物地理优化作为探索阶段,海鸥群优化作为开发阶段,因此平衡了优化算法的探索与开发能力,改善了其求解高维可靠性问题的收敛性能。
本文提出了基于生物地理-海鸥群优化的一阶可靠度方法BBSOA-FORM,其执行步骤为, (1) 基于一阶可靠性理论,建立可靠性分析中的优化列式; (2) 设置改进海鸥优化算法的种群数量和迭代步等参数; (3) 生成初始种群,并计算其适应度; (4) 更新算法参数与随机数,并根据概率值判断位置更新策略; (5) 若rand 图1 所提方法总流程 采用两个高维工程算例来验证本文算法,以蒙特卡洛(106次抽样)结果为参考解,将本文算法与粒子群算法、哈里斯鹰算法和海鸥群算法进行比较。所有算法种群数量取100,最大迭代步为500,本文算法Rd=0.6;统计20次独立运行结果。 考虑一个7层3跨的框架结构,各层层高3 m,柱间距7 m,如图2所示。该高维问题包含106个随机变量,其分布信息列入表1。极限状态函数为g(x)=0.08-d(x)。 图2 框架结构 表1 框架结构随机变量分布 蒙特卡洛抽样得到的可靠度指标和失效概率分别为β=3.8787和Pf=5.36×10-5。各群智能优化算法经过20次独立运行的可靠度指标统计结果列入表2,其中CPU时间为算法单次运行时间的平均值。结果表明,在500次迭代后,粒子群算法和海鸥群算法误差很大,没有收敛到正确结果;哈里斯鹰算法和本文方法与蒙特卡洛法结果相比,误差分别为5.78%和4.59%,而且本文方法的计算时间仅为哈里斯鹰算法的43%。另外,本文方法的标准差仅为0.0134,体现了良好的鲁棒性。需要指出的是,由于功能函数为高维非线性,随机变量正态分布,会导致蒙特卡洛法得到的失效概率转化为可靠指标有一定误差。各方法的收敛曲线如 图3 所示,本文所提方法的收敛速度与性能明显优于海鸥群算法和粒子群算法。 表2 框架结构可靠指标计算结果 图3 框架结构可靠指标收敛曲线 考虑一个空间网架穹顶结构,包含31个节点,75个杆件,平面最大跨度为20.153 m,高度为 6.448 m,其结构如图4所示。 图4 网架结构 结构受到作用于顶部1号节点的竖向垂直荷载。该高维结构可靠度问题包括77个随机变量,其随机变量统计信息列入表3。极限状态函数定义为g(x)=0.035-d(x)。 表3 网架结构随机变量分布 基于蒙特卡洛抽样的可靠度指标和失效概率计算结果为β=2.3325和Pf=0.00984。表4统计了各算法独立运行20次的计算结果。结果表明,其余三种算法的可靠度指标误差较大,本文方法的计算误差最小,且标准差仅为0.0109,具有良好的精度和鲁棒性。图5的收敛曲线表明,本文方法的收敛速度优于其他三种群智能优化算法。 表4 网架结构可靠指标计算结果 图5 网架结构可靠指标收敛曲线 本文提出了一种基于生物地理-海鸥群优化的改进一阶可靠度方法,并采用高维框架结构和网架结构算例验证了所提方法的有效性。其结论如下。 (1) 在混合优化算法中,采用生物地理优化中物种迁移和变异策略对海鸥优化算法进行改善,平衡了优化算法的探索与开发能力,使海鸥优化算法寻找全局最优可靠度指标的能力增强,可有效避免局部最优。计算过程中不需要梯度信息,便于处理梯度信息难以获取的可靠性问题。 (2) 本文方法在高维框架结构和高维网架结构可靠度算例进行了验证,并与粒子群算法、哈里斯鹰算法和海鸥群算法进行了比较。结果表明,本文提出的生物地理-海鸥群优化算法在计算精度、计算效率和鲁棒性等方面均有优势。3 算例验证
3.1 高维框架结构
3.2 高维空间网架结构
4 结 论