CN110779743B - 基于车载悬浮控制器的磁浮列车轨排监测方法 - Google Patents

基于车载悬浮控制器的磁浮列车轨排监测方法 Download PDF

Info

Publication number
CN110779743B
CN110779743B CN201911098711.6A CN201911098711A CN110779743B CN 110779743 B CN110779743 B CN 110779743B CN 201911098711 A CN201911098711 A CN 201911098711A CN 110779743 B CN110779743 B CN 110779743B
Authority
CN
China
Prior art keywords
track panel
suspension
levitation
electromagnet
track
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911098711.6A
Other languages
English (en)
Other versions
CN110779743A (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.)
National University of Defense Technology
CRRC Tangshan Co Ltd
Original Assignee
National University of Defense Technology
CRRC Tangshan Co Ltd
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 National University of Defense Technology, CRRC Tangshan Co Ltd filed Critical National University of Defense Technology
Priority to CN201911098711.6A priority Critical patent/CN110779743B/zh
Publication of CN110779743A publication Critical patent/CN110779743A/zh
Application granted granted Critical
Publication of CN110779743B publication Critical patent/CN110779743B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/08Railway vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/005Testing of electric installations on transport means
    • G01R31/008Testing of electric installations on transport means on air- or spacecraft, railway rolling stock or sea-going vessels

Abstract

本发明公开了一种基于车载悬浮控制器的磁浮列车轨排状况监测方法,完全嵌入在磁浮列车悬浮控制器中,根据磁浮列车自身配置的传感器获得悬浮间隙、电磁铁垂向加速度、电磁铁电流信息,利用实时估计方法估计出悬浮单元所处位置的轨排垂向位移,然后采用自适应辨识器实时估计轨排的动力学参数,由此构建出轨排的动力学模型传递函数,得出轨排的模态频率、阻尼系数、柔度系数等刻画轨排结构状况的关键参数。这些参数既可以作为悬浮系统自适应振动控制算法的参考输入,又可以为轨排的日常监测和维护提供参考。相比现有技术,该方法无需专用的轨检设备,也不依赖轨排振动幅值检测,具有结构简单、成本低廉、应用范围广且检测准确的优点。

Description

基于车载悬浮控制器的磁浮列车轨排监测方法
技术领域
本发明涉及轨道交通的轨道检测技术领域,特别是涉及一种基于车载悬浮控制器的磁浮列车轨排监测方法。
背景技术
电磁吸力型磁浮列车通过通电的电磁铁产生磁场吸引轨道来产生悬浮力,然而这种悬浮力的大小与电磁铁和轨道之间的距离成平方反比关系,因而使得吸力型悬浮结构自身是不稳定的,需要外加主动控制来实时调节电磁铁的电流方能使悬浮系统稳定工作。然而这种主动控制的引入使得磁浮车的悬浮系统和轨道之间的耦合问题变得复杂,当轨道的刚度不足、或者出现结构松动、部件疲劳等问题时,易引发电磁铁和轨道之间的耦合振动问题。在实际运营的中低速磁浮线路中发现,由轨排的支撑结构松动、部件老化等原因容易引发高频自激振动问题,而这种振动又会加剧轨排结构的松动、甚至引发螺栓断裂等故障,影响行车安全。因此,对轨排结构的监测是商业磁浮线路维护中的重要内容。然而,当前的轨道检测设备多以轨道几何不平顺为检测对象,无法检测螺栓松动及可能引发轨排自激振动的潜在因素。
中国专利CN201810763591.6公开了一种带定位功能的中低速磁浮F型轨道检测仪,该检测仪由支架单元、滑动单元、测量单元、定位和数据传输单元构成,是一种独立的轨道检测设备,用于检测磁浮轨排中左右两侧的F型钢轨的几何形状。但是这种方法仅对F型钢轨的几何平顺度进行检测,无法检测螺栓松动以及由此可能引发的振动等涉及轨排状况的隐性问题。中国专利CN201610548735.7公开了一种基于四传感器的中低速磁浮轨道不平顺检测方法,该方法在悬浮架沿轨道运动方向上设置四个间隙传感器,通过四个间隙测点拟合直线,再根据四个间隙测点与第一拟合直线的偏差获取所处位置的四点直线度。这种方法可以在磁浮车的悬浮架上实现,但其同样只能检测轨道的直线度、错台、折角等几何不平顺参数。中国专利CN201811253700.6公开了一种用于中低速磁浮列车的轨排振动检测装置,该装置包括设于中低速磁浮列车上的车载机箱、嵌设于中低速磁浮列车的悬浮控制器中的轨道检测模块、以及设于悬浮控制器和车载机箱的第一和第二通信模块,用来检测轨排安装误差以及高频振动,并对超差或振动轨排进行定位,以指导轨排维护和保障行车安全。但该装置需要在磁浮车的电器柜中设置专用的车载机箱,其检测需要在磁浮车和轨排发生显著的耦合振动时才能对振动频率和悬浮电磁铁加速度幅值进行估计,而振动的幅值同悬浮控制器的控制参数依赖度很大,且加速度幅值的估计易受车辆行进时电磁铁的正常波动、传感器噪声干扰的影响。
故此,如何研发一种结构简单、易于实现、节约能源且适用范围广的基于车载悬浮控制器的磁浮列车轨排监测方法,成了本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提出了一种基于车载悬浮控制器的磁浮列车轨排监测方法,该方法完全不依赖轨排振动幅值检测,具有结构简单、成本低廉、应用范围广且检测准确的优点。
一方面,本发明提供了一种基于车载悬浮控制器的磁浮列车轨排监测方法,包括如下步骤:
S100、根据磁浮列车自身配置的悬浮传感器采集电磁铁和轨排之间的悬浮间隙、电磁铁垂向加速度以及通过电流传感器采集电磁铁电流信息;
S200、利用实时估计方法估计出悬浮单元所处位置的轨排垂向位移和悬浮力;
S300、采用自适应辨识器实时估计轨排的动力学参数,并构建轨排的动力学模型传递函数,得出用于刻画轨排结构状况的关键参数;
S400、通过所述关键参数对磁浮列车轨排进行日常监测和维护。
进一步地,所述步骤S200中轨排垂向位移通过如下公式进行估计:
Figure BDA0002269169800000031
式中,
Figure BDA0002269169800000032
为估计得到的轨排垂向位移,x1、x2、x3和x4均为状态变量,
Figure BDA0002269169800000033
Figure BDA0002269169800000034
Figure BDA0002269169800000035
分别为状态变量x1、x2、x3和x4的一阶导数,σ1为采集的加速度和间隙数据低频特征参数,σ2为采集的加速度和悬浮间隙数据高频特征参数,δ1为电磁铁和轨排之间的间隙,c1为电磁铁的垂向加速度。
进一步地,所述步骤S200中悬浮单元的悬浮力的表达公式如下:
Figure BDA0002269169800000036
式中,i1表示悬浮电磁铁的电流,i10表示悬浮单元的稳态悬浮电流,z0是稳态悬浮间隙,kf和kz是与悬浮系统的悬浮间隙、电磁铁尺寸相关的参数,可取稳态悬浮间隙情况下的标称值,其中:
Figure BDA0002269169800000037
式中,A是电磁铁的磁极面积,N是电磁铁的线圈匝数,μ0是空气磁导率。
进一步地,所述步骤S300采用自适应辨识器实时估计轨排的动力学参数,并构建轨排的动力学模型传递函数,得出用于刻画轨排结构状况的关键参数,具体表现为:
S301、更新观测向量:
Figure BDA0002269169800000038
式中,
Figure BDA0002269169800000039
表示第n步的观测向量,上标T为向量的转置,
Figure BDA00022691698000000310
表示悬浮单元第n步的悬浮力、
Figure BDA00022691698000000311
表示悬浮单元第n–1步的悬浮力、
Figure BDA00022691698000000312
表示悬浮单元第n–2步的悬浮力,
Figure BDA0002269169800000041
表示轨排第n–1步估计得到的垂向位移、
Figure BDA0002269169800000042
表示轨排第n–2步估计得到的垂向位移,n表示当前迭代的步数,可取任意一个整数;
S302、更新增益向量:
Figure BDA0002269169800000043
式中,K(n)是5×1的增益向量,P(n-1)是第n-1步的误差协方差矩阵,γ是遗忘因子,取0.9<γ<1;
S303、更新状态估计向量:
Figure BDA0002269169800000044
式中,
Figure BDA0002269169800000045
为第n+1步的状态估计向量,
Figure BDA0002269169800000046
为第n步的状态估计向量,K(n+1)是5×1的增益向量,
Figure BDA0002269169800000047
表示轨排第n步估计得到的电磁铁垂向位移;
S304、更新协方差矩阵:
Figure BDA0002269169800000048
式中,P(n+1)是第n+1步的误差协方差矩阵,P(n)是第n步的误差协方差矩阵,I是5×5的单位矩阵,
Figure BDA0002269169800000049
表示第n+1步的观测向量,上标T为向量的转置;
S305,得到轨排的模型传递函数的实时估计参数:
Figure BDA00022691698000000410
式中,
Figure BDA00022691698000000411
表示估计的模型传递函数,
Figure BDA00022691698000000412
Figure BDA00022691698000000413
分别表示状态估计向量
Figure BDA00022691698000000414
的第1个、第2个、第3个、第4个、第5个分量。
进一步地,所述轨排的关键参数包括模态频率、阻尼系数和柔度系数。
进一步地,所述步骤S100之前还包括如下步骤:
S00A、变量初始化:创建初始观测向量
Figure BDA00022691698000000415
状态估计向量
Figure BDA00022691698000000416
以及估计误差协方差矩阵P(0),其中
Figure BDA0002269169800000051
P(0)为5×5的矩阵,且矩阵P(0)的元素初始值取较大的值。
进一步地,所述步骤S00A和步骤S100之间还包括如下步骤:
S00B、判断悬浮系统的悬浮状态:若悬浮系统处于悬浮状态,则进入步骤S100,否则,停止磁浮列车轨排监测流程。
进一步地,所述磁浮列车轨排检测结果包括如下类别:
对于辨识出的轨排的模型传递函数的幅频特性峰值超出设定的第一阈值,则判断轨排在该位置的柔度系数较大,存在支撑不良情况,用于提示该处轨排需要进一步进行维护保养;
对于辨识出的轨排的模型传递函数的幅频特性峰值对应的频率小于设定的第二阈值,且幅频特性峰值对应的阻尼系数小于第三阈值,则判断该轨排存在螺栓松动、支撑刚度下降等问题,用于提示该处轨排需要进一步进行维护保养;
进一步地,得到的轨排实时估计参数同时还可以作为悬浮控制自适应振动控制算法的参考输入,以提高悬浮控制性能。
本发明提供的基于车载悬浮控制器的磁浮列车轨排监测方法,相比现有技术,不依赖轨排振动幅值检测,可以完全嵌入在磁浮列车悬浮控制器中,并根据磁浮列车自身配置的传感器获得悬浮间隙、电磁铁垂向加速度、电磁铁电流信息,利用实时估计方法估计出悬浮单元所处位置的轨排垂向位移,然后采用自适应辨识器实时估计轨排的动力学参数,由此构建出轨排的动力学模型传递函数,得出轨排的模态频率、阻尼系数、柔度系数等刻画轨排结构状况的关键参数,这些参数的获得既可以作为悬浮系统自适应振动控制算法的参考输入,又可以为轨排的日常监测和维护提供参考,具有结构简单、成本低廉、应用范围广且检测准确的优点。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为简化的悬浮模块和轨排结构关系示意图;
图2为轨排、悬浮系统构成的闭环耦合系统;
图3为轨排、悬浮系统构成的闭环耦合系统简化框图;
图4为本发明一实施例提供的基于车载悬浮控制器的磁浮列车轨排监测方法的流程图;
图5为中低速磁浮轨排结构;
图6为图5所示轨排单侧F轨的模态振型;
图7为图5所示轨排振动发生时两个悬浮单元的悬浮间隙曲线图;
图8为图5所示轨排参数辨识过程;
图9为设定的轨排模型的频率响应特性图;
图10为采用本发明的估计方法得到的轨排频率响应特性图。
附图标记说明
1-轨排 2-电磁铁
3-第一悬浮传感器 4-第二悬浮传感器
5-第一控制器 6-第二控制器
11-F轨 12-轨枕。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
为更好地理解本发明,在阐述本发明基于车载悬浮控制器的磁浮列车轨排监测方法的具体流程之前,先对本发明的基本原理展开论述:
图1为磁浮列车悬浮电磁铁和轨排的典型结构关系示意图。图中只画出了一个悬浮模块(包含左右两个悬浮单元)的示意图,为了清晰起见,隐去了悬浮架上连接悬浮电磁铁的托臂、直线电机、空气弹簧、车体等部件。悬浮电磁铁位于轨排1的下方,在电磁铁2的端部分别安装有第一悬浮传感器3和第二悬浮传感器4用于测量电磁铁2距离轨排1的垂向距离,以及电磁铁2的垂向加速度,并将测量的上述信息分别反馈至第一控制器5和第二控制器6处。
构成轨排1的主要部件,如两侧的导轨、中间的轨枕、以及轨枕下面的支撑垫通常都是有弹性的,因此轨排1存在自身对应的模态振型和模态频率。对于轨排的第k阶模态振型来说,其受迫振动的响应可以表示为:
Figure BDA0002269169800000071
式中,qk(t)是轨排的第k阶振动模态的广义位移,
Figure BDA0002269169800000072
为qk(t)的一阶导数,
Figure BDA0002269169800000073
为qk(t)的二阶导数,ζk为第k阶模态阻尼系数,ωk表示第k阶模态频率,Qk(t)是作用在该阶模态上的广义力,由作用在同一段轨排上的所有电磁力产生,Mk是该阶模态的广义质量,且
Figure BDA0002269169800000074
式中,ρ是轨排的线密度,φk(x)为轨排的k阶模态振型。
由于轨排的振动位移通常很小,因此可以认为作用在轨排上的电磁力是均布力,于是有
Qk(t)=α1F1(t)+α2F2(t)+…+αmFm(t) (3)
其中,αm表示第m个电磁力对考察的轨段的广义力的贡献程度,其表达式为:
Figure BDA0002269169800000075
式中,xm表示电磁力所处的起始位置,lm表示电磁力所处的覆盖长度。
对于第m个悬浮控制单元来说,假定其悬浮传感器所处的位置为x0m,则该处的轨排垂向位移可以表示为:
y0m(t)=qk(t)βm (5)
其中βm=φk(x0m),轨排振动产生的垂向变形位移会叠加在悬浮传感器的检测信号中,从而对悬浮控制系统产生激励,使得悬浮电磁力中包含振动激励的响应分量:
Fm(s)=Gm(s)y0m(s) (6)
式中,Gm(s)是第m个悬浮控制单元从轨排位移激励到电磁力输出的传递函数,s是拉普拉斯算子。在磁浮列车中,由于所有的悬浮控制器参数和电磁铁结构、标称悬浮间隙都是一样的,因此可以认为G1(s)=G2(s)=…=Gm(s)。
而由式(1)可得在广义力作用下轨排第k阶模态的垂向位移响应:
qk(s)=H(s)Fk(s) (7)
其中
Figure BDA0002269169800000081
公式(5)~(8)表征的轨排、悬浮系统相互耦合组成一个闭环系统,具体如图2所示。
下面以第一个悬浮单元为例,介绍在磁浮列车上基于单个悬浮单元对轨排参数进行实时辨识的方法。在图2中,由于所有的悬浮单元的参数是一致的,因此可以将图中除第一个悬浮单元外的其它悬浮单元进行合并,得到图3所示的简化框图。图中
Figure BDA0002269169800000082
它们代表这些悬浮单元的电磁力综合起来的增益大小。
对于式(8)的表述的轨排二阶系统来说,其离散化的形式可以表示为
Figure BDA0002269169800000083
写成差分方程的形式,得
qk(n)=b0Qk(n)+b1Qk(n-1)+b2Qk(n-2)-a1qk(n-1)-a2qk(n-2) (10)
由式(3)、(5),可得
y01(n)=α1β1b0F1(n)+α1β1b1F1(n-1)+α1β1b2F1(n-2)
0β0b0F0(n)+α0β0b1F0(n-1)+α0β0b2F0(n-2)-a1y01(n-1)-a2y01(n-2) (11)
由于α0、β0皆为实数,因此F0和F1的相位是相同的,仅幅值不同,因此可以把式(11)改写为
y01(n)=ηb0F1(n)+ηb1F1(n-1)+ηb2F1(n-2)-a1y01(n-1)-a2y01(n-2) (12)
式中,η=α0β01β1
定义向量
Figure BDA0002269169800000091
θ=[ηb0 ηb1 ηb2 -a1 -a2]T (14)
其中,上标T为向量的转置。则
Figure BDA0002269169800000092
式中,e(n)表示在实际的磁浮列车中因测量误差、噪声干扰等因素引起的综合误差,为不可测量。
利用式(15)就可以辨识出对应的系数,进而得到式(9)描述的轨排模型。需要指出的是,由于η是个未知的参数,其大小与作用在轨排上的所有悬浮单元所处的位置有关,因此该参数无法辨识。但是由式(15)辨识得到的结果与式(9)描述的轨排模型仅相差一个系数,不影响对轨排基本参数(如不稳定模态频率、模态阻尼系数等)的辨识。
有鉴于此,当磁浮列车处于悬浮状态时,选择其中的一台悬浮控制器,通过如下步骤进行基于车载悬浮控制器的磁浮列车轨排监测:
S000、实施准备,
S00A、变量初始化。创建初始观测向量
Figure BDA0002269169800000095
状态估计向量
Figure BDA0002269169800000096
以及估计误差协方差矩阵P(0),其中
Figure BDA0002269169800000093
Figure BDA0002269169800000094
P(0)为5×5的矩阵,矩阵P(0)的元素初始值可以取较大的值,并设置遗忘因子γ,可取0.9<γ<1。
S00B、判断悬浮系统的悬浮状态:若悬浮系统处于悬浮状态,则进入步骤S100,否则,停止磁浮列车轨排监测流程。
S100、数据采集。通过磁浮列车自身配置的悬浮传感器采集电磁铁和轨排之间的悬浮间隙数据δ1、悬浮电磁铁的垂向加速度a1,以及通过电流传感器采集的电磁铁电流i1
S200、状态估计。通过采集到的数据利用实时估计该悬浮单元的轨排垂向位移
Figure BDA0002269169800000101
和悬浮力
Figure BDA0002269169800000102
Figure BDA0002269169800000103
Figure BDA0002269169800000104
式中,
Figure BDA0002269169800000105
为估计得到的轨排垂向位移,x1、x2、x3和x4均为状态变量,
Figure BDA0002269169800000106
Figure BDA0002269169800000107
Figure BDA0002269169800000108
分别为状态变量x1、x2、x3和x4的一阶导数,σ1为采集的加速度和悬浮间隙数据低频特征参数,σ2为采集的加速度和悬浮间隙数据高频特征参数,δ1为电磁铁和轨排之间的间隙,c1为悬浮电磁铁的垂向加速度,i1表示悬浮电磁铁的电流,i10表示悬浮单元的稳态悬浮电流,z0是稳态悬浮间隙,kf和kz是与悬浮系统的悬浮间隙、电磁铁尺寸相关的参数,可取稳态悬浮间隙情况下的标称值,其中:
Figure BDA0002269169800000109
式中,A是电磁铁的磁极面积,N是电磁铁的线圈匝数,μ0是空气磁导率。
需要说明的是,上述公式中,σ1可选择远小于振动圆频率的正值,σ2则应选择两倍于振动圆频率以上的值,μ0=4π×10-7H/m。
S301、更新观测向量:
Figure BDA0002269169800000111
式中,
Figure BDA0002269169800000112
表示第n步的观测向量,上标T为向量的转置,
Figure BDA0002269169800000113
表示悬浮单元第n步的悬浮力、
Figure BDA0002269169800000114
表示悬浮单元第n–1步的悬浮力、
Figure BDA0002269169800000115
表示悬浮单元第n–2步的悬浮力,
Figure BDA0002269169800000116
表示轨排第n–1步的估计得到垂向位移、
Figure BDA0002269169800000117
表示轨排第n–2步估计得到的垂向位移,n表示当前迭代的步数,可取任意一个整数。
S302、更新增益向量:
Figure BDA0002269169800000118
式中,K(n)是5×1的增益向量,P(n-1)是第n-1步的误差协方差矩阵,γ是遗忘因子,取0.9<γ<1;
S303、更新状态估计向量:
Figure BDA0002269169800000119
式中,
Figure BDA00022691698000001110
为第n+1步的状态估计向量,
Figure BDA00022691698000001111
为第n步的状态估计向量,K(n+1)是5×1的增益向量,
Figure BDA00022691698000001112
表示轨排第n步估计得到的电磁铁垂向位移;
S304、更新协方差矩阵:
Figure BDA00022691698000001113
式中,P(n+1)是第n+1步的误差协方差矩阵,P(n)是第n步的误差协方差矩阵,I是5×5的单位矩阵,
Figure BDA00022691698000001114
表示第n+1步的观测向量,上标T为向量的转置;
S305,得到轨排的模型传递函数的实时估计参数:
Figure BDA00022691698000001115
式中,
Figure BDA00022691698000001116
表示估计得到的模型传递函数,
Figure BDA00022691698000001117
Figure BDA00022691698000001118
分别表示状态估计向量
Figure BDA00022691698000001119
的第1个、第2个、第3个、第4个、第5个分量。
进一步需要说明的是,得到磁浮列车轨排的关键参数后需要对其进行判断,以便确定轨排的状况,步骤S400即为该步骤。本发明中磁浮列车轨排检测异常情况包括如下两种:
(1)对于辨识出的轨排模型传递函数的幅频特性峰值超出设定的第一阈值,则判断轨排在该位置的柔度系数较大,存在支撑不良情况,用于提示该处轨排需要进一步进行维护保养,以图10为例,前述幅频特性峰值即为图中纵坐标的尖峰幅值;
(2)对于辨识出的轨排的模型传递函数的幅频特性峰值对应的频率小于设定的第二阈值,且幅频特性峰值对应的阻尼系数小于第三阈值,则判断该轨排存在螺栓松动、支撑刚度下降等问题,用于提示该处轨排需要进一步进行维护保养,同样,以图10为例,前述幅频特性峰值对应的频率为图中的横坐标,具体为图中尖峰幅值对应的横坐标。
判断完毕后,重返步骤S00B进行下一轮检测。
以上过程可以在磁浮列车的悬浮控制器中通过软件算法实现,估计得到的轨排实时参数既可以作为悬浮控制自适应振动控制算法的参考输入,又可以作为轨道检测的辅助手段,为轨排的检修和维护提供参考。
进一步地,为了更好地理解本发明,下面以实际的中低速磁浮列车为例,按照以上实施步骤,演示本发明对轨排的模态参数识别的效果。
假定磁浮列车的一个悬浮模块(包含左右两个悬浮单元)位于一段如图5所示的轨排1上,该轨排包括两根平行间隔设置的F轨11和设置于两根F轨之间的轨枕12。悬浮系统的主要参数如下:
悬浮电磁铁匝数:360匝;
单个线圈对应的磁极面积:0.028m×0.66m;
线圈形式:2个串联。
轨排的某阶振动模态形状如图6所示,其模态频率为70Hz,模态的阻尼比设置为0.005。则当悬浮模块位于轨排的中央位置时,悬浮系统将发生轨排自激振动,对应的两个悬浮单元的悬浮间隙波形如图7所示。假设在左侧的第一个悬浮单元中实施上述监测流程,则对应的轨排参数辨识过程如图8所示(图中以轨排模型式(9)中的a1和a2为例,其它参数的收敛过程与此类似)。从图中可以看出,轨排的模型参数在1s之内便可收敛到稳定值,且收敛的数值和真值差别很小。
为了便于对比说明本发明的辨识精度,图9所示是给定的真实轨排模型的频率响应特性曲线,图10是在辨识稳定后(2s时)利用辨识的结果绘制的轨排频率响应曲线。从二者的对比上可以看出,估计的参数(图10)和真实的参数(图9)在曲线形状上几乎一样,在轨排的模态频率(70Hz)位置均有一个明显的尖峰。通过对该尖峰的幅频特性分析,可以得出轨排的模态阻尼比约为0.005,与设定的值一致。这说明本发明提出的方法能够有效地对轨排的模态参数进行辨识。对于模态频率低、且模态阻尼系数小的轨排,提示轨排结构中可能存在螺栓松动、结构疲劳等问题,需要进一步进行检修。
图9和图10唯一的区别是二者的幅频响应曲线存在整体上的增益差别。这个原因是由式(12)中的参数η引起的,该参数可以看作是其它悬浮单元对轨排的整体幅值响应造成的影响。由于η无法估计,且η和其它悬浮单元所处的位置有关,因此获得的轨排模型在幅频增益上存在误差。但是该误差不影响对轨排特性的判断。
综上所述,本发明与现有技术相比,具有如下优点:
1)本发明完全嵌入在悬浮控制器的控制计算机中,无需单独构造轨检设备,也无需设置专用的车载机箱或外部显示系统,无需更改现有磁浮车辆的硬件配置,结构简单,易于实现;
2)本发明无需使用车载电源或外接电源,不额外消耗能量;
3)本发明提出的轨排状况检测方法收敛速度快,可以适应较宽的车速范围,应用范围广;
4)相对振动幅值检测的方法,本方法可以得到更多能够刻画轨排状况的参数,包括轨排的模态频率、阻尼系数、柔度系数等;而且本方法对于车辆行驶时的间隙波动、测量噪声不敏感,在轨排几乎不发生显著振动的情况下也能辨识出轨排参数;
5)本发明不影响磁浮列车的正常运行。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,包括如下步骤:
S100、根据磁浮列车自身配置的悬浮传感器采集电磁铁和轨排之间的悬浮间隙、电磁铁垂向加速度以及通过电流传感器采集电磁铁电流信息;
S200、利用实时估计方法估计出悬浮单元所处位置的轨排垂向位移和悬浮力;
S300、采用自适应辨识器实时估计轨排的动力学参数,并构建轨排的动力学模型传递函数,得出用于刻画轨排结构状况的关键参数,具体表现为:
S301、更新观测向量:
Figure FDA0003015479720000011
式中,
Figure FDA0003015479720000012
表示第n步的观测向量,上标T为向量的转置,
Figure FDA0003015479720000013
表示悬浮单元第n步的悬浮力、
Figure FDA0003015479720000014
表示悬浮单元第n–1步的悬浮力、
Figure FDA0003015479720000015
表示悬浮单元第n–2步的悬浮力,
Figure FDA0003015479720000016
表示轨排第n–1步估计得到的垂向位移、
Figure FDA0003015479720000017
表示轨排第n–2步估计得到的垂向位移,n表示当前迭代的步数,可取任意一个整数;
S302、更新增益向量:
Figure FDA0003015479720000018
式中,K(n)是5×1的增益向量,P(n-1)是第n-1步的误差协方差矩阵,γ是遗忘因子,取0.9<γ<1;
S303、更新状态估计向量:
Figure FDA0003015479720000019
式中,
Figure FDA00030154797200000110
为第n+1步的状态估计向量,
Figure FDA00030154797200000111
为第n步的状态估计向量,K(n+1)是5×1的增益向量,
Figure FDA00030154797200000112
表示轨排第n步估计得到的电磁铁垂向位移;
S304、更新协方差矩阵:
Figure FDA0003015479720000021
式中,P(n+1)是第n+1步的误差协方差矩阵,P(n)是第n步的误差协方差矩阵,I是5×5的单位矩阵,
Figure FDA0003015479720000022
表示第n+1步的观测向量,上标T为向量的转置;
S305,得到轨排的模型传递函数的实时估计参数:
Figure FDA0003015479720000023
式中,
Figure FDA0003015479720000024
表示估计的轨排模型传递函数,
Figure FDA0003015479720000025
Figure FDA0003015479720000026
分别表示状态估计向量
Figure FDA0003015479720000027
的第1个、第2个、第3个、第4个、第5个分量;
S400、通过所述关键参数对磁浮列车轨排进行日常监测和维护。
2.根据权利要求1所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述步骤S200中轨排垂向位移通过如下公式进行估计:
Figure FDA0003015479720000028
式中,
Figure FDA0003015479720000029
为估计得到的轨排垂向位移,x1、x2、x3和x4均为状态变量,
Figure FDA00030154797200000210
Figure FDA00030154797200000211
Figure FDA00030154797200000212
分别为状态变量x1、x2、x3和x4的一阶导数,σ1为采集的电磁铁和轨排之间的悬浮间隙数据低频特征参数,σ2为采集电磁铁和轨排之间的悬浮间隙数据高频特征参数,δ1为电磁铁和轨排之间的间隙,c1为电磁铁的垂向加速度。
3.根据权利要求2所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述步骤S200中悬浮单元的悬浮力的表达公式如下:
Figure FDA00030154797200000213
式中,i1表示悬浮电磁铁的电流,i10表示悬浮单元的稳态悬浮电流,z0是稳态悬浮间隙,kf和kz是与悬浮系统的悬浮间隙、电磁铁尺寸相关的参数,可取稳态悬浮间隙情况下的标称值,其中:
Figure FDA0003015479720000031
式中,A是电磁铁的磁极面积,N是电磁铁的线圈匝数,μ0是空气磁导率。
4.根据权利要求3所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述轨排的关键参数包括模态频率、阻尼系数和柔度系数。
5.根据权利要求4所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述步骤S100之前还包括如下步骤:
S00A、变量初始化:创建初始观测向量
Figure FDA0003015479720000032
状态估计向量
Figure FDA0003015479720000033
以及估计误差协方差矩阵P(0),其中
Figure FDA0003015479720000034
P(0)为5×5的矩阵,且矩阵P(0)的元素初始值取较大的值。
6.根据权利要求5所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述步骤S00A和步骤S100之间还包括如下步骤:
S00B、判断悬浮系统的悬浮状态:若悬浮系统处于悬浮状态,则进入步骤S100,否则,停止磁浮列车轨排监测流程。
7.根据权利要求6所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,所述磁浮列车轨排检测异常情况包括如下两种:
(1)对于辨识出的轨排模型传递函数的幅频特性峰值超出设定的第一阈值,则判断轨排在该位置的柔度系数较大,存在支撑不良情况,用于提示该位置轨排需要进一步进行维护保养;
(2)对于辨识出的轨排的模型传递函数的幅频特性峰值对应的频率小于设定的第二阈值,且幅频特性峰值对应的阻尼系数小于第三阈值,则判断该轨排存在螺栓松动、支撑刚度下降等问题,提示该处轨排需要进一步进行维护保养。
8.根据权利要求1至7中任一项所述的基于车载悬浮控制器的磁浮列车轨排监测方法,其特征在于,得到的轨排实时估计参数同时还可以作为悬浮控制自适应振动控制算法的参考输入。
CN201911098711.6A 2019-11-12 2019-11-12 基于车载悬浮控制器的磁浮列车轨排监测方法 Active CN110779743B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911098711.6A CN110779743B (zh) 2019-11-12 2019-11-12 基于车载悬浮控制器的磁浮列车轨排监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911098711.6A CN110779743B (zh) 2019-11-12 2019-11-12 基于车载悬浮控制器的磁浮列车轨排监测方法

