大型起重机械在服役过程中受到各种载荷作用, 同时由于材料老化、恶劣环境、自然灾害等多种因素的影响, 其结构性能不再完全等同于设计初期.随着损伤积累, 其承载能力也会逐渐下降.为了能动态在线监测结构的健康状况, 一个经过试验验证并能准确反映结构实际状态的有限元模型是极其必要的[1].
塔机的有限元模型修正原则上是一个多目标优化问题, 以基于模态参数的模型修正为例, 修正后的有限元模型仅重现某一阶试验模态参数是不够的, 一般应至少可以重现感兴趣频段内的各阶模态参数[2-4].
非支配排序遗传算法Ⅱ(non-dominated sorting genetic algorithm-Ⅱ, NSGA-Ⅱ)作为一种多目标进化算法, 不需要人为选择各子目标权重, 避免了人为选择主观影响.同时它采用快速非支配排序算法降低了计算复杂度, 引进精英策略, 保证某些优良的种群个体在进化过程中不会被丢弃, 从而提高了优化结果的精度, 采用拥挤度和拥挤度比较算子保证了种群的多样性.
本文提出适用于起重机械模型修正的NSGA-Ⅱ方法, 为减小基于遗传算法的模型修正的计算量, 提出用响应面模型代替原有的有限元模型, 将修正迭代中计算量巨大的有限元分析转移到计算速度更快的响应面分析, 以某型塔式起重机为工程实例完成了模型修正.
1 基于非支配排序遗传算法的多目标优化问题 1.1 优化问题数学建模用二次多项式响应面模型代替原有的有限元模型进行修正.二次多项式响应面模型包含多个待修正参数, 含有多个目标函数, 对模型修正过程实际是一个多目标优化问题.
多目标优化问题的数学模型一般可以表示为
(1) |
式中:x=[x1, x2, …, xn]为设计变量, xi, min≤xi≤xi, max(i=1, 2, …, n), xi, max和xi, min为第i个设计变量xi的上下限, n为设计变量的个数; r为子目标的个数; p为等式约束的个数; q为不等式约束的个数.
由于固有频率对结构的刚度变化比较敏感, 所以在有限元模型修正中, 采用固有频率残差来构建目标函数[5].目标函数由实际结构动力学实验得到的模态频率和由响应面模型计算得到的频率构成.
通过实验测试得到实际模型Mtrue, 其对应的前5阶模态频率值为ytrue=[f1, true, f2, true, f3, true, f4, true, f5, true].响应面模型为Mres, 由响应面计算出的前4阶模态频率为yres=[f1, res, f2, res, f3, res, f4, res].为了减小响应面模型计算所得的前4阶模态频率yres与实际模型计算值ytrue之间的差值, 就需要优化响应面模型Mres中的设计变量x.
根据修正模型构造的优化数学模型如下:
(2) |
这里, xi, min≤xi≤xi, max.
设修正完成后的有限元模型为Mopt, 由修正后的设计变量xtrue建立的有限元模型计算结果为yopt=[f1, opt, f2, opt, f3, opt, f4, opt, f5, opt].有限元模型修正的效果可以由修正后计算得到的前4阶模态频率yopt与实际模型测试所得模态频率ytrue之间的接近程度来衡量.为验证修正后有限元模型Mopt在测试频段外的预测能力, 将Mopt计算所得第5阶模态频率与实际模型Mtrue第5阶频率进行比较.整个模型的修正过程如图 1所示.
NSGA-Ⅱ[6]算法是Deb等于2000年在NSGA的基础上提出的, 该算法降低了非劣排序遗传算法的复杂性, 具有运行速度快、解集的收敛性好等优点, 其算法框图如图 2所示.
基于响应面方法[7]的有限元模型修正的效率和精度主要取决于试验设计、参数筛选及响应面函数形式的选择和拟合3个方面.本文采用的响应面函数是二次多项式[8].二次多项式是在响应面建模过程中广泛采用的一种函数形式, 但因为含交叉项的二次多项式系数多, 需要的样本点多, 相应的计算效率低[9], 致使目前使用比较广泛的是不考虑交叉项影响的二次多项式.不考虑交叉项影响的二次多项式形式如下:
(3) |
式中, 待定系数β0, βi和βii分别为响应面模型的常数项、一次项和二次项系数; ε为替代模型的误差.对于不考虑交叉项的二次型, 待定系数的个数为s=2n+1.
为了确定待定系数, 进行m次(m>s)独立试验, 每次试验中各设计变量的取值不同, 第j个样本点对应的响应值yj, 上标j表示第j次试验.
(4) |
将式(4)改写为矩阵形式:
(5) |
其中:y为m维系统输出值矢量; A为m×s维矩阵,
(6) |
采用最小二乘法求解式(4), 得到待定系数向量为
(7) |
根据式(3)拟合得到二次响应面有限元模型.
3 工程实例 3.1 初始有限元模型建立以某型塔机为工程对象, 根据塔机施工图纸中的材料参数、几何参数及边界条件等建立该塔机的初始有限元模型, 如图 3所示.
修正参数的选取对提高模型修正的质量和效率至关重要, 一个最基本的原则就是尽量选择那些对响应面影响比较显著的参数[10].在有限元模型修正中, 采用了固有频率残差来构建目标函数.对于该型塔机, 通过比较不同参数对各阶频率的灵敏度高低, 找出最灵敏的7个参数, 列于表 1.
以上述7个设计变量来构造响应面.应用最小二乘回归分析技术拟合的二次多项式响应面函数如式(8)~式(10)所示.
(8) |
(9) |
(10) |
(11) |
响应面的有效性用相对均方根误差RMSE和决定系数R2进行评价.RMSE值越接近0, R2越接近1, 响应面越有效.计算塔机前4阶频率二次响应面的RMSE和R2, 结果如表 2所示.
从表 2中可以看出, 塔机前4阶频率的二次响应面相对均方根误差RMSE最大仅为0.001 5, 而决定系数R2的最小值为0.993 1.由响应面的有效性评价准则可知, 在设计参数空间内, 拟合的响应面函数能以较高的精度反映结构响应和输入参数之间的关系, 回归的响应面模型可以作为有限元模型的替代模型进行模型修正.
为了验证不考虑二次交叉项的响应面模型不会影响计算精度, 同时计算了考虑交叉项后的相对均方根误差RMSE和决定系数R2, 结果如图 4所示.
图 4显示的结果说明:是否考虑二次交叉项对模型拟合精度影响不大, 不考虑交叉项的模型完全满足精度要求.但是由式(4)所示的响应面函数可以看出, 如果考虑二次交叉项会使目标函数项数大大增加, 从而增加了修正优化的计算时间, 降低修正优化的效率.
3.3 参数优化与模型修正采用基于NSGA-Ⅱ的多目标优化方法, 对二次响应面模型进行迭代优化,同时与实测风振响应数据得到的固有频率进行比较.其中实测数据通过布置在起重机上的传感器得到, 对结构进行激励, 采用模态识别技术识别出结构动力学特性.得到的低阶固有频率如表 3所示:
经过50次优化迭代, 得到了该多目标优化问题的Pareto最优解, 将所有的Pareto最优解向第4个子目标所在的平面投影, 即可得到如图 5所示的Pareto前沿的三维图.
为了更清晰地表现Pareto前沿, 再将三维的Pareto前沿向f1-f3和f2-f3两个平面投影, 分别得到二维的Pareto前沿如图 6所示.
每一个Pareto最优解对应一个修正后的有限元模型, 分别计算它们的响应值, 并求其统计值.结果如表 4和表 5所示, 从响应值中选取一个较优解列于最后一列.
将基于NSGA-Ⅱ算法的多目标优化结果与实际模型进行对比, 误差如表 6和表 7所示, 从表中可以看出, 使用多目标函数修正后的模型频率误差较小, 能够反映结构力学特性, 证明有限元模型修正是有效的.
1) 本文提出基于非支配排序遗传算法的机械结构有限元模型修正方法.将有限元模型修正建模为一个多目标优化问题, 无需人为主观选择各子目标的权重.
2) 利用二次响应面来代替遗传算法计算中的有限元模型, 大大提高了修正效率.
3) 完成了某型塔式起重机工程实例的有限元模型修正, 说明了本文所提算法的可行性和有效性.
[1] |
秦仙蓉, 徐俭, 赵坤, 等.
基于分层进化寻优的塔机结构有限元模型修正[J]. 同济大学学报(自然科学版), 2015, 43(6): 900–903.
( Qin Xian-rong, Xu Jian, Zhao Kun, et al. Dynamic model hierarchy updating of tower cranes using evolutionary computation[J]. Journal of Tongji University (Natural Science), 2015, 43(6): 900–903. ) |
[2] |
Mottershead J E, Friswell M I.
Model updating in structural dynamics:a survey[J]. Journal of Sound and Vibration, 1993, 167(2): 347–375.
DOI:10.1006/jsvi.1993.1340 |
[3] |
Brownjohn J M W, Dumanoglu A A, Severn R T.
Ambient vibration survey of fatih sultan mehmet (second Bosporus) suspension bridge[J]. Earthquake Engineering and Structural Dynamics, 1992, 21(10): 907–924.
DOI:10.1002/(ISSN)1096-9845 |
[4] |
李辉, 丁桦.
结构动力模型修正方法研究进展[J]. 力学进展, 2005, 35(2): 170–180.
( Li Hui, Ding Hua. Progress in model updating for structural dynamic[J]. Advances in Mechanics, 2005, 35(2): 170–180. DOI:10.6052/1000-0992-2005-2-J2003-124 ) |
[5] |
Arora V, Singh S P, Kundra T K.
Finite element model updating with damping identification[J]. Journal of Sound and Vibration, 2009, 324(3/4/5): 1111–1123.
|
[6] |
Deb K, Agrawal S, Pratap A, et al. A fast elitist non-dominated sorting genetic algorithm for multi-objective optimization: NSGA-Ⅱ[C]//International Conference on Parallel Problem Solving from Nature. Berlin: Springer-Verlag, 2000: 849-858.
|
[7] |
Ghrib F, Li L, Wilbur P.
Damage identification of Euler-Bernoulli beams using static responses[J]. Journal of Engineering Mechanics, 2012, 138(5): 405–415.
DOI:10.1061/(ASCE)EM.1943-7889.0000345 |
[8] |
Ma Y, Zhang J, Wang L, et al.
Probabilistic prediction with Bayesian updating for strength degradation of RC bridge beams[J]. Structural Safety, 2013, 44(2334): 102–109.
|
[9] |
Khuri A I.
Response surface methodology[M]. Berlin: Springer, 2011.
|
[10] |
Govers Y, Link M.
Stochastic model updating-covariance matrix adjustment from uncertain experimental modal data[J]. Mechanical Systems & Signal Processing, 2010, 24(3): 696–706.
|