CN111336871A - Vertical attack guidance method based on circuitous flight - Google Patents

Vertical attack guidance method based on circuitous flight Download PDF

Info

Publication number
CN111336871A
CN111336871A CN202010212599.0A CN202010212599A CN111336871A CN 111336871 A CN111336871 A CN 111336871A CN 202010212599 A CN202010212599 A CN 202010212599A CN 111336871 A CN111336871 A CN 111336871A
Authority
CN
China
Prior art keywords
target point
flight
guidance
point
arc
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
Application number
CN202010212599.0A
Other languages
Chinese (zh)
Other versions
CN111336871B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202010212599.0A priority Critical patent/CN111336871B/en
Publication of CN111336871A publication Critical patent/CN111336871A/en
Application granted granted Critical
Publication of CN111336871B publication Critical patent/CN111336871B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B15/00Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
    • F42B15/01Arrangements thereon for guidance or control

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The invention provides a vertical attack guidance method based on circuitous flight, which is beneficial to reducing the requirement of the guidance ammunition on overload capacity in the terminal guidance stage, can bypass barriers and prevent air-raid fire and plays a role in tactical deception. Setting a virtual target point of an arc raising section by applying a drop point prediction technology, and sharing maneuvering pressure of the guided munition to the front section and the middle section of a trajectory; and then, combining the arc-dropping section virtual target point, lifting the terminal point of the proportional guidance law constrained by the falling angle to be right above the actual target point, so that the guided ammunition can continuously complete the vertical attack task, and a certain terminal guidance allowance is reserved. On one hand, the maneuvering pressure of the final guidance section is reduced, the improvement of the final guidance precision and the improvement of the reliability of engineering are facilitated, on the other hand, the influence of the fall of a launching point and a target point on a guidance scheme in the actual battlefield application can be effectively reduced, and the improvement of the robustness of the scheme is facilitated.

Description

