文章快速检索  
  高级检索
基于加速度频响函数小波分解的模型修正方法
彭珍瑞, 曹明明, 刘满东     
兰州交通大学 机电工程学院, 兰州 730070
摘要: 为提高模型修正效率,满足修正方法对于实测环境噪声的鲁棒性,将Kriging模型和小波分解引入加速度频响函数模型修正。首先,将加速度频响函数进行小波分解,用得到的第1层幅值较大的小波系数来表征原频响函数。其次,采用拉丁超立方抽样对初选待修正参数进行设计,根据设计结果对各参数进行灵敏度分析,从而确定模型修正的待修正参数,以待修正参数作为Kriging模型输入,所对应的小波系数作为Kriging模型输出,通过混合灰狼算法寻得最优Kriging模型相关系数,建立精确有效的Kriging模型。最后,以目标加速度频响函数小波分解的小波系数与Kriging模型输出的小波系数差值最小为目标,通过水循环算法求解模型待修正参数。数值算例表明,所提模型修正方法具有良好的修正效果,即使在加速度频响函数中加入信噪比为5 dB的高斯白噪声时,修正误差也低于4%,证明了该方法对于随机噪声的鲁棒性。
关键词: 模型修正    Kriging模型    加速度频响函数    小波分解    相关系数    
Model updating method based on wavelet decomposition of acceleration frequency response function
PENG Zhenrui, CAO Mingming, LIU Mandong     
School of Mechanical Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
Abstract: To improve the efficiency of model updating and satisfy its robustness for the measured environmental noise, the Kriging model and wavelet decomposition are introduced into the model updating of the acceleration frequency response function. Firstly, the acceleration frequency response function is decomposed by wavelets, and the obtained wavelet coefficients with large amplitudes in the first layer are used to represent the original frequency response function. Secondly, the Latin hypercube sampling is utilized to design the primary parameters to be updated, and the sensitivity analysis of each parameter carried out according to the design results to determine the parameters to be modified, which are then used as the inputs of the Kriging model, while the corresponding wavelet coefficients as the outputs of the model. The optimal correlation coefficients of the Kriging model are found through the mixed grey wolf algorithm, with which an accurate and effective Kriging model is established. Finally, with the error between the wavelet coefficients calculated from the Kriging model and those from tests as the objective function, a minimization problem is solved by the water cycle algorithm for parameter updating. Numerical examples show the effectiveness of the proposed model updating method. When Gaussian white noise with a signal-to-noise ratio of 5 dB is added to the acceleration frequency response function, the updating error is smaller than 4%, proving the robustness of the method against random noise.
Keywords: model updating    Kriging model    acceleration frequency response function    wavelet decomposition    correlation coefficients    

在航空、航天、机械、土木结构和武器系统的研制过程中,有限元仿真和试验分析已成为2种不可或缺的重要手段[1]。但是在建模过程中的许多不确定性因素会导致有限元模型与试验模型存在一定的差异,因此需要利用试验模型的测试数据对有限元模型进行修正,使修正后的有限元模型能较准确地反映结构的动态特性[2]。其中,基于模态参数和基于频响函数(FRFs)的模型修正方法较为常见[3]。相比较而言,基于频响函数的模型修正方法直接利用测量得到的频响函数进行模型修正不需要模态识别,从而避免了识别误差的引入,具有一定的优越性[4]。文献[5-6]分别从理论上介绍了设计参数型频响函数修正法在工程应用中的优势及局限。殷红等[7]利用位移频响函数结合结构损伤识别,得到了与结构实际参数相符的模型修正结果。Esfandiari等[8]采用基于频响函数的模型修正方法,实现了对一标准桁架结构的损伤识别。杨海峰等[9]使用加速度频响函数对一简支梁进行了有限元模型修正。虽然基于频响函数的模型修正方法有一定的优势,但因其包含大量频率点信息,所以频率点的选择对于修正结果有重要影响,目前尚无统一方法。Kwon和Lin[10]引入测量依赖性指标和参数依赖性指标来选择最优频率点。张勇等[11]利用频响函数截断奇异值分解进行模型修正避开了频率点的选择难题。

目前,作为一种可以替代有限元模型进行快速分析的代理模型技术,引起了众多学者的关注。代理模型通过建立一个显式的函数来代替模型参数与结构响应之间复杂的隐式关系,近年来逐渐被应用到结构模型修正与确认中[12]。常用的代理模型有响应面法(RSM)、径向基函数(RBF)、神经网络(NN)、支持向量回归机(SVR)、Kriging模型等[13]。其中Kriging模型是基于Kriging插值技术的一种等效模型,与其他代理模型不同的是,Kriging模型除了能给出对未知函数的预估值,还能得到预估值的误差估计[14],只用少量的样本就可以较准确地描述结构输入与输出间的关系,在航天飞行器设计和结构优化设计等领域应用较为广泛。文献[15-17]将Kriging方法应用于频响函数模型修正中,减少了求解时间,提高了频响函数模型修正的效率。

