CN106228014A - A kind of acquisition methods of missile aerodynamic coefficient - Google Patents
A kind of acquisition methods of missile aerodynamic coefficient Download PDFInfo
- Publication number
- CN106228014A CN106228014A CN201610600360.4A CN201610600360A CN106228014A CN 106228014 A CN106228014 A CN 106228014A CN 201610600360 A CN201610600360 A CN 201610600360A CN 106228014 A CN106228014 A CN 106228014A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- missile
- angle
- aerodynamic
- coordinate system
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
The present invention relates to the acquisition methods of a kind of missile aerodynamic coefficient, the method solves the problem that polar coordinate system Aerodynamic Coefficient can not use under rectangular coordinate system.The method obtains missile trajectory parameter after first passing through the step that wind tunnel test obtains Aerodynamic Coefficient under the polar coordinate system of guided missile, and obtains missile aerodynamic parameter according to missile trajectory parameter calculation;Wherein, described trajectory parameter is flight speedWith flight attitude angle, described missile aerodynamic parameter is flying drilling angle α and flight sideslip angle beta;The total angle of attack under polar coordinate system is obtained further according to flying drilling angle α, flight sideslip angle beta and conversion formula∑With roll angle Φ;Afterwards according to described total angle of attack∑With roll angle Φ, interpolation obtains interpolation polar coordinate Aerodynamic Coefficient;Rectangular coordinate system Aerodynamic Coefficient is obtained further according to described interpolation polar coordinate Aerodynamic Coefficient and conversion formula;And then obtain acceleration and the angular acceleration of guided missile;Finally obtain the trajectory parameter flight speed of guided missile subsequent time
Description
Technical field
The invention belongs to Aero-Space wind tunnel test data processing field, be specifically related to the acquisition of a kind of missile aerodynamic coefficient
Method.
Background technology
According to body and airflow space geometric angle relation, the Aerodynamic Coefficient that Missile Design uses is divided into two kinds of systems: straight
Angular coordinate system and polar coordinate system.Rectangular coordinate system Aerodynamic Coefficient changes with angle of attack and sideslip angle beta bidimensional, and data point exists
Spatially present matrix distribution.Polar coordinate system Aerodynamic Coefficient is with total angle of attackΣChanging with roll angle Φ bidimensional, data point is at sky
Presenting circular distribution between, closeer near the data of circle centre position, the data away from circle centre position are sparse.Rectangular coordinate system is pneumatic
Coefficient and polar coordinate system Aerodynamic Coefficient, in addition to a small amount of data point overlaps, most of data point is misaligned, point
Cloth rule differs greatly.
There are differences between rectangular coordinate system Aerodynamic Coefficient and polar coordinate system Aerodynamic Coefficient, cause Ballistic Simulation of Underwater and control
Larger difference is there is in loop processed designs etc. in the Aerodynamic Coefficient mode using two kinds of coordinate-systems.Therefore, Ballistic Simulation of Underwater and control
Loop processed design can only use the Aerodynamic Coefficient under same coordinate system, and the Aerodynamic Coefficient of two kinds of coordinate-systems can not be general,
Change Ballistic Simulation of Underwater and control the coordinate-system of loop design, it is necessary to obtaining identical seat by modes such as wind tunnel test, CFD calculating
Aerodynamic Coefficient under mark system.
Summary of the invention
The purpose of the present invention: the main object of the present invention is, finds out the acquisition methods of a kind of missile aerodynamic coefficient in other words
It is the test method obtaining missile aerodynamic coefficient, solves what polar coordinate system Aerodynamic Coefficient can not use under rectangular coordinate system
Problem.
Technical scheme:
Technical solution of the present invention:
The acquisition methods of a kind of missile aerodynamic coefficient is provided,
The first step, obtains the Aerodynamic Coefficient under the polar coordinate system of guided missile by wind tunnel test;
Second step, it is thus achieved that missile trajectory parameter, and obtain missile aerodynamic parameter according to missile trajectory parameter calculation;Wherein,
Described trajectory parameter is flight speedWith flight attitude angle, described missile aerodynamic parameter is flying drilling angle α and flight side
Sliding angle beta;
3rd step, obtains the total angle of attack under polar coordinate system according to flying drilling angle α, flight sideslip angle beta by equation below
α∑With roll angle Φ, described formula is:
CosA=cos β .cos α
4th step, according to described total angle of attack∑With roll angle Φ, interpolation obtains interpolation polar coordinate Aerodynamic Coefficient, described
Interpolation polar coordinate Aerodynamic Coefficient includes axial force coefficientNormal force coefficient Cξ(αΣ, Φ), cross force Cη(αΣ,Φ)、
Rolling moment coefficientYawing moment coefficient mξ(αΣ, Φ), pitching moment coefficient mη(αΣ,Φ);
5th step, obtains the pneumatic system of rectangular coordinate system according to described interpolation polar coordinate Aerodynamic Coefficient and conversion formula
Number, described rectangular coordinate system Aerodynamic Coefficient is axial force coefficientNormal force coefficientCross force
CoefficientRolling moment coefficientYawing moment coefficientAnd pitching moment coefficientConversion formula therein is as follows:
With
6th step, the rectangular coordinate system Aerodynamic Coefficient using the 5th step to obtain carries out body kinetics resolving, is led
The acceleration of bullet and angular acceleration;
7th step, the acceleration of the guided missile that use the 6th step obtains and angular acceleration, carry out kinesiology resolving, obtain guided missile
The trajectory parameter flight speed of subsequent timeWith flight attitude angle;
8th step, resolves the aerodynamic parameter obtaining guided missile subsequent time according to the trajectory parameter of subsequent time.
Advantages of the present invention:
The present invention is converted into the Aerodynamic Coefficient under polar coordinate system under right angle mark system permissible by conversion using method
The Aerodynamic Coefficient used, saves development cost and the time cost of the consumings such as CFD calculating, wind tunnel test.
Detailed description of the invention:
Below the present invention is described in further detail.
The acquisition methods of a kind of missile aerodynamic coefficient is provided,
The first step, obtains the Aerodynamic Coefficient under the polar coordinate system of guided missile by wind tunnel test;
Second step, it is thus achieved that missile trajectory parameter, and obtain missile aerodynamic parameter according to missile trajectory parameter calculation;Wherein,
Described trajectory parameter is flight speedWith flight attitude angle, described missile aerodynamic parameter is flying drilling angle α and flight side
Sliding angle beta;
3rd step, obtains the total angle of attack under polar coordinate system according to flying drilling angle α, flight sideslip angle beta by equation below
αΣWith roll angle Φ, described formula is:
CosA=cos β .cos α
4th step, according to described total angle of attackΣWith roll angle Φ, interpolation obtains interpolation polar coordinate Aerodynamic Coefficient, described
Interpolation polar coordinate Aerodynamic Coefficient includes axial force coefficientNormal force coefficient Cξ(αΣ, Φ), cross force Cη(αΣ,Φ)、
Rolling moment coefficientYawing moment coefficient mξ(αΣ, Φ), pitching moment coefficient mη(αΣ,Φ);
5th step, obtains the pneumatic system of rectangular coordinate system according to described interpolation polar coordinate Aerodynamic Coefficient and conversion formula
Number, described rectangular coordinate system Aerodynamic Coefficient is axial force coefficientNormal force coefficientCross force
CoefficientRolling moment coefficientYawing moment coefficientAnd pitching moment coefficientConversion formula therein is as follows:
With
6th step, the rectangular coordinate system Aerodynamic Coefficient using the 5th step to obtain carries out body kinetics resolving, is led
The acceleration of bullet and angular acceleration;
7th step, the acceleration of the guided missile that use the 6th step obtains and angular acceleration, carry out kinesiology resolving, obtain guided missile
The trajectory parameter flight speed of subsequent timeWith flight attitude angle etc.;
8th step, resolves the aerodynamic parameter obtaining guided missile subsequent time according to the trajectory parameter of subsequent time.
Claims (1)
1. an acquisition methods for missile aerodynamic coefficient, the method comprises the steps:
The first step, obtains the Aerodynamic Coefficient under the polar coordinate system of guided missile by wind tunnel test;
Second step, it is thus achieved that missile trajectory parameter, and obtain missile aerodynamic parameter according to missile trajectory parameter calculation;Wherein, described
Trajectory parameter be flight speedWith flight attitude angle, described missile aerodynamic parameter is flying drilling angle α and flight yaw angle
β;
3rd step, obtains the total angle of attack under polar coordinate system according to flying drilling angle α, flight sideslip angle beta by equation below∑With
Roll angle Φ, described formula is:
CosA=cos β .cos α
4th step, according to described total angle of attack∑With roll angle Φ, interpolation obtains interpolation polar coordinate Aerodynamic Coefficient, described interpolation
Polar coordinate Aerodynamic Coefficient includes axial force coefficientNormal force coefficient Cξ(αΣ, Φ), cross force Cη(αΣ, Φ), rolling
Moment coefficientYawing moment coefficient mξ(αΣ, Φ), pitching moment coefficient mη(αΣ,Φ);
5th step, obtains rectangular coordinate system Aerodynamic Coefficient, institute according to described interpolation polar coordinate Aerodynamic Coefficient and conversion formula
The rectangular coordinate system Aerodynamic Coefficient stated is axial force coefficientNormal force coefficientCornering ratioRolling moment coefficientYawing moment coefficientAnd pitching moment coefficientWherein
Conversion formula as follows:
With
6th step, the rectangular coordinate system Aerodynamic Coefficient using the 5th step to obtain carries out body kinetics resolving, obtains guided missile
Acceleration and angular acceleration;
7th step, uses the acceleration of guided missile and angular acceleration that the 6th step obtains, carries out kinesiology resolving, obtain guided missile next
The trajectory parameter flight speed in momentWith flight attitude angle;8th step, resolves according to the trajectory parameter of subsequent time and is led
Play the aerodynamic parameter of subsequent time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600360.4A CN106228014B (en) | 2016-07-27 | 2016-07-27 | A kind of acquisition methods of missile aerodynamic coefficient |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600360.4A CN106228014B (en) | 2016-07-27 | 2016-07-27 | A kind of acquisition methods of missile aerodynamic coefficient |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106228014A true CN106228014A (en) | 2016-12-14 |
CN106228014B CN106228014B (en) | 2019-01-29 |
Family
ID=57533464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610600360.4A Active CN106228014B (en) | 2016-07-27 | 2016-07-27 | A kind of acquisition methods of missile aerodynamic coefficient |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106228014B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845507A (en) * | 2018-06-15 | 2018-11-20 | 上海航天控制技术研究所 | The appraisal procedure of strong pneumatic nonlinearity characteristic flight control system rapidity index |
CN109492237A (en) * | 2017-09-12 | 2019-03-19 | 江西洪都航空工业集团有限责任公司 | A kind of Aerodynamic Coefficient preparation method |
CN109612676A (en) * | 2018-12-12 | 2019-04-12 | 北京空天技术研究所 | Aerodynamic parameter reverse calculation algorithms based on test flight data |
CN109827472A (en) * | 2018-12-12 | 2019-05-31 | 中国人民解放军陆军工程大学 | Method and device for processing downlink signal and electronic equipment |
CN112487731A (en) * | 2020-11-13 | 2021-03-12 | 北京电子工程总体研究所 | Air-defense missile trajectory modeling method based on pneumatic fixed connection coordinate system |
CN112762776A (en) * | 2021-01-22 | 2021-05-07 | 北京理工大学 | Rocket projectile tail end speed estimation method |
CN114237295A (en) * | 2021-12-20 | 2022-03-25 | 北京航空航天大学 | Unconventional flight control technology for high-agility air-to-air missile at large angle of attack |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103307938A (en) * | 2013-04-23 | 2013-09-18 | 北京电子工程总体研究所 | Method for acquiring aerodynamic parameters of spinning projectile |
CN105466660A (en) * | 2015-12-29 | 2016-04-06 | 清华大学 | Airplane aerodynamic parameter partition multi-step identification method adopting acceleration observer |
CN105629725A (en) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | Elastic motion modeling method of trailing edge rudder gliding aircraft |
CN105716826A (en) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | Reynolds number effect correcting method of zero lift-drag force coefficient |
-
2016
- 2016-07-27 CN CN201610600360.4A patent/CN106228014B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103307938A (en) * | 2013-04-23 | 2013-09-18 | 北京电子工程总体研究所 | Method for acquiring aerodynamic parameters of spinning projectile |
CN105629725A (en) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | Elastic motion modeling method of trailing edge rudder gliding aircraft |
CN105466660A (en) * | 2015-12-29 | 2016-04-06 | 清华大学 | Airplane aerodynamic parameter partition multi-step identification method adopting acceleration observer |
CN105716826A (en) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | Reynolds number effect correcting method of zero lift-drag force coefficient |
Non-Patent Citations (1)
Title |
---|
倪原等: "某飞行器的气动特性计算及仿真", 《西安工业大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492237A (en) * | 2017-09-12 | 2019-03-19 | 江西洪都航空工业集团有限责任公司 | A kind of Aerodynamic Coefficient preparation method |
CN108845507A (en) * | 2018-06-15 | 2018-11-20 | 上海航天控制技术研究所 | The appraisal procedure of strong pneumatic nonlinearity characteristic flight control system rapidity index |
CN109612676A (en) * | 2018-12-12 | 2019-04-12 | 北京空天技术研究所 | Aerodynamic parameter reverse calculation algorithms based on test flight data |
CN109827472A (en) * | 2018-12-12 | 2019-05-31 | 中国人民解放军陆军工程大学 | Method and device for processing downlink signal and electronic equipment |
CN112487731A (en) * | 2020-11-13 | 2021-03-12 | 北京电子工程总体研究所 | Air-defense missile trajectory modeling method based on pneumatic fixed connection coordinate system |
CN112487731B (en) * | 2020-11-13 | 2024-03-29 | 北京电子工程总体研究所 | Air defense missile trajectory modeling method based on pneumatic fixedly-linked coordinate system |
CN112762776A (en) * | 2021-01-22 | 2021-05-07 | 北京理工大学 | Rocket projectile tail end speed estimation method |
CN114237295A (en) * | 2021-12-20 | 2022-03-25 | 北京航空航天大学 | Unconventional flight control technology for high-agility air-to-air missile at large angle of attack |
Also Published As
Publication number | Publication date |
---|---|
CN106228014B (en) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228014A (en) | A kind of acquisition methods of missile aerodynamic coefficient | |
CN104792340B (en) | A kind of star sensor installation error matrix and navigation system star ground combined calibrating and the method for correction | |
CN107688295A (en) | A kind of quadrotor finite time self-adaptation control method based on fast terminal sliding formwork | |
CN109612676B (en) | Pneumatic parameter back calculation method based on flight test data | |
CN105509750B (en) | A kind of astronomy test the speed combined with terrestrial radio Mars capture section air navigation aid | |
CN104950901B (en) | Depopulated helicopter attitude error finite time convergence control nonlinear robust control method | |
CN109446582B (en) | High-precision order-reduction steady gliding dynamics modeling method considering earth rotation | |
CN111695193B (en) | Modeling method and system of globally relevant three-dimensional aerodynamic mathematical model | |
Baranowski | Equations of motion of a spin-stabilized projectile for flight stability testing | |
CN109492237A (en) | A kind of Aerodynamic Coefficient preparation method | |
CN103453907B (en) | Based on the planet approach section Navigation method of stratified atmosphere model | |
CN109059914B (en) | Projectile roll angle estimation method based on GPS and least square filtering | |
CN110320927A (en) | Flight control method and system of intelligent deformable aircraft | |
CN108225323B (en) | Method, medium, and apparatus for determining landing zone boundaries based on deviation influence direction combinations | |
Hui et al. | Dynamic attack zone of air-to-air missile after being launched in random wind field | |
CN106403934A (en) | Shell-borne geomagnetic attitude measurement and processing algorithm | |
CN106570242B (en) | The big dynamic pressure monoblock type radome fairing high-speed separation fluid structurecoupling emulation mode in low latitude | |
CN102566446A (en) | Method for establishing full-envelope mathematical model of unmanned helicopter based on linear model group | |
Desai et al. | Six-degree-of-freedom trajectory optimization using a two-timescale collocation architecture | |
Shen et al. | A 6DOF mathematical model of parachute in Mars EDL | |
Condaminet et al. | Identification of aerodynamic coefficients of a projectile and reconstruction of its trajectory from partial flight data | |
CN105987695A (en) | Interval quartering lagrangian method used for high-speed spinning projectile attitude algorithm | |
CN108414185B (en) | Zero error processing method for wind tunnel test data of symmetric aircraft | |
Harish et al. | Store separation dynamics using grid-free Euler solver | |
CN111649738B (en) | Method for calculating initial attitude of accelerometer under microgravity field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |