用谐波能量平衡法求解单摆方程-动力学与控制学报.PDF

第１４卷第３期２０１６年６月 动力 学 与控 制 学报 Ｖｏｌ．１４Ｎｏ．３ １６７２６５５３／２０１６／１４ ／１９７８ Ｊｕｎ．２０１６ ⑶ ＪＯＵＲＮＡＬＯＦＤＹＮＡＭＩＣＳＡＮＤＣＯＮＴＲＯＬ 用谐波能量平衡法求解单摆方程 １ ２ １ １ 李银山 　李彤　韦炳威　李欣业 （１．河北工业大学机械工程学院力学系，天津　３００１３０） （２．华东理工大学承压系统与安全教育部重点实验室，上海　２００２３７） 摘要　应用谐波能量平衡法求解了强非线性单摆方程，谐波能量平衡法与经典的摄动法和谐波平衡法不 同，不是把微分方程和初始条件分离处理；而是把微分方程和初始条件同时处理．用谐波平衡，将描述动力 系统的二阶常微分方程，化为以角频率、振幅为变量的非线性代数方程组，考虑能量平衡，构成角频率、振幅 为变量的封闭方程组求得解析解．谐波能量平衡法将谐波平衡与能量平衡相结合，克服了二者的缺点吸取 了二者的优点．实例表明，谐波能量平衡法方法简单，取较少谐波就可以达到较高的精度． 关键词　强非线性，　单摆，　谐波能量平衡法 ＤＯＩ：　１０．６０５２／１６７２６５５３２０１５０４７ 献［１４－１５］研究了采用谐波能量平衡法求解对称 引言 强非线性动力系统问题．文献［１６］研究了采用谐波 角频率是描述周期振动的最主要因素，采用通 能量平衡法求解非对称强非线性动力系统问题． ［１－２］ 常的摄动法 不能求解强非线性振动问题．近三 本文采用谐波能量平衡法对强非线性单摆方 十多年来，强非线性振动研究所取得的一系列成 程进行求解研究，并与ＫＢＭ法进行了对比． 果，其突破点一般最终都可导出振动频率的瞬变 １　谐波能量平衡法 性．比如时间变换法，椭圆函数法，频闪法，推广Ｌ Ｐ法，等效线性化方法，改进的多尺度法 ，ＦＦＴ快速 谐波能量平衡法的基本思想是把非线性微分 Ｇａｌｅｒｋｉｎ法，增量谐波平衡法和摄动增量法等［３－６］． 方程组的解，用等效的线性微分方程组的解来解析 ［７］ 逼近．首先采用谐波平衡，得到以振幅，角频率为未 张琪昌 等将待定固有频率法与规范性方法 相结合研究强非线性振动问题的求解． 知数的不完备非线性代数方程组（方程数小于未知 数）；然后利用能量守恒原理，增加关于初始条件、 文献［８］将Ａｒｎｏｌｄ舌头的结构作了推广，给出 振幅，角频率之间协调的补充方程，从而构成了关 了强非线性强迫振动解的完整结构．提出了求解强 于振幅，角频率为未知数的完备非线性代数方程 非线性简谐强迫振动的叠加迭代谐波平衡法．将 组；对这个非线性代数方程组进行求解，就可以得 描述动力系统的二阶常微分方程，化为基本解的基 到近似解析解． 本微分方程和分岔解的增量微分方程．利用叠加 研究形如 １ １