计算n阶固有频率的一个简便公式.docVIP

第23卷　第3期　2003年9月 西安科技学院学报 Vol.23　No.3Sep.2003 　JOURNALOFXI’ANUNIVERSITYOFSCIENCEANDTECHNOLOGY 　　文章编号:1671-1912(2003)03-0354-03 计算n阶固有频率的一个简便公式 郭志勇 (西安科技大学基础课部,陕西西安　710054) Ξ 摘　要:对多自由度无阻线性正定系统固有频率的特征方程利用根与系数的关系,给出了一个计算多自由度系统n阶固有频率的简便公式,并结合实例说明该式精度较高,具有一定的实用价值。关键词:振动系统;固有频率;计算方法中图分类号:O321　　　文献标识码:A 动系统,(Rayleigh)法、邓克莱(Dunkerley)法、李兹(Ritz)法、子空间迭代法、[1],,但瑞雷法需要准确的估计主振型列阵A,]],导出了一个计算多自由度系统基频的近似公式。 本文将文献[3],给出了一个计算多自由度系统各阶固有频率的近似公式。并用实例说明了应用本公式其计算过程较瑞雷法简便,精度与瑞雷法相当。因此,本文所推出的公式在理论和实际中均有一定的实用价值。 1　振动系统的基本方程 设有n个自由度的无阻尼线性正定振动系统,其振动微分方程为 MX¨+KX=0 (1) 若以各质量的静平衡位置作为广义坐标的原点,则质量阵M为一对角阵,刚度阵K为一对称阵,X,X¨均为列阵。 设其振型列阵为A,则振型方程为 -MpA+KA=0 2 (2)(3) 式中p为系统的固有频率。频率方程为 | K-Mp|=0 2 2　计算固有频率的近似公式 将(3)式展开,则有 n p 2n - i=1 2(n-1) p+mi n i=1j=2ij kmi 2(n-2) -p+…=0 mimjmimj 2 (4) 22 设(4)式有n个不相等的实根,由小到大其依次为　p2则由根与系数的关系有1,p2,…,pn。 n 2 p12p1 + 2p2 2p2 +…+pn= 2p1 2 2 i=1 (5) n +…+pn+…+pn-1pn= 22 i=1j=2ij -mimjm imj 2 (6) Ξ收稿日期:2002-06-06 作者简介:郭志勇(1958-),男,陕西扶风人,副教授,主要从事工程力学教学与研究. ? 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. 第3期　　　　　　　　　　　　　　郭志勇　计算n阶固有频率的一个简便公式(5)式两端平方后减去(6)式两端的2倍有 n 355 4 p1 + 4p2 +…+pn=( 4 i=1 )-2mi n i=1j=2ijn -mimjmimj n = i=1 ∑ ()2+2mi i=1j=2ij ∑ ()mimj (7) 由于 p1p2…pn n 因而可在(7)式左端依次略去前n-1项小量,得到 pn= 4 i=1 ∑ ()2+2mi n i=1j=2ij ∑ ()mimj (8) 为了能一次求得最高两阶固有频率,保留(5)式和(7)式等号左边的最大两项,写为 n pn-1+pn= 4 4 22 i=1 m ni=1 i n (9) pn-1+pn=( (9)式和(10)式联立消去p4n-1,得到 n k2 )-2mi i=1j=2ij kkkmiji j (10) pn- 42 imi n i=1j=2ij -mijmimj =0(11) 解(11)式得到 2n　 pn -1 =2 n i=1 ±mi n i=1 )-4mi n i=1j =2ij -mimjmimj (12) (12)式即为计算最高二阶固有频率的近似公式,(12)式根号前取正号为n阶,根号前取负号为n-1阶。当 当n=3时,由(12)式可计算出二、三阶固有n=2时,(12)式即为二自由度振动系统固有频率的精确计算公式。 频率,其近似程度很高。随着系统自由度数的增大,其更高阶固有频率的误差也随之增大,但始终与瑞雷法求基频的精度相当,而比其它近似方法的精确度高。 3　应用实例 如图1所示,简支梁的质量不计,沿梁长均匀地作用3个集中质量,梁的跨度为l,抗弯刚度为EJ。求梁作横向振动时的各阶固有频率[1]。 解:系统的质量阵为 m M= 02m 系统的柔度阵为 δ= 系统的刚度阵为 -1 K=δ= 3 9117 111611 711923 -22 768EJ 图1　梁的横向振动 Fig.1　Transversevibrationofbeam -22 28l3 9 92232 -22 有关各量计算如下 3 　i=1 α=3(23+32+23)=mi28ml28ml3 ? 1994-2009 C