基于无源化的非线性弹簧悬架主动控制方法

摘要

本发明公开了一种基于无源化的非线性弹簧悬架主动控制方法，在考虑悬架弹簧非线性刚度特性的基础上，对1/4车辆模型描述的悬架系统进行了主动控制研究，通过分析非线性弹簧悬架系统数学模型的特点，选定合适的期望互联矩阵、阻尼矩阵及能量函数，基于无源化控制理论，获得了非线性弹簧悬架主动控制器，该主动控制器适用于恶劣地况，可以有效降低车身加速度幅值及控制量幅值，具有更快的收敛速度，并且在4‑8赫兹这个人体敏感的频率区域内，车身加速度的振动分贝数也有明显下降。

权利要求

1.基于无源化的非线性弹簧悬架主动控制方法，其特征在于：包括以下步骤，构建非线性弹簧悬架模型；

分析非线性弹簧悬架模型的特点，选定期望互联矩阵、阻尼矩阵及能量函数，基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架主动控制器；

非线性弹簧悬架模型为，

其中， 为状态变量导数，状态变量 z1为车轮的位移，z2为车身的位移， 为车轮的速度， 为车身的速度，x1＝z2，x2＝z1，p为速度分量，p＝[z2,z1]T，q为位移分量，m1为车轮的重量，m2为车身的重量，c2为非线性弹簧悬架阻尼系数，k1为轮胎等价刚度系数，u＝fa，fa为

非线性弹簧悬架主动控制力，k2为非线性弹簧悬架刚度系数，ε为弹簧非线性程度。

2.根据权利要求1所述的基于无源化的非线性弹簧悬架控制方法，其特征在于：非线性弹簧悬架模型的特点为，A12≠-A21， A22≠02×2。

3.根据权利要求2所述的基于无源化的非线性弹簧悬架控制方法，其特征在于：选定的期望互联矩阵Jd(x)为，q12＝q22P22P11-1z1z2-1，q13＝P33-1，q24＝P44-1，Pii＞0为可设计参数，qjj＞0为可设计参数，i＝1,2,3,4，j＝2,3,4；

选定的阻尼矩阵Rd(x)为，

q11＝q22(P22P11-1z1z2-1)2；

选定的能量函数Hd(x)为，

P＝diag{P11,P22,P33,P44}。

4.根据权利要求3所述的基于无源化的非线性弹簧悬架控制方法，其特征在于：基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架控制器为，即在非线性弹簧悬架控制力的作用下，非线性弹簧悬架对于平衡点x*＝0渐近稳定的，并且可变成端口受控的哈密尔顿系统结构

5.根据权利要求4所述的基于无源化的非线性弹簧悬架控制方法，其特征在于：基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架控制器的具体过程为：由选定的能量函数Hd(x)可知，T

P＝P＞0

对于平衡点x*＝0，有 为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的一阶梯度， 为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的二阶梯度；

根据选定的期望互联矩阵Jd(x)和阻尼矩阵Rd(x)，可得，T

Jd(x)+Jd(x)＝0

Rd(x)＝RdT(x)≥0

由非线性弹簧悬架模型的哈密尔顿系统结构公式、选定的期望互联矩阵Jd(x)以及选定的阻尼矩阵Rd(x)可知，其中，F1(x)＝[-q11 q12 q13 0]，F2(x)＝[-q12 -q22 0 q24]，F3(x)＝[-q13 0 -q33 0]，F4(x)＝[0 -q24 0 -q44]，假设g(x)的左零矩阵为，

其中，γ1＞0，γ2＞0，γ1和γ2为可设计参数，则，

其中， 为非线性弹簧悬架系统非理想能量函数的一阶梯度；

由于

基于选定的期望互联矩阵Jd(x)可知，因此 成立，

进而根据无源化的互联和阻尼分配控制理论可知，存在使得非线性弹簧悬架对于平衡点x*＝0渐近稳定的，并且可变成端口受控的哈密尔顿系统结构；

由于[gT(x)g(x)]-1gT(x)＝[0,0,m2,-m1]，所以基于选定的期望互联矩阵Jd(x)可得，则根据 可得非线性弹簧悬架控制器。

6.根据权利要求5所述的基于无源化的非线性弹簧悬架控制方法，其特征在于：通过选定可设计参数可对非线性弹簧悬架控制器进行简化，具体为，选定 q33＝(m2P33)-1c2，q44＝(m1P44)-1c2时，非线性弹簧悬架控制器可简化为，

此时有

进一步选定P33＝P44，则

此时，非线性弹簧悬架可变成端口受控的哈密尔顿系统，可知，q22和P44的选择决定了非线性弹簧悬架控制器特征值。

说明书

基于无源化的非线性弹簧悬架主动控制方法

技术领域

[0001] 本发明涉及一种基于无源化的非线性弹簧悬架主动控制方法，属于非线性弹簧悬架控制技术领域。

背景技术

[0002] 由于工程车辆、农用拖拉机和军用车辆等等非路面车辆的行驶道路条件差，激振源复杂(发动机振动、道路不平等)，振动问题尤为突出。近几年，悬架主动控制研究受到了国内外诸多学者的广泛关注，在理论分析和物理试验方面均开展了一些研究工作。主动悬架系统对于减小车辆振动以及提高车辆行驶平顺性和驾驶舒适性，都具有重要意义。Elnashai AS(Analysis of the damage potential of the Kocaeli(Turkey)earthquake of 17August 1999[J].Engineering Structures,2000,22(7):746―754.)针对带有磁流变减振器的二自由度座椅悬架系统，设计了基于积分滑模的半主动控制器。施仲衡(地下铁道设计与施工[M].西安:陕西科学技术出版社,1997.)针对1/4车辆模型设计了悬架主动控制器，保证车身加速度在合理频带内，满足驾驶舒适度要求。然而，以上文献均将悬架振动运动视为线性系统，事实上汽车悬架是典型的非线性系统，如弹簧的非线性刚度特性、减振器的不对称阻尼特性等，因此，研究悬架的非线性特性更符合悬架的实际情况。Wang J N(Seismic design of tunnels:A simple state-of-the–art design approach[M].New York:Parsons Brinckerhoff Quade and Douglas Incorporated,1993.)针对变刚度弹簧非线性悬架，采用基于线性二次型最优控制方法(LQR)，设计了主动控制器，但是仅就B级路面进行仿真。川島一彦(地下構造物の耐震設計[M].日本:鹿島出版会,1994)为非线性悬架设计了backstepping控制器，但是仅针对小幅冲击和正弦形式描述的地面激励给予了仿真验证。因此急需一种控制方法能保证悬架系统适用于大幅冲击激励和E级恶劣地面。

发明内容

[0003] 为了解决上述技术问题，本发明提供了基于无源化的非线性弹簧悬架主动控制方法。

[0004] 为了达到上述目的，本发明所采用的技术方案是：

[0005] 基于无源化的非线性弹簧悬架主动控制方法，包括以下步骤，

[0006] 构建非线性弹簧悬架模型；

[0007] 分析非线性弹簧悬架模型的特点，选定期望互联矩阵、阻尼矩阵及能量函数，基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架主动控制器。

[0008] 非线性弹簧悬架模型为，

[0009]

[0010] 其中，为状态变量导数，状态变量z1为车轮的位移，z2为车身的位移， 为车
轮的速度， 为车身的速度，x1＝z2，x2＝z1， p为速度分量，p＝[z2,z1]T，q为位移分量，

[0011]

[0012]

[0013] m1为车轮的重量，m2为车身的重量，c2为非线性弹簧悬架阻尼系数，k1为轮胎等价刚度系数；，u＝fa，fa为
非线性弹簧悬架主动控制力，k2为非线性弹簧悬架刚度系数，ε为弹簧非线性程度。

[0014] 非线性弹簧悬架模型的特点为，A12≠-A21， A22≠02×2。

[0015] 选定的期望互联矩阵Jd(x)为，

[0016] q12＝q22P22P11-1z1z2-1，q13＝P33-1，

[0017] q24＝P44-1，Pii＞0为可设计参数，qjj＞0为可设计参数，i＝1,2,3,4，j＝2,3,4；

[0018] 选定的阻尼矩阵Rd(x)为，

[0019] q11＝q22(P22P11-1z1z2-1)2；

[0020] 选定的能量函数Hd(x)为，

[0021] P＝diag{P11,P22,P33,P44}。

[0022] 基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架控制器为，[0023]

[0024] 即在非线性弹簧悬架控制器力的作用下，非线性弹簧悬架对于平衡点x*＝0渐近稳定的，并且可变成端口受控的哈密尔顿系统结构

[0025]

[0026] 基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架控制器的具体过程为：

[0027] 由选定的能量函数Hd(x)可知，

[0028] ▽Hd(x)＝Px

[0029] ▽2Hd(x)＝P

[0030] P＝PT＞0

[0031] 对于平衡点x*＝0，有▽H(x*)＝0，▽2H(x*)＞0，▽H(x*)为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的一阶梯度，▽2H(x*)为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的二阶梯度；

