CN108108524A - 一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 - Google Patents
一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 Download PDFInfo
- Publication number
- CN108108524A CN108108524A CN201711231763.7A CN201711231763A CN108108524A CN 108108524 A CN108108524 A CN 108108524A CN 201711231763 A CN201711231763 A CN 201711231763A CN 108108524 A CN108108524 A CN 108108524A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- aircraft
- omega
- mfrac
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,首先对在轨滑行、脱架至火箭脱落和上升到安全高度、速度飞行三个阶段的运动机理进行分析,以确定各飞行阶段作用在飞机上的力和力矩,再利用三状态机体运动方程组建立数学模型,并运用Matlab仿真环境以完成数模仿真,进而获得可靠的无人机安全发射包线,为无人机安全发射起飞提供有力技术支撑。
Description
技术领域
本发明涉及无人机起飞仿真技术领域,尤其涉及一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法。
背景技术
无人机发射系统是无人机的一个重要功能系统,而无人机安全发射包线则属于无人机发射系统的核心参数,旨在保证无人机顺利升空以达到安全的高度与速度飞行。
对于采用火箭助推发射方式起飞的无人机,发射起飞阶段故障率高,风险大,因此为确保飞行安全,需提供无人机安全发射包线,一套成熟的用于火箭助推无人机起飞发射段仿真方法显得十分重要。由于助推火箭在很短的时间(一般为2s左右)内向无人机提供大量能量,使其速度迅速增加,以达到安全速度,并达到一定的安全高度,时间很短但划分为三个连续阶段:在轨滑行、脱架至火箭脱落和上升到安全高度、速度飞行阶段;目前采用的发射起飞段仿真方法,虽然也是基于三状态机体动力学方程组,能够得到具有一定意义的发射起飞段特性,但仍然存在以下不足:
1)针对单个或部分阶段进行分离式研究,虽能够获得某个或部分阶段的力和力矩特性,但难以全局掌握整个发射起飞阶段的飞行特性,而实际飞行中,既希望能够安全脱轨、脱架,也能够安全爬升达到安全高度、速度,每一个阶段都很重要,连续考虑才能全局了解飞机飞行特性,掌握安全包线,显然仅考虑某个或部分阶段的飞行特性是不够;
2)对三状态机体动力学方程组进行了较多简化假设,如忽略了作用在飞机上的导轨摩擦力、前后支架支反力、质量的变化等,严重降低模型的逼真度,从参数敏感性分析,敏感性高的不容简化,而忽略敏感性低也易导致精准度不足;
3)对在轨滑行阶段的处理,忽略了脱钩或解锁前飞机受发射架牵引约束,飞行轨迹实际是弧形运动,简单以平轨滑行的方式处理,在俯仰角速度等参数上过于保守,进而造成模型失真;此外,对于以俯仰角速度作为输入的控制律设计,忽略牵引约束,既无法模拟真实飞行情况,也不利于指导飞行。
发明内容
本发明所解决的技术问题在于提供一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,以解决上述背景技术中的缺点。
本发明所解决的技术问题采用以下技术方案来实现:
一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,首先对在轨滑行、脱架至火箭脱落和上升到安全高度、速度飞行三个阶段的运动机理进行分析,以确定各飞行阶段作用在飞机上的力和力矩,再利用三状态机体运动方程组建立数学模型,并运用Matlab仿真环境以完成数模仿真,具体步骤如下:
1)在轨滑行阶段
①此阶段飞机系统由飞行器、发射架和火箭助推器组成,当发射架受力达到倒架启动力后,发射架牵引带来飞机弧形运动提供向下的俯仰角速度,在数学模型中引入俯仰角速度运动特性,利用三状态机体运动方程组建立数学模型:
其中,Fx为合力在机体X轴的分量,Fy为合力在机体Y轴的分量,Fz为合力在机体Z轴的分量,飞机所受的合力为:发动机推力、火箭推力、导轨摩擦力、前支架支反力、后支架支反力、倒架力、自身重力、升力、阻力;
∑Mx为合力矩在机体X轴的分量,∑My为合力矩在机体Y轴的分量,∑Mz为合力矩在机体Z轴的分量,合力矩由推力矩、摩擦力矩、支反力矩、倒架力矩、气动力矩组成;
ωx为围绕机体X轴转动的角速度,ωy为围绕机体Y轴转动的角速度,ωz为围绕机体Z轴转动的角速度;
仿真过程,根据平尾预偏、发射架安装参数、助推器安装参数等初始输入条件,获得发射起飞的初始高度、姿态角等初始飞行参数;插值发动机推力曲线和火箭推力曲线,获得不同时刻的推力;根据姿态角、舵偏角、马赫数、高度等获得飞机气动力特性;再根据机体Y轴方向力的平衡,计算得到作用在飞机上的前支架支反力与后支架支反力,从而获得飞机受到的摩擦力;
②在Fx达到启动力后,飞机前支点受到发射架前支架的牵引而作围绕前支架底端支点的弧形运动,此时飞机重心位置的俯仰角速度为:
ωz=Vsj×sin(arcsin(L/Hq))/R
其中,Hq为飞机前支点高度,L为飞机的水平飞行距离,R为飞机实际重心到飞机前支点的距离,Vsj为飞机实际飞行速度;
③根据此阶段仿真结束条件,即发射架前支架旋转角度达到倒架角后,结束此阶段仿真,并将此脱轨时刻获得的姿态角、舵偏角、速度等参数作为下一阶段仿真的初始参数输入;
2)脱架至火箭脱落阶段
①飞机脱架,不再受导轨和飞机前支架的作用力,此阶段飞机系统由飞行器和火箭助推器组成,飞机所受的力包括发动机推力、火箭推力、自身重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据此阶段仿真结束条件,即火箭提供的推力在地面坐标系Y轴的分量不足以平衡自身重力为止,完成此阶段仿真;在此阶段火箭脱落时刻,分析此时重量、重心、惯量矩、实际飞机重心距发动机推力线垂直距离等参数的瞬态变化,并与姿态角、舵偏角、高度、速度等参数一起作为下一阶段仿真的初始参数输入;
3)上升达到安全高度、速度飞行阶段
①在此阶段,火箭脱落,飞机系统仅含飞行器,此时飞机受到的力包括发动机推力、自身重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据阶段仿真结束条件,即预先给定的飞行安全高度、速度,完成仿真;
通过对上述三个阶段的仿真,逼真模拟给定初始输入条件下的发射起飞段飞行轨迹、高度、速度等特性,以获得可靠的无人机安全发射包线,为无人机安全发射起飞提供有力技术支撑。
在本发明中,为获得安全发射起飞包线,需进行安全性检查,基于安全性约束条件,以循环方式优化搜索初始输入条件。
有益效果:本发明通过建立系统、完备的数学模型,逼真模拟给定初始输入条件下的发射起飞段飞行轨迹、高度、速度等特性,进而获得可靠的无人机安全发射包线,为无人机安全发射起飞提供有力技术支撑。
附图说明
图1为本发明的较佳实施例的流程示意图。
图2~图3为本发明的较佳实施例中的无人机起飞阶段仿真结果示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
参见图1的一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,首先对在轨滑行、脱架至火箭脱落和上升到安全高度、速度飞行三个阶段的运动机理进行分析,以确定各飞行阶段作用在飞机上的力和力矩,再利用三状态机体运动方程组建立数学模型,并运用Matlab仿真环境以完成数模仿真,具体步骤如下:
1)在轨滑行阶段
①此阶段飞机系统由飞行器、发射架和火箭助推器组成,考虑倒架启动力、平尾预偏、发射架安装参数、助推器安装参数等初始输入条件,通过计算飞机的气动特性,结合飞机系统的质量特性,建立适用于在轨滑行阶段的数学模型;当发射架受力达到倒架启动力后,发射架牵引带来飞机弧形运动,飞机弧形运动本质是发射架牵引飞机围绕前支架底端支点进行的旋转运动,在该运动过程中,提供一个向下的俯仰角速度,俯仰角速度导致飞机低头,对发射起飞不利,只有在数学模型中引入该运动特性,才能提高安全性,获得更为可靠的发射起飞包线;基于飞行器、发射架和助推器组成的系统,分析飞机所受到的力和力矩,利用三状态机体运动方程组建立数学模型:
其中,Fx为合力在机体X轴的分量,Fy为合力在机体Y轴的分量,Fz为合力在机体Z轴的分量,飞机所受的合力为:发动机推力、火箭推力、导轨摩擦力、前支架支反力、后支架支反力、倒架力、自身重力、升力、阻力;
∑Mx为合力矩在机体X轴的分量,∑My为合力矩在机体Y轴的分量,∑Mz为合力矩在机体Z轴的分量,合力矩由推力矩、摩擦力矩、支反力矩、倒架力矩、气动力矩组成;
ωx为围绕机体X轴转动的角速度,ωy为围绕机体Y轴转动的角速度,ωz为围绕机体Z轴转动的角速度;
仿真过程,根据平尾预偏、发射架安装参数、助推器安装参数等初始输入条件,获得发射起飞的初始高度、姿态角等初始飞行参数;插值发动机推力曲线和火箭推力曲线,以获得不同时刻的推力;根据姿态角、舵偏角、马赫数、高度等获得飞机气动力特性;再根据机体Y轴方向力的平衡,计算得到作用在飞机上的前支架支反力与后支架支反力,从而获得飞机受到的摩擦力;
②在Fx达到启动力后,飞机前支点受到发射架前支架的牵引而作围绕前支架底端支点的弧形运动,此时飞机重心位置的俯仰角速度为:
ωz=Vsj×sin(arcsin(L/Hq))/R
其中,Hq为飞机前支点高度,L为飞机的水平飞行距离,R为飞机实际重心到飞机前支点的距离,Vsj为飞机实际飞行速度;
③根据此阶段仿真结束条件,即发射架前支架旋转角度达到倒架角后,结束此阶段仿真,并将此脱轨时刻获得的姿态角、舵偏角、速度等参数存储作为下一阶段仿真的初始参数输入;
2)脱架至火箭脱落阶段
①飞机脱架,不再受导轨和飞机前支架的作用力,此阶段飞机系统由飞行器和火箭助推器组成,飞机所受的力包括发动机推力、火箭推力、重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据此阶段仿真结束条件,即火箭提供的推力在地面坐标系Y轴的分量不足以平衡自身重力为止,完成此阶段仿真;在此阶段火箭脱落时刻,分析此时重量、重心、惯量矩、实际飞机重心距发动机推力线垂直距离等参数的瞬态变化,并与姿态角、舵偏角、高度、速度等参数一起作为下一阶段仿真的初始参数输入;
3)上升达到安全高度、速度飞行阶段
①在此阶段,火箭脱落,飞机系统仅含飞行器,此时飞机受到的力包括发动机推力、重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据阶段仿真结束条件,即预先给定的飞行安全高度、速度,完成仿真;
通过对上述三状态仿真,逼真模拟给定初始输入条件下的发射起飞段飞行轨迹、高度、速度等特性,以获得可靠的无人机安全发射包线,为无人机安全发射起飞提供有力技术支撑。
在本实施例中,为获得安全发射起飞包线,需进行安全性检查,基于安全性约束条件,以循环方式优化搜索初始输入条件,以某型飞机为例,安全性约束条件假定为:发射起飞段不允许掉高,迎角不大于20°,航迹角不小于6.5°,俯仰角不大于30°,仿真时间4s,仿真结果如图2~图3所示,从图中可知,该初始状态在安全发射包线内,能够满足安全发射起飞要求。
以上显示和描述了本发明的基本原理和主要特征以及本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (2)
1.一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,其特征在于,首先对在轨滑行、脱架至火箭脱落和上升到安全高度、速度飞行三个阶段的运动机理进行分析,以确定各飞行阶段作用在飞机上的力和力矩,再利用三状态机体运动方程组建立数学模型,并运用Matlab仿真环境以完成数模仿真,具体步骤如下:
1)在轨滑行阶段
①此阶段飞机系统由飞行器、发射架和火箭助推器组成,当发射架受力达到倒架启动力后,发射架牵引带来飞机弧形运动提供向下的俯仰角速度,在数学模型中引入俯仰角速度运动特性,利用三状态机体运动方程组建立数学模型:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>dv</mi>
<mrow>
<mi>x</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>z</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>dv</mi>
<mrow>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>x</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>z</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>dv</mi>
<mrow>
<mi>z</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>x</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&Sigma;M</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&Sigma;M</mi>
<mi>y</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mi>z</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mi>z</mi>
</msub>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&Sigma;M</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Fx为合力在机体X轴的分量,Fy为合力在机体Y轴的分量,Fz为合力在机体Z轴的分量,飞机所受的合力为:发动机推力、火箭推力、导轨摩擦力、前支架支反力、后支架支反力、倒架力、自身重力、升力、阻力;
∑Mx为合力矩在机体X轴的分量,∑My为合力矩在机体Y轴的分量,∑Mz为合力矩在机体Z轴的分量,合力矩由推力矩、摩擦力矩、支反力矩、倒架力矩、气动力矩组成;
ωx为围绕机体X轴转动的角速度,ωy为围绕机体Y轴转动的角速度,ωz为围绕机体Z轴转动的角速度;
②在Fx达到启动力后,飞机前支点受到发射架前支架的牵引而作围绕前支架底端支点的弧形运动,此时飞机重心位置的俯仰角速度为:
ωz=Vsj×sin(arcsin(L/Hq))/R
其中,Hq为飞机前支点高度,L为飞机的水平飞行距离,R为飞机实际重心到飞机前支点的距离,Vsj为飞机实际飞行速度;
③根据此阶段仿真结束条件,即发射架前支架旋转角度达到倒架角后,结束此阶段仿真,并将此脱轨时刻获得的飞行参数作为下一阶段仿真的初始参数输入;
2)脱架至火箭脱落阶段
①飞机脱架,不再受导轨和飞机前支架的作用力,此阶段飞机系统由飞行器和火箭助推器组成,飞机所受的力包括发动机推力、火箭推力、重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据此阶段仿真结束条件,即火箭提供的推力在地面坐标系Y轴的分量不足以平衡自身重力为止,完成此阶段仿真;在此阶段火箭脱落时刻,分析此时重量、重心、惯量矩、实际飞机重心距发动机推力线垂直距离等参数的瞬态变化,获得火箭脱落时刻的飞行参数,并将此飞行参数作为下一阶段仿真的初始参数输入;
3)上升达到安全高度、速度飞行阶段
①在此阶段,火箭脱落,飞机系统仅含飞行器,此时飞机受到的力包括发动机推力、重力、升力、阻力,引入控制律设计的舵偏,通过计算飞机的气动特性,结合飞机系统的质量特性,进行实时仿真;
②根据阶段仿真结束条件,即预先给定的飞行安全高度、速度,完成仿真;
通过对上述三状态仿真,逼真模拟给定初始输入条件下的发射起飞段飞行轨迹、高度、速度等特性,进而获得无人机安全发射包线。
2.根据权利要求1所述的一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法,其特征在于,为获得安全发射起飞包线,需进行安全性检查,基于安全性约束条件,以循环方式优化搜索初始输入条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711231763.7A CN108108524A (zh) | 2017-11-30 | 2017-11-30 | 一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711231763.7A CN108108524A (zh) | 2017-11-30 | 2017-11-30 | 一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108108524A true CN108108524A (zh) | 2018-06-01 |
Family
ID=62207964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711231763.7A Pending CN108108524A (zh) | 2017-11-30 | 2017-11-30 | 一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108108524A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583041A (zh) * | 2018-11-07 | 2019-04-05 | 湖北航天技术研究院总体设计所 | 一种运载火箭弹道设计方法 |
CN110990947A (zh) * | 2019-11-19 | 2020-04-10 | 中国人民解放军总参谋部第六十研究所 | 一种火箭助推无人机发射过程多场耦合仿真分析方法 |
CN111290299A (zh) * | 2020-03-20 | 2020-06-16 | 湖南云顶智能科技有限公司 | 一种高超声速试飞器助推段控制系统仿真平台 |
CN113094817A (zh) * | 2021-04-20 | 2021-07-09 | 中国直升机设计研究所 | 一种直升机武器发射响应的计算方法 |
CN113120250A (zh) * | 2021-04-27 | 2021-07-16 | 成都飞机工业(集团)有限责任公司 | 一种零长发射无人机发射参数的快速选取方法 |
CN113204196A (zh) * | 2021-04-20 | 2021-08-03 | 西安爱生技术集团公司 | 一种中小型无人机火箭助推发射仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592007A (zh) * | 2011-12-30 | 2012-07-18 | 成都飞机工业(集团)有限责任公司 | 用于飞行控制律设计调参的无人机对象模型建模方法 |
CN106970531A (zh) * | 2017-05-02 | 2017-07-21 | 西北工业大学 | 倾转翼垂直起降无人机模态转换控制策略确定方法 |
-
2017
- 2017-11-30 CN CN201711231763.7A patent/CN108108524A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592007A (zh) * | 2011-12-30 | 2012-07-18 | 成都飞机工业(集团)有限责任公司 | 用于飞行控制律设计调参的无人机对象模型建模方法 |
CN106970531A (zh) * | 2017-05-02 | 2017-07-21 | 西北工业大学 | 倾转翼垂直起降无人机模态转换控制策略确定方法 |
Non-Patent Citations (6)
Title |
---|
刘建新著: "《空天科学与工程系列教材 飞行器总体技术与应用 导弹总体设计导论》", 31 January 2017, 国防工业出版社 * |
刘文珽主编: "《结构可靠性设计手册》", 29 February 2008, 国防工业出版社 * |
姜椿芳总编辑,中国大百科全书总编辑委员会航空、航天编辑委员会编: "《中国大百科全书 航空、航天》", 30 April 1992, 中国大百科全书出版社 * |
方辉熠主编: "《防空导弹武器系统仿真》", 31 August 1995, 宇航出版社 * |
杨广 等: "单发夹角式火箭助推发射起飞仿真技术研究", 《飞行力学》 * |
陶于金: "小型固定翼无人机零长发射参数安全边界研究", 《海军航空工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583041A (zh) * | 2018-11-07 | 2019-04-05 | 湖北航天技术研究院总体设计所 | 一种运载火箭弹道设计方法 |
CN110990947A (zh) * | 2019-11-19 | 2020-04-10 | 中国人民解放军总参谋部第六十研究所 | 一种火箭助推无人机发射过程多场耦合仿真分析方法 |
CN111290299A (zh) * | 2020-03-20 | 2020-06-16 | 湖南云顶智能科技有限公司 | 一种高超声速试飞器助推段控制系统仿真平台 |
CN113094817A (zh) * | 2021-04-20 | 2021-07-09 | 中国直升机设计研究所 | 一种直升机武器发射响应的计算方法 |
CN113204196A (zh) * | 2021-04-20 | 2021-08-03 | 西安爱生技术集团公司 | 一种中小型无人机火箭助推发射仿真方法 |
CN113094817B (zh) * | 2021-04-20 | 2022-09-06 | 中国直升机设计研究所 | 一种直升机武器发射响应的计算方法 |
CN113120250A (zh) * | 2021-04-27 | 2021-07-16 | 成都飞机工业(集团)有限责任公司 | 一种零长发射无人机发射参数的快速选取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108108524A (zh) | 一种单发夹角式火箭助推发射无人机起飞阶段的仿真方法 | |
Selig | Modeling full-envelope aerodynamics of small UAVs in realtime | |
Selig | Real-time flight simulation of highly maneuverable unmanned aerial vehicles | |
CN104777844B (zh) | 一种高超声速临近空间飞行器航迹跟踪方法 | |
CN109614633A (zh) | 一种复合式旋翼飞行器非线性建模方法及配平方法 | |
CN103631153B (zh) | 可视化双模制导航空时敏炸弹半实物仿真方法 | |
CN100576124C (zh) | 一种空天飞机高空爬升纵向控制方法 | |
CN102592007A (zh) | 用于飞行控制律设计调参的无人机对象模型建模方法 | |
CN105242679A (zh) | 一种四旋翼飞行器的控制系统设计方法 | |
Chiappinelli et al. | Modeling and Control of a Tailsitter UAV | |
Young et al. | Design considerations for a Mars highland helicopter | |
CN106570242B (zh) | 低空大动压整体式整流罩高速分离流固耦合仿真方法 | |
CN110162818A (zh) | 伞弹系统弹道计算方法 | |
Xue et al. | Derivation and verification of a similarity law for wind-tunnel free-flight tests of heavy-store separation | |
CN113204196A (zh) | 一种中小型无人机火箭助推发射仿真方法 | |
CN211685678U (zh) | 一种多旋翼无人机实时尾迹的仿真分析系统 | |
CN110826189A (zh) | 一种飞行器缩比模型实验系统的确定方法 | |
Parada | Conceptual and preliminary design of a long endurance electric UAV | |
Husseyin et al. | Design considerations for a stopped-rotor cyclocopter for Venus exploration | |
Park et al. | Hybrid analysis for quadrotor type uav and modified blade element momentum theory considering gust and flight condition | |
CN112441253A (zh) | 一种多旋翼无人机实时尾迹的仿真分析系统及方法 | |
Zaw et al. | Development of mathematical model and stability analysis for UAH | |
Black | High-speed store separation-Correlation between wind-tunnel and flight-test data. | |
Rogers | Aerodynamics—retrospect and prospect | |
Liu et al. | Development of a Modularized Virtual Flight Simulator based on Multiple Discipline Coupled Method |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180601 |
|
RJ01 | Rejection of invention patent application after publication |