振型截断法振动力学演示文稿.ppt

振型截断法振动力学演示文稿
第一页，共二十页。
振型叠加法中，需要求出各个阶的固有频率 和与之对应的主振型 ，然后分析响应x(t)。
若系统自由度数n很大时， 及 不便于也率求出同样精度的响应。式（）中的 可以用积分号下的微分法算出为：
第六页，共二十页。
利用分部积分，上式也可写为(备后用)：
当考虑阻尼时，式（）成为：
将式（）代入上式，近似地得：
结合第四章公式：
故而由式（）及主振型的正交性，上式右端第二项为：
2）考虑阻尼时，系统的响应：
第七页，共二十页。
于是系统的响应近似地为：
第八页，共二十页。
。
：如下图所示，四层楼建筑，简化为刚性楼板和弹性支柱。其余四张为不同的振型图。
已知：顶层楼板上作用有简谐激振力： ；
若激振频率分别为：
1） ； 2） ； 3）
分别用振型位移法和振型加速度法计算顶层楼板的响应 。
第九页，共二十页。
其中各主振型的归一化是使最大的元素为1。
系统刚度矩阵、质量矩阵、固有频率及振型矩阵已知如下。
第十页，共二十页。
解： 由公式 求出主质量、主刚度：
已知激振力向量为：
由第4章知：假设 简谐激振力P(t) 与响应同频率，即：
其中 是激振力幅的常数列向量；
则系统在主坐标下对该激振力的稳态响应幅值为：
故激振力幅为：
第十一页，共二十页。
又由第4章知，此时主坐标的稳态响应为：
（1）当采用阵型位移法时，系统的的响应近似为：
其中顶层楼板的响应为：
因为振型叠加法有n项，下面只截取前4项，将 写出；
并指出当所截取振型个数为s=1,s=2及s=3时的响应部分，即：
其中
第十二页，共二十页。
此时，激振频率分别取：
将上述顶层楼板的响应表示为：
下表列出了不同频率下系数 的值 ：
第十三页，共二十页。
可以看出：
当振型个数取s=1时，振型位移法得到的响应对三种激振频率的任何一种都存在较大的误差；
取s=3时，响应在 时是相当精确的，但在
时，响应的误差任较大。这是因为 接近于 （前），第四阶主坐标的响应在 中占重要成分，而振型截断法却没有包括它。
（2）当采用振型加速度法计算响应时，先算出柔度矩阵：
第十四页，共二十页。
由式（），顶层楼板的响应近似为：
将式(a)代入上式，得：
为与精确解比较，仍将上式按（c）的形式写为：
第十五页，共二十页。
将上述顶层楼板的响应表示为：
下表列出了不同频率下系数 的值 ：
从上表可以看出：
对于 的静态载荷，振型加速度法得到精确解，实际上由式(d)得知，这个精确解是由伪静态响应给出的；
对于 的低频情况，振型个数取s=1时已经得到相当好的近似解；取s=2时，响应的精度相当于振型位移法中取s=3时的精度；
而 时，出于与振型位移法相同的原因，振型加速度法同样得不到精度较好的解。
第十六页，共二十页。
根据上例可知，在使用振型截断法求系统响应时，必须把分布在激振频率 附近的固有频率 所对应的主振型都包括在内。
工程实践当中，当计入激振频率值±20%范围内的固有频率对应的主振型时，一般已能得到较好的近似解。
另，有结论：对于低频激振力，振型加速度法求出的响应比振型位移法所得到的更好一些。
下面以无阻尼系统为例说明原因：
第4章有公式 ：
从而得
将上式代入()，得到：
第十七页，共二十页。
根据第4章柔度矩阵的模态展开式可知，上式右端第二项圆括号中的部分可以写为：
于是式()可表示为：
称为剩余柔度矩阵，上式右端第二项正是振型加速度法比振型位移法多出的部分。
先考虑振型位移法中撇去的高阶振型部分：
第十八页，共二十页。
对于低频激振力，当 (i=s+1,s+2,…，n)时，上式近似为
正因为振型加速度法中增添了对高阶振型部分的近似项 ，因而求出的响应比振型位移法求出的要好。
振型截断法中包括的主振型个数不仅与激振频