“天棚”阻尼控制.docx

天棚”阻尼控制“天棚”阻尼是 D.Karnopp 利用最优控制理论在1974 年提出来的一种悬架系 统主动控制策略，其控制性能优越，具有一定的鲁棒性，但由于它是基于悬架速度 的负反馈主动控制，对于移动的车辆来说无法实现但将“天棚”阻尼悬架系统作 为控制的参考模型，即把“天棚”系统作为实际系统控制的动态目标得到广泛的应 用但由于可调参数只有“天棚”阻尼系数，系统性能无法进一步提高本文采用天棚阻尼悬架作为研究对象，将分数阶微积分引入到“天棚”阻尼控 制系统中，取代原来的整数阶导数以 B 级路面为输入信号，根据优化理论找到最 优的阶数和阻尼系数最终，通过分析比较分数阶“天棚”阻尼悬架、整数阶“天 棚”阻尼悬架和被动悬架，得出分数阶“天棚”阻尼悬架能够全面提高整数阶“天 棚”阻尼悬架的性能1 车辆半主动悬架模型车辆悬架按振动控制的方法分为被动、半主动、主动 3 个类型，其中主动悬架 可很好地提高车辆的乘坐舒适性和操纵稳定性，但因其价格昂贵、能耗高、结构复 杂、可靠性差，限制了它的推广;被动悬架系统减震器的阻尼特性不能根据路面状 况和车辆运行状态进行实时的调节，因而控制效果有限;半主动悬架相比于主动悬 架，结构相对简单，能量消耗少，价格低廉，而性能接近主动悬架，特别是磁流变 材料的出现，其应用前景非常良好。

