多自由度非线性系统周期解的同伦分析方法.pdfVIP

中国力学学会学术大会2009 （CCTAM2009 ） 1) 多自由度非线性系统周期解的同伦分析方法 *，**，2) * 钱有华 张伟 * （北京工业大学机电学院，北京 100124 ） ** （浙江师范大学数理信息学院，浙江金华 321004 ） 摘要 大部分工程实际问题可以用多自由度非线性系统来描述，这些系统的数学模型是许多个耦合的两阶常微分 方程。一般地，要精确求解这些方程非常困难，因此可以考虑它们的解析近似解。同伦分析方法是解非线性系统响 应的有用工具，本文将它应用于多自由度非线性系统周期解的求解中。利用多自由度耦合 van der Pol 振子和非自治 薄板系统的实例，展示了同伦分析方法的有效性和巨大潜力。同时，把得到的解析近似解与系统的 Runge-Kutta 数值 解作了比较，结果表明同伦分析方法是求多自由度非线性系统周期解的有效方法。 关键词 同伦分析法，非线性系统，多自由度，周期解 引 言 文献[1]系统地描述了求解非线性常微分方程和偏微分方程的一种解析近似方法——同伦分析方 法。该方法自 1992 年廖世俊[2]在他的博士学位论文中提出后，经过十几年的发展，在非线性振动、 非线性波浪、多孔介质中的黏性流动、非线性热传导、Navier-Stokes方程、Thomas-Fermi方程和Gelfand 方程等方面得到了应用。1998 年，廖世俊[3]将该方法首次应用于非线性振动的研究中；2003 年，廖 [4] [5] 世俊 用该方法研究了正阻尼系统的自由振动；2004 年，廖世俊 用该方法研究了自激系统的自由 振动；温建明等[6-7]用该方法研究了参数激励非线性振动系统周期解的解析近似；2008 年，廖世俊[8] 对同伦分析方法作了综述。同伦分析方法在解决非线性问题的过程中发挥了越来越大的作用，它不 仅克服了摄动方法依赖于小参数的局限性，而且在逻辑上包含了其他一些非摄动方法，如Lyapunov 人工小参数法、Adomian分解法和δ展开法。因此，该方法应用范围广，具有一般性，与其他解析近 似方法不同，该方法提供了一个简便的途径来调节和控制所得级数解的收敛区域和收敛速度。所以 该方法既适用于弱非线性系统又适用于强非线性系统。 由于工程实际中很多问题发生的都是大挠度、大变形的振动，当发生振幅较大的振动时，不能 应用线性理论描述结构的大变形。非线性的大振幅振动可以引起飞机机翼折断，卫星星体的大挠度 振动可能引起卫星系统的内能耗散，最后引起卫星姿态失稳而翻滚。因此用非线性理论研究工程中 的振动模型具有重要的理论和实际意义。两自由度耦合van der Pol振子是典型的非线性振动模型，很 多实际工程中的动力学问题最终可简化为该模型的求解。随着科学技术的发展，薄板被广泛应用于 航空航天等工程领域。近十几年来，关于薄板的非线性振动、分叉和混沌动力学的研究取得了一些 进展。2001 年，张伟[9]研究了参数激励下的矩阵薄板的全局分叉和混沌动力学；2007 年，姚明辉等 [10]研究了参数激励和外激励联合作用下的矩阵薄板的多脉冲Shilnikov型轨道和混沌动力学；2008 年， 张君华等[11]利用广义Melnikov方法研究了参数激励下的矩阵薄板的多脉冲Shilnikov型混沌动力学。 本文将同伦分析方法应用于两自由度耦合van der Pol振子和非自治薄板系统周期解的求解中。同时， 把得到的解析近似解与系统的Runge-Kutta数值解作了比较，结果表明同伦分析方法是求多自由度非 线性系统周期解的有效方法。 1 耦合van der Pol 振子的同伦分析方法 1) 国家自然科学基金项目 资助 2) E-mail: zjjhqyh@ 1