首页机械工程机械设计模态分析理论与试验 03 模态参数估计.pptx

03 模态参数估计.pptx

页数:25页文件大小:96.7 KB来源: 模态分析理论与试验ppt课件

本课件先给出用模态参数表示的频响函数矩阵,说明矩阵维数限制与上、下剩余项的作用;随后按单自由度与多自由度、局部与整体划分估计方法,讲解峰值检测、圆拟合,以及 ITD、LSCE、LSFD 与 ISSPA 的原理和步骤。

第三章模态参数估计封面页,列出含义与基本方程、方法分类、单自由度与多自由度估计及方法比较
模态参数估计含义页,说明利用频响函数或时间历程估计极点、振型与模态参预因子
模态模型基本方程页,给出由模态参数表示的频响函数矩阵及各矩阵含义与极点构成
查看本课件的文字版(共 25 页,供检索与复制)

第1页:第三章 模态参数估计

第三章 模态参数估计

主要研究内容

  • 模态参数估计含义及模态模型基本方程
  • 模态参数估计方法分类简述
  • 单自由度模态参数估计
  • 多自由度时域法模态参数估计
  • 多自由度频域法模态参数估计
  • 模态参数估计方法比较

第2页:模态参数估计的含义

第一节 模态参数估计及模态模型方程 一、模态参数估计:利用测得的频响函数或时间历程来估计模态参数,包括系统极点(阻尼固有频率和阻尼因子)、振型以及模态参预因子等。

第3页:频响函数矩阵的模态表达式

二、模态模型的基本方程 1. 利用模态参数表示的频响函数矩阵:

[H(jω)]=Σ(Q_r{ψ}_r{ψ}_r^T/(jω−λ_r)+Q_r*{ψ}_r*{ψ}_r*^T/(jω−λ_r*)),或写为 [H(jω)]=[V](jω[I]−[Λ])⁻¹[L]。

  • Q_r:模态比例换算因子
  • {ψ}_r:第 r 阶模态向量,[V] 为模态向量矩阵
  • λ_r=σ_r+jω_r:系统极点,σ_r 为阻尼因子,ω_r 为阻尼固有频率
  • [Λ]:由系统极点构成的对角矩阵,[L]=Q[V]^T 为模态参预因子矩阵

第4页:矩阵维数限制与剩余项

2. 实际中各矩阵维数受限制的频响函数矩阵:[H(jω)] 为 N₀×N_i 维,[V] 为 N₀×2N_m,[L] 为 2N_m×N_i;其中 N₀ 为输出点数,N_i 为输入点数,N_m 为可观测模态数。

3. 需要近似考虑所关心频带范围之外模态影响的频响函数矩阵:引入上剩余项矩阵 [U_R](远高于所关心频率范围的模态的影响)与下剩余项矩阵 [L_R](远低于所关心频率范围的模态的影响)。

第5页:时域脉冲响应函数表达式

注意:时域中不能采用剩余项。

4. 时域中用模态参数表达的脉冲响应函数:

[h(t)]=Σ(Q_r{ψ}_r{ψ}_r^T e^(λ_r t)+Q_r*{ψ}_r*{ψ}_r*^T e^(λ_r* t)),或写为 [h(t)]=[V]e^[Λ]t[L]。

第6页:单自由度法与多自由度法

第二节 模态参数估计方法分类简述 一、单自由度法 VS. 多自由度法

假设:一个系统的动态响应是由其若干阶模态振型叠加而成。

单自由度法(SDOF)
假定在给定的频带内只有一个模态是重要的,该模态的参数可以单独确定;只有当系统的各阶模态能够很好解耦时才正确。
多自由度法(MDOF)
需要包含若干个模态的模型对测得的数据进行近似,同时识别这些模态参数。

第7页:局部识别与整体识别的结论

二、局部参数识别 VS. 整体参数识别:由单个脉冲响应函数 h_ij(t)=Σ(ψ_ir L_rj e^(λ_r t)+ψ_ir* L_rj* e^(λ_r* t)) 可得结论:

  • 极点 λ_r 与输出点 i、输入点 j 位置无关
  • 模态振型系数 ψ_ir 与输入点 j 位置无关
  • 模态参预因子 L_rj 与输出点 i 位置无关

第8页:局部参数识别与整体参数识别

局部参数识别(局部估计)
一次使用一个脉冲响应函数估计,对同一参数各次估计的结果可能不相同;最终估计结果必须从一组不同的结果中加以选择,例如计算其平均值,或通过某一指标估计精度。
整体参数识别(整体估计)
利用单个脉冲响应函数的特性对模态参数进行整体估计;可能遭受一些典型的测量误差(如移动加速度计造成的附加质量)侵害,可用数据组(如多个输入)的一致性至关重要。

第9页:单输入与多输入

三、单输入 VS. 多输入

单输入
利用单输入数据进行模态参数估计。
多输入
利用多输入数据进行模态参数估计,能够识别耦合紧密的极点甚至双极点。

第10页:模态模型与直接模型

四、模态模型 VS. 直接模型

模态模型识别法
把输出信号描述为机械结构运动微分方程的特征解(称系统模态)的线性组合,识别问题简化为估计已知模型中的未知参数即模态参数。
直接模型识别法
估计输入—输出微分方程系数,再从这些系数来估计模态参数。

第11页:低阶完整模型与高阶不完整模型

五、低阶完整模型 VS. 高阶不完整模型(N_m:模态数,N₀:响应点数)

低阶完整模型
响应描述为测量数据中全部可观测模态的一个线性组合,N₀≤N_m。
高阶不完整模型
响应点数多于可观测的模态数,N₀>N_m。

第12页:实模态与复模态、时域与频域

六、实模态振型 VS. 复模态振型(根据系统的阻尼形式划分)

  • 实模态振型:比例阻尼或无阻尼系统
  • 复模态振型:非比例阻尼系统

七、时域估计 VS. 频域估计

  • 时域法:直接以输入—输出信号的时间历程为根据的参数估计
  • 频域法:直接以输入—输出信号的频率为根据的参数估计
  • 可由傅立叶变换实现两者变换

第13页:峰值检测法

第三节 单自由度模态参数估计法介绍 一、峰值检测

  • 在固有频率附近,频响函数出现极值,此时其实部为 0,而虚部的幅值最大
  • 出现极值的该点频率就是阻尼固有频率的良好估计
  • 阻尼固有频率处的频响函数关系:|H(jω₁)|=|H(jω₂)|=|H(jω_r)|₂/√2(ω₁<ω_r<ω₂ 为固有频率两侧的频率)
  • 根据半功率点法,可得阻尼比的估计:ζ_r=(ω₂−ω₁)/ω_r

第14页:模态检测法

二、模态检测

  • 根据频域中模态模型,对复模态或实模态向量进行局部估计的一种单自由度方法
  • 根据单个频响函数在阻尼固有频率处的近似值 H_ij(jω_r)≈Q_r ψ_ir ψ_jr/(−σ_r),可得留数 A_ijr 的估计=H_ij(jω_r)×σ_r
  • 利用该方法,需先估计 ω_r

第15页:圆拟合法

三、圆拟合

  • 理论:单自由度系统的速度频响函数(速度对力)在奈奎斯特图(实部对虚部图)上呈现为一个圆
  • 频响函数基本公式:H_ij(jω)=u+jv/(−σ_r+j(ω−ω_r));ω_r 为相位角与圆心的相位角最接近的数据点频率
  • 阻尼比估计:ζ_r=(ω₂−ω₁)/(ω_r[tan(θ₁/2)+tan(θ₂/2)]),ω₁、ω₂ 分据 ω_r 两侧,θ₁、θ₂ 为对应半径与 ω_r 半径之间的夹角

第16页:圆拟合的直径与复留数

4. 圆的直径:复留数 u+jv 与 σ_r 的关系,|u+jv|=√(u²+v²),圆心与固有频率点连线跟虚轴之间的夹角 α 满足 tanα=u/v。

第17页:Ibrahim时域法(ITD)

第四节 多自由度时域模态参数估计法介绍 一、Ibrahim 时域法(ITD)

特点:多自由度、整体估计、多输入、模态模型。

1. 将模态时域方程就采样次数扩展后得到 Hankel 矩阵 [h(t)]_ik,其元素由 [h(t)]、[h(t+Δt)]、[h(t+2Δt)]…按行列时移排列而成,i、k 分别为列向和行向时移计数符号,并有 [h(t)]=[V][L] 的分解形式。

