双物理摆的混沌分析.pdfVIP

第37卷第15期 数学的实践与认识 voI．37No．15 2007年8月MATHEMATICSINPRACTICEAND August，2007 THEORY 双物理摆的混沌分析 陈汉军， 黄东卫， 苏永福 (天津工业大学理学院．天津30016D) 摘要： 为了给一种新型振动机的参数设计提供理论依据．将其抽象为双物理摆的力学、数学模型．借助 Matlab软件对其振动系统发生的动力学特性进行仿真分析，得到了在给定参数n=o．01．6=o．ozl8．，= 0．04．”l；2．16．”2=3．2，厶=o．1772．以=3．4798×10—3下的力学行为的五个混沌特征． 关键词： MatIab 双物理摆，混沌，Runge．Kutta数值积分I 过去人们往往认为搜集到的许多复杂的随机信息是一种偶然现象的反映，甚至认为是 “噪声”，是实验的失败，因而弃置一旁，不予理睬．现在，人们意识到以前可能错了，这些信息 里有相当一部分或许可以归人混沌一类，而可以另辟蹊径，用混沌的方法去进行研究，或者 能从混沌中找到出路，甚至取得令人意想不到的结果．另外，鉴于混沌的复杂性，要用好它往 往也不是一件轻而易举的事，对此我们应当心中有数．本文在文献[1]的物理模型和数学分 析的基础上，利用Matlab软件对该系统的振动特性进行了给定参数下的仿真分析，并得到 了一些有用的结论． 1 实际问题 多摩擦片式混沌振动发生机构的主、从动片均为偏心盘，在 轴向压力作用下，主、从动片相互接触，当接触表面之间有相对 滑动趋势或相对滑动时，彼此间产生干摩擦力矩．其简化力学模 型为双物理摆，如图1所示． 。 设物体I、I的质量分别为，，l。和优：，质心在C和B点，物 体I相对于0轴和物体I相对于质心B的转动惯量分别为以 与JH． 图1双物理摆 2 数学模型 依力学模型及动力学方程，文献[1]有 』亏1～2．． ， △ 。 【屯4瓦万丽焉—刁习瓦驴 其中： 收稿日期；2006一03—14 50 数学的实践与认识 37卷 sin △={，，12口6 sin硼￡一研269sin(而+硼f)} zl[(z2+叫)2+训2]+优29矾一(优1+研2)口g 口=O．01，6=O．0218，，．一O．04，优l一2．16，优2=3．2，^=O．1772， 以一3．4798×10一，九=砒硼=筹 3解决方法 本问题为两个自由度的质点系，用数值计算证明在恰当的参数下，双物理摆将出现混沌 现象，为此从以下几个方面考察混沌现象： 1)对初值的敏感依赖性，即在经过较长时间以后其运动是不可预测的． 2)通过相图的演变，观察奇怪吸引子的形成． 3)通过快速傅里叶变换，作出其功率谱，混沌现象的功率谱应为连续谱． 4)通过庞加莱截面，观察奇怪吸引子在此截面上的形状． 5)通过I。yapunov指数来描述相空间相邻轨道的发散程度． 4‘理论知识 4．1庞加莱截面 为了分析相空间中的复杂轨线，庞加莱在相空问分析中 发展了一种截面方法．该方法是在相空间里取某一坐标为常 数的截面，通过研究相轨线与该截面的交点来分析系统动力 学的复杂行为．设相空间是村维的，原则上可以取出一个(以 一1)维的相平面，称为庞加莱截面．人们就可通过相轨线留 在截面上的点所组成的图象来掌握复杂的轨线情况，如图2 所示．人们将这种时间上连续的运动转变为离散的图象处理