CN106021679B - 一种末端机动小过载制导指令的奇异处理方法 - Google Patents
一种末端机动小过载制导指令的奇异处理方法 Download PDFInfo
- Publication number
- CN106021679B CN106021679B CN201610319674.7A CN201610319674A CN106021679B CN 106021679 B CN106021679 B CN 106021679B CN 201610319674 A CN201610319674 A CN 201610319674A CN 106021679 B CN106021679 B CN 106021679B
- Authority
- CN
- China
- Prior art keywords
- guidance
- force
- calculating
- longitudinal
- overload
- 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.)
- Expired - Fee Related
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 9
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种末端机动小过载制导指令的奇异处理方法,该方法包括:根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx;对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;根据奇异处理后的制导力计算制导所需的倾侧角;根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。通过使用本发明所提供的方法,可以解决高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。
Description
技术领域
本发明涉及高超声速飞行器制导技术领域,特别涉及一种末端机动小过载制导指令的奇异处理方法。
背景技术
升力式高超声速飞行器由于升力大,当飞行器接近目标时,该飞行器在末端机动段飞行时极易出现小过载制导指令的需求,此时制导律给出的纵向和横向制导力需求均比较小,可能会出现纵向制导力|Fycx|≤0.1g0,横向制导力|Fzcx|≤0.1g0的情况,此时如果还是按照γcx=tg-1(Fzcx/Fycx)计算倾侧角指令γcx,将会导致倾侧角大幅滚转,甚至到90°的情况,从而对姿态控制和制导尤为不利。
小过载制导指令引起的倾侧实际上是不必要的机动,其产生的侧向加速度并不大,虽然倾侧大,但弹体法向过载小,其对机动飞行的贡献基本可以忽略,不必要的滚转反而对控制***造成不利,甚至影响制导精度。
因此,制导***设计需要考虑消除由于小过载制导指令导致的奇异问题,提高制导的有效性。
发明内容
有鉴于此,本发明提供一种末端机动小过载制导指令的奇异处理方法,从而可以解决高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。
本发明的技术方案具体是这样实现的:
一种末端机动小过载制导指令的奇异处理方法,该方法包括:
根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx;
对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量;g0为地球引力加速度;
根据奇异处理后的制导力计算制导所需的倾侧角;
根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
较佳的,所述Kc的取值大于1.0。
较佳的,使用如下的公式计算得到制导所需的倾侧角
较佳的,使用如下的公式计算得到制导需用的法向力Fn:
如上可见,在本发明中的末端机动小过载制导指令的奇异处理方法中,由于先根据预设的末端机动段制导律和方案,计算得到纵向制导力和横向制导力;然后对制导力进行奇异处理,并根据奇异处理后的制导力计算制导所需的倾侧角,最后可以根据纵向制导力和横向制导力计算制导需用的法向力,因此可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转,从而可以从制导力需求和倾侧机动的效率入手,当弹体法向为小过载的条件下时(即出现横法向小过载需求时),消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,因而可以避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性,完全可避免高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。与现有技术中的采用倾侧角进行直接限幅的方式相比,尤其当飞行器真正需要大倾侧转弯时,采用本发明中的末端机动小过载制导指令的奇异处理方法,将不会因为倾侧角别限制而影响制导精度,也不会因为过载变化而引起指令跳跃。
附图说明
图1为本发明实施例中的末端机动小过载制导指令的奇异处理方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种末端机动小过载制导指令的奇异处理方法,该方法适用于以气动力飞行为主的升力式飞行器。
图1为本发明实施例中的末端机动小过载制导指令的奇异处理方法的流程示意图。如图1所示,本发明实施例中的末端机动小过载制导指令的奇异处理方法包括:
步骤101,根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx。
在本发明的技术方案中,在本步骤中,将首先按照预设的末端机动段制导律和方案,进行纵向制导力Fycx和横向制导力Fzcx需求的计算。
另外,在本发明的技术方案中,末端机动制导段可以采用最优比例导引进行制导。因此,在本发明的一个较佳实施例中,可以按照最优比例导引进行制导计算,得到纵向和横向所需要的制导力,即纵向制导力Fycx和横向制导力Fzcx。其中,可以使用本领域中的常用方法按照最优比例导引进行制导计算,因此在此不再赘述。
步骤102,对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx。
其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量,是一个常数;g0为地球引力加速度,较佳的,在本发明的一个具体实施例中,所述g0取值为9.8055。
在本步骤中,将对制导力进行奇异处理,即当弹体法向为小过载的条件下时(即出现横法向小过载需求时,例如,Fycx<Kc·m·g0),则需要消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,从而避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性。
较佳的,在本发明的一个具体实施例中,所述Kc的取值大于1.0。Kc值的选取原则是能够消除倾侧角指令的奇异性,但又不至于影响制导精度。Kc值过大,虽然能够消除倾侧角指令的奇异性,但无法满足飞行器真正需要的倾侧指令,最终导致制导精度差。
步骤103,根据奇异处理后的制导力计算制导所需的倾侧角。
较佳的,在本发明的一个具体实施例中,所述倾侧角可以使用如下所述的公式计算得到:
其中,为制导所需的倾侧角。
步骤104,根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
在本发明的技术方案中,由于飞行器运动所需的力由弹体提供,因此,较佳的,在本发明的一个具体实施例中,所述制导需用的法向力可以使用如下所述的公式计算得到:
通过上述的步骤101~104,即可得到制导所需的倾侧角和法向力,从而可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转。
综上可知,在本发明中的末端机动小过载制导指令的奇异处理方法中,由于先根据预设的末端机动段制导律和方案,计算得到纵向制导力和横向制导力;然后对制导力进行奇异处理,并根据奇异处理后的制导力计算制导所需的倾侧角,最后可以根据纵向制导力和横向制导力计算制导需用的法向力,因此可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转,从而可以从制导力需求和倾侧机动的效率入手,当弹体法向为小过载的条件下时(即出现横法向小过载需求时),消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,因而可以避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性,完全可避免高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。与现有技术中的采用倾侧角进行直接限幅的方式相比,尤其当飞行器真正需要大倾侧转弯时,采用本发明中的末端机动小过载制导指令的奇异处理方法,将不会因为倾侧角别限制而影响制导精度,也不会因为过载变化而引起指令跳跃。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (4)
1.一种末端机动小过载制导指令的奇异处理方法,其特征在于,该方法包括:
根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx;
对纵向制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量;g0为地球引力加速度;
根据横向制导力Fzcx和限幅后的纵向制导力Fyn计算制导所需的倾侧角;
根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
2.根据权利要求1所述的方法,其特征在于:
所述Kc的取值大于1.0。
3.根据权利要求2所述的方法,其特征在于,使用如下的公式计算得到制导所需的倾侧角
4.根据权利要求3所述的方法,其特征在于,使用如下的公式计算得到制导需用的法向力Fn:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610319674.7A CN106021679B (zh) | 2016-05-13 | 2016-05-13 | 一种末端机动小过载制导指令的奇异处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610319674.7A CN106021679B (zh) | 2016-05-13 | 2016-05-13 | 一种末端机动小过载制导指令的奇异处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021679A CN106021679A (zh) | 2016-10-12 |
CN106021679B true CN106021679B (zh) | 2017-04-19 |
Family
ID=57100806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610319674.7A Expired - Fee Related CN106021679B (zh) | 2016-05-13 | 2016-05-13 | 一种末端机动小过载制导指令的奇异处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021679B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108801081B (zh) * | 2018-06-15 | 2020-12-18 | 上海航天控制技术研究所 | 一种针对btt导弹的非奇异滚转指令生成算法 |
CN110209197B (zh) * | 2019-06-25 | 2022-03-25 | 湖北航天技术研究院总体设计所 | 一种飞行器控制***设计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7014150B2 (en) * | 2004-07-30 | 2006-03-21 | Honeywell International Inc. | Method and system for optimizing torque in a CMG array |
CN104266546B (zh) * | 2014-09-22 | 2015-10-07 | 哈尔滨工业大学 | 一种基于视线的有限时间收敛主动防御制导控制方法 |
CN105022403B (zh) * | 2015-05-11 | 2016-03-23 | 北京航天自动控制研究所 | 滑翔飞行器的纵向轨迹控制增益的确定方法 |
CN105466660B (zh) * | 2015-12-29 | 2017-11-28 | 清华大学 | 一种采用加速度观测器的飞机气动参数分区分步辨识方法 |
-
2016
- 2016-05-13 CN CN201610319674.7A patent/CN106021679B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106021679A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11286065B2 (en) | Method for designing reentry trajectory based on flight path angle planning | |
CN104615144B (zh) | 基于目标规划的高超声速飞行器再入轨迹在线优化方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN106643341B (zh) | 基于准平衡滑翔原理的力热控制耦合设计方法 | |
CN103587681A (zh) | 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法 | |
CN106021679B (zh) | 一种末端机动小过载制导指令的奇异处理方法 | |
Zhu et al. | Impact time and angle control guidance independent of time-to-go prediction | |
CN110471456A (zh) | 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法 | |
CN104085539B (zh) | 成像定标的姿态控制方法 | |
CN104019701B (zh) | 一种利用直接力气动力复合控制的前向拦截制导方法 | |
CN105759827B (zh) | 一种抑制不期望柔性振动的航天器姿态控制*** | |
CN109398763A (zh) | 一种基于有限推力有限工作时间情况下的航天器精确入轨控制方法 | |
CN114485672B (zh) | 小天体附着的面状探测器耦合约束轨迹规划方法 | |
CN111591470A (zh) | 一种适应推力可调模式的飞行器精确软着陆闭环制导方法 | |
CN103486905A (zh) | 一种再入飞行器末制导交班条件确定方法 | |
CN109542111B (zh) | 基于分段式的无人机超低空飞行控制方法 | |
CN109343551A (zh) | 一种旋翼机协调转弯控制方法及*** | |
CN104656659B (zh) | 舰载机滑跃起飞自动飞行控制方法 | |
CN109460049A (zh) | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 | |
CN102880049B (zh) | 一种基于帆板挠性形变测量的自适应振动控制方法 | |
CN105759830B (zh) | 一种升力式飞行器高动态下压段制导方法 | |
CN115615261B (zh) | 大长细比火箭弹弹性辨识与视线角速率提取融合方法 | |
CN110879604A (zh) | 一种带落角控制的飞行器航向导引方法 | |
CN105836160B (zh) | 一种升力式飞行器在稀薄流区飞行的倾侧制导方法 | |
CN105923172B (zh) | 一种升力式飞行器的倾斜转弯翻转制导方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170419 Termination date: 20200513 |
|
CF01 | Termination of patent right due to non-payment of annual fee |