机械振动大作业——简支梁的各情况分析2
- 格式:doc
- 大小:576.50 KB
- 文档页数:18
机械振动大作业姓名:徐强学号:SX1302106专业:航空宇航推进理论与工程能源与动力学院2013年12月简支梁的振动特性分析题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型,计算其固有频率、固有振型。
单、双、三自由度模型要求理论解;十自由度模型要求使用李兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频;连续体要求推导理论解,并通过有限元软件进行数值计算。
解答:一、 单自由度简支梁的振动特性如图1,正方形截面(取5mm ×5mm )的简支梁,跨长为l =1m ,质量m 沿杆长均匀分布,将其简化为单自由度模型,忽略阻尼,则运动微分方程为0=+••kx x m ,固有频率ωn =eqeq m k ,其中k 为等效刚度,eq m 为等效质量。
因此,求出上述两项即可知单自由度简支梁的固有频率。
根据材料力学的结果,由于横向载荷F 作用在简支梁中间位置而引起的变形为)(224348EI F -)(x l x x y -=(20l x ≤≤), 48EI F -3max l y =为最大挠度,则: eq k =δF=348EIl梁本身的最大动能为:)(224348EI F -)(x l xx y -==)(223max43x l l x y -T max =2×dx x y l m l 220)(21⎭⎬⎫⎩⎨⎧•⎰=2max 351721•y m )(如果用eq m 表示简支梁的质量等效到中间位置时的大小,它的最大动能可表示为:T max =2max21•y m eq所以质量为m 的简支梁,等效到中间位置的全部质量为: m m eq 3517=故单自由度简支梁横向振动的固有频率为:ωn =eqeq m k =3171680mlEImk图1 简支梁的单自由度模型二、 双自由度简支梁的振动特性如图2,将简支梁简化为双自由度模型,仍假设在简支梁中间位置作用载荷,根据对称性,等效质量相等,因此只要求出在3/l 处的等效质量即可。
2.1 弹簧下悬挂一物体,弹簧静伸长为δ。
设将物体向下拉,使弹簧有静伸长3δ,然后无初速度地释放,求此后的运动方程。
解:设物体质量为m ,弹簧刚度为k ,则:mg k δ=,即:n ω==取系统静平衡位置为原点0x =,系统运动方程为: δ⎧+=⎪=⎨⎪=⎩&&&00020mx kx x x (参考教材P14)解得:δω=()2cos n x t t2.2 弹簧不受力时长度为65cm ,下端挂上1kg 物体后弹簧长85cm 。
设用手托住物体使弹簧回到原长后无初速度地释放,试求物体的运动方程、振幅、周期及弹簧力的最大值。
解:由题可知:弹簧的静伸长0.850.650.2()m =-=V所以:7(/)n rad s ω=== 取系统的平衡位置为原点,得到:系统的运动微分方程为:20n x x ω+=&& 其中,初始条件:(0)0.2(0)0x x=-⎧⎨=⎩& (参考教材P14) 所以系统的响应为:()0.2cos ()n x t t m ω=- 弹簧力为:()()cos ()k n mg F kx t x t t N ω===-V因此:振幅为0.2m 、周期为2()7s π、弹簧力最大值为1N 。
2.3 重物1m 悬挂在刚度为k 的弹簧上并处于静平衡位置,另一重物2m 从高度为h 处自由落到1m 上而无弹跳,如图所示,求其后的运动。
解:取系统的上下运动x 为坐标,向上为正,静平衡位置为原点0x =,则当m 有x 位移时,系统有: 2121()2T E m m x =+& 212U kx =由()0T d E U +=可知:12()0m m x kx ++=&& 即:12/()n k m m ω=+系统的初始条件为:⎧=⎪⎨=-⎪+⎩&2020122m gx k m x gh m m (能量守恒得:221201()2m gh m m x =+&) 因此系统的响应为:01()cos sin n n x t A t A t ωω=+其中:ω⎧==⎪⎨==-⎪+⎩&200021122n m g A x k x m g ghk A k m m即:ωω=-2()(cos )n n m g x t t t k2.4 一质量为m 、转动惯量为I 的圆柱体作自由纯滚动,圆心受到一弹簧k 约束,如图所示,求系统的固有频率。
单自由度系统自由振动(简支梁)一、 实验目的 1、测定简支梁的等效弹簧常数k ; 2、记录简支梁的自由振动曲线,用分析仪测定系统的有阻尼时的固有频率d ω及相对阻尼系数ζ; 3、用附加质量法测定简支梁的等效质量m ; 4、初步了解振动测试的一些仪器设备及测试方法。
二、 实验装置及原理 1、 实验装置 一根均匀的、截面为矩形的简支梁,其简图如图1所示。
这个系统可看作如图2所示的,有阻尼的单自由度弹簧质量系统,有阻尼时的振动微分方程为: 0=++kx x c xm &&& (1) 令m c n =2,mk n =2ω (2) 则(1)式为:022=++x x n x n ω&&& (3) 再令nn ωζ= (4) 则式(3)为:022=++x x x n n ωςω&&& (5) 其中: m :为简支梁系统的等效质量; k :为简支梁系统对于跨度中点的等效弹簧常数; c :为简支梁下的阻尼常数,n 称为衰减系数,ζ称为相对阻尼系数; n ω:为简支梁系统固有频率,n n f πω2=,d ω为系统的有阻尼固有频率,d d f πω2=。
2、 实验原理 (1) 等效弹簧常数的测定 由于梁在弹性范围内的挠度与梁所受载荷成正比,因此只要在简支梁的跨中点加载,同时图2用百分表读出该点的挠度值,即可测出等效弹簧常数。
(2)记录简支梁系统的自由振动曲线 在简支梁跨度中点贴应变片作用是使梁在振动时的应变量变化转化成电阻量的变化,再将应变片按半桥接法接到动态应变仪上,把电阻量的变化信号放大,并转化成电压量的变化信号,输出到示波器或分析仪,这样即可观察和记录波形。
测试系统框图如图3所示。
(3)附加质量法测等效质量 根据式(2),因为()222n n f m k πω==,21ζωω−=n d ,d d f πω2=要测出简支梁的等效质量m ,只要在原来的简支梁上附加一个已知质量∆,再次求得带有附加质量∆时的固有频率2∆n ω,然后通过下式计算得到m : ()()()()22222222∆∆∆==∆+=n n n n n n f f f f m m ππωω (6) ()()1111222222−∆=−−−∆=∆∆∆d d d d f f f f m &ζζ (7) 三、 实验步骤 1、 测定简支梁系统的等效弹簧常数 在简支梁跨中点处用砝码加载(i=1,2, …., 5),同时用百分表读出该点相对应的挠度值,并记录表1中,按公式算出。
简支梁的变形与振动分析简支梁是一种常见的结构形式,广泛应用于桥梁、楼板等工程中。
在实际工程项目中,我们需要对简支梁的变形和振动进行分析,以确保结构的安全性和稳定性。
本文将从数学模型到应用实例,全面深入地探讨简支梁的变形与振动分析。
一、简支梁的基本理论简支梁是在两端支座的约束下,承受集中力或均布力作用下的一种结构形式。
为了研究其变形和振动特性,我们需要建立数学模型。
1. 简支梁的受力分析在进行简支梁的变形和振动分析前,首先需要了解其受力情况。
在两端支座的约束下,简支梁主要受到弯矩和剪力的作用。
通过弯矩和剪力的分析,可以得出简支梁的受力公式,进而计算结构在承受力作用下的变形。
2. 简支梁的变形分析简支梁在受力作用下会发生一定的变形。
根据梁的假设和力学原理,可以建立简支梁的弹性变形方程。
通过求解弹性变形方程,可以得到简支梁在各个位置的变形情况。
3. 简支梁的振动分析在实际工程中,简支梁还可能受到外力的激励,导致振动现象的发生。
为了分析简支梁的振动特性,我们可以建立简支梁的振动微分方程,并求解得到简支梁的振动模态。
二、简支梁的应用实例1. 桥梁工程简支梁在桥梁工程中得到广泛应用。
为了确保桥梁在运行过程中的安全性和稳定性,需要进行简支梁的变形与振动分析。
通过分析得到的变形和振动数据,可以对桥梁的结构参数进行优化,提高桥梁的工作性能。
2. 建筑结构在楼板、屋顶等建筑结构中,简支梁也扮演着重要的角色。
在设计建筑结构时,需要对简支梁进行变形与振动分析,以确保结构的稳定性和安全性。
通过合理调整支座位置或增加梁的截面尺寸,可以改善简支梁的变形和振动特性。
三、总结简支梁的变形与振动分析对于工程项目的设计和施工至关重要。
通过建立数学模型,进行受力分析和变形分析,可以预测结构在实际工况下的变形情况。
同时,通过振动分析,可以了解简支梁的振动特性,为结构的稳定性提供参考。
在实际工程中,我们还可以利用现代软件进行简支梁的有限元分析,获得更加准确的变形和振动数据。
简支梁模态参数测定之一—测定固有频率与振型一、实验目的1、加深对系统固有频率和主振型的理解;2、掌握振动系统固有频率及主振型的一种测量方法(共振法);3.了解压电式传感器及与它相配的测量系统的工作原理,掌握正确使用的方法;4、了解激振系统的工作原理。
二、实验装置框图图1 表示实验装置系统框图图1 实验装置系统框图三、实验原理试验模态分析法是确定结构固有频率的有效方法,在结构分析中应用广泛,而简支梁也是桥梁结构中一种常见的模型,现代桥梁中依然存在不少采用简支梁模型的桥梁结构。
所以本事通过试验模态法得到简支梁的固有频率和振型,也是桥梁结构分析中一种常用的方法很有实际意义,实验所用的均质等截面简支梁模型,属于小阻尼和连续的无限自由度的振动系统。
本实验模型是一矩形截面简支梁,它是一无限自由度系统。
理论上说,它应有无限个固有频率和主振型,在一般情况下,梁的振动是无穷多个主振型的迭加。
如果给梁施加一个合适大小的激振力,且该力的频率正好等于梁的某阶固有频率,就会产生共振,对应于这一阶固有频率而确定的振动形态叫做这一阶主振型,这时其它各阶振型的影响小得可以忽略不计。
用共振法确定梁的各阶固有频率及振型,具体步骤是首先得找到梁的各阶固有频率,并让激扰力频率等于某阶固有频率,使梁产生共振,然后,测定共振状态下梁上各测点的振动加速度值,从而确定前三阶主振型。
振型:即振动形态,即梁上各个测量点和振幅的关系图。
如图所示为一阶,二阶和三阶的振型图。
在正弦激励下振幅的比值等于加速度的比值。
所以本次试验测量加速度与位置之间的关系就能正确画出振型,大致如图2所示。
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.80.9 10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.80.9 1图2 前三阶振型图根据材料力学理论下简支梁固有频率的计算:2012f l ππ⎛⎫= ⎪⎝⎭E 为材料的弹性模量,查表取E=210Gpa 测量得简支梁b=0.05m h=0.15m l=1m312bh I =s 为梁的横截面积37850kgm ρ=2201135.1622f Hzl l ππππ⎛⎫⎛⎫=== ⎪ ⎪⎝⎭⎝⎭10f f =214140.6f f Hz == 319316.4f f Hz==四、实验方法1、 激振器安装把激振器安装在支架上,将激振器和支架固定在实验台基座上,并保证激振器顶杆对简支梁有一定的预压力(不要超过激振杆上的红线标识),用专用连接线连接激振器和DH1301输出接口。
振动实验报告讲解振动与控制系列实验姓名:李⽅⽴学号:201520000111电⼦科技⼤学机械电⼦⼯程学院实验1 简⽀梁强迫振动幅频特性和阻尼的测量⼀、实验⽬的1、学会测量单⾃由度系统强迫振动的幅频特性曲线。
2、学会根据幅频特性曲线确定系统的固有频率f 0和阻尼⽐。
⼆、实验装置框图图3.1表⽰实验装置的框图图3-1 实验装置框图KCX图3-2 单⾃由度系统⼒学模型三、实验原理单⾃由度系统的⼒学模型如图3-2所⽰。
在正弦激振⼒的作⽤下系统作简谐强迫振动,设激振⼒F 的幅值B 、圆频率ωo(频率f=ω/2π),系统的运动微分⽅程式为:扫频信号源动态分析仪计算机系统及分析软件打印机或绘图仪简⽀梁振动传感器激振器⼒传感器质量块M或 M F x dt dxdt x d M F x dt dx n dtx d FKx dt dx C dtx d M /2/222222222=++=++=++ωξωω(3-1)式中:ω—系统固有圆频率ω =K/Mn ---衰减系数 2n=C/M ξ---相对阻尼系数ξ=n/ωF ——激振⼒ )2sin(sin 0ft B t B F πω== ⽅程①的特解,即强迫振动为:)2sin()sin(0?π?ω-=-=f A A x (3-2)式中:A ——强迫振动振幅--初相位20222024)(/ωωωn M B A +-=(3-3)式(3-3)叫做系统的幅频特性。
将式(3-3)所表⽰的振动幅值与激振频率的关系⽤图形表⽰,称为幅频特性曲线(如图3-3所⽰):3-2 单⾃由度系统⼒学模型 3-3 单⾃由度系统振动的幅频特性曲线图3-3中,Amax 为系统共振时的振幅;f 0为系统固有频率,1f 、2f 为半功率点频率。
振幅为Amax 时的频率叫共振频率f 0。
在有阻尼的情况下,共振频率为:221ξ-=f f a (3-4) 当阻尼较⼩时,0f f a =故以固有频率0f 作为共振频率a f 。
简支梁受力组合变形-概述说明以及解释1.引言概述部分的内容可以如下所示:1.1 概述简支梁是一种常见的结构形式,由于其结构简单、使用方便,广泛应用于建筑、桥梁、机械等领域。
简支梁在受到外力作用时,会发生变形,这种变形对于梁的安全性和使用寿命至关重要。
因此,研究简支梁受力组合变形是提高梁的设计和使用性能的重要方面。
本文将深入探讨简支梁受力组合变形的原因、特点以及对梁结构的影响。
首先,我们将介绍简支梁的定义和特点,包括它的基本结构和建筑原理。
接着,我们将通过对简支梁的受力分析,揭示不同受力组合对梁的变形产生的原因。
随后,我们将对梁的变形进行详细的分析,包括弯曲变形、剪切变形和挠度等。
最后,我们将研究受力组合在简支梁上的影响,探讨其对梁的变形程度和安全性的影响。
通过本文的研究,我们将对简支梁受力组合变形的机理有更深入的了解,同时也能为简支梁的设计和使用提供有用的指导。
这对于提高梁的结构性能、延长梁的使用寿命具有重要意义。
此外,对于简支梁受力组合变形的应用前景,本文也将进行展望,探讨其在未来工程领域中的可能应用和发展方向。
总之,通过本文的研究和分析,我们将为读者提供一个全面的简支梁受力组合变形的概述,从而增进对该领域的理解和应用。
相信本文的内容将对相关领域的研究人员和工程师具有一定的参考价值。
1.2 文章结构文章结构部分的内容可以参考以下示例:2.文章结构本文将按照以下结构进行叙述和分析简支梁受力组合变形的相关内容:2.1 简支梁的定义和特点首先,我们将介绍简支梁的定义和特点。
简支梁是一种常见的结构形式,其特点是两端支座可以自由转动,同时梁自身在受力作用下会发生弯曲变形。
我们将详细探讨简支梁的定义、结构特征以及其在工程实践中的应用。
2.2 受力分析在本节中,我们将进行简支梁的受力分析。
通过分析简支梁在不同荷载作用下的受力情况,我们可以了解到梁的内力分布以及受力大小。
我们将介绍常见的荷载类型,并利用力学原理进行受力计算和分析。
某简支梁桥的抖振响应分析简支梁桥是一种常见的桥梁结构,受到地震、风力、车辆行驶等外部载荷作用时,其会有弹性变形和抖振响应。
抖振又称为疲劳波,是桥梁结构中比较常见的一种现象,本文将从抖振响应的分析角度来探讨某简支梁桥的抖振问题。
简支梁桥的结构形式比较简单,它是由两根简支梁和一个跨中的钢混凝土桥面板构成的。
在受到外部载荷作用时,桥梁会产生弯曲和横向位移,从而导致桥面板上的车辆产生颠簸和震动。
因此,在设计简支梁桥的时候,抖振响应分析是非常重要的。
在进行抖振响应分析时,需要先进行模态分析,即分析桥梁的自振频率、振型和振幅。
简支梁桥的振型可以通过梁的自振方程来求解,简支梁桥的自振频率为:f = (1/2L) * sqrt(EI/m)其中,L为梁的跨度,E为梁的弹性模量,I为梁的截面惯性矩,m为梁的质量。
在确定了桥梁的自振频率和振型后,可以进行抖振响应分析。
抖振响应分析的核心是求出桥梁在外部载荷作用下的响应函数,然后根据响应函数求解桥梁的位移、速度和加速度等参数。
响应函数是指桥梁受到单位冲击力时的响应值,之所以使用单位冲击力,是因为它是一种突然的、瞬间的外部载荷,可以激发出桥梁的最大响应。
在计算抖振响应时,需要考虑桥梁的阻尼效应。
弹性结构的阻尼主要由材料的内耗和结构自身的剪切阻尼来提供,阻尼的大小直接影响着结构的动力响应和稳定性。
构造采用柔性滞回减震器以消除桥梁的抖振响应,从而提高桥梁的抗震性能。
在某简支梁桥的抖振响应分析中,需要考虑到桥梁的自身特性、外部载荷作用、阻尼效应等因素。
通过对桥梁的模态分析和响应函数计算,可以获得桥梁在不同外部载荷下的位移、速度和加速度等参数,从而评估桥梁在抗震、防风、防震等方面的性能,并采取相应的措施对其加固和改造。
机械振动大作业姓名:徐强学号: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 重复上述过程进行第二次迭代。