基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法技术

本发明专利技术涉及一种基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法，属于模型修正方法技术领域。本发明专利技术选取损伤特征参数作为模型修正的输入参数，以对应损伤状态下的物理参数作为模型修正的输出参数，损伤特征参数运用了频率和振型变化进行组合，选择待修正的物理参数时综合考虑可能对结构模态信息及损伤识别产生影响的结构边界条件及材料性能参数，本发明专利技术应用支持向量机统计学习的算法，并采用分布修正的策略，不仅减少了一次修正样本点数，使训练样本大大减少，而且提高了模型修正的准确性,使修正后的有限元模型在火灾过程中及火灾后能够更加真实的反应结构振动特性及其他结构响应。

【技术实现步骤摘要】
基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法
本专利技术创造涉及一种模型修正方法，具体地说，涉及一种基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法。
技术介绍
结构健康监测早已成为结构防震减灾领域的研究热点，但对于高温下结构动力特性研究因受试验条件及有限元模型建立等诸多困难的困扰，在国内外还十分罕见。桥梁方面的模型修正研究较多，但针对钢筋混凝土简支梁的模型修正鲜有研究，大部分基于动力测试的模型修正方法较繁琐且未能考虑对损伤更敏感的振型的影响，传统方法极易导致“数据爆炸”及映射能力降低等现象，并且仅以材料参数作为待修正参数，对边界条件考虑不足，无法使有限元模型全面、正确地反映结构的真实性。针对有限元模型修正算法，部分学者采用神经网络方法作为回归机制，进行有限元模型，但从理论上来说，神经网络利用经验风险最小化选取模型，但经验风险最小化替代期望风险小化选取模型是存在错误的。
技术实现思路
本专利技术的目的在于克服现有模型修正方法存在的上述缺陷，提出了一种基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法，采用分步修正的方式，不仅减少了一次修正样本点数，使训练样本大大减少，而且提高了模型修正的准确性,使修正后的有限元模型在火灾过程中及火灾后能够更加真实的反应结构振动特性及其他结构响应。本专利技术具体通过以下步骤实现技术方案：步骤一、建立ANSYS精细化模型；步骤二、选定模型中不确定物理参数作为待修正参数；为充分考虑可能对结构模态信息及损伤识别产生影响的结构边界条件及材料性能参数，分别选取简支梁支座偏移(D1、D2)、支座刚度(K1、K2)、混凝土弹性模量(E)、混凝土密度(ρ)等可能存在误差的物理参数作为有限元模型修正的待修正物理参数参数库。火灾后的待修正物理参数包括混凝土弹性模量及混凝土密度。步骤三、确定待修正物理参数基准值及取值范围；对于物理参数的取值范围，参考经验值，并基于实测模态信息进行预估。步骤四、进行参数灵敏度分析，确定模型的修正物理参数；步骤五、依据均匀设计表构造样本点，并利用ANSYS精细化模型计算样本点对应的结构响应，依据结构响应构造损伤特征参数；火灾前的损伤特征参数为频率及振型的变化组合：{VCIP}＝{FCR1，FCR2…，FCRm，DF1，DF2…DFn}其中：式中：FCRi-第i阶频率的变化率；fui-结构损伤前第i阶频率；fdi-结构损伤后第i阶频率；n-使用的振型阶数；m-使用的频率阶数；DFi-i阶模态对应于p个测试自由度的归一化后的振型向量；Φij-第i阶模态第j测试自由度的振型分量。火灾后的损伤特征参数为频率及振型的变化组合叠加位移信息：{VCIP+}＝{FCR1，FCR2...，FCRm，DF1，DF2...DFn，L1，L2...，LZ}式中Z为拾取的位移数量，Lz为对应第Z个位移值。步骤六、将样本点对应的物理参数值及相应的损伤特征参数输入到修正界面中，建立对应的支持向量机修正模型；步骤七、将实测模态信息进行组合，带入到建立的模型中，输出对应的物理参数回归值；本专利技术对传统的支持向量机GUI操作界面进行改进，改进后界面能够独立完成组合参数的计算，能够显示平方相关系数随核函数的改变，并且能够考虑多种核函数，本专利技术选用高斯径向基核函数(RBF)进行回归；支持向量机回归机参数的选取对于回归精度也至关重要，选取合适的惩罚参数c和函数的参数g的取值范围，在此取值范围内，运用交叉验证方法，综合考虑平均平方误差(MSE)和平方相关系数(r2)选择合理的输入参数。具体的支持向量机回归流程如图1所示，指标表达式如下：平均平方误差(MSE)平方相关系数(r2)步骤八、依据修正结果重新设置模型物理参数，获得模型结构响应，并与实测对比；步骤九、评价模型修正收敛性，若不收敛，重复步骤一到八；火灾前模型修正的收敛准则为同时满足以下两个条件：结构响应修正终止：式中为第i阶n次修正后结构频率，δ为理想收敛界限值；物理参数收敛终止：式中Yi、Xi为第i个修正物理参数的初始修正区间上下限，a’ij为第i个物理参数j次修正后结果，δ及ε依据实际工程经验进行取值。火灾后模型修正的收敛准则为：其中J1、J2分别为绝对收敛终止指标和相对收敛终止指标，J1、J2上限依据工程经验进行取值，当两个收敛指标有一个满足时，就终止修正，分别为修正前结构响应及n次修正后结构响应；m表示选取的结构响应总数。步骤十、对达到收敛指标的结构模型，进行修正评价。模型修正的准确性的评价准则为：a、频率相关性评价指标令结构第i阶频率的预测值和实测值分别为那么他们之间的相关性表达式可以如下表示：当Ew越趋近于0时，表示模型修正结果越好。b、模态置信准则评价模态置信准侧同时考察了结构固有频率和振型的相关性，对结构修正前后振型的相关性做出了较好的评价，计算公式如下：上式中，表示有限元计算的第i阶各测点模态振型，[ψ]i表示实验测试中获取的第i阶各测点模态振型，当Ew趋近于零，MAC趋近于1时，表示有限元模型和实际模型的模态相关性越好。通常，这种趋势达到一定范围内时，我们就认为模型计算值与实验计算值误差在允许范围内，不再对其进一步修正。本专利技术的有益效果是：(1)该专利技术应用支持向量机统计学习的算法，并采用分布修正的策略，避免了传统方法“数据爆炸”及映射能力降低等现象，减少了一次修正的样本点数，使训练样本大大减少，解决多物理参数、多样本修正问题，提高了模型修正的准确性；(2)该专利技术应用支持向量机统计学习的算法构建模型修正网络，在进行模型修正时综合考虑固有频率与振型二者的相关性，克服单独运用其中一种参数的不足，是一种更全面、更精确的模型修正方法；(3)该专利技术利用提出的支持向量机分步有限元模型修正方法将边界条件纳入模型修正对象之中，使有限元模型全面、正确地反映结构的真实性；(4)该专利技术利用支持向量机统计学习的算法,可最终转化为凸优化问题,因而可保证算法的全局最优性,避免了神经网络无法解决的局部最小问题，另外支持向量机的最终决策函数只由少数的支持向量所确定，计算的复杂性取决于支持向量的数目，而不是样本空间的维数，这在某种意义上避免了高维变量引起的灾害；(5)该专利技术对传统的支持向量机GUI操作截面进行改进，改进后界面能够独立完成组合参数的计算，能够显示平方相关系数随核函数的改变，并且能够考虑多种核函数。附图说明图1是本专利技术所述支持向量机的回归流程。图2是本专利技术所述修正方法的流程图。图3是本专利技术实施方式所述混凝土简支梁的配筋详图。图4是本专利技术实施方式所述结构响应灵敏度分析结果。具体实施方式本实施方式设计制作了4根尺寸相同的钢筋混凝土简支梁，分别用于研究简支梁在60min、90min、120min及150min火灾作用下结构振动特性。考虑到实验室条件限制，简支梁设计总长度3000mm，有效长度2700mm(两端各留有150mm支撑长度)，截面尺寸为250mm×400mm，试件采用C35商品混凝土浇筑，简支梁保护层厚度取30mm。配筋详图如图3所示，试验工况表如表1所示。表1试验工况表为充分考虑可能对结构模态信息及损伤识别产生影响的结构边界条件及材料性能参数，分别选取简支梁支座偏移(D1、D2)、支座刚度(K1、K2)、混凝土弹性模量(E)、混凝土密度(ρ)作本文档来自技高网...

