[发明专利]基于比例积分观测器的非线性主动悬架容错追踪控制方法有效
申请号: | 201910351648.6 | 申请日: | 2019-04-28 |
公开(公告)号: | CN110435377B | 公开(公告)日: | 2022-06-07 |
发明(设计)人: | 庞辉;杨军杰;尚玉婷;姚睿;刘楠;王鹏 | 申请(专利权)人: | 西安理工大学 |
主分类号: | B60G17/018 | 分类号: | B60G17/018;B60G17/0185 |
代理公司: | 北京国昊天诚知识产权代理有限公司 11315 | 代理人: | 杨洲 |
地址: | 710048*** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: |
基于比例积分观测器的非线性主动悬架容错追踪控制方法,首先,为描述悬架系统的模型不确定性,建立基于Takagi‑Sugeno(T‑S)模糊模型的1/2车辆悬架模型;其次,基于该模型,开发鲁棒H |
||
搜索关键词: | 基于 比例 积分 观测器 非线性 主动 悬架 容错 追踪 控制 方法 | ||
【主权项】:
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∞性能指标
即:![]()
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安理工大学,未经西安理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201910351648.6/,转载请声明来源钻瓜专利网。