在上述背景下,为保留利用频响函数进行模型修正的优点,同时避免频率点的选择难题和提高模型修正效率,本文提出了一种基于加速度频响函数小波分解的模型修正方法。将计算得到的加速度频响函数进行小波分解,用小波分解得到的第1层幅值较大的小波系数来表征频响函数作为Kriging模型的输出,以模型待修正参数作为Kriging模型的输入,并通过混合灰狼算法对Kriging模型相关系数进行寻优来构建精确的Kriging代理模型。然后,以实测的加速度频响函数小波分解得到的第1层小波系数与代理模型输出小波系数的差值最小作为目标函数,通过水循环算法进行待修正参数的求解。最后,通过一空间桁架结构证明了本文方法的有效性,特别是对于随机噪声的鲁棒性。

1 加速度频响函数小波分解 1.1 加速度频响函数

频域内,一个n自由度阻尼系统在简谐激励作用下的运动方程可以表示为

$ ( - {\omega ^2}\mathit{\boldsymbol{M}} + {\rm{i}}\omega \mathit{\boldsymbol{C}} + \mathit{\boldsymbol{K}})\mathit{\boldsymbol{X}}(\omega ) = \mathit{\boldsymbol{F}}(\omega ) $ (1)
 

式中:MKC分别表示质量矩阵、刚度矩阵和阻尼矩阵;ω为激励频率;F(ω)为简谐激励;稳态响应X(ω)可表示为

$ \mathit{\boldsymbol{X}}(\omega ) = \mathit{\boldsymbol{H}}(\omega )\mathit{\boldsymbol{F}}(\omega ) $ (2)
 

其中:H(ω)为加速度频响函数,可表示为

$ \mathit{\boldsymbol{H}}(\omega ) = {\left( {\mathit{\boldsymbol{M}} - \frac{{{\rm{i}}\mathit{\boldsymbol{C}}}}{\omega } - \frac{\mathit{\boldsymbol{K}}}{{{\omega ^2}}}} \right)^{ - 1}} $ (3)
 
1.2 小波变换

小波变换是信号处理领域的一种重要研究方法,被用于图像压缩、信号滤波、数据融合和特征提取等很多方面。Donoho[18]提出了小波软、硬阈值消噪法。该方法计算量小且容易实现,得到了广泛的关注和应用。

一维含噪信号可以表示为

$ y(t) = s(t) + o(t)\;\;\;{\kern 1pt} t = 0,1, \cdots ,T - 1 $ (4)
 

式中:y(t)为含噪信号;s(t)为有用信号;o(t)为噪声信号;t为采样点;T为采样点数,其小波变换为

$ {\lambda _{jk}} = {2^{ - \frac{j}{2}}}\sum\limits_{n = 0}^{T - 1} y (n)\psi ({2^{ - j}} - k)\;\;\;{\kern 1pt} j,k \in {\rm{Z}} $ (5)
 

其中:λjk为含噪信号的小波系数,ψ(·)为小波函数。因为小波变换是线性变换,所以式(5)可以变为

$ {\lambda _{jk}} = {\lambda _{jk}}(s) + {\lambda _{jk}}(o) $ (6)
 

其中:λjk(s)和λjk(o)分别为有用信号和噪声信号的小波系数。

小波变换特别是正交小波变换具有很强的去数据相关性,它能够在小波域中使信号的能量集中在一些大的小波系数中;而噪声的能量却分布于整个小波域内。因此,经小波分解后,有用信号的小波系数幅值要大于噪声信号的系数幅值。可以认为,幅值比较大的小波系数一般以有用信号为主,而幅值比较小的系数在很大程度上是噪声信号[19]。对比多次试验效果,设定小波基为db4,当无噪声或者弱噪声(信噪比为30 dB、20 dB)干扰时分解层数设为5,当强噪声(信噪比为10 dB、5 dB)干扰时分解层数设为4,用加速度频响函数经小波分解后的第1层幅值较大的小波系数作为加速度频响函数的特征量进行模型修正。

2 基于Kriging模型的模型修正理论 2.1 Kriging模型的建立

Kriging模型是一个基于随机过程的代理模型,最初在地质统计学中提出。模型包含了线性回归部分和非参数部分:

$ y(\mathit{\boldsymbol{x}}) = \sum\limits_{i = 1}^p {{\beta _i}} {f_i}(\mathit{\boldsymbol{x}}) + z(\mathit{\boldsymbol{x}}) = {\mathit{\boldsymbol{f}}^{\rm{T}}}(\mathit{\boldsymbol{x}})\mathit{\boldsymbol{\beta }} + z(\mathit{\boldsymbol{x}}) $ (7)
 

