CN107599778A - 一种电磁混合悬架的控制方法 - Google Patents

一种电磁混合悬架的控制方法 Download PDF

Info

Publication number
CN107599778A
CN107599778A CN201710637762.6A CN201710637762A CN107599778A CN 107599778 A CN107599778 A CN 107599778A CN 201710637762 A CN201710637762 A CN 201710637762A CN 107599778 A CN107599778 A CN 107599778A
Authority
CN
China
Prior art keywords
suspension
control
lqg
linear electric
electromagnetism
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.)
Granted
Application number
CN201710637762.6A
Other languages
English (en)
Other versions
CN107599778B (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 CN201710637762.6A priority Critical patent/CN107599778B/zh
Priority to PCT/CN2017/100109 priority patent/WO2019024164A1/zh
Priority to US16/491,879 priority patent/US11584184B2/en
Publication of CN107599778A publication Critical patent/CN107599778A/zh
Application granted granted Critical
Publication of CN107599778B publication Critical patent/CN107599778B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G17/00—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load
    • B60G17/015—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G17/00—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load
    • B60G17/015—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements
    • B60G17/018—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements characterised by the use of a specific signal treatment or control method
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G17/00—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load
    • B60G17/015—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements
    • B60G17/0195—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements characterised by the regulation being combined with other vehicle control systems
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G17/00—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load
    • B60G17/06—Characteristics of dampers, e.g. mechanical dampers
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2202/00—Indexing codes relating to the type of spring, damper or actuator
    • B60G2202/10—Type of spring
    • B60G2202/16—Magnetic spring
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2202/00—Indexing codes relating to the type of spring, damper or actuator
    • B60G2202/20—Type of damper
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2400/00—Indexing codes relating to detected, measured or calculated conditions or factors
    • B60G2400/10—Acceleration; Deceleration
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2400/00—Indexing codes relating to detected, measured or calculated conditions or factors
    • B60G2400/50—Pressure
    • B60G2400/52—Pressure in tyre
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2400/00—Indexing codes relating to detected, measured or calculated conditions or factors
    • B60G2400/60—Load
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60G—VEHICLE SUSPENSION ARRANGEMENTS
    • B60G2600/00—Indexing codes relating to particular elements, systems or processes used on suspension systems or suspension control systems
    • B60G2600/18—Automatic control means
    • B60G2600/184—Semi-Active control means

Abstract

本发明公开了一种电磁混合悬架的控制方法，将电磁混合悬架分为四种模式：舒适模式、运动模式、综合模式、馈能模式，四种模式间的切换由驾驶员决定。在舒适、运动、综合模式下，采用混合控制，分为两种子模式：主动控制模式、半主动控制模式，两种子模式间的切换条件为引入新参数Cact，与直线电机最大等效电磁阻尼系数Ceqmax比较判断。本发明解决了悬架舒适性与轮胎接地性性的协调问题，由人为决定模式切换满足了不同工况及不同用户的需求。同时混合控制方法解决了主动悬架能耗过高、半主动悬架控制效果有限的问题，在保证汽车悬架动力学性能的同时，能够减少能量消耗。同时馈能模式的设计使得悬架能够进行能量回收，符合节能减排要求。

Description

一种电磁混合悬架的控制方法
技术领域
本发明涉及汽车底盘控制领域，特指一种电磁混合悬架的控制方法。
背景技术
我国汽车需求量和保有量都在不断增加，由此所带来的能源紧张与环境问题日益突出，节能减排成为当今汽车发展的主题之一。汽车行驶过程中存在大量的能量损耗，其中振动能量损耗通过阻尼器转化为热能损耗掉，占所有能量损耗约20％。
电磁混合悬架可以通过直线电机回收部分振动能量提高汽车的燃油经济性，也可以通过直线电机主动输出控制力来改善汽车悬架的动力学性能。在电磁混合悬架进行主动控制时，悬架将消耗能量，且悬架的隔振性与轮胎接地性存在相互矛盾的关系。
因此，设计合适的控制策略来对电磁混合悬架馈能模式与动力学性能模式进行切换，并且在主动模式下协调悬架的能耗特性、隔振性、轮胎接地性具有很高的研究价值。
发明内容
本发明所要解决的技术问题是：设计一种电磁混合悬架的控制方法，能够对电磁混合悬架馈能模式与动力学性能模式进行切换，并且在悬架处于动力学性能模式的情况下，对悬架的能耗特性、隔振性以及轮胎接地性进行协调控制。
本发明解决该技术问题所采用的技术方案是：
一种电磁混合悬架的控制方法，按以下步骤实现：
A建立四分之一电磁混合悬架系统模型，线性数学模型可用动力学微分方程表示：
其中，mb为簧载质量，mw为非簧载质量，ks为弹簧刚度，Cs为被动可调阻尼器阻尼系数，kt为轮胎等效刚度。Z0为路面输入位移，数学模型可表示为：
其中，G0为路面不平度；P为π，u为车速；f0为下截止频率；w(t)为高斯白噪声且均值为0，用来生成随机路面。本文设置路面等级为C级，路面不平度G0为256·10-6m3，车速设置为20m/s，
B将电磁混合悬架分为四种工作模式，分别为舒适、运动、综合、馈能模式，四种模式间的切换由人为决定。当电磁混合悬架处于舒适模式时，以车身加速度均方根值作为悬架的主要改善目标；处于运动模式时，以轮胎动载荷均方根值作为悬架的主要改善目标；处于综合模式时，要求悬架能兼顾车身加速度和轮胎动载荷均方根值这两个性能指标；处于馈能模式时，要求悬架在保证基本的动力学性能的前提下，实现馈能最大化。
C在电磁混合悬架处于舒适、运动、综合模式时，分为两种子模式：主动控制模式、半主动控制模式。
在主动控制模式下，ECU通过控制策略计算出所需控制力Fdes，除以电机推力系数kf得到直线电机所需电流Ides，输入至直线电机控制器，车载电源与直线电机控制器相连接，直线电机控制器与直线电机相连接，输入电流I至直线电机，直线电机向悬架输出主动控制力Fact。此时悬架消耗能量；
在半主动控制模式下，ECU通过控制策略计算出所需控制力Fdes，除以电机推力系数kf得到直线电机所需电流Ides，输入至半主动回路控制器，半主动回路控制器与半主动控制回路相连接，向半主动控制回路输入占空比，调节半主动控制回路中的电流I，从而调节直线电机的等效阻尼系数Ceq，半主动控制回路与直线电机相连接，直线电机向悬架系统输出电磁阻尼力，直线电机用作电磁阻尼器即发电机，直线电机作为发电机回收的电能储存在超级电容中。引入参数：
式中，Fdes表示控制策略所计算出的直线电机所需提供的控制力，将Cact与直线电机最大等效阻尼系数Ceqmax比较，当0<Cact<Ceqmax时，处于半主动控制模式，其他情况处于主动控制模式。
D在电磁混合悬架处于馈能模式时，直线电机用作发电机，相当于一个不可调电磁阻尼器在悬架中作用，电磁阻尼系数为最大等效电磁阻尼系数Ceqmax，以保证能回收尽可能多的电能，电能储存在超级电容中。
进一步，步骤C中所述的直线电机所需提供的控制力Fdes由LQG控制策略计算得到，通过改变LQG加权系数的选择，分别对应舒适、运动、综合模式。将LQG性能指标函数设计为：
其中，Zw-Z0为轮胎动位移，将其与轮胎等效刚度kt相乘则为轮胎动载荷，Zb-Zw为悬架动行程，为车身加速度，所以与之对应的，q1、q2、q3为轮胎动载荷、悬架动行程、车身加速度的加权系数。
电磁混合悬架的状态空间方程为：
其中，选取为状态变量，U＝[Fact]为输入控制变量，可得：
选取输出变量输出方程如下：
Y＝CX+DU (6)
对于LQG控制算法，目标为求得输入控制变量U＝[Fact]，使性能指标函数取得最小值，属于二次型最优控制问题中的一种，故将性能指标函数写作标准二次型形式如下：
由式(4)、式(6)可得：
其中，Q0为LQG加权系数矩阵，表达式如下：
因此可以求得状态变量的加权矩阵Q、控制变量的加权矩阵R以及交叉项的权重矩阵N：
可以看出状态变量矩阵Q的大小与q1、q2、q3的取值大小有关，所以加权系数对于LQG控制策略的控制效果起决定性作用。
要得到LQG控制策略的最优控制力Fact，需要求得反馈状态增益矩阵K，即：
Fact＝-KX (9)
也可表示为：
根据LQG控制理论，增益矩阵K可表示为：
K＝R-1(SB+N)T (11)
其中，矩阵S为Ricatti方程的解，Ricatti方程的形式如下：
(SA)T+SA-(SB+N)R-1(SB+N)T+Q＝0 (12)
在软件Matlab中，可利用其自带的LQR函数求得矩阵K、S：
[K S E]＝lqr(A,B,Q,R,N) (13)
进一步，通过采集车身加速度、轮胎动载荷、悬架动行程三个状态变量输入至LQG控制器中，在不同混合模式下具有不同的LQG加权系数组合，由状态变量及加权系数计算出Fdes。
进一步，LQG控制策略加权系数由遗传算法优化得到，优化目标分别为车身加速度和轮胎动载荷，代表了车辆隔振性与轮胎接地性，而悬架动行程这一性能参数用来作为优化算法的约束条件。优化的设计变量为LQG加权系数：
x＝(q1 q2 q3)T (14)
遗传算法优化LQG加权系数的流程中，需要通过目标函数计算得到适应度函数。电磁混合悬架运用LQG控制策略时具有三种模式：舒适、运动、综合，其中舒适模式要求车身加速度最小，运动模式要求轮胎动载荷最小，而综合模式要求兼顾两者。所以，三种模式下进行遗传算法优化时，适应度函数的计算方法均不同。
舒适模式下，目标函数为车身加速度均方根值，由于要求车身加速度均方根值取最小值，而适应度函数应与遗传概率成正比，要求取最大值，所以适应度函数设计为：
其中，AP为相同仿真条件下传统被动悬架的车身加速度均方根值，AH为电磁混合悬架LQG控制策略下的车身加速度均方根值。
运动模式下，适应度函数设计为：
其中，TP为相同仿真条件下传统被动悬架的轮胎动载荷均方根值，TH为电磁混合悬架LQG控制策略下的轮胎动载荷均方根值。
综合模式下，目标函数为车身加速度均方根值和轮胎动载荷均方根值两个，所以属于多目标遗传算法优化，适应度函数设计为：
具体步骤为：
(1)参数初始化。将q1、q2、q3分别进行二进制编码，然后串联构成染色体串，形成一个个体。进行解码时将二进制编码转换为十进制数值，进行悬架系统控制。
(2)生成种群。首先要通过随机方法来产生种群大小M＝80的初始种群，而新生种群则通过由上一代种群经过遗传算子操作后生成。
(3)适应度函数评估。将解码后的加权系数所控制的悬架系统的目标值函数AH、TH提取，转换为适应度函数F(x)。
(4)选择操作。采用比例选择方法，设第i个个体的适应度为F(i)，则遗传概率psi为：
(5)交叉操作。采用单点交叉，依照交叉概率pc＝0.6对每一对相互配对的个体，选定单个交叉点，交换其部分染色体。
(6)变异操作。采用基本位变异，以变异概率pm＝0.001对个体染色体的编码串中某一部分的基因值进行取反运算。
(7)终止条件判断。设定终止代数T＝300，将第300进化代数群体中的最优个体作为最优解输出。若连续15代个体的平均适应度差异小于0.5％则提前结束运算，当前群体中的最优个体作为最优解输出。
进一步，步骤C中所述主动控制模式、半主动控制模式之间的电路切换通过开关器件IGBT实现：车载电源与逆变器相连接，逆变器另一端与开关器件VT1相连接，开关器件VT1另一端连接开关器件VT2和直线电机M，开关器件VT2与整流桥、DC/DC变换器相连接，DC/DC变换器输出与二极管正端连接，二极管负端连接超级电容SC。
进一步，步骤D中所述的馈能模式中，设计控制电路，直线电机M与整流桥相连，整流桥一端与电感L一端相连，电感L另一端分别与开关S1一端和开关S2一端相连，开关S2另一端与二极管正端相连，二极管负端与超级电容SC一端相连，超级电容SC另一端分别和开关S1的另一端、整流桥的另一端相连。
进一步，可调减振阻尼器分为四档可调，分别对应舒适、运动、综合、馈能四种模式。
本技术具有如下优点：
(1)将电磁混合悬架分为舒适、运动、综合、馈能四种模式，满足了汽车在不同工况下的需求，由人为决定切换，满足用户自身需求。
(2)在舒适、运动、综合模式下分为两种子模式：主动控制模式、半主动控制模式，并且设计了两种模式下的切换条件，该方法能够结合主动控制控制效果好、半主动控制能量消耗低的优点，能够兼顾悬架的动力学性能以及能耗特性。
(3)在馈能模式下，悬架能够进行能量回收，且悬架此时相当于传统被动悬架运作，保证基本的悬架功能的同时，符合节能减排要求。
附图说明
图1为本发明总体流程图；
图2为遗传算法优化流程图；
图3为主动控制子模式系统框图；
图4为半主动控制子模式系统框图；
图5为半主动控制回路拓扑图；
图6为切换电路框图；
图7为一般主动LQG控制与混合LQG控制电机顺时功率时域仿真对比；
图8为主动悬架与电磁混合悬架10s内能量消耗；
图9为馈能模式充电电路图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明的一种电磁混合悬架的控制方法，总体流程如图1所示。
首先建立四分之一电磁混合悬架系统模型，线性数学模型可用动力学微分方程表示：
其中，mb为簧载质量，mw为非簧载质量，ks为弹簧刚度，Cs为被动可调阻尼器阻尼系数，kt为轮胎等效刚度。Z0为路面输入位移，数学模型可表示为：
其中，G0为路面不平度；u为车速；f0为下截止频率；w(t)为高斯白噪声且均值为0，用来生成随机路面。本文设置路面等级为C级，路面不平度G0为256·10-6m3，车速设置为20m/s。
本发明将电磁混合悬架分为四种模式：舒适、运动、综合、馈能模式，四种模式间的切换由驾驶员决定。每一种模式下可调减振器的阻尼系数不同，将可调减振器设计为四档可调，在驾驶员进行模式切换时，控制器控制步进电机调节减振器节流口大小，从而调节减振器阻尼系数Cs。
在电磁混合悬架处于舒适、运动、综合模式时，通过LQG控制策略来对悬架进行控制，LQG性能指标函数设计为：
式中，q1、q2、q3分别为轮胎动载荷、悬架动行程、车身加速度的加权系数，不同的加权系数选取对悬架动力学性能会产生影响，本发明要求在舒适模式下侧重车身加速度的改善，运动模式下侧重轮胎动载荷的改善，综合模式下兼顾车身加速度与轮胎动载荷，悬架动行程指标起限位作用。经过遗传算法优化，优化目标为车身加速度、轮胎动载荷，悬架动行程作为优化约束条件，优化流程如图2所示，三种模式下加权系数优化后如下表：
由于目前悬架的主动控制存在消耗能量大的缺点，而半主动控制存在控制效果有限的缺点，所以本发明提出了一种混合控制方法，在电磁混合悬架处于舒适、运动、综合模式时，设计两种子模式：主动控制模式、半主动控制模式。
在主动控制模式下，由车载电源向直线电机供电，主动输出控制力，实际是控制直线电机的输入电流I，本文设计双环控制系统，如图3所示，从悬架系统中测得的状态反馈变量输入至外控制器，由控制器计算得到电机所需提供的控制力Fdes，进而得到所需电流Ides输入至内环控制，内环控制为电流跟踪控制，使得直线电机的输入电流I能够跟踪所需电流Ides，从而输出控制力Fact，外部能量来源于车载电源的直接供电。
在半主动控制模式下，控制系统如图4所示，通过从悬架系统中反馈的状态变量计算出所需要的等效电磁阻尼力Fdes，计算得到此时直线电机绕组中流过的所需电流Ides，半主动回路控制器通过调节半主动控制回路中的占空比来使直线电机绕组中的实际电流Ireal跟踪Ides，从而使直线电机提供的等效电磁阻尼力Fact跟踪Fdes。图5为半主动控制回路拓扑结构图，直线电机M与整流桥连接，整流桥连接开关S1的漏极，开关S1的源极分别连接电感L和二极管D1的负端，电感L的另一端连接开关S2的漏极和二极管D2的正端，二极管D2的负端连接超级电容SC，超级电容SC的另一端、二极管D1的正端、开关S2的源级相连，当直线电机端电压小于超级电容端电压时，S1导通，S2断开，使得回路中产生电流，避免出现死区现象，其他情况S1断开，S2导通，电能向超级电容中储存；由于悬架的上下运动会使直线电机中电流的方向不断发生变化，所以需要在回路中添加整流桥来保证电流方向一致性。通过DC-DC变换器起到升降压的作用，半主动回路控制器经过运算向S1、S2输出两路脉冲信号在调节占空比。在半主动控制中，直线电机产生的电能能够通过半主动控制回路后储存至超级电容中，起到回收能量的作用。
在这里，引入参数：
式中，Fdes表示控制策略所计算出的直线电机所需提供的控制力，将Cact与直线电机最大等效阻尼系数Ceqmax比较，当0<Cact<Ceqmax时，处于半主动控制模式，其他情况处于主动控制模式。两种子模式间的切换由开关器件IGBT实现，开关器件的频率远大于悬架振动频率，如图6所示。
上述Fdes由LQG控制策略计算的具体过程为：
电磁混合悬架的状态空间方程为：
其中，选取为状态变量，U＝[Fact]为输入控制变量，可得：
选取输出变量输出方程如下：
Y＝CX+DU (6)
对于LQG控制算法，目标为求得输入控制变量U＝[Fact]，使性能指标函数取得最小值，属于二次型最优控制问题中的一种，故将性能指标函数写作标准二次型形式如下：
由式(4)、式(6)可得：
其中，Q0为LQG加权系数矩阵，表达式如下：
因此可以求得状态变量的加权矩阵Q、控制变量的加权矩阵R以及交叉项的权重矩阵N：
可以看出状态变量矩阵Q的大小与q1、q2、q3的取值大小有关，所以加权系数对于LQG控制策略的控制效果起决定性作用；
要得到LQG控制策略的最优控制力Fact，需要求得反馈状态增益矩阵K，即：
Fact＝-KX (9)
也可表示为：
根据LQG控制理论，增益矩阵K可表示为：
K＝R-1(SB+N)T (11)
其中，矩阵S为Ricatti方程的解，Ricatti方程的形式如下：
(SA)T+SA-(SB+N)R-1(SB+N)T+Q＝0 (12)
在软件Matlab中，可利用其自带的LQR函数求得矩阵K、S：
[K S E]＝lqr(A,B,Q,R,N) (13)
LQG控制策略加权系数由遗传算法优化得到，优化目标分别为车身加速度和轮胎动载荷，代表了车辆隔振性与轮胎接地性，而悬架动行程这一性能参数用来作为优化算法的约束条件，优化的设计变量为LQG加权系数：
x＝(q1 q2 q3)T (14)
遗传算法优化LQG加权系数的流程中，需要通过目标函数计算得到适应度函数，电磁混合悬架运用LQG控制策略时具有三种模式：舒适、运动、综合，其中舒适模式要求车身加速度最小，运动模式要求轮胎动载荷最小，而综合模式要求兼顾两者，所以，三种模式下进行遗传算法优化时，适应度函数的计算方法均不同；
舒适模式下，目标函数为车身加速度均方根值，由于要求车身加速度均方根值取最小值，而适应度函数应与遗传概率成正比，要求取最大值，所以适应度函数设计为：
其中，AP为相同仿真条件下传统被动悬架的车身加速度均方根值，AH为电磁混合悬架LQG控制策略下的车身加速度均方根值；
运动模式下，适应度函数设计为：
其中，TP为相同仿真条件下传统被动悬架的轮胎动载荷均方根值，TH为电磁混合悬架LQG控制策略下的轮胎动载荷均方根值；
综合模式下，目标函数为车身加速度均方根值和轮胎动载荷均方根值两个，所以属于多目标遗传算法优化，适应度函数设计为：
上述遗传算法优化LQG控制策略加权系数的具体步骤为：
(1)参数初始化：将q1、q2、q3分别进行二进制编码，然后串联构成染色体串，形成一个个体，进行解码时将二进制编码转换为十进制数值，进行悬架系统控制；
(2)生成种群：首先要通过随机方法来产生种群大小M＝80的初始种群，而新生种群则通过由上一代种群经过遗传算子操作后生成；
(3)适应度函数评估：将解码后的加权系数所控制的悬架系统的目标值函数AH、TH提取，转换为适应度函数F(x)；
(4)选择操作：采用比例选择方法，设第i个个体的适应度为F(i)，则遗传概率psi为：
(5)交叉操作：采用单点交叉，依照交叉概率pc＝0.6对每一对相互配对的个体，选定单个交叉点，交换其部分染色体；
(6)变异操作：采用基本位变异，以变异概率pm＝0.001对个体染色体的编码串中某一部分的基因值进行取反运算；
(7)终止条件判断：设定终止代数T＝300，将第300进化代数群体中的最优个体作为最优解输出，若连续15代个体的平均适应度差异小于0.5％则提前结束运算，当前群体中的最优个体作为最优解输出。
此种混合控制方法的优点在于：能量消耗相比于主动控制的悬架有明显降低，且控制效果与主动悬架相差不大。如图7所示，为一般主动LQG控制与混合LQG控制电机顺时功率时域仿真对比，可见在大部分时间下，两者的电机瞬时功率相同，且主动控制所有功率皆为电机消耗功率，而混合控制在0线以下部分为半主动控制，这些功率为电机回收能量功率。如图8所示，为10s内能量消耗情况，可见混合控制相比与主动控制能量消耗减少了约20％。
在电磁混合悬架处于馈能模式时，要求直线电机能够回收尽可能多的振动能量，所以直线电机用作电磁阻尼器，且等效阻尼系数为直线电机的最大等效阻尼系数Ceqmax。馈能模式能量回收充电电路如图9所示。当直线电机输出电压Um小于超级电容端电压Uc时，回路中将没有感应电流产生，出现死区现象，影响悬架性能，所以此时控制开关S1导通，S2断开，此时回收的振动能量储存在电感L中，当直线电机输出电压Um大于超级电容端电压Uc时，控制S1断开，S2导通，向超级电容充电。
应理解上述施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

Claims (9)

1.一种电磁混合悬架的控制方法，其特征在于：包括以下步骤：
A：建立四分之一电磁混合悬架系统模型，线性数学模型可用动力学微分方程表示：
其中，mb为簧载质量，mw为非簧载质量，ks为弹簧刚度，Cs为被动可调阻尼器阻尼系数，kt为轮胎等效刚度，Z0为路面输入位移，数学模型可表示为：
其中，G0为路面不平度；P为π，u为车速；f0为下截止频率；w(t)为高斯白噪声且均值为0，用来生成随机路面；
B：将电磁混合悬架分为四种工作模式，分别为舒适、运动、综合、馈能模式，当电磁混合悬架处于舒适模式时，以车身加速度均方根值作为悬架的主要改善目标；处于运动模式时，以轮胎动载荷均方根值作为悬架的主要改善目标；处于综合模式时，要求悬架能兼顾车身加速度和轮胎动载荷均方根值这两个性能指标；处于馈能模式时，要求悬架在保证基本的动力学性能的前提下，实现馈能最大化；
C：在电磁混合悬架处于舒适、运动、综合模式时，分为两种子模式：主动控制模式、半主动控制模式：
在主动控制模式下，ECU通过控制策略计算出所需控制力Fdes，除以电机推力系数kf得到直线电机所需电流Ides，输入至直线电机控制器，车载电源与直线电机控制器相连接，直线电机控制器与直线电机相连接，输入电流I至直线电机，直线电机向悬架输出主动控制力Fact，此时悬架消耗能量；
在半主动控制模式下，ECU通过控制策略计算出所需控制力Fdes，除以电机推力系数kf得到直线电机所需电流Ides，输入至半主动回路控制器，半主动回路控制器与半主动控制回路相连接，向半主动控制回路输入占空比，调节半主动控制回路中的电流I，从而调节直线电机的等效阻尼系数Ceq，半主动控制回路与直线电机相连接，直线电机向悬架系统输出电磁阻尼力，直线电机用作电磁阻尼器即发电机，直线电机作为发电机回收的电能储存在超级电容中，引入参数：
式中，Fdes表示控制策略所计算出的直线电机所需提供的控制力，将Cact与直线电机最大等效阻尼系数Ceqmax比较，当0<Cact<Ceqmax时，处于半主动控制模式，其他情况处于主动控制模式；
D：在电磁混合悬架处于馈能模式时，直线电机用作发电机，相当于一个不可调电磁阻尼器在悬架中作用，电磁阻尼系数为最大等效电磁阻尼系数Ceqmax，以保证能回收尽可能多的电能，电能储存在超级电容中。
2.根据权利要求1所述的一种电磁混合悬架的控制方法，其特征在于，步骤C中所述控制力Fdes由LQG控制策略计算得到，通过改变LQG加权系数的选择，分别对应舒适、运动、综合模式；将LQG性能指标函数设计为：
其中，Zw-Z0为轮胎动位移，将其与轮胎等效刚度kt相乘则为轮胎动载荷，Zb-Zw为悬架动行程，为车身加速度，所以与之对应的q1、q2、q3为轮胎动载荷、悬架动行程、车身加速度的加权系数。
3.根据权利要求2所述的一种电磁混合悬架的控制方法，其特征在于，所述Fdes由LQG控制策略计算的具体过程为：
电磁混合悬架的状态空间方程为：
其中，选取为状态变量，U＝[Fact]为输入控制变量，可得：
选取输出变量输出方程如下：
Y＝CX+DU (6)
对于LQG控制算法，目标为求得输入控制变量U＝[Fact]，使性能指标函数取得最小值，属于二次型最优控制问题中的一种，故将性能指标函数写作标准二次型形式如下：
由式(4)、式(6)可得：
其中，Q0为LQG加权系数矩阵，表达式如下：
因此可以求得状态变量的加权矩阵Q、控制变量的加权矩阵R以及交叉项的权重矩阵N：
可以看出状态变量矩阵Q的大小与q1、q2、q3的取值大小有关，所以加权系数对于LQG 控制策略的控制效果起决定性作用；
要得到LQG控制策略的最优控制力Fact，需要求得反馈状态增益矩阵K，即：
Fact＝-KX (9)
也可表示为：
根据LQG控制理论，增益矩阵K可表示为：
K＝R-1(SB+N)T (11)
其中，矩阵S为Ricatti方程的解，Ricatti方程的形式如下：
(SA)T+SA-(SB+N)R-1(SB+N)T+Q＝0 (12)
在软件Matlab中，可利用其自带的LQR函数求得矩阵K、S：
[K S E]＝lqr(A,B,Q,R,N) (13)。
4.根据权利要求3所述的一种电磁混合悬架的控制方法，其特征在于，通过采集车身加速度、轮胎动载荷、悬架动行程三个状态变量输入至LQG控制器中，在不同混合模式下具有不同的LQG加权系数组合，由状态变量及加权系数计算出Fdes。
5.根据权利要求2所述的一种电磁混合悬架的控制方法，其特征在于，LQG控制策略加权系数由遗传算法优化得到，优化目标分别为车身加速度和轮胎动载荷，代表了车辆隔振性与轮胎接地性，而悬架动行程这一性能参数用来作为优化算法的约束条件，优化的设计变量为LQG加权系数：
x＝(q1 q2 q3)T (14)
遗传算法优化LQG加权系数的流程中，需要通过目标函数计算得到适应度函数，电磁混合悬架运用LQG控制策略时具有三种模式：舒适、运动、综合，其中舒适模式要求车身加速度最小，运动模式要求轮胎动载荷最小，而综合模式要求兼顾两者，所以，三种模式下进行遗传算法优化时，适应度函数的计算方法均不同；
舒适模式下，目标函数为车身加速度均方根值，由于要求车身加速度均方根值取最小值，而适应度函数应与遗传概率成正比，要求取最大值，所以适应度函数设计为：
其中，AP为相同仿真条件下传统被动悬架的车身加速度均方根值，AH为电磁混合悬架LQG控制策略下的车身加速度均方根值；
运动模式下，适应度函数设计为：
其中，TP为相同仿真条件下传统被动悬架的轮胎动载荷均方根值，TH为电磁混合悬架LQG控制策略下的轮胎动载荷均方根值；
综合模式下，目标函数为车身加速度均方根值和轮胎动载荷均方根值两个，所以属于多目标遗传算法优化，适应度函数设计为：
。
6.根据权利要求5所述的一种电磁混合悬架的控制方法，其特征在于，遗传算法优化LQG控制策略加权系数的具体步骤为：
(1)参数初始化：将q1、q2、q3分别进行二进制编码，然后串联构成染色体串，形成一个个体，进行解码时将二进制编码转换为十进制数值，进行悬架系统控制；
(2)生成种群：首先要通过随机方法来产生种群大小M＝80的初始种群，而新生种群则通过由上一代种群经过遗传算子操作后生成；
(3)适应度函数评估：将解码后的加权系数所控制的悬架系统的目标值函数AH、TH提取，转换为适应度函数F(x)；
(4)选择操作：采用比例选择方法，设第i个个体的适应度为F(i)，则遗传概率psi为：
(5)交叉操作：采用单点交叉，依照交叉概率pc＝0.6对每一对相互配对的个体，选定单个交叉点，交换其部分染色体；
(6)变异操作：采用基本位变异，以变异概率pm＝0.001对个体染色体的编码串中某一部分的基因值进行取反运算；
(7)终止条件判断：设定终止代数T＝300，将第300进化代数群体中的最优个体作为最优解输出，若连续15代个体的平均适应度差异小于0.5％则提前结束运算，当前群体中的最优个体作为最优解输出。
7.根据权利要求1所述的一种电磁混合悬架的控制方法，其特征在于，步骤C中所述主动控制模式、半主动控制模式之间的电路切换通过开关器件IGBT实现：车载电源与逆变器相连接，逆变器另一端与开关器件VT1相连接，开关器件VT1另一端连接开关器件VT2和直线电机M，开关器件VT2与整流桥、DC/DC变换器相连接，DC/DC变换器输出与二极管正端连接，二极管负端连接超级电容SC。
8.根据权利要求1所述的一种电磁混合悬架的控制方法，其特征在于，步骤D中所述馈能模式中采用的控制电路为：直线电机M与整流桥相连，整流桥一端与电感L一端相连，电感L另一端分别与开关S1一端和开关S2一端相连，开关S2另一端与二极管正端相连，二极管负端与超级电容SC一端相连，超级电容SC另一端分别和开关S1的另一端、整流桥的另一端相连。
9.根据权利要求1所述的一种电磁混合悬架的控制方法，其特征在于，所述的舒适、运动、综合、馈能四种模式下其可调减振阻尼器分别对应不同的四个档位。
CN201710637762.6A 2017-07-31 2017-07-31 一种电磁混合悬架的控制方法 Active CN107599778B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201710637762.6A CN107599778B (zh) 2017-07-31 2017-07-31 一种电磁混合悬架的控制方法
PCT/CN2017/100109 WO2019024164A1 (zh) 2017-07-31 2017-09-01 一种电磁混合悬架的控制方法
US16/491,879 US11584184B2 (en) 2017-07-31 2017-09-01 Control method for hybrid electromagnetic suspension

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710637762.6A CN107599778B (zh) 2017-07-31 2017-07-31 一种电磁混合悬架的控制方法

Publications (2)

Publication Number Publication Date
CN107599778A true CN107599778A (zh) 2018-01-19
CN107599778B CN107599778B (zh) 2020-09-25

Family

ID=61063962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710637762.6A Active CN107599778B (zh) 2017-07-31 2017-07-31 一种电磁混合悬架的控制方法

Country Status (3)

Country Link
US (1) US11584184B2 (zh)
CN (1) CN107599778B (zh)
WO (1) WO2019024164A1 (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108399298A (zh) * 2018-03-02 2018-08-14 中船第九设计研究院工程有限公司 一种管材切割的套料算法
CN108482199A (zh) * 2018-04-28 2018-09-04 北京理工大学 一种主动与半主动混合控制的座椅悬架及车辆
CN108674118A (zh) * 2018-05-07 2018-10-19 青岛理工大学 一种满足卡车用户悬架性能需求变更的智能悬架系统
CN109080399A (zh) * 2018-07-30 2018-12-25 江苏大学 一种可实现自供能的混合电磁悬架及其控制方法
CN109334377A (zh) * 2018-09-11 2019-02-15 江苏大学 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制系统及其构造方法
CN110001336A (zh) * 2019-03-12 2019-07-12 江苏大学 一种基于add正实网络优化的车辆isd悬架一阶理想模型
CN110155103A (zh) * 2019-05-27 2019-08-23 西南交通大学 一种基于lqg控制的高速列车半主动悬挂控制系统及方法
CN110281725A (zh) * 2019-05-24 2019-09-27 江西科技学院 存在作动器死区的主动悬架的优化控制方法及系统
CN111125837A (zh) * 2019-12-31 2020-05-08 北京理工大学 用于主动悬架动力学性能和能耗优化的控制方法
CN111546850A (zh) * 2020-03-31 2020-08-18 重庆交通大学 基于混合逻辑动态模型的车身高度与车辆姿态协调控制方法
CN112776688A (zh) * 2019-11-06 2021-05-11 操纵技术Ip控股公司 用于座椅振动消除的系统和方法
CN113352832A (zh) * 2021-07-06 2021-09-07 南昌大学 一种基于路面等级识别的多目标动态最优主动悬架控制方法
CN113635725A (zh) * 2021-07-23 2021-11-12 中国北方车辆研究所 一种机电悬架多模式控制装置
CN115360912A (zh) * 2022-08-17 2022-11-18 燕山大学 基于pi前馈模型预测控制算法的新型馈能型悬架控制系统

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112364431B (zh) * 2020-10-12 2024-04-05 宁波思明汽车科技股份有限公司 一种馈能悬架馈能装置的优化方法
CN112895833B (zh) * 2021-02-05 2022-10-25 奇瑞新能源汽车股份有限公司 悬架作动器及悬架作动器的控制方法
CN112948971B (zh) * 2021-03-04 2023-09-29 北京交通大学 一种磁悬浮列车速度曲线的节能优化方法
CN115214812B (zh) * 2022-07-25 2023-08-22 北京理工大学 一种高速履带车辆悬挂动行程跟踪主动控制方法
CN117708997B (zh) * 2024-02-05 2024-05-03 中汽研汽车检验中心（天津）有限公司 基于遗传算法的车辆悬架设计参数确定方法、设备和介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101417596A (zh) * 2008-11-18 2009-04-29 江苏大学 一种汽车空气悬架的控制方法
WO2012139406A1 (zh) * 2011-04-12 2012-10-18 Zhao Wenzhi 一种汽车系统综合改进装置
CN103419588A (zh) * 2013-07-30 2013-12-04 江苏大学 三级阻尼可调的主动馈能悬架分层控制器及其构造方法
CN105172510A (zh) * 2015-09-30 2015-12-23 江苏大学 一种基于虚拟惯容的车辆悬架控制系统及控制方法
RU2603703C1 (ru) * 2015-08-28 2016-11-27 Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") Способ снижения динамической нагруженности транспортного средства при движении по поверхности и транспортное средство

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6598885B2 (en) 2001-10-23 2003-07-29 Liquidspring Technologies, Inc. Single valve control of damping and stiffness in a liquid spring system
JP5212015B2 (ja) * 2008-10-27 2013-06-19 アイシン精機株式会社 減衰力制御装置
JP5668872B2 (ja) 2011-12-28 2015-02-12 日産自動車株式会社 車両の制御装置
JP6396414B2 (ja) * 2013-03-15 2018-09-26 クリアモーション，インコーポレイテッド 多経路流体ダイバータバルブ
CN104309438B (zh) 2014-10-20 2016-07-20 江苏大学 一种多工况车辆悬架
JP6347248B2 (ja) * 2015-10-07 2018-06-27 トヨタ自動車株式会社 車両用サスペンション制御装置
CN105480043B (zh) * 2015-12-17 2016-08-24 西安科技大学 一种多模式电磁馈能式车辆主动悬架作动器及其控制方法
CN105751845B (zh) 2016-04-08 2018-11-09 江苏大学 一种馈能型半主动悬架系统的半主动控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101417596A (zh) * 2008-11-18 2009-04-29 江苏大学 一种汽车空气悬架的控制方法
WO2012139406A1 (zh) * 2011-04-12 2012-10-18 Zhao Wenzhi 一种汽车系统综合改进装置
CN103419588A (zh) * 2013-07-30 2013-12-04 江苏大学 三级阻尼可调的主动馈能悬架分层控制器及其构造方法
RU2603703C1 (ru) * 2015-08-28 2016-11-27 Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") Способ снижения динамической нагруженности транспортного средства при движении по поверхности и транспортное средство
CN105172510A (zh) * 2015-09-30 2015-12-23 江苏大学 一种基于虚拟惯容的车辆悬架控制系统及控制方法

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108399298B (zh) * 2018-03-02 2022-02-15 中船第九设计研究院工程有限公司 一种管材切割的套料算法
CN108399298A (zh) * 2018-03-02 2018-08-14 中船第九设计研究院工程有限公司 一种管材切割的套料算法
CN108482199A (zh) * 2018-04-28 2018-09-04 北京理工大学 一种主动与半主动混合控制的座椅悬架及车辆
CN108674118A (zh) * 2018-05-07 2018-10-19 青岛理工大学 一种满足卡车用户悬架性能需求变更的智能悬架系统
CN108674118B (zh) * 2018-05-07 2021-10-01 青岛理工大学 一种满足卡车用户悬架性能需求变更的智能悬架系统
CN109080399A (zh) * 2018-07-30 2018-12-25 江苏大学 一种可实现自供能的混合电磁悬架及其控制方法
CN109334377B (zh) * 2018-09-11 2021-07-20 江苏大学 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制系统及其构造方法
CN109334377A (zh) * 2018-09-11 2019-02-15 江苏大学 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制系统及其构造方法
CN110001336A (zh) * 2019-03-12 2019-07-12 江苏大学 一种基于add正实网络优化的车辆isd悬架一阶理想模型
CN110001336B (zh) * 2019-03-12 2021-11-23 江苏大学 一种基于add正实网络优化的车辆isd悬架一阶理想模型
CN110281725A (zh) * 2019-05-24 2019-09-27 江西科技学院 存在作动器死区的主动悬架的优化控制方法及系统
CN110155103B (zh) * 2019-05-27 2020-10-23 西南交通大学 一种基于lqg控制的高速列车半主动悬挂控制系统及方法
CN110155103A (zh) * 2019-05-27 2019-08-23 西南交通大学 一种基于lqg控制的高速列车半主动悬挂控制系统及方法
CN112776688A (zh) * 2019-11-06 2021-05-11 操纵技术Ip控股公司 用于座椅振动消除的系统和方法
CN111125837A (zh) * 2019-12-31 2020-05-08 北京理工大学 用于主动悬架动力学性能和能耗优化的控制方法
CN111125837B (zh) * 2019-12-31 2022-12-09 北京理工大学 用于主动悬架动力学性能和能耗优化的控制方法
CN111546850A (zh) * 2020-03-31 2020-08-18 重庆交通大学 基于混合逻辑动态模型的车身高度与车辆姿态协调控制方法
CN113352832A (zh) * 2021-07-06 2021-09-07 南昌大学 一种基于路面等级识别的多目标动态最优主动悬架控制方法
CN113635725A (zh) * 2021-07-23 2021-11-12 中国北方车辆研究所 一种机电悬架多模式控制装置
CN113635725B (zh) * 2021-07-23 2024-03-08 中国北方车辆研究所 一种机电悬架多模式控制装置
CN115360912A (zh) * 2022-08-17 2022-11-18 燕山大学 基于pi前馈模型预测控制算法的新型馈能型悬架控制系统

Also Published As

Publication number Publication date
US20200031188A1 (en) 2020-01-30
WO2019024164A1 (zh) 2019-02-07
CN107599778B (zh) 2020-09-25
US11584184B2 (en) 2023-02-21

Similar Documents

Publication Publication Date Title
CN107599778A (zh) 一种电磁混合悬架的控制方法
CN103419588B (zh) 三级阻尼可调的主动馈能悬架分层控制器及其构造方法
CN110171261B (zh) 一种磁流变馈能悬架减振与发电半主动协调控制方法
Long et al. Regenerative active suspension system with residual energy for in-wheel motor driven electric vehicle
Zand et al. Using adaptive fuzzy logic for intelligent energy management in hybrid vehicles
CN109080399A (zh) 一种可实现自供能的混合电磁悬架及其控制方法
CN101844498B (zh) 无外界动力源的半主动/主动复合控制悬架及控制方法
Liu et al. Modeling and simulation of energy-regenerative active suspension based on BP neural network PID control
Yu et al. A loss-minimization port-controlled Hamilton scheme of induction motor for electric vehicles
CN108454343A (zh) 基于姿态补偿的混合电磁悬架节能型半主动控制方法
CN111709163A (zh) 一种基于电磁主动馈能悬架的提升轮毂驱动汽车性能的方法
CN115793445A (zh) 一种基于多智能体深度强化学习的混合动力汽车控制方法
CN103625236B (zh) 确定基于分级变压充电的esasre悬架充电电压方法
CN113177267B (zh) 一种基于改进模糊pid的全流程多学科建模方法
CN114619824A (zh) 车身控制方法、控制装置、电子设备及存储介质
CN116201890B (zh) 一种多挡amt纯电动城市客车自适应换挡规律设计方法
yin Tang et al. Research on suspension system based on genetic algorithm and neural network control
Mohamed et al. LQR and PID controllers performance on a half car active suspension system
CN208914919U (zh) 一种悬架系统双向馈能控制系统
CN112947087B (zh) 半主动悬架增强型多重模糊pid控制系统及方法
CN109130758B (zh) 一种悬架系统双向馈能控制系统及方法
CN112440974B (zh) 基于分布式深度确定性策略梯度的hev能量管理方法
Zhang et al. An online energy management strategy for fuel cell hybrid vehicles
CN106168800A (zh) 电磁悬架调节控制器
Liu et al. Driver-in-the-loop handling stability control of 4WID-EV

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Sun Zeyu

Inventor after: Wang Ruochen

Inventor after: Qian Yuchen

Inventor after: Meng Xiangpeng

Inventor after: Xie Jian

Inventor after: Chen Long

Inventor before: Wang Ruochen

Inventor before: Qian Yuchen

Inventor before: Meng Xiangpeng

Inventor before: Sun Zeyu

Inventor before: Xie Jian

Inventor before: Chen Long

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant