CN107357976B - Method for calculating dynamic derivative of aircraft - Google Patents
Method for calculating dynamic derivative of aircraft Download PDFInfo
- Publication number
- CN107357976B CN107357976B CN201710499236.8A CN201710499236A CN107357976B CN 107357976 B CN107357976 B CN 107357976B CN 201710499236 A CN201710499236 A CN 201710499236A CN 107357976 B CN107357976 B CN 107357976B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- aerodynamic force
- residual error
- flow field
- boundary condition
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention relates to the technical field of fluid dynamics and discloses a method for calculating a dynamic derivative of an aircraft. Determining static flow field boundary conditions b for flight states of an aircraft1The governing equation of the flow field is converged to 0, i.e. R ═ Ax + b1When 0, the flow field is for the boundary condition b1Solution vector x of1Wherein A is a matrix of the governing equation, and when the residual error does not converge to absolute zero, it is set to converge to the allowable residual error R1=Ax1+b1At this time, the aerodynamic force isSolving equationsWherein F is aerodynamic force, R is residual error, Lambda is adjoint matrix, x is solution vector, and solution vector x is obtained1Conditional corresponding adjoint matrix lambda1(ii) a Increasing a rotation angular speed delta omega around a reference center by the aircraft, generating a normal motion speed on the surface of the aircraft, and calculating a residual error increment delta R generated by the normal speed; the aerodynamic force increment generated by the normal speed is delta F ═ Lambda1Δ R, then the dynamic derivative of the aerodynamic force with respect to the angular velocity can be obtained asThe solution scheme of the dynamic derivative is simple and accurate.
Description
Technical Field
The invention relates to the technical field of fluid dynamics, in particular to a method for calculating a dynamic derivative of an aircraft.
Background
The calculation of aircraft aerodynamic forces by the CFD method is usually solved in the form of
Ax+b=0 (1)
The matrix A is a matrix representing a control equation, b represents a boundary condition, and x is a vector to be solved.
The method for solving the flow field control equation (1) is mainly a time correlation method, and an initial value x is set for a vector x to be solved0Over a time-varying course
Finally, the solution of the formula (1) is obtained. Wherein R is the residual, and the process of solving is the process of gradually converging the residual to 0. The existing methods comprise tests and calculations, and the methods for solving the dynamic derivative are all that the aircraft is subjected to periodic forced vibration (rotation around the gravity center when the derivative of angular velocity is solved, translation when the derivative of the change rate of the angle of attack/sideslip angle is solved), the time history of the aerodynamic force change is obtained through the calculations or tests, and the corresponding dynamic derivative is obtained according to the phase relation between the aerodynamic force change and the change rate of the angular velocity/angle of attack.
The following is an example of the derivative of the pitch moment with respect to the pitch angle rate and the rate of change of the angle of attack. Setting the aircraft (model) at equilibrium position alpha0The nearby part is in pitching sinusoidal vibration, and the state equation is
Where A is the angular amplitude, T is the period of vibration, and T is the time. Then there is
Where Δ α is the angle of attack deviation amplitude, ω is the angular velocity,is the rate of change of angle of attack. The aerodynamic force at small amplitude rotation/oscillation satisfies the following relation:
F=F(α0)+ΔαFα+ω(Fω+Fα) (4)
after the change course of the aerodynamic force F is obtained, the integral is carried out on a complete period
The sum of the angular velocity and the derivative of the angle of attack change rate of the aerodynamic force can be obtained
These two terms cannot be separately obtained with forced rotation only. If the aircraft (model) is subjected to up-and-down sinusoidal vibration, the state equation is
Wherein Y is the amplitude. Then there is
Where is the v flight velocity. Only influencing the aerodynamic force during this movementThe magnitude of the deviation, ω is the angular velocity,is the rate of change of angle of attack. Formula (A), (B) and4) the following steps are changed:
after the change course of the aerodynamic force F is obtained, the integral is carried out on a complete period
The derivative of the aerodynamic force to the change rate of the angle of attack can be obtained
When the dynamic derivative is calculated by adopting a numerical calculation method, the unsteady aerodynamic force change course of the forced oscillation process needs to be calculated, and then the dynamic derivative term is distinguished by using the phase component of the aerodynamic force change. The calculation in the unsteady process has high precision requirement, long time consumption and great resource consumption, and if the amplitude of vibration is overlarge, the calculation result contains excessive high-order terms, so that the precision is poor; if the amplitude is too small, the amount of variation may be too small, and the calculation error may introduce a large deviation. Wind tunnel tests are subject to test accuracy, as are conflicts in amplitude selection.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: in view of the existing problems, a method for calculating the dynamic derivative of the aircraft is provided.
The technical scheme adopted by the invention is as follows: a method for calculating the dynamic derivative of an aircraft specifically comprises the following processes:
step 1, determining static flow field boundary condition b aiming at flight state of aircraft1The governing equation of the flow field is converged to 0, i.e. R ═ Ax + b1When 0, the flow field is for the boundary condition b1Solution vector x of1Wherein A is a matrix of the governing equation, and when the residual error does not converge to absolute zero, it is set to converge to the allowable residual error R1=Ax1+b1At this time, it is pneumaticForce is
Step 2, solving equationWherein F is aerodynamic force, R is residual error, Lambda is adjoint matrix, x is solution vector, and solution vector x is obtained1Conditional corresponding adjoint matrix lambda1;
Step 3, increasing a rotation angular velocity delta omega around the reference center by the aircraft, generating a normal motion velocity on the surface of the aircraft, and calculating a residual error increment delta R generated by the normal velocity;
and 4, increasing the aerodynamic force generated by the normal speed to delta F ═ Lambda1Δ R, then the dynamic derivative of the aerodynamic force with respect to the angular velocity can be obtained as
Further, the calculation method of the dynamic derivative of the line former further comprises the following processes: (a) a tiny rotation angular velocity delta omega is added around the reference center by the aircraft, the normal motion velocity is generated on the surface of the aircraft, and a corresponding boundary condition b is obtained2(ii) a (b) With x1Is the initial field, corresponds to b2As a boundary condition, there is a residual R2=Ax1+b2(ii) a (c) Boundary condition becomes b2The aerodynamic force at the initial moment of the flow field is still the aerodynamic force of the static flow fieldIn the process that the residual R converges to 0, satisfy Then(d) When the residual error converges to 0, it can be obtainedIs corresponding to the boundary condition b2The aerodynamic force of (a); (f) similarly, if the solution in step 1 is continuously calculated and converged to the residual error equal to 0, the solution can also be obtained(g) The amount of change Δ F in aerodynamic force is represented by a adjoint matrix as F2-F1=Λ1Δ R; wherein Δ R ═ R2-R1=(Ax1+b2)-(Ax1+b1)=b2-bx; (h) obtaining the dynamic derivative of the aerodynamic force to the angular velocity as
Compared with the prior art, the beneficial effects of adopting the technical scheme are as follows: the technical scheme of the invention can obtain the dynamic derivative of the aerodynamic force to the angular velocity only by requiring a solution field and an adjoint matrix thereof without calculating unsteady motion. According to the principle of the method, the obtained dynamic derivative is only related to residual variation delta R caused by the change of the boundary condition, and the residual existing in the initial flow field has no influence, so that the degree of residual convergence in the calculation of the initial flow field can be properly relaxed. The method greatly reduces resource consumption, obtains a strict dynamic derivative, does not contain the influence of other factors such as high-order components, incidence angle change rate and the like, improves the accuracy, and does not have the limitation of vibration amplitude and frequency.
Drawings
Fig. 1 is a schematic flow chart of the method for calculating the dynamic derivative of the aircraft according to the invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
As shown in fig. 1, a method for calculating a dynamic derivative of an aircraft specifically includes the following processes:
step 1, determining static flow field boundary condition b aiming at flight state of aircraft1The governing equation of the flow field is converged to 0, i.e. R ═Ax+b1When 0, the flow field is for the boundary condition b1Solution vector x of1Wherein A is a matrix of the governing equation, and when the residual error does not converge to absolute zero, it is set to converge to the allowable residual error R1=Ax1+b1At this time, the aerodynamic force is
Step 2, solving equationWherein F is aerodynamic force, R is residual error, Lambda is adjoint matrix, x is solution vector, and solution vector x is obtained1Conditional corresponding adjoint matrix lambda1;
Step 3, increasing a rotation angular velocity delta omega around the reference center by the aircraft, generating a normal motion velocity on the surface of the aircraft, and calculating a residual error increment delta R generated by the normal velocity;
and 4, increasing the aerodynamic force generated by the normal speed to delta F ═ Lambda1Δ R, then the dynamic derivative of the aerodynamic force with respect to the angular velocity can be obtained as
When the aircraft increases a rotation angular speed delta omega around the reference center, a corresponding boundary condition b is introduced2. The calculation method of the dynamic derivative of the line former further comprises the following processes: (a) a tiny rotation angular velocity delta omega is added around the reference center by the aircraft, the normal motion velocity is generated on the surface of the aircraft, and a corresponding boundary condition b is obtained2(ii) a (b) With x1Is the initial field, corresponds to b2As a boundary condition, there is a residual R2=Ax1+b2(ii) a (c) Boundary condition becomes b2The aerodynamic force at the initial moment of the flow field is still the aerodynamic force of the static flow fieldIn the process that the residual R converges to 0, satisfyThen(d) When the residual error converges to 0, it can be obtainedIs corresponding to the boundary condition b2The aerodynamic force of (a); (f) similarly, if the solution in step 1 is continuously calculated and converged to the residual error equal to 0, the solution can also be obtained(g) The amount of change Δ F in aerodynamic force is represented by a adjoint matrix as F2-F1=Λ1Δ R; wherein Δ R ═ R2-R1=(Ax1+b2)-(Ax1+b1)=b2-b1(ii) a (h) Obtaining the dynamic derivative of the aerodynamic force to the angular velocity asThe invention is not limited to the foregoing embodiments. The invention extends to any novel feature or any novel combination of features disclosed in this specification and any novel method or process steps or any novel combination of features disclosed. Those skilled in the art to which the invention pertains will appreciate that insubstantial changes or modifications can be made without departing from the spirit of the invention as defined by the appended claims.
Claims (2)
1. A method for calculating a dynamic derivative of an aircraft is characterized by specifically comprising the following processes:
step 1, determining static flow field boundary condition b aiming at flight state of aircraft1The governing equation of the flow field is converged to 0, i.e. R ═ Ax + b1When 0, the flow field is for the boundary condition b1Solution vector x of1Wherein A is a matrix of the governing equation, and when the residual error does not converge to absolute zero, it is set to converge to the allowable residual error R1=Ax1+b1At this time, the aerodynamic force is
Step 2, solving equationWherein F is aerodynamic force, R is residual error, Lambda is adjoint matrix, x is solution vector, and solution vector x is obtained1Conditional corresponding adjoint matrix lambda1;
Step 3, increasing a rotation angular velocity delta omega around the reference center by the aircraft, generating a normal motion velocity on the surface of the aircraft, and calculating a residual error increment delta R generated by the normal velocity;
2. The method of calculating the dynamic derivatives of an aircraft according to claim 1, characterized in that it further comprises the following processes: (a) a tiny rotation angular velocity delta omega is added around the reference center by the aircraft, the normal motion velocity is generated on the surface of the aircraft, and a corresponding boundary condition b is obtained2(ii) a (b) With x1Is the initial field, corresponds to b2As a boundary condition, there is a residual R2=Ax1+b2(ii) a (c) Boundary condition becomes b2The aerodynamic force at the initial moment of the flow field is still the aerodynamic force of the static flow fieldIn the process that the residual R converges to 0, satisfy Then(d) When the residual error converges to 0, it can be obtainedIs corresponding to the boundary condition b2The aerodynamic force of (a); (f) similarly, if the solution in step 1 is continuously calculated and converged to the residual error equal to 0, the solution can also be obtained(g) The amount of change Δ F in aerodynamic force is represented by a adjoint matrix as F2-F1=Λ1Δ R; wherein Δ R ═ R2-R1=(Ax1+b2)-(Ax1+b1)=b2-b1(ii) a (h) Obtaining the dynamic derivative of the aerodynamic force to the angular velocity as
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710499236.8A CN107357976B (en) | 2017-06-27 | 2017-06-27 | Method for calculating dynamic derivative of aircraft |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710499236.8A CN107357976B (en) | 2017-06-27 | 2017-06-27 | Method for calculating dynamic derivative of aircraft |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107357976A CN107357976A (en) | 2017-11-17 |
CN107357976B true CN107357976B (en) | 2020-12-11 |
Family
ID=60272577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710499236.8A Active CN107357976B (en) | 2017-06-27 | 2017-06-27 | Method for calculating dynamic derivative of aircraft |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107357976B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063391B (en) * | 2018-09-30 | 2020-10-09 | 上海机电工程研究所 | Dynamic derivative calculation detection method under rotation condition and dynamic derivative wind tunnel test method |
CN112393876B (en) * | 2019-08-16 | 2022-04-12 | 北京空天技术研究所 | Dynamic pneumatic derivative prediction method suitable for internal and external flow integrated appearance |
CN110619160B (en) * | 2019-09-02 | 2022-12-02 | 四川腾盾科技有限公司 | Implicit solution method based on accompanying residual sorting |
CN110674607B (en) * | 2019-09-02 | 2022-11-18 | 四川腾盾科技有限公司 | Implicit solution method based on residual magnitude ordering |
CN114088330B (en) * | 2021-11-19 | 2023-03-17 | 中国航空工业集团公司哈尔滨空气动力研究所 | Longitudinal dynamic derivative low-speed wind tunnel continuous measurement test method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246515A (en) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | Digital ship model planar motion mechanism experimental method based on CFD software |
CN106650095A (en) * | 2016-12-21 | 2017-05-10 | 中国航天空气动力技术研究院 | Method for correcting unmanned aerial vehicle control matrix based on wind tunnel test data and CFD calculation |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8010326B2 (en) * | 2005-12-28 | 2011-08-30 | Caterpillar Inc. | Method and apparatus for automated grid formation in multi-cell system dynamics models |
US9223754B2 (en) * | 2012-06-29 | 2015-12-29 | Dassault Systèmes, S.A. | Co-simulation procedures using full derivatives of output variables |
CN103400035B (en) * | 2013-07-29 | 2016-12-28 | 中国空气动力研究与发展中心高速空气动力研究所 | A kind of method of high credibility fast prediction aircraft rolling dynamic derivative |
CN106055791B (en) * | 2016-05-31 | 2019-06-04 | 西北工业大学 | Aircraft overall situation Aerodynamic optimization method based on Predictor-Correcting Algorithm |
-
2017
- 2017-06-27 CN CN201710499236.8A patent/CN107357976B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246515A (en) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | Digital ship model planar motion mechanism experimental method based on CFD software |
CN106650095A (en) * | 2016-12-21 | 2017-05-10 | 中国航天空气动力技术研究院 | Method for correcting unmanned aerial vehicle control matrix based on wind tunnel test data and CFD calculation |
Non-Patent Citations (3)
Title |
---|
CFD simulations of turbulent nonreacting and reacting flows for forrocket engine applications;Martin J. Seidla et.al;《The Journal of Supercritical Fluids》;20161120;第121卷;63–77 * |
流体力学有限元分析中的边界条件处理;赵永凯等;《东北重型机械学院学报》;19970930;第21卷(第3期);189-195 * |
高超声速飞行器飞行性能评估***设计与实现;廉波;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑(月刊)》;20160315(第 03 期);C031-388 * |
Also Published As
Publication number | Publication date |
---|---|
CN107357976A (en) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107357976B (en) | Method for calculating dynamic derivative of aircraft | |
CN106078742B (en) | A kind of vibration control method for being directed to the flexible mechanical arm with output constraint | |
AU2015213424B2 (en) | Methods and apparatus to control aircraft horizontal stabilizers | |
CN108256264B (en) | Pneumatic servo elastic stability prediction method based on ground frequency response test | |
CN111006835B (en) | Rocket projectile pitching moment coefficient and pressure center coefficient correction method and storage medium | |
KR102021498B1 (en) | Design method of attitude control system for flight vehicle and computer program | |
CN111258216B (en) | Sliding mode repetitive controller suitable for four-rotor aircraft | |
CN106896722B (en) | The hypersonic vehicle composite control method of adoption status feedback and neural network | |
CN109063391B (en) | Dynamic derivative calculation detection method under rotation condition and dynamic derivative wind tunnel test method | |
CN112180965A (en) | High-precision overload control method | |
Anton et al. | Stability derivatives for a delta-wing X-31 aircraft validated using wind tunnel test data | |
CN108469734A (en) | Consider the motor servo system Auto-disturbance-rejection Control of state constraint | |
CN107976296A (en) | A kind of aerodynamic characteristics of vehicle on-line identification method based on backtracking adaptive algorithm | |
CN107256290B (en) | Method for calculating influence difference of boundary change on flow field based on residual perturbation method | |
CN111427267A (en) | High-speed aircraft attack angle tracking method adopting force and moment adaptive estimation | |
Zou et al. | Three-dimensional ventilated supercavity on a maneuvering trajectory | |
Voß et al. | A ROM based flutter prediction process and its validation with a new reference model | |
CN112904898B (en) | Method and system for evaluating unsteady pneumatic response characteristic of rotary rocket | |
CN111015661A (en) | Active vibration control method and system for flexible load of robot | |
KR102114051B1 (en) | Method of non-linear control for aircraft considering center of gravity movement | |
Farcy et al. | Sensitivity of spin parameters to uncertainties of the aerodynamic model of a light aircraft | |
Embacher et al. | Coupled CFD-simulation of a helicopter in free-flight trim | |
Beyers | Interpretation of experimental high-alpha aerodynamics-Implications for flight prediction | |
Kuzmina et al. | Dynamic response of an airplane elastic structure in transonic flow | |
Kasumov | Modeling of helicopter behavior dynamics under the controlling action of the stabilization and control system |
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 |