基于多响应频响函数的加筋壁板结构模型修正

2020-06-17展铭郭勤涛岳林张保强

中南大学学报（自然科学版） 2020年5期

展铭，郭勤涛，岳林，张保强

(1.南京航空航天大学机电学院，江苏南京，210016；2.厦门大学航空航天学院，福建厦门，361005)

在土木、机械、航空航天、武器装备等领域的工程结构中，有限元建模和分析技术得到越来越多的应用，并逐渐成为理论分析和试验测试后的另一大支柱[1-2]。虽然有限元技术在工程实践中十分普及，但是随着工程结构的复杂化，在有限元建模过程中，局部细节简化和边界条件处理等不可避免地会使分析结果与结构真实测试结果之间存在差异。将有限元分析结果用于结构设计或工程决策的前提是有限元模型具有足够的精度，即所建立的模型能够准确反映结构的输入和输出之间的关系。因此，在诸如上述的工程领域中，根据结构测试的试验数据对有限元模型进行修正，可以提高有限元模型反映结构特性的能力，并可以将修正后的有限元模型用于结构优化设计、结构损伤识别、状态监测和结构响应预示等[3]。大量学者对有限元模型修正的关键问题进行了研究，并取得了丰硕的研究成果。在模型修正前，首先需要对模型中的参数进行筛选，ZANG等[4]采用等效单元模态应变能和等效单元模态动能选择模型修正中误差的敏感参数。万华平等[5]基于全局灵敏度分析方法，提出了有限元模型修正待修正参数选择的新方法。经过实践验证，待修正响应信息对有限元模型修正过程、修正后模型的精度以及修正后模型的用途具有重要的影响。REN等[6]基于结构静态试验测试结果修正了结构中的未知参数。THONON等[7]采用模态频率和模态置信度作为修正目标，实现了梁结构的模型修正。XIAO等[8]采用模态频率和影响线作为目标修正了斜拉桥的多尺度有限元模型。张保强等[9]使用模态频率和有效模态质量作为结构动力学模型修正的目标，实现了复杂边界梁模型的修正。ESFANDIARI[10]利用共振区附近的应变传递函数修正了二维桁架结构。JAISHI等[11]采用模态频率及模态应变能残差，基于多目标优化算法修正了简支梁及箱型梁结构的有限元模型。GUO等[12]通过定义应变频响函数的相关性，修正了结构的材料属性。姜东等[13]基于模态试验参数修正了蜂窝夹层复合材料有限元模型中的不确定性参数。郭铁能等[14]研究了应用频响函数辨识结合部参数的方法，并基于三自由度弹簧质量系统验证了方法的有效性。LINK等[15]将模态参数和时域响应作为目标响应，基于模型修正方法对复合板结构中的损伤进行了识别。针对复杂结构模型修正过程计算成本过高的问题，曹文斌等[16]以某型航天结构为例，基于代理模型以及二级优化策略实现了有限元模型修正过程。JUNG等[17]通过分层模型校准方法，提高了电机有限元模型的预测能力。XIONG等[18]对有限元模型修正方法进行了总结，通过最大似然估计确定模型中随机参数的分布特征。MODAK等[19]对模型修正方法进行了比较，认为不同类型的模型修正方法各具优势，在工程应用时需要酌情选取。从上述的文献可知，目前常用的有限元模型修正方法中，采用的响应主要包括静力学响应、模态参数、加速度时频域响应等。静力学响应、结构整体模态参数等对于结构中局部连接参数的灵敏度较低。此外，采用模态参数作为修正目标，要求试验能够精确测量多阶模态，并且对阻尼参数的修正效果通常不理想。本文作者采用复杂加筋壁板结构为研究对象，联合加速度和应变频响函数对模型中的材料参数、连接参数以及模态阻尼比进行修正，并对修正后的模型进行检验。

1 基于频响函数的模型修正

1.1 模型修正方法

有限元模型修正的基本思想是通过调整模型中的参数达到缩小有限元分析与结构试验测试之间的差距的目的，本质上是一个力学反问题，并可以表示为如下数学形式：

式中：X为待修正参数；Xlb和Xub分别为其下限和上限；R(X)为响应残差；RA(X)和RE分别为有限元分析响应以及试验测试响应；WR为残差的加权系数。

1.2 多响应频响函数分析方法

在有限元分析时，通常将结构离散为多自由度系统，其动力学方程为

式中：M，C和K分别为质量、阻尼和刚度矩阵；X(t)，˙(t)和¨(t)分别为位移、速度及加速度向量；F(t)为激励力向量。

