CN113687096B - Crosswind estimation method based on embedded atmospheric data system - Google Patents
Crosswind estimation method based on embedded atmospheric data system Download PDFInfo
- Publication number
- CN113687096B CN113687096B CN202111104238.5A CN202111104238A CN113687096B CN 113687096 B CN113687096 B CN 113687096B CN 202111104238 A CN202111104238 A CN 202111104238A CN 113687096 B CN113687096 B CN 113687096B
- Authority
- CN
- China
- Prior art keywords
- beta
- bomb
- bullet
- pressure
- mother
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
- G01P13/025—Indicating direction only, e.g. by weather vane indicating air data, i.e. flight variables of an aircraft, e.g. angle of attack, side slip, shear, yaw
Abstract
The invention discloses a crosswind estimation method based on an embedded atmospheric data system. The crosswind estimation method is used for a female bullet with an inertial navigation system, wherein an embedded atmospheric data system is arranged on the female bullet and consists of a plurality of pressure measuring points which are arranged on the left side and the right side of the head of the female bullet. The crosswind estimation method includes the steps of selecting a pressure measuring point position on the surface of a mother bullet model, constructing a mathematical model of a mother bullet sideslip angle beta and pressure of the pressure measuring point, calculating the mother bullet sideslip angle beta, obtaining flight parameters of the mother bullet during the flat flight measured by an inertial navigation system, and calculating crosswind speed. The crosswind estimation method is suitable for the requirements of a new generation of aircrafts on stealth and high-precision performance, can obtain a high-precision sideslip angle, can estimate the wind speed and the wind direction of local crosswind more accurately, and has engineering popularization value.
Description
Technical Field
The invention belongs to the technical field of flight power, and particularly relates to a crosswind estimation method based on an embedded atmospheric data system.
Background
The primary and secondary bomb comprises a primary bomb and a plurality of bullets loaded in the primary bomb, and the primary and secondary bomb is mainly used for attacking the surface targets of cluster armored vehicles, technical weapon places, command communication centers, airport runways, parking ramps, highway and railway transportation hubs, large-scale electric facilities, ports, bridges, oil depots, heavy equipment factories and the like.
The actuation process of the snap bomb is as follows:
1. the shrapnel is buried in a carrier and flies to a predetermined combat airspace;
2. after the primary bomb is launched and departed from the aircraft, the combined navigation system of the primary bomb fuses the inertial navigation system and satellite positioning information, the position, speed, attitude and other information of the bomb body are measured in real time, a guidance computer of the primary bomb combines the combined target point parameters to carry out guidance law resolving, a steering engine is controlled to drive a control surface to deflect to generate aerodynamic force, and the primary bomb and the secondary bomb are enabled to accurately fly to a target according to a preset guidance rule;
3. before the mother bullet leaves the cabin and approaches a target, the mother bullet deflects through a deflection rudder to rotate (also called to be rotated) the mother bullet in order to enlarge a damage area. Before the mother bounce rotates, the mother is in a guided flight state without rolling and with an attack angle of zero basically, namely a plane flight state. In a flat flight state, measuring and pre-estimating an open cabin point region wind field through a device on a mother bomb, and correcting the open cabin point position by applying a wind correction technology (WCMD) to improve the hit precision;
4. and after the condition of opening the cabin is met, giving a cabin opening and scattering instruction by the combined navigation system of the mother bullet, scattering the bullet out of the cabin by the scattering system, and damaging the target by the action of the bullet according to a preset time sequence after the bullet is out of the cabin.
In order to ensure the accuracy of the bullet drop point, besides the integrated navigation system which needs high accuracy, the interference of crosswind during the bullet drop process must be considered. If the disturbance of the crosswind is not reasonably considered before the bullet reaches the drop point, the drop accuracy of the bullet is greatly reduced and even deviates from the predetermined target point, so that the crosswind correction technology needs to be vigorously developed, and the accurate measurement of the crosswind is the basis of the crosswind correction technology.
Currently, a traditional atmospheric data system on a parent bomb is an exposed airspeed tube and wind vane manufactured based on a pitot tube principle, and various crosswind correction technologies are developed on the basis. However, because the airspeed head and the vane occupy a large volume and have a complex structure, the requirements of the new generation of aircraft on stealth and high-precision performance cannot be met, a novel embedded air data system (FADS system) needs to be sought, and a crosswind estimation method is developed on the basis of the embedded air data system to provide basic data for a crosswind correction technology.
Disclosure of Invention
The invention aims to provide a crosswind estimation method based on an embedded atmospheric data system.
The crosswind estimation method based on the embedded atmospheric data system is used for a parent bomb with an inertial navigation system, is based on the embedded atmospheric data system arranged on the parent bomb, and consists of a plurality of pressure measuring points arranged on the left side and the right side of the head of the parent bomb; the crosswind estimation method comprises the following steps:
a. arranging pressure measuring points on the mother bullet model
a1. Determining a horizontal symmetry plane of the female warhead;
a2. determining the intersecting line of the horizontal symmetry plane and the head of the female bullet, namely a left side line L0 and a right side line R0 of the female bullet;
a3. determining a series of vertical sections perpendicular to the central axis at the head of the female bullet from front to back along the central axis of the female bullet, drawing section lines S1, S2, … … and SN, wherein the intersection points of a left side line L0 and the section lines are L01, L02, … … and L0N in sequence, the intersection points of a right side line R0 and the section lines are R01, R02, … … and R0N in sequence, and N is more than or equal to 3 and less than or equal to 6;
a4. on a section line S1, determining L11 above L01 and L21 below L01 according to experience to obtain a group of L11, L01 and L21, and obtaining L12, L02, L22, … …, L1N, L0N, L N by the same method;
a5. on a section line S1, according to experience, determining R11 above R01, determining R21 below R01 to obtain a group of R11, R01 and R21, and obtaining R12, R02, R22, … …, R1N, R0N, R N by the same method;
a6. arranging pressure measuring points on the parent bomb model according to the marking points in the steps a4 and a 5;
b. on the mother bullet model, a pressure measuring point is selected
b1. Determining flat flight state parameters of the grenade, including flight speedv e Altitude H, and speed of flightv e Step amount Δ ofv e Step amount delta alpha of an attack angle alpha and step amount delta beta of a sideslip angle beta;
b2. determining a flight speed sequence:v e -Δv e 、v e 、v e +Δv e (ii) a Determining an angle of attack sequence: -2 Δ α, - Δ α, 0 °, Δ α,2 Δ α; determining a sideslip angle sequence: -2 Δ β, - Δ β, 0 °, Δ β,2 Δ β;
b3. obtaining pressure values of each point in the steps a4 and a5 under each parameter sequence in the step b2 through wind tunnel tests or numerical simulation;
b4. calculating the pressure P of each pressure measuring point on the left side of each section line of the parent bomb under each parameter sequence in the step b2 L Pressure P of the pressure measurement points on the corresponding right side R Calculating the respective sideslip angle pressure coefficients K β by:
b5. at the same flying speedv e Fly height H and angle of attack alphaFitting a K beta-beta data curve by adopting least square normative;
b6. according to Kbeta pair of flight speedv e Selecting at least one group of left pressure measuring points and at least one group of right pressure measuring points from the points of the steps a4 and a5 as the pressure measuring points of the embedded atmospheric data system according to the principle that the change of the attack angle alpha is insensitive, namely the root mean square of the K beta is small and the linearity of a K beta-beta data curve is good;
c. mathematical model for constructing lateral slip angle beta of mother projectile and pressure of pressure measuring point
c1. Installing a pressure sensor on the parent bomb according to the pressure measuring point position of the embedded atmospheric data system determined by the parent bomb model, and meanwhile, adding a sideslip angle calculation module in the combined navigation system of the parent bomb;
c2. introducing the following formula into a sideslip angle calculation module to calculate the sideslip angle beta of the parent projectile:
in the formula, K is a proportionality coefficient, the reciprocal of K is a derivative of K beta relative to beta, after a pressure measuring point of the embedded atmospheric data system is determined, a K value and a b value of the selected pressure measuring point can be obtained by linearly fitting a K beta-beta data curve of the pressure measuring point; b is an intercept, and for a mother bullet shape which is symmetrical left and right, b =0;
d. calculating the sideslip angle beta of the parent projectile
In the flat flight stage before the opening of the cabin of the mother bullet, the sideslip angle calculation module acquires the pressure value of the pressure sensor in the delta t time period, and the pressure mean value of the pressure sensor on the left side is taken as P L Pressure mean value of the pressure sensor on the right side is P R (ii) a The sideslip angle calculation module calculates according to the formula (2) to obtain a sideslip angle beta;
e. obtaining flight parameters of the primary bomb during the flat flight measured by the inertial navigation system
Obtaining ground speed by using inertial navigation system of mother bombv e Yaw angleψDeviation angle of trajectoryψ v ；
f. Calculating crosswind velocity
Calculating crosswind by using a guidance computer of a combined navigation system of the grenade;
under the ground coordinate system, the ground coordinate system is set,x 0 taking a connection line between the projection of the parent bomb on the ground and a bomb throwing target point as a bomb target direction;x b is in the axial direction of the projectile body;v e the flying speed of the mother bullet relative to the ground, namely the ground speed,v a the moving speed of the mother bullet relative to the air is the airspeed,v w the wind speed is the plane wind;ψin order to determine the yaw angle,ψ v is the ballistic declination; the inertial navigation system of the female projectile gives the flying speed of the female projectile relative to the groundv e Yaw angleψDeviation angle of trajectoryψ v If the parent projectile flies in the direction of the bullet, the trajectory deflection angle is determinedψ v =0；
The combined navigation system of the parent bomb gives the following calculation results:
ground speedv e The component in the lateral direction isv ez =│v e │sin(ψ v )；
Airspeedv a The lateral component is |v a │sin(ψ-β)；
v a =v e - v w Due to wind speedv w Relative ground speedv e Smaller, approximately |v a │=│v e L, then airspeedv a The component in the lateral direction isv az =│v e │sin(ψ-β) ；
Wind speedv w =v e - v a Velocity of windv w The component in the lateral direction, i.e. the side wind size, is:
the crosswind estimation method based on the embedded atmospheric data system is characterized in that one or more pressure sensors are respectively arranged at the left and right positions of the surface of a grenade on the basis of grenade hardware with an inertial navigation system to form the embedded atmospheric data system, and the sideslip angle of the grenade is calculated by measuring the pressure of the embedded atmospheric data system. After the sideslip angle of the parent projectile is obtained, the wind speed and the wind direction of local crosswind are estimated by combining the movement speed, the yaw angle and the trajectory deflection angle of the parent projectile relative to the ground, which are provided by the inertial navigation system, so that basic data are provided for a crosswind correction technology.
The pressure measuring point pressure selected by the crosswind estimation method based on the embedded atmospheric data system is sensitive to change along with the sideslip angle, the linearity is good, the correlation degree along with the flight speed and the attack angle is small, the physical representation is visual, and the sideslip angle obtained through calculation is ensured to have higher precision.
The crosswind estimation method based on the embedded atmospheric data system meets the requirements of a new generation of aircrafts on stealth and high-precision performance, can obtain a high-precision sideslip angle, can estimate the wind speed and the wind direction of local crosswind more accurately, and has engineering popularization value.
Drawings
FIG. 1 is a flow chart of a crosswind estimation method based on an embedded atmospheric data system according to the present invention;
FIG. 2 is a schematic diagram (left side) of the distribution of pressure points of the crosswind estimation method based on the embedded atmospheric data system of the present invention;
FIG. 3 is a K β - β curve of example 1;
fig. 4 is a diagram of the relationship between the velocity vectors in the horizontal plane.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and examples.
The crosswind estimation method based on the embedded atmospheric data system is used for a parent bomb with an inertial navigation system, is based on the embedded atmospheric data system arranged on the parent bomb, and consists of a plurality of pressure measuring points arranged on the left side and the right side of the head of the parent bomb; the crosswind estimation method comprises the following steps:
a. arranging pressure measuring points on the mother bullet model
a1. Determining a horizontal symmetry plane of the head of the female bullet;
a2. determining the intersecting line of the horizontal symmetry plane and the head of the female bullet, namely a left side line L0 and a right side line R0 of the female bullet;
a3. determining a series of vertical sections perpendicular to the central axis at the head of the female bullet from front to back along the central axis of the female bullet, drawing section lines S1, S2, … … and SN, wherein the intersection points of a left side line L0 and the section lines are L01, L02, … … and L0N in sequence, the intersection points of a right side line R0 and the section lines are R01, R02, … … and R0N in sequence, and N is more than or equal to 3 and less than or equal to 6;
a4. on a cross-sectional line S1, determining L11 above L01 and L21 below L01 according to experience to obtain a group of L11, L01 and L21, and obtaining L12, L02, L22, … …, L1N, L0N, L N by the same method;
a5. on a section line S1, according to experience, determining R11 above R01, determining R21 below R01 to obtain a group of R11, R01 and R21, and obtaining R12, R02, R22, … …, R1N, R0N, R N by the same method;
a6. arranging pressure measuring points on the parent bomb model according to the marking points in the steps a4 and a 5;
b. on the mother bullet model, a pressure measuring point is selected
b1. Determining flat flight state parameters of the grenade, including flight speedv e Altitude H, and speed of flightv e Step amount ofv e Step amount delta alpha of an attack angle alpha and step amount delta beta of a sideslip angle beta;
b2. determining a flight speed sequence:v e -Δv e 、v e 、v e +Δv e (ii) a Determining an attack angle sequence: -2 Δ α, - Δ α, 0 °, Δ α,2 Δ α; determining a sideslip angle sequence: -2 Δ β, - Δ β, 0 °, Δ β,2 Δ β;
b3. obtaining pressure values of each point in the steps a4 and a5 under each parameter sequence in the step b2 through wind tunnel tests or numerical simulation;
b4. calculating the pressure P of each pressure measuring point on the left side of each section line of the parent bomb under each parameter sequence in the step b2 L Pressure P of the pressure measurement points on the corresponding right side R Calculating the respective sideslip angle pressure coefficients K β by:
b5. at the same flying speedv e Drawing a Kbeta-beta data curve under the flying height H and the attack angle alpha, and fitting the Kbeta-beta data curve by adopting least square normativity;
b6. according to Kbeta pair of flight speedv e Selecting at least one group of left pressure measuring points and at least one group of right pressure measuring points from the points of the steps a4 and a5 as the pressure measuring points of the embedded atmospheric data system according to the principle that the variation of the attack angle alpha is insensitive, namely the root mean square of the K beta is small and the linearity of a K beta-beta data curve is good;
c. mathematical model for constructing lateral slip angle beta of mother projectile and pressure of pressure measuring point
c1. Installing a pressure sensor on the parent bomb according to the pressure measuring point position of the embedded atmospheric data system determined by the parent bomb model, and meanwhile, adding a sideslip angle calculation module in the combined navigation system of the parent bomb;
c2. introducing the following formula into a sideslip angle calculation module to calculate the sideslip angle beta of the parent projectile:
in the formula, K is a proportionality coefficient, the reciprocal of K is a derivative of K beta relative to beta, and after a pressure measuring point of the embedded atmospheric data system is determined, the K value and the b value of the selected pressure measuring point can be obtained by linearly fitting a K beta-beta data curve of the pressure measuring point; b is an intercept, and for a mother bullet shape with left and right symmetry, b =0;
d. calculating the sideslip angle beta of the parent projectile
In the flat flight stage before the opening of the cabin of the mother bullet, the sideslip angle calculation module acquires the pressure value of the pressure sensor in the delta t time period, and the pressure mean value of the pressure sensor on the left side is taken as P L Pressure mean value of the pressure sensor on the right side is P R (ii) a The sideslip angle calculation module calculates according to the formula (2) to obtain a sideslip angle beta;
e. acquiring flight parameters of a mother bomb during flat flight measured by an inertial navigation system
Obtaining ground speed by using inertial navigation system of mother projectilev e Yaw angleψDeviation angle of trajectoryψ v ；
f. Calculating crosswind velocity
Calculating crosswind by using a guidance computer of a combined navigation system of the parent bomb;
under the condition of a ground coordinate system, the ground coordinate system,x 0 taking a connection line between the projection of the parent bomb on the ground and a bomb throwing target point as a bomb target direction;x b is in the axial direction of the elastomer;v e the flying speed of the mother bullet relative to the ground, namely the ground speed,v a the moving speed of the bullet relative to the air is the airspeed,v w the wind speed is the plane wind;ψin order to determine the yaw angle,ψ v is the ballistic declination; the inertial navigation system of the mother bomb gives the flying speed of the mother bomb relative to the groundv e Yaw angleψDeviation angle of trajectoryψ v If the parent projectile flies in the direction of the bullet, the trajectory deflection angle is determinedψ v =0；
The combined navigation system of the parent bomb gives the following calculation results:
ground speedv e The component in the lateral direction isv ez =│v e │sin(ψ v )；
Airspeedv a The component in the lateral direction is |v a │sin(ψ-β)；
v a =v e - v w Due to wind speedv w Relative ground speedv e Smaller, approximately |v a │=│v e L, then airspeedv a The component in the lateral direction isv az =│v e │sin(ψ-β) ；
Wind speedv w =v e - v a Velocity of windv w The component in the lateral direction, i.e. the side wind size, is:
example 1
The flowchart of the crosswind estimation method based on the embedded atmospheric data system of the embodiment is shown in fig. 1. The shape of the bullet is symmetrical left and right, and b =0, n =5, Δ α =2 °, Δ β =2 °. The implementation process of this embodiment is as follows:
1. selecting pressure point position on surface of mother bullet model
Drawing the pressure measuring point position shown in figure 2 on the surface of the mother bullet model; the leveling flight parameters are as follows: the flight speed V is 270m/s, 285m/s, 300 m/s, the height H =1000 meters, the attack angle alpha is-4 degrees, -2 degrees, 0 degrees, 2 degrees and 4 degrees, the sideslip angle beta is 0 degree, 2 degrees and 4 degrees, and the pressure of the pressure measuring point is obtained through a wind tunnel test; calculate the average pressure P of the left point L Average pressure P of left side point R Calculating the sideslip angle pressure coefficient K beta under each level flight parameter by the following formula:
by comparison, L12, L02 and L22 of the section S2 and R12, R02 and R22 are selected as pressure measuring points of the embedded atmospheric data system, and the pressure values are shown in a table 1;when the pressure measurement point of the S2 section has the sideslip angle beta of 0 degrees, 2 degrees and 4 degrees, the K beta values are-5 e-6, 6.75e-2 and 1.324e-1, a linear fitting K beta-beta data curve is shown in figure 3, and the linear fitting expression is K beta=0.0331 β +0.0004, correlation coefficient R 2 =0.9999, close to 1, indicating good conformity and good linearity of K β with sideslip angle variation; the root mean square values of K beta at 0 degree, 2 degree and 4 degree of sideslip are respectively 5e-6, 8e-4 and 2.3e-3, the root mean square value is far smaller than the K beta value with sideslip, the K beta value is insensitive to the change of flight and attack angle, and the K beta value accords with the standard for constructing an embedded atmospheric data system;
2. mathematical model for constructing lateral slip angle beta of mother projectile and pressure of pressure measuring point
Pressure measurement point pressures at the sideslip angle β of 0 °,2 °, and 4 ° were obtained with the flight speed V =285m/s and the angle of attack α =0, and the pressure mean value of L12, L02, and L22 was taken as P L The mean value of the pressures of R12, R02 and R22 is taken as P R (ii) a Calculating K beta, and obtaining K through least square normal fitting of the beta and the K beta; calculating the parent slip angle β by the following formula:
where the reciprocal of the proportionality coefficient K is the derivative of K β with respect to β, the derivative of K β with respect to β is 0.0331, the proportionality coefficient K =1/0.0331=30.2, b =0; the mathematical model of the sideslip angle and the gauge pressure of this embodiment is as follows:
3. calculating the sideslip angle beta of the main bullet
Pressure data of left and right pressure measuring points when the mother bomb flies flatly, P of the embodiment R =84000Pa、P L =88000 Pa, and the sideslip angle β =0.7 ° is calculated according to equation (4);
4. acquiring flight parameters of a mother bomb during flat flight measured by an inertial navigation system
The embodiment utilizes the ground speed measured by the inertial navigation systemv e =280m/s, yaw angleψ=0.1 degree ballistic declinationψ v =1.5°。
5. Calculating crosswind velocity
The crosswind velocity calculated using equation (3) is:
v wz =v ez - v az =│v e │sin(ψ v )-│v e │sin(ψ-β)=3.4m/s
the coordinate axis system is shown in figure 4,v wz the direction is from left to right (bullet back view).
Claims (1)
1. A crosswind estimation method based on an embedded atmospheric data system is characterized in that the crosswind estimation method is used for a mother bomb with an inertial navigation system, the crosswind estimation method is based on the embedded atmospheric data system installed on the mother bomb, the embedded atmospheric data system is composed of a plurality of pressure measuring points arranged on the left side and the right side of the head of the mother bomb, and the crosswind estimation method comprises the following steps:
a. arranging pressure measuring points on the mother bullet model
a1. Determining a horizontal symmetry plane of the head of the female bullet;
a2. determining the intersecting line of the horizontal symmetry plane and the head of the female bullet, namely a left side line L0 and a right side line R0 of the female bullet;
a3. determining a series of vertical sections perpendicular to the central axis at the head of the female bullet from front to back along the central axis of the female bullet, drawing section lines S1, S2, … … and SN, wherein the intersection points of a left side line L0 and the section lines are L01, L02, … … and L0N in sequence, the intersection points of a right side line R0 and the section lines are R01, R02, … … and R0N in sequence, and N is more than or equal to 3 and less than or equal to 6;
a4. on a section line S1, determining L11 above L01 and L21 below L01 according to experience to obtain a group of L11, L01 and L21, and obtaining L12, L02, L22, … …, L1N, L0N, L N by the same method;
a5. on a section line S1, according to experience, determining R11 above R01, determining R21 below R01 to obtain a group of R11, R01 and R21, and obtaining R12, R02, R22, … …, R1N, R0N, R N by the same method;
a6. arranging pressure measuring points on the parent bomb model according to the marking points in the steps a4 and a 5;
b. on the mother bullet model, a pressure measuring point is selected
b1. Determining flat flight state parameters of the grenade, including flight speedv e Altitude H, and speed of flightv e Step amount ofv e Step quantity delta alpha of an attack angle alpha and step quantity delta beta of a sideslip angle beta;
b2. determining a flight speed sequence:v e -Δv e 、v e 、v e +Δv e (ii) a Determining an attack angle sequence: -2 Δ α, - Δ α, 0 °, Δ α,2 Δ α; determining a sideslip angle sequence: -2 Δ β, - Δ β, 0 °, Δ β,2 Δ β;
b3. obtaining pressure values of each point in the steps a4 and a5 under each parameter sequence in the step b2 through wind tunnel tests or numerical simulation;
b4. calculating the pressure P of each pressure measuring point on the left side of each section line of the parent bomb under each parameter sequence in the step b2 L Pressure P of the pressure measurement points on the corresponding right side R Calculating the respective sideslip angle pressure coefficients K β by:
b5. at the same flying speedv e Drawing a K beta-beta data curve under the flying height H and the attack angle alpha, and fitting the K beta-beta data curve by adopting least square normative;
b6. according to Kbeta pair of flight speedv e Selecting at least one group of left pressure measuring points and at least one group of right pressure measuring points from the points of the steps a4 and a5 as the pressure measuring points of the embedded atmospheric data system according to the principle that the change of the attack angle alpha is insensitive, namely the root mean square of the K beta is small and the linearity of a K beta-beta data curve is good;
c. mathematical model for constructing lateral slip angle beta of female projectile and pressure of pressure measuring point
c1. Installing a pressure sensor on the parent bomb according to the pressure measuring point position of the embedded atmospheric data system determined by the parent bomb model, and meanwhile, adding a sideslip angle calculation module in the combined navigation system of the parent bomb;
c2. introducing the following formula into a sideslip angle calculation module to calculate the sideslip angle beta of the parent projectile:
in the formula, K is a proportionality coefficient, the reciprocal of K is a derivative of K beta relative to beta, after a pressure measuring point of the embedded atmospheric data system is determined, a K value and a b value of the selected pressure measuring point can be obtained by linearly fitting a K beta-beta data curve of the pressure measuring point; b is an intercept, and for a mother bullet shape with left and right symmetry, b =0;
d. calculating the sideslip angle beta of the parent projectile
In the flat flight stage before the opening of the cabin of the mother bullet, the sideslip angle calculation module acquires the pressure value of the pressure sensor in the delta t time period, and the pressure mean value of the pressure sensor on the left side is taken as P L Pressure mean value of the pressure sensor on the right side is P R (ii) a The sideslip angle calculation module calculates according to the formula (2) to obtain a sideslip angle beta;
e. acquiring flight parameters of a mother bomb during flat flight measured by an inertial navigation system
Obtaining ground speed by using inertial navigation system of mother projectilev e Yaw angleψDeviation angle of trajectoryψ v ；
f. Calculating the crosswind velocity
Calculating crosswind by using a guidance computer of a combined navigation system of the parent bomb;
under the ground coordinate system, the ground coordinate system is set,x 0 taking a connection line between the projection of the parent bomb on the ground and a bomb throwing target point as a bomb target direction;x b is in the axial direction of the elastomer;v e the flying speed of the mother bullet relative to the ground, namely the ground speed,v a is a female bulletRelative to the velocity of movement of the air i.e. the airspeed,v w the wind speed is the plane wind;ψin order to be the angle of yaw,ψ v is the ballistic declination; the inertial navigation system of the mother bomb gives the flying speed of the mother bomb relative to the groundv e Yaw angleψDeviation angle of trajectoryψ v If the parent bullet flies in the direction of the bullet, the deviation angle of the trajectoryψ v =0；
The combined navigation system of the parent bomb gives the following calculation results:
ground speedv e The component in the lateral direction isv ez =│v e │sin(ψ v )；
Airspeedv a The component in the lateral direction is |v a │sin(ψ-β)；
v a =v e - v w Due to wind speedv w Relative ground speedv e Smaller, approximately |)v a │=│v e L, then airspeedv a The component in the lateral direction isv az =│v e │sin(ψ-β) ；
Wind speedv w =v e - v a Velocity of windv w The component in the lateral direction, i.e. the side wind size, is:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111104238.5A CN113687096B (en) | 2021-09-22 | 2021-09-22 | Crosswind estimation method based on embedded atmospheric data system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111104238.5A CN113687096B (en) | 2021-09-22 | 2021-09-22 | Crosswind estimation method based on embedded atmospheric data system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113687096A CN113687096A (en) | 2021-11-23 |
| CN113687096B true CN113687096B (en) | 2023-03-21 |
Family
ID=78586720
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111104238.5A Active CN113687096B (en) | 2021-09-22 | 2021-09-22 | Crosswind estimation method based on embedded atmospheric data system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113687096B (en) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE59204569D1 (en) * | 1991-08-22 | 1996-01-18 | Porsche Ag | Transducer for recording cross wind influences on a vehicle. |
| JP2003237696A (en) * | 2002-02-12 | 2003-08-27 | Tech Res & Dev Inst Of Japan Def Agency | System for detecting angle of elevation and angle of sideslip |
| US6640173B1 (en) * | 2003-02-11 | 2003-10-28 | Visteon Global Technologiee, Inc. | System and method of controlling a vehicle having yaw stability control |
| JP2006022792A (en) * | 2004-07-09 | 2006-01-26 | Topy Ind Ltd | Wind power generating system capable of increasing energy efficiency and method of controlling wind power generating system |
| CN101832851A (en) * | 2010-05-31 | 2010-09-15 | 中国航空工业空气动力研究院 | Air flotation force measuring device of wind tunnel model |
| CN103625637A (en) * | 2013-12-04 | 2014-03-12 | 中国航空工业第六一八研究所 | Large aircraft lateral gust moderating method |
| CN105424975A (en) * | 2015-11-03 | 2016-03-23 | 中国科学院电子学研究所 | Wind sensing part applied to wind speed measuring instrument |
| CN105628051A (en) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | Embedded atmosphere measuring equipment performance evaluation method |
| CN108061813A (en) * | 2017-12-26 | 2018-05-22 | 北京航天益森风洞工程技术有限公司 | A kind of bullet train crosswind data measurin system based on five-hole probe |
| JP2018159587A (en) * | 2017-03-22 | 2018-10-11 | 株式会社豊田中央研究所 | Mobile entity motion controller, mobile entity motion control method, and mobile entity motion control program |
| CN108663536A (en) * | 2018-05-02 | 2018-10-16 | 中山大学 | A kind of air speed data modification method measuring three-dimensional velocity and ultrasonic temperature instrument |
| CN109063256A (en) * | 2017-06-30 | 2018-12-21 | 北京航空航天大学 | It is a kind of for assessing the airplane digital virtual flight simulation computing system of passenger plane airworthiness |
| CN109871628A (en) * | 2019-02-27 | 2019-06-11 | 北京航空航天大学 | It is a kind of for assessing the simulation computing system and method for amphibious aircraft seaworthiness accordance |
| CN212275200U (en) * | 2020-07-27 | 2021-01-01 | 中国空气动力研究与发展中心高速空气动力研究所 | Large-stroke single-point displacement measuring device for wind tunnel supersonic flow correction |
| CN112560191A (en) * | 2021-01-08 | 2021-03-26 | 中国空气动力研究与发展中心高速空气动力研究所 | Turboprop power airplane thrust resistance system with slipstream influence correction and performance analysis method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021524013A (en) * | 2018-05-22 | 2021-09-09 | インディアン・スペース・リサーチ・オーガニゼイションＩｎｄｉａｎ Ｓｐａｃｅ Ｒｅｓｅａｒｃｈ Ｏｒｇａｎｉｓａｔｉｏｎ | Systems and methods for detecting erroneous pressure measurements in flash air data systems using pressure patterns between adjacent ports |
-
2021
- 2021-09-22 CN CN202111104238.5A patent/CN113687096B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE59204569D1 (en) * | 1991-08-22 | 1996-01-18 | Porsche Ag | Transducer for recording cross wind influences on a vehicle. |
| JP2003237696A (en) * | 2002-02-12 | 2003-08-27 | Tech Res & Dev Inst Of Japan Def Agency | System for detecting angle of elevation and angle of sideslip |
| US6640173B1 (en) * | 2003-02-11 | 2003-10-28 | Visteon Global Technologiee, Inc. | System and method of controlling a vehicle having yaw stability control |
| JP2006022792A (en) * | 2004-07-09 | 2006-01-26 | Topy Ind Ltd | Wind power generating system capable of increasing energy efficiency and method of controlling wind power generating system |
| CN101832851A (en) * | 2010-05-31 | 2010-09-15 | 中国航空工业空气动力研究院 | Air flotation force measuring device of wind tunnel model |
| CN103625637A (en) * | 2013-12-04 | 2014-03-12 | 中国航空工业第六一八研究所 | Large aircraft lateral gust moderating method |
| CN105628051A (en) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | Embedded atmosphere measuring equipment performance evaluation method |
| CN105424975A (en) * | 2015-11-03 | 2016-03-23 | 中国科学院电子学研究所 | Wind sensing part applied to wind speed measuring instrument |
| JP2018159587A (en) * | 2017-03-22 | 2018-10-11 | 株式会社豊田中央研究所 | Mobile entity motion controller, mobile entity motion control method, and mobile entity motion control program |
| CN109063256A (en) * | 2017-06-30 | 2018-12-21 | 北京航空航天大学 | It is a kind of for assessing the airplane digital virtual flight simulation computing system of passenger plane airworthiness |
| CN108061813A (en) * | 2017-12-26 | 2018-05-22 | 北京航天益森风洞工程技术有限公司 | A kind of bullet train crosswind data measurin system based on five-hole probe |
| CN108663536A (en) * | 2018-05-02 | 2018-10-16 | 中山大学 | A kind of air speed data modification method measuring three-dimensional velocity and ultrasonic temperature instrument |
| CN109871628A (en) * | 2019-02-27 | 2019-06-11 | 北京航空航天大学 | It is a kind of for assessing the simulation computing system and method for amphibious aircraft seaworthiness accordance |
| CN212275200U (en) * | 2020-07-27 | 2021-01-01 | 中国空气动力研究与发展中心高速空气动力研究所 | Large-stroke single-point displacement measuring device for wind tunnel supersonic flow correction |
| CN112560191A (en) * | 2021-01-08 | 2021-03-26 | 中国空气动力研究与发展中心高速空气动力研究所 | Turboprop power airplane thrust resistance system with slipstream influence correction and performance analysis method |
Non-Patent Citations (4)
| Title |
|---|
| Laurence,R.J. 等.A low-cost system for wind field estimation through sensor networks and aircraft design.2015,第12页. * |
| 屈展文 ; 赵志强 ; 张惠中 ; 揭裕文 ; .民用运输类飞机侧风飞行试验研究.2019,(03),第 101-105页. * |
| 李林 ; 王立新 ; .小展弦比飞翼布局作战飞机偏航轴飞行品质评定.2009,(06),8-14页. * |
| 武涛 等.基于SINS/GPS的子母弹侧风解算方法.2011,第23卷(第1期),第31-34+39页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113687096A (en) | 2021-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109597427B (en) | Bomb random attack planning method and system based on unmanned aerial vehicle | |
| CN102607639A (en) | BP (Back Propagation) neural network-based method for measuring air data in flight state with high angle of attack | |
| US8606435B2 (en) | Method and a system for estimating a trajectory of a moving body | |
| CN106681348A (en) | Guidance and control integrated design method considering all-strapdown seeker view field constraint | |
| CN105180728B (en) | Front data based rapid air alignment method of rotary guided projectiles | |
| CN105115508B (en) | Alignment methods in rotation guided cartridge Quick air based on rear data | |
| CN108153323B (en) | A kind of high-altitude unmanned vehicle high-precision reentry guidance method | |
| CN109724624B (en) | Airborne self-adaptive transfer alignment method suitable for wing deflection deformation | |
| CN108073742B (en) | Method for estimating flight state of intercepted missile tail section based on improved particle filter algorithm | |
| CN107943079B (en) | Online estimation method for residual flight time | |
| CN105486307A (en) | Line-of-sight angular rate estimating method of maneuvering target | |
| CN113359856B (en) | Unmanned aerial vehicle designated course target point guiding method and system | |
| CN113687096B (en) | Crosswind estimation method based on embedded atmospheric data system | |
| Tian et al. | Wake encounter simulation and flight validation with UAV close formation flight | |
| CN109084772B (en) | Unscented Kalman based sight line conversion rate extraction method and system | |
| RU2498342C1 (en) | Method of intercepting aerial targets with aircraft | |
| CN111649734B (en) | Particle swarm algorithm-based strapdown seeker target positioning method | |
| CN109376364A (en) | High speed rotation bullet Aerodynamic Parameter Identification method based on Extended Kalman filter | |
| Zhang et al. | An improved MCS/INS integrated navigation algorithm for multi-rotor UAV in indoor flight | |
| CN111090830A (en) | On-orbit light pressure identification method for high-orbit non-cooperative target | |
| Zhang et al. | The ballistic missile tracking method using dynamic model | |
| CN115235475B (en) | MCC-based EKF-SLAM back-end navigation path optimization method | |
| CN114967749B (en) | Maneuvering trajectory design method for low-cost altimeter | |
| Ahmad et al. | Aerodynamic performances of cruise missile flying above local terrain | |
| CN113064443B (en) | Gain online adjustment method and damping loop control method using same |
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 |