式中:β=[β1 β2 … βp]T为回归模型系数;f(x)=[f1(xf2(x) … fp(x)]T为样本点x的多项式函数, p为回归多项式的数量;z(x)为协方差非零的随机分布,服从正态分布N(0, σ2),其协方差为

$ {\rm{ Cov}} [z({x_i}),z({x_j})] = {\sigma ^2}R({x_i},{x_j}) $ (8)
 

其中:σ2z(x)的方差;xixj为样本中任意2个样本点;R(xi, xj)是相关函数,表示样本点间的空间相关性。选计算效率较好的高斯函数为相关函数:

$ R({x_i},{x_j}) = {\rm{exp}}\left( { - \sum\limits_{k = 1}^m {{\theta _k}} |x_i^k - x_j^k{|^2}} \right) $ (9)
 

式中:$x_i^k$$x_j^k$分别为样本点xixj的第k个分量; m表示设计变量的个数;θk为相关系数。利用最小二乘估计法,可得β估计值为

$ \mathit{\boldsymbol{\hat \beta }} = {({\mathit{\boldsymbol{\hat F}}^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}\mathit{\boldsymbol{\hat F}})^{ - 1}}{\mathit{\boldsymbol{\hat F}}^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}\mathit{\boldsymbol{Y}} $ (10)
 

式中:${\mathit{\boldsymbol{\hat F}}}$为样本点的向量组成的矩阵;Y为样本点响应组成的列向量;R为空间相关矩阵,其中元素Rij=R(xi, xj); i, j=1, 2, …, n; σ2的估计值为

$ {\hat \sigma ^2} = \frac{1}{n}{(\mathit{\boldsymbol{Y}} - \mathit{\boldsymbol{\hat F\beta }})^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}(\mathit{\boldsymbol{Y}} - \mathit{\boldsymbol{\hat F\beta }}) $ (11)
 

其中:n为试验点数;σ2β均为θk的函数,那么Kriging模型中唯一未知数即为相关系数θk,其决定着预测响应精度,本文采用混合灰狼算法对相关系数θk进行寻优。

Kriging模型构造完成后,下一步就是预测待测点的响应值。对于任意待测点x0,其响应值${\hat y}$(x0)为

$ \hat y({x_0}) = {\mathit{\boldsymbol{f}}^{\rm{T}}}({x_0})\mathit{\boldsymbol{\hat \beta }} + {\mathit{\boldsymbol{r}}^{\rm{T}}}({x_0}){\mathit{\boldsymbol{R}}^{ - 1}}(\mathit{\boldsymbol{Y}} - \mathit{\boldsymbol{\hat F\hat \beta }}) $ (12)
 

式中:rT(x0)=[R(x0, x1R(x0, x2) …R(x0, xn)]。

2.2 混合灰狼算法

混合灰狼算法(DE-GWO)是由Zhu等[20]提出,该算法将差分进化算法(DE)强大的搜索能力融入灰狼算法(GWO)中,使GWO能够跳出停滞区,降低了GWO寻得局部最优值的可能性,加快了收敛速度,优于粒子群(PSO)算法、差分进化算法(DE)和灰狼算法(GWO)。算法流程如下:

1) 参数初始化。设置种群规模Npop、最大迭代次数MaxIt、自变量维数nVar、缩放因子下上界beta_min和beta_max,交叉概率pCR,参数下上界lb和ub。

2) 对种群个体进行DE变异操作,产生中间体;并进行竞争选择操作形成父代种群、子代种群和变异种群个体。

3) 社会等级。计算种群中每个灰狼个体的适应度值,并依据适应度值的大小进行排序,确定灰狼种群中的社会等级,最优解Xα为头狼,第2和第3最优解XβXδ为狼群中第2和第3等级狼,其余的候选解X为狼群中最低等级狼。

4) 搜索包围猎物。灰狼搜索猎物过程中,灰狼接近并包围猎物行为可以表示为

$ {\mathit{\boldsymbol{D}} = |\mathit{\boldsymbol{C}} \cdot {\mathit{\boldsymbol{X}}_p}(t) - \mathit{\boldsymbol{X}}(t)|} $ (13)
 
$ {\mathit{\boldsymbol{X}}(t + 1) = {\mathit{\boldsymbol{X}}_p}(t) - \mathit{\boldsymbol{A}} \cdot \mathit{\boldsymbol{D}}} $ (14)
 

式中:D为灰狼与猎物之间的距离;Xp(t)和X(t)分别是第t次迭代猎物和灰狼的位置向量;CA为系数向量,且C=2r1A=2ar2a,其中r1r2为[0, 1]之间的随机向量,a是一个动态向量,随着迭代次数的增加从2线性递减至0。

