振型截断法----振动力学.pptVIP

振动力学 ————振型截断法 主讲人：易之振型叠加法中，需要求出各个阶的固有频率和与之对应的主振型 ，然后分析响应xt。若系统自由度数n很大时， 及 不便于也不可能全部求出。若激励频率主要包含低频成分，可撇去高阶振型及固有频率对响应的贡献，只利用较低的前面若干项 及主振型近似分析系统的响应，这就是工程上常采用的振型截断法。振型叠加法振型截断法 撇去高阶 及 对响应的贡献 利用较低的前面若干项 及 振型截断法 振型位移法 振型加速度法主要知识点：1）以上两类方法的介绍及对比；2）如何进行截断，即阶数s的确定。 若记 由前知，有坐标变换公式：i1,2，……，s假设已求得系统较低的前s阶固有频率i1,2，……，s及 相应的主振型i1,2，……，s，由第4章知系统在第 i 个主坐标的响应为： 1.振型位移法 则有： 撇去高阶振型部分，就可以得到下列近似的系统响应： 其中由于在上式中响应是由主振型及主坐标的位移叠加组成的，因而这种振型截断法称为振型位移法。如果考虑系统的阻尼，并且假定其主阻尼矩阵是对角阵，那么只需要确定前s阶的振型阻尼比 i1,2,…,s，而将高阶的阵型阻尼比is+1,s+2,…,n都假定为零，即有： 这时，第i i1,2，……，s个主坐标的响应式为： 1）考虑阻尼时，系统的响应： 而表示系统的阻尼矩阵的表达式为： 上式可变为： 已知强迫振动的振动方程为： 2.振型加速度法 将式（5.89）代入上式，并结合下式： 可得到系统的响应近似地为：上式右端第一项是伪静态响应，第二项是由前s阶主振型及主坐标的加速度叠加而成的，因而这种方法称为振型加速度法。由于第二项有存在，比较起振型位移法，振型加速度法改善了收敛性，即可用更少的主振型和固有频率求出同样精度的响应。式（5.94）中的 可以用积分号下的微分法算出为： 利用分部积分，上式也可写为备后用：当考虑阻尼时，式（5.93）成为： 将式（5.89）代入上式，近似地得： 结合第四章公式： 故而由式（5.92）及主振型的正交性，上式右端第二项为： 2）考虑阻尼时，系统的响应： 于是系统的响应近似地为： 下面通过例5.8来观察使用振型截断法时如何选取阵型的个数s。 例5.8：如下图所示，四层楼建筑，简化为刚性楼板和弹性支柱。其余四张为不同的振型图。已知：顶层楼板上作用有简谐激振力：；若激振频率分别为：1）；2）； 3）分别用振型位移法和振型加速度法计算顶层楼板的响应。 其中各主振型的归一化是使最大的元素为1。 系统刚度矩阵、质量矩阵、固有频率及振型矩阵已知如下。 解： 由公式求出主质量、主刚度： 已知激振力向量为： 由第4章知：假设 简谐激振力Pt 与响应同频率，即： 其中是激振力幅的常数列向量； 则系统在主坐标下对该激振力的稳态响应幅值为： 故激振力幅为： 又由第4章知，此时主坐标的稳态响应为： （1）当采用阵型位移法时，系统的的响应近似为： 其中顶层楼板的响应为： 因为振型叠加法有n项，下面只截取前4项，将写出； 并指出当所截取振型个数为s1,s2及s3时的响应部分，即： 其中 此时，激振频率分别取： 将上述顶层楼板的响应表示为： 下表列出了不同频率下系数的值 ： 可以看出：当振型个数取s1时，振型位移法得到的响应对三种激振频率的任何一种都存在较大的误差；取s3时，响应在时是相当精确的，但在 时，响应的误差任较大。这是因为接近于 （前），第四阶主坐标的响应在中占重要成分，而振型截断法却没有包括它。 （2）当采用振型加速度法计算响应时，先算出柔度矩阵： 由式（5.94），顶层楼板的响应近似为： 将式a代入上式，得： 为与精确解比较，仍将上式按（c）的形式写为： 将上述顶层楼板的响应表示为： 下表列出了不同频率下系数的值 ： 从上表可以看出：对于的静态载荷，振型加速度法得到精确解，实际上由式d得知，这个精确解是由伪静态响应给出的； 对于的低频情况，振型个数取s1时已经得到相当好的近似解；取s2时，响应的精度相当于振型位移法中取s3时的精度； 而时，出于与振型位移法相同的原因，振型加速度法同样得不到精度较好的解。根据上例可知，在使用振型截断法求系统响应时，必须把分布在激振频率 附近的固有频率 所对应的主振型都包括在内。工程实践当中，当计入激振频率值±20%范围内的固有频率对应的主振型时，一般已能得到较好的近似解。另，有结论：对于低频激振力，振型加速度法求出的响应比振型位移法所得到的更好一些。下面以无阻尼系统为例说明原因： 第4章有公式 ：从而得将上式代入5.94，得到：根据第4章柔度矩阵的模态展开式可知，上式右端第二项圆括号中的部分可以写为：于是式5.100可表示为：称为剩余柔度矩阵，上式右端