多自由度转子系非线性动力学数值分析及混沌控制.pdf

摘要 本文继承线性常微分方程精细积分法高精度等特点，提出适合于多自由 度非线性动力学方程的数值求解方法。给出了求解非线性动力学方程的精细 积分表达式，利用插值近似该方程的非线性项。得到一个具有四阶精度并且 是单步显式、自起步、预测．校正的Lagrange(或Hermite)插值精细积分算 法，适于强非线性、非保守系统。对于多自由度非线性动力学方程组分块地 应用插值精细积分算法。从而降低转换矩阵阶数，避免高阶矩阵、向量运算， 有效地提高了计算效率，得到多自由度非线性动力学方程组的分块精细积分 法。另外提出了求解动力学方程的一个新型的逐步积分法，基于线性动力学 方程的解析齐次解及Duhamel积分，构造出适用于非线性动力学方程解的一 般积分表达式，对包含非线性项的非齐次项采用插值近似的方法，得到一个 单步显式、自起步、预测校正具有四阶精度的解析逐步积分算法。借助分块 积分的概念，得到了一个求解线性动力学方程的积分型(相对于差分而言) 直接积分法。算例表明解析逐步积分法和积分型直接积分法比中心差分法、 Newmark、Wilson-0、Houbok法等有较高的幸膏度。 本文采用动载短轴承非稳态非线性油膜力模型，给出了多自由度柔性 Jeffcott转子系统稳定性及其Hopf分岔的适用数值方法，利用李雅普诺夫稳 定性理论和区间二分法代替R-H颤振失稳准则，求解系统的失稳角速度国． (即Hopf分岔点)。数值求解系统动力学方程，在国，处，转子轴心轨迹为极 小幅值的周期涡动，此为Hopf分岔，分又后，系统运动为拟周期运动。该方 法较R-H准则更适于多自由度系统。 研究了非稳态动载短轴承支撑的Jeffcott柔性转子系统的动力特性，基于 可用三个函数表示动态万油膜的非稳态非线性油膜力模型，将转速比、不平 衡量、阻尼比、粘度作为控制参数，利用Floquet乘子预测周期解的局部稳定 性，发现倍周期分叉是由一定量的不平衡引起的，而Hopf分叉是由油膜失稳 造成的。通过Lagrangc插值精细积分法数值给出系统的轴心轨迹图、Poincar 6映射图、分叉图，检验Floquet理论预测结果并预测系统的长期性态，显示 系统在四个参数组合的某些范围内具有丰富的菲线性特性，还存在多形式次 谐波解，以及由倍周期分叉、二次Hopf分叉通往混沌的现象。利用理论与数 值方法，在较宽的参数范围内研究了不平衡柔性转子系统运动的局部稳定性、 分叉和长期性态，结果显示二者在求周期解的局部稳定性是一致的，但数值 方法可分析周期解分叉以后系统随参数变化的动力特性，这是F|oque=t理论的 局限性。计算结果表明不平衡参数、阻尼比及粘度越大，对系统的非同步运 动就越具有抑制作用，能提高油膜失稳转速。结合这两种方法为油膜轴承柔 性转子系统相关参数的适当设计提供了定性分析的依据，将动态油膜力模型 和稳态油膜力模型的数值结果进行比较，表明动态非线性油膜力模型的合理 性。 考虑具有横向裂纹柔性轴Jeffcott转子在裂纹开闭、转轴剐度和转轴位 移三者耦舍的影响。构造了裂纹转子，动态池膜轻承非线性躅期对变系绕模 型，研究了转子系统在非线性涡动影响下的分叉现象及其分叉后的力学行为。 通过数值计算表明，随着裂纹轴剐度变化比的增大，系统在低转速区域内具 有丰富的菲线性动力行为，出现倍周期分又及混沌现象。涡动振幅随转速升 高而减小，直到非稳态非线性油膜失稳。在无裂纹转子油膜临界失稳点处发 现了类Hopf分叉现象，即系统运动由平衡变为拟周期运动{裂纹转子在油膜 临界失稳时的系统由周期运动变为拟周期运动，此为第二Hopf分叉。裂纹 转子轴刚度变化对油膜失稳点及油膜失稳之后转子的运动影响不大，转子系 统作拟周期运动。 本文利用混沌系统的敏感性、共存性及各态历经性，提出控制混沌的两 个新方法，并应用控制方法控制油膜轴承-柔性转子系统的混沌运动。一个方 法是小参数干扰反馈控制动力系统中的混沌运动，在混沌吸引子内不稳定周 期轨的邻域，通过一个小的参数干扰反馈，将混沌轨逼近并稳定在该不稳定 周期轨上。另一个方法是压缩映射．参数微扰控制漫沌，在漫沌吸引子内的不 稳周期轨的邻域，利用一个压缩映射和线性逼近，得到一个参数微扰反馈序 列，使混沌轨快速靠近并稳定在周期轨上，压缩映射参数可在f．1。1)区间内 灵活选取，用本方法控制动态稳膜．柔性转子系统和Duffing振子系统的混沌 运动，数值结果袭明该控制方法不仅能控制低维动力系统的混沌，也能有效 控制高维动力系统的混沌运动