利用改进交叉模型交叉模态的随机模型修正方法_王炎.pdf
《利用改进交叉模型交叉模态的随机模型修正方法_王炎.pdf》由会员分享,可在线阅读,更多相关《利用改进交叉模型交叉模态的随机模型修正方法_王炎.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 36 卷第 2 期2023 年 4 月振 动 工 程 学 报Journal of Vibration EngineeringVol.36 No.2Apr.2023利用改进交叉模型交叉模态的随机模型修正方法王炎1,陈辉1,2,黄斌1,柴满1(1.武汉理工大学土木工程与建筑学院,湖北 武汉 430070;2.武汉工程大学邮电与信息工程学院,湖北 武汉 430073)摘要:将混合摄动伽辽金方法和改进的交叉模型交叉模态技术相结合,提出了一种随机模型修正方法。该方法有效缓解了模型修正过程中测量数据有限和测量误差不确定的影响。考虑到实测模态数据具有不确定性,基于改进的交叉模型交叉模态方法,建立了一个新
2、的描述结构随机参数和随机响应关系的模型修正方程。利用混合摄动伽辽金方法求解该随机修正方程,进而得到结构随机修正参数的统计特征。简支梁的数值结果表明,该方法在测量数据不确定性较大时仍能保持很高的修正精度,同时计算效率比蒙特卡罗模拟方法高出一个数量级。在测量模态数据较少的情况下,该方法比单独的混合摄动伽辽金修正方法修正效果好,且比交叉模型交叉模态法的修正精度更高。框架试验的结果表明,该方法可以同时修正结构的刚度和质量,修正后的结构参数与预设工况基本吻合,同时能复现结构的测量模态,从而验证了所提方法的有效性。关键词:随机模型修正;随机混合摄动伽辽金方法;改进的交叉模型交叉模态方法中图分类号:O324
3、;TU311.4 文献标志码:A 文章编号:1004-4523(2023)02-0498-09 DOI:10.16385/ki.issn.1004-4523.2023.02.021引 言近几十年来,基于动力测量数据的有限元模型修正方法越来越受到关注。许多研究人员在这一领域进行了广泛的研究,并取得了大量研究成果14。在动力有限元模型修正中,修正参数的选择对修正结果有很大影响。如果修正参数过多,在修正过程中往往会出现病态问题,所以在修正模型之前首先要排除不敏感的修正参数5。关于动力有限元模型修正方法,Hu等6提出了一种基于交叉模型交叉模态(CMCM)方法的模型修正技术。与传统的模型修正方法不同,该
4、方法可以同时修正结构的刚度矩阵、质量矩阵和阻尼矩阵。此外,该方法不用迭代计算,计算效率较高。在 CMCM 方法中,通过将结构的实测模态和计算模态相乘,就可以仅用少量的测量模态构建多个模型修正方程。李世龙等7利用 CMCM 方法,提出了一种有效识别子结构边界约束状态的模型修正方法。Wang 等8使用了 CMCM方法对海上平台进行了试验研究,证明了当结构的实际测量模态不完整且只有低阶测量模态可用时CMCM 方法的有效性。在已有的 CMCM 方法的基础上,Liu等9提出了一种基于改进的交叉模型交叉模态(ICMCM)的模型修正方法,该方法充分利用实测数据,进一步增加了修正方程的个数。然而,这些方法仅涉
5、及确定性有限元模型修正,当结构参数的 不 确 定 性 或 者 测 量 噪 声 无 法 避 免 时,现 有 的CMCM 方法将不适用。因此,充分利用 CMCM 方法的优点,并将它融入随机模型修正中,是一项非常有意义的工作。在随机模型修正领域中,蒙特卡罗方法、摄动法以及贝叶斯方法被广泛使用。Schuller等10使用了具有大样本的蒙特卡罗模拟来计算模型修正的统计特性。宗周红等11在对下白石连续刚构桥进行模型修正的过程中,利用蒙特卡罗模拟方法和有限元方法进行不确定性量化分析,并评价模型的预测精度,实现对于连续刚构桥的有限元模型确认。但是对于大型结构而言,这种方法的计算效率太低,耗时过长。与蒙特卡罗方
6、法不同的是,摄动法具有推导简单、计算效率高的特点。Hua 等12使用一种改进的摄动法,利用随机实测模态数据对桁架桥有限元模型进行修正,并估计了结构参数的均值和均方差。尽管摄动法的计算效率比较高,但其对测量误差的变异性要求比较苛刻。随机模型修正方法中,另一个具有代表性的方法是基于马尔可夫链与蒙特卡罗抽样的贝叶斯方法1315,但是基于此种抽样的贝叶斯方法会面临较大的挑战,即需要非常耗时的重复有限元计算。为了提高计算效率,Wan等16和 Fang收稿日期:2021-07-26;修订日期:2021-12-05基金项目:国家自然科学基金面上项目(51978545)。第 2 期王炎,等:利用改进交叉模型交
7、叉模态的随机模型修正方法等17分别采用了高斯代理模型和随机响应面模型对原有的贝叶斯方法进行了改进。与上述方法不同,Huang 等18提出了一种基于混合摄动伽辽金方法(HPG)的随机模型修正方法(HPGSMUM),该方法在测量变异性较大情况下具有比较高的修正精度和效率,此方法也为确定性模型修正方法扩展到随机领域提出了一个新的思路和完整的框架。本文将随机摄动伽辽金方法与改进的交叉模型交叉模态方法结合,提出一种随机模型修正方法。该方法可利用含测量误差的少量模态测量数据实现结构有限元模型的有效修正。文中用一个简支梁的数值算例来验证该方法的有效性和不同模态组合的稳定性,并利用七层框架的模态试验来验证所提
8、方法在较少测量模态情况下仍能同时有效地修正结构刚度和质量。1基于 ICMCM 的随机模型修正方程考虑具有 N 个自由度的无阻尼结构,该结构初始模型满足以下特征值方程:Kai=iMai;i=1,nc(1)式中 Ka和Ma分别为初始结构模型的整体刚度矩阵和质量矩阵;i和i分别为初始模型的第 i阶特征值和特征向量;nc为初始模型的计算模态的个数。类似地,实际结构的特征值方程可以表示为:Kdj=jMdj;j=1,nm(2)式中 Kd和Md分别为实际结构模型的整体刚度矩阵和质量矩阵;j和j分别为实际模型的第 j阶特征值和特征向量;nm为实际模型的计算模态的个数。初始结构与实际结构的质量矩阵、刚度矩阵存在
9、以下关系:Md=Ma+n=1NenMn(3)Kd=Ka+n=1NenKn(4)式中 Ne为结构的单元个数;Kn和Mn分别为结构第 n 个单元的N N单元组装矩阵;n和n分别为结构第 n 个单元的刚度和质量的修正系数,表示实际结构的单元刚度和质量相对于初始矩阵的变化率。通过文献 6 可以得到确定性的基于交叉模型交叉模态的模型修正方程为:1+n=1NenTiKnjTiKaj=ji(1+n=1NenTiMnjTiMaj)(5)对式(5)进行因式变换可以得到:n=1NenTiKnjTiKaj-n=1NenjiTiMnjTiMaj=ji-1(6)通过求解式(6)所示的方程组可以得到刚度和质量的修正系数n
10、和n。但是由于在实际的模态测量中只能精确测量出前几阶模态,使得修正系数方程组的方程个数比较少,导致求解结果不正确且不稳定。因此 Liu等9对传统的 CMCM 方法进行改进,充分利用测量模态数据,在式(2)方程两边同时乘Tj,得到如下所示的基于 ICMCM 的模型修正方程:Ti(Ka+n=1NenKn)j=TjTi(Ma+n=1NenMn)j(7)对式(7)进行因式变换,可以得到:n=1NenTiKnj-n=1NenTiMnj=TjTiMaj-TiKaj(8)显然,式(8)也含有与方程(6)相同的修正系数,结合式(6)和(8),就可以得到更多的修正方程,确保修正方程的适定性。在实际结构的模态试验
11、过程中,不可避免地会遇到测量误差。假定第 j阶的特征值和特征向量可以表示为:j=0j+j1j(9)j=0j+j1j(10)式中 0j和0j分别为测量的第 j阶特征值和特征向量均值;1j和1j分别为第 j阶测量误差的确定性部分;j为与测量误差相关的随机变量,且随机变量的分布类型由实测数据的统计特征或者是工程经验决定。假设所有随机变量j完全相关,并且表示为随机变量,则第 n个刚度单元和质量单元的修正系数可以分别用下式表示:n()=n0+n1+n22+(11)n()=n0+n1+n22+(12)式中 ni(i=1,2,)和ni(i=1,2,)分别为展开式中的确定性系数。同样,将式(9)(12)代入到
12、式(6)和(8)中,可以得到基于改进交叉模型交叉模态的随机模型修正方程为:n=1Ne(n0+n1+n22+)TiKn(0j+1j)+n=1Ne-(n0+n1+n22+)jTiMn(0j+1j)=(0j+1j)TiMa(0j+1j)-TiKa(0j+1j)(13)499振 动 工 程 学 报第 36 卷n=1Ne(n0+n1+n22+)(0j+1j)TKn(0j+1j)+n=1Ne-(n0+n1+n22+)(0j+1j)(0j+1j)TMn(0j+1j)=(0j+1j)(0j+1j)TMa(0j+1j)-(0j+1j)TKa(0j+1j)(14)这里,式(13)和(14)也是关于ni(i=0,1
13、,2,)和ni(i=0,1,2,)的随机代数方程。2随 机 模 型 修 正 方 程 的 HPG 方 法求解首先采用高阶摄动法,递推求解对应多项式基0,1,2,的方程。考虑与式(13)和(14)中与0对应的项,有:n=1Nen0TiKn0j+n=1Ne(-n0jTiMn0j)=0jTiMa0j-TiKa0j(15)n=1Nen0T0jKn0j+n=1Ne(-n00jT0jMn0j)=0jT0jMa0j-T0jKa0j(16)将式(15)和(16)写为以下矩阵形式:C(0)E(0)(0)=f(0)(17)C(0)E(0)(0)=f(0)(18)式中(0)=(0)(0)T。式(17)和(18)中矩阵
14、C(0),E(0),C(0),E(0)以及向量f(0),f(0),(0)和(0)中的元素表示如下:C(0)pn=TiKn0j,E(0)pn=-0jTiMn0j,C(0)qn=T0jKn0j,E(0)qn=-0jT0jMn0j,f(0)p=0jTiMa0j-TiKa0j,(0)n=n0,f(0)q=0jT0jMa0j-T0jKa0j,(0)n=n0式中 下标 p 或 q 代表矩阵的第 p 或第 q 行;下标 n代表矩阵的第 n 列。假定测量的模态数为 s,则p=j+(i-1)s,q=j+(j-1)s。合并式(17)和(18),可得:C(0)I E(0)I(0)=f(0)I(19)式中 C(0)I
15、=|C(0)C(0),E(0)I=|E(0)E(0),f(0)I=|f(0)f(0)。类似地,考虑与式(13)和(14)中与1对应的项,有:C(0)E(0)(1)+C(1)E(1)(0)=f(1)(20)C(0)E(0)(1)+C(1)E(1)(0)=f(1)(21)式中(1)=(1)(1)T。式(20)和(21)中矩阵C(1),E(1),C(1),E(1)以及向量f(1),f(1),(1)和(1)中的元素表示如下:C(1)pn=TiKn1j,E(1)pn=-0jTiMn1j,C(1)qn=T0jKn1j+T1jKn0j,E(1)qn=-0jT0jMn1j-0jT1jMn0j-1jT0jMn0
16、j,f(1)p=(0j-i)TiMa1j+1jTiMa0j,f(1)q=1jT0jMa0j+0jT0jKa1j+0jT1jMa0j-T0jKa1j-T1jKa0j,(1)n=n1,(1)n=n1。将式(20)和(21)合并,可得:C(0)I E(0)I(1)+C(1)I E(1)I(0)=f(1)I(22)式中C(1)I=|C(1)C(1),E(1)I=|E(1)E(1),f(1)I=|f(1)f(1)。然后,考虑二阶多项式基2,可以得到:C(0)E(0)(2)+C(1)E(1)(1)=f(2)(23)C(0)E(0)(2)+C(1)E(1)(1)+C(2)E(2)(0)=f(2)(24)式中
17、(2)=(2)(2)T。式(23)和(24)中矩阵C(2),E(2)以及向量f(2),f(2),(2)和(2)中的元素表示如下:C(2)qn=T1jKn1j,f(2)p=1jTiMa1j,E(2)pn=-1jT0jMn1j-1jT1jMn0j-0jT1jMn1j,f(2)q=1jT0jMa1j+1jT1jM0j+0jT1jM1j-T1jK1j,(2)n=n2,(2)n=n2。综合式(23)和(24),有:C(0)IE(0)I(2)+C(1)IE(1)I(1)+C(2)IE(2)I(0)=f(2)I(25)式中 C(2)I=|0C(2),E(2)I=|0E(2),f(2)I=|0f(2)。通过求
18、解方程(19),(22)和(25),可以得到向量(0),(1)和(2)。同理可以求出更高阶的系数向量。之后,采用伽辽金投影法提高含有ni(i=0,1,2,)和ni(i=0,1,2,)的摄动解精度。基于式(19),(22)和(25),将方程(13)和(14)合并,可以有:C(0)IE(0)I(0)+(1)+(2)2+)+C(1)IE(1)I(0)+(1)+(2)2+)+C(2)IE(2)I(0)+(1)+(2)2+)2=f(0)I+f(1)I+f(2)I2(26)进一步,假定结构的修正因子向量为:=i=0mii(27)式中 i(i=0,1,m)为新展开基i的系数;基500第 2 期王炎,等:利用
19、改进交叉模型交叉模态的随机模型修正方法向量i分别为(0),(1),(2)2,;m为展开阶数。为了求出系数i,可以将式(27)代入到式(26)中,并且在式(26)的左右两边左乘(C(0)IE(0)I k)T(k=0,1,m),则可得确定性的代数方程为:i=0m-CE(0)kii+i=0m-1-CE(1)kii+i=0m-2-CE(2)kii=fk(28)式中-CE(0)ki=,-CE(1)ki=,-CE(2)ki=,fk=i=02式 中 是 一 个 取 数 学 期 望 的 符 号。显 然,-CE(0)ki,-CE(1)ki,-CE(2)ki和fk为 标 量,方 程 组(28)包 含m+1 个待定
20、系数i。通过求解(m+1)(m+1)维线性方程组,就可以得到i。本文所提出的随机模型修正方法的流程图如图 1所示。上述方法就是本文所提出的结合 HPG 和 ICMCM 的随机模型修正方法(HPGICMCM)。假设用 C(0)E(0)和 C(1)E(1)代替式(19),(22)和(25)中的 C(0)IE(0)I和 C(1)IE(1)I,相应的向量(0),(1),(2)也可以通过上述方式递推得到。此时,HPGICMCM 方法退化为 HPGCMCM 方法。需要注意的是,实际结构的转角模态往往难以测量。此外,由于测量条件的限制,仅能测量包括部分测点的振型。因此,本文使用文献 19 的模态扩阶方法得到
21、完备振型的均值和标准差。同时,在求解方程组(19),(22)和(25)的过程中,采用截断奇异值分解或者 L1正则化技术20避免方程病态的问题。3数值算例考虑一个简支梁,如图 2 所示。简支梁的跨度为 6 m,截 面 为0.2 m 0.25 m,弹 性 模 量 为2.8 1010 Pa,密 度 为2.5 103 kg m3。将 该 Euler Bernoulli梁的有限元模型沿梁长度方向划分为 15个相同的单元,每个节点包含竖向位移和转角两个自由度。根据测量经验,可以假设实测模态数据服从某种概率分布,如正态分布或者分布。由于实测数据是有界的,本文假设实测模态数据服从分布。依据工程经验,动力特性测
22、试数据的变异系数一般在 0.010.02之间,这里假设变异系数为 0.02。首先,考虑结构质量不发生变化,单元 1,3,7,9和 11 的 刚 度 分 别 减 小 30%,15%,20%,20%和30%,其余单元的刚度和初始模型相同。取前六阶初始模型的计算模态和前六阶测量模态,使用 HPGICMCM 方法和 HPGCMCM 方法对模型进行修正,同时利用与本文所提出的方法对应的蒙特卡罗模 拟 方 法(MC ICMCM)和 Huang 等18的 HPG SMUM 方法求解上述方法中的修正系数的统计特性,在求解过程中使用奇异值分解正则化技术降低矩阵求逆的不适定性,以提高计算精度。修正结果如图 3和
23、4所示。图 1 基于 HPGICMCM 方法的随机模型修正流程图Fig.1 Flowchart of stochastic model updating by means of the HPGICMCM method图 2 简支梁的有限元模型Fig.2 Finite element model of the simply supported beam图 3 单元刚度修正系数均值Fig.3 Mean of updated coefficient of element stiffness图 4 单元刚度修正系数标准差Fig.4 Standard deviation of updated coeff
24、icient of element stiffness501振 动 工 程 学 报第 36 卷观察图 3 和 4,不难看出,当结构的自由度比较多但测量模态有限时,通过 HPGSMUM 方法得到的修正系数均值和预设的真实值差别比较大。例如,HPGSMUM 方法得到的单元 1 刚度修正系数均值为0.02,和 MCICMCM 方法结果相比,相对误差接近 90%。同时,单元 2,15的刚度修正系数均小于0.1,出现了明显误判。对于单元 8,10,13 和15,HPGSMUM 方法的修正系数标准差结果和MCICMCM 方法结果最大相对误差达到 400%,说明在这种情况下 HPGSMUM 方法修正效果不能
25、令人 满 意。而 通 过 HPG ICMCM 方 法 和 HPG CMCM 方法得到的各单元修正系数与 MCICMCM相比较,均值的绝对误差均未超过 0.03,标准差的相对误差基本小于 20%。可以说明统计结果和仿真试验预设的结果基本吻合,并且 HPGICMCM 方法的均值结果吻合更好。为了验证同时修正质量和刚度时本文方法的有效性,假设单元 3,5,6,8,9,11和 13的实际质量分别增加 10%,20%,20%,20%,20%,20%和 10%,同时,单元 1,3,5,7,9,11 和 13 的弹性模量分别降低30%,15%,20%,20%,20%,30%和 30%,其余单元的质量和刚度和
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 利用 改进 交叉 模型 随机 修正 方法 王炎
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。