振动分析矩阵迭代法课件

1、振动分析矩阵迭代法1振动分析的矩阵迭代法振动分析的矩阵迭代法求解结构自由振动的频率和阵型 学号：2014200040 姓名：吴龙旺振动分析矩阵迭代法2一、基本频率和主阵型的计算基本频率和主阵型的计算 在工程中对结构影响较大的体系的最低的几个频率和阵型，它们对结构影响较大。 计算自由振动的频率常常是利用柔度法的公式进行迭代。 由公式 ，可以得到： 用迭代法反复进行迭代，就可以得到体系的最低频率和相应的阵型。0)1(2xmxxm21振动分析矩阵迭代法3迭代收敛性的证明： stodola法和holzer法stodola法是先假定初始振型并不断迭代调整获得一个实际振型的适当近似形式为止，然后从运动方程

2、确定振动频率。holzer法，恰好与之相反，先假定振动频率并迭代调整直到满足边界条件为止，频率和振型就在满足边界条件的过程中确定了。stodola法（迭代分析法）不断迭代就能得到相应的振型（第一振型）下面是其理论证明，即振型收敛性证明。振动分析矩阵迭代法4.*.231)0(3333221)0(2322)0(1211)0()1(1231)0(3333221)0(2322)0(1211)0()0(21)0(121)0()0(33)0(22)0(11)0()0(1yyymff fvyyymfymmvfyyyyviii是这些惯性力产生的挠度)0(11)0(21)()(1)1(1)1(1)0(21)1(

3、)1(21)1(1)()(yyyvyvyyymfvnsnssnnnnnnnnnnn振动分析矩阵迭代法5例题e13-100. 950.1650.2211122245 . 7545 . 71136001d)0(1v)1(1v80. 5400. 010.12733. 010.18000. 1)1(1v)2(1v287. 5320. 0296.11669. 0296.17000. 1)2(1v)3(1v182. 5306. 0121.11653. 0121.17000. 1)3(1v)4(1v)4(1v159. 5303. 0082.11650. 0082.17000. 1)5(1v振动分析矩阵迭代法

4、6(s 1)(4)21111(s)(5)111max()210.77max()vvvv(1)(0)1211(1)(1)11()(22.5,24.5,18)(1,1,1) 3600218(22.5,24.5,18)(22.5,16.5,9)()tttvmvvmv振动分析矩阵迭代法7二、高阶阵型分析高阶阵型分析第二振型以及更高振型的分析可以通过引入淘汰矩阵的方法实现更高一阶收敛，依次就能求得相关振型。如果第一阵型对假设阵型的贡献为零，那么占支配地位的就是第二阵型；类似的如果前两接阵型贡献均为零，那么占支配地位的阵型就是第三阵型；利用阵型的正交性可以从任意的第二阵型中消去第一阵型分量的影响振动分析矩

5、阵迭代法8 1)(1)o(2)o(21)o(2t1)(1)(21t1)(11t1)o(2t1t1)()o(2moooooymyymymmmy消去由阵型的正交性可知：前乘此式两边，得到：用：以任意的二阶阵型为例振动分析矩阵迭代法9逆迭代法-首选的方法 直接迭代法：优点：简单容易实现，且同时由于该方法是基于动力矩阵的柔度形式，所以向最低阶振动形状收敛。缺点：柔度矩阵是满阵，与窄带刚度矩阵相比的运算相比会使计算效率降低。补充：基于刚度矩阵的动力矩阵：由于收敛于最高阶阵型，所以实际应用不大。mkd1 -k1me振动分析矩阵迭代法10逆迭代法基本思路:逆迭代法是利用是利用刚度矩阵窄带性质的首选方法，因为其是使用刚度矩阵的逆矩阵，所以此法向最低阵型收敛。为保证刚度矩阵k的窄带性质，不形成动力矩阵e，而是把假定的位移向量和质量举证结合，然后