主权项: |
1.基于比例积分观测器的非线性主动悬架容错追踪控制方法,其特征在于,包括以下步骤: 步骤1,建立线性主动悬架模型 根据牛顿第二定律,1/2车辆主动悬架系统动力学方程为: 式中,ms和Iy分别为车身质量和车身转动惯量;zc,分别表示车身的垂向位移和俯仰角;muf,mur为前后悬架车轮质量;zuf,zur分别为前后车轮的垂向位移;F1,F2为前、后悬架组件中的弹性力;Ftf,Ftr分别为前、后轮胎组件中的弹性力;其表达式分别为: 其中:cf,cr为前后悬架阻尼系数;kf,kr为前后悬架的刚度;zrf,zrr为前后轮的路面输入位移;ktf,ktr为前后轮胎刚度系数,uFf,uFr为前后悬架作动器产生的主动控制力,Δyf=zsf-zuf为前悬架动行程,Δyr=zsr-zur为后悬架动行程,且 考虑选取系统状态变量为干扰输入w(t)=[zrf,zrr]T,测量输出以及控制输出则主动悬架状态空间方程可描述为: 其中,A、B、B1、C、D、C1、D1、E1均为状态空间的系数矩阵;具体可写为: D=[0], 系数矩阵A和C1中的相应元素如下所示: 步骤2,基于T-S模糊模型建立非线性主动悬架模型 由于不同工况下乘客的数量和车身载荷在一定范围内变化,因此本文考虑车身质量ms、前悬架车轮质量muf和后悬架车轮质量mur存在不确定性,即ms∈[msmin,msmax],muf∈[mufmin,mufmax],mur∈[murmin,murmax],基于上述情况,结合T-S模糊模型建立悬架系统的非线性模型,为此定义 其中,ξ1(t)=1/ms,ξ2(t)=1/muf,ξ3(t)=1/mur,此外,M1(ξ1(t)),M2(ξ1(t)),N1(ξ2(t)),N2(ξ2(t)),O1(ξ3(t)),O2(ξ3(t))为隶属函数,且可描述为: 其中隶属函数M1(ξ1(t))和M2(ξ1(t))分别定义为Heavy和Light,隶属函数N1(ξ2(t))和N2(ξ2(t))分别定义为Heavy和Light,同时隶属函数O1(ξ3(t))和O2(ξ3(t))分别定义为Heavy和Light,且满足M1(ξ1(t))+M2(ξ1(t))=1,N1(ξ2(t))+N2(ξ2(t))=1,O1(ξ3(t))+O2(ξ3(t))=1, 此时,考虑悬架系统质量不确定性的悬架系统(4)可用如下T-S模型来描述: 模型规则i:如果ξ1(t)为Mr,ξ2(t)为Nj,ξ3(t)为Ol,那么 其中,r=1,2(j=1,2;l=1,2;i=1,2,3,4,5,6,7,8);将矩阵A、B、B1、C、D、C1、D1、E1中的ms、muf、mur分别替换为ms min(或ms max)、muf min(或muf max)、mur min(或mur max)即可得到矩阵Ai、Bi、B1i、Ci、Di、C1i、D1i、E1i; 由上述模型规则可得,考虑悬架系统质量不确定性的车辆非线性主动悬架系统可通过T-S模糊模型表示为: 其中: h1(ξ(t))=M1(ξ1(t))×N1(ξ2(t))×O1(ξ3(t)), h2(ξ(t))=M1(ξ1(t))×N2(ξ2(t))×O1(ξ3(t)), h3(ξ(t))=M2(ξ1(t))×N1(ξ2(t))×O1(ξ3(t)), h4(ξ(t))=M2(ξ1(t))×N2(ξ2(t))×O1(ξ3(t)), h5(ξ(t))=M2(ξ1(t))×N2(ξ2(t))×O2(ξ3(t)), h6(ξ(t))=M2(ξ1(t))×N1(ξ2(t))×O2(ξ3(t)), h7(ξ(t))=M1(ξ1(t))×N2(ξ2(t))×O2(ξ3(t)), h8(ξ(t))=M1(ξ1(t))×N1(ξ2(t))×O2(ξ3(t)). 其中hi(ξ(t))为模糊权重函数,且hi(ξ(t))≥0,满足 步骤3,设计正常模式下的鲁棒H∞输出反馈控制器 在悬架系统无故障状态下时,主动悬架系统采用鲁棒H∞输出反馈控制器: u(t)=K(s)y(t) (9) 且待设计的控制律u(t)满足如下状态空间表达形式 式中,xc(t)为H∞输出反馈控制器的状态向量,Ac、Bc、Cc、Dc是待设计的控制器参数矩阵, 将待设计的控制律(10)应用于无故障系统(8),可得1/2车辆主动悬架闭环控制系统 其中,xcl(t)=[x(t),xc(t)]T为系统的状态向量,其它各系数矩阵分别为: Cci=[C1i+D1iDcCi D1iCc],Dci=[E1i]. 因此,将设计鲁棒H∞输出反馈控制器转化为求解参数矩阵Aci、Bci、Cci、Dci使得闭环系统(11)保持渐近稳定,且满足H∞性能指标||Tzw||∞<β,即干扰输入w(t)到控制输出z(t)传递函数的H∞范数||Tzw||∞小于给定上界β,根据有界实引理,满足这一设计要求的充要条件是存在一个对称正定矩阵Q,使得: 在Matlab中利用hinflmi求解器确定参数矩阵Ac、Bc、Cc、Dc,得形如式(10)的控制律u(t),并将其代入矩阵Aci、Bci、Cci、Dci,可得基于鲁棒H∞输出反馈控制器的1/2车辆主动悬架闭环控制系统; 步骤4,建立悬架故障模型 当作动器发生偏差故障,同时也将悬架系统质量不确定性考虑在内时,故障悬架模型为: 其中,f(t)表示因作动器偏差故障导致作动器产生的控制力偏差,即故障信号,作动器的偏差故障是车辆悬架系统中的最常见故障情况,因此作动器故障受到了广泛的关注,自然磨损、未知的外部严重干扰、参数变化、系统老化和部分堵塞都是造成作动器偏差故障的原因; 步骤5,设计基于比例积分观测器的容错追踪控制器 在考虑作动器偏差故障时,基于比例积分观测器设计主动悬架系统的容错追踪控制器为: 其中,代表容错追踪控制器的补偿力,K1i为待确定的状态反馈增益矩阵,和为故障及状态向量的估计值, 为了得到上述容错追踪控制律,建立如下的比例积分观测器: 其中,H1i,H2i是观测器的增益矩阵,在容错追踪控制策略中,通过获得故障系统的信息,来实现对故障信号的检测和隔离, 由式(13)和式(15)可得故障系统与观测器系统的输出误差为: 式中 定义动态跟踪误差e(t)=x(t)-xf(t),则 其中 由式(13)和式(15)可得故障系统与观测器系统的状态向量误差为 其中 综合式(17)和式(18),可得如下动态增广误差系统 其中 如果存在H1i、H2i、K1j和正定矩阵X1>0、P2>0,且有γ>0、P3=I,使得下面的凸优化问题成立: 其中,Qij=[BiK1j 0 0]T,W=[X1 0 0]T,Γ=[X1 I I]T,*表示矩阵中的对称转置项, 则有以下的设计指标成立: 1)闭环系统(19)是渐进稳定的; 2)在零初始条件下,对假设从干扰输入v(t)到误差输出的闭环传递函数为那么闭环系统(19)满足H∞性能指标即: |