Publications (2)

Publication Number Publication Date
CN110779743A CN110779743A (zh) 2020-02-11
CN110779743B true CN110779743B (zh) 2021-06-01

Family

ID=69390452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911098711.6A Active CN110779743B (zh) 2019-11-12 2019-11-12 基于车载悬浮控制器的磁浮列车轨排监测方法

Country Status (1)

Country Link
CN (1) CN110779743B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113525098B (zh) * 2020-04-20 2023-02-03 株洲中车时代电气股份有限公司 一种磁浮车辆的悬浮控制方法及其装置
CN113619400B (zh) * 2020-05-07 2023-05-23 株洲中车时代电气股份有限公司 一种悬浮系统的反馈线性化控制方法及其装置
CN113865691A (zh) * 2020-06-30 2021-12-31 株洲中车时代电气股份有限公司 一种磁浮列车车轨共振辨识的方法及系统
CN113212180B (zh) * 2021-06-09 2022-04-08 中车株洲电力机车有限公司 一种磁浮列车、悬浮控制系统和垂向阻尼信号的计算方法
CN113836637B (zh) * 2021-09-15 2023-10-03 四川发展磁浮科技有限公司 一种测定磁浮列车重量的方法、装置以及可读存储介质
CN113619402B (zh) * 2021-09-17 2023-04-07 中车株洲电力机车有限公司 一种磁浮列车、悬浮控制系统及提高运行平稳性的方法
CN113997795B (zh) * 2021-11-04 2024-04-26 同济大学 抑制磁浮列车悬浮系统横向冲击干扰信号处理方法和系统
CN115817185B (zh) * 2023-02-15 2023-05-05 华东交通大学 一种中低速磁浮列车速度控制方法及系统
CN116358401B (zh) * 2023-05-25 2023-08-08 宁波中车时代传感技术有限公司 一种中低速磁浮间隙传感器及其检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110411595A (zh) * 2019-08-15 2019-11-05 中国人民解放军国防科技大学 一种无温度传感器的磁浮列车电磁铁温度估计方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3704630B2 (ja) * 1996-10-24 2005-10-12 株式会社日立製作所 超電導磁石のねじり振動検出コイルを備えた磁気浮上列車
JP6184926B2 (ja) * 2014-09-30 2017-08-23 富士フイルム株式会社 椎骨セグメンテーション装置、方法およびプログラム
CN104477049B (zh) * 2014-11-28 2016-08-17 中国人民解放军国防科学技术大学 基于虚拟能量俘获器的磁浮列车-桥梁自激振动抑制方法
CN207352143U (zh) * 2017-11-08 2018-05-11 西南交通大学 一种高温超导磁悬浮的磁-力测试装置
CN109229138B (zh) * 2018-10-25 2023-05-16 湖南根轨迹智能科技有限公司 一种用于中低速磁浮列车的轨排振动检测装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110411595A (zh) * 2019-08-15 2019-11-05 中国人民解放军国防科技大学 一种无温度传感器的磁浮列车电磁铁温度估计方法

