本发明涉及模型修正技术领域，尤其涉及一种考虑相关性的区间模型修正方法。
背景技术：
随着有限元法的广泛运用，基于实测或试验结果的有限元模型修正技术也得到了很大程度的发展。传统的模型修正方法都假设参数是确定性的，然而，工程实际中却普遍存在着不确定性，这些不确定性常常存在于材料的本构模型参数、结构几何参数、边界条件、初始条件、测量信息、认知判断等方面。若仍然采用确定性的方法和理论对有限元模型进行修正，必然导致修正结果的不可靠，与实际情况有较大出入。因此，考虑参数不确定性的随机模型修正方法越来越得到研究者们的关注。
结构不确定性的量化及传播分析方法有很多种，但应用较多的主要是其中的三种，分别是概率方法、模糊方法和区间分析方法。其中，概率方法是最为普遍使用的方法，即把结构参数视为随机变量来对问题进行建模和分析。但是，应用概率手段处理实际问题必须满足事件明确定义、大量样本存在、样本具有可重复性和不受人为因素影响等四个前提条件，在没有足量信息的情形下勉为其难地应用随机不确定理论及相应的可靠性分析方法，不确定参数概率分布的较小误差可能导致很大的分析偏差。而采用模糊方法来描述不确定性时，需要参数的隶属度函数。但在很多情况下，确定隶属度函数甚至比概率分布函数更为困难，即研究人员往往不得不带有很大的主观性来选取相应的隶属函数，使得分析结果的可靠性也值得怀疑。因此，为了反映客观实际、减少主观因素的干扰，在结构测试信息不够充分的情况下，将工程问题中的不确定性参数视为有界的区间变量，采用区间分析方法来修正有限元模型是十分必要的。
目前，基于区间分析的模型修正方法有很多。其中《mechanicalsystemsandsignalprocessing》2011年第25卷第4期中，公开了“intervalmodelupdatingwithirreducibleuncertaintyusingthekrigingpredictor”，基于全局优化理论，运用kriging模型替代有限元模型，采用确定性的模型修正技术，优化得到与每组实验数据相对应的结构参数，进而得到结构参数的区间。上述过程实际是多次的确定性模型修正过程，具有一定的局限性。《振动工程学报》2015年第28卷第1期中，公开了“不确定性参数识别的区间响应面模型修正方法”，该方法先将响应面表达式转化为完全平方项，再引入区间变量，以避免区间运算过程中的区间扩张问题。利用区间响应的上下界建立优化目标函数，并构建区间优化反问题，然后直接基于区间响应面模型实现修正过程。该方法避免了复杂的区间参数灵敏度计算，使得区间模型修正问题得到简化，并大幅提高修正效率。但是，该方法只能获取模型参数的区间，并不能求解模型参数之间的相关性矩阵。
由此可见，现有的模态更新方法大多数都是利用某些算法更新模态参数区间，在计算效率方面进行研究，而没有考虑过模型响应之间相关性对模型参数之间相关性的影响与传播规律。因此，亟需一种考虑相关性的区间模型修正方法，不仅能准确地修正模型参数的区间，而且能修正模型参数之间的相关性矩阵，从而提高模型的精度及可信度水平。
技术实现要素：
本发明的目的是提供一种考虑相关性的区间模型修正方法，既能准确地修正模型参数的区间，还能准确地描述响应之间相关性对模型参数之间相关性的影响与传播规律，提高模型的精度及可信度水平。
本发明的技术方案是提供一种考虑相关性的区间模型修正方法，包括如下步骤：
步骤1：建立结构或系统的仿真模型；
步骤2：多次测量结构或系统的多个响应，对响应信息的不确定性及相关性进行统计分析并建立椭球凸模型；
步骤3：建立区间模型修正的目标函数，并通过区间优化方法获取模型参数的修正区间；
步骤4：设定模型参数的初始相关性矩阵并利用椭球特征矩阵传播公式获取计算响应的椭球特征矩阵；
步骤5：将测量响应的椭球特征矩阵与计算响应的椭球特征矩阵作最小二乘优化，从而实现模型参数相关性的修正；
进一步地，所述步骤1还包括：
仿真模型包括有限元模型、响应面模型、kriging模型、神经网络模型。
进一步地，所述步骤2还包括：
测量响应
式中，
式中，λii＝1,2,…,n、λjj＝1,2,…,n是n维椭球的第i维和第j维的半主轴长，θ称为椭圆的姿态角，并满足-45°≤θ≤45°，
进一步地，所述步骤3还包括：
建立区间模型修正的目标函数如下：
式中，
区间优化方法即采用单目标优化方式，通过建立线性约束来设置参数变化的界限，并利用二次规划算法求得目标函数的最优解。
进一步地，所述步骤4还包括：
椭球特征矩阵的传播公式如下：
式中，ωy表示计算响应的椭球特征矩阵，fi表示第i个系统gi(·)对模型参数x＝(x1,x2,…,xn)的一阶导系数列矩阵，covx表示模型参数的n维协方差方阵，cij表示矩阵的迹，其计算公式如下：
式中，tr表示求迹运算符，fi和fj分别表示第i个系统gi(·)和第j个系统gj(·)对模型参数x＝(x1,x2,…,xn)的二阶导系数方阵。
进一步地，所述步骤5还包括：
将测量响应的椭球特征矩阵与计算响应的椭球特征矩阵作最小二乘优化，其过程如下：
本发明创新性地提出一种考虑相关性的区间模型修正方法。该方法不仅能准确地修正模型参数的不确定性区间，而且能修正模型参数之间的相关性。
根据上述技术方案，本发明的有益效果包括：
(1)本发明提出一种考虑相关性的区间模型修正方法，该方法能够准确高效地修正模型参数的区间。
(2)本发明考虑了测量响应之间的相关性矩阵，通过椭球特征矩阵传播公式能准确地修正模型参数之间的相关性矩阵。
(3)本发明通过结合模型参数的区间与相关性矩阵，运用椭球凸模型的建模方法，提高了模型的精度及可信度水平。
附图说明
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
图1是考虑相关性的模型修正流程图；
图2是齿轮装配体的有限元模型；
图3是修正模型参数的椭球凸模型；
图4是第一阶和第二阶模态频率的椭球凸模型。
具体实施方式
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理。
如图2所示，本实施例提供了一种考虑相关性的区间模型修正方法，具体包括以下步骤：
步骤1：建立结构或系统的仿真模型；
仿真模型包括有限元模型、响应面模型、kriging模型、神经网络模型。
该实施例中，根据齿轮装配体建立如图2所示的有限元模型，其中齿轮装配体中大小齿轮的弹性模量e和泊松比μ未知，且e和μ为的区间变量。由于本算例为仿真算例，故假定e、μ的真实区间分别为[2.0350,2.0901]×1011pa、[0.2507,0.3392]，相关系数ρ＝0.3671，并在对应真实区间内抽取100个样本。
步骤2：多次测量结构或系统的多个响应，对响应信息的不确定性及相关性进行统计分析并建立椭球凸模型；
该步骤中，测量响应
式中，
式中，λii＝1,2,…,n、λjj＝1,2,…,n是n维椭球的第i维和第j维的半主轴长，θ称为椭圆的姿态角，并满足-45°≤θ≤45°，
将步骤1抽取的100个待修正模型参数e、μ样本通过有限元模型计算得到100个第一、二阶自由模态频率，选取其作为实验测量响应，并对响应信息的不确定性及相关性进行统计分析，根据椭球凸模型公式，建立齿轮装配体的第一、二阶自由模态频率的椭球凸模型；
步骤3：建立区间模型修正的目标函数，并通过区间优化方法获取模型参数的修正区间；
该步骤中，建立如下区间模型修正的目标函数，并设定待修正模型参数e和μ的初始搜索区间分别为[1.8000,2.3000]×1011pa、[0.2000,0.4000]：
式中，
区间优化方法即采用单目标优化方式，通过建立线性约束来设置参数变化的界限，并利用二次规划算法求得目标函数的最优解，所求模型参数的修正区间如表1所示。
步骤4：设定模型参数的初始相关性矩阵并利用椭球的特征矩阵传播公式获取计算响应的椭球特征矩阵；
设定待模型参数e、μ的相关性矩阵为
式中，ωy表示计算响应的椭球特征矩阵，fi表示第i个系统gi(·)对模型参数x＝(x1,x2,…,xn)的一阶导系数列矩阵，covx表示模型参数的n维协方差方阵，cij表示矩阵的迹，其计算公式如下：
式中，tr表示求迹运算符，fi和fj分别表示第i个系统gi(·)和第j个系统gj(·)对模型参数x＝(x1,x2,…,xn)的二阶导系数方阵。
步骤5：将测量响应的相关性矩阵与计算响应的相关性矩阵作最小二乘优化，从而实现模型参数相关性的修正；
将测量响应的椭球特征矩阵与计算响应的椭球特征矩阵作最小二乘优化，其过程如下：
直到模型参数的相关信息满足精度，输出的结果如表1所示。
表1齿轮装配体系统参数区间估计值
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，未公开的内容均是本领域的公知常识或惯用技术手段，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。