5) 追捕猎物。当灰狼搜索到猎物所在位置时,首先,头狼Xα带领狼群对猎物进行包围;然后,头狼Xα带领Xβ狼和Xδ狼对猎物进行攻击捕捉。在灰狼群体中,XαXβXδ狼距离猎物位置最近,因此可以通过这三者的位置来计算灰狼个体向猎物移动的位置。并根据式(15)计算出种群中最优个体XαXβXδ与其他灰狼个体之间的距离,再根据式(16)计算出其余灰狼的移动方向。最后由式(17)更新灰狼位置。

$ \left\{ \begin{array}{l} {\mathit{\boldsymbol{D}}_\alpha } = |{\mathit{\boldsymbol{C}}_1} \cdot {\mathit{\boldsymbol{X}}_\alpha }(t) - \mathit{\boldsymbol{X}}(t)|\\ {\mathit{\boldsymbol{D}}_\beta } = |{\mathit{\boldsymbol{C}}_2} \cdot {\mathit{\boldsymbol{X}}_\beta }(t) - \mathit{\boldsymbol{X}}(t)|\\ {\mathit{\boldsymbol{D}}_\delta } = |{\mathit{\boldsymbol{C}}_3} \cdot {\mathit{\boldsymbol{X}}_\delta }(t) - \mathit{\boldsymbol{X}}(t)| \end{array} \right. $ (15)
 
$ \left\{ \begin{array}{l} {\mathit{\boldsymbol{X}}_1}(t) = {\mathit{\boldsymbol{X}}_\alpha }(t) - {\mathit{\boldsymbol{A}}_1} \cdot {\mathit{\boldsymbol{D}}_\alpha }\\ {\mathit{\boldsymbol{X}}_2}(t) = {\mathit{\boldsymbol{X}}_\beta }(t) - {\mathit{\boldsymbol{A}}_2} \cdot {\mathit{\boldsymbol{D}}_\beta }\\ {\mathit{\boldsymbol{X}}_3}(t) = {\mathit{\boldsymbol{X}}_\delta }(t) - {\mathit{\boldsymbol{A}}_3} \cdot {\mathit{\boldsymbol{D}}_\delta } \end{array} \right. $ (16)
 
$ \mathit{\boldsymbol{X}}(t + 1) = \frac{{{\mathit{\boldsymbol{X}}_1}(t) + {\mathit{\boldsymbol{X}}_2}(t) + {\mathit{\boldsymbol{X}}_3}(t)}}{3} $ (17)
 

式中:X1(t)、X2(t)和X3(t)代表最优个体XαXβXδ与其他灰狼个体X指导后更新的位置;X(t+1)为子代灰狼的最终寻优位置。

6) 对种群个体进行交叉、选择操作保留优良成分并产生新子代个体,计算个体的适应度值。

7) 更新灰狼XαXβXδ及其他对应的位置信息。

8) 判断是否达到最大迭代次数MaxIt,如果是则停止并输出当前最优解,否则返回3)继续执行。

2.3 混合灰狼算法优化Kriging模型

Kriging模型中的相关系数θk决定着预测响应值精度,只有当所建立Kriging模型精度足够高时,代理模型对修正结果的误差影响才会最小。所以相关系数的确定对于构建代理模型至关重要。本文采用拉丁超立方抽样方法抽取一定数量的样本点,然后把抽取的样本点分成训练集和测试集。以训练集作为Kriging模型的输入变量,以加速度频响函数小波分解得到的第1层小波系数作为响应值来构建Kriging模型,以测试集Kriging模型的均方误差(MSE)均值平均值作为混合灰狼算法的目标函数,寻得最优相关系数。最后以训练集作为样本点建立Kriging模型。

2.4 模型修正

当Kriging模型构造完成时,模型修正归根结底是一优化问题,在目标函数最小的条件下,通过求解该优化问题,得到设计参数的修正值。当激励点和响应点确定时,目标函数为

$ {\rm{obj}} = \sum\limits_{i = 1}^K {{{({{\hat \lambda }_i} - {\lambda _i})}^2}} $ (18)
 

式中:K为加速度频响函数小波分解得到的第1层小波系数个数;${{\hat \lambda }_i}$为Kriging模型的第i个预测响应小波系数;λi为试验模型的第i个真实响应小波系数。

在优化求解时,除了灵敏度方法外,智能算法也是很简便的。水循环算法(WCA)[21]的基本概念和思想来源于自然,并基于现实世界中的水循环过程,具有良好的随机搜索能力、鲁棒性和快速收敛等优点。故本文使用水循环算法进行优化求解待修正参数,模型修正流程如图 1所示。

图 1 模型修正流程图 Fig. 1 Flow chart of model updating
3 数值算例

选择一空间桁架结构来验证本文方法。该结构包括28个节点,66个单元和48个自由度。桁架约束条件为4个支座固定(节点编号:1、8、9、16),节点铰接,如图 2所示,桁架结构参数如表 1所示。

