含非线性摩擦系数的单自由度颤振系统全局动力学行为数值研究.pdfVIP

第28卷第3期 计 算 力 学 学 报 Vo1．28，NO．3 2O11年 6月 ChineseJournalofComputationalMechanics June2011 文章编号 ：1007—4708(20l1)03—0417—06 含非线性摩擦系数的单 自由度颤振系统 全局动力学行为数值研究 唐进 元 ， 熊兴波 ， 陈思雨 (现代复杂装备设计与极端制造教育部重点实验室 (中南大学)，中南大学 机电工程学院，长沙 410083) 摘 要 ：摩擦系数模型取更具普适性的Stribeck非线性模型 ，基于事件驱动理论 ，利用 C与 Matlab联合仿真 的方 法开发 了干摩擦颤振 问题 的快速求解程序。给 出改进 的胞 映射算法 ，对含非线性摩擦的单 自由度摩擦颤振系统 的演化过程及其全局性态进行数值分析和研 究，得到系统在任意的初始状态下的响应特性、系统收敛域 的数值计 算分析结果 。分析结果表 明：在激励速度 不断减小 的过程 中，系统 的稳定解 由稳定 的平衡 点转变为稳定 的极 限环 ，并依次出现颤振 、粘滑等非线性动力学现象；在整个相平面 内，系统的收敛特性不发生改变，总是收敛于 同 一 稳定 的奇 点或 同一极 限环。 关键词 ：摩擦颤振 ；胞映射 ；非线性；全局性态 中图分类号 ：O322 文献标志码 ：A 段的近似解析分析和整个运动阶段的数值模拟 ；颜 1 引 言 海燕等n 也对转 向架的摩擦颤振分岔现象进行 了 业 已知道摩擦是产生颤振 的主要原 因口 ，而 分析 ；梁爽等[】用实验 的方法对摩擦颤振系统进 各种类型的机械系统 中，摩擦广泛存在，引起 了一 行 了研究。对于这类非光滑系统来说 ，精确的仿真 系列的颤振问题 ，对机械系统的性能产生了破坏性 计算是非常困难 的。主要表现为，粘滑计算时间非 的影响。如铁道机车车辆运行 中所产生的颤振，对 常长 ；同时由于干摩擦 导致的粘滑运动 ，使数值求 车体、转向架和悬挂系统产生较大的破坏作用 ，颤 解计算时，精确判断粘滑进人和退出的时间非常困 振也对提高车辆 的运行平稳性 、减小铁道沿线的噪 难D7,183。Piiroinen和 Kuznetsov=[1叼提出了一种基 音污染产生了巨大的负面影响。 于事件驱动的求解 Filippov系统的数值计算方法 ， 胞映射方法是求解非线性复杂系统全局性态 可以判断粘滑进入退出的时间，但是计算 的时间非 的一种数值计算技术 ，结合多线程并行处理的可视 常长。尽管对摩擦颤振有较多的研究 ，但对全 化技术，能够实时地对非线性动力学体系的复杂演 局动力学行为的研究并不多见，尤其对含非线性摩 化过程、奇异吸引域和混沌性态等进行有效、动态 擦系数的单 自由度颤振系统全局动力学行为的研 和准确地解析Ll¨]。胞 映射理论 由 HsuCS创 究尚没见有文献发表。 建 引，其理论基础是庞加莱 (Poincar6)点映射原 基于上述问题 ，本文主要研究内容和特色 ：(1) 理 、马尔科夫链 (MarkovChains)理论和相空 间重 引进非线性 Stribeck干摩擦模型，建立了更加合理 构原理 。 的干摩擦颤振模型。(2)基于事件驱动理论 ，利用 对摩擦颤振现象，国内外学者进行过大量的研 C与 Matlab的联