第3章多自由度系统的振动机械工程学院王文瑞博士,副教授gmbitwrw@ustb.edu.cnchap3多自由度系统振动2习题课1建立系统的动力学微分方程牛顿第二运动定律)()(0021213222212121tPtPxxkkkkkkxxmmm1m2k3k1k2x1x2P1(t)P2(t)xMk/2mk/2l拉格朗日方程法2()00Mmxmlkxmlmlxmgl2多自由度线性系统的振动求解两自由度线性系统自由振动求解步骤:Step1:建立运动微分方程:质量矩阵[M]、刚度矩阵[K]Step2:解频率方程,求固有频率12;ppStep3:求振幅比Step4:根据初始条件求12;vv21)2(1)1(1,,,AAStep5:自由振动响应(1)(2)(1)(2)111111122(1)(2)(1)(2)22211112122()sin()sin()()sin()sin()xtxxAptAptxtxxAptApt【例题】:例:图示二自由度系统:k1=k2=k,k3=2k,m1=m,m2=2m。求:在初始条件下系统的响应。102010201.2;0;0xxxx当t=0时初始条件m1m2k3k1k2x1x2解:20[][]302kkmKMkkm系统微分方程令系数行列式为零,得到特征方程:假设方程解1122sin()sin()xAptxApt将(2)代入(1),得到1221112232220000kkkmxxkkkmxx2111122212220kpmkkkpm解:求解特征方程系统固有频率为1kpm252kpm求振幅比22222222112AapcbdpA1221121111AapcbdpA1v11v2节点解:由初始条件求振幅102010201.2;0;0xxxx2(1)22102012102022112(2)2110202110202122()1()0.4()1()0.8xxAxxpxxAxxp121020121020211020211020()arctan2()arctan2pxxxxpxxxx解:(1)(2)(1)(2)111111122(1)(2)(1)(2)22211112122()sin()sin()50.4cos0.8cos2()sin()sin()1510.4cos0.8cos2250.4cos0.4cos2xtxxAptAptkkttmmxtxxAptAptkkttmmkkttmm系统的自由响应为102010201;0xxxx讨论:当初始条件为的系统响应?解:102010201;0xxxx由初始条件(1)21(2)2111(1)0112121(111)00112AA12系统响应(1)(2)(1)(2)111111122(1)(2)(1)(2)22211112122()sin()sin()cos()sin()sin()coskxtxxAptApttmkxtxxAptApttm【课堂练习】:对例所示系统中,已知m1=m2=m,k1=k3=k,k2=4k,求该系统对以下两组初始条件的响应:102010201;0;0xxxx(1)初始条件102010201;1;0xxxx(2)初始条件4.2自由振动的求解第4章两自由度线性系统的振动【内容回顾】─振幅比与主振型20KpMA1221111111210ApmkkA2222111221210ApmkkA•第一阶主振型:质量m1与m2沿相同方向运动;•第二阶主振型:质量m1与m2沿相反方向运动。1v11v2节点4.2自由振动的求解第4章两自由度线性系统的振动【内容回顾】─初始条件的响应2(1)22102012102022112(2)2110201110202122()1()()1()xxAxxpxxAxxp121020121020211020211020()arctan()arctanpxxxxpxxxx─自由振动响应是两种主振动的叠加(1)(2)(1)(2)111111122(1)(2)(1)(2)22211112122()sin()sin()()sin()sin()xtxxAptAptxtxxAptApt•一般初始条件下,两种不同频率主振动的叠加,不一定是简谐振动;•某种特殊条件下,以某一阶固有频率按相应的主振型振动,为简谐振动。•某种特殊条件下,会发生一种称为拍振的特殊振动形态。3拍振现象两个摆长、质量相同的单摆,中间以弹簧相连,形成两自由度系统。得到摆做微小振动的微分方程:221121222221()()mlmglkamlmglkalgp12222mlkalgp求解频率方程:求系统的振幅比:1112第一主振动)sin(11)1()1(1tp)sin(11)1()1(2tp)sin(22)2()2(1tp)sin(22)2()2(2tp第二主振动系统的自由响应振动(1)(2)111(1)(2)1122sin()sin()ptpt(1)(2)222(1)(2)1122sin()sin()ptpt若初始条件(t=0时):10212(0);(0)0;(0)(0)00)2()1(21122得到双摆作自由振动的规律0112(coscos)2ptpt0212(coscos)2ptpt若弹簧的刚度k很小,即222kamllg<<212110coscos22pppptt212120sinsin22pppptt12Δppp拍频率p1、p2相差很少。拍振周期322π2πΔBmTglpka10Δcoscos2aptpt20Δsinsin2aptpt212pppa•拍振是一种比较普遍的现象,凡是由两个频率相近的简谐振动合成的振动,都可能产生拍的现象;•某种特殊初始条件下,使得两个质量的振幅相等或相近的,可能导致拍振的出现•拍振的周期取决于系统参数。一般刚度越小,周期越长;质量越大,周期越长。【小结】1122sinsinxBtxBt4受迫振动求解相应非齐次方程的特解(考虑简谐激励)1112122122212322()sin()sinmxkkxkxFtmxkxkkxFt引入121kkam21kbm22kcm232kkdm111Ffm222Ffm11212122sinsinxaxbxftxcxdxft原方程的改写为设特解为11212122sinsinxaxbxftxcxdxft得到关于振幅B1、B2的非齐次代数方程组将特解带入原方程21212122()()aBbBfcBdBf由此解出强迫振动的振幅2121221222()()()()dfbfBcfafB进而得到系统的稳态响应其中222222212()()()()()adbcpp1122sinsinxBtxBt讨论:1、响应为简谐振动,振动频率为干扰力频率。2、振幅大小不仅和干扰力的幅值大小F1、F2有关,且和干扰力的频率有关。3、当=p1或=p2时,振幅B1、B2将会增大,发生共振,可能出现两次共振,共振时的振型就是相应的主振型。4、实践中经常用共振法测定系统的固有频率,并根据测出的振型判定固有频率的阶次。22122112()()BcfafBdfbf122112121112()()pBcfapfvBdpfbf222122221212()()pBcfapfvBdpfbf2121221222()()()()dfbfBcfafB1122sinsinxBtxBt【例题】:已知:k1=k2=k,k3=2k,m1=m,m2=2m;F1(t)=F1sinωt,F2(t)=0。求:系统的响应。m1m2k3k1k2x1x2F1(t)F2(t)解:1kpm252kpm11212前面求解得到11Ffm20f根据已知带入公式求得振幅221211222222122121222222212()(32)()()()(52)()()()()(52)dfbfkmFBppkmkmcfafkFBppkmkm1122sinsinxBtxBt系统响应1、11122sin0sinsinxBtFxBttk2、22132BkBkm121211132pBkvBkmp2222121322pBkvBkmp当2112211222(32)0()(52)()(52)kmFBkmkmkFFBkmkmk232km反共振思考题已知机器质量m1=90kg,减振器质量m2=2.25kg,机器上有一偏心质量m'=0.5kg,偏心矩e=1cm,机器转速n=1800r/min。试问∶(1)减振器的弹簧刚度k2多大,才能使机器振幅为零?(2)此时减振器的振幅B2为多大?(3)若使减振器的振幅B2不超过2mm,应如何改变减振器的参数?