CN106970524B - 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法 - Google Patents

用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法 Download PDF

Info

Publication number
CN106970524B
CN106970524B CN201710034592.2A CN201710034592A CN106970524B CN 106970524 B CN106970524 B CN 106970524B CN 201710034592 A CN201710034592 A CN 201710034592A CN 106970524 B CN106970524 B CN 106970524B
Authority
CN
China
Prior art keywords
roll
vehicle
vector
matrix
safety
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710034592.2A
Other languages
English (en)
Other versions
CN106970524A (zh
Inventor
陈士安
蔡宇萌
王建
苑磊
孙文强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University
Original Assignee
Jiangsu University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu University filed Critical Jiangsu University
Priority to CN201710034592.2A priority Critical patent/CN106970524B/zh
Publication of CN106970524A publication Critical patent/CN106970524A/zh
Application granted granted Critical
Publication of CN106970524B publication Critical patent/CN106970524B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明公开了一种用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，基于车辆三自由度转向‑侧倾运动模型的状态方程，选取横向载荷转移率、侧倾角加速度和侧倾角构建车辆侧倾运动安全综合性能指标；为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角进行满足最小相位系统的微分变形与原系统状态方程组成增广系统方程，并以新的状态变量改写车辆侧倾运动安全综合性能指标；在新的综合性能评价指标中引入包含控制项的无穷小量，以满足LQG控制器设计条件；最后，针对新的状态方程和车辆侧倾运动安全综合性能指标，依据传统LQG控制器设计方法计算出车辆主动悬架的主动控制力。

Description

用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制 器设计方法
技术领域
本发明属于车辆悬架控制领域，尤其涉及一种用于车辆主动悬架的LQG控制器参数设计方法。
背景技术
悬架是汽车重要的结构与功能部件，主要任务是将车轮受到的力与力矩传递给车身，保证汽车能稳定行驶。悬架性能对汽车行驶安全性有重要影响。传统被动悬架没有外界能量输入，参数不可调，无法较好的适应路面激励和行驶工况的变化。即使经过结构优化，悬架性能提高有限。主动悬架则弥补了被动悬架的缺陷。主动悬架具有悬架控制力发生器，可以依据路面激励和行驶工况提供主动控制力，同时兼顾悬架的行驶安全性和舒适性，提高悬架的综合性能。
线性二次型最优控制(Linear Quadratic Gaussian，LQG)可以在名义工况下，通过建立综合性能评价指标，合理设计控制器，可使线性系统获得最优性能，是较常使用的控制方法。LQG控制多应用于车辆平顺性及轮胎接地性的优化与控制。目前，侧翻事故发生几率逐年上升，提高汽车侧倾运动安全性优先级应高于平顺性，主动悬架控制器侧倾安全控制迫在眉睫。在主动悬架系统中，悬架主动控制力是LQG控制器依据二次性能指标来确定。但是存在以下几项问题：
1)在目前的车辆侧倾运动安全的相关研究中，横向载荷转移率(Lateral-loadTransfer Ratio，LTR)是最主要的安全评价指标。LTR与车身侧倾角及侧倾角加速度相互影响，仅使用LTR作为侧倾安全控制指标会导致车身侧倾角和侧倾角加速度异常变化。同时，车身侧倾加剧不仅和车身姿态有关，还和驾驶员的驾驶操作相关，在选取车辆侧倾安全评价指标时，应综合考虑车辆状态、车身姿态和驾驶员的影响。
2)求取汽车系统的三自由度转向-侧倾动力学模型和侧倾安全综合性能评价指标。将侧倾安全综合性能评价指标写成二次型标准形式
Figure GDA0001306771170000011
其中X是系统的状态向量；U是系统的控制向量；W是系统的干扰向量；Q是系统状态向量加权矩阵；R是系统控制向量加权矩阵；N是系统状态向量与系统控制向量交叉加权矩阵；P是系统干扰向量加权矩阵；L是系统状态向量与系统干扰向量交叉加权矩阵；M是系统干扰向量与系统控制向量交叉加权矩阵。T是总工作时间；t是时间变量。
依据最优控制理论，左、右悬架理想主动控制力
Figure GDA0001306771170000021
其中K为反馈增益向量矩阵，K＝R-1(SB+N)T，通过(SA)T+SA-(SB+N)R-1(SB+N)T+Q＝0可以求出矩阵S。但是发现按照传统的LQG控制器设计方法，与系统干扰量相关的P、M和L矩阵未被考虑。这意味着按照传统方法设计的控制器只能跟踪一个不准确的性能指标进行反馈控制。
3)在上述设计方法中，会导致det(R)＝0，不满足LQG控制器的设计条件det(R)>0。
由此可知：根据上述方法设计出的LQG控制器无法使主动悬架取得最佳的抗侧倾/侧翻安全性能，且有可能LQG控制器不工作。
发明内容
为解决常规LQG控制器设计方法中追踪错误性能指标和det(R)＝0导致控制器无法工作，以及侧倾安全评价指标选取的难题。本发明提出一种用于主动悬架的LQG控制器设计方法。
本发明所述的用于车辆主动悬架的LQG控制器涉及：在垂直方向上，车轮位于车身的下方外侧，车轮与簧载质量之间并联有悬架弹簧，固定阻尼减振器和悬架控制力发生器，车辆转向运动时由于惯性力产生的横向加速度，车身侧倾而使悬架产生变形；在簧载质量上设有簧载质量运动状态传感器，在车轮质量上设有车轮质量运动状态传感器，簧载质量运动状态传感器和车轮质量运动状态传感器分别通过信号线连接于LQG控制器，阻尼减振器也通过信号线连接于LQG控制器；前车轮转角以干扰向量W的形式输入主动悬架系统中；将簧载质量运动状态传感器测得的侧倾角、质心侧偏角、侧倾角速度和横摆角速度组合成状态向量X的形式；在LQG控制器中，状态向量X与干扰向量W组成扩展状态向量X1,LQG控制器根据扩展状态向量X1和侧倾安全综合性能评价指标J1求解出控制器输出信号(主动悬架主动力作动器输入信号)U1，最后将U1输送至悬架控制力发生器产生控制力FAL、FAR来控制主动悬架。
本发明采用的技术方案是：用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，包括以下步骤：步骤1，针对转向工况，对整车受力进行分析，求取车辆三自由度转向-侧倾运动系统状态方程；步骤2，综合考虑车辆行驶状态、车身姿态和驾驶员的影响，选取横向载荷转移率、侧倾角加速度和侧倾角作为LQG控制器侧倾安全评价指标，建立侧倾安全综合性能指标，基于步骤1提出的状态方程，求取侧倾安全综合性能指标二次型标准形式；步骤3，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角进行满足最小相位系统的微分变形与步骤1中的系统状态方程组成增广系统方程，获得新的综合性能评价指标；步骤4，在步骤三中新的综合性能评价指标中引入包含控制项的无穷小量，以满足LQG控制器设计条件；步骤5，针对新的状态方程和车辆侧倾运动安全综合性能指标，依据传统LQG控制器设计方法计算出车辆主动悬架的主动控制力。
进一步，所述步骤1具体过程为：
步骤1.1，车辆转向-侧倾运动动力学模型求取：依据牛顿力学，分析整车在横轴线、横摆轴线和侧倾轴线上的受力，建立三自由度汽车转向-侧倾运动动力学模型；
步骤1.2，转向-侧倾运动系统的状态方程建立：选择侧倾角θ、车身侧偏角β、侧倾角速度
Figure GDA0001306771170000031
和横摆角速度
Figure GDA00013067711700000311
组成转向-侧倾运动系统状态向量X，标记为
Figure GDA0001306771170000032
并且选择前轮转向角δv作为干扰向量W，标记为W＝[δv]，左、右悬架主动控制力FAL和FAR为控制向量U，标记为
Figure GDA0001306771170000033
建立转向-侧倾运动系统的状态方程
Figure GDA0001306771170000034
其中A是系统状态向量矩阵；B是系统控制向量矩阵；G是系统干扰向量矩阵。
进一步，所述步骤2具体过程为：
步骤2.1，侧倾安全综合性能指标的建立：选取侧倾角加速度
Figure GDA0001306771170000035
横向载荷转移率LTR和侧倾角θ作为侧倾控制评价指标，建立侧倾安全综合性能指标
Figure GDA0001306771170000036
其中：T是总工作时间；t是时间变量；δ1、δ2和δ3分别是侧倾角加速度
Figure GDA0001306771170000037
横向载荷转移率LTR和侧倾角θ的加权系数，δ1可默认为1；
步骤2.2，转向-侧倾运动系统输出方程的建立：选取侧倾角加速度
Figure GDA0001306771170000038
横向载荷转移率LTR和侧倾角θ组成系统输出状态向量Y，标记为
Figure GDA0001306771170000039
依据步骤1.2中的X、U和W矩阵，建立系统输出方程Y＝CX+DU+EW，其中C是系统输出状态向量矩阵；D是系统输出控制向量矩阵；E是系统输出干扰向量矩阵；
步骤2.3，侧倾安全综合性能指标二次型标准形式的求取：将步骤2.1中的侧倾安全综合性能指标写为二次型标准形式
Figure GDA00013067711700000310
其中Q是状态向量加权矩阵；R是控制向量加权矩阵；N是状态向量与控制向量交叉加权矩阵；P是干扰向量加权矩阵；L是状态向量与干扰向量交叉加权矩阵；M是干扰向量与控制向量交叉加权矩阵，进一步建立加权系数矩阵
Figure GDA0001306771170000041
联合步骤2.2提出的转向-侧倾运动系统输出方程，可以求得二次型标准形式中的Q、R、N、P、L和M矩阵。
进一步，所述步骤3具体过程为：
首先，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角δV，将前轮转向角 δV进行满足最小相位系统的微分变形然后，将微分变形后的前轮转向角与步骤 1.2提出的车辆三自由度转向-侧倾运动系统状态方程组成增广系统状态方程其中增广系统状态向量X1标记为：增广系统控制向量U1标 记为：增广系统干扰向量W1标记为：A1是增广系统状态向量矩阵、 B1是增广系统控制向量矩阵和G1是增广系统干扰向量矩阵，由此增广系统状态方程可以将 侧倾安全综合性能指标二次型标准形式改写为
Figure GDA0001306771170000047
进一步，所述步骤4具体过程为：
首先，为了满足LQG控制器的工作条件，在侧倾安全综合性能指标函数上增加两个和控制向量相关的无穷小量
Figure GDA0001306771170000048
和
Figure GDA0001306771170000049
然后新的侧倾安全综合性能指标二次型标准形式可以写为
Figure GDA00013067711700000410
其中Q1是增广系统状态向量加权矩阵，
Figure GDA00013067711700000411
R1是增广系统控制向量加权矩阵，
Figure GDA00013067711700000412
和N1是增广系统状态向量与增广系统控制向量交叉加权矩阵，
Figure GDA00013067711700000413
当
Figure GDA00013067711700000414
和
Figure GDA00013067711700000415
相对侧倾安全综合性能指标J为无穷小时，J1略等于J。
进一步，所述步骤5具体过程为：
依据传统的LQG控制的设计方法，可以得到主动悬架理想控制力
Figure GDA00013067711700000416
其中K1为反馈增益向量矩阵，
Figure GDA00013067711700000417
矩阵S1为黎卡提方程的解，可以通过(S1A1)T+S1A1-(S1B1+N1)R1 -1(S1B1+N1)T+Q1＝0计算得出。
本发明的有益效果是：
与由传统LQG控制方法设计的主动悬架控制器相比，通过原系统方程的增广，可使LQG控制器能跟踪前轮转向角，从而保证了控制参考的准确性；通过在综合性能指标中引入两个和控制相关的无穷小量，在保证控制目标基本不变时避免det(R)＝0，顺利设计出控制器；在设计LQG控制器时综合考虑车辆行驶状态、车身姿态和驾驶员的影响，避免了侧倾角加速度和侧倾角异常变化。通过该方法设计出的LQG控制器可使主动悬架控制系统获得最优的侧倾运动安全性。
附图说明
图1是车辆主动悬架控制原理示意图。
图2是主动悬架LQG控制器的工作原理图。
图3是车辆转向模型示意图。
图4是半车侧倾模型示意图。
图5是鱼钩工况，前轮转向角。
图6是根据本发明LQG控制器设计方法设计LQG控制器控制的主动悬架主动控制力随时间变化曲线图。
图7是根据本发明LQG控制器设计方法设计LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的横向载荷转移率LTR比较图。
图8是根据本发明LQG控制器设计方法设计LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的侧倾角加速度比较图。
图9是根据本发明LQG控制器设计方法设计LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的侧倾角比较图。
图中：1左侧悬架控制力发生器；2.左侧悬架阻尼；3.左侧悬架弹簧；4.前车轮转角传感器；5.侧倾角传感器；6.车身姿态传感器；7.簧载质量；8.右侧悬架控制力发生器；9.主动悬架LQG控制器；10.右侧车轮变形传感器；11.右侧车轮；12.右侧悬架阻尼；13.右侧悬架弹簧；14.左侧轮胎变形传感器；15.左侧车轮。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
图1为：本发明所运用于的1/2车三自由度车辆主动悬架系统为：在垂直方向上，左侧车轮15和右侧车轮11位于簧载质量7的下方。左侧车轮15与簧载质量7之间并联有左侧悬架弹簧3，左侧悬架阻尼2和左侧悬架控制力发生器1。右侧车轮11与簧载质量7之间并联有右侧悬架弹簧13，右侧悬架阻尼12和右侧悬架控制力发生器8。车辆转向运动时由于惯性力产生的横向加速度，车身侧倾而使悬架产生变形。在簧载质量7上固定设有前车轮转角传感器4、侧倾角传感器5和车身姿态传感器6。在左侧车轮15上固定设有左侧轮胎变形传感器14。在右侧车轮11上固定设有右侧轮胎变形传感器10。侧倾角加速度传感器4、侧倾角传感器5、车身姿态传感器6、左侧轮胎变形传感器14和右侧轮胎变形传感器10各自通过信号线连接于主动悬架LQG控制器9，左侧悬架控制力发生器1和右侧悬架控制力发生器8也通过信号线连接于主动悬架LQG控制器9；前车轮转角传感器4测得的前车轮转角输入到LQG控制器中，使之与侧倾角传感器5和车身姿态传感器6分别测得包括侧倾角速度、侧倾角、车身侧偏角和横摆角的状态向量X组成扩展状态向量X1，其中
Figure GDA0001306771170000061
主动悬架LQG控制器9依据增广系统方程和侧倾安全综合性能指标J1得到LQG控制器控制向量U1。
步骤1，针对转向工况，对整车受力进行分析，求取车辆三自由度转向-侧倾运动系统的状态方程。
(1)车辆转向-侧倾运动动力学模型求取。
针对图3和图4，整车质量m，车辆行驶速度v，车辆转向时产生的横向加速度ay，前轴转向角δV，车辆质心到前轴的距离lV，车辆质心到后轴的距离lH，前轴处所受侧偏力FyV，后轴处所受侧偏力FyH，簧载质量ms，左、右悬架主动控制力FAL和FAR，车身质侧倾角θ，车身质心到侧倾中心的垂直距离hs，车身绕质心纵轴的转动惯量Ix，车身绕质心横轴的转动惯量Iz。
依据牛顿第二定律，对整车在横轴线、横摆轴线和侧倾轴线上的受力分析，可以得到：
Figure GDA0001306771170000062
Figure GDA0001306771170000063
式中：β为SP处侧偏角，φ为车辆横摆角，cV和cH分别为前轴、后轴轮胎的侧偏刚度，αV和αH分别为前轴、后轴的侧偏角。αV和αH经过线性化表达如下：
Figure GDA0001306771170000071
将式(2)-(3)带入式(1)中，略去无穷小量
Figure GDA0001306771170000072
三自由度汽车转向-侧倾运动动力学模型可表达成：
Figure GDA0001306771170000073
(2)转向-侧倾运动系统的状态方程建立。
基于三自由度汽车转向-侧倾运动动力学模型，选择侧倾角θ、车身侧偏角β、侧倾角速度
Figure GDA0001306771170000074
和横摆角速度
Figure GDA0001306771170000075
组成转向-侧倾运动系统状态向量X，标记为
Figure GDA0001306771170000076
并且选择前轮转向角δv作为干扰向量W，标记为W＝[δv]，左、右悬架主动控制力FAL和FAR为控制向量U，标记为
Figure GDA0001306771170000077
车辆转向-侧倾的状态方程如下：
Figure GDA0001306771170000078
式中A是系统状态向量矩阵；B是系统控制向量矩阵；G是系统干扰向量矩阵，如下所示：
A＝[A01 A02]
Figure GDA0001306771170000079
Figure GDA00013067711700000710
Figure GDA0001306771170000081
步骤2，综合考虑车辆行驶状态、车身姿态和驾驶员的影响，选取横向载荷转移率、侧倾角加速度和侧倾角作为LQG控制器侧倾安全评价指标，建立侧倾安全综合性能指标，基于步骤1提出的状态方程求取侧倾安全综合性能指标二次型标准形式。
(1)侧倾安全综合性能指标的建立。
车身侧倾角θ与车身姿态相关。侧倾角加速度
Figure GDA0001306771170000085
与驾驶员身体上的受力相关，会影响驾驶员的心理，进而影响他对车辆的操纵。LTR是目前较常使用的车辆侧倾安全评价指标。仅使用LTR作为侧倾安全控制指标会导致车身侧倾角和侧倾角加速度异常变化。综合考虑车辆行驶状态、车身姿态和驾驶员的影响，选取侧倾角加速度、LTR和侧倾角作为侧倾控制评价指标，建立侧倾安全综合性能指标函数：
Figure GDA0001306771170000082
式中：T是总工作时间；t是时间变量；δ1、δ2和δ3分别是侧倾角加速度
Figure GDA0001306771170000086
LTR和侧倾角θ的加权系数，其中δ1可默认为1。
(2)转向-侧倾运动系统输出方程的建立。
选取侧倾角加速度、LTR和侧倾角组成系统输出状态向量Y，标记为
Figure GDA0001306771170000083
基于步骤1中提出的转向-侧倾运动系统的状态方程(5)，系统输出方程如下：
Y＝CX+DU+EW (7)
式中：
Figure GDA0001306771170000084
Figure GDA0001306771170000091
Figure GDA0001306771170000092
其中C是系统输出状态矩阵；D是系统输出控制矩阵；E是系统输出干扰矩阵。
(3)侧倾安全综合性能指标二次型标准形式的求取
基于步骤2.1中提出侧倾安全综合性能指标(6)写为二次型标准形式：
Figure GDA0001306771170000093
式中：Q是状态向量加权矩阵；R是控制向量加权矩阵；N是状态向量与控制向量交叉加权矩阵；P是干扰向量(本专利中为前轮转向角)加权矩阵；L是状态向量与干扰向量交叉加权矩阵；M是干扰向量与控制向量交叉加权矩阵。
设加权系数矩阵Z为：
Figure GDA0001306771170000094
则可以推导出：
Figure GDA0001306771170000095
步骤3，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角进行满足最小相位系统的微分变形与步骤1中的系统状态方程组成增广系统方程，获得新的综合性能评价指标。
首先，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角δV进行满足最小相位系统的微分变形
Figure GDA0001306771170000096
保证系统的可控性，式(11)必须为最小相位系统，因此必须使ξ＜0，可取ξ＝-0.001。
然后，将微分变形后的前轮转向角式(11)与步骤1.2提出的车辆三自由度转向-侧倾运动系统状态方程式(5)组成增广系统状态方程
Figure GDA0001306771170000101
式中：X1为增广系统状态向量；U1为增广系统控制向量；W1为增广系统干扰向量；A1是增广系统状态向量矩阵、B1是增广系统控制向量矩阵和G1是增广系统干扰向量矩阵。
Figure GDA0001306771170000102
Figure GDA0001306771170000103
用X1表达新的侧倾安全综合性能指标二次型标准形式有：
Figure GDA0001306771170000104
步骤4，在步骤三中新的的综合性能评价指标中引入包含控制项的无穷小量，以满足LQG控制器设计条件。
为了满足LQG控制器的工作条件，在侧倾安全综合性能指标函数上增加两个和控制向量相关的无穷小量
Figure GDA0001306771170000105
和
Figure GDA0001306771170000106
新的侧倾安全综合性能指标二次型标准形式可以写为：
Figure GDA0001306771170000107
其中Q1是增广系统状态向量加权矩阵；R1是增广系统控制向量加权矩阵；N1是增广系统状态向量与增广系统控制向量交叉加权矩阵。
Figure GDA0001306771170000108
当
Figure GDA0001306771170000109
和
Figure GDA00013067711700001010
相对J为无穷小时，J1略等于J，因此可以选取
Figure GDA00013067711700001011
步骤5，针对新的状态方程和车辆侧倾运动安全综合性能指标，依据传统LQG控制器设计方法计算出车辆主动悬架的主动控制力。
依据传统的LQG控制的设计方法，LQG控制器输出为：
U1＝-K1X1 (15)
其中：K1为反馈增益向量矩阵
Figure GDA00013067711700001012
矩阵S1为黎卡提方程的解，可以通过下式求出：
(S1A1)T+S1A1-(S1B1+N1)R1 -1(S1B1+N1)T+Q1＝0 (17)
优选实施例：
本发明的一个最优具体实施方法：
实际应用时所需的参数：整车质量m＝1146.6kg，车辆行驶速度v＝60km/h，车辆质心到前轴的距离lv＝0.88m，车辆质心到后轴的距离lh＝1.32m，簧载质量ms＝984.6kg，车身质心到侧倾中心的垂直距离hm＝0.51m，车身绕质心纵轴的转动惯量Ix＝442kg·m2，车身绕质心横轴的转动惯量Iz＝1302kg·m2，侧倾阻尼cθ＝1000N·s·m/rad；侧倾角刚度kθ＝60000N·m/rad；前轴轮胎的侧偏刚度cV＝27000N/rad；后轴轮胎的侧偏刚度cH＝20000N/rad。
在车辆上安装前轮转角传感器、侧倾角传感器、车身姿态传感器、左侧悬架控制力发生器、右侧悬架控制力发生器、主动悬架LQG控制器、右侧车轮变形传感器和左侧轮胎变形传感器。在前车轮转角传感器测得的前车轮转角输入到LQG控制器中，使之与侧倾角传感器和车身姿态传感器分别测得包括侧倾角速度、侧倾角、车身侧偏角和横摆角的状态向量X组成扩展状态向量X1，其中
Figure GDA0001306771170000111
和增广系统方程
Figure GDA0001306771170000112
主动悬架LQG控制器依据增广系统方程和侧倾安全综合性能指标J1得到LQG控制器控制向量U1。
进一步在计算侧倾安全综合性能指标J1时：
监测LTR对车身侧倾严重程度进行评估。当LTR数值过大时，提示驾驶员，车身处于不安全状态。侧倾角加速度会影响驾驶员操纵，驾驶员通过对方向盘改变车辆原有转向运动，调节车身侧倾程度。侧倾角θ与车辆的运动姿态相关。所以各侧倾评价指标之间的重要性为LTR比
Figure GDA0001306771170000113
很重要，LTR比θ重要；θ比
Figure GDA0001306771170000114
略重要。
根据层次分析法、归一法和侧倾评价指标之间的重要性确定侧倾角加速度
Figure GDA0001306771170000115
的加权系数δ1、LTR的加权系数δ2和侧倾角θ的加权系数δ3，可以得到δ1＝1、δ2＝1.2963和δ3＝7.1188。可以求得加权系数矩阵Z，可以求出侧倾安全综合性能指标J1。
如图6所示，显示了LQG控制器控制的主动悬架主动控制力随时间的变化曲线图。使用LQG控制器主动悬架的控制力左右对称，且最大值约为3200N,出现在前轮转向角速度最大的时候。主动悬架控制力的合力矩与LTR变化趋势相反，说明该控制力可以有效抑制车身侧倾趋势，提高汽车侧倾安全性。
如图7所示，显示了LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的LTR指标比较。被动悬架的LTR在仿真模拟期间最大值约为0.47，而LQG控制器控制的主动悬架的LTR相对较小，最大值约为0.4。说明LQG控制器控制的主动悬架在鱼钩工况时，可以有效抑制车身横向载荷转移率，提高车身侧倾安全性。
如图8所示，显示了LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的侧倾角加速度指标比较。LQG控制器控制的主动悬架的侧倾角加速度最大值为0.075rad/s2，小于被动悬架侧倾角加速度最大值0.48rad/s2，而且LQG控制器控制的主动悬架的侧倾角加速度曲线较为平滑。同时，LQG控制器控制的主动悬架的侧倾角加速度曲线方向与被动悬架侧倾角加速度曲线方向相反。以上分析说明LQG控制器控制的主动悬架在鱼钩工况时，可以有效减少车身侧倾变化程度，提高车身侧倾安全性。
如图9所示，显示了LQG控制器控制的主动悬架与被动悬架在鱼钩工况下的侧倾角指标比较。LQG控制器控制的主动悬架的侧倾角最大值为0.036rad，小于被动悬架侧倾角最大值0.086rad。同时，LQG控制器控制的主动悬架的侧倾角曲线方向与被动悬架侧倾角曲线方向相反
图7-9显示：与被动悬架相比，根据本发明提供的用于主动悬架的LQG控制器设计方法，可以使汽车在鱼钩工况时获得较好的侧倾/侧翻安全性能。
综上所述：本发明的一种提高车辆侧倾运动安全的主动悬架LQG控制器设计方法，属于车辆悬架控制领域。首先，基于车辆三自由度转向-侧倾运动模型的状态方程，选取横向载荷转移率、侧倾角加速度和侧倾角构建车辆侧倾运动安全综合性能指标；其次，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角进行满足最小相位系统的微分变形与原系统状态方程组成增广系统方程，并以新的状态变量改写车辆侧倾运动安全综合性能指标；进一步，在新的综合性能评价指标中引入包含控制项的无穷小量，以满足LQG控制器设计条件；最后，针对新的状态方程和车辆侧倾运动安全综合性能指标，依据传统LQG控制器设计方法计算出车辆主动悬架的主动控制力。本发明为车辆侧倾运动安全控制提供了新思路。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (4)

