CN112389450A - 一种基于板簧形变的罐车整车质量动态估计方法 - Google Patents

一种基于板簧形变的罐车整车质量动态估计方法 Download PDF

Info

Publication number
CN112389450A
CN112389450A CN202011310790.5A CN202011310790A CN112389450A CN 112389450 A CN112389450 A CN 112389450A CN 202011310790 A CN202011310790 A CN 202011310790A CN 112389450 A CN112389450 A CN 112389450A
Authority
CN
China
Prior art keywords
tank car
plate spring
resistance
mass
resistance strain
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.)
Pending
Application number
CN202011310790.5A
Other languages
English (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.)
Research Institute of Highway Ministry of Transport
Original Assignee
Research Institute of Highway Ministry of Transport
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 Research Institute of Highway Ministry of Transport filed Critical Research Institute of Highway Ministry of Transport
Priority to CN202011310790.5A priority Critical patent/CN112389450A/zh
Publication of CN112389450A publication Critical patent/CN112389450A/zh
Pending legal-status Critical Current

Links

Images

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
    • 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/12Estimation 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 parameters of the vehicle itself, e.g. tyre models
    • B60W40/13Load or weight

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明提出了一种基于板簧形变的罐车整车质量动态估计方法,首先明确用于估计整车质量的参数及其获取方式,然后设计罐车典型行驶场景并开展实车试验获得样本信息,然后利用LS‑SVM拟合电阻应变片电阻与整车质量的函数表达式,最后基于拟合出的函数表达式实现罐车整车质量的动态估计。该方法利用电阻感应板簧形变来估计罐车整车质量,操作方便,便于工程使用。同时通过两侧板簧形变程度的不同体现出罐车质心位置的影响,减小晃动对整车质量估计产生的误差。

Description

一种基于板簧形变的罐车整车质量动态估计方法
技术领域
本发明涉及一种罐车质量动态估计方法,具体的是一种基于板簧形变的罐车整车质量动态估计方法,属于车载估计技术领域。
背景技术
近年来,由于装载量大、运输效率高,罐车的使用越来越多,人们对于罐车的需求不再简单局限于运输工具的范畴,对其运行效率、行驶稳定性和主动安全性等的追求越来越高,故驾驶辅助电子系统广泛运用于罐车上,如制动防抱死系统(ABS)、车身电子稳定系统(ESP)等。这些电子系统大多是根据罐车行驶状态参数进行工作的,其中整车质量是一个关键参数,然而罐车容积大,整车质量随着装载货物种类和体积的不同而经常变化,因此准确实时估计整车质量对上述电子系统自适应调整其控制算法具有重要意义。同时罐车装载状态异常是导致车辆发生事故主要原因之一,我国每年由于此类原因引发的交通事故占罐车交通事故的60%以上,而对罐车整车质量的动态估计是避免此类事故的有效方法。
现有的整车质量估计方法主要是基于车辆纵向动力学模型实现整车质量的动态估计,这些方法将车辆的速度、加速度和发动机扭矩等信息带入模型从而估算出整车质量,如青岛大学臧政等人利用车辆发动机扭矩作为动力学模型输入估算整车质量(臧政,霍炜,王玉海,李兴坤,郑旭光,李圆圆.重型商用车辆质量估计算法研究[J].中国机械工程)。然而建立的动力学模型复杂,工程化困难,并且大多数方法是通过多维仿真试验来验证其算法的有效性,未对算法进行实车道路测试;并且,在估计罐车整车质量时,由于受液体扰动的影响,罐车质心位置运动多变,给模型估算带来一定的误差。
发明内容
基于现状,针对使用钢板弹簧悬架的罐车,本发明提出了一种基于板簧形变的罐车整车质量动态估计方法,该方法首先明确用于估计整车质量的参数及其获取方式,然后设计罐车典型行驶场景并开展实车试验获得样本信息,接着利用LS-SVM拟合电阻应变片电阻与整车质量的函数表达式,最后基于拟合出的函数表达式实现罐车整车质量的动态估计。该方法利用板簧受力形变来估计罐车整车质量,操作方便,便于工程使用;通过两侧板簧形变程度不同体现出的液体晃动对罐车质心位置的影响,就此弥补液体晃动对整车质量估计产生的误差;开展罐车实车试验获得样本数据,建立板簧形变与整车质量之间的联系,从而实现整车质量的动态估计。
本发明方案的具体步骤如下:一种基于板簧形变的罐车整车质量动态估计方法,其特征在于:通过在罐车每个钢板弹簧悬架最上端的钢板表面固定安装电阻应变片,根据预先拟合建立的电阻应变片电阻与整车质量的函数关系,应用电阻应变片感应钢板弹簧受力形变所输出的相应电阻值体现整车质量变化。
进一步地,利用LS-SVM模型拟合电阻应变片电阻与整车质量的函数关系。
具体拟合方法为:
采用“驾驶行为”、“载荷”和“纵向坡度”作为行驶场景基本元素,几种元素遍历组合,形成一定数量的车辆行驶试验场景,在每种试验场景下经过一定次数的试验,获得电阻应变片电阻和整车质量的关联样本数据,然后基于LS-SVM模型拟合成函数关系。
所述“驾驶行为”、“载荷”和“纵向坡度”量化如下:
Figure BDA0002789688700000021
总计有3×5×2=30种车辆行驶试验场景。
进一步地,设罐车左、右两侧车轴上方钢板弹簧悬架上均设置电阻应变片,设左侧电阻应变片电阻总和为
Figure BDA0002789688700000022
右侧电阻应变片电阻总和为
Figure BDA0002789688700000023
其中s表示罐车车轴编号,n表示罐车车轴总数,
Figure BDA0002789688700000024
表示罐车第s根车轴左侧钢板弹簧悬架上电阻应变片的电阻值,
Figure BDA0002789688700000025
表示罐车第s根车轴右侧钢板弹簧悬架上电阻应变片的电阻值;
应用LS-SVM模型,对于N组试验数据,建立训练集(xi,yi),其中
Figure BDA0002789688700000031
yi=Mi,i=1,2,…,N;
进行回归输出,函数表达式表示为:
Figure BDA0002789688700000032
式中,f(x)即整车质量,μi和b由式(7)、(8)、(9)求得,
Figure BDA0002789688700000033
式(7)中,1N为N×N维的全1向量,IN为N×N维的单位矩阵,μN=[μ12,…,μN]T,Y=[y1,y2,…,yN]T,Ω为N×N维的核矩阵,其定义为:
Ωij=φ(xi)Tφ(xj)=K(xi,xj),i,j=1,2,…,N (8)
式(8)中,K(xi,xj)为核函数,使用高斯核函数,表达式为:
Figure BDA0002789688700000034
Figure BDA0002789688700000035
是样本数据,预先存储,基于式(10),根据每次测得的x,x=[Rl,Rr]T就可输出相应的f(x)。
与现有技术相比,本发明的优点体现在:
1.本发明利用电阻应变片感应板簧形变来估计罐车整车质量,操作方便,便于工程使用;
2.本发明通过在车轴两侧安装电阻应变片,通过感应两侧板簧形变程度的不同体现出液体晃动对罐车质心位置的影响,减小了液体晃动对整车质量估计产生的误差;
3.本发明通过实车试验获得关于电阻应变片电阻值与质量的样本数据,基于LS-SVM模型拟合出电阻与整车质量的函数表达式,实现整车质量的动态估计。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制。
图1为罐车整车质量动态估计方法的总体实施流程图;
图2为电阻应变片在车辆上的安装位置示意图;
图3为用训练好的样本,应用LS-SVM模型进行实际整车质量估计的过程图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述,但本领域的技术人员应该知道,以下实施例并不是对本发明技术方案作的唯一限定,凡是在本发明技术方案精神实质下所做的任何等同变换或改动,均应视为属于本发明的保护范围。
本发明罐车整车质量动态估计方法的总体实施流程,如图1所示,具体步骤包括如下:
步骤一:确定用于整车质量估计的参数及获取方式
罐车整车质量分布作用在各个车轮上,通过车轮所受垂直载荷可以间接估算整车质量,然而高速运转状态下车轮的载荷测量比较困难。针对使用钢板弹簧悬架的罐车,当钢板弹簧所承受的垂直载荷为正向时,各弹簧片都受力变形,有向上拱弯的趋势,车轮处钢板弹簧形变与车轮载荷存在对应关系,因此可以通过钢板弹簧形变估计整车质量。
电阻应变片是用于测量物体应变的元件,它能将机械构件上应变的变化转换为电阻变化,通过电信号的方式呈现出来。因此,如图2所示,将电阻应变片1固定在罐车每个钢板弹簧悬架2最上端的钢板表面,当钢板弹簧受力发生形变时,电阻应变片1就能感应形变,通过输出相应的电阻值体现该变化,也就是体现车辆受载情况。
通过电阻测量仪测得罐车左侧电阻应变片电阻和为
Figure BDA0002789688700000041
右侧电阻应变片电阻和为
Figure BDA0002789688700000042
其中s表示罐车车轴编号,n表示罐车车轴总数,
Figure BDA0002789688700000043
表示罐车第s根车轴左侧钢板弹簧悬架上电阻应变片的电阻值,
Figure BDA0002789688700000044
表示罐车第s根车轴右侧钢板弹簧悬架上电阻应变片的电阻值。
步骤二:设计罐车典型行驶场景并开展实车试验获得样本信息
设计罐车行驶场景并开展实车试验,获得在不同整车质量下的电阻应变片电阻值,通过拟合建立电阻应变片电阻与整车质量的函数关系。
设罐车行驶场景基本元素包括“驾驶行为”、“载荷”和“纵向坡度”,各基元划分如下表所示。
Figure BDA0002789688700000045
Figure BDA0002789688700000051
J转向试验轨迹设置参考GB/T6323-2014《汽车操纵稳定性试验方法》和JT/T1094-2016《营运客车安全技术条件》中转向试验规定,双移线试验轨迹设置参考ISO 3888-2《乘用车--急剧变换车道操纵用试验车道--第2部分:障碍物规避》中规定要求。
场景基元排列组合后有3(驾驶行为)×5(载荷)×2(纵向坡度)=30种测试场景。设罐车的整车质量为M,则计算公式为:
Figure BDA0002789688700000052
式(1)中,m为罐车空载时整车质量,V为罐车容积,ρ为装载物的密度,m和V由罐车产品手册获得。
设置测试场景后开展试验,在保证安全的前提下,驾驶员可自由改变车速并换挡,30种场景试验结束后,保存罐车左侧电阻应变片电阻和、右侧电阻应变片电阻和、整车质量的信息,得到
Figure BDA0002789688700000053
组样本数据,tα为第α种场景试验中关于罐车电阻应变片电阻和整车质量的数据量,数据集如式(2)所示,在单一场景试验中,整车质量为定值,所以有M1=M2=...M
Figure BDA0002789688700000054
步骤三:利用LS-SVM拟合电阻应变片电阻与整车质量的函数表达式
获得电阻应变片电阻与整车质量的样本数据后,采用泛化性能好的最小二乘支持向量机(LS-SVM)建立非线性回归模型。
对于N组数据,建立训练集(xi,yi),
Figure BDA0002789688700000061
yi=Mi,i=1,2,…,N,上角标T表示对矩阵转置,LS-SVM模型在特征空间内的表达形式为:
f(x)=wTφ(x)+b (3)
式(3)中,f(x)即yi,w是可调节的权重向量;φ(·)表示非线性映射,把输入量映射到更高维的特征空间;b为标量阈值。LS-SVM函数表达式求解的目标函数为:
Figure BDA0002789688700000062
式(4)中,ei是误差变量,γ是正则化常数,可以在训练过程中调整优化,较小的γ可以在数据噪声较大时避免过拟合现象的发生。利用拉格朗日乘数法将式(4)改写成无约束形式可得:
Figure BDA0002789688700000063
进一步推导:
Figure BDA0002789688700000064
由式(6)可以得到一个线性系统,从而避免了二次规划问题,如式(7)所示:
Figure BDA0002789688700000065
式(7)中,1N为N×N维的全1向量,IN为N×N维的单位矩阵,μN=[μ12,…,μN]T,Y=[y1,y2,…,yN]T,Ω为N×N维的核矩阵,其定义为:
Ωij=φ(xi)Tφ(xj)=K(xi,xj),i,j=1,2,…,N (8)
式(8)中,K(xi,xj)为核函数,使用高斯核函数,表达式为:
Figure BDA0002789688700000071
式(9)中,σ为高斯核的带宽,可以在训练过程中进行调整优化。因此,利用LS-SVM模型进行回归输出的函数表达式可以写成:
Figure BDA0002789688700000072
式(10)中,μi和b由式(7)、(8)、(9)求得,
Figure BDA0002789688700000073
是样本数据,预先已存储。这样,关于LS-SVM模型的表达式就由式(3)转变为式(10)表达。
步骤四:基于拟合出的函数表达式实现罐车整车质量的动态估计
如图3所示,罐车装载货物行驶时,实时读取罐车左侧电阻应变片的电阻和Rl、右侧电阻应变片的电阻和Rr,从而有x=[Rl,Rr]T,应用样本数据,加载到LS-SVM模型中,就可以估计出罐车整车质量f(x)。

Claims (5)

1.一种基于板簧形变的罐车整车质量动态估计方法,其特征在于:
通过在罐车每个钢板弹簧悬架最上端的钢板表面固定安装电阻应变片,根据预先拟合建立的电阻应变片电阻与整车质量的函数关系,应用电阻应变片感应钢板弹簧受力形变所输出的相应电阻值体现整车质量变化。
2.根据权利要求1所述基于板簧形变的罐车整车质量动态估计方法,其特征在于:利用LS-SVM模型拟合电阻应变片电阻与整车质量的函数关系。
3.根据权利要求2所述基于板簧形变的罐车整车质量动态估计方法,其特征在于:具体拟合方法为:
采用“驾驶行为”、“载荷”和“纵向坡度”作为行驶场景基本元素,几种元素遍历组合,形成一定数量的车辆行驶试验场景,在每种试验场景下经过一定次数的试验,获得电阻应变片电阻和整车质量的关联样本数据,然后基于LS-SVM模型拟合成函数关系。
4.根据权利要求3所述基于板簧形变的罐车整车质量动态估计方法,其特征在于:所述“驾驶行为”、“载荷”和“纵向坡度”量化如下:
Figure FDA0002789688690000011
总计有3×5×2=30种车辆行驶试验场景。
5.根据权利要求1或2或3或4所述基于板簧形变的罐车整车质量动态估计方法,其特征在于:设罐车左、右两侧车轴上方钢板弹簧悬架上均设置电阻应变片,设左侧电阻应变片电阻总和为
Figure FDA0002789688690000012
右侧电阻应变片电阻总和为
Figure FDA0002789688690000013
其中s表示罐车车轴编号,n表示罐车车轴总数,
Figure FDA0002789688690000014
表示罐车第s根车轴左侧钢板弹簧悬架上电阻应变片的电阻值,
Figure FDA0002789688690000015
表示罐车第s根车轴右侧钢板弹簧悬架上电阻应变片的电阻值;
应用LS-SVM模型,对于N组试验数据,建立训练集(xi,yi),其中
Figure FDA0002789688690000021
yi=Mi,i=1,2,…,N;
进行回归输出,函数表达式表示为:
Figure FDA0002789688690000022
式中,f(x)即整车质量,μi和b由式(7)、(8)、(9)求得,
Figure FDA0002789688690000023
式(7)中,1N为N×N维的全1向量,IN为N×N维的单位矩阵,μN=[μ12,…,μN]T,Y=[y1,y2,…,yN]T,Ω为N×N维的核矩阵,其定义为:
Ωij=φ(xi)Tφ(xj)=K(xi,xj),i,j=1,2,…,N (8)
式(8)中,K(xi,xj)为核函数,使用高斯核函数,表达式为:
Figure FDA0002789688690000024
Figure FDA0002789688690000025
是样本数据,预先存储,基于式(10),根据每次测得的x,x=[Rl,Rr]T就可输出相应的f(x)。
CN202011310790.5A 2020-11-20 2020-11-20 一种基于板簧形变的罐车整车质量动态估计方法 Pending CN112389450A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011310790.5A CN112389450A (zh) 2020-11-20 2020-11-20 一种基于板簧形变的罐车整车质量动态估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011310790.5A CN112389450A (zh) 2020-11-20 2020-11-20 一种基于板簧形变的罐车整车质量动态估计方法

Publications (1)

Publication Number Publication Date
CN112389450A true CN112389450A (zh) 2021-02-23

Family

ID=74605964

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011310790.5A Pending CN112389450A (zh) 2020-11-20 2020-11-20 一种基于板簧形变的罐车整车质量动态估计方法

Country Status (1)

Country Link
CN (1) CN112389450A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114877980A (zh) * 2022-04-22 2022-08-09 武汉工程大学 一种混凝土罐车实时称重系统及方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202454077U (zh) * 2012-01-04 2012-09-26 华北电力大学 一种汽车超载实时监控和报警装置
CN104568096A (zh) * 2014-12-23 2015-04-29 潍柴动力股份有限公司 一种具有板簧悬架的车辆及其重量测量系统、方法
CN105910692A (zh) * 2016-06-20 2016-08-31 河北科技大学 基于加速度传感器的车辆载重测量方法
CN207379598U (zh) * 2017-10-17 2018-05-18 山东交通学院 一种动静两用货车称重系统
WO2018214229A1 (zh) * 2017-05-24 2018-11-29 深圳市靖洲科技有限公司 一种三孔悬臂梁应变片的动态电子小型无人车衡
CN109916487A (zh) * 2017-12-13 2019-06-21 北京万集科技股份有限公司 行车重量智能监控系统及方法
CN111639715A (zh) * 2020-06-01 2020-09-08 重庆大学 基于ls-svm的汽车仪表组装质量预测方法与系统
CN111695197A (zh) * 2020-05-23 2020-09-22 东南大学 一种罐车侧翻阈值高可靠动态估计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202454077U (zh) * 2012-01-04 2012-09-26 华北电力大学 一种汽车超载实时监控和报警装置
CN104568096A (zh) * 2014-12-23 2015-04-29 潍柴动力股份有限公司 一种具有板簧悬架的车辆及其重量测量系统、方法
CN105910692A (zh) * 2016-06-20 2016-08-31 河北科技大学 基于加速度传感器的车辆载重测量方法
WO2018214229A1 (zh) * 2017-05-24 2018-11-29 深圳市靖洲科技有限公司 一种三孔悬臂梁应变片的动态电子小型无人车衡
CN207379598U (zh) * 2017-10-17 2018-05-18 山东交通学院 一种动静两用货车称重系统
CN109916487A (zh) * 2017-12-13 2019-06-21 北京万集科技股份有限公司 行车重量智能监控系统及方法
CN111695197A (zh) * 2020-05-23 2020-09-22 东南大学 一种罐车侧翻阈值高可靠动态估计方法
CN111639715A (zh) * 2020-06-01 2020-09-08 重庆大学 基于ls-svm的汽车仪表组装质量预测方法与系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
魏洋: "基于最小二乘支持向量机的称重系统偏载优化方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114877980A (zh) * 2022-04-22 2022-08-09 武汉工程大学 一种混凝土罐车实时称重系统及方法

Similar Documents

Publication Publication Date Title
Cheng et al. Parameter and state estimation for articulated heavy vehicles
Doumiati et al. Lateral load transfer and normal forces estimation for vehicle safety: experimental test
Doumiati et al. Observers for vehicle tyre/road forces estimation: experimental validation
US11008014B2 (en) Methods and apparatus to determine vehicle weight information based on ride height
Jin et al. Stability and optimised H∞ control of tripped and untripped vehicle rollover
CN103476657A (zh) 用于路面摩擦系数的实时估计的方法和装置
CN113359457B (zh) 智能车底盘域控制器的高维动力学模型解算装置及方法
Beal Rapid road friction estimation using independent left/right steering torque measurements
CN108556850B (zh) 一种多轴重型车辆车轮垂向载荷通用计算方法
CN108819950A (zh) 汽车稳定性控制系统的车速估计方法及系统
Mashadi et al. Vehicle lift-off modelling and a new rollover detection criterion
Zhang et al. Investigation into untripped rollover of light vehicles in the modified fishhook and the sine maneuvers. Part I: Vehicle modelling, roll and yaw instability
WO2022134929A1 (zh) 一种确定车辆质量方法、装置、设备及介质
Sharma et al. Vehicle dynamic state observation and rolling resistance estimation via unknown input adaptive high gain observer
CN112389450A (zh) 一种基于板簧形变的罐车整车质量动态估计方法
Widner et al. Framework for vehicle dynamics model validation
US10836386B2 (en) Determination of roll angle and bank angle with suspension displacement data
Abdulwahab et al. Estimation of LTR rollover index for a high-sided tractor semitrailer vehicle under extreme crosswind conditions through dynamic simulation
CN109131344B (zh) 用于实时确定轮胎法向力的系统和方法
CN115406669A (zh) 一种多轴特种车辆侧翻指标优化方法
CN111695197B (zh) 一种罐车侧翻阈值高可靠动态估计方法
CN112526880B (zh) 一种车辆行驶过程中的路面高度实时估计方法
Marotta et al. Estimation of the Tire-Road Interaction Forces by using Pacejka’s Formulas with Combined Slips and Camber Angles
Van et al. Dynamic load determination of a multi-purpose forest fire fighting vehicle when operating on forestry road
Marimuthu et al. A study on SUV parameters sensitivity on rollover propensity

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210223