图 2 三维桁架模型结构图 Fig. 2 Structure drawing of three-dimensional truss model
表 1 三维桁架结构参数 Table 1 Structure parameters of three-dimensional truss
杆类型 弹性模量/GPa 横截面积/mm2 密度/(kg·m-3)
上弦杆 190 85.5 7 800
下弦杆 190 85.5 7 800
直腹杆 190 141.0 7 800
斜腹杆 190 45.0 7 800

根据样本点参数,运用动力学方法计算其质量矩阵、刚度矩阵并且添加一定的比例阻尼,将其代入式(3)中计算出其加速度频响函数矩阵。由于本文所用方法间接利用频响函数,避开了频率点选择问题,所以在频率点的选择问题上,只需保证频响函数在感兴趣频带内有足够数量的共振峰。所选结构的激励点和响应点如图 2所示。

3.1 待修正参数选取与试验设计

根据表 1将各杆的弹性模量E和密度ρ、上弦杆横截面积A1、下弦杆横截面积A2、直腹杆横截面积A3、斜腹杆横截面积A4,这6个参数作为初选待修正参数。以表 1中的各参数数值作为“试验模型”参数,然后对“试验模型”参数进行偏移得到“有限元模型”参数。其中,弹性模量E增加10%,密度ρ减小10%,上弦杆横截面积A1增加10%,下弦杆横截面积A2减小10%,直腹杆横截面积A3增加10%,斜腹杆横截面积A4减小10%,得到对应的有限元模型参数值如表 2所示

表 2 试验模型和有限元模型参数 Table 2 Parameters of test model and finite element
参数 试验值 初始误差/% 有限元值
E/GPa 190 10 209
ρ/(kg·m-3) 7 800 -10 7 020
A1/mm2 85. 5 10 95
A2/mm2 85. 5 -10 76
A3/mm2 141 10 156
A4/mm2 45 -10 40

对于全局优化问题,比较好的方法是通过试验设计(DoE)来选取一组样本点。对于Kriging模型,初始样本点数理论上不受设计空间维数的限制,且优化效率对初始样本点数的依赖也不明显[14]。为此,采用拉丁超立方抽样在上述有限元模型6个参数的±20%区间内抽取200个样本点,使所构建Kriging模型满足精度要求。根据样本点运用全局灵敏度中的Sobol法计算各参数对于响应的灵敏度值,将各参数灵敏度值归一化并转换成百分比形式,如图 3所示。

图 3 各参数灵敏度 Fig. 3 Sensitivity of each parameter

由图 3可以看出,各杆的弹性模量E、密度ρ和上弦杆横截面积A1灵敏度最高,而其他杆的横截面积灵敏度相对较小,所以选择EρA1作为待修正参数进行模型修正。

3.2 Kriging模型的构造及评估

基于拉丁超立方抽样抽出的200个样本构造Kriging模型。用前150个样本数据作为训练集,后50个数据作为测试集。以训练集作为Kriging模型的输入,以训练集所对应的加速度频响函数小波分解得到的第1层的小波系数作为输出,构造初始Kriging模型。以测试集Kriging模型的均方误差(MSE)均值平均值最小作为混合灰狼算法的个体目标函数,寻得最优相关系数θk。混合灰狼算法参数设置:种群规模Npop=40,最大迭代次数MaxIt=300,自变量维数nVar=1,缩放因子上下界beta_max和beta_min分别为0.2和0.8,交叉概率pCR=0.2,参数上下界ub和lb分别为0.01和100。最终寻得最优相关系数为0.875 8,迭代曲线如图 4所示。

图 4 混合灰狼算法迭代曲线 Fig. 4 Iteration curve of DE-GWO

图 4可以看出,当迭代次数增加至250次时,迭代曲线基本收敛,寻得最优值。相关系数确定以后,就可以构造Kriging模型。为评估所构建Kriging模型的精度,训练集样本第5个小波系数真实值和Kriging模型预测值的拟合曲线如图 5所示。为进一步评估Kriging模型的预测精度,测试集样本第10个小波系数真实值和Kriging模型预测值的拟合曲线如图 6所示,并根据式(19)计算出小波系数真实值和Kriging模型预测值间的均方根误差(RMSE),当RMSE越接近0时Kriging模型预测精度越高。

图 5 训练集样本第5个小波系数预测值 Fig. 5 Predicted value of the fifth wavelet coefficient of training set sample
图 6 测试集样本第10个小波系数预测值 Fig. 6 Predicted value of the 10th wavelet coefficient of test set sample
$ {\rm{RMSE}} = \frac{1}{{k\bar \lambda }}\sqrt {\sum\limits_{i = 1}^k {{{({\lambda _i} - {{\hat \lambda }_i})}^2}} } $ (19)
 

式中:k为每个样本点对应的小波系数数量;λ为小波系数真实值均值;λi为样本点对应的第i个小波系数;${{\hat \lambda }_i}$为Kriging模型对样本点对应的第i个小波系数的预测值。

