机械振动大作业-简支梁的各情况分析2

./机械振动大作业姓名：徐强学号：SX1302106专业：航空宇航推进理论与工程能源与动力学院20XX12月简支梁的振动特性分析题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型,计算其固有频率、固有振型。单、双、三自由度模型要求理论解；十自由度模型要求使用兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频；连续体要求推导理论解,并通过有限元软件进行数值计算。解答：单自由度简支梁的振动特性如图1,正方形截面〔取5mm×5mm的简支梁,跨长为=1m,质量m沿杆长均匀分布,将其简化为单自由度模型,忽略阻尼,则运动微分方程为,固有频率ωn=,其中k为等效刚度,为等效质量。因此,求出上述两项即可知单自由度简支梁的固有频率。根据材料力学的结果,由于横向载荷F作用在简支梁中间位置而引起的变形为〔,为最大挠度,则：==梁本身的最大动能为：=Tmax=2×=如果用表示简支梁的质量等效到中间位置时的大小,它的最大动能可表示为：Tmax=所以质量为m的简支梁,等效到中间位置的全部质量为：故单自由度简支梁横向振动的固有频率为：ωn==图1简支梁的单自由度模型双自由度简支梁的振动特性如图2,将简支梁简化为双自由度模型,仍假设在简支梁中间位置作用载荷,根据对称性,等效质量相等,因此只要求出在处的等效质量即可。在至之间积分,利用最大动能进行质量等效,略去小量得：所以,质量矩阵为：双自由度简支梁的柔度矩阵：在b=处作用单位力,挠曲线方程为：则处的变形为：,同理可求：,,其中。所以,柔度矩阵为：动力矩阵：令特征行列式为零,得到频率方程为：其中,,将上式整理得：其中,。解上述方程的根为：,,由式,其中,分别将、代入上式,得第一、二阶主振型分别为：,图2简支梁的双自由度模型三自由度简支梁的振动特性如图3,将简支梁简化为三自由度模型,按照双自由度类似的等效思想,可得等效质量：因此,质量矩阵为：由机械振动中文教材例6.6可知,系统的柔度矩阵为：其中,。动力矩阵：令特征行列式为零,得到频率方程为：其中,,将上式整理得：其中,。利用Matlab软件,求解上述方程的根为：,,,由式,其中,分别将、、代入上式,得第一、二、三阶主振型分别为：,,图3简支梁的三自由度模型十自由度简支梁的数值方法将简支梁简化为十自由度模型〔如图4。图4简支梁的十自由度模型通过在一点施加单位力,计算其余点的挠度,可得柔度矩阵：0.01370.02400.03060.03390.03440.03240.02840.02270.01580.00810.02400.04430.05790.06500.06640.06280.05520.04430.03090.01580.03060.05790.07870.09040.09340.08910.07870.06330.04430.02270.03390.06500.09040.10710.11310.10930.09730.07870.05520.02840.03440.06640.09340.11310.12290.12120.10930.08910.06280.03240.03240.06280.08910.10930.12120.12290.11310.09340.06640.03440.02840.05520.07870.09730.10930.11310.10710.09040.06500.03390.02270.04430.06330.07870.08910.09340.09040.07870.05790.03060.01580.03090.04430.05520.06280.06640.06500.05790.04430.02400.00810.01580.02270.02840.03240.03440.03390.03060.02400.0137表1十自由度挠度变形矩阵十自由度简支梁为十个集中质量的振动模型,每个质量都近似等于,因此,质量矩阵为：动力矩阵为：下面,用如下几种方法计算十自由度简支梁的固有频率与振型。1、邓克莱法利用邓克莱法求基频〔比准确值小：因此,将柔度矩阵主对角线上各元素相加并乘以,可求得：2、瑞利法〔1瑞利第一商柔度矩阵求逆得刚度矩阵：,其中,矩阵见表2。2.0433-1.90030.7778-0.26490.16470.0356-0.28170.2834-0.0446-0.0926-1.90032.9228-2.46751.4375-0.68620.02990.5193-0.62260.3193-0.04460.7778-2.46753.8387-3.34681.6331-0.1417-0.66840.8588-0.62260.2834-0.26491.4375-3.34684.2339-2.91230.7780.4309-0.66840.5193-0.28170.1647-0.68621.6331-2.91233.4193-2.29320.778-0.14170.02990.03560.03560.0299-0.14170.778-2.29323.4193-2.91231.6331-0.68620.1647-0.28170.5193-0.66840.43090.778-2.91234.2339-3.34681.4375-0.26490.2834-0.62260.8588-0.6684-0.14171.6331-3.34683.8387-2.46750.7778-0.04460.3193-0.62260.51930.0299-0.68621.4375-2.46752.9228-1.9003-0.0926-0.04460.2834-0.28170.03560.1647-0.26490.7778-1.90032.0433表2矩阵各元素假设力作用在简支梁中间位置而得到各点的静变形,可以表示为：其中,。因此,可以假设振型：则由瑞利第一商公式：,可得：〔2瑞利第二商同样假设力作用在简支梁中间位置,由瑞利第二商公式：可得:瑞利法中,代表质量矩阵,代表刚度矩阵,代表柔度矩阵,为模态向量。3、兹法将十自由度简支梁缩减为三自由度,假设振型为：则可求出：72.960-0.6023.061.2-0.61.2180.127900.135404.3927-1.23220.1354-1.23224.2842由式,得：,,其中,,因此可得：以及：所以系统的前三阶主振型的近似为：-1.04×-1.39×-0.28×4、矩阵迭代法单位力作用在简支梁中间位置得到各点的挠度变形,将首项化一,得：其中,。因此可假设振型：利用矩阵迭代求第一阶固有频率和主振型：由上,仅3次矩阵迭代后,与基本相等,因此可以认为系统的第一阶主振型为：第一阶固有频率为：5、雅可比法根据雅可比法原理,依次找出上三角非对角线上〔考虑对称性的最大元素,利用公式得到值,代入旋转矩阵,可得：0.349-0.056-0.4170.2150.396-0.394-0.2170.4150.058-0.3470.3960.323-0.134-0.420-0.2120.2120.4200.134-0.324-0.3950.398-0.369-0.0320.334-0.185-0.2150.3400.002-0.4360.4550.3210.4240.230-0.122-0.387-0.387-0.1220.2310.4230.3220.161-0.3240.419-0.3460.0870.244-0.4600.454-0.2820.1020.2320.3880.4220.3220.120-0.120-0.322-0.422-0.388-0.2320.104-0.2310.354-0.4420.466-0.4240.351-0.2630.155-0.0530.1200.2310.3220.3880.4220.4220.3880.3220.2310.1200.4390.101-0.383-0.2190.3040.337-0.236-0.3910.1270.415-0.4040.453-0.158-0.1840.324-0.2420.0440.217-0.4480.399则：0.000500000000000.002700000000000.000300000000000.008400000000000.000100000000000.042400000000000.000100000000000.677500000000000.001100000000000.0002由此得到系统的固有频率：对应各特征值的特征向量即为振型,即的列向量为各阶振型。6、子空间迭代法取假设振型：由动力矩阵迭代得到：0.6941-0.0460.0341-0.01761.3333-0.07680.0706-0.02871.8642-0.08340.108-0.03122.2447-0.06290.138-0.03042.4421-0.02340.1537-0.03182.44210.02340.1506-0.03612.24470.06290.1282-0.03881.86420.08340.0932-0.03471.33330.07680.0568-0.02370.69410.0460.0258-0.011将各列分别归一化得：求得和分别为：再由兹法得特征值问题为：解出：,,,其中,,相应的主振型为：所以：0.2850.0527-0.0249-0.53970.54770.0679-0.0244-0.90330.76610.03320.0048-0.9850.9229-0.01980.0311-0.74631.0042-0.05890.0222-0.27851.0042-0.0607-0.01470.2790.9229-0.0196-0.03670.74770.76610.037-0.01460.98550.54770.06650.02770.9010.2850.04950.03680.5365各列分别归一化后,得：0.2838080.776141-0.67663-0.547640.5454091-0.66304-0.916590.7628960.4889540.130435-0.999490.91904-0.291610.845109-0.757281-0.867450.603261-0.28261-0.89396-0.399460.2831050.91904-0.28866-0.997280.7587010.7628960.544919-0.3967410.5454090.9793810.7527170.9142570.2838080.72901310.544394重复上述过程进行第二次迭代。由：0.19260.0065-0.0021-0.02320.36990.0085-0.0017-0.03880.51710.00480.0008-0.04220.6225-0.0020.0025-0.03220.6771-0.0070.0015-0.0120.6771-0.007-0.0010.01210.6225-0.0019-0.00240.03220.51710.0049-0.00070.04220.36990.00860.00220.03880.19260.00640.00260.0232归一化后得：0.2844480.755814-0.80769-0.549760.54630.988372-0.65385-0.919430.7636980.558140.307692-10.919362-0.232560.961538-0.763031-0.813950.576923-0.284361-0.81395-0.384620.286730.919362-0.22093-0.923080.7630330.7636980.569767-0.2692310.546310.8461540.9194310.2844480.74418610.549763则有：5.61560.32950.41680.00240.32955.1660.17580.02290.41680.17585.22040.08370.00240.02290.08375.62270.00830.00090.000600.00090.61070.01570.00040.00060.01571.91860.002300.00040.00230.1328由：解出：,,,其中,,相应的主振型为：0.284-0.7390.8411-0.54950.5457-0.95780.7085-0.91930.7633-0.5176-0.2488-1.00020.91950.2806-0.9104-0.76371.00050.8685-0.5323-0.28511.00050.87260.43450.28620.91950.27890.99070.76280.7633-0.51970.35560.99990.5457-0.9641-0.75970.91930.284-0.7258-0.94530.5496各列分别归一化后,得：0.2838580.7665180.8489960.549390.5454270.9934650.7151510.9191160.7629190.536874-0.2511410.91904-0.29105-0.918950.7635471-0.90084-0.53730.2850431-0.905090.438579-0.286140.91904-0.289291-0.762650.7629190.5390520.358938-0.99970.5454271-0.76683-0.919120.2838580.752826-0.95417-0.54949进行第三次迭代：0.19260.00650.00240.02330.370.00860.0020.03890.51710.0049-0.00080.04240.6225-0.002-0.00260.03230.6772-0.0072-0.00140.01210.6772-0.00720.0013-0.01210.6225-0.0020.0026-0.03230.51710.00490.0009-0.04240.370.0086-0.002-0.03890.19260.0065-0.0025-0.0233归一化得：0.2844480.755814-0.80769-0.549760.54630.988372-0.65385-0.919430.7636980.558140.307692-10.919