1.用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，其特征在于：包括以下步骤：步骤1，针对转向工况，对整车受力进行分析，求取车辆三自由度转向-侧倾运动系统状态方程；步骤2，综合考虑车辆行驶状态、车身姿态和驾驶员的影响，选取横向载荷转移率LTR、侧倾角加速度和侧倾角作为LQG控制器侧倾安全评价指标，建立侧倾安全综合性能指标，基于步骤1提出的状态方程，求取侧倾安全综合性能指标二次型标准形式；步骤3，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角，将前轮转向角进行满足最小相位系统的微分变形与步骤1中的系统状态方程组成增广系统方程，获得新的综合性能评价指标；步骤4，在步骤3中新的综合性能评价指标中引入包含控制项的无穷小量
Figure FDA0002526489060000011
和
Figure FDA0002526489060000012
当
Figure FDA0002526489060000013
和
Figure FDA0002526489060000014
相对J为无穷小时，选取
Figure FDA0002526489060000015
以满足LQG控制器设计条件；步骤5，针对新的状态方程和车辆侧倾运动安全综合性能指标，依据传统LQG控制器设计方法计算出车辆主动悬架的主动控制力；
所述步骤1具体过程为：
步骤1.1，车辆转向-侧倾运动动力学模型求取：依据牛顿力学，分析整车在横轴线、横摆轴线和侧倾轴线上的受力，建立三自由度汽车转向-侧倾运动动力学模型；
步骤1.2，转向-侧倾运动系统的状态方程建立：选择侧倾角θ、车身侧偏角β、侧倾角速度
Figure FDA0002526489060000016
和横摆角速度
Figure FDA0002526489060000017
组成转向-侧倾运动系统状态向量X，标记为
Figure FDA0002526489060000018
并且选择前轮转向角δv作为干扰向量W，标记为W＝[δv]，左、右悬架主动控制力FAL和FAR为控制向量U，标记为
Figure FDA0002526489060000019
建立转向-侧倾运动系统的状态方程
Figure FDA00025264890600000110
其中A是系统状态向量矩阵；B是系统控制向量矩阵；G是系统干扰向量矩阵；
所述步骤2具体过程为：
步骤2.1，侧倾安全综合性能指标的建立：选取侧倾角加速度
Figure FDA00025264890600000111
横向载荷转移率LTR和侧倾角θ作为侧倾控制评价指标，建立侧倾安全综合性能指标
Figure FDA00025264890600000112
其中：T是总工作时间；t是时间变量；δ1、δ2和δ3分别是侧倾角加速度
Figure FDA00025264890600000113
横向载荷转移率LTR和侧倾角θ的加权系数，δ1默认为1；
步骤2.2，转向-侧倾运动系统输出方程的建立：选取侧倾角加速度
Figure FDA0002526489060000021
横向载荷转移率LTR和侧倾角θ组成系统输出状态向量Y，标记为
Figure FDA0002526489060000022
依据步骤1.2中的X、U和W矩阵，建立系统输出方程Y＝CX+DU+EW，其中C是系统输出状态向量矩阵；D是系统输出控制向量矩阵；E是系统输出干扰向量矩阵；
步骤2.3，侧倾安全综合性能指标二次型标准形式的求取：将步骤2.1中的侧倾安全综合性能指标写为二次型标准形式
Figure FDA0002526489060000023
其中Q是状态向量加权矩阵；R是控制向量加权矩阵；N是状态向量与控制向量交叉加权矩阵；P是干扰向量加权矩阵；L是状态向量与干扰向量交叉加权矩阵；M是干扰向量与控制向量交叉加权矩阵，进一步建立加权系数矩阵
Figure FDA0002526489060000024
联合步骤2.2提出的转向-侧倾运动系统输出方程，可以求得二次型标准形式中的Q、R、N、P、L和M矩阵。
2.根据权利要求1所述的用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，其特征在于：所述步骤3具体过程为：
首先，为克服传统LQG控制器不能跟踪作为干扰量的前轮转向角δV，将前轮转向角δV进行满足最小相位系统的微分变形
Figure FDA0002526489060000025
其中，ξ＜0，取ξ＝-0.001；然后，将微分变形后的前轮转向角与步骤1.2提出的车辆三自由度转向-侧倾运动系统状态方程组成增广系统状态方程
Figure FDA0002526489060000026
其中增广系统状态向量X1标记为：
Figure FDA0002526489060000027
增广系统控制向量U1标记为：
Figure FDA0002526489060000028
增广系统干扰向量W1标记为：
Figure FDA0002526489060000029
A1是增广系统状态向量矩阵、B1是增广系统控制向量矩阵和G1是增广系统干扰向量矩阵，由此增广系统状态方程可以将侧倾安全综合性能指标二次型标准形式改写为
Figure FDA00025264890600000210
3.根据权利要求2所述的用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，其特征在于：所述步骤4具体过程为：
首先，为了满足LQG控制器的工作条件，在侧倾安全综合性能指标函数上增加两个和控制向量相关的无穷小量
Figure FDA00025264890600000211
和
Figure FDA00025264890600000212
然后新的侧倾安全综合性能指标二次型标准形式可以写为
Figure FDA0002526489060000031
其中Q1是增广系统状态向量加权矩阵，
Figure FDA0002526489060000032
R1是增广系统控制向量加权矩阵，
Figure FDA0002526489060000033
和N1是增广系统状态向量与增广系统控制向量交叉加权矩阵，
Figure FDA0002526489060000034
当
Figure FDA0002526489060000035
和
Figure FDA0002526489060000036
相对侧倾安全综合性能指标J为无穷小时，J1略等于J。
4.根据权利要求3所述的用于主动悬架的车辆侧倾运动安全线性二次型最优LQG控制器设计方法，其特征在于：所述步骤5具体过程为：
依据传统的LQG控制的设计方法，可以得到主动悬架理想控制力
Figure FDA0002526489060000037
其中K1为反馈增益向量矩阵，
Figure FDA0002526489060000038
矩阵S1为黎卡提方程的解，可以通过(S1A1)T+S1A1-(S1B1+N1)R1 -1(S1B1+N1)T+Q1＝0计算得出。
CN201710034592.2A 2017-01-17 2017-01-17 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法 Active CN106970524B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710034592.2A CN106970524B (zh) 2017-01-17 2017-01-17 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710034592.2A CN106970524B (zh) 2017-01-17 2017-01-17 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法

