随机子空间模态参数识别
模态参数辨识改进算法及其在振动台试验中的应用
( ; a . C o l l e eo fC i v i lE n i n e e r i n b . C o l l e eo fC o n s t r u c t i o nM a n a e m e n t a n dR e a l g g g g g , , , ) h o n i n n i v e r s i t h o n i n 0 0 0 4 5 C h i n a E s t a t e C g q gU y C g q g4 : ) , , 犃 犫 狊 狋 狉 犪 犮 狋 B a s e do nt h ed a t a d r i v e ns t o c h a s t i cs u b s a c ei d e n t i f i c a t i o n( S S I a ni m r o v e dI T Dm e t h o d a n p p i m r o v e dS T Dm e t h o da n da n i m r o v e dc o m l e x i n d e xm e t h o da r ep u t f o r w a r d . S S Ih a sh i h i d e n t i f i c a t i o n p p p g , a c c u r a c a n dH a n k e lm a t r i xo r t h o o n a lp r o e c t i o nc a l c u l a t i o no fd a t ac a nr e t a i na l l t h eu s e f u l i n f o r m a t i o n y g j i nt h eo r i i n a ld a t aa n ds i m u l t a n e o u s l e m o v en o i s e .T h ed a t ao b t a i n e di nP m a t r i xi st a k e na s i n u tf o r g yr p , t h e I T D, S T Da n dc o m l e xi n d e xm e t h o d s a n dt h e nt h et h r e em e t h o d sw o u l dn ol o n e rn e e dr a n d o m p g , t h e r e b l i m i n a t i n h ee r r o rc a u s e db h et w o d e c r e m e n t t e c h n i u eo rN E x tm e t h o df o rp r e r o c e s s i n ye gt yt q p g r e t r e a t m e n tm e t h o d s . T h ei m r o v e dm e t h o d sa r ea l i e dt ot h em o d a lp a r a m e t e ri d e n t i f i c a t i o nf o r p p p p s t r u c t u r es h a k i n t a b l e t e s tm o d e l b s i n t h e s i m u l a t e de a r t h u a k e r e s o n s ed a t a f r o mT o n iU n i v e r s i t . g yu g q p g j y T h er e c o n i t i o nr e s u l t s f u r t h e rp r o v e t h ec o r r e c t n e s sa n df e a s i b i l i t f t h e i m r o v e dm e t h o d s . g yo p : ; ; ; 犓 犲 狅 狉 犱 狊f r a m es t r u c t u r e m o d a lp a r a m e t e r si d e n t i f i c a t i o n i m r o v e dI T D m e t h o d i m r o v e dS T D p p 狔狑 ; m e t h o d i m r o v e dc o m l e x i n d e xm e t h o d p p 不 对土木结 构 进 行 实 时 安 全 监 测 或 振 动 控 制 , 更需要充分了解其动 仅需要掌握结 构 的 静 力 特 性 , 力特性 , 以避 免 灾 难 事 故 的 发 生 。 而 结 构 的 模 态 参 数( 频率 、 阻尼 、 振 型) 是 结 构 动 力 特 性 的 主 要 参 数, 所以如何准确识别结构的模态参数 是研 究 结构 动力
一种分段平稳随机信号的参数识别方法
中图分类号 :U 4 41 文献标 识码 :A
f + 1=Ax +
参数识别是结构健康监测领域 中的重点¨ 。随机 子空间方法是近年来发展起来 的一种行之有效 的结构 模态 参数 识 别 方 法 。 由于 它 具 有 如 下 优 点 : 只需 一 个 给定 的 “ 数 ” 即 系统 阶 次 ; 参 , 由于 该 方 法 中所 用 到 的 Q R分解和奇异值分解 (V ) S D 不需要进行迭代运算 , 因 此不 存在 是 否 收 敛 的 问题 等 , 来 越 受 到 国 内外 业 内 越 人士 的广 泛关 注 。 随机 子 空 间方 法 的 基 本 假 定 是 输 入 为 白噪 声 , 输 出为平稳随机信号 , 方法 全部算法都 是 以这个假 定 该 为前 提 的 。实 际应 用 中很 难 完 全 满 足 这 个 假 定 , 了 为 尽可 能满 足 白噪声 假 定 , 测 时 需 要 记 录 尽 可 能 长 的 实 时间, 因为有 时认 为在 一 个 较 长 的时 间 内 , 梁 环 境 激 桥 励近似为一个稳态 的白噪声过 程。但这种假设并不是 绝对 的 , 时 恰好 相 反 , 一 个 比较 长 的时 间段 内认 为 有 在 信号是非平稳的 , 而在 比较短的时 间段 内认为信 号是 平稳 的 。这就 是所 谓 的“ 分段 平 稳 ” 。 为了能将随机子 空间方法 的应用前提 拓展 , 其 使 能够识别 比较实际 的非平稳信 号 , 文提 出一种 分析 本 分段平稳随机信号两阶段稳定图的方法 。
维普资讯
振 第2 6卷第 6 期 B A I N A D S C OU N F VI R T O N HO K
一
大型工程结构模态识别的理论和方法研究
引言
齿轮箱是机械设备中的重要组成部分,其性能和可靠性对整个设备的运行有 着重要影响。在齿轮箱的设计和制造过程中,模态分析和结构优化是提高其性能 和可靠性的一种重要手段。本次演示旨在研究齿轮箱模态分析和结构优化方法, 以提高其性能和可靠性,为机械设备的安全稳定运行提供保障。
文献综述
齿轮箱模态分析主要是通过分析齿轮箱的振动特性,预测其在实际运行中的 振动响应,从而为结构优化提供依据。结构优化则是在模态分析的基础上,对齿 轮箱的结构进行改进,以降低其振动和噪声。传统的模态分析方法主要包括有限 元法和实验法,其中有限元法是通过计算机模拟来进行分析,实验法则通过实验 测试来进行分析。
实验设计与实现
实验设计是模态识别的重要环节,包括数据采集、数据预处理、特征提取和 模式识别等步骤。首先,针对大型工程结构的特点,选择合适的测点布置和数据 采集设备,获取结构在振动过程中的动态响应。然后,对采集的数据进行预处理, 如去除噪声、填充缺失值等,以便后续分析。接下来,采用适当的特征提取方法, 如傅里叶变换、小波变换等,提取数据的时频特征。
2、进行模态分析:利用有限元软件对齿轮箱进行模态分析,得到其固有频 率和振型。
3、设计实验方案:根据模态分析的结果,设计实验方案,包括实验测试系 统和测试数据处理方法。
4、进行实验测试:利用实验测试系统对齿轮箱进行振动测试,并记录测试 数据。
5、分析实验结果:将实验测试数据与有限元分析结果进行对比,分析误差 和影响因素,优化有限元模型和实验方案。
最后,利用模式识别技术,如支持向量机、随机森林等,对提取的特征进行 分类和识别,得到结构的模态参数。
在实验过程中,需要注意以下几点:
1、数据采集设备的性能和精度对模态识别的结果有很大影响,应选择适合 大型工程结构特点的高效、可靠的设备。
基于Welch法的协方差随机子空间方法的模态参数识别
基于Welch法的协方差随机子空间方法的模态参数识别李雪艳;官宇航;罗铭涛;吴博宇【期刊名称】《力学学报》【年(卷),期】2022(54)10【摘要】对工程结构进行环境激励下的模态参数识别具有重要意义,而随机子空间法作为适合环境激励下模态参数识别的时域方法,由于噪声和复杂激励的原因,会产生虚假模态、真实模态遗漏、系统自动定阶难和计算效率等问题,这些问题阻碍了该方法在实际工程中的广泛应用.本文提出了基于Welch法的随机子空间方法,通过Welch法对振动响应在频域进行去噪、降低环境激励和其他不确定性因素影响的处理,把结构固有模态从噪声和激励频率中突显出来,形成富含更多结构模态的Toeplitz矩阵,然后进行奇异值分解和状态矩阵计算,最后进行特征值分析.为了实现自动定阶,对不同奇异值分量构建的状态矩阵得到的特征参数,进行模糊C均值聚类分析和模态的平均相位偏移分析,剔除虚假模态,实现结构模态参数的自动识别.并把本文所提出方法应用于一座大跨悬索桥的实测加速度响应分析,和一座七十层的高层建筑的加速度响应分析,跟频域分解法、传统随机子空间法和基于相关分析的随机子空间法的计算结果进行了比较,发现基于Welch方法的随机子空间法相比于传统随机子空间法和基于相关分析的随机子空间法,在避免模态遗漏和计算效率方面有显著提高,而相对于频域分解法则在自动识别和剔除虚假模态方面有明显优势.【总页数】12页(P2850-2860)【作者】李雪艳;官宇航;罗铭涛;吴博宇【作者单位】暨南大学力学与建筑工程学院【正文语种】中文【中图分类】TU375.4【相关文献】1.协方差驱动随机子空间算法的梁桥模态参数识别2.基于数据驱动的应变模态参数随机子空间识别法3.基于锤击激励的随机子空间法模态参数识别4.基于随机子空间法的城市桥梁模态参数识别研究与程序开发5.基于随机子空间法的风机模态参数识别因版权原因,仅展示原文概要,查看原文内容请购买。
协方差驱动随机子空间算法的梁桥模态参数识别
【 关键词 】 环境振动 ; 随机子空间识别 ; 协 方差; 模 态参数 【 中图分类号】 U4 . 41 3
不同于传统的结构模 态参数识别 , 于输 出的模 态参数 基
识别 是 在 激 励 未 知 的情 况 下 对 结 构 系 统 进 行 辨 识 。针 对 此
【 文献标识 码】 A
参 数 识 别 中的 应 用 进 行 了研 究 。首 先 由 Q R分 解计 算 由 系统 输 出组 成 的 H ne 矩 阵 , 构 建 基 于输 出的 自 akl 并 协 方 差 Tp lz 阵 , 对 该矩 阵进 行 奇 异值 分 解 ( V , 取 系统 的 可观 测 矩 阵及 可控 制 矩 阵 , 据 协 方 差 oet 矩 i 再 S D) 获 依
出协 方 差 驱 动 随 机 子 空 间 模 态 参 数 识 别 流 程 图 , 文 仅 从 本 H n e 矩 阵 和 系 统 阶 数 的 确 定 以 及 模 态 参 数 的 提 取 等 方 面 a kl
进行讨论 。
2 试验 验证
按上述 SI 法 , MA L S算 以 T B为开发平 台 , 编制 了相 应的
驱动的随机子空 间算法构建结构 离散状 态矩 阵并进行 特征值 分解 , 而识 别结构 的模 态参 数。编制 了协 方 从 差驱动随机子空间算法识别桥梁结构模 态参数的 Ma a tb程序 , l 以某三跨连续 梁桥 梁的模 态参数识 别为例 , 验 证 了方法的可靠性和可行性。结果证明 , 方差驱动 随机 子空间算法可应用于实际桥 梁的模态参数识别 中。 协
理想的效果 。
由于系统状态 向量不一定有特定 的物理意义 , 因此状态 空间模型的特征向量 必 须转化 为结 构的振 型。可以计算 出结构的模态参数第 i 阶固有频率, 、 i阻尼 比 s 以及振型 ;
模态参数识别原理
模态参数识别原理
模态参数识别是一种结构动力学分析技术,它是通过对结构系统进行激励和响应的测量,来估计结构系统的振动特性。
模态参数识别的目的是确定结构体系的固有频率、阻尼和振动模态(模态形状),这些参数可以用来评估结构的稳定性、安全性和可靠性。
模态参数识别的原理是通过结构系统的振动响应,采用最小二乘法、奇异值分解法、支持向量机、神经网络等数学方法,来计算结构系统的固有频率、阻尼和振动模态。
在实际应用中,结构系统的振动响应可以通过传感器、激励器和信号分析仪等设备来获取,这些设备可以分别安装在结构系统的不同位置,通过测量响应信号的时程和频谱特征,来计算结构系统的模态参数。
模态参数识别的应用领域非常广泛,包括工程结构的监测、损伤诊断、结构优化设计等方面。
在实际应用中,由于结构系统的复杂性和多变性,模态参数识别存在一定的难度和挑战,因此需要结合实际情况选用合适的方法和技术,来保证识别结果的准确性和可靠性。
基于谱系聚类的随机子空间模态参数自动识别
: 入矩 阵 , ∈ 输 C
对矩阵 A进行特征值分解 :
A = (2 1) 其 中 : ∈C 州是 一 个 对 角矩 阵 , A 由离 散 时 间 系统 的极
Rh 输 出矩 阵 , ∈R : : D 直馈 矩 阵 , 输入 个数 ,: m: z输 、 () 2
点 A组成 。系统参数可 以由状态矩阵 A 输 出矩 阵 c 、
() 5
A =B~ = ( ×T ) T 2
1
(5 1)
可以看出通过式( 1 和式 (5 得到的控制矩阵 A 1) 1)
2 一1 2 一2 … z
具有 相 同 的系 统 特 征值 , 两 种方 法 所 得 结 果 中物 理 即 模态 将会 成 对 出现 。而 对 于 由噪 声 、 阶次 过 估 计 等 因
实现结果的 自动拾取。通过数值仿真 和实例分析验证该方法 的有效性 。 关键词 :参数识别 ; 随机子空间法 ; 系聚类 ; 谱 稳定 图
中 图分 类 号 :T 12 B 2 文 献 标 识 码 :A
Au o a i t c a tc s bs a e i e tfc to fm o a t m tc s o h s i u p c d n i a i n o d l i p r m e e s ba e n hir r h c lc u t rng m e h d a a tr s d o e a c ia l se i t o T ANG opig - Ba - n ’,ZHANG o we ,CHEN u Gu — n Zh o
A l +
J
() 1
由式 (0 可 以看 出 : 1)
A = T = ( × ) T ( 1 1)
环境激励下大型桥梁模态参数识别的一种方法
第一作者 秦世 强 男, 博士生 ,97年生 18
振 动 与 冲 击
21 第 3 0 2年 1卷
梁 自振 频 率 低 、 态 密 集 , 之 测 试 过 程 中 的 噪声 影 模 加
为此 , 以通过 如下 的算 法 实 现 : 可
响, E 在 MD过 程 中往 往 会 出现 模 态 混 叠 ; 纯 的 利 用 单 SI S 方法 识别 模 态参数 又存 在 虚假 模 态 的 问题 , 因此 尝 试 从解 决这 两 个 问题 的角 度 出发 , 引入 屏 蔽 信 号 以消 除E MD中的模 态混 叠 , 将原 始 信 号分 解 成 一 系列 只 包 含 结构某 一 阶 固有振 动信 息 的 I , MF 然后 再 利用 SI S识
w spooe .T em to a b sd o ad it— s i e mp i l o ed cm oio B E D)a d s c at a rp sd h eh dw s ae nb n w d r tc d e ic d eo p s i h ert ra m t n( R M n t h s c o i
模 态参 数 识 别 是 桥 梁 健 康 监 测 的重 要 组 成 部 分 , 准 确 的识 别模 态参数 , 是进 行 有 限元 模 型修 正 、 构 损 结 伤 识别 以及 性 能评 定
通 过 测得 的动 力 响应 , 别 桥 梁 的 固有 振 动 频率 、 识 阻尼 比和 振型 的过 程 。环境 激 励下 的模 态参 数 识 别 由于不 需 要外 部 激励 设 备 且 不 影 响 正 常交 通 , 已经 成 为 大 型
振
第 3 第 2期 1卷
动
与
冲
击
J OURNAL OF VI BRAT ON AND HOC I S K
结构应变模态识别的随机子空间方法
振
动
与
冲
击
第 2 卷第 6 7 期
J OURNAL B OF VI RATI AND S ON H0CK
结构应 变 模 态识 别 的 随机 子 空 间方 法
彭细荣 路新瀛 陈肇 元 , ,
(. 1哈尔滨工业 大学 深圳研究生 院 , 深圳 5 85 ; 105 2清华大学 深圳研究生 院 , . 深圳 585 ; 华大学 土木 工程 系 , 105清 北京 108 ) 004
( )=A ()+B () Y t =C ()( ) t t 。 t , () t 7
其 中 , , R 。 K∈ 为质量 、 尼及 刚度 矩 阵 , 离 阻 n为
对应 的离散状 态空 间描 述 为 :
+ l=A。 +B , Y = C 。 () 8 其 中 = ( a )=( £ ) x kt £ 为离 散 状态 向量 ; , Y
散系统的 自由度数 ( 0 )q t ∈R D F ; () 是位移 向量 , 为
时 间 t 函数 ; () q( ) 的 q f及 t 为速 度 及 加速 度 向量 , 时
间函数上一 点 及 二 点 表示 对 时 间 的一 阶及 二 阶 导 数 ;
为输 入及 输 出采样 数 据 ; 为离散 状态 矩 阵 ; 为离 A B
的应用…。应变模态对结构局部损伤 比位移模态更敏 感 , 但传统的应变测试技术在可靠性及耐久性方面 不能 满足 土木 工 程 结 构 长 期 监 测 的需 求 。近 年 来 , 随 着光纤应变测试技术的发展 , J应用应变模态测试技 术进 行结 构损伤诊 断 的方 法备受 关 注 _ ] 4 。
维普资讯
进行 了验证 。此 方法 可应 用 于土 木 工程 结 构 在 线应 变 模态参 数识别 。
模态参数辨识方法
模态参数辨识方法
一、基于离散时间数据的方法:
1.自相关法:基于自相关函数的方法,通过自相关函数的峰值位置估计模态参数。
2.频率法:通过频率域上的峰值提取方法,估计模态参数。
3.时域法:通过观察结构的动态响应曲线,提取相关的信息计算模态参数。
二、基于连续时间数据的方法:
1.基于有限元模型的方法:通过有限元模型与观测数据拟合,利用最小二乘法估计模态参数。
2.系统辨识方法:利用系统辨识理论,将结构动力学模型视为线性时不变系统,通过观测数据建立结构的状态空间模型,再通过参数辨识算法估计模态参数。
3.压缩感知方法:利用稀疏表示理论,将结构动力学模型表示为稀疏信号,通过压缩感知算法估计模态参数。
在实际应用中,以上方法可以相互结合以提高模态参数辨识的准确性和可靠性。
此外,值得一提的是,模态参数辨识方法的选择也需要根据具体的实验条件和数据特点进行合理的选择。
总之,模态参数辨识方法是结构动力学领域中常用的方法,可以通过使用合适的辨识方法和合理的实验设计,从实验数据中准确地获取结构的模态参数,为结构动力学分析和结构设计提供有力支持。
结构模态参数识别的随机子空间法
~
一 一
一 口
1Y 2
结 构 的参数 识 别是 结构 健康 监测 领 域研 究 的重点 。随机子 空 间方法 是 近年 来发 展起 来 的一 种线 性 系统 辨识 方法 , 该法 不 需要 人 工激 励 , 接从 环 境 激励 的 相应输 出信 号 中提 取结 构 的模 态 参数 , 些 模 态 参数 将 直 这 作 为结 构 健康 监测 、 限元模 型修 正 、 有 受损 结 构物 的评 定 和结 构控 制 的基 础 或输 入 lj I。因此 , 方 法愈 来 愈 . 2 该
法 可用 于有 较 强噪声 情 形 下 的模 态识 别 。那 些 满足 条件 的稳 定极 点被 认 为是 系 统 的真 实极 点 , 在应 用 中 但
需要 改 进,比如 尽 管一 些 极点 也 可 以作 为稳 定 的极 点 , 它 们对 系统 的影 响 很 小 , 以作 为 噪声 模 态 , 现 但 可 用 在 的稳定 图无 法将 其剔 除 。但 通 过考 察其 对 系统 的影响 , 以对稳 定 图方 法进 行改 进 , 可 使其 能 剔除 这 一类 噪 声模态 。
Vo .1 1 9 No 3 .
20 0 6年 9月
( n ier ga dT c n lg ) E gn ei n eh oo n y
S p. 20 e 06
结构模态参数识别 的随机子空间法
常 军 , 张启伟 , 孙 利民
(. 1 苏州科 技学 院 管理科学 与工程系 , 江苏 苏州 25 1 ; . 10 1 2同济大学 土木工程防灾国家重点实验室 , 上海 2 0 9 ) 0 02 摘 要: 参数识别是结构健康监测领域研究 中的重点 随机子空间法是近些 年发展 起来 的一种线性系统辩 识方法,
随机子空间法参数识别与有限元分析
E
㈦ : ) T ( 】 妻
式中 , E是数学期望 ; 是 Koekr rnce 函数 。
12 随机子空间系统参数识别理论 。
在环境 震动测试过程 中, 需要很多个测点才 能量测到结
构的信息 , 然而优化后布置 的传 感器数量 是有限 的, 需要 分
P 鼬 吲 【 lI A Ⅱ ) I EM 1FCA 0N D 珊 I Ⅱ 低 N ' NA 口1 Ⅱ AN E 1A I
BAS 同D ON 【 ) =j sI C S Bs CE 1 (I 【 . I U 】 M匮1 [ 1 oD I
L i Z IHu, HOU Jn , DU Y n ig o g—fn , Z NG B eg HA o
wan ie t c t n me o a xrc tu tr lmo a a a tr ae n te s t p c o e ,a odn i n i a o td C e t ts cu a d lp rmees b s d o t e s a e m d l v i ig d f i i h n a r h a c n e t n lma u l ie t c t n a d i mf e p oe s T e man te r ft i a p o c s it d c d a o v n o a n a d n i a o n t i rc s . h i h oy o h s p ra h i n r u e t i i f i e v o i tUs N YS p o rm u t o po e t d l rmeeso a d l si n i d b t h s f s . i g A S r ga sr s l t rv e mo a a a tr ftef memo e e t e y s c a — r n e h p h r i d f i o
结构模态参数识别-随机子空间法详解
s ⎡=-⎣C()()()()()()()()()11exp exp 1k tc c c k tk t t k t k t d τττ+∆∆+∆=∆∆++∆-⎰x A x A B u 0()()() T T T k kkk t =∆=x x q q&表示由采样开始时刻位移和速度向量组成的状态向量1k +x 表示k +1时刻的状态()1k t ττ'=+∆-令()τu 在一个采样间隔内是常数假设1k k k+=+x Ax Bu ()exp c t =∆A A离散随机状态空间模型kw 1k +x kx kv ky A∆C1i Ti -=R CA G协方差驱动随机子空间识别方法建立的核心表达式。
)协方差驱动随机子空间法(3)随机状态空间模型(t)和对应协方差序列输出信号y150060070080090010002x 10-3时间 (s))channel5 channel6 channel7 channel10 channel11 channel12 channel17 channel19 channel21 channel23 channel26 channel27 channel296.2 随机子空间法2)数据驱动随机子空间法由系统状态空间出发,直接利用实测数据进行分析,无需进行协方差计算。
基于协方差的SSI方法基于数据驱动的SSI方法系统阶次n+1系统阶次n f n+1,ξn+1,Φn+1f n, ξn, Φn|f n+1‐f n|<f的评判标准是否不稳定是否不稳定|Φn+1‐Φn|<Φ的评判标准是否不稳定|ξn+1‐ξn|<ξ的评判标准完全稳定。
模态参数辨识方法——综述
模态参数辨识方法——综述模态参数辨识方法综述摘要:本文对模态分析和模态参数识别进行了综述,对当前识别方法的原理、识别精度及适用条件进行阐述和比较,提出环境激励下模态参数识别方法需解决的关键问题及模态分析在缺陷检测和结构优化中作用。
关键词:模态分析模态参数识别模态分析与缺陷检测结构工作模态0引言模态分析是将线性时不变系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,坐标变换的变换矩阵为振型矩阵,其每列即为各阶振型。
模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。
这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。
振动模态是弹性结构固有的、整体的特性。
如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内,各阶主要模态的特性,就可能预知结构在此频段内,在外部或内部各种振源作用下实际振动响应,而且一旦通过模态分析知道模态参数并给予验证,就可以把这些参数用于(重)设计过程,优化系统动态特性,或者研究把该结构连接到其他结构上时所产生的影响。
模态分析的最终目标是识别出系统的模态参数,为结构系统的振动分析、振动故障诊断和预报、结构动力特性的优化设计提供依据。
解析模态分析可用有限元计算实现,而实验模态分析则是对结构进行可测可控的动力学激励,由激振力和响应的信号求得系统的频响函数矩阵,再在频域或转到时域采用多种识别方法求出模态参数,得到结构固有的动态特性,这些特性包括固有频率、振型和阻尼比等。
有限元法是当前分析机械结构模态的主要方法,很多学者研究了单裂缝和多裂缝缺陷对不同结构动态特性的影响,但这些研究仅局限于出现缺陷结构的当前状态,考虑到缺陷在机械结构使用过程中的扩展,提出了模态分析与缺陷扩展理论相结合的方法分析缺陷的发展趋势,便于机械结构剩余寿命的评估,使已达到设计寿命的结构在失效前仍然发挥其功能,节约了经济成本。
结构模态参数识别的随机子空间法
结构模态参数识别的随机子空间法随机子空间法是一种用于结构模态参数识别的有效方法。
结构模态参数识别是指通过分析结构的振动响应数据,确定结构的固有频率、阻尼比和模态形态等参数,以揭示结构的动力特性和健康状况。
在实际工程中,准确的结构模态参数识别对于结构的安全评估、健康监测和结构优化设计等方面具有重要意义。
随机子空间法是一种基于结构的振动响应数据的子空间分析方法。
它通过将结构的振动响应数据矩阵进行奇异值分解,提取其中的主要特征向量,进而得到结构的模态参数。
与传统的频域法相比,随机子空间法具有以下优势:首先,它不需要事先假设结构的模态数量,能够自动识别结构的模态数量;其次,它可以在较短的时间内完成模态参数识别,适用于长周期结构的分析;此外,随机子空间法对于信号中的噪声和干扰具有较好的抗干扰能力。
随机子空间法的基本步骤如下:1. 收集结构的振动响应数据,包括加速度或位移信号;2. 构建振动响应数据矩阵,将采集到的振动响应数据按照时间序列排列;3. 对振动响应数据矩阵进行奇异值分解,得到奇异值和奇异向量;4. 根据奇异值的大小选择主要特征向量,确定结构的模态数量;5. 根据选择的主要特征向量重构结构的振动响应数据矩阵;6. 对重构的振动响应数据矩阵进行模态参数识别,包括计算固有频率、阻尼比和模态形态等参数。
为了提高随机子空间法的识别精度和稳定性,还可以采取以下措施:1. 在数据采集过程中,应该选择合适的采样频率和采样点数,以充分反映结构的振动特性;2. 在进行奇异值分解时,可以采用截断奇异值的方法,去除奇异值中的噪声和干扰;3. 在选择主要特征向量时,可以采用奇异值比值法或百分比能量法,以确定结构的模态数量;4. 在模态参数识别过程中,可以采用最小二乘法或最大似然法,对模态参数进行估计,提高参数的准确性和稳定性。
随机子空间法在结构模态参数识别中得到了广泛应用。
例如,在桥梁结构的健康监测中,可以通过振动响应数据的采集和分析,实时监测桥梁的振动特性和结构健康状况,及时发现结构的异常变化和损伤;在建筑结构的优化设计中,可以通过对不同结构方案进行振动模态参数识别,评估结构的动力特性,为结构优化设计提供参考。
随机子空间方法结构模态识别的数值仿真
2 简 支 梁 环 境激 励 下 S I S 的数 值 实 现
模 拟 一 单 跨 简 支 梁在 环境 激 励 下 的动 态 响 应 。
() 1 () 2
该 梁 的 几 何 尺 寸 为 : 度 2 T, 面 尺 寸 为 跨 O1 截 I
作者 简 介 : 伟 宜 (9 O )男 , 徽 安 庆 人 , 庆 市 城 乡 规划 设 计 院 工 程 师 章 17 一 , 安 安 1 2 《 程 与建 设 》 2 1 年 第 2 6 工 01 5卷第 2期
方 法 实 现该 梁 的 模态 参 数 识 别 。
关键 词 :S 方 法 ; SI 环境 激 励 ; 白噪 声信 号
中 图分 类 号 : 4 . U4 1 3
文 献标 识码 : A
文章 编 号 :6 35 8 (0 1 0 1 20 1 7 71 2 1 )20 6 —2
0 引
言
阶和第 2阶 阻尼 比均 为 一0 O , . 2 则 一2 z/ + w ∞ (
0 ) 一2 ( ) , 2 c + ) U 。环 境 随机 激 励 信 号作 用 于 梁
节点 上 , 了与 真实 情 况 相 符 , 构 模 态 分 析 与 时 程 为 结 分析 均考虑 了梁 的 自重 恒载 的影 响 。
ftr i e 估计 卫 +在 系统矩 阵 未知 时就 可 由该 状态 前 k l 。 个观察 数据 序列 { Y Y … Y 一 } 定: 1确
将估 计 出的 代 人状态 方程 , 即有
机 子 空 间方法 在 内的 子空 间方 法 的数 学 算 法 , 在 并
前 人 的基 础 上 提 出 了基 于 数 据 的 子 空 间识 别 的 概
( ) () () 一 . 针
基于随机子空间识别的模态参数不确定性量化方法
基于随机子空间识别的模态参数不确定性量化方法罗杰;康杰;孙嘉宝;曾舒洪【期刊名称】《振动与冲击》【年(卷),期】2024(43)8【摘要】基于协方差驱动随机子空间识别(covariance-driven stochastic subspace identification,SSI-COV)方法的模态参数识别具有强鲁棒性、高精度的优势,在结构工作模态分析中应用广泛。
为保证模态参数识别的准确性,新近提出的基于随机子空间(stochastic subspace identification,SSI)的模态参数不确定性量化方法,可有效估计模态参数的方差,但由于其计算各中间变量时,需显式表示出Jacobian矩阵,导致矩阵运算维度高、计算效率低。
为此,提出一种基于SSI-COV的模态参数不确定度高效计算方法。
首先,计算振动响应相关函数的方差,通过奇异值分解(singular value decomposition,SVD),选取合适的奇异值截断阶数,由每阶奇异向量组装出多组Hankel矩阵的扰动。
其次,根据一阶矩阵摄动理论,隐式计算SSI-COV算法各中间变量的一阶扰动,最终,由多组模态参数的扰动叠加计算出方差。
最后,以桁架结构模型为例,采用所提方法辨识结构模态参数并计算模态参数方差,分析了Hankel矩阵维度及相关函数奇异值截断阶数对辨识结果的影响,结果表明计算得到的模态参数方差与蒙特卡洛仿真(Monte Carlo simulation,MCS)结果非常接近,且模态参数不确定度可作为剔除虚假模态的有效依据。
【总页数】8页(P272-279)【作者】罗杰;康杰;孙嘉宝;曾舒洪【作者单位】南京航空航天大学航天学院【正文语种】中文【中图分类】O32【相关文献】1.基于随机子空间结合稳定图的拱桥模态参数识别方法2.基于延时随机子空间方法的非白噪声环境激励结构模态参数识别3.基于随机子空间方法的模态参数识别研究(英文)4.基于GNSS信号的随机子空间模态参数识别方法5.基于随机子空间法的结构模态参数自动识别方法因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
随机状态空间模型:
xk 1 yk
Axk Cxk
wk vk
A:状态矩阵,表示系统内部状态的联系 C:输出矩阵,表示输出与系统内部状态的关系 Wk:环境干扰和建模的不精确而引起的过程噪声 Vk:传感器的不精确或环境对传感器的影响引起的测量噪声
子空间基本原理
正交投影
ΠB BT (BBT ) B A/B AΠB ABT (BBT ) B
Xˆ i
S1/2 1
V1T
同理可得Oi 1
Y f
Y p
Γi1Xˆ i1
卡尔曼滤波估计的状态方程
Xˆ i1 Yi/i
A C
Xˆ
i
wi vi
A C
Xˆ i1 Yi/i
Xˆ
i
特征值分解
A 1, diag(i ),i 1,2,
从统计概念而言,即选取了出现频率较高的数 值,出现频率较高的数值其可信度较大,稳定 性较好。
稳定图
系统阶次n+1 fn+1,ξn+1,Φn+1
系统阶次n fn, ξn, Φn
|fn+1-fn|<f的评判标准 是
|ξn+1-ξn|<ξ的评判标准 是
|Φn+1-Φn|<Φ的评判标准 是
完全稳定
否 不稳定
不同的系统阶次对应一个余弦值 在系统阶次顺序增长的过程中,该余弦值会有一个
大的跳跃 该处跳跃所对应的系统阶次认为是真正系统阶次 由于存在噪声的影响,高于真正系统阶次的余弦值
并不为零,而是跟零比较接近的一个值 研究表明该方法对低阶系统有较好的判断,对高阶
系统效果不佳,奇异值与阶次关系曲线比较平滑, 没有明显的跳跃,难以判定系统阶次
否 不稳定
否 不稳定
试验验证
试验模型全长26.625m 跨径组成
5.1875+16.25+5.1875m 主塔高5.1816m,边墩
高1.7755m,桥面宽 2.9375m 模型除承台采用现浇 混凝土,塔柱和箱梁 均采用钢材。
试验验证
一次性采集
•足够的传感器 •7个采集点 (Case I)
连续与离散的关系
A e Act
ci ,ci
ai
ci
ln i
t
bij
ci ,ci ii i
1 i2 j
i ai2 bi2 (rad / s)
fi
ai2 bi2 / 2 (HZ) i
ai ai2 bi2
振型:
C
e B u( (k1)t Ac ((k1)t )
0
c
)d
xk1 Axk Buk
同理,输出方程即也可以写成
yk Cxk Du k 离散的空间状态模型
xykk1CAxxkkDBuukk
随机状态空间模型
测试必然存在噪声干扰
过程噪声wk 测量噪声vk
离散的空间状态模型可以写成
xk 1 yk
Axk Cxk
Buk Duk
wk vk
随机状态空间模型
uk
B
wk
xk+1
xk
+
△
A
wk
yk
C
+
D
wk
滞后
xk+1
xk
+
△
C
A
vk yk
+
结构振动的状态空间模型
••
•
M q(t) C q(t) Kq(t) F(t) B2u(t)
Q1T QT2
T
R11
1
0 T
R11
R 21
R 22 I R11
0 T R11
0 I R11
0 T 1 R11
0
Q1T QT2
R
RT
21 11
(R11R
T
11
)
1
R
QT
11 1
Y
Oi 1
f
Y
R 41
p
R 42
Q1T
R Q R Q T
43
2
T [l(i1) 1:2li,1:l(i1)] [1:l(i1),:]
Q3T
卡尔曼滤波
SSI方法是通过卡尔曼滤波估计,由输出矩阵得到系 统状态的。
将状态向量������������的卡尔曼滤波估计写成���ො���������
方 程
q(t)
引入状态向量
x(t) • q(t)
y(t) Cc x(t) Dcu(t)
状态方程-离散
实际采样是离散的,连续状态空间模型不 符合实际情况
应将连续状态空间模型转换成离散状态空 间模型。
假设连续状态空间模型存在初始条件:当 时t=t0,x(t)=x(t0) 。那么状态变量 的通解可 以求得:
卡尔曼滤波
Oi
Yf
Yp
Φ Φ Y [Yf,Yp ] [Yp,Yp ] p
CiLi1Yp
Oi
Ci Li1Yp
Γ
i
Δ
c i
Li1Yp
Γi Xˆ i
Oi USVT U1
U2
S1 0
0 0
V1T V2T
U1S1V1T
Γi U1S11/2
环境激励、完好状态
试验环境 采集方法 试验目的
多次采集
•以某些点为参考点 •移动其余传感器至不 同位置 •多次测量到所有测点 数据 •四次采集(Case II) •五次采集(Case III)
卡尔曼滤波
以上几式写成如下形式,即可得到输出数据与状态 向量之间的关系
���ො��������� = Δ������������L−������1
������0 ������1 … ������������−1
������������ = Δ������������L−������ 1
稳定图
稳定图确定系统阶次的方法假定模型有不同的 阶次,得到多个不同阶次的状态空间模型,对 每个模型进行模态参数识别,将得到的所有模 态参数绘制在稳定图上,以确定模态阶次和模 态参数
稳定图包含统计思想,在未知准确系统阶次前 提下,选择一定范围进行多次重复计算,以计 算结果稳定性选取所得模态参数
yi j2 yi j1
1 j
Y0,i1 Yi,2i1
1 j
Yp Yf
yi 1
yi2
yi j
y2i 1
y 2i
y2i j2
y0 y1 y j1
R 22
Q1T QT2
B R11
0
Q1T
Q
T 2
A/B ABT (BBT ) B
R 21
R 22
Q1T
Q
T 2
Q1T
Q
T 2
T
R11
0
T
R11
0
Q1T QT2
稳定图定阶
随机子空间方法模态参数识别过程中,唯 一需要确定的参数为系统阶次。
准确的系统阶次往往是未知的,在识别系 统的过程中估计系统阶次是关键的一环。
系统阶次的估计主要有两种方法:
奇异值的跳跃性 稳定图
奇异值跳跃
根据子空间的原理,矩阵奇异值分解以后所得的奇 异值可认为是将来输出行空间对过去输出行空间投 影过程中的主对角余弦值
R
21R
T 11
(R
T
11
)
1
(R
11)
1
R
QT
11 1
R
QT
21 1
C
B A
R11 R 21
0 Q1T
R 22
Q
T 2
0
Q1T QT2
投影矩阵
l(i 1) R11 0
H 0,2i 1 QR分解
l R 21 R 22
x(t) e Ac (tt0 ) x(t0 )
t
e
0
Ac
(t
)
Bc
u
(
)d
t t0
状态方程-离散
设时间间隔为Δt,离散的时间序列为0、Δt、2Δt,∙∙∙,
(k+1)Δt,∙∙∙,把t=(k+1)Δt,t0=kΔt代入上式可得:
x((k 1)t) e Act x(kt)
yi 1
yi2
yi j
y
2i
1
y 2i
y 2i
j2
同一矩阵的投影
C
B A
R R
11 21
0 R 22
Q1T QT2
,