【技术保护点】
一种基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法，其特征在于：包括如下步骤：步骤一、建立ANSYS精细化模型；步骤二、选定模型中不确定物理参数作为待修正参数；步骤三、确定待修正物理参数基准值及取值范围；步骤四、进行参数灵敏度分析，确定模型的修正物理参数；步骤五、依据均匀设计表构造样本点，并利用ANSYS精细化模型计算样本点对应的结构响应，依据结构响应构造损伤特征参数；步骤六、将样本点对应的物理参数值及相应的损伤特征参数输入到修正界面中，建立对应的支持向量机修正模型；步骤七、将实测模态信息进行组合，带入到建立的模型中，输出对应的物理参数回归值；步骤八、依据修正结果重新设置模型物理参数，获得模型结构响应，并与实测对比；步骤九、评价模型修正收敛性，若不收敛，重复步骤一到八；步骤十、对达到收敛指标的结构模型，进行修正评价。

【技术特征摘要】
1.一种基于动力或静动力的钢筋混凝土简支梁火灾模型修正方法，其特征在于：包括如下步骤：步骤一、建立ANSYS精细化模型；步骤二、选定模型中不确定物理参数作为待修正参数；步骤三、确定待修正物理参数基准值及取值范围；步骤四、进行参数灵敏度分析，确定模型的修正物理参数；步骤五、依据均匀设计表构造样本点，并利用ANSYS精细化模型计算样本点对应的结构响应，依据结构响应构造损伤特征参数；步骤六、将样本点对应的物理参数值及相应的损伤特征参数输入到修正界面中，建立对应的支持向量机修正模型；步骤七、将实测模态信息进行组合，带入到建立的模型中，输出对应的物理参数回归值；步骤八、依据修正结果重新设置模型物理参数，获得模型结构响应，并与实测对比；步骤九、评价模型修正收敛性，若不收敛，重复步骤一到八；步骤十、对达到收敛指标的结构模型，进行修正评价。2.根据权利要求1所述的模型修正方法，其特征在于：步骤二所述的火灾前的待修正物理参数为包括简支梁支座偏移、支座刚度、混凝土弹性模量、混凝土密度。3.根据权利要求1所述的模型修正方法，其特征在于：步骤二所述的火灾后的待修正物理参数包括混凝土弹性模量及混凝土密度。4.根据权利要求1所述的模型修正方法，其特征在于：步骤五所述的火灾前的损伤特征参数为频率及振型的变化组合：{VCIP}＝{FCR1，FCR2...，FCRm，DF1，DF2...DFn}其中：式中：FCRi—第i阶频率的变化率；fui—结构损伤前第i阶频率；fdi—结构损伤后第i阶频率；n—使用的振型阶数；m—使用的频率阶数；DFi—i阶模态对应于p个测试自由度的归一化后的振型向量；Φij—第i阶模态第j测试自由度的振型分...

【专利技术属性】
技术研发人员：刘才玮，巴光忠，苗吉军，
申请(专利权)人：青岛理工大学，
类型：发明
国别省市：山东,37