CN104802826B - 用于估计轨道车辆的横向力的装置 - Google Patents
用于估计轨道车辆的横向力的装置 Download PDFInfo
- Publication number
- CN104802826B CN104802826B CN201510088778.7A CN201510088778A CN104802826B CN 104802826 B CN104802826 B CN 104802826B CN 201510088778 A CN201510088778 A CN 201510088778A CN 104802826 B CN104802826 B CN 104802826B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- estimate
- msub
- rail vehicle
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000001133 acceleration Effects 0.000 claims abstract description 20
- 238000005259 measurement Methods 0.000 claims description 14
- 238000011156 evaluation Methods 0.000 claims description 4
- 238000005452 bending Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/042—Track changes detection
- B61L23/047—Track or rail movements
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61F—RAIL VEHICLE SUSPENSIONS, e.g. UNDERFRAMES, BOGIES OR ARRANGEMENTS OF WHEEL AXLES; RAIL VEHICLES FOR USE ON TRACKS OF DIFFERENT WIDTH; PREVENTING DERAILING OF RAIL VEHICLES; WHEEL GUARDS, OBSTRUCTION REMOVERS OR THE LIKE FOR RAIL VEHICLES
- B61F9/00—Rail vehicles characterised by means for preventing derailing, e.g. by use of guide wheels
- B61F9/005—Rail vehicles characterised by means for preventing derailing, e.g. by use of guide wheels by use of non-mechanical means, e.g. acoustic or electromagnetic devices
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L15/00—Indicators provided on the vehicle or train for signalling purposes
- B61L15/0081—On-board diagnosis or maintenance
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Electromagnetism (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本公开涉及一种用于估计当轨道车辆行驶于弯曲路段时由于车轮与轨道之间的接触而导致的施加到转向架的横向力的装置及方法,该装置包括:横向速度估计观测器,其配置为通过基于轨道车辆的纵向加速度、横向加速度、横摆角速度、以及车轮角速度估计横向速度,来计算横向速度估计值;以及横向力估计观测器,其配置为通过基于轨道车辆的转向角、施加到轨道车辆的纵向力、以及由横向速度估计观测器计算出的横向速度估计值估计施加到轨道车辆的转向架的横向力,来计算横向力估计值。
Description
技术领域
本公开涉及一种用于估计轨道车辆的横向力的装置及方法。更具体地,本公开涉及一种用于估计当轨道车辆行驶于弯曲路段时由于车轮与轨道之间的接触而导致的施加到转向架的横向力的装置及方法。
背景技术
关于施加到轨道车辆的转向架的横向力的信息是确定火车脱轨可能性的因素。由于该原因,横向力是在行驶在弯曲路段中时表示火车运动的一个关键因素。
此外,关于横向力的信息被用作轨道车辆主动转向控制的关键控制因素。
用于测量在弯曲路段中时的横向力的相关技术公开在韩国专利公开号No.10-2013-0055110(“Tire lateral force estimation method and device(轮胎横向力估计方法和设备)”,下文中称作“参考文献1”)和美国专利号No.7,853,412(“Estimation ofwheel rail interaction forces(车轮轨道相互作用力的估计)”,下文中称作“参考文献2”)中。
参考文献1公开了一种用于检测施加到机动车辆的轮胎的横向力的设备。它涉及一种用于检测施加到轮胎的横向力的方法,通过进行配置有多个传感器的车辆的实际驾驶测试,收集关于车辆运动的数据,并且应用该数据到参考车辆模型和卡尔曼(Kalman)估计以计算轮胎模型的参数。
参考文献2公开了一种用于检测在轨道车辆的车轮与轨道之间的横向力和纵向力的设备。它涉及一种用于检测横向力的方法,通过将轨道车辆构造为十三个自由度的动态模型,并且使用从安装在车辆中的加速度传感器获得的信息和由于在轨道与车轮之间的接触而产生的横向力和纵向力的模型计算横向力和纵向力。
参考文献1公开了一种用于检测施加到机动车辆的轮胎的横向力的方法。然而,该方法难以直接应用到轨道车辆,并且具有需要复杂轮胎模型的缺点。
此外,使用轮胎模型检测横向力的技术需要该轮胎模型的精度。因此,估计的值依赖于轮胎模型的精度。
此外,参考文献2公开了一种检测轨道车辆的横向力和纵向力的方法。然而,该方法是基于与横向力相关的数学模型的。因此,该方法的缺点是,估计的横向力依赖于该数学模型的精度。
发明内容
为了克服常规技术的问题,本公开提供一种用于估计施加于轨道车辆的前和后转向架的横向力的装置和方法,通过使用轨道车辆主体的动态模型和来自传感器的数据测量值,而不需要任何复杂的横向力数学模型。
在本公开的总体方面,提供了一种用于估计轨道车辆的横向力的装置,该装置包括:横向速度估计观测器,其配置为通过基于轨道车辆的纵向加速度、横向加速度、横摆角速度以及车轮角速度估计横向速度的方式,来计算横向速度估计值;以及横向力估计观测器,其配置为通过基于轨道车辆的转向角、施加到轨道车辆的纵向力以及由横向速度估计观测器计算出的横向速度估计值估计施加到轨道车辆的转向架的横向力的方式,来计算横向力估计值。
在本公开的一些示例性实施例中,横向速度估计观测器可以包括:纵向速度计算器,其配置为基于由车轮传感器测量的前轮角速度和后轮角速度计算轨道车辆的纵向速度;以及横向速度估计器,其配置为基于纵向加速度、横向加速度以及由主体传感器测量的横摆角速度,并且基于由纵向速度计算器计算出的纵向速度,来计算横向速度估计值。
在本公开的一些示例性实施例中,横向速度估计观测器可以使用卡尔曼滤波器来计算横向速度估计值,并且横向力估计观测器可以使用扩展卡尔曼滤波器来计算横向力估计值。
在本公开的另一个总体方面中,提供了一种用于估计轨道车辆的横向力的方法,该方法包括:通过使用轨道车辆的前轮角速度和后轮角速度计算纵向速度;通过将轨道车辆的纵向加速度、横向加速度以及横摆角速度和上述纵向速度应用于卡尔曼滤波器,来计算横向速度估计值;以及通过将轨道车辆的转向角、施加到轨道车辆的纵向力、以及横向速度估计值应用于扩展卡尔曼滤波器估计施加到轨道车辆的转向架的横向力来计算横向力估计值。
根据本公开的一个示例性实施例,可以通过使用轨道车辆主体的动态模型和来自传感器的数据测量值来估计施加到轨道车辆的前和后转向架的横向力,而不需要任何复杂的横向力数学模型。
附图说明
图1为图示说明根据本公开的示例性实施例的用于估计轨道车辆的横向力的装置的方框图。
图2为图示说明根据本公开的示例性实施例的用于估计轨道车辆的横向力的装置的横向速度估计器的方框图。
图3为图示说明轨道车辆行驶在弯曲路段中的车辆模型的视图。
图4为图示说明轨道车辆的横向模型的自行车模型的视图。
具体实施方式
下文中,将参考附图详细地描述本公开的示例性实施例,以便本领域的技术人员可以实现和使用同样的实施例。为了清楚和方便描述,在此放大了在附图中图示说明的线条宽度和组件尺寸。此外,在下文提及的术语定义为考虑了在本公开中的功能,其可以根据用户或操作者、或实际消费者的意图改变。因此,应该基于本公开的整体内容作出术语的定义。
图1为图示说明根据本公开的示例性实施例的用于估计轨道车辆的横向力的装置的方框图;图2为图示说明根据本公开的示例性实施例的用于估计轨道车辆的横向力的装置的横向速度估计器的方框图;图3为图示说明轨道车辆行驶在弯曲路段中的车辆模型的视图;以及图4为图示说明轨道车辆的横向模型的自行车模型的视图。
参考图1,根据本公开的示例性实施例的一种用于估计轨道车辆的横向力的装置,可以包括横向速度估计观测器100和横向力估计观测器200。
横向速度估计观测器100可以通过基于轨道车辆的纵向加速度(ax)、横向加速度(ay)、横摆角速度(r)、以及车轮角速度(ωf,ωr)估计横向速度,来计算横向速度估计值。
这里,参考图2,横向速度估计观测器100可以包括纵向速度计算器110,其配置为基于由车轮传感器S1测量的前轮角速度(ωf)和后轮角速度(ωr)计算轨道车辆的纵向速度;以及横向速度估计器120,其配置为基于纵向加速度、横向加速度、以及由主体传感器S2测量的横摆角速度,并且基于由纵向速度计算器110计算出的纵向速度来计算横向速度估计值。
同时,横向力估计观测器200可以通过基于转向角(δ)、施加到车轮的纵向力(Fx1,Fx2,Fx3,Fx4)、以及由横向速度估计观测器100计算的横向速度估计值估计施加到转向架的横向力,来计算横向力估计值。
如上所述,横向速度估计观测器100计算横向速度估计值。下文中,将详细描述用于计算横向速度估计值的方法。
在图3所示的轨道车辆中心处的动力学动态特性可以由如下等式1表示。
[等式1]
其中vx和vy分别是在轨道车辆的质心处的纵向速度和横向速度,r是横摆角速度,并且ax和ay分别是纵向加速度和横向加速度。
上述等式1可以表示为如下面等式2的一种形式。
[等式2]
此外,假使***中存在扰动而将等式2表示为离散化等式时,等式2可以表示为下面的等式3。
[等式3]
x(k)=A(k-1)·x(k-1)+B(k-1)·u(k-1)+wd(k-1)
y(k)=C(k)·x(k)+wv(k)
其中
B(k-1)=ΔT
ΔT是测量间隔(步长),wd(k-1)和wv(k)分别表示在第k-1步中施加到***的扰动和在第k步中施加到输出的传感器噪声。
此外,假使可以测量在轨道车辆质心处的纵向速度,那么等式2可以表示为如下的等式4。
[等式4]
y(k)=vx(k)
C(k)=[1 0]
可以从前轮角速度和后轮角速度测量轨道车辆质心处的纵向速度。也就是,轨道车辆的纵向速度(vx(k))可以如下面等式5中作为前轮角速度和后轮角速度的平均值而被计算出。
[等式5]
其中ωf(k)和ωr(k)分别表示在第k步中的前轮角速度和后轮角速度,并且D表示车轮的直径。
因此,纵向速度计算器110可以基于上述等式5,使用由车轮传感器S2测量的前轮角速度和后轮角速度计算轨道车辆的纵向速度。
使用线性观测器来估计轨道车辆质心处的横向速度,并且具有各种类型的观测器用以估计线性***中的状态变量。在当前示例性实施例中,将横向速度估计器120设计为使用卡尔曼滤波器。
可以如下文所述对估计横向速度的线性卡尔曼滤波器进行设计。
首先,根据下面等式6估计状态变量的估计值。
[等式6]
其中是在第k-1步中的状态变量估计值,u(k-1)是在第k-1步中的输入估计值,并且是通过使用在第k-1步中的状态估计值、在第k-1步中的输入测量值等而预测出的第k个状态变量值。
继续地,使用下面等式7估计误差协方差。
[等式7]
P(k|k-1)=A(k-1)P(k-1|k-1)AT(k-1)+Q(k-1),
其中P(k-1|k-1)是误差协方差估计值,其中估计误差定义为实际状态变量与估计状态变量之间的差值。此外,Q(k-1)是施加到***的扰动wd(k-1)的协方差。P(k|k-1)是在第k步中通过使用***矩阵和扰动的协方差、以及在之前步骤中的估计误差协方差而预测出的状态变量的估计误差协方差。
接下来,使用下面等式8计算卡尔曼滤波器增益。
[等式8]
K(k)=P(k|k-1)CT(k)(C(k)P(k|k-1)CT(k)+R(k))-1,
其中K(k)是在第k步中的卡尔曼滤波器增益,并且R(k)是在第k步中的传感器测量噪声的协方差。
接下来,使用下面等式9校准状态变量。
[等式9]
其中y(k)是在第k步中的传感器测量值,并且是在第k步中的状态变量估计值。
当考虑它时,通过使用相对于来自在第k步中测量的值的输出变量的估计误差校准在第k-1步中预测的第k个状态变量值,来估计在第k步中的状态变量。
从而使用估计的状态变量,可以根据下面等式10计算在轨道车辆质心处的横向速度。
[等式10]
其中是在第k步中估计的轨道车辆的横向速度。
根据本公开的示例性实施例的横向力估计观测器200计算横向力估计值。下文中,将具体描述用于计算横向力估计值的方法。
图4为将图3的轨道车辆模型图示为自行车模型的视图。轨道车辆模型可以简化为自行车模型;因为可以假设当轨道车辆行驶在弯曲路段中时,施加到轨道车辆的左轮和右轮的力几乎相同。将描述其中具有四个轨道车辆的示例性示例。
沿着纵向方向、横向方向以及横摆方向的如图4所示自行车模型的轨道车辆动态模型分别如下面等式11至13。
[等式11]
[等式12]
[等式13]
其中∑Fx是施加到每个轨道车辆的纵向方向的力的总和,∑Fy是施加到每个轨道车辆的横向方向的力的总和,∑Fz是施加到每个轨道车辆的横摆方向的力的总和,并且每个力的总和(∑Fx,∑Fy,∑Fz)可以根据下面等式14计算。
[等式14]
当轨道车辆行驶在弯曲路段中,轨道车辆行驶在曲率为常数的轨道上。因此,可以假设每个轨道车辆的前轮以相同角度转向,并且后轮沿着相反方向以相同角度转向。因此,转向角可以如下面的等式15假设。
[等式15]
δ1=δ2=δ
δ3=δ4=-δ
因此,当应用等式14和15到等式11至13时,可以获得下面的等式16。
[等式16]
施加到轨道车辆的前转向架的横向力是施加到两个前轮的横向力的总和,并且施加到轨道车辆的后转向架的横向力是施加到两个后轮的横向力的总和。因此,施加到前和后转向架的横向力可以如下面等式17中那样定义。
[等式17]
lfFyf=l1Fy1+l2Fy2
lrFyr=l3Fy3+l4Fy4,
其中,lf是沿着纵向方向从轨道车辆中心到前轮转向架的长度,lr是沿着纵向方向从轨道车辆中心到后轮转向架的长度,Fyf是施加到前轮转向架的横向力,并且Fyr是施加到后轮转向架的横向力。此外,l1是沿着纵向方向从轨道车辆中心到第一前轮的长度,l2是沿着纵向方向从轨道车辆中心到第二前轮的长度,Fy1是施加到第一前轮的横向力,并且Fy2是施加到第二前轮的横向力。同样,l3是沿着纵向方向从轨道车辆中心到第一后轮的长度,l4是沿着纵向方向从轨道车辆中心到第二后轮的长度,Fy3是施加到第一后轮的横向力,并且Fy4是施加到第二后轮的横向力。
当用上述等式17替代等式15时,可以推导出下面等式18。
[等式18]
为了将等式18表示为状态等式,状态变量可以如在等式19中那样定义。
[等式19]
X1=vx
X2=vy
X3=r
X4=Fyf
X5=Fyr
附加地,当假设施加到前轮和后轮转向架的横向力的值缓慢改变时,可以假设横向力几乎为恒定常数。因此,可以假设横向力的差值为零(0)。
[等式20]
当再次使用等式19和20表示等式18时,可以推导出下面等式21。
[等式21]
当离散化等式21时,它可以如下面等式22中那样表示。
[等式22]
假设***中存在扰动并且在测量时产生传感器噪声,当将等式22重新定义为状态等式时,它可以如在等式23中那样表示。
[等式23]
X(k)=f(X(k-1),U(k-1))+wd(k-1)
Y(k)=h(X(k))+wv(k),
其中
wd(k-1)是施加到***的扰动,并且wv(k)是测量的噪声。
如在上述等式23中确定,将纵向速度、横向速度、沿着轨道车辆中心施加的横摆角速度、以及施加到前和后轮转向架的横向力定义为状态变量。此外,将在轨道车辆的质心处的纵向速度、在轨道车辆的质心处估计的横向速度、以及在轨道车辆的质心处的横摆角速度定义为测量变量。
在本公开的示例性实施例中将扩展卡尔曼滤波器使用作为横向力估计观测器200。然而,这仅是为描述本公开的一个示例。因此,本公开所属领域的技术人员可以清楚,也可以使用其他类型的观测器来估计施加到轨道车辆的转向架的横向力。
可以通过下面等式24计算使用扩展卡尔曼滤波器估计施加到转向架的横向力的状态变量值。
[等式24]
其中是在第k-1步中的状态变量估计值,U(k-1)是在第k-1步中的输入测量值。此外,是通过使用在第k-1步中的状态变量估计值、在第k-1步中的输入测量值等预测的第k个状态变量值。
同时,可以通过下面等式25获得在第k步中预测的状态变量的估计误差协方差(P(k|k-1))。
[等式25]
P(k|k-1)=F(k-1)P(k-1|k-1)F(k-1)T+Q(k-1),
其中其定义为函数f(X(k),U(k))的关于X(k)的雅克比(Jacobian)矩阵。
此外,P(k-1|k-1)是在第k-1步中的已估计的误差协方差估计值,并且该已估计的误差定义为在实际状态变量与已估计的状态变量之间的差值。此外,Q(k-1)是作为施加到***的扰动wd(k-1)的协方差,并且P(k|k-1)是在第k步中通过使用***矩阵、扰动的协方差、以及在之前步骤中预测出的状态变量的已估计的误差协方差值,来预测的状态变量的已估计的误差协方差。
同时,根据下面等式26,可以基于由等式24计算的状态变量值等式估计测量变量值。
[等式26]
此外,可以通过下面等式27计算在第k步中的卡尔曼滤波器增益(L(k))。
[等式27]
L(k)=P(k|k-1)H(k)T(H(k)P(k|k-1)H(k)T+R(k))-1,
其中R(k)是在第k步中的传感器测量噪声的协方差。
此外,可以通过下面等式28计算状态变量估计值。
[等式28]
其中Y(k)是在第k步中的传感器测量值,并且是在第k步中的状态变量估计值。
当考虑它时,通过使用关于来自在第k步中测量的值的输出变量的估计误差校准在第k-1步中预测的第k个状态变量值,来估计在第k步中的状态变量。
此外,使用通过等式25预测的状态变量值的已估计的误差协方差以及通过等式27计算的卡尔曼滤波器增益更新的已估计协方差(P(k|k))可根据下面等式29计算。
[等式29]
P(k|k)=(I-L(k)H(k))P(k|k-1),
其中其定义为函数h(X(k))的关于X(k)的雅克比(Jacobian)矩阵。
同样,通过使用在等式24至29中定义的扩展卡尔曼滤波器来估计状态变量。此外,可以如在等式30中那样,使用在第k步中估计的状态变量值来估计施加到轨道车辆的前和后轮转向架的横向力。
[等式30]
其中是在第k步中的状态变量估计值,是在第k步中施加到前轮转向架的横向力的估计值,并且是在第k步中施加到后轮转向架的横向力的估计值。
同时,虽然一种根据本公开示例性实施例的用于估计轨道车辆的横向力的装置及方法已经在上文中描述,然而,本公开的范围不受上述实施例的限制。因此,本公开可以在该限制范围内替代地以各种变形或修改来实施,以便对于本公开所属领域的普通技术人员来说那些差异是显而易见的。
因此,上述示例性实施例和附图旨在进行图示说明,而不是限制权利要求的范围。本公开的保护范围由下面的权利要求解释,并且应该包括在本公开范围的等同范围内的所有技术思想。
Claims (5)
1.一种用于估计轨道车辆的横向力的装置,所述装置包括:
横向速度估计观测器,其配置为通过基于所述轨道车辆的纵向加速度、横向加速度、横摆角速度、以及车轮角速度估计横向速度,来计算横向速度估计值;以及
横向力估计观测器,其配置为通过基于所述轨道车辆的转向角、施加到所述轨道车辆的纵向力、以及由所述横向速度估计观测器计算出的横向速度估计值估计施加到所述轨道车辆的转向架的横向力,来计算横向力估计值,
其中所述横向速度估计观测器包括:
纵向速度计算器,其配置为基于由车轮传感器测量的前轮角速度和后轮角速度,来计算所述轨道车辆的纵向速度;以及
横向速度估计器,其配置为基于由主体传感器测量的所述纵向加速度、所述横向加速度和所述横摆角速度,并且基于由所述纵向速度计算器计算出的纵向速度,来计算所述横向速度估计值,
其中所述纵向速度通过下面等式作为所述前轮角速度和所述后轮角速度的平均值而被计算出:
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mi>D</mi>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
其中ωf(k)和ωr(k)分别表示在第k步中的前轮角速度和后轮角速度,并且D表示车轮的直径,其中通过下面等式计算所述横向速度估计值
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中是在第k步估计的所述轨道车辆的横向速度估计值,并且是在第k步中的状态变量估计值。
2.如权利要求1所述的装置,其中
所述横向速度估计观测器使用卡尔曼滤波器来计算所述横向速度估计值,以及
所述横向力估计观测器使用扩展卡尔曼滤波器来计算所述横向力估计值。
3.如权利要求1所述的装置,其中通过下面等式计算所述横向力估计值
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中是在第k步中的状态变量估计值,是在第k步中施加到前轮转向架的横向力的估计值,以及是在第k步中的施加到后轮转向架的横向力的估计值。
4.一种用于估计轨道车辆的横向力的方法,所述方法包括:
通过计算所述轨道车辆的前轮角速度和后轮角速度的平均值来计算纵向速度;
通过将所述轨道车辆的纵向加速度、横向加速度、横摆角速度、以及计算出的所述纵向速度应用到卡尔曼滤波器来计算横向速度估计值;以及
通过将所述轨道车辆的转向角、施加到所述轨道车辆车轮的纵向力、以及所述横向速度估计值应用到扩展卡尔曼滤波器来估计施加到所述轨道车辆的转向架的横向力,来计算横向力估计值,
其中所述纵向速度通过下面等式作为所述前轮角速度和所述后轮角速度的平均值而被计算出:
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mi>D</mi>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
其中ωf(k)和ωr(k)分别表示在第k步中的前轮角速度和后轮角速度,并且D表示车轮的直径,
其中,通过如下等式,使用在第k步估计的状态变量计算所述横向速度估计值,通过使用相对于第k-1步中预测的第k个状态变量估计值与在第k步中测量的值之间的输出变量的估计误差来校准所述状态变量:
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中是在第k步估计的所述轨道车辆的横向速度估计值;并且是在第k步中的状态变量估计值。
5.如权利要求4所述的方法,其中通过计算在第k步中施加到前轮转向架的横向力的估计值和在第k步中施加到后轮转向架的横向力的估计值,通过将在第k步中的状态变量估计值应用到下面等式,来计算所述横向力估计值:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中是在第k步中的状态变量估计值,是在第k步中施加到前轮转向架的横向力的估计值,以及是在第k步中的施加到后轮转向架的横向力的估计值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR10-2014-0009343 | 2014-01-27 | ||
KR1020140009343A KR101870482B1 (ko) | 2014-01-27 | 2014-01-27 | 철도 차량의 횡력 추정장치 및 추정방법 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104802826A CN104802826A (zh) | 2015-07-29 |
CN104802826B true CN104802826B (zh) | 2018-02-13 |
Family
ID=53678297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510088778.7A Expired - Fee Related CN104802826B (zh) | 2014-01-27 | 2015-01-26 | 用于估计轨道车辆的横向力的装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9637148B2 (zh) |
KR (1) | KR101870482B1 (zh) |
CN (1) | CN104802826B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3434552B1 (en) * | 2016-03-23 | 2021-05-05 | Nippon Steel Corporation | Inspection system, inspection method and program |
ITUA20162317A1 (it) * | 2016-04-05 | 2017-10-05 | Corghi Spa | Apparato e metodo per la valutazione dell'assetto di un veicolo. |
CN105946858B (zh) * | 2016-06-08 | 2019-02-15 | 吉林大学 | 基于遗传算法的四驱电动汽车状态观测器参数优化方法 |
CN111065567B (zh) * | 2017-08-31 | 2021-04-13 | 日本制铁株式会社 | 检查***、检查方法及计算机可读取存储介质 |
CN107966307B (zh) * | 2017-12-14 | 2020-06-02 | 常州路航轨道交通科技有限公司 | 一种基于振动加速度的列车运行安全在线监测方法 |
CN109855770B (zh) * | 2018-12-31 | 2021-08-20 | 华东交通大学 | 一种基于钢轨横向加速度功率谱密度检测温度力的方法 |
CN112722003B (zh) * | 2019-10-14 | 2022-12-13 | 北京昊鹏智能技术有限公司 | 列车脱轨风险的监测方法与设备 |
CN111060335B (zh) * | 2020-01-02 | 2021-07-02 | 中车青岛四方车辆研究所有限公司 | 一种主动轮对连接机构的试验方法及试验*** |
KR102465763B1 (ko) * | 2020-06-29 | 2022-11-15 | 한국철도기술연구원 | 철도차량 탈선계수 측정장치 및 측정방법 |
CN112061178B (zh) * | 2020-09-18 | 2022-07-05 | 北京世纪东方智汇科技股份有限公司 | 一种列尾设备脱落状态检测方法及装置 |
CN114212104B (zh) * | 2021-12-14 | 2024-06-18 | 京东鲲鹏(江苏)科技有限公司 | 车辆控制方法、装置、车辆和存储介质 |
JP7332008B1 (ja) | 2022-09-09 | 2023-08-23 | いすゞ自動車株式会社 | 車両制御装置及び車両制御方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100352703C (zh) * | 2002-03-19 | 2007-12-05 | 汽车***实验室公司 | 车辆倾翻检测*** |
CN101417654A (zh) * | 2007-05-07 | 2009-04-29 | 通用汽车环球科技运作公司 | 用于估计车辆状态以便避免倾翻的*** |
CN102009653A (zh) * | 2010-11-10 | 2011-04-13 | 清华大学 | 融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法 |
CN102556075A (zh) * | 2011-12-15 | 2012-07-11 | 东南大学 | 一种基于改进扩展卡尔曼滤波的车辆运行状态估计方法 |
CN102673569A (zh) * | 2012-05-25 | 2012-09-19 | 同济大学 | 车辆状态测算装置、方法及使用该装置的车辆 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9233696B2 (en) * | 2006-03-20 | 2016-01-12 | General Electric Company | Trip optimizer method, system and computer software code for operating a railroad train to minimize wheel and track wear |
KR100521169B1 (ko) * | 2002-12-27 | 2005-10-12 | 현대자동차주식회사 | 롤 오버 제어 방법 |
US7451032B2 (en) | 2004-06-02 | 2008-11-11 | Ford Global Technologies, Llc | System and method for determining desired yaw rate and lateral velocity for use in a vehicle dynamic control system |
JP4285548B2 (ja) | 2007-02-05 | 2009-06-24 | エプソントヨコム株式会社 | ジャイロセンサモジュールおよび角速度検出方法 |
US7835846B2 (en) * | 2007-02-28 | 2010-11-16 | Gm Global Technology Operations, Inc. | Nonlinear vehicle yaw/roll/sideslip command interpreter |
US8050838B2 (en) | 2008-11-24 | 2011-11-01 | GM Global Technology Operations LLC | Kinematic estimator for vehicle lateral velocity using force tables |
DE102010003205B4 (de) | 2010-03-24 | 2020-01-16 | Ford Global Technologies, Llc | Verfahren zur Bestimmung der vertikalen Beschleunigung, der longitudinalen Winkelbeschleunigung und der transversalen Winkelbeschleunigung eines Körpers, insbesondere eines Kraftfahrzeugs |
-
2014
- 2014-01-27 KR KR1020140009343A patent/KR101870482B1/ko active IP Right Grant
-
2015
- 2015-01-20 US US14/601,024 patent/US9637148B2/en not_active Expired - Fee Related
- 2015-01-26 CN CN201510088778.7A patent/CN104802826B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100352703C (zh) * | 2002-03-19 | 2007-12-05 | 汽车***实验室公司 | 车辆倾翻检测*** |
CN101417654A (zh) * | 2007-05-07 | 2009-04-29 | 通用汽车环球科技运作公司 | 用于估计车辆状态以便避免倾翻的*** |
CN102009653A (zh) * | 2010-11-10 | 2011-04-13 | 清华大学 | 融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法 |
CN102556075A (zh) * | 2011-12-15 | 2012-07-11 | 东南大学 | 一种基于改进扩展卡尔曼滤波的车辆运行状态估计方法 |
CN102673569A (zh) * | 2012-05-25 | 2012-09-19 | 同济大学 | 车辆状态测算装置、方法及使用该装置的车辆 |
Also Published As
Publication number | Publication date |
---|---|
US20150210300A1 (en) | 2015-07-30 |
KR101870482B1 (ko) | 2018-06-22 |
KR20150089127A (ko) | 2015-08-05 |
US9637148B2 (en) | 2017-05-02 |
CN104802826A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104802826B (zh) | 用于估计轨道车辆的横向力的装置 | |
Reina et al. | Vehicle parameter estimation using a model-based estimator | |
Singh et al. | Literature review and fundamental approaches for vehicle and tire state estimation | |
Zhao et al. | Design of a nonlinear observer for vehicle velocity estimation and experiments | |
Stephant et al. | Virtual sensor: Application to vehicle sideslip angle and transversal forces | |
Boada et al. | Vehicle sideslip angle measurement based on sensor data fusion using an integrated ANFIS and an Unscented Kalman Filter algorithm | |
Doumiati et al. | A method to estimate the lateral tire force and the sideslip angle of a vehicle: Experimental validation | |
Doumiati et al. | Observers for vehicle tyre/road forces estimation: experimental validation | |
CN102112354B (zh) | 路面摩擦系数估计装置和路面摩擦系数估计方法 | |
Doumiati et al. | Lateral load transfer and normal forces estimation for vehicle safety: experimental test | |
US20140371990A1 (en) | Sensor system comprising a vehicle model unit | |
US20190368878A1 (en) | Method for determining an orientation of a vehicle | |
CN103909933B (zh) | 一种分布式电驱动车辆的前轮侧向力估算方法 | |
CN105835889B (zh) | 一种基于二阶滑模观测器的车辆质心侧偏角的估计方法 | |
Ahn | Robust Estimation of Road Friction Coefficient for Vehicle Active Safety Systems. | |
EP3309024A1 (en) | Method and system for determining friction between the ground and a tire of a vehicle | |
Madhusudhanan et al. | Load sensing bearing based road-tyre friction estimation considering combined tyre slip | |
JP2011128093A5 (zh) | ||
Han et al. | Monitoring system design for lateral vehicle motion | |
Wielitzka et al. | State estimation of vehicle's lateral dynamics using unscented Kalman filter | |
KR20130133414A (ko) | 열차 주행속도 보정장치 | |
CN103279675A (zh) | 轮胎-路面附着系数与轮胎侧偏角的估计方法 | |
JP2009073466A (ja) | 車両姿勢角推定装置及びプログラム | |
Ding et al. | An extended Luenberger observer for estimation of vehicle sideslip angle and road friction | |
Jiang et al. | Real-time estimation and prediction of tire forces using digital map for driving risk assessment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180213 Termination date: 20210126 |