CN106681337A - 基于奇次滑模的平流层飞艇定高飞行控制方法 - Google Patents
基于奇次滑模的平流层飞艇定高飞行控制方法 Download PDFInfo
- Publication number
- CN106681337A CN106681337A CN201710008391.5A CN201710008391A CN106681337A CN 106681337 A CN106681337 A CN 106681337A CN 201710008391 A CN201710008391 A CN 201710008391A CN 106681337 A CN106681337 A CN 106681337A
- Authority
- CN
- China
- Prior art keywords
- dirigible
- height
- pitch
- angle
- sin
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 title abstract 2
- 230000006641 stabilisation Effects 0.000 claims abstract description 13
- 238000011105 stabilization Methods 0.000 claims abstract description 13
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 238000009415 formwork Methods 0.000 claims description 25
- 238000013461 design Methods 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000000205 computational method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 206010016256 fatigue Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/04—Control of altitude or depth
- G05D1/042—Control of altitude or depth specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0825—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了基于奇次滑模的平流层飞艇定高飞行控制方法,生成飞艇高度与垂向速度的测量以及高度误差;基于误差与微分的奇次函数的非线性滑模构建;基于线性与非线性增益组合的期望俯仰角指令生成;根据飞艇俯仰通道的微分方程所建立的数学模型,近似模拟飞艇俯仰通道的特性;将期望俯仰角指令输入给俯仰角姿态稳定跟踪控制器,并由该控制器生成俯仰舵偏角信号,将该舵偏角代入步骤四所建立的数学模型,最终确定一组飞艇定高飞行的控制方案参数,使得飞艇定高飞行具有满意的动态响应与稳态响应性能。
Description
技术领域
本发明属于飞行器控制技术领域,涉及一种基于奇次滑模的平流层飞艇定高飞行控制方法。
背景技术
早在20世纪60年代美国空军就已经开始尝试平流层平台的研制和试飞。进入20世纪90年代美国的Sky Station公司Skysat公司和Lockheed Martin公司都提出了实用飞艇平台的设计方案。目前国内外针对飞艇的控制主要集中于姿态控制,位置控制,温度控制。而其中位置控制包含高度控制与定点悬停。目前文献大部分集中于定点悬停,而且无论是姿态控制还是位置控制主要采用PID方法为多。
飞艇高度上的大空域快速机动与定高飞行控制还具有很多问题,动态响应与稳态响应性能不好。
发明内容
为实现上述目的,本发明提供一种基于奇次滑模的平流层飞艇定高飞行控制方法,解决飞艇高度上的大空域快速机动与定高飞行控制问题,实现了飞艇俯仰通道按照给定高度飞行的控制目的,飞艇定高飞行具有满意的动态响应与稳态响应性能。
本发明所采用的技术方案是,基于奇次滑模的平流层飞艇定高飞行控制方法,具体按照以下步骤进行:
步骤一:飞艇高度与垂向速度的测量以及高度误差的生成;
步骤二:基于误差与微分的奇次函数的非线性滑模构建;
步骤三:基于线性与非线性增益组合的期望俯仰角指令生成;
步骤四:利用计算机,根据飞艇俯仰通道的微分方程所建立的数学模型,近似模拟飞艇俯仰通道的特性;
步骤五:将期望俯仰角指令输入给俯仰角姿态稳定跟踪控制器,并由该控制器生成俯仰舵偏角信号,将该舵偏角代入步骤四所建立的数学模型,通过不断调整控制参数,并观察飞艇各状态的数据并画图,观测飞艇高度变化的数据曲线,分析定高飞行的动态响应,从而最终确定一组飞艇定高飞行的控制方案参数,使得飞艇定高飞行具有满意的动态响应与稳态响应性能。
进一步的,所述步骤一具体按照以下步骤进行:
有飞艇上高度表测量飞艇的实时高度,记为z,并通过A/D转换传递给艇上计算机,通过和期望高度信号的比较,生成高度误差信号,其中期望高度信号记为zd,高度误差信号记为ez,其满足ez=z-zd,采用垂直速度传感器测量飞艇的垂向速度,记为w,并通过A/D转换传递给艇上计算机,为第二步做准备。
进一步的,所述步骤二具体按照以下步骤进行:
首先由高度误差信号生成非线性其次误差项Π,其定义如下:
由上述高度误差信号,在艇上计算机中生成误差积分信号Ω,其定义如下:
Ω=∫ezdt
而误差微分信息其中定高飞行时期望高度为常值,故(为期望高度zd的导数,由于其一般为常值,因此导数为0),故有即可由测量的垂向速度信号w代替误差微分信号;
最终构成如下奇次非线性滑模信息Sz如下:
Cz、Czs的含义是滑模面中的控制参数,为正数,具体值可在设计调整中调整。
进一步的,所述步骤三具体按照以下步骤进行:
基于奇次滑模信息,构造俯仰角的期望值θd:
其中kz1sz为滑模的线性增益项,为滑模的非线性增益项;kz1、kz2、kz3、ξ1与ξ2为控制参数,选取为正常数。
进一步的,所述步骤四具体按照以下步骤进行:
飞艇俯仰通道的数学模型如下:
其中,u1为飞艇俯仰舵偏角,用于稳定与控制飞艇的俯仰姿态角;u2为飞艇的发动机推力,用于提高飞艇向前的飞行速度;f1-f6仅为变量,无物理含义,表达数为:
而a11,a13,a22,a31,a33为飞艇质量分布与转动惯量相关的参数,其计算方法通过下面M的逆阵获得,即满足
而M矩阵有飞艇的质量与转动惯量所决定,其求取方法如下:
I3为3阶单位矩阵;
其中,M3为M的子矩阵,用于计算M;m为飞艇的质量,az为飞艇结构参数,az=16.8;m11、m33、m55分别飞艇在不同方向的质量分布系数决定的参数,由飞艇质量分布与转动惯量所决定:m11=km1Mr,m33=km2Mr,m55=km3Iy,其中km1=0.1053;km2=0.8260;km3=0.1256;km1、km2、km3含义是飞艇的x,y,z三个方向的质量分布系数;Iy为飞艇沿y轴方向的转动惯量,Mr是飞艇所排开气体的质量,Mr=ρV,其中ρ为大气密度,V为飞艇的体积;
Q为动压头,其计算方法为Q=0.5ρVf 2;Vf为飞艇的运动速度;
为飞艇的前向飞行加速度;u为艇体坐标系中飞艇的前向飞行速度;
为飞艇的垂向飞行加速度;w为艇体坐标系中飞艇的垂向飞行速度;
为飞艇的俯仰角加速度;q为飞艇的俯仰角速度;
为飞艇的俯仰角速度,θ为飞艇的俯仰角;
为发射坐标系中飞艇的前向飞行速度;x为飞艇的前向飞行距离;
为发射坐标系中飞艇的垂向飞行速度;z为飞艇的飞行高度;
α含义是为飞艇向前与向上速度所形成的夹角;
kg1与kg2为舵效常数,为空气动力学系数,其数据来自于飞艇风洞试验;
CX1、CX2、Cz1、Cz2与Cz3为飞艇受力相关的空气动力系数,CM1、CM2、CM1为飞艇受力矩相关的空气动力系数;
针对上述复杂模型的分析,可以简化为如下一阶模型:
其中飞艇俯仰通道姿态稳定控制的设计是通过设计飞艇俯仰舵偏角u1来控制飞艇的俯仰角θ跟踪期望的姿态角指令θd。
本发明的有益效果是:本发明主要针对飞艇俯仰通道的高度控制问题,提出一类新型非线性奇次滑模控制,能够解决飞艇高度上的大空域快速机动与定高飞行控制问题,与现有方向相比,其能避免控制量的饱和,而且由于奇次滑模的引入,使得飞艇机动控制具有更好的鲁棒性与快速性。因此本发明不仅具有较好的创新性,而且具有较高的工程实用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的原理图;
图2是本发明实施例提供的500米定高飞行情况下的飞艇的前向运动速度曲线;
图3是本发明实施例提供的500米定高飞行情况下的飞艇的垂向运动速度曲线;
图4是本发明实施例提供的500米定高飞行情况下的俯仰角小幅波动曲线;
图5是本发明实施例提供的500米定高飞行情况下的飞艇的俯仰角速率曲线;
图6是本发明实施例提供的500米定高飞行情况下的飞艇的水平飞行距离曲线;
图7为本发明实施例提供的500米定高飞行情况下的飞艇的飞行高度曲线;
图8为本发明实施例提供的500米定高飞行情况下的飞艇的俯仰舵偏曲线;
图9是本发明实施例提供的500米定高飞行情况下的飞艇的俯仰角期望值。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明一种基于奇次滑模的平流层飞艇定高飞行控制方法,思路是:通过测量飞艇的高度与垂向速度,采用艇上计算机得出飞艇实时高度与期望高度的误差信号,再由该误差信号的奇次项、误差积分信号以及垂向速度信号的奇次项构成一类新的奇次非线性滑模信息,基于该新的滑模面,生成期望的飞艇俯仰角指令信号,再由飞艇俯仰角跟踪控制器实现对该指令跟踪,当俯仰角跟踪上该期望指令信号时,即实现了飞艇俯仰通道按照给定高度飞行,从而消除高度误差。
因此本发明中默认飞艇姿态稳定控制器已完成设计。本发明主要针对飞艇俯仰通道的高度控制问题,提出一类新的奇次非线性滑模控制,主要解决飞艇高度上的大空域快速机动与定高飞行控制问题,并尽量避免控制量的饱和。因此本发明不仅具有较好的创新性,而且具有较高的工程实用价值。
本发明采用的实施例原理如图1所示,具体按照以下步骤进行:
步骤一:飞艇高度与垂向速度的测量以及高度误差的生成;
有飞艇上高度表测量飞艇的实时高度,记为z,并通过A/D转换传递给艇上计算机,通过和期望高度信号的比较,生成高度误差信号,其中期望高度信号记为zd,高度误差信号记为ez,其满足ez=z-zd,采用垂直速度传感器测量飞艇的垂向速度,记为w,并通过A/D转换传递给艇上计算机,为第二步做准备。
步骤二:基于误差与微分的奇次函数的非线性滑模构建;
首先由高度误差信号生成非线性其次误差项Π,其定义如下:
由上述高度误差信号,在艇上计算机中生成误差积分信号Ω,其定义如下:
Ω=∫ezdt
而误差微分信息其中定高飞行时期望高度为常值,故故有即可由测量的垂向速度信号w代替误差微分信号。
最终构成如下奇次非线性滑模信息Sz如下:
Cz、Czs的含义是滑模面中的控制参数,为正数,具体值可在设计调整中调整;
步骤三:基于线性与非线性增益组合的期望俯仰角指令生成;
基于奇次滑模信息,构造俯仰角的期望值θd:
其中kz1sz为滑模的线性增益项,为滑模的非线性增益项。
kz1、kz2、kz3、ξ1与ξ2为控制参数,可选取为正常数。其具体选取见示例与仿真。
步骤四:利用计算机,根据下面的飞艇俯仰通道的微分方程所建立的数学模型,近似模拟飞艇俯仰通道的特性。
为了确保上述步骤一至步骤三中控制器的参数选取合理,可用通过计算机模拟仿真的手段进行编程,从而进行参数调整。其中飞艇俯仰通道的数学模型如下:
其中,u1为飞艇俯仰舵偏角,用于稳定与控制飞艇的俯仰姿态角;u2为飞艇的发动机推力,用于提高飞艇向前的飞行速度;f1-f6仅为变量,无物理含义,表达数为:
而a11,a13,a22,a31,a33为飞艇质量分布与转动惯量相关的参数,其计算方法通过下面M的逆阵获得,即满足
而M矩阵有飞艇的质量与转动惯量所决定,其求取方法如下:
I3为3阶单位矩阵。
其中,M3是为M的子矩阵,用于计算M;m为飞艇的质量,az为飞艇结构参数,az=16.8;m11、m33、m55分别飞艇在不同方向的质量分布系数决定的参数,由飞艇质量分布与转动惯量所决定:m11=km1Mr,m33=km2Mr,m55=km3Iy,其中km1=0.1053;km2=0.8260;km3=0.1256;km1、km2、km3含义是飞艇的x,y,z三个方向的质量分布系数;Iy为飞艇沿y轴方向的转动惯量,Mr是飞艇所排开气体的质量,Mr=ρV,其中ρ为大气密度,V为飞艇的体积。
Q为动压头,其计算方法为Q=0.5ρVf 2;Vf为飞艇的运动速度。
为飞艇的前向飞行加速度;u为艇体坐标系中飞艇的前向飞行速度;
为飞艇的垂向飞行加速度;w为艇体坐标系中飞艇的垂向飞行速度;
为飞艇的俯仰角加速度;q为飞艇的俯仰角速度;
为飞艇的俯仰角速度,θ为飞艇的俯仰角;
为发射坐标系中飞艇的前向飞行速度;x为飞艇的前向飞行距离;
为发射坐标系中飞艇的垂向飞行速度;z为飞艇的飞行高度;
α含义是为飞艇向前与向上速度所形成的夹角;
kg1与kg2为舵效常数,为空气动力学系数,其数据来自于飞艇风洞试验。
CX1、CX2、Cz1、Cz2与Cz3为飞艇受力相关的空气动力系数,CM1、CM2、CM1为飞艇受力矩相关的空气动力系数,各型飞艇的计算方式略有不同,其数据来自于飞艇的风洞实验数据,非本发明所保护与所讨论的内容,故不详细累述。
针对上述复杂模型的分析,可以简化为如下一阶模型:
其中飞艇俯仰通道姿态稳定控制的设计是通过设计飞艇俯仰舵偏角u1来控制飞艇的俯仰角θ跟踪期望的姿态角指令θd,有关姿态稳定控制器的设计在此不再详细讨论,非本专利的核心内容,本发明是在上述姿态稳定控制器设计完成的基础上进行的。
根据上述简化的一阶模型,本发明的基本思想是根据高度误差信号生成期望的姿态角指令θd。
步骤五:将步骤一至步骤三所得的期望俯仰角指令输入给俯仰角姿态稳定跟踪控制器,并由该控制器生成俯仰舵偏角信号,将该舵偏角代入步骤四所建立的模型,通过不断调整控制参数,并观察飞艇各状态的数据并画图,尤其是观测飞艇高度变化的数据曲线,分析定高飞行的动态响应,从而最终确定一组飞艇定高飞行的控制方案参数,使得飞艇定高飞行具有满意的动态响应与稳态响应性能。
实施例:
首先采用PID控制规律,设定飞艇的俯仰角稳定控制器,也可采用其它控制规律设计俯仰角稳定控制器,实现飞艇姿态稳定的控制功能,在此不再详细阐述姿态稳定控制器的设计,主要是由于本发明的核心技术在于高度指令的生成,因此本发明的讨论是在默认姿态稳定控制器设计完好的基础上进行的,故此处仅以PID姿态稳定控制为例,以完成高度控制的示例。
在上述姿态稳定控制器设计完好的基础上,设定飞艇初始高度为0米,初始速度为0m/s,假设期望的给定高度为500米。并设置发动机推力为常值u2=8000,飞艇最终稳定的前向飞行速度为25m/s左右。
则按照上述步骤一至步骤三,选取方案参数为cz=0.06,czs=0.003,kz1=0.1,kz2=5,kz3=2,ξ1=10,ξ2=6。
最终得到的期望俯仰角信号如下:
代入步骤四中模型,得到的结果如图2至图8所示。
通过以上仿真结果与曲线可以看出,在定高飞行阶段,俯仰角最终能够稳定在0附近附近,有波动但周期比较长,频率比较低,符合飞艇定高飞行运动特点与物理特性。由图7可以看出飞艇可以稳定在期望的500米高度,因此本发明是完全可行的,且具有较高的工程应用价值。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (5)
1.基于奇次滑模的平流层飞艇定高飞行控制方法,其特征在于,具体按照以下步骤进行:
步骤一:飞艇高度与垂向速度的测量以及高度误差的生成;
步骤二:基于误差与微分的奇次函数的非线性滑模构建;
步骤三:基于线性与非线性增益组合的期望俯仰角指令生成;
步骤四:利用计算机,根据飞艇俯仰通道的微分方程所建立的数学模型,近似模拟飞艇俯仰通道的特性;
步骤五:将期望俯仰角指令输入给俯仰角姿态稳定跟踪控制器,并由该控制器生成俯仰舵偏角信号,将该舵偏角代入步骤四所建立的数学模型,通过不断调整控制参数,并观察飞艇各状态的数据并画图,观测飞艇高度变化的数据曲线,分析定高飞行的动态响应,从而最终确定一组飞艇定高飞行的控制方案参数,使得飞艇定高飞行具有满意的动态响应与稳态响应性能。
2.根据权利要求1所述的基于奇次滑模的平流层飞艇定高飞行控制方法,其特征在于,所述步骤一具体按照以下步骤进行:
有飞艇上高度表测量飞艇的实时高度,记为z,并通过A/D转换传递给艇上计算机,通过和期望高度信号的比较,生成高度误差信号,其中期望高度信号记为zd,高度误差信号记为ez,其满足ez=z-zd,采用垂直速度传感器测量飞艇的垂向速度,记为w,并通过A/D转换传递给艇上计算机,为第二步做准备。
3.根据权利要求2所述的基于奇次滑模的平流层飞艇定高飞行控制方法,其特征在于,所述步骤二具体按照以下步骤进行:
首先由高度误差信号生成非线性其次误差项Π,其定义如下:
由上述高度误差信号,在艇上计算机中生成误差积分信号Ω,其定义如下:
Ω=∫ezdt
而误差微分信息其中定高飞行时期望高度为常值,故期望高度信号zd的导数故有即可由测量的垂向速度信号w代替误差微分信号;
最终构成如下奇次非线性滑模信息Sz如下:
Cz、Czs的含义是滑模面中的控制参数,为正数。
4.根据权利要求3所述的基于奇次滑模的平流层飞艇定高飞行控制方法,其特征在于,所述步骤三具体按照以下步骤进行:
基于奇次滑模信息,构造俯仰角的期望值θd:
其中kz1sz为滑模的线性增益项,为滑模的非线性增益项;kz1、kz2、kz3、ξ1与ξ2为控制参数,选取为正常数。
5.根据权利要求4所述的基于奇次滑模的平流层飞艇定高飞行控制方法,其特征在于,所述步骤四具体按照以下步骤进行:
飞艇俯仰通道的数学模型如下:
其中,u1为飞艇俯仰舵偏角,用于稳定与控制飞艇的俯仰姿态角;u2为飞艇的发动机推力,用于提高飞艇向前的飞行速度;f1-f6仅为变量,无物理含义,表达数为:
而a11,a13,a22,a31,a33为飞艇质量分布与转动惯量相关的参数,其计算方法通过下面M的逆阵获得,即满足
而M矩阵有飞艇的质量与转动惯量所决定,其求取方法如下:
I3为3阶单位矩阵;
其中,M3为M的子矩阵,用于计算M;m为飞艇的质量,az为飞艇结构参数,az=16.8;m11、m33、m55分别飞艇在不同方向的质量分布系数决定的参数,由飞艇质量分布与转动惯量所决定:m11=km1Mr,m33=km2Mr,m55=km3Iy,其中
km1=0.1053;km2=0.8260;km3=0.1256;km1、km2、km3含义是飞艇的x,y,z三个方向的质量分布系数;Iy为飞艇沿y轴方向的转动惯量,Mr是飞艇所排开气体的质量,Mr=ρV,其中ρ为大气密度,V为飞艇的体积;
Q为动压头,其计算方法为Q=0.5ρVf 2;Vf为飞艇的运动速度;
为飞艇的前向飞行加速度;u为艇体坐标系中飞艇的前向飞行速度;
为飞艇的垂向飞行加速度;w为艇体坐标系中飞艇的垂向飞行速度;
为飞艇的俯仰角加速度;q为飞艇的俯仰角速度;
为飞艇的俯仰角速度,θ为飞艇的俯仰角;
为发射坐标系中飞艇的前向飞行速度;x为飞艇的前向飞行距离;
为发射坐标系中飞艇的垂向飞行速度;z为飞艇的飞行高度;
α含义是为飞艇向前与向上速度所形成的夹角;
kg1与kg2为舵效常数,为空气动力学系数,其数据来自于飞艇风洞试验;
CX1、CX2、Cz1、Cz2与Cz3为飞艇受力相关的空气动力系数,CM1、CM2、CM1为飞艇受力矩相关的空气动力系数;
针对上述复杂模型的分析,可以简化为如下一阶模型:
其中飞艇俯仰通道姿态稳定控制的设计是通过设计飞艇俯仰舵偏角u1来控制飞艇的俯仰角θ跟踪期望的姿态角指令θd。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710008391.5A CN106681337B (zh) | 2017-01-05 | 2017-01-05 | 基于奇次滑模的平流层飞艇定高飞行控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710008391.5A CN106681337B (zh) | 2017-01-05 | 2017-01-05 | 基于奇次滑模的平流层飞艇定高飞行控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106681337A true CN106681337A (zh) | 2017-05-17 |
CN106681337B CN106681337B (zh) | 2019-03-15 |
Family
ID=58850242
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710008391.5A Active CN106681337B (zh) | 2017-01-05 | 2017-01-05 | 基于奇次滑模的平流层飞艇定高飞行控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106681337B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109625248A (zh) * | 2018-11-23 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | 给定高度控制方法 |
CN110244754A (zh) * | 2019-06-28 | 2019-09-17 | 上海工程技术大学 | 一种用于平流层浮空器定点驻空的控制***及方法 |
CN111650947A (zh) * | 2020-07-06 | 2020-09-11 | 上海交通大学 | 一种平流层飞艇高度非线性控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5431084A (en) * | 1992-05-19 | 1995-07-11 | United Technologies Corporation | Helicopter integrated fire and flight control having azimuth and pitch control |
CN102205876A (zh) * | 2010-12-22 | 2011-10-05 | 中国特种飞行器研究所 | 平流层飞艇定点高度保持综合控制方法 |
CN105173051A (zh) * | 2015-09-15 | 2015-12-23 | 北京天航华创科技股份有限公司 | 一种平流层飞艇的制导控制一体化及控制分配方法 |
CN105718727A (zh) * | 2016-01-19 | 2016-06-29 | 中国科学院光电研究院 | 一种估计平流层飞艇飞行性能参数的方法和*** |
-
2017
- 2017-01-05 CN CN201710008391.5A patent/CN106681337B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5431084A (en) * | 1992-05-19 | 1995-07-11 | United Technologies Corporation | Helicopter integrated fire and flight control having azimuth and pitch control |
CN102205876A (zh) * | 2010-12-22 | 2011-10-05 | 中国特种飞行器研究所 | 平流层飞艇定点高度保持综合控制方法 |
CN105173051A (zh) * | 2015-09-15 | 2015-12-23 | 北京天航华创科技股份有限公司 | 一种平流层飞艇的制导控制一体化及控制分配方法 |
CN105718727A (zh) * | 2016-01-19 | 2016-06-29 | 中国科学院光电研究院 | 一种估计平流层飞艇飞行性能参数的方法和*** |
Non-Patent Citations (3)
Title |
---|
柯艺杰等: "平流层飞艇俯仰附加转动惯量的估算与辨识", 《飞行力学》 * |
王延等: "基于变结构动态逆控制的平流层飞艇时间最优航迹设计", 《测控技术》 * |
肖鸣超等: "一种平流层飞艇高度控制方法", 《电光与控制》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109625248A (zh) * | 2018-11-23 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | 给定高度控制方法 |
CN109625248B (zh) * | 2018-11-23 | 2022-04-01 | 中国航空工业集团公司沈阳飞机设计研究所 | 给定高度控制方法 |
CN110244754A (zh) * | 2019-06-28 | 2019-09-17 | 上海工程技术大学 | 一种用于平流层浮空器定点驻空的控制***及方法 |
CN110244754B (zh) * | 2019-06-28 | 2021-12-10 | 上海工程技术大学 | 一种用于平流层浮空器定点驻空的控制***及方法 |
CN111650947A (zh) * | 2020-07-06 | 2020-09-11 | 上海交通大学 | 一种平流层飞艇高度非线性控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106681337B (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106527122A (zh) | 平流层飞艇定高飞行非线性pid控制方法 | |
CN108490788B (zh) | 一种基于双干扰观测的飞行器俯仰通道反演控制方法 | |
CN104950901B (zh) | 无人直升机姿态误差有限时间收敛非线性鲁棒控制方法 | |
CN102073755B (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
CN104133479B (zh) | 一种采用单轴气浮台模拟挠性卫星三轴姿态耦合运动的测试方法 | |
CN109614633A (zh) | 一种复合式旋翼飞行器非线性建模方法及配平方法 | |
Lu et al. | Real-time simulation system for UAV based on Matlab/Simulink | |
CN107608367A (zh) | 多变量干扰补偿四旋翼无人机轨迹与姿态协同控制方法 | |
CN107390545A (zh) | 一种无人机及其载荷的仿真训练*** | |
CN109635494A (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN105759609B (zh) | 一种基于显式非线性模型预测控制的舰载机自主着舰方法 | |
CN105182989A (zh) | 一种风场影响下的飞机姿态控制方法 | |
CN106681337A (zh) | 基于奇次滑模的平流层飞艇定高飞行控制方法 | |
CN109358634A (zh) | 一种高超声速飞行器鲁棒自适应控制方法 | |
CN109507890A (zh) | 一种基于eso的无人机动态逆广义预测控制器 | |
CN105652880B (zh) | 用于飞行器大空域飞行的非线性抗饱和高度指令生成方法 | |
CN106406333B (zh) | 一种基于积分型终端滑模的平流层飞艇俯仰角跟踪方法 | |
CN111208844B (zh) | 基于过载角速度的飞行器纵向过载非线性控制方法 | |
CN107505846A (zh) | 一种空间机械臂***抗干扰姿态协调验证装置及控制方法 | |
CN107065544A (zh) | 基于攻角幂函数的高超飞行器神经网络控制方法 | |
CN102654772A (zh) | 一种基于控制力受限情况下飞行器航迹倾角反演控制方法 | |
CN103400035A (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN107063244A (zh) | 一种飞行器飞行过程模拟方法 | |
CN106774385B (zh) | 一种采用自适应变结构的飞艇定点悬停控制方法 | |
CN106484969A (zh) | 一种大包线强机动飞行器动力学高精度仿真方法 |
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 |