图 5看出,通过训练集构建的Kriging模型预测值和真实值几乎重合,表明训练集的拟合精度很高。训练集的拟合精度并不能说明Kriging模型预测精度,必须用训练集以外的样本点来验证Kriging模型的预测精度,本文用测试集样本来验证。由图 6看出,除个别样本点外,测试集样本点对应的Kriging模型预测值几乎和真实值重合,并且通过式(19)计算出的RMSE值为8.164 1×10-4,表明所构建的Kriging模型预测精度很高,可以代替理论模型进行模型修正。

3.3 模型修正

当完成Kriging模型构造以后,就可以根据式(12)通过水循环算法以式(18)作为目标函数迭代求解待修正参数。假定试验参数值在有限元参数值的±20%区间内,使用水循环算法进行迭代寻优。水循环参数设置:雨滴层数为50,河流和大海个数总和为4,极小值为1×10-5,最大迭代次数为100。修正结果如表 3所示,水循环迭代曲线如图 7所示。

表 3 模型修正参数及其误差 Table 3 Parameters and their errors of model updating
修正参数 试验值 有限元值 初始误差/% 修正值 平均误差/%
E/GPa 190 209 10 190.011 856 6.24×10-3
ρ/(kg·m-3) 7 800 7 020 -10 7 802.893 8 3.71×10-2
A1/mm2 85.5 95 10 85.510 858 5 1.27×10-2
图 7 水循环迭代曲线 Fig. 7 Iteration curve of WCA

图 7可以看出水循环算法具有良好的搜索能力和快速收敛性能,算法在10次迭代以前就已收敛。多次运行寻优程序,迭代曲线基本不变化,也表明了该算法寻优的鲁棒性。由表 3可以看出,用本文所提方法进行模型修正取得了很好的修正效果。使用表 3中的各参数试验值和有限元值计算出试验模型和有限元模型对应的加速度频响函数。使用所得修正参数的平均修正值计算得到修正后的加速度频响函数,称其为Kriging模型加速度频响函数。Kriging模型、试验模型及有限元模型加速度频响函数曲线如图 8所示。

图 8 Kriging模型、试验模型及有限元模型频响函数曲线 Fig. 8 FRF curves of Kriging model, test model and finite element model

图 8可以看出修正过后的结构加速度频响函数曲线与试验加速度频响函数曲线几乎重合,与有限元模型加速度频响函数曲线形状相似,只是沿坐标轴有一定的移动。

为进一步检验本文修正方法的修正效果,比较频响函数修正前后实、虚部曲线,如图 9所示。由图 9可以看出修正后的有限元模型无论是实部还是虚部,都能很好地与试验模型的频响函数曲线重合。这进一步证明了本文所提方法的有效性。

图 9 修正前后频响函数曲线 Fig. 9 FRF curves before and after updating

实测试验时目标频响函数不可避免会受到随机噪声的干扰,在一些特殊环境下外界噪声干扰还很大。所以在计算的试验频响函数中加入信噪比为30、20、10和5 dB的高斯白噪声,来验证本文方法对于噪声的鲁棒性。图 10为加入信噪比为20 dB和5 dB时加速度频响函数图。

图 10 20 dB和5 dB信噪比下加速度频响函数 Fig. 10 FRF under 20 dB and 5 dB SNRs

将加入噪声的加速度频响函数按前文所提小波参数设置方法进行小波分解,得到不同信噪比下小波系数曲线和局部放大曲线,如图 11所示。由图 11可以看出,将加速度频响函数进行小波分解以后得到的序号较小的小波系数幅值较大,而绝大部分小波系数趋近于0。当在信号中加入随机噪声时,随着加入的高斯白噪声信噪比逐渐减小,对于幅值较大的小波系数影响较小,而对于幅值小的小波系数干扰较大。为了更直观地看出不同信噪比下噪声对于小波系数的影响,将小波系数随小波系数序号的变化分成3个部分,如各图中的3个放大曲线所示。放大曲线中左侧为加速度频响函数小波分解得到的第1层小波系数,个数较少,幅值较大;右上为加速度频响函数小波分解得到的中间靠前的小波系数,个数中等,幅值有大有小;右下为加速度频响函数小波分解得到的中间靠后的小波系数,个数较多,幅值趋近于0。

图 11 不同信噪比下小波系数曲线和局部放大图 Fig. 11 Wavelet coefficient curves and local enlarged curves under different SNRs

根据图 11(a)~图 11(d),由其中的左侧放大曲线可以看出第1层小波系数曲线随着加入的噪声加大,加噪小波系数曲线与未加噪小波系数曲线相似,大部分重合,只是当噪声很大时略有波动。由右上放大曲线可以看出,随着噪声的加大幅值大的小波系数曲线波动较小,而幅值小的小波系数曲线波动开始加大。由右下放大曲线可以看出,随着噪声的加大对于幅值趋近于0的小波系数影响较大,噪声越大小波系数曲线波动越大。综上所述,加速度频响函数经小波分解后得到较少数目幅值较大的小波系数,绝大部分小波系数趋近于0。当加速度频响函数受到噪声干扰时,对于幅值较大的小波系数干扰较小,而对于幅值较小的小波系数干扰较大。

