CN114889446A - An off-road vehicle two-direction torque vector distribution method, device and storage medium - Google Patents

An off-road vehicle two-direction torque vector distribution method, device and storage medium Download PDF

Info

Publication number
CN114889446A
CN114889446A CN202210499593.5A CN202210499593A CN114889446A CN 114889446 A CN114889446 A CN 114889446A CN 202210499593 A CN202210499593 A CN 202210499593A CN 114889446 A CN114889446 A CN 114889446A
Authority
CN
China
Prior art keywords
vehicle
wheel
yaw moment
torque
tire
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
CN202210499593.5A
Other languages
Chinese (zh)
Other versions
CN114889446B (en
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.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
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 Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN202210499593.5A priority Critical patent/CN114889446B/en
Publication of CN114889446A publication Critical patent/CN114889446A/en
Application granted granted Critical
Publication of CN114889446B publication Critical patent/CN114889446B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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/00Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles
    • B60L15/20Methods, 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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/00Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles
    • B60L15/32Control or regulation of multiple-unit electrically-propelled vehicles
    • B60L15/38Control or regulation of multiple-unit electrically-propelled vehicles with automatic control
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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
    • B60L2220/00Electrical machine types; Structures or applications thereof
    • B60L2220/40Electrical machine applications
    • B60L2220/42Electrical machine applications with use of more than one motor
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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
    • B60L2220/00Electrical machine types; Structures or applications thereof
    • B60L2220/40Electrical machine applications
    • B60L2220/44Wheel Hub motors, i.e. integrated in the wheel hub
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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/00Control parameters of input or output; Target parameters
    • B60L2240/10Vehicle control parameters
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION 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/00Control parameters of input or output; Target parameters
    • B60L2240/40Drive Train control parameters
    • B60L2240/42Drive Train control parameters related to electric machines
    • B60L2240/423Torque
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/72Electric energy management in electromobility

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Steering Control In Accordance With Driving Conditions (AREA)

Abstract

本发明公开了一种越野车辆两方向转矩矢量分配方法、设备及存储介质,方法包括纵向运动控制方法和横向稳定性控制方法。本发明公开的一种越野车辆两方向转矩矢量分配方法、设备和存储介质,以机动性、通过性、稳定性以及转向灵活性为优化目标,开发了轮毂电脑及驱动层重构优化矢量控制器。本发明提出的纵向控制方法能够有效提升越野车辆在极端起伏地形下的行驶稳定性和通过性;横向控制方法在提高转向灵活性和横摆稳定性的前提下,解决了滑模控制器常见的抖震问题,通过趋近率的自适应调节实现了横摆运动误差跟踪的全局优化。

Figure 202210499593

The invention discloses a two-direction torque vector distribution method, equipment and storage medium for off-road vehicles. The method includes a longitudinal motion control method and a lateral stability control method. The invention discloses a method, device and storage medium for distributing torque vector in two directions of off-road vehicles. Taking maneuverability, passability, stability and steering flexibility as optimization goals, a wheel hub computer and drive layer reconstruction optimization vector control are developed. device. The longitudinal control method proposed by the invention can effectively improve the driving stability and passability of off-road vehicles under extreme undulating terrain; the lateral control method solves the common problems of sliding mode controllers on the premise of improving steering flexibility and yaw stability. For the chattering problem, the global optimization of the yaw motion error tracking is realized through the adaptive adjustment of the approach rate.

Figure 202210499593

Description

一种越野车辆两方向转矩矢量分配方法、设备及存储介质An off-road vehicle two-direction torque vector distribution method, device and storage medium

技术领域technical field

本发明涉及轮毂电机车辆驱动领域,尤其是一种越野车辆两方向转矩矢量分配方法、设备及存储介质。The invention relates to the field of in-wheel motor vehicle driving, in particular to a method, device and storage medium for distributing torque vector in two directions of off-road vehicles.

背景技术Background technique

近年来,得益于集成技术的逐渐成熟,微型控制器等电子器件得到了长足发展。将其应用到电驱动系统开发中,可显著提升车辆的控制精度和动力性能。配备鲁棒性较好的算法于微型控制器,并结合集成度高、传动链短、机械损耗低、转矩响应速度快、动力强劲等诸多优势的轮毂电机车辆具有良好的发展前景。In recent years, thanks to the gradual maturity of integration technology, electronic devices such as microcontrollers have developed by leaps and bounds. Applying it to the development of electric drive system can significantly improve the control precision and dynamic performance of the vehicle. Equipped with a robust algorithm in the microcontroller, and combined with high integration, short transmission chain, low mechanical loss, fast torque response, strong power and many other advantages, the in-wheel motor vehicle has a good development prospect.

越野车辆因其动力性强、通过性高越来越受广大消费者青睐。然而越野车的行驶环境较为复杂,需要具备在各种等级公路、泥泞路和沙地草地等恶劣路面的通过能力,这样对以轮毂电机驱动的越野车机动性和稳定性受到了极大挑战。Off-road vehicles are increasingly favored by consumers because of their strong power and high passability. However, the driving environment of off-road vehicles is complex, and it needs to have the ability to pass on various grades of roads, muddy roads, sandy grass and other harsh roads, which greatly challenges the mobility and stability of off-road vehicles driven by in-wheel motors.

发明内容SUMMARY OF THE INVENTION

有鉴于此,本发明实施例提供一种越野车辆两方向转矩矢量分配方法、设备及存储介质。In view of this, embodiments of the present invention provide a method, device, and storage medium for distributing torque vectors in two directions of an off-road vehicle.

本发明的第一方面提供了一种越野车辆两方向转矩矢量分配方法,包括纵向运动控制方法和横向稳定性控制方法;A first aspect of the present invention provides a two-directional torque vector distribution method for off-road vehicles, including a longitudinal motion control method and a lateral stability control method;

所述纵向运动控制方法具体包括:The longitudinal motion control method specifically includes:

建立轮胎力学模型,根据轮胎力学模型得出车辆的驱动力矩传递方程;Establish a tire mechanics model, and obtain the driving torque transfer equation of the vehicle according to the tire mechanics model;

计算车辆的实时轮端输出力矩,通过力矩矢量预分配,调节车辆的车轮附着利用率;Calculate the real-time wheel end output torque of the vehicle, and adjust the wheel adhesion utilization rate of the vehicle through the torque vector pre-distribution;

设计滑膜控制器,通过滑膜控制器输出的驱动力矩指令跟踪车轮的滑转率误差,以提高纵向运动控制效果;A synovial controller is designed to track the slip rate error of the wheel through the driving torque command output by the synovial controller, so as to improve the longitudinal motion control effect;

所述横向稳定性控制方法具体包括:The lateral stability control method specifically includes:

设计车辆前馈控制方法,得到前馈横摆力矩;Design the vehicle feedforward control method to obtain the feedforward yaw moment;

设计车辆反馈控制方法,得到反馈横摆力矩;Design the vehicle feedback control method to get the feedback yaw moment;

根据前馈横摆力矩和反馈横摆力矩,通过力矩矢量再分配,调节车辆的车轮附着利用率。According to the feedforward yaw moment and the feedback yaw moment, the vehicle's wheel adhesion utilization is adjusted through the redistribution of the torque vector.

进一步地,所述建立轮胎力学模型,根据轮胎力学模型得出车辆的驱动力矩传递方程,具体包括:Further, the tire mechanics model is established, and the driving torque transmission equation of the vehicle is obtained according to the tire mechanics model, which specifically includes:

建立轮胎力学模型,自轮胎模型中得到车辆轮胎的接地长度、纵向力、滑动摩擦和车轮输出力矩;Establish a tire mechanics model, and obtain the ground contact length, longitudinal force, sliding friction and wheel output torque of the vehicle tire from the tire model;

其中接地长度通过如下公式表示:The grounding length is expressed by the following formula:

Figure BDA0003634912840000021
Figure BDA0003634912840000021

纵向力通过如下公式表示:The longitudinal force is expressed by the following formula:

Figure BDA0003634912840000022
Figure BDA0003634912840000022

滑动摩擦通过如下公式表示:Sliding friction is expressed by the following formula:

Figure BDA0003634912840000023
Figure BDA0003634912840000023

车轮输出力矩通过如下公式表示:The wheel output torque is expressed by the following formula:

Figure BDA0003634912840000024
Figure BDA0003634912840000024

式中:z为轮胎平均弹性变形量,vr为轮胎与接触面之间的相对速度,ε为路面附着条件的动态变化,σ0、σ1、σ2分别为轮胎接触面归一化集总刚度、集总阻尼以及粘性相对阻尼,m为整车质量,vx为实际纵向车速,Fz为车辆垂向载荷,g(vr)为滑动摩擦,μc为归一化库仑摩擦系数,μs为静态摩擦系数,vs为描述低速滑动阶段的Stribeck相对速度,Tw为车轮输出力矩,J为车轮转动惯量,w为车轮角速度,σw为黏性转动摩擦系数,rw为车轮有效滚动半径;In the formula: z is the average elastic deformation of the tire, v r is the relative velocity between the tire and the contact surface, ε is the dynamic change of the road adhesion condition, σ 0 , σ 1 , and σ 2 are the normalized set of tire contact surfaces, respectively Total stiffness, lumped damping and viscous relative damping, m is vehicle mass, v x is actual longitudinal vehicle speed, F z is vehicle vertical load, g(v r ) is sliding friction, μ c is normalized Coulomb friction coefficient , μ s is the static friction coefficient, v s is the Stribeck relative velocity describing the low-speed slip stage, Tw is the wheel output torque, J is the wheel rotational inertia, w is the wheel angular velocity, σ w is the viscous rotational friction coefficient, and r w is The effective rolling radius of the wheel;

根据轮胎力学模型结合车辆四轮滑转变化率,得到驱动力矩传递方程;According to the tire mechanics model combined with the four-wheel slip change rate of the vehicle, the driving torque transmission equation is obtained;

Figure BDA0003634912840000025
Figure BDA0003634912840000025

式中:fi为附着系数拟合函数;μi为轮胎附着系数;λi为轮胎滑转率;ax为车辆纵向加速度;wi为轮胎角速度;rw为车轮有效滚动半径;Fxi为轮胎纵向力;Fzi为垂向载荷;Ji为车轮转动惯量;im为减速器传动比,

Figure BDA0003634912840000026
为四轮滑转率变化率,L为拉氏变换。where f i is the fitting function of the adhesion coefficient; μ i is the tire adhesion coefficient; λ i is the tire slip rate; a x is the longitudinal acceleration of the vehicle; w i is the tire angular velocity; r w is the effective rolling radius of the wheel; F xi is the longitudinal force of the tire; F zi is the vertical load; J i is the moment of inertia of the wheel; i m is the transmission ratio of the reducer,
Figure BDA0003634912840000026
is the change rate of the four-wheel slip rate, and L is the Laplace transform.

进一步地,所述计算车辆的实时轮端输出力矩,通过力矩矢量预分配,调节车辆的车轮附着利用率,具体包括:Further, the calculation of the real-time wheel end output torque of the vehicle, through the torque vector pre-allocation, adjusts the wheel adhesion utilization rate of the vehicle, specifically including:

通过以下公式计算实时轮端输出力矩:Calculate the real-time wheel end output torque by the following formula:

Figure BDA0003634912840000031
Figure BDA0003634912840000031

Figure BDA0003634912840000032
Figure BDA0003634912840000032

式中:Ti(k)为实时轮端输出力矩;Tmaxi(k)为实时轮端输出最大输出力矩;ΔTi(k)为驱动防滑调节指令;

Figure BDA0003634912840000033
为加速踏板开度变化率;OSi为稳定性判定标志位;
Figure BDA0003634912840000034
Figure BDA0003634912840000035
的一阶惯性延时系统;τacc为电子踏板、电机峰值电流采样周期,t为采样时间,m为整车质量,g为重力加速度。Where: T i (k) is the real-time wheel end output torque; T maxi (k) is the real-time wheel end output maximum output torque; ΔT i (k) is the drive anti-skid adjustment command;
Figure BDA0003634912840000033
is the rate of change of the accelerator pedal opening; OS i is the stability determination flag;
Figure BDA0003634912840000034
for
Figure BDA0003634912840000035
τ acc is the sampling period of the electronic pedal and motor peak current, t is the sampling time, m is the mass of the vehicle, and g is the acceleration of gravity.

进一步地,所述设计滑膜控制器,通过滑膜控制器输出的驱动力矩指令跟踪车轮的滑转率误差,以提高纵向运动控制效果,具体包括:Further, the designed synovial controller tracks the slip rate error of the wheel through the driving torque command output by the synovial controller, so as to improve the longitudinal motion control effect, specifically including:

设计自适应滑模控制器的滑动面和趋近律:Design the sliding surface and reaching law of the adaptive sliding mode controller:

s(k)=1-DPURi(k)s(k)=1-DPUR i (k)

Figure BDA0003634912840000036
Figure BDA0003634912840000036

式中:s表示滑模面函数;q为趋近速率增益用以调节λi的收敛速度;In the formula: s represents the sliding mode surface function; q is the approach rate gain to adjust the convergence speed of λ i ;

设计增益q的计算式:The formula for calculating the design gain q:

Figure BDA0003634912840000037
Figure BDA0003634912840000037

得到驱动力矩指令Ti(k):Obtain the driving torque command T i (k):

Figure BDA0003634912840000038
Figure BDA0003634912840000038

式中驱动力矩指令Ti(k)受电机最小力矩(Tmin=0)约束。In the formula, the driving torque command T i (k) is constrained by the minimum torque of the motor (T min =0).

进一步地,所述设计车辆前馈控制方法,得到前馈横摆力矩,具体包括:Further, the design of the vehicle feedforward control method to obtain the feedforward yaw moment specifically includes:

根据车辆前轮转向角、车身质心侧偏角、车辆质心到前轴距离、车辆横摆角速度和纵向车速,计算车辆前轮侧偏角;Calculate the front wheel slip angle of the vehicle according to the steering angle of the front wheel of the vehicle, the side slip angle of the center of mass of the vehicle, the distance from the center of mass of the vehicle to the front axle, the yaw rate of the vehicle and the longitudinal vehicle speed;

Figure BDA0003634912840000041
Figure BDA0003634912840000041

根据车辆质心到后轴距离、车辆横摆角速度、纵向车速和车身质心侧偏角,计算车辆的后轮侧偏角;Calculate the rear wheel slip angle of the vehicle according to the distance from the center of mass of the vehicle to the rear axle, the yaw rate of the vehicle, the longitudinal vehicle speed and the side slip angle of the center of mass of the vehicle;

Figure BDA0003634912840000042
Figure BDA0003634912840000042

建立侧偏力方程:Establish the cornering force equation:

Fyi=Ciαi F yi =C i α i

展开得到:Expand to get:

Figure BDA0003634912840000043
Figure BDA0003634912840000043

建立前馈横摆力矩方程:Establish the feedforward yaw moment equation:

Figure BDA0003634912840000044
Figure BDA0003634912840000044

式中:αF、αR为前、后轮侧偏角;δF为前轮转角;LF、LR为车辆质心到前后轴的距离;β为车身质心侧偏角;γ为车辆横摆角速度;Ci为前后轮侧偏刚度,下标i使用FL、FR、RL、RR分别指代左前轮、右前轮、左后轮、右后轮;Fyi为侧偏力;m为整车质量,vx为纵向车速,vy为横向车速;Iz为车辆绕z轴转动惯量;MzFF为前馈横摆力矩;In the formula: α F , α R are the front and rear wheel side slip angles; δ F is the front wheel rotation angle; LF , LR are the distances from the center of mass of the vehicle to the front and rear axles; β is the side slip angle of the body center of mass; γ is the vehicle transverse Swing angular velocity; C i is the cornering stiffness of the front and rear wheels, and the subscript i uses FL, FR, RL, and RR to refer to the left front wheel, right front wheel, left rear wheel, and right rear wheel respectively; F yi is the cornering force; m is the vehicle mass, v x is the longitudinal vehicle speed, v y is the lateral vehicle speed; I z is the moment of inertia of the vehicle around the z-axis; M zFF is the feedforward yaw moment;

稳态工况下

Figure BDA0003634912840000045
近似等于质心侧偏角速度
Figure BDA0003634912840000046
将侧偏力方程与前馈横摆力矩方程联立化简得到前馈横摆力矩表达式:under steady state
Figure BDA0003634912840000045
approximately equal to the centroid sideslip angular velocity
Figure BDA0003634912840000046
Simultaneously simplify the cornering force equation and the feedforward yaw moment equation to obtain the feedforward yaw moment expression:

Figure BDA0003634912840000047
Figure BDA0003634912840000047

所述前馈横摆力矩表达式与车辆前轮转角正相关。The feedforward yaw moment expression is positively correlated with the vehicle front wheel rotation angle.

进一步地,在得到前馈横摆力矩后,还包括:Further, after the feedforward yaw moment is obtained, the method further includes:

使用期望前轮转角替代前轮转角进行前馈横摆力矩计算、以施加差动转向横摆力矩Feedforward yaw moment calculation using desired front wheel angle instead of front wheel angle to apply differential steering yaw moment

δFd=δF(1+K)δ Fd = δ F (1+K)

Figure BDA0003634912840000048
Figure BDA0003634912840000048

式中δFd为车辆期望前轮转角,K为表征车辆转向特性的侧偏稳定性参数。where δFd is the expected front wheel turning angle of the vehicle, and K is the cornering stability parameter that characterizes the steering characteristics of the vehicle.

进一步地,所述设计车辆反馈控制方法,得到反馈横摆力矩,具体包括:Further, the designed vehicle feedback control method to obtain the feedback yaw moment specifically includes:

计算质心侧偏角最大值:Compute the maximum value of the centroid sideslip angle:

Figure BDA0003634912840000051
Figure BDA0003634912840000051

式中,g为重力加速度;μ为路面附着系数;In the formula, g is the acceleration of gravity; μ is the adhesion coefficient of the road surface;

设计协调控制权重:Design coordination control weights:

Figure BDA0003634912840000052
Figure BDA0003634912840000052

式中,Fxmax表示实时的轮端纵向驱动力峰值;m为整车质量;ε1与ε2为ξ的约束系数;In the formula, F xmax represents the real-time peak value of the longitudinal driving force at the wheel end; m is the mass of the vehicle; ε 1 and ε 2 are the constraint coefficients of ξ;

计算跟踪误差:Calculate tracking error:

Figure BDA0003634912840000053
Figure BDA0003634912840000053

其中:βd为期望质心侧偏角,β为质心侧偏角,γcorrect为附加横摆力矩;Where: β d is the expected center of mass side slip angle, β is the side slip angle of the center of mass, and γ correct is the additional yaw moment;

Figure BDA0003634912840000054
Figure BDA0003634912840000054

p(t)为有界不确定扰动;p(t) is a bounded uncertain perturbation;

得到反馈横摆力矩:Get the feedback yaw moment:

Figure BDA0003634912840000055
Figure BDA0003634912840000055

式中,特征参数k1与k2均为正数,sgn为符号函数,k4为控制权重因子。In the formula, the characteristic parameters k 1 and k 2 are both positive numbers, sgn is the sign function, and k 4 is the control weight factor.

进一步地,所述根据前馈横摆力矩和反馈横摆力矩,通过力矩矢量再分配,调节车辆的车轮附着利用率,具体包括:Further, according to the feedforward yaw moment and the feedback yaw moment, the redistribution of the torque vector is used to adjust the wheel adhesion utilization rate of the vehicle, which specifically includes:

将前馈横摆力矩与反馈横摆力矩相加输入至执行层;Add the feedforward yaw moment and the feedback yaw moment to the executive layer;

计算车辆各个车轮中轮毂电机的调节量,完成力矩矢量分配。Calculate the adjustment amount of the in-wheel motor in each wheel of the vehicle to complete the torque vector distribution.

本发明的第二方面公开了一种电子设备,包括处理器以及存储器;A second aspect of the present invention discloses an electronic device, including a processor and a memory;

所述存储器用于存储程序;the memory is used to store programs;

所述处理器执行所述程序实现一种越野车辆两方向转矩矢量分配方法。The processor executes the program to implement a two-direction torque vector distribution method for off-road vehicles.

本发明的第三方面公开了一种计算机可读存储介质,所述存储介质存储有程序,所述程序被处理器执行实现一种越野车辆两方向转矩矢量分配方法。A third aspect of the present invention discloses a computer-readable storage medium, the storage medium stores a program, and the program is executed by a processor to implement a two-direction torque vector distribution method for an off-road vehicle.

本发明的具有如下有益效果:本发明公开的一种越野车辆两方向转矩矢量分配方法、设备和存储介质,以机动性、通过性、稳定性以及转向灵活性为优化目标,开发了轮毂电脑及驱动层重构优化矢量控制器。本发明提出的纵向控制方法能够有效提升越野车辆在极端起伏地形下的行驶稳定性和通过性;横向控制方法在提高转向灵活性和横摆稳定性的前提下,解决了滑模控制器常见的抖震问题,通过趋近率的自适应调节实现了横摆运动误差跟踪的全局优化。The invention has the following beneficial effects: the invention discloses a method, device and storage medium for distributing torque vector in two directions of an off-road vehicle. Taking maneuverability, passability, stability and steering flexibility as optimization goals, a hub computer is developed. And the driver layer reconfigures the optimized vector controller. The longitudinal control method proposed by the invention can effectively improve the driving stability and passability of off-road vehicles under extreme undulating terrain; the lateral control method solves the common problems of sliding mode controllers on the premise of improving steering flexibility and yaw stability. For the chattering problem, the global optimization of the yaw motion error tracking is realized through the adaptive adjustment of the approach rate.

本发明的附加方面和优点将在下面的描述部分中给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the present invention will be presented in the description which follows, in part, which will become apparent from the following description, or may be learned by practice of the invention.

附图说明Description of drawings

为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to illustrate the technical solutions in the embodiments of the present application more clearly, the following briefly introduces the drawings that are used in the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative effort.

图1为本发明一种越野车辆两方向转矩矢量分配方法的主要流程图;1 is a main flow chart of a method for distributing torque vectors in two directions of an off-road vehicle according to the present invention;

图2为本发明一种越野车辆两方向转矩矢量分配方法的车辆参数示意图。FIG. 2 is a schematic diagram of vehicle parameters of a method for distributing torque vectors in two directions of an off-road vehicle according to the present invention.

具体实施方式Detailed ways

为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。In order to make the purpose, technical solutions and advantages of the present application more clearly understood, the present application will be described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application, but not to limit the present application.

本实施例公开了一种越野车辆两方向转矩矢量分配方法,主要流程如图1所示,包括SA.纵向运动控制方法和SB.横向稳定性控制方法;This embodiment discloses a two-directional torque vector distribution method for off-road vehicles. The main process is shown in Figure 1, including SA. longitudinal motion control method and SB. lateral stability control method;

纵向运动控制方法具体包括:The longitudinal motion control methods specifically include:

SA-1.建立轮胎力学模型,根据轮胎力学模型得出车辆的驱动力矩传递方程;SA-1. Establish a tire mechanics model, and obtain the driving torque transfer equation of the vehicle according to the tire mechanics model;

SA-2.计算车辆的实时轮端输出力矩,通过力矩矢量预分配,调节车辆的车轮附着利用率;SA-2. Calculate the real-time wheel end output torque of the vehicle, and adjust the wheel adhesion utilization rate of the vehicle through the torque vector pre-distribution;

SA-3.设计滑膜控制器,通过滑膜控制器输出的驱动力矩指令跟踪车轮的滑转率误差,以提高纵向运动控制效果;SA-3. Design a synovial controller to track the slip rate error of the wheel through the driving torque command output by the synovial controller to improve the longitudinal motion control effect;

所述横向稳定性控制方法具体包括:The lateral stability control method specifically includes:

SB-1.设计车辆前馈控制方法,得到前馈横摆力矩;SB-1. Design the vehicle feedforward control method to obtain the feedforward yaw moment;

SB-2.设计车辆反馈控制方法,得到反馈横摆力矩;SB-2. Design the vehicle feedback control method to obtain the feedback yaw moment;

SB-3.根据前馈横摆力矩和反馈横摆力矩,通过力矩矢量再分配,调节车辆的车轮附着利用率。SB-3. According to the feedforward yaw moment and the feedback yaw moment, through the moment vector redistribution, adjust the wheel adhesion utilization of the vehicle.

下面将具体论述方法各部分的执行方式:The implementation of each part of the method will be discussed in detail below:

SA-1.建立轮胎力学模型,根据轮胎力学模型得出车辆的驱动力矩传递方程:SA-1. Establish a tire mechanics model, and obtain the driving torque transfer equation of the vehicle according to the tire mechanics model:

搭建LuGre轮胎力学模型,自轮胎模型中得到车辆轮胎的接地长度、纵向力、滑动摩擦和车轮输出力矩。The LuGre tire mechanics model was built, and the contact length, longitudinal force, sliding friction and wheel output torque of the vehicle tire were obtained from the tire model.

其中接地长度通过如下公式表示:The grounding length is expressed by the following formula:

Figure BDA0003634912840000071
Figure BDA0003634912840000071

纵向力通过如下公式表示:The longitudinal force is expressed by the following formula:

Figure BDA0003634912840000072
Figure BDA0003634912840000072

滑动摩擦通过如下公式表示:Sliding friction is expressed by the following formula:

Figure BDA0003634912840000073
Figure BDA0003634912840000073

车轮输出力矩通过如下公式表示:The wheel output torque is expressed by the following formula:

Figure BDA0003634912840000074
Figure BDA0003634912840000074

式中:

Figure BDA0003634912840000075
表示轮胎接地长度;z表示轮胎平均弹性变形量;vr表示轮胎与接触面之间的相对速度;ε表征路面附着条件的动态变化;σ0、σ1、σ2分别为轮胎接触面归一化集总刚度、集总阻尼以及粘性相对阻尼;m为整车质量;vx为实际纵向车速;Fz为车辆垂向载荷;g(vr)为滑动摩擦,描述了摩擦力在不同相对速度、滑动条件下的变化;μc为归一化库仑摩擦系数;μs为静态摩擦系数;vs为描述低速滑动阶段的Stribeck相对速度;Tw为车轮输出力矩;J为车轮转动惯量;w为车轮角速度;σw为黏性转动摩擦系数;rw为车轮有效滚动半径。where:
Figure BDA0003634912840000075
Represents the tire contact length; z represents the average elastic deformation of the tire; v r represents the relative speed between the tire and the contact surface ; ε represents the dynamic change of the road adhesion condition ; lumped stiffness, lumped damping and viscous relative damping; m is the mass of the vehicle; v x is the actual longitudinal speed; F z is the vertical load of the vehicle; g(v r ) is the sliding friction, which describes the friction force in different relative Changes under speed and sliding conditions; μ c is the normalized Coulomb friction coefficient; μ s is the static friction coefficient; v s is the Stribeck relative velocity describing the low-speed sliding stage; Tw is the wheel output torque; J is the wheel moment of inertia; w is the wheel angular velocity; σ w is the viscous rotational friction coefficient; r w is the effective rolling radius of the wheel.

SA-2.计算车辆的实时轮端输出力矩,通过力矩矢量预分配,调节车辆的车轮附着利用率:SA-2. Calculate the real-time wheel-end output torque of the vehicle, and adjust the wheel adhesion utilization rate of the vehicle through the torque vector pre-distribution:

将四轮滑转率变化率

Figure BDA0003634912840000076
与LuGre轮胎力学模型联立得到驱动时力矩的传递方程:The four-wheel slip rate change rate
Figure BDA0003634912840000076
Combined with the LuGre tire mechanics model, the torque transfer equation during driving is obtained:

Figure BDA0003634912840000077
Figure BDA0003634912840000077

Figure BDA0003634912840000081
Figure BDA0003634912840000081

式中:fi为附着系数拟合函数;μi为轮胎附着系数;λi为轮胎滑转率;ax为车辆纵向加速度;wi为轮胎角速度;rw为车轮有效滚动半径;Fxi为轮胎纵向力;Fzi为垂向载荷;Ji为车轮转动惯量;im为减速器传动比;where f i is the fitting function of the adhesion coefficient; μ i is the tire adhesion coefficient; λ i is the tire slip rate; a x is the longitudinal acceleration of the vehicle; w i is the tire angular velocity; r w is the effective rolling radius of the wheel; F xi is the longitudinal force of the tire; F zi is the vertical load; J i is the moment of inertia of the wheel; i m is the transmission ratio of the reducer;

上述公式联立,利用拉式变换可得驱动力矩的传递函数:The above formulas are combined, and the transfer function of the driving torque can be obtained by using the pull transformation:

Figure BDA0003634912840000082
Figure BDA0003634912840000082

为了更直观地分析转矩控制系统的性能,可基于式(7)将传递函数表达式变换为多项式以便通过计算传递函数的零极点分析系统稳定性、收敛能力等:In order to analyze the performance of the torque control system more intuitively, the transfer function expression can be transformed into a polynomial based on Equation (7) to analyze the system stability, convergence ability, etc. by calculating the zeros and poles of the transfer function:

Figure BDA0003634912840000083
Figure BDA0003634912840000083

式中:Kdi、τdi均为公式(8)内部参数,L表示拉氏变换(拉普拉斯变换)。In the formula: K di and τ di are both internal parameters of formula (8), and L represents Laplace transform (Laplace transform).

据式(1)、(8)分析可知,由于滑转率满足条件1-λi>0,因此τdi(k)为一动态非负数,即传递函数

Figure BDA0003634912840000084
有负实根,证明力矩传递为稳定系统。According to the analysis of equations (1) and (8), since the slip rate satisfies the condition 1-λ i >0, τ di (k) is a dynamic non-negative number, that is, the transfer function
Figure BDA0003634912840000084
There are negative real roots, which proves that the torque transfer is a stable system.

在实际行车过程中,各轮毂电机驱动力矩是由加速踏板开度、电机峰值电流、驱动防滑调节指令ΔTi等多个离散变量决定和约束的,参见式(9),其更新周期与电子踏板采样周期τacc相等。而考虑到实时轮端输出力矩变化

Figure BDA0003634912840000085
受CAN通讯时滞与功率逆变器响应时间的影响,此处可将
Figure BDA0003634912840000086
简化成延时时长为τw的一阶惯性系统:In the actual driving process, the driving torque of each in-wheel motor is determined and constrained by multiple discrete variables such as the accelerator pedal opening, the motor peak current, and the driving anti-skid adjustment command ΔT i . See equation (9), and its update period is related to the electronic pedal The sampling periods τ acc are equal. Considering the real-time wheel end output torque change
Figure BDA0003634912840000085
Affected by the time delay of CAN communication and the response time of the power inverter, the
Figure BDA0003634912840000086
Simplify to a first-order inertial system with delay time τw :

Figure BDA0003634912840000091
Figure BDA0003634912840000091

Figure BDA0003634912840000092
Figure BDA0003634912840000092

式中:Ti(k)为实时轮端输出力矩;Tmaxi(k)为实时轮端输出最大输出力矩;ΔTi(k)为驱动防滑调节指令;

Figure BDA0003634912840000093
为加速踏板开度变化率;OSi为稳定性判定标志位;
Figure BDA0003634912840000094
Figure BDA0003634912840000095
的一阶惯性延时系统;τacc为电子踏板、电机峰值电流采样周期,t为采样时间,m为整车质量,g为重力加速度。Where: T i (k) is the real-time wheel end output torque; T maxi (k) is the real-time wheel end output maximum output torque; ΔT i (k) is the drive anti-skid adjustment command;
Figure BDA0003634912840000093
is the rate of change of the accelerator pedal opening; OS i is the stability determination flag;
Figure BDA0003634912840000094
for
Figure BDA0003634912840000095
τ acc is the sampling period of the electronic pedal and motor peak current, t is the sampling time, m is the mass of the vehicle, and g is the acceleration of gravity.

基于式(10),可将轮毂电机驱动功率利用率(DPURi)表示为各轮纵向估计力

Figure BDA0003634912840000096
与轮毂电机传递至轮端的驱动力(Twi(k)/rW)的比值,并得到其离散化表达式,并基于各轮稳定性标志位设计了式(11)中的误差校正与积分计算的触发机制。Based on equation (10), the in-wheel motor drive power utilization ratio (DPUR i ) can be expressed as the longitudinal estimated force of each wheel
Figure BDA0003634912840000096
and the ratio of the driving force (T wi (k)/r W ) transmitted by the in-wheel motor to the wheel end, and its discretized expression is obtained, and the error correction and integration in Eq. (11) are designed based on the stability flags of each wheel Calculation trigger mechanism.

Figure BDA0003634912840000097
Figure BDA0003634912840000097

Figure BDA0003634912840000098
Figure BDA0003634912840000098

其中:i和j旨在区分k时刻和k+1不同时刻下的驱动功率利用率公式。Among them: i and j are designed to distinguish the driving power utilization formula at different moments k and k+1.

SA-3.设计滑膜控制器,通过滑膜控制器输出的驱动力矩指令跟踪车轮的滑转率误差,以提高纵向运动控制效果:SA-3. Design a synovial controller to track the slip rate error of the wheel through the driving torque command output by the synovial controller to improve the longitudinal motion control effect:

基于式(11)设计自适应滑模控制器的滑动面和趋近律:The sliding surface and reaching law of the adaptive sliding mode controller are designed based on equation (11):

s(k)=1-DPURi(k) (13)s(k)=1-DPUR i (k) (13)

Figure BDA0003634912840000099
Figure BDA0003634912840000099

式(14)中:s表示滑模面函数;q为趋近速率增益用以调节λi的收敛速度;同时,将符号函数sgn(s(k))的增益设为时变变量(|s(k)|/2),以减小系统在滑动面周围的抖振。In formula (14): s represents the sliding mode surface function; q is the approach rate gain to adjust the convergence speed of λ i ; at the same time, the gain of the sign function sgn(s(k)) is set as a time-varying variable (|s (k)|/2) to reduce the buffeting of the system around the sliding surface.

为满足稳定收敛条件,使滑动模态区存在,设计了增益q的计算式:In order to satisfy the stable convergence condition and make the sliding mode region exist, the calculation formula of the gain q is designed:

Figure BDA0003634912840000101
Figure BDA0003634912840000101

将满足设计要求的离散趋近率代入式(12)可得受电机最小力矩(Tmin=0)约束的驱动力矩指令Ti(k):Substitute the discrete approach rate that meets the design requirements into equation (12) to obtain the driving torque command T i (k) constrained by the minimum torque of the motor (T min =0):

Figure BDA0003634912840000102
Figure BDA0003634912840000102

SB-1.设计车辆前馈控制方法,得到前馈横摆力矩:SB-1. Design the vehicle feedforward control method to obtain the feedforward yaw moment:

前馈环节控制目的是通过施加前馈横摆力矩MzFF缩短车辆达到理想稳态的响应时间,其控制目标包括期望横摆角加速度

Figure BDA0003634912840000103
期望质心侧偏角与角速度
Figure BDA0003634912840000104
The purpose of the feedforward link control is to shorten the response time for the vehicle to reach the ideal steady state by applying the feedforward yaw moment M zFF , and its control objectives include the desired yaw angular acceleration
Figure BDA0003634912840000103
Desired centroid slip angle and angular velocity
Figure BDA0003634912840000104

Figure BDA0003634912840000105
Figure BDA0003634912840000105

式(17)表明理想稳态模型下车辆横摆和侧滑运动应收敛为稳态运动,参考图2所示的简化二自由度车辆动力学模型可求得前馈环节车辆理想的侧向运动、轮胎侧偏角以及横摆运动的计算表达式:Equation (17) shows that the vehicle yaw and sideslip motions should converge to steady-state motions under the ideal steady state model. Referring to the simplified two-degree-of-freedom vehicle dynamics model shown in Figure 2, the ideal sideways motion of the vehicle in the feedforward link can be obtained. , the tire slip angle and the calculation expressions of the yaw motion:

Figure BDA0003634912840000106
Figure BDA0003634912840000106

Figure BDA0003634912840000107
Figure BDA0003634912840000107

Fyi=Ciαi (20)F yi =C i α i (20)

Figure BDA0003634912840000108
Figure BDA0003634912840000108

Figure BDA0003634912840000109
Figure BDA0003634912840000109

式中:αF、αR为前、后轮侧偏角;δF为前轮转向角;LF、LR为车辆质心到前后轴的距离;β为车身质心侧偏角;γ为车辆横摆角速度;Ci为前后轮侧偏刚度;Fyi为侧偏力;m为整车质量,vx为纵向车速,vy为横向车速;Iz为车辆绕z轴转动惯量;下标i可用L、R指代左、右;MzFF为前馈横摆力矩。In the formula: α F , α R are the front and rear wheel side slip angles; δ F is the steering angle of the front wheels; LF , LR are the distances from the vehicle center of mass to the front and rear axles; β is the body center of mass side slip angle; γ is the vehicle yaw rate; C i is the cornering stiffness of the front and rear wheels; F yi is the cornering force; m is the vehicle mass, v x is the longitudinal vehicle speed, v y is the lateral vehicle speed; I z is the moment of inertia of the vehicle around the z-axis; subscript i can use L and R to refer to the left and right; M zFF is the feedforward yaw moment.

稳态工况下,质心侧偏角速度

Figure BDA0003634912840000111
可近似为
Figure BDA0003634912840000112
代入式(21),并与式(22)联立可得横摆和侧滑运动的计算式:Under steady-state conditions, the side-slip angular velocity of the center of mass
Figure BDA0003634912840000111
can be approximated as
Figure BDA0003634912840000112
Substitute into Equation (21) and combine with Equation (22) to obtain the calculation formulas of the yaw and sideslip motions:

Figure BDA0003634912840000113
Figure BDA0003634912840000113

Figure BDA0003634912840000114
Figure BDA0003634912840000114

将式(17)中的目标值代入上式,可得到期望横摆角速度γd的表达式:Substituting the target value in equation (17) into the above equation, the expression for the desired yaw rate γ d can be obtained:

Figure BDA0003634912840000115
Figure BDA0003634912840000115

接着将式(24)与(25)联立得与前轮转角δF正相关的前馈横摆力矩:Then, formulas (24) and (25) are combined to obtain the feedforward yaw moment that is positively related to the front wheel angle δF :

Figure BDA0003634912840000116
Figure BDA0003634912840000116

此外,考虑到低速转向阻力矩较大,而基于轮毂电机驱动系统的差动转向也能起到调整转向增益的作用,因此本发明引入了基于操纵灵活性优化的期望前轮转角δFd替代δF进行前馈力矩MZFF的计算、施加差动转向横摆力矩。In addition, considering that the low-speed steering resistance torque is relatively large, and the differential steering based on the in-wheel motor drive system can also play the role of adjusting the steering gain, the present invention introduces the desired front wheel rotation angle δFd based on the optimization of handling flexibility instead of δ F calculates the feedforward moment M ZFF and applies the differential steering yaw moment.

δFd=δF(1+K) (27)δ Fd = δ F (1+K) (27)

Figure BDA0003634912840000117
Figure BDA0003634912840000117

上式中K为表征车辆转向特性的侧偏稳定性参数。In the above formula, K is the cornering stability parameter that characterizes the steering characteristics of the vehicle.

SB-2.设计车辆反馈控制方法,得到反馈横摆力矩:SB-2. Design the vehicle feedback control method to obtain the feedback yaw moment:

本发明采用离散控制技术实现横向运动的反馈跟踪。首先选择横向运动状态参数为横摆角速度γ与质心侧偏角β,前者可用于调整期望的转向特性,后者表征车辆侧滑运动的剧烈程度与横向稳定裕度。进一步设计了反馈控制的状态跟踪误差以实现横摆运动与横向稳定性的协调控制:The invention adopts discrete control technology to realize feedback tracking of lateral motion. Firstly, the lateral motion state parameters are selected as the yaw rate γ and the center of mass sideslip angle β. The former can be used to adjust the desired steering characteristics, and the latter represents the severity of the vehicle's sideslip motion and lateral stability margin. The state tracking error of feedback control is further designed to realize the coordinated control of yaw motion and lateral stability:

Figure BDA0003634912840000121
Figure BDA0003634912840000121

Figure BDA0003634912840000122
Figure BDA0003634912840000122

式中:g为重力加速度;μ为路面附着系数。In the formula: g is the acceleration of gravity; μ is the adhesion coefficient of the road surface.

将路面附着条件(μ=∑μi/4)作为约束以提高稳定性控制的环境适应性。进一步基于β的相平面分析法设计了协调控制权重ξ:The road surface adhesion condition (μ=Σμ i /4) is used as a constraint to improve the environmental adaptability of the stability control. Further, based on the phase plane analysis method of β, the coordinated control weight ξ is designed:

Figure BDA0003634912840000123
Figure BDA0003634912840000123

式(31)中Fxmax表示实时的轮端纵向驱动力峰值;m为整车质量;ε1与ε2为ξ的约束系数。In formula (31), F xmax represents the real-time peak value of the wheel-end longitudinal driving force; m is the vehicle mass; ε 1 and ε 2 are the constraint coefficients of ξ.

为了消除离散滑模控制固有的抖震效应、优化控制平顺性并提升误差的收敛速度,本发明采用了由连续和非连续两种分量组成的二阶滑模控制器(SOSMC)。In order to eliminate the inherent chattering effect of discrete sliding mode control, optimize the control smoothness and improve the convergence speed of the error, the present invention adopts a second-order sliding mode controller (SOSMC) composed of continuous and discontinuous components.

Figure BDA0003634912840000124
Figure BDA0003634912840000124

式中:特征参数k1与k2均为正数,p(t)表征有界不确定扰动。由轮毂电机驱动系统输出的附加横摆力矩可用于修正车辆姿态:In the formula: the characteristic parameters k 1 and k 2 are both positive numbers, and p(t) represents the bounded uncertain disturbance. The additional yaw moment output by the in-wheel motor drive system can be used to correct the vehicle attitude:

Figure BDA0003634912840000125
Figure BDA0003634912840000125

将式(33)代入跟踪误差e的微分方程可得到

Figure BDA0003634912840000126
的表达式:Substitute equation (33) into the differential equation of tracking error e to get
Figure BDA0003634912840000126
expression:

Figure BDA0003634912840000127
Figure BDA0003634912840000127

其中:in:

Figure BDA0003634912840000128
Figure BDA0003634912840000128

由于

Figure BDA0003634912840000131
均为零,因此包含不确定扰动项的p(t)可简化表示为:because
Figure BDA0003634912840000131
are all zero, so p(t) including the uncertain disturbance term can be simplified as:

Figure BDA0003634912840000132
Figure BDA0003634912840000132

p(t)为有界扰动,因此可引入正常数k3表示|P(t)|的极值:p(t) is a bounded disturbance, so a constant k 3 can be introduced to represent the extreme value of |P(t)|:

|P(t)|≤P≤k3 (37)|P(t)|≤P≤k 3 (37)

与滑模控制器因符号函数的作用在滑模面附近出现振动不同,由于二阶滑模控制趋近率中非连续部分的符号函数增益

Figure BDA0003634912840000133
在逼近滑模面过程中逐渐收敛,即滑模面宽度(SW)逐渐减小,可有效减小抖振:Unlike the sliding mode controller which vibrates near the sliding mode surface due to the action of the sign function, the gain of the sign function due to the discontinuous part of the approach rate of the second-order sliding mode control
Figure BDA0003634912840000133
Gradually converge in the process of approaching the sliding mode surface, that is, the sliding mode surface width (SW) gradually decreases, which can effectively reduce the chattering:

Figure BDA0003634912840000134
Figure BDA0003634912840000134

由式(38)可知,SOSMC的滑模面宽度可收敛于零。然而,在直线行驶或以较大半径转向时,微小的环境噪声和传感器测量噪声使得e(t)在零点附近出现高频窄幅的振荡,使驾驶平顺性发生恶化。另一方面,考虑到越野车辆质心位置较高,在中高速转向工况下车身侧倾、轮荷转移较为明显,内侧轮胎失稳风险将显著提高,而在低速转向时较大的车重将显著提高转向阻力矩。因此,本文引入了控制权重因子k4,在直行或轻微转向、低速转向、中高速转向等不同工况下动态调整驾驶平顺性、转向灵活性、横向稳定性等三种控制目标的优化权重:It can be seen from equation (38) that the width of the sliding mode surface of SOSMC can converge to zero. However, when driving in a straight line or turning with a large radius, the small ambient noise and sensor measurement noise cause e(t) to oscillate with a high frequency and a narrow amplitude near the zero point, which deteriorates driving comfort. On the other hand, considering the high position of the center of mass of off-road vehicles, the body roll and wheel load transfer are more obvious under medium and high-speed steering conditions, and the risk of inner tire instability will be significantly increased. Significantly increases steering resistance torque. Therefore, this paper introduces a control weight factor k 4 to dynamically adjust the optimization weights of three control objectives such as driving smoothness, steering flexibility, and lateral stability under different working conditions such as straight or slight steering, low-speed steering, and medium-high-speed steering:

Figure BDA0003634912840000135
Figure BDA0003634912840000135

Figure BDA0003634912840000136
Figure BDA0003634912840000136

如式(39)所示,与方向盘角度、纵向车速相关的两个分段函数描述了车辆的侧向运动状态。其中车速被分为了低、中、高三个连续的速域,在低速和高速工况下,SOSMC的趋近律相对较快;另一方面,在横向稳定裕度较高的工况下,k4设置为零,以确保横摆力矩不会干扰主观操作意图,控制趋近率也将随方向盘角度幅值|δ|的增加而加快;计算得到的k4为一不大于3/2的非负数。As shown in equation (39), two piecewise functions related to steering wheel angle and longitudinal vehicle speed describe the lateral motion state of the vehicle. Among them, the vehicle speed is divided into three continuous speed domains: low, medium and high. Under the low speed and high speed conditions, the reaching law of SOSMC is relatively fast; on the other hand, under the conditions with high lateral stability margin, k 4 is set to zero to ensure that the yaw moment will not interfere with the subjective operation intention, and the control approach rate will also speed up with the increase of the steering wheel angle amplitude |δ| ; negative number.

SB-3.根据前馈横摆力矩和反馈横摆力矩,通过力矩矢量再分配,调节车辆的车轮附着利用率:SB-3. According to the feedforward yaw moment and the feedback yaw moment, through the moment vector redistribution, adjust the wheel adhesion utilization of the vehicle:

将前馈与反馈横摆力矩相加输入至执行层,并以各轮附着利用率协调优化为控制目标,根据各轮附着极限的观测值得到式(41)中各个轮毂电机的调节量ΔTmi,实现力矩矢量再分配、响应直接横摆力矩指令。The feedforward and feedback yaw moment are added and input to the execution layer, and the coordination and optimization of the adhesion utilization rate of each wheel is taken as the control objective, and the adjustment amount ΔT mi of each in-wheel motor in Eq. , to realize the redistribution of the torque vector and to respond to the direct yaw moment command.

Figure BDA0003634912840000141
Figure BDA0003634912840000141

本发明实施例还公开了一种计算机程序产品或计算机程序,该计算机程序产品或计算机程序包括计算机指令,该计算机指令存储在计算机可读存储介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令,处理器执行该计算机指令,使得该计算机设备执行图1所示的方法。The embodiment of the present invention also discloses a computer program product or computer program, where the computer program product or computer program includes computer instructions, and the computer instructions are stored in a computer-readable storage medium. A processor of the computer device can read the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions to cause the computer device to perform the method shown in FIG. 1 .

在一些可选择的实施例中,在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如,取决于所涉及的功能/操作,连续示出的两个方框实际上可以被大体上同时地执行或所述方框有时能以相反顺序被执行。此外,在本发明的流程图中所呈现和描述的实施例以示例的方式被提供,目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的,其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。In some alternative implementations, the functions/operations noted in the block diagrams may occur out of the order noted in the operational diagrams. For example, two blocks shown in succession may, in fact, be executed substantially concurrently or the blocks may sometimes be executed in the reverse order, depending upon the functionality/operations involved. Furthermore, the embodiments presented and described in the flowcharts of the present invention are provided by way of example in order to provide a more comprehensive understanding of the technology. The disclosed methods are not limited to the operations and logic flows presented herein. Alternative embodiments are contemplated in which the order of the various operations are altered and in which sub-operations described as part of larger operations are performed independently.

在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。In the description of this specification, description with reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples", etc., mean specific features described in connection with the embodiment or example , structure, material or feature is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。Although embodiments of the present invention have been shown and described, it will be understood by those of ordinary skill in the art that various changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, The scope of the invention is defined by the claims and their equivalents.

以上是对本发明的较佳实施进行了具体说明,但本发明并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the described embodiments, and those skilled in the art can also make various equivalent deformations or replacements on the premise of not violating the spirit of the present invention, These equivalent modifications or substitutions are all included within the scope defined by the claims of the present application.

Claims (10)

1. A two-direction torque vector distribution method for an off-road vehicle is characterized by comprising a longitudinal motion control method and a transverse stability control method;
the longitudinal motion control method specifically comprises the following steps:
building a tire mechanics model, and obtaining a driving torque transfer equation of the vehicle according to the tire mechanics model;
calculating the real-time wheel end output torque of the vehicle, and adjusting the wheel attachment utilization rate of the vehicle through torque vector pre-distribution;
designing a sliding mode controller, and tracking the slip rate error of the wheel through a driving torque command output by the sliding mode controller so as to improve the control effect of longitudinal motion;
the lateral stability control method specifically comprises the following steps:
designing a vehicle feedforward control method to obtain a feedforward yaw moment;
designing a vehicle feedback control method to obtain a feedback yaw moment;
and adjusting the wheel attachment utilization rate of the vehicle through moment vector redistribution according to the feedforward yaw moment and the feedback yaw moment.
2. A method of distributing a two-directional torque vector for an off-road vehicle as in claim 1, wherein the establishing a tire mechanics model from which to derive a drive torque transfer equation for the vehicle comprises:
building a tire mechanical model, and obtaining the ground contact length, longitudinal force, sliding friction and wheel output torque of the vehicle tire from the tire model;
wherein the grounding length is expressed by the following formula:
Figure FDA0003634912830000011
the longitudinal force is expressed by the following formula:
Figure FDA0003634912830000012
the sliding friction is expressed by the following formula:
Figure FDA0003634912830000013
the wheel output torque is expressed by the following equation:
Figure FDA0003634912830000014
in the formula: z is the average elastic deformation of the tire, v r Is the relative speed between the tyre and the contact surface,. epsilon.is the dynamic variation of the road adhesion conditions,. sigma 0 、σ 1 、σ 2 Normalized lumped stiffness, lumped damping and viscous relative damping of a tire contact surface respectively, m is the mass of the whole vehicle, v x For actual longitudinal vehicle speed, F z For vertical loading of the vehicle, g (v) r ) For sliding friction, mu c To normalize the Coulomb coefficient of friction, μ s Is the static coefficient of friction, v s To describe the Stribeck relative velocity, T, for the low speed slip phase w Is the wheel output torque, J is the wheel moment of inertia, w is the wheel angular velocity, σ w Is a viscous coefficient of rotational friction, r w Is the effective rolling radius of the wheel;
obtaining a driving torque transfer equation according to a tire mechanics model and the vehicle four-wheel slip change rate;
Figure FDA0003634912830000021
in the formula: f. of i Fitting a function to the adhesion coefficient; mu.s i Is the tire adhesion coefficient; lambda i Is the tire slip ratio; a is a x Is the vehicle longitudinal acceleration; w is a i Is the tire angular velocity; r is w Is the effective rolling radius of the wheel; f xi For longitudinal forces of the tyre;F zi Is a vertical load; j. the design is a square i Is the rotational inertia of the wheel; i.e. i m In order to realize the transmission ratio of the speed reducer,
Figure FDA0003634912830000022
the four wheel slip rate change rate, L is the Laplace transform.
3. The method as claimed in claim 2, wherein the calculating of the real-time wheel-end output torque of the vehicle, the adjusting of the wheel attachment utilization of the vehicle by the pre-distribution of the torque vector, specifically comprises:
calculating the real-time wheel end output torque by the following formula:
Figure FDA0003634912830000023
Figure FDA0003634912830000024
in the formula: t is i (k) Outputting torque for the real-time wheel end; t is maxi (k) Outputting the maximum output torque for the real-time wheel end; delta T i (k) Adjusting instructions for driving the antiskid devices;
Figure FDA0003634912830000025
is the accelerator pedal opening change rate; OS i Determining a flag bit for stability;
Figure FDA0003634912830000026
is composed of
Figure FDA0003634912830000027
The first-order inertial delay system of (1); tau is acc The method is characterized in that the sampling period of peak current of an electronic pedal and a motor is shown, t is sampling time, m is the mass of the whole vehicle, and g is gravity acceleration.
4. A method as claimed in claim 3, wherein the slip controller is configured to track slip error of the wheels via drive torque commands output from the slip controller to improve longitudinal motion control, and comprises:
designing a sliding surface and an approach law of the adaptive sliding mode controller:
s(k)=1-DPUR i (k)
Figure FDA0003634912830000031
in the formula: s represents a sliding mode surface function; q is the approximate rate gain to adjust lambda i The convergence rate of (2);
designing a calculation formula of the gain q:
Figure FDA0003634912830000032
obtaining a driving torque command T i (k):
Figure FDA0003634912830000033
In which the drive torque command T is i (k) Minimum torque (T) of motor min 0) constraint.
5. The method for allocating two-directional torque vectors of an off-road vehicle according to claim 1, wherein the method for designing the vehicle feed-forward control to obtain the feed-forward yaw moment specifically comprises:
calculating the side slip angle of the front wheel of the vehicle according to the steering angle of the front wheel of the vehicle, the side slip angle of the mass center of the vehicle body, the distance from the mass center of the vehicle to a front shaft, the yaw velocity of the vehicle and the longitudinal speed of the vehicle;
Figure FDA0003634912830000034
calculating the side slip angle of the rear wheel of the vehicle according to the distance from the mass center of the vehicle to the rear axle, the yaw angular velocity of the vehicle, the longitudinal speed and the mass center side slip angle of the vehicle body;
Figure FDA0003634912830000035
establishing a lateral deviation force equation:
F yi =C i α i
unfolding to obtain:
Figure FDA0003634912830000036
establishing a feedforward yaw moment equation:
Figure FDA0003634912830000041
in the formula: alpha is alpha F 、α R Is the side deflection angle of the front wheel and the rear wheel; delta F Is a front wheel corner; l is F 、L R The distance from the center of mass of the vehicle to the front and rear axes; beta is the vehicle body mass center slip angle; gamma is the yaw velocity of the vehicle; c i For front and rear wheel side deflection rigidity, subscript i respectively refers to a left front wheel, a right front wheel, a left rear wheel and a right rear wheel by using FL, FR, RL and RR; f yi Is the lateral bias force; m is the total vehicle mass, v x For longitudinal vehicle speed, v y Is the transverse vehicle speed; I.C. A z Rotating inertia around a z-axis for the vehicle; m zFF A feedforward yaw moment;
under steady state operating conditions
Figure FDA0003634912830000042
/v x Approximately equal to centroid slip angular velocity
Figure FDA0003634912830000043
And (3) simultaneously simplifying a lateral deviation force equation and a feedforward yaw moment equation to obtain a feedforward yaw moment expression:
Figure FDA0003634912830000044
the feed-forward yaw moment expression is positively correlated with the vehicle front wheel steering angle.
6. An off-road vehicle two-way torque vector allocation method as claimed in claim 5, further comprising, after obtaining a feed forward yaw moment:
feedforward yaw moment calculation using desired front wheel turning angle instead of front wheel turning angle to apply differential steering yaw moment
δ Fd =δ F (1+K)
Figure FDA0003634912830000045
In the formula of Fd The desired front wheel angle for the vehicle, and K is a yaw stability parameter that characterizes the steering characteristics of the vehicle.
7. An off-road vehicle two-way torque vector distribution method as claimed in claim 6, wherein said design vehicle feedback control method, resulting in a feedback yaw moment, specifically comprises:
calculating the maximum value of the centroid slip angle:
Figure FDA0003634912830000046
wherein g is the acceleration of gravity; mu is the road surface adhesion coefficient;
designing a coordination control weight:
Figure FDA0003634912830000051
in the formula, F xmax Representing a real-time wheel end longitudinal drive force peak; m is the mass of the whole vehicle; epsilon 1 And epsilon 2 A constraint coefficient of ξ;
calculating a tracking error:
Figure FDA0003634912830000052
wherein: beta is a d To obtain a desired centroid slip angle, β is the centroid slip angle, γ correct An additional yaw moment;
Figure FDA0003634912830000053
p (t) is a bounded uncertain disturbance;
obtaining a feedback yaw moment:
Figure FDA0003634912830000054
in the formula, the characteristic parameter k 1 And k is 2 Are all positive numbers, sgn is a sign function, k 4 To control the weighting factors.
8. An off-road vehicle two-way torque vector distribution method as claimed in claim 7, wherein the adjusting of the wheel attachment utilization of the vehicle by torque vector redistribution according to the feed-forward yaw moment and the feedback yaw moment comprises:
adding the feedforward yaw moment and the feedback yaw moment and inputting the sum to an execution layer;
and calculating the adjustment quantity of the hub motor in each wheel of the vehicle to complete the distribution of the torque vector.
9. An electronic device comprising a processor and a memory;
the memory is used for storing programs;
the processor executing the program realizes the method according to any one of claims 1-8.
10. A computer-readable storage medium, characterized in that the storage medium stores a program, which is executed by a processor to implement the method according to any one of claims 1-8.
CN202210499593.5A 2022-05-09 2022-05-09 A method, device and storage medium for allocating torque vectors in two directions of off-road vehicles Active CN114889446B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210499593.5A CN114889446B (en) 2022-05-09 2022-05-09 A method, device and storage medium for allocating torque vectors in two directions of off-road vehicles

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210499593.5A CN114889446B (en) 2022-05-09 2022-05-09 A method, device and storage medium for allocating torque vectors in two directions of off-road vehicles

