[0023] 下面结合附图和实施实例对本发明作进一步描述。
[0024] 本发明提出的基于多智能体的分布式驱动汽车集成控制系统,由信号采集层agent、决策层agent、协调层agent、区域层agent、执行层agent和应用知识库组成,对应的硬件装置如图1所示,图中,决策控制器为决策层agent的硬件装置,协调控制器为协调层agent的硬件装置,EPS控制器和转矩分配控制器为区域层agent的硬件装置,助力电机和驱动电机为执行层agent的硬件装置,信号采集装置为信号采集层agent的硬件装置,信号采集装置包括车速传感器、转向盘转角传感器、转向盘转矩传感器、助力电机电枢电流传感器、油门踏板信号处理模块、制动踏板信号处理模块、驱动电机状态监测模块、陀螺仪;其中,决策控制器的输出端与协调控制器相连,应用知识库分别与EPS控制器和转矩分配控制器双向连接,协调控制器、EPS控制器、转矩分配控制器、助力电机、驱动电机和信号采集装置通过CAN总线相连。
[0025] 应用知识库用于存储初始化数据,包括:不同车辆行驶工况所对应的EPS系统和转矩分配系统所需要实际补偿的力矩大小关系表,以及用于实时计算无法从表中直接获取数据的近似拟合函数。
[0026] 信号采集层agent用于采集车辆当前动力学参数信号,并发送给决策层agent和区域层agent;决策层agent实时分析计算来自于信号采集agent的车辆行驶参数数据,得出分布式驱动汽车的当前所处的行驶工况,并向各协调层agent发送指令;协调层agent根据当前的行驶工况和应用知识库存储的预先设计好的经验数据,协调EPS系统与转矩分配系统,计算得出两个系统所占权重,并向区域层agent发送指令;区域层agent根据EPS系统与转矩分配系统的权重,结合信号采集层agent传送过来的数据信息,与应用知识库存储的数据进行分析和对比,得出当前EPS系统与转矩分配系统各自所需要得输出量,并向各执行层agent发送指令进行控制;执行层agent用于将区域层agent的输出命令转化成物理量实现输出,同时通过助力电机电枢电流传感器和驱动电机状态监测模块检测输出的执行情况,并通过CAN通信模块发送给区域层agent。
[0027] 本发明还提供基于多智能体的分布式驱动汽车集成控制方法,包括以下步骤:
[0028] 1)系统初始化,并将初始化信息存储在应用知识库中;其中,初始化过程为:首先,针对不同车辆行驶工况建立所对应的EPS系统和转矩分配系统所需要实际补偿的力矩大小关系表;其次,对EPS系统和转矩分配系统所需要实际补偿的力矩大小关系表进行拟合,获取近似拟合函数,用于实时计算出无法从表中直接获取的数据。
[0029] 2)信号采集层agent采集EPS系统信号、转矩分配系统信号和车况信号,通过CAN总线发送给决策层agent和区域层agent,其中,EPS系统信号包括转向盘转角传感器信号、转向盘转矩传感器信号和助力电机电枢电流传感器信号;转矩分配系统信号包括驱动电机状态监测模块提供的驱动电机转速信号、转矩信号和陀螺仪信号;车况信号包括车速传感器信号、油门踏板信号和制动踏板信号。
[0030] 3)决策层agent根据信号采集层agent传递过来的信息,判断车辆当前状态,并向协调层agent发送指令;其中,车辆当前状态分为常规转向工况、转向制动/驱动工况、常规行驶工况与紧急避让工况四种,判断方法如图2所示,包括如下步骤:
[0031] 3.1)判断当前转向盘转角δsw是否等于0,即判断驾驶员是否有转向的意图,如果转向盘转角δsw等于0,转步骤3.2),否则转步骤3.3);
[0032] 3.2)判断转矩传感器是否检测到突变阶跃脉冲信号,若有突变阶跃脉冲信号则车辆当前状态为紧急避让工况,否则为常规行驶工况;
[0033] 3.3)启动计时器,对δsw不等于0的时间进行计时,转步骤3.4);
[0034] 3.4)判断计时器时间是否大于2s,如果不大于2s,转步骤3.1),否则可判定驾驶员正在转向过程中,转步骤3.5);
[0035] 3.5)判断驱动/制动踏板信号是否为0,如果为0则车辆当前状态为常规转向工况,否则车辆当前状态为转向制动/驱动工况。
[0036] 4)协调层agent计算当前工况下,EPS系统和转矩分配系统各自所占的权重系数和实际需要补偿的力矩,使EPS系统和转矩分配系统实现相对最优,并向区域agent发送指令;各工况下EPS系统和转矩分配系统各自所占的权重系数和实际需要补偿的力矩的计算方法如下:
[0037] 对于常规转向工况,转矩分配控制器采取差速控制,即保证左前轮和右前轮的转矩,以及左后轮和右后轮转矩保持相等,同时满足
[0038] ωlf+ωrf=2ωfo
[0039] ωlr+ωrr=2ωro
[0040] 其中ωlf为左前轮角速度,ωrf为右前轮角速度,ωfo为车辆前轴中点处角速度,ωlr为左后轮角速度,ωrr为右后轮角速度,ωro为车辆等效后轴中点处角速度。对于各个轮毂驱动电机而言,可以将角速度用电机转速n来表示,根据公式 计算,即
[0041] nlf+nrf=2nfo
[0042] nlr+nrr=2nro
[0043] 同时对于EPS系统,当前状态属于正常转向工况,只需要EPS控制器提供常规助力控制即可。
[0044]
[0045] 其中Tmr为助力电机助力力矩,eh为目标助力电机电枢电流和实际助力电机电枢电流的偏差,Kp为比例系数,Ki为积分系数,Kd为微分系数。通过Ziegler-Nichols整定方法确定比例系数Kp为1.695,积分系数Ki为3.673,微分系数Kd为0.26。
[0046] 对于转向驱动/制动工况,通过陀螺仪检测出横摆角加速度的变化值Δωr′,如下式所示:
[0047] Mz′=IZ·Δωr′
[0048] 其中IZ为车辆绕Z轴的转动惯量,根据车辆参数可知IZ=1157.1kg·m2,通过上式就可以求得系统需要补偿的总力矩Mz′,经过大量仿真试验可求出下表,即得出经过协调后的EPS系统和转矩分配系统实际需要补偿的力矩,而查表只是给出了某些特定的情况下的数据,如只有Δωr′=0.2或0.4等情况的数据,而Δωr′取其他数据则无法查表,而应是通过对已有数据进行拟合,获取近似拟合函数,用来实时计算出无法从表中直接获取的数据。
[0049] 表1转向驱动/制动工况(规定向左转向为正)
[0050]Δωr′ TEPS_real TTD_real
0 0 0
0.2 23.14 208.28
0.4 21.53 441.31
0.6 20.89 673.37
0.8 20.45 905.23
1.0 19.87 1137.23
1.2 17.56 1370.96
1.4 16.86 1603.08
1.6 15.87 1835.49
1.8 15.38 2067.40
2.0 14.29 2299.91
2.2 14.13 2531.49
2.4 13.79 2763.25
2.5 13.75 2879.00
[0051] 对于紧急避让工况,通过陀螺仪检测出侧向加速度的变化值Δay,如下式所示[0052] Mz′=m·Δay·lr
[0053] 其中m为车辆总质量,根据车辆参数可知m=830kg,lr为车辆质心到后轴的距离,根据车辆参数可知lr=1.244m。通过上式就可以求得系统需要补偿的总力矩Mz′,经过大量仿真试验可求出下表,即得出经过协调后的EPS系统和转矩分配系统需要提供的实际力矩。
[0054] 表2紧急避让工况(规定向左转向为正)
[0055]Δay TEPS_real TTD_real
0 0 0
0.3 15.49 294.27
0.6 15.35 604.16
0.9 15.21 914.06
1.2 14.97 1224.05
1.5 14.85 1533.93
1.8 14.78 1843.76
2.1 14.55 2153.74
2.4 14.32 2463.73
2.7 13.88 2773.92
3.0 13.54 3084.02
3.3 13.27 3394.05
3.6 13.18 3703.89
3.9 13.10 4013.73
[0056] 5)区域层agent接收协调层agent传递过来的信息,送EPS控制器和转矩分配控制器,分别计算得到EPS控制器的控制策略和转矩分配控制器的转矩分配算法,并向各执行层agent发送指令进行控制;
[0057] EPS控制器根据当前车速、转向盘转矩以及可变助力特性曲线计算得出当前系统应该提供的助力转矩,同时根据车辆不同的运动状态给出相应的控制策略,其中,控制策略包括常规助力控制和补偿控制;
[0058] 转矩分配控制器根据车辆的当前运动状态,实时给出相应的转矩分配算法,其中,转矩分配算法包括常规差速控制和左、右轮驱动力矩分配控制;不同工况下的计算方法如下:
[0059] 对于常规转向工况,转矩分配控制器采取差速控制,即保证左前轮和右前轮的转矩,以及左后轮和右后轮转矩保持相等,同时满足
[0060] ωlf+ωrf=2ωfo
[0061] ωlr+ωrr=2ωro
[0062] 其中ωlf为左前轮角速度,ωrf为右前轮角速度,ωfo为车辆前轴中点处角速度,ωlr为左后轮角速度,ωrr为右后轮角速度,ωro为车辆等效后轴中点处角速度。对于各个轮毂驱动电机而言,可以将角速度用电机转速n来表示,根据公式 计算,即
[0063] nlf+nrf=2nfo
[0064] nlr+nrr=2nro
[0065] 同时对于EPS系统,当前状态属于正常转向工况,只需要EPS控制器提供常规助力控制即可。
[0066]
[0067] 其中Tmr为助力电机助力力矩,eh为目标助力电机电枢电流和实际助力电机电枢电流的偏差,Kp为比例系数,Ki为积分系数,Kd为微分系数。通过Ziegler-Nichols整定方法确定比例系数Kp为1.695,积分系数Ki为3.673,微分系数Kd为0.26。
[0068] 对于转向制动/驱动工况
[0069] EPS系统:
[0070] TEPS_real=ΔTmr1=-(Kpeh+Ki∫ehdt),
[0071] TEPS_real通过查表获得,其中ΔTmr1为助力电机附加的正向补偿力矩,eh为目标转向盘转角和实际转向盘转角的偏差,Kp为比例系数,Ki为积分系数。通过Ziegler-Nichols整定方法确定比例系数Kp为0.5,积分系数Ki为0.083。
[0072] 转矩分配系统:
[0073]
[0074] 其中ΔFxrf,ΔFxlf,ΔFxrr,ΔFxlr分别为右前轮、左前轮、右后轮和左后轮纵向附加驱动力,B为车轮轮距。
[0075] 轮毂驱动电机转矩与车轮纵向力的关系为:
[0076]
[0077] 其中ΔTmij为各个轮毂驱动电机附加转矩,R为车轮半径。根据横摆力矩与纵向力的关系,并带入上式可求得各个轮毂驱动电机附加转矩与转矩分配系统提供的附加横摆力矩分量的关系:
[0078]
[0079] 上式就是右前轮驱动电机控制器、左前轮驱动电机控制器、右后轮驱动电机控制器和左后轮驱动电机控制器控制各自驱动电机输出相应附加转矩的依据。具体而言,此时车辆以向左转向为例,此时左、右轮驱动力矩分配即增加左侧两个车轮的纵向驱动力矩,同时降低右侧两个车轮的纵向驱动力矩。
[0080] 将各轮转矩用TTD_real表示,如下式:
[0081]
[0082] TTD_real通过查表获得。
[0083] 对于紧急避让工况,与转向制动/驱动工况同理,不同之处只是查询表格的不同。
[0084] 6)执行层agent根据区域agent层传递过来的信息,负责将区域层agent的输出命令转化成物理量实现输出,同时通过助力电机电枢电流传感器和驱动电机状态监测模块分别对EPS系统和转矩分配系统进行检测,对各自输出的执行情况进行实时反馈,并通过CAN总线发送给区域层agent。