第四章结构试验模态参数识别.ppt

从（3）式中我们可以看出（ 都是已知了 ）， 是 的线性函数。那么 可将写成矩阵形式，即： （4） （5） （6） 式中： 的元素为： 由向量最小二乘法定理，方程（4）式的最小二乘解为： ③ 先 将在 及 附近展开成台劳级数，并略去高阶项，则有： 由（7）式写成矩阵形式（上式是元素，共有M个（7）式,因 ） (8) (7) 式中: (9) 的元素为 及 组成。（是 个阶矩阵） 而方程（8）的最小二乘解为： （10） 那么我们所求参数为: （11） ④检查所求参数的精度 , ；如不满足，则将 和 作为初始值进行迭代，直到满足精度。有了 和 ，就可求得 即 ，从而由（6）式得 ，经标准化，则可得主模态 。 全部参数识别完成。 （二）克罗斯特曼法 我们先写出N个自由度的粘性阻尼系统频响函数的实部、虚部公式为 （12） （13） 式中 设将要识别的参数为: （14） 那么（12）式与（13）式可写成 （15） （16） 式中系数： （17） （18） 说明：①此方法中固有频率 不是待识别量，是从图解法中得到（即单自由度识别方法得到） ②阻尼比 取上一次迭代值；初始值 由图解识别方法给出。故所以 都可看作已知量。 同前面一样，针对 的采样值 ，可将方程（15）式、（16）式写成矩阵形式，即 （20） （19） 式中： 和 的元素为： 方程（19），（20）式的最小二乘解为 （22） （21） 当解得 和 后，可得各阶阻尼比为 若 的精度不够时，可将 作为初始值代入 中再次求解，直到精度满足为止。 （23） 为求结构的主振型，必须有多个测点， ，这样（19）式与（20）式为相应的增广矩阵方程： （24） （25） 这里 为 阶矩阵。 同理其最小二乘解为 （26） （27） 这样就可将按行标准化，即可得模态矩阵的转换矩阵 最后可由下列公式解得模态参数 所有模态参数都识别了 （28） 我们这里就介绍两种，其它还有很多方法，如： 1、改进的克罗斯特曼方法；2、多项式拟合法；3、利用优化原理识别法，等等。 这些方法都没有跳出最小二乘的方法，基本上是最小二乘法基础上的改进，或具体系统的特殊识别方法。 以上所讨论的是频域模态识别方法，另外还有时域识别方法。这种方法现在也很热门（但太专门化了）但要补上时间序列分析的专门知识，在这里我们就不讨论了，有兴趣的同学可看些专门的书。 试验模态分析 第四章引言 第4章 模态参数的获取方法（频域识别法） 4.1引言 前面已经讨论了模态试验的第一阶段——测得原始数据，并用以推导出所要求的数学模型。也就是说通过介绍模态分析理论，将频响函数与模态参数的关系搞清了，下面我们就可以进行模态参数识别了。 一般情况下工程中遇到的问题都不太可能是简化为一个自由度的系统。其实，实际结构一般都很复杂的。从数学上来讲，其推出的运动方程为一组变量之间相互藕联的方程组。学过振动分析的同学知道：在一定条件下，我们总可以将这一组藕连的方程组解藕，也就是说，可以转换为N个互不相关的单自由度方程来进行研究。这是从理论上或从数学上来说。但实际测试问题是复杂的。从实际测量到的导纳函数数据来看，我们可以将工程中结构的导纳函数分成两大类：一类是模态密集型的，如图所示。另一类是模态稀疏型的，如图所示。 H (a) 模态密集型 H (b) 模态稀疏型 这两种类型，从物理上来说，取决于被测系统的阻尼性质（当然与系统的固有频率也有关）。对小阻尼或比例阻尼的被测系统就是（b）图所示的形式。工程中这类结构是较多的。 从（b）图中我们可以看出，其各阶模态的相互影响较小。故我们由此可以联想到，将各阶模态作为单自由度系统来识别，那就简单多了。也就是，如果我们将单自由度系统的参数识别方法了解了。可将此方法用于多自由度系统的参数识别。 4.2单自由度系统参数频域识别方法和原理 单自由度系