一种实时主动重构容错控制方法技术

一种实时主动重构容错控制方法，该方法基于滑模控制、有限时间控制技术和切比雪夫神经网络，能满足挠性飞行器系统在执行器出现故障和饱和等情况下的实时主动容错控制。引入只需要期望信号的切比雪夫神经网络来估计包含故障和饱和的系统总扰动，设计标称控制律和补偿控制律，补偿故障和饱和造成的影响，削弱滑模系统的固有抖振，提高姿态跟踪系统的精度。

【技术实现步骤摘要】
一种实时主动重构容错控制方法
本专利技术涉及一种实时主动重构容错控制方法，属于飞行器姿态控制领域。
技术介绍
随着科学技术的不断发展和应用领域的不断拓展，各类飞行器的结构和任务日益复杂和庞大，飞行环境恶劣多变，其可靠性已经成为飞行控制系统设计中的重要问题。通过控制系统的重构或者重组实现飞行器最低安全性要求，这对于保证飞行器顺利完成任务或避免坠毁具有重要意义。如何研发具有较强容错能力的飞行控制系统以满足高可靠性要求具有重要的学术价值和应用前景。容错控制技术是一种能适应环境的明显变化，可容许控制系统中的一个或多个部件失效的控制系统。容错控制的指导思想是一个控制系统一旦发生故障，系统依然能够维持其自身运行在安全状态，并在条件允许下满足一定的性能指标。滑模变结构控制是一种特殊的非线性不连续控制方法，这种控制方法与其他控制不同在于系统的结构在动态过程中，会根据系统当前的状态，使得系统按照预定滑动模态的状态轨迹运行。其设计与模型参数及扰动无关，使得变结构控制具有反应速度快、对参数变化不敏感、对扰动不敏感、物理实现简单等优点。切比雪夫神经网络在一定条件下可以任意精度逼近非线性函数，且具有较强的自学习、自适应和自组织能力。将切比雪夫神经网络与滑模变结构控制相结合，对模型不确定性和非线性部分估计补偿，一定程度上能够消除滑模控制的抖振问题。有限时间控制方法是一种非线性控制方法，是时间最优的控制方法，与渐近稳定的系统相比，有限时间稳定的系统在存在外部扰动和内部不确定的情况下不仅具有更快的收敛速度，还有更好的鲁棒性和抗扰能力。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，以挠性飞行器的主动容错控制为背景，提出一种基于滑模控制技术、有限时间控制技术和切比雪夫神经网络的飞行器实时主动重构容错控制方法。实现了挠性飞行器实时容错控制，在执行器出现乘性或加性故障下，最大程度满足飞行器姿态跟踪控制需求。本专利技术的技术解决方案是：一种实时主动重构容错控制方法，步骤如下：(1)建立挠性飞行器系统模型；(2)利用步骤(1)得到的所述挠性飞行器系统模型，基于四元数建立挠性飞行器运动学误差方程和动力学误差方程；(3)根据步骤(2)中的挠性飞行器运动学误差方程和动力学误差方程，建立有限时间非奇异终端滑模面；(4)根据切比雪夫神经网络以及步骤(3)中建立的有限时间非奇异终端滑模面，确定标称控制律un和补偿控制律ua，从而得到完整的主动重构容错控制器，进而实现了实时主动重构容错控制。本专利技术与现有技术相比的有益效果是：1、本专利技术控制方法能实现挠性飞行器实时主动重构容错控制。2、本专利技术控制方法将非奇快速终端滑模控制运用到挠性飞行器姿态跟踪控制领域，使系统在有限时间内快速稳定且避免了奇异问题3、本专利技术控制方法将神经网络与滑模控制相结合，提出用基函数仅依赖于期望信号的切比雪夫神经网络来高效地逼近系统未知总扰动。附图说明图1为本专利技术主动重构容错控制器结构框图；图2为本专利技术主动重构容错控制姿态误差及角速度误差；图3为本专利技术PID控制姿态误差及角速度误差；图4为本专利技术姿态角及角速度曲线；图5为本专利技术滑模面及控制力矩曲线；图6为本专利技术切比雪夫神经网络辅助控制器的输出曲线；图7为本专利技术挠性模态频率衰减曲线。具体实施方式下面结合附图对本专利技术的具体实施方式进行进一步的详细描述。如图1所示，本专利技术提出了一种实时主动重构容错控制方法，具体步骤如下：(1)考虑飞行器挠性特性、转动惯量不确定、外部扰动、执行器故障与饱和等因素的影响，建立如下挠性飞行器系统模型：其中：d∈R3是外部扰动，δ∈R4×3为刚体与挠性附件的耦合矩阵，δT是δ的转置，η为挠性模态，和分别为η的一阶导数和二阶导数；J0∈R3×3为已知的标称惯量矩阵，且为正定矩阵；ΔJ为惯量矩阵中的不确定部分，Ω＝[Ω1,Ω2,Ω3]T是飞行器在本体坐标系中的角速度分量，是Ω的一阶导数；×是运算符号，将×用于向量b＝[b1,b2,b3]T可得到：L＝diag{2ζiωni,i＝1,2,...,N}和分别为阻尼矩阵和刚度矩阵，N为模态阶数，ωni,i＝1,2,...,N为振动模态频率矩阵，ζi,i＝1,2,...,N为振动模态阻尼比；u＝[u1,u2,u3]T是主动重构容错控制器，sat(u)＝[sat(u1),sat(u2),sat(u3)]T是执行器产生的实际控制向量，sat(ui),i＝1,2,3表示执行器的非线性饱和特性且满足sat(ui)＝sign(ui)·min{umi,|ui|},i＝1,2,3，sat(ui)表述为sat(ui)＝θoi+ui,i＝1,2,3，其中θoi,i＝1,2,3为：umi,i＝1,2,3是执行器饱和值，超出执行器饱和值部分为θo＝[θo1,θo2,θo3]T，且满足‖θo‖≤lδθ，lδθ是正实数，Gδ＝[Gδ1,Gδ2,Gδ3]T是加性故障，即故障以加性方式影响系统且满足‖Gδ‖≤lδf，lδf是正实数；D＝diag{δo1,δo2,δo3}是执行器效率指标值且满足0<ετi≤δoi≤1,i＝1,2,3；0<ετi≤1,i＝1,2,3表示执行器最低执行能力，δoi＝1,i＝1,2,3表示第i个执行器工作正常；0<ετi≤δoi≤1,i＝1,2,3表示第i个执行器部分失效，但是该执行器仍能提供部分执行能力。(2)利用步骤(1)得到的挠性飞行器系统模型，基于四元数建立挠性飞行器运动学误差方程和动力学误差方程如下：挠性飞行器运动学误差方程：其中:(ev,e4)∈R3×R，ev＝[e1,e2,e3]T是当前飞行器姿态与期望姿态的误差四元数矢量部分，e4是标量部分，且满足和分别是ev、e4的一阶导数；(qv,q4)∈R3×R，qv＝[q1,q2,q3]T是描述飞行器姿态的单位四元数矢量部分，q4是标量部分，且满足qdv＝[qd1,qd2,qd3]T是描述期望姿态的单位四元数矢量部分，qd4是标量部分，且满足Ωe＝Ω-CΩd＝[Ωe1Ωe2Ωe3]T是建立在本体坐标系和目标坐标系之间的角速度误差向量，Ωd∈R3是期望角速度向量，是转换矩阵，且满足‖C‖＝1，是C的一阶导数，I3是3×3单位矩阵；挠性飞行器动力学误差方程为：其中，是Ωe的一阶导数，Ωd是期望角速度，是Ωd的一阶导数；挠性飞行器动力学误差方程改写为：其中：F为模型确定部分，R为未知总扰动；(3)根据步骤(2)中的挠性飞行器运动学误差方程和动力学误差方程，建立有限时间非奇异终端滑模面：S＝Ωe+K1ev+K2Sc(9)其中S＝[S1,S2,S3]T∈R3，Kj＝diag{kji}>0,i＝1,2,3,j＝1,2，diag(a1,a2,…,an)表示对角线元素为a1,a2,…,an的对角矩阵；且定义Sc＝[Sc1,Sc2,Sc3]T如下：其中r1,r2是正奇数，且0<r<1，l1i、l2i，i＝1,2,3是参数；εi,i＝1,2,3、ι1、ι2为设计参数，sign(a)是符号函数，定义如下：基于有限时间非奇异终端滑模面,如公式(9)所示，设计适当参数，当满足时，在有限时间内可实现控制目标{ev≡0,e4≡1,Ωe≡0}。(4)根据切比雪夫神经网络以及步骤(3)中建立的有限时间非奇异终端滑模本文档来自技高网...

