CN107544243A - 基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 - Google Patents
基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 Download PDFInfo
- Publication number
- CN107544243A CN107544243A CN201710609924.5A CN201710609924A CN107544243A CN 107544243 A CN107544243 A CN 107544243A CN 201710609924 A CN201710609924 A CN 201710609924A CN 107544243 A CN107544243 A CN 107544243A
- Authority
- CN
- China
- Prior art keywords
- derivation
- formula
- control
- controller
- lateral
- 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
Links
Landscapes
- Vehicle Body Suspensions (AREA)
Abstract
本发明公开了基于H∞控制的高速列车横向半主动悬挂控制系统及控制方法,包括分别设置在列车前后的两套结构相同系统构成,每套系统包括用于采集车身垂向振动加速度的加速度传感器、基于H∞的悬挂控制器、电流控制器和磁流变阻尼器;基于H∞的悬挂控制器连接加速度传感器和电流控制器;电流控制器连接磁流变阻尼器,磁流变阻尼器连接转向架;本发明可根据路面状况,计算出所需要的主动控制力,通过磁流变阻尼器输出适当的阻尼来改善列车的舒适度和行驶品质。
Description
技术领域
本发明属于高速列车悬挂领域,具体涉及一种基于H∞控制的高速列车横向半主动悬挂控制系统及其控制方法。
背景技术
最近几年高速列车取得了很大的进步,随着越来越多的列车速度和轻量化设计的应用,由轨道不规则、隧道、桥梁引起的振动,侧风将极大地影响到稳定性、安全性和乘坐质量;因此主动/半主动的铁路车辆悬挂系统已广泛应用于减少列车车身不必要的振动;而我们发现,现有的技术很难同时提高列车的稳定性和乘坐质量。由此本发明设计基于H∞控制器的半主动悬挂系统通过优化控制输出指标来同时改善稳定性和乘坐质量。
发明内容
本发明公开了一种可改善列车半主动悬挂性能的基于H∞控制的高速列车横向半主动悬挂控制系统及其控制方法。
本发明采用的技术方案是:一种基于H∞控制的高速列车横向半主动悬挂控制系统,包括分别设置在列车前后的两套结构相同系统构成,每套系统包括用于采集车身垂向振动加速度的加速度传感器、基于H∞的悬挂控制器、电流控制器和磁流变阻尼器;基于H∞的悬挂控制器连接加速度传感器和电流控制器;电流控制器连接磁流变阻尼器,磁流变阻尼器连接转向架。
一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,包括以下步骤:
A、构建列车半主动悬挂控制系统的被控对象模型,获得模型参数;
B、根据牛顿定律,建立关于车身、转向架和轮对的动力学方程,根据上述方程得到状态空间模型;
C、构建基于H∞的悬挂控制器,根据步骤A构建的模型和H∞理论计算控制器增益,根据步骤B得到的状态空间模型和控制器增益得到列车半主动悬挂的控制所需要的主动力;
D、电流控制器根据主动力控制电流大小,磁流变阻尼器根据电流大小输出相应阻尼到转向架上实现控制。
进一步的,所述车身动力学方程为:
式中:mc为车身质量,ksx为二级悬挂双倍主纵向刚度,ksy为二级悬挂双倍侧向刚度,ksz为二级悬挂双倍垂向刚度,yc为车身横向位移,hcs为车身重心到二级悬挂的垂直距离,hts为车架中心到二级悬挂的垂直距离,csx为二级悬挂双倍主纵向阻尼,csy为二级悬挂双倍侧向阻尼,csz为二级悬挂双倍垂向阻尼,yc为17-自由度车身的横向位移,为yc的一次求导,为yc的二次求导,θc为17-自由度车身的侧倾角,为θc的一次求导,为θc的二次求导,yt1、yt2为17-自由度转向架的横向位移,为yt1的一次求导,为yt2的一次求导,θt1、θt2为17-自由度转向架的侧倾角,为θt1的一次求导,为θt2的一次求导,f1、f2分别为安装在前后转向架上的致动器产生的前后作用力,Jcx为车身横摆惯性矩,Jcz为车身旋转惯性矩,为ψc的二次求导,为ψc的一次求导,ψc为17-自由度车身摆横角,ψt1、ψt2为17-自由度转向架的横摆角,ds为,为ψt1的一次求导,为ψt2的一次求导,l为车身半间距。
进一步的,所述转向架动力方程如下,式中i=1,2:
式中:mt为转向架质量,htp为车架重心到主悬挂的垂直距离,kpx为主悬挂双倍纵向刚度,kpy为主悬挂双倍侧向刚度,kpz为主悬挂双倍垂向刚度,cpx为主悬挂双倍主纵向阻尼,cpy为主悬挂双倍侧向阻尼,cpz为主悬挂双倍垂向阻尼,htp为车架重心到主悬挂的垂直距离,yw(2i-1)、yw(2i)为17-自由度轮对的横向位移,为yw(2i-1)的一次求导,为yw(2i)的一次求导,Jtx为转向架旋转惯性矩,Jtz为转向架横摆惯性矩,ψw(2i-1)、ψw(2i)为17-自由度轮对的横摆角,为ψw(2i-1)的一次求导,为ψw(2i)的一次求导,dp为主悬挂半间距,b为半轴距,为17-自由度转向架的侧倾角。
进一步的,所述轮对动力学方程为如下,式中,当i=1时,j=1,2;当i=2时,j=3,4:
式中:mw为轮对质量,ywi为17-自由度轮对的横向位移,为ywi的一次求导,为ywi的二次求导,ψtj为17-自由度转向架的横摆角,为ψtj的一次求导,ytj为17-自由度转向架的横向位移,为ytj的一次求导,θtj为17-自由度转向架的侧倾角,为θtj的一次求导,f22为轮轨侧向蠕变系数,V为车速,σ为轮轨车轮卷轴,r0为车轮滚动半径,θcli为四轮的对到不平顺横向校正,为θcli的一次求导,ψwi为17-自由度轮对的横摆角,为ψwi的一次求导,为ψwi的二次求导,Jwz为轮对横摆惯性力矩,f11为轮轨纵向蠕变系数,a为半轮对接触距离,yaj为四轮的轨道不平顺的横向对齐,为yaj一次求导,λe为有效轮锥度,Kgw为横向引力刚度,Kgψ为侧向重力刚度,由下式给出
进一步的,所述状态空间模型建立方法如下:
令q=[yc,ψc,θc,yt1,ψt1,θt1,yt2,ψt2,θt2,yw1,ψw1,yw2,ψw2,yw3,ψw3,yw4,ψw4]
d1=[ya1,ya2,ya3,ya4,θcl1,θcl2,θcl3,θcl4]T,d=[d1,d1]T,f=[f1,f2]T
则方程(1)-(8)可转换为:
式中:M为列车系统质量矩阵,C为列车系统阻尼矩阵,K为列车系统刚度矩阵,f为由执行器产生的控制力,d为路面不平的干扰,Ff为控制系数矩阵,Fd为干扰系数矩阵,为q的一次求导,为q的二次求导;
式中,yc1,yc2分别为前、后转向架上方的车身横向位移,为yc1的二次求导,为yc2的二次求导;
根据铁路车辆标准,是用于评价驾驶质量,用于评价行驶稳定性。那么控制输出zv1为:
在实际情况中,执行器、悬挂横向位移和干扰都有所限制,主要有以下几点:
1)执行器饱和现象,控制力满足|fi|≤fmax,i=1,2,其中fmax是执行器能输出最大值;
2)悬挂横向位移,相对位移yc1-yt1和yc2-yt2被用于表示悬挂横向行程;由于执行器的限制,这里存在一个正数使|yci-yti|≤smax,i=1,2,其中smax表示最大悬挂横向位移;
3)干扰限制,每个分量di满足|di|≤dmax,其中dmax表示路面不平度幅值的估计值。
由此,我们定义:
式中,u,w分别表示归一化控制输入和干扰,zv2是相对于约束的控制输出。通过定义状态分量和控制输出系统动力学方程在状态空间中可表示为:
式中:A、B1、B2、C1、D11、D12、C2、D12、D22为式(9)-(12)导出的动力学矩阵,zv为被控输出,xv为,为,。
进一步的,所述基于H∞的悬挂控制器构建过程如下:
考虑测量噪声v后,
式中:u为控制输入,w为干扰;
上式干扰w和测量噪声v为白噪声,可得
式中:E为期望,为干扰w的期望,为白噪声的期望;
H∞控制器如下:
式中:为xv的估计,L为基于H∞性能的观测器增益矩阵,Kh为控制器增益;;
式中:为P的期望;
其中,P为下列方程的解:
进一步的,所述控制器增益Kh计算方法如下:
控制器u=Khxv
状态空间模型可表示为:
状态空间模型(20)在w=0是渐进稳定的;
闭环传递函数Twz(s)满足干扰衰减条件:
||Twz(s)||∞=||(C1+D12Kh)[sI-(A+B2Kh)]-1B1+D11||∞<γ (21)
式中:s为复数变量,I为单位矩阵,γ为表示干扰衰减等级的正数;
考虑连续传递函数T(s)=C(sI-A)-1B+D
其中:||C(sI-A)-1B+D||∞<γ和A在连续时间意义Re(λi(A))<0下稳定;
式中:X为式(22)的对称正定矩阵解;
设置权重矩阵Q=diag{q1,q2I2,q3I4},式中,q1为行驶品质的权重,q2为稳定性的权重,q3为约束zv2的权重;I2为、I4为下标所表示维度的单位矩阵;
给定一个正数γ,满足下式:
则状态空间模型(20)和控制器(17)是渐进稳定的,则控制器增益的干扰衰减等级为γ;
通过(23)可计算出控制器增益。
本发明的有益效果是:
(1)本发明通过设计控制器使控制输出提高车辆的操纵性,可保证其稳定性;
(2)本发明由于线性矩阵不等式计算的方便性,可用于17自由度全尺寸车辆模型。
(3)相比于频域的H∞控制,基于时域的H∞控制器更适合该系统的控制器设计;利用这一特点,可将基于时域的H∞控制器应用该全尺寸车辆模型。
附图说明
图1为本发明中控制方法流程示意图。
图2为本发明系统结构示意图。
图3为本发明中列车全尺寸模型简图。
图4位本发明中关于轨道不平度定义的横向对齐图。
图5为本发明中关于轨道不平度定义的斜度等级图。
图6为本发明与现有技术相比稳定性、驾驶质量和平衡性幅值对比图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
一种基于H∞控制的高速列车横向半主动悬挂控制系统,包括分别设置在列车前后的两套结构相同系统构成,每套系统包括用于采集车身垂向振动加速度的加速度传感器、基于H∞的悬挂控制器、电流控制器和磁流变阻尼器;基于H∞的悬挂控制器连接加速度传感器和电流控制器;电流控制器连接磁流变阻尼器,磁流变阻尼器连接转向架;基于H∞的悬挂控制器根据列车全尺寸模型和H∞控制理论计算所需要的增益,根据加速度传感器所采集到的车身垂向振动加速度信号,计算出悬挂半主动控制所需要的主动力,并发出控制指令,具体计算方法如下详细介绍;电流控制器根据悬挂控制器的控制指令实时控制执行器所需要的电流,磁流变阻尼器在电流控制的作用下输出所需的阻尼到转向架上,提高列车舒适度和稳定性。
一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,包括以下步骤:
A、构建列车半主动悬挂控制系统的被控对象模型,获得模型参数;
B、根据牛顿定律,建立关于车身、转向架和轮对的动力学方程,根据上述方程得到状态空间模型;
所述车身动力学方程为:
式中:mc为车身质量,ksx为二级悬挂双倍主纵向刚度,ksy为二级悬挂双倍侧向刚度,ksz为二级悬挂双倍垂向刚度,yc为17-自由度车身横向位移,hcs为车身重心到二级悬挂的垂直距离,hts为车架中心到二级悬挂的垂直距离,csx为二级悬挂双倍主纵向阻尼,csy为二级悬挂双倍侧向阻尼,csz为二级悬挂双倍垂向阻尼,yc为17-自由度车身的横向位移,为yc的一次求导,为yc的二次求导,θc为17-自由度车身的侧倾角,为θc的一次求导,为θc的二次求导,yt1、yt2为17-自由度转向架的横向位移,为yt1的一次求导,为yt2的一次求导,θt1、θt2为17-自由度转向架的侧倾角,为θt1的一次求导,为θt2的一次求导,f1、f2分别为安装在前后转向架上的致动器产生的前后作用力,Jcx为车身横摆惯性矩,Jcz为车身旋转惯性矩,为ψc的二次求导,为ψc的一次求导,ψc为17-自由度车身摆横角,ψt1、ψt2为17-自由度转向架的横摆角,ds为,为ψt1的一次求导,为ψt2的一次求导,l为车身半间距。
所述转向架动力方程如下,式中i=1,2:
式中:mt为转向架质量,htp为车架重心到主悬挂的垂直距离,kpx为主悬挂双倍纵向刚度,kpy为主悬挂双倍侧向刚度,kpz为主悬挂双倍垂向刚度,cpx为主悬挂双倍主纵向阻尼,cpy为主悬挂双倍侧向阻尼,cpz为主悬挂双倍垂向阻尼,htp为车架重心到主悬挂的垂直距离,yw(2i-1)、yw(2i)为17-自由度轮对的横向位移,为yw(2i-1)的一次求导,为yw(2i)的一次求导,Jtx为转向架旋转惯性矩,Jtz为转向架横摆惯性矩,ψw(2i-1)、ψw(2i)为17-自由度轮对的横摆角,为ψw(2i-1)的一次求导,为ψw(2i)的一次求导,dp为主悬挂半间距,b为半轴距,为17-自由度转向架的侧倾角。
所述轮对动力学方程为如下,式中,当i=1时,j=1,2;当i=2时,j=3,4:
式中:mw为轮对质量,ywi为17-自由度轮对的横向位移,为ywi的一次求导,为ywi的二次求导,ψtj为17-自由度转向架的横摆角,为ψtj的一次求导,ytj为17-自由度转向架的横向位移,为ytj的一次求导,θtj为17-自由度转向架的侧倾角,为θtj的一次求导,f22为轮轨侧向蠕变系数,V为车速,σ为轮轨车轮卷轴,r0为车轮滚动半径,θcli为四轮的对到不平顺横向校正,为θcli的一次求导,ψwi为17-自由度轮对的横摆角,为ψwi的一次求导,为ψwi的二次求导,Jwz为轮对横摆惯性力矩,f11为轮轨纵向蠕变系数,a为半轮对接触距离,yaj为四轮的轨道不平顺的横向对齐,为yaj一次求导,λe为有效轮锥度,Kgw为横向引力刚度,Kgψ为侧向重力刚度,由下式给出
状态空间模型建立方法如下:
令q=[yc,ψc,θc,yt1,ψt1,θt1,yt2,ψt2,θt2,yw1,ψw1,yw2,ψw2,yw3,ψw3,yw4,ψw4]
d1=[ya1,ya2,ya3,ya4,θcl1,θcl2,θcl3,θcl4]T,d=[d1,d1]T,f=[f1,f2]T
则方程(1)-(8)可转换为:
式中:M为列车系统质量矩阵,C为列车系统阻尼矩阵,K为列车系统刚度矩阵,其中M、C、K均是17×17的矩阵;f为由执行器产生的控制力,d为路面不平的干扰,Ff为控制系数矩阵为17×2的矩阵,Fd为干扰系数矩阵为17×16的矩阵,为q的一次求导,为q的二次求导;;
为了实现该半悬挂控制系统,系统的测量输出yv为:
式中,yc1,yc2分别为前、后转向架上方的车身横向位移,为yc1的二次求导,为yc2的二次求导;
根据铁路车辆标准,是用于评价驾驶质量,用于评价行驶稳定性。那么控制输出zv1为:
在实际情况中,执行器、悬挂横向位移和干扰都有所限制,主要有以下几点:
1)执行器饱和现象,控制力满足|fi|≤fmax,i=1,2,其中fmax是执行器能输出最大值;
2)悬挂横向位移,相对位移yc1-yt1和yc2-yt2被用于表示悬挂横向行程;由于执行器的限制,这里存在一个正数使|yci-yti|≤smax,i=1,2,其中smax表示最大悬挂横向位移;
3)干扰限制,每个分量di满足|di|≤dmax,其中dmax表示路面不平度幅值的估计值。
由此,我们定义:
式中,u,w分别表示归一化控制输入和干扰,zv2是相对于约束的控制输出。
通过定义状态分量和控制输出系统动力学方程在状态空间中可表示为:
式中:A、B1、B2、C1、D11、D12、C2、D12、D22为式(9)-(12)导出的动力学矩阵,zv为被控输出,xv为状态变量,为xv的求导。
列车系统的横向振动主要由横向对齐ya和斜度等级θcl引入的,
其中yl,yr分别表示左右轮轨的不对齐程度;zl,zr分别表示左右轮轨的垂直不对齐程度。
考虑轨道不平度,由图4可得如下方程:
式中:ya1(t)、ya2(t)、ya3(t)、ya4(t)为四轮的轨道不平顺的横向对齐,V为列车速度;
令和那么
N(s)=Gwn(s)W(s) (26)
其中,N(s)和W(s)为n和w的拉普拉斯变换,Gwn(s)为其传递函数;
随机横向对准和轨道不平顺常用他们的功率密度表示:
式中:Sa(ω)为,Sc(ω)为,Λa、Λv为路面不平的标量,Ωc、Ωr、Ωs为截断波数,为Ωc的估计,为Ωr的估计,为Ωs的估计;
ω(rad/s),Ω(rad/m)分别是角频率和空间频率;L∈[Lmin,Lmax]是路径的空间波长范围;
引进周期不平度:
式中:Λvp是修正余弦波的幅值,Λap是周期不平度的幅值。
C、构建基于H∞的悬挂控制器,根据步骤A构建的模型和H∞理论计算控制器增益,根据步骤B得到的状态空间模型和控制器增益得到列车半主动悬挂的控制所需要的主动力;
根据分离原理,将线性系统的最优输出反馈控制器设计成两个独立的部分:1)一个优化观测器;2)一个优化反馈控制器;结合观测器和反馈控制器可以保持单独设计的观测器和控制器的稳定性;为了方便起见,分离原则在这里被用来设计H∞控制器。
H∞控制器构建过程如下:
首先设计观测器估计不可测的状态xv,考虑测量噪声v后,状态空间模型可转换为:
式中:u为控制输入,w为干扰;
上式干扰w和测量噪声v为白噪声,可得
式中:E为期望,为干扰w的期望,为白噪声的期望;
H∞控制器如下:
式中:为xv的估计,L为基于H∞性能的观测器增益矩阵,Kh为控制器增益;
式中:为P的期望;
其中,P为下列方程的解:
特别的,上述黎卡提方程有一个特殊解;1)(C2,A)是可测的;2) 在虚轴上没有不可测的模式。
假设状态可测,那么在控制器u=Khxv的情况下系统(14)被写为
状态空间模型(20)在w=0是渐进稳定的;
闭环传递函数Twz(s)满足干扰衰减条件:
||Twz(s)||∞=||(C1+D12Kh)[sI-(A+B2Kh)]-1B1+D11||∞<γ (21)
式中:s为复数变量,I为单位矩阵,γ为表示干扰衰减等级的正数;
为了求出增益Kh,有下面的计算过程;考虑连续传递函数T(s)=C(sI-A)-1B+D
其中:||C(sI-A)-1B+D||∞<γ和A在连续时间意义Re(λi(A))<0下稳定;
这里存在一个线性矩阵不等式(LMI)的对称正定矩阵解X:
式中:X为式(22)的对称正定矩阵解;
设置权重矩阵Q=diag{q1,q2I2,q3I4},式中,q1为行驶品质的权重,q2为稳定性的权重,q3为约束zv2的权重;I2为、I4为下标所表示维度的单位矩阵。
给定一个正数γ,满足下式:
则状态空间模型(20)和控制器(17)是渐进稳定的,则控制器增益的干扰衰减等级为γ;
通过(23)可计算出控制器增益。
计算出控制器增益后,由H∞控制器求得列车半主动悬挂的控制输入。
D、电流控制器根据主动力控制电流大小,磁流变阻尼器根据电流大小输出相应阻尼到转向架上实现控制;提高列车的舒适度和稳定新。
该半主动悬挂系统及方法根据路面状况,计算出所需要的主动控制力,通过执行器也就是磁流变阻尼器输出适当的阻尼来改善列车的舒适度与行驶品质。
仿真分析,以闭环系统的频谱分析为例,图6是从dmaxya和dmaxθcl到zv1的闭环系统幅值响应和表示稳定性、驾驶质量、平衡性的3个H∞控制系统;由图可知,本发明设计的H∞控制器对于列车的稳定性、驾驶质量、平衡性都有一定程度的提高。
本发明中半主动悬挂控制系统主要由前后两套独立的半主动悬挂系统组成,悬挂控制器根据传感器检测的垂向振动信号计算出不同的控制力,通过执行器输出不同的阻尼来调节列车不同路况下的舒适度和稳定性,使半主动悬挂始终工作在最佳状态;本发明的半主动悬挂控制系统用于在17自由度全尺寸模型,并且可提高车辆的操纵性,保证其稳定性。
Claims (8)
1.一种基于H∞控制的高速列车横向半主动悬挂控制系统,其特征在于,包括分别设置在列车前后的两套结构相同系统构成,每套系统包括用于采集车身垂向振动加速度的加速度传感器、基于H∞的悬挂控制器、电流控制器和磁流变阻尼器;基于H∞的悬挂控制器连接加速度传感器和电流控制器;电流控制器连接磁流变阻尼器,磁流变阻尼器连接转向架。
2.如权利要求1所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,包括以下步骤:
A、构建列车半主动悬挂控制系统的被控对象模型,获得模型参数;
B、根据牛顿定律,建立关于车身、转向架和轮对的动力学方程,根据上述方程得到状态空间模型;
C、构建基于H∞的悬挂控制器,根据步骤A构建的模型和H∞理论计算控制器增益,根据步骤B得到的状态空间模型和控制器增益得到列车半主动悬挂的控制所需要的主动力;
D、电流控制器根据主动力控制电流大小,磁流变阻尼器根据电流大小输出相应阻尼到转向架上实现控制。
3.根据权利要求2所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述车身动力学方程为:
式中:mc为车身质量,ksx为二级悬挂双倍主纵向刚度,ksy为二级悬挂双倍侧向刚度,ksz为二级悬挂双倍垂向刚度,yc为17-自由度车身横向位移,hcs为车身重心到二级悬挂的垂直距离,hts为车架中心到二级悬挂的垂直距离,csx为二级悬挂双倍主纵向阻尼,csy为二级悬挂双倍侧向阻尼,csz为二级悬挂双倍垂向阻尼,yc为17-自由度车身的横向位移,为yc的一次求导,为yc的二次求导,θc为17-自由度车身的侧倾角,为θc的一次求导,为θc的二次求导,yt1、yt2为17-自由度转向架的横向位移,为yt1的一次求导,为yt2的一次求导,θt1、θt2为17-自由度转向架的侧倾角,为θt1的一次求导,为θt2的一次求导,f1、f2分别为安装在前后转向架上的致动器产生的前后作用力,Jcx为车身横摆惯性矩,Jcz为车身旋转惯性矩,为ψc的二次求导,为ψc的一次求导,ψc为17-自由度车身摆横角,ψt1、ψt2为17-自由度转向架的横摆角,ds为,为ψt1的一次求导,为ψt2的一次求导,l为车身半间距。
4.根据权利要求3所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述转向架动力方程如下,式中i=1,2:
式中:mt为转向架质量,htp为车架重心到主悬挂的垂直距离,kpx为主悬挂双倍纵向刚度,kpy为主悬挂双倍侧向刚度,kpz为主悬挂双倍垂向刚度,cpx为主悬挂双倍主纵向阻尼,cpy为主悬挂双倍侧向阻尼,cpz为主悬挂双倍垂向阻尼,htp为车架重心到主悬挂的垂直距离,yw(2i-1)、yw(2i)为17-自由度轮对的横向位移,为yw(2i-1)的一次求导,为yw(2i)的一次求导,Jtx为转向架旋转惯性矩,Jtz为转向架横摆惯性矩,ψw(2i-1)、ψw(2i)为17-自由度轮对的横摆角,为ψw(2i-1)的一次求导,为ψw(2i)的一次求导,dp为主悬挂半间距,b为半轴距,为17-自由度转向架的侧倾角。
5.根据权利要求4所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述轮对动力学方程为如下,式中,当i=1时,j=1,2;当i=2时,j=3,4:
式中:mw为轮对质量,ywi为17-自由度轮对的横向位移,为ywi的一次求导,为ywi的二次求导,ψtj为17-自由度转向架的横摆角,为ψtj的一次求导,ytj为17-自由度转向架的横向位移,为ytj的一次求导,θtj为17-自由度转向架的侧倾角,为θtj的一次求导,f22为轮轨侧向蠕变系数,V为车速,σ为轮轨车轮卷轴,r0为车轮滚动半径,θcli为四轮的对到不平顺横向校正,为θcli的一次求导,ψwi为17-自由度轮对的横摆角,为ψwi的一次求导,为ψwi的二次求导,Jwz为轮对横摆惯性力矩,f11为轮轨纵向蠕变系数,a为半轮对接触距离,yaj为四轮的轨道不平顺的横向对齐,为yaj一次求导,λe为有效轮锥度,Kgw为横向引力刚度,Kgψ为侧向重力刚度,由下式给出
。
6.根据权利要求5所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述状态空间模型建立方法如下:
令q=[yc,ψc,θc,yt1,ψt1,θt1,yt2,ψt2,θt2,yw1,ψw1,yw2,ψw2,yw3,ψw3,yw4,ψw4]
d1=[ya1,ya2,ya3,ya4,θcl1,θcl2,θcl3,θcl4]T,d=[d1,d1]T,f=[f1,f2]T
则方程(1)-(8)可转换为:
式中:M为列车系统质量矩阵,C为列车系统阻尼矩阵,K为列车系统刚度矩阵,f为由执行器产生的控制力,d为路面不平的干扰,Ff为控制系数矩阵,Fd为干扰系数矩阵,为q的一次求导,为q的二次求导;
系统的测量输出yv为:
式中,yc1,yc2分别为前、后转向架上方的车身横向位移,为yc1的二次求导,为yc2的二次求导;
平顺性指标的控制输出zv1为:
由上述可得,系统动力学方程在状态空间中可表示为:
式中:A、B1、B2、C1、D11、D12、C2、D12、D22为式(9)-(12)导出的动力学矩阵,zv为被控输出,xv为 为xv的一次求导。
7.根据权利要求6所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述基于H∞悬挂控制器构建过程如下:
考虑测量噪声v后,状态空间模型变为:
式中:u为控制输入,w为干扰;
上式干扰w和测量噪声v为白噪声,可得
式中:E为期望,为干扰w的期望,为白噪声的期望;
H∞控制器如下:
式中:为xv的估计,L为基于H∞性能的观测器增益矩阵,Kh为控制器增益;
式中:为P的期望;
其中,P为下列方程的解:
8.根据权利要求7所述的一种基于H∞控制的高速列车横向半主动悬挂控制系统的控制方法,其特征在于,所述控制器增益Kh计算方法如下:
控制器u=Khxv
状态空间模型可表示为:
状态空间模型(20)在w=0是渐进稳定的;
闭环传递函数Twz(s)满足干扰衰减条件:
||Twz(s)||∞=||(C1+D12Kh)[sI-(A+B2Kh)]-1B1+D11||∞<γ (21)
式中:s为复变量,I为单位矩阵,γ为表示干扰衰减等级的正数;
考虑连续传递函数T(s)=C(sI-A)-1B+D
其中:||C(sI-A)-1B+D||∞<γ和A在连续时间意义Re(λi(A))<0下稳定;
式中:X为式(22)的对称正定矩阵解;
设置权重矩阵Q=diag{q1,q2I2,q3I4},式中,q1为行驶品质的权重,q2为稳定性的权重,q3为约束zv2的权重;I2为、I4为下标所表示维度的单位矩阵;
给定一个正数γ,满足下式:
则状态空间模型(20)和控制器(17)是渐进稳定的,则控制器增益的干扰衰减等级为γ;
通过(23)可计算出控制器增益。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710609924.5A CN107544243A (zh) | 2017-07-25 | 2017-07-25 | 基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710609924.5A CN107544243A (zh) | 2017-07-25 | 2017-07-25 | 基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107544243A true CN107544243A (zh) | 2018-01-05 |
Family
ID=60970977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710609924.5A Pending CN107544243A (zh) | 2017-07-25 | 2017-07-25 | 基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107544243A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110155103A (zh) * | 2019-05-27 | 2019-08-23 | 西南交通大学 | 一种基于lqg控制的高速列车半主动悬挂控制系统及方法 |
CN111123706A (zh) * | 2019-12-26 | 2020-05-08 | 湖南工业大学 | 高速列车半主动悬挂系统控制方法 |
CN113110186A (zh) * | 2021-04-19 | 2021-07-13 | 华东交通大学 | 可随路面激励调整的汽车悬架系统用磁流变阻尼器控制器 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003320931A (ja) * | 2002-05-07 | 2003-11-11 | Kayaba Ind Co Ltd | 鉄道車両振動抑制装置 |
CN102518737A (zh) * | 2011-12-28 | 2012-06-27 | 哈尔滨工业大学 | 一体化智能磁流变阻尼器及其控制方法 |
CN102616245A (zh) * | 2012-03-23 | 2012-08-01 | 株洲南车时代电气股份有限公司 | 列车转向架半主动减振器控制装置和系统及其控制方法 |
CN102797202A (zh) * | 2012-08-29 | 2012-11-28 | 北京交通大学 | 基于观测器的轨道横向不平顺检测方法 |
CN105015572A (zh) * | 2015-07-07 | 2015-11-04 | 北京交通大学 | 基于磁流液变阻尼器的列车半主动控制方法 |
CN106647256A (zh) * | 2016-10-08 | 2017-05-10 | 西南交通大学 | 一种基于h∞ pid的主动悬架控制系统及控制方法 |
-
2017
- 2017-07-25 CN CN201710609924.5A patent/CN107544243A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003320931A (ja) * | 2002-05-07 | 2003-11-11 | Kayaba Ind Co Ltd | 鉄道車両振動抑制装置 |
CN102518737A (zh) * | 2011-12-28 | 2012-06-27 | 哈尔滨工业大学 | 一体化智能磁流变阻尼器及其控制方法 |
CN102616245A (zh) * | 2012-03-23 | 2012-08-01 | 株洲南车时代电气股份有限公司 | 列车转向架半主动减振器控制装置和系统及其控制方法 |
CN102797202A (zh) * | 2012-08-29 | 2012-11-28 | 北京交通大学 | 基于观测器的轨道横向不平顺检测方法 |
CN105015572A (zh) * | 2015-07-07 | 2015-11-04 | 北京交通大学 | 基于磁流液变阻尼器的列车半主动控制方法 |
CN106647256A (zh) * | 2016-10-08 | 2017-05-10 | 西南交通大学 | 一种基于h∞ pid的主动悬架控制系统及控制方法 |
Non-Patent Citations (1)
Title |
---|
宗路航: "磁流变阻尼器的动力学模型及其在车辆悬架中的应用研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110155103A (zh) * | 2019-05-27 | 2019-08-23 | 西南交通大学 | 一种基于lqg控制的高速列车半主动悬挂控制系统及方法 |
CN110155103B (zh) * | 2019-05-27 | 2020-10-23 | 西南交通大学 | 一种基于lqg控制的高速列车半主动悬挂控制系统及方法 |
CN111123706A (zh) * | 2019-12-26 | 2020-05-08 | 湖南工业大学 | 高速列车半主动悬挂系统控制方法 |
CN111123706B (zh) * | 2019-12-26 | 2022-05-27 | 湖南工业大学 | 高速列车半主动悬挂系统控制方法 |
CN113110186A (zh) * | 2021-04-19 | 2021-07-13 | 华东交通大学 | 可随路面激励调整的汽车悬架系统用磁流变阻尼器控制器 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA3112415C (en) | Inertial regulation active suspension system based on vehicle posture deviation, and control method thereof | |
CN110155103B (zh) | 一种基于lqg控制的高速列车半主动悬挂控制系统及方法 | |
CN106647256B (zh) | 一种基于h∞ pid的主动悬架控制系统及控制方法 | |
CN107539332B (zh) | 基于共振控制的高速列车横向半主动悬挂控制系统及控制方法 | |
CN107544243A (zh) | 基于h∞控制的高速列车横向半主动悬挂控制系统及控制方法 | |
CN105824238B (zh) | 一种高速动车组车体弹性模态振动控制方法 | |
CN105117556B (zh) | 高铁一系和二系及端部减振器阻尼系数的协同优化方法 | |
IT202100015170A1 (it) | Metodo di controllo degli ammortizzatori attivi di un veicolo stradale che prevede l'abbassamento del baricentro | |
CN108681257B (zh) | 一种主动防侧倾系统的控制器的设计方法 | |
CN102975587A (zh) | 基于双可控阻尼器的车辆半主动悬架及其控制方法 | |
CN105138783B (zh) | 高速轨道车辆车体端部横向减振器阻尼系数的设计方法 | |
CN107839425B (zh) | 一种履带车辆垂直振动与俯仰振动协同控制方法 | |
CN110032807A (zh) | 一种高速列车车体下部悬吊设备的垂向悬挂频率设计方法 | |
CN110443003B (zh) | 一种主动稳定杆系统的控制及优化设计方法 | |
CN105573116B (zh) | 一种基于自抗扰控制的坦克底盘主动悬挂系统的控制方法 | |
CN114633597A (zh) | 一种整车协调控制方法及系统 | |
CN112765728B (zh) | 单轴转向架跨座式单轨车辆的摆振优化控制方法 | |
Li et al. | Driver steering control and full vehicle dynamics study based on a nonlinear three-directional coupled heavy-duty vehicle model | |
Hirao et al. | A semi-active suspension system using ride control based on bi-linear optimal control theory and handling control considering roll feeling | |
Nguyen et al. | A motion-scheduled LPV control of full car vertical dynamics | |
Trikande et al. | Evaluation of semi-active suspension control strategies for 8x8 armoured vehicle using stochastic road profile inputs | |
CN115195799B (zh) | 一种车辆蛇行运动稳定性在线识别与主动控制仿真方法及系统 | |
WO2023108699A1 (zh) | 一种新型lqg控制方法 | |
CN109109598A (zh) | 一种基于逆模型pid的主动悬架控制系统及方法 | |
CN115570924A (zh) | 一种基于迭代学习改进的h∞控制方法 |
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: 20180105 |