CN106228014A - A kind of acquisition methods of missile aerodynamic coefficient - Google Patents

A kind of acquisition methods of missile aerodynamic coefficient Download PDF

Info

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
Application number
CN201610600360.4A
Other languages
Chinese (zh)
Other versions
CN106228014B (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.)
Jiangxi Hongdu Aviation Industry Group Co Ltd
Original Assignee
Jiangxi Hongdu Aviation Industry Group Co Ltd
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 Jiangxi Hongdu Aviation Industry Group Co Ltd filed Critical Jiangxi Hongdu Aviation Industry Group Co Ltd
Priority to CN201610600360.4A priority Critical patent/CN106228014B/en
Publication of CN106228014A publication Critical patent/CN106228014A/en
Application granted granted Critical
Publication of CN106228014B publication Critical patent/CN106228014B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/02Wind tunnels
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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 formulaWith roll angle Φ;Afterwards according to described total angle of attackWith 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

A kind of acquisition methods of missile aerodynamic coefficient
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 α
t a n Φ = t a n β sin α ;
4th step, according to described total angle of attackWith 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 α
t a n Φ = t a n β sin α ;
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 belowWith Roll angle Φ, described formula is:
CosA=cos β .cos α
t a n Φ t a n β sin α ;
4th step, according to described total angle of attackWith 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.
CN201610600360.4A 2016-07-27 2016-07-27 A kind of acquisition methods of missile aerodynamic coefficient Active CN106228014B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* 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
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