[发明专利]一种导弹多驱动状态执行机构故障的容错控制方法有效
申请号: | 201911096907.1 | 申请日: | 2019-11-11 |
公开(公告)号: | CN110727198B | 公开(公告)日: | 2022-06-14 |
发明(设计)人: | 李彤;陈小庆;孟志鹏;王波;吕良 | 申请(专利权)人: | 中国人民解放军军事科学院国防科技创新研究院 |
主分类号: | G05B13/04 | 分类号: | G05B13/04 |
代理公司: | 中国兵器工业集团公司专利中心 11011 | 代理人: | 刘瑞东 |
地址: | 100071 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 导弹 驱动 状态 执行机构 故障 容错 控制 方法 | ||
1.一种导弹多驱动状态执行机构故障的容错控制方法,其特征在于,包括如下步骤:
步骤一、导弹模型建立
S1针对一般无翼导弹,在满足冻结系数和小扰动假设的条件下构建如下表示的导弹动力学模型:
其中,θ、σ、α、β和γ分别为导弹速度倾角、航迹偏航角、攻角、侧滑角和滚转角(rad);P、Y和Z分别代表推力、升力和侧向力(N);m和V分别代表导弹质量(kg)和速度(m/s);ωx1、ωy1和ωz1分别为弹体滚转、偏航和俯仰角速度(rad/s);Jx1、Jy1和Jz1分别为弹体体系x、y和z轴转动惯量(kg·m);mx1 ,my1 和mz1 分别为滚转、偏航和俯仰方向力矩(N·m);δ1,δ2,δ3和δ4依次为#1舵至#4舵的舵偏角(O);Sm和lk分别代表弹体特征面积(m2)和特征长度(m);和分别为滚转、偏航和俯仰方向无因次角速度;代表动压(Pa);偏导数上标代表相应的气动力系数和气动力矩系数;
S2、定义{σ1,σ2,σ3,σ4}-rj,(j=1,2,3...11)形式表示导弹执行机构工作模式,其中σi(i=1,2,3,4)代表故障状态,σi=1,舵机工作状态良好,若σi接近于0,则舵机发生故障,其模式主要包括以下几种:
(1)无执行机构完全失效的过驱动状态:
{1,1,1,1}-r1
(2)单个执行机构完全失效的全驱动状态:
{1,1,1,0}-r2,{1,1,0,1}-r3,{1,0,1,1}-r4,{0,1,1,1}-r5
(3)两个执行机构完全失效的欠驱动状态:
{1,1,0,0}-r6,{1,0,1,0}-r7,{1,0,0,1}-r8,{0,1,1,0}-r9,{0,1,0,1}-r10,{0,0,1,1}-r11
步骤二、针对导弹执行机构的不同驱动状态重构控制策略
S1、模式r1至r5的标准控制策略构建
设计控制系统使导弹纵向和侧向跟踪设计弹道,滚转通道保持稳定,即控制目标变量为和γc=0;定义导弹系统状态为x=[α,β,γ,ωx1,ωy1,ωz1]T;系统输入u=[δ1,δ2,δ3,δ4]T;系统输出y=[qH,ql,γ]T;
则系统状态方程描述如下:
其中,
对系统进行基于Lie导数的微分同胚映射可得
系统动态特性由反馈线性化状态进行描述
其中,
为在rj故障模式下无完全失效故障舵机序号的集合;
定义系统跟踪指令为yc=[yc1,yc2,yc3]T=[qHc,qlc,γc]T,采用反步控制方法设计如下
其中,和为系统跟踪误差,和为反步控制稳定函数,通过设计令相应状态二次型形式Lyapunov函数导数为负;
其中,和为稳定函数设计增益,且需满足跟踪性能特征多项式为Hurwitz多项式,为不同通道跟踪性能设计带宽;
根据系统输入矩阵在标准控制策略下是非奇异的,通过控制分配得到系统输入为:
其中,表示矩阵的广义逆矩阵
S2、模式r6至r11的欠驱动控制策略构建
针对故障模式所导致的欠驱动输入状态,其控制输入作如下变换
其中,δa和δb为故障模式rj下无完全失效故障舵机舵偏输入,因此,系统输入矩阵相应变换为系统状态方程表示如下
针对不同欠驱动故障模式,设计成形函数为
成形函数包含唯一变量γ,即成形变量,且存在关于γ的唯一积分形式;应用成形函数对系统输出进行变换如下
对系统进行基于Lie导数的微分同胚映射可得
则系统动态特性描述如下
其中,
并满足
由此,通过应用成形函数进行变换将系统转变为串联系统;
对于和子系统应用反步控制如下
其中,符号定义同上不变,反步控制设计中和被替换成qH和ql来保证弹道跟踪精度;
根据设计跟踪指令为
则对于子系统反步控制设计如下
根据反步控制设计可得等效控制输入为
从而有
由此,欠驱动控制策略舵偏输入为有效舵偏δa和δb,以及故障舵偏零输入;
步骤三、自适应故障识别因子重构策略
根据步骤二所设计的控制策略,控制系统采用自适应故障模式识别因子,对不同故障模式进行识别和控制策略切换;
S1重构机制设计
设置重构机制为:其中,为故障模式识别因子κj的估计值,定义故障识别因子κj满足以下条件
S2识别因子设计
设计系统Lyapunov函数如下
其中,Vx和Vκ分别对应系统跟踪状态和识别因子估计误差Lyapunov函数,为κ=[κ1,κ2,...,κ11]T的估计值,Λ=diag{Λ1,Λ2,...,Λ11}为自适应趋近律;
令Lyapunov函数导数负定,同时其中故障因子估计误差部分导数为0,则得到故障识别因子自适应律如下
其中,
以及
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国人民解放军军事科学院国防科技创新研究院,未经中国人民解放军军事科学院国防科技创新研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201911096907.1/1.html,转载请声明来源钻瓜专利网。