Vertical attack guidance method based on circuitous flight
Technical Field
The invention relates to the technical field of ammunition guidance control, in particular to a vertical attack guidance method based on circuitous flight, which is not limited to guidance rocket projectiles and can be used for guidance cannonballs, guidance aeronautical projectiles, missiles and the like.
Background
With the gradual development of radar detection technology and defense system in modern battlefield, the battle mission puts higher requirements on the penetration capability and survival capability of guided munitions on the battlefield. Taking axe cruise missiles as an example, the missiles fly roundly above the land, bypass barriers and prevent air fire, even further utilize the barriers to shield the other party to detect own equipment such as radars and the like, and can protrude into and attack targets from the optimal direction. Meanwhile, the adoption of circuitous flight is also beneficial to confusing the prejudgment of enemy on own tactical intention, and plays a role in tactical deception.
Meanwhile, for some warheads, the falling angle is one of the great factors determining the warhead power. Taking the example of killing and blasting the warhead, because the fragment flying direction of the warhead is mainly concentrated near the normal plane of the bullet axis, most fragments of the warhead fly to the near ground and the sky when the falling angle is small, and the exertion of the power of the warhead is greatly influenced. And taking the penetration warhead as an example, the large falling angle is beneficial to the penetration effect of the warhead. Under the condition of a small falling angle, the penetration warhead is easy to slip on the surface of a hard-spot target to cause slippage, so that the penetration into the target body is impossible, and the attack efficiency is completely lost.
A related research field of guidance law with a falling angle constraint and the like is also carried out, and a thesis of optimal polynomial guidance law with an impact angle constraint (Zhao, Yang book wave, Mo jin He, et. optimal polynomial guidance law with an impact angle constraint [ J ]. the university of Beijing Physician university, 2018, v.38; No.280 (06): 91-94 +100.) aims at the problem of hitting of an air-to-ground missile on a static target on the ground, designs a polynomial guidance law, meets the constraints of miss distance and the impact angle, and provides an analytic solution of the state of a guidance system. However, the designed optimal polynomial trajectory needs larger maneuvering capacity at the end of the trajectory, and once the maneuvering capacity is saturated due to disturbance in practical engineering application, the final miss distance is larger, and the optimal polynomial trajectory can only fly roundly in a small range at the end of the trajectory, so that the optimal polynomial trajectory cannot play a role in tactical fraud to a great extent. According to the 'limited time convergence sliding mode guidance law with falling angle constraint' (Lepengcheng, limited time convergence sliding mode guidance law with falling angle constraint [ J ]. modern defense technology, 2017 (06): 70-77 +162.), the limited time convergence sliding mode guidance law with falling angle constraint is improved and designed, and the condition that the guidance system state converges to the sliding mode surface, the missile eye line angle rate converges to zero and the trajectory angle meets the requirement of the terminal falling angle in the limited time is proved. However, in an actual guidance control system, the sliding mode variable structure controls the switching characteristic which is discontinuous in nature, and the system shakes. Buffeting can affect the stability of the control system, resulting in a reduction in the accuracy of the missile hit.
Disclosure of Invention
In view of the above, the invention provides a vertical attack guidance method based on circuitous flight, which is beneficial to reducing the requirement of the guidance ammunition on overload capacity in the terminal guidance stage, can bypass obstacles and prevent air-fire power, and plays a role in tactical deception.
The invention is realized by the following technical scheme:
in the stage from the stop of a guided ammunition engine to the top point of a trajectory, taking a virtual target point of an arc raising section as a predicted drop point, and adjusting the trajectory by adopting a zero-effect miss distance control strategy; the arc-rising section virtual target point is arranged behind the task target point and is 0.1 to 1 time of the range from the task target point; the zero-effect miss distance is the horizontal coordinate difference value of the uncontrolled flight landing point and the arc raising section virtual target point;
guiding the guided ammunition to the virtual target point of the arc-falling section by using a proportional guidance law with a falling angle constraint at the stage from the trajectory vertex to the virtual target point of the arc-falling section, wherein the constraint falling angle is selected to form a vertical trajectory when the projectile flies to the virtual target point of the arc-falling section; entering a final guidance stage when the flight height meets the set threshold height; the virtual target point of the arc-reducing section is arranged above the plumb of the task target point, and the horizontal projections of the virtual target point and the task target point are overlapped.
The threshold height is the sum of the altitude of the virtual target point of the arc-dropping section and the distance reserved for avoiding the divergence of the proportional guidance instruction.
The elevation of the virtual target point of the arc-dropping section can be selected within the range of 0.1-0.3 times of the height of the full trajectory.
In the second stage, according to ammunition speed and position information provided by a guided ammunition navigation system, calculating the horizontal distance between the current point and the predicted landing point of the uncontrolled flight; obtaining the horizontal distance between the current point and the virtual target point of the arc raising section; and obtaining the horizontal coordinate difference value of the uncontrolled flight landing point and the arc raising section virtual target point according to the difference value of the two horizontal distances.
The method comprises the following steps of obtaining an uncontrolled predicted drop point by carrying out numerical integration on a rigid six-degree-of-freedom ballistic kinetic equation of a guided ammunition through a missile-borne computer.
The method comprises the steps of taking a residual estimated trajectory numerical value integral initial value as a table head in a flight envelope range, taking an estimated falling speed output by an integral terminal point as a table value, generating a data table in advance through offline calculation, writing the data table into an missile-borne computer, and obtaining an uncontrolled predicted falling point through interpolation operation in a real-time flight process.
Wherein, the integration initial value meshing density in the flight envelope reaches simultaneously: the length grid width is less than 400m, the velocity grid width is less than 10m/s, and the angle grid width is less than 1 deg.
The method comprises the steps of adopting a proxy model to carry out regression and fitting on a data table generated offline to generate an algebraic formula, and obtaining an estimated falling speed through algebraic formula operation in the real-time flight process.
And in the final guidance stage, the attack target is guided according to the proportion to complete guidance.
Wherein the guided munition comprises a guided rocket projectile, a guided cannonball, a guided aerobomb or a guided missile.
Has the advantages that:
1. the invention sets the virtual target point of the arc raising section by applying the drop point prediction technology, shares the maneuvering pressure of the guided ammunition to the front section and the middle section of the trajectory, can bypass barriers and air defense fire, and can intrude and attack the target from the optimal direction. Meanwhile, the zero-effect predicted landing point of the trajectory of the arc raising section is far away from the real target, so that the prejudgment of the enemy on own tactical intention can be effectively confused, and the tactical deception effect is played; and then, combining the arc-dropping section virtual target point, lifting the terminal point of the proportional guidance law of the falling angle constraint right above the actual target point, so that the projection of the center of mass of the projectile body in the horizontal plane is close to the target at the end of the adjustment section, and the trajectory is nearly vertical, so that the guided ammunition can continuously complete the vertical attack task, and a certain final guidance allowance is reserved. On one hand, the maneuvering pressure of the final guidance section is reduced, the improvement of the final guidance precision and the improvement of the reliability of engineering are facilitated, on the other hand, the influence of the fall of a launching point and a target point on a guidance scheme in the actual battlefield application can be effectively reduced, and the improvement of the robustness of the scheme is facilitated.
2. In the third stage, the proportional guidance law with the falling angle constraint is adopted, so that the requirement on the maneuvering capability of the guided ammunition is greatly reduced, the requirement on the on-ammunition hardware is not high, the additional cost is low, and the implementation is easy.
Drawings
FIG. 1 is a flow chart of a vertical attack guidance method based on circuitous flight.
FIG. 2 is a diagram illustrating second-stage null-effect prediction according to the present invention.
FIG. 3 is a schematic diagram of a virtual target point of a third-stage arc-dropping segment according to the present invention.
FIG. 4 is a graph showing the results of simulation experiments of the present invention, FIG. 4(a) is a diagram of ballistic curves ① - ③ in simulation situations, and FIG. 4(b) is a diagram of ballistic inclination angles ① - ③ in simulation situations.
Detailed Description
The invention is described in detail below by way of example with reference to the accompanying drawings.
The invention provides a vertical attack guidance method based on circuitous flight, which is divided into four stages according to the flight trajectory of guided munitions for convenient expression, wherein the first stage is from the launching time of the guided munition to t1Time of day t1=tmax+ δ t, where tmaxIs the maximum operating time of the ammunition engine under different temperature conditions, deltat is the design margin taken to ensure that the first stage covers the operating time of the engine,generally, the time is more than 1 s; the second stage is from t1Moment to ballistic vertex; the third stage is that the altitude from the trajectory vertex to the center of mass of the guided munition is reduced to the altitude close to the virtual target point of the arc-dropping section, and the altitude is mathematically described as H being less than or equal to Hd+ δ H. Wherein Hd+ δ H is called the threshold height, H is the current flight altitude of the guided munition, HdThe elevation of the virtual target point of the arc-dropping section designed according to the tactical use requirement is selected within the range of about 0.1-0.3 times of the total trajectory height, and the delta H leaves a distance for avoiding the divergence of the proportional guidance instruction and can be 500m +/-200 m; and the fourth stage is a final guidance stage, and the final guidance is carried out by adopting classical proportion guidance or a proportion guidance law with falling angle constraint, so that the guidance error is further reduced under the condition of ensuring vertical attack.
The flow chart of the invention is shown in fig. 1, and specifically comprises the following steps:
the first stage, the flight track control is carried out;
in the second stage, an arc-rising section virtual target point is arranged at a distance of 0.1 to 1 times of the range from the target point behind the task target point; obtaining the coordinates of the predicted uncontrolled flight landing point through a rigid six-degree-of-freedom trajectory kinetic equation of the guided munition; obtaining a deflection instruction of a steering engine of the guided ammunition according to a horizontal coordinate difference value between the uncontrolled flight falling point and the virtual target point of the arc raising section, deflecting the steering engine under the action of the instruction to enable the guided ammunition to maneuver by taking the target point of the arc raising section as a predicted falling point, and entering a third stage when the trajectory inclination angle is less than 0 degrees;
in order to realize circuitous flight, an arc raising section virtual target point is set according to a tactical mission in the stage; and adjusting the trajectory according to the relation between the virtual target point of the arc raising section and the predicted uncontrolled flight landing point (the difference value of the two is the zero-effect miss distance).
The arc-lifting section virtual target point can be arranged behind the task target point by 0.1 to 1 time of range from the target point according to requirements. If the virtual target point of the arc-lifting section is in front of or behind the task target point but the distance is less than 0.1 time of the range, the circuitous flying effect is not obvious by adopting the invention. Meanwhile, the roundabout flight turning distance is larger as the distance of the arc raising section virtual target point behind the task target point is larger. However, when the virtual target point of the arc-lifting section is behind the task target point and the distance is greater than 1 time of the range, on one hand, the effect of further increasing the turn-back distance of roundabout flight is not obvious, on the other hand, a more rigorous requirement is provided for the maneuvering capability of the guided rocket, so that the waste of missile-borne resources can be caused, and the cost is greatly increased.
For a rocket projectile, the coordinates of the predicted uncontrolled flight landing point P can be obtained by numerical integration of the missile-borne computer, and the kinetic equation is as follows:
Figure BDA0002423331020000051
the formula (2) is a rigid six-degree-of-freedom ballistic kinetic equation of a conventional rocket projectile, and can accurately describe the mass center motion and the attitude motion of the rocket projectile. For the rocket projectiles at medium and long distances, inertial navigation is mounted on the common projectiles, the initial value of the state variable in the formula (2) can be obtained through the current value output by the inertial navigation, the kinetic equation is subjected to numerical integration, and the integral termination condition is that y is less than or equal to ytAnd obtaining the uncontrolled predicted drop point P of the rocket projectile.
Numerical integration is carried out on the formula (2), and the obtained predicted drop point precision can meet the guidance requirement, but the time-consuming problem caused by the numerical integration cannot be completely avoided. Therefore, a residual estimated trajectory numerical value integral initial value is used as a table head in a flight envelope range, an estimated falling speed output by an integral terminal point is used as a table value, a data table is generated in advance through offline calculation and written into an missile-borne computer, and a predicted falling point is obtained through interpolation operation in the real-time flight process. Trial calculation shows that when the grid division density in the flight envelope simultaneously reaches the length grid width below 400m, the speed grid width below 10m/s and the angle grid width below 1 degree, the requirement of the guidance method can be met.
Further, in the case where there is mesh correspondence data for estimating the amount of miss, the speed of obtaining the predicted amount of miss by interpolation is related to the size of the data amount of the stored mesh. If the flight envelope range is large or the grid division is dense, the speed of interpolation calculation is very significantly slow, even slower than the speed of direct online prediction by adopting rigid six-degree-of-freedom model integration. Therefore, a polynomial function or a proxy model such as a neural network can be adopted to carry out regression and fitting on the off-line generated prior data to generate an algebraic formula, and algebraic calculation is adopted to replace interpolation calculation, so that the calculation amount is further reduced.
As shown in FIG. 2, the task target point T has the coordinate of (x)T,yT) The set virtual target point T' of the arc rising section has the coordinates of
Figure BDA0002423331020000061
The predicted uncontrolled flight landing point P coordinate is (x)P,yP) Wherein y isP=yTZero effective miss distance
Figure BDA0002423331020000062
The control amount output at this time is the angle of attack αzExpressed as:
Figure BDA0002423331020000063
wherein K1Is a proportionality coefficient, is selected by numerical simulation according to specific guided ammunition products, the y coordinates of a virtual target point T' of an arc rising section and a task target point T are equal, the x coordinate is different by a target offset L,
Figure BDA0002423331020000064
is an important variable for realizing circuitous flight.
The third stage, firstly, selecting a virtual target point of the arc descending section above the plumb of the task target point, superposing the horizontal projections of the virtual target point and the task target point, then guiding the guided ammunition to the virtual target point of the arc descending section by using a proportional guidance law with a falling angle constraint, wherein the constraint falling angle is selected to enable the missile to be close to forming a vertical trajectory when flying to the vicinity of the virtual target point of the arc descending section, and when the missile approaches the virtual target point and the flying height of the missile meets the set threshold height, switching to the fourth stage;
in order to lift the terminal point of the falling angle constraint proportional guidance law to be right above the actual target point, the projection of the center of mass of the projectile body in the horizontal plane is close to the target at the final stage of the adjustment section, and meanwhile, the terminal guidance allowance is reserved, after the flying posture of the missile is judged to enter the third stage, the virtual target point T' of the arc-reducing section is selected to be vertically above the task target point T, and the horizontal projections of the two are overlapped, as shown in fig. 3. And guiding the guided munition to the virtual target point of the arc-descending segment by using a proportional guidance law with a falling angle constraint, wherein the constraint falling angle is selected to be close to forming a vertical trajectory when the projectile body flies to the vicinity of the virtual target point of the arc-descending segment, and the projection of the center of mass of the guided munition on the horizontal plane is close to and coincident with the mission target point.
And fourthly, guiding the attack target according to the proportion to finish guidance.
The terminal guidance can be carried out by adopting classical proportion guidance or proportion guidance law with falling angle constraint, the guidance error is further reduced under the condition of ensuring vertical attack, and the terminal guidance is not detailed because the prior art is adopted.
Furthermore, in the second stage, according to ammunition speed and position information provided by a guided ammunition navigation system, the horizontal distance between a task target point and an uncontrolled flight predicted landing point is calculated by means of numerical integration and the like, the horizontal distance between the current task target point and a virtual target point is solved by a space distance solving method, a deflection instruction of a guided ammunition steering engine is obtained by multiplying a certain scale factor according to the difference value between the task target point and the virtual target point, and the steering engine deflects under the action of the instruction, so that the guided ammunition is maneuvered by taking the virtual target point as the predicted landing point.
The navigation system related in the invention is a navigation system which adopts means such as inertial navigation, satellite navigation, ground radio positioning navigation and the like, or adopts other technical measures to obtain motion state information such as speed, position and the like in flight.
The invention relates to a vertical attack guidance method of a guided rocket projectile based on miss distance estimation, which comprises the following specific implementation steps:
step 1, adopting flight path control guidance in the first stageLaw, trajectory inclination of control rocket projectile trajectory inclination tracking theory according to formula αtrack=Kt×(θ-θtrack) And calculating to obtain the control quantity of the attack angle. Wherein, KtAnd the scale factor is selected through numerical simulation according to specific rocket projectile products. Theta, thetatrackThe real-time trajectory inclination angle and the theoretically set trajectory inclination angle of the rocket projectile are respectively used.
Step 2, establishing a ground coordinate system for describing missile motion in the second stage, wherein an X axis is positioned in a horizontal plane and points to the direction of a gun mesh connecting line, a Y axis is vertical to the horizontal plane and upwards, an X, Y, Z axis forms a launching coordinate system, and a motion model for establishing the second stage of a kinetic equation is given by an equation (3), wherein v, theta, psi in the equation (3)VRespectively the speed, trajectory inclination angle and trajectory deflection angle of the guided rocket projectile, m is mass, g is gravity acceleration,
Figure BDA0002423331020000081
psi, gamma are pitch angle, yaw angle, roll angle α*、β*For the angle of attack and the angle of sideslip,
Figure BDA0002423331020000082
is the velocity ramp angle. x, y and z are X, Y, Z axis coordinates of the position of the guided rocket projectile respectively, and X, Y, Z are forces in three directions obtained by decomposing aerodynamic force applied to the guided rocket projectile according to a speed coordinate system and are called as resistance, lift force and lateral force respectively. Its expression can be written as:
Figure BDA0002423331020000083
wherein,
Figure BDA0002423331020000084
is the zero lift drag coefficient of the projectile,
Figure BDA0002423331020000085
the coefficients of lift and lateral forces for the projectile are dimensionless scaling coefficients, in particular, for the sake of simplicity of calculationThe consideration of the process is that,
Figure BDA0002423331020000086
may all take zero.
The aerodynamic moment and the aerodynamic force are similar in calculation method, and the expression is as follows:
Figure BDA0002423331020000087
wherein M isx、My、MzThe rolling force, the yawing moment and the pitching moment which are respectively applied to the projectile body; m isx、my、mzThe rolling moment coefficient, the yawing moment coefficient and the pitching moment coefficient of the projectile body are obtained; q is dynamic pressure, SrefIs the maximum cross-sectional area. The formula for q is:
Figure BDA0002423331020000088
wherein rho is the atmospheric density, the atmospheric density rho is the function of the flight height y, the example adopts the international standard atmospheric condition, and the atmospheric density is fitted into a 6 th-order polynomial function of the flight height y, and the fitting expression is as follows:
ρ=λ1y62y53y44y35y26y+λ7(7)
values of fitting coefficients of various items are shown in table 1:
TABLE 1 atmospheric Density fitting coefficients
Figure BDA0002423331020000089
Figure BDA0002423331020000091
C in formula (4)d0The coefficient of zero lift resistance is related to the Mach number Ma of the guided rocket and is obtained through a wind tunnel force measurement testGiven in the form of a table function, there are different results for the types of projectiles of a missile, for example table 2 is a table function of zero-liter resistance coefficients of a certain projectile, which is solved in the calculation by linear interpolation.
TABLE 2 coefficient of zero lift drag Cd0Tabular function
Ma 0.4 0.9 1.15 1.5 2 3 4
Cd0 0.4095 0.4843 0.6720 0.5846 0.4894 0.3707 0.3057
Wherein the Mach number Ma of the missile is the speed and the local sound velocity vsThe ratio of (A) to (B):
Figure BDA0002423331020000092
in the invention, the sound velocity is a function of the flying height y, and the sound velocity is fitted into a 6 th-order polynomial function of the flying height y according to the international standard atmospheric condition, wherein the fitting expression is as follows:
vs=η1y62y53y44y35y26y+η7(9)
values of fitting coefficients of the terms are shown in table 3:
TABLE 3 Sound velocity fitting coefficient Table
Figure BDA0002423331020000093
Figure BDA0002423331020000101
Step 3, current time t0Current position v of guided rocket received from navigation system0、θ0
Figure BDA0002423331020000104
x0、y0、z0
Figure BDA0002423331020000102
ψ0、γ0And (3) performing numerical integration solution on the formula (2) by adopting a fixed step length numerical integration algorithm, and taking an initial value of the integral as the current time t of the missile motion information measured by the navigation system0Of motion information values, i.e.
Figure BDA0002423331020000103
Taking integral termination condition as y<yt
Step 4, calculating the virtual target point (x) from the launching point to the arc ascending sectiont+L,yt) Is horizontally thrownShadow distance R ═ xt+ L. Wherein (x)t,yt) The horizontal projection coordinates of the task target points are pre-stored in an missile-borne computer of the guided rocket before launching; l is the target point offset;
using the formula Δ R-R' xT-xP+ L obtaining the zero-effect miss distance delta R to the virtual target point of the arc-raising section, and then α according to the control formula of the arc-raising sectionz=K1× Delta R calculates an attack angle instruction to obtain a deflection instruction of a steering engine of the guided rocket projectile, and the steering engine deflects under the action of the instruction to make the projectile maneuver by taking a virtual target point of an ascending arc section as a predicted drop pointPIs the X-axis coordinate of the integral endpoint of equation (3).
And 5, entering a third stage of the guidance scheme when the trajectory inclination angle is less than 0 DEG, and setting a virtual target point of the arc-dropping section. And selecting the virtual target point T 'of the arc-reducing section to be vertically above the actual target point T, and superposing the horizontal projections of the virtual target point T' and the actual target point T. Guiding the rocket projectile to a virtual target point of a descending arc section by using a proportional guidance law with a falling angle constraint, and expressing an instruction formula by using an attack angle as follows:
Figure BDA0002423331020000111
wherein, N, K2The specific values can be calculated by solving the optimal control problem off line according to specific rocket projectile products, and can also refer to the value range N ∈ [2,20 ]],K2∈[-50,-2]And (4) selecting according to trajectory simulation data. V is the relative flight speed of the rocket projectile and the virtual target of the descending arc section,
Figure BDA0002423331020000112
angular velocity of the projectile in line of sight, theta, for the rocket to reach the virtual targetdThe trajectory inclination angle of the rocket projectile at the current moment and the trajectory inclination angle required by the battle mission are theta in the vertical attack guidance missiondTake 90 degrees. t is tgoThe remaining flight time of the projectile to the virtual target may be estimated by dividing the distance of the projectile from the virtual target by the relative flight speed of the projectile.
When the missile approaches the virtual target point of the arc-dropping section, the trajectory inclination angle of the missile is very close to the trajectory inclination angle required by the battle mission according to the characteristic of the falling angle constraint proportion guidance law. When the flying altitude is lower than the set threshold height HdAfter + δ H, the fourth stage is entered. Wherein HdThe elevation of the virtual target in the arc-dropping section designed according to the tactical use requirement is selected within the range of about 0.1-0.3 times of the height of the whole trajectory, and the delta H leaves a distance for avoiding the divergence of the proportional guidance instruction and can be 500m +/-200 m.
And 6, after the fourth stage is carried out, the guided rocket projectile flies to the target according to the control instruction calculated by the proportional guidance law, and because the guided rocket projectile forms a vertical trajectory inclination angle and the projection of the center of mass of the guided rocket projectile in the horizontal plane of the target point almost coincides with the task target, a vertical falling angle is naturally formed when the guided rocket projectile hits the target.
The method is characterized in that different radio angles are uniformly adopted for emission, different offsets L, target positions (20, 0, 0) (unit: km) and threshold heights (5 km) are set, and numerical simulation is carried out on different target point positions shown in a table 4 by adopting a guidance method introduced by the patent. The simulation results are shown in fig. 4.
TABLE 4 simulation scenarios and Main index parameters
Simulation scenario Angle of emanation/° L/m Miss amount/m Falling angle/°
30 5000 0.1289 -89.4625
25 2000 0.2624 -90.5693
40 2000 0.0396 -90.5165
As can be seen from the graph, the guidance method of the invention realizes the circuitous flight of the guided rocket, ensures the vertical attack on the task target and the miss distance within 0.3 m, and has higher precision. The simulation conditions and the results can be compared to obtain the theta required by the first stage of the schemetrackThe size of the offset L determines the circuitous flight degree of the rocket projectile, in the simulation scene ②, if the value of L is too small, the circuitous flight is not obvious, and in the simulation scene ①, the value of L is relatively large, and the obvious circuitous flight occurs.
In summary, the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. A vertical attack guidance method based on circuitous flight is characterized in that in the stage from the stop of a guided ammunition engine to the top of a trajectory, a zero-effect miss control strategy is adopted to adjust the trajectory by taking a virtual target point of an arc raising section as a predicted drop point; the arc-rising section virtual target point is arranged behind the task target point and is 0.1 to 1 time of the range from the task target point; the zero-effect miss distance is the horizontal coordinate difference value of the uncontrolled flight landing point and the arc raising section virtual target point;
guiding the guided ammunition to the virtual target point of the arc-falling section by using a proportional guidance law with a falling angle constraint at the stage from the trajectory vertex to the virtual target point of the arc-falling section, wherein the constraint falling angle is selected to form a vertical trajectory when the projectile flies to the virtual target point of the arc-falling section; entering a final guidance stage when the flight height meets the set threshold height; the virtual target point of the arc-reducing section is arranged above the plumb of the task target point, and the horizontal projections of the virtual target point and the task target point are overlapped.
2. The vertical attack guidance method based on detour flight according to claim 1, wherein the threshold height is the sum of the altitude of the virtual target point of the arc-dropping segment and the distance set aside for avoiding the divergence of the proportional guidance command.
3. The vertical attack guidance method based on roundabout flight according to claim 1, wherein the elevation of the arc-dropping virtual target point is selectable within the range of 0.1-0.3 times of the full ballistic height.
4. The vertical attack guidance method based on roundabout flight according to claim 1, characterized in that in the second phase, the horizontal distance between the current point and the predicted landing point of uncontrolled flight is calculated according to ammunition speed and position information provided by a guided ammunition navigation system; obtaining the horizontal distance between the current point and the virtual target point of the arc raising section; and obtaining the horizontal coordinate difference value of the uncontrolled flight landing point and the arc raising section virtual target point according to the difference value of the two horizontal distances.
5. The vertical attack guidance method based on circuitous flight according to claim 1, characterized in that the uncontrolled predicted drop point is obtained by numerical integration of rigid six-degree-of-freedom ballistic kinetic equations of guided munitions by an onboard computer.
6. The roundabout-type-flight-based vertical attack guidance method according to claim 5, wherein a residual estimated ballistic numerical value integral initial value is used as a header in a flight envelope range, an estimated falling speed output by an integral end point is used as a table value, a data table is generated in advance through offline calculation and written into a missile-borne computer, and an uncontrolled predicted falling point is obtained through interpolation operation in a real-time flight process.
7. The vertical attack guidance method based on roundabout flight according to claim 6, wherein the integration initial value meshing density in the flight envelope simultaneously achieves: the length grid width is less than 400m, the velocity grid width is less than 10m/s, and the angle grid width is less than 1 deg.
8. The vertical attack guidance method based on detour flight as claimed in claim 6, wherein the proxy model is adopted to carry out regression and fitting on the data table generated off-line to generate an algebraic formula, and the estimated falling speed is obtained through algebraic formula operation in the real-time flight process.
9. The vertical attack guidance method based on roundabout flight according to claim 1, characterized in that the terminal guidance phase guides the attack targets according to the proportion to complete guidance.
10. The vertical attack guidance method based on roundabout flight according to claim 1, wherein the guided munition comprises a guided rocket projectile, a guided cannonball, a guided missile or a guided missile.
CN202010212599.0A 2020-03-24 2020-03-24 Vertical attack guidance method based on circuitous flight Active CN111336871B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010212599.0A CN111336871B (en) 2020-03-24 2020-03-24 Vertical attack guidance method based on circuitous flight

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010212599.0A CN111336871B (en) 2020-03-24 2020-03-24 Vertical attack guidance method based on circuitous flight

Publications (2)

Publication Number Publication Date
CN111336871A true CN111336871A (en) 2020-06-26
CN111336871B CN111336871B (en) 2021-04-02

Family

ID=71180393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010212599.0A Active CN111336871B (en) 2020-03-24 2020-03-24 Vertical attack guidance method based on circuitous flight

Country Status (1)

Country Link
CN (1) CN111336871B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856924A (en) * 2020-08-06 2020-10-30 西安睿高测控技术有限公司 Control method of guided ammunition using relay type steering engine
CN112034702A (en) * 2020-08-06 2020-12-04 西安睿高测控技术有限公司 Intelligent control method for guided ammunition by using relay type steering engine
CN114754628A (en) * 2022-03-31 2022-07-15 南京理工大学 Flight body trajectory control method based on drop point prediction and virtual tracking
CN116362163A (en) * 2023-06-01 2023-06-30 西安现代控制技术研究所 Nonsingular multi-constraint trajectory rapid optimization method
CN117113619A (en) * 2023-06-09 2023-11-24 中国人民解放军92493部队试验训练总体研究所 Simulation method for coastal target area environment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007089243A2 (en) * 2005-02-07 2007-08-09 Bae Systems Information And Electronic Systems Integration Inc. Optically guided munition control system and method
CN102353301A (en) * 2011-09-15 2012-02-15 北京理工大学 Guidance method with terminal restraint based on virtual target point
CN104792232A (en) * 2015-04-28 2015-07-22 北京理工大学 Minimum overload terminal guiding method with terminal angular constraint
CN110319736A (en) * 2019-06-10 2019-10-11 西北工业大学 A kind of STT missile method based on vertical strike guidance law over the ground

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007089243A2 (en) * 2005-02-07 2007-08-09 Bae Systems Information And Electronic Systems Integration Inc. Optically guided munition control system and method
CN102353301A (en) * 2011-09-15 2012-02-15 北京理工大学 Guidance method with terminal restraint based on virtual target point
CN104792232A (en) * 2015-04-28 2015-07-22 北京理工大学 Minimum overload terminal guiding method with terminal angular constraint
CN110319736A (en) * 2019-06-10 2019-10-11 西北工业大学 A kind of STT missile method based on vertical strike guidance law over the ground

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王蒙一,江涌: "一种基于零效脱靶量的末制导律", 《现代防御技术》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856924A (en) * 2020-08-06 2020-10-30 西安睿高测控技术有限公司 Control method of guided ammunition using relay type steering engine
CN112034702A (en) * 2020-08-06 2020-12-04 西安睿高测控技术有限公司 Intelligent control method for guided ammunition by using relay type steering engine
CN114754628A (en) * 2022-03-31 2022-07-15 南京理工大学 Flight body trajectory control method based on drop point prediction and virtual tracking
CN114754628B (en) * 2022-03-31 2023-08-04 南京理工大学 Flying body trajectory control method based on drop point prediction and virtual tracking
CN116362163A (en) * 2023-06-01 2023-06-30 西安现代控制技术研究所 Nonsingular multi-constraint trajectory rapid optimization method
CN116362163B (en) * 2023-06-01 2023-09-01 西安现代控制技术研究所 Multi-constraint trajectory rapid optimization method
CN117113619A (en) * 2023-06-09 2023-11-24 中国人民解放军92493部队试验训练总体研究所 Simulation method for coastal target area environment
CN117113619B (en) * 2023-06-09 2024-02-06 中国人民解放军92493部队试验训练总体研究所 Simulation method for coastal target area environment

