Movatterモバイル変換


[0]ホーム

URL:


CN112013726B - An integrated design method for fully strapdown guidance and control based on a third-order model - Google Patents

An integrated design method for fully strapdown guidance and control based on a third-order model
Download PDF

Info

Publication number
CN112013726B
CN112013726BCN202010866153.XACN202010866153ACN112013726BCN 112013726 BCN112013726 BCN 112013726BCN 202010866153 ACN202010866153 ACN 202010866153ACN 112013726 BCN112013726 BCN 112013726B
Authority
CN
China
Prior art keywords
missile
design
guidance
angle
control
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
Application number
CN202010866153.XA
Other languages
Chinese (zh)
Other versions
CN112013726A (en
Inventor
侯明哲
吴爱国
石文锐
段广仁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology ShenzhenfiledCriticalHarbin Institute of Technology Shenzhen
Priority to CN202010866153.XApriorityCriticalpatent/CN112013726B/en
Publication of CN112013726ApublicationCriticalpatent/CN112013726A/en
Application grantedgrantedCritical
Publication of CN112013726BpublicationCriticalpatent/CN112013726B/en
Expired - Fee Relatedlegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Images

Classifications

Landscapes

Abstract

Translated fromChinese

本发明公开了一种基于三阶模型的全捷联制导控制一体化设计方法,所述方法包括如下步骤:第一步、建立三阶制导控制一体化设计模型;第二步、明确考虑全捷联导引头视场约束的制导控制一体化算法的设计任务;第三步、构造辅助系统,设计第一层期望虚拟控制量ηd,并将其通过近似饱和函数处理后得到第一层虚拟控制量ηc;第四步:利用Barrier Lyapunov函数,设计第二层虚拟控制量ωzc;第五步、设计实际舵偏角指令δz;第六步、综合第三至第五步,得到考虑视场约束的制导控制一体化算法;第七步、检验制导控制一体化算法的性能。本发明的方法能够实现对目标的精确打击,并确保全捷联导引头视场约束得以满足。

Figure 202010866153

The invention discloses a fully strap-down guidance and control integrated design method based on a third-order model. The method includes the following steps: first, establishing a third-order guidance and control integrated design model; The design task of the guidance and control integration algorithm with the constraints of the combined seeker's field of view; the third step is to construct the auxiliary system, design the first layer of the desired virtual control variable ηd , and process it through the approximate saturation function to obtain the first layer of virtual control. Control variable ηc ; fourth step: use the Barrier Lyapunov function to design the second-layer virtual control variable ωzc ; fifth step, design the actual rudder deflection angle command δz ; sixth step, integrate the third to fifth steps, get The guidance and control integration algorithm considering the field of view constraints; the seventh step, check the performance of the guidance and control integration algorithm. The method of the invention can achieve precise strike on the target and ensure that the constraints of the field of view of the full strapdown seeker are satisfied.

Figure 202010866153

Description

Translated fromChinese
一种基于三阶模型的全捷联制导控制一体化设计方法An integrated design method for fully strapdown guidance and control based on a third-order model

技术领域technical field

本发明属于航空航天领域,涉及一种制导控制一体化设计方法,具体涉及一种基于三阶模型的考虑全捷联导引头视场约束的制导控制一体化设计方法。The invention belongs to the field of aerospace, and relates to a guidance and control integrated design method, in particular to a guidance and control integrated design method based on a third-order model considering the constraints of a full strapdown seeker field of view.

背景技术Background technique

开展全捷联探测制导体制下的制导控制一体化设计研究对发展低成本、高性能战术导弹具有十分重要的意义。全捷联制导控制一体化设计中的一个重要问题是处理导引头视场约束。从所建立的制导控制一体化设计模型角度看,目前该领域的主要进展包括四阶模型方法和二阶模型方法。四阶模型方法以文献1“考虑全捷联导引头视场约束的制导控制一体化设计方法,中国,2017-01-13,CN201710023831.4”和文献2“考虑全捷联导引头视场约束的三维制导控制一体化设计方法,中国,2019-01-29,CN2019100834420.3”为代表。这类方法以体视线角、视线角速率、攻角(侧滑角、滚转角)和俯仰(偏航、滚转)角速率为作为四阶模型的状态变量。这类方法主要存在两个问题:一是将俯仰(偏航)角速率这些状态变量当作有界干扰处理,缺乏严格的系统稳定性保障;二是模型阶次较高,导致所设计的一体化算法相对复杂。二阶模型方法以文献2“Integrated guidance and control formissile with narrow field-of-view strapdown seeker,ISA Transactions,2020,https://doi.org/10.1016/j.isatra.2020.06.012”为代表。该文献建立了以体视线角及其变化率为状态的二阶制导控制一体化设计模型,通过设计俯仰舵偏角指令使体视线角变化率跟踪负的导弹俯仰角速率信号,从而实现对目标的精确打击。但是,由于导弹俯仰角速率信号并不是一个独立的外部信号,而是直接受俯仰舵偏角指令所控制的,因此,这种设计方法本质上存在控制器的循环设计问题。It is of great significance to carry out the research on the integrated design of guidance and control under the full strapdown detection guidance system for the development of low-cost, high-performance tactical missiles. An important issue in the integrated design of full strapdown guidance and control is to deal with the constraints of the seeker field of view. From the perspective of the established guidance and control integrated design model, the main progress in this field includes the fourth-order model method and the second-order model method. The fourth-order model method is based on document 1 "Guidance and control integrated design method considering the full strapdown seeker field of view constraints, China, 2017-01-13, CN201710023831.4" anddocument 2 "Considering the full strapdown seeker view Field Constrained 3D Guidance and Control Integrated Design Method, China, 2019-01-29, CN2019100834420.3" as the representative. This kind of method takes the body line of sight angle, line of sight angular rate, angle of attack (sideslip angle, roll angle) and pitch (yaw, roll) angular rate as the state variables of the fourth-order model. There are two main problems with this type of method: one is that the state variables such as pitch (yaw) angular rate are treated as bounded disturbances, and there is no strict system stability guarantee; The algorithm is relatively complex. The second-order model method is represented bydocument 2 "Integrated guidance and control formissile with narrow field-of-view strapdown seeker, ISA Transactions, 2020, https://doi.org/10.1016/j.isatra.2020.06.012". This document establishes a second-order guidance and control integrated design model based on the body line of sight angle and its rate of change. By designing the pitch rudder declination command, the rate of change of the body line of sight angle tracks the negative missile pitch angle rate signal, so as to realize the target detection. precision strike. However, since the missile pitch rate signal is not an independent external signal, but is directly controlled by the pitch angle command, this design method essentially has the problem of loop design of the controller.

发明内容SUMMARY OF THE INVENTION

为了克服现有基于四阶模型和二阶模型的考虑全捷联导引头视场约束的制导控制一体化设计方法的不足,本发明提供了一种基于三阶模型的全捷联制导控制一体化设计方法。该方法能够实现对目标的精确打击,并确保全捷联导引头视场约束得以满足。In order to overcome the shortcomings of the existing guidance and control integration design methods based on the fourth-order model and the second-order model considering the constraints of the full strapdown seeker field of view, the present invention provides a third-order model-based full strapdown guidance and control integration. Design method. This method can achieve precise strikes on targets and ensure that the constraints of the full strapdown seeker field of view are satisfied.

本发明的目的是通过以下技术方案实现的:The purpose of this invention is to realize through the following technical solutions:

一种基于三阶模型的全捷联制导控制一体化设计方法,包括如下步骤:An integrated design method for full strapdown guidance and control based on a third-order model, comprising the following steps:

第一步、建立三阶制导控制一体化设计模型:The first step is to establish a third-order guidance and control integrated design model:

Figure BDA0002649808600000021
Figure BDA0002649808600000021

式中,

Figure BDA0002649808600000022
Figure BDA0002649808600000023
λ=θ-ε为导弹速度追踪误差角,θ为导弹弹道倾角,ε为视线倾角,
Figure BDA0002649808600000024
为导弹攻角,
Figure BDA0002649808600000025
为导弹俯仰角,ωz为导弹俯仰角速率,δz为导弹俯仰舵偏角,m为导弹质量,V为导弹速度,P为导弹发动机推力,q为导弹动压头,S为导弹特征面积,g为重力加速度,
Figure BDA0002649808600000031
为导弹升力系数cy对α的偏导数,L为导弹特征长度,Jz为导弹绕弹体z轴的转动惯量,
Figure BDA0002649808600000032
分别为导弹俯仰力矩系数mz相对于α,δz的偏导数;In the formula,
Figure BDA0002649808600000022
Figure BDA0002649808600000023
λ=θ-ε is the missile speed tracking error angle, θ is the missile ballistic inclination angle, ε is the line of sight inclination,
Figure BDA0002649808600000024
is the missile angle of attack,
Figure BDA0002649808600000025
is the missile pitch angle, ωz is the missile pitch angle rate, δz is the missile pitch rudder declination angle, m is the missile mass, V is the missile speed, P is the missile engine thrust, q is the missile dynamic pressure head, and S is the missile characteristic area , g is the acceleration of gravity,
Figure BDA0002649808600000031
is the partial derivative of the missile lift coefficient cy to α, L is the characteristic length of the missile, Jz is the moment of inertia of the missile around the z-axis of the missile body,
Figure BDA0002649808600000032
are the partial derivatives of the missile pitching moment coefficient mz with respect to α, δz , respectively;

第二步、明确考虑全捷联导引头视场约束的制导控制一体化算法的设计任务,具体要求如下:The second step is to clearly consider the design task of the guidance and control integration algorithm considering the constraints of the full strapdown seeker field of view. The specific requirements are as follows:

根据第一步建立的三阶制导控制一体化设计模型,设计导弹的舵偏角指令使导弹速度追踪误差角λ尽快收敛到零,同时保证全捷联导引头视场约束得以满足,即满足:

Figure BDA0002649808600000033
其中,
Figure BDA0002649808600000034
表示体视线角的最大允许值,η=λ+α表示导弹体视线角;According to the third-order guidance and control integrated design model established in the first step, the rudder deflection angle command of the missile is designed to make the missile speed tracking error angle λ converge to zero as soon as possible, and at the same time to ensure that the full strapdown seeker field of view constraints are satisfied, that is :
Figure BDA0002649808600000033
in,
Figure BDA0002649808600000034
Represents the maximum allowable value of the body line of sight angle, η=λ+α represents the missile body line of sight angle;

第三步、构造辅助系统,设计第一层期望虚拟控制量ηd,并将其通过近似饱和函数处理后得到第一层虚拟控制量ηcThe third step is to construct an auxiliary system, design the desired virtual control variable ηd of the first layer, and process it through the approximate saturation function to obtain the virtual control variable ηc of the first layer;

第四步、利用BarrierLyapunov函数,设计第二层虚拟控制量ωzcThe fourth step is to use the BarrierLyapunov function to design the second-layer virtual control quantity ωzc ;

第五步、设计实际舵偏角指令δzThe fifth step, design the actual rudder deflection angle command δz ;

第六步、综合第三至第五步,得到如下考虑全捷联导引头视场约束的制导控制一体化算法:The sixth step, synthesizing the third to fifth steps, the following integrated algorithm of guidance and control considering the constraints of the full strapdown seeker field of view is obtained:

Figure BDA0002649808600000041
Figure BDA0002649808600000041

其中,设计参数的取值范围为:ke>0,ε>0,

Figure BDA0002649808600000042
ki>0,i=1,2,3;Among them, the value range of the design parameters is:ke > 0, ε > 0,
Figure BDA0002649808600000042
ki >0, i=1, 2, 3;

第七步、检验制导控制一体化算法的性能:The seventh step, check the performance of the guidance and control integration algorithm:

在允许的范围内选择好设计参数后,借助计算机数值计算和仿真软件进行仿真计算并进行性能检验,如果制导控制一体化算法的性能满足要求,则设计结束;否则,需要调整设计参数,重新进行仿真计算并进行性能检验。After selecting the design parameters within the allowable range, use computer numerical calculation and simulation software to carry out simulation calculation and perform performance inspection. If the performance of the integrated guidance and control algorithm meets the requirements, the design ends; otherwise, the design parameters need to be adjusted and re-run Simulation calculation and performance verification.

相比于现有技术,本发明具有如下优点:Compared with the prior art, the present invention has the following advantages:

本发明既克服了现有基于四阶模型的考虑全捷联导引头视场约束的制导控制一体化设计方法由于模型阶次较高导致的控制器复杂性问题以及系统的稳定性问题,同时也避免了基于二阶模型的考虑全捷联导引头视场约束的制导控制一体化设计方法存在的控制器循环设计问题。The present invention not only overcomes the complexity of the controller and the stability of the system caused by the higher model order of the existing guidance and control integrated design method based on the fourth-order model considering the constraints of the full strapdown seeker field of view, and at the same time It also avoids the controller loop design problem existing in the guidance and control integrated design method based on the second-order model considering the full strapdown seeker field of view constraints.

附图说明Description of drawings

图1为本发明基于三阶模型的全捷联制导控制一体化设计流程图;Fig. 1 is the integrated design flow chart of the fully strapdown guidance control based on the third-order model of the present invention;

图2为纵向平面拦截几何关系;Figure 2 shows the geometric relationship of longitudinal plane interception;

图3为弹-目相对距离变化曲线;Figure 3 is the change curve of the relative distance between the projectile and the eye;

图4为导弹速度追踪误差角变化曲线;Figure 4 is the change curve of the missile velocity tracking error angle;

图5为导弹体视线角变化曲线;Figure 5 is the change curve of the line-of-sight angle of the missile body;

图6为导弹俯仰舵偏角变化曲线。Figure 6 shows the change curve of the pitch angle of the missile.

具体实施方式Detailed ways

下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。The technical solutions of the present invention will be further described below in conjunction with the accompanying drawings, but are not limited thereto. Any modification or equivalent replacement of the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention shall be included in the present invention. within the scope of protection.

本发明提供了一种基于三阶模型的全捷联制导控制一体化设计方法,如图1所示,其设计步骤如下:The present invention provides a fully strap-down guidance and control integrated design method based on a third-order model, as shown in FIG. 1 , and the design steps are as follows:

第一步:建立三阶制导控制一体化设计模型。The first step is to establish a third-order guidance and control integrated design model.

纵向平面内的拦截几何关系如图2所示,其中,M表示导弹,T表示目标,LOS表示视线,ε表示视线倾角,r表示弹-目相对距离,V表示导弹的速度,xb表示导弹弹体纵轴,

Figure BDA0002649808600000051
表示导弹俯仰角,θ表示导弹弹道倾角,
Figure BDA0002649808600000052
表示导弹攻角,λ=θ-ε表示导弹速度追踪误差角,η=λ+α表示导弹体视线角,Vt表示目标速度,θt表示目标弹道倾角。The interception geometric relationship in the longitudinal plane is shown in Figure 2, where M represents the missile, T represents the target, LOS represents the line of sight, ε represents the line of sight inclination, r represents the relative distance between the missile and the target, V represents the speed of the missile, and xb represents the missile projectile longitudinal axis,
Figure BDA0002649808600000051
represents the missile pitch angle, θ represents the missile ballistic inclination angle,
Figure BDA0002649808600000052
represents the missile attack angle, λ=θ-ε represents the missile velocity tracking error angle, η=λ+α represents the missile body sight angle, Vt represents the target speed, and θt represents the target ballistic inclination angle.

根据速度追踪导引原理,为了实现对目标的精确打击,只需使导弹速度方向指向目标飞行,也就是使导弹速度追踪误差角λ尽快收敛到零。因此,取λ作为第一个状态变量,有:According to the principle of speed tracking and guidance, in order to achieve a precise attack on the target, it is only necessary to make the missile speed direction point to the target flight, that is, to make the missile speed tracking error angle λ converge to zero as soon as possible. Therefore, taking λ as the first state variable, we have:

Figure BDA0002649808600000061
Figure BDA0002649808600000061

导弹的弹道倾角θ满足如下动态方程:The ballistic inclination angle θ of the missile satisfies the following dynamic equation:

Figure BDA0002649808600000062
Figure BDA0002649808600000062

式中,m为导弹质量,P为导弹发动机推力,g为重力加速度,Y为导弹所受升力,

Figure BDA0002649808600000063
其中,q为动压头(q=0.5ρV2,ρ为导弹所处高度的空气密度),S为导弹特征面积,δz为导弹俯仰舵偏角,
Figure BDA0002649808600000064
分别为导弹升力系数cy对α,δz的偏导数。因为导弹所受升力主要由弹体产生,即
Figure BDA0002649808600000065
同时,当攻角不大时,有sinα≈α,故可得:In the formula, m is the mass of the missile, P is the thrust of the missile engine, g is the acceleration of gravity, Y is the lift of the missile,
Figure BDA0002649808600000063
Among them, q is the dynamic pressure head (q=0.5ρV2 , ρ is the air density at the height of the missile), S is the characteristic area of the missile, δz is the pitch angle of the missile,
Figure BDA0002649808600000064
are the partial derivatives of the missile lift coefficient cy to α, δz , respectively. Because the lift force on the missile is mainly generated by the projectile body, namely
Figure BDA0002649808600000065
At the same time, when the angle of attack is not large, there is sinα≈α, so it can be obtained:

Figure BDA0002649808600000066
Figure BDA0002649808600000066

选取导弹的攻角α与俯仰角速率ωz作为状态变量,可以建立导弹的控制系统模型如下:Selecting the attack angle α and the pitch rate ωz of the missile as the state variables, the control system model of the missile can be established as follows:

Figure BDA0002649808600000067
Figure BDA0002649808600000067

Figure BDA0002649808600000068
Figure BDA0002649808600000068

其中,L表示导弹特征长度,Jz表示导弹绕弹体z轴的转动惯量,

Figure BDA0002649808600000069
分别表示导弹俯仰力矩系数mz相对于α,δz的偏导数。Among them, L represents the characteristic length of the missile, Jz represents the moment of inertia of the missile around the z-axis of the missile body,
Figure BDA0002649808600000069
respectively represent the partial derivatives of the missile pitching moment coefficient mz with respect to α, δz .

定义

Figure BDA00026498086000000610
Figure BDA00026498086000000611
则有:definition
Figure BDA00026498086000000610
Figure BDA00026498086000000611
Then there are:

Figure BDA0002649808600000071
Figure BDA0002649808600000071

式(6)即为三阶制导控制一体化设计模型。Equation (6) is the integrated design model of the third-order guidance and control.

第二步:明确考虑全捷联导引头视场约束的制导控制一体化算法的设计任务。The second step: the design task of the guidance and control integration algorithm considering the full strapdown seeker field of view constraints is clearly defined.

考虑全捷联导引头视场约束的制导控制一体化算法的设计任务可以描述为:根据制导控制一体化设计模型(6),设计导弹的舵偏角指令使导弹速度追踪误差角λ尽快收敛到零,同时保证全捷联导引头视场约束得以满足,即满足:

Figure BDA0002649808600000072
其中
Figure BDA0002649808600000073
表示体视线角的最大允许值。The design task of the guidance and control integration algorithm considering the full strapdown seeker field of view constraints can be described as follows: According to the guidance and control integration design model (6), design the missile's rudder deflection angle command to make the missile speed tracking error angle λ converge as soon as possible to zero, while ensuring that the full strapdown seeker field of view constraints are satisfied, that is:
Figure BDA0002649808600000072
in
Figure BDA0002649808600000073
Indicates the maximum allowable value of the stereoscopic line-of-sight angle.

第三步:构造辅助系统,设计第一层期望虚拟控制量,并将其通过近似饱和函数处理后得到第一层虚拟控制量。The third step: construct the auxiliary system, design the desired virtual control variable of the first layer, and process it through the approximate saturation function to obtain the virtual control variable of the first layer.

构造辅助系统:Construction assistance system:

Figure BDA0002649808600000074
Figure BDA0002649808600000074

其中,ke>0为设计参数,ηd为待求解的第一层期望虚拟控制量,ηc为ηd通过如下饱和近似函数后的输出,即第一层虚拟控制量:Among them,ke > 0 is the design parameter, ηd is the expected virtual control variable of the first layer to be solved, and ηc is the output of ηd after passing through the following saturation approximation function, that is, the virtual control variable of the first layer:

Figure BDA0002649808600000075
Figure BDA0002649808600000075

其中,设计参数ε>0,

Figure BDA0002649808600000076
满足:
Figure BDA0002649808600000077
由上式可知,
Figure BDA0002649808600000078
Among them, the design parameter ε>0,
Figure BDA0002649808600000076
Satisfy:
Figure BDA0002649808600000077
From the above formula, it can be seen that
Figure BDA0002649808600000078

定义z1=λ-e,则有:Define z1 =λ-e, then we have:

Figure BDA0002649808600000081
Figure BDA0002649808600000081

构造第一层期望虚拟控制量为:The expected virtual control quantity of the first layer is constructed as:

Figure BDA0002649808600000082
Figure BDA0002649808600000082

其中,k1>0为设计参数;并定义z2=η-ηc,则有:Among them, k1 >0 is the design parameter; and z2 =η-ηc is defined, there are:

Figure BDA0002649808600000083
Figure BDA0002649808600000083

定义

Figure BDA0002649808600000084
则有:definition
Figure BDA0002649808600000084
Then there are:

Figure BDA0002649808600000085
Figure BDA0002649808600000085

第四步:利用BarrierLyapunov函数,设计第二层虚拟控制量。Step 4: Use the BarrierLyapunov function to design the second layer of virtual control quantities.

考虑consider

Figure BDA0002649808600000086
Figure BDA0002649808600000086

其中,ωzc为待求的第二层虚拟控制量。Among them, ωzc is the virtual control quantity of the second layer to be obtained.

构造第二层虚拟控制量为:The construction of the second-level virtual control quantity is:

Figure BDA0002649808600000087
Figure BDA0002649808600000087

其中,k2>0为设计参数;并定义z3=ωzzc,则有:Among them, k2 >0 is a design parameter; and z3zzc is defined, there are:

Figure BDA0002649808600000088
Figure BDA0002649808600000088

定义

Figure BDA0002649808600000089
则有:definition
Figure BDA0002649808600000089
Then there are:

Figure BDA00026498086000000810
Figure BDA00026498086000000810

第五步:设计实际舵偏角指令。Step 5: Design the actual rudder declination command.

考虑consider

Figure BDA0002649808600000091
Figure BDA0002649808600000091

构造舵偏角控制指令为:The rudder deflection angle control command is constructed as:

Figure BDA0002649808600000092
Figure BDA0002649808600000092

其中,k3>0为设计参数,则有:Among them, k3 >0 is the design parameter, then:

Figure BDA0002649808600000093
Figure BDA0002649808600000093

定义Lyapunov函数为

Figure BDA0002649808600000094
则有:Define the Lyapunov function as
Figure BDA0002649808600000094
Then there are:

Figure BDA0002649808600000095
Figure BDA0002649808600000095

由上式可知,若初值z2(0)<kb,则有z2<kb且zi,i=1,2,3渐近收敛到零。因而有:

Figure BDA0002649808600000096
即导引头视场约束得以满足。进一步地,当ηc的生成环节严格退饱和后,e将渐近收敛到零,于是,λ也将渐近收敛到零。It can be seen from the above formula that if the initial value z2 (0)<kb , then z2 <kb andzi ,i=1,2,3 asymptotically converge to zero. Hence:
Figure BDA0002649808600000096
That is, the seeker field of view constraints are satisfied. Further, when the generation link of ηc is strictly desaturated, e will converge to zero asymptotically, so λ will also converge to zero asymptotically.

第六步:综合前述第三至第五步,便得到考虑全捷联导引头视场约束的制导控制一体化算法。最终的制导控制一体化算法如下所示:Step 6: Combining the above-mentioned third to fifth steps, an integrated guidance and control algorithm considering the constraints of the full strapdown seeker field of view is obtained. The final guidance and control integration algorithm is as follows:

Figure BDA0002649808600000101
Figure BDA0002649808600000101

其中,设计参数的取值范围为:ke>0,ε>0,

Figure BDA0002649808600000102
ki>0,i=1,2,3,设计参数的具体取值需要结合具体应用场景开展非线性数值仿真进行。Among them, the value range of the design parameters is:ke > 0, ε > 0,
Figure BDA0002649808600000102
ki > 0, i=1, 2, 3, the specific values of design parameters need to be carried out in combination with specific application scenarios to carry out nonlinear numerical simulation.

第七步:检验制导控制一体化算法的性能。Step 7: Check the performance of the guidance and control integration algorithm.

为了检验所设计的考虑全捷联导引头视场约束的制导控制一体化算法的性能,需将其应用于导弹在纵向平面内的非线性制导与控制系统,借助常用的计算机数值计算和仿真软件如Matlab/Simulink等来进行。在允许的范围内选择好设计参数后,进行仿真计算并进行性能检验。如果制导控制一体化算法的性能满足要求,则设计结束;否则,需要调整设计参数,重新进行仿真计算并进行性能检验。In order to test the performance of the guidance and control integration algorithm designed considering the constraints of the full strapdown seeker field of view, it needs to be applied to the nonlinear guidance and control system of the missile in the longitudinal plane. With the help of commonly used computer numerical calculation and simulation Software such as Matlab/Simulink etc. After selecting the design parameters within the allowable range, carry out simulation calculation and perform performance inspection. If the performance of the guidance and control integration algorithm meets the requirements, the design is over; otherwise, the design parameters need to be adjusted, and the simulation calculation and performance inspection should be performed again.

实施例:Example:

这里通过介绍一个具有一定代表性的实施例,来进一步说明本发明技术方案中的相关设计。Here, a representative embodiment is introduced to further illustrate the related designs in the technical solution of the present invention.

将所设计的考虑全捷联导引头视场约束的制导控制一体化算法应用于如下所示的导弹在纵向平面内的非线性制导与控制系统:The designed guidance and control integration algorithm considering the constraints of the full strapdown seeker field of view is applied to the nonlinear guidance and control system of the missile in the longitudinal plane as shown below:

Figure BDA0002649808600000111
Figure BDA0002649808600000111

式中,导弹所受阻力X,升力Y和俯仰力矩Mz的计算公式为:

Figure BDA0002649808600000112
其中,cx0为零升阻力系数,
Figure BDA0002649808600000113
分别为阻力系数相对于α,δz的偏导数。In the formula, the calculation formulas of the resistance X, lift Y and pitching moment Mz on the missile are:
Figure BDA0002649808600000112
where cx0 is zero-lift drag coefficient,
Figure BDA0002649808600000113
are the partial derivatives of the drag coefficient with respect to α, δz , respectively.

在仿真中,设导弹的结构和气动参数分别为:S=0.42m2,L=0.68m,m=1200Kg,Jz=5600Kg·m2,P=5000N,

Figure BDA0002649808600000114
Figure BDA0002649808600000115
设目标为地面慢速运动目标,速度大小为Vt=10m/s,目标弹道倾角为θt=0°,设导弹的初始速度为V=250m/s,导弹的俯仰角初始值为
Figure BDA0002649808600000116
导弹的俯仰角速度初始值为ωz=6°/s,导弹的弹道倾角初始值为θ0=-27°,弹-目相对距离初始值为R0=3000m,视线倾角初始值为ε0=-30°,视场约束设为
Figure BDA0002649808600000117
设计参数选取如下:ke=2,ε=0.1,
Figure BDA0002649808600000118
k1=1,k2=1,k3=1。考虑到实际物理限制,设导弹的最大允许舵偏角为30°。同时,设弹-目相对距离小于50m时,导弹导引头进入盲区,此后导弹舵偏角保持不变,进入无控飞行状态,直至仿真结束。弹-目相对距离小于1m时,停止仿真。In the simulation, the structure and aerodynamic parameters of the missile are set as: S=0.42m2 , L=0.68m, m=1200Kg, Jz =5600Kg·m2 , P=5000N,
Figure BDA0002649808600000114
Figure BDA0002649808600000115
Let the target be the ground slow moving target, the speed is Vt = 10m/s, the target ballistic inclination angle is θt = 0°, the initial speed of the missile is V = 250m/s, the initial value of the pitch angle of the missile is
Figure BDA0002649808600000116
The initial value of the pitch angular velocity of the missile is ωz =6°/s, the initial value of the missile's ballistic inclination angle is θ0 =-27°, the initial value of the relative distance between the missile and the target is R0 =3000m, and the initial value of the line of sight inclination is ε0 = -30°, the field of view constraint is set to
Figure BDA0002649808600000117
The design parameters are selected as follows:ke = 2, ε = 0.1,
Figure BDA0002649808600000118
k1 =1, k2 =1, k3 =1. Considering the actual physical limitations, the maximum allowable rudder deflection angle of the missile is set to 30°. At the same time, when the relative distance between the missile and the target is less than 50m, the missile seeker enters the blind area, and then the missile rudder declination angle remains unchanged and enters the uncontrolled flight state until the end of the simulation. When the relative distance between the projectile and the target is less than 1m, the simulation is stopped.

弹-目相对距离变化曲线如图3所示,导弹的脱靶量小于1m,导弹能够精确命中目标。在导弹的受控飞行阶段(即弹-目相对距离不小于50m阶段),导弹的速度追踪误差角变化曲线如图4所示,导弹的速度追踪误差角λ逐渐收敛并保持在零附近。导弹的体视线角η变化曲线如图5所示,体视线角满足|η|≤20°,即导引头视场约束得到满足。导弹的俯仰舵偏角变化曲线如图6所示。仿真结果表明了所提的考虑全捷联导引头视场约束的制导控制一体化算法的有效性。The change curve of the relative distance between the missile and the target is shown in Figure 3. The missile misses less than 1m, and the missile can accurately hit the target. In the controlled flight stage of the missile (that is, the relative distance between the missile and the target is not less than 50m), the change curve of the missile's velocity tracking error angle is shown in Figure 4. The missile's velocity tracking error angle λ gradually converges and remains near zero. The change curve of the body line of sight angle η of the missile is shown in Figure 5. The body line of sight angle satisfies |η|≤20°, that is, the seeker field of view constraint is satisfied. The change curve of the pitch rudder deflection angle of the missile is shown in Figure 6. The simulation results show the effectiveness of the proposed integrated guidance and control algorithm considering the full strapdown seeker field of view constraints.

Claims (1)

1. A three-order model-based full strapdown guidance control integrated design method is characterized by comprising the following steps:
firstly, establishing a three-order guidance control integrated design model:
Figure FDA0002950855270000011
in the formula,
Figure FDA0002950855270000012
Figure FDA0002950855270000013
lambda is theta-epsilon as missile speed tracking error angle, theta is missile trajectory inclination angle, epsilon is sight inclination angle,
Figure FDA0002950855270000014
is a missile angle of attack,
Figure FDA0002950855270000015
is missile pitch angle, omegazFor missile pitch angle rate, deltazIs missile pitching rudder deflection angle, m is missile mass, V is missile speed, P is missile engine thrust, q is missile dynamic pressure head, S is missile characteristic area, g is gravity acceleration,
Figure FDA0002950855270000016
coefficient of lift of missileyPartial derivative of alpha, L being the missile characteristic length, JzIs the rotational inertia of the missile around the z-axis of the missile,
Figure FDA0002950855270000017
respectively the missile pitching moment coefficient mzRelative to alpha, deltazPartial derivatives of (d);
and step two, specifically considering the design task of a guidance control integrated algorithm of the full strapdown seeker view field constraint, wherein the specific requirements are as follows:
according to the three-order guidance control integrated design model established in the first stepThe missile speed tracking error angle lambda is converged to zero as soon as possible by counting the rudder deflection angle instruction of the missile, and meanwhile, the requirement that the field of view constraint of a full strapdown seeker is met:
Figure FDA0002950855270000018
wherein,
Figure FDA0002950855270000019
represents the maximum allowable value of the body line-of-sight angle, and η ═ λ + α represents the projectile line-of-sight angle;
thirdly, constructing an auxiliary system, designing a first layer expected virtual control quantity etadAnd processing the control signal by an approximate saturation function to obtain a first layer of virtual control quantity etacThe auxiliary system is as follows:
Figure FDA0002950855270000021
wherein k iseGreater than 0 as a design parameter, ηdExpectation of a virtual control quantity, eta, for the first layer to be solvedcIs etadAnd the output after the following saturation approximation function is carried out, namely the first layer virtual control quantity:
Figure FDA0002950855270000022
wherein the design parameter epsilon is more than 0,
Figure FDA0002950855270000023
satisfies the following conditions:
Figure FDA0002950855270000024
fourthly, designing a second layer virtual control quantity omega by utilizing a Barrier Lyapunov functionzc
Fifthly, designing an actual rudder deflection angle instruction deltaz
And sixthly, integrating the third step to the fifth step to obtain a guidance control integrated algorithm considering the field restriction of the full strapdown seeker:
Figure FDA0002950855270000025
wherein, the value range of the design parameter is as follows: k is a radical ofe>0,ε>0,
Figure FDA0002950855270000026
ki>0,i=1,2,3;
Seventhly, checking the performance of the guidance control integrated algorithm: after the design parameters are selected in the allowed range, carrying out simulation calculation and performance inspection by means of computer numerical calculation and simulation software, and finishing the design if the performance of the guidance control integrated algorithm meets the requirements; otherwise, the design parameters need to be adjusted, and the simulation calculation is carried out again and the performance inspection is carried out.
CN202010866153.XA2020-08-252020-08-25 An integrated design method for fully strapdown guidance and control based on a third-order modelExpired - Fee RelatedCN112013726B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN202010866153.XACN112013726B (en)2020-08-252020-08-25 An integrated design method for fully strapdown guidance and control based on a third-order model

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN202010866153.XACN112013726B (en)2020-08-252020-08-25 An integrated design method for fully strapdown guidance and control based on a third-order model

Publications (2)

Publication NumberPublication Date
CN112013726A CN112013726A (en)2020-12-01
CN112013726Btrue CN112013726B (en)2021-05-18

Family

ID=73502258

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN202010866153.XAExpired - Fee RelatedCN112013726B (en)2020-08-252020-08-25 An integrated design method for fully strapdown guidance and control based on a third-order model

Country Status (1)

CountryLink
CN (1)CN112013726B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN114636352B (en)*2022-02-222023-08-04西安航天三沃机电设备有限责任公司Strapdown laser guide head angle output calibration method
CN115421501B (en)*2022-09-142023-04-25中国科学院声学研究所 Method and system for finite-time signal estimation of underwater vehicle
CN115406312B (en)*2022-09-232023-06-23北京航空航天大学 An Integrated Method of Missile Guidance and Control Considering Field of View and Steering Gear Delay Constraints
CN116483109B (en)*2023-06-252023-09-05西安羚控电子科技有限公司Sliding mode control-based integrated aircraft guidance control method and system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN106681348B (en)*2017-01-132019-04-19西北工业大学 Guidance and control integrated design method considering full strapdown seeker field of view constraints
CN108106597B (en)*2017-11-302020-07-07中国人民解放军国防科技大学Method for measuring angle of full strapdown laser seeker under condition of target linear field of view
CN109254594A (en)*2018-09-212019-01-22北京大翔航空科技有限公司A kind of full strapdown terminal guidance method for unmanned plane
CN109581892B (en)*2019-01-292021-10-12西北工业大学Double-turntable semi-physical simulation system and method for full strapdown missile guidance control system
CN110008502B (en)*2019-01-292022-07-15西北工业大学 An integrated design method for 3D guidance and control considering the full strapdown seeker field of view constraints

Also Published As

Publication numberPublication date
CN112013726A (en)2020-12-01

Similar Documents

PublicationPublication DateTitle
CN112013726B (en) An integrated design method for fully strapdown guidance and control based on a third-order model
CN110008502B (en) An integrated design method for 3D guidance and control considering the full strapdown seeker field of view constraints
CN106681348B (en) Guidance and control integrated design method considering full strapdown seeker field of view constraints
CN102706217B (en)Method for controlling attack angle and attack time of multiple missiles
CN113758380B (en)Guidance control integration method considering steering engine dynamic performance and attack angle constraint
CN115406312B (en) An Integrated Method of Missile Guidance and Control Considering Field of View and Steering Gear Delay Constraints
Bao et al.Integrated guidance and control for hypersonic morphing missile based on variable span auxiliary control
CN114153143B (en)Design method of non-singular fixed time sliding mode guidance law of missile
CN112486023B (en)Simulation design method for flight control system of cruise missile flight path and control system
CN106091816B (en)A kind of half strapdown air-to-air missile method of guidance based on sliding mode variable structure theory
CN116839429B (en) An integrated method of guidance and control considering seeker field of view angle constraints
CN113359819B (en) An Optimal Guidance Law with Collision Angle Constraint and Acceleration Limitation
CN103486904B (en)A kind of plan Velocity Pursuit method of guidance of simple and easy guided munition
CN111176315B (en)Variant cruise missile control method based on L1 adaptive control
CN111240204B (en) A cruise missile control method based on model reference sliding mode variable structure control
CN112050693B (en)Semi-strapdown guidance control integrated design method considering attack angle and view field constraint
Xu et al.Six‐Degree‐of‐Freedom Digital Simulations for Missile Guidance and Control
He et al.Coning motion stability of spinning missiles with strapdown seekers
CN118625654A (en) Integrated method of multi-constraint adaptive guidance and control for hypersonic vehicle in dive phase
CN119902538A (en) A fast fitting method for nonlinear state of aircraft
CN118533002A (en) A design method for trajectory-adjustable terminal guidance law with attack angle constraint
Qian et al.Integrated guidance and control for missiles with three-dimensional impact angle constrained
CN115615261B (en) Elasticity identification and line-of-sight angular rate extraction fusion method for rockets with large slenderness ratio
CN115685756A (en) An Optimization Method Based on Improved Gray Wolf&#39;s Glide Projectile Controller
CN115421388B (en) An Online Trajectory Planning Method Based on Convex Optimization with Multiple Attitude Constraints for the Last Stage of Long-range Missile

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant
CF01Termination of patent right due to non-payment of annual fee

Granted publication date:20210518

CF01Termination of patent right due to non-payment of annual fee

[8]ページ先頭

©2009-2025 Movatter.jp