东北大学机械学院机电系统与非线性振动控制课程大作业

1、机电系统及其控制过程中的非线性振动课程报告姓名：学号：指导教师：姚红良时间：2016年7月对非线性振动课程的理解和认识10前言11含立方项非线性方程的椭圆函数解12具体算例22. 1线性谐振子32. 2立方非线性振动32. 3杜芬系统的自由振动33解析解与数值解结果的比较43. 1程序设计43. 2数值实验结果5参考文献:56对非线性振动课程的理解和认识0前言立方项的非线性振动是物理学及工程应用中出现较多的一类非线性振动。立方项系数远 小于线性项系数的杜芬方程,作为弱非线性的典型代表,得到了非常广泛的应用。近年来, 有关含立方项的强作线性实际振动系统的研究越來越多，如双弹簧振子的横向振动、新材

2、料 中的纳米机械共振子的振动及悬索的振动等。描述这些振动系统的微分方程中，线牲项常常 小于立方项，甚至仅存在立方项。为了解决强非线性振动在工程设计中的实际应用，出现了诸如能量法、广义谐波函数平 均袪、范式理论方法、同伦摄动法及迭代摄动法等多种强非线性振动系统周期解的近似求解 方法。这些方法原则卜.都町以用来求解强立方非线性振动方程周期解，但只能得到近似结果。 本文将依据弹性力作用下系统机械能守恒的原理，求解出一类含有线性项和立方项的非线性 微分方程的楕确解析解。1含立方项非线性方程的椭圆函数解常见的含立方项非线性口由振动微分方程町表示为：策+比必+=0,式中，kl 0, k2N 0是由振动系统

3、性质决定的非负常数。为了方便，设初始条件为:x(0) - A x (0) = ()可化为:xx + kxx + k2x 兀二(),上式表明：方程表示的系统在振动过程中总机械能守恒。设总机械能为E,则在式中 初始条件h E =斗仏/1 +十尼屮，从而；y-.v + kx +尼.丫 = kiA +他/。匕式表明：系统的相图为闭合凸曲线，则式(1)的解町设为：x = 4cos 3o将上式求导得:代入得:f) = k、+ A：242 (1 + cos%) o0 = A* | 4- -yA*2/l2 (1 +丄2 cosT),则上式可化为:两端枳分得:+ k.A2 , A2J -入Ldt o()Jl -

4、入si0F (弘入)=Lt,Ldt 90其中，F(0,是模为入，参数为o的勒让徳第一类椭圆积分，其反函数为:0 二 am(Ls A)。代入可得:x = Acos am (Ls A) = .4cn (Ag A) rx = 一 L4sn (/, A)(ln (Lz, A) 0其中sn、cn和dn是雅可比椭圆函数。椭圆函数是双周期的亚纯函数。因此，由LT=4K( X)可以得振动系统的周期为：T _ 4K (入)/r02具体算例为了考察方程的解析解式的合理性，具休分析几个算例。2. 1线性谐振子当k2=0时,化为： X + kX = Oo代入可得：x = /lcn ( yTptjO)二 Acos y/

5、kC,卞4/(O)2tt1= =o可以看出结果与众所周知的结果完全相同。2. 2立方非线性振动当ki=o时，化为:X + A*2%3= Oo代入可得:Acn4K （丄匕式表明立方振子的周期与振幅成反比，比例系数由k2确定。这与己有的报道完全相 同。2. 3杜芬系统的自由振动此时方程可以化为：I 9 3 cX + 3&咒 += Uo上式是杜芬系统自由振动微分方程，町采用多种近似方法进行求解。(1 + gA1) 2 ;八一八A =/2（1+ 7）。2（1 +胡）由于（）而.4冇限，囚而一般仃0 胡2 I.拿 A * yl. * y3);hold onplot(y(: , 1), y(: , 2),

6、 4r);xlabel( 4x( m) ), ylabel( v(msAAl)* )3- 2数值实验结果采用以上程序进行数值实验发现：无论kl、k2和A如何取值，由方程得到的振动曲线 和相图与数值解得到的振动曲线和相图完全觅合。图1是kl=ls-2, k2=lm-2s-2, A=lm时方程的数值解（虚线）与解析解（实线）结果的 比牧。图la足振动曲线的比牧，t足时间，x是位移。图lb是相图的比牧，其中，x农示位 移，v表示速度。由图1可见:无论是振动曲线还是相图，解析解与数值解的曲线完全重介。 这充分说明椭圆函数型解式确实有效。（a）援动曲线1.51.0 卞5 ?。-0.5-1.0-1.5（b

7、）相图图1 解析解结果（实线与数佰祥结果（虚线）比较综上所述，本文利用仅受弹性力作用的系统，其乱械能守恒的原理，桔确求解了一类含 线性项和立方项的非线性振动微分方程，并通过与常见的算例及一般情况下的数值解结果进 行对比，说明解析解是有效的.丰富了含立方项非线性振动的研究。参考文献：1 席德勋，席沁.非线性物理学M.南京：南京大学出版社，2007.2 Tomasz K.面向工程的混沌学M.施引,译.北京:国防工业出版社,2008:1-5.3 刘延柱，陈立群.非线性振动M.北京：高等教育出版社，2001:8-109.4 何松林，黄燄，戴祖诚.对称双弹赞振子横向振动的复杂性研究J.昆明学院学报， 2

8、010. 32( 3):86-88.Karabalin R B. Cross M C, Roukes M L. Nonlinear Dynamics and Chaos in Two CoupledNanomechanical Resonators j Physical Review B 2008, 79( 16).6 Che n Shengmi ng, Chen Zili, Luo Yingshe. Non-U near Dynamics A nalysis of inpla ne Motion for Suspended Cable UnderConcentrated Load j. J

9、Cent South UnivTechnol, 2008, 15( si).7_ Yuriy A R Marina V S. New Approach for the Analysis of Damped Vibratio no f Fractionl Oscilltors J. ShockandVibration, 2009( 16).8 丁强生，蒋威，朱小进.一类二阶强迫非线性FDE解的振动性和渐近性J河南科 技大学学报:自然科学版，2011, 32(4):74-78.9 李骊，叶红玲.强非线性系统周期解的能量法M.北京:科学出版社.10 徐兆，詹杰民.强非线性振子的受迫振动J.中山大学学报：自然科学版，1995, 34( 2).11 王永岗，宋慧芳，胥学世.双层旋转扁壳非线性振动分析的同伦摄动法J.清华 人学学报：自然科学版，2007, 47( 8).12 He J H. Some Asymptotic Methods for Strongly Nonlinear Equations J. In