CN113547928B - 一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 - Google Patents
一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 Download PDFInfo
- Publication number
- CN113547928B CN113547928B CN202110795840.1A CN202110795840A CN113547928B CN 113547928 B CN113547928 B CN 113547928B CN 202110795840 A CN202110795840 A CN 202110795840A CN 113547928 B CN113547928 B CN 113547928B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- motor
- torque
- efficiency
- drive system
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012546 transfer Methods 0.000 claims abstract description 28
- 238000011217 control strategy Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 17
- 238000005265 energy consumption Methods 0.000 claims abstract description 16
- 230000004044 response Effects 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000000875 corresponding effect Effects 0.000 description 13
- 238000012360 testing method Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000007599 discharging Methods 0.000 description 3
- 238000009472 formulation Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000003638 chemical reducing agent Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- QYAPHLRPFNSDNH-MRFRVZCGSA-N (4s,4as,5as,6s,12ar)-7-chloro-4-(dimethylamino)-1,6,10,11,12a-pentahydroxy-6-methyl-3,12-dioxo-4,4a,5,5a-tetrahydrotetracene-2-carboxamide;hydrochloride Chemical compound Cl.C1=CC(Cl)=C2[C@](O)(C)[C@H]3C[C@H]4[C@H](N(C)C)C(=O)C(C(N)=O)=C(O)[C@@]4(O)C(=O)C3=C(O)C2=C1O QYAPHLRPFNSDNH-MRFRVZCGSA-N 0.000 description 1
- 102100026127 Clathrin heavy chain 1 Human genes 0.000 description 1
- 101000912851 Homo sapiens Clathrin heavy chain 1 Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L15/00—Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles
- B60L15/20—Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles for control of the vehicle or its driving motor to achieve a desired performance, e.g. speed, torque, programmed variation of speed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L2240/00—Control parameters of input or output; Target parameters
- B60L2240/40—Drive Train control parameters
- B60L2240/42—Drive Train control parameters related to electric machines
- B60L2240/423—Torque
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Arrangement And Driving Of Transmission Devices (AREA)
Abstract
本发明涉及一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法,属于新能源汽车底盘智能控制领域。该方法包括:S1:搭建整车动力学及各部件模型,包括车辆动力学模型、电池模型、电机模型、轮胎模型和电驱动系统损耗模型;S2:分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,得出主动防滑控制策略;S3:综合考虑由于工况的变化导致的工作模式频繁切换而带来的电驱动系统损耗,得到整车驱动系统最优效率曲面;S4:构建考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略。本发明在保障车辆安全稳定性的前提下,显著降低能耗和模式切换频率,改善了车辆的经济性,提高了整车的舒适性。
Description
技术领域
本发明属于新能源汽车底盘智能控制领域,涉及一种考虑轮胎滑移的整车效率最优的双电机四驱电动汽车转矩分配方法。
背景技术
随着社会环保需求的日益提高,从节能的角度对电动汽车提出了更高的要求。在现有电池技术未取得重大突破的前提下,为进一步改善纯电动汽车的经济性,通过轻量化设计或高效的电驱动系统构型等方法使整个系统效率达到最优,从而实现最长的续驶里程。在驱动系统构型的选择上,前后轴双电机驱动方式相比传统的单电机四驱驱动方式,省略前后轴传动系统,降低了机械损耗,减少整备质量,同时通过直接调节各驱动电机驱动力可控制车轮附着特性,使整车获得最大驱动力;而相比分布式轮毂电机驱动控制复杂,空间布置难度大,成本高等特点,前后轴双电机驱动电动汽车逐渐成为汽车行业关注的焦点之一。前后轴双电机驱动电动汽车的转矩分配控制是新能源汽车底盘智能控制的关键技术,对整车性能起到决定性作用。因此开展新能源汽车底盘智能控制的关键技术研究,对汽车产业的转型和提升新能源汽车的研发水平有着十分重要的意义。
车辆安全稳定性和经济性是整车性能的重要指标,车辆驱动转矩分配控制直接决定了整车的性能表现。现阶段以经济性为研究目标,提出电机或系统效率优化的转矩分配策略,并未提及该策略下转矩分配是否满足当前驱动车轮附着条件。如果转矩分配不合理,过大的转矩输出使得部分车轮的驱动力大于地面提供的最大附着力,会使车轮发生滑转,从而使车辆触发ABS或ASR防滑系统来保证车辆的行驶稳定性。没有考虑动态工况下前后轴荷变化引起的整车附着条件变化而带来的整车能耗变化。在良好路面下,以经济性为优化目标所提出的转矩分配策略如果频繁触发ABS或ASR防滑系统,在这些被动防滑策略的控制下,虽然抑制了车轮打滑,但过大的转矩输出会增加整车能耗,说明该优化策略是不合理的。同时,现有技术以经济性为目标所提出的驱动系统最优效率曲面,划分出在不同工况点下对应的不同工作模式,通过不同工作模式的选择来改善车辆经济性;但忽略了前后电机的拖拽损耗、励磁损耗、开关损耗等由于工况的变化导致的工作模式频繁切换而带来的损耗。因此,现有技术所使用的驱动系统效率曲面不一定是最优效率,这也进一步影响了转矩分配策略的制定与控制,从而限制了车辆经济性能的改善空间。
发明内容
有鉴于此,本发明的目的在于提供一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法,结合整车行驶的附着条件,考虑前后轴荷转移的动态约束变化,综合考虑由于工况的变化导致的工作模式频繁切换而带来的前后电机拖拽损耗、励磁损耗、开关损耗等损耗问题,所得到的整车驱动系统最优效率曲面更能反映出实车的真实情况,可为转矩分配策略的制定与控制提供更好的决策依据;同时,在不触发ABS或ASR防滑系统控制的前提下,尽可能的减小车轮过度滑转现象,起到主动防滑作用,保障车辆的安全稳定性,从而避免过大的转矩输出带来的整车能耗增加,进一步改善车辆经济性。
为达到上述目的,本发明提供如下技术方案:
一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法,具体包括以下步骤:
S1:搭建整车动力学及各部件模型,包括车辆动力学模型、电池模型、电机模型、轮胎模型和电驱动系统损耗模型;
S2:分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,得出主动防滑控制策略;
S3:综合考虑由于工况的变化导致的工作模式频繁切换而带来的电驱动系统损耗(前后电机拖拽损耗、励磁损耗、阶跃转矩瞬态响应实际功率损耗等损耗问题),得到整车驱动系统最优效率曲面;
S4:构建考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略。
进一步,步骤S2中,得出主动防滑控制策略,具体包括:在可得到实时路面附着系数的情况下,通过考虑前后轴荷转移的车辆动力性模型,得到前/后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe;再与前/后电机最大峰值转矩Tmf_max、Tmr_max做对比,较小的一方作为外特性约束边界,从而对车辆动力性进行动态约束,以起到主动防滑的作用。
进一步,步骤S2中,前后轴荷转移的约束条件和整车行驶的附着条件,具体包括:
其中,Fzf为前轴垂直载荷,Fzr为后轴垂直载荷,Tmf_safe、Tmr_safe分别为前/后轴所能提供的最大附着力对应的电机轴端输出转矩,ηtran为传动系统效率,为路面附着系数,if为前轴转动比,ir为后轴转动比;Tmf、Tmr为前/后电机轴端实际输出转矩,Tmf_max、Tmr_max为前/后电机最大峰值转矩。
进一步,步骤S3中,得到整车驱动系统最优效率曲面,具体包括:根据工况,综合考虑电驱动系统损耗模型与边界区域效率差值模型,以经济性和切换频率为目标,分析模式边界带宽的影响因素,进行边界寻优;然后分析不同边界带宽下,每一段车速内,所有工况点的能耗及切换频率进行多目标边界寻优,得到最优的整车驱动系统综合效率曲面。
进一步,所述的边界区域效率差值模型的表达式为:Δη=f(T,n,η1,η2),其中,T为当前总需求转矩,n为转速,η1为前电机效率,η2为后电机效率。
进一步,综合考虑整车动力性约束的动态变化,结合电池、电机和传动系统的效率,得到不同(前驱、后驱、四驱)驱动模式下整车驱动系统综合效率;
整车效率取决于车速、加速度、电池效率、电机效率和传动系统效率;
以整车效率为目标函数maxηsys,得到各驱动模式下整车行驶效率最优的前后电机的最佳转矩分配,进而得到整车的模式划分规则;其中,以整车效率为目标函数的约束条件为:
其中,SOC为电荷状态,SOCmin为最小截止电荷状态,Pbat为电池功率,Pbat_max为最大电池放电功率,Tmf、Tmr为前/后电机轴端实际输出转矩,Tmf_max、Tmr_max为前/后电机最大峰值转矩,Tmf_safe、Tmr_safe分别为前/后轴所能提供的最大附着力对应的电机轴端输出转矩,nf、nr分别为前/后电机转速,nf_max,nr_max分别为前/后电机最大转速。
进一步,模式切换的控制策略为:根据不同模式下的整车驱动系统综合效率,得到在不同车速和不同需求转矩下的不同模式划分规则,即各模式的运行区间;然后根据速度和总需求转矩的变化自动选择工作模式。
进一步,步骤S4中,构建的多目标双电机四驱转矩分配控制策略,具体包括:根据步骤1搭建的整车动力学及各部件模型,以及步骤2得到的主动防滑控制策略,驾驶员模型通过驱动/制动踏板位置对总需求转矩进行解析,根据当前的车速、需求转矩,通过步骤3得到的整车驱动系统最优效率曲面进行模式划分;考虑前后轴荷转移的动态变化,对当前模式下的总需求转矩进行前/后转矩分配Tf_req、Tr_req,再与前/后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe以及前/后电机外特性所能提供的最大转矩Tmf_max、Tmr_max做对比,得到前后轴实际输出转矩Tmf、Tmr;通过车辆动力学模型输出当前状态下的车速、前后车轮角速度,然后进行滑移率判断;若当前滑移率小于路面附着系数对应的最优滑移率,则进入下一个时刻,从而起到主动防滑的作用。
本发明的有益效果在于:
1)本发明根据所搭建的整车及前后轴系统动力学模型,分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,提出主动防滑控制策略,在不触发ABS(防抱死刹车系统,控制所有车轮)或ASR(驱动防滑系统,只控制驱动轮)防滑系统控制的前提下,尽可能的减小车轮过度滑转现象,保障车辆的安全稳定性。
2)本发明综合考虑了由于工况的变化导致的工作模式频繁切换而带来的前后电机拖拽损耗、励磁损耗、开关损耗等损耗问题,所得到的整车驱动系统最优效率曲面更能反映出实车的真实情况,可为转矩分配策略的制定与控制提供更好的决策依据。
3)本发明根据车辆车速、需求转矩、SOC等状态变量,着重分析模式切换带来的电驱动系统损耗以及模式之间的边界区域效率差值对模式边界带宽确定的影响规律;以经济性和切换频率为目标,分析不同边界带宽下,每一段车速内,所有工况点的能耗及切换频率进行多目标边界寻优。
4)本发明所提出的考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略方法,驾驶员模型通过驱动/制动踏板位置对总需求转矩进行解析,根据当前的车速、需求转矩通过最优整车驱动系统综合效率曲面进行模式划分;考虑前后轴荷转移的动态变化,对当前模式下的总需求转矩进行前后转矩分配,再与主动防滑控制策略下的输出转矩作对比,得到前后轴实际输出转矩。在保障车辆安全稳定性的前提下,能耗显著降低,改善了车辆的经济性,并极大的降低了模式切换频率,提高了整车的舒适性。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为双电机四驱电动汽车结构示意图;
图2为整车纵/垂向动力学模型示意图;
图3为前驱与四驱转矩分配响应过程;
图4为考虑轴荷转移及附着条件动态约束的主动防滑控制策略;
图5为前轴驱动模式电机轴端提供最大输出转矩;
图6为后轴驱动模式电机轴端提供最大输出转矩;
图7为四轮驱动模式电机轴端提供最大输出转矩;
图8为60%SOC各模式下整车最优综合效率;
图9为60%SOC模式划分;
图10为各工况下的模式切换频率;
图11为60%SOC带有模式边界的模式划分;
图12为工况点对应的模式选择;
图13为WLTC(550s-900s)无边界与有边界的模式切换对比;
图14为整车驱动系统最优效率曲面;
图15为考虑轴荷转移及附着条件动态约束的转矩分配控制策略。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
请参阅图1~图15,本实施例针对前后独立驱动四驱电动汽车的转矩分配及驱动防滑协调控制器设计,具体的转矩分配方法主要包括以下步骤:
步骤1:搭建整车动力学及各部件模型;
本实施例所研究的双电机四驱电动汽车为前后轴各一个驱动电机,电机轴端驱动转矩经减速器传到差速器、半轴至驱动前后轴左右两边车轮,动力源由电池组提供。为满足整车驱动系统的驱动转矩需求,通过MCU直接调节各驱动电机驱动转矩,实现前后电机独立工作及双电机同时工作的3种驱动模式。图1所示为双电机四驱电动汽车结构示意图。
1)车辆动力学模型
分析双电机四驱电动汽车位于地面坐标系上的平面运动,忽略侧向力和横摆力矩作用,视左右两边垂直载荷相等,如图2所示,O为整车质心位置,b、c分别为整车质心到前后轴的纵向距离。由于车辆在运动过程中存在轴荷转移现象,将对该车辆前后轴受力情况进行分别研究,并综合分析该车辆在纵向、垂向运动时产生的惯性力和力矩。
整车纵向动力学方程为:
Fx-FR=mδa (1)
其中,Fxf为前轴轮端驱动力,Fxr为后轴轮端驱动力,m为整车整备质量,f为滚动阻力系数,CD为空气阻力系数,A为迎风面积,m为整车整备质量,a为整车加速度,δ为整车旋转质量换算系数,θ为坡度,v整车车速。
车辆在行驶过程中,前后轴垂直载荷和行驶阻力会受到升力变化的影响,同时纵向加速度也会造成载荷纵向转移,而前后轴荷的变化直接影响前后电机如何进行最优转矩分配,这对进一步研究整车能耗和车辆稳定性有着重要意义。下面分别对车辆前后轴车轮及整车垂直载荷进行分析。
前轴垂直载荷及动力学方程表达式为:
其中,Fzf为前轴垂直载荷,FRf为前轴行驶阻力,c为质心到后轴的距离,L为轴距,h为质心高度,CLf为前轴升力系数;
后轴垂直载荷及动力学方程表达式为:
其中,Fzr为后轴垂直载荷,FRr为后轴行驶阻力,Jer后轴等效转动惯量,b为质心到前轴的距离,CLr为后轴升力系数;
2)搭建驾驶员模型
本实施例采用PID控制来模拟驾驶员操作,以目标车速和实际车速作为控制输入,踏板开度作为输出信号。PID控制算式为:
其中,Kp为比例系数,Ki为积分系数,Kd为微分系数,Kaw为抗饱和系数,ev为目标车速误差,eout为踏板开度误差,y表示踏板开度,踏板开度为正时表示加速,踏板开度为负时表示减速,由于踏板开度在[-1,1]范围内变化,所以引入ysat函数作为输出少的限制边界。
由于车辆总需求转矩与踏板开度成正相关,则定义踏板开度与总需求转矩的变化关系如下:
Treq=y·Tmax(nm)=y·[Tmax_mf(nm)+Tmax_mr(nm)] (6)
其中,Tmax(nm)为转速nm时前后电机所能提供的最大转矩之和;Tmax mf(nm)为转速nm时前电机所能发出的最大转矩,Tmax_mr(nm)为转速nm时后电机能发出的最大转矩。
3)搭建电池模型
电池组的充放电特性受到温度、充放电电流以及温度的影响,为简化电池组模型,假设电池组环境温度为25℃;
电池组功率为:
Pb=Pm_in=EI-I2R0 (7)
其中,Pb为电池组的充放电功率,E为电池组的电动势,I为充放电电流,R0为等效内阻;
对式(8)求解,可得充放电电流为:
电池组荷电状态为:
其中,SOCinitial为电池组初始电荷状态,SOC(t)为时间t时的电荷状态,Cb为电池组额定容量。
4)搭建电机模型
考虑到电机响应具有一定延迟,采用惯性环节表示电机响应特性:
其中,Tm表示电机输出转矩,单位为Nm;Tcmd表示电机的目标转矩,单位为Nm;Tc为时间常数;s复变量。
电机、减速器、车轮之间通过齿轮啮合传递动力,则对于单个电机而言,其转速、转矩和功率的转换关系为:
其中,nm为电机转速,Pm_out为电机输出功率,Pm_in为电机输入功率,ηm(Tm,nm)为电机在相应转矩和转速下的工作效率。
5)车辆轮胎模型
假设4个车轮规格相同,前后驱动轮受力平衡图如图2所示。在车轮坐标系下驱动轮转动力矩平衡方程为
Jwaij=Tdij-Fxijr-Fzijfr (14)
其中,Jw为车轮转动惯量(kg·m2),Td为车轮驱动力矩(N·m),r为车轮滚动半径(m),ij代表fl,fr,rl,rr;
轮胎纵向滑移率的计算公式为
其中,ω为轮胎旋转的角速度(rad/s),rω为轮胎半径,Sd为轮胎驱动滑转率,Sb为轮胎制动滑移率。
国内外对轮胎的特性从各个角度进行了详细的研究。其中,魔术公式因为模拟准确性较高而影响力较大。因此本发明采用非线性轮胎模型来较为精确地描述轮胎力学特性。
其中,s为纵向滑移率,B为刚度因子,C为曲线形状因子,D为曲线峰值因子,E为曲线曲率因子,b0-b8为魔术公式主要参数系数。
6)电驱动系统损耗模型
①拖拽损耗
由于该整车电驱动系统具有3种驱动模式,前驱、后驱、四驱,当VCU根据当前工况决策进行后轴驱动时,后电机需求转矩为总需求转矩,前电机需求转矩为0;当前电机进入0转矩模式时,由于整车电驱动系统处于上电运行状态,故前轴驱动系统存在一定的功率损失,且随着转速的变化而变化,通过转鼓试验台架所测得的试验数据拟合得到前轴拖拽损耗:
Ploss_fkz=f(n)=a1n+b1 (17)
其中,a1为拟合系数,b1为常数。
当VCU根据当前工况决策进行前轴驱动时,前电机需求转矩为总需求转矩,后电机需求转矩为0;当后电机进入0转矩模式时,由于整车电驱动系统处于上电运行状态,故后轴驱动系统也存在一定的功率损失,且随着转速的变化而变化,通过转鼓试验台架所测得的试验数据拟合得到后轴拖拽损耗:
Ploss_rkz=f(n)=a2n+b2 (18)
其中,a2为拟合系数,b2为常数。
②励磁损耗
由于本发明所研究的双电机前后轴驱动构型中,前电机采用异步电机,后电机采用永磁同步电机。所以当前电机接到指令进入转矩模式时,异步电机需要一个励磁的过程,而这一过程中异步电机所产生的励磁损耗不可忽略不计。通过转鼓试验台架所测得的试验数据拟合得到前电机励磁损耗:
Ploss_flc=f(n)=a3n+b3 (19)
其中,a3为拟合系数,b3为常数。
③电驱动系统阶跃转矩瞬态响应功率损耗
车辆驱动模式的选择根据工况变化而发生变化。当工况的变化频繁的改变驱动模式时,即前驱、后驱、四驱3种模式会进行频繁切换。当前驱切为后驱时,如图3所示,在t0~t2时间内,由前电机提供的总需求转矩,在t1时刻转为了后电机提供;无论电机退出还是进入,电机响应是由一个变化过程的;而在这一变化过程中,前电机在t0~t2时间内的功率消耗,与前电机在t0~t1时间内的功率消耗和后电机在t1~t2时间内功率消耗之和是不相等的,并且由于模式频繁切换引起的电驱动系统阶跃转矩瞬态响应功率损耗与当前的需求转矩,转速均有关系,不可忽略不计;故通过转鼓试验台架所测得的试验数据同样可以拟合得到电驱动系统阶跃转矩瞬态响应功率损耗关系表达式,其他模式切换对应的电驱动系统阶跃转矩瞬态响应功率损耗同理可得。
Ploss_onoff=f(Tq,n) (20)
步骤2:分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,提出主动防滑控制策略;
由于车辆在行驶过程中,存在纵向轴荷转移,车辆的动力性在满足整车动力学方程的同时,也受前后车轮的路面附着条件约束和前后电机最大峰值转矩特性的约束。因此,在电池组所提供的需求功率能够满足的情况下,车辆的实际输出转矩T是否能够满足车辆需求转矩T_need,其约束条件表达式如下:
其中,Fzf为前轴垂直载荷,Fzr为后轴垂直载荷,Tmf_safe、Tmr_safe分别为前/后轴所能提供的最大附着力对应的电机轴端输出转矩,ηtran为传动系统效率,为路面附着系数,if为前轴转动比,ir为后轴转动比;Tmf、Tmr为前/后电机轴端实际输出转矩,Tmf_max、Tmr_max为前/后电机最大峰值转矩。
根据步骤1所搭建的整车及前后轴系统动力学模型,分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,提出一种主动防滑策略,如图4所示。在可得到实时路面附着系数的情况下,通过考虑前后轴荷转移的车辆动力性模型,得到前后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe;再与前后电机外特性所能提供的最大转矩Tmf_max、Tmr_max做对比,较小的一方作为外特性约束边界,从而对车辆动力性进行动态约束,以起到主动防滑的作用。
由上述提出的主动防滑控制策略,得到车辆在前驱、后驱和四驱三种驱动模式下,所能提供的最大驱动转矩边界。当车辆进入模式1前轴单独驱动时,前轴载荷向后转移,前轴载荷减小的同时又受到地面附着条件的约束,因此在保证前轴驱动不出现过度滑转的情况下,前电机所能提供的最大驱动转矩曲线Tmf_real,如图5所示。
当车辆进入模式2后轴单独制动时,后轴载荷向前转移,后轴载荷减小的同时又受到地面附着条件的约束,后电机所能提供的最大驱动转矩曲线Tmr_real,如图6所示;若转矩超过Tmr_safe后轴会出现过度滑移,造成能耗损失。
当车辆进行四轮驱动时,综合考虑各轴垂直载荷受纵向轴荷转移及路面附着条件的影响,得到某一特定路面附着系数下前后电机所能提供的有效输出转矩与加速度-速度的关系,当路面附着系数发生变化时,前后电机所能提供的有效输出转矩也会随之改变。由图7可以看出,车辆在各车速下电机所能提供的最大有效输出转矩(实线)并不等于理想状态的前后电机最大输出转矩之和(虚线)。因此,如果在理想状态下(虚线内)进行转矩分配,当转矩超过最大有效输出转矩(实线)时,车轮一定会出现过度打滑现象,造成大量的能耗损失。
步骤3:综合考虑由于工况的变化导致的工作模式频繁切换而带来的电驱动系统损耗问题,得到整车驱动系统最优效率曲面;
对于双电机四驱电动汽车的前后轴电机之间的转矩如何分配直接影响整车系统效率。整车效率取决于车速、加速度、电池效率、电机效率和传动系统效率。
当车辆驱动时,整车综合效率为:
其中,ηmf、ηmr为前后电机效率,ηm_dis为前后电机总驱动效率,ηbat_dis为电池组放电效率,ηtran为传动系统效率,ηsys_dis为整车驱动系统总效率。
当车辆制动时,整车综合效率为:
其中,ηm_chg为前后电机等效制动能量回收总效率,ηbat_chg为电池组充电效率,η_sys_chg为整车等效制动能量回收总效率。
因此,以整车效率为目标函数maxηsys,可以得到各驱动模式下整车行驶效率最优的前后电机的最佳转矩分配,进而得到整车的模式划分规则。公式(25)为以整车效率为目标函数的约束条件。
约束条件为:
其中,SOC为电荷状态,SOCmin为最小截止电荷状态,Pbat为电池功率,Pbat_max为最大电池放电功率,nf、nr分别为前/后电机转速,nf_max,nr_max分别为前/后电机最大转速。
综合考虑整车动力性约束的动态变化,结合电池、电机和传动系统的效率,可得到前驱、后驱、四驱三种驱动模式下整车驱动系统综合效率。如图8所示,当SOC为60%时,图中三个曲面分别是前驱(FWD)、后驱(RWD)、四驱(AWD)模式的驱动系统效率,用颜色和箭头表示。通过比较各运行点不同模式下的驱动系统效率,选择最佳运行点,可以得到速度和总需求转矩全范围内的整车驱动系统综合效率。
为了实现更好的整车系统效率,结合车辆动力性约束的动态变化,满足考虑前后轴荷变化的驱动转矩需求,有必要制定可行的模式切换控制策略。根据不同模式下的整车驱动系统综合效率,可以得到在不同车速和不同需求转矩下的三种模式划分规则,即各模式的运行区间。根据速度和总需求转矩的变化自动选择工作模式。
图9显示了当SOC为60%时3种模式中每一种模式在速度和总需求转矩状态下的运行区间。该车辆在低速、低加速区域选择mode2后驱模式,在中速/高速、低加速区域选择mode1前驱模式,可以看出在低加速区域单电机运行模式以最佳的整车效率满足车辆驱动需求;为了获得更高的动力性,前后电机同时工作进入mode3四轮驱动模式,以最佳的整车效率满足车辆动力性需求。
然而根据上述已有的模式划分进行NEDC、UDDS、CLTC、WLTC等经济性循环工况试验分析,发现在4种不同的循环工况下存在非常高频的模式切换,如图10所示;这会引起高频的冲击,极大的降低车辆乘坐的舒适性,且频繁的模式切换不仅增加能耗,还对电驱动控制技术要求极高。因此,本实施例提出模式边界作为模式切换的过度区,如图11所示,从而降低切换频率,提高车辆乘坐的舒适性,同时避免能耗的增加。
图12为图11的局部放大,从图12中可以看出,当工作点模式2走到边界0模式时,维持上一个模式;当工况点继续走到黄色区域模式3时,才真正的进入模式3状态。这样可以极大的避免工作点仅在边界处徘徊导致的模式2与模式3频繁切换的问题。同时,电机响应需要一个时间,如果模式维持的时间小于电机响应时间以及can总线信号传输时间,那么当前模式的选择与执行是不必要的,也会增加能量的消耗。如图13所示,模式3维持时间延长,降低了切换频率,也避免了极短时间的模式切换。
对于如何定义模式边界带宽大小,其影响因素与车速和转矩波动误差有关;不同的边界带宽会得到不同的整车驱动系统综合效率曲面。边界带宽区域越大,模式切换频率越低,电驱动损耗越小,能耗越低;但模式边界是维持上一个模式,区域越大,说明两模式之间的效率差值越大,经济性越差;效率差值的大小与两电机本体的效率有关;因此,模式边界带宽区域越大,模式切换频率固然降低,但对于能耗而言,模式切换引起的电驱动损耗与两模式之间的效率差值两者相互矛盾。
为解决这一问题,根据工况,综合考虑步骤1所建立电驱动系统损耗模型与边界区域效率差值模型公式(26),以经济性和切换频率为目标,分析模式边界带宽的影响因素,进行边界寻优,如表1所示。分析不同边界带宽下,每一段车速内,所有工况点的能耗及切换频率进行多目标边界寻优,如图14所示,得到最优的整车驱动系统综合效率曲面。
边界区域效率差值
Δη=f(T,n,η1,η2) (26)
其中,T为当前总需求转矩,n为转速,η1为前电机效率,η2为后电机效率。
表1多目标优化中具体的控制变量与影响因素
步骤4:构建考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略;
根据步骤1所搭建的整车及前后轴系统动力学模型,以及步骤2所提出一种主动防滑控制策略,如图15所示,驾驶员模型通过驱动/制动踏板位置对总需求转矩进行解析,根据当前的车速、需求转矩通过由步骤3所得到的最优整车驱动系统综合效率曲面进行模式划分;考虑前后轴荷转移的动态变化,对当前模式下的总需求转矩进行前后转矩分配Tf_req、Tr_req,再与前后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe以及前后电机外特性所能提供的最大转矩Tmf_max、Tmr_max做对比,得到前后轴实际输出转矩Tmf、Tmr;通过车辆动力学模型输出当前状态下的车速、前后车轮角速度,然后进行滑移率判断;若当前滑移率小于路面附着系数对应的最优滑移率,则进入下一个时刻,从而起到主动防滑的作用。在这一策略下车轮一般不会出现打滑的现象,除非在路面附着系数极低且动力需求极高的变极端工况下,为了保障车辆的安全稳定性,在极端工况若出现打滑现象,才去触发ABS或ASR防滑系统作为最后保障。
综上可以看出,本实施例所提出的考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略方法,在保障车辆安全稳定性的前提下,能耗显著降低,改善了车辆的经济性,并极大的降低了模式切换频率,提高了整车的舒适性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法,其特征在于,该方法具体包括以下步骤:
S1:搭建整车动力学及各部件模型,包括车辆动力学模型、电池模型、电机模型、轮胎模型和电驱动系统损耗模型;
S2:分析前后轴荷转移的动态约束变化,结合整车行驶的附着条件,得出主动防滑控制策略;其中,前后轴荷转移的约束条件和整车行驶的附着条件,具体包括:
其中,Fzf为前轴垂直载荷,Fzr为后轴垂直载荷,Tmf_safe、Tmr_safe分别为前/后轴所能提供的最大附着力对应的电机轴端输出转矩,ηtran为传动系统效率,为路面附着系数,if为前轴转动比,ir为后轴转动比,rω为轮胎半径;Tmf、Tmr为前/后电机轴端实际输出转矩,Tmf_max、Tmr_max为前/后电机最大峰值转矩;
S3:综合考虑电驱动系统的拖拽损耗、异步电机所产生的励磁损耗以及由于工况的变化导致的工作模式频繁切换而带来的电驱动系统阶跃转矩瞬态响应实际功率损耗,得到整车驱动系统最优效率曲面;
S4:构建考虑前后轴荷转移和附着特性动态变化的整车效率最优的多目标双电机四驱转矩分配控制策略,具体包括:根据步骤1搭建的整车动力学及各部件模型,以及步骤2得到的主动防滑控制策略,驾驶员模型通过驱动/制动踏板位置对总需求转矩进行解析,根据当前的车速、需求转矩,通过步骤3得到的整车驱动系统最优效率曲面进行模式划分;考虑前后轴荷转移的动态变化,对当前模式下的总需求转矩进行前/后转矩分配Tf_req、Tr_req,再与前/后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe以及前/后电机外特性所能提供的最大转矩Tmf_max、Tmr_max做对比,得到前后轴实际输出转矩Tmf、Tmr;通过车辆动力学模型输出当前状态下的车速、前后车轮角速度,然后进行滑移率判断;若当前滑移率小于路面附着系数对应的最优滑移率,则进入下一个时刻,从而起到主动防滑的作用。
2.根据权利要求1所述的双电机四驱电动汽车转矩分配方法,其特征在于,步骤S2中,得出主动防滑控制策略,具体包括:在得到实时路面附着系数的情况下,通过考虑前后轴荷转移的车辆动力性模型,得到前/后轴所能提供的最大附着力对应的输出转矩Tmf_safe、Tmr_safe;再与前/后电机最大峰值转矩Tmf_max、Tmr_max做对比,小的一方作为外特性约束边界,从而对车辆动力性进行动态约束,以起到主动防滑的作用。
3.根据权利要求1所述的双电机四驱电动汽车转矩分配方法,其特征在于,步骤S3中,得到整车驱动系统最优效率曲面,具体包括:根据工况,综合考虑电驱动系统的拖拽损耗、异步电机所产生的励磁损耗以及由于工况的变化导致的工作模式频繁切换而带来的电驱动系统阶跃转矩瞬态响应实际功率损耗模型与边界区域效率差值模型,以经济性和切换频率为目标,分析模式边界带宽的影响因素,进行边界寻优;然后分析不同边界带宽下,每一段车速内,所有工况点的能耗及切换频率进行多目标边界寻优,得到最优的整车驱动系统综合效率曲面。
4.根据权利要求3所述的双电机四驱电动汽车转矩分配方法,其特征在于,所述的边界区域效率差值模型的表达式为:Δη=f(T,n,η1,η2),其中,T为当前总需求转矩,n为转速,η1为前电机效率,η2为后电机效率。
5.根据权利要求4所述的双电机四驱电动汽车转矩分配方法,其特征在于,综合考虑整车动力性约束的动态变化,结合电池、电机和传动系统的效率,得到不同驱动模式下整车驱动系统综合效率;
整车效率取决于车速、加速度、电池效率、电机效率和传动系统效率;
以整车效率为目标函数max ηsys,得到各驱动模式下整车行驶效率最优的前后电机的最佳转矩分配,进而得到整车的模式划分规则;其中,以整车效率为目标函数的约束条件为:
其中,SOC为电荷状态,SOCmin为最小截止电荷状态,Pbat为电池功率,Pbat_max为最大电池放电功率,Tmf、Tmr为前/后电机轴端实际输出转矩,Tmf_max、Tmr_max为前/后电机最大峰值转矩,Tmf_safe、Tmr_safe分别为前/后轴所能提供的最大附着力对应的电机轴端输出转矩,nf、nr分别为前/后电机转速,nf_max、nr_max分别为前/后电机最大转速。
6.根据权利要求5所述的双电机四驱电动汽车转矩分配方法,其特征在于,模式切换的控制策略为:根据不同模式下的整车驱动系统综合效率,得到在不同车速和不同需求转矩下的不同模式划分规则,即各模式的运行区间;然后根据速度和总需求转矩的变化自动选择工作模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110795840.1A CN113547928B (zh) | 2021-07-14 | 2021-07-14 | 一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110795840.1A CN113547928B (zh) | 2021-07-14 | 2021-07-14 | 一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113547928A CN113547928A (zh) | 2021-10-26 |
CN113547928B true CN113547928B (zh) | 2022-11-25 |
Family
ID=78103095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110795840.1A Active CN113547928B (zh) | 2021-07-14 | 2021-07-14 | 一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113547928B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114428494B (zh) * | 2022-01-18 | 2024-02-13 | 上汽通用五菱汽车股份有限公司 | 电动车辆整车控制器pid参数标定方法、设备及介质 |
CN114347805B (zh) * | 2022-01-27 | 2022-12-27 | 奇瑞商用车(安徽)有限公司 | 电动汽车集成车辆控制的双电机驱动系统及其控制方法 |
CN115056656B (zh) * | 2022-08-03 | 2024-12-27 | 中国第一汽车股份有限公司 | 电动车辆的控制方法、装置、电动车辆和存储介质 |
CN115214382A (zh) * | 2022-08-10 | 2022-10-21 | 中国第一汽车股份有限公司 | 一种基于车路协感知融合的纯电动车能耗优化方法、系统和装置 |
CN118907107B (zh) * | 2024-10-09 | 2025-01-07 | 张家港长城汽车研发有限公司 | 车辆扭矩调节方法、装置、电子设备及车辆 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2587666A1 (fr) * | 1985-09-20 | 1987-03-27 | Mitsubishi Electric Corp | Procede pour controler le mode d'entrainement de vehicules electriques |
CN102275580A (zh) * | 2010-06-10 | 2011-12-14 | 福特全球技术公司 | 机动车辆及其控制方法 |
CN103879305A (zh) * | 2014-03-10 | 2014-06-25 | 清华大学 | 用于四轮独立驱动电动车的最大转矩估计驱动防滑算法 |
CN109552312A (zh) * | 2018-11-14 | 2019-04-02 | 吉林大学 | 车辆稳定性模型预测控制方法 |
CN112356685A (zh) * | 2020-11-25 | 2021-02-12 | 重庆大学 | 一种前后独立驱动四驱电动汽车的转矩分配及驱动防滑协调控制方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3589208B2 (ja) * | 2001-08-13 | 2004-11-17 | 日産自動車株式会社 | ハイブリッド車両の駆動装置 |
US8504238B2 (en) * | 2010-04-16 | 2013-08-06 | Ford Global Technologies | Vehicle stability and steerability control via electronic torque distribution |
US8965609B2 (en) * | 2011-12-29 | 2015-02-24 | Kawasaki Jukogyo Kabushiki Kaisha | Electric vehicle |
US20160090005A1 (en) * | 2014-03-10 | 2016-03-31 | Dean Drako | Distributed Torque Generation System and Method of Control |
US11052757B2 (en) * | 2019-04-05 | 2021-07-06 | GM Global Technology Operations LLC | Methods of controlling axle torque distribution |
-
2021
- 2021-07-14 CN CN202110795840.1A patent/CN113547928B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2587666A1 (fr) * | 1985-09-20 | 1987-03-27 | Mitsubishi Electric Corp | Procede pour controler le mode d'entrainement de vehicules electriques |
CN102275580A (zh) * | 2010-06-10 | 2011-12-14 | 福特全球技术公司 | 机动车辆及其控制方法 |
CN103879305A (zh) * | 2014-03-10 | 2014-06-25 | 清华大学 | 用于四轮独立驱动电动车的最大转矩估计驱动防滑算法 |
CN109552312A (zh) * | 2018-11-14 | 2019-04-02 | 吉林大学 | 车辆稳定性模型预测控制方法 |
CN112356685A (zh) * | 2020-11-25 | 2021-02-12 | 重庆大学 | 一种前后独立驱动四驱电动汽车的转矩分配及驱动防滑协调控制方法 |
Non-Patent Citations (1)
Title |
---|
姜涛.前后轴独立驱动电动汽车转矩控制策略研究.《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》.2020, * |
Also Published As
Publication number | Publication date |
---|---|
CN113547928A (zh) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113547928B (zh) | 一种考虑轮胎滑移的双电机四驱电动汽车转矩分配方法 | |
CN104786804B (zh) | 一种车辆及其轮边驱动系统和轮边驱动扭矩分配方法 | |
CN101814876B (zh) | 一种电动汽车后桥双电机驱动控制系统 | |
CN103879307B (zh) | 一种用于电动汽车的后轮独立驱动控制系统及方法 | |
CN103419676B (zh) | 一种用于电动汽车的轮式驱动分层控制系统及方法 | |
CN106864306A (zh) | 一种分布式驱动电动汽车多模式电子差速控制系统 | |
CN107953801A (zh) | 一种全轮毂电机驱动车辆的驱动力控制方法 | |
CN110422053B (zh) | 四轮轮毂电机驱动电动汽车节能控制方法 | |
CN108790940A (zh) | 轮边驱动转向差速控制方法、控制装置、设备及汽车 | |
Cao et al. | All-wheel-drive torque distribution strategy for electric vehicle optimal efficiency considering tire slip | |
CN107117073A (zh) | 一种四轮轮毂电动汽车牵引力控制方法 | |
CN104210383A (zh) | 一种四轮独立驱动电动汽车转矩分配控制方法及系统 | |
CN111746305B (zh) | 线控四轮驱动轮毂电机电动汽车节能控制方法及系统 | |
CN103303157A (zh) | 一种四轮驱动电动汽车转矩分配方法 | |
CN112356685B (zh) | 一种前后独立驱动四驱电动汽车的转矩分配及驱动防滑协调控制方法 | |
CN104477051A (zh) | 一种双驱轴双电机纯电动车驱动电机功率差异化匹配方法 | |
CN105620310B (zh) | 一种三电机混合驱动载货汽车及动力系统参数匹配方法 | |
CN108544913A (zh) | 三电机传动结构纯电动汽车及其转矩分配优化算法 | |
CN105667346A (zh) | 一种三电机混合驱动载货汽车及动力系统参数匹配方法 | |
CN116512934A (zh) | 实现三电机四驱电动汽车能耗优化的扭矩分配控制方法 | |
CN114083995B (zh) | 一种轮毂电机汽车的转矩分配的方法、系统和介质 | |
CN115723590A (zh) | 一种轮毂电机驱动汽车的节能转矩矢量控制方法 | |
CN110341498A (zh) | 一种轮毂电机分布式分时四驱电动汽车底盘构型、四驱电动汽车及控制方法 | |
CN207374100U (zh) | 一种半分布式前后多电机驱动系统 | |
CN116691371A (zh) | 一种分布式驱动车辆动力系统能耗优化方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |