CN110395263A - 基于mfac算法的液罐车侧倾稳定性控制方法、系统 - Google Patents

基于mfac算法的液罐车侧倾稳定性控制方法、系统 Download PDF

Info

Publication number
CN110395263A
CN110395263A CN201910756801.3A CN201910756801A CN110395263A CN 110395263 A CN110395263 A CN 110395263A CN 201910756801 A CN201910756801 A CN 201910756801A CN 110395263 A CN110395263 A CN 110395263A
Authority
CN
China
Prior art keywords
tank truck
moment
tank
yaw velocity
control
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
CN201910756801.3A
Other languages
English (en)
Other versions
CN110395263B (zh
Inventor
郑雪莲
任园园
李显生
赵兰
崔晓彤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN201910756801.3A priority Critical patent/CN110395263B/zh
Publication of CN110395263A publication Critical patent/CN110395263A/zh
Application granted granted Critical
Publication of CN110395263B publication Critical patent/CN110395263B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W30/00Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
    • B60W30/02Control of vehicle driving stability
    • B60W30/04Control of vehicle driving stability related to roll-over prevention
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2300/00Indexing codes relating to the type of vehicle
    • B60W2300/12Trucks; Load vehicles
    • B60W2300/125Heavy duty trucks
    • B60W2300/128Silo or fluid transporting vehicles
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2520/00Input parameters relating to overall vehicle dynamics
    • B60W2520/12Lateral speed
    • B60W2520/125Lateral acceleration

Landscapes

  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Vehicle Body Suspensions (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及一种基于MFAC算法的液罐车侧倾稳定性控制方法、系统,基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定液罐车侧倾稳定性的影响因素包括横摆角速度;针对液罐车难以建模、模型通用性差、模型准确度低的问题,基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量,省去了建模过程,而仅利用控制输入和输出进行车辆侧倾控制,适用于各种形状罐体的液罐车,控制效果好,对于保障液罐车的道路运输安全具有十分重要的意义。

Description

基于MFAC算法的液罐车侧倾稳定性控制方法、系统
技术领域
本发明涉及机动车辆行驶安全技术领域,尤其涉及一种基于MFAC算法的液罐车侧倾稳定性控制方法、系统。
背景技术
液罐车是流体货物道路运输的主要载体,在生产生活中有着广泛应用。
与固体货物不同,流体货物具有易于流动的性质,在车辆运行状态发生改变时(如加减速、转弯等),货物在外力作用下冲击罐体壁面,增加了车辆惯性力,降低了车辆的制动和行驶稳定性。
罐体内多安装横向防波板以抑制纵向液体晃动的影响。由于纵向防波板的设计和安装难度大,其较为少见,致使液罐车侧向稳定性问题突出。在侧向液体晃动的影响下,罐车相比普通载货汽车更易发生侧翻事故。因此,采用主动控制的方法提高液罐车侧倾稳定性很有必要,以此来保障道路运输安全。
液罐车是一个具有无限多自由度的流固耦合动力学系统,其动力学特性受到液体冲击的影响,而液体冲击又受到罐体形状、罐内分仓数量及其结构、罐内货物充装率以及货物物理属性等的影响。因此,液罐车的动力学建模困难、模型通用性差、模型准确性低。对液罐车的主动控制不适宜采用基于模型的控制方法。
发明内容
为了解决上述技术问题,本发明的目的在于提供一种基于MFAC算法的液罐车侧倾稳定性控制方法、系统。
根据本发明的一个方面,提供了一种基于MFAC算法的液罐车侧倾稳定性控制方法,包括以下步骤:
基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;
基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
进一步的,所述MFAC算法包括基于全格式动态线性化的MFAC算法。
进一步的,液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系为:
其中,
αy为液罐车的车身侧向加速度;
V为液罐车行驶速度;
r为液罐车横摆角速度;
为液罐车质心侧偏角速度。
进一步的,基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,包括:
基于MFAC算法通过液罐车动力学系统的横摆角速度输出变化量、液罐车动力学系统的伪偏导数矩阵、滑动窗口内的系统输出和理论控制输入变化量矩阵之间的第三预设关系一,液罐车动力学系统的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二,液罐车动力学系统的理论控制输入、第二步长因子、液罐车动力学系统的伪偏导数矩阵、车辆理想横摆角速度、车辆实际横摆角速度、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三,液罐车动力学系统的实际控制输入(附加横摆力矩)、液罐车动力学系统的理论控制输入之间的第三预设关系四,对行驶状态下的液罐车进行侧倾稳定性控制。
其中,液罐车动力学系统(k+1)时刻的横摆角速度输出变化量、液罐车动力学系统k时刻的伪偏导数矩阵、滑动时间窗口[k-ny+1,k]内的系统横摆角速度输出变化量和滑动时间窗口[k-nu+1,k]内的理论控制输入变化量所组成的矩阵之间的第三预设关系一为:
其中,
△r(k+1)为液罐车动力学系统在(k+1)时刻的横摆角速度变化量;
是液罐车动力学系统的伪偏导数矩阵。其中,的维度为1×1,的维度为1×1。的维度为(ny+nu)×1;
是滑动窗口内的系统输出和理论控制输入变化量矩阵,其维度为(ny+nu)×1;
ny为控制输入线性化长度,nu为输出线性化长度。
其中,液罐车动力学系统k时刻的伪偏导数估计值、液罐车动力学系统(k-1)时刻的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二为:
其中,
是液罐车动力学系统k时刻的伪偏导数估计值;
是液罐车动力学系统(k-1)时刻的伪偏导数估计值;
η是第一步长因子;
μ是惩罚因子;
r(k)是液罐车动力学系统在k时刻的横摆角速度输出;
r(k-1)是液罐车动力学系统在(k-1)时刻的横摆角速度输出;
是滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵。
其中,液罐车动力学系统k时刻的理论控制输入、液罐车动力学系统(k-1)时刻的理论控制输入、第二步长因子、液罐车动力学系统k时刻的伪偏导数矩阵、液罐车动力学系统(k+1)时刻的理想横摆角速度输出、液罐车动力学系统在k时刻的实际横摆角速度输出、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三为:
其中,
u(k)为k时刻的系统理论控制输入;
u(k-1)为(k-1)时刻的系统理论控制输入;
λ为权重因子;
ρ为第二步长因子;
r*(k+1)为液罐车动力学系统在(k+1)时刻的理想横摆角速度输出;
r(k)为液罐车动力学系统在k时刻的实际横摆角速度输出;
△r为滑动窗口内的系统横摆角速度输出变化量;
△u为滑动窗口内的系统理论控制输入变化量;
为伪偏导数矩阵中的第ny+1个子矩阵;
为系统伪偏导数矩阵中的第i个子矩阵。
其中,液罐车动力学系统k时刻的实际控制输入(附加横摆力矩)与理论控制输入之间的第三预设关系四为:
U(k)=amp×u(k),
其中:
U(k)为液罐车动力学系统在k时刻的实际控制输入量(附加横摆力矩);
amp为扩充系数。
根据本发明的另一个方面,提供了一种基于MFAC算法的液罐车侧倾稳定性控制系统,包括:
液罐车侧倾稳定性影响因素确定单元,配置用于基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;
液罐车侧倾稳定性控制单元,配置用于基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
进一步的,所述MFAC算法包括基于全格式动态线性化的MFAC算法。
与现有技术相比,本发明具有以下有益效果:
1、本发明示例的基于MFAC算法的液罐车侧倾稳定性控制方法,基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;针对液罐车难以建模、模型通用性差、模型准确度低的问题,基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量,省去了建模过程,而仅利用系统的控制输入和输出进行车辆侧倾控制,适用于装配各种形状罐体的液罐车,控制效果好,这对于保障液罐车的道路运输安全具有十分重要的意义。
2、本发明示例的基于MFAC算法的液罐车侧倾稳定性控制系统,组成简单,通过各个组成系统及单元之间相互配合,基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,省去了建模过程,而仅利用系统的控制输入和输出进行车辆侧倾控制,控制效果好,适用于装配各种形状罐体的液罐车,可有效保障液罐车的道路运输安全。
附图说明
图1为实施例一基于MFAC算法的液罐车侧倾控制方法流程图;
图2为罐体为圆柱、罐内充液比为0.6时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图3为罐体为圆柱、罐内充液比为0.6时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的质心侧偏角示意图;
图4为罐体为圆柱、罐内充液比为0.6时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的横摆角速度示意图;
图5为罐体为圆柱、罐内充液比为0.6时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的车身侧倾角示意图;
图6为罐体为圆柱、罐内充液比为0.6时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的液体振荡摆角示意图;
图7为罐体为圆柱、罐内充液比为0.6时,液罐车CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的附加横摆力矩的示意图;
图8为罐体为椭圆柱,长短轴之比为1.5时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图9为罐体为椭圆柱,长短轴之比为1.5时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的质心侧偏角示意图;
图10为罐体为椭圆柱,长短轴之比为1.5时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的横摆角速度示意图;
图11为罐体为椭圆柱,长短轴之比为1.5时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的车身侧倾角示意图;
图12为罐体为椭圆柱,长短轴之比为1.5时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的液体振荡摆角示意图;
图13为罐体为椭圆柱,长短轴之比为1.5时,液罐车CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的附加横摆力矩的示意图;
图14为罐体为椭圆柱,长短轴之比为2.0时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图15为罐体为椭圆柱,长短轴之比为2.0时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的质心侧偏角示意图;
图16为罐体为椭圆柱,长短轴之比为2.0时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的横摆角速度示意图;
图17为罐体为椭圆柱,长短轴之比为2.0时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的车身侧倾角示意图;
图18为罐体为椭圆柱,长短轴之比为2.0时,液罐车在无控制与CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的液体振荡摆角示意图;
图19为罐体为椭圆柱,长短轴之比为2.0时,液罐车CFDL-MFAC、PFDL-MFAC和FFDL-MAFC控制下的附加横摆力矩的示意图;
图20为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图21为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的质心侧偏角示意图;
图22为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的横摆角速度示意图;
图23为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的车身侧倾角示意图;
图24为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的液体振荡摆角示意图;
图25为罐体为圆柱、罐内充液比为0.6,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的附加横摆力矩的示意图;
图26为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图27为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的质心侧偏角示意图;
图28为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的横摆角速度示意图;
图29为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的车身侧倾角示意图;
图30为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的液体振荡摆角示意图;
图31为罐体为椭圆柱,长短轴之比为1.5,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的附加横摆力矩的示意图;
图32为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的后轴横向载荷转移率示意图;
图33为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的质心侧偏角示意图;
图34为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的横摆角速度示意图;
图35为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的车身侧倾角示意图;
图36为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的液体振荡摆角示意图;
图37为罐体为椭圆柱,长短轴之比为2.0,ny=nu=3、ny=nu=4、ny=nu=5时,液罐车在FFDL-MAFC控制下的附加横摆力矩的示意图。
具体实施方式
为了更好的了解本发明的技术方案,下面结合具体实施例、说明书附图对本发明作进一步说明。
实施例一:
本实施例提供了一种基于MFAC算法的液罐车侧倾稳定性控制方法,如图1所示,包括以下步骤:
S1、基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度。
液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系为:
式中,αy为液罐车的车身侧向加速度;V为液罐车行驶速度;r为液罐车横摆角速度;为液罐车质心侧偏角速度。
一般地,货车车身侧向加速度不超过0.3g时,其侧倾稳定性可以得到保证。当车辆做稳态圆周运动时,质心侧偏角速度为0,此时车辆的侧向加速度由其横摆角速度决定。
因此,可以通过控制车辆的横摆角速度来达到侧倾稳定性控制的目的,通过差动制动控制策略实现液罐车的侧倾控制。控制输入量为附加横摆力矩U,被控量为罐车横摆角速度r。控制系统为单输入(附加横摆力矩U)、单输出(横摆角速度r)系统。
S2、基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
MFAC有三种控制器的计算方法:基于紧格式动态线性化(CFDL)的MFAC控制,基于偏格式动态线性化(PFDL)的MFAC控制,基于全格式动态线性化(FFDL)的MFAC控制。
(1)基于CFDL的MFAC控制
设定液罐车动力学系统的横摆角速度变化仅与理论控制输入的当前变化量相关,其可以表示为理论控制输入变化量的线性函数,即为:
△r(k+1)=φc(k)△u(k)
式中,△r(k+1)为液罐车动力学系统在(k+1)时刻的横摆角速度变化量;Δu为理论控制输入变化量;φc为液罐车动力学系统的伪偏导数,其维度为1*1。
此时,理论控制输入量u可表示为:
式中,u(k)为k时刻的系统理论控制输入;u(k-1)为(k-1)时刻的系统理论控制输入;ρ为步长因子,计算中取值为0.5;λ为权重因子,计算中取值为1;r*(k+1)为液罐车动力学系统在(k+1)时刻的理想横摆角速度输出;r(k)为液罐车动力学系统在k时刻的实际横摆角速度输出。
液罐车动力学系统的伪偏导数估计值为:
式中,为系统在k时刻的伪偏导数估计值;为系统在k-1时刻的伪偏导数估计值;△u为系统理论控制输入的变化量;η为步长因子,计算中取值为0.8;μ为惩罚因子,计算中取值为0.5;r(k)是液罐车动力学系统在k时刻的横摆角速度输出;r(k-1)是液罐车动力学系统在(k-1)时刻的横摆角速度输出。
此时,液罐车实际的控制输入量U为:
U(k)=amp×u(k)
式中,U(k)为液罐车动力学系统实际控制输入量(附加横摆力矩);amp为扩充系数,计算中取值为2e+5。
(2)基于PFDL的MFAC控制
设定液罐车动力学系统的横摆角速度变化与滑动时间窗口[k-L+1,k]内的理论控制输入变化量相关,其可表示为:
式中,△r(k+1)为液罐车动力学系统在(k+1)时刻的横摆角速度变化量,φL(k)=[φ1(k),φ2(k),...,φL(k)]T,是液罐车动力学系统的伪偏导数矩阵,其维度为L*1;△UL(k)=[△u(k),....,△u(k-L+1)]T是滑动时间窗口内的系统理论控制输入变化量矩阵,其维度为L*1;L为控制输出线性化长度。
此时,理论控制输入量u表示为:
式中,u(k)为k时刻的系统理论控制输入;u(k-1)为(k-1)时刻的系统理论控制输入;ρ为步长因子,计算中取值为0.5,λ为权重因子,计算中取值为1;r*(k+1)为液罐车动力学系统在(k+1)时刻的理想横摆角速度输出;r(k)为液罐车动力学系统在k时刻的实际横摆角速度输出;为伪偏导数矩阵φL(k)中的第1个子矩阵;为伪偏导数矩阵φL(k)中的第i个子矩阵;△u为系统理论控制输入的变化量。
液罐车动力学系统的伪偏导数估计值为:
式中,为系统在k时刻的伪偏导数估计值;为系统在k-1时刻的伪偏导数估计值;η为步长因子,计算中取值为0.8,μ为惩罚因子,计算中取值为0.5;△UL是滑动时间窗口内的系统理论控制输入变化量矩阵;r(k)是液罐车动力学系统在k时刻的横摆角速度输出;r(k-1)是液罐车动力学系统在(k-1)时刻的横摆角速度输出。
此时,液罐车的实际控制输入量U为:
U(k)=amp×u(k)
式中,U(k)为液罐车动力学系统实际控制输入量(附加横摆力矩);amp为扩充系数,计算中取值为2e+5。
(3)基于FFDL的MFAC控制
设定液罐车动力学系统的横摆角速度变化不仅与滑动时间窗口内的理论控制输入变化量相关,还与滑动时间窗口内的输出变化量相关,其可表示为:
式中,φny,nu(k)=[φ1(k),φ2(k),...,φny(k),φny+1(k),...,φny+nu(k)]T为液罐车动力学系统的伪偏导数矩阵,其维度为(ny+nu)*1;
为滑动时间窗口[k-nu+1,k]内的控制输出变化量与滑动时间时间窗口[k-ny+1,k]内的理论控制输入变化量所组成的矩阵,其维度为(ny+nu)*1;ny为控制输出线性化长度;nu为理论控制输入线性化长度。
此时,理论控制输入量u表示为:
式中,u(k)为系统k时刻的理想控制输入;u(k-1)为系统k-1时刻的理想控制输入;ρ为步长因子,计算中取值为0.5,λ为权重因子,计算中取值为1;r*(k+1)为液罐车在k+1时刻的理想横摆角速度输出;r(k)为液罐车在k时刻的实际横摆角速度输出;△r为滑动窗口内的系统横摆角速度输出变化量;△u为滑动窗口内的系统理论控制输入变化量;为伪偏导数矩阵中的第ny+1个子矩阵;为系统伪偏导数矩阵中的第i个子矩阵。
液罐车动力学系统的伪偏导数估计值为:
式中,为系统k时刻的伪偏导数估计值;为系统k-1时刻的伪偏导数估计值;η为步长因子,计算中取值为0.8,μ惩罚因子,计算中取值为0.5;为滑动时间窗口内的系统控制输出和理论控制输入变化量矩阵;r(k)为系统在k时刻的控制输出;r(k-1)为系统在k-1时刻的控制输出。
此时,液罐车的实际控制输入量U为:
U(k)=amp×u(k)
式中,U(k)为液罐车动力学系统实际控制输入量(附加横摆力矩);amp为扩充系数,计算中取值为2e+5。
而本实施例MFAC算法为基于全格式动态线性化(FFDL)的MFAC控制算法。故基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,具体包括:
基于MFAC算法通过液罐车动力学系统的横摆角速度输出变化量、液罐车动力学系统的伪偏导数矩阵、滑动窗口内的系统输出和理论控制输入变化量矩阵之间的第三预设关系一,液罐车动力学系统的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二,液罐车动力学系统的理论控制输入、第二步长因子、液罐车动力学系统的伪偏导数矩阵、理想横摆角速度、车辆实际横摆角速度、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三,液罐车动力学系统的实际控制输入(附加横摆力矩)、液罐车动力学系统的理论控制输入之间的第三预设关系四,对行驶状态下的液罐车进行侧倾稳定性控制。
故,液罐车动力学系统(k+1)时刻的横摆角速度输出变化量、液罐车动力学系统k时刻的伪偏导数矩阵、滑动时间窗口[k-ny+1,k]内的系统横摆角速度输出变化量和滑动时间窗口[k-nu+1,k]内的理论控制输入变化量所组成的矩阵之间的第三预设关系一为:
其中,
△r(k+1)为液罐车动力学系统在(k+1)时刻的横摆角速度变化量;
是液罐车动力学系统的伪偏导数矩阵。其中,的维度为1×1,的维度为1×1,的维度为(ny+nu)×1;
是滑动窗口内的系统输出和理论控制输入变化量矩阵,其维度为(ny+nu)×1;
ny为控制输入线性化长度,nu为控制输出线性化长度。
液罐车动力学系统k时刻的伪偏导数估计值、液罐车动力学系统(k-1)时刻的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二为:
其中,
是液罐车动力学系统k时刻的伪偏导数估计值;
是液罐车动力学系统(k-1)时刻的伪偏导数估计值;
η是第一步长因子;
μ是惩罚因子;
r(k)是液罐车动力学系统在k时刻的横摆角速度输出;
r(k-1)是液罐车动力学系统在(k-1)时刻的横摆角速度输出;
是滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵。
液罐车动力学系统k时刻的理论控制输入、液罐车动力学系统(k-1)时刻的理论控制输入、第二步长因子、液罐车动力学系统k时刻的伪偏导数矩阵、液罐车动力学系统(k+1)时刻的理想横摆角速度输出、液罐车动力学系统在k时刻的实际横摆角速度输出、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三为:
其中,
u(k)为k时刻的系统理论控制输入;
u(k-1)为(k-1)时刻的系统理论控制输入;
λ为权重因子;
ρ为第二步长因子;
r*(k+1)为液罐车动力学系统在(k+1)时刻的理想横摆角速度输出;
r(k)为液罐车动力学系统在k时刻的实际横摆角速度输出;
△r为滑动窗口内的系统横摆角速度输出变化量;
△u为滑动窗口内的系统理论控制输入变化量;
为伪偏导数矩阵中的第ny+1个子矩阵;
为系统伪偏导数矩阵中的第i个子矩阵。
液罐车动力学系统k时刻的实际控制输入(附加横摆力矩)与理论控制输入之间的第三预设关系四为:
U(k)=amp×u(k),
其中:
U(k)为液罐车动力学系统在k时刻的实际控制输入量(附加横摆力矩);
amp为扩充系数。
分别使用CFDL-MAFC、PFDL-MAFC、FFDL-MAFC三种方法求解附加横摆力矩u、实施液罐车侧倾稳定性控制。其中,PFDL-MAFC中的L取值为3,FFDL-MAFC中的ny及nu取值均为3。
(1)当罐车装载圆柱罐体(罐体横截面半径为0.89m)、罐内充液比为0.6时,设定车辆行驶速度为15m/s,并通过前轮输入0.01rad的阶跃车轮转角。此时,无控制输入下的罐车发生侧倾失稳(液罐车后轴横向载荷转移率超过1.0)。三种控制方法下的液罐车动力学响应如图2-7所示。可知,CFDL-MFAC控制下车辆后轴横向载荷转移率的第一个峰值仍会超过侧翻极限值1.0,车辆依然会发生侧翻事故。PFDL-MFAC和FFDL-MAFC控制下的罐车后轴横向载荷转移率均在1.0以下,车辆侧倾稳定性得到控制。与PFDL-MFAC相比,FFDL-MAFC的控制效果更佳。
(2)更改罐体形状为椭圆柱罐体(罐体长轴半径1.09m,短轴半径0.73m,长短轴之比为1.5;罐体长轴半径1.26m,短轴半径0.63m,长短轴之比为2),车辆运行条件不变。无控制输入下的液罐车发生侧倾失稳。此时,三种控制方法下的液罐车动力学响应分别如图8-19所示。同样地,三种控制策略中,FFDL-MFAC的控制效果最佳。
(3)进一步探究滑动窗口大小对FFDL-MFAC控制方法的影响,在ny=nu=3的基础上,对比ny=nu=4以及ny=nu=5时的液罐车控制效果。装载圆柱罐体液罐车的控制效果如图20-25所示(充液比0.6,行驶速度15m/s,前轮角阶跃输入0.01rad)。可知,ny=nu=3及ny=nu=4时,FFDL-MFAC对液罐车的控制效果基本相同。当ny=nu=5时,车辆横摆角速度的超调量得到了显著下降,车辆后轴横向载荷转移率的第一个峰值也得到了很好的控制,车辆侧倾稳定性要高于另外两组参数下的控制效果。
(4)罐体装载椭圆柱罐体(罐体长轴半径1.09m,短轴半径0.73m,长短轴之比为1.5;罐体长轴半径1.26m,短轴半径0.63m,长短轴之比为2)、车辆运行条件不变时,三种FFDL-MFAC控制方法对液罐车侧倾稳定性的控制效果如图26-37所示。同样地,ny=nu=5时FFDL-MFAC对罐车侧倾稳定性的控制效果最佳。
通过上述分析,针对液罐车难以建模、模型通用性差、模型准确度低的问题,MFAC控制方法能略去建模过程,而仅利用系统的控制输入和输出进行车辆侧倾控制,且FFDL-MFAC控制方法能够取得很好的控制效果。这对于保障液罐车的道路运输安全具有十分重要的意义。
本实施例还提供了一种基于MFAC算法的液罐车侧倾稳定性控制系统,包括:
液罐车侧倾稳定性影响因素确定单元,配置用于基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;
液罐车侧倾稳定性控制单元,配置用于基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
所述MFAC算法包括基于全格式动态线性化的MFAC算法。
应当理解,上述基于MFAC算法的液罐车侧倾控制系统记载的诸子系统或单元与上述基于MFAC算法的液罐车侧倾控制方法描述的各个步骤相对应。由此,上文针对方法描述的操作和特征同样适用于上述基于MFAC算法的液罐车侧倾控制系统的诸子系统及其中包含的单元,在此不再赘述。
以上描述仅为本申请的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本申请中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离所述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能。

Claims (10)

1.一种基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,包括以下步骤:
基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;
基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
2.根据权利要求1所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,所述MFAC算法包括基于全格式动态线性化的MFAC算法。
3.根据权利要求1所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系为:
其中,
αy为液罐车的车身侧向加速度;
V为液罐车行驶速度;
r为液罐车横摆角速度;
为液罐车质心侧偏角速度。
4.根据权利要求1所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,包括:
基于MFAC算法通过液罐车动力学系统的横摆角速度输出变化量、液罐车动力学系统的伪偏导数矩阵、滑动窗口内的系统输出和理论控制输入变化量矩阵之间的第三预设关系一,液罐车动力学系统的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二,液罐车动力学系统的理论控制输入、第二步长因子、液罐车动力学系统的伪偏导数矩阵、理想横摆角速度、车辆实际横摆角速度、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三,液罐车动力学系统的实际控制输入(附加横摆力矩)、液罐车动力学系统的理论控制输入之间的第三预设关系四,对行驶状态下的液罐车进行侧倾稳定性控制。
5.根据权利要求4所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,液罐车动力学系统(k+1)时刻的横摆角速度输出变化量、液罐车动力学系统k时刻的伪偏导数矩阵、滑动时间窗口[k-ny+1,k]内的系统横摆角速度输出变化量和滑动时间窗口[k-nu+1,k]内的理论控制输入变化量所组成的矩阵之间的第三预设关系一为:
其中,
△r(k+1)为液罐车动力学系统在(k+1)时刻的横摆角速度变化量;
是液罐车动力学系统的伪偏导数矩阵,其中,的维度为1×1,的维度为1×1,的维度为(ny+nu)×1;
是滑动窗口内的系统输出和理论控制输入变化量矩阵,其维度为(ny+nu)×1;
ny为控制输入线性化长度,nu为控制输出线性化长度。
6.根据权利要求4所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,液罐车动力学系统k时刻的伪偏导数估计值、液罐车动力学系统(k-1)时刻的伪偏导数估计值、第一步长因子、滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵、惩罚因子之间的第三预设关系二为:
其中,
是液罐车动力学系统k时刻的伪偏导数估计值;
是液罐车动力学系统(k-1)时刻的伪偏导数估计值;
η是第一步长因子;
μ是惩罚因子;
r(k)是液罐车动力学系统在k时刻的横摆角速度输出;
r(k-1)是液罐车动力学系统在(k-1)时刻的横摆角速度输出;
是滑动时间窗口内的系统横摆角速度输出和理论控制输入变化量矩阵。
7.根据权利要求4所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,液罐车动力学系统k时刻的理论控制输入、液罐车动力学系统(k-1)时刻的理论控制输入、第二步长因子、液罐车动力学系统k时刻的伪偏导数矩阵、液罐车动力学系统(k+1)时刻的理想横摆角速度输出、液罐车动力学系统在k时刻的实际横摆角速度输出、权重因子、系统理论控制输入变化量、系统控制输出变化量之间的第三预设关系三为:
其中,
u(k)为k时刻的系统理论控制输入;
u(k-1)为(k-1)时刻的系统理论控制输入;
λ为权重因子;
ρ为第二步长因子;
r*(k+1)为液罐车动力学系统在(k+1)时刻的理想横摆角速度输出;
r(k)为液罐车动力学系统在k时刻的实际横摆角速度输出;
△r为滑动窗口内的系统横摆角速度输出变化量;
△u为滑动窗口内的系统理论控制输入变化量;
为伪偏导数矩阵中的第ny+1个子矩阵;
为系统伪偏导数矩阵中的第i个子矩阵。
8.根据权利要求4所述的基于MFAC算法的液罐车侧倾稳定性控制方法,其特征是,液罐车动力学系统k时刻的实际控制输入(附加横摆力矩)与理论控制输入之间的第三预设关系四为:
U(k)=amp×u(k),
其中:
U(k)为液罐车动力学系统在k时刻的实际控制输入量(附加横摆力矩);
amp为扩充系数。
9.一种基于MFAC算法的液罐车侧倾稳定性控制系统,其特征是,包括:
液罐车侧倾稳定性影响因素确定单元,配置用于基于液罐车的车身侧向加速度、质心侧偏角速度、横摆角速度之间的第一预设关系及液罐车的车身侧向加速度与侧倾稳定性之间的第二预设关系,确定所述液罐车侧倾稳定性的影响因素包括液罐车的横摆角速度;
液罐车侧倾稳定性控制单元,配置用于基于MFAC算法通过附加横摆力矩与液罐车横摆角速度之间的第三预设关系对行驶状态下的液罐车进行侧倾稳定性控制,其中,附加横摆力矩为控制输入量,罐车横摆角速度为被控量。
10.根据权利要求7所述的基于MFAC算法的液罐车侧倾稳定性控制系统,其特征是,所述MFAC算法包括基于全格式动态线性化的MFAC算法。
CN201910756801.3A 2019-08-16 2019-08-16 基于mfac算法的液罐车侧倾稳定性控制方法、系统 Active CN110395263B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910756801.3A CN110395263B (zh) 2019-08-16 2019-08-16 基于mfac算法的液罐车侧倾稳定性控制方法、系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910756801.3A CN110395263B (zh) 2019-08-16 2019-08-16 基于mfac算法的液罐车侧倾稳定性控制方法、系统

Publications (2)

Publication Number Publication Date
CN110395263A true CN110395263A (zh) 2019-11-01
CN110395263B CN110395263B (zh) 2020-08-11

Family

ID=68328460

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910756801.3A Active CN110395263B (zh) 2019-08-16 2019-08-16 基于mfac算法的液罐车侧倾稳定性控制方法、系统

Country Status (1)

Country Link
CN (1) CN110395263B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114454199A (zh) * 2022-02-11 2022-05-10 吉林大学 一种基于数据控制算法的末端执行器及其设计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050015191A1 (en) * 2003-07-17 2005-01-20 Masaki Banno Anti-rollover device for vehicle
KR20070076675A (ko) * 2006-01-19 2007-07-25 주식회사 만도 차량 안정성 제어 시스템의 횡가속도 추정 방법
CN105946863A (zh) * 2016-06-23 2016-09-21 吉林大学 一种车辆行驶稳定性区域的确定方法
CN108657175A (zh) * 2018-05-04 2018-10-16 合肥工业大学 一种基于gis系统的液罐车侧翻预警系统
CN109101697A (zh) * 2018-07-18 2018-12-28 吉林大学 液罐汽车流-固耦合系统的整车动力学模型的建模方法
US20190217854A1 (en) * 2018-01-12 2019-07-18 Hyundai Motor Company Vehicle stability control system and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050015191A1 (en) * 2003-07-17 2005-01-20 Masaki Banno Anti-rollover device for vehicle
KR20070076675A (ko) * 2006-01-19 2007-07-25 주식회사 만도 차량 안정성 제어 시스템의 횡가속도 추정 방법
CN105946863A (zh) * 2016-06-23 2016-09-21 吉林大学 一种车辆行驶稳定性区域的确定方法
US20190217854A1 (en) * 2018-01-12 2019-07-18 Hyundai Motor Company Vehicle stability control system and method
CN108657175A (zh) * 2018-05-04 2018-10-16 合肥工业大学 一种基于gis系统的液罐车侧翻预警系统
CN109101697A (zh) * 2018-07-18 2018-12-28 吉林大学 液罐汽车流-固耦合系统的整车动力学模型的建模方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑雪莲: "基于液体冲击等效机械模型的汽车罐车行驶稳定性研究", 《工程科技Ⅱ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114454199A (zh) * 2022-02-11 2022-05-10 吉林大学 一种基于数据控制算法的末端执行器及其设计方法
CN114454199B (zh) * 2022-02-11 2024-02-23 吉林大学 一种基于数据控制算法的末端执行器及其设计方法

Also Published As

Publication number Publication date
CN110395263B (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
US11364895B2 (en) Yaw motion control method for four-wheel distributed vehicle
CN111267835B (zh) 基于模型预测算法的四轮独立驱动汽车稳定性控制方法
US10967870B2 (en) Hill descent system for vehicle and control method thereof
Chen et al. Electronic stability control for electric vehicle with four in-wheel motors
CN103303367B (zh) 一种四轮驱动电动车车身稳定控制方法
US20090108671A1 (en) Vehicle drive apparatus
WO2012023162A1 (en) A vehicle
CN108237950A (zh) 车辆的控制方法、系统及车辆
US11318804B2 (en) Vehicle state estimation device, control device, suspension control device, and suspension device
CN104925054A (zh) 一种基于微分平坦的车辆稳定转向集成控制方法
CN110422053A (zh) 四轮轮毂电机驱动电动汽车节能控制方法
Ni et al. Relaxed static stability based on tyre cornering stiffness estimation for all-wheel-drive electric vehicle
CN110395263A (zh) 基于mfac算法的液罐车侧倾稳定性控制方法、系统
Peng et al. Integrated vehicle control based on tire force reserve optimization concept
Hu et al. MTTE-based motion stabilization control for in-wheel motor electric vehicles
Oke et al. H∞ dynamic output feedback control for independently driven four-wheel electric vehicles with differential speed steering
CN115703363A (zh) 用于确定道路车辆的驱动车轮的优化的扭矩分配的方法和相关的道路车辆
Jacobsen Potential of electric wheel motors as new chassis actuators for vehicle manoeuvring
Zhao et al. Traction control system of electric vehicle with 4 in-wheel motors using lyapunov stability analysis algorithm
Ueno et al. Direct yaw moment control for electric vehicles with variable-rate-slip-ratio-limiter based driving force control
Mutoh et al. Cornering performance of front and rear wheel independent drive type electric vehicles (FRID EVs) under severe road conditions
Tao et al. Research on Control Strategy of the Side Forklift Steering Movement Stability and Its Simulation Analysis
Mihály et al. Integrated vehicle control of in-wheel electric vehicle
Jing et al. Effect of the tire lateral force saturation on stability control of 4WDEV steering
CN112758081B (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