Publications (2)

Publication Number Publication Date
CN106970524A CN106970524A (zh) 2017-07-21
CN106970524B true CN106970524B (zh) 2020-08-28

Family

ID=59334807

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710034592.2A Active CN106970524B (zh) 2017-01-17 2017-01-17 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法

Country Status (1)

Country Link
CN (1) CN106970524B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107662468B (zh) * 2017-07-28 2020-07-31 江苏大学 用于主动悬架的车辆侧倾运动安全h2/h∞控制器设计方法
CN107976904B (zh) * 2017-11-09 2021-01-15 江苏大学 磁流变半主动悬架泰勒级数-二重h2时滞补偿控制方法
CN109849610B (zh) * 2019-02-16 2021-10-12 江苏大学 基于差动制动和主动悬架的车辆侧倾和横摆运动控制系统及其设计方法
CN110175349B (zh) * 2019-04-04 2023-07-18 江苏大学 一种基于mpc算法的独立悬架结构车辆控制方法
CN110568758B (zh) * 2019-09-12 2022-07-01 中汽研（天津）汽车工程研究院有限公司 一种自动驾驶汽车的参数自适应横向运动lqr控制方法
CN110775046B (zh) * 2019-10-22 2023-02-17 浙江万安科技股份有限公司 一种车辆防侧翻控制系统和防侧翻控制方法
CN111399380B (zh) * 2020-03-24 2021-10-22 湖南大学 一种基于Fiala刷子轮胎模型的多点预瞄LQR横向控制方法
CN112526876B (zh) * 2020-08-12 2022-09-13 江苏大学 一种基于数据驱动的lpv系统的lqg控制器的设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102756686A (zh) * 2011-04-26 2012-10-31 贵州大学 预警车辆侧翻的方法和装置
CN104865079A (zh) * 2015-06-14 2015-08-26 山东交通学院 一种汽车侧倾综合评价指标测量装置及其计算方法
CN106004870A (zh) * 2016-06-23 2016-10-12 吉林大学 一种基于变权重模型预测算法的车辆稳定性集成控制方法
CN106080553A (zh) * 2016-07-13 2016-11-09 南京航空航天大学 一种融合车速变化的四轮转向汽车防侧翻控制系统及方法
CN106197416A (zh) * 2016-07-04 2016-12-07 南京航空航天大学 一种多工况汽车侧翻指标计算装置及其计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102756686A (zh) * 2011-04-26 2012-10-31 贵州大学 预警车辆侧翻的方法和装置
CN104865079A (zh) * 2015-06-14 2015-08-26 山东交通学院 一种汽车侧倾综合评价指标测量装置及其计算方法
CN106004870A (zh) * 2016-06-23 2016-10-12 吉林大学 一种基于变权重模型预测算法的车辆稳定性集成控制方法
CN106197416A (zh) * 2016-07-04 2016-12-07 南京航空航天大学 一种多工况汽车侧翻指标计算装置及其计算方法
CN106080553A (zh) * 2016-07-13 2016-11-09 南京航空航天大学 一种融合车速变化的四轮转向汽车防侧翻控制系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Design of a robust controller for rollover prevention with active suspension and differential braking;Seongjin Yim;《Journal of Mechanical Science and Technology》;20121231;第26卷(第1期);213-222 *
二次型最优控制的半挂汽车列车主动侧倾控制算法研究;宗长富 等;《中国机械工程》;20080430;第19卷(第07期);872-877 *
基于混合H_2/H_∞控制的电动助力系统转向路感;赵万忠 等;《机械工程学报》;20090430;第45卷(第04期);142-147 *