[0032] 根据选定的期望互联矩阵Jd(x)和阻尼矩阵Rd(x)，可得，

[0033] Jd(x)+JdT(x)＝0

[0034] Rd(x)＝RdT(x)≥0

[0035] 由非线性弹簧悬架模型的哈密尔顿系统结构公式、选定的期望互联矩阵Jd(x)以及选定的阻尼矩阵Rd(x)可知，

[0036]

[0037] 其中，F1(x)＝[-q11 q12 q13 0]，F2(x)＝[-q12 -q22 0q24]，F3(x)＝[-q13 0-q33 0]，F4(x)＝[0-q24 0-q44]，

[0038] 假设g(x)的左零矩阵为，

[0039]

[0040] 其中，γ1＞0，γ2＞0，γ1和γ2为可设计参数，

[0041] 则，

[0042]

[0043] 其中， ▽H(x)为非线性弹簧悬架系统非理想能量函数的一阶梯度；

[0044] 由于

[0045]

[0046]

[0047] 基于选定的期望互联矩阵Jd(x)可知，

[0048]

[0049]

[0050] 因此g⊥(x)[f(x)-F(x)▽H(x)]＝0成立，

[0051] 进而根据无源化的互联和阻尼分配控制理论可知，存在

[0052]

[0053] 使得非线性弹簧悬架对于平衡点x*＝0渐近稳定的，并且可变成端口受控的哈密尔顿系统结构；

[0054] 由于[gT(x)g(x)]-1gT(x)＝[0,0,m2,-m1]，所以

[0055]

[0056] 基于选定的期望互联矩阵Jd(x)可得，

[0057]

[0058] 则根据 可得非线性弹簧悬架控制器。

[0059] 通过选定可设计参数可对非线性弹簧悬架控制器进行简化，具体为，[0060] 选定 q33＝(m2P33)-1c2，q44＝(m1P44)-1c2时，

[0061] 非线性弹簧悬架控制器可简化为，

[0062]

[0063] 此时有

[0064]

[0065] 进一步选定P33＝P44，则

[0066]

[0067] 此时，非线性弹簧悬架可变成端口受控的哈密尔顿系统，

[0068]

[0069] 可知，q22和P44的选择决定了非线性弹簧悬架控制器特征值。

[0070] 本发明所达到的有益效果：1、本发明在考虑悬架弹簧非线性刚度特性的基础上，对1/4车辆模型描述的悬架系统进行了主动控制研究，通过分析非线性弹簧悬架系统数学模型的特点，选定合适的期望互联矩阵、阻尼矩阵及能量函数，基于无源化控制理论，获得了非线性弹簧悬架主动控制器，该主动控制器适用于恶劣地况，可以有效降低车身加速度幅值及控制量幅值，具有更快的收敛速度，并且在4-8赫兹这个人体敏感的频率区域内，车身加速度的振动分贝数也有明显下降；2、本发明通过合适选择可设计参数形式，给出了所研究悬架的无源化控制器简化形式，降低了需要整定的可设计参数个数。

附图说明

[0071] 图1为本发明的流程图；

[0072] 图2为非线性弹簧悬架模型；

[0073] 图3为弹簧的回复力和位移关系；

[0074] 图4为第一种情况地面激励位移量；

[0075] 图5为第二种情况地面激励位移量；

[0076] 图6为第一种情况车轮位移量；

[0077] 图7为第一种情况车身加速度；

[0078] 图8为第一种情况车身加速度谱密度；

[0079] 图9为第一种情况输入控制力；

[0080] 图10为第二种情况车轮位移量；

[0081] 图11为第二种情况车身加速度；

[0082] 图12为第二种情况输入控制力。

具体实施方式

[0083] 下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

[0084] 如图1所示，基于无源化的非线性弹簧悬架主动控制方法，包括以下步骤：

[0085] 步骤1，构建非线性弹簧悬架模型。

[0086] 这里采用1/4车辆模型进行座椅悬架系统的主动控制研究，同时考虑悬架弹簧的非线性刚度特性，具体模型如图2所示。其中，m1为车轮的重量，m2为车身的重量，c2为非线性弹簧悬架阻尼系数，k1为轮胎等价刚度系数，fa为非线性弹簧悬架主动控制力，k2为非线性弹簧悬架刚度系数，z1为车轮的位移，z2为车身的位移，z0表示地面位移激励。

[0087] 非线性弹簧的回复力-位移关系表示为，

[0088] fs＝k2Δz+εk2Δz3

[0089] 其中，fs为非线性弹簧的回复力，位移Δz＝z2-z1，ε为弹簧非线性程度。

