机械振动大作业姓名：徐强学号：SX1302106专业：航空宇航推进理论与工程能源与动力学院2013年12月简支梁的振动特性分析题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型，计算其固有频率、固有振型。
单、双、三自由度模型要求理论解；十自由度模型要求使用兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频；连续体要求推导理论解，并通过有限元软件进行数值计算。
解答：一、单自由度简支梁的振动特性如图1，正方形截面（取5mm×5mm，质量m沿杆长均匀分布，将其简化为单自由度模型，忽略阻尼，则运动固有频率ωn其中k等效质量。
因此，求出上述两项即可知单自由度简支梁的固有频率。
根据材料力学的结果，由于横向载荷F作用在简支梁中间位置而度，则：梁本身的最大动能为：T max=2大动能可表示为：T max所以质量为m的简支梁，等效到中间位置的全部质量为： 故单自由度简支梁横向振动的固有频率为：ωn图1 简支梁的单自由度模型二、 双自由度简支梁的振动特性如图2，将简支梁简化为双自由度模型，仍假设在简支梁中间位置作用载荷，根据对称性，等效质量相等,去小量得：所以，质量矩阵为：双自由度简支梁的柔度矩阵：在所以，柔度矩阵为：动力矩阵：令特征行列式为零，得到频率方程为：解上述方程的根为：第一、二阶主振型分别为：图2 简支梁的双自由度模型三、 三自由度简支梁的振动特性如图3，将简支梁简化为三自由度模型，按照双自由度类似的等效思想，可得等效质量：因此，质量矩阵为：由机械振动中文教材例6.6可知，系统的柔度矩阵为：动力矩阵：令特征行列式为零，得到频率方程为：利用Matlab 软件,求解上述方程的根为：第一、二、三阶主振型分别为：图3 简支梁的三自由度模型四、十自由度简支梁的数值方法将简支梁简化为十自由度模型（如图4）。
图4 简支梁的十自由度模型通过在一点施加单位力，计算其余点的挠度，可得柔度矩阵：0.0137 0.0240 0.0306 0.0339 0.0344 0.0324 0.0284 0.0227 0.0158 0.0081 0.0240 0.0443 0.0579 0.0650 0.0664 0.0628 0.0552 0.0443 0.0309 0.0158 0.0306 0.0579 0.0787 0.0904 0.0934 0.0891 0.0787 0.0633 0.0443 0.0227 0.0339 0.0650 0.0904 0.1071 0.1131 0.1093 0.0973 0.0787 0.0552 0.0284 0.0344 0.0664 0.0934 0.1131 0.1229 0.1212 0.1093 0.0891 0.0628 0.0324 0.0324 0.0628 0.0891 0.1093 0.1212 0.1229 0.1131 0.0934 0.0664 0.0344 0.0284 0.0552 0.0787 0.0973 0.1093 0.1131 0.1071 0.0904 0.0650 0.0339 0.0227 0.0443 0.0633 0.0787 0.0891 0.0934 0.0904 0.0787 0.0579 0.0306 0.0158 0.0309 0.0443 0.0552 0.0628 0.0664 0.0650 0.0579 0.0443 0.0240 0.0081 0.0158 0.0227 0.0284 0.0324 0.0344 0.0339 0.0306 0.0240 0.0137表1十自由度简支梁为十个集中质量的振动模型，每个质量都近似等动力矩阵为：下面，用如下几种方法计算十自由度简支梁的固有频率与振型。
1、邓克莱法利用邓克莱法求基频（比准确值小）：2、瑞利法（1）瑞利第一商柔度矩阵求逆得刚度矩阵：2。
2.0433 -1.9003 0.7778 -0.2649 0.1647 0.0356 -0.2817 0.2834 -0.0446 -0.0926-1.9003 2.9228 -2.4675 1.4375 -0.6862 0.0299 0.5193 -0.6226 0.3193 -0.0446 0.7778 -2.4675 3.8387 -3.3468 1.6331 -0.1417 -0.6684 0.8588 -0.6226 0.2834-0.2649 1.4375 -3.3468 4.2339 -2.9123 0.778 0.4309 -0.6684 0.5193 -0.2817 0.1647 -0.6862 1.6331 -2.9123 3.4193 -2.2932 0.778 -0.1417 0.0299 0.0356 0.0356 0.0299 -0.1417 0.778 -2.2932 3.4193 -2.9123 1.6331 -0.6862 0.1647 -0.2817 0.5193 -0.6684 0.4309 0.778 -2.9123 4.2339 -3.3468 1.4375 -0.2649 0.2834 -0.6226 0.8588 -0.6684 -0.1417 1.6331 -3.3468 3.8387 -2.4675 0.7778 -0.0446 0.3193 -0.6226 0.5193 0.0299 -0.6862 1.4375 -2.4675 2.9228 -1.9003 -0.0926 -0.0446 0.2834 -0.2817 0.0356 0.1647-0.2649 0.7778 -1.9003 2.0433表2假设力作用在简支梁中间位置而得到各点的静变形，可以表示为：因此，可以假设振型：（2）瑞利第二商同样假设力作用在简支梁中间位置，由瑞利第二商公式：可得3、兹法将十自由度简支梁缩减为三自由度，假设振型为：则可求出：由式，得：以及：所以系统的前三阶主振型的近似为：4、矩阵迭代法单位力作用在简支梁中间位置得到各点的挠度变形，将首项化72.96 0 -0.6 0 23.06 1.2 -0.6 1.2 180.1279 0 0.13540 4.3927 -1.2322 0.1354 -1.2322 4.2842一，得：因此可假设振型：利用矩阵迭代求第一阶固有频率和主振型：由上，仅3的第一阶主振型为：第一阶固有频率为：5、雅可比法根据雅可比法原理，依次找出上三角非对角线上（考虑对称性）的最大元素，代入旋转矩阵，可得：0.349 -0.056 -0.417 0.215 0.396 -0.394 -0.217 0.415 0.058 -0.3470.396 0.323 -0.134 -0.420 -0.212 0.212 0.420 0.134 -0.324 -0.395 0.398 -0.369 -0.032 0.334 -0.185 -0.215 0.340 0.002 -0.436 0.455 0.321 0.424 0.230 -0.122 -0.387 -0.387 -0.122 0.231 0.423 0.3220.161 -0.324 0.419 -0.346 0.087 0.244 -0.460 0.454 -0.282 0.102 0.232 0.388 0.422 0.322 0.120 -0.120 -0.322 -0.422 -0.388 -0.2320.104 -0.231 0.354 -0.442 0.466 -0.424 0.351 -0.263 0.155 -0.053 0.120 0.231 0.322 0.388 0.422 0.422 0.388 0.322 0.231 0.120 0.439 0.101 -0.383 -0.219 0.304 0.337 -0.236 -0.391 0.127 0.415 -0.404 0.453 -0.158 -0.184 0.324 -0.242 0.044 0.217 -0.448 0.399则：0.0005 0 0 0 0 0 0 0 0 00 0.0027 0 0 0 0 0 0 0 0 0 0 0.0003 0 0 0 0 0 0 0 0 0 0 0.0084 0 0 0 0 0 00 0 0 0 0.0001 0 0 0 0 00 0 0 0 0 0.0424 0 0 0 0 0 0 0 0 0 0 0.0001 0 0 0 0 0 0 0 0 0 0 0.6775 0 0 0 0 0 0 0 0 0 0 0.0011 0 0 0 0 0 0 0 0 0 0 0.0002由此得到系统的固有频率：6、子空间迭代法取假设振型：由动力矩阵迭代得到：0.6941 -0.046 0.0341 -0.01761.3333 -0.0768 0.0706 -0.02871.8642 -0.0834 0.108 -0.03122.2447 -0.0629 0.138 -0.03042.4421 -0.0234 0.1537 -0.03182.4421 0.0234 0.1506 -0.03612.2447 0.0629 0.1282 -0.03881.8642 0.0834 0.0932 -0.03471.3333 0.0768 0.0568 -0.02370.6941 0.046 0.0258 -0.011 将各列分别归一化得：再由兹法得特征值问题为：解出：所以：0.2850.0527 -0.0249 -0.5397 0.5477 0.0679 -0.0244 -0.9033 0.7661 0.0332 0.0048 -0.985 0.9229 -0.0198 0.0311 -0.7463 1.0042 -0.0589 0.0222 -0.2785 1.0042 -0.0607 -0.0147 0.279 0.9229 -0.0196 -0.0367 0.7477 0.7661 0.037 -0.0146 0.9855 0.5477 0.0665 0.0277 0.901 0.2850.04950.0368 0.5365各列分别归一化后，得：0.283808 0.776141-0.67663 -0.54764 0.545409 1-0.66304 -0.91659 0.762896 0.488954 0.130435 -0.99949 0.91904-0.29161 0.845109 -0.75728 1 -0.86745 0.603261 -0.2826 1 -0.89396 -0.39946 0.283105 0.91904 -0.28866 -0.99728 0.7587010.7628960.544919-0.3967410.545409 0.979381 0.752717 0.9142570.283808 0.729013 1 0.544394 重复上述过程进行第二次迭代。