以具有两自由度的 1/4 车辆悬架模型作为研究对象，具有磁流变阻尼器的半主 动悬架模型如图 1 所示，其动力学方程:圏1平主创悬嗯模熨式中，m 簧载质量，m 簧下质量；s 悬架结构阻尼；k 悬架stcs刚度，k ――轮胎刚度;x ――车身位移，x ――轮胎位移，x――路面位移； tstgF ――半主动控制力，Fb――磁流变阻尼器的可调阻尼系数d半主动悬架是1974年由美国加州大学戴维斯分校机械工程系D. E. Karnopp教授等提出的，并利用天棚阻尼控制理论给出半主动悬架的控制策略，近十多 年来，基于各种控制理论和磁流变阻尼器技术的半主动悬架控制策略相继发表，例 LQR/LQG 控制、滑模变结构控制、自适应控制、人式神经网络控制、模糊控制、鲁 棒控制等，相比较优这些较复杂的控制理论，天棚阻尼控制方法以其简单有效一直 在半主动振动控制方面占有重要的一席之地2(分数阶天棚阻尼悬架模型2(1 天棚阻尼悬架模型天棚阻尼控制实际是形象化的最优输出反馈控制方法，见图(2)，其动力学方 程为图2天棚阻尼悬架模型这儿的b为天棚阻尼系数由于这种悬架只能是一种存在于实验室中的模型 实际的汽车上并没有“天棚”作为惯性参考，该模型不能在实际悬架上应用，必须 通过其它控制策略来间接实现，因此它被广泛地作为悬架控制的理想动力学模型应 用于其它控制方法中。

根据D.KARNOPP的天棚阻尼半主动悬架理论，半主动悬架阻尼力的控制律为:其中，f为天棚阻尼力，f为电流（或电压）最大时阻尼器的最大阻尼力；skyMRmaxf为电流（或电压）最小（0 A）时阻尼器的最小阻尼力，以美国LORD公司生MRmin产的RD-1005型磁流变阻尼器为例，此减振器允许的最大输入电流为2A,电压 12V,最大阻尼力为4448N即通过阻尼器的电流（或电压）调节部分实现天棚阻尼 控制2.2 分数阶天棚阻尼及半主动悬架控制策略天棚阻尼半主动悬架，控制策略简单且效果理想，但它毕竟只是简单的一阶系 统，控制效果无法进一步提高，而其它控制方法比较复杂，对控制器设计要求高， 下面引入新的天棚阻尼控制方法，它不仅在性能上比经典的天棚阻尼算法的进一步 的提高，而且其控制器的设计也不复杂，只是比原控制算法要求增加一些数据存贮 和乘法和加法运算原天棚阻尼算法是对悬架速度的反馈，这种算法主要通过调整阻尼系数的大小 改变控制力的值，无法改变控制力的相位，如果能实现任意调整控制力的大小和相 位，则天棚阻尼控制的效果有望进一步提高为了实现这一想法，本文引入分数微 积分思想，用对悬架位移的非整数阶导数代替原一阶导数即速度，即天棚阻尼控制 力F变为：da这儿的D表示为分数阶微分算子，其数学定义式常用的有Riemann-Liouville定义円"m-丁n列说罷m.h. I ■- iiT.niOi^Eifc』•—O -j 片也T 5-业加丁 'a'式中［?］代表取整数，当h为小量时，(6)式就是(5)式的离散表达，因此分数 阶导数的计算可以通过对函数的离散值进行加权和计算，这在控制器设计中是易实 现的。

从优化角度看，式(3)中b是唯一的可优化值，而式(5)中，b和a值都是可 优化参数也就意味着，当我们优化参数b取得最优值之后，可进一步优化参数 a,使分数阶系统比整数阶系统有更好的性能考虑到悬架位移无法测量，可以利 用指数定律，将(5)式变为类似一般文献中半主动悬架的控制方法，按照能耗原则和分数阶天棚阻尼力的计算值，半主动控制力F的控制策略为d2.3 控制参数的选取为了确定式(6 )最优的控制参数b和a,选择标准的B级随机路面激励路面谱 为路面输入，路面谱在标准ISO,TClO&SC2N67中的定义为Gg(的=曲(叫)山 如)一耳 廿＜%-1-1式中n为空间频率(m)，n=0.1m为参考空间频率，Gq(n)为参考空间频率 n0002-l所对应的路面谱值，即路面不平度系数(m,m)；3为频率指数，决定路面谱 的频-62率结构，分级路面的频率指数3,2对应于ISO的B级路面标准， Gq(n)=64X10m0-l,m，构造出时域内的随机激励路面模型本文研究的某型号轿车的悬架参数如表1后1II悬托部讣质皐■IR昭30t> kg悬牡冋阻尼蔡数LOO N/ (in.'s)轮胎垂n刚屢200000 S悬架剛厘 孔17000 y/m取控制泛函指标其中q, i=l,2,3为加权因子。

在B级公路路面行驶条件下对悬架进行仿真计算，i得到（10）式中各项和综合指标J在不同b和a时的曲线（图3）,当a=0.21, b=4500时指标J取得极小值从图3a可以看出，当a=0.31时，加 速度最小；当a=0.1时，悬架动挠度最小；当a=0.17时，轮胎动载荷最小同时 从图3a中还可以看出，当a ? （0,0.1 ）时，四条曲线都单调递减，说明如果a 在此范围内取值，那么上述悬架状态都将变得更优异当选取a=0.21时，加速度特性将更好，但悬架动挠度并非最优值3（分数阶系统的仿真目前分数阶微积分的计算方法可分为三类，一类是解析法，其主要有基于双参 数指数函数Mittag—Leffler函数的级数法、Adomian分解法和平均法等;第二类 为数值法，主要有差分法，样条函数法，预估校正法等；第三类为滤波器算法，例 如Oustaloup算法，FIR算法，IIR算法和Pade逼近等其中滤波器算法主动应用于仿真系统中的分数阶微积分的处理，其原理是用高阶的传递函数近似替代分数 阶传递函数Oustaloup算法是目前国际上最流行的近似法，首先按系统的要求选择需近似处理的频率范围[3, 3],考虑频域边界的影响，应选择更大的频率范围[3, bhA3]，要求3<3，3>3，然后在选择的范围内对分数阶微积分的拉氏变化进BAbBh行近似处理，即对于的计算采用零点和极点递归计算法其中零点3’和极点3的递归计算公式"（切厂>=（劝匚其中N为逼近函数的阶数，通常取5?10足够,G满足下式：05.结束语汽车行驶的平顺性和操纵稳定性之间是相互矛盾的，悬架的天棚阻尼控制策略 提高了汽车的平顺性，让乘客感到更舒适，但却牺牲了部分操纵稳定性。

为寻找更 优异的控制效果我们将分数阶微积分理论引入到天棚阻尼控制策略中经过仿真比 较，我们发现分数阶天棚阻尼控制策略比整数阶控制策略有更好的行驶平顺性，而 且并没有影响其操纵稳定性，相反地，若取得恰当的分数阶导数，其操纵稳定性还 会略微提高。