【技术保护点】
一种实时主动重构容错控制方法，其特征在于步骤如下：(1)建立挠性飞行器系统模型；(2)利用步骤(1)得到的所述挠性飞行器系统模型，基于四元数建立挠性飞行器运动学误差方程和动力学误差方程；(3)根据步骤(2)中的挠性飞行器运动学误差方程和动力学误差方程，建立有限时间非奇异终端滑模面；(4)根据切比雪夫神经网络以及步骤(3)中建立的有限时间非奇异终端滑模面，确定标称控制律u

【技术特征摘要】
1.一种实时主动重构容错控制方法，其特征在于步骤如下：(1)建立挠性飞行器系统模型；(2)利用步骤(1)得到的所述挠性飞行器系统模型，基于四元数建立挠性飞行器运动学误差方程和动力学误差方程；(3)根据步骤(2)中的挠性飞行器运动学误差方程和动力学误差方程，建立有限时间非奇异终端滑模面；(4)根据切比雪夫神经网络以及步骤(3)中建立的有限时间非奇异终端滑模面，确定标称控制律un和补偿控制律ua，从而得到完整的主动重构容错控制器，进而实现了实时主动重构容错控制。2.根据权利要求1所述的一种实时主动重构容错控制方法，其特征在于：所述建立挠性飞行器系统模型，具体为：其中：d∈R3是外部扰动，δ∈R3×3为刚体与挠性附件的耦合矩阵，δT是δ的转置，η为挠性模态，和分别为η的一阶导数和二阶导数；J0∈R3×3为已知的标称惯量矩阵，且为正定矩阵；ΔJ为惯量矩阵中的不确定部分，Ω＝[Ω1,Ω2,Ω3]T是飞行器在本体坐标系中的角速度分量，是Ω的一阶导数；×是运算符号，将×用于向量b＝[b1,b2,b3]T得到：L＝diag{2ζiωni,i＝1,2,...,N}和分别为阻尼矩阵和刚度矩阵，N为模态阶数，ωni,i＝1,2,...,N为振动模态频率矩阵，ζi,i＝1,2,...,N为振动模态阻尼比；u＝[u1,u2,u3]T是主动重构容错控制器，sat(u)＝[sat(u1),sat(u2),sat(u3)]T是执行器产生的实际控制向量，sat(ui),i＝1,2,3表示执行器的非线性饱和特性且满足sat(ui)＝sign(ui)·min{umi,|ui|},i＝1,2,3，sat(ui)表述为sat(ui)＝θoi+ui,i＝1,2,3，其中θoi为：umi,i＝1,2,3是执行器饱和值，超出执行器饱和值部分为θo＝[θo1,θo2,θo3]T，且满足‖θo‖≤lδθ，lδθ是正实数，Gδ＝[Gδ1,Gδ2,Gδ3]T是加性故障，即故障以加性方式影响系统且满足‖Gδ‖≤lδf，lδf是正实数；D＝diag{δo1,δo2,δo3}是执行器效率指标值且满足0<ετi≤δoi≤1,i＝1,2,3；0<ετi≤1,i＝1,2,3表示执行器最低执行能力，δoi＝1,i＝1,2,3表示第i个执行器工作正常；0<ετi≤δoi≤1,i＝1,2,3表示第i个执行器部分失效，但是该执行器仍能提供部分执行能力。3.根据权利要求2所述的一种实时主动重构容错控制方法，其特征在于：建立挠性飞行器运动学误差方程和动力学误差方程具体为：挠性飞行器运动学误差方程：其中:(ev,e4)∈R3×R，ev＝[e1,e2,e3]T是当前飞行器姿态与期望姿态的误差四元数矢量部分，e4是标量部分，且满足和分别是ev、e4的一阶导数；(qv,q4)∈R3×R，qv＝[q1,q2,q3]T是描述飞行器姿态的单位四元数矢量部分，q4是标量部分，且满足qdv＝[qd1,qd2,qd3]T是描述期望...

【专利技术属性】
技术研发人员：路坤锋，李天涯，刘海亮，郑春胜，张隽，高磊，王辉，杨广慧，李德强，张喆，
申请(专利权)人：北京航天自动控制研究所，中国运载火箭技术研究院，
类型：发明
国别省市：北京,11