[发明专利]一种基于响应面建模和改进粒子群算法的有限元模型修正方法

说明书

技术领域

本发明涉及一种基于响应面建模和改进粒子群算法的有限元模型修正方法。

背景技术

现代航天器大多采用平台化、模块化设计，新型航天器往往采用对成熟平台的改进就能完成设计，尽量缩短研制周期，节约研究经费，以适应空间市场的快速发展。随着有限元计算软件的不断研发和完善，有限元计算结果已能够在一定程度上代替结构试验，得到具有一定精度的分析数据，从而对航天器的结构强度和在轨运行状态进行快速分析。

有限元模型修正技术在航空航天工程中应用广泛，高精度的有限元模型是准确进行结构力学特性分析的基础。由于建模误差等因素影响，有限元计算结果与实测结果之间总存在着难以消除的差异，利用结构实测响应修正有限元计算响应，使得修正后有限元模型计算响应值与试验测量值一致的过程即为有限元模型修正。有限元模型修正利用实际结构现场实测响应数据修正有限元理论模型计算数据，从而缩小理论模型与实际结构之间的误差，使得修正后模型计算的响应值与试验值趋于一致。运用模型修正技术可以省掉一些大型结构试验，根据需要仅作零部件级试验进行验证和修正，以使修正后的有限元模型能够代表真实的结构模型。修正过程充分利用结构试验和有限元分析两者优点，用少量的试验数据对有限元模型进行修正，修正后模型可以进行代替实际结构进行多种分析，节约试验经费和时间。

发明内容

本发明的目的是为了解决现有有限元模型修正方法效率低、精度有限的问题，本发明提供了一种基于响应面建模和改进粒子群算法的有限元模型修正方法。

本发明的基于响应面建模和改进粒子群算法的有限元模型修正方法是通过以下步骤实现的：

一、采用正交试验设计方法设计试验，将结构参数x＝[x₁，x₂，...x_n]按不同水平分为k组，则第k组的结构参数表示为k取整数；

二、选定用于构造响应面的结构参数x＝[x₁，x₂，...x_n]的初始迭代点，记为初始迭代点取k个组处结构参数的平均值；

三、用式(1)计算第j组处结构的响应值如下：

y^j=a0+Σi=1naixi+Σi=1nbiexp(-||xi-xi‾||22σ2)---(1)]]>

其中，x_i为相应第j组的第i个结构参数，为k个组中结构参数x_i的平均值；为第j组结构参数对应的响应面模型计算值；1≤j≤k，j取整数；σ表示响应面模型中的控制参数；