[发明专利]一种基于二次响应面反演的随机模型修正方法

摘要

本发明提供一种基于二次响应面反演的随机模型修正方法，该方法针对工程系统参数识别或模型修正中存在的不确定性问题，通过单点逐次反演迭代，有效回避了常规随机模型修正方法中，因不确定性分析和灵敏度矩阵求解带来的误差。通过实例验证和方法对比说明，本方法可有效进行含明显不确定性结构的不确定参数识别，且识别效率和精度均很高。

主权项

1.一种三自由度弹簧质量系统基于二次响应面反演的随机模型修正方法，其特征在于：包括以下步骤：步骤1：根据结构设计图纸建立初步有限元模型，并确定结构中存在不确定性的所有参数，结合实验模型，以结构响应输出构建优化目标函数，在不确定参数初始区间下进行优化修正，得到不确定参数名义值和与实验模型匹配度最高的初步修正模型；步骤2：以步骤1中得到的不确定参数名义值为中心，构造包络涉及不确定性的初始超立方体；利用拉丁超立方抽样实验设计得到二次响应面构造所需原数据，并完成不确定参数初始超立方下的系统输入输出的二次响应面函数拟合；步骤3：根据步骤2中得到的二次响应面函数，利用最小二乘法建立基于二次响应面反演的不确定参数修正迭代算法，结合实验模型拉丁超立方抽样响应测量结果，进行不确定参数初始超立方修正；步骤4：有效性判别：如果步骤3中所得修正后不确定参数超立方小于初始超立方体，则算法有效，转入下一步，否则，返回步骤2，重新确定初始超立方和二次响应面函数；步骤5：计算迭代终止后不确定参数超立方中的二次响应面模型输出，并与实验模型响应拉丁超立方抽样结果对比，如果两者一致，则输出修正后不确定参数超立方体，否则，返回步骤2继续修正；在步骤2中，包含以下步骤：步骤2.1：综合考虑结构可能涉及的不确定性范畴，以不确定性参数名义值为中心，构建不确定参数初始超立方体；步骤2.2：在超立方体范围内，利用拉丁超立方实验设计确定实验方案，并计算结构响应输出，构建结构输入输出二次多项式响应面函数：其中，a,b_i,d_ij为响应面系数，为响应面输出，X_i,X_j为输入参数，ε为残差；在步骤3中，包含以下步骤：步骤3.1：首先将步骤2.2中的公式写为如下形式：其矩阵形式为：其中，ε＝(ε₁,ε₂,...,ε_K)^T为残差向量；为响应测量向量；C＝(a₁,a₂,...,a_K)^T为二次响应面常数向量；B＝[b_ki](k＝1,...,K；i＝1,...,n)为二次响应面一次项系数矩阵；为二次响应面二次项系数矩阵；为待修正参数二次齐次矩阵；X＝(X₁,X₂,...,X_n)^T为待修正系数向量；通过系统输出修正不确定参数区间转化为如下优化问题：根据最小二乘法，▽(εTε)＝0其中，根据上式构造基于二次响应面反演的不确定参数修正迭代公式：其中，W＝λI为加权矩阵；步骤3.2：将实验模型响应拉丁超立方抽样结果依次代入上式，进行不确定参数反演修正；系统参数取值为：k2,k4,k5为不确定参数，且处于如下区间：k2＝[7.5 8.5]N/m,k4＝[1.8 2.2]N/m,k5＝[1.8 2.2]N/m不确定参数初始区间假设如下：k2＝[6.5 9.5]N/m,k4＝[1.6 2.4]N/m,k5＝[1.6 2.4]N/m系统输出定义为结构前三阶特征值，即：z＝[eig1 eig2 eig3]T其中，z为系统输出，eig1为第一阶特征值，eig2为第二阶特征值，eig3为第三阶特征值；1)系统输出二次响应面拟合：在参数初始区间内，二次响应面待定形式如下：根据最小二乘法，得到响应面计算结果如下：一阶特征值：二阶特征值：三阶特征值：2)不确定反向传播模型修正：基于二次响应面推导不确定参数修正迭代公式和拉丁超立方抽样测量响应，修正三自由度弹簧质量系统不确定参数；在步骤1中，包含以下步骤：步骤1.1：根据结构设计图纸建立初步有限元模型，对结构中可能存在不确定性的参数或连接，利用不确定参数变量标记并用区间变量描述其可能变化范围；步骤1.2：以结构前n阶固有频和模态振型为系统输出，构建如下目标函数：其中，w_fre,i、w_shap,ij分别为第i阶频率和第i阶振型的第j自由度位移权重系数；f_a,i(X)为第i阶频率理论计算值；f_m,i为第i阶频率实验测量值；为第i阶振型的第j自由度位移理论计算值；为第i阶振型的第j个自由度位移实验测量值；步骤1.3：利用优化算法求取使得目标函数最小的不确定参数名义值；在步骤4中，操作如下：将修正后不确定参数超立方体与初始超立方体比较，如果修正后超立方体大于初始超立方体，则退回步骤2，重新构造不确定参数初始超立方体及二次响应面；如果修正后参数超立方体小于不确定参数初始超立方体，则进行下一步；在步骤5中，操作如下：判断迭代收敛所得参数超立方体范围内，二次响应面输出范围与实验模型拉丁超立方抽样是否一致；如果两者一致，则输出最终修正所得不确定参数超立方体，即为修正结果；如果两者响应分布超出误差容限，则退回步骤2重新进行修正；三自由度弹簧质量系统参数取值为：k2,k4,k5为不确定参数，且处于如下区间：k2＝[7.5 8.5]N/m,k4＝[1.8 2.2]N/m,k5＝[1.8 2.2]N/m不确定参数初始区间假设如下：k2＝[6.5 9.5]N/m,k4＝[1.6 2.4]N/m,k5＝[1.6 2.4]N/m系统输出定义为结构前三阶特征值，即：z＝[eig1 eig2 eig3]T其中，z为系统输出，eig1为第一阶特征值，eig2为第二阶特征值，eig3为第三阶特征值；1)系统输出二次响应面拟合：在参数初始区间内，二次响应面待定形式如下：根据最小二乘法，得到响应面计算结果如下：一阶特征值：二阶特征值：三阶特征值：2)不确定反向传播模型修正：基于二次响应面推导不确定参数修正迭代公式和拉丁超立方抽样测量响应，修正三自由度弹簧质量系统不确定参数，修正结果列于表1，表1三自由度弹簧质量系统不同方法修正结果与基于Kriging模型和径向基模型修正结果对比，本方法所得结果区间与真实区间吻合度最高。