1、 年第 卷 月第 期机 械 科 学 与 技 术 ：收稿日期：基金项目：国家自然科学基金项目（）、安徽省高校优秀拔尖人才培育项目（）及安徽省教育厅自然科学重点项目（）作者简介：王娜（），副教授，硕士，研究方向为汽车动力学，王娜，毛忠民，任翠锋，等汽车摆振系统 分岔及参数灵敏度分析机械科学与技术，（）：汽车摆振系统 分岔及参数灵敏度分析王娜，毛忠民，任翠锋，高大威（安徽三联学院 机械工程学院，合肥；上海理工大学 机械工程学院，上海）摘要：为了分析汽车摆振系统的失稳机制，以某样车的右前轮为例，基于 车辆动力学模型，通过第二类拉格朗日方程建立了三自由度单轮摆振系统动力学模型。在系统状态方程的基础上，应
2、用 准则求解得到了摆振系统的临界失稳车速，进而获得了摆振系统在临界失稳车速处的特征值。结果表明，临界车速工况下的系统平衡点是非双曲平衡点，因此，以车速为分岔参数，通过中心流形理论得到了系统在临界车速处的二维中心流形。在此基础上深入分析了摆振系统的 分岔特性，并得到了摆振系统的近似周期解。最后，通过四水平五参数正交实验设计，获得了转向系和悬架系的结构参数对系统临界失稳车速的灵敏度。研究发现，系统临界失稳车速对主销后倾角最为敏感，其次是转向系和悬架系的阻尼。因此，在工程实践中，应该优先考虑通过调整主销后倾角来抑制汽车摆振。关键词：机摆振；中心流形；分岔；灵敏度分析中图分类号：文献标志码：文章编号：
3、（），（，；，）：，：；摆振是一种复杂的自激振动现象，广泛存在于汽车、摩托车和飞机起落架等动力学系统中。摆振现象的发生会对汽车的行驶安全性和操纵稳定性产生不利影响，因此，为了揭示摆振的发生机理，进而完善汽车防摆振设计，许多学者围绕摆振动力机 械 科 学 与 技 术第 卷：学开展了大量研究。由于汽车摆振主要表现为前轮绕主销的角振动，因此，轮胎力学特性与摆振的发生密切相关，一些学者结合轮胎力学特性深入研究了汽车摆振系统的失稳机理、能量传递规律及分岔特性。此外，部分学者也关注到了转向系统间隙和干摩擦等强非线性因素对汽车摆振的影响规律，研究发现间隙主导下的汽车摆振系统动力学响应将表现为混沌运动，而干摩
4、擦会诱发摆振系统的多极限环运动。随着对汽车摆振的研究不断深入，非线性振动理论也被用于摆振动力学机理的研究之中。其中，学者围绕经典的三自由非独立悬架汽车摆振动力学模型深入探究了摆振系统 分岔特性，指出摆振现象是系统在发生 分岔后出现的极限环运动。此外，一些学者也针对不同类型的车辆摆振系统开展了研究工作。肖闯和殷智宏建立了轮椅车轮摆振系统动力学模型，应用 分岔理论分析了系统的稳定性。等基于一个十二自由度摆振模型讨论了双前桥重型汽车的多极限环振动特性，并分析了路面附着系数对极限环的影响。等建立了五自由度重型车辆摆振模型，讨论了油气悬架参数对系统稳定性的影响。上述研究成果为进一步研究汽车摆振提供了良好
5、的理论基础。然而，为了简化动力学分析，大多数学者忽略了汽车轮胎、转向系和悬架系统之间的动力学耦合作用，从而使得汽车垂向运动引起的动载荷对汽车摆振的影响尚未明确。需要指出的是，轮胎垂向载荷的变化对汽车动力学稳定性有着重要的影响。因此，本文基于单轮摆振系统动力学模型，研究考虑轮胎垂向载荷影响的摆振系统分岔特性。此外，由于影响汽车摆振的因素较多，目前尚不清楚汽车摆振对哪些因素的变化较为敏感，所以本文将对转向系和悬架系统参数进行敏感性分析，相关结论可为工程中的摆振抑制提供有力的技术支持。摆振系统的动力学模型由于独立悬架汽车的左轮和右轮之间的相互作用较小，因此仅以某样车的右前轮为研究对象，其简化动力学模
6、型如图 所示，基于以下假设建立摆振系统的动力学模型。）转向系统的阻尼和刚度分别等效为车轮绕主销旋转的角阻尼和角刚度。）仅考虑轮胎的垂向刚度，忽略轮胎的垂向阻尼。）非簧载质量主要集中在轮心。图 汽车摆振系统动力学模型显然，动力学模型共有 个自由度，包括前轮摆角、簧载质量的垂向位移和非簧载质量的垂向位移，相关车辆结构参数如表 所示。表 某样车结构参数 参数数值车轮围绕主销的转动惯量 车轮绕其旋转轴的转动惯量 非簧载质量 簧载质量 转向系统等效角阻尼 悬架阻尼 转向系统等效角刚度 轮胎的垂向刚度 悬架刚度 轮胎滚动半径量 轮胎拖距 从车轮中心到主销的横向距离 车轮中心到下控制臂和副车架铰链中心距离
7、轮胎的松弛长度 接地印迹长度 后倾角 轮胎公式拟合参数 轮胎公式拟合参数 轮胎公式拟合参数 轮胎垂向静载荷 第 期 王娜，等：汽车摆振系统 分岔及参数灵敏度分析：动力学方程基于第二类拉格朗日方程，图 所示的摆振系统的动能、势能、耗散能分别为：（）（）（）（）（）各自由度对应的广义力为：（）（）（）（）式中：为轮胎的侧偏力；为车轮中心的垂向合力，为作用于车轮中心的地面垂向力，（），为作用于车轮中心的悬架垂向力，（）（）。将式（）式（）代入第二类拉格朗日方程，可得到摆振系统的微分方程：（）（）（）（）（）（）（）（）轮胎模型由于轮胎的垂向载荷对侧向力有着至关重要的影响，考虑垂向载荷影响的侧向力表达
8、式为 （）（）（）式中：为侧偏角；。汽车右前轮侧偏角和摆角的约束方程为 （）分岔特性为了分析摆振系统的稳定性，引入状态变量 （，），。以车速为分岔参数，将摆振系统的微分方程转化为状态矩阵。（，）（）临界车速令（，），可求得 是摆振系统的一个平衡点。因此，系统的状态矩阵可转换为（）（，）（）式中：（）为摆振系统在零平衡点处的雅可比矩阵，可得（），；（，）为摆振系统的非线性项。因此，系统特征方程可由 （）求解得到。（）式中：；。根据 准则，特征方程（）具有一对纯虚根，且其余特征根具有负实部的判别条件为：），）将表 中的参数值代入 行列式，便可求得汽车临界失稳车速，如表 所示。显然，仅有.和.具有实
9、际物理意义。将车速.代入其余的 行式中，得到 .，.机 械 科 学 与 技 术第 卷：。因此，摆振系统在此车速工况下处于临界稳定状态。同样地，当车速 .时，.，.，摆振系统也处于临界稳定状态。表 临界车速的求解（）.中心流形为了进一步研究摆振系统的分岔特性，将上述两个临界车速代入方程（），即可得到摆振系统的特征值，如表 所示。显然，摆振系统的零平衡点在两个临界失稳车速工况下是非双曲的。因此，以车速 为例，应用中心流形理论对原摆振系统进行降维。表 摆振系统的特征值特征值.引入新变量，令，且，则可以将状态矩阵转换为 （，）（，）（）式中：可重新组合系统特征向量的实部和虚部获得；为一个对角矩阵。则中
10、心函数可设为（，）（）式中：，；初始条件（，），（，），（，），（，）；为高阶无穷小。此外，对式（）求导可得（，）（，）（，）（）将式（）和式（）代入式（），可获得中心流形为：|（）上述中心流形可在极坐标系下转换为：（）式中 可分别从式（）中获得。对于临界失稳车速 ，根据上述步骤，可获得中心流形如下：|（）同样，在极坐标系中，式（）转化为：.（）第 期 王娜，等：汽车摆振系统 分岔及参数灵敏度分析：式中 可结合式（）求得。令，可结合式（）和式（）中得到摆振系统的分岔图，如图 所示。从图 中可以看出，当车速在 时，摆振系统零平衡点是不稳定的，系统将会由于 分岔而出现极限环振动现象，即汽车摆振。当
11、车速 或车速 时，系统零平衡点将成为稳定的焦点，从而摆振现象消失。图 摆振系统的分岔特性 数值验证基于龙格库塔算法，对不同车速下的摆振幅值和垂向载荷进行了数值求解，结果如图 所示。图 摆振幅值与轮胎垂向载荷最值随车速的变化由图 可以看出，数值结果与上述求得的临界车速结果一致。此外，摆振幅值将随着垂向载荷的增加而增加。因此，垂向载荷对摆振幅值有重要的影响。在工程中，可以通过减小垂向载荷来抑制汽车摆振。通过将式（）和式（）的变量转换为原摆振系统变量，可分别获得图 所示的摆角相图和表 所列的摆振幅值。显然，用中心流形方法得到的结果与数值结果基本一致。因此，验证了动力学模型的准确性，也表明中心流形理论
12、可以有效求解摆振系统微分方程。图 两种方法求得的摆角相图（）表 不同速度下的摆振幅值车辆速度（）数值结果 解析结果.灵敏度分析汽车作为一个复杂的机械动力学系统，各部件及子系统之间存在着重要的耦合作用，因此影响汽车摆振的因素较多，工程中也存在多种途径来抑制汽车摆振。然而，轮胎半径、侧偏刚度以及整车质量等参数在工程实际中很难进行调整优化，鉴于此，本文以易于调整的转向系和悬架系参数为研究对象，基于正交实验设计对其进行灵敏度分析，从而为工程实际中更好地进行防摆振设计提供技术参考。首先，选择 个因素（即转向系统的刚度和阻尼、悬架的刚度和阻尼以及主销后倾角）进行灵敏度分析，以分析各因素对临界车速的影响程度
13、。对于不同的参数，利用 准则计算出两个临界车速（即 和），并以临界车速之差（即 ）作为实验结果。各因素的水平如表 所示。机 械 科 学 与 技 术第 卷：表 各因素的水平表水平因素（）（）（）（）.选择标准的四水平五参数正交表（），并根据（）的实验安排来求解临界车速。正交实验的安排和结果如表 所示。表 正交实验的安排和结果编号因素水平（）.由正交实验结果可得各因素的主效应图，如图 所示。图 各因素主效应图由图可以看出，增加悬架阻尼、转向系阻尼和刚度可以使得摆振系统更加稳定，从而抑制汽车摆振。主销后倾角和悬架刚度的减小也可达到抑制摆振的效果。对于四水平正交实验，各因素的灵敏度计算公式为 （）式中
14、：?，?，?，?，?，?，?，?为因素 在 水平下的所有实验结果的平均值；为因素 在 水平处的值；为因素 在 水平处的值。各因素的灵敏度如表 所示。表 各因素的灵敏度因素灵敏度 由表 可以看出，汽车摆振系统失稳车速区间对主销后倾角最为敏感，其次是转向系阻尼。各因素的敏感度排序为：。因此，在工程实践中，应优先考虑通过调整主销后倾角来抑制汽车摆振。此外，调整转向系和悬架系的阻尼也是抑制汽车摆振的有效方式。值得注意的是，汽车失稳车速区间对转向系和悬架系统的刚度变化不敏感，在工程实际中，不建议通过调整系统刚度来抑制汽车摆振，这是由于尽管大范围的调整刚度，汽车摆振不仅得不到抑制，而汽车其他动力学性能将会
15、被恶化。结论基于汽车单轮摆振系统动力学模型，采用 准则求得了系统的临界分岔车速，进而根据系统在临界车速处的特征值，分析发现摆振系统的平衡点是非双曲的。因此，利用中心流形理论分析了系统的分岔特性。研究表明，在一定的速度范围内，由于 分岔的发生使得系统失稳而出现摆振现象，且摆振幅值与轮胎的垂向载荷密切相关。值得注意的是，较大的垂向载荷将导致更大的摆振幅第 期 王娜，等：汽车摆振系统 分岔及参数灵敏度分析：值，因此可以通过减小垂向载荷来抑制汽车摆振。此外，选取转向系和悬架系统的结构参数作为 个因素，并将 准则求得的汽车失稳车速区间作为实验结果。结合灵敏度分析可知，汽车失稳的车速区间对主销后倾角最为敏
16、感。因此，在工程实践中，优先考虑通过减小后倾角来抑制汽车摆振。此外，增加转向系和悬架系统的阻尼来减小汽车摆振也是直接有效的方法。然而，汽车失稳车速区间对转向系和悬架系统的刚度不敏感，工程中不建议通过调整系统刚度来抑制汽车摆振。参考文献 张宁，殷国栋，陈南，等车辆动力学中的摆振问题研究现状综述机械工程学报，（）：，（）：（），（）：谢帅，舒成辉，贾天娇，等起落架摆振多体动力学建模及其在摆振飞行试验中的应用机械科学与技术，（）：，（）：（），（）：，（）：，（）：，：，（）：，（）：，（）：王威，陈军，宋玉玲含双间隙模型轿车转向系非线性动力学特性研究农业机械学报，（）：，（）：（），（）：张琪昌，李小涛，田瑞兰汽车转向轮摆振的稳定性及分岔行为分析振动与冲击，（）：，（）：（）刘献栋，张紫广，何田，等基于增量谐波平衡法的汽车转向轮非线性摆振的研究汽车工程，（）：，（）：（）肖闯，殷智宏单转向轮摆振非线性动力学研究中国机械工程，（）：，（）：（），（）：，（）：，（）：，：，：，刘文卿实验设计北京：清华大学出版社，：，（）