Also Published As

Publication number Publication date
CN110779743A (zh) 2020-02-11

Similar Documents

Publication Publication Date Title
CN110779743B (zh) 基于车载悬浮控制器的磁浮列车轨排监测方法
CN102607867B (zh) 基于glrt列车悬挂系在途故障检测系统及检测方法
ES2354147T3 (es) Procedimiento y dispositivo para el análisis de rodamientos en máquinas.
Esteban et al. Model-based approach for elevator performance estimation
JP5139163B2 (ja) 移動体の異常検出方法
CN112590559B (zh) 一种磁悬浮轨道平顺性偏点弦测的配置方法
JP6550654B2 (ja) 軌道制御装置
JPWO2014126177A1 (ja) 軌道制御装置
Aceituno et al. On the design of a scaled railroad vehicle for the validation of computational models
Meskell et al. Investigation of the nonlinear behaviour of damping controlled fluidelastic instability in a normal triangular tube array
Jin et al. Modal parameters identification method based on symplectic geometry model decomposition
CN113844274B (zh) 悬浮电磁铁横向位移修正系统、方法、设备及存储介质
Koch et al. Design and modeling of a quarter-vehicle test rig for active suspension control
CN113653734A (zh) 一种磁悬浮轴承定子及基座振动抑制试验系统
CN109932152A (zh) 一种汽车喇叭共振检测装置及检测方法
CN109902414A (zh) 超低频高阻尼隔振器及参数确定方法、装置和轨道
CN111103809B (zh) 一种用于高速和中低速磁悬浮列车的悬浮控制仿真平台
CN111144037B (zh) 一种轨道交通变流器柜体与车体联接刚度确定的方法
CN105035099A (zh) 引入桥梁振动速度的磁浮列车—桥梁自激振动抑制方法
CN112729734B (zh) 一种测量串联式隔振器传递特性的方法
JP2012011860A (ja) 鉄道車両の振動制御装置及び鉄道車両の振動制御方法
US10890499B2 (en) System and method for predicting strain power spectral densities of light machine structure
CN114166443B (zh) 一种仪器安装结构组件动包络干涉测试与改进方法
Datta et al. On the optimal location of sensors in multi-storeyed buildings
Diana et al. High speed railways: pantograph and overhead lines modelling and simulation

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