第18页:ITD的特征矩阵与振型求解

2. [h(t)]_ik 移动一个采样时间后:[h(t+Δt)]_ik=[w][h(t)]_ik,式中 [w]=[v]e^[Λ]Δt[v]⁺。

测得 [h(t+Δt)] 和 [h(t)] 后即可估计矩阵 [w];由 [w][v]=[v]e^[Λ]Δt 求得特征值 e^(λ_r Δt),进而求得系统极点 λ_r 和模态振型向量 [v]。

第19页:多参考点最小二乘复指数法

二、多参考点最小二乘复指数法(LSCE)

特点:多自由度、整体估计、多输入、模态模型。

  • 与响应点 O 相对应的脉冲响应矩阵:h(nΔt)=Σ(ψ_or L_r z_r^n+ψ_or* L_r* z_r*^n),其中 z_r=e^(λ_r Δt)
  • 复指数 z_r 与模态参预因子 L_r 的组合满足 p 阶有限差分方程
  • 矩阵 [w]_1,…,[w]_p 的维数为 N_i×N_i,即所考虑的输入点数;为求得 2N 个特征解,阶数 p 需满足 p≥2N/N_i

第20页:LSCE的差分方程与极点求解

脉冲响应是特征解 z_r 或 z_r* 的线性组合,满足:

h(nΔt)[I]+h((n−1)Δt)[w]_1+…+h((n−p)Δt)[w]_p=[0]

  • 由测得的脉冲响应,即可对矩阵系数 [w]_1,…,[w]_p 进行整体最小二乘估计
  • 由矩阵系数可求出 pN_i 个特征值 z_r,进而得到系统极点 λ_r 的估计,以及模态参预因子 L_r

第21页:最小二乘频域法(LSFD)

第五节 多自由度频域模态参数估计法介绍 一、(非线性)最小二乘频域法(LSFD)

特点:多自由度、整体估计、多输入、模态模型。

  • 实测频响函数 H_ij(jω) 包含留数、极点以及上、下剩余项 U_Rij、L_Rij/ω² 的贡献
  • 估计频响函数 G_ij(jω)=G(jω,λ,ψ,L,U_R,L_R),误差 e_ij(jω)=H_ij(jω)−G_ij(jω)
  • 关心频带(ω₀,…,ω_Nf)上的总方差 E_ij=Σ e_ij* e_ij;再对 N_i 个输入点和 N₀ 个输出点求和得到总方差 E

第22页:LSFD参数估计与ISSPA简介

2. 通过使总误差最小来估计模态参数 λ_r、ψ_ir、L_rj、U_Rij、L_Rij,即满足 ∂E/∂u_r=0(r=1,…,N)。

二、结构系统参数识别(ISSPA):特点为多自由度、整体估计、多输入、直接模型。

1. 频域系统方程左乘 [M]⁻¹,写成含 (−ω²+jω[M]⁻¹[C]+[M]⁻¹[K]) 的形式。

第23页:ISSPA的系统方程分解

2. 把 ω_i 集中写成矩阵形式:[−ω²][x]+j[M]⁻¹[C][x][Ω]+[M]⁻¹[K][x]=[M]⁻¹[F]。

令 [x]=[x_R]+j[x_I],可分解为两个分别含实部和虚部的系统方程,从而把未知矩阵 [M]⁻¹[K] 与 [M]⁻¹[C] 分离出来求解。

第24页:ISSPA的矩阵求解与模态参数

3. 按某种最小二乘意义,可解出 [M]⁻¹[K] 和 [M]⁻¹[C];前提是激励频率的数目超过输出数目(即质量矩阵、刚度矩阵和阻尼矩阵的维数)。

4. 根据 [M]⁻¹[K] 可估计模态参数:[M]⁻¹[K]{ψ}_r=ω_r²{ψ}_r。

第25页:模态参数估计法比较

第六节 模态参数估计法比较

  • 时域法适用于处理噪声大的数据
  • 进行宽带分析适于采用时域法,带外模态分析适于频域法
  • 估计大阻尼系统极点宜采用频域法
顶部目录