Also Published As

Publication number Publication date
CN111336871B (en) 2021-04-02

Similar Documents

Publication Publication Date Title
CN111336871B (en) Vertical attack guidance method based on circuitous flight
CN111306998B (en) Parameter perturbation self-adaptive guided rocket projectile vertical attack guidance method
CN111473696B (en) Guidance rocket vertical attack guidance method based on drop point estimation
CN103090728B (en) Tail angle restraining guidance method based on sliding mode control
CN111580547B (en) Hypersonic aircraft formation control method
CN110032206B (en) Long-distance guidance aircraft large-falling-angle top attack control method and control system
CN111649624B (en) Space miniature precise guided weapon control method
CN111692919B (en) Precise guidance control method for aircraft with ultra-close range
CN114035616B (en) Method and system for controlling striking of aircraft to moving target
CN112526872B (en) Method for processing guidance and terminal guidance handover and guidance information in constraint with large falling angle
CN113739635B (en) Guidance method for realizing missile large-sector-angle launching
CN116360489A (en) Collaborative middle guidance law design method based on forming guidance law
CN114610057B (en) Design method for maneuver burst prevention strategy of high Mach aircraft
Kim et al. Guidance and control system design for impact angle control of guided bombs
CN114610077A (en) Multi-hypersonic aircraft trajectory planning method and system
Baba et al. Air combat guidance law for an UCAV
CN111486755B (en) Falling speed control method of variable pneumatic profile guidance equipment
CN117688727B (en) Short-distance air combat maneuver planning method considering over-stall maneuver of airplane
CN116579118B (en) Online quick calculation method for three-dimensional dynamic emission area for combat in short-distance air combat
CN114608391B (en) Cannonball guidance method and system with stealth effect
CN118089485B (en) Method for recovering projectile body after thermal ablation
Morrison A lesson learned about cannon-launched guided projectiles
CN118034337A (en) Aircraft guidance control method for reverse incline target
CN118550306A (en) Aircraft reentry track planning control method with altitude control
CN115826608A (en) Multi-missile online collaborative track planning method based on waypoints

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