Publications (2)

Publication Number Publication Date
CN114889446A true CN114889446A (en) 2022-08-12
CN114889446B CN114889446B (en) 2024-06-25

Family

ID=82721503

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210499593.5A Active CN114889446B (en) 2022-05-09 2022-05-09 A method, device and storage medium for allocating torque vectors in two directions of off-road vehicles

Country Status (1)

Country Link
CN (1) CN114889446B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230044869A1 (en) * 2021-07-22 2023-02-09 GM Global Technology Operations LLC Vehicle actuation commands to affect transient handling
CN115723590A (en) * 2022-12-13 2023-03-03 吉林大学 Energy-saving torque vector control method for hub motor driven automobile

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2290884A (en) * 1994-06-27 1996-01-10 Fuji Heavy Ind Ltd Torque distribution control system
EP1207090A2 (en) * 2000-11-20 2002-05-22 Toyota Jidosha Kabushiki Kaisha Vehicle motion control device and method
US20110218700A1 (en) * 2010-03-05 2011-09-08 Honda Motor Co., Ltd. Vehicle motion control apparatus
CN208544105U (en) * 2018-07-12 2019-02-26 上海汽车集团股份有限公司 The automobile-used electric drive axle of integrated double-rotor machine
CN113147422A (en) * 2021-03-30 2021-07-23 江苏大学 Direct yaw moment control system and method for distributed driving electric automobile

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2290884A (en) * 1994-06-27 1996-01-10 Fuji Heavy Ind Ltd Torque distribution control system
EP1207090A2 (en) * 2000-11-20 2002-05-22 Toyota Jidosha Kabushiki Kaisha Vehicle motion control device and method
US20110218700A1 (en) * 2010-03-05 2011-09-08 Honda Motor Co., Ltd. Vehicle motion control apparatus
CN208544105U (en) * 2018-07-12 2019-02-26 上海汽车集团股份有限公司 The automobile-used electric drive axle of integrated double-rotor machine
CN113147422A (en) * 2021-03-30 2021-07-23 江苏大学 Direct yaw moment control system and method for distributed driving electric automobile

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张一西;马建;唐自强;刘晓东;贺伊琳;: "分布式驱动电动汽车直接横摆力矩控制策略与仿真验证", 甘肃农业大学学报, no. 04, 15 August 2018 (2018-08-15) *
皮大伟;陈南;张丙军;: "基于主动制动的车辆稳定性系统最优控制策略", 农业机械学报, no. 11, 25 November 2009 (2009-11-25) *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230044869A1 (en) * 2021-07-22 2023-02-09 GM Global Technology Operations LLC Vehicle actuation commands to affect transient handling
US11724739B2 (en) * 2021-07-22 2023-08-15 GM Global Technology Operations LLC Vehicle actuation commands to affect transient handling
CN115723590A (en) * 2022-12-13 2023-03-03 吉林大学 Energy-saving torque vector control method for hub motor driven automobile

