CN109334377B - 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 - Google Patents
磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 Download PDFInfo
- Publication number
- CN109334377B CN109334377B CN201811058385.1A CN201811058385A CN109334377B CN 109334377 B CN109334377 B CN 109334377B CN 201811058385 A CN201811058385 A CN 201811058385A CN 109334377 B CN109334377 B CN 109334377B
- Authority
- CN
- China
- Prior art keywords
- magneto
- semi
- controller
- rheological
- time
- 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
Links
- 239000000725 suspension Substances 0.000 title claims abstract description 65
- 238000010276 construction Methods 0.000 title claims abstract description 6
- 239000006096 absorbing agent Substances 0.000 claims abstract description 47
- 230000035939 shock Effects 0.000 claims abstract description 47
- 230000001133 acceleration Effects 0.000 claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims description 33
- 238000006073 displacement reaction Methods 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 24
- 239000002131 composite material Substances 0.000 claims description 20
- 238000013016 damping Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000004069 differentiation Effects 0.000 claims description 3
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 3
- 230000032683 aging Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000005294 ferromagnetic effect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60G—VEHICLE SUSPENSION ARRANGEMENTS
- B60G17/00—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load
- B60G17/015—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements
- B60G17/018—Resilient suspensions having means for adjusting the spring or vibration-damper characteristics, for regulating the distance between a supporting surface and a sprung part of vehicle or for locking suspension during use to meet varying vehicular or surface conditions, e.g. due to speed or load the regulating means comprising electric or electronic elements characterised by the use of a specific signal treatment or control method
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60G—VEHICLE SUSPENSION ARRANGEMENTS
- B60G2600/00—Indexing codes relating to particular elements, systems or processes used on suspension systems or suspension control systems
- B60G2600/18—Automatic control means
- B60G2600/187—Digital Controller Details and Signal Treatment
- B60G2600/1871—Optimal control; Kalman Filters
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
本发明公开了一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***及其构造方法,通过簧载质量加速度传感器和车轮质量加速度传感器的输出端连接扩展卡尔曼滤波器,扩展卡尔曼滤波器的连接等效替换泰勒级数LQG时滞补偿控制器,等效替换泰勒级数LQG时滞补偿控制器分别连接时滞补偿预测控制器与半主动力求取控制器,时滞补偿预测控制器连接等效替换泰勒级数LQG时滞补偿控制器的输入,半主动力求取控制器分别连接扩展卡尔曼滤波器与控制电流求解控制器,将控制电流求解控制器连接数控电流源,进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制,通过等效替换泰勒级数LQG时滞补偿控制,可以实现精确跟踪悬架综合性能指标,进而能取得更理想的时滞补偿控制效果。
Description
技术领域
本发明属于车辆悬架控制领域,尤其涉及一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***及其构造方法。
背景技术
悬架是汽车重要的结构与功能部件,对汽车的乘坐舒适性和行驶安全性有重要的影响。磁流变半主动悬架不需要输出主动力,能根据车辆行驶工况的变化输出具有正阻尼特性的半主动控制力,无需为磁流变减振器专门配备动力源,接近主动悬架的综合性能。
但目前磁流变半主动悬架存在的问题是响应时滞,该响应时滞产生原因主要有:传感器传输时滞、控制器计算时滞、电流响应时滞、磁流变液响应时滞等影响车辆行驶的平顺性和安全性。因此,在磁流变半主动悬架的控制***设计时,必须要考虑磁流变减振器的响应时滞问题。
目前,关于磁流变减振器的响应时滞补偿,采用Smith预测器方法、基于Lyapunov-Krasovskii稳定性定理的鲁棒控制方法。虽然关于磁流变减振器的力学特性时变与响应时滞补偿控制研究现在取得了一定的成果,但与理想半主动悬架的性能还有较大的差距。
发明内容
本发明根据现有技术的不足与缺陷,提出了一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***及其构造方法,目的在于获得理想的时滞补偿控制效果。
磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***,包括在垂直方向上,车辆的车轮质量与等效成弹簧的轮胎组成车轮,车轮位于簧载质量的下方,簧载质量与车轮质量之间并联有悬架弹簧和磁流变减振器;在簧载质量上设有簧载质量加速度传感器,在车轮质量上设有车轮质量加速度传感器,所述簧载质量加速度传感器和车轮质量加速度传感器分别连接磁流变半主动悬架控制器,所述磁流变半主动悬架控制器连接数控电流源,所述数控电流源连接磁流变减振器。
进一步,所述磁流变半主动悬架控制器包括扩展卡尔曼滤波器、时滞补偿预测控制器、等效替换泰勒级数LQG时滞补偿控制器、理想半主动力求解控制器与控制电流求解控制器;所述簧载质量加速度传感器和车轮质量加速度传感器的输出端连接扩展卡尔曼滤波器,将所接收的当前时刻的簧载质量加速度与当前时刻的车轮质量加速度传输给扩展卡尔曼滤波器,所述扩展卡尔曼滤波器的输出端连接等效替换泰勒级数LQG时滞补偿控制器,等效替换泰勒级数LQG时滞补偿控制器的输出分别连接时滞补偿预测控制器与理想半主动力求解控制器,所述时滞补偿预测控制器的输出连接等效替换泰勒级数LQG时滞补偿控制器的输入,所述理想半主动力求解控制器的输出端分别连接扩展卡尔曼滤波器与控制电流求解控制器,所述控制电流求解控制器的输出端连接数控电流源。
磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***构造方法,包括以下步骤:
步骤2,等效替换泰勒级数LQG时滞补偿控制器对当前时刻的复合输出力变化量的观测量和上一工作循环理想时滞主动力Fb(t-τ)进行处理,输出给上一工作循环的理想主动控制预测力以及当前时刻的理想主动控制预测力Fp;
步骤4,理想半主动力求解控制器输出当前时刻的理想半主动控制力Fps;
步骤5,控制电流求解控制器根据当前时刻的理想半主动控制力Fps得到理想控制电流Ii;
步骤6,将理想控制电流值Ii输入至数控电流源,数控电流源输出实际控制电流Ia至磁流变减振器,进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制。
步骤1.1,针对车辆垂向方向上的运动进行动力学分析,得到磁流变半主动悬架运动状态方程;
其中,m1是车轮质量,m2是簧载质量,k1是轮胎刚度,k2是悬架刚度,z1是车轮质量的垂直位移,z2是簧载质量的垂直位移,分别为z1、z2的一阶微分,分别为z1、z2的二阶微分,q是路面不平度对车轮的位移输入,t是当前时刻,τ为响应时滞时间,c0是磁流变减振器的基值黏性阻尼系数,FMR-I是为稳态磁流变减振器通过电流控制的阻尼力部分,简称稳态控制力,磁流变减振器复合输出力时变量Fd;
步骤1.2,通过扩展卡尔曼滤波器观测出磁流变减振器复合输出力时变量Fd:
根据:Fd=FMR_r-FMR+ks+{fy0tanh[b1v+b2sgn(a)]+f0}
将Fd写成干扰形式:
其中,FMR_r为磁流变减振器的实时输出力,FMR为磁流变减振器的稳态输出力,k为磁流变减振器的稳态刚度系数,s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度,b1、b2、fy0为拟合参数,tanh()表示双曲正切函数,f0为稳态偏置力,qd为干扰量,λ为正实数,是Fd的一阶导数;
其中,为状态方程的输出,为的一阶微分,是X1的观测值,X=(x1,x2,x3,x4)T,x1=z1-q,x2=z2-z1,z1是车轮质量的垂直位移,z2是簧载质量的垂直位移,分别为z1、z2的一阶微分 λ为正实数,m1与m2分别是车轮质量和簧载质量,k1与k2分别是轮胎刚度和悬架刚度,c0为磁流变减振器的基值黏性阻尼系数,U1=[F(t-τ)], G=[-1 0 0 0]T; 为X的一阶微分;G=[-1 0 0 0]T,U=[FMR-I(t-τ)+Fd],q是路面不平度对车轮的位移输入,其一阶微分用下式表示为:
进一步,Pk、Nk、Rk通过以下方法求得:
Pk是以下黎卡提方程的唯一解;
进一步,获得当前时刻的理想主动控制预测力Fp的方法是:
Kp为控制计算向量,Kp按LQR函数(Kp,Sp,Ep)=LQR(A2,B2,Q2s,R2s,N2)求取,Sp和Ep分别是求解Kp用黎卡提方程的和特征向量, N2=[01×6],R2s=ψ[0 0 0 0 A(4,5) 0]T[0 0 0 0 A(4,5) 0], ψ为小于1的正等效替换系数,必须满足正定,则是悬架综合性能指标的等效替换,
进一步,求得的上一工作循环理想时滞主动力Fb(t-τ),具体过程为:
进一步,当前时刻的理想半主动控制力Fps的方法为:
进一步,得到理想控制电流Ii的方法为:
其中,fMR_I为稳态控制力的励磁电流系数,fMR_I=c1v+fy1tanh[b1v+b2sgn(a)],c1为受电流影响的屈服后的稳态黏性阻尼系数,v、a分别减振器活塞相对于减振器筒的速度、加速度,b1、b2为拟合参数,tanh()表示双曲正切函数。
进一步,得到磁流变减振器的稳态输出力FMR的方法为:
FMR_I=fMR_I×Ia
fMR_I=c1v+fy1tanh[b1v+b2sgn(a)]
FMR=FMR_I+FMR_0
FMR_0=ks+fy0tanh[b1v+b2sgn(a)]+f0+c0v
式中:fMR_I为稳态控制力的励磁电流系数;FMR_0为不受电流影响的稳态输出力;c0为磁流变减振器的稳态基值黏性阻尼系数;c1为受电流影响的屈服后的稳态黏性阻尼系数;k为磁流变减振器的稳态刚度系数;tanh()表示双曲正切函数,f0为稳态偏置力;s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度;b1、b2、fy0、fy1为拟合参数。
本发明的有益效果:
1、不需要进行受温度变化、部件老化、磁流变液中铁磁颗粒沉积影响的磁流变减振器力学特性准确建模,控制简单精准。
2、通过扩展卡尔曼滤波器可以直接将磁流变减振器输出力的时变量实时观测出来,为磁流变半主动悬架的时滞补偿控制奠定基础;通过等效替换泰勒级数LQG时滞补偿控制,可以实现精确跟踪悬架综合性能指标,进而能取得更理想的时滞补偿控制效果
附图说明
图1是本发明控制***与磁流变半主动悬架的连接结构示意图;
图2是图1中所示的本发明控制***的控制结构框图;
图中,1、悬架弹簧;2、簧载质量;3、簧载质量加速度传感器;4、磁流变半主动悬架控制器;5、车轮质量加速度传感器;6、磁流变减振器;7、车轮质量;8、等效成弹簧的轮胎;9、数控电流源;10、扩展卡尔曼滤波器;11、时滞补偿预测控制器;12、等效替换泰勒级数LQG时滞补偿控制器;13、理想半主动力求解控制器;14、控制电流求解控制器。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用于解释本发明,并不用于限定本发明。
如图1所示,采用1/4车2自由度车辆的磁流变半主动悬架为例,一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***,包括在垂直方向上,车辆的车轮质量7与等效成弹簧的轮胎8组成车轮,车轮位于簧载质量2的下方,簧载质量2与车轮质量7之间并联有悬架弹簧1和磁流变减振器6;在簧载质量2上设有簧载质量加速度传感器3,在车轮质量7上设有车轮质量加速度传感器5,所述簧载质量加速度传感器3和车轮质量加速度传感器5分别连接磁流变半主动悬架控制器4,所述磁流变半主动悬架控制器4连接数控电流源9,所述数控电流源9连接磁流变减振器6;
如图2所示,磁流变半主动悬架控制器4包括扩展卡尔曼滤波器10、时滞补偿预测控制器11、等效替换泰勒级数LQG时滞补偿控制器12、半主动力求取控制器13与控制电流求解控制器14;簧载质量加速度传感器3和车轮质量加速度传感器5的输出端连接扩展卡尔曼滤波器10,将所接收的当前时刻的簧载质量加速度与当前时刻的车轮质量加速度传输给扩展卡尔曼滤波器10,同时,半主动力求取控制器13将上一个工作控制循环的上一工作循环理想半主动力给扩展卡尔曼滤波器10,扩展卡尔曼滤波器10进行处理后输出当前时刻的车辆运动状态及磁流变减振器的复合输出力变化量的观测量扩展卡尔曼滤波器10的输出端连接等效替换泰勒级数LQG时滞补偿控制器12,将当前时刻复合输出力变化量的观测量传送给等效替换泰勒级数LQG时滞补偿控制器12,等效替换泰勒级数LQG时滞补偿控制器12的输出分别连接时滞补偿预测控制器11与半主动力求取控制器13,等效替换泰勒级数LQG时滞补偿控制器12将处理后得到的将输出的上一个工作循环的理想主动控制预测力传送给时滞补偿预测控制器11,将当前时刻的理想主动控制预测力Fp传送给理想半主动力求解控制器13;所述时滞补偿预测控制器11的输出连接等效替换泰勒级数LQG时滞补偿控制器12的输入,时滞补偿预测控制器11将输出的上一个工作循环求得的上一工作循环理想时滞主动力Fb(t-τ)传输给等效替换泰勒级数LQG时滞补偿控制器12;半主动力求取控制器13的输出端分别连接扩展卡尔曼滤波器10与控制电流求解控制器14,。理想半主动力求解控制器13以当前时刻的理想主动控制预测力Fp为输入,输出上一个工作循环求得的理想半主动力和当前时刻的理想半主动控制力Fps,将当前时刻的理想半主动控制力Fps控制电流求解控制器14;控制电流求解控制器14的输出端串接数控电流源9,控制电流求解控制器14以当前时刻的理想半主动控制力Fps为输入,输出当前时刻的理想控制电流值Ii给数控电流源9。数控电流源9以当前时刻的理想控制电流值Ii为输入,输出当前时刻的实际控制电流Ia至磁流变减振器6,进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制。
一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***构造方法,
步骤1.1,针对车辆垂向方向上的运动进行动力学分析,得到磁流变半主动悬架运动状态方程;
其中,m1是车轮质量,m2是簧载质量,k1是轮胎刚度,k2是悬架刚度,z1是车轮质量的垂直位移,z2是簧载质量的垂直位移,分别为z1、z2的一阶微分,分别为z1、z2的二阶微分,q是路面不平度对车轮的位移输入,t是当前时刻,τ为响应时滞时间,c0是磁流变减振器的基值黏性阻尼系数,FMR-I是为稳态磁流变减振器通过电流控制的阻尼力部分,简称稳态控制力,磁流变减振器复合输出力时变量Fd;
步骤1.2,通过扩展卡尔曼滤波器10观测出磁流变减振器复合输出力时变量Fd:
根据公式:
Fd=FMR_r-FMR+ks+{fy0tanh[b1v+b2sgn(a)]+f0} (2)
将Fd写成干扰形式:
其中,FMR_r为磁流变减振器的实时输出力,FMR为磁流变减振器的稳态输出力,k为磁流变减振器的稳态刚度系数,s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度,b1、b2、fy0为拟合参数,tanh()表示双曲正切函数,f0为稳态偏置力,qd为干扰量,λ为正实数,是Fd的一阶导数;
其中,为状态方程的输出,为的一阶微分,是X1的观测值,当前时刻的车辆运动状态向量X=(x1,x2,x3,x4)T,x1=z1-q,x2=z2-z1,向量向量向量向量λ为正实数,m1与m2分别是车轮质量和簧载质量,k1与k2分别是轮胎刚度和悬架刚度,c0为磁流变减振器的基值黏性阻尼系数,向量U1=[F(t-τ)], G=[-1 0 0 0]T; 为X的一阶微分;G=[-1 0 0 0]T,U=[FMR-I(t-τ)+Fd],q是路面不平度对车轮的位移输入,其一阶微分用下式表示为:
进一步,Pk、Nk、Rk通过以下方法求得:
Pk是以下黎卡提方程的唯一解;
步骤2,等效替换泰勒级数LQG时滞补偿控制器12对当前时刻的复合输出力变化量的观测量和上一工作循环理想时滞主动力Fb(t-τ)进行处理,输出给上一工作循环的理想主动控制预测力以及当前时刻的理想主动控制预测力Fp的方法是:
Kp按LQR函数(Kp,Sp,Ep)=LQR(A2,B2,Q2s,R2s,N2)求取,Sp和Ep分别是求解Kp用黎卡提方程的和特征向量, N2=[01×6],R2s=ψ[0 0 0 0 A(4,5) 0]T[0 0 0 0 A(4,5) 0], ψ为小于1的正等效替换系数,必须满足正定,则是悬架综合性能指标的等效替换,
因此,本发明提供的等效替换泰勒级数LQG时滞补偿控制可以解决因R2=[0]而导致LQR()函数不能工作,从而实现对悬架综合性能指标进行精确的时滞补偿控制。
步骤4,理想半主动力求解控制器13输出当前时刻的理想半主动控制力Fps的方法为:
步骤5,控制电流求解控制器14根据当前时刻的理想半主动控制力得到理想控制电流Ii的方法为:
其中,fMR_I为稳态控制力的励磁电流系数,fMR_I=c1v+fy1tanh[b1v+b2sgn(a)],c1为受电流影响的屈服后的稳态黏性阻尼系数,v、a分别减振器活塞相对于减振器筒的速度、加速度,b1、b2为拟合参数,tanh()表示双曲正切函数。
步骤6,将理想控制电流值Ii输入至数控电流源9,数控电流源9输出实际控制电流Ia至磁流变减振器6,进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制;
得到磁流变减振器6的稳态输出力FMR的方法为:
FMR_I=fMR_I×Ia (11)
fMR_I=c1v+fy1tanh[b1v+b2sgn(a)] (12)
FMR=FMR_I+FMR_0 (13)
FMR_0=ks+fy0tanh[b1v+b2sgn(a)]+f0+c0v (14)
式中:fMR_I为稳态控制力的励磁电流系数;FMR_0为不受电流影响的稳态输出力;c0为磁流变减振器的稳态基值黏性阻尼系数;c1为受电流影响的屈服后的稳态黏性阻尼系数;k为磁流变减振器的稳态刚度系数;tanh()表示双曲正切函数,f0为稳态偏置力;s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度;b1、b2、fy0、fy1为拟合参数。
将理想控制电流值Ii输入至数控电流源9,数控电流源9输出实际控制电流Ia至磁流变减振器6,进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制。
本发明中的扩展卡尔曼滤波器10能观测磁流变减振器复合输出力时变量,等效替换泰勒级数LQG时滞补偿控制器12进行磁流变半主动悬架力学特性时变与响应时滞的补偿控制,利用扩展卡尔曼滤波器10与等效替换泰勒级数LQG时滞补偿控制器12能有效改善车辆磁流变半主动悬架存在的减振器力学特性时变与响应时滞问题,实现车辆磁流变半主动悬架控制***的减振器力学特性时变与响应时滞补偿控制。
以上实施例仅用于说明本发明的设计思想和特点,其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施,本发明的保护范围不限于上述实施例。所以,凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰,均在本发明的保护范围之内。
Claims (9)
1.磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***,其特征在于,包括在垂直方向上,车辆的车轮质量(7)与等效成弹簧的轮胎(8)组成车轮,车轮位于簧载质量(2)的下方,簧载质量(2)与车轮质量(7)之间并联有悬架弹簧(1)和磁流变减振器(6);在簧载质量(2)上设有簧载质量加速度传感器(3),在车轮质量(7)上设有车轮质量加速度传感器(5),所述簧载质量加速度传感器(3)和车轮质量加速度传感器(5)分别连接磁流变半主动悬架控制器(4),所述磁流变半主动悬架控制器(4)连接数控电流源(9),所述数控电流源(9)连接磁流变减振器(6);
所述磁流变半主动悬架控制器(4)包括扩展卡尔曼滤波器(10)、时滞补偿预测控制器(11)、等效替换泰勒级数LQG时滞补偿控制器(12)、半主动力求取控制器(13)与控制电流求解控制器(14);所述簧载质量加速度传感器(3)和车轮质量加速度传感器(5)的输出端连接扩展卡尔曼滤波器(10),所述扩展卡尔曼滤波器(10)的输出端连接等效替换泰勒级数LQG时滞补偿控制器(12),等效替换泰勒级数LQG时滞补偿控制器(12)的输出分别连接时滞补偿预测控制器(11)与半主动力求取控制器(13),所述时滞补偿预测控制器(11)的输出连接等效替换泰勒级数LQG时滞补偿控制器(12)的输入,所述半主动力求取控制器(13)的输出端分别连接扩展卡尔曼滤波器(10)与控制电流求解控制器(14),所述控制电流求解控制器(14)的输出端连接数控电流源(9)。
2.一种磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***的构造方法,其特征在于,应用于如权利要求1所述的磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***;包括以下步骤:步骤1,扩展卡尔曼滤波器(10)对所接收的当前时刻两个加速度和上一个工作循环求得的理想半主动力进行处理,输出当前时刻的车辆运动状态向量X及磁流变减振器的复合输出力变化量Fd的观测量
步骤2,等效替换泰勒级数LQG时滞补偿控制器(12)对当前时刻的复合输出力变化量的观测量和上一工作循环理想时滞主动力Fb(t-τ)进行处理,输出给上一工作循环的理想主动控制预测力以及当前时刻的理想主动控制预测力Fp;
步骤4,理想半主动力求解控制器(13)输出当前时刻的理想半主动控制力Fps;
步骤5,控制电流求解控制器(14)根据当前时刻的理想半主动控制力Fps得到理想控制电流Ii;
步骤6,将理想控制电流值Ii输入至数控电流源(9),数控电流源(9)输出实际控制电流Ia至磁流变减振器(6),进而实现磁流变半主动悬架的减振器响应时滞与力学特性时变的补偿控制。
3.根据权利要求2所述的磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***的构造方法,其特征在于,步骤1中获得当前时刻的车辆运动状态向量X及磁流变减振器的复合输出力变化量Fd的观测量具体过程如下:
步骤1.1,针对车辆垂向方向上的运动进行动力学分析,得到磁流变半主动悬架运动状态方程;
其中,m1是车轮质量,m2是簧载质量,k1是轮胎刚度,k2是悬架刚度,z1是车轮质量的垂直位移,z2是簧载质量的垂直位移,分别为z1、z2的一阶微分,分别为z1、z2的二阶微分,q是路面不平度对车轮的位移输入,t是当前时刻,τ为响应时滞时间,c0是磁流变减振器的基值黏性阻尼系数,FMR-I是为稳态磁流变减振器通过电流控制的阻尼力部分,简称稳态控制力,磁流变减振器复合输出力时变量Fd;
步骤1.2,通过扩展卡尔曼滤波器10观测出磁流变减振器复合输出力时变量Fd:
根据:Fd=FMR_r-FMR+ks+{fy0tanh[b1v+b2sgn(a)]+f0}
将Fd写成干扰形式:
其中,FMR_r为磁流变减振器的实时输出力,FMR为磁流变减振器的稳态输出力,k为磁流变减振器的稳态刚度系数,s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度,b1、b2、fy0为拟合参数,tanh()表示双曲正切函数,f0为稳态偏置力,qd为干扰量,λ为正实数,是Fd的一阶导数;
其中,为状态方程的输出,为的一阶微分,是X1的观测值,矩阵X当前时刻的车辆运动状态向量,X=(x1,x2,x3,x4)T,x1=z1-q,x2=z2-z1,z1是车轮质量的垂直位移,z2是簧载质量的垂直位移,分别为z1、z2的一阶微分λ为正实数,m1与m2分别是车轮质量和簧载质量,k1与k2分别是轮胎刚度和悬架刚度,c0为磁流变减振器的基值黏性阻尼系数,U1=[F(t-τ)], G=[-1 0 0 0]T; 为X的一阶微分;G=[-1 0 0 0]T,U=[FMR-I(t-τ)+Fd],q是路面不平度对车轮的位移输入,为q的一阶微分;Pk是以下黎卡提方程的唯一解,
式中,V1是簧载质量加速度与车轮质量加速度测量信号噪声向量;E{W1 TW1}为干扰量W1的自谱;E{V1 TV1}为测量信号噪声向量V1的自谱;E{W1 TV1}为W1与V1的互谱。
9.根据权利要求2所述的磁流变半主动悬架等效替换泰勒级数LQG时滞补偿控制***的构造方法,其特征在于,得到磁流变减振器(6)的稳态输出力FMR的方法为:
FMR_I=fMR_I×Ia
fMR_I=c1v+fy1tanh[b1v+b2sgn(a)]
FMR=FMR_I+FMR_0
FMR_0=ks+fy0tanh[b1v+b2sgn(a)]+f0+c0v
式中:fMR_I为稳态控制力的励磁电流系数;FMR_0为不受电流影响的稳态输出力;c0为磁流变减振器的稳态基值黏性阻尼系数;c1为受电流影响的屈服后的稳态黏性阻尼系数;k为磁流变减振器的稳态刚度系数;tanh()表示双曲正切函数,f0为稳态偏置力;s、v、a分别减振器活塞相对于减振器筒的位移、速度、加速度;b1、b2、fy0、fy1为拟合参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811058385.1A CN109334377B (zh) | 2018-09-11 | 2018-09-11 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811058385.1A CN109334377B (zh) | 2018-09-11 | 2018-09-11 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109334377A CN109334377A (zh) | 2019-02-15 |
CN109334377B true CN109334377B (zh) | 2021-07-20 |
Family
ID=65304997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811058385.1A Active CN109334377B (zh) | 2018-09-11 | 2018-09-11 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109334377B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109334377B (zh) * | 2018-09-11 | 2021-07-20 | 江苏大学 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
CN110155103B (zh) * | 2019-05-27 | 2020-10-23 | 西南交通大学 | 一种基于lqg控制的高速列车半主动悬挂控制***及方法 |
CN111055650B (zh) * | 2019-12-27 | 2021-11-23 | 江苏大学 | 磁流变半主动悬架粒子群-时滞依赖h∞鲁棒控制方法 |
CN112487553B (zh) * | 2020-11-18 | 2024-04-09 | 江苏大学 | 一种用于可控悬架***的时滞补偿控制器的设计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003045718A1 (fr) * | 2001-11-28 | 2003-06-05 | Peugeot Citroen Automobiles S.A. | Procede pour evaluer la frequence instantanee d'une excitation mecanique exercee sur une roue d'un vehicule automobile, et applications |
CN105172511A (zh) * | 2015-07-15 | 2015-12-23 | 西南交通大学 | 一种汽车悬架减振器控制***及方法 |
JP2016002844A (ja) * | 2014-06-16 | 2016-01-12 | トヨタ自動車株式会社 | 車両のばね上−ばね下相対速度推定装置 |
CN107599778A (zh) * | 2017-07-31 | 2018-01-19 | 江苏大学 | 一种电磁混合悬架的控制方法 |
CN107976904A (zh) * | 2017-11-09 | 2018-05-01 | 江苏大学 | 磁流变半主动悬架泰勒级数-二重h2时滞补偿控制方法 |
CN108032698A (zh) * | 2017-11-09 | 2018-05-15 | 江苏大学 | 磁流变半主动悬架泰勒级数-三重h2时滞补偿控制方法 |
CN109334377A (zh) * | 2018-09-11 | 2019-02-15 | 江苏大学 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
-
2018
- 2018-09-11 CN CN201811058385.1A patent/CN109334377B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003045718A1 (fr) * | 2001-11-28 | 2003-06-05 | Peugeot Citroen Automobiles S.A. | Procede pour evaluer la frequence instantanee d'une excitation mecanique exercee sur une roue d'un vehicule automobile, et applications |
JP2016002844A (ja) * | 2014-06-16 | 2016-01-12 | トヨタ自動車株式会社 | 車両のばね上−ばね下相対速度推定装置 |
CN105172511A (zh) * | 2015-07-15 | 2015-12-23 | 西南交通大学 | 一种汽车悬架减振器控制***及方法 |
CN107599778A (zh) * | 2017-07-31 | 2018-01-19 | 江苏大学 | 一种电磁混合悬架的控制方法 |
CN107976904A (zh) * | 2017-11-09 | 2018-05-01 | 江苏大学 | 磁流变半主动悬架泰勒级数-二重h2时滞补偿控制方法 |
CN108032698A (zh) * | 2017-11-09 | 2018-05-15 | 江苏大学 | 磁流变半主动悬架泰勒级数-三重h2时滞补偿控制方法 |
CN109334377A (zh) * | 2018-09-11 | 2019-02-15 | 江苏大学 | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 |
Non-Patent Citations (3)
Title |
---|
An optimal design for suspension based on LQG control;陈士安等;《2011 International Conference on Electric Information and Control Engineering》;IEEE;20111231;第1987页-1992页 * |
一种用于车辆半主动悬架控制的磁流变阻尼器模型;沈延等;《汽车工程》;20090531;第31卷(第5期);第462页-第466页 * |
磁流变半主动悬架的泰勒级数-LQG时滞补偿控制方法;陈士安等;《振动与冲击》;20170831;第36卷(第8期);第190页-第196页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109334377A (zh) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109334377B (zh) | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制***及其构造方法 | |
CN109334376B (zh) | 磁流变半主动悬架的减振器力学特性时变补偿控制***及其构造方法 | |
Sun et al. | Performance investigation of vehicle suspension system with nonlinear ball-screw inerter | |
CN106647256B (zh) | 一种基于h∞ pid的主动悬架控制***及控制方法 | |
CN106828005B (zh) | 一种预规定瞬态性能的汽车主动悬架自适应控制方法 | |
Wu et al. | A load-dependent PWA-H∞ controller for semi-active suspensions to exploit the performance of MR dampers | |
CN112009577B (zh) | 一种用于重卡驾驶室半主动悬置的控制方法 | |
CN111055650B (zh) | 磁流变半主动悬架粒子群-时滞依赖h∞鲁棒控制方法 | |
CN110597064A (zh) | 基于非线性和不确定模型的主动悬挂输出反馈控制方法 | |
Zhang et al. | A new SSUKF observer for sliding mode force tracking H∞ control of electrohydraulic active suspension | |
Carratù et al. | Semi-active suspension system for motorcycles: From the idea to the industrial product | |
CN109398020B (zh) | 一种基于非线性模型的车辆液电耦合式isd悬架的预测控制方法 | |
CN114488805A (zh) | 滑模干扰观测补偿磁流变阻尼器的悬架***和控制方法 | |
Feng et al. | Mu-synthesis PID control of full-car with parallel active link suspension under variable payload | |
Liang et al. | Optimal vibration control for nonlinear systems of tracked vehicle half-car suspensions | |
Chen et al. | Nonlinear analysis of a quasi-zero stiffness air suspension based on the cell-mapping method | |
Mahala et al. | Mathematical models for designing vehicles for ride comfort | |
CN109849932B (zh) | 一种路面自适应的车轮动载荷估计方法 | |
CN110077191B (zh) | 一种汽车主动悬挂***控制方法 | |
Zhao et al. | A linear variable parameter observer-based road profile height estimation for suspension nonlinear dynamics improvements | |
CN112526880B (zh) | 一种车辆行驶过程中的路面高度实时估计方法 | |
Canale et al. | Semi-active suspension control using" fast" model predictive control | |
Hanafi et al. | Speed effect to a quarter car ARX model based on system identification | |
Junyao et al. | Vibration control for vehicle active suspension based on ANFIS method | |
Zhu et al. | Optimal damping matching of semi-active air suspension based on vehicle comfort and handling |
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 |