[0090] 当k2＝27358，ε＝10时，弹簧的回复力和位移关系如图3所示。

[0091] 非线性弹簧悬架的动力学微分方程为，

[0092]

[0093] 取状态变量 u＝fa，x1＝z2，x2＝z1，T
p为速度分量，p＝[z2,z1] ，q为位移分量，

[0094] 假设车身加速度 车身的速度 车轮的速度 z1为车轮的位移，z2为车身的位移均可测，引入辅助控制量，

[0095]

[0096] 则式(1)可改写为，

[0097]

[0098] 其中，为状态变量导数，

[0099] 式(3)对应的标称系统模型，即非线性弹簧悬架模型为，

[0100]

[0101] 其中，

[0102]

[0103]

[0104]

[0105]

[0106] 步骤2：分析非线性弹簧悬架模型的特点，选定期望互联矩阵、阻尼矩阵及能量函数，基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架主动控制器。

[0107] 无源化控制理论在最初主要是基于无源化的能量平衡控制，然而这种方法中存在耗散性问题，因此基于无源化的互联和阻尼分配控制(IDA-PBC)被用于克服该问题，并随后得到了广泛关注。IDA-PBC方法主要是通过构造控制器将系统转变成理想的哈密尔顿系统，在设计过程中最关键的一步是求解带有可设计能量函数的偏微分方程。

[0108] 通过分析非线性弹簧悬架模型，可以发现该数学模型具有以下特点：

[0109] A12≠-A21， A22≠02×2  (6)

[0110] 与典型机械系统模型要求A12＝-A21， A22＝02×2，两者是不同的，因此，根据非线性弹簧悬架模型的特点，选取合适的期望互联矩阵、阻尼矩阵及能量函数，保证含有能量函数的偏微分方程有解，基于IDA-PBC理论，设计出新颖的非线性弹簧悬架主动控制器，实现悬架系统的振动抑制和获得良好的动态性能。

[0111] 基于无源化的互联和阻尼分配控制理论，获得了非线性弹簧悬架主动控制器的具体过程如下：

[0112] 通过构造控制器将非线性弹簧悬架系统转变成理想的哈密尔顿系统，[0113] 即

[0114] 选定的能量函数Hd(x)为：

[0115]

[0116] 其中，Pii＞0为可设计参数，i＝1,2,3,4；

[0117] 选定的期望互联矩阵Jd(x)为：

[0118] q12＝q22P22P11-1z1z2-1，q13＝P33-1，

[0119] q24＝P44-1；  (9)

[0120] 其中，qjj＞0为可设计参数，j＝2,3,4；

[0121] 选定的阻尼矩阵Rd(x)为：

[0122] q11＝q22(P22P11-1z1z2-1)2。  (10)

[0123] 证明非线性弹簧悬架控制器公式的过程如下：

[0124] 由选定的能量函数Hd(x)，即式(8)可知，

[0125] ▽Hd(x)＝Px＝[P11x1 P22x2 P33x3 P44x4]T   (11)

[0126] ▽2Hd(x)＝P，P＝PT＞0 (12)

[0127] 对于平衡点x*＝0，有

[0128] ▽H(x*)＝0，▽2H(x*)＞0  (13)

[0129] 其中，▽H(x*)为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的一阶梯度，2
▽H(x*)为平衡点x*＝0处非线性弹簧悬架系统非理想能量函数的二阶梯度；

[0130] 根据式(9)和(10)可得，

[0131] Jd(x)+JdT(x)＝0，Rd(x)＝RdT(x)≥0  (14)

[0132] 由式(7)、(9)、(10)可知，

[0133]

[0134] 其中，F1(x)＝[-q11 q12 q13 0]，F2(x)＝[-q12 -q22 0 q24]，F3(x)＝[-q13 0-q33 0]，F4(x)＝[0-q24 0-q44]，

[0135] 假设g(x)的左零矩阵为，

[0136]

[0137] 其中，γ1＞0，γ2＞0，γ1和γ2为可设计参数，

[0138] 由式(5)、(11)、(15)、(16)，可推出，

[0139]

[0140] 其中， ▽H(x)为非线性弹簧悬架系统非理想能量函数的一阶梯度；

[0141] 由于

[0142]

[0143]

[0144] 基于式(9)可知，

[0145]

[0146]

[0147] 因此，式(17)成立

[0148] g⊥(x)[f(x)-F(x)▽H(x)]＝0  (17)

[0149] 进而根据无源化的互联和阻尼分配控制理论可知，在式(12)、(13)、(14)、(17)成立时，存在

[0150]