在不同噪声水平下连续运行10次程序的模型修正结果均值如表 4所示。由表可以看出/有无噪声对于模型修正结果影响很大,而当在有噪声情况下,随着噪声的加大用本文所述方法仍然能够得出较满意的模型修正误差,即使当信噪比低至5 dB时,各参数模型修正误差仍然低于4%。由此证明了本文所提模型修正方法对于噪声具有较强的鲁棒性。

表 4 不同信噪比下修正结果 Table 4 Updating results under different SNRs
修正参数 初始误差/% 无噪声误差/% 平均误差/%
信噪比30 dB 信噪比20 dB 信噪比10 dB 信噪比5 dB
E/GPa 10 6.24×10-3 0.286 0.533 1.430 3.416
ρ/(kg·m-3) -10 3.71×10-2 0.485 0.453 1.718 3.718
A1/mm2 10 1.27×10-2 0.326 0.595 1.067 2.079
4 结 论

1) 将小波分解引入频响函数模型修正中,既保留了利用频响函数进行模型修正无需进行模态识别的优点,同时也避开了频响函数中频率点的选择难题。

2) 在构造Kriging模型时,用混合灰狼算法对Kriging模型相关系数进行寻优,使所构造的Kriging模型具有良好的拟合精度和预测能力,能代替有限元模型进行迭代计算,提高了模型修正效率。

3) 将加速度频响函数进行小波分解,用幅值较大的小波系数来表征频响函数进行模型修正,修正效果较好。算例表明,在不同噪声水平下,使用本文所提方法仍然可以得出较满意的修正效果,即使当信噪比低至5 dB时,各参数修正误差仍然低于4%,证明了本文所提方法对于强噪声的鲁棒性。

