CN109190232A - A kind of aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure - Google Patents
A kind of aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure Download PDFInfo
- Publication number
- CN109190232A CN109190232A CN201810977464.6A CN201810977464A CN109190232A CN 109190232 A CN109190232 A CN 109190232A CN 201810977464 A CN201810977464 A CN 201810977464A CN 109190232 A CN109190232 A CN 109190232A
- Authority
- CN
- China
- Prior art keywords
- kinetic energy
- flow
- horizontal tail
- grid
- aircraft
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of aircraft horizontal tail area's kinetic energy rejections to calculate appraisal procedure, mainly by generate calculate grid and Kinetic Energy Calculation grid surface, numerical solution flow field, calculate horizontal tail area by grid surface come the energy of flow and calculate horizontal tail area kinetic energy rejection amount and etc. form, wherein the landing configuration model aircraft used includes the components such as fuselage, high lift device, engine nacelle, horizontal tail and vertical fin；Purpose is to the flow field data containing model aircraft being calculated based on numerical value, it obtains flowing through by post-processing interpolation computing method all come the energy of flow on the certain position Kinetic Energy Calculation grid surface of aircraft horizontal tail leading edge upstream, the difference for carrying out the energy of flow carrying out the energy of flow and being collected into Kinetic Energy Calculation grid surface in the flow field containing model aircraft that Kinetic Energy Calculation grid surface in the free Inlet-air Velocity Distribution without model aircraft is collected into is calculated under the conditions ofs identical incoming flow angle of attack etc., realizes that aircraft horizontal tail area's kinetic energy rejection calculates assessment.
Description
Technical field
The present invention relates to aerodynamic computation fields, calculate assessment side more particularly, to a kind of aircraft horizontal tail area's kinetic energy rejection
Method.
Background technique
In the big aircraft stall characteristic research process of landing configuration, fuselage, machine under the conditions of the different incoming flow angles of attack are needed to be grasped
One of an important factor for aeroperformance of the components such as the wing, horizontal tail, wherein horizontal tail aeroperformance is influence big aircraft stall characteristic.It is logical
Horizontal tail aeroperformance is often characterized with horizontal tail torque factor, and horizontal tail torque factor can be calculated by numerical value or test measurement hand
Section directly obtains.In fact, for the big aircraft of Yi Diaodi horizontal tail aerodynamic arrangement form, when aircraft is flown with large angle of attack
When, horizontal tail is easy the interference effect by the low energy separation vortex carried from wing tail, this, which will lead at horizontal tail, flows
It can reduce, to influence horizontal tail torque factor, and then influence rise and fall performance and the safety of aircraft.Aircraft horizontal tail upstream is come
The measurement for flowing tail, generallys use the measuring techniques such as PIV method and seven-hole probe both at home and abroad, the deficiency of this methods and techniques
Be in:
PIV method is also known as particle image velocimetry method, the main velocity profile information for measuring space flow field, but is unable to get incoming flow
The information such as the density of air and kinetic energy can not assess aircraft horizontal tail area kinetic energy rejection situation, and equipment is expensive, at
This height；
Seven-hole probe is mounted on tail rake velocity magnitude, direction and the stagnation pressure and static pressure that can measure local airflow near horizontal tail
Etc. information, can not but measure the information such as kinetic energy of incoming wake in a certain range of horizontal tail upstream, aircraft horizontal tail area can not be assessed
Kinetic energy rejection situation, and the design and difficulty of processing of seven-hole probe are larger；
It is the sum of all Incoming gas kinetic energy within the scope of certain space due to carrying out the energy of flow at aircraft horizontal tail, is tried currently without correlation
It tests equipment or numerical computation method and directly measures and calculate kinetic energy, therefore aircraft horizontal tail area kinetic energy change situation is not
It is concerned by people, and there are no carry out the quantitative method for calculating assessment to aircraft horizontal tail area kinetic energy rejection in engineering.
Summary of the invention
It is an object of the invention to be inserted by post-processing to the flow field data containing model aircraft being calculated based on numerical value
Value calculating method obtain flowing through on the certain position Kinetic Energy Calculation grid surface of aircraft horizontal tail leading edge upstream it is all come the energy of flow, in phase
The incoming flow being collected into calculating under the conditions of incoming flow angle of attack etc. without Kinetic Energy Calculation grid surface in the free Inlet-air Velocity Distribution of model aircraft
The difference for carrying out the energy of flow that Kinetic Energy Calculation grid surface is collected into kinetic energy and the flow field containing model aircraft realizes that aircraft horizontal tail area is dynamic
The loss of energy calculates assessment.
To achieve the goals above, the present invention adopts the following technical scheme:
A kind of aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure, comprising the following steps:
Step 1: the landing configuration model aircraft with high lift device and horizontal tail is generated using Grid Generation Software and contains model aircraft
Calculating grid, to without model aircraft same physical domain space generate be free of model aircraft calculating grid；
Step 2: using Grid Generation Software between aircraft components high lift device and horizontal tail, one is generated for calculating incoming flow
The two-dimensional grid face of kinetic energy, i.e. Kinetic Energy Calculation grid surface；
Step 3: under identical Mach number and Reynolds number, for certain incoming flow angle of attack condition, using Flow Field Calculation solver point
It is other that flow field numerical solution is carried out to the calculating grid containing model aircraft and the calculating grid without model aircraft, it obtains the incoming flow and attacks
Flow field data under corner condition containing model aircraft and the flow field data without model aircraft；
Step 4: using post-processing interpolation computing method to the net on Kinetic Energy Calculation grid surface in the flow field data containing model aircraft
Lattice point carries out the flow field data that interpolation obtains grid infinitesimal on grid surface, and the institute by Kinetic Energy Calculation grid surface is then calculated
Have and carry out the energy of flow, i.e., certain incoming flow angle of attack condition horizontal tail area that gets off the plane carrys out the energy of flow；
Step 5: using post-processing interpolation computing method on Kinetic Energy Calculation grid surface in the flow field data without model aircraft
Mesh point carries out the flow field data that interpolation obtains grid infinitesimal on grid surface, is then calculated through Kinetic Energy Calculation grid surface
It is all come the energy of flow；
Step 6: calculating carrys out the energy of flow without what Kinetic Energy Calculation grid surface in the free Inlet-air Velocity Distribution data of model aircraft was collected into
With the difference for carrying out the energy of flow that Kinetic Energy Calculation grid surface is collected into the flow field data containing model aircraft, certain incoming flow angle of attack item is obtained
Part is got off the plane horizontal tail area kinetic energy rejection amount；
Step 7: changing the incoming flow angle of attack, repeats step 3 to six, and the available difference incoming flow angle of attack conditions horizontal tail area that gets off the plane is dynamic
Loss of energy situation.
In the above-mentioned technical solutions, the Kinetic Energy Calculation grid surface is close between aircraft components high lift device and horizontal tail
The up-front position of horizontal tail.
In the above-mentioned technical solutions, it is disposed with several Grid dimensions in the grid surface, is arrived for collecting the angle of attack at 0 °
Stream information is come to all upstreams that horizontal tail has an impact within the scope of 30 °.
In the above-mentioned technical solutions, in the step 4 step 5, the post-processing interpolation computing method, which refers to, completes number
Arbitrary point position carries out interpolation and obtains the flow field data on the aspect in the space flow field that value resolves, and calculates acquisition and pass through kinetic energy meter
Calculate grid surface it is all come the energy of flow.
In the above-mentioned technical solutions, the post-processing interpolation computing method specifically includes the following steps:
Firstly for mesh point, that is, interpolation point on Kinetic Energy Calculation grid surface, positions and look into the calculating grid for completing Flow Field Calculation
Look for the lattice heart point for surrounding eight space lattice units of the interpolation point；
Then the flow field data of the interpolation point are calculated by Tri linear interpolation relationship again, have thus obtained Kinetic Energy Calculation net
The flow field data of mesh point on lattice face；
Secondly the face vector of each grid cell on Kinetic Energy Calculation grid surface is calculated, and with the stream of four angle points of the grid cell
Flow field data of the average value of field data as the grid cell；
Finally according to face vector sum flow field data calculate Kinetic Energy Calculation grid surface on each grid cell kinetic energy and sum to get
To by all come the energy of flow on the Kinetic Energy Calculation grid surface.
In the above-mentioned technical solutions, the Grid Generation Software be comparative maturity grid generate business software, as ICEM,
Pointwise etc.；
In the above-mentioned technical solutions, the grid that calculates is the structured grid that engineer application is more in Fluid Mechanics Computation field,
Including traditional muti-piece docking structure grid, the structured grid based on overlapping grid configuration and the Structure Network based on splicing operator form
Lattice；
In the above-mentioned technical solutions, the Flow Field Calculation solver can be business hydrodynamics method popular in the market
Software or fluid calculation open source program, such as FLUENT, CFL3D, are also possible to the fluid calculation software of independent research.
In conclusion by adopting the above-described technical solution, the beneficial effects of the present invention are:
It is different from traditional measurement method and technological means, get off the plane horizontal tail area by calculating certain incoming flow angle of attack condition come the energy of flow,
The investigation parameter object different from physical quantitys such as horizontal tail torques is provided for the horizontal tail aeroperformance research of the big aircraft of landing configuration；
It is got off the plane horizontal tail area kinetic energy rejection situation by the way that certain incoming flow angle of attack condition is quantitatively evaluated, for research Elevator stall characteristic peace
Tail aeroperformance and its optimization design provide reference；
By investigating aircraft horizontal tail area's kinetic energy rejection with the variation tendency of the incoming flow angle of attack, ginseng is provided to study big aircraft stall characteristic
It examines；
It, can be from pitching moment before horizontal tail area kinetic energy rejection interpretation aircraft stall by analyzing horizontal tail area kinetic energy rejection situation
There is the reason of inflection point, reference can be provided to promote the maneuvering performance of aircraft.
Detailed description of the invention
Examples of the present invention will be described by way of reference to the accompanying drawings, in which:
Fig. 1 is that aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure flow chart；
Fig. 2 is the calculating grid schematic diagram containing model aircraft；
Fig. 3 is free from the calculating grid schematic diagram of model aircraft；
Fig. 4 is to calculate space grid cell schematic diagram in grid；
Fig. 5 is grid cell schematic diagram in Kinetic Energy Calculation grid surface
Fig. 6 is the aircraft horizontal tail area incoming flow Kinetic Energy Calculation assessment result of embodiment one to three；
Fig. 7 is that aircraft horizontal tail area's kinetic energy rejection of embodiment one to three calculates assessment result；
Fig. 8 is that aircraft horizontal tail area's kinetic energy rejection of example IV calculates assessment result and relevant aeroplane pitching moment coefficient compares；
In figure: 1 is fuselage, and 2 be the calculating grid containing model aircraft, and 3 are free from the calculating grid of model aircraft, and 4 be kinetic energy meter
Grid surface is calculated, 5 be vertical fin, and 6 be horizontal tail, and 7 be engine nacelle, and 8 be high lift device, and 9 be to calculate space grid cell in grid,
10 be to calculate mesh point in grid, and 11 be lattice heart point, and 12 be grid cell in Kinetic Energy Calculation grid surface, and 13 be Kinetic Energy Calculation grid
Mesh point in face.
Specific embodiment
All features disclosed in this specification or disclosed all methods or in the process the step of, in addition to mutually exclusive
Feature and/or step other than, can combine in any way.
Aircraft horizontal tail area's kinetic energy rejection of the invention calculates appraisal procedure as shown in Figure 1, aircraft horizontal tail area kinetic energy rejection meter
It calculates appraisal procedure and mainly passes through net by generating calculating grid and Kinetic Energy Calculation grid surface, numerical solution flow field, calculating horizontal tail area
The energy of flow and calculating horizontal tail area kinetic energy rejection amount composition are come in lattice face, need to use Grid Generation Software, Flow Field Calculation
Solver, the calculating grid containing model aircraft, the flow field data containing model aircraft, is free of aircraft mould at landing configuration model aircraft
Calculating grid, the free Inlet-air Velocity Distribution data without model aircraft, Kinetic Energy Calculation grid surface and the post-processing interpolation calculation of type
Method etc., wherein landing configuration model aircraft includes the components such as fuselage, high lift device, engine nacelle, horizontal tail and vertical fin.
The present invention is to the flow field data containing model aircraft being calculated based on numerical value, by post-processing interpolation computing method
Obtain flowing through on the certain position Kinetic Energy Calculation grid surface of aircraft horizontal tail leading edge upstream it is all come the energy of flow, in the identical incoming flow angle of attack
Calculating carrys out the energy of flow and containing winged without what Kinetic Energy Calculation grid surface in the free Inlet-air Velocity Distribution of model aircraft was collected into Deng under the conditions of
What Kinetic Energy Calculation grid surface was collected into carrys out the difference of the energy of flow in the flow field of machine model, realizes that aircraft horizontal tail area's kinetic energy rejection calculates
Assessment:
Certain incoming flow angle of attack condition horizontal tail area that gets off the plane is calculated and carrys out the energy of flow, increases the investigation parameter of research horizontal tail aeroperformance
Object；
Certain incoming flow angle of attack condition is quantitatively evaluated to get off the plane horizontal tail area kinetic energy rejection situation, for research Elevator stall characteristic and pneumatic property
Reference can be provided；
Aircraft horizontal tail area's kinetic energy rejection is obtained with the variation tendency of the incoming flow angle of attack, provides reference to study big aircraft stall characteristic；
There is the reason of inflection point from pitching moment before horizontal tail area kinetic energy rejection interpretation aircraft stall.
Embodiment one
1) Pointwise Grid Generation Software is respectively adopted to certain landing configuration people containing components such as high lift device 8 and horizontal tails 6
Machine model generates the calculating grid 2 containing model aircraft, as shown in Fig. 2, the calculating grid is the structured grid for being overlapped grid configuration,
The calculating grid 3 for being free of model aircraft is generated to the same physical domain space without model aircraft, as shown in figure 3, the calculating net
Lattice are traditional muti-piece docking structure grid；
2) up-front close to horizontal tail using Pointwise Grid Generation Software between aircraft components high lift device 8 and horizontal tail 6
Position generates one for calculating the two-dimensional grid face for carrying out the energy of flow, i.e. Kinetic Energy Calculation grid surface 4, which is disposed with 101
× 161 mesh points, the gross area are about 38.3368m2, can be collected into all upstream incoming flows letter having an impact to horizontal tail
Breath；
3) in Mach 2 ship 0.2, the Reynolds number based on average air power chord length is 19,800,000, the incoming flow angle of attack be 0 ° it is identical come
Under the conditions of stream, to the calculating grid 2 containing model aircraft and be free of respectively using the fluid calculation software PMB3D of independent research
The calculating grid 3 of model aircraft carries out flow field numerical solution, obtain after restraining under the conditions of 0 ° of incoming flow angle of attack containing model aircraft
Flow field data and free Inlet-air Velocity Distribution data without model aircraft, which includes speed, density, pressure etc., is deposited
Storage is calculating grid cell centers, i.e. lattice heart point 11, as shown in Figure 4；
4) passed through on Kinetic Energy Calculation grid surface 4 using post-processing interpolation computing method in the flow field data containing model aircraft
It is all come the energy of flow calculated, specific calculating process are as follows: firstly for the mesh point 13 on Kinetic Energy Calculation grid surface 4 i.e. insert
Be worth point, as shown in figure 5, complete Flow Field Calculation calculating grid in (i.e. comprising model aircraft calculating grid 2 or do not include
The calculating grid 3 of model aircraft) positioning searching surround the interpolation point 8 space lattice units 9 lattice heart point 11, such as Fig. 4 institute
Show；Acquisition is calculated since the related streams field data of the lattice heart point 11 of these space lattice units 9 has passed through Flow Field Calculation solver
(i.e. then the flow field data of the interpolation point are calculated, thus in step 3) in technical solution by Tri linear interpolation relationship again
The flow field data of the mesh point 13 on Kinetic Energy Calculation grid surface 4 are obtained；Secondly each grid on Kinetic Energy Calculation grid surface 4 is calculated
The face vector of unit 12, and with four angle points of the grid cell 12 (i.e. mesh point 13, that is, above-mentioned interpolation point)
Flow field data flow field data of the average value as the grid cell 12；Finally kinetic energy is calculated according to face vector sum flow field data
It calculates on grid surface 4 kinetic energy of each grid cell 12 and sums and get off the plane horizontal tail area incoming flow to get to 0 ° of incoming flow angle of attack condition
Kinetic energy is about 7300695.31J/s；
5) similarly, using post-processing interpolation computing method to without model aircraft free Inlet-air Velocity Distribution data in pass through kinetic energy
Calculate grid surface 4 on it is all come the energy of flow calculated, specific calculating process is consistent with step 4), obtains 0 ° of incoming flow angle of attack item
It by Kinetic Energy Calculation grid surface 4 is about 7389592.53J/s come the energy of flow under part；
6) calculate without model aircraft free Inlet-air Velocity Distribution data in Kinetic Energy Calculation grid surface 4 be collected into come the energy of flow with
To the difference for carrying out the energy of flow that Kinetic Energy Calculation grid surface 4 in the flow field data containing model aircraft is collected into, 0 ° of incoming flow angle of attack is obtained
Condition is got off the plane horizontal tail area kinetic energy rejection amount are as follows: 7389592.53J/s -7300695.31J/s=88897.22J/s.
Embodiment two
Embodiment two and embodiment one are almost the same, and the main distinction is, in step 3), the incoming flow angle of attack is 10 °, keep Mach
Several and Reynolds number is constant, continues the fluid calculation software PMB3D using independent research respectively to the calculating containing model aircraft
Grid 2 and calculating grid 3 without model aircraft carry out flow field numerical solution, obtain after restraining under the conditions of 10 ° of incoming flow angles of attack
Flow field data containing model aircraft and the free Inlet-air Velocity Distribution data without model aircraft；In step 4), after continuing to use
Manage interpolation computing method in the flow field data containing model aircraft by Kinetic Energy Calculation grid surface 4 it is all come the energy of flow progress
It calculates, obtains 10 ° of incoming flow angle of attack conditions and get off the plane that carry out the energy of flow be about 6951836.67J/s in horizontal tail area；In step 5), continue
Using post-processing interpolation computing method in the free Inlet-air Velocity Distribution data without model aircraft by Kinetic Energy Calculation grid surface 4
It is all come the energy of flow calculated, obtain be about come the energy of flow by Kinetic Energy Calculation grid surface 4 under the conditions of 10 ° of incoming flow angles of attack
7278507.28 J/s；In step 6), Kinetic Energy Calculation grid surface 4 in the free Inlet-air Velocity Distribution data without model aircraft is calculated
The difference for carrying out the energy of flow carried out the energy of flow Yu Kinetic Energy Calculation grid surface 4 in the flow field data containing model aircraft is collected into being collected into
Value obtains 10 ° of incoming flow angle of attack conditions and gets off the plane horizontal tail area kinetic energy rejection amount are as follows: 7278507.28J/s -6951836.67J/s
= 326670.61J/s。
Embodiment three
Embodiment three and embodiment one and embodiment two are almost the same, and the main distinction is, change the incoming flow angle of attack respectively, repeat
Step 3) -6), it has obtained different incoming flow angle of attack conditions and has got off the plane horizontal tail area kinetic energy rejection situation.In the present embodiment, incoming flow is attacked
Angle takes 4 °, 8 °, 12 °, 14 °, 16 °, 17 °, 18 °, 19 °, 20 °, 21 ° etc. respectively, has obtained different incoming flow angle of attack conditions and has got off the plane
Horizontal tail area comes the energy of flow and kinetic energy rejection situation.The calculating assessment result of embodiment three and the meter of embodiment one and embodiment two
Assessment result is calculated to show together in figure 6 and figure 7.Fig. 7 shows that aircraft horizontal tail area's kinetic energy rejection becomes with the variation of the incoming flow angle of attack
Gesture.
Example IV
Example IV and the embodiment of front three are almost the same, and the main distinction is, in step 1), containing 8 peace of high lift device
Certain landing configuration model aircraft shape of the components such as tail 6 is changed, not with the model aircraft shape in embodiment one to three
Unanimously, therefore using Pointwise Grid Generation Software the calculating grid 2 containing model aircraft that generates and without model aircraft
It calculates grid 3 and the calculating grid in embodiment one to three is inconsistent；In step 3), using the hydrodynamics meter of independent research
Calculate software PMB3D carry out flow field numerical solution obtain the flow field data containing model aircraft after restraining under the conditions of 0 ° of incoming flow angle of attack and
Free Inlet-air Velocity Distribution data without model aircraft are inconsistent；In step 4), using post-processing interpolation computing method to containing winged
It is calculated in the flow field data of machine model by all on Kinetic Energy Calculation grid surface 4 come the energy of flow, obtains 0 ° of incoming flow angle of attack
Condition is got off the plane, and to carry out the energy of flow be about 7274893.86J/s in horizontal tail area；In step 5), post-processing interpolation computing method pair is used
It is calculated in free Inlet-air Velocity Distribution data without model aircraft by all on Kinetic Energy Calculation grid surface 4 come the energy of flow,
It obtains under the conditions of 0 ° of incoming flow angle of attack through Kinetic Energy Calculation grid surface 4 being about 7389439.53J/s come the energy of flow；In step 6)
In, calculate that Kinetic Energy Calculation grid surface 4 in the free Inlet-air Velocity Distribution data without model aircraft is collected into come the energy of flow with to containing
The difference for carrying out the energy of flow that Kinetic Energy Calculation grid surface 4 is collected into the flow field data of model aircraft obtains 0 ° of incoming flow angle of attack condition
Horizontal tail of getting off the plane area kinetic energy rejection amount are as follows: 7389439.53J/s -7274893.86J/s=114545.67J/s；
Different incoming flow angle of attack conditions are got off the plane horizontal tail area kinetic energy rejection situation in order to obtain, are changed the incoming flow angle of attack respectively, are repeated
Step 3) -6).In the present embodiment, the incoming flow angle of attack take respectively 0 °, 4 °, 8 °, 10 °, 12 °, 14 °, 16 °, 18 °, 19 °, 20 °,
21 °, 22 °, 24 ° etc., different incoming flow angle of attack conditions are calculated and have got off the plane horizontal tail area kinetic energy rejection amount.The aircraft of example IV
It is as shown in Figure 8 that horizontal tail area kinetic energy rejection calculates assessment result.Fig. 8 is also shown based on example IV model aircraft shape numerical value
Obtained aircraft pitch torque coefficient with the incoming flow angle of attack variation tendency.As can be seen from Figure 8, after 16 ° of angles of attack, aircraft
Pitching moment starts to face upward, and since stalling angle is about 20 °, inflection point occurs in advance in aircraft pitch torque, and from aircraft horizontal tail
Area's kinetic energy rejection is seen with the variation tendency of the incoming flow angle of attack, after 16 ° of angles of attack, aircraft horizontal tail area kinetic energy rejection start it is increased compared with
It is more, according to the analysis, aircraft horizontal tail area kinetic energy occur greater loss may be before leading to aircraft stall pitching moment there is inflection point
One of reason.
The invention is not limited to specific embodiments above-mentioned.The present invention, which expands to, any in the present specification to be disclosed
New feature or any new combination, and disclose any new method or process the step of or any new combination.
Claims (5)
1. a kind of aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure, it is characterised in that the following steps are included:
Step 1: the landing configuration model aircraft with high lift device and horizontal tail is generated using Grid Generation Software and contains model aircraft
Calculating grid, to without model aircraft same physical domain space generate be free of model aircraft calculating grid；
Step 2: using Grid Generation Software between aircraft components high lift device and horizontal tail, one is generated for calculating incoming flow
The two-dimensional grid face of kinetic energy, i.e. Kinetic Energy Calculation grid surface；
Step 3: under identical Mach number and Reynolds number, for certain incoming flow angle of attack condition, using Flow Field Calculation solver point
It is other that flow field numerical solution is carried out to the calculating grid containing model aircraft and the calculating grid without model aircraft, it obtains the incoming flow and attacks
Flow field data under corner condition containing model aircraft and the flow field data without model aircraft；
Step 4: using post-processing interpolation computing method to the net on Kinetic Energy Calculation grid surface in the flow field data containing model aircraft
Lattice point carries out the flow field data that interpolation obtains grid infinitesimal on grid surface, and the institute by Kinetic Energy Calculation grid surface is then calculated
Have and carry out the energy of flow, i.e., certain incoming flow angle of attack condition horizontal tail area that gets off the plane carrys out the energy of flow；
Step 5: using post-processing interpolation computing method on Kinetic Energy Calculation grid surface in the flow field data without model aircraft
Mesh point carries out the flow field data that interpolation obtains grid infinitesimal on grid surface, is then calculated through Kinetic Energy Calculation grid surface
It is all come the energy of flow；
Step 6: calculating carrys out the energy of flow without what Kinetic Energy Calculation grid surface in the free Inlet-air Velocity Distribution data of model aircraft was collected into
With the difference for carrying out the energy of flow that Kinetic Energy Calculation grid surface is collected into the flow field data containing model aircraft, certain incoming flow angle of attack item is obtained
Part is got off the plane horizontal tail area kinetic energy rejection amount；
Step 7: changing the incoming flow angle of attack, repeats step 3 to six, and the available difference incoming flow angle of attack conditions horizontal tail area that gets off the plane is dynamic
Loss of energy situation.
2. a kind of aircraft horizontal tail area's kinetic energy rejection according to claim 1 calculates appraisal procedure, it is characterised in that described dynamic
Grid surface can be calculated between aircraft components high lift device and horizontal tail close to the up-front position of horizontal tail.
3. a kind of aircraft horizontal tail area's kinetic energy rejection according to claim 2 calculates appraisal procedure, it is characterised in that the net
It is disposed with several Grid dimensions in lattice face, horizontal tail is had an impact for collecting the angle of attack within the scope of 0 ° to 30 ° all
Trip carrys out stream information.
4. a kind of aircraft horizontal tail area's kinetic energy rejection according to claim 1 calculates appraisal procedure, it is characterised in that the step
In rapid four steps five, the post-processing interpolation computing method refers to that completing arbitrary point position in the space flow field that numerical value resolves carries out
Interpolation obtains the flow field data on the aspect, and calculates acquisition by all come the energy of flow of Kinetic Energy Calculation grid surface.
5. a kind of aircraft horizontal tail area's kinetic energy rejection according to claim 4 calculates appraisal procedure, it is characterised in that after described
Handle interpolation computing method specifically includes the following steps:
Firstly for mesh point, that is, interpolation point on Kinetic Energy Calculation grid surface, positions and look into the calculating grid for completing Flow Field Calculation
Look for the lattice heart point for surrounding eight space lattice units of the interpolation point；
Then the flow field data of the interpolation point are calculated by Tri linear interpolation relationship again, have thus obtained Kinetic Energy Calculation net
The flow field data of mesh point on lattice face；
Secondly the face vector of each grid cell on Kinetic Energy Calculation grid surface is calculated, and with the stream of four angle points of the grid cell
Flow field data of the average value of field data as the grid cell；
Finally according to face vector sum flow field data calculate Kinetic Energy Calculation grid surface on each grid cell kinetic energy and sum to get
To by all come the energy of flow on the Kinetic Energy Calculation grid surface.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810977464.6A CN109190232B (en) | 2018-08-27 | 2018-08-27 | Method for calculating and evaluating kinetic energy loss of plane stern area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810977464.6A CN109190232B (en) | 2018-08-27 | 2018-08-27 | Method for calculating and evaluating kinetic energy loss of plane stern area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109190232A true CN109190232A (en) | 2019-01-11 |
| CN109190232B CN109190232B (en) | 2022-03-25 |
Family
ID=64919877
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810977464.6A Active CN109190232B (en) | 2018-08-27 | 2018-08-27 | Method for calculating and evaluating kinetic energy loss of plane stern area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109190232B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113420379A (en) * | 2021-06-29 | 2021-09-21 | 西北工业大学 | Method for extracting surface-average pressure distribution from CFL3D calculation result |
| CN113868779A (en) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | Propeller aircraft blade flying-off safety assessment method based on multidisciplinary coupling simulation |
| CN116227244A (en) * | 2023-05-09 | 2023-06-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for designing relative position of plane taking Mach number change and optimal drag reduction effect into consideration |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1856654A (en) * | 2004-05-08 | 2006-11-01 | 全力 | Vortex producing method and apparatus and aircraft therewith |
| US20100133377A1 (en) * | 2008-10-17 | 2010-06-03 | Airbus Operations (Societe Par Actions Simplifiee) | Airplane with flat rear fuselage said queue-de-morue empennage |
| CN105183975A (en) * | 2015-09-01 | 2015-12-23 | 中国空气动力研究与发展中心高速空气动力研究所 | Trailing vortex field transfer based numerical simulation method for multi-aircraft flight formation |
| CN106158230A (en) * | 2015-04-05 | 2016-11-23 | 田景华 | The firearms rocket jet engine steam turbine internal combustion engine of additional magnetic power apparatus and aircraft |
-
2018
- 2018-08-27 CN CN201810977464.6A patent/CN109190232B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1856654A (en) * | 2004-05-08 | 2006-11-01 | 全力 | Vortex producing method and apparatus and aircraft therewith |
| US20100133377A1 (en) * | 2008-10-17 | 2010-06-03 | Airbus Operations (Societe Par Actions Simplifiee) | Airplane with flat rear fuselage said queue-de-morue empennage |
| CN106158230A (en) * | 2015-04-05 | 2016-11-23 | 田景华 | The firearms rocket jet engine steam turbine internal combustion engine of additional magnetic power apparatus and aircraft |
| CN105183975A (en) * | 2015-09-01 | 2015-12-23 | 中国空气动力研究与发展中心高速空气动力研究所 | Trailing vortex field transfer based numerical simulation method for multi-aircraft flight formation |
Non-Patent Citations (7)
| Title |
|---|
| LIU PENGYUAN等: "《WTA model study of air defense missile system based on particle algorithm》", 《PROCEEDINGS OF 2014 IEEE CHINESE GUIDANCE, NAVIGATION AND CONTROL CONFERENCE》 * |
| XIANG ZHENG等: "《Robust finite-time attitude tracking control for nonlinear quadrotor with uncertainties and delays》", 《2017 IEEE INTERNATIONAL CONFERENCE ON ROBOTICS AND BIOMIMETICS (ROBIO)》 * |
| 崔兴达等: "《超声速气流中低频燃烧振荡实验研究》", 《高等学校工程热物理第二十届全国学术会议论文集——燃烧学专辑》 * |
| 肖中云等: "《船尾形状对旋成体马格努斯效应的影响》", 《航空学报》 * |
| 胡孟权等: "《自适应模糊—滑模控制在重构飞行控制中的应用》", 《航空学报》 * |
| 邱名等: "《超声速压气机叶栅研究进展》", 《2017年（第三届）中国航空科学技术大会论文集（上册）》 * |
| 黄江涛等: "《大型民用飞机气动外形典型综合设计方法》", 《航空学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113420379A (en) * | 2021-06-29 | 2021-09-21 | 西北工业大学 | Method for extracting surface-average pressure distribution from CFL3D calculation result |
| CN113868779A (en) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | Propeller aircraft blade flying-off safety assessment method based on multidisciplinary coupling simulation |
| CN113868779B (en) * | 2021-12-02 | 2022-02-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Propeller aircraft blade flying-off safety assessment method based on multidisciplinary coupling simulation |
| CN116227244A (en) * | 2023-05-09 | 2023-06-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for designing relative position of plane taking Mach number change and optimal drag reduction effect into consideration |
| CN116227244B (en) * | 2023-05-09 | 2023-07-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for designing relative position of plane taking Mach number change and optimal drag reduction effect into consideration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109190232B (en) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Visbal et al. | Three-dimensional vortex formation on a heaving low-aspect-ratio wing: Computations and experiments | |
| Hartloper et al. | On the competition between leading-edge and tip-vortex growth for a pitching plate | |
| CN109190232A (en) | A kind of aircraft horizontal tail area's kinetic energy rejection calculates appraisal procedure | |
| Chitsaz et al. | An experimental study of the aerodynamics of micro corrugated wings at low Reynolds number | |
| Ma et al. | Time-resolved topology of turbulent boundary layer separation over the trailing edge of an airfoil | |
| Hahn et al. | Implicit large-eddy simulation of swept-wing flow using high-resolution methods | |
| Shetty et al. | CFD Simulation of Flow Past MAV Wings. | |
| CN115879231A (en) | Method and device for acquiring flow transition information of micro aircraft | |
| Chitsaz et al. | Numerical and experimental analysis of three-dimensional microcorrugated wing in gliding flight | |
| Khalil et al. | Computational analyses of aerodynamic characteristics of naca653218airfoil | |
| Garmann et al. | Implicit les computations for a rapidly pitching plate | |
| Jia et al. | Force and flow field measurement system for tandem flapping wings | |
| Baik et al. | Unsteady force generation and vortex dynamics of pitching and plunging flat plates at low Reynolds number | |
| Newsome et al. | Numerical simulation of vortical flow over an elliptical-body missile at high angles of attack | |
| Ramesh | Numerical and experimental investigation of the effectof geometry modification on the aerodynamic characteristics of a NACA 64 (2)-415wing | |
| Gordnier et al. | Aeroelastic Simulations of an Aspect Ratio Two Flexible Membrane Wing | |
| Hanke et al. | Numerical Simulation of DLR-F11 High Lift Configuration from HiLiftPW-2 using STAR-CCM+ | |
| Uthra et al. | Comparative Investigation of Laminar Separation Bubble on a Wing at Low Reynolds Number | |
| Kashitani et al. | Study on Busemann Biplane Airfoil in Low-Speed Smoke Wind Tunnel | |
| Kojima et al. | Computational study of flow characteristics of thick and thin airfoil with implicit large-eddy simulation at low reynolds number | |
| Wibowo et al. | Vortex dynamic investigation of wing slotted gap of saab jas gripen c-like fighter | |
| Smith | Large-Eddy-Simulation-Based Wind Tunnel Assessments | |
| Shin et al. | Parametric study on efficient formation flying for a blended-wing uav | |
| Verhaagen | Investigation of the vortex flow over a 76/60-deg double delta wing at 20 deg incidence | |
| Breuer et al. | High-performance computing for the investigation of the flow past an airfoil with trailing-edge stall |
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 |