对上述动力学方程进行变换，引入模态坐标对方程组进行解耦，根据模态叠加原理可得在任意坐标下结构的位移响应为

式中：φr为第r阶模态振型向量；N为模态截断的总阶次；Mr，Cr和Kr分别为模态质量、模态阻尼及模态刚度；为第r阶模态频率；ζr=为第r阶模态阻尼比；Hd(ω)为位移的频率响应函数。

根据有限元理论中的应变和位移之间的关系以及上述位移响应表达式，通过坐标变换等过程可得结构中任一节点处的应变响应为

式中：ψr称为第r阶应变模态振型；Hε(ω)为应变频响函数。

1.3 频响函数的相关性

假设有限元分析获得的结构关键点处的频响函数矩阵为HA，试验的频响函数矩阵为HE。借鉴评价模态振型相关性的概念，定义频响函数的相关系数描述仿真分析与试验测试结果之间相关性[20]：

式中：HiA(X)和HiE分别为仿真分析和试验测试频响函数矩阵在频率点i处的列向量；RSAC和RMAC分别为形状相关系数和幅值相关系数。

由两类相关系数的定义式可以看出，RSAC和RMAC分别是对频响函数矩阵中各频率点处的频响函数向量的形状相似程度和幅值一致程度的评价指标。基于频响函数相关系数的有限元模型修正目标函数可表示为

2 加筋壁板结构建模及试验

2.1 加筋壁板结构建模

加筋壁板结构在航空航天等领域广泛使用，其动力学特性决定着结构在各个方面的研究进展。加筋壁板主要由基础板、夹持边以及筋条等零部件组成，各零部件之间通过螺栓或铆钉相连接。在有限元建模时，考虑到各零部件厚度相对于其余方向尺寸较小，采用板单元建立加筋壁板各零部件的模型。对于零部件之间的连接，采用弹簧阻尼单元描述。由于各零部件之间的连接方向不一，因此在建模时对弹簧阻尼单元按照连接方向进行分组，并分别赋予不同的属性值。将螺栓视为弹性杆件，根据材料力学对其轴向刚度进行近似计算，并将计算结果圆整后作为模型中参数的初值。此外，由于结构中的连接数量较多，为了减小其质量对结构响应计算的影响，在建模时采用集中质量单元描述螺栓的质量。加筋壁板结构的有限元模型如图1所示。

2.2 加筋壁板试验测试

图1 加筋壁板结构有限元模型Fig.1 Finite elementmodelof stiffened structure

在对加筋壁板进行试验前，首先对各零部件进行组装，并通过扭矩扳手控制螺栓的拧紧力矩，保证同一方向的螺栓具有相同的拧紧力。试验时采用弹性绳将加筋壁板悬挂于刚性支架上，通过激振器产生0～200 Hz的连续扫频信号对结构基础板施加垂向激励。采用动态信号采集仪，同时记录激励力信号以及基础板上测点处的加速度和应变信号。随后基于信号处理方法获得各测点处加速度及应变频响函数信息。加筋壁板试验测试激励点及响应测点分布如图2所示，其中，C1为激励点，激励方向垂直于基础板；A1～A5为加速度测点，方向垂直于基础板；S1～S6为应变测点，其中S1～S2以及S5～S6测试方向沿Y向，S3～S4测试方向沿X向。

图2 加筋壁板结构试验测点Fig.2 Testpointof stiffened structure

表1 加筋壁板模型修正参数Table1 Updating parametersof stiffened structure

3 加筋壁板结构模型修正

正如建模部分所述，加筋壁板模型中存在多组未知参数，本文采用多响应频响函数对模型中参数进行修正，一方面，可以避免模态参数识别引起的误差，另一方面，频响函数上可用于修正的信息更多。在对加筋壁板结构模型修正前，首先采用灵敏度分析方法对模型中参数进行筛选，确定结构基础板与夹持边连接的Z向平动刚度Kn1以及X和Y向转动刚度Krx1和Kry1、基础板与筋条连接的Z向平动刚度Kn2以及Y向转动刚度Kry2、结构前6阶模态的阻尼比ξ1～ξ6为待修正参数。在模型修正时，选取频响函数各共振峰的半功率带宽点以及相邻共振峰的峰值间距的中点作为目标频率点，分别计算加速度和应变频响函数的两类相关系数，基于式(7)构造用于修正的目标函数。采用遗传算法对所构造的目标函数进行优化，从而获得待修正参数合理空间内的最优组合。加筋壁板模型修正各待修正的初值、修正值如表1所示。分别基于修正前后的模型参数计算加筋壁板结构的模态频率，并与试验测试结果比较，其各阶频率的误差如表2所示。