参考文献
[1] 张令弥. 计算仿真与模型确认及在结构环境与强度中的应用[J]. 强度与环境, 2002, 29(2): 42-47, 56.
ZHANG L M. Computer simulation & model validation with application to strength and environment engineering[J]. Structure & Environment Engineering, 2002, 29(2): 42-47, 56. (in Chinese)
Cited By in Cnki (51) | Click to display the text
[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.
Click to display the text
[3] 郭勤涛, 张令弥, 费庆国. 结构动力学有限元模型修正的发展-模型确认[J]. 力学进展, 2006, 36(1): 36-42.
GUO Q T, ZHANG L M, FEI Q G. From FE model updating to model validation:Advances in modeling of dynamic structures[J]. Advances in Mechanics, 2006, 36(1): 36-42. (in Chinese)
Cited By in Cnki (175) | Click to display the text
[4] 李伟明, 洪嘉振. 基于频响函数的模型修正方法[J]. 上海交通大学学报, 2011, 45(10): 1455-1459.
LI W M, HONG J Z. Research on model updating method based on frequency response functions[J]. Journal of Shanghai Jiao Tong University, 2011, 45(10): 1455-1459. (in Chinese)
Cited By in Cnki (30) | Click to display the text
[5] 朱凼凼, 冯咬齐, 向树红. 基于频率响应函数的动力学模型修正方法研究[J]. 中国工程科学, 2005, 7(8): 89-94.
ZHU D D, FENG Y Q, XIANG S H. Research of refinement methods of dynamic model based on frequency response functions[J]. Engineering Science, 2005, 7(8): 89-94. (in Chinese)
Cited By in Cnki (37) | Click to display the text
[6] HEMEZ F M, BROWN G W. Improving structural dynamics models by correlating simulated to measured frequency response functions[C]//39th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference and Exhibit. Reston: AIAA, 1998: 772-782.
Click to display the text
[7] 殷红, 马静静, 董小圆, 等. 基于频响函数的结构损伤识别模型修正方法[J]. 噪声与振动控制, 2019, 39(2): 7-13.
YIN H, MA J J, DONG X Y, et al. Structural damage identification using model updating method[J]. Noise and Vibration Control, 2019, 39(2): 7-13. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[8] ESFANDIARI A, SANAYEI M, BAKHTIARI-NEJAD F, et al. Finite element model updating using frequency response function of incomplete strain data[J]. AIAA Journal, 2010, 48(7): 1420-1433.
Click to display the text
[9] 杨海峰, 吴子燕, 吴丹. 基于加速度频率响应函数的结构损伤测量方法研究[J]. 振动与冲击, 2007, 26(2): 90-92, 97, 178.
YANG H F, WU Z Y, WU D. Structural damage detection method based on acceleration frequency response function[J]. Journal of Vibration and Shock, 2007, 26(2): 90-92, 97, 178. (in Chinese)
Cited By in Cnki (15) | Click to display the text
[10] KWON K S, LIN R M. Frequency selection method for FRF-based model updating[J]. Journal of Sound and vibration, 2004, 278(1-2): 285-306.
Click to display the text
[11] 张勇, 侯之超, 赵永玲. 基于频响函数截断奇异值响应面的有限元模型修正[J]. 振动工程学报, 2017, 30(3): 341-348.
ZHANG Y, HOU Z C, ZHAO Y L. Finite element model updating based on response surface of the truncated singular values of frequency response functions[J]. Journal of Vibration Engineering, 2017, 30(3): 341-348. (in Chinese)
Cited By in Cnki (8) | Click to display the text
[12] 张冬冬, 郭勤涛. Kriging响应面代理模型在有限元模型确认中的应用[J]. 振动与冲击, 2013, 32(9): 187-191.
ZHANG D D, GUO Q T. Application of Kriging response surface in finite element model validation[J]. Journal of Vibration and Shock, 2013, 32(9): 187-191. (in Chinese)
Cited By in Cnki (28) | Click to display the text
[13] 郭俊龙, 马立元, 李永军, 等. 大型复杂钢管焊接结构损伤识别新方法[J]. 振动与冲击, 2016, 35(18): 186-189, 195.
GUO J L, MA L Y, LI Y J, et al. New damage identification method for large complex steel tube welded structure[J]. Journal of Vibration and Shock, 2016, 35(18): 186-189, 195. (in Chinese)
Cited By in Cnki (4) | Click to display the text
[14] 韩忠华. Kriging模型及代理优化算法研究进展[J]. 航空学报, 2016, 37(11): 3197-3225.
HAN Z H. Kriging surrogate model and its application to design optimization:A review of recent progress[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(11): 3197-3225. (in Chinese)
Cited By in Cnki (162) | Click to display the text
[15] 王巨涛, 王春洁, 宋顺广. 结构加速度频响函数模型修正的Kriging方法[J]. 工程科学学报, 2017, 39(7): 1087-1093.
WANG J T, WANG C J, SONG S G. Structural model updating of acceleration frequency response function based on Kriging method[J]. Chinese Journal of Engineering, 2017, 39(7): 1087-1093. (in Chinese)
Cited By in Cnki (5) | Click to display the text
[16] WANG J T, WANG C J, ZHAO J P. Frequency response function-based model updating using Kriging Model[J]. Mechanical Systems & Signal Processing, 2017, 87(A): 218-228.
Click to display the text
[17] 杨修铭, 郭杏林, 李东升. 基于Kriging模型的频响函数有限元模型修正方法[J]. 计算力学学报, 2018, 35(4): 487-493.
YANG X M, GUO X L, LI D S. Kriging model based finite element model updating method using frequency response function[J]. Chinese Journal of Computational Mechanics, 2018, 35(4): 487-493. (in Chinese)
Cited By in Cnki (7) | Click to display the text
[18] DONOHO D L. De-noising by soft-thresholding[J]. IEEE Transactions on Information Theory, 1995, 41(3): 613-627.
Click to display the text
[19] 郭晓霞, 杨慧中. 基于粒子群优化算法的小波阈值去噪方法研究[J]. 计算机工程与应用, 2009, 45(31): 122-124.
GUO X X, YANG H Z. Research of wavelet threshold de-noising method based on particle swarm optimization[J]. Computer Engineering and Applications, 2009, 45(31): 122-124. (in Chinese)
Cited By in Cnki (17) | Click to display the text
[20] ZHU A J, XU C P, LI Z, et al. Hybridizing grey wolf optimization with differential evolution for global optimization and test scheduling for 3D stacked SoC[J]. Journal of Systems Engineering and Electronics, 2015, 26(2): 317-328.
Click to display the text
[21] ESKANDAR H, SADOLLAH A, BAHREININEJAD A, et al. Water cycle algorithm-a novel metaheuristic optimization method for solving constrained engineering optimization problems[J]. Computers and Structures, 2012, 110-111: 151-166.
Click to display the text
http://dx.doi.org/10.7527/S1000-6893.2020.23548
中国航空学会和北京航空航天大学主办。
0

文章信息

彭珍瑞, 曹明明, 刘满东
PENG Zhenrui, CAO Mingming, LIU Mandong
基于加速度频响函数小波分解的模型修正方法
Model updating method based on wavelet decomposition of acceleration frequency response function
航空学报, 2020, 41(7): 223548.
Acta Aeronautica et Astronautica Sinica, 2020, 41(7): 223548.
http://dx.doi.org/10.7527/S1000-6893.2020.23548

文章历史

收稿日期: 2019-10-08
退修日期: 2019-12-23
录用日期: 2020-01-19
网络出版时间: 2020-02-17

相关文章

工作空间