1、第 卷 第期 年 月广西大学学报（自然科学版）（）收稿日期：；修订日期：基金资助：国家重点研发计划项目（）；四川省科技计划项目（）通讯作者：杨蔡进（），男，江苏淮安人，西南交通大学副研究员，博士；：。引文格式：陈俊，李艳，杨蔡进，等多编组铰接车辆的平顺性时域分析与优化广西大学学报（自然科学版），（）：多编组铰接车辆的平顺性时域分析与优化陈俊，李艳，杨蔡进，张德华，张卫华（西南交通大学 牵引动力国家重点实验室，四川 成都 ；成都西南交通大学设计研究院有限公司，四川 成都 ）摘要：为了分析优化多编组铰接车辆的平顺性问题，建立了考虑悬架非线性阻尼特性的节编组铰接车辆的振动模型，采用粒子群算法对影响车
2、辆平顺性的关键参数（悬架非线性阻尼参数）进行了优化，并提出了一种基于车间铰接约束力补偿和线性二次高斯（，）控制相结合的主动悬架控制方法。以车体的垂向、侧倾和俯仰振动加速度均方根作为车辆平顺性关键评价指标，采用级路面以不同速度对节编组车辆实施仿真分析。结果表明：随车速增加车辆平顺性变差，且存在首车至尾车车体垂向振动加速度均方根依次变大等振动特点，且通过悬架参数优化与采用主动悬架控制后，车辆平顺性得到明显提高 车辆在垂向、侧倾和俯仰方向的振动加速度均方根分别减少了 ，和 。关键词：多编组铰接车辆；悬架参数优化；平顺性；粒子群算法；控制方法中图分类号：文献标识码：文章编号：（），（，；，）：，（）（
3、），广西大学学报（自然科学版）第 卷 ，：；引言随着城市化进程的发展，城市人口规模不断扩大，居民的出行需求也显著增多。由公交车、地铁和轻轨等组成的城市公共交通运输体系的运输能力逐渐趋于饱和。多编组铰接车辆通过多车体模块铰接实现中等运量，可形成对现有城市公共交通运输体系的有益补充，缓解当前日趋严重的交通运输压力。多编组铰接车辆在行驶中受路面不平顺激励作用而发生振动。由于车身长大，车体振动经过传播和叠加，形成复杂振动效应，影响车辆行驶平顺性和操纵稳定性，因此对于城市客运车辆，具有良好的平顺性非常重要。通常，悬架参数优化是改善车辆平顺性的重要手段，如 等提出了一种计算效率较高的平面双横臂悬架全范围动
4、力学模型，实现悬架几何形状的快速优化设计。等基于粒子群算法与二次规划算法，建立了悬架最优参数选取原则。等采用 算法优化被动悬架参数。车辆基于上述的被动悬架参数优化方法难以匹配不同路面条件，不能达到最佳的平顺性状态。围绕悬架控制技术，国内外学者开展了有效研究，并取得一定成果。为改善悬架振动性能，提出了一种延时加速度反馈控制方法。对于四角互联空气悬架难以发挥多可控结构优势的问题，孙丽琴等设计了 控制器，实现了互联空气悬架之间的协同控制，改善了整车性能。等提出了一种自适应谐波控制算法，并与多种控制方案进行对比，该方法具有更好的性能。等提出了一种用于主动悬架系统振动控制的最优滑模控制器，并采用粒子群算
5、法求得最优控制器的参数，具有结构简单易调节的特点。等通过比对基于 控制的主动悬架使用情况，确定了两自由度汽车模型的被动悬架参数。等考虑了车辆主动悬架模型，设计了一种高阶终端滑模控制器。秦东晨等 采用系统建模和机器视觉的方法研究了一种基于完全预瞄的主动悬架控制系统。除上述研究外，基于悬架控制技术也可以改善诸如多轴特种车辆和半挂车等特殊铰接车辆的平顺性性能。孙帮成 建立了节浮车式虚拟轨道列车模型，分析线性悬架参数对列车平顺性影响。朱亚夫等 建立三轴车辆的系统模型，并设计一种基于滑模控制器的主动悬架系统。在上述研究中，尽管采用了多种方法对不同车辆进行了平顺性的分析与优化，但研究对象很少涉及多节编组车
6、辆，且研究的方向没有从车辆的编组连接关系着手，很少体现不同编组之间振动特性的差异，在主动悬架控制的研究中，针对多编组车辆提出的主动控制悬架较少。多编组铰接车辆，特别是长编组铰接车辆可用于公共运输，提升城市客运能力，是城市轨道交通未来发展的一个重要方向。目前，针对这类新制式运输系统的研究还不够充分，本文针对上述系统的平顺性问题开展前期研究，主要工作包括：建立具有节编组铰接车辆的振动模型，通过仿真分析了节编组车辆系统的平顺性特点，研究基于悬架系统的参数优化和主动控制改善平顺性的方法。相关研究可为系统研制与应用提供理论依据，具有一定理论和工程意义。多编组铰接车辆振动模型图所示为简化的节编组铰接车辆系
7、统。考虑车体垂向、侧倾、俯仰和车轮垂向运动，并考虑第期陈俊，等：多编组铰接车辆的平顺性时域分析与优化悬架非线性阻尼特性，本节建立车辆系统的多自由度振动模型，过程如下。图多编组铰接车辆示意图 基本描述假设表示第节车体的垂向位移，表示该车体的侧倾角，表示该车体的俯仰角，表示该车体上位置车轮的垂向位移，这里指标，表示编组序号，、和 分别表示左前、右前、左后和右后处车轮或悬架位置，且为了统一性，下文将采用相同记号规定。系统具有个自由度，且描述系统振动的广义坐标列阵可定义为，（）其中，。（）图悬架阻尼特性 悬架系统本文中悬架系统的力学特性可表征为线性刚度和具有分段线性函数 的非线性阻尼，如图所示，其中：
8、阻尼系数可根据附录 中参数计算获得。悬架弹簧作用力与等效弹簧变形量成正比。下文仅以第节车体左前轮处悬架作用力为例，说明计算过程。其他悬架作用力分析过程类似，不再赘述。假设车体发生小角度转动，则左前轮处悬架两端垂向相对位移，可表示为，。（）式中：为半倍轮距；为半倍轴距。将式（）对时间求导一次，可得到悬架两端点的相对速度为，。（）利用式（）和式（），该悬架作用力，表示为，。（）当系统使用主动悬架时，其作用力只需在方程（）等号右端中增加控制力项即可。车间铰接系统车间铰接系统实现车体间的机械连接。本节采用约束模型来描述车间铰接处力学行为，在不考虑转向的情况下，将铰接式公共交通车辆常采用的盘式铰广西大学
9、学报（自然科学版）第 卷接器等效为转动副。根据转动副约束特点，相邻车体的铰接位置的垂向位移相等，在不考虑车体变形的情况下，相邻车体的侧倾角相等，以第（，）个车间铰接为对象的约束方程为 ，（）式中和 表示前、后铰接点的垂向位移，计算公式为（），（），（）其中，为半倍车间轴距，因此，铰接系统的等效约束方程表示为（），。（）路面激励和轮胎力路面通常具有随机不平度，形成车体振动的激励源。假设，表示第节车体左前轮胎处路面激励，且基于滤波白噪声方法可定义为 ，（），（）式中：表示下截止空间频率，取 ；表示车速；（）表示与路面等级相关的不平度系数。本文考虑级路面条件，故该系数为 ；，表示白噪声输入。为匹配不
10、同车速工况，本文白噪声的采样时间间隔设置为 。另一方面，车轮在同一区域受到路面激励具有一定的相干性。因此，参考文献 ，设第节车体右前轮处路面激励白噪声输入为，。（）代入式（），得右前轮处路面激励，。此外，当车辆匀速直线行驶时，车辆后轮与前轮处路面激励存在时间滞后性，故后轴左右轮的路面激励可表示为，（），（），（），（），（）式中表示延迟时间，且有。最后根据前后路面输入相关性的状态空间方程得到后轮处路面激励输入 。设轮胎垂向刚度为，则轮胎力可表示为，（，），（）图路面不平顺功率谱密度图 图 路面不平顺激励时域图：（）左轮；（）后轮；（）右轮 ：（）；（）；（）图表示车速 时左、后与右轮路面输入的
11、功率谱密度图。图分别为左轮路面、后轮路第期陈俊，等：多编组铰接车辆的平顺性时域分析与优化面和右轮路面的垂向激励时域图，取第 节车的四轮激励为第节车对应的四轮滞后激励，滞后时间为（），以此得到铰接车辆多轮路面输入模型。车辆系统振动方程基于上文的推导，车辆系统的振动方程可表示为，（）式中：表示约束雅可比矩阵，由约束方程（）对广义坐标求偏导而得；表示与约束相关的拉格朗日乘子向量；由 约束稳定方法计算得到；为车辆系统的作用力的广义力阵，包括重力、轮胎力和悬架作用力等，其具体形式为，（）式中：为刚度矩阵；为阻尼矩阵；主要为重量和路面激励力向量，其中刚度矩阵具体形式为 ，（）式中表示第（，）节编组的刚度矩
12、阵，其具体形式为 （）（）（）（）。（）阻尼矩阵的形式与刚度矩阵相似；但由于采用非线性阻尼，因此各悬架处的实际阻尼值可能各不相同。主要包含路面激励作用产生的垂向作用力和车体所受重力，计算公式为，（）式中表示第（，）节编组的作用力，由下式计算：，（）表示车辆系统的质量矩阵，且有 ，（）式中 ，其中表示车体质量，表示簧下质量，和 分别为俯仰和侧倾方向的转动惯量。多编组铰接车辆平顺性时域评价指标与分析采用四阶龙格库塔法求解车辆振动方程（），获得系统振动响应，并结合相关评价指标，分析系统的平顺性特点。平顺性评价指标本文选取车体的垂向、侧倾角和俯仰角加速度均方根作为车辆平顺性的评价指标，其相关计算公式为
13、（）（），（）（），（）（）。（）广西大学学报（自然科学版）第 卷车轮相对动载荷均方根能较好反映出轮胎与地面的附着特性，当相对动载荷均方根过大时车轮附着力减小，车轮可能会脱离地面，故本文采用该指标来评价车轮与路面接触性能，其计算公式为（，）（，）（），（）式中：为重力加速度；，为车轮无量纲的相对动载荷均方根。多编组铰接车辆平顺性特点分析不失一般性，本节以节编组铰接车辆为例，结合附录中表 的车辆参数，仿真分析系统在车速、工况下的振动响应，并基于上文相关指标计算得到车体的垂向加速度、俯仰角加速度、侧倾角加速度和车轮相对动载荷均方根，研究车辆系统的平顺性，主要结果见图。图垂向振动加速度均方值 图侧倾
14、角加速度均方根 图俯仰角加速度均方根 图车轮相对动载荷均方根 由图可见，在一定车速下各节车体的垂向振动加速度均方根值呈现出依次增大的趋势，这是由车体振动传播、叠加以及刚体运动造成，此外，车速增加会导致垂向激励作用增加，从而使各节车体的垂向振动加速度均方根值变大，但首车垂向振动加速度均方根与单车系统的垂向振动加速度均方根相接近，且两者接近程度随车速增加并无显著变化。图、分别为车辆的侧倾角和俯仰角的加速度均方根值随车速增加的变化图。图表明各节车体的侧倾角加速度均方根值随车速增加先增大后减小，单车系统的侧倾角加速度均方根值相比与各编组在各个速度下都更大，由于受车间铰接作用，因此各节车体侧倾振动具有一
15、致性。如图所示，各节车第期陈俊，等：多编组铰接车辆的平顺性时域分析与优化图时速 不同车轴车轮相对动载荷均方根 体的俯仰角加速度均方根值随速度增加而增大，从车间分布来看，该评价指标在首车到第二节编组之间减小，随后依次增加，且单车系统的俯仰角加速度均方根分布情况与尾车类似。由图可见，各编组车轮的相对动载荷均方根平均值随车速增加而增大，而单车系统车轮的相对动载荷均方根值比各编组都要小。图表明从车轴方向比较，铰接车辆车轮的相对动载荷均方根平均值近似地按线性方式增长。车辆平顺性优化由上文仿真结果可知，随着车速增加，车辆的垂向振动与俯仰振动加剧，且相比于侧倾运动，车辆的垂向与俯仰振动对平顺性具有更重要的影
16、响，为抑制车辆振动，提升车辆整体的平顺性，本文分别通过悬架参数优化与主动悬架控制来提升车辆的平顺性。基于粒子群算法的悬架参数优化本节基于自适应惯性权重粒子群算法，针对悬架减震器阻尼特性参数进行优化设计，考虑到悬架动扰度对悬架刚度约束，仅对悬架阻尼系数进行优化，包括减震器压缩低速段系数、压缩中速段系数、减震器拉伸低速段系数 和拉伸中速段系数，以及相应低速过渡节点速度 和，并以最小化各节车体的垂向、俯仰角和侧倾角加速度均方根为优化目标，此外，本节在优化过程中还考虑下列约束条件：为保证阻尼大小，悬架阻尼比，满足：，（）式中和分别取 和 ；为避免悬架行程过大而触及限位块，悬架动挠度均方根应小于限位行程的，即（，），（）式中 为；轮胎应具有足够的附着力，满足车轮脱离地面的概率小于 ，保证操纵稳定性，其相对动载荷均方根应小于，即（，）。（）基于优化目标和约束条件，取粒子群规模为 ，经过 次迭代后，得到优化后的减震器阻尼特性参 数 为：，和 。表给出了在悬架参数优化前后，以车速 工况下所有编组平顺性指标的平均值大小。表时速 铰接车辆平顺性优化效果 均方根平均值优化前优化后优化效果垂向振动加速度（）侧