一种基于二次响应面反演的随机模型修正方法技术

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

【技术实现步骤摘要】

本专利技术涉及不确定性模型修正和系统参数识别领域，具体涉及一种基于二次响应 面反演的随机模型修正方法。
技术介绍
现行结构设计趋向精细化和大型化发展，以有限元分析为主的CAE分析技术在结 构产品设计中扮演越来越重要的地位。在结构力学性能分析、失效预警、损伤识别等方面， 准确而有效的有限元模型必不可少。然而由于结构本身的复杂性和建模过程中采取的不恰 当简化策略等多方面原因，根据设计图纸建立的有限元模型并不能反应真实结构的力学性 能。为了保证有限元模型和实验模型静动力响应一致性，可以利用模型修正技术进行有限 元模型参数识别。 近年来，国内外学者在模型修正理论研究方面取得很大进展，发展了很多模型修 正方法。但大多数模型修正方法都属于确定性模型修正方法，忽略了参数不确定性和测量 响应不确定性。在实际问题中，不确定性是普遍存在的，如传感器配置或信号处理方面的不 足；材料分散或设计制造误差；对复杂连接的不准确模拟等。确定性模型修正方法为寻找 参数准确名义值努力，忽略了客观存在的不确定性，这极大限制了其在实际中的应用，可靠 性大大降低。不确定性模型修正方法综合考虑系统存在的参数不确定性和响应测量不确定 性，以对不确定参数进行可靠范围估计为目标，应用更广，尤其适用于含明显不确定性的复 杂工程结构。 不确定性模型修正技术通过考虑实际结构存在的不确定性，并对不确定性量化分 析，扩展了模型修正的适用范围，提高了修正可靠性。对于高频冲击、高度非线性大变形和 随机耦合等存在明显不确定性的复杂结构，尤其适用。因此，开展相关理论研究并在实际结 构上检验理论方法的有效性具有重大意义。
技术实现思路
本专利技术针对工程系统参数识别或模型修正中存在的不确定性问题，提供一种基于 二次响应面反演的随机模型修正方法，该方法提高随机模型修正效率并保证了修正精度。 该方法可有效进行含明显不确定性复杂结构不确定参数识别，且识别效率和精度均较高。 本专利技术采用如下技术方案实现：， 其实施步骤如下： 步骤1 :根据结构设计图纸建立初步有限元模型，并确定结构中存在不确定性的 所有参数；结合实验模型，以结构响应输出（频率、位移等）构建优化目标函数，在不确定参 数初始区间下进行优化修正，得到不确定参数名义值和与实验模型匹配度最高的初步修正 模型； 步骤2 :以步骤1中得到的不确定参数名义值为中心，构造包络涉及不确定性的初 始超立方体；利用拉丁超立方抽样实验设计得到二次响应面构造所需原数据，并完成不确 定参数初始超立方下的系统输入输出二次响应面函数拟合； 步骤3 :根据步骤2中得到的二次响应面模型，利用最小二乘法建立基于二次响应 面反演的不确定参数修正迭代算法，结合实验模型拉丁超立方抽样响应结果，进行不确定 参数初始超立方修正； 步骤4 :有效性判别：如果步骤3中所得修正后不确定参数超立方小于初始超立方 体，则算法有效，转入下一步，否则，返回步骤2,重新确定初始超立方和二次响应面函数； 步骤5:计算迭代终止后不确定参数超立方中的二次响应面模型输出，并与实验 模型响应拉丁超立方抽样结果对比，如果两者范围一致，则输出修正后不确定参数超立方 体，否则，返回步骤2继续修正。 本专利技术的有益效果是在模型修正中综合考虑参数不确定性和响应测量不确定性， 并利用非不确定量分析方法完成随机模型修正，提出一种基于二次响应面反演的随机模型 修正方法。采用多项式响应面模拟系统输入输出关系并构建不确定参数反演迭代算法，不 仅回避了运算复杂且精度较低的不确定性分析方法（概率方法、模糊方法和凸模型方法 等），提高了运算效率，同时，无需构建灵敏度矩阵，简化问题分析且回避了病态灵敏度矩阵 问题，提高了修正精度。本专利技术可用于含明显不确定性的结构损伤识别、失效预警等方面， 为结构静动力响应提供可靠估计，具有重要理论和现实意义。【附图说明】 图1是本专利技术实现流程图； 图2是三自由度弹簧质量系统模型； 图3是三自由度弹簧质量系统一阶特征值二次响应面； 图4是三自由度弹簧质量系统二阶特征值二次响应面； 图5是三自由度弹簧质量系统三阶特征值二次响应面； 图6是三自由度弹簧质量系统修正前后eigl~eig2二维域分布； 图7是三自由度弹簧质量系统修正前后eigl~eig3二维域分布； 图8是三自由度弹簧质量系统修正前后eigl~|modalll|二维域分布； 图9是三自由度弹簧质量系统修正前后eig3~|modalll|二维域分布。【具体实施方式】 下面结合附图以及具体实施例进一步说明本专利技术。 如图1所示，本专利技术为，包括以下 步骤： 步骤1 :根据结构设计图纸建立初步有限元模型，并确定结构中存在不确定性的 所有参数。结合实验模型，以结构响应输出（频率、位移等）构建优化目标函数，在不确定 参数初始区间下进行优化修正，得到不确定参数名义值和与实验模型匹配度最高的初步修 正模型； 步骤2 :以步骤1中得到的不确定参数名义值为中心，构造包络涉及不确定性的初 始超立方体。利用拉丁超立方抽样实验设计得到二次响应面构造所需原数据，并完成不确 定参数初始超立方下的系统输入输出二次响应面函数拟合； 步骤3 :根据步骤2中得到的二次响应面模型，利用最小二乘法建立基于二次响应 面反演的不确定参数修正迭代算法，结合实验模型拉丁超立方抽样响应结果，进行不确定 参数初始超立方修正； 步骤4 :有效性判别：如果步骤3中所得修正后不确定参数超立方小于初始超立方 体，则算法有效，转入下一步，否则，返回步骤2,重新确定初始超立方和二次响应面函数； 步骤5:计算迭代终止后不确定参数超立方中的二次响应面模型输出，并与实验 模型响应拉丁超立方抽样结果对比，如果两者范围一致，则输出修正后不确定参数超立方 体，否则，返回步骤2继续修正。 在步骤1中，确定性模型修正包含以下步骤： 步骤1. 1 :根据结构设计图纸建立初步有限元模型，对结构中可能存在不确定性 的参数或连接，如材料弹性模量、构件几何尺寸、复杂连接阻尼系数等，利用不确定参数变 量标记并用区间变量描述其可能变化范围。 步骤1. 2 :以结构前n阶固有频和模态振型为系统输出，构建如下目标函数： 其中，wfl^、wshaPil>v别为第i阶频率和第i阶振型的第j自由度位移权重系数； fail(X)为第i阶频率理论计算值；为第i阶频率实验测量值；％jX)为第i阶振型的 第j自由度位移理论计算值；为第i阶振型的第j个自由度位移实验测量值。 步骤1. 3 :利用优化算法求取使得目标函数最小的不确定参数名义值。 在步骤2中，拟合二次响应面包含以下步骤： 步骤2. 1 :综合考虑结构可能涉及的不确定性范畴，以不确定性参数名义值为中 心，构建不确定参数初始超立方体。 步骤2. 2 :在超立方体范围内，利用拉丁超立方实验设计确定实验方案，并计算结 构响应输出，构建结构输入输出二次多项式响应面函数。 其中，为响应面系数，以^)为响应面输出，为输入参数，e为残差。 在步骤3中，基于二次响应面反演进行不确定参数修正包含以下步骤： 步骤3. 1 :首先将上式写为如下形式： 其中，e= (ei，e2, . ..，eK)T为残差向量；Y 为响应测量向量；C =(a:,a2, . ? .，aK)T为二次本文档来自技高网...