Also Published As

Publication number Publication date
CN106970524A (zh) 2017-07-21

Similar Documents

Publication Publication Date Title
CN106970524B (zh) 用于主动悬架的车辆侧倾运动安全线性二次型最优lqg控制器设计方法
CN104385873B (zh) 一种汽车悬架系统的多目标优化方法
Mastinu et al. Integrated controls of lateral vehicle dynamics
CN107662468B (zh) 用于主动悬架的车辆侧倾运动安全h2/h∞控制器设计方法
Mahmoodi-Kaleibar et al. Optimization of suspension system of off-road vehicle for vehicle performance improvement
CN108891221A (zh) 一种基于模态能量分配法的主动悬架系统及其工作方法
CN111391595B (zh) 车辆防侧翻主动倾摆模型预测控制方法
CN112339517A (zh) 一种半主动悬架控制方法及控制系统
Chen et al. A novel LQG controller of active suspension system for vehicle roll safety
Javanshir et al. Optimization of suspension system of heavy off-road vehicle for stability enhancement using integrated anti-roll bar and coiling spring mechanism
JP6543393B1 (ja) ステアリング制御装置及びステアリング装置
CN108984888A (zh) 基于敏感性分析的麦弗逊悬架多目标优化函数构建方法
CN113320525A (zh) 在沿弯道行驶时控制具有可变刚度和转向式后轮的道路车辆的方法
Liang et al. Integration of active tilting control and full-wheel steering control system on vehicle lateral performance
Her et al. Development of integrated control of electronic stability control, continuous damping control and active anti-roll bar for vehicle yaw stability
Hyvärinen The improvement of full vehicle semi-active suspension through kinematical model
CN113071282B (zh) 一种车辆调节方法及装置
CN109849610B (zh) 基于差动制动和主动悬架的车辆侧倾和横摆运动控制系统及其设计方法
Lee et al. Enhancement of vehicle handling characteristics by suspension kinematic control
CN112590483B (zh) 一种基于观测器的汽车横向稳定性与主动悬架多目标集成控制方法
CN111338208B (zh) 一种车辆侧倾和平顺性协调控制方法
Ukamnal et al. Design of Trailing Arm Suspension
Shim et al. Suspension design and dynamic analysis of a lightweight vehicle
Biglarbegian et al. Design of a novel fuzzy controller to enhance stability of vehicles
Mihon et al. The behaviour of a vehicle’s suspension system on dynamic testing conditions

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant