两自由度碰撞振动系统的Lyapunov指数谱分析.pdf

振动与冲击 第28卷第1期 JOURNAl．OFVIBRATIONANDSHOCK 两自由度碰撞振动系统的Lyapunov指数谱分析 魏艳辉1，。徐洁琼2，黄龙生1 (1．浙江林学院天目学院，l临安311300；2．北京航空航天大学理学院，北京100083) 摘 要：针对一类两自由度碰撞振动系统。取碰撞前瞬时的定相位面为Poincar6截面，引入局部映射，构造 敛性很好．为了验证该计算方法的正确性和有效性，分析当系统参数在大范围内变化时，相应的最大Lyapunov指数的变 化规律。 关键词：碰撞振动；Lyapunov指数谱；Poincar6映射；局部映射 中图分类号：0322 文献标识码：A 碰撞振动问题产生于机械、车辆和核反应堆等多 成功地应用到单自由度碰撞振动系统。文[10]讨论了 个工程领域，其研究具有重要的意义．许多学者对碰振 系统的动力学行为进行了广泛的研究．文[1—2]推导征乘子之间的解析关系，提出了非光滑碰振系统动力 两自由度碰振系统存在Hopf分岔和次谐分岔的条件，学分析的Lyapunov指数判据。 系统地分析了系统的两参数动力学性态，研究了Hopf． 本文应用文[9]的方法，研究了一类两自由度碰撞 Flip和Hopf-Hopf分岔等余维二分岔，通过数值模拟进 一步讨论了系统的，r2环面分岔。文[3]针对两自由度用数值模拟，计算了系统周期吸引子和混沌吸引子的 碰振系统，提出了系统参数依据系统输出的变化自动 调节连续反馈控制的策略．文[4]研究了一类双自由度 性和有效性，讨论了当参数在大范围内变化时，相应的 碰振系统，给出了次谐周期运动的存在性、稳定性和共 最大Lyapunov指数谱变化情况。 存性的结果。 1力学模型及其运动方程 Lyapunov指数谱是一种判定动力系统运动稳定性 和混沌特性的重要工具。它表示相空间内邻近轨迹的 两自由度碰振系统的力学模型如图1所示，质量 平均指数发散率的数值特征。光滑动力系统Lyapunov为M，，M：的质块分别由刚度为K．，恐的线性弹簧和阻 指数谱的研究目前已比较成熟。在非光滑动力系统方 尼系数为C，，C：的线性阻尼器相接，两质块只作水平 面，由于非光滑性使得光滑动力系统Lyapunov指数谱 的计算方法不能直接应用到非光滑动力系统上。近年 来，一些学者针对某些特定的非光滑动力系统，研究了 块肘．将与刚性约束A碰撞．假设力学模型中的阻尼是 其Lyapunov指数谱的计算方法。文[5一]针对两个 相同的具有不连续性的离散动力系统，运用同步现象 碰撞恢复系数尺确定，碰撞持续时间忽略不计。 计算该系统的最大Lyapunov指数。文[7]针对不连续 的非线性动力系统，在不连续的瞬时，附加一定的转换 条件，将光滑动力系统Lyapunov指数谱的计算方法应 用到了不连续的非线性动力系统。文[8]针对两个碰 撞振动系统，在连续两次碰撞之间，引入超越映射，将 光滑离散动力系统Lyapunov指数谱的计算方法应用到 该系统。文[9]引入局部映射，利用Poincar6映射分析 方法及Gram．Schmidt正交化和范数归一化方法得到了 图l两自由度碰振系统的力学模型 非光滑动力系统Lyapunov指数谱的计算方法，该方法 在任意连续两次碰撞之间，系统的运动微分方 程为 基金项目：广西自然科学基金(0640002)资助项目 收稿日期：2008—02—29修改稿收到日期：2008—05—15 第一作者魏艳辉女，硕士，助教，1980年11月生 万方数据 第1期 魏艳辉等：两自由度碰撞振动系统的Lyapunov指数谱分析 6l Xl+=一麟．．(X，=B)， (2) 式中一，+分别表示碰撞前后瞬时．