Also Published As

Publication number Publication date
CN114889446B (en) 2024-06-25

Similar Documents

Publication Publication Date Title
CN112572411B (en) A vehicle chassis coordinated control method and system considering tire cornering characteristics
CN112677992B (en) Path tracking optimization control method for distributed driving unmanned vehicle
CN111391822B (en) A collaborative control method for vehicle lateral and longitudinal stability under extreme working conditions
CN111055921A (en) A data-driven model predictive control method for four-wheel steering
CN109552312A (en) Intact stability model predictive control method
CN107215329B (en) A distributed drive electric vehicle lateral stability control method based on ATSM
CN109017759B (en) Desired path vehicle yaw control method
CN108107732A (en) Active front wheel steering and the united Vehicle Stability Control method of direct yaw moment
CN114889446B (en) A method, device and storage medium for allocating torque vectors in two directions of off-road vehicles
CN114801781B (en) A Fuzzy-Sliding Mode Composite Control System of Four-Wheel Drive AFS/DYC Integrated Control System
CN113221257B (en) Vehicle lateral and longitudinal stability control method under extreme conditions considering control area
CN111483467A (en) Vehicle control method and device
CN107600176B (en) A kind of intelligent vehicle active steering control method based on extension control theory
CN113650621B (en) State Parameter Estimation Method for Distributed Drive Electric Vehicles for Complex Operating Conditions
Lai et al. Analysis of vehicle driving stability based on longitudinal-lateral and vertical unified dynamics model
CN114771503A (en) Automobile transverse stability control method based on nonsingular terminal sliding mode control
CN113147422A (en) Direct yaw moment control system and method for distributed driving electric automobile
CN114572191A (en) Independently-driven electric automobile trajectory tracking and stability integrated control method
Li et al. Adaptive sliding mode control of lateral stability of four wheel hub electric vehicles
CN116552547A (en) Additional yaw moment real-time control method for distributed driving electric automobile
CN113830094A (en) An adaptive fusion and compensation method for vehicle centroid sideslip angle considering multi-source input information
CN114734983A (en) Distributed electric vehicle stability control method based on stable domain
CN115837843A (en) Direct yaw moment coordinated steering control method for distributed driving electric automobile
CN118323258A (en) On-the-spot steering control method, device and vehicle based on driving intention
CN114684109A (en) Electric automobile yaw stability control method and device

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