关于增量谐波平衡法的一个注解.pdfVIP

第5卷增刊 非线性动力学学报 Vol5SI inSciencemld Jun．1998 1998年6月 JounmlofNonlinearDyAl，unics Teclmology 关于增量谐波平衡法的一个注解+ 陈泳斌杜建伟陈树辉 (中山大学应用力学与工程系广州，510275， 摘要本文对采用增量谐波平衡法时如何选取谐波项的项数加以注解 关键词 非线性振动，增量谐波平衡法 中图分类号0322 1引言 HarmonicBalance 增量谐波平衡法(Incremental 树辉121将其应用于求解具有二次、三次非线性系统的联合共振问题．IHB法最基本的思想 波项的项数的问题，初学者对此会感到有些困难，本文旨在对此问题加以注释． 为叙述方便，考虑～个自由度的振动方程 m量+砖+．b+如x2+k3x3=Fcos09t (1) IHB法的第一步就是采用N．R方法，称增量过程．设Xo．oo为振动过程的某一已 知状态，则其邻近状态可表示为 x=Xo+缸，09=COo+Ato (2，3) ’ 把(2、3)-式代入方程(1)，井略去高阶项可得到一个关于△r和A09的增量方程 (4) 其中 (5) 再=Fcos09t一(，打印：x孑+棚ocx；+缸o+七2工；+||}3霹) 万称为修正项，或称为不平衡力．若x0，国。为方程(1)的准确解，则R=o． IHB法的第二步是采用谐波平衡法(IIB法)，称谐波平衡过程。即把‰和△r取为 Fourier级数 Ⅳ x0=ao+∑(口，cosi09t+岛sini09t) (6) Ⅳ Ax=觇+∑(衄cosi09t+她sinicot) (7) i=l ·国家自然科学基金，广东省自然科学基金(960029)资助项目 本文子1998年3月20日收到 非线性动力学学报 J998年 把(6)、(7)二式代八上述得到的关于Ar、Aao的增量方程“)，并应用谐波平衡原理．令 方程两边的谐波项系数相等，晟后可得到一个线性代数方程组，共有zⅣ十1个方程。2N +2个未知数，即△(T．、△6．和△珊．当指定其中某一增量(例如A：o)后．其它增量就可 以唯一求解了．{旨波项^r取多少是一个重要的问题，也是初学者较难掌握的问题Ⅳ愈大． 方程纽阶数愈大，对于个自由度系统．则方程组的阶数将是月(2Ⅳ+l卜1．如果是工程 结构经有限元法离敬后化成大量自由度系统，H是一个很大的数目，而n(2N+l卜l将是 更大的数目了．这时，方程组将占据很大的计算机内存，求解方程组则要花费很多计算 时问．所以Ⅳ不是愈大愈好．但是．若^，取得太小．则会影响解的精确度，对某些问题． 这种误差还是不能忽略的．下