结构恢复力非参数化模型识别的改进容积卡尔曼滤波方法_杜义邦.pdf
《结构恢复力非参数化模型识别的改进容积卡尔曼滤波方法_杜义邦.pdf》由会员分享,可在线阅读,更多相关《结构恢复力非参数化模型识别的改进容积卡尔曼滤波方法_杜义邦.pdf(11页珍藏版)》请在咨信网上搜索。
1、第 36 卷第 2 期2023 年 4 月振 动 工 程 学 报Journal of Vibration EngineeringVol.36 No.2Apr.2023结构恢复力非参数化模型识别的改进容积卡尔曼滤波方法杜义邦1,许斌1,2,赵冶1,邓百川3(1.华侨大学土木工程学院,福建 厦门 361021;2.福建省智慧基础设施与监测重点实验室(华侨大学),福建 厦门 361021;3.纽约州立大学布法罗分校土木、结构与环境工程系,纽约 布法罗 14260)摘要:对地震等强动力荷载作用过程中结构损伤的发生发展过程进行识别,必须考虑结构行为的非线性。本文运用相对位移和相对速度的幂级数多项式表征结
2、构恢复力模型,提出一种基于改进的容积卡尔曼滤波算法(Updated Cubature Kalman Filter,UCKF)和结构部分自由度上加速度响应时程的结构参数、未知响应及恢复力非参数化模型识别方法。以一个含磁流变阻尼器的多自由度数值模型为例,考虑 20%的加速度响应测量噪声影响,识别出模型的结构参数、未知响应及阻尼力。并将本文方法所得结果分别与基于扩展卡尔曼滤波算法、传统容积卡尔曼滤波算法及含记忆衰退的扩展卡尔曼滤波算法所得结果进行比较。对一个带磁流变阻尼器的四层剪切型框架模型进行激振试验,基于部分自由度上的加速度响应时程实测值,识别出结构参数、未知动力响应以及阻尼器阻尼力的非参数化模
3、型,通过与实测结果的比较,验证了本文方法的可行性。关键词:非线性恢复力;改进的容积卡尔曼滤波;非参数化识别;幂级数多项式;磁流变阻尼器中图分类号:TU311.3 文献标志码:A 文章编号:1004-4523(2023)02-0389-11 DOI:10.16385/ki.issn.1004-4523.2023.02.010引言土木工程结构在服役过程中,除了因材料劣化等因素会出现结构损伤和性能退化外,也往往会因地震等强动力荷载的作用出现不同程度损伤甚至破坏。对强动力荷载作用过程中结构损伤的发生发展过程进行识别,并据此对结构剩余承载力和剩余寿命进行预测是亟待解决的课题。严格来讲,基于结构特征值和特
4、征向量抽取的结构识别方法仅适用于线性系统,可近似应用于材料劣化等引起的结构刚度缓慢变化情况下的识别问题。但结构在地震等强动力荷载作用过程中损伤的发生发展过程是一个典型的非线性过程,不同结构构件在不同时刻进入非线性阶段,构件和结构层次的刚度并非保持不变,运用特征值或特征向量来识别刚度,并用其描述不断发展的损伤存在不合理性。不同于刚度,结构或构件在振动过程中的恢复力是其非线性行为的最直接描述,结合结构动力响应识别结果,不仅可以评价不同构件在不同时刻和不同变形情况下的损伤状态,还可以定量计算振动过程中的耗能1-2。此外,由于实际土木工程结构材料和类型的多样性,结构恢复力往往难以通过某一事先假定的统一
5、参数化数学模型来准确描述。因此,针对强动力荷载作用后或灾后结构损伤的识别问题,通过建立恢复力的非参数化模型并识别结构恢复力具有重要意义。相较于线性结构,由于结构材料和类型的多样性,非线性结构的识别问题难度更大,因而得到国内外学者的重视。针对非线性结构的识别问题,Masri等3提出了恢复力曲面法并将其推广应用到了多自由度非线性动力系统。随后,Smyth 等4提出了一种基于最小二乘的结构参数的在线自适应识别方法,并估计了结构非线性恢复力。Xu等5-6基于线性等效思想,运用最小二乘法识别了结构恢复力,并运用一个具有磁流变(Magneto Rheological,MR)阻尼器的多自由度剪切型框架模型的
6、动力试验实测数据验证了该方法的可行性。考虑实际结构恢复力参数化模型难以预先确定的问题,许斌等7-8分别运用幂级数多项式和切比雪夫多项式作为结构非线性恢复力的非参数化模型,提出了恢复力的非参数化识别方法,并分别通过含有 MR 阻尼器和形状记忆合金阻尼器(Shape Memory Alloy,SMA)的非线性多自由度系统的数值模拟和模型试验实测数据,验证了所提出方法的有效性。收稿日期:2021-09-27;修订日期:2022-02-02基金项目:国家自然科学基金资助项目(50978092);华侨大学科研基金资助项目(605-50Y18016)。振 动 工 程 学 报第 36 卷基于状态空间模型递推
7、的卡尔曼滤波类算法被广泛用于解决实际工程结构中激励与系统动力响应信息不完全已知情况下结构的识别问题。Jaz-winski9利用扩展卡尔曼滤波方法(Extended Kal-man Filter,EKF)对加速度测量进行滤波,并通过预测与估计识别得到了结构参数。Hoshiyam 等10提 出 了 加 权 全 局 迭 代 扩 展 卡 尔 曼 滤 波 算 法(Ex-tended Kalman Filter-Weighted Global Iteration,EKF-WGI)。王祥建等11又引入记忆衰退技术提高了 EKF 算法的稳定性。Lei 等12-13和 Liu 等14基于等效线性方法,实现了部分
8、观测下结构参数及非线性恢复力的识别。张肖雄等15-16则通过引入投影矩阵,推导了改进的观测方程,实现了外激励的实时识别。Xu 等17基于勒让德多项式和位移与加速度的数据融合方法,对结构恢复力和未知激励进行识别。EKF 通过一阶泰勒级数将非线性函数线性化,该方法存在精度不高、易于发散以及只适用于弱非线性系统的缺点。由于对概率分布进行近似要比对非线性函数进行近似容易得多,基于该思想,Julier等18-19提出了无迹变换(Unscented Transform,UT)与无迹卡尔曼滤波算法(Unscented Kalman Filter,UKF)。UKF对非线性函数的概率密度分布进行近似,用一系列确
9、定样本来逼近状态的后验概率密度。Wan等20验证了 UKF 的有效性。Wu等21-22的研究结果表明,在噪声更大的情况下对高维系统的参数进行估计时,UKF 相较于 EKF 具有更高的识别精度。Xie 等23运 用 迭 代 无 迹 卡 尔 曼 滤 波(Iterated Unscented Kalman Filter,IUKF)对 Bouc-Wen滞回系统进行识别,结果表明,相较于 UKF,IUKF 的识别精度更高和算法鲁棒性更好。针对地震作用下结构刚度突变的识别问题,Bisht等24提出了一种具有跟踪结构参数突变能力的自适应无迹卡尔曼滤波算法(Adaptive Unscented Kalman
10、Filter,AUKF),数值模拟结果表明,该方法能对多个结构构件在不同时刻的参数突变进行有效识别。实际应用卡尔曼滤波时,模型误差、噪声误差、计算误差均可能会造成预测误差协方差矩阵和增益矩阵随迭代次数增加而减弱修正状态估计的情况(即数据饱和或者观测老化),并进而导致滤波发散。针对此问题,渐消记忆滤波的扩展卡尔曼滤波算法(Extended Kalman Fil-ter-Memory Fading,EKF-MF)被提出,该算法通过增大新数据的作用而降低旧数据的负面影响25-26。为准确估计系统的噪声统计特性,Sage等27提出了一种可实时估计系统及测量噪声的自适应滤波算法,该方法在滤波的同时利用观
11、测量信息,对模型参数、噪声特性进行实时调整,进一步提高了滤波精度。为克服 UFK 在高维系统中出现滤波精度低和数值计算不稳定的问题,加拿大学者 Arasaratnam等28于 2009 年首次提出了基于 Cubature 变换的容积 卡 尔 曼 滤 波 算 法(Cubature Kalman Filter,CKF)。基于三阶球面径向容积准则(Cubature 准则),使用 2n 个(n 为扩展状态向量的维数)权值相同的容积点(Cubature 点),经非线性系统方程转换后进行加权处理,来逼近具有附加高斯噪声的非线性系统的状态均值和协方差。常宇健等29对比了EKF,UKF 和 CKF 的滤波性能
12、,结果表明 EKF 的滤波性能最差,而 UKF 对于高维非线性系统不仅调解参数困难,且在滤波过程中可能会出现协方差非正定的情况,使滤波结果不稳定甚至发散,而CKF 的滤波结果优于 EKF 与 UKF。孙枫等30比较了 CKF 与 UKF 的滤波精度,并通过数值模拟验证了对于高维(n3)的非线性系统,CKF 的精度、稳定性均高于 UKF,建议选择 CKF 作为滤波方法。CFK 是理论上当前最接近贝叶斯滤波的近似算法,是解决非线性系统状态估计的强有力工具。其中,将积分形式变换成球面径向积分形式和三阶球面径向准则是最为重要的两个步骤。CKF 方法一经提出便在姿态估计、导航、连续系统和混合滤波等领域得
13、到应用30。本文将记忆衰退权重、奇异值分解及 Sage-Hu-sa自适应滤波算法引入 CKF,提出一种运用改进的容 积 卡 尔 曼 滤 波 算 法(Updated Cubature Kalman Filter,UCKF)及部分加速度测量,适用于结构质量、刚度和阻尼参数均未知情况,不依赖于结构恢复力参数化模型的非线性结构识别方法。该方法基于UCKF,利用幂级数多项式表征结构非线性恢复力的非参数化模型,利用部分加速度测量响应,预测结构质量、刚度、阻尼系数,识别未知速度和位移响应以及幂级数多项式的系数,建立恢复力的非参数化模型,进而可以得到结构在振动过程中恢复力时程曲线及滞回曲线。为验证所提出方法的
14、可行性,首先建立了一个含 MR 阻尼器的四自由度的集中质量非线性数值模型,考虑 20%的加速度测量噪声的影响,通过本方法的识别结果与其理论值的比较,验证本文所提出方法的有效性与精确性。其次,将基于本 文 UCKF 方 法 的 识 别 结 果 与 基 于 EKF,CKF,EKF-MF 方法的识别结果进行比较,验证了本文方法的优点。最后,通过一个带 MR 阻尼器的四自由度非线性钢框架模型的动力试验,利用部分自由度上的加速度响应观测值,对模型的质量、刚度、阻尼系数、未测量动力响应时程及 MR 阻尼器阻尼力的非参数化模型进行识别,通过比较识别所得 MR 阻尼力与与试验实测值,进一步验证了本文方法的识别
15、结果的准确性。390第 2 期杜义邦,等:结构恢复力非参数化模型识别的改进容积卡尔曼滤波方法1基于 UCKF 与幂级数多项式的恢复力非参数化模型识别方法1.1等效线性化在外激励作用下,一个含非线性元件的多自由度动力系统的运动平衡方程可写为:Mx?(t)+Cx?(t)+Kx(t)+fnon(t)=f(t)(1)式中M,K和C分别为结构质量、刚度和阻尼矩阵;x?(t),x(t)和x?(t)分别为加速度、位移和速度响应向量;fnon(t)为非线性元件的恢复力;f(t)为外激励。将式(1)进行等效线性化31,可得到如下方程:MEx?(t)+CEx?(t)+KEx(t)=f(t)(2)式中ME,CE和K
16、E分别表示等效线性结构的质量、阻尼和刚度矩阵。在结构非线性的发展过程中,可合理认为结构质量是不变的,故可将ME视为结构质量的识别值。一般而言,式(1)中结构非线性恢复力会体现在等效线性结构参数CE与KE的变化上。1.2恢复力非参数化模型的幂级数多项式表示式(2)可进一步表达为:Mx?(t)+Rnon(t)=f(t)(3)式中Rnon(t)为结构总非线性恢复力,可表示为:Rnon(t)=CEx?(t)+KEx(t)(4)对于链式结构,两个自由度间的非线性恢复力可通过一组相对速度及相对位移的幂级数多项式来表达31,如下式所示:Rnoni,i-1(t)k=1Kj=0J krnoni,i-1,k,jv
17、ji,i-1sK-ji,i-1(5)式中Rnoni,i-1(t)表示结构第i个与第i-1个自由度间总非线性恢复力;vi,i-1与si,i-1分别表示结构第i个与第i-1个自由度间的相对速度及相对位移;rnoni,i-1,k,j为幂级数多项式的系数。J取整数,K 的取值与结构的非线性程度相关,本文中取 K=3。根据式(5),结构第i个自由度的运动平衡方程可离散为:mix?i+k=1Kj=0J krnoni,i-1,k,jvji,i-1sK-ji,i-1+k=1Ki=0J krnoni,i+1,k,jvji,i+1sK-ji,i+1=fi(t)(6)式中mi为第 i自由度的层间质量。式(6)中结构
18、质量、部分未知加速度响应、速度和位移响应以及恢复力非参数化模型中的参数有待识别。识别所得幂级数多项式即为结构非线性恢复力的非参数化模型。将识别所得速度和位移响应代入该模型即为恢复力的识别结果,可分别与数值模拟中的理论值和实验验证中的实测值进行比较。需要指出的是,本文的数值模拟及试验验证均选择在线弹性结构中引入 MR 阻尼器来模拟结构的非线性行为。此时,MR 阻尼器的阻尼力Fnon(t),需要从所识别的Rnon(t)中减去线性结构本身所提供的线弹性恢复力与阻尼力,即Fnon(t)=Rnon(t)-Cx?(t)-Kx(t)(7)在本文的数值模拟和试验验证中,Fnon(t)将分别与数值模型中阻尼力的
19、理论值及试验中阻尼力的实测值进行比较,以验证本文方法的准确性。1.3基于 UCKF的非线性结构识别方法非线性结构在动力作用下的系统状态方程和观测方程可表示为:Xk=f()Xk-1+wk-1(8)Zk=h()Xk+vk(9)式中k 为时间步数;Xk为第 k 步的状态值;f()表示系统的状态转移函数,服从一阶马尔可夫假设;wk-1表示协方差为Qk-1的系统过程噪声;Zk为第 k步的观测值;h()表示系统的观测函数,服从观测独立假设;vk表示协方差为Rk的系统观测噪声。若k-1时刻的后验概率p()Xk-1Z1:k-1 N()Xk-1;X?k-1,Pk-1,其 中,X?k-1为k-1时刻状态估计值,P
20、k-1为k-1时刻误差协方差矩阵,则 UCKF的递推算法如下:(1)基本容积点及对应权值的计算根据三阶容积原则,有:j=l 2 1 j,j=1/l(10)l=2nx(11)式中j表示容积点序号,l表示容积点总数,取j=1,2,l;j和j分别为基于三阶容积规则获得的第 j个基本容积点和相应权值;nx表示系统维数。记nx维单位向量为e=1,0,0,0 T,1表示对e的元素进行全排列和改变元素符号产生的点集,称为完整全对称点集,1 j表示点集中的第 j个点32-33。(2)时间更新计算 k1时刻第 j个容积点Xj,k-1:Xj,k-1=Sk-1j+Xk-1(12)式中Sk-1为Pk-1的平方根形式,
21、可根据Pk-1=Sk-1STk-1得到。但在实际情况中,由于矩阵Pk-1往往非半正定,使用传统的 Cholesky 变换会导致滤波发散,故在 本 文 UCKF 算 法 中,对 矩 阵Pk采 用 奇 异 值 分解34,即Sk-1=svd()Pk-1(13)式中svd()为求解的奇异值分解矩阵。计 算 k-1 时 刻 经 系 统 方 程 传 递 后 的 容391振 动 工 程 学 报第 36 卷积点X*j,k:X*j,k=f()Xj,k-1(14)预测 k时刻状态值Xk:Xk=j=1ljX*j,k(15)预测 k时刻状态协方差矩阵P(0)k:P(0)k=j=1ljX*j,kX*Tj,k-XkXTk
22、+Qk-1(16)式中P(0)k为引入记忆衰退权重系数 G 之前的状态协方差矩阵。对于 UCKF,考虑到滤波发散的原因及卡尔曼滤波记忆无限增长的特点,本文使用记忆衰退技术在滤波过程中增大新数据在观测中的作用,降低旧数据的负面影响。于是,在式(16)的P(0)k中引入记忆衰退权重系数 G,写为Pk,有:Pk=()j=1ljX*j,kX*Tj,k-XkXTkG+Qk-1(17)式中Pk为引入记忆衰退权重系数 G 后的状态协方差矩阵。G 根据文献 25-26 的建议取1 G 1.05。(3)量测更新计算 k时刻容积点Xj,k:Xj,k=Skj+Xk(18)计算 k时刻经量测方程传递后的容积点Zj,k
23、:Zj,k=h()Xj,k(19)估计 k时刻观测值Zk:Zk=j=1ljZj,k(20)估计 k 时刻状态协方差矩阵Pzz,k、互协方差矩阵Pxz,k:Pzz,k=j=1jZj,kZTj,k-ZkZTk+Rk(21)Pxz,k=j=1ljXj,kZTj,k-XkZTk(22)计算 k时刻滤波增益Wk:Wk=Pxz,kP-1zz,k(23)计 算 k 时 刻 状 态 估 计 值X?k和 误 差 协 方 差矩阵Pk:X?k=Xk+Wk()zk-z k(24)Pk=Pk-WkPzz,kWTk(25)在本文 UCKF 中,为增加滤波抗噪性能并提高滤波精度,引入 Sage-Husa 算法在估计当前噪声
24、的同时对观测噪声方差进行实时修正。根据简化的Sage-Husa 次优无偏极大后验估计器,可对Rk和Qk作出估计35-36:Rk=()1-dkRk-1+dk()kTk-j=1ljZj,kZTj,k-ZkZTk(26)Qk=()1-dkQk-1+dk|WkkTkWTk+Pk-()j=1ljX*j,kX*Tj,k-XkXTk(27)其中:k=Zk-Zk(28)dk=1-b1-bk(29)式中b为遗忘因子,根据文献 35 可确定 b的取值范围为 0.95b0.99。b 越大,则新量测数据对实时估计的作用越明显,故当噪声统计特性变化速度较快时,b的值应向 0.99靠近;当噪声统计特性变化速度较慢时,b的
25、值应向 0.95靠近。在 k 时刻至 k+1 时刻的计算中,根据式(26)计算所得到的Rk代入式(21)中,进行后续计算,并按上述步骤对Rk进行更新。1.4算法实现本文提出的识别方法的具体实现步骤如下:(1)假设结构质量、刚度、阻尼系数初始值;(2)运用等效线性化及 UCKF 算法,在已知部分加速度观测量的条件下,对结构各物理参数及未知结构速度、位移响应时程进行识别;(3)收敛判断:本文采用一种弱化的收敛判断条件。结构各层前后两次质量识别的差与前一次识别质量的比值,其最大值小于 1%,或上述各层质量前后误差的绝对值之和小于 3%,停止迭代,否则将上一次识别值当做下一次识别初始值,并循环以上步骤
- 配套讲稿:
如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。