[0151] 使得非线性弹簧悬架对于平衡点x*＝0渐近稳定的，并且可变成端口受控的哈密尔顿系统结构；

[0152] 根据式(5)，有

[0153] [gT(x)g(x)]-1gT(x)＝[0,0,m2,-m1]

[0154] 所以

[0155]

[0156] 基于式(9)可得，

[0157]

[0158] 则根据 可得非线性弹簧悬架控制器为，

[0159]

[0160] 从式(19)中可以看出，过选定可设计参数可对非线性弹簧悬架控制器进行简化，降低了需要整定的可设计参数个数，如下：

[0161] 选定 q33＝(m2P33)-1c2，q44＝(m1P44)-1c2时，

[0162] 非线性弹簧悬架控制器可简化为，

[0163]

[0164] 此时有

[0165]

[0166] 进一步选定P33＝P44，则

[0167]

[0168] 此时，非线性弹簧悬架可变成端口受控的哈密尔顿系统，

[0169]

[0170] 不难发现式(21)的特征值

[0171] 其中，

[0172] 从上述可知，仅需要确定q22和P44这两个可设计参数。

[0173] 由(21)可知，q22和P44的选择决定了非线性弹簧悬架的闭环控制系统特征值，因此，在整定q22和P44时，可以根据经典控制理论的极点配置方法来使得闭环系统具有满意的动态性能。

[0174] 为了验证上述方法，做以下仿真试验。

[0175] 某型车的相关计算参数：m1＝70，k1＝309.511，m2＝310，k2＝27.358，c2＝0.984，其中质量单位为kg，刚度系数单位为KN/m，阻尼系数单位为KN s/m，假设ε＝10。

[0176] 针对以下两种情况进行仿真：

[0177] 第一种(Case 1)：假设地面激励为某种冲击形式，地面激励位移量z0如式(21)及图4所示；

[0178]

[0179] 其中，地面质量块高度a＝0.5m，长度为l＝5m，前行速度v0＝30km/h，t表示时间。

[0180] 第二种(Case 2)：假设行驶的地面为E级地面，地面不平度功率谱密度满足Gq(n0)＝1024×10-6,n0＝0.1，前行速度v0＝7km/h，地面激励位移量图5所示。

[0181] 控制器参数选用q22＝1，q33＝0.5，q44＝20，P11＝1，P22＝1，P33＝5，P44＝0.08。

[0182] 为了体现本发明的优势，同时设计了LQR主动控制器，相应的控制器状态加权参数矩阵选择为QLQR＝diag{1,1,1,1}，控制加权参数标量选择为RLQR＝0.000001。在两种情况下的仿真结果如图6～图12所示，其中PBC表示采用的是本文设计的控制器，用实线表示，虚线表示采用LQR控制器。

[0183] 从图6～9中可以发现，针对地面冲击形式的激励，本文设计的控制器可以显著降低车身加速度幅值及控制量幅值，具有更快的收敛速度；在4-8赫兹这个人体敏感的频率区域内，车身加速度的振动分贝数也有明显下降。图10～12表明，针对E级地面这种恶劣地况，相对于LQR控制器，本文所设计的控制器可以使得轮胎位移下降，车身加速度和控制量变化更为平滑。

[0184] 表1和表2中对两种控制器作用下的轮胎位移及车身加速度的RMS值(均方根值)、振幅峰值及控制力总量给予了对比，其中，

[0185] RMS值

[0186] 振幅峰值

[0187] 控制力总量

[0188] τ为时间常数；

[0189] 表1第一种情况两种控制器作用效果对比表

[0190]

[0191] 表2第二种情况两种控制器作用效果对比表

[0192]

[0193]

[0194] 从表1和表2不难发现，除了第一种情况中，轮胎位移的最大幅值 在LQR控制器的作用下略小于PBC控制器作用时，其他各个RMS值，峰值(max值)及控制量总和值均在本文设计的PBC控制器作用下具有更小值。

[0195] 因此，本文设计的PBC控制器具有更优的总体性能。

[0196] 综上所述，上述方法在考虑悬架弹簧非线性刚度特性的基础上，对1/4车辆模型描述的悬架系统进行了主动控制研究，通过分析非线性弹簧悬架系统数学模型的特点，选定合适的期望互联矩阵、阻尼矩阵及能量函数，基于无源化控制理论，获得了非线性弹簧悬架主动控制器，该主动控制器适用于恶劣地况，可以有效降低车身加速度幅值及控制量幅值，具有更快的收敛速度，并且在4-8赫兹这个人体敏感的频率区域内，车身加速度的振动分贝数也有明显下降。

[0197] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。