【技术保护点】
一种基于二次响应面反演的随机模型修正方法，其特征在于：包括以下步骤：步骤1：根据结构设计图纸建立初步有限元模型，并确定结构中存在不确定性的所有参数，结合实验模型，以结构响应输出构建优化目标函数，在不确定参数初始区间下进行优化修正，得到不确定参数名义值和与实验模型匹配度最高的初步修正模型；步骤2：以步骤1中得到的不确定参数名义值为中心，构造包络涉及不确定性的初始超立方体；利用拉丁超立方抽样实验设计得到二次响应面构造所需原数据，并完成不确定参数初始超立方下的系统输入输出二次响应面函数拟合；步骤3：根据步骤2中得到的二次响应面模型，利用最小二乘法建立基于二次响应面反演的不确定参数修正迭代算法，结合实验模型拉丁超立方抽样响应测量结果，进行不确定参数初始超立方修正；步骤4：有效性判别：如果步骤3中所得修正后不确定参数超立方小于初始超立方体，则算法有效，转入下一步，否则，返回步骤2，重新确定初始超立方和二次响应面函数；步骤5：计算迭代终止后不确定参数超立方中的二次响应面模型输出，并与实验模型响应拉丁超立方抽样结果对比，如果两者范围一致，则输出修正后不确定参数超立方体，否则，返回步骤2继续修正。

【技术特征摘要】

【专利技术属性】
技术研发人员：丁和武，邱